5.1 Slope Stability Analysis & Limit Equilibrium Methods

Introduction to Slope Stability Analysis

Slope stability analysis is a fundamental discipline in geotechnical engineering, essential for evaluating the stability of natural slopes, engineered embankments, earth dams, deep excavations, and soil cuts. The primary objective is to determine the Factor of Safety ($FS$) against shear failure along a potential slip surface.

The general definition of the Factor of Safety against slope failure is expressed as the ratio of the available shear strength ($\tau_f$) along the failure surface to the mobilized shear stress ($\tau_m$) required to maintain equilibrium:

$$FS = \frac{\tau_f}{\tau_m} = \frac{c' + \sigma'_n \tan\phi'}{\tau_m}$$

where:

• $c'$ is the effective cohesion intercept of the soil
• $\sigma'_n = (\sigma_n - u)$ is the effective normal stress on the failure plane
• $u$ is the pore water pressure
• $\phi'$ is the effective internal friction angle

Slope Failure Mechanisms & Surface Geometries

Depending on soil stratigraphy, strength anisotropy, and geological boundary conditions, slope failures exhibit distinct kinematic modes:

1. Circular Rotational Failures: Prevalent in homogeneous clay cuts, soft ground embankments, and uniform fills. The slip surface forms a continuous curved arc that can be approximated by a circular segment.
2. Planar (Translational) Failures: Common in stratified deposits where a weak soil layer or bedding plane runs parallel to the slope face. The failure mass slides along a flat plane.
3. Non-Circular Arbitrary Slip Surfaces: Occur in non-homogeneous, multi-layered soils with variable pore pressure regimes or structural weak seams.
4. Infinite Slopes: Shallow failure surfaces where the depth of the slip surface ($z$) is significantly smaller than the length of the slope ($L \gg z$).

Total Stress Analysis (TSA) vs. Effective Stress Analysis (ESA)

Geotechnical engineers must select the appropriate stress analysis framework based on drainage conditions and time frames:

• Total Stress Analysis (TSA / Short-Term / Undrained): Used for short-term evaluation immediately after construction or excavation in saturated cohesive soils. The shear strength is governed by the undrained shear strength ($s_u$ or $c_u$, with $\phi_u = 0$). Excess pore water pressure generated by construction load changes does not have time to dissipate.
• Effective Stress Analysis (ESA / Long-Term / Drained): Used for long-term stability analysis under steady-state seepage conditions, or rapid drawdown conditions in free-draining materials. Soil strength is defined by effective strength parameters ($c', \phi'$) and the equilibrium pore water pressure grid ($u$) derived from seepage flow nets.

Infinite Slope Mechanics & Mathematical Formulations

An infinite slope represents an idealized condition where soil properties, groundwater depth, and slope inclination are uniform over an infinite extent. This model is highly accurate for shallow translational slides parallel to the slope face.

Case 1: Dry or Fully Drained Cohesionless Soil ($c' = 0, u = 0$)

Consider an element of soil at depth $z$ on a slope inclined at angle $\beta$ to the horizontal. The weight of the soil block of unit width is $W = \gamma z \cos\beta$. Resolving $W$ into normal and shear components along the base plane yields:

$$\sigma_n = \gamma z \cos^2\beta$$ $$\tau_m = \gamma z \sin\beta \cos\beta$$

The available frictional strength is $\tau_f = \sigma_n \tan\phi' = \gamma z \cos^2\beta \tan\phi'$. The Factor of Safety is:

$$FS = \frac{\tau_f}{\tau_m} = \frac{\gamma z \cos^2\beta \tan\phi'}{\gamma z \sin\beta \cos\beta} = \frac{\tan\phi'}{\tan\beta}$$

Key Concept: A dry, cohesionless slope is stable ($FS \ge 1.0$) as long as the slope angle $\beta$ does not exceed the soil's effective friction angle $\phi'$. Depth $z$ and unit weight $\gamma$ do not affect stability.

Case 2: Cohesive-Frictional Soil without Seepage ($c' > 0, \phi' > 0, u = 0$)

Including effective cohesion $c'$ introduces a depth-dependent term:

$$FS = \frac{c' + \gamma z \cos^2\beta \tan\phi'}{\gamma z \sin\beta \cos\beta} = \frac{c'}{\gamma z \sin\beta \cos\beta} + \frac{\tan\phi'}{\tan\beta}$$

Here, cohesion increases the factor of safety, particularly at shallow depths ($z$). The critical depth $z_{cr}$ at which failure occurs ($FS = 1.0$) can be calculated directly by rearranging the equation.

Case 3: Cohesionless Soil with Parallel Seepage to Ground Surface ($c' = 0$)

When seepage flows parallel to the slope face with the phreatic surface at the ground surface, the pore water pressure at depth $z$ is $u = \gamma_w z \cos^2\beta$. The effective normal stress becomes:

$$\sigma'_n = \sigma_n - u = \gamma_{sat} z \cos^2\beta - \gamma_w z \cos^2\beta = (\gamma_{sat} - \gamma_w) z \cos^2\beta = \gamma' z \cos^2\beta$$

where $\gamma'$ is the submerged (buoyant) unit weight of the soil ($\gamma' = \gamma_{sat} - \gamma_w$). Substituting $\sigma'_n$ into the factor of safety equation yields:

$$FS = \frac{\gamma' z \cos^2\beta \tan\phi'}{\gamma_{sat} z \sin\beta \cos\beta} = \left( \frac{\gamma'}{\gamma_{sat}} \right) \frac{\tan\phi'}{\tan\beta}$$

Since $\gamma' \approx 0.5 \gamma_{sat}$ for typical soils (e.g., $\gamma_{sat} \approx 19-20\text{ kN/m}^3, \gamma' \approx 9-10\text{ kN/m}^3$), seepage parallel to the slope cuts the Factor of Safety in half compared to the dry condition!

Limit Equilibrium Method of Slices

For non-infinite slopes with curved or arbitrary failure surfaces, the sliding mass is discretized into $n$ vertical slices to account for varying soil properties, pore water pressures, and slope geometry along the slip surface.

Force Equilibrium on an Individual Slice

Each vertical slice $i$ of width $b_i$ experiences the following system of forces:

• Slice weight $W_i = \gamma_i b_i h_i$
• Base normal force $N_i$ and base shear force $T_i$
• Pore water pressure force at the base $U_i = u_i l_i$, where $l_i = b_i / \cos\alpha_i$
• Interslice normal forces $E_{L,i}$ and $E_{R,i}$ acting on the left and right vertical boundaries
• Interslice shear forces $X_{L,i}$ and $X_{R,i}$ acting along the vertical boundaries
• Pseudostatic seismic force $k_h W_i$ (if seismic loading is analyzed)

Method Comparison and Equilibrium Assumptions

Different limit equilibrium methods make specific simplifying assumptions regarding interslice forces to resolve the statically indeterminate system:

1. Ordinary Method of Slices (Fellenius / Swedish Circle)

The Ordinary Method of Slices assumes that interslice forces $E$ and $X$ are parallel to the base of each slice and cancel each other out ($E_L = E_R, X_L = X_R$). This simplifies base normal force resolution:

$$N'_i = N_i - U_i = W_i \cos\alpha_i - u_i l_i$$

The global factor of safety is derived solely from overall moment equilibrium about the center of the circular slip surface:

$$FS = \frac{\sum \left[ c' l_i + (W_i \cos\alpha_i - u_i l_i) \tan\phi' \right]}{\sum W_i \sin\alpha_i}$$

• Pros: Simple to calculate manually without iterative software routines.
• Cons: Neglects interslice forces; can be significantly conservative, underestimating $FS$ by 5% to 15% (or overestimating error in deep circles with high pore pressures).

2. Bishop's Simplified Method

Bishop's Simplified Method assumes interslice shear forces are zero ($X_L = X_R = 0$), but includes horizontal interslice forces ($E_L, E_R$). It satisfies vertical force equilibrium for each slice and overall moment equilibrium about the center of rotation.

The resulting equation for $FS$ is:

$$FS = \frac{\sum \frac{c' b_i + (W_i - u_i b_i) \tan\phi'}{m\alpha}}{\sum W_i \sin\alpha_i}$$

where the term $m_\alpha$ is defined as:

$$m\alpha = \cos\alpha_i \left( 1 + \frac{\tan\alpha_i \tan\phi'}{FS} \right)$$

• Pros: Highly accurate for circular failure surfaces (typically within 1-2% of rigorous methods).
• Cons: Requires iterative solution because $FS$ appears on both sides of the equation (inside $m_\alpha$).

3. Janbu's Simplified Method

Janbu's Simplified Method is designed for non-circular (composite or planar) slip surfaces. It assumes horizontal interslice forces ($X_i = 0$) and satisfies force equilibrium (horizontal and vertical) across all slices. It applies an empirical correction factor $f_0$ to account for interslice shear forces:

$$FS_{corrected} = f_0 \cdot FS_{uncorrected}$$

4. Spencer's Method & Morgenstern-Price Method

• Spencer's Method: Assumes interslice forces are inclined at a constant angle $\theta$ across all slices ($X/E = \tan\theta$). Solves for both $\theta$ and $FS$, fully satisfying both force equilibrium and moment equilibrium for arbitrary slip surfaces.
• Morgenstern-Price Method: Assumes interslice force inclination varies according to a prescribed function $f(x)$ scaled by $\lambda$ ($X/E = \lambda f(x)$). Satisfies both force and moment equilibrium and is considered the standard for modern 2D and 3D limit equilibrium software.

Comparison of Limit Equilibrium Slope Stability Methods

---

Taylor's Stability Charts & Mass Procedures

For homogeneous, purely cohesive slopes ($\phi_u = 0$) under total stress conditions, Taylor developed chart solutions based on the dimensionless Stability Number ($N_s$):

$$N_s = \frac{c_u}{\gamma H \cdot FS} \quad \implies \quad FS = \frac{c_u}{\gamma H N_s}$$

where:

• $c_u$ is the undrained shear strength
• $\gamma$ is the total unit weight of the soil
• $H$ is the height of the slope cut
• $N_s$ is Taylor's stability number, obtained from charts based on slope angle $\beta$ and depth factor $D_d = D/H$ (where $D$ is the depth from top of slope to a hard rock stratum).

---

Worked Engineering Calculation: Slope Stability Analysis

Problem Statement

A permanent road cut slope is to be excavated in a saturated clay deposit at an angle of $\beta = 30^\circ$ to a total vertical height $H = 10.0\text{ m}$. Field and laboratory testing indicate:

• Undrained shear strength $c_u = 45.0\text{ kPa}$
• Undrained friction angle $\phi_u = 0^\circ$
• Total saturated unit weight $\gamma_{sat} = 19.0\text{ kN/m}^3$
• A hard bedrock stratum lies at a depth of $15.0\text{ m}$ below the ground crest ($D = 15.0\text{ m}$).

Determine:

1. The depth factor $D_d$.
2. The Factor of Safety ($FS$) against short-term undrained failure using Taylor's stability chart value of $N_s = 0.160$.
3. The maximum allowable slope height $H_{max}$ if a minimum $FS = 1.50$ is required by design specifications.

Step-by-Step Solution

Step 1: Calculate the Depth Factor ($D_d$)

$$D_d = \frac{D}{H} = \frac{15.0\text{ m}}{10.0\text{ m}} = 1.50$$

Step 2: Calculate the Short-Term Factor of Safety ($FS$)

Using Taylor's equation $FS = \frac{c_u}{\gamma_{sat} H N_s}$:

$$FS = \frac{45.0\text{ kPa}}{(19.0\text{ kN/m}^3) \times (10.0\text{ m}) \times 0.160}$$ $$FS = \frac{45.0}{30.40} = 1.48$$

Step 3: Calculate Maximum Allowable Slope Height ($H_{max}$) for $FS = 1.50$

Rearranging Taylor's formula for $H$ with $FS_{req} = 1.50$:

$$H_{max} = \frac{c_u}{\gamma_{sat} \cdot FS_{req} \cdot N_s}$$ $$H_{max} = \frac{45.0\text{ kPa}}{(19.0\text{ kN/m}^3) \times 1.50 \times 0.160} = \frac{45.0}{4.56} = 9.87\text{ m}$$

Design Recommendation

Since the computed $FS = 1.48$ is slightly below the required $FS_{req} = 1.50$, the slope height must be trimmed to $9.87\text{ m}$ or the slope angle flattened to $\beta = 26^\circ$ to lower $N_s$ and satisfy civil infrastructure safety criteria.