Vertical Stresses

Learning Objectives

  • Differentiate between geostatic and induced vertical stresses.
  • Calculate geostatic stresses based on soil unit weight and depth.
  • Calculate induced stresses from point loads using Boussinesq's and Westergaard's theories.
  • Compute stresses under area loads using Fadum's Chart and Newmark's Influence Chart.
  • Apply the Approximate 2:1 Method for preliminary stress estimation.

In geotechnical engineering, we must calculate the vertical stress (σv\sigma_v) at any depth to predict settlement and bearing capacity. Stresses arise from two primary sources: geostatic stresses from soil weight, and induced stresses from external loads like foundations.

Geostatic Stresses

In-situ vertical stresses arising from the self-weight of the overlying soil.

Induced Stresses

Additional vertical stresses caused by applied external loads, such as buildings, embankments, or vehicles.

Stresses due to Self-Weight

The geostatic vertical stress increases linearly with depth in a homogeneous soil profile.

Geostatic Stress Equation

The vertical stress is calculated by summing the products of the unit weight and thickness of each overlying soil layer.

Geostatic Vertical Stress

In-situ vertical stress at any depth due to the self-weight of all overlying soil layers; increases linearly with depth.

σv=γz\sigma_v = \sum \gamma z

Variables

SymbolDescriptionUnit
σv\sigma_vTotal vertical stress-
γ\gammaUnit weight of the soil layer-
zzThickness of the layer-

Effect of Groundwater

Below the water table, total stress uses γsat\gamma_{sat}, while effective stress uses γ\gamma'.

Stresses due to Point Loads (Boussinesq)

Boussinesq's Theory

An analytical solution that calculates stress distribution in a semi-infinite, homogeneous, isotropic, elastic medium due to a point load at the surface.

Boussinesq (1885) provided a foundational solution for point loads that is still widely used in modern geotechnical analysis.

Boussinesq's Equation

The vertical stress increase Δσz\Delta \sigma_z at depth zz and radial distance rr is:

Boussinesq Point Load

Vertical stress increase in an elastic half-space due to a concentrated surface point load; valid for isotropic, homogeneous soils.

Δσz=3P2πz2[11+(r/z)2]5/2\Delta \sigma_z = \frac{3P}{2\pi z^2} \left[ \frac{1}{1 + (r/z)^2} \right]^{5/2}

Variables

SymbolDescriptionUnit
Δσz\Delta \sigma_zVertical stress increase-
PPPoint load applied at the surface-
zzDepth below the surface-
rrRadial distance from the point load-
- Directly under the load (r=0r=0):

Stress Directly Under Load

Simplified Boussinesq formula for vertical stress directly on the axis of the applied point load (r = 0).

Δσz=3P2πz2=0.4775Pz2\Delta \sigma_z = \frac{3P}{2\pi z^2} = \frac{0.4775 P}{z^2}

Variables

SymbolDescriptionUnit
Δσz\Delta \sigma_zVertical stress increase-
PPPoint load-
zzDepth-

Stress Dissipation

Stress decreases rapidly with depth (proportional to 1/z21/z^2) and with radial distance from the point of application.

Westergaard's Theory

Westergaard's Theory

An alternative stress distribution theory assuming the soil mass is reinforced by inextensible horizontal layers, limiting horizontal deformation.

Westergaard's approach is more appropriate when the soil is highly stratified (e.g., layered sedimentary soils like varved clays) and restricts horizontal expansion.

Westergaard's Equation

An alternative stress distribution formula for horizontally layered (anisotropic) soils.

Westergaard Point Load

Alternative stress distribution formula for horizontally layered (anisotropic) soils; generally yields lower stresses than Boussinesq.

Δσz=Pπz21[1+2(r/z)2]3/2\Delta \sigma_z = \frac{P}{\pi z^2} \frac{1}{[1 + 2(r/z)^2]^{3/2}}

Variables

SymbolDescriptionUnit
Δσz\Delta \sigma_zVertical stress increase-
PPPoint load-
zzDepth-
rrRadial distance-

Comparison with Boussinesq

  • Westergaard's theory typically yields lower vertical stresses directly beneath the load (2/3\approx 2/3 of Boussinesq).
  • It is considered more appropriate for highly stratified or horizontally layered soils.

Stresses due to Area Loads

Foundations rarely apply true point loads. Instead, they distribute loads over an area (such as a strip, square, rectangle, or circle). The stress from area loads is found by integrating the point load solution over the foundation's footprint.

Fadum's Chart (Corner of a Rectangular Area)

To find the vertical stress increase (Δσz\Delta \sigma_z) exactly under the corner of a flexible rectangular area (B×LB \times L) loaded with a uniform pressure (qq), engineers use the mathematical integration provided by Fadum (1948).

Fadum's Stress Influence

Vertical stress increase under the corner of a uniformly loaded rectangular area using the Fadum influence chart factor.

Δσz=qIz\Delta \sigma_z = q \cdot I_z

Variables

SymbolDescriptionUnit
Δσz\Delta \sigma_zVertical stress increase-
qqUniform pressure-
IzI_zInfluence factor obtained from Fadum's Chart-

Using Fadum's Chart

  • The influence factor IzI_z is based on the dimensionless ratios m=B/zm = B/z and n=L/zn = L/z.
  • Superposition: To find the stress under the center (or any other point), divide the total area into smaller rectangles such that the point of interest is a corner to all of them, then sum their individual contributions.

Newmark's Influence Chart

A graphical method developed by Nathan Newmark (1942) to determine the vertical stress under any irregularly shaped loaded area. The chart consists of concentric circles and radial lines that divide the area into smaller sectors, each contributing an equal fraction of stress to the center point.

How to Use Newmark's Influence Chart

  1. Draw the foundation footprint to scale based on the depth zz where the stress is required (the scale is usually indicated on the chart, where a specific length represents zz).
  2. Superimpose the scaled drawing over the center of the influence chart such that the point where you want to find the stress is exactly at the origin (center) of the chart.
  3. Count the total number of sectors (NN) covered by the footprint of the foundation.

Newmark's Equation

Graphical method for determining vertical stress increase under any irregularly shaped foundation using Newmark's influence chart.

Δσz=INq\Delta \sigma_z = I \cdot N \cdot q

Variables

SymbolDescriptionUnit
Δσz\Delta \sigma_zVertical stress increase-
IIInfluence value per sector (typically 0.005)-
NNNumber of sectors covered by the foundation plan-
qqUniform contact pressure-

Approximate 2:1 Method

A simple and practical approximation used to estimate the vertical stress increase at a depth zz below a loaded rectangular foundation.

2:1 Method Calculation

The 2:1 Method assumes the applied load spreads linearly outwards at a slope of 2 vertical to 1 horizontal (an angle of 26.6\approx 26.6^{\circ}).

It is widely used because it is simple and generally conservative for preliminary design, meaning it tends to overestimate the stress increase compared to exact elastic solutions.

Approximate 2:1 Method

Practical approximation of vertical stress increase assuming loads spread at a 2:1 (vertical:horizontal) ratio; conservative for preliminary design.

Δσz=P(B+z)(L+z)\Delta \sigma_z = \frac{P}{(B + z)(L + z)}

Variables

SymbolDescriptionUnit
Δσz\Delta \sigma_zVertical stress increase-
PPTotal applied load-
BBWidth of foundation-
LLLength of foundation-
zzDepth below foundation base-

2:1 Method Characteristics

  • The area distributing the load effectively becomes (B+z)×(L+z)(B+z) \times (L+z).
  • This method is reasonably accurate for shallow depths but tends to diverge from precise elastic solutions at greater depths.
Key Takeaways
  • Vertical stress (σv\sigma_v) is the sum of geostatic stress (weight of soil) and induced stress (external loads).
  • Boussinesq's Theory assumes a homogeneous, isotropic, elastic half-space and is the standard for point-load stress distribution.
  • Westergaard's Theory is used for highly stratified soils and predicts lower stresses than Boussinesq.
  • Fadum's Chart provides the exact elastic solution for the stress increase beneath the corner of a uniformly loaded rectangular area.
  • The 2:1 Method is a simple, conservative approximation where the stress is assumed to spread over an expanding area (B+z)(L+z)(B+z)(L+z).