4.2 Engineering Properties of Soils

3.2 Engineering Properties of Soils

Permeability and Hydraulic Conductivity

The rate at which water flows through a soil mass is governed by its permeability (or hydraulic conductivity, k). This property is crucial for assessing seepage beneath dams, calculating flow rates into excavations, and designing drainage layers in pavements and retaining structures.

Darcy's Law

Water flow through soil is assumed to be laminar and is described by Darcy's Law, which states that the discharge velocity (v) of water through a soil cross-section is directly proportional to the hydraulic gradient (i): v = k × i The hydraulic gradient (i) is defined as the head loss (Δh) over the flow path distance (L): i = Δh / L The total volume flow rate (Q) through a cross-sectional area (A) perpendicular to the flow direction over a time (t) is given by: Q = k × i × A × t which can also be written as the flow rate per unit time: q = Q / t = k × i × A Where q is the volumetric discharge rate.

Laboratory Determination of Permeability

Two standard laboratory tests are utilized to measure k:

1. Constant Head Boundary Test (ASTM D2434): Ideal for coarse-grained soils (sands and gravels) with hydraulic conductivity values greater than 0.001 cm/s. The hydraulic head difference (h) across a soil specimen of length $L$ and area $A$ is held constant. The volume of water Q collected over time t is measured: k = (Q × L) / (A × h × t)
2. Falling Head Boundary Test: Utilized for fine-grained soils (silts and clays) where flow rates are very small. The specimen of length L and area A is connected to a standpipe of small cross-sectional area a. The water head in the standpipe drops from h1 to h2 over a time interval t: k = [a × L / (A × t)] × ln(h1 / h2) = 2.303 × [a × L / (A × t)] × log10(h1 / h2)

Seepage in Stratified Soil Deposits

Natural soils are deposited in horizontal layers, resulting in anisotropic permeability. Geotechnical engineers determine equivalent permeability coefficients (k_eq) for two flow directions:

• Flow Parallel to Stratification (Horizontal Flow): k_eq,h = (k1 × H1 + k2 × H2 + ... + kn × Hn) / H
• Flow Perpendicular to Stratification (Vertical Flow): k_eq,v = H / (H1/k1 + H2/k2 + ... + Hn/kn) Where Hi is the thickness of layer i, and H is the total thickness of the deposit (H = Σ Hi). The equivalent horizontal permeability k_eq,h is always larger than the vertical permeability k_eq,v due to the restrictive nature of lower-permeability layers in vertical flow.

---

Soil Compaction and Proctor Curves

Compaction is the process of using mechanical energy to pack soil particles closely together, reducing the volume of air voids, which increases shear strength and decreases future settlement.

Moisture-Density Relationship

Laboratory compaction tests (ASTM D698 for Standard Proctor and ASTM D1557 for Modified Proctor) establish the relationship between moisture content (w) and dry unit weight (γ_d).

• Standard Proctor: Compacts soil using a 24.4 N (5.5 lb) hammer dropped from a height of 304.8 mm (12 inches). The soil is compacted in 3 layers in a mold, receiving 25 blows per layer. Compaction energy = 592.7 kJ/m³ (12,400 ft-lbf/ft³).
• Modified Proctor: Uses a heavier 44.5 N (10.0 lb) hammer dropped from a height of 457.2 mm (18 inches). Soil is compacted in 5 layers with 25 blows per layer. Compaction energy = 2693 kJ/m³ (56,250 ft-lbf/ft³).

The Proctor curve is plotted as dry unit weight (γ_d) versus water content (w). As moisture is added, it acts as a lubricant, allowing particles to slide into a denser configuration. However, past a certain point, water begins to occupy space that would otherwise be occupied by solid particles, causing the dry density to decrease. This peak is the optimum moisture content (OMC) and the corresponding maximum dry unit weight (γ_d,max).

The Zero Air Voids (ZAV) Curve represents the theoretical dry unit weight of a soil containing zero air voids (degree of saturation S = 100%): γ_zav = (G_s × γ_w) / (1 + w × G_s) Because it is impossible to expel 100% of air voids through mechanical compaction, the compaction curve will always plot to the left of and below the ZAV curve.

Relative Compaction and Relative Density

For field quality control, construction specifications define the minimum relative compaction (RC): RC = (γ_d,field / γ_d,max,lab) × 100% For granular soils where laboratory compaction curves are difficult to define, relative density (Dr) is used to describe the state of packing: Dr = [(e_max − e) / (e_max − e_min)] × 100% Where:

• e_max = Void ratio of the soil in its loosest state.
• e_min = Void ratio of the soil in its densest state.
• e = In-situ void ratio of the soil.

---

One-Dimensional Soil Consolidation

Consolidation is a time-dependent process by which saturated, fine-grained soils (clays and silts) undergo settlement under sustained vertical loads. Because clay has low permeability, the pore water cannot escape instantaneously, transferring the load to the water phase as excess pore water pressure before slowly dissipating.

Terzaghi’s Spring Analogy and Time Factor

Karl Terzaghi modeled this using a water-filled cylinder containing a spring and a valve. When a load is applied:

1. Initial State (t = 0): Valve is closed. The load is carried entirely by the water pressure (excess pore pressure Δu = Δσ). Spring carries zero stress.
2. Intermediate State (0 < t < ∞): Valve is opened. Water begins to drain. Pore pressure decreases, and the spring (representing the soil skeleton) starts compressing, carrying a portion of the load.
3. Final State (t = ∞): Pore pressure has completely dissipated (Δu = 0). The spring carries the entire load (Δσ' = Δσ).

The progress of consolidation is defined by the dimensionless Time Factor (Tv): Tv = (cv × t) / d² Where:

• cv = Coefficient of consolidation (L²/T).
• t = Time.
• d = Drainage path length. For a clay layer draining from both the top and bottom boundary (double drainage), d = H_clay / 2. For single drainage (draining from only one boundary, e.g. bedrock at bottom), d = H_clay.

Primary Consolidation Settlement Calculations

The magnitude of consolidation settlement depends on the soil's stress history:

• Pre-consolidation Pressure (σ'c): The maximum historical effective vertical overburden stress the soil has ever experienced.
• Normally Consolidated (NC) Clay: The current effective overburden stress σ'v0 is equal to the pre-consolidation pressure σ'c (OCR = σ'c / σ'v0 = 1). Sc = [Cc × H / (1 + e0)] × log10[(σ'v0 + Δσ') / σ'v0]
• Overconsolidated (OC) Clay: The current effective overburden stress σ'v0 is less than the pre-consolidation pressure σ'c (OCR > 1).
  • Case 1: σ'v0 + Δσ' ≤ σ'c (Stress stays on the recompression curve): Sc = [Cr × H / (1 + e0)] × log10[(σ'v0 + Δσ') / σ'v0]
  • Case 2: σ'v0 + Δσ' > σ'c (Stress exceeds pre-consolidation pressure, moving past the bend): Sc = [Cr × H / (1 + e0)] × log10(σ'c / σ'v0) + [Cc × H / (1 + e0)] × log10[(σ'v0 + Δσ') / σ'c] Where:
• Cc = Compression index.
• Cr = Recompression (swell) index.
• e0 = Initial void ratio.
• H = Thickness of clay layer.
• σ'v0 = Initial effective vertical stress at the midpoint of the clay layer.
• Δσ' = Stress increase at the midpoint of the clay layer (calculated using Boussinesq or 2:1 distribution methods).

---

Shear Strength of Soils

The shear strength of a soil is its capacity to resist shearing stresses on a failure plane. It determines the bearing capacity of foundations, the stability of slopes, and lateral earth pressures on retaining walls.

Mohr-Coulomb Failure Criterion

The Mohr-Coulomb failure criterion states that shear strength (τ_f) along a plane is a linear function of the effective normal stress (σ') acting on that plane: τ_f = c' + σ' × tan(φ') Where:

• c' = Effective cohesion (attributable to intermolecular forces in clays).
• φ' = Effective angle of internal friction (attributable to particle interlocking and friction in sands/silts).
• σ' = Effective normal stress = σ − u (total normal stress minus pore water pressure).

Mohr's Circle and Stress Transformation

The stress state at failure can be plotted on a τ − σ diagram. The failure envelope is tangent to Mohr's circle at failure. The relationship between the major principal stress (σ'1) and minor principal stress (σ'3) at failure is: σ'1 = σ'3 × tan²(45° + φ'/2) + 2c' × tan(45° + φ'/2)

Laboratory Shear Strength Testing

1. Direct Shear Test (ASTM D3080): Soil is placed in a split metal box, consolidated under a normal force, and then sheared horizontally. The failure plane is forced along the split plane of the box, which is a major limitation because it may not represent the weakest plane of the soil.
2. Triaxial Shear Test (ASTM D2850/D4767): A cylindrical specimen is placed in a cell where it is subjected to an isotropic chamber confining pressure (σ3) and then loaded axially by a deviator stress (Δσd = σ1 − σ3) until failure. Three standard configurations exist:
   • Unconsolidated-Undrained (UU) Test: No drainage is permitted during consolidation or shear. Used to evaluate short-term stability (e.g., end-of-construction for an embankment on soft clay).
   • Consolidated-Undrained (CU) Test: Drainage is allowed under the confining pressure to consolidate the sample, but the drainage valves are closed during shear. Pore water pressures are measured to determine effective strength parameters (c' and φ').
   • Consolidated-Drained (CD) Test: Drainage is allowed throughout both stages. The loading rate must be slow enough to prevent pore pressure build-up. Used to model long-term stability.
3. Unconfined Compression Test (ASTM D2166): A rapid laboratory test performed on cohesive soils where the confining pressure σ3 = 0. The peak axial stress is the unconfined compressive strength (q_u). The undrained shear strength (s_u or c_u) is half of q_u: s_u = c_u = q_u / 2