6.2 Soil Mechanics

Stress Distribution in Soils

Soil mechanics is the study of how soil behaves under loads and changing conditions. It forms the core of geotechnical engineering, and you'll rely on it whenever you design foundations, retaining walls, or slopes. This section covers four major areas: stress distribution, settlement calculations, slope and retaining wall stability, and effective stress with pore water pressure.

Effective Stress Principle and Stress Calculation Methods

All stress distribution in soils starts with the effective stress principle, which relates three quantities: total stress, pore water pressure, and effective stress. (More on this relationship in the final section.) The practical question is: when you place a load on the ground surface, how does the added stress spread through the soil below?

Several methods answer that question:

• Boussinesq theory is the classic approach. It calculates the vertical stress increase beneath point loads, line loads, and distributed loads, assuming the soil is a homogeneous, elastic half-space. For a point load:

$\Delta\sigma_z = \frac{3Q}{2\pi z^2} \cdot \frac{1}{(1 + (r/z)^2)^{5/2}}$

where $Q$ is the point load, $z$ is the depth below the surface, and $r$ is the horizontal distance from the load.

• Newmark's influence chart is a graphical tool built on Boussinesq theory. You overlay a scaled drawing of the loaded area onto the chart and count influence blocks to find the vertical stress increase at a given depth. It's especially useful for irregularly shaped foundations.
• The 2:1 method is a quick approximation. It assumes the loaded area spreads outward at a slope of 2 vertical to 1 horizontal as you go deeper. It's less accurate than Boussinesq but handy for preliminary estimates.
• Westergaard's solution applies when the soil is layered rather than uniform. It accounts for differences in stiffness between layers, which Boussinesq does not.

Visualization and Analysis of Stress Distribution

Engineers use two related tools to visualize how stress spreads beneath a load:

• Isobars are contour lines connecting points of equal vertical stress. They give you a map of the stress distribution pattern through the soil mass.
• Stress bulbs (or pressure bulbs) are the zones enclosed by isobars where the stress increase is significant. Their size and shape depend on the load configuration and soil properties.

Some useful rules of thumb for stress bulb depth:

• A strip footing's stress bulb extends to roughly 1.5 times the footing width.
• A square footing's stress bulb extends to about 2 times the footing width.

When adjacent footings are close together, their stress bulbs overlap. That overlap increases the total stress in the soil between them, which can lead to greater settlement than either footing would cause alone. This is why foundation spacing matters in design.

Soil Settlement Calculation

Components and Calculations of Soil Settlement

Settlement has three components, and they happen in sequence:

1. Immediate (elastic) settlement occurs right as the load is applied. The soil deforms elastically, with no change in water content. You calculate it using:

$S_e = \frac{qB(1-\nu^2)}{E_s} I_s$

where $q$ is the applied pressure, $B$ is the foundation width, $\nu$ is Poisson's ratio, $E_s$ is the soil's elastic modulus, and $I_s$ is an influence factor that depends on the shape and rigidity of the foundation.

1. Primary consolidation settlement is the gradual squeezing out of pore water from saturated fine-grained soils (clays and silts) under sustained load. This is usually the largest component in clay soils. For a normally consolidated clay:

$S_c = \frac{C_c H}{1+e_0} \log_{10}\frac{\sigma'_0 + \Delta\sigma'}{\sigma'_0}$

where $C_c$ is the compression index, $H$ is the thickness of the compressible layer, $e_0$ is the initial void ratio, $\sigma'_0$ is the initial effective stress, and $\Delta\sigma'$ is the stress increase from the new load.

1. Secondary compression (creep) continues after primary consolidation is complete. It's estimated using the secondary compression index $C_\alpha$ and tends to be smaller than primary consolidation, though it can be significant in highly organic soils.

The coefficient of consolidation ($c_v$) controls how fast primary consolidation happens. It's determined from laboratory oedometer (consolidation) tests.

Time-Dependent Settlement Analysis and Mitigation Techniques

Predicting when settlement occurs matters just as much as predicting how much. Terzaghi's one-dimensional consolidation theory provides the framework.

The time factor $T$ relates to the degree of consolidation $U$ (the fraction of total primary consolidation that has occurred). For early stages of consolidation ($U < 60\%$):

$U = \sqrt{\frac{4T}{\pi}}$

The time factor itself is defined as $T = \frac{c_v t}{H_{dr}^2}$, where $t$ is elapsed time and $H_{dr}$ is the length of the longest drainage path. This means consolidation time is proportional to the square of the drainage path length, which is why thicker clay layers take dramatically longer to consolidate.

When consolidation would take too long on its own, engineers use mitigation techniques:

• Preloading places a temporary surcharge on the site to force consolidation before construction begins. Once the target settlement is reached, the surcharge is removed.
• Vertical drains (sand drains or prefabricated vertical drains, often called "wick drains") shorten the drainage path by giving pore water a shorter route to escape. Since consolidation time depends on drainage path length squared, even modest reductions in path length speed things up significantly.
• Preloading + vertical drains combined can reduce consolidation time from years to months.
• Stone columns are vibro-compacted gravel columns installed through soft soil. They both stiffen the ground and provide drainage, improving bearing capacity while reducing settlement.

Slope and Retaining Structure Stability

Slope Stability Analysis Methods

Slope stability analysis determines the factor of safety ($F_S$) against a slope failing by sliding or rotating along a slip surface. A factor of safety of 1.0 means the slope is on the verge of failure. In practice, engineers target $F_S$ values of 1.3 to 1.5 for long-term stability, with the exact requirement depending on the consequences of failure and the confidence in the soil data.

Infinite slope method applies to long, shallow slopes where the potential slip surface runs roughly parallel to the ground surface. The factor of safety is:

$F_S = \frac{c'}{\gamma H \sin\beta \cos\beta} + \frac{\tan\phi'}{\tan\beta}$

where $c'$ is effective cohesion, $\gamma$ is the soil unit weight, $H$ is the depth to the slip surface, $\beta$ is the slope angle, and $\phi'$ is the effective friction angle. Notice that if the soil is cohesionless ($c' = 0$), the factor of safety depends only on the ratio of friction angle to slope angle.

Method of slices handles more complex slopes with curved (usually circular) failure surfaces. The potential sliding mass is divided into vertical slices, and equilibrium is checked for each. Two common versions:

• Fellenius (Ordinary) method satisfies moment equilibrium but ignores interslice forces, making it simpler but less accurate.
• Bishop's simplified method accounts for interslice normal forces and satisfies moment equilibrium. It's more accurate and widely used in practice.

Both are limit equilibrium methods, meaning they balance the driving forces (gravity pulling the soil downslope) against the resisting forces (shear strength along the slip surface).

Retaining Wall Stability and Earth Pressure Theories

A retaining wall must resist three modes of failure:

1. Overturning about the toe
2. Sliding along the base
3. Bearing capacity failure of the soil beneath the base

To check these, you first need to know the lateral earth pressure acting on the wall. Two theories provide this:

• Rankine's theory assumes a smooth wall (no friction between wall and soil) and calculates the active earth pressure coefficient as:

$K_a = \tan^2\left(45^\circ - \frac{\phi'}{2}\right)$

This is the simpler approach and works well when wall friction is small.

• Coulomb's theory accounts for wall friction ($\delta$) and a sloping backfill ($\beta$), giving a more general but more complex expression:

$K_a = \frac{\cos^2(\phi' - \alpha)}{\cos^2\alpha \cos(\delta + \alpha)\left[1 + \sqrt{\frac{\sin(\phi' + \delta)\sin(\phi' - \beta)}{\cos(\delta + \alpha)\cos(\beta - \alpha)}}\right]^2}$

where $\alpha$ is the inclination of the wall face from vertical. For a vertical wall with no backfill slope and no wall friction, Coulomb's equation reduces to Rankine's.

Soil reinforcement techniques can improve stability for both slopes and retaining structures:

• Geogrid reinforcement consists of polymer grids placed in horizontal layers within the soil. The grids increase the soil's effective shear strength and allow construction of steeper slopes or mechanically stabilized earth (MSE) walls.
• Soil nailing involves drilling steel bars (nails) into the face of a slope or excavation and grouting them in place. The nails act as passive reinforcement, providing additional resisting forces in both natural slopes and cut slopes.

Effective Stress and Pore Water Pressure

Effective Stress Principle and Applications

Effective stress is the stress carried by the soil skeleton (the grain-to-grain contact forces). It's the single most important concept in soil mechanics because it controls soil strength, compressibility, and volume change. The relationship is straightforward:

$\sigma' = \sigma - u$

where $\sigma$ is total stress and $u$ is pore water pressure.

Soil shear strength is expressed in terms of effective stress through the Mohr-Coulomb failure criterion:

$\tau_f = c' + \sigma'_n \tan\phi'$

where $\tau_f$ is the shear strength at failure, $c'$ is effective cohesion, $\sigma'_n$ is the effective normal stress on the failure plane, and $\phi'$ is the effective friction angle.

Liquefaction is a dramatic example of why effective stress matters. In saturated, loose sandy soils, earthquake shaking can generate excess pore water pressure rapidly. If the excess pore pressure rises to equal the total stress, effective stress drops to zero. With no grain-to-grain contact force, the soil loses all shear strength and behaves like a liquid.

Pore Water Pressure and Seepage Analysis

Pore water pressure comes in two forms:

• Hydrostatic pore pressure exists under static (no-flow) conditions and increases linearly with depth below the water table:

$u = \gamma_w h$

where $\gamma_w$ is the unit weight of water (approximately 9.81 kN/m³) and $h$ is the depth below the water table.

• Excess pore pressure develops when loads are applied (or removed) faster than the water can drain. It's what drives consolidation and dissipates over time.

Seepage analysis determines how water flows through soil under a hydraulic gradient. The primary tool is the flow net, a graphical construction of two families of curves:

• Flow lines show the paths water follows.
• Equipotential lines connect points of equal hydraulic head.

Flow lines and equipotential lines always cross at right angles. From a completed flow net, you can calculate both the seepage quantity and the pore pressure distribution throughout the soil.

A critical condition arises when the upward seepage force equals the buoyant weight of the soil, reducing effective stress to zero. This is the quicksand condition, and it occurs at the critical hydraulic gradient:

$i_c = \frac{\gamma_{sat} - \gamma_w}{\gamma_w}$

where $\gamma_{sat}$ is the saturated unit weight of the soil. For most soils, $i_c$ is approximately 1.0.

Above the water table, capillary rise pulls water upward through small pore spaces. The height of capillary rise is inversely proportional to particle size: fine-grained soils (clays, silts) can have capillary rise of several meters, while coarse sands may only see a few centimeters. This capillary water creates negative pore pressures (suction) above the water table, which temporarily increases effective stress in that zone.