Plastic Soil Structure Interaction: Foundation Design Guide (2026)

1. Fundamentals of Inelastic Soil-Structure Interaction (SSI)

In conventional structural design, foundations are frequently assumed to rest on perfectly rigid, unyielding supports. However, during strong earthquakes, high wind overturning, or heavy industrial loadings, foundation soils undergo significant plastic deformation, yielding, and interface debonding. Accounting for inelastic foundation-soil interaction allows engineers to capture the true coupled compliance and energy dissipation of the combined superstructure-substructure system.

                    COUPLED PLASTIC SSI EQUILIBRIUM
          Superstructure
             |      |
             |  M   |  Overturning Moment
             +--+---+
                |
                v  V (Vertical Gravity)
          +------------+  <-- Concrete Footing
          |            |
          +------------+ ---> H (Lateral Shear)
          / / / / / / /
        /  Zone I    /        /  Active    /   \  Zone II (Radial Shear Fan)
      /   Wedge    /-----     +------------+       \ Zone III (Passive Rankine Zone)
                           +---------------------------------

When soil reaches its plastic capacity, localized yielding and foundation uplift act as a natural seismic base isolation mechanism. Evaluating bearing capacity plasticity alongside structural frame hinging prevents overestimating base shears and ensures realistic plastic redistribution between the building and the ground.

By implementing inelastic foundation-soil interaction principles, geotechnical and structural engineers can determine the true foundation collapse load and prevent catastrophic bearing failures under complex multiaxial loads. Modeling plastic soil structure interaction accurately produces resilient foundation systems.

2. Geotechnical Plasticity Formulations and Failure Criteria

2.1 The Mohr-Coulomb and Drucker-Prager Yield Surfaces

Unlike metals whose yielding is pressure-independent ($J_2$ von Mises criterion), geomaterials (soils and rocks) exhibit shear strength that increases with hydrostatic confinement pressure.

The classical Mohr-Coulomb yield criterion is defined in terms of principal effective stresses $\sigma_1′ \ge \sigma_2′ \ge \sigma_3’$:

$$f(\boldsymbol{\sigma}’) = (\sigma_1′ – \sigma_3′) – (\sigma_1′ + \sigma_3′) \sin \phi’ – 2 c’ \cos \phi’ = 0$$

where $\phi’$ is the effective friction angle and $c’$ is effective cohesion. In stress invariant space ($p’ = \frac{1}{3}\text{ tr}(\boldsymbol{\sigma}’)$, $q = \sqrt{3 J_2}$), the smooth Drucker-Prager approximation is expressed as:

$$f(p’, q) = q – M p’ – k \le 0, \quad M = \frac{6 \sin \phi’}{3 – \sin \phi’}, \quad k = \frac{6 c’ \cos \phi’}{3 – \sin \phi’}$$

2.2 Non-Associated Plastic Flow Rules and Soil Dilatancy

If an associated flow rule ($g = f$) is applied to frictional geomaterials, the predicted volumetric plastic expansion (dilatancy) significantly exceeds experimental observations. A non-associated flow rule ($g e f$) with dilatancy angle $\psi < \phi'$ is utilized:

$$g(\boldsymbol{\sigma}’) = (\sigma_1′ – \sigma_3′) – (\sigma_1′ + \sigma_3′) \sin \psi – 2 c’ \cos \psi$$

2.3 Limit Analysis Theorems in Geotechnical Media

• Lower Bound Geotechnical Theorem: If a statically admissible stress field satisfies equilibrium, matches boundary tractions, and nowhere violates Mohr-Coulomb yield ($f \le 0$), the calculated load is a safe lower bound.

• Upper Bound Geotechnical Theorem: For any kinematically admissible velocity field compatible with plastic flow, equating external work rate to internal soil dissipation yields an upper bound collapse load.

Applying these bounds within plastic soil structure interaction models provides rigorous safety margins.

3. Shallow Foundation Plasticity Under Combined V-H-M Actions

3.1 The 3D Interaction Yield Envelope ($V-H-M$)

In building frameworks, footings experience simultaneous vertical forces ($V$), horizontal shears ($H$), and moments ($M$). The ultimate capacity is defined by a convex 3D yield surface in $(V, H, M)$ force space:

$$f(V, H, M) = \left( \frac{H}{\mu V_{ult}} ight)^2 + \left( \frac{M}{\psi B V_{ult}} ight)^2 – \left[ \frac{V}{V_{ult}} \left( 1 – \frac{V}{V_{ult}} ight) ight]^{2\alpha} \le 0$$

3.2 Foundation Uplift, Gap Debonding, and Soil Crushing

Under intense overturning moments, footings undergo partial interface separation (gap formation). Contact width reduces to $B’ = B – 2e$ where $e = M/V$. When $e \ge B/6$, footing edge uplift caps seismic demands transmitted to the superstructure.

3.3 Macro-Element Formulations for Computational SSI

Modern non-linear software models the entire soil-foundation interface via macro-elements. A macro-element collapses boundary degrees of freedom into $\mathbf{u} = [u_z, u_x, \theta]^T$ governed by generalized force vectors $\mathbf{Q} = [V, H, M]^T$.

4. Deep Foundation Inelastic Mechanics: Pile Groups and Nonlinear Springs

4.1 Nonlinear $p-y$, $t-z$, and $q-z$ Plastic Spring Idealizations

For pile foundations, non-linear SSI is modeled via distributed Winkler springs:

• Lateral Soil Resistance ($p-y$ curves): Models lateral pile-soil interaction yielding at $p_{ult}$.

• Axial Skin Friction ($t-z$ curves): Captures shear stress transfer along the pile shaft.

• Tip End Bearing ($q-z$ curves): Models compressive plastic punch-through at the pile base.

4.2 Kinematic vs Inertial Inelastic Interactions

Kinematic SSI filters high-frequency seismic waves due to foundation stiffness, while inertial SSI transmits superstructure dynamic overturning shears to the soil, triggering localized plastic yielding. Incorporating plastic soil structure interaction links both mechanisms smoothly.

3.4 Hardening Soil Models and Strain-Dependent Soil Stiffness

In cohesive and cohesionless soils, soil stiffness degrades non-linearly with shear strain amplitude well before ultimate plastic collapse. Modern geotechnical plasticity models, such as the Hardening Soil (HS) and Hardening Soil with Small-Strain Stiffness (HSsmall) models, incorporate hyperbolic stress-strain curves with distinct loading stiffness ($E_{50}$), unloading/reloading stiffness ($E_{ur}$), and very small-strain shear modulus ($G_0$). Under cyclic lateral loads, these formulations capture hysteretic energy dissipation, cyclic soil degradation, and progressive permanent plastic settlements beneath foundation edges.

4.3 Seismic Kinematic Bending in End-Bearing and Friction Piles

During seismic shaking, seismic shear waves propagating upward through soil strata impose lateral curvature on embedded piles even in the absence of superstructure inertial forces (kinematic interaction). At interfaces between soil layers with contrasting shear wave velocities ($V_{s1} e V_{s2}$), high kinematic bending moments develop. When combined with inertial moments from superstructure lateral drift, pile sections can form localized plastic hinges beneath the ground surface. Structural design codes mandate spiral transverse reinforcement throughout the upper pile length and at geotechnical layer boundaries to provide rotational ductility $\mu_\phi \ge 6.0$.

4.4 Liquefaction-Induced Lateral Spreading and Foundation Kinematics

In loose, saturated sandy soils, cyclic shear stresses generate excess pore water pressure leading to soil liquefaction and lateral ground spreading. During liquefaction, effective overburden stress drops to zero, depriving shallow foundations of bearing capacity and exposing pile shafts to severe lateral kinematic pressures from sliding non-liquefied crust layers. Plastic limit analysis models liquefied soil as a purely cohesive residual medium with undrained residual shear strength $S_r$, allowing engineers to verify whether pile groups can undergo plastic flexural redistribution without catastrophic structural collapse.

5. Comprehensive Step-by-Step Worked Numerical Example

Let us evaluate the ultimate capacity and foundation collapse load factor for a strip footing under combined vertical, lateral, and overturning actions.

5.1 Strip Footing Geometry, Soil Parameters, and Ultimate Vertical Capacity

Consider a continuous strip footing:

• Footing Width $B = 3.0\text{ m}$, Embedment $D = 1.0\text{ m}$.

• Soil properties: $c’ = 25.0\text{ kPa}$, $\phi’ = 30.0^\circ$, $\gamma = 18.5\text{ kN/m}^3$.

• Applied loads: $V = 600.0\text{ kN/m}$, $H = 120.0\text{ kN/m}$, $M = 240.0\text{ kNm/m}$.

5.2 Prandtl Bearing Capacity Factors

For $\phi’ = 30^\circ$: $N_q = 18.40$, $N_c = 30.14$, $N_\gamma = 22.40$.
Surcharge $q = \gamma D = 18.5\text{ kPa}$.
Pure vertical bearing capacity:
$$q_{ult,0} = c’ N_c + q N_q + 0.5 \gamma B N_\gamma = 753.50 + 340.40 + 621.60 = 1715.50\text{ kPa}$$
$$V_{ult,0} = 1715.50 \times 3.0 = 5146.50\text{ kN/m}$$

5.3 3D Interaction Function Evaluation Under Inclined Eccentric Load

Eccentricity: $e = \frac{M}{V} = \frac{240.0}{600.0} = 0.40\text{ m} < \frac{B}{6} = 0.50\text{ m}$ (No uplift gap).
Effective footing width: $B' = B – 2e = 3.0 – 2(0.40) = 2.20\text{ m}$.
Load inclination: $\tan \theta_L = \frac{120.0}{600.0} = 0.20 \implies \theta_L = 11.31^\circ$.

Meyerhof inclination factors:
$$i_c = i_q = (1 – 11.31/90)^2 = 0.7644, \quad i_\gamma = (1 – 11.31/30)^2 = 0.3881$$

Reduced bearing pressure:
$$q_{ult, inc} = (25.0 \times 30.14 \times 0.7644) + (18.5 \times 18.40 \times 0.7644) + (0.5 \times 18.5 \times 2.20 \times 22.40 \times 0.3881) = 1013.18\text{ kPa}$$

5.4 Foundation Plastic Collapse Load Factor Determination

Ultimate vertical capacity under combined $(V-H-M)$:
$$V_{ult, combined} = 1013.18 \times 2.20 = 2228.99\text{ kN/m}$$

Factor of safety against plastic collapse:
$$FoS_{collapse} = \frac{2228.99}{600.0} = 3.715$$

Combined shear and moment reduce bearing capacity by 56.7%, proving the need for coupled plastic soil structure interaction evaluations.

6. Code Provisions, Eurocode 7 / ASCE 7 Guidelines, and Settlement Tolerances

1. Eurocode 7 (EN 1997-1): Enforces partial safety factors ($\gamma_{\phi’} = 1.25, \gamma_{c’} = 1.25$); limits angular distortion to $\beta \le 1/500$.
2. ASCE 7-22 Chapter 19: Allows base shear reductions via foundation damping ($\beta_{SSI} \le 20\%$) and kinematic filtering.

7. Geotechnical-Structural Integration Synthesis

Incorporating plastic soil structure interaction eliminates the artificial barrier between geotechnical soil mechanics and structural frame analysis. By treating the foundation-soil interface as a yielding, energy-dissipating plastic mechanism, engineers achieve safer, more cost-effective foundation designs that realistically capture ultimate collapse limits.

References & Standards Cited

1. Meyerhof, G. G. (1953). “Bearing Capacity of Foundations under Eccentric and Inclined Loads.” Proc. 3rd ICSMFE, Zurich, 1, 440-445.
2. European Committee for Standardization (CEN). (2004). Eurocode 7: Geotechnical Design (EN 1997-1). Brussels, Belgium.
3. ASCE. (2022). Minimum Design Loads for Buildings (ASCE/SEI 7-22). Reston, VA.
4. Salençon, J. (2002). Yield Design. Springer-Verlag, Vienna.
5. Nova, R., & Montrasio, L. (1991). “Settlements of Shallow Foundations on Sand.” Géotechnique, 41(2), 243-256.

Frequently Asked Questions (FAQ)

How does plastic soil-structure interaction improve seismic building performance?

Localized soil yielding and foundation rocking dissipate seismic energy, acting as a natural base isolation mechanism that reduces superstructure ductility demands.

What is the fundamental difference between associated and non-associated soil plasticity?

Associated plasticity ($g = f$) overpredicts volumetric dilatancy in sands. Non-associated plasticity uses a distinct potential function ($g
e f$) with $psi < phi'$ to model realistic soil volume changes.

What is a macro-element in foundation engineering?

A macro-element represents the non-linear foundation behavior (bearing yield, uplift, sliding) as a single element with coupled $V-H-M$ degrees of freedom.

Why does overturning moment dramatically reduce vertical bearing capacity?

Moments shift the resultant load toward the footing edge, reducing effective contact width ($B' = B - 2e$) and concentrating shear stresses.

When is foundation uplift beneficial rather than hazardous?

Controlled rocking and uplift during earthquakes limit maximum overturning moments transmitted to columns, provided settlements remain acceptable.