10.1 Soil Compaction and Soil Compressibility/Consolidation

10.1 Soil Compaction and Soil Compressibility/Consolidation

Geotechnical analysis for civil infrastructure requires understanding how soil deposits deform and densify under artificial compaction and sustained foundation loading. In Philippine civil engineering practice, candidates frequently encounter problems distinguishing between the rapid mechanical expulsion of pore air (compaction) and the gradual, time-dependent expulsion of pore water from saturated void spaces (consolidation). Both processes govern structural serviceability, settlement limits, and subgrade load-bearing capacity.

---

Principles of Soil Compaction

Compaction is the artificial densification of soil by mechanical manipulation (rolling, tamping, or vibrating) to expel air from the void space without significant alteration of moisture content during the process. Densification increases soil shear strength, reduces future settlement, decreases hydraulic conductivity, and mitigates frost susceptibility or swelling potential.

Standard Proctor vs. Modified Proctor Tests

Laboratory compaction tests establish the relationship between soil moisture content ($w$) and dry unit weight ($\gamma_d$). Two primary standards define this procedure:

The compactive energy per unit volume is derived algebraically: $E = \frac{(\text{Number of Blows/Layer}) \times (\text{Number of Layers}) \times (\text{Hammer Weight}) \times (\text{Drop Height})}{\text{Mold Volume}}$

Compaction Curve & Zero Air Voids Line

Plotting dry unit weight $\gamma_d$ against gravimetric moisture content $w$ reveals an inverted parabolic compaction curve. At low water contents, water acts as a lubricating film between soil particles, allowing them to pack closer together under mechanical blows. Dry density increases until reaching the Optimum Moisture Content (OMC), corresponding to the Maximum Dry Unit Weight ($\gamma_{d,\max}$).

Beyond OMC, additional water displaces solid mineral particles because water has a significantly lower unit weight than solid mineral grains ($\gamma_w = 9.81\text{ kN/m}^3$ vs $\gamma_s \approx 26.5\text{ kN/m}^3$). Consequently, dry density declines.

   Dry Unit Weight (γ_d)
       ^
       |               ZAV Curve (S = 100%)
       |              / 
       |      Peak   /   Modified Proctor (Higher Energy)
       |       *--- / --.  
       |      / \  /     \
       |     /   *------- \ --. Standard Proctor (Lower Energy)
       |    /   / \        \
       |   /   /   \        \
       +-----------------------------------> Moisture Content (w)
              OMC_mod  OMC_std

The theoretical upper limit of compaction corresponds to complete saturation ($S = 100\% = 1.0$), termed the Zero Air Voids (ZAV) line: $\gamma_{zav} = \frac{G_s \gamma_w}{1 + e} = \frac{G_s \gamma_w}{1 + w G_s}$

Where $G_s$ is soil specific gravity, $w$ is gravimetric water content, and $\gamma_w = 9.81\text{ kN/m}^3$. If a soil contains a specified percentage of air voids $A$ (where $A = V_a / V$): $\gamma_d = \frac{(1 - A) G_s \gamma_w}{1 + w G_s}$

Field Compaction Control & Relative Compaction

Field specifications dictate the required quality of earth fill through the Relative Compaction ($RC$): $RC = \frac{\gamma_{d,\text{field}}}{\gamma_{d,\max(\text{lab})}} \times 100\%$ Most structural earthworks in highway and structural engineering require $RC \ge 95\%$ of Standard or Modified Proctor maximum dry unit weight.

Two classical testing methods verify field density:

1. Sand Cone Method (ASTM D1556): A test hole is excavated in the compacted lift. All excavated soil is preserved, weighed, and dried to obtain wet mass $M_{\text{wet}}$ and moisture content $w$. The hole is filled with pre-calibrated, dry Ottawa sand of known bulk density $\rho_{\text{sand}}$. The volume of the hole is calculated: $V_{\text{hole}} = \frac{M_{\text{sand, total}} - M_{\text{sand, cone}}}{\rho_{\text{sand}}}$ $\gamma_{d,\text{field}} = \frac{M_{\text{wet}} / (1 + w)}{V_{\text{hole}}} \times g$
2. Nuclear Density Gauge (ASTM D6938): Uses gamma radiation attenuation (Cesium-137) to measure total wet density and neutron thermalization (Americium-241/Beryllium) to determine volumetric moisture content instantly.

---

Terzaghi's One-Dimensional Consolidation Theory

Consolidation is the gradual reduction in volume of a saturated cohesive soil deposit resulting from the dissipation of excess pore water pressure under sustained static load. Unlike cohesionless soils (sands/gravels), whose high hydraulic conductivity allows instantaneous drainage, fine-grained saturated clays have extremely low permeability ($k \sim 10^{-7}$ to $10^{-10}\text{ m/s}$), causing settlement to extend over months, years, or decades.

The Spring-Piston Analogy

Karl Terzaghi conceptualized one-dimensional consolidation using a water-filled cylinder containing a spring and a perforated piston with a drainage valve:

• Initial State ($t = 0$): Valve closed. Applied vertical stress increment $\Delta \sigma$ is transferred entirely to the incompressible water. Excess pore water pressure $\Delta u = \Delta \sigma$, while the effective stress increase in the spring is zero: $\Delta \sigma' = 0$.
• Transient State ($0 < t < \infty$): Valve opened. Water escapes through the orifice under hydraulic gradient. Pore pressure $\Delta u$ dissipates, and load progressively transfers to the compressive spring. Effective stress increases: $\Delta \sigma' = \Delta \sigma - \Delta u$.
• Final Equilibrium ($t \to \infty$): Water outflow ceases. Excess pore pressure completely dissipates ($\Delta u = 0$). The entire stress increment is carried by the soil skeleton: $\Delta \sigma' = \Delta \sigma$.

$\sigma = \sigma' + u \implies \Delta \sigma' = \Delta \sigma - \Delta u$

Preconsolidation Pressure & Casagrande Construction

In the laboratory, an undisturbed soil specimen is tested in an oedometer (consolidometer) under incremental vertical loading. Plotting void ratio $e$ versus the logarithm of effective vertical stress $\log \sigma'$ yields the consolidation curve.

The preconsolidation pressure ($\sigma'_c$ or $p_c$) is the maximum past effective overburden pressure the soil has experienced in its geological history. Arthur Casagrande's graphical method locates $\sigma'_c$:

1. Identify the point of minimum radius of curvature (maximum curvature) on the $e - \log \sigma'$ curve.
2. Draw a horizontal reference line from this point.
3. Draw a tangent line to the curve through this point.
4. Bisect the angle between the horizontal line and the tangent line.
5. Project the straight-line portion of the virgin compression curve backward. The intersection of this projection with the angle bisector defines $\sigma'_c$.

Consolidation Stress History Classification

Soil is classified based on the Overconsolidation Ratio ($OCR$): $OCR = \frac{\sigma'_c}{\sigma'_0}$ Where $\sigma'_0$ is the current in-situ effective vertical overburden stress ($=\sum \gamma' z$).

• Normally Consolidated ($NC$) Clay ($OCR = 1$): Current overburden is the highest effective stress the deposit has ever experienced ($\sigma'_0 = \sigma'_c$). The soil operates entirely on the steep virgin compression curve governed by the compression index ($C_c$).
• Overconsolidated ($OC$) Clay ($OCR > 1$): Past preconsolidation pressure exceeds current overburden ($\sigma'_c > \sigma'_0$), caused by past glacial loads, erosion of overlying strata, or groundwater fluctuations. The soil deforms along the flatter recompression/swell curve ($C_s$ or $C_r$) until $\sigma'$ exceeds $\sigma'_c$.
• Underconsolidated Clay ($OCR < 1$): Deposit is still actively consolidating under its own weight or recently placed fill; excess pore water pressure has not fully dissipated.

---

Primary Consolidation Settlement Calculations

The general formula for one-dimensional primary consolidation settlement of a clay stratum of initial thickness $H$ and initial void ratio $e_0$ is: $S_c = H \frac{\Delta e}{1 + e_0}$

Where $\Delta e$ is the change in void ratio evaluated across three possible stress cases:

Case 1: Normally Consolidated Clay ($OCR = 1$)

Because $\sigma'_0 = \sigma'_c$, any applied stress increment $\Delta \sigma'$ immediately pushes the soil along the virgin compression branch: $S_c = \frac{C_c H}{1 + e_0} \log_{10} \left( \frac{\sigma'_0 + \Delta \sigma'}{\sigma'_0} \right)$

Where Skempton's empirical correlation estimates $C_c$ for undisturbed clays of normal sensitivity: $C_c \approx 0.009 (LL - 10)$ (with Liquid Limit $LL$ entered as an integer percentage, e.g., $LL = 50 \implies C_c = 0.36$).

Case 2: Overconsolidated Clay — Condition A ($\sigma'_0 + \Delta \sigma' \le \sigma'_c$)

The final effective stress remains below the preconsolidation pressure. Deformation occurs entirely along the recompression curve governed by the swell/recompression index ($C_s$ or $C_r$), where typically $C_s \approx \frac{1}{5} C_c \text{ to } \frac{1}{10} C_c$: $S_c = \frac{C_s H}{1 + e_0} \log_{10} \left( \frac{\sigma'_0 + \Delta \sigma'}{\sigma'_0} \right)$

Case 3: Overconsolidated Clay — Condition B ($\sigma'_0 + \Delta \sigma' > \sigma'_c$)

The applied stress exceeds the preconsolidation pressure. Settlement splits into two distinct segments: recompression from $\sigma'_0$ up to $\sigma'_c$, followed by virgin compression from $\sigma'_c$ to final stress $\sigma'_f = \sigma'_0 + \Delta \sigma'$: $S_c = \frac{C_s H}{1 + e_0} \log_{10} \left( \frac{\sigma'_c}{\sigma'_0} \right) + \frac{C_c H}{1 + e_0} \log_{10} \left( \frac{\sigma'_0 + \Delta \sigma'}{\sigma'_c} \right)$

---

Time Rate of Consolidation

Terzaghi's fundamental governing one-dimensional differential equation relates excess pore water pressure $u$, depth $z$, and elapsed time $t$: $\frac{\partial u}{\partial t} = c_v \frac{\partial^2 u}{\partial z^2}$

Where $c_v$ is the coefficient of consolidation ($m^2/\text{s}$ or $\text{cm}^2/\text{s}$), defined by soil permeability $k$ and coefficient of volume compressibility $m_v$: $c_v = \frac{k}{m_v \gamma_w} = \frac{k (1 + e_0)}{a_v \gamma_w}$ Here, $a_v = -\Delta e / \Delta \sigma'$ is the coefficient of compressibility, and $m_v = a_v / (1 + e_0)$.

Dimensionless Time Factor ($T_v$)

The progress of consolidation is quantified by the dimensionless Time Factor ($T_v$): $T_v = \frac{c_v t}{H_{dr}^2}$

Where $H_{dr}$ is the maximum length of the drainage path:

• Double Drainage (Two-Way): Pervious sand or gravel strata exist at both top and bottom boundaries. Water at the mid-depth travels a maximum distance of half the stratum thickness: $H_{dr} = \frac{H}{2}$
• Single Drainage (One-Way): Pervious material on one face, impermeable rock or hardpan on the opposite face. Water must travel the entire stratum thickness: $H_{dr} = H$

Average Degree of Consolidation ($U\%$)

The average degree of consolidation $U$ across the stratum is empirically approximated from $T_v$:

• For $U \le 60\%$ ($T_v \le 0.283$): $T_v = \frac{\pi}{4} \left( \frac{U\%}{100} \right)^2 \implies U\% = 100 \sqrt{\frac{4 T_v}{\pi}}$
• For $U > 60\%$ ($T_v > 0.283$): $T_v = 1.781 - 0.933 \log_{10} (100 - U\%)$

Standard benchmark values tested on the CELE include:

• At $U = 50\%$: $T_v = \frac{\pi}{4}(0.50)^2 = 0.197$
• At $U = 90\%$: $T_v = 0.848$

Drainage Comparison Scaling Law

For the same soil deposit ($c_v = \text{const}$) reaching the same degree of consolidation ($T_v = \text{const}$), the time required is proportional to $H_{dr}^2$: $\frac{t_1}{t_2} = \left( \frac{H_{dr,1}}{H_{dr,2}} \right)^2$ Thus, a clay layer with single drainage requires four times longer to achieve the same degree of consolidation as an identical layer with two-way drainage ($t_{\text{single}} = 4 \times t_{\text{double}}$).

---

Step-by-Step Worked Problem Examples

Worked Example 1: Consolidation Settlement of Overconsolidated Clay

Problem: A $4.0\text{ m}$ thick saturated clay layer is situated between an overlying dense sand deposit and underlying permeable gravel bed. The groundwater table is located at the top of the clay layer. The clay properties are: initial void ratio $e_0 = 0.85$, compression index $C_c = 0.32$, recompression index $C_s = 0.05$. The existing effective overburden pressure at the mid-depth of the clay stratum is $\sigma'_0 = 120\text{ kPa}$, and its preconsolidation pressure is $\sigma'_c = 160\text{ kPa}$. A surface building load produces a uniform vertical stress increase of $\Delta \sigma' = 90\text{ kPa}$ at the clay mid-depth.

1. Determine the primary consolidation settlement $S_c$.
2. If the laboratory-derived coefficient of consolidation is $c_v = 2.8 \times 10^{-4}\text{ cm}^2/\text{s}$, calculate the time in days required for the clay stratum to reach $50\%$ consolidation.

Solution:

Step 1: Check Consolidation Stress Category $\sigma'_0 = 120\text{ kPa}, \quad \sigma'_c = 160\text{ kPa} \implies OCR = \frac{160}{120} = 1.33 > 1.0 \quad (\text{Overconsolidated})$ $\sigma'_f = \sigma'_0 + \Delta \sigma' = 120 + 90 = 210\text{ kPa}$ Since $\sigma'_0 < \sigma'_c < \sigma'_f$ ($120 < 160 < 210\text{ kPa}$), this problem falls under Condition B (Case 3), traversing both the recompression curve and the virgin compression curve.

Step 2: Calculate Consolidation Settlement ($S_c$) $S_c = \frac{C_s H}{1 + e_0} \log_{10} \left( \frac{\sigma'_c}{\sigma'_0} \right) + \frac{C_c H}{1 + e_0} \log_{10} \left( \frac{\sigma'_f}{\sigma'_c} \right)$ Evaluate each term separately: $S_{c1} = \frac{0.05 \times 4.0}{1 + 0.85} \log_{10} \left( \frac{160}{120} \right) = \frac{0.20}{1.85} \log_{10} (1.3333) = 0.10811 \times 0.12494 = 0.01351\text{ m} = 13.51\text{ mm}$ $S_{c2} = \frac{0.32 \times 4.0}{1 + 0.85} \log_{10} \left( \frac{210}{160} \right) = \frac{1.28}{1.85} \log_{10} (1.3125) = 0.69189 \times 0.11810 = 0.08171\text{ m} = 81.71\text{ mm}$ $S_c = S_{c1} + S_{c2} = 13.51 + 81.71 = 95.22\text{ mm} \approx 95.2\text{ mm}$

Step 3: Determine Time to 50% Consolidation ($t_{50}$) The clay is bounded by sand on top and gravel on bottom, providing double drainage: $H_{dr} = \frac{H}{2} = \frac{4.0\text{ m}}{2} = 2.0\text{ m}$ Convert $c_v$ to consistent SI units ($\text{m}^2/\text{s}$): $c_v = 2.8 \times 10^{-4}\text{ cm}^2/\text{s} \times \left( \frac{1\text{ m}}{100\text{ cm}} \right)^2 = 2.8 \times 10^{-8}\text{ m}^2/\text{s}$ For $U = 50\%$: $T_v = \frac{\pi}{4}(0.50)^2 = 0.19635 \approx 0.197$ Solve for time $t$ from $T_v = \frac{c_v t}{H_{dr}^2}$: $t = \frac{T_v H_{dr}^2}{c_v} = \frac{0.197 \times (2.0\text{ m})^2}{2.8 \times 10^{-8}\text{ m}^2/\text{s}} = \frac{0.788}{2.8 \times 10^{-8}} = 28,142,857\text{ seconds}$ Convert seconds to days: $t = \frac{28,142,857}{86,400\text{ s/day}} = 325.7\text{ days}$

---

CELE Board Exam Traps & Strategic Checklists