Element Formulation & Architecture

The Pipe Stress FEA engine is built on a robust 3D structural solver. The core architecture utilizes a three-phase solving sequence:

1. Matrix Assembly & Preprocessing: Discretizing curves, generating soil-spring meshes, and building the global stiffness and mass matrices.
2. State Resolution: Applying loads and solving for displacements (utilizing an iterative penalty method for non-linear boundary conditions like friction and gaps).
3. Code Compliance: Extracting localized forces to evaluate ASME B31 stress ratios and flange leakage.

1.1 The 12-DOF Beam Element

The fundamental building block of the pipeline model is a 2-node, 12-Degree-of-Freedom (DOF) 3D Euler-Bernoulli beam element. Each node possesses 6 DOFs: three translational ($x, y, z$) and three rotational ($\theta_x, \theta_y, \theta_z$).

The local element stiffness matrix $[k]$ is a $12 \times 12$ matrix formulated using the pipe's cross-sectional properties:

• $A$: Cross-sectional area
• $E$: Elastic Modulus
• $G$: Shear Modulus
• $I$: Moment of Inertia
• $J$: Polar Moment of Inertia
• $L$: Element Length

ASME B31J Flexibility Integration

A distinct feature of this engine is the direct integration of ASME B31J Flexibility Factors ($k$) into the displacement solution. Rather than treating bends and tees as standard rigid pipes and applying stress multipliers post-solution, the engine inherently softens the stiffness matrix for these components.

The effective inertia and torsional constants are modified prior to matrix assembly:

$I_{in, eff} = \frac{I}{k_{in}}$ $I_{out, eff} = \frac{I}{k_{out}}$ $J_{eff} = \frac{J}{k_{tor}}$

This ensures that the global displacement geometry accurately reflects the localized ovalization and flexibility of piping fittings.

1.2 Coordinate Transformation

Because the local stiffness matrix $[k]$ is aligned with the element's longitudinal axis, it must be rotated to align with the Global $(X, Y, Z)$ coordinate system before assembly.

The engine calculates a $3 \times 3$ direction cosine matrix based on the element's spatial vector. A $12 \times 12$ transformation matrix $[T]$ is then constructed to rotate both translations and rotations.

The global stiffness matrix for the element $[K_e]$ is obtained via:

$[K_e] = [T]^T [k] [T]$

The Vertical Axis

The engine dynamically adjusts the local transverse and vertical axes based on the user's Global Settings. If the user defines Z as the vertical elevation axis, the horizontal plane is evaluated as X-Y, and the transformation matrix is rotated accordingly.

1.3 The Consistent Mass Matrix

For dynamic load cases (such as Modal Analysis or Seismic Occasional loads), the engine eschews simplified lumped-mass models in favor of a Consistent Mass Matrix $[M]$.

This $12 \times 12$ matrix distributes the mass of the steel pipe, internal fluid, external coating, and inline valves continuously along the beam's shape functions. It accounts for both translational mass ($\mu$) and rotational inertia ($J_{\rho}$):

$J_{\rho} = \mu \left( \frac{r_o^2 + r_i^2}{2} \right)$

This continuous mass distribution guarantees highly accurate fundamental frequency extraction during eigenvalue modal analysis, avoiding the "stair-stepping" inaccuracies common in lumped-mass models.