Mathematical Problems in Engineering
Volume 2017 (2017), Article ID 3569826, 6 pages
https://doi.org/10.1155/2017/3569826

Research Article

A New Calculation Method for the Soil Slope Safety Factor

Fu Zhu,1,2 Wanxi Zhang,¹ and Mingzhi Sun¹

¹College of Material Science and Engineering, Jilin University, Changchun 130022, China
²School of Transportation Science and Engineering, Jilin Jianzhu University, Changchun 130118, China

Correspondence should be addressed to Mingzhi Sun; mzsun14@mails.jlu.edu.cn

Received 27 February 2017; Revised 25 July 2017; Accepted 6 September 2017; Published 9 October 2017

Academic Editor: Dane Quinn

Copyright © 2017 Fu Zhu 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.

Abstract.

Based on the unified strength theory, a new method to calculate the plane soil slope safety factor was derived that considers the effect of intermediate principal stress and at-rest lateral pressure coefficient . Calculation examples from the literature were used to compare the new calculation method and the current slice method; the results showed that both provided good consistency. The new method can provide a reference for slope stability evaluation. The new method was used to calculate the soil slope safety factors for different values of intermediate principal stress parameter , double shear stress parameters , and static lateral pressure coefficient . The results showed that the safety factor increased when was increased; first increased and then decreased when was increased; and increased when was increased. These results show that the intermediate principal stress as well as the stress state and its changes cannot be ignored during soil slope stability analysis. The slope soil characteristics and stress state should be considered to determine the unified strength theoretical parameters and static lateral pressure coefficient, maximize the potential of slope soil strength, and effectively reduce the costs of soil slope engineering.

1. Introduction

In roads, bridges, and construction projects, slope stability problems are often encountered during cutting or foundation pit excavation. Slope instability is due to the destruction of the original soil stress state of equilibrium caused by external forces—such as cutting or foundation pit excavation—and the reduction of soil antishear strength by the influence of various external factors—such as rainwater intrusion and soil freeze-thaw. In practical engineering, slope stability is analyzed to test the reasonableness of the soil slope section design. If the slope is too steep, it will slump easily; if the slope is too gentle, it will increase the amount of earthwork needed.

The characteristics of the Ordinary Method of Slices [1], Bishop’s Modified Method [2], Force Equilibrium Methods (e.g., Lowe and Karafiath [3]), Janbu’s Genralized Procedure of Slices [4], Morgenstern and Price’s Method [5], and Spencer’s Method [6] were summarized in most textbooks. Fall et al. [7] have conducted study on the stability analysis of landslides by the finite element method. Cheng and Yip [8] have shown that rigorous method is necessary to give a reliable estimation of stability of landslides in 3D analyses. Zhu and Lee [9] conducted study on the factor of safety based on Bell’s assumption. Bell’s method was improved by Zheng and Tham [10]. Zheng and Tham’ method can be regarded subsequently as the enhancement of Fellenius’ method.

The slope stability safety factor refers to the ratio of the soil shear strength to the shear stress of a possible sliding surface in the slope. The soil stress state and its changes are prerequisites of slope stability; the existing slope circular slipping method (Petterson 1916) and slice method (Fellenius 1927) ignore the impact of the stress state. In reality, the slope stability changes with changes in the stress state. Researchers [11–14] are currently looking for the sliding center and slip surface, supplementing and modifying the basic assumptions of the slice method, and providing a fundamental basis in engineering applications for the slice method. However, defects of slice method and statically indeterminate problem of this method [15] have produced challenges in practical engineering applications.

Based on multislip mechanism and the model of multishear element, Yu established the unified strength theory which takes into consideration the different contribution of all stress components on the yield of failure of materials [16, 17]. The unified strength theory encompasses the twin shear strength theory [18–20] and single strength theory. The excellent agreement between the predicted results by the unified strength theory and the experiment results indicates that the unified strength theory is applicable for a wide range of stress states in many materials (Ma et al. 1985 [21]).

The coefficient of the earth pressure at rest () is defined as the ratio of the in situ horizontal effective stress to the in situ vertical effective stress. The parameter is needed in the interpretation of laboratory and in situ tests and design of retaining structures and excavation support systems. Schnaid and Yu [22] believe that is an important input parameter for the numerical analyses of geotechnical boundary value problems.

In this study, the perspective of a total stress state was considered to derive a new calculation method for the soil slope safety factor based on the unified strength theory. The slope safety factor was defined by considering the effects of the intermediate principal stress and at-rest lateral pressure coefficient. The method was compared and verified with the current slice method and can provide a reference for stability evaluation in soil slope engineering.

2. Basic Theory and Formula Derivation

2.1. Unified Strength Theory

The Mohr–Coulomb strength theory is simple and practical. It is conducive to engineering applications but does not reflect the effect of intermediate principal stress, and the calculated results are relatively conservative. In 1991, Yu suggested the unified strength theory to compensate for the shortcomings of the Mohr–Coulomb strength theory. The unified strength theory can consider the intermediate principal stress effect of the material and can simulate nearly all materials on the partial plane to develop the material strength potentials. There are two equations with a conditional formula for both the mathematical model and the theoretical expression of the unified strength theory, which considers the different contributions of various stress components to the material yield and destruction, reduces the number of material parameters, and makes the limiting surface to reach the outer boundary of the convex criteria; these cannot be achieved by other criteria. The mathematical modeling methods from the two equations can also be used to solve issues with determining the intermediate principal shear stress. Yu derived a mathematical expression of the unified strength theory using the unified twin shear model and a new mathematical model [17]:

where and are yield functions; is the ratio of the tensile to compression strengths of the material; is the tensile strength; and are the cohesion and internal friction angle, respectively, of rock and soil; and is the selected failure criterion introduced in the unified strength theory that also reflects the destructive effects of the intermediate principal shear stress and normal stress of the corresponding surface on the material.

The unified strength theory was converted into a formula similar to the Mohr–Coulomb strength theory in order to obtain the friction angle and unified cohesive force ; these are expressed by the internal friction angle and cohesion as follows [23].

When ,

When ,

where is the twin shear stress parameter.

The unified internal friction angle and unified cohesion can be used to express the Mohr–Coulomb strength theory:

2.2. Basic Assumptions

(1)The excavated soil is simplified as plane soil slopes.(2)Soil is homogeneous.(3)The stress state can be represented by formula (4).(4)The soil antishear strength satisfies formula (3).(5)The horizontal stress of inner points along the depth causes soil slope instability.(6)During excavation of the soil, the static lateral pressure coefficient remains unchanged.(7)The impact of pore water and groundwater is not considered.

2.3. Formula Derivation

The extending direction of the soil slope is taken as the plane stress state, and elastic half-space plane stress analysis is performed for soil in a steady state under gravity stress; the main stress expression of any point is as follows:

where is the maximum principal stress, is the minimum principal stress, is the soil gravity, is the static lateral pressure coefficient of the soil, is the distance from the ground surface to any point, is the horizontal stress at any point, and is the vertical stress at any point.

As shown in Figure 1, assuming that the cutting or foundation soil is not excavated, when slope angle , then . In vertical excavation of the soil, when slope angle , then . The horizontal stress of the soil excavation slope angle satisfies the following formula:

Figure 1: Simplified soil slope excavation.

According to the conventional definition of the factor of safety for some point within a soil mass, the safety factor is the ratio between the shear strength and shear stress at that point [24, 25]:

In line with the Mohr–Coulomb criterion for the shear strength of soils, for stresses at some point within a soil mass, differences in the magnitude of shear stress in an arbitrary direction will result in shear strength differences. In other words, the factor of safety at a point in a soil mass, defined in (6), will vary with the direction. This leads to complexities and difficulties in the methods of slope stability analysis and to a variety of assumptions in calculational theories. To ensure the uniqueness of safety factors computed at each point within a soil mass, the factor of safety was defined as described below.

Given a point with a determined stress state within some soil mass, its margin of safety is the ratio between the shear strength corresponding to the maximum shear strength at that point and the total maximum shear strength as described in Figure 2.

Figure 2: Schematic diagram of accumulated safety in soil.

Then, the safety margin of a slope is the ratio between the cumulative shear strength and the cumulative maximum shear stress within the slope’s height; thus

since

Therefore, (9) are substituted into (8) to obtain

3. Safety Factor Calculation and Analysis

The new approach was derived from cutting or foundation pit soil excavation. Slope engineering calculation examples in textbooks [26, 27] and literature [11] were used to verify the general application of the calculation method with formula (10).

Example 1. Slope height  m, slope angle , soil gravity kN/m³, angle of soil internal friction , and cohesion  kPa were known. The Fellenius’ slice method and Bishop formula were used to calculate the slope safety factor; the results were 1.18 and 1.19, respectively [9].
Using and , which were calculated using Jaky’s formula, formula (10) was used to calculate the slope safety factor, which was 0.98.

Example 2. Question EX1 (c) from Australian Association of Computer Applications (ACADS) assessment in 1987: a heterogeneous soil slope of material properties was shown in Table 1, and slope shape was shown in Figure 3. The problem was simplified into one with a homogeneous soil slope with a height of  m, slope gradient = 1 : 2, soil gravity  kN/m³, soil internal friction angle , and cohesion  kPa. The Bishop formula and genetic algorithm were used in the calculation, and the ten safety factors of the slide facing were in the range of 1.398~1.40 [11]. Using and , which were calculated using Jaky’s formula, the slope safety factor was calculated using formula (10) to be 1.34. The reference value was 1.39.

Table 1: Material properties.

Soil number
Gravity (kN/m³)
Cohesion (kN/m²)
Internal friction angle (°)

#1 soil
19.5
0.0
38.0

#2 soil
19.5
5.3
23.0

#3 soil
19.5
7.2
20.0

Figure 3: Contour map of ACADS assessment question EX1 (c).

The examples above show that the new method is simpler than the current slice method, has a clearer theoretical basis and concept, does not require programming, and effectively reduces the computational workload. It can be used as a general basis for safety evaluation of slope stability.

In Example 1, the data substituted into formulas (2a) and (2b) were used to analyze the relationship among the unified internal friction angle , cohesion , , and . The calculation results are shown in Figures 4 and 5.

Figure 4: Relationship among , , and .

Figure 5: Relationship among , , and .

Figures 4 and 5 show that the unified internal friction angle and cohesion increased when was increased. They firstly increased and then decreased when was increased and reached their maximum when .

and were obtained with different and values and substituted into formula (10) to calculate various soil slope safety factors , as shown in Figure 6. Figure 6 shows that when was a fixed value, the safety factor increased with increasing . When was a fixed value, changed from 0 to 1, and the safety factor first increased and then decreased. When , the safety factor was at its maximum; when , or , the safety factor was at its minimum value. When , formula (10) degrades to a formula under the Mohr–Coulomb strength theory, and the obtained safety factor will be the smallest fixed value. This indicates that the potential of the soil is far from being developed and may result in significant waste.

Figure 6: Relationship among , , and .

Using , and obtained with different and were substituted into formula (10), and various soil slope safety factors were obtained as shown in Figure 7. Figure 7 shows that when was a fixed value, increased with , which is consistent with the above conclusions. When was a fixed value, varied from 0.39 to 0.79, and the safety factor increased gradually. This indicates that and are combined factors that determine slope stability. Formula (10) considers the effects of both and can provide a reference for safety evaluation of a soil slope.

Figure 7: Relationship among , , and .

4. Conclusions

The new method was verified by comparing results with the existing slice method using examples in the literature. The new method can be combined with slope engineering to obtain parameters of the unified strength theory and the static lateral pressure coefficient. The slope stability and safety were evaluated with the new method.

The effects of different factors on soil slope safety and stability were analyzed, including the intermediate principal stress parameter, twin shear stress state parameters, and static lateral pressure coefficient. These results indicate that the intermediate principal stress and static lateral pressure coefficient cannot be ignored in slope stability analysis.

This study only examined the effects of unified strength theory parameters and the static lateral pressure coefficient on the slope safety factor. In order to determine the parameters and practical applications of the new method, further research and verification are needed.

The theoretical formula was derived, calculated, and analyzed from the perspective of the total stress state. The effects on pore water pressure and groundwater should be further examined from the perspective of effective stress.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research is financially supported by the Project of Education Department of Jilin Province (JJKH20170260KJ), the Project of Ministry of Housing and Urban-Rural Development (2017-K4-004), and the Plan Projects of Transportation Science and Technology in Jilin Province of China (2011103).

References

W. Fellenius, “Calculation of the stability of earth dams,” in In Proceedings of the 2nd Congress on Large Dams, vol. 4, pp. 445–463, Wash, USA, 1936.

A. W. Bishop, “The use of the slip circle in the stability analysis of slopes,” Geotechnique, vol. 5, pp. 7–17, 1955.

J. Lowe and L. Karafiath, “Stability of earth dams upon drawdown,” in Proceedings of the 1st Pan-Am Conf. On Soil Mech. And Founf. Engrg., pp. 537–552, Mexico City, Mexico, 1960.

N. Janbu, “Slope stability computations,” in Soil Mechanics and Foundation Engineering Journal, Technical University of Norway, Trondheim, Norway, 1968.

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.

E. Spencer, “A method of analysis of the stability of embankments assuming parallel inter-slice forces,” Geotechnique, vol. 17, no. 1, pp. 11–26, 1967.

M. Fall, R. Azzam, and C. Noubactep, “A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping,” Engineering Geology, vol. 82, no. 4, pp. 241–263, 2006.

Y. M. Cheng and C. J. Yip, “Three dimensional asymmetrical slope stability analysis extension of bishop's, Janbu's, and Morgenstern-Price's Techniques,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 133, no. 12, pp. 1544–1555, 2007.

D. Y. Zhu and C. F. Lee, “Explicit limit equilibrium solution for slope stability,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 26, no. 15, pp. 1573–1590, 2002.

H. Zheng and L. G. Tham, “Improved Bell's method for the stability analysis of slopes,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 33, no. 14, pp. 1673–1689, 2009.

W. J. Lv, X. J. Li, and H. H. Zhu, “GA-based generalized slope stability analysis method,” Chinese Journal of Geotechnical Engineering, vol. 27, no. 5, pp. 595–599, 2005.

A. I. H. Malkawi, W. F. Hassan, and S. K. Sarma, “Global search method for locating general slip surface using Monte Carlo techniques,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 127, no. 8, pp. 688–698, 2001.

G.-D. Zou, “A global optimization method of the slice method for slope stability analysis,” Chinese Journal of Geotechnical Engineering, vol. 24, no. 3, pp. 309–312, 2002.

G. Zou and Y. Chen, “Coupling algorithm of simulated annealing algorithm and random search method for slope stability analysis,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, no. 12, pp. 2032–2037, 2004.

Z. Y. Chen, X. G. Wang, Y. C. Xing et al., “Theoretical study and physical modeling on “Principle of Maximum” in slope stability analysis,” Chinese Journal of Geotechnical Engineering, vol. 27, no. 5, pp. 495–499, 2005.

M. H. Yu, New System of Strength Theory, Xian Jiaotong University Press, 1992.

M.-H. Yu, “Unified strength theory for geomaterials and its application,” Chinese Journal of Geotechnical Engineering, vol. 16, no. 2, pp. 1–9, 1994.

M.-H. Yu, “Twin shear stress yield criterion,” International Journal of Mechanical Sciences, vol. 25, no. 1, pp. 71–74, 1983.

M. H. Yu, L. He, and L. Song, “Twin shear stress theory and its generalization,” Scientia Sinica, Series A-Mathematical, Physical, Astronomical and Technical Sciences, vol. 28, pp. 1174–1183, 1985.

M. H. Yu and L. N. He, “A new model and theory on yield and failure of materials under complex stress state,” Mechanical Behaviour of Materials, vol. 6, no. 3, pp. 841–846, 1991.

M. Yu, J. Li, and Y. Zhang, “Unified characteristics line theory of spacial axisymmetric plastic problem,” Science in China, Series E: Technological Sciences, vol. 44, no. 2, pp. 207–215, 2001.

F. Schnaid and H. S. Yu, “Interpretation of the seismic cone test in granular soils,” Geotechnique, vol. 57, no. 3, pp. 265–272, 2007.

W. Fan, L. S. Deng, X. Y. Bai et al., “Application of unified strength theories to the stability of slope analysis,” Coal Geology & Exploration, vol. 35, no. 1, pp. 63–66, 2007.

G. T. Wang, “The method and formula for stability factor of slope with state of soil original stress,” Chinese Journal of Engineering Science, vol. 8, no. 12, pp. 80–84, 2006.

G. T. Wang, “Illuminating the new method process of slope stability analysis by soil stress state,” Chinese Journal of Engineering Science, vol. 12, no. 1, pp. 52–55, 2010.

Z. Y. Chen, Soil Slope Stability Analysis-Theory, Methods and Programs, China Water Power Press, Beijing, China, 2003.

J. Y. Yuan, J. G. Qian, H. M. Zhang et al., Soil Properties and Soil Mechanics, China Communication Press, Beijing, china, 2009.