Loads Part I — Environmental Loads
Comprehensive treatment of wave, wind, and current loads in offshore structural design, with emphasis on modern spectral analysis and API/ISO code methods.
1. Overview of Environmental Loads
Environmental loads—principally wave, wind, and current—are the dominant drivers of offshore structural behavior. Design codes (API RP 2A-WSD, ISO 19902) prescribe methodologies to estimate these loads based on water depth, return period, exposure category, and site geography. Modern practice increasingly incorporates site-specific data from hindcast models and metocean measurements rather than generic code assumptions, yielding more cost-effective designs.
2. Wave Loads — Fundamental Theory
Airy Linear Wave Theory
Airy (small-amplitude) wave theory provides the foundation for estimating water particle motion and applied forces. Key parameters include:
- H (Wave Height): Vertical distance from trough to crest. Typically quantified as significant wave height (Hs), the average height of the highest one-third of waves in a sea state.
- T (Period): Time for one complete wave cycle. Peak spectral period (Tp) and mean period (Tz) are commonly used.
- λ (Wavelength): Spatial period of the wave; related to T via dispersion relation: λ = gT²/(2π)·tanh(2πd/λ), where d is water depth and g = 9.81 m/s².
- d (Water Depth): Controls shallow vs. deep-water behavior. Deep water: d > λ/2; Shallow water: d < λ/20.
Water particle velocity u and acceleration u̇ at a given depth z are sinusoidal functions of wave phase:
where k = 2π/λ (wavenumber), ω = 2π/T (angular frequency), z = elevation above seafloor.
Morison Equation and Drag/Inertia Forces
The Morison equation, developed in the 1950s and still widely used in API RP 2A and ISO 19902, separates wave loading on slender cylindrical members (D/T > 50) into two components: drag (velocity-squared) and inertia (acceleration):
F = 0.5·ρ·Cd·D·u|u| + ρ·Cm·(π·D²/4)·u̇
where:
ρ = seawater density (typically 1025 kg/m³)
Cd = drag coefficient (typically 1.0–1.2 for smooth cylinders)
D = member outside diameter (m)
u = particle velocity (m/s)
Cm = inertia coefficient (typically 1.9–2.1)
u̇ = particle acceleration (m/s²)
Physical Interpretation: The drag term dominates when wave periods are short (high-frequency, small-scale waves). The inertia term dominates in longer-period waves and larger-diameter members. Combined, they produce cyclic loading that drives fatigue in jacket members.
Force Coefficients and Marine Growth
Drag and inertia coefficients are empirically determined functions of:
- Keulegan–Carpenter Number (KC): KC = u_max·T / D, which compares inertial to viscous effects. Low KC (<1): inertia-dominated; High KC (>10): drag-dominated.
- Reynolds Number (Re): Re = u_max·D / ν, which characterizes turbulence; ν = kinematic viscosity (~1.3e-6 m²/s for seawater at 10°C).
Marine growth (algae, shells, barnacles) increases effective diameter by 50–100 mm in tropical regions, up to 250+ mm in productive North Sea waters. Modern codes require explicit accounting: effective diameter = D + 2·growth thickness. Growth also roughens the surface, increasing Cd by 20–40%. API RP 2A Section 6 and ISO 19902 Appendix J provide design guidance on growth assumptions.
3. Wave Spectra and Return Period Concept
Spectral Representation (JONSWAP)
Real ocean waves are irregular, comprising many frequency components. The spectrum S(f) represents energy as a function of frequency. The JONSWAP (Joint North Sea Wave Observation Project) spectrum, adopted by API RP 2A and ISO 19902, is parameterized by:
where:
Hs = significant wave height (m)
fp = peak frequency (Hz)
γ = peakedness parameter (typically 1–7, default 3.3)
σ = frequency spreading parameter (0.07 for f < fp, 0.09 for f > fp)
The spectrum is integrated numerically to obtain force time-histories, which are then analyzed for fatigue damage using rainflow cycle counting. Modern fatigue analysis couples spectral analysis with S-N (stress-life) curves to estimate cumulative damage over the design life.
Return Period and Extreme Events
Design wave height is not the maximum observed wave, but a "return period" event—the statistically largest wave expected in a 100-year period (or other design event). Extreme value analysis fits historical wave data to distributions (Weibull, Gumbel) and extrapolates to rare events. For example, a 100-year return period Hs in the Gulf of Mexico might be 18–22 m, whereas the North Sea could see 15–18 m.
- Operational Limit: Highest sea state for safe operations; typically 10–15 year return period.
- Ultimate Limit State (ULS): Design extreme wave; typically 100-year return period for normal structures, 10,000-year for high-consequence facilities.
- Serviceability Limit State (SLS): Design for fatigue; cumulative effect of all sea states over design life (20–30 years).
4. Wind Loads
Wind Speed Profile and Exposure
Wind speed increases with elevation due to surface roughness. API RP 2A prescribes a power-law wind profile:
where:
V(z) = wind speed at elevation z above mean water level
V_ref = reference speed at z_ref = 10 m (or specified reference height)
Exponent 1/7 represents "normal" exposure (open water)
Rougher surfaces (coastal structures, industrial parks) use exponent 0.25–0.35; smooth open ocean uses exponent 0.10–0.15. Typical design wind speeds (10-minute average, 100-year return) are 45–55 m/s in U.S. Gulf of Mexico, 40–50 m/s in North Sea, and 30–45 m/s in Southeast Asia.
Drag Coefficients for Structural Members
Wind drag on a cylindrical element is analogous to wave drag:
where:
ρ_air = air density (1.225 kg/m³ at sea level)
Cs = shape drag coefficient (0.4–1.2 for cylinders, depending on surface finish)
D = projected width (m)
V = wind speed at height z (m/s)
API RP 2A Section 6 provides shape factors for cylinders, angles, deck areas, and superstructure components. Modern codes also account for dynamic wind-structure interaction and vortex-induced vibration (VIV), especially for slender topside structures.
5. Current Loads
Current Profiles and Origins
Offshore currents arise from tidal motion, ocean circulation (Gulf Stream, etc.), density-driven flows, and storm surge. Tidal currents are predictable and cyclic; ocean currents are more persistent but vary seasonally. Design codes typically combine worst-case tidal current with a percentage of the design wave's associated current.
- Surface Current: Strongest at mean water level; driven by wind and external pressure gradients.
- Tidal Current: Varies from surface to near-seafloor; typically symmetric about slack water.
- Subsurface/Gradient Current: May reverse direction at depth; common in continental shelves and fjords.
Current-Induced Drag
Current load is calculated similarly to steady-state wind, using the Morison equation without the inertia term (current is quasi-steady):
where:
ρ = seawater density
Cd = drag coefficient (1.0–1.2)
D = member diameter
V_c = current speed (m/s)
API RP 2A Table 6-3 and ISO 19902 Table 10 provide current drag coefficients. For grouped members (e.g., multiple legs of a jacket), blockage factors reduce drag on shielded members: F_blockage = (1 - solid fraction)·F_nominal.
6. Combined and Directional Loads
Wave, Wind, and Current Directionality
In reality, waves, wind, and current act simultaneously and may approach from different directions. Modern analysis uses directional spreading functions (e.g., cos²(θ/2)) to split energy among multiple approach directions. Code provisions often simplify to a few worst-case combinations:
- Aligned Loading: Wave, wind, and current all from the same direction (most severe).
- Opposing Wind and Current: Wind from one direction, current from opposite (tests lateral stability).
- Orthogonal Loading: Wave from one axis, wind/current from perpendicular axis (tests torsional and bi-directional effects).
Load Combination Factors
When combining different environmental loads, reduction factors account for the low probability of simultaneous extremes:
| Load Combination Scenario | Wave Height Factor | Wind Speed Factor | Current Speed Factor |
|---|---|---|---|
| Wave governs (design case) | 1.0 (design Hs) | 0.5–0.7 | 0.7–0.9 |
| Wind governs | 0.6–0.8 | 1.0 (design V) | 0.5–0.7 |
| Current governs (shallow water) | 0.6–0.8 | 0.6–0.8 | 1.0 (design Vc) |
| Fatigue (operational sea states) | Variable spectrum | Variable profile | Seasonal variation |
7. Ice Loads (Arctic Structures)
Platforms in Arctic waters (Canada, Alaska, Siberia) experience ice load from seasonal fast ice and multiyear ice floes. Ice loads are highly uncertain and difficult to model; typical design approaches include:
- Quasi-static approach: Estimate ice crushing stress and floe contact area; apply as steady distributed load.
- Dynamic approach: Model ice breaking as a time-dependent process; account for peak loads from buckling and crushing phases.
- Probabilistic approach: Fit historical ice impact data to extreme value distributions.
Typical ice loads for Arctic structures range from 500 kN to 10,000+ kN per impact, depending on floe thickness, velocity, and hull geometry. Modern designs use sloping or conical leg geometry to reduce ice pile-up and peak contact forces. Regulatory agencies (e.g., Canada-Newfoundland & Labrador Board) mandate strict ice load analysis and full-scale field testing.
8. Worked Example: Morison Equation Application
Example: Calculate Drag and Inertia Forces on a Jacket Leg
Given:
- Jacket leg diameter D = 0.762 m (30 in) with marine growth = 50 mm → D_eff = 0.862 m
- Water depth d = 150 m; at mid-depth z = -75 m
- Design wave: H = 14 m, T = 12 s (Tp = 12.5 s)
- Wavelength λ ≈ 224 m (deep-water limit for this period)
- Seawater density ρ = 1025 kg/m³
- Cd = 1.1, Cm = 1.95 (from API RP 2A for typical KC range)
Step 1: Calculate particle velocity and acceleration at z = -75 m
Amplitude of motion from Airy theory at the wave crest (worst case):
u_max ≈ (πH/T) · cosh(2π(z+d)/λ) / sinh(2πd/λ) ≈ (π·14/12) · e^(2π·75/224) / sinh(2π·150/224)
Simplified (deep-water approximation, d ≫ λ): u_max ≈ (πH/T) · e^(2π(z+d)/λ) ≈ 1.32 m/s
u̇_max ≈ (2π²H/T²) · e^(2π(z+d)/λ) ≈ 0.69 m/s²
Step 2: Calculate drag force per unit length
F_d = 0.5 · 1025 · 1.1 · 0.862 · 1.32 · 1.32 ≈ 800 N/m at peak velocity
Step 3: Calculate inertia force per unit length
F_i = 1025 · 1.95 · (π·0.862²/4) · 0.69 ≈ 380 N/m at peak acceleration
Step 4: Resultant force envelope
Total inertia-dominated force: F_max ≈ sqrt(800² + 380²) ≈ 890 N/m
For a 30 m segment of this leg, total force ≈ 26.7 kN. This cyclic load is applied to every wave cycle, driving fatigue damage assessment.
Summary: Environmental Load Design Principles
- Spectral Analysis: Use JONSWAP or site-specific spectra for fatigue assessment; single design waves for ultimate limit states.
- Morison Equation: Still the industry standard for wave loading on slender members, despite its 1950s origins; non-linear drag term u|u| is critical.
- Return Period and Risk: 100-year events are probabilistic; actual extremes may exceed or fall short. Risk-based design allows for economic optimization.
- Marine Growth: Cannot be ignored in fatigue analysis; increases both drag and effective inertia. Conservative assumption: 100–150 mm growth in tropics, 250+ mm in productive temperate zones.
- Combined Loading: Simultaneous waves, wind, and current; reduction factors prevent over-conservative design. Directional spreading is increasingly important in modern risk assessments.
- Code Compliance: API RP 2A and ISO 19902 provide tabulated coefficients and load combinations; always cross-check with site-specific studies.
Advertisement
Educational use only. All tools and course content on SmartUtilz are for informational and study purposes. Results must be independently verified by a qualified engineer before use in any design or safety-critical application. Read full disclaimer →