9.1 Ultimate & Allowable Bearing Capacity Theories

Key Takeaways

  • Terzaghi's ultimate bearing capacity equation models general shear failure beneath continuous strip footings under drained conditions using three components: cohesion, surcharge, and soil unit weight.
  • Shape factors (s_c, s_q, s_gamma), depth factors (d_c, d_q, d_gamma), and inclination factors (i_c, i_q, i_gamma) modify Terzaghi's base equation in Meyerhof, Hansen, and Vesic general bearing capacity frameworks.
  • Soil shear failure modes depend on density and compressibility: general shear failure occurs in dense sand and stiff clay, local shear failure in medium dense soil, and punching shear failure in loose sand and soft clay.
  • High groundwater tables reduce effective vertical stress, requiring buoyant unit weight adjustments depending on whether the water table is above the footing base, within depth B below the base, or below depth B.
  • Net allowable bearing capacity incorporates a factor of safety (typically FS = 3.0) applied to net ultimate bearing capacity (q_net,ult = q_ult - gamma * D_f), safeguarding against shear failure while accounting for excavation relief.
Last updated: July 2026

Introduction to Shallow Foundations

A shallow foundation is defined geotechnical engineering-wise as a structural element whose depth of embedment $D_f$ is less than or equal to its minimum width $B$ ($D_f / B \le 1$), though some modern codes extend this boundary to $D_f / B \le 4 \text{ to } 8$. Shallow foundations transfer structural column and wall loads to upper soil strata through direct base contact pressure. The fundamental design requirements demand that:

  1. The ultimate shear strength of the supporting soil mass is not exceeded (ultimate limit state for bearing capacity with an adequate Factor of Safety, $FS \ge 3.0$).
  2. Immediate elastic and long-term consolidation settlements remain within tolerable structural limits (serviceability limit state).

Shear Failure Modes Beneath Shallow Foundations

Depending on soil relative density $D_r$, soil compressibility, and foundation embedment depth, three distinct shear failure modes can develop beneath loaded footings:

Failure ModeSoil ConditionsPhysical CharacteristicsStress-Settlement Response
General Shear FailureDense sand ($D_r > 70%$), stiff to hard cohesive clayWell-defined continuous slip surfaces extend from footing edges to ground surface. Significant surface heaving on both sides.Sudden, catastrophic collapse with a clear peak ultimate load $q_{ult}$.
Local Shear FailureMedium dense sand ($30% \le D_r \le 70%$), medium stiff claySlip surfaces extend outward into radial shear zones but do not reach ground surface. Slight surface heave.Significant settlement before failure; well-defined peak missing ($q_{ult}$ designated at specific settlement).
Punching Shear FailureLoose sand ($D_r < 30%$), soft clay, weak upper crustSlip surfaces are restricted to vertical perimeter boundaries directly beneath footing. Zero surface heave; soil around footing remains flat.Continuous downward penetration with high settlements; curve shows no distinct yield point.
                      General Shear Failure Geometry
                      
                             Applied Load P
                                   │
                       ┌───────────┴───────────┐
                       │   Footing (B x L)     │
         Ground Surface └───┬───────────────┬───┘ Ground Surface
        ───┬───────────────┘               └───────────────┬───
          ╱ │ Zone I:                      Zone I: │ ╲
        ╱   │ Active Wedge               Active Wedge │   ╲  Zone III: Passive
      ╱     │ (Elastic)                     (Elastic) │     ╲ Rankine Zone
    ╱       └─────────────────┬───────────────────────┘       ╲
  ╱     Zone II: Radial       │      Zone II: Radial            ╲
 └────── Shear Zone ──────────┴────── Shear Zone ────────────────┘

For local and punching shear failures, Vesic and Terzaghi recommended adjusting soil shear strength parameters before applying standard bearing capacity equations: c=23cc'^* = \frac{2}{3} c' tanϕ=23tanϕ\tan\phi'^* = \frac{2}{3} \tan\phi'


Terzaghi's Ultimate Bearing Capacity Theory

Karl Terzaghi (1943) developed the first comprehensive theory for evaluating ultimate bearing capacity beneath a continuous strip footing of infinite length ($L \gg B$) under drained conditions. Terzaghi assumed a rigid active wedge (Zone I) directly beneath the footing base inclined at angle $\alpha = \phi'$ to the horizontal, flanked by radial shear zones (Zone II) and passive Rankine zones (Zone III).

Terzaghi's Equation for Continuous Strip Footings

qult=cNc+qNq+12γBNgammaq_{ult} = c' N_c + q N_q + \frac{1}{2} \gamma B N_gamma

Where:

  • $q_{ult} =$ Gross ultimate bearing capacity (kPa or psf)
  • $c' =$ Effective soil cohesion (kPa or psf)
  • $q = \gamma D_f =$ Effective vertical surcharge stress at footing base level (kPa or psf)
  • $\gamma =$ Total/moist unit weight of soil supporting the wedge (kN/m³ or pcf)
  • $B =$ Footing width (m or ft)
  • $N_c, N_q, N_gamma =$ Terzaghi non-dimensional bearing capacity factors (functions strictly of friction angle $\phi'$)

Mathematical Formulation of Bearing Capacity Factors

Nq=e2(3π/4ϕ/2)tanϕ2cos2(45+ϕ/2)N_q = \frac{e^{2(3\pi/4 - \phi'/2)\tan\phi'}}{2 \cos^2(45^\circ + \phi'/2)} Nc=(Nq1)cotϕN_c = (N_q - 1) \cot\phi' Ngamma=12(Kpγcos2ϕ1)tanϕN_gamma = \frac{1}{2} \left( \frac{K_{p\gamma}}{\cos^2\phi'} - 1 \right) \tan\phi'

Undrained Cohesive Soils (Saturated Clay, $\phi_u = 0$)

Under quick, undrained loading in saturated clay ($\phi_u = 0, c = c_u = s_u$): Nq=1.0,Nc=5.7 (Terzaghi) or 5.14 (Prandtl/Vesic),Ngamma=0N_q = 1.0, \quad N_c = 5.7 \text{ (Terzaghi)} \text{ or } 5.14 \text{ (Prandtl/Vesic)}, \quad N_gamma = 0 qult=cuNc+q=5.14su+γDfq_{ult} = c_u N_c + q = 5.14 s_u + \gamma D_f


General Bearing Capacity Equation (Meyerhof, Hansen, Vesic)

Terzaghi's original formulation was restricted to continuous strip footings loaded vertically. Meyerhof (1963), Hansen (1970), and Vesic (1973) expanded the equation to accommodate non-strip footprints, shallow-to-deep embedment, and inclined loadings:

qult=cNcscdcic+qNqsqdqiq+12γBNgammasgammadgammaigammaq_{ult} = c' N_c s_c d_c i_c + q N_q s_q d_q i_q + \frac{1}{2} \gamma B N_gamma s_gamma d_gamma i_gamma

Where $s_i$ are shape factors, $d_i$ are depth factors, and $i_i$ are load inclination factors.

Terzaghi Empirical Shape Modifications

For square and circular footings, Terzaghi simplified shape effects as follows:

  • Square Footing ($B \times B$): qult=1.3cNc+qNq+0.4γBNgammaq_{ult} = 1.3 c' N_c + q N_q + 0.4 \gamma B N_gamma
  • Circular Footing (Diameter $B$): qult=1.3cNc+qNq+0.3γBNgammaq_{ult} = 1.3 c' N_c + q N_q + 0.3 \gamma B N_gamma

Meyerhof / Vesic General Shape Factors

For rectangular footings ($B \times L$ where $B \le L$): sc=1+(BL)(NqNc)s_c = 1 + \left(\frac{B}{L}\right)\left(\frac{N_q}{N_c}\right) sq=1+(BL)tanϕs_q = 1 + \left(\frac{B}{L}\right) \tan\phi' sgamma=10.4(BL)s_gamma = 1 - 0.4 \left(\frac{B}{L}\right)

Vesic Bearing Capacity Factor $N_gamma$

In modern PE exam specifications, Vesic's $N_gamma$ equation is widely accepted: Ngamma=2(Nq+1)tanϕN_gamma = 2 (N_q + 1) \tan\phi'


Groundwater Table (GWT) Corrections

The presence of a groundwater table alters effective vertical stresses and unit weights within the shear failure zone. The correction depends on the depth $d_w$ of the water table relative to ground surface and footing base $D_f$.

               Case 1: GWT Above Footing Base (0 <= dw <= Df)
               
               Ground Surface ───────────────────────────── dw = 0 to Df
                              │                        │
                              ▼                        ▼
               GWT ═════════════════════════════════════════
                              │ (Df - dw)
               Footing Base   └───► ┌───────────┐
                                    │ B x L     │
                                    └───────────┘
GWT Depth PositionSurcharge Term $q$ AdjustmentUnit Weight Term $\gamma$ Adjustment
Case 1: $0 \le d_w \le D_f$<br/>Water table is above footing base$q = \gamma d_w + (\gamma_{sat} - \gamma_w)(D_f - d_w)$Use effective buoyant unit weight $\gamma' = \gamma_{sat} - \gamma_w$ in the third term.
Case 2: $D_f < d_w \le D_f + B$<br/>Water table is within failure wedge depth $B$ below baseUnmodified overburden: $q = \gamma D_f$Use weighted average unit weight $\gamma_{eff} = \gamma' + \frac{d}{B}(\gamma - \gamma')$, where $d = d_w - D_f$.
Case 3: $d_w > D_f + B$<br/>Water table is below failure wedge zoneUnmodified overburden: $q = \gamma D_f$Unmodified moist/total unit weight $\gamma$ in third term.

Factors of Safety and Allowable Bearing Capacity

To prevent shear failure with an adequate margin of safety, ultimate bearing capacity values are converted to allowable design pressures:

Gross Allowable Bearing Capacity ($q_{all}$)

qall=qultFSq_{all} = \frac{q_{ult}}{FS} Where $FS$ is typically $3.0$ for primary static building loads ($FS = 2.0 \text{ to } 2.5$ for temporary or seismic loading).

Net Ultimate Bearing Capacity ($q_{net,ult}$)

Net ultimate bearing capacity represents the net stress increase above existing overburden pressure: qnet,ult=qultq=qultγDfq_{net,ult} = q_{ult} - q = q_{ult} - \gamma D_f

Net Allowable Bearing Capacity ($q_{net,all}$)

qnet,all=qnet,ultFS=qultγDfFSq_{net,all} = \frac{q_{net,ult}}{FS} = \frac{q_{ult} - \gamma D_f}{FS}

Maximum allowable column force $P_{all}$ for a footing area $A = B \times L$: Pall=qnet,allAP_{all} = q_{net,all} \cdot A


Comprehensive Worked Calculation Example

Problem Statement

A square footing measuring $B = 2.0\text{ m} \times 2.0\text{ m}$ is founded at depth $D_f = 1.5\text{ m}$ in a sandy soil. The soil properties are:

  • Soil friction angle $\phi' = 32^\circ$
  • Effective cohesion $c' = 0\text{ kPa}$
  • Moist soil unit weight above GWT $\gamma = 18.0\text{ kN/m}^3$
  • Saturated soil unit weight $\gamma_{sat} = 19.5\text{ kN/m}^3$
  • Water unit weight $\gamma_w = 9.81\text{ kN/m}^3$
  • Groundwater table is located at $d_w = 0.5\text{ m}$ below ground surface.

Using Terzaghi bearing capacity factors for $\phi' = 32^\circ$ ($N_q = 23.2, N_gamma = 22.0$), compute:

  1. The gross ultimate bearing capacity $q_{ult}$
  2. The net ultimate bearing capacity $q_{net,ult}$
  3. The allowable net column load $P_{all}$ using $FS = 3.0$

Solution Steps

Step 1: Calculate Effective Surcharge $q$ (Case 1 GWT)

Since $d_w = 0.5\text{ m} < D_f = 1.5\text{ m}$, Case 1 applies: γ=γsatγw=19.59.81=9.69 kN/m3\gamma' = \gamma_{sat} - \gamma_w = 19.5 - 9.81 = 9.69\text{ kN/m}^3 q=γdw+γ(Dfdw)=(18.0×0.5)+(9.69×(1.50.5))=9.0+9.69=18.69 kPaq = \gamma d_w + \gamma' (D_f - d_w) = (18.0 \times 0.5) + (9.69 \times (1.5 - 0.5)) = 9.0 + 9.69 = 18.69\text{ kPa}

Step 2: Compute Terzaghi Gross Ultimate Capacity $q_{ult}$ for Square Footing

For square footing with $c' = 0$: qult=1.3cNc+qNq+0.4γBNgammaq_{ult} = 1.3 c' N_c + q N_q + 0.4 \gamma' B N_gamma qult=0+(18.69×23.2)+(0.4×9.69×2.0×22.0)q_{ult} = 0 + (18.69 \times 23.2) + (0.4 \times 9.69 \times 2.0 \times 22.0) qult=433.61+170.54=604.15 kPaq_{ult} = 433.61 + 170.54 = 604.15\text{ kPa}

Step 3: Compute Net Ultimate Bearing Capacity $q_{net,ult}$

qnet,ult=qultq=604.1518.69=585.46 kPaq_{net,ult} = q_{ult} - q = 604.15 - 18.69 = 585.46\text{ kPa}

Step 4: Compute Net Allowable Bearing Pressure and Column Load $P_{all}$

qnet,all=qnet,ultFS=585.463.0=195.15 kPaq_{net,all} = \frac{q_{net,ult}}{FS} = \frac{585.46}{3.0} = 195.15\text{ kPa} Pall=qnet,all(B×B)=195.15×(2.0×2.0)=780.6 kNP_{all} = q_{net,all} \cdot (B \times B) = 195.15 \times (2.0 \times 2.0) = 780.6\text{ kN}

Loading diagram...
Bearing Capacity Calculation Workflow and GWT Decision Tree
Loading diagram...
Terzaghi Tri-Zone Failure Mechanism Beneath Continuous Footing
Test Your Knowledge

A continuous strip footing of width B = 2.0 m is founded at a depth of D_f = 1.5 m in a saturated cohesive clay layer under undrained conditions (phi_u = 0, c_u = 80 kPa, gamma = 18.0 kN/m³). Using Vesic undrained bearing capacity theory with N_c = 5.14, what is the gross ultimate bearing capacity q_ult of the footing?

A
B
C
D
Test Your Knowledge

How does a high groundwater table located at the ground surface (d_w = 0) affect the ultimate bearing capacity of a shallow foundation compared to a completely dry soil condition?

A
B
C
D
Test Your Knowledge

A square footing (2.0 m x 2.0 m) and a continuous strip footing of width B = 2.0 m are constructed at D_f = 1.0 m in dry sand (phi' = 35°, c' = 0, gamma = 18.0 kN/m³). According to Terzaghi bearing capacity equations, why does the square footing have a different unit-weight term coefficient than the strip footing?

A
B
C
D
Test Your Knowledge

A column footing has a gross ultimate bearing capacity q_ult = 900 kPa at depth D_f = 1.5 m in soil with unit weight gamma = 18.0 kN/m³. If a factor of safety FS = 3.0 is specified, what is the net allowable bearing pressure q_net,all?

A
B
C
D