Compaction, Consolidation and Poroelasticity

Compaction vs Consolidation & Biot Poroelasticity

A complete interactive guide from classical soil mechanics to modern coupled theory

Intro Compaction vs Consolidation Terzaghi Theory Biot Theory Mandel-Cryer Poroelasticity Equations Solved Examples Interactive Calculator

1. Introduction

Compaction and consolidation are two fundamentally different processes that change the volume of soil. Understanding the difference is essential for foundation design, embankment construction and long-term settlement prediction. This article progresses from the classical one-dimensional theory of Terzaghi to the fully coupled three-dimensional poroelasticity of Biot, including the famous Mandel-Cryer effect.

2. Compaction vs Consolidation

Aspect	Compaction	Consolidation
Mechanism	Rolling, tamping, vibration	Sustained static load
What is expelled	Air	Water
Time	Immediate	Time-dependent (days–years)
Soils	Mostly granular	Fine-grained (clays)
Purpose	Increase density & strength	Causes settlement

Key Takeaway: Compaction is an active construction process; consolidation is a natural response of saturated clay under load.

3. Terzaghi’s One-Dimensional Consolidation Theory

Karl Terzaghi (1923–1925) developed the first rigorous theory of consolidation. It remains the workhorse of practical settlement calculations.

Assumptions

• Soil is homogeneous and fully saturated
• Solid grains and water are incompressible
• Flow and deformation are strictly one-dimensional (vertical)
• Darcy’s law is valid
• Coefficients \(k\) and \(m_v\) are constant
• Small strains

Governing Equation

\[ \frac{\partial u}{\partial t} = c_v \frac{\partial^2 u}{\partial z^2} \] where \( c_v = \dfrac{k}{m_v \gamma_w} \) is the coefficient of consolidation.

Time Factor & Degree of Consolidation

\[ T_v = \frac{c_v t}{H_{dr}^2}, \qquad U = \frac{s_t}{s_\infty} \]

Approximate relations:

• For \(U \le 60\%\): \( T_v \approx \dfrac{\pi}{4}\left(\dfrac{U}{100}\right)^2 \)
• For \(U > 60\%\): \( T_v \approx -0.933\log(1-U/100)-0.085 \)

4. Biot’s Three-Dimensional Consolidation Theory

Maurice Biot (1941) generalised consolidation to three dimensions and introduced true solid–fluid coupling. The theory is the foundation of modern poroelasticity.

Solid Skeleton
Deformation

↔

Pore Fluid
Pressure & Flow

→

Coupled
Response

Feature	Terzaghi	Biot
Dimensionality	1-D only	Fully 3-D
Coupling	Uncoupled	Fully coupled
Compressibility	Incompressible grains & fluid	Both compressible
Effective stress	\(\sigma'=\sigma-u\)	\(\sigma'=\sigma-\alpha p\)
Key phenomenon	Classic isochrones	Mandel-Cryer effect

5. The Mandel-Cryer Effect

In multi-dimensional consolidation the pore pressure at interior points can temporarily rise above the applied load before dissipating. This non-monotonic behaviour is called the Mandel-Cryer effect.

Physical reason: Drainage starts at the boundaries → outer zones contract → they squeeze the still-undrained interior → interior pore pressure rises temporarily.

6. Biot’s Poroelasticity Equations (Detailed)

Constitutive Law

\[ \sigma_{ij} = 2G\varepsilon_{ij} + \lambda\varepsilon\delta_{ij} - \alpha p\delta_{ij} \] \[ \zeta = \alpha\varepsilon + \frac{p}{M} \]

Equilibrium

\[ G\nabla^2\mathbf{u} + (K+\tfrac13G)\nabla(\nabla\cdot\mathbf{u}) - \alpha\nabla p + \mathbf{f} = 0 \]

Fluid Continuity + Darcy

\[ \alpha\frac{\partial}{\partial t}(\nabla\cdot\mathbf{u}) + \frac{1}{M}\frac{\partial p}{\partial t} = \nabla\cdot\left(\frac{k}{\mu}\nabla p\right) \]

Important parameters
• \(\alpha\) = Biot–Willis coefficient (\(1-K/K_s\))
• \(M\) = Biot modulus
• \(K_u = K + \alpha^2 M\) = undrained bulk modulus

7. Solved Examples

Example 1 – Terzaghi Time Calculation

A 4 m thick clay layer (double drainage) has \(c_v = 2.5\times10^{-7}\) m²/s. How long will it take to reach 90 % consolidation?

\(H_{dr} = 2\) m
For \(U=90\%\), \(T_v \approx 0.848\)
\[ t = \frac{T_v H_{dr}^2}{c_v} = \frac{0.848\times(2)^2}{2.5\times10^{-7}} = 1.357\times10^7\ \text{s} \approx 157\ \text{days} \]

Example 2 – Undrained vs Drained Response

A poroelastic material has \(K=20\) MPa, \(\alpha=0.9\), \(M=100\) MPa. Calculate the undrained bulk modulus.

\[ K_u = K + \alpha^2 M = 20 + (0.9)^2\times100 = 20 + 81 = 101\ \text{MPa} \]

8. Interactive Degree-of-Consolidation Calculator

Time Factor \(T_v\): 0.2 Drainage path \(H_{dr}\) (m): 2 Coefficient of consolidation \(c_v\) (m²/s × 10⁻⁷): 2.5

Average Degree of Consolidation U ≈ 50.4 %
Real time t ≈ 37.0 days

9. Quick Reference Flowchart

Load Applied

→

Excess Pore Pressure

→

Drainage Begins

→

Effective Stress ↑

→

Settlement

In multi-dimensional problems the intermediate step can produce the Mandel-Cryer peak.

Educational interactive article • Compaction, Terzaghi, Biot & Mandel-Cryer • For geotechnical learning

Yogendra Gopal Borse

Civil Engineer | Assistant Engineer Grade-I, Maharashtra PWD

B.Tech (Civil) from VJTI Mumbai. Experienced in bridge design, road works, estimation, project monitoring and digital engineering tools. Creator of YogiPWD – practical technical resources for civil engineers.

Read full About →