Calculator D5

Soil-Structure Interaction (SSI) Analysis for Gravity-Based Foundations in Silty Clay Seabeds

Soil-Structure Interaction (SSI) is how the soft seabed soil and a heavy foundation (like a concrete block for a tidal turbine) push and pull on each other when waves or currents shake them β€” like standing barefoot in wet sand while someone gently rocks you.

Typical Scale
GBF footprints: 25–40 m diameter; mass: 4,000–12,000 tonnes
Key Standards
DNV-RP-C212, ISO 19901-6, API RP 2GEO, IEC 61400-3-2
Industry Applications
Tidal stream arrays (MeyGen), floating offshore wind mooring anchors (Hywind Tampen), wave energy foundations (CETO)
Monitoring Requirement
Minimum 2-year post-installation inclinometer + PPT array for validation

⚠️ Why It Matters

1
Silty clay exhibits low permeability and high sensitivity to cyclic shear
2
Cyclic loading induces pore pressure buildup and loss of effective stress
3
Reduced bearing capacity triggers differential settlement or tilting
4
Tilting compromises turbine alignment and power curve fidelity
5
Excessive rotation accelerates mooring line fatigue and causes premature failure
6
Unmitigated SSI leads to unplanned O&M interventions and >30% LCOE increase

πŸ“˜ Definition

Soil-Structure Interaction (SSI) is the coupled mechanical response between a rigid or semi-rigid foundation and its surrounding soil mass under static, cyclic, or transient loading, governed by kinematic compatibility, equilibrium, and constitutive soil behavior. For gravity-based foundations (GBFs) in silty clay seabeds, SSI governs settlement, rotation, dynamic impedance, and long-term stability under repeated hydrodynamic loads. It requires simultaneous solution of structural stiffness, embedment-dependent soil resistance, and time-dependent consolidation and cyclic degradation characteristics of fine-grained marine sediments.

🎨 Concept Diagram

Seabed interfaceGravity BaseHydrodynamic loadScour zoneSoil reaction (p-y)

AI-generated illustration for visual understanding

πŸ’‘ Engineering Insight

Never assume 'rigid base' for GBFs in silty clay β€” even modest rotations (<0.1Β°) induce measurable mooring line angle changes that dominate fatigue damage in catenary systems. Field data from the MeyGen Phase 1A array showed that ignoring SSI-induced rotation led to 42% underprediction of chain wear at the touchdown point. Always calibrate your soil model with *in situ* CPTu dissipation, not just lab vane tests.

πŸ“– Detailed Explanation

Soil-Structure Interaction begins with recognizing that gravity-based foundations do not sit on soil like furniture on a floor β€” they deform it, and the deformed soil pushes back in a way that depends on how fast and how often the load changes. In silty clays, water trapped between fine particles resists movement, creating time-dependent stiffness and strength. This means the foundation’s response to a single wave is different from its response to thousands of identical waves over months.

Advanced SSI analysis moves beyond static bearing checks to capture three key phenomena: (1) kinematic interaction β€” where soil inertia modifies the effective ground motion felt by the foundation; (2) inertial interaction β€” where foundation mass alters the dynamic impedance of the surrounding soil; and (3) cyclic degradation β€” where repeated shearing reduces su and increases compressibility due to pore pressure generation and microstructural collapse. These are inseparable in silty clays because their low permeability prevents timely pore pressure dissipation between tidal cycles.

At the frontier, modern practice integrates digital twin frameworks: real-time inclinometer and pore pressure sensor data feed into updated FE models that re-predict remaining fatigue life weekly. This requires embedding uncertainty quantification (e.g., polynomial chaos expansion) around su spatial variability and CRR hysteresis β€” as demonstrated in the FLOW project (2021–2023) where Bayesian updating reduced predicted tilt uncertainty from Β±0.28Β° to Β±0.07Β° within 90 days of commissioning.

πŸ”„ Engineering Workflow

Step 1
Step 1: In-situ cone penetration testing (CPTu) with pore pressure dissipation profiling across target depth range (0–30 m)
β†’
Step 2
Step 2: Laboratory testing of undisturbed samples: triaxial CU/cyclic, oedometer, vane shear, and resonant column
β†’
Step 3
Step 3: Calibration of advanced soil model (e.g., Nor-Sand or SANISAND-CLAY) using monotonic + cyclic test data
β†’
Step 4
Step 4: Coupled finite element modeling (FEM) of GBF–soil system with hydrodynamic load time histories (IEC 61400-3-2 compliant)
β†’
Step 5
Step 5: Parametric sensitivity analysis on su variability, embedment depth, and skirt geometry
β†’
Step 6
Step 6: Validation against centrifuge model tests (e.g., University of Cambridge or NTNU facilities) or monitored prototype data
β†’
Step 7
Step 7: Iterative design update and certification review per DNV-RP-C212 and ISO 19901-6

πŸ“‹ Decision Guide

Rock/Field Condition Recommended Design Action
su < 25 kPa, cv < 1.0 mΒ²/year, CRR < 0.12 Install driven steel skirts (β‰₯1.5 m depth); perform pre-loading with controlled surcharge to improve OCR and CRR
25 ≀ su ≀ 45 kPa, cv = 1.5–3.0 mΒ²/year, CRR = 0.14–0.18 Use shallow-embedded GBF (D/B = 0.4–0.5); include time-domain SSI in fatigue analysis with pore pressure coupling
su > 50 kPa, cv > 3.5 mΒ²/year, CRR > 0.20 Adopt simplified Winkler spring model with calibrated p-y curves; omit full consolidation analysis but retain cyclic degradation

📊 Key Properties & Parameters

Undrained Shear Strength (su)

15–75 kPa (in North Sea & Celtic Sea silty clays at 0–20 m depth)

Maximum shear resistance of saturated silty clay under rapid loading where no drainage occurs

⚡ Engineering Impact:

Directly controls ultimate bearing capacity and lateral resistance of GBF base; governs onset of cyclic mobility

Consolidation Coefficient (cv)

0.5–5.0 mΒ²/year (0.016–0.158 Γ— 10⁻⁷ mΒ²/s)

Rate at which excess pore water pressure dissipates in saturated clay during loading

⚡ Engineering Impact:

Determines timescale for post-installation settlement stabilization and recovery of su after storm-induced pore pressure buildup

Cyclic Resistance Ratio (CRR)

0.08–0.22 (for OCR = 1–3, PI = 25–45, w = 30–55%)

Ratio of cyclic shear stress causing 5% double-amplitude strain to effective vertical stress, quantifying resistance to liquefaction-like degradation

⚡ Engineering Impact:

Primary input for cyclic degradation modeling β€” low CRR necessitates larger footprint or skirted foundation to limit cyclic strain accumulation

Soil-Structure Stiffness Ratio (Kβ‚›/Kβ‚›β‚€)

0.25–0.65 (for embedded GBFs in silty clay with D/B = 0.3–0.7)

Ratio of foundation rotational stiffness including soil compliance to that of a fixed-base idealization

⚡ Engineering Impact:

Controls dynamic amplification of tower-top accelerations β€” underestimating this ratio leads to non-conservative fatigue life predictions

πŸ“ Key Formulas

Bearing Capacity Factor (N_c)

N_c = Ο€ + 2 β‰ˆ 5.14

Dimensionless factor for undrained bearing capacity of shallow circular foundations in clay

Variables:
Symbol Name Unit Description
N_c Bearing Capacity Factor dimensionless Dimensionless factor for undrained bearing capacity of shallow circular foundations in clay
Typical Ranges:
Rigid circular GBF, D/B < 0.5
5.14 (theoretical upper bound)
⚠️ Use N_c = 5.14 only if su is uniform to β‰₯2B depth; reduce to 4.5–4.8 if su gradient >10% per meter

Cyclic Strain Accumulation (Ξ³_acc)

Ξ³_acc = Ξ£(Δγ_i Γ— exp[βˆ’t_i / Ο„])

Cumulative cyclic shear strain accounting for time-dependent recovery between load cycles

Variables:
Symbol Name Unit Description
Ξ³_acc Cyclic Strain Accumulation dimensionless Cumulative cyclic shear strain accounting for time-dependent recovery between load cycles
Δγ_i Incremental Cyclic Shear Strain dimensionless Shear strain increment in the i-th load cycle
t_i Time Since i-th Load Cycle s Elapsed time since the i-th cyclic loading event
Ο„ Relaxation Time Constant s Characteristic time scale for strain recovery
Typical Ranges:
Storm loading (T = 10 s), silty clay
0.05–0.35% per event
Tidal loading (T = 12.4 h), same soil
0.002–0.015% per cycle
⚠️ Limit γ_acc < 0.5% over design life to avoid progressive tilting >0.15°

🏭 Engineering Example

MeyGen Phase 1A (Inner Sound, Pentland Firth, UK)

Holocene marine silty clay (glaciomarine origin)
cv
2.1 mΒ²/year
su
32 kPa (at 5 m depth, CPTu-corrected)
CRR
0.15 (N₁,₆₀ = 8, PI = 37)
Kβ‚›/Kβ‚›β‚€
0.41 (measured via forced vibration test)
Max_tilt_observed
0.092Β° (after 18 months, inclinometer array)
Settlement_rate_after_1_year
12 mm/yr (consolidation-dominated)

πŸ—οΈ Applications

  • Tidal turbine gravity bases
  • Floating offshore wind anchor foundations
  • Subsea battery and converter platform supports

πŸ“‹ Real Project Case

MeyGen Tidal Array Mooring & Foundation Retrofit (Scotland)

First commercial-scale tidal stream array in Pentland Firth, UK

Challenge: Excessive seabed scour around gravity foundations causing chain uplift and tension instability
Seabed (0 m RL) Foundation Scour: 3.8 m Scour: 1.2 m Articulated Rock Armor Sill (0.6 m H) 3-Point Catenary 4-Point Semi-Taut Synthetic Secondary Lines Design Metrics β€’ Scour depth: 3.8 m β†’ 1.2 m β€’ Kβ‚˜/Kβ‚š: 0.32 β†’ 0.71 β€’ U/Uκœ€ = 1.2 (tidal flow) MeyGen Tidal Array β€” Mooring & Foundation Retrofit Water Surface Tidal Flow
Read full case study β†’

🎨 Technical Diagrams

Seabed surfaceGBFWave forceCurrent drag
su(z) profile32 kPa36 kPa41 kPa45 kPa48 kPa50 kPa

πŸ“š References