Skip to main content

Environmental Occasional Loads

Occasional loads represent short-duration, high-intensity dynamic events—most commonly Wind and Seismic activity.

Rather than requiring the engineer to calculate and manually apply wind pressures to every individual pipe, the engine automatically generates 3D spatial load profiles based on the selected structural design codes.

Dynamic Restraint Behavior

During Occasional load cases, the engine performs a "Locked Linear" state resolution. All dynamic snubbers instantly lock and act as rigid supports to resist the shock loads, while standard friction drag is ignored due to the dynamic nature of the event.

1. Wind Loading

Wind forces are modeled as uniform distributed loads (UDLs) projected onto the exposed area of the pipe and its insulation. The velocity pressure increases non-linearly with the elevation (zz) of the element.

The engine calculates the wind magnitude and then decomposes the resulting UDL into global XX, YY, and ZZ vectors based on the user-defined wind attack angle.

ASCE 7 (United States)

For ASCE 7, the engine calculates the velocity pressure (qzq_z) in N/m2N/m^2 using the standard formulation:

qz=0.613KzKztKdV2q_z = 0.613 K_z K_{zt} K_d V^2

The velocity pressure exposure coefficient (KzK_z) is derived dynamically based on the element's centroid elevation (zz) and the user-selected Exposure Category (B, C, or D):

Kz=2.01(zzg)2/αK_z = 2.01 \left( \frac{z}{z_g} \right)^{2/\alpha}

The final wind UDL (wwindw_{wind}) applied to the pipe is a function of the pressure, the coated outer diameter (DoD_o), and the cylindrical force coefficient (Cf=0.8C_f = 0.8):

wwind=qzDoCfw_{wind} = q_z D_o C_f

AS/NZS 1170.2 (Australia / New Zealand)

For AS 1170.2, the design wind speed (VdesV_{des}) incorporates regional multipliers:

Vdes=VRMdMz,catMsMtV_{des} = V_R M_d M_{z,cat} M_s M_t

The terrain/height multiplier (Mz,catM_{z,cat}) is dynamically interpolated from Table 4.1 based on the element's elevation and the selected Terrain Category (1, 2, 3, or 4). The resulting pressure and UDL are calculated as:

p=0.5ρairVdes2Cfigp = 0.5 \rho_{air} V_{des}^2 C_{fig}

wwind=pDow_{wind} = p D_o

(Note: ρair\rho_{air} is taken as 1.2 kg/m31.2 \text{ kg/m}^3 and the aerodynamic shape factor CfigC_{fig} for cylinders is taken as 1.21.2).

2. Seismic Loading

Seismic events are modeled using the Equivalent Static Force method. The engine extracts the total operating mass (mopem_{ope}) of each element—including steel, coating, fluid, and inline components—and subjects it to a uniform lateral acceleration.

AS 1170.4 Calculation

If the AS 1170.4 code is selected, the engine calculates the seismic coefficient (CdC_d) based on the site hazard factor (ZZ), probability factor (kpk_p), and the site sub-soil class multiplier (Ch(0)C_h(0)):

Cd=kpZCh(0)(acRc)C_d = k_p Z C_h(0) \left( \frac{a_c}{R_c} \right)

The resulting seismic UDL (wseisw_{seis}) applied along the user-defined earthquake axis is:

wseis=mopeCdgw_{seis} = m_{ope} C_d g

(If the "Generic" seismic code is selected, the user provides a raw G-factor which is directly multiplied by the element mass and gravity).

3. Valve Actuator Torsion (Advanced)

A critical feature of the FEA engine is its automated handling of center-of-mass offsets during seismic events.

Heavy valve actuators are rarely perfectly centered on the pipe axis. When subjected to a lateral seismic acceleration (FseisF_{seis}), this offset mass acts as a lever arm, inducing severe torsional twisting into the pipeline.

If a valve component contains an Actuator Mass (mactm_{act}) with vertical (dvd_v) or transverse (dtd_t) offsets, the engine automatically calculates the induced torsional moment (MtorM_{tor}) and injects it as a fixed-end rotational force into the local element matrix:

Mtor=(Fseis×dv)(Fseis×dt)M_{tor} = (F_{seis} \times d_{v}) - (F_{seis} \times d_{t})

This guarantees that torsional fatigue at valve stations is accurately captured without requiring the engineer to model rigid dummy-legs.