Ice-Structure Interaction
Ice-Structure Interaction: Theoretical Foundations
Overview of the Phenomenon
What kind of problem is ice-structure interaction?
It's the load evaluation when Arctic offshore oil platforms or icebreakers collide with sea ice. Ice exhibits complex failure behaviors such as crushing, buckling, and creep, so simple load models are sometimes insufficient.
Governing Equations
What kind of mechanical models are there for ice?
Ice is often modeled as a visco-elasto-plastic body. In Sinha's (1978) model,
the strain rate is decomposed into three components: elastic strain, delayed elastic strain (grain boundary slip), and viscous creep strain. For high strain rates (during impact), Maxwell-Brittle type models are also used.
The structure side uses standard elasto-plastic FEM. The ice-structure contact problem is handled as an interface using the penalty method or Augmented Lagrangian method. Ice crushing is represented by element erosion or CZM (Cohesive Zone Model).
Are there empirical formulas for ice load?
ISO 19906 (Arctic offshore structures) gives the ice pressure relative to the contact area $A$ as,
where $C_R$ is the Ice Reference Strength and $h$ is the ice thickness. Comparison with this empirical formula is useful for validating numerical simulations.
Ice is "More Complex Than Metal" โ The Mechanical Properties of Ice Change with Temperature, Speed, and Salinity
When dealing with ice in structural analysis, the first wall you face is the problem of "what material model should I use for ice?" Ice's strength more than doubles between -2ยฐC and -20ยฐC, it creeps (viscous flow) at low strain rates, and fails in a brittle manner at high strain rates. Furthermore, sea ice is weaker than pure ice because it contains salt, so even the same "ice" can have vastly different properties. In engineering, the boundary where "the failure mode transitions with strain rate" is said to be around 10โปยณ/s. Brittle fracture dominates in ship collisions (high speed), while creep dominates in ice pressure on bridge piers (low speed). Representing these two failure modes with a single material model is the theoretical core of ice-structure interaction simulation.
Computational Methods for Ice-Structure Interaction
Discretization Methods
How is ice crushing handled numerically?
There are three main approaches.
| Method | Characteristics | Application |
|---|---|---|
| FEM + Element Erosion | Element deletion upon reaching failure criteria | LS-DYNA, Abaqus/Explicit |
| DEM (Discrete Element Method) | Represents ice as an aggregate of particles | PFC, YADE |
| SPH | Mesh-free. Easy to track crushing. | LS-DYNA SPH |
| Peridynamics | Non-local model. Cracks occur naturally. | Peridigm |
Using DEM for ice sounds interesting.
In DEM, ice is represented as an aggregate of many disks (2D) or spheres (3D), and bonds between particles break when they exceed failure criteria. It has the advantage of naturally reproducing ice crushing patterns (radial crack, circumferential crack).
Contact Algorithm
How is contact between ice and structure handled?
In LS-DYNA, frictional contact is defined using *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE. The friction coefficient of ice strongly depends on temperature, varying in the range $\mu = 0.01$ to $0.3$.
The explicit method time step follows the Courant condition determined by the minimum element size and sound speed.
The speed of sound in ice is about 3,000 m/s, so for an element size of 0.01 m, $\Delta t \approx 3 \times 10^{-6}$ s, which is very short.
Calculating "Ice Floes" with DEM (Discrete Element Method) โ Large-Scale Icebreaking Impossible with FEM
In actual icebreaker navigation, countless broken ice fragments continue to flow around the hull. Modeling each of these "ice fragment floes" individually with FEM (Finite Element Method) would result in millions to tens of millions of elements, making computation impossible. This is where DEM (Discrete Element Method) is used. In DEM, each ice fragment is represented by a simple "rigid body + spring + dashpot" model, efficiently calculating contact forces between ice fragments and between fragments and the hull. Norway's SINTEF research institute has conducted simulations containing over 1 million ice fragments, matching the resistance of icebreakers with measurements within ยฑ15%. DEM is widely used not only for ice but also for analyzing sand, rock, and granular materials, and its development for ice-structure problems contributes to the overall technological advancement of granular mechanics.
Ice-Structure Interaction in Practice
Model Construction Procedure
Please tell me the steps to start an ice-structure interaction simulation.
1. Create a 3D FE model of the structure (for steel structures, use shell elements)
2. Create an ice plate model (solid elements. Add erosion settings if handling failure)
3. Define contact (surface-to-surface contact, set friction coefficient)
4. Set initial velocity / drift velocity of ice
5. Define material models (ice: e.g., Tsai-Wu failure criterion; structure: elasto-plastic)
6. Run with explicit method
Ice Material Parameters
How are the material properties of ice determined?
Sea ice properties strongly depend on temperature, salinity, and strain rate.
| Parameter | First-year ice (-10ยฐC) | Multi-year ice (-10ยฐC) |
|---|---|---|
| Young's modulus | 3โ9 GPa | 5โ10 GPa |
| Compressive strength | 2โ10 MPa | 5โ15 MPa |
| Tensile strength | 0.5โ2 MPa | 1โ3 MPa |
| Poisson's ratio | 0.33 | 0.33 |
| Density | 900 kg/mยณ | 910 kg/mยณ |
The property variation is large.
Therefore, parametric studies are essential. ISO 19906 specifies using characteristic values with a 50-year return period. Probabilistic evaluation using Monte Carlo simulation is also sometimes performed.
Icebreaker Design โ Optimizing the "Ice Breaking Process" with Simulation
Icebreakers don't simply "collide and break" ice; they induce "bending failure" in the ice through the inclination angle of the bow to break it efficiently. If the angle is too shallow, the ice doesn't break and slides under the hull; if it's too steep, it pushes the ice in compression, generating extremely high loads. The optimal angle varies with ice thickness and strength, so in practice, simulations are run under design conditions like "continuously breaking sea ice of 1.5m thickness and 2MPa compressive strength" to optimize the combination of bow shape and propulsion. In the design of Russia's Arktika-class (nuclear-powered icebreakers), simulation results directly led to changes in bow shape, reportedly improving icebreaking capability by 20% over the initial design.
Related Topics
Experience the theory firsthand with the interactive simulator for this field
All Simulators