1. Fredlund DG, Krahn J. Comparison of slope stability methods of analysis. Canadian Geotechnical Journal. 1977;14(3):429–439.
2. Li KS. A unified solution scheme for slope stability analysis. Proceedings of the 5th International Symposium on Landslides, New Zealand; 1992; Rotterdam, The Netherlands. Balkema; pp. 481–486.
3. Niu W-J, Ye W-M, Liu S-G, Yu H-T. Limit analysis of a soil slope considering saturated-unsaturated seepage. Rock and Soil Mechanics. 2009;30(8):2477–2482.
4. Chen WF. Limit Analysis and Soil Plasticity. Amsterdam, The Netherlands: Elsevier; 1975.
5. Michalowski RL. Slope stability analysis: a kinematical approach. Geotechnique. 1995;45(2):283–293.
6. Donald IB, Chen ZY. Slope stability analysis by the upper bound approach: fundamentals and methods. Canadian Geotechnical Journal. 1997;34(6):853–862.
7. Sun JP, Li JC, Liu QQ. Search for critical slip surface in slope stability analysis by spline-based GA method. Journal of Geotechnical and Geoenvironmental Engineering. 2008;134(2):252–256.
8. Javadzadeh R, Javadzadeh E. Locating critical failure surface in rock slope stability with hybrid model based on artificial immune system and Cellular Learning Automata (CLA-AIS) World Academy of Science, Engineering and Technology. 2011;59:662–665.
9. Gitirana FNG, Jr., Fredlund DG. Analysis of transient embankment stability using the dynamic programming method. Proceedings of the 56th Canadian Geotechnical Conference; 2003.
10. Pham HTV, Fredlund DG. Dynamic programming method in slope stability computations. Proceedings of the 12th Asian Regional Conference on Soil Mechanics and Geotechnica Engineering; August 2003; Singapore.
11. Pham HTV, Fredlund DG, Gitirana FNG., Jr. Slope stability analysis using dynamic prograamming combined with finite element stress analysis. Proceedings of the International Conference on the Management of the Land and Water Resources (MLWR '01); October 2001; Hanoi, Vietnam.
12. Pham HTV, Fredlund DG. The application of dynamic programming to slope stability analysis. Canadian Geotechnical Journal. 2003;40(4):830–847.
13. Taha MR, Khajehzadeh M, El-Shafie A. Slope stability assessment using optimization techniques: an overview. Electronic Journal of Geotechnical Engineering. 2010;15:1901–1915.
14. Rao AU, Sabhahit N. Genetic algorithm in stability of non-homogeneous slopes. Proceedings of the 12th International Conference of International Association for Computer Methods and Advances in Geomechanics (IACMAG '08); October 2008; Goa, India.
15. Nie Y-G, Liu W-Q, Shi J-Y, Zhao W-B. Application of accelerating genetic algorithm for embankment slope stability analysis. China Journal of Highway and Transport. 2003;16(4):p. 16.
16. Zhou Y-P, Li S-J, Liu Y-X, Sun H-L, Jiang F. Searching critical failure surface in slope stability analysis with genetic algorithm. Chinese Journal of Rock Mechanics and Engineering. 2005;24:5226–5230.
17. McCombie P. Critical failure surface location using Simple Genetic Algorithm and multiple wedge slope stability. Proceedings of the 1st International Symposium on Computational Geomechanics; May 2009; Juan-les-Pins, France. pp. 900–909.
18. Gavin K, Xue JF. Use of a genetic algorithm to perform reliability analysis of unsaturated soil slopes. Geotechnique. 2009;59(6):545–549.
19. Wikipedia. Genetic algorithm. 2012, http://en.wikipedia.org/wiki/Genetic_algorithm.
20. China University of Geosciences. Cohesive soil slope stability analysis. http://course.cug.edu.cn/cug/soil_mechanics/COURSE/CHAPTER7/Chap7_3_5.htm.
21. Thomaz JE. JTRP Technical Reports, Purdue University; 1986. A general method for three dimensional slope stability analysis: informational report.
22. Zhang TB. Chengdu University of Science and Technology Press; 1987. Earth slope stability analysis and geotechnology slope design.
23. He BG. Application Examples in Civil Engineering with Ansys. Beijing, China: China Water Power Press; 2011.
24. Liu Y, Yang Q, Zhu L. Abutment stability analysis of arch dam based on 3D nonlinear finite element method. Chinese Journal of Rock Mechanics and Engineering. 2008;27(supplement 1):3222–3228.
25. Liu YR, He Z, Li B, Yang Q. Slope stability analysis based on a multigrid method using a nonlinear 3D finite element model. Frontiers of Structural and Civil Engineering. 2013;7(1):24–31.
26. Wikipedia. Drucker–Prager yield criterion. 2013, http://en.wikipedia.org/wiki/Drucker%E2%80%93Prager_yield_criterion.
27. Doran B, Koksal HO, Polat Z, Karakoc C. The use of “Drucker-Prager Criterion” in the analysis of reinforced concrete members by finite elements. Teknik Dergi. 1998;9(2):1617–1625.
28. Zoran B, Verka P, Biljana M. Mathematical modeling of materially nonlinear problems in structural analyses (part I—theoretical fundamentals) Facta Universitatis. Architecture and Civil Engineering. 2010;8(1):67–78.
29. Yang Q, Chen X, Zhou W-Y. A practical 3D elasto-plastic incremental method in FEM based on D-P yield criteria. Chinese Journal of Geotechnical Engineering. 2002;24(1):16–20.
30. Zhang GD. Soil Constitutive Model and Its Application in Engineering. Beijing, China: Science Press; 1995.
31. Rakić D, Živković M, Slavković R, Kojić M. Stress integration for the Drucker-Prager material model without hardening using the Incremental Plasticity Theory. Journal of the Serbian Society For Computational Mechanics. 2008;2(1):80–89.
32. Zhang XY. Soil and Rock Plastic Mechanics. Beijing, China: China Communications Press; 1993.
33. Zheng YR, Shen ZJ, Gong XN. The Principles of Geotechnical Plastic Mechanics. Beijing, China: China Architecture and Building Press; 2002.
34. Akers SA. Two-dimensional finite element analysis of porous geomaterials at multikilobar stress levels [PhD thesis] Virginia Polytechnic Institute and State University; 2001.
35. Tan X-H, Wang J-G, Wang Y. Nonlinear finite element analysis of slope stability. Rock and Soil Mechanics. 2008;29(8):2047–2050.
For a broader coverage related to this topic, see Slope stability.
Slope stability analysis is performed to assess the safe design of a human-made or natural slopes (e.g. embankments, road cuts, open-pit mining, excavations, landfills etc.) and the equilibrium conditions.Slope stability is the resistance of inclined surface to failure by sliding or collapsing. The main objectives of slope stability analysis are finding endangered areas, investigation of potential failure mechanisms, determination of the slope sensitivity to different triggering mechanisms, designing of optimal slopes with regard to safety, reliability and economics, designing possible remedial measures, e.g. barriers and stabilization.
Successful design of the slope requires geological information and site characteristics, e.g. properties of soil/rock mass, slope geometry, groundwater conditions, alternation of materials by faulting, joint or discontinuity systems, movements and tension in joints, earthquake activity etc. The presence of water has a detrimental effect on slope stability. Water pressure acting in the pore spaces, fractures or other discontinuities in the materials that make up the pit slope will reduce the strength of those materials. Choice of correct analysis technique depends on both site conditions and the potential mode of failure, with careful consideration being given to the varying strengths, weaknesses and limitations inherent in each methodology.
Before the computer age stability analysis was performed graphically or by using a hand-held calculator. Today engineers have a lot of possibilities to use analysis software, ranges from simple limit equilibrium techniques through to computational limit analysis approaches (e.g. Finite element limit analysis, Discontinuity layout optimization) to complex and sophisticated numerical solutions (finite-/distinct-element codes). The engineer must fully understand limitations of each technique. For example, limit equilibrium is most commonly used and simple solution method, but it can become inadequate if the slope fails by complex mechanisms (e.g. internal deformation and brittle fracture, progressive creep, liquefaction of weaker soil layers, etc.). In these cases more sophisticated numerical modelling techniques should be utilised. Also, even for very simple slopes, the results obtained with typical limit equilibrium methods currently in use (Bishop, Spencer, etc.) may differ considerably. In addition, the use of the risk assessment concept is increasing today. Risk assessment is concerned with both the consequence of slope failure and the probability of failure (both require an understanding of the failure mechanism).
Within the last decade (2003) Slope Stability Radar has been developed to remotely scan a rock slope to monitor the spatial deformation of the face. Small movements of a rough wall can be detected with sub-millimeter accuracy by using interferometry techniques.
Limit equilibrium analysis
Conventional methods of slope stability analysis can be divided into three groups: kinematic analysis, limit equilibrium analysis, and rockfall simulators. Most slope stability analysis computer programs are based on the limit equilibrium concept for a two- or three-dimensional model. Two-dimensional sections are analyzed assuming plane strain conditions. Stability analyses of two-dimensional slope geometries using simple analytical approaches can provide important insights into the initial design and risk assessment of slopes.
Limit equilibrium methods investigate the equilibrium of a soil mass tending to slide down under the influence of gravity. Translational or rotational movement is considered on an assumed or known potential slip surface below the soil or rock mass. In rock slope engineering, methods may be highly significant to simple block failure along distinct discontinuities. All these methods are based on the comparison of forces, moments, or stresses resisting movement of the mass with those that can cause unstable motion (disturbing forces). The output of the analysis is a factor of safety, defined as the ratio of the shear strength (or, alternatively, an equivalent measure of shear resistance or capacity) to the shear stress (or other equivalent measure) required for equilibrium. If the value of factor of safety is less than 1.0, the slope is unstable.
All limit equilibrium methods assume that the shear strengths of the materials along the potential failure surface are governed by linear (Mohr-Coulomb) or non-linear relationships between shear strength and the normal stress on the failure surface. The most commonly used variation is Terzaghi's theory of shear strength which states that
where is the shear strength of the interface, is the effective stress ( is the total stress normal to the interface and is the pore water pressure on the interface), is the effective friction angle, and is the effective cohesion.
The methods of slices is the most popular limit equilibrium technique. In this approach, the soil mass is discretized into vertical slices. Several versions of the method are in use. These variations can produce different results (factor of safety) because of different assumptions and inter-slice boundary conditions.
The location of the interface is typically unknown but can be found using numerical optimization methods. For example, functional slope design considers the critical slip surface to be the location where that has the lowest value of factor of safety from a range of possible surfaces. A wide variety of slope stability software use the limit equilibrium concept with automatic critical slip surface determination.
Typical slope stability software can analyze the stability of generally layered soil slopes, embankments, earth cuts, and anchored sheeting structures. Earthquake effects, external loading, groundwater conditions, stabilization forces (i.e., anchors, geo-reinforcements etc.) can also be included.
Analytical techniques: Method of slices
Many slope stability analysis tools use various versions of the methods of slices such as Bishop simplified, Ordinary method of slices (Swedish circle method/Petterson/Fellenius), Spencer, Sarma etc. Sarma and Spencer are called rigorous methods because they satisfy all three conditions of equilibrium: force equilibrium in horizontal and vertical direction and moment equilibrium condition. Rigorous methods can provide more accurate results than non-rigorous methods. Bishop simplified or Fellenius are non-rigorous methods satisfying only some of the equilibrium conditions and making some simplifying assumptions. Some of these approaches are discussed below.
Swedish Slip Circle Method of Analysis
The Swedish Slip Circle method assumes that the friction angle of the soil or rock is equal to zero, i.e., . In other words, when friction angle is considered to be zero, the effective stress term goes to zero, thus equating the shear strength to the cohesion parameter of the given soil. The Swedish slip circle method assumes a circular failure interface, and analyzes stress and strength parameters using circular geometry and statics. The moment caused by the internal driving forces of a slope is compared to the moment caused by forces resisting slope failure. If resisting forces are greater than driving forces, the slope is assumed stable.
Ordinary Method of Slices
In the method of slices, also called OMS or the Fellenius method, the sliding mass above the failure surface is divided into a number of slices. The forces acting on each slice are obtained by considering the mechanical (force and moment) equilibrium for the slices. Each slice is considered on its own and interactions between slices are neglected because the resultant forces are parallel to the base of each slice. However, Newton's third law is not satisfied by this method because, in general, the resultants on the left and right of a slice do not have the same magnitude and are not collinear.
This allows for a simple static equilibrium calculation, considering only soil weight, along with shear and normal stresses along the failure plane. Both the friction angle and cohesion can be considered for each slice. In the general case of the method of slices, the forces acting on a slice are shown in the figure below. The normal () and shear () forces between adjacent slices constrain each slice and make the problem statically indeterminate when they are included in the computation.
For the ordinary method of slices, the resultant vertical and horizontal forces are
where represents a linear factor that determines the increase in horizontal force with the depth of the slice. Solving for gives
Next, the method assumes that each slice can rotate about a center of rotation and that moment balance about this point is also needed for equilibrium. A balance of moments for all the slices taken together gives
where is the slice index, are the moment arms, and loads on the surface have been ignored. The moment equation can be used to solve for the shear forces at the interface after substituting the expression for the normal force:
Using Terzaghi's strength theory and converting the stresses into moments, we have
where is the pore pressure. The factor of safety is the ratio of the maximum moment from Terzaghi's theory to the estimated moment,
Modified Bishop’s Method of Analysis
The Modified Bishop’s method is slightly different from the ordinary method of slices in that normal interaction forces between adjacent slices are assumed to be collinear and the resultant interslice shear force is zero. The approach was proposed by Alan W. Bishop of Imperial College. The constraint introduced by the normal forces between slices makes the problem statically indeterminate. As a result, iterative methods have to be used to solve for the factor of safety. The method has been shown to produce factor of safety values within a few percent of the "correct" values.
The factor of safety for moment equilibrium in Bishop's method can be expressed as
where, as before, is the slice index, is the effective cohesion, is the effective internal angle of internal friction, is the width of each slice, is the weight of each slice, and is the water pressure at the base of each slice. An iterative method has to be used to solve for because the factor of safety appears both on the left and right hand sides of the equation.
Lorimer's Method is a technique for evaluating slope stability in cohesive soils. It differs from Bishop's Method in that it uses a clothoid slip surface in place of a circle. This mode of failure was determined experimentally to account for effects of particle cementation. The method was developed in the 1930s by Gerhardt Lorimer (Dec 20, 1894-Oct 19, 1961), a student of geotechnical pioneer Karl von Terzaghi.
Spencer’s Method of analysis requires a computer program capable of cyclic algorithms, but makes slope stability analysis easier. It is not as accurate as the Modified Bishop’s method, but is acceptably accurate in engineering practices.
Main article: Sarma method
The Sarma method, proposed by Sarada K. Sarma of Imperial College is a Limit equilibrium technique used to assess the stability of slopes under seismic conditions. It may also be used for static conditions if the value of the horizontal load is taken as zero. The method can analyse a wide range of slope failures as it may accommodate a multi-wedge failure mechanism and therefore it is not restricted to planar or circular failure surfaces. It may provide information about the factor of safety or about the critical acceleration required to cause collapse.
The assumptions made by a number of limit equilibrium methods are listed in the table below.
|Ordinary method of cells||Interslice forces are neglected|
|Bishop's simplified/modified ||Resultant interslice forces are horizontal. There are no interslice shear forces.|
|Janbu's simplified||Resultant interslice forces are horizontal. An empirical correction factor is used to account for interslice shear forces.|
|Janbu's generalized||An assumed line of thrust is used to define the location of the interslice normal force.|
|Spencer ||The resultant interslice forces have constant slope throughout the sliding mass.|
|Chugh||Same as Spencer's method but with a constant acceleration force on each slice.|
|Morgenstern-Price||The direction of the resultant interslice forces is defined using an arbitrary function. The fractions of the function value needed for force and moment balance is computed.|
|Fredlund-Krahn (GLE) ||Similar to Morgenstern-Price.|
|Corps of Engineers ||The resultant interslice force is either parallel to the ground surface or equal to the average slope from the beginning to the end of the slip surface..|
|Lowe and Karafiath ||The direction of the resultant interslice force is equal to the average of the ground surface and the slope of the base of each slice.|
|Sarma ||The shear strength criterion is applied to the shears on the sides and bottom of each slice. The inclinations of the slice interfaces are varied until a critical criterion is met.|
The table below shows the statical equilibrium conditions satisfied by some of the popular limit equilibrium methods.
|Method||Force balance (vertical)||Force balance (horizontal)||Moment balance|
|Janbu's generalized||Yes||Yes||Used to compute interslice shear forces|
|Corps of Engineers||Yes||Yes||No|
|Lowe and Karafiath||Yes||Yes||No|
Rock slope stability analysis
Rock slope stability analysis based on limit equilibrium techniques may consider following modes of failures:
- Planar failure -> case of rock mass sliding on a single surface (special case of general wedge type of failure); two-dimensional analysis may be used according to the concept of a block resisting on an inclined plane at limit equilibrium
- Polygonal failure -> sliding of a nature rock usually takes place on polygonally-shaped surfaces; calculation is based on a certain assumptions (e.g. sliding on a polygonal surface which is composed from N parts is kinematically possible only in case of development at least (N - 1) internal shear surfaces; rock mass is divided into blocks by internal shear surfaces; blocks are considered to be rigid; no tensile strength is permitted etc.)
- Wedge failure -> three-dimensional analysis enables modelling of the wedge sliding on two planes in a direction along the line of intersection
- Toppling failure -> long thin rock columns formed by the steeply dipping discontinuities may rotate about a pivot point located at the lowest corner of the block; the sum of the moments causing toppling of a block (i.e. horizontal weight component of the block and the sum of the driving forces from adjacent blocks behind the block under consideration) is compared to the sum of the moments resisting toppling (i.e. vertical weight component of the block and the sum of the resisting forces from adjacent blocks in front of the block under consideration); toppling occur if driving moments exceed resisting moments
Limit equilibrium analysis software
- SLIDE provides 2D stability calculations in rocks or soils using these rigorous analysis methods: Spencer,Sarma and Morgenstern-Price/General limit equilibrium; and non-rigorous methods: Bishop simplified, Corps of Engineers, Janbu simplified/corrected, Lowe-Karafiath and Ordinary/Fellenius. Searching of the critical slip surface is realized with the help of a grid or as a slope search in user-defined area. Program includes also probabilistic analysis using Monte Carlo or Latin Hypercube simulation techniques where any input parameter can be defined as a random variable. Probabilistic analysis determine the probability of failure and reliability index, which gives better representation of the level of safety. Back analysis serves for calculation of a reinforcement load with a given required factor of safety. Program enables finite element groundwater seepage analysis.
- SLOPE/W is formulated in terms of moment and force equilibrium factor of safety equations. Limit equilibrium methods include Morgenstern-Price, General limit equilibrium, Spencer, Bishop, Ordinary, Janbu etc. This program allows integration with other applications. For example, finite element computed stresses from SIGMA/W or QUAKE/W can be used to calculate a stability factor by computing total shear resistance and mobilized shear stress along the entire slip surface. Then a local stability factor for each slice is obtained. Using a Monte Carlo approach, program computes the probability of failure in addition to the conventional factor of safety.STABL WV is a limit equilibrium-based, Windows software based on the stabl family of algorithms. It allows analysis using Bishop's, Spencer's and Janbu's method. Regular slopes as well as slopes with various types of inclusions may be analyzed.
- HYDRUS add-on modules can check the stability of embankments, dams, earth cuts and anchored sheeting structures with the influence of the water. The values of the pore pressure in transport domain are imported automatically for the selected time to Stability module. The analysis can be repeated for all time shots of the water movement simulated by basic program. The common method of slices (the Bishop, Fellenius/Petterson, Morgenstern-Price or the Spencer) can be set as well as the different type of Geo-reinforcement or Earthquake effects.
- SVSlope is formulated in terms of moment and force equilibrium factor of safety equations. Limit equilibrium methods include Morgenstern-Price, General limit equilibrium, Spencer, Bishop, Ordinary, Kulhawy and others This program allows integration with other applications in the geotechnical software suite. For example, finite element computed stresses from SVSolid or pore-water pressures from SVFlux can be used to calculate the factor of safety by computing total shear resistance and mobilized shear stress along the entire slip surface. The software also utilizes Monte Carlo, Latin Hypercube, and the APEM probabilistic approaches. Spatial variability through random fields computations may also be included in the analysis.
- dotSlope provides limit equilibrium analyses through the methods of Fellenius, Bishop simplified, Janbu simplified/corrected, Corps of Engineers, Lowe & Karafiath, Spencer, GLE/Morgenstern & Price. The slope can have multiple soils, impenetrable layers, cuts and embankments, multiple groundwater conditions, ponded water, dry and water filled tension cracks, soil reinforcements (anchors, nails, piles and geo-synthetics). Slip surfaces can be defined through six surface generators in order to find the critical case. dotSlope can run deterministic analyses, surface optimization, sensitivity analyses along with probabilistic analyses using the Monte Carlo method. The program runs on Microsoft Windows, Mac OS X and Android.
- GALENA - includes stability analysis, back analysis, and probability analysis, using the Bishop, Spencer-Wright and Sarma methods.
- GSLOPE - provides limit equilibrium slope stability analysis of existing natural slopes, unreinforced man-made slopes, or slopes with soil reinforcement, using Bishop’s Modified method and Janbu’s Simplified method applied to circular, composite or non-circular surfaces.
- CLARA-W - three-dimensional slope stability program includes calculation with the help of these methods: Bishop simplified, Janbu simplified, Spencer and Morgenstern-Price. Problem configurations can involve rotational or non-rotational sliding surfaces, ellipsoids, wedges, compound surfaces, fully specified surfaces and searches.
- TSLOPE3 - two- or three-dimensional analyses of soil and rock slopes using Spencer method.
- A program specific for rock slope analysis isAutoBlock. It overcome two problems inherently complicating the analysis in engineering practice: firstly determining volumes and magnitudes of sliding areas of potentially unstable blocks based on the real topography, and secondly finding the critical blocks which are formed by an intersection of various discontinuities. It allows importing arbitrarily complex terrain surfaces which have been digitized beforehand using a topographic map. These surfaces are then extruded to a 3D solid which may be intersected by various sets of discontinuities. By combining all possible locations of all discontinuities potentially unstable blocks are determined. For each block, the factor of safety against sliding is computed using the limit equilibrium method. AutoBlock is an add-on to the popular program "AutoCAD" and exploits its possibilities and its power (e.g. for 3D-visualizations).
See also: Slope stability
A more rigorous approach to slope stability analysis is limit analysis. Unlike limit equilibrium analysis which makes ad-hoc though often reasonable assumptions, limit analysis is based on rigorous plasticity theory. This enables, among other things, the computation of upper and lower bounds on the true factor of safety.
Programs based on limit analysis include:
- OptumG2 (2014-) General purpose software for geotechnical applications including slope stability.
- LimitState:GEO (2008-) General purpose geotechnical software application based on Discontinuity layout optimization for plane strain problems including slope stability.
- GEO5 Slope Stability (1989-) Program is used to perform slope stability analysis of embankments, earth cuts, anchored retaining structures and MSE walls.
Stereographic and kinematic analysis
See also: Stereonet § Geology
Kinematic analysis examines which modes of failure can possibly occur in the rock mass. Analysis requires the detailed evaluation of rock mass structure and the geometry of existing discontinuities contributing to block instability.Stereographic representation (stereonets) of the planes and lines is used. Stereonets are useful for analyzing discontinuous rock blocks. Program DIPS allows for visualization structural data using stereonets, determination of the kinematic feasibility of rock mass and statistical analysis of the discontinuity properties.
Rock slope stability analysis may design protective measures near or around structures endangered by the falling blocks. Rockfall simulators determine travel paths and trajectories of unstable blocks separated from a rock slope face. Analytical solution method described by Hungr & Evans assumes rock block as a point with mass and velocity moving on a ballistic trajectory with regard to potential contact with slope surface. Calculation requires two restitution coefficients that depend on fragment shape, slope surface roughness, momentum and deformational properties and on the chance of certain conditions in a given impact.