Research Article  Open Access
Yu Zhao, ZhiYi Tong, Qing Lü, "Slope Stability Analysis Using SliceWise Factor of Safety", Mathematical Problems in Engineering, vol. 2014, Article ID 712145, 6 pages, 2014. https://doi.org/10.1155/2014/712145
Slope Stability Analysis Using SliceWise Factor of Safety
Abstract
The concept of slicewise factor of safety is introduced to investigate the state of both the whole slope and each slice. The assumption that the interslice force ratio is the same between any two slices is made and the eccentric moment of slice weight is also taken into account. Then four variables equations are formulated based on the equilibrium of forces and moment and the assumption of interslice forces, and then the slicewise factor of safety along the slip surface can be obtained. The active and passive sections of the slope can be determined based on the distribution of factor of safety. The factor of safety of the whole slope is also defined as the ratio of the sum of antisliding force to the sum of sliding force on the slip surface. Two examples with different slip surface shapes are analysed to demonstrate the usage of the proposed method. The slicewise factor of safety enables us to determine the sliding mechanism and pattern of a slope. The reliability is verified by comparing the overall factor of safety with that calculated by conventional methods.
1. Introduction
Among various methods currently available for slope stability analysis, conventional methods based on the concept of limit equilibrium have been most widely used in engineering practice. Though finite element analysis is becoming an attractive alternative [1–3], the limit equilibrium technique will probably continue to play an important role in the further slope engineering due to its simplicity and its ease of use.
The key procedure and main purpose of slope stability analysis using the limit equilibrium technique are the calculation of the factor of safety. Given a predefined slip surface, the factor of safety is determined with these methods from the equilibrium of force and/or momentum of the mass contained between the slip surface and the free ground surface. Some of the proposed methods are only for circular slip surfaces [4, 5], while more recent ones are for any shape of slip surfaces [6–9]. The list is not exhaustive. In addition, these methods are also different in the equilibrium conditions that they satisfy. The ordinary methods of slices [4] satisfy only the moment equilibrium; Bishop’s modified method [5] satisfies the moment equilibrium and vertical force equilibrium. Morgenstern and Price’s method [6], Janbu’s generalized procedure of slice [10], Sarma’s method [11], and slope stability charts [9] satisfy all conditions of equilibrium and differ from each other in the assumptions about interslice forces. Recently, researches focus on methods to find the critical slip surfaces instead of the limit equilibrium technique itself. Sarma and Tan [12] used the stress acceptability criterion to locate the critical slip surface; Li et al. [13] employed a realcoded genetic algorithm to develop a search approach for locating the noncircular critical slip surface.
However, we realized that the limit equilibrium technique for slope analyses still has some aspects to improve. The equation in each method mentioned above is an implicit equation of the factor of safety and there is a connotative assumption that the factor of safety of each slice is equal to the factor of safety of the slope. The assumption is just a simple treatment of the factor of safety in order to solve the implicit equation more conveniently. Wright et al. [14], Tavenas et al. [15], and others noted that the factor of safety varies from place to place along the slip surface. Furthermore, the slicewise factor of safety is helpful for engineers to determine the active or passive section of a slope and the countermeasures for slope stabilization. Misleading measures based on incorrectly determined passive sections resulted in disastrous consequences in many cases [16, 17].
In this paper, the slicewise factor of safety is first proposed (the definition of factor of safety is based on the conception of shear strength reduction), and the influence of its distribution along the slip surface on the stability of the whole slope is discussed. For each slice, four variables, including the horizontal and vertical interslice forces, the normal force on the slip surface, and the slicewise factor of safety, are used to establish the slicewise equations, and then the group of equations are solved iteratively using the boundary conditions at the first and the last slice.
2. Basic Assumption
2.1. Definition of SliceWise Factor of Safety
For the limit equilibrium method, the mass of a slope between the slip surface and the free ground surface is evenly divided into slices as shown in Figure 1(a). The main forces acted on the th slice (Figure 1(b)) are the weight (), the horizontal earthquake force (), the horizontal and vertical load (), the water pressure , the interslice forces (), the normal force (), and the sliding force () of the slip surface.
The slicewise factor of safety is defined as the ratio of the antisliding force to the sliding force, which is the same as the conception of shear strength reduction [1], where is the shear strength parameters of the slip surface and the width of the slice.
In the conventional strength reduction method, the factor of safety () means that the shear strength of all the slices’ bottom surface should be reduced by times, which is obviously not the realistic state of the slope. By defining the slicewise factor of safety, an individual reduction coefficient is given to each slice and all the slip surfaces can reach the limit equilibrium state simultaneously.
2.2. Discussion on the Height of Thrust Line
The rigorous Janbu’s method [10] suggests that the height of the thrust line is between 1/3 and 1/2 of the height of the slice’s profile. Actually the height of the thrust line is dependent on the property of the slope mass. For a loose deposit slope, the distribution of horizontal soil pressure would be triangular and thus the height of the thrust line is onethird of the height of the slice’s profile. For a hard soil or block rock slope, the distribution of horizontal soil would be a parallelogram, and the height of the thrust line is half of the height of the slice’s profile. Therefore, the determination of the thrust line height should be based on the geological investigation of the slope and differ from slope to slope.
2.3. Equilibrium of Forces and Boundary Condition
From the horizontal and vertical equilibrium of the forces
According to (1), the sliding force in Figure 1 can be replaced by . Substituting into (2a) and (2b) yields
At the boundaries, the horizontal force acting on the left side of the first slice and the force on the right side of the last slice are usually assigned certain values, 0 in this case.
2.4. Assumption of Interslice Forces
The major difference among all the limit equilibrium methods is the assumption about the interslice force. No matter when a slope is sliding or stable, the interaction pattern between slices is generally similar to one another and seldom influenced by the states of the slices. Therefore, we assume the interslice force ratio (the horizontal force to the vertical force) equal to the same value : where is the height of the profile of the slice. The zero horizontal force indicates there is no interaction between the two slices, and thus must be zero too.
As a sign of global mobilization of the slope, the value of the interslice force ratio must be between 0 and 1. The larger the is, the better the antisliding forces are exerted. When the factor of safety of a slope is less than 1 and the slope is unstable, the antisliding force between the slices is exerted to the maximum at that time and the value of should be larger. In contrast, when the slope is stable, the antisliding force between the slices will not be exerted sufficiently and the value of should be much smaller.
2.5. Discussion on Balance of Moment
From the equilibrium of the moment at the midpoint of slice bottom, we can obtain
An assumption connoted in (5) is that the weight , the vertical load , and the normal force all pass through the midpoint of slice bottom. According to (5), the boundary condition at the first slice will be when the horizontal load is zero regardless of the earthquake force . Consequently, a simple calculation using (4) suggests that all the interslice forces are zero, which is not consistent with actual stress state of the slope.
In order to avoid this erroneous situation, we recommended that the eccentric moment of slice weight is to be taken into account. For a parallelogramshape slice, the mass center passes through the midpoint of slice bottom, whereas the eccentric moment is no longer zero in the case of a trapezoidal slice. A straightforward way is to divide a trapezoid into a triangle and a parallelogram, and the eccentric moment is only caused by the weight of triangle part. Hence, we can define the eccentric moment as where is the height of the upper profile of the slice and is the height of the lower profile of the slice. We neglect the negative eccentric moment and assign zero to it when is lower than .
After introducing the eccentric moment , the balance of moment is expressed as follows:
3. Solution of the SliceWise and the Overall Factor of Safety
3.1. Solution of the SliceWise Factor of Safety
From (3a), (3b), (4), and (6) of the four variables , , , , we can express as the following formula: where
We can attain by substituting the boundary condition at the first slice into (8). The expression can be obtained successively, and then the value of can be solved with the boundary condition . The interslice forces and can be obtained using the value of and (4) and (8). After substituting into (3a) and (3b), we can solve the simultaneous equations for the normal force and the slicewise factor of safety .
After obtaining all the slicewise factors of safety, the sliding mechanism of the slope also can be discussed accordingly.
3.2. Solution of the Overall Factor of Safety
The factor of safety of the whole slope is also defined by the ratio of the sum of antisliding force to the sum of sliding force on the slip surface: where
The method for slope stability analysis proposed in this paper can satisfy all the equilibrium conditions of forces and moment. Compared with the conventional limit equilibrium methods, the method can obtain not only the factor of the whole slope, but also the slicewise factor of safety which can enable us to determine the mechanism and sliding pattern of a slope. In addition, the interslice force ratio also has the special physical meaning that reflects the exertion degree of the antisliding forces of slope mass.
4. Two Examples
In this section, we select two classical examples from the tests published by ACADS [18] to verify the method proposed in this paper. One slope has circular slip surface and the other has polylines slip surface.
4.1. Slope with Circular Slip Surface
The calculation example of circular slip surface is the first one of the ten tests published by ACADS EX1(a). Figure 2 shows the profile of the slope and the parameter of the soil mass.
According to the method proposed, firstly we should solve the value of by the equation .
The slope is composed of soil mass, and thus the height of thrust line of every slice is equal to 1/3 of the height of the slice’s profile, . There are two solutions for equation : and (2) . For the first one, the factor of safety of the whole slope is 0.968; for the second, the corresponding factor of safety is 0.996. The two factors of safety that are less than 1.0 indicate the slope is unstable. Hence, the antisliding forces between the slices should be exerted to the maximum and we must select the bigger and the corresponding factor of safety of the slope as 0.996.
Figure 3 plots the distribution of slicewise factor of safety along the slip surface when .
Figure 3 shows that the slicewise factors of safety of the slices from 22 to 30 at the lower region of the slope mass are all less than 1.0, which indicates that the region is the active section and will slide before the other slices. Hence, the sliding mechanism of the slope can be judged as traction sliding.
Table 1 shows the factor of safety calculated by a number of methods. The difference of the factor of safety between the proposed method and the standard answer introduced by Donald is only 0.4%, which verifies the reliability of the method.

4.2. Slope with Polylines Slip Surface
The calculation example of polylines slip surface is the first one of the ten tests published by ACADS EX3(b). Figure 4 and Table 2 show the profile of the slope and the parameters of the soil mass, respectively.

The slope is composed of soil mass, and thus the height of thrust line of every slice is 1/3 of the height of the slice’s profile, . There is only one solution for equation , , and the factor of safety of the whole slope 1.247. The nearzero indicates that the slope mass is integrated tightly, and the vertical forces are small enough to be ignored. In this case, the result is very close to the result given by Bishop’s method which ignores the influence of vertical forces.
Figure 5 shows that the slicewise factors of safety of slices from 11 to 29 in the middle of the slope are smaller than 1 and distribute reposefully. Hence, this region is probably the active section, and the region from slices 1 to 10 and the region from slices 30 to 32 are the passive sections of the slope.
Table 3 shows the calculation results of different methods. From Table 3 we can conclude that the error of the factors of slope safety between the proposed method and Morgenstern and Price’s methods is only 1.4%. With an accuracy of about ±6%, the factor of safety calculated using methods that satisfy all conditions of equilibrium can be considered to be the correct answer [19]. Because the error is small enough to be ignored, the proposed method is quite reliable.

4.3. Influence of the Height of Thrust Line
In order to evaluate the influence of thrust line height on the factor of safety, we use 1/2 of the height of the slice’s profile as the thrust line height and calculate the two examples again.
For the example of circular slip surface, the single solution of equation is , and the factor of safety is 0.972. For the example of polylines slip surface, the single solution of equation is , and the corresponding factor of safety is 1.263. Though the height of thrust line has little effect on the factor of safety of the slope, it has great effect on the interslice forces ratio . Hence the height of thrust line should be determined based on the property of the slope mass. Usually, we use for soft soil and for rigid rock.
5. Conclusion
Considering the limitation of the definition of factor of safety in traditional slice methods, the authors introduce the concept of slicewise factor of safety and propose a new limit equilibrium method based on it. In the method, an assumption that the interslice force ratio is the same between any two slices is first made to ease the solution process. The ratio reflects the exertion degree of the antisliding forces of slip mass. In addition, the eccentric moment is considered in the analysis. A fourvariable implicit equation is established based on the equilibrium of forces and moment and the assumption of interslice forces. And then the interslice force ratio can be solved from the equation using the boundary condition. With the ratio, the slice factors of safety along the slip surface can be obtained straightforwardly. The slicewise factor of safety is useful for engineers to determine the sliding mechanism and passive section of a slope, which can help the engineers with the design of stabilizing piles. The factor of safety of a slope is also defined by the ratio of the sum of antisliding force to the sum of sliding force on the slip surface. The results of the two calculation examples verified the reliability of the proposed method. The height of thrust line has little effect on the factor of slope safety but has great effect on the state of interslice forces.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The study is financially supported by the NSFC projects (no. 51208461 and no. 41202216), and the Fundamental Research Funds for the Central Universities (no. 2014QNA4016 and no. 2014QNA4020).
References
 E. M. Dawson, W. H. Roth, and A. Drescher, “Slope stability analysis by strength reduction,” Géotechnique, vol. 49, no. 6, pp. 835–840, 1999. View at: Publisher Site  Google Scholar
 H. Lin and P. Cao, “A dimensionless parameter determining slip surfaces in homogeneous slopes,” KSCE Journal of Civil Engineering, vol. 18, no. 2, pp. 470–474, 2014. View at: Google Scholar
 H. Lin, W. Zhong, P. Cao, and T. Liu, “Variational safety factors and slip surfaces of slope using threedimensional strength reduction analysis,” Journal of the Geological Society of India, vol. 82, no. 5, pp. 545–552, 2013. View at: Google Scholar
 W. Fellenius, “Calculation of the stability of earth dams,” in Proceedings of the 2nd International Congress on Large Dams, vol. 4, p. 445, Washington, DC, USA, 1936. View at: Google Scholar
 A. W. Bishop, “The use of the slip circle in the stability analysis of slopes,” Géotechnique, vol. 5, no. 1, pp. 7–17, 1955. View at: Google Scholar
 N. R. Morgenstern and V. E. Price, “The analysis of the stability of general slip surfaces,” Géotechnique, vol. 15, no. 1, pp. 79–93, 1965. View at: Google Scholar
 E. Spencer, “A method of analysis of the stability of embankments assuming parallel interslice forces,” Géotechnique, vol. 17, no. 1, pp. 11–26, 1967. View at: Google Scholar
 S. K. Sarma, “Stability analysis of embankments and slopes,” Géotechnique, vol. 23, no. 3, pp. 423–433, 1973. View at: Google Scholar
 J. M. Duncan and A. L. Buchignani, An Engineering Manual for Slope Stability Studies, Virginia Tech, Blacksburg, Va, USA, 1987.
 N. Janbu, “Slope stability computations,” in Embankment Dam Engineering, R. C. Hirschfeld and S. J. P. John, Eds., pp. 47–86, John Wiley & Sons, New York, NY, USA, 1973. View at: Google Scholar
 S. K. Sarma, “Stability analysis of embankments and slopes,” Journal of the Geotechnical Engineering Division, vol. 105, no. 12, pp. 1511–1524, 1979. View at: Google Scholar
 S. K. Sarma and D. Tan, “Determination of critical slip surface in slope analysis,” Géotechnique, vol. 56, no. 8, pp. 539–550, 2006. View at: Publisher Site  Google Scholar
 Y. C. Li, Y. M. Chen, T. L. T. Zhan, D. Ling, and P. J. Cleall, “An efficient approach for locating the critical slip surface in slope stability analyses using a realcoded genetic algorithm,” Canadian Geotechnical Journal, vol. 47, no. 7, pp. 806–820, 2010. View at: Publisher Site  Google Scholar
 S. G. Wright, F. H. Kulhawy, and J. M. Duncan, “Accuracy of equilibrium slope stability analysis,” Journal of the Soil Mechanics and Foundations Division, vol. 99, no. 10, pp. 783–791, 1973. View at: Google Scholar
 F. Tavenas, B. Trak, and S. Leroueil, “Remarks on the validity of stability analyses,” Canadian Geotechnical Journal, vol. 17, no. 1, pp. 61–73, 1980. View at: Publisher Site  Google Scholar
 H. Y. Sun, Y. Zhao, Y. Q. Shang, Y. Yu, and Q. Zhao, “Deepseated slope failures induced by inappropriate cutting in china,” Rock Mechanics and Rock Engineering, vol. 45, no. 6, pp. 1103–1111, 2012. View at: Publisher Site  Google Scholar
 H. Y. Sun, Y. Zhao, Y. Q. Shang, and J. Zhong, “Field measurement and failure forecast during the remediation of a failed cut slope,” Environmental Earth Sciences, vol. 69, no. 7, pp. 2179–2187, 2013. View at: Publisher Site  Google Scholar
 I. B. Donald and P. Giam, “The ACADS slope stability programs review,” in Proceedings of the 6th International Symposium on Landslides, vol. 3, pp. 1665–1670, Christchurch, New Zealand, February 1992. View at: Google Scholar
 J. M. Duncan, “State of the art: limit equilibrium and finiteelement analysis of slopes,” Journal of Geotechnical Engineering, vol. 122, no. 7, pp. 577–596, 1996. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Yu Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.