# Slope stability analysis

This article may be too technical for most readers to understand. (December 2012) |

**Slope stability analysis** is a static or dynamic, analytical or empirical method to evaluate the stability of earth and rock-fill dams, embankments, excavated slopes, and natural slopes in soil and rock. Slope stability refers to the condition of inclined soil or rock slopes to withstand or undergo movement. The stability condition of slopes is a subject of study and research in soil mechanics, geotechnical engineering and engineering geology. Analyses are generally aimed at understanding the causes of an occurred slope failure, or the factors that can potentially trigger a slope movement, resulting in a landslide, as well as at preventing the initiation of such movement, slowing it down or arresting it through mitigation countermeasures.

The stability of a slope is essentially controlled by the ratio between the available shear strength and the acting shear stress, which can be expressed in terms of a safety factor if these quantities are integrated over a potential (or actual) sliding surface. A slope can be globally stable if the safety factor, computed along any potential sliding surface running from the top of the slope to its toe, is always larger than 1. The smallest value of the safety factor will be taken as representing the global stability condition of the slope. Similarly, a slope can be locally stable if a safety factor larger than 1 is computed along any potential sliding surface running through a limited portion of the slope (for instance only within its toe). Values of the global or local safety factors close to 1 (typically comprised between 1 and 1.3, depending on regulations) indicate marginally stable slopes that require attention, monitoring and/or an engineering intervention (slope stabilization) to increase the safety factor and reduce the probability of a slope movement.

A previously stable slope can be affected by a number of predisposing factors or processes that make the safety factor decrease - either by increasing the shear stress or by decreasing the shear strength - and can ultimately result in slope failure. Factors that can trigger slope failure include hydrologic events (such as intense or prolonged rainfall, rapid snowmelt, progressive soil saturation, increase of water pressure within the slope), earthquakes (including aftershocks), internal erosion (piping), surface or toe erosion, artificial slope loading (for instance due to the construction of a building), slope cutting (for instance to make space for roadways, railways or buildings), or slope flooding (for instance by filling an artificial lake after damming a river).

## Examples[edit]

As seen in Figure 1, earthen slopes can develop a cut-spherical weakness area. The probability of this happening can be calculated in advance using a simple 2-D circular analysis package.^{[1]} A primary difficulty with analysis is locating the most-probable slip plane for any given situation.^{[2]} Many landslides have only been analyzed after the fact. More recently slope stability radar technology has been employed, particularly in the mining industry, to gather real time data and assist in determining the likelihood of slope failure.

Real life failures in naturally deposited mixed soils are not necessarily circular, but prior to computers, it was far easier to analyze such a simplified geometry. Nevertheless, failures in 'pure' clay can be quite close to circular. Such slips often occur after a period of heavy rain, when the pore water pressure at the slip surface increases, reducing the effective normal stress and thus diminishing the restraining friction along the slip line. This is combined with increased soil weight due to the added groundwater. A 'shrinkage' crack (formed during prior dry weather) at the top of the slip may also fill with rain water, pushing the slip forward. At the other extreme, slab-shaped slips on hillsides can remove a layer of soil from the top of the underlying bedrock. Again, this is usually initiated by heavy rain, sometimes combined with increased loading from new buildings or removal of support at the toe (resulting from road widening or other construction work). Stability can thus be significantly improved by installing drainage paths to reduce the destabilizing forces. Once the slip has occurred, however, a weakness along the slip circle remains, which may then recur at the next monsoon.

Slope stability issues can be seen with almost any walk down a ravine in an urban setting. An example is shown in Figure 3, where a river is eroding the toe of a slope, and there is a swimming pool near the top of the slope. If the toe is eroded too far, or the swimming pool begins to leak, the forces driving a slope failure will exceed those resisting failure, and a landslide will develop, possibly quite suddenly.

## Measuring the Angle of Repose[edit]

The angle of repose is defined as the steepest angle of granular unconfined material measured from the horizontal plane on which the granular material can be heaped on without collapsing, ranging between 0-90°.^{[3]} For granular material, the angle of repose is the main factor that influences the slope's stability under different conditions in relation to the cohesiveness/friction of the material, the grain size, and the particle shape.^{[4]}

### Theoretical Measurement[edit]

A simple free body diagram can be used to understand the relationship between the angle of repose and the stability of the material on the slope. For the heaped material to collapse, the frictional forces must be equivalent to the horizontal component of the gravitational force , where is the mass of the material, is the gravitational acceleration and is the slope angle:

The frictional force is equivalent to the multiplication product of the coefficient of static friction and the Normal Force or :

Where is the angle of repose, or the angle at which the slope fails under regular conditions, and is the coefficient of static friction of the material on the slope.

### Experimental Measurement[edit]

*Tilting Box Method*[edit]

This method is particularly well suited for finer grained material < 10mm in diameter with relatively no cohesion. The material is placed in a base of a box, which is gradually tilted at a rate of 18°/min. The angle of repose is then measured to be the angle at which the material begins to slide.^{[3]}

*Fixed Funnel Method*[edit]

In this method, the material is poured down a funnel from a certain height into a horizontal base. The material is then allowed to heap until the pile reaches a predetermined height and width. The angle of repose is then measured by observing the height and radius of the cone and applying the arctangent rule.^{[3]}

## Angle of Repose and Slope Stability[edit]

The angle of repose and the stability of a slope are impacted by climatic and non-climatic factors.

*Water Content*[edit]

Water content is an important parameter that could change the angle of repose. reportedly, a higher water content can stabilize a slope and increase the angle of repose.^{[3]} However, water saturation can result in a decrease in the slope's stability since it acts as a lubricant and creates a detachment where mass wasting can occur.^{[5]}

Water content is dependent on soil properties such as grain size, which can impact infiltration rate, runoff, and water retention. Generally, finer-grained soils rich in clay and silt retain more water than coarser sandy soils. This effect is mainly due to capillary action, where the adhesive forces between the fluid, particle, and the cohesive forces of the fluid itself counteract gravitational pull. Therefore, smaller grain size results in a smaller surface area on which gravitational forces can act. Smaller surface area also leads to more capillary action, more water retention, more infiltration, and less runoff.^{[6]}

*Vegetation*[edit]

The presence of vegetation does not directly impact the angle of repose, but it acts as a stabilizing factor in a hillslope, where the tree roots anchor into deeper soil layers and form a fiber‐reinforced soil composite with a higher shear resistance (mechanical cohesion).^{[7]}

*Roundness of Grains*[edit]

The shape of the grain can have an impact on the angle of repose and the stability of the slope. The more rounded the grain is, the lower the angle of repose. A decrease in roundness, or an increase in angularity, results in interlocking via particle contact. This linear relationship between the angle of repose and the roundness of grain can also be used as a predictor of the angle of repose if the roundness of the grain is measured.^{[3]}

## Applications of Angle of Repose in Science and Engineering[edit]

The angle of repose is related to the shear strength of geologic materials, which is relevant in construction and engineering contexts.^{[8]} For granular materials, the size and shape of grains can impact angle of repose significantly. As the roundness of materials increases, the angle of repose decreases since there is less friction between the soil grains.^{[9]}

When the angle of repose is exceeded, mass wasting and rockfall can occur. It is important for many civil and geotechnical engineers to know the angle of repose to avoid structural and natural disasters. As a result, the application of retaining walls can help to retain soil so that the angle of repose is not exceeded.^{[3]}

## Slope stabilization[edit]

Since the stability of the slope can be impacted by external events such as precipitation, an important concern in civil/geotechnical engineering is the stabilization of slopes.

### Application of Vegetation[edit]

The application of vegetation to increase the slope stability against erosion and landslide is a form of bioengineering that is widely used in areas where the landslide depth is shallow. Vegetation increases the stability of the slope mechanically, by reinforcing the soils through plant roots, which stabilize the upper part of the soil. Vegetation also stabilizes the slope via hydrologic processes, by the reduction of soil moisture content through the interception of precipitation and transpiration. This results in a drier soil that is less susceptible to mass wasting.^{[10]}

Stability of slopes can also be improved by:

- Flattening of slopes results in reduction in weight which makes the slope more stable
- Soil stabilization
- Providing lateral supports by piles or retaining walls
- Grouting or cement injections into special places
- Consolidation by surcharging or electro osmosis increases the stability of slope

## Analysis methods[edit]

If the forces available to resist movement are greater than the forces driving movement, the slope is considered stable. A factor of safety is calculated by dividing the forces resisting movement by the forces driving movement. In earthquake-prone areas, the analysis is typically run for static conditions and pseudo-static conditions, where the seismic forces from an earthquake are assumed to add static loads to the analysis.

**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.^{[11]}^{[12]} Slope stability is the resistance of inclined surface to failure by sliding or collapsing.^{[13]} 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.^{[11]}^{[12]}

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.^{[14]}^{[15]} 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.^{[16]}
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.^{[17]}

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).^{[11]} 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).^{[18]}^{[19]}

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[edit]

Conventional methods of slope stability analysis can be divided into three groups: kinematic analysis, limit equilibrium analysis, and rock fall simulators.^{[18]}
Most slope stability analysis computer programs are based on the limit equilibrium concept for a two- or three-dimensional model.^{[20]}^{[21]} 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.^{[22]} In rock slope engineering, methods may be highly significant to simple block failure along distinct discontinuities.^{[18]} 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.^{[22]} 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.^{[21]}^{[23]} 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.^{[22]}^{[24]}

The location of the interface is typically unknown but can be found using numerical optimization methods.^{[25]} 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[edit]

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.^{[23]}^{[24]} Some of these approaches are discussed below.

#### Swedish Slip Circle Method of Analysis[edit]

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[edit]

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.^{[26]}

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[edit]

The Modified Bishop's method^{[27]} 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

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[edit]

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[edit]

Spencer's Method of analysis^{[28]} requires a computer program capable of cyclic algorithms, but makes slope stability analysis easier. Spencer's algorithm satisfies all equilibria (horizontal, vertical and driving moment) on each slice. The method allows for unconstrained slip plains and can therefore determine the factor of safety along any slip surface. The rigid equilibrium and unconstrained slip surface result in more precise safety factors than, for example, Bishop's Method or the Ordinary Method of Slices.^{[28]}

#### Sarma method[edit]

The Sarma method,^{[29]} 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.

### Comparisons[edit]

The assumptions made by a number of limit equilibrium methods are listed in the table below.^{[30]}

Method | Assumption |
---|---|

Ordinary method of slices | Interslice forces are neglected |

Bishop's simplified/modified ^{[27]} |
Resultant interslice forces are horizontal. There are no interslice shear forces. |

Janbu's simplified^{[31]} |
Resultant interslice forces are horizontal. An empirical correction factor is used to account for interslice shear forces. |

Janbu's generalized^{[31]} |
An assumed line of thrust is used to define the location of the interslice normal force. |

Spencer ^{[28]} |
The resultant interslice forces have constant slope throughout the sliding mass. The line of thrust is a degree of freedom. |

Chugh^{[32]} |
Same as Spencer's method but with a constant acceleration force on each slice. |

Morgenstern-Price^{[33]} |
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) ^{[26]} |
Similar to Morgenstern-Price. |

Corps of Engineers ^{[34]} |
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 ^{[35]} |
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 ^{[29]} |
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.^{[30]}

Method | Force balance (vertical) | Force balance (horizontal) | Moment balance |
---|---|---|---|

Ordinary MS | Yes | No | Yes |

Bishop's simplified | Yes | No | Yes |

Janbu's simplified | Yes | Yes | No |

Janbu's generalized | Yes | Yes | Used to compute interslice shear forces |

Spencer | Yes | Yes | Yes |

Chugh | Yes | Yes | Yes |

Morgenstern-Price | Yes | Yes | Yes |

Fredlund-Krahn | Yes | Yes | Yes |

Corps of Engineers | Yes | Yes | No |

Lowe and Karafiath | Yes | Yes | No |

Sarma | Yes | Yes | Yes |

### Rock slope stability analysis[edit]

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^{[36]}^{[37]}**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.)^{[37]}**Wedge failure**-> three-dimensional analysis enables modelling of the wedge sliding on two planes in a direction along the line of intersection^{[37]}^{[38]}**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^{[39]}^{[40]}

## Limit analysis[edit]

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 (also includes elastoplasticity, seepage, consolidation, staged construction, tunneling, and other relevant geotechnical analysis types).
- LimitState:GEO (2008-) General purpose geotechnical software application based on Discontinuity layout optimization for plane strain problems including slope stability.

## Stereographic and kinematic analysis[edit]

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.^{[41]}^{[42]} Stereographic representation (stereonets) of the planes and lines is used.^{[43]} Stereonets are useful for analyzing discontinuous rock blocks.^{[44]} Program **DIPS**^{[45]} allows for visualization structural data using stereonets, determination of the kinematic feasibility of rock mass and statistical analysis of the discontinuity properties.^{[41]}^{[45]}

## Rockfall simulators[edit]

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.^{[46]} Analytical solution method described by Hungr & Evans^{[47]} 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.^{[48]}

## Numerical methods of analysis[edit]

Numerical modelling techniques provide an approximate solution to problems which otherwise cannot be solved by conventional methods, e.g. complex geometry, material anisotropy, non-linear behavior, in situ stresses. Numerical analysis allows for material deformation and failure, modelling of pore pressures, creep deformation, dynamic loading, assessing effects of parameter variations etc. However, numerical modelling is restricted by some limitations. For example, input parameters are not usually measured and availability of these data is generally poor. User also should be aware of boundary effects, meshing errors, hardware memory and time restrictions. Numerical methods used for slope stability analysis can be divided into three main groups: continuum, discontinuum and hybrid modelling.^{[49]}

### Continuum modelling[edit]

Modelling of the continuum is suitable for the analysis of soil slopes, massive intact rock or heavily jointed rock masses. This approach includes the *finite-difference* and *finite element* methods that discretize the whole mass to finite number of elements with the help of generated mesh (Fig. 3). In *finite-difference* method (FDM) differential equilibrium equations (i.e. strain-displacement and stress-strain relations) are solved. *finite element* method (FEM) uses the approximations to the connectivity of elements, continuity of displacements and stresses between elements. ^{[50]} Most of numerical codes allows modelling of discrete fractures, e.g. bedding planes, faults. Several constitutive models are usually available, e.g. elasticity, elasto-plasticity, strain-softening, elasto-viscoplasticity etc.^{[49]}

### Discontinuum modelling[edit]

Discontinuum approach is useful for rock slopes controlled by discontinuity behaviour. Rock mass is considered as an aggregation of distinct, interacting blocks subjected to external loads and assumed to undergo motion with time. This methodology is collectively called the *discrete-element* method (DEM). Discontinuum modelling allows for sliding between the blocks or particles. The DEM is based on solution of dynamic equation of equilibrium for each block repeatedly until the boundary conditions and laws of contact and motion are satisfied. Discontinuum modelling belongs to the most commonly applied numerical approach to rock slope analysis and following variations of the DEM exist:^{[49]}

- distinct-element method
- discontinuous deformation analysis (DDA)
- particle flow codes

The *distinct-element* approach describes mechanical behaviour of both, the discontinuities and the solid material. This methodology is based on a force-displacement law (specifying the interaction between the deformable rock blocks) and a law of motion (determining displacements caused in the blocks by out-of-balance forces). Joints are treated as [boundary conditions. Deformable blocks are discretized into internal constant-strain elements.^{[49]}

Discontinuum program **UDEC**^{[51]} (Universal distinct element code) is suitable for high jointed rock slopes subjected to static or dynamic loading. Two-dimensional analysis of translational failure mechanism allows for simulating large displacements, modelling deformation or material yielding.^{[51]} Three-dimensional discontinuum code **3DEC**^{[52]} contains modelling of multiple intersecting discontinuities and therefore it is suitable for analysis of wedge instabilities or influence of rock support (e.g. rockbolts, cables).^{[49]}

In *discontinuous deformation analysis* (DDA) displacements are unknowns and equilibrium equations are then solved analogous to *finite element* method. Each unit of *finite element* type mesh represents an isolated block bounded by discontinuities. Advantage of this methodology is possibility to model large deformations, rigid body movements, coupling or failure states between rock blocks.^{[49]}

Discontinuous rock mass can be modelled with the help of *distinct-element* methodology in the form of *particle flow* code, e.g.
program **PFC2D/3D**.^{[53]}^{[54]} Spherical particles interact through frictional sliding contacts. Simulation of joint bounded blocks may be realized through specified bond strengths. Law of motion is repeatedly applied to each particle and force-displacement law to each contact. *Particle flow* methodology enables modelling of granular flow, fracture of intact rock, transitional block movements, dynamic response to blasting or seismicity, deformation between particles caused by shear or tensile forces. These codes also allow to model subsequent failure processes of rock slope, e.g. simulation of rock^{[49]}

### Hybrid/coupled modelling[edit]

Hybrid codes involve the coupling of various methodologies to maximize their key advantages, e.g. *limit equilibrium* analysis combined with *finite element* groundwater flow and stress analysis ; coupled *particle flow* and *finite-difference* analyses. Hybrid techniques allows investigation of piping slope failures and the influence of high groundwater pressures on the failure of weak rock slope. Coupled *finite-*/*distinct-element* codes provide for the modelling of both intact rock behavior and the development and behavior of fractures.^{[49]}

^{[55]}

## Rock mass classification[edit]

Various rock mass classification systems exist for the design of slopes and to assess the stability of slopes. The systems are based on empirical relations between rock mass parameters and various slope parameters such as height and slope dip.

The Q-slope method for rock slope engineering and rock mass classification developed by Barton and Bar^{[56]} expresses the quality of the rock mass for assessing slope stability using the Q-slope value, from which long-term stable, reinforcement-free slope angles can be derived.

## Probability classification[edit]

The **slope stability probability classification** (SSPC)^{[57]}^{[58]} system is a rock mass classification system for slope engineering and slope stability assessment. The system is a three-step classification: *‘exposure’*, *‘reference’*, and *‘slope’* rock mass classification with conversion factors between the three steps depending on existing and future weathering and damage due to method of excavation. The stability of a slope is expressed as probability for different failure mechanisms.

A rock mass is classified following a standardized set of criteria in one or more exposures (*‘exposure’* classification). These values are converted per exposure to a *‘reference’* rock mass by compensating for the degree of weathering in the exposure and the method of excavation that was used to make the exposure, i.e. the *‘reference’* rock mass values are not influenced by local influences such as weathering and method of excavation. A new slope can then be designed in the *‘reference’* rock mass with compensation for the damage due to the method of excavation to be used for making the new slope and compensation for deterioration of the rock mass due to future weathering (the *‘slope’* rock mass). If the stability of an already existing slope is assessed the *‘exposure’* and *‘slope’* rock mass values are the same.

The failure mechanisms are divided in *orientation dependent* and *orientation independent*. Orientation dependent failure mechanisms depend on the orientation of the slope with respect to the orientation of the discontinuities in the rock mass, i.e. sliding (plane and wedge sliding) and toppling failure. Orientation independent relates to the possibility that a slope fails independently from its orientation, e.g. circular failure completely through newly formed discontinuities in intact rock blocks, or failing partially following existing discontinuities and partially new discontinuities.

In addition the shear strength along a discontinuity ('sliding criterion')^{[57]}^{[58]}^{[59]} and 'rock mass cohesion' and 'rock mass friction' can be determined. The system has been used directly or modified in various geology and climate environments throughout the world.^{[60]}^{[61]}^{[62]} The system has been modified for slope stability assessment in open pit coal mining.^{[63]}

## See also[edit]

- Rock mass rating
- SMR classification
- Plaxis
- Angle of repose
- Retaining wall
- Discontinuous Deformation Analysis
- Discontinuity layout optimization
- Discrete element method
- Finite difference method
- Finite element limit analysis
- Finite element method
- Stereonet
- Q-slope
- Discontinuity layout optimization
- Mohr-Coulomb theory
- SMR classification

## References[edit]

**^**"Slope Stability Calculator". Retrieved 2006-12-14.**^**Chugh, Ashok K. (2002). "A method for locating critical slip surfaces in slope stability analysis: Discussion".*Canadian Geotechnical Journal*.**39**(3): 765–770. doi:10.1139/t02-042.- ^
^{a}^{b}^{c}^{d}^{e}^{f}Beakawi Al-Hashemi, Hamzah M.; Baghabra Al-Amoudi, Omar S. (May 2018). "A review on the angle of repose of granular materials".*Powder Technology*.**330**: 397–417. doi:10.1016/j.powtec.2018.02.003. **^**Fu, Jian-Jun; Chen, Cheng; Ferellec, Jean-Francois; Yang, Juan (2020-08-17). "Effect of Particle Shape on Repose Angle Based on Hopper Flow Test and Discrete Element Method".*Advances in Civil Engineering*.**2020**: 1–10. doi:10.1155/2020/8811063. ISSN 1687-8086.**^**Prof A Balasubramanian (2011). "MASS-WASTING". doi:10.13140/RG.2.2.10405.50407. Cite journal requires`|journal=`

(help)**^**Kozicki, J.; Donzé, F.V. (2009-10-09). "YADE‐OPEN DEM: an open‐source software using a discrete element method to simulate granular material".*Engineering Computations*.**26**(7): 786–805. doi:10.1108/02644400910985170. ISSN 0264-4401.**^**Kim, John H.; Fourcaud, Thierry; Jourdan, Christophe; Maeght, Jean-Luc; Mao, Zhun; Metayer, James; Meylan, Louise; Pierret, Alain; Rapidel, Bruno; Roupsard, Olivier; de Rouw, Anneke (2017-05-28). "Vegetation as a driver of temporal variations in slope stability: The impact of hydrological processes: Variable Stability of Vegetated Slopes".*Geophysical Research Letters*.**44**(10): 4897–4907. doi:10.1002/2017GL073174.**^**Kim, Donghwi; Nam, Boo Hyun; Youn, Heejung (December 2018). "Effect of clay content on the shear strength of clay–sand mixture".*International Journal of Geo-Engineering*.**9**(1): 19. doi:10.1186/s40703-018-0087-x. ISSN 2092-9196.**^**Santamarina, J. Carlos (2003-01-13). "Soil Behavior at the Microscale: Particle Forces".*Soil Behavior and Soft Ground Construction*. Reston, VA: American Society of Civil Engineers: 25–56. doi:10.1061/40659(2003)2. ISBN 978-0-7844-0659-5.**^**Mulyono, A; Subardja, A; Ekasari, I; Lailati, M; Sudirja, R; Ningrum, W (February 2018). "The Hydromechanics of Vegetation for Slope Stabilization".*IOP Conference Series: Earth and Environmental Science*.**118**: 012038. doi:10.1088/1755-1315/118/1/012038. ISSN 1755-1307.- ^
^{a}^{b}^{c}Eberhardt 2003, p. 4 - ^
^{a}^{b}Abramson 2002, p. 2 **^**Kliche 1999, p. 2**^**USArmyCorps 2003, pp. 1–2**^**Abramson 2002, p. 1**^**Beale, Geoff; Read, John, eds. (2014).*Guidelines for Evaluating Water in Pit Slope Stability*. CSIRO Publishing. ISBN 9780643108356.**^**Stead 2001, p. 615- ^
^{a}^{b}^{c}Eberhardt 2003, p. 6 **^**Cardenas, IC (2019). "On the use of Bayesian networks as a meta-modelling approach to analyse uncertainties in slope stability analysis".*Georisk: Assessment and Management of Risk for Engineered Systems and Geohazards*.**13**(1): 53–65. doi:10.1080/17499518.2018.1498524. S2CID 216590427.**^**Abramson 2002, p. 329- ^
^{a}^{b}Abramson 2002, p. 363 - ^
^{a}^{b}^{c}USArmyCorps 2003, p. 2 - ^
^{a}^{b}Zhu 2003, pp. 377–395 - ^
^{a}^{b}Abramson 2002, pp. 363–367 **^**USArmyCorps 2003, p. 5- ^
^{a}^{b}Fredlund, DG; Krahn, J (1977), "Comparison of slope stability methods of analysis",*Canadian Geotechnical Journal*,**14**(3): 429–439, doi:10.1139/t77-045 - ^
^{a}^{b}Bishop, A. W. (1955). "The use of the Slip Circle in the Stability Analysis of Slopes".*Géotechnique*.**5**: 7–17. doi:10.1680/geot.1955.5.1.7. - ^
^{a}^{b}^{c}Spencer, E. (1967). "A Method of analysis of the Stability of Embankments Assuming Parallel Inter-Slice Forces".*Géotechnique*.**17**: 11–26. doi:10.1680/geot.1967.17.1.11. - ^
^{a}^{b}Sarma, S. K. (1975). "Seismic stability of earth dams and embankments".*Géotechnique*.**25**(4): 743–761. doi:10.1680/geot.1975.25.4.743. - ^
^{a}^{b}Fredlund, D.G. (1984), "Analytical methods for slope stability analysis" (PDF),*Proceedings of the Fourth International Symposium on Landslides, State-of-the-Art*: 229–250 - ^
^{a}^{b}Janbu, Nilmar (1973), RC Hirschfeld; SJ Poulos (eds.), "Slope stability computations",*In Embankment-dam Engineering*, Jon Wiley and Sons Inc., NY: 40P **^**Chugh, Ashok K (1982), "Slope stability analysis for earthquakes",*International Journal for Numerical and Analytical Methods in Geomechanics*,**6**(3): 307–322, Bibcode:1982IJNAM...6..307C, doi:10.1002/nag.1610060304**^**Morgenstern, N. R.; Price, V. Eo (1965), "The analysis of the stability of general slip surfaces",*Geotechnique*,**15**(1): 79–93, doi:10.1680/geot.1965.15.1.79**^**"Slope Stability" (PDF). US Army Corps of Engineers. Retrieved 15 April 2015.**^**Lowe, John; Karafiath, Leslie (1960), "Stability of earth dams upon drawdown",*In Proc. 1st. Pan American Conference on Soil Mechanics and Foundation Engineering, México*,**2**: 537–552**^**Kliche 1999, pp. 125–137- ^
^{a}^{b}^{c}Kovari 1978, pp. 103–124 **^**Kliche 1999, pp. 153–169**^**Kliche 1999, p. 15**^**Kliche 1999, pp. 139–152- ^
^{a}^{b}Eberhardt 2003, p. 7 **^**Kliche 1999, p. 111**^**Kliche 1999, pp. 111–123**^**Kliche 1999, pp. 43–65- ^
^{a}^{b}"DIPS – Graphical and Statistical Analysis of Orientation Data",*Rocscience.com*, Toronto, Canada: Rocscience, retrieved 21 July 2009 **^**Bourrier, Franck; Hungr, Oldrich (2013-03-06), Lambert, Stéphane; Nicot, François (eds.), "Rockfall Dynamics: A Critical Review of Collision and Rebound Models",*Rockfall Engineering*, Hoboken, NJ, USA: John Wiley & Sons, Inc., pp. 175–209, doi:10.1002/9781118601532.ch6, ISBN 978-1-118-60153-2, retrieved 2021-01-18**^**Hungr 1988, pp. 685–690**^**Eberhardt 2003, pp. 15–17- ^
^{a}^{b}^{c}^{d}^{e}^{f}^{g}^{h}Eberhardt 2003, pp. 17–38 **^**Kaczmarek, Łukasz D.; Popielski, Paweł (1 January 2019). "Selected components of geological structures and numerical modelling of slope stability".*Open Geosciences*.**11**(1): 208–218. doi:10.1515/geo-2019-0017.- ^
^{a}^{b}"UDEC - Universal Distinct Element Code",*Itascacg.com*, Minneapolis, USA: Itasca, retrieved 27 July 2009 **^**"3DEC - Three Dimensional Distinct Element Code",*Itascacg.com*, Minneapolis, USA: Itasca, retrieved 27 July 2009**^**"PFC2D - Particle Flow Code in Two Dimensions",*Itascacg.com*, Minneapolis, USA: Itasca, retrieved 27 July 2009**^**"PFC3D - Particle Flow Code in Three Dimensions",*Itascacg.com*, Minneapolis, USA: Itasca, retrieved 27 July 2009**^**Van der Meij, R. (25 May 2010). "A Genetic Algorithm for Solving Slope Stability Problems: From Bishop to a Free Slip Plane".*Numerical Methods in Geotechnical Engineering*.**4**: 345–350. ISBN 9780203842362.**^**Bar, N.; Barton, N.R. (2017). "The Q-slope Method for Rock Slope Engineering".*Rock Mechanics & Rock Engineering*, Vol 50, Springer, Vienna, https://doi.org/10.1007/s00603-017-1305-0.- ^
^{a}^{b}Hack, R. (1996 and 1998).*Slope Stability Probability Classification (SSPC)*(PDF). ITC publication 43. Technical University Delft & Twente University - International Institute for Aerospace Survey and Earth Sciences (ITC Enschede), Netherlands. p. 258. ISBN 978-90-6164-154-4. Check date values in:`|date=`

(help) - ^
^{a}^{b}Hack, R.; Price, D.; Rengers, N. (2003). "A new approach to rock slope stability – a probability classification (SSPC)".*Bulletin of Engineering Geology and the Environment*.**62**(2): 167–184. doi:10.1007/s10064-002-0155-4. S2CID 140693335. **^**Andrade, P.S.; Saraiva, A.A. (2008). "Estimating the joint roughness coefficient of discontinuities found in metamorphic rocks" (PDF).*Bulletin of Engineering Geology and the Environment*.**67**(3, number 3): 425–434. doi:10.1007/s10064-008-0151-4. hdl:10316/7611. S2CID 129119508.**^**Filipello, A.; Giuliani, A.; Mandrone, G. (2010). "Rock Slopes Failure Susceptibility Analysis: From Remote Sensing Measurements to Geographic Information System Raster Modules".*American Journal of Environmental Sciences*.**6**(6): 489–494. doi:10.3844/ajessp.2010.489.494.**^**Hailemariam, G.T.; Schneider, J.F. (May 2–7, 2010). "Rock Mass Classification of Karstic Terrain in the Reservoir Slopes of Tekeze Hydropower Project" (PDF).*EGU General Assembly 2010*. EGU2010-831, 2010.**12**. Vienna, Austria. p. 831.**^**Dhakal, S.; Upreti, B.N.; Yoshida, M.; Bhattarai, T.N.; Rai, S.M.; Gajurel, A.P.; Ulak, P.D.; Dahal, R.K. (2005). "Application of the SSPC system in some of the selected slopes along the trekking route from Jomsom to Kagbeni, central-west Nepal". In Yoshida, M.; Upreti, B.N.; Bhattarai, T.N.; Dhakal, S. (eds.).*Natural disaster mitigation and issues on technology transfer in South and Southeast Asia; proceedings of the JICA Regional Seminar*. Kathmandu, Nepal: Department of Geology, Tri-Chandra Campus, Tribhuvan University, Kathmandu, Nepal. pp. 79–82.**^**Lindsay, P.; Campbellc, R.N.; Fergussonc, D.A.; Gillarda, G.R.; Moore, T.A. (2001). "Slope stability probability classification, Waikato Coal Measures, New Zealand".*International Journal of Coal Geology*.**45**(2–3): 127–145. doi:10.1016/S0166-5162(00)00028-8.

## Further reading[edit]

- Devoto, S.; Castelli, E. (September 2007). "Slope stability in an old limestone quarry interested by a tourist project".
*15th Meeting of the Association of European Geological Societies: Georesources Policy, Management, Environment*. Tallinn. - Douw, W. (2009).
*Entwicklung einer Anordnung zur Nutzung von Massenschwerebewegungen beim Quarzitabbau im Rheinischen Schiefergebirge*. Hackenheim, Germany: ConchBooks. p. 358. ISBN 978-3-939767-10-7. - Hack, H.R.G.K. (25–28 November 2002). "An evaluation of slope stability classification. Keynote Lecture.". In Dinis da Gama, C.; Ribeira e Sousa, L. (eds.).
*Proc. ISRM EUROCK’2002*. Funchal, Madeira, Portugal: Sociedade Portuguesa de Geotecnia, Lisboa, Portugal. pp. 3–32. ISBN 972-98781-2-9. - Liu, Y.-C.; Chen, C.-S. (2005). "A new approach for application of rock mass classification on rock slope stability assessment".
*Engineering Geology*.**89**(1–2): 129–143. doi:10.1016/j.enggeo.2006.09.017. - Pantelidis, L. (2009). "Rock slope stability assessment through rock mass classification systems".
*International Journal of Rock Mechanics and Mining Sciences*.**46**(2, number 2): 315–325. doi:10.1016/j.ijrmms.2008.06.003. - Rupke, J.; Huisman, M.; Kruse, H.M.G. (2007). "Stability of man-made slopes".
*Engineering Geology*.**91**(1): 16–24. doi:10.1016/j.enggeo.2006.12.009. - Singh, B.; Goel, R.K. (2002).
*Software for engineering control of landslide and tunnelling hazards*.**1**. Taylor & Francis. p. 358. ISBN 978-90-5809-360-8.

- Coduto, Donald P. (1998).
*Geotechnical Engineering: Principles and Practices.*Prentice-Hall. ISBN 0-13-576380-0 - Fredlund, D. G., H. Rahardjo, M. D. Fredlund (2014).
*Unsaturated Soil Mechanics in Engineering Practice.*Wiley-Interscience. ISBN 978-1118133590 - Kliche, Charles A. (1999),
*Rock Slope Stability*, Colorado, USA: Society for Mining, Metallurgy, and Exploration, ISBN 0-87335-171-1CS1 maint: ref duplicates default (link) - Eberhardt, Erik (2003),
*Rock Slope Stability Analysis - Utilization of Advanced Numerical Techniques*(PDF), Vancouver, Canada: Earth and Ocean Sciences, University of British ColumbiaCS1 maint: ref duplicates default (link) - US Army Corps of Engineers (2003),
*Engineering and Design - Slope Stability*(PDF), Washington DC, USA: US Army Corps of Engineers - Stead, Doug; Eberhardt, E.; Coggan, J.; Benko, B. (2001), M. Kühne; H.H. Einstein; E. Krauter; H. Klapperich; R. Pöttler (eds.),
*Advanced numerical techniques in rock slope stability analysis - Applications and limitations*(PDF), Davos, Switzerland: Verlag Glückauf GmbH, pp. 615–624 Unknown parameter`|conference=`

ignored (help) - Abramson, Lee W.; Lee, Thomas S.; Sharma, Sunil; Boyce, Glenn M. (2002),
*Slope Stability and Stabilization Methods*(2nd ed.), New York, USA: John Wiley & Sons, ISBN 0-471-38493-3 - Zhu, D.Y.; Lee, C.F.; Jiang, H.D. (2003), "Generalised framework of limit equilibrium methods for slope stability analysis",
*Geotechnique*, Telford, London, Great Britain,**53**(4): 377–395, doi:10.1680/geot.2003.53.4.377, hdl:10722/71758, ISSN 0016-8505 - Kovári, Kalman; Fritz, P. (1978),
*Slope Stability with Plane, Wedge and Polygonal Sliding Surfaces*, Rio de Janeiro, Brazil, pp. 103–124 Unknown parameter`|conference=`

ignored (help) - Yang, Xiao-Li; Li, L.; Yin, J.H. (2004), "Stability analysis of rock slopes with a modified Hoek-Brown failure criterion",
*International Journal for Numerical and Analytical Methods in Geomechanics*, Chichester, Great Britain: John Wiley & Sons,**28**(2): 181–190, Bibcode:2004IJNAM..28..181Y, doi:10.1002/nag.330, ISSN 0363-9061 - Barton, N.R.; Bandis, S.C. (1990), Barton, Nick (ed.),
*Review of predictive capabilities of JRC-JCS model in engineering practice*, Rotterdam: Balkema, pp. 603–610, ISBN 978-90-6191-109-8 Unknown parameter`|book-title=`

ignored (help); Unknown parameter`|conference=`

ignored (help) - Hungr, O.; Evans, S.G. (1988), Bonnard, C. (ed.),
*Engineering evaluation of fragmental rockfall hazards*, Rotterdam: Balkema, pp. 685–690 Unknown parameter`|book-title=`

ignored (help); Unknown parameter`|conference=`

ignored (help)