Soil-Structure Interaction

Standard plant piping is supported by discrete, localized restraints (shoes, guides, anchors). Buried pipelines, however, are continuously supported and restrained by the surrounding soil.

To accurately capture this soil-structure interaction (SSI), the Pipe Stress FEA engine utilizes an automated meshing algorithm coupled with the empirical soil spring formulations defined by the American Society of Civil Engineers (ASCE) and the American Lifelines Alliance (ALA).

4.1 Automated Discretization

When an element is assigned a valid soil_id, the engine's preprocessor automatically subdivides the element into a dense mesh of shorter pipe segments.

The maximum mesh length ($L_{max}$) is calculated to ensure numerical stability and an accurate distribution of soil forces:

$L_{max} = \max(0.5\text{ m}, \min(2.0\text{ m}, 5 \times D_o))$

At each newly generated internal node, the engine attaches a 3D non-linear spring restraint (springRateX, springRateY, springRateZ) oriented to the local axis of the pipe.

4.2 ASCE/ALA Soil Spring Formulations

The engine evaluates the resistance of the soil in four distinct directions. For standard linear analysis, the ultimate yield force is divided by the characteristic yield displacement ($\Delta$) to determine the equivalent linear spring stiffness ($k$) per unit length, which is then multiplied by the node's tributary length ($L_{trib}$).

Parameters

Symbol	Definition
$\gamma$	Effective soil density ($N/m^3$)
$H$	Depth to the pipe centerline ($m$)
$\phi$	Soil friction angle (degrees)
$c$	Soil cohesion ($Pa$)
$\mu$	Pipe coating friction coefficient
$D_o$	Pipe outer diameter ($m$)

---

1. Axial Soil Spring ($T_u$)

Axial resistance is generated by the friction between the pipe coating and the surrounding soil.

Ultimate Force: $T_u = \pi D_o H \gamma \mu$

Yield Displacement: $\Delta_t = 0.005 \text{ m}$

---

2. Lateral Soil Spring ($P_u$)

Lateral resistance occurs when the pipe expands sideways into the trench wall (e.g., at horizontal bends). It is heavily dependent on the horizontal bearing capacity factor ($N_{qh}$).

Ultimate Force: $N_{qh} = \tan^2(45^\circ + \frac{\phi}{2})$ $P_u = D_o H \gamma N_{qh}$

Yield Displacement: $\Delta_p = 0.04 \left(H + \frac{D_o}{2}\right)$

---

3. Vertical Downward Bearing ($Q_d$)

Downward resistance acts as the primary foundational support for the pipeline's dead weight. It utilizes standard Terzaghi bearing capacity equations.

Ultimate Force: $N_q = e^{\pi \tan\phi} \tan^2\left(45^\circ + \frac{\phi}{2}\right)$ $N_\gamma = e^{0.18\phi - 2.5}$ $Q_d = (c D_o \times 5.14) + (\gamma H D_o N_q) + (0.5 \gamma D_o^2 N_\gamma)$

Yield Displacement: $\Delta_d = 0.1 D_o$

---

4. Vertical Upward Uplift ($Q_u$)

Upward resistance is engaged when thermal expansion or seismic forces attempt to push the pipe out of the trench. It is generally the weakest plane of resistance.

Ultimate Force: $N_{qv} = \frac{\phi H}{D_o}$ $Q_u = (c D_o) + (\gamma H D_o N_{qv})$

Yield Displacement: $\Delta_q = 0.015 H$

Equivalent Vertical Stiffness

Because a standard linear 6-DOF matrix cannot inherently differentiate between upward and downward stiffness on the same axis without iterating, the solver averages the upward ($Q_u$) and downward ($Q_d$) spring rates to define the local vertical stiffness ($k_y$) prior to matrix injection.