During icebreaker navigation in ice-covered waters, icebreaking resistance and dynamic ice loads acting on the bow critically determine the vessel’s icebreaking performance. Quantitative characterization of the icebreaking resistance behavior and ice load distribution on the bow is essential for elucidating ship-ice interaction mechanisms, assessing icebreaking capability, and optimizing structural design. This study conducted comparative icebreaking tests on two icebreaker bow models with distinct geometries in the small ice model basin of China Ship Scientific Research Center (CSSRC SIMB). Systematic measurements were performed to quantify icebreaking resistance, capture spatiotemporal ice load distributions, and document ice failure patterns under level ice conditions. The analysis reveals that bow geometry profoundly influences icebreaking efficiency: the stem angle governs the proportion of bending failure during vertical ice penetration, while the flare angle modulates circumferential failure modes along the hull-ice interface. Notably, the sunken keel configuration enhances ice clearance by mechanically expelling fractured ice blocks. Ice load distributions exhibit pronounced nonlinearity, with localized pressure concentrations and stochastic load center migration driven by ice fracture dynamics. Furthermore, icebreaking patterns—such as fractured ice dimensions and kinematic behavior during ship-ice interaction—are quantitatively correlated with the bow designs. These experimentally validated findings provide critical insights into ice-structure interaction physics, offering an empirical foundation for performance prediction and bow-form optimization in the preliminary design of icebreakers.

Wave propagation is impacted by the presence of ice floes. The influence of waves, on the other hand, causes ice floes to overlap and accumulate. In this paper, the interaction of ice floes and regular waves was simulated using the Finite Element Method. Firstly, natural ice floe fields were generated using the Python 3.10 programing language, with floe size distribution and randomness taken into consideration. Then, using the velocity inlet boundary wave generation method, regular simple harmonic waves were produced. The process where ice floes couple with waves was simulated with the Coupled Eulerian–Lagrangian (CEL) approach. Variations in wave height after passing through the ice floe field were investigated, and further research was conducted on the movement and fragmentation characteristics of ice floes. Simulations employing the Coupled Eulerian–Lagrangian (CEL) approach reveal that (1) ice floe motion exhibits periodic characteristics synchronized with incident wave periods; (2) wave height attenuation increases by 62–80% with rising ice concentration (70–90%); and (3) fragmentation predominantly occurs at wave trough phases due to flexural stress concentration. These findings quantitatively characterize wave–ice energy transfer mechanisms critical for polar navigation safety assessments.

The crushing failure of sea ice is a critical design issue for polar offshore structures and ship structures because ice-induced loads may generate pronounced local damage and dynamic responses. Accurately modelling this process remains challenging because ice crushing involves localized fragmentation, crack propagation, rubble accumulation, and repeated contact release. This paper presents a controlled numerical sensitivity study of level-ice crushing against a vertical structure using a coupled LS-DYNA framework that combines the Structured Arbitrary Lagrangian–Eulerian (S-ALE) formulation with the Cohesive Element Method (CEM). The study focuses on a benchmark-scale indentation configuration and examines how mesh topology, mesh size, and imposed indentation velocity affect the predicted fracture morphology and load-time histories. The results show that random triangular meshes better reproduce stochastic fragmentation and lateral flaking than regular triangular or quadrilateral meshes, while finer meshes reduce excessive load oscillations and provide more stable force histories. The velocity study indicates a transition from gradual crushing and fragment retention at lower velocities to more rapid brittle chipping and stronger dynamic fluctuations at higher velocities. A benchmark-level comparison with published ice-indentation simulations shows that the predicted peak line load is of the same order of magnitude as reference results. The proposed framework is therefore useful for investigating numerical sensitivities and failure-mode trends in ice-crushing simulations, although final design-load application requires further calibration and formal mesh-independence assessment.