Calculation Methods in Soil Bearing Capacity Analysis - Civil Engineering Knowledge & Design Tools

⚠️ Why It Matters

1

Inaccurate bearing capacity estimate

↓

2

Excessive settlement or differential movement

↓

3

Structural cracking in foundations or superstructure

↓

4

Premature failure of retaining systems or embankments

↓

5

Costly post-construction remediation or underpinning

↓

6

Regulatory non-compliance and project delay

📘 Definition

Soil bearing capacity is the maximum average contact pressure between a foundation and the supporting soil mass at which the soil fails in shear or undergoes excessive, non-recoverable deformation. It is determined by soil strength parameters (cohesion c, friction angle φ, unit weight γ), foundation geometry (width B, depth D), and loading conditions. Ultimate bearing capacity (qᵤ) must be reduced by an appropriate factor of safety (typically 2.5–3.0) to obtain the allowable bearing capacity (qₐ).

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Never default to Terzaghi for all shallow foundations — his method assumes rough, continuous, strip footings on homogeneous, weightless soil. Real-world applications demand Vesic for eccentric/inclined loads, Meyerhof for layered or sloping ground, and always cross-check with field correlations (e.g., SPT-based qₐ = 0.12×N₁₆₀×B⁰·⁵ for sands). When settlement governs (not strength), bearing capacity becomes secondary to modulus and compressibility calibration.

📖 Detailed Explanation

Bearing capacity analysis begins with recognizing that soil fails in shear, not compression — so the problem is fundamentally about mobilizing shear resistance along potential rupture surfaces. Terzaghi’s 1943 solution was the first rigorous limit equilibrium model for strip footings, introducing three dimensionless bearing capacity factors (N_c, N_q, N_γ) tied to φ', and assuming a rigid, smooth, weightless soil wedge beneath the footing.

Meyerhof extended this in 1951 to account for foundation shape (square, circular), depth (embedment ratio D/B), and load inclination — critical for real footings with finite width and non-vertical loads. His method introduced shape (s_c, s_q, s_γ) and depth (d_c, d_q, d_γ) correction factors, enabling more accurate estimates for isolated footings and shallow mats.

Vesic (1973) refined the framework further by deriving generalized forms for all factors based on plasticity theory and cavity expansion analogies, incorporating ground slope, foundation tilt, and base roughness. Modern practice combines these analytical methods with empirical correlations (e.g., DMT, CPT-based qₐ models), numerical validation (PLAXIS, FLAC), and probabilistic assessment where parameter uncertainty exceeds ±20% — especially in residual or weathered soils where c and φ' are spatially variable and scale-dependent.

🔄 Engineering Workflow

Step 1

Step 1: Site investigation — CPT/SPT, boreholes, lab testing (triaxial, CU/CD, consolidation)

→

Step 2

Step 2: Soil classification & parameter selection — identify c, φ', γ, k₀, Eₛ, and consolidation characteristics

→

Step 3

Step 3: Select bearing capacity theory — match soil type, foundation geometry, and loading (static/dynamic, drained/undrained)

→

Step 4

Step 4: Compute ultimate capacity (qᵤ) using selected method + shape/depth/inclination factors

→

Step 5

Step 5: Apply factor of safety (FS = 2.5–3.0 for permanent loads; ≥1.5 for temporary) → derive qₐ

→

Step 6

Step 6: Verify serviceability — compute immediate & consolidation settlement (Schmertmann, Skempton-Bjerrum) against tolerances

→

Step 7

Step 7: Integrate with structural design — iterate footing dimensions, reinforcement, and connection details

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
Saturated soft clay (cᵤ < 25 kPa, φ' ≈ 0°, OCR < 1.2)	Use undrained analysis (Terzaghi with N_c = 5.14); limit qₐ ≤ 1.5×cᵤ; consider preloading or wick drains.
Dense, well-graded sand (φ' ≥ 38°, γ = 20 kN/m³, B > 2 m)	Apply Vesic’s method with shape/depth corrections; verify against SPT-N₁₆₀ correlation (qₐ ≈ 0.15×N₁₆₀×B¹ᐟ² MPa).
Layered profile: 2 m stiff clay over loose sand (risk of punching shear)	Perform two-layer analysis (Hansen or Meyerhof); use equivalent φ' and weighted γ; verify interface shear transfer and embed foundation below weak layer.

📊 Key Properties & Parameters

Cohesion (c)

0–100 kPa (clays); 0 kPa (clean sands/gravel)

Shear strength intercept representing the inherent bonding between soil particles under zero normal stress.

⚡ Engineering Impact:

Dominates bearing capacity in fine-grained soils; governs short-term stability of excavations and footings on clays.

Effective Friction Angle (φ')

25°–45° (sands & gravels); <20° (soft clays)

Angle quantifying inter-particle resistance to sliding under effective stress conditions.

⚡ Engineering Impact:

Primary driver of bearing capacity in cohesionless soils; directly influences foundation embedment efficiency and slope stability.

Unit Weight (γ)

15–22 kN/m³ (dry to saturated soils)

Weight per unit volume of soil, including solids and pore fluids.

⚡ Engineering Impact:

Affects overburden pressure and self-weight contribution to bearing resistance; critical for depth factor calculations in Terzaghi/Vesic methods.

Foundation Width (B)

0.6–6.0 m (typical spread footings); up to 30 m (rafts)

Smaller plan dimension of a shallow foundation (e.g., strip or square footing).

⚡ Engineering Impact:

Directly scales shape and depth factors; wider footings increase qᵤ nonlinearly but may induce problematic settlement if compressible layers exist.

📐 Key Formulas

Terzaghi Ultimate Bearing Capacity (Strip Footing)

qᵤ = cN_c + qN_q + 0.5γBN_γ

Classic limit equilibrium solution for continuous footings on homogeneous soil.

Variables:

Symbol	Name	Unit	Description
q_u	Ultimate Bearing Capacity	kPa	Maximum pressure the soil can support without failure
c	Cohesion	kPa	Shear strength of soil at zero effective normal stress
N_c	Bearing Capacity Factor for Cohesion	dimensionless	Dimensionless factor dependent on soil friction angle
q	Effective Overburden Pressure	kPa	Vertical effective stress at the base of the footing
N_q	Bearing Capacity Factor for Surcharge	dimensionless	Dimensionless factor dependent on soil friction angle
γ	Unit Weight of Soil	kN/m3	Weight per unit volume of soil
B	Width of Footing	m	Width of the continuous (strip) footing
N_γ	Bearing Capacity Factor for Unit Weight	dimensionless	Dimensionless factor dependent on soil friction angle

Typical Ranges:

Clay (φ'=0°)

qᵤ = 5.14c (c = 10–70 kPa → qᵤ = 50–360 kPa)

Dense sand (φ'=38°)

N_q ≈ 49, N_γ ≈ 82 → qᵤ ≈ 1,200–2,800 kPa for B=1–3 m, q=30 kPa

⚠️ qₐ = qᵤ / FS; FS ≥ 3.0 for sustained dead loads

🔧 Open Calculator

Vesic Shape Factors (Square Footing)

s_q = 1 + (B/L)(N_q/N_c), s_γ = 1 − 0.4(B/L)

Corrects Terzaghi/Meyerhof for footing geometry; L = longer plan dimension.

Variables:

Symbol	Name	Unit	Description
s_q	Vesic shape factor for surcharge	dimensionless	Shape factor for square footing accounting for surcharge component
s_γ	Vesic shape factor for unit weight	dimensionless	Shape factor for square footing accounting for soil unit weight component
B	Shorter plan dimension of footing	m	Width of footing in the shorter plan direction
L	Longer plan dimension of footing	m	Length of footing in the longer plan direction
N_q	Bearing capacity factor for surcharge	dimensionless	Dimensionless bearing capacity factor dependent on soil friction angle
N_c	Bearing capacity factor for cohesion	dimensionless	Dimensionless bearing capacity factor dependent on soil friction angle