# Slope Stability Analysis Using Limit Equilibrium Method in Nonlinear Criterion

Authors: Hang Lin; Wenwen Zhong; Wei Xiong; Wenyu Tang
Journal: The Scientific World Journal
 (2014)

---

## Abstract

In slope stability analysis, the limit equilibrium method is usually used to calculate the safety factor of slope based on Mohr-Coulomb criterion. However, Mohr-Coulomb criterion is restricted to the description of rock mass. To overcome its shortcomings, this paper combined Hoek-Brown criterion and limit equilibrium method and proposed an equation for calculating the safety factor of slope with limit equilibrium method in Hoek-Brown criterion through equivalent cohesive strength and the friction angle. Moreover, this paper investigates the impact of Hoek-Brown parameters on the safety factor of slope, which reveals that there is linear relation between equivalent cohesive strength and weakening factorD. However, there are nonlinear relations between equivalent cohesive strength and Geological Strength Index (GSI), the uniaxial compressive strength of intact rock σci, and the parameter of intact rock mi. There is nonlinear relation between the friction angle and all Hoek-Brown parameters. With the increase of D, the safety factor of slope F decreases linearly; with the increase of GSI, F increases nonlinearly; when σci is relatively small, the relation between F and σci is nonlinear, but when σci is relatively large, the relation is linear; with the increase of mi, F decreases first and then increases.


---

## Body

## 1. Introduction

The analysis of slope stability is an important aspect of research in disaster reduction engineering. In the analysis of slope stability, the primary task is to calculate the safety factor of the slope [1–5]. Because the analysis of slope stability is statically indeterminate, the factor of safety cannot be obtained directly from the static equilibrium condition. In order to avoid the complex relation between stress and strain and the difficulty in transforming the statically indeterminate problem into the statically determinate problem, reasonable approximations for the analysis of slope stability are required. Thus the problem can be solved in statically determinate form with limit equilibrium method (LEM) [6–8]. This method is at the cost of the rigor but has little impact on the precision. It is applied widely in engineering projects with many distinctive advantages; for example, it simplifies the analytical calculation and shortens the computing time. The fundamental feature of limit equilibrium method is only considering the static equilibrium condition and Mohr-Coulomb criterion; that is, the solution to the problem is the analysis of the equilibrium of forces at the moment of the failure of soil mass. However, there is a restriction of Mohr-Coulomb criterion to the description of rock mass; for example, it fails to explain the impact of high stress areas but only reflects the feature of linear failure in rock mass [9, 10]. To overcome the above shortcomings, based on large amount of rock tests and related materials, Hoek and Brown proposed Hoek-Brown criterion which can reflect intrinsic characteristics including the rock strength [10–12], the number of rock joints, and the impact of stress state on rock strength as well as the feature of nonlinear failure. The criterion matches characteristics of deformation and failure well [13–16] which won the high reputation and wide application in geological engineering. Zhang et al. [17] proposed a three-dimensional (3D) version of the generalized Hoek-Brown strength criterion. Halakatevakis and Sofianos [9] investigated Hoek-Brown criterion analytically through an extended plane of weakness theory and validated it with experimental evidence on physical specimens. Besides, some scholars applied numerical analytical methods to slope stability on the basis of Hoek-Brown criterion. While numerical methods consume long time to calculate the safety factor of the slope [18], for instance, when taking the CPU time required by each of the models on the same computer into consideration, limit equilibrium method ran considerably faster, taking about 2 minutes per run compared with 30 minutes and 7 hours using the coarse (about 2000 elements) and fine grids (about 7000 elements), respectively, by numerical calculation method. Though the traditional method for the safety factor, limit equilibrium method, owns the advantage of simplicity, it fails to combine Hoek-Brown criterion. Taking all the above factors into account, this paper proposed an equation for the safety factor with limit equilibrium method in Hoek-Brown criterion through theoretical deduction. Moreover, this paper investigates the impact of Hoek-Brown parameters on the safety factor of slope.


## 2. Definition of Safety Factor

So far, there are three types of safety factors [2, 3, 19]: (1) the strength reserve factor of safety obtained through lowering the strength of rock and soil mass; (2) the overload reserve factor of safety obtained through increasing the exterior load; (3) the driving force overloading reserve factor, a design value of landslide thrust calculated through amplifying the driving force along the slope while keeping the corresponding resisting force constant. The limit equilibrium method mainly adopts the concept of the strength reserve factor.In 1955, Bishop proposed the modified Bishop method [20], an effective and simple method to analyze circular slip surfaces. In this method, the safety factor Fs1 is defined as the potential sliding along a failure surface. After reducing the shear strength index by Fs1 times, that is, c0/Fs1 and tanϕ0/Fs1, the rock and soil mass along the failure surface are in a state of limiting equilibrium; that is,
(1)τ=ccr+σtan⁡ϕcr,
where τ is the shear strength and σ is the normal stress at failure surfaces; ccr=c0/Fs1; tanϕcr=tanϕ0/Fs1, where c0 and ϕ0 are the cohesive strength and the friction angle, respectively, in original state; ccr and ϕcr are the cohesive strength and the friction angle, respectively, in the limit equilibrium state of the slope.The mentioned method that takes the reserve of strength index to define the factor of safety has clear physical meanings. The factor of safety is the ratio of resisting to deriving forces acting at the slope along the failure surface. This type of safety factor is only reducing the resisting forces while keeping deriving forces constant. Comparatively, the definition of safety factor with strength reserve matches the practical engineering projects more perfectly, because frequently they are external factors that contribute to failures of slopes and the decrease of strength of rock and soil mass. Thus, this definition has been widely accepted by engineers around the world with the examination of practical projects for many years.


## 3. Hoek-Brown Equivalent Limit Equilibrium Method

The modified Bishop method is applied widely in practical engineering projects [21–23]. As a limit equilibrium method, this method is not confined to the shape of failure surface and possesses relatively high precision. The result with this method is very close to the one with the strict slice method. Therefore, this paper adopts the modified Bishop method to calculate the factor of safety. The modified Bishop method sets three hypothetical conditions: (1) the actual tangential force at the failure surface is the part that the shear strength works; (2) the vertical shear strength between slices is ignored; (3) the yield of soil mass matches the Mohr-Coulomb criterion. The calculation model for safety factor of slope is illustrated in Figure 1. Through deduction, the factor of safety can be expressed as
(2)Fs=∑mai[cibi+(Wi-uibi)tan⁡ϕi]∑Wisin⁡αi,
where
(3)mai=1cos⁡⁡αi+(tan⁡ϕisin⁡αi/Fs),
where Fs is the safety factor of slope and ci and ϕi are the effective cohesive strength and the friction angle of soil layer. If the slope is homogeneous, both the cohesive strength and the friction angle are constant; ui is the pore water pressure loading at the failure surface of the soil slice i; bi is the width of soil slice; αi is the inclination of the bottom of soil slice from the horizontal; Wi is the gravity of soil slice i, Wi=γ1bih1i+γ2bih2i+⋯+γjbihji+⋯+γkbihki; k is the number of soil layer where the failure surface is; γj is the bulk density of the soil at the number j layer; hji is the height of slice i in the number j soil layer.Figure 1
Theoretical model for deducing the safety factor of slope.On the basis of (1), if the relation between Hoek-Brown strength parameters and Mohr-Coulomb strength parameters can be set up, substituting parameters gives the safety factor of slope through the modified Bishop method in Hoek-Brown criterion. Therefore, the first step is to deduce the relation between Hoek-Brown criterion and Mohr-Coulomb criterion.Hoek-Brown criterion is expressed as [10–12]
(4)σ1′=σ3′+(mbσciσ3′+sσci2)a,
where σ1′ is the effective maximum principle stress when the rock mass is destroyed; σ3′ is the effective minimum principle stress loaded at the rock mass; σci is the uniaxial compressive strength of intact rock; mb is the constant of rock mass and is related to mi of intact rock; s and a are the coefficient depending on features of rock mass. All of these parameters can be described with a function of the Geological Strength Index (GSI).Specific functions are as follows [12, 17]:
(5)mb=mi·exp⁡[(GSI-100)(28-14D)],s=exp⁡[(GSI-100)(9-3D)],a=12+e-GSI/15-e-20/36,
where D is the weakening factor of rock which is related to the excavation way and the disturbing degree. The value of D is in the range of 0~1, where 0 indicates the undisturbed state.Equation (4) is differentiated to obtain
(6)dσ1′dσ3′=1+ambσci(mbσ3′σci+sσci2)a-1.Based on (6), the shear stress and the normal stress are expressed with the principle stress [10]:
(7)τ=(σ1′-σ3′)dσ1′/dσ3′(dσ1′/dσ3′+1),σn′=σ1′+σ3′2-(σ1′-σ3′)2·(dσ1′/dσ3′)-1(dσ1′/dσ3′)+1.Since Mohr-Coulomb failure criterion is still widely used in geotechnical engineering analysis, it is necessary to determine equivalent angles of friction and cohesive strengths for each rock mass and stress range. In Mohr-Coulomb criterion,tan⁡ϕ=dτ/dσn, c=τ-σntan⁡ϕ. By fitting an average linear relationship to the curve generated by solving (4) for a range of minor principal stress values, as illustrated in Figure 2, the equivalent angle of internal friction and the cohesive strength can be obtained as [12]
(8)ϕ=arcsin(fbfc2fa+fbfc),(9)c=σcifc[s(1+2a)+(1-a)mbσ3n]fa1+fbfc/fa,
where fa=(1+a)(2+a); fb=6amb; fc=(s+mbσ3n′)a-1; σ3n′=σ3max⁡′/σci; σ3max⁡′ is the upper limit of confining stress over which the relationship between the Hoek-Brown and the Mohr-Coulomb criteria is considered, has to be determined for each individual case [12]. In the studies of slopes using Bishop’s circular failure analysis for a wide range of slope geometries and rock mass properties, σ3max⁡′=0.72(σcm/γH)-0.91; γ is the bulk density of rock mass; H is the height of slope; σcm indicates the strength of rock mass:
(10)σcm=σci[mb+4s-a(mb-8s)(mb/4+s)a-1]2(1+a)(2+a).Figure 2
The relation between the maximum principle stress and the minimum principle stress in Hoek-Brown criterion and Mohr-Coulomb criterion.The fitting process involves balancing the areas above and below the Mohr-Coulomb plot. Transferring Figure2 into the strength envelope diagram Figure 3 shows that the Hoek-Brown envelope diagram is a curve, while the Mohr-Coulomb envelope diagram is a straight line [12]. Figure 3 is divided into three parts and marked as regions 1, 2, and 3, respectively. When the normal stress is in region 1 or region 3, the equivalent Mohr-Coulomb strength parameter will overestimate the shear strength of rock mass. When most of normal stress at the bottom of failure surface concentrates in region 1 or region 3, the factor of safety calculated with equivalent Mohr-Coulomb strength parameters may slightly be insecure. Through the above deduction, substituting (8) and (9) into (2) gives the safety factor of slope calculated by the modified Bishop method in Hoek-Brown criterion.Figure 3
The Hoek-Brown strength envelope diagram and equivalent Mohr-Coulomb strength envelope diagram.


## 4. Parametric Analysis

The height of a slope is 20 m and the inclination is 1 : 0.5. For the convenience of discussion, the homogenous model of slope is set up. Basic calculating parameters are bulk densityγ=25.0 kN/m3, mi=17, D=0.5, GSI=50, and σci=95.0 MPa.Changing parameters, respectively,mi=2~32, D=0~1, GSI=5~95, and σci=5~185 MPa, in (8) and (9), gives the corresponding Mohr-Coulomb parameters, that is, the cohesive strength c and the friction angle ϕ. Figures 4, 5, 6, and 7 present relations between Hoek-Brown parameters, that is, D, GSI, σci and mi, and the cohesive strength c as well as the friction angle ϕ. These figures illustrate that the equivalent friction angle and the cohesive strength of rock material decrease gradually with the increase of weakening factor D. Furthermore, the relation between the cohesive strength and D is linear but the relation between the friction angle and D is nonlinear. The relation between GSI and the cohesive strength is nonlinear as well as the friction angle. Moreover, with the increase of GSI, the cohesive strength increases, but the friction angle increases first and then decreases. With the increase of σci, the cohesive strength presents a tendency of linear increase, while the friction angle presents a tendency of nonlinear increase. What is more, the slope of relation curve of the friction angle and σci decreases gradually as σci increases, which reveals that the friction angle is very sensitive to the impact of σci when the value of σci is relatively small. As mi increases, the cohesive strength shows the nonlinear decreasing characteristic, while the friction angle shows the nonlinear increasing characteristic. That is to say, mi has the opposite impact on the cohesive strength and the friction angle.Figure 4
The relation betweenD and the cohesive strength/angle of internal friction.Figure 5
The relation between GSI and the cohesive strength/angle of internal friction.Figure 6
The relation betweenσci and the cohesive strength/angle of internal friction.Figure 7
The relation betweenmi and the cohesive strength/angle of internal friction.Figure8 displays the relation between the safety factor of slope F and other parameters. It can be concluded from the figure that as for the D, the linear characteristic of F decreases with the increase of D and the relation matches the function F=-4.918·(D-1.920). As for GSI, the nonlinear characteristic of F increases with the increase of GSI and the relation matches the exponential equation F=0.186·1.068GSI. When GSI>60, the amplitude of increase of F is obviously large which is much higher than F with the change of other parameters. As for σci, when σci is relatively small, the relation between F and σci is nonlinear, while it turns out to be linear when σci is relatively large and it can be fitted with the function F=15.275-13.156·0.995σci with high correlation. As for mi, F decreases first and then increases with the increase of mi which matches the function F=5.222·exp⁡(-mi/2.698)+6.533+0.0267·mi well.The relation between the safety factor of slopeF and all Hoek-Brown parameters.
(a)
The relation betweenD and F
(b)
The relation between GSI andF
(c)
The relation betweenσci and F
(d)
The relation betweenmi and F


## 5. Conclusions

(1)
Relations between Hoek-Brown parameters, that is,D, GSI, σci and mi, and the cohesive strength as well as the friction angle are set up based on theoretical calculation and the corresponding equivalent cohesive strength and the friction angle with Hoek-Brown parameters.
(2)
The equation for calculating the safety factor of slope with limit equilibrium method in Hoek-Brown criterion is set up through equivalent cohesive strength and the friction angle. This method is applicable with practical value.
(3)
There is linear relation between equivalent cohesive strength and weakening factorD. However, there are nonlinear relations between equivalent cohesive strength and Geological Strength Index (GSI), the uniaxial compressive strength of intact rock σci, and the parameter of intact rock mi. There is nonlinear relation between the friction angle and all Hoek-Brown parameters.
(4)
With the increase ofD, the safety factor of slope F decreases linearly; with the increase of GSI, F increases nonlinearly; when σci is relatively small, the relation between F and σci is nonlinear, but when σci is relatively large, the relation is linear; with the increase of mi, F decreases first and then increases.



---

*Source: 206062-2014-07-23.xml*