Stress Distribution and Soil Response

Introduction

When a foundation is constructed, its load is transferred to the ground through a contact area. The resulting increase in stress spreads downward and outward through the soil. This stress increase governs settlement, bearing-capacity mobilization, lateral deformation, and the interaction between the foundation and the supporting soil.

The stress state in soil is three-dimensional. In practice, engineers estimate stress increases using idealized elastic solutions, influence factors, and approximate distribution methods. The choice of method depends on foundation geometry, soil stratification, loading conditions, and the accuracy required.

Learning Objectives

By the end of this chapter, a student should be able to:

1. Distinguish total, effective, initial, and induced stresses.
2. Explain how foundation loads spread through soil and form stress bulbs.
3. Calculate vertical stress increases using Boussinesq and Westergaard solutions.
4. Apply influence-factor and approximate methods to loaded areas.
5. Evaluate contact-pressure distributions for flexible and rigid foundations.
6. Relate stress increases to elastic deformation and settlement.
7. Assess the effects of layered soils and soil–foundation interaction.
8. Recognize the assumptions and limitations of common stress-distribution methods.

1. Soil Stress and Stress Increase

The total vertical stress at depth is the weight of the soil and water above that depth. For a saturated soil under hydrostatic conditions, the effective vertical stress is

where σ′_v is effective vertical stress, σ_v is total vertical stress, and u is pore-water pressure. All stresses must be expressed in consistent units, such as kPa.

A foundation produces an additional vertical stress, commonly written as Δσ_z. The final effective stress, after drainage and pore-pressure dissipation, is approximately

for a drained loading condition. During undrained loading, part of the applied stress may initially appear as excess pore-water pressure rather than effective stress. Consequently, the time scale of loading is important, especially in clay.

Stress increase is not identical to settlement. Settlement depends on stress increase, compressibility, drainage, soil structure, stress history, and the thickness of compressible layers.

2. Concept of Stress Distribution

A surface load does not remain confined to the foundation footprint. In an idealized homogeneous, isotropic, elastic half-space, the load spreads both vertically and laterally. The region within which the induced stress is significant is called the stress bulb.

Stress bulbs are usually drawn using contours of equal vertical stress increase, such as 10%, 20%, or 30% of the applied foundation pressure. They are useful for identifying:

• the depth of soil significantly affected by a foundation;
• the potential interaction of neighboring foundations;
• whether a weak layer lies within the influence zone;
• the likely extent of settlement-related ground movement.

There is no universal depth at which stress increase becomes exactly zero. The practical depth of influence depends on the selected stress criterion and the sensitivity of the problem.

3. Boussinesq Solution

3.1 Point-load solution

Boussinesq developed an elastic solution for a vertical point load acting on the surface of a homogeneous, isotropic, linearly elastic, semi-infinite soil mass. The vertical stress increase at a point below the surface is

where:

• Q = vertical point load, in kN;
• z = depth below the ground surface, in m;
• r = horizontal radial distance from the load axis, in m;
• Δσ_z = induced vertical stress, in kPa or kN/m².

The expression is dimensionally consistent because Q/z² has units of kN/m².

Directly below the load, r = 0, so

The solution assumes a point load, linear elasticity, homogeneity, isotropy, and a semi-infinite soil mass. It is most useful for theoretical analysis and for constructing solutions for finite loaded areas by integration.

3.2 Circular and rectangular loaded areas

A uniformly loaded area can be represented as a collection of infinitesimal point loads. The vertical stress beneath the area is obtained by integrating the Boussinesq point-load solution. In routine practice, engineers commonly use influence factors rather than performing the full integration.

For a uniformly loaded flexible circular area of radius a, the stress increase below the center at depth z is

where q is the uniform contact pressure. The equation applies to a circular area on the surface of an elastic half-space and is evaluated below its center.

As z approaches zero, the expression approaches q. At great depth, the stress increase diminishes toward zero.

4. Westergaard Solution

The Westergaard solution idealizes the soil as horizontally reinforced by closely spaced, frictionless, perfectly flexible layers. This assumption restricts lateral deformation and produces a different stress distribution from Boussinesq’s isotropic model.

For a vertical point load, the vertical stress increase is commonly expressed as

At the axis of loading, r = 0, and therefore

The Westergaard and Boussinesq equations are not interchangeable. Boussinesq is generally associated with an isotropic elastic half-space, whereas Westergaard is more suitable as an idealized representation of a soil mass with restricted lateral strain, such as a heavily stratified deposit. Real soils fall between these idealizations, and field judgment is required.

5. Influence Factors

For a uniformly loaded foundation, stress increase is often written in the general form

where I_z is a dimensionless vertical-stress influence factor. The factor depends on foundation shape, location, depth ratio, and whether the point is below the center, corner, edge, or another position.

Influence-factor charts and tables are derived from elastic solutions. Their use requires careful matching of geometry and location. A factor for a corner point must not be used for a center point without a valid superposition procedure.

Rectangular areas

For a rectangular uniformly loaded area, the vertical stress beneath a corner can be obtained from an influence factor based on the ratios m = B/z and n = L/z, where B and L are the rectangle dimensions. Stress beneath the center can be calculated by dividing the rectangle into four equal rectangles and using superposition, provided the tabulated factor corresponds to a corner solution.

The principle of superposition is valid for linear elastic analysis:

It is not strictly valid when soil behavior is strongly nonlinear, when yielding occurs, or when the geometry and boundary conditions change substantially during loading.

6. Approximate Stress-Distribution Methods

6.1 The 2:1 method

A simple approximate method assumes that the load spreads downward at a ratio of two horizontal units to one vertical unit. For a uniformly loaded rectangular foundation of width B and length L, the average stress increase at depth z is estimated as

The expanded area at depth z is (B+z)(L+z). The method assumes uniform stress over this expanded area and is therefore an average-stress method, not a point-by-point elastic solution.

For a square foundation, B = L, and

The 2:1 method is quick and useful for preliminary settlement calculations, but it does not reproduce the exact stress variation near edges and corners.

6.2 Equivalent-area and strip methods

For irregular foundations, the loaded area may be divided into rectangles or strips. The stress contribution from each part is calculated separately and then combined. This is particularly useful for combined footings, mats, embankments, and staged construction.

Approximate methods should be checked against an elastic influence-factor method when the foundation is large, the soil profile is complex, or a weak layer is close to the foundation.

7. Worked Example 1: Point Load

A vertical point load of 600 kN acts at the ground surface. Determine the Boussinesq vertical stress increase at a point 4 m below the load and 2 m horizontally from its axis.

Given Q = 600 kN, z = 4 m, and r = 2 m,

Therefore,

The first factor is

and

Thus,

The estimated vertical stress increase is approximately 10.2 kPa.

8. Worked Example 2: 2:1 Distribution

A flexible rectangular footing is 3 m wide and 5 m long and carries a uniform pressure of 150 kPa. Estimate the average vertical stress increase 4 m below the footing using the 2:1 method.

The estimated average stress increase is 35.7 kPa. This result represents the average stress over the enlarged area at that depth; it is not the exact stress at the center or edge.

9. Contact Pressure Beneath Foundations

Contact pressure is the stress transmitted from the foundation base to the soil. It is not always equal to the applied structural load divided by foundation area.

Flexible foundations

A flexible foundation bends sufficiently to conform approximately to soil deformation. Under uniform loading, contact pressure may be relatively uniform, while settlement varies according to the soil profile and the foundation’s position.

Rigid foundations

A rigid foundation maintains approximately the same settlement across its base. The contact pressure therefore adjusts to the stiffness and support conditions of the soil. For a centrally applied vertical load on a uniform soil, pressure may be approximately uniform. An eccentric load produces a pressure gradient.

For a rectangular base carrying vertical load P and moments M_x and M_y, a simplified linear pressure distribution is

where A is base area, x and y are coordinates, and I_x and I_y are the second moments of area of the contact base about the relevant axes. The signs depend on the selected coordinate convention.

For one-axis eccentricity e = M/P, the maximum and minimum pressures for a rectangular base of width B are

These expressions assume full contact and a linear pressure distribution. Full compression requires e ≤ B/6; otherwise the calculated minimum pressure is negative, indicating loss of contact and the need for a compressed-contact analysis.

10. Soil–Foundation Interaction

The foundation and soil form a coupled system. Soil stiffness affects contact pressure, while foundation stiffness affects the distribution of settlement and load. Important interaction variables include:

• foundation rigidity;
• soil modulus and its variation with depth;
• load eccentricity and inclination;
• neighboring foundations;
• construction sequence;
• groundwater and drainage conditions;
• nonlinearity and yielding.

A perfectly rigid foundation does not necessarily produce a uniform soil pressure. Similarly, a flexible foundation does not necessarily experience uniform settlement. The pressure and deformation patterns must be considered together.

Elastic settlement is often estimated using a form of

where s is elastic settlement, q is representative net foundation pressure, B is foundation width, ν is Poisson’s ratio, E_s is soil Young’s modulus, and I_s is a dimensionless influence factor accounting for foundation shape, rigidity, depth, and the point or average settlement being evaluated.

This equation is an elastic estimate. It is most appropriate for small strains and soils that can reasonably be represented by an elastic modulus. The modulus must be compatible with the stress range and drainage condition. A modulus measured in one loading direction or strain range should not be applied uncritically to another condition.

11. Layered Soils

Elastic half-space equations assume homogeneous soil, but natural deposits are layered. A stiff layer over soft clay, a compressible stratum beneath dense sand, or alternating sand and clay can substantially change stress transmission and settlement.

Stress increase is affected by stiffness contrast. A stiff layer may redistribute load and restrict deformation, while a soft layer may undergo significant settlement even when the stress increase at the foundation level appears moderate. The following procedure is useful for preliminary assessment:

1. Divide the profile into layers with approximately uniform properties.
2. Calculate the stress increase at the top and bottom of each layer.
3. Use the average stress increase for settlement calculations, commonly the average of the top and bottom values for a sufficiently thin layer.
4. Select compressibility parameters appropriate to each layer.
5. Sum the layer settlements, while considering drainage and stress history.

For a layer of thickness H with constrained modulus M, a simple one-dimensional elastic estimate is

where Δσ′_avg is the average effective stress increase in the layer. This equation is appropriate only when the modulus definition and boundary conditions are compatible with one-dimensional compression.

Layered-soil analysis may require elastic-layer solutions, finite-element analysis, or numerical methods when stiffness contrasts are large or when the foundation is adjacent to slopes, excavations, or structural loads.

12. Case Application: Warehouse on a Layered Deposit

Background

This illustrative case concerns a 12 m × 18 m warehouse mat founded at ground level. The net contact pressure is 90 kPa. Beneath the mat is 3 m of dense sand, followed by 6 m of normally consolidated clay and then dense gravel.

Problem

The design team must determine whether the clay layer is significantly affected by the warehouse load and whether settlement may be excessive.

Analysis

A preliminary 2:1 estimate at the top of clay, 3 m below the mat, gives

At the bottom of the clay, 9 m below the mat,

The clay therefore receives a substantial stress increase over its full thickness. The average preliminary increase is approximately

The project team should not conclude that the dense sand alone controls performance. Consolidation settlement of the clay must be evaluated using its initial effective stress, compression index or constrained modulus, drainage path, and construction time.

Approach and solution

The recommended approach is to:

• verify the mat pressure and construction sequence;
• obtain representative consolidation and stiffness parameters for the clay;
• calculate stress increases using an elastic influence method as a check on the 2:1 estimate;
• evaluate immediate and consolidation settlement separately;
• assess differential settlement caused by variable clay thickness;
• consider staged construction or ground improvement if predicted settlement is unacceptable.

Outcome and learning points

The approximate analysis identifies the clay as a controlling layer. It also demonstrates that the stress bulb extends well below the footing dimensions. A foundation may appear to rest on competent material while a deeper compressible stratum governs long-term performance.

13. Limitations and Professional Judgment

Stress-distribution formulas are models, not direct measurements. Their limitations include:

• soil nonlinearity and stress-dependent stiffness;
• anisotropy and fabric effects;
• layering and interfaces;
• finite soil thickness or rigid boundaries;
• groundwater and coupled pore-pressure response;
• construction-induced changes in stress history;
• cyclic, dynamic, or seismic loading;
• footing embedment and excavation relief;
• adjacent structures and overlapping stress bulbs.

The engineer should compare methods rather than rely mechanically on one formula. A large difference between Boussinesq, Westergaard, and 2:1 results is a signal to examine the soil model, geometry, and design sensitivity.

14. Technology and Contemporary Practice

Spreadsheet calculations, geotechnical databases, and finite-element software allow stress distributions to be evaluated for irregular geometry and layered profiles. Numerical models can incorporate nonlinear constitutive laws, staged construction, groundwater flow, and soil–structure interaction. However, software output is not automatically reliable. Model selection, parameter calibration, mesh quality, boundary placement, and verification remain engineering responsibilities.

Digital site investigation systems and three-dimensional geological models can improve interpretation of variable strata. Machine-learning tools may assist with pattern recognition and parameter estimation, but they should not replace laboratory testing, professional judgment, or traceable design assumptions. Safety-critical calculations require independent checking and compliance with applicable building codes and geotechnical standards.

Best-Practice Checklist

Before adopting a stress-distribution result, confirm that:

1. The load, area, depth, and coordinate system are clearly defined.
2. Units are consistent.
3. The selected solution matches the foundation geometry.
4. The soil assumptions are stated and defensible.
5. Stress increases are distinguished from effective-stress increases.
6. Important layers lie within the stress bulb and are assessed.
7. Settlement and bearing-capacity checks use compatible parameters.
8. Eccentricity, contact loss, and neighboring foundations are considered.
9. Approximate results are checked where consequences are significant.
10. Calculations are independently reviewed and documented.

Chapter Summary

Foundation loads produce stress increases that spread through the soil and may affect deep compressible layers. Boussinesq’s solution models an isotropic elastic half-space, while Westergaard’s solution represents restricted lateral deformation. Influence factors provide practical solutions for finite loaded areas, and the 2:1 method offers a rapid estimate of average stress increase. Contact pressure depends on foundation rigidity, load position, and soil stiffness. Stress increase must be interpreted together with effective stress, drainage, compressibility, and settlement. Layered soils, eccentric loading, adjacent foundations, and nonlinear response can make simple formulas unreliable. Sound practice combines appropriate analytical methods, realistic soil parameters, engineering judgment, and independent verification.

Key Terms

Stress increase: Additional stress caused by an applied load.

Effective stress: Total stress minus pore-water pressure.

Boussinesq solution: Elastic solution for stress from a point load in a homogeneous isotropic half-space.

Westergaard solution: Stress solution based on restricted lateral deformation in a layered idealization.

Stress bulb: Region beneath a loaded area where induced stress is significant.

Influence factor: Dimensionless factor relating applied pressure to induced stress.

Superposition: Addition of responses from separate loads in a linear system.

Contact pressure: Stress transmitted between a foundation base and the soil.

Flexible foundation: Foundation that bends and conforms substantially to soil deformation.

Rigid foundation: Foundation that maintains approximately uniform settlement across its base.

Eccentricity: Distance between the load line and the centroidal axis of a foundation.

Elastic settlement: Deformation estimated using an elastic stiffness model.

Constrained modulus: Stiffness parameter for one-dimensional compression.

Layered soil: Soil profile consisting of strata with differing properties.

Consolidation: Time-dependent compression associated with drainage of pore water.

Homogeneous soil: Soil assumed to have uniform properties throughout the analysis domain.

Poisson’s ratio: Ratio describing lateral strain relative to axial strain in an elastic material.

Soil–foundation interaction: Mutual influence between foundation deformation, contact pressure, and soil response.

Review Questions

Short Answer Questions

1. What is meant by a stress bulb?
2. State the principal assumptions of the Boussinesq solution.
3. Why is stress increase not equivalent to settlement?
4. What does a vertical-stress influence factor represent?
5. What condition indicates possible loss of contact beneath an eccentrically loaded footing?

Long Answer Questions

1. Explain and compare the Boussinesq and Westergaard models, including their assumptions and practical limitations.
2. Describe how contact pressure differs beneath flexible and rigid foundations.
3. Explain how the 2:1 method estimates stress increase and discuss its advantages and disadvantages.
4. Discuss the effects of layered soils on stress transmission and settlement.
5. Explain the importance of soil–foundation interaction in foundation design.

Analytical / Application Questions

1. A 4 m × 6 m footing carries a uniform pressure of 120 kPa. Using the 2:1 method, estimate the average stress increase at 5 m depth.
2. A 500 kN point load acts at the ground surface. Calculate the Boussinesq stress increase 3 m below the load axis.
3. A foundation is supported by dense sand overlying soft clay. Propose a method sequence for deciding whether the clay controls settlement, and identify the field and laboratory data required.

Multiple Choice Questions with Answers

Question 1

Which expression gives the Boussinesq vertical stress increase from a surface point load?

A. Δσ_z = 3Q/(2πz²)[1+(r/z)²]^(−5/2)

B. Δσ_z = Q/(πz²)[1+2(r/z)²]^(−3/2)

C. Δσ_z = qBL/[(B+z)(L+z)

D. Δσ_z = qI_z

Correct Answer: A. Δσ_z = 3Q/(2πz²)[1+(r/z)²]^(−5/2)

Question 2

At the axis of a Boussinesq point load, the radial distance r is:

A. z

B. 0

C. B/2

D. infinite

Correct Answer: B. 0

Question 3

A 2:1 method calculation for a rectangular footing uses which denominator?

A. (B−z)(L−z)

B. BL

C. (B+z)(L+z)

D. B+L+z

Correct Answer: C. (B+z)(L+z)

Question 4

The Westergaard model primarily represents soil with:

A. No vertical stress

B. Restricted lateral deformation

C. Infinite tensile strength

D. Uniform pore pressure increase only

Correct Answer: B. Restricted lateral deformation

Question 5

For full compression beneath a rectangular footing under one-axis eccentricity, the usual condition is:

A. e ≤ B/6

B. e ≥ B/2

C. e = B

D. e > B/6

Correct Answer: A. e ≤ B/6

Question 6

A 2 m × 3 m footing carries 100 kPa. At 1 m depth, the 2:1 average stress increase is:

A. 20 kPa

B. 30 kPa

C. 50 kPa

D. 67 kPa

Correct Answer: C. 50 kPa

Question 7

Effective vertical stress is calculated as:

A. σ′_v = σ_v + u

B. σ′_v = σ_v − u

C. σ′_v = u − σ_v

D. σ′_v = σ_v/u

Correct Answer: B. σ′_v = σ_v − u

Question 8

Which factor is most directly associated with long-term consolidation settlement in saturated clay?

A. Drainage and compressibility

B. Paint thickness

C. Foundation color

D. Air temperature alone

Correct Answer: A. Drainage and compressibility

Question 9

The equation Δσ_z = qI_z requires I_z to be:

A. A length

B. A force

C. Dimensionless

D. A unit weight

Correct Answer: C. Dimensionless

Question 10

If a calculated minimum contact pressure is negative under a full-contact formula, the likely interpretation is:

A. The soil has zero stiffness

B. The foundation may lose contact over part of its base

C. The load is necessarily tensile

D. The stress bulb disappears

Correct Answer: B. The foundation may lose contact over part of its base