Soil Mechanics Simulator Back
Geotechnical Engineering

Soil Mechanics Simulator

Visualize Mohr circles and the Coulomb failure envelope in real time. Adjust cohesion and internal friction angle to calculate Terzaghi bearing capacity and consolidation settlement for shallow foundations interactively.

Soil Parameters
Cohesion c
kPa
Friction Angle φ
°
Unit Weight γ
kN/m³
Stress Deviator (σ₁-σ₃)/2
kPa
Confining Stress σ₃
kPa
Foundation Design
Foundation Type
Footing Width B
m
Embedment Depth Df
m
Consolidation
Coeff. of Consolidation cv
m²/yr
Clay Layer Thickness H
m
Final Settlement Sf
mm
Results
Results
480
Ult. Bearing Cap. qu (kPa)
3.2
Factor of Safety Fs
50
Undrained Su (kPa)
6.8
T₉₀ Consol. Time (yr)
Mohr Circle and Failure Envelope
Mohr
Main

What is Soil Failure & Foundation Settlement?

🙋
What exactly is a "Mohr circle" in soil mechanics? I see it drawn in the simulator, but what does it represent?
🎓
Basically, it's a clever graphical tool that shows all the possible combinations of normal stress ($\sigma$) and shear stress ($\tau$) on a soil element at a point. The circle's center and radius change based on the soil's state. In practice, when you adjust the "Vertical Stress" and "Horizontal Stress" sliders above, you're directly changing the size and position of that circle. A common case is analyzing the stress under a building's foundation.
🙋
Wait, really? So the sliders for "Cohesion" and "Friction Angle" are the soil's strength. How does that relate to the circle?
🎓
Great question! Those two parameters define the "Coulomb Failure Envelope," the diagonal line on the plot. The soil fails when the Mohr circle just touches that line. Try it: increase the vertical load (stress) with the slider. The circle grows until it touches the line—that's the failure condition you're simulating! For instance, this is how we determine if a slope will collapse or a foundation will punch into the ground.
🙋
So the "Bearing Capacity" and "Settlement" results below are calculated from this? Why does changing the "Foundation Width" affect both so much?
🎓
Exactly! They're two sides of the same coin. Bearing capacity (how much load before sudden shear failure) uses the failure criterion you just visualized. Settlement (slow, long-term sinking) depends on soil compressibility and layer thickness. A wider foundation spreads the load over more soil, increasing capacity but also involving a thicker compressible layer, which affects settlement. Play with the width slider and watch how the two results change in opposite directions—it's a key design trade-off.

Physical Model & Key Equations

The Mohr-Coulomb failure criterion defines the shear strength of the soil. Failure occurs when the shear stress on a plane reaches a value that depends linearly on the normal stress on that same plane.

$$\tau_f = c + \sigma \tan(\phi)$$

Where $\tau_f$ is the shear stress at failure, $c$ is the cohesion (soil's "stickiness"), $\sigma$ is the normal stress, and $\phi$ is the internal friction angle. In the simulator, $c$ and $\phi$ are set by their respective sliders.

Terzaghi's bearing capacity equation calculates the ultimate pressure a shallow foundation can support before a shear failure mechanism develops in the underlying soil.

$$q_u = c N_c + \gamma D_f N_q + 0.5 \gamma B N_\gamma$$

Here, $q_u$ is the ultimate bearing capacity, $\gamma$ is the soil unit weight, $D_f$ is the foundation depth, $B$ is the foundation width, and $N_c$, $N_q$, $N_\gamma$ are bearing capacity factors that depend only on the friction angle $\phi$. The simulator solves this when you adjust the geometry and soil strength parameters.

Real-World Applications

Shallow Foundation Design: This is the direct application. Every house, building, or bridge pier resting on footings or rafts requires a bearing capacity check to prevent catastrophic sinking and a settlement analysis to ensure doors won't jam and pipes won't crack from differential movement.

Slope Stability Analysis: The same Mohr-Coulomb failure criterion is used to assess the risk of landslides. Engineers analyze potential failure planes within a slope to calculate a factor of safety, helping to design safe road cuts, embankments, and open-pit mines.

Retaining Wall Design: The lateral earth pressure exerted on a retaining wall, which determines how thick and deep the wall must be, is derived directly from the soil's shear strength parameters (cohesion and friction angle) visualized in the simulator.

Construction Planning: Predicting consolidation settlement over time is critical for large infrastructure projects like embankments for highways or pre-loading of port terminals. It informs how much surcharge to apply and how long to wait before building the final structure.

Common Misconceptions and Points to Note

First, note that cohesion and internal friction angle are not determined independently. While the simulator lets you adjust them with separate sliders, in real soil—especially cohesive soil (clay)—the internal friction angle φ can be an apparent value. There's the concept of the "effective internal friction angle φ'" used in effective stress analysis and the "φ=0" approach used in total stress analysis. For example, a simplified method evaluating stability of soft clay using only cohesion c with "φ=0" is often used. After learning about the two extremes of sandy soil (high φ, low c) and clay layers (low φ, high c) with the tool, consider how both parameters interact in intermediate soils like "silt" or "sandy clay".

Next, Terzaghi's formula is not a "universal solution". It assumes a continuous footing (strip foundation). For square or circular isolated footings, you need to apply shape factors. For instance, the bearing capacity of a square footing with width B=2m is about 1.3 times that of a continuous footing with the same width. Also, the formula assumes ideal conditions of homogeneous, horizontal soil. In reality, soil is layered, and inclined loads or seismic forces may act. Use the "ultimate bearing capacity" obtained from the tool for understanding the basic concept and for parameter sensitivity analysis (to see which factors are influential).

Finally, understand the limitation that the "consolidation time" calculation is based on a one-dimensional model. The tool uses Terzaghi's one-dimensional consolidation theory; once the drainage distance H (e.g., half the clay layer thickness for double drainage) is determined, the settlement time history can be calculated. However, in the field, vertical drains like sand drains are installed to shorten the drainage path and accelerate consolidation. You can confirm in the tool that reducing the "drainage distance" dramatically shortens consolidation time—this is precisely one of the theoretical bases for ground improvement.

How to Use

  1. Enter cohesion (c) in kPa using the slider—typical values: 0 kPa for sand, 25–50 kPa for clay
  2. Set friction angle (φ) in degrees—range 25° for soft clay to 45° for dense sand
  3. Input unit weight (γ) in kN/m³—18 kN/m³ for loose sand, 20 kN/m³ for clay, 22 kN/m³ for dense material
  4. The simulator automatically generates the Mohr circle and calculates Terzaghi bearing capacity factor Nc, Nq, Nγ in real time
  5. Observe the failure envelope intersecting the stress circles to visualize the Coulomb criterion: τ = c + σ tan(φ)

Worked Example

For a shallow foundation on clay: c = 30 kPa, φ = 28°, γ = 19 kN/m³, foundation depth Df = 1.2 m, width B = 2 m. The simulator yields Nc ≈ 30.1, Nq ≈ 17.7, Nγ ≈ 15.1. Ultimate bearing capacity qu = cNc + γDfNq + 0.5γBNγ = 30(30.1) + 19(1.2)(17.7) + 0.5(19)(2)(15.1) = 903 + 402 + 287 = 1,592 kPa. Allowable capacity at factor of safety 3.0 = 531 kPa.

Practical Notes

  1. Undrained clay exhibits φu ≈ 0° with elevated c value; drained conditions show realistic φ angles—toggle between conditions to compare Mohr circles
  2. Sandy soils consolidate rapidly (hours to days); clayey soils require months or years—adjust γ and settlement rate predictions accordingly
  3. Bearing capacity factors increase nonlinearly with φ; a 5° increase in friction angle can raise Nγ by 30–40%, significantly affecting shallow foundation design
  4. The Mohr circle radius expands with applied principal stress difference; witness failure initiation when the circle touches the Coulomb envelope