Deep Foundation Geotechnics: High-Strain Dynamic Pile Testing and CAPWAP Analysis
Wave equation propagation, Case Method formulation, and signal matching CAPWAP numerical models for cast-in-place drilled shaft capacity verification.
High-rise skyscrapers, suspension bridges, and offshore wind turbines transfer mega-Newton axial foundation loads through deep drilled shafts and driven steel pipe piles. Verifying ultimate pile capacity without executing destructive static load tests relies on High-Strain Dynamic Pile Testing governed by one-dimensional stress wave mechanics.
1. 1D Stress Wave Propagation in Elastic Rods
When a heavy drop hammer ($5 \text{ to } 40\text{ tons}$) strikes a foundation pile head, an axial compression stress wave propagates downward at wave velocity $c = \sqrt{E / \rho}$, where $E$ is Young's modulus and $\rho$ is concrete/steel density.
The wave propagation satisfies the classic second-order hyperbolic wave equation:
$$\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2}$$
Instrumented strain transducers and piezoresistive accelerometers bolted opposite each other near the pile head record force $F(t) = A \cdot E \cdot \varepsilon(t)$ and particle velocity $v(t) = \int a(t) dt$.
Dynamic Pile Head Instrumentation Rig
[ Impact Hammer: Mass M ]
|
V Impact Strike
+---------------+ Pile Top
| [G1] [G2] | <-- Paired Strain Gauges (F)
| [A1] [A2] | <-- Paired Accelerometers (v)
| |
| Pile Shaft |
| (E, A, c) | Wave Speed c = sqrt(E/rho)
| |
| | Soil Shaft Resistance (Rs)
| |
| |
+---------------+
[ Toe ] Toe End Bearing (Rt)
2. The Case Method Total Resistance Formulation
By decoupling downward-traveling waves ($W_d$) and upward-traveling reflection waves ($W_u$) via acoustic impedance $Z = \frac{E A}{c}$:
$$W_d = \frac{1}{2} [F(t) + Z \cdot v(t)], \quad W_u = \frac{1}{2} [F(t) - Z \cdot v(t)]$$
The total static + dynamic soil resistance $R_{\text{total}}$ mobilized during impact is derived at wave turnaround time $t_1 + 2L/c$:
$$R_{\text{total}} = \frac{1}{2} [F(t_1) + Z \cdot v(t_1)] + \frac{1}{2} [F(t_1 + 2L/c) - Z \cdot v(t_1 + 2L/c)]$$
Subtracting velocity-dependent dynamic damping yields static capacity $R_s$:
$$R_s = R_{\text{total}} - J_c \left[ Z \cdot v(t_1) + F(t_1) - R_{\text{total}} \right]$$
3. Numerical Signal Matching: The CAPWAP Algorithm
While the Case Method yields instantaneous head approximations, CAPWAP (Case Pile Wave Analysis Program) utilizes continuous iterative non-linear optimization:
- The pile is discretized into 50-100 lumped mass-spring-damper elements.
- Soil surrounding each element is modeled with elasto-plastic springs (quake $q$) and Smith-type linear viscous dashpots.
- The measured velocity $v(t)$ is applied as boundary excitation. The numerical solver computes an output force curve $F_{\text{calc}}(t)$ and adjusts soil resistance distribution until signal matching error between $F_{\text{calc}}(t)$ and $F_{\text{meas}}(t)$ converges below $5%$.
import numpy as np
def compute_case_capacity(force_arr, velocity_arr, time_arr, pile_length, wave_speed, impedance):
dt = time_arr[1] - time_arr[0]
two_l_over_c = 2.0 * pile_length / wave_speed
samples_delay = int(round(two_l_over_c / dt))
# Find initial impact peak
peak_idx = np.argmax(force_arr)
t1 = peak_idx
t2 = t1 + samples_delay
if t2 >= len(force_arr):
raise ValueError("Record duration insufficient for 2L/c turnaround")
r_total = 0.5 * (force_arr[t1] + impedance * velocity_arr[t1]) + 0.5 * (force_arr[t2] - impedance * velocity_arr[t2])
return r_total
Dynamic wave analysis empowers civil engineers to verify multi-thousand-ton foundation capacities with millisecond non-destructive acoustic telemetry, safeguarding mega-structures against differential settlement failure.