With the shrinking of Arctic sea ice due to global climate change, potential access to Arctic waters has increased for non-typical icebreaking or strengthened ships. Numerous studies have been conducted on hull form designs and ice resistance predictions for ships with typical icebreaking bows, but published research for ships with non-icebreaking bows in ice is still rare. The objective of this study was to investigate the ice resistance of a naval surface ship with a non-icebreaking bow through model tests in an ice tank. The naval surface combatant concept DTMB 5415 was used as the ship model. The tests were conducted under different levels of ice thicknesses and speeds. During the tests, the total resistance of the model ship was measured, accompanied by monitoring of the ice load at the stem area with a flexible tactile sensor sheet. Compared with the test results of icebreaker models in former studies, the total ice resistance, as well as the stem ice load, of the present ship was significantly higher. The ice crushing resistance component in the stem area accounted for more than 60% of the total resistance in the ice. Discussions on the applicability of a semi-empirical formula for predicting the ice resistance of the present ship are also presented. Keinonen’s formula was found to be relatively more consistent with the predictions produced by model tests, and a preliminary modification was proposed to obtain more accurate predictions.

This study presents a 2D numerical simulation method for predicting the resistance of ships navigating in pack ice. The key contribution of this study lies in the derivation of analytical closed-form solutions for calculating the flexural deformation and stress distribution in an elastic plate using Symplectic Mechanics and Hooke’s laws. These solutions are used to determine the failure mode of ice floes. Linear Elastic Fracture Mechanics (LEFMs) and the weight function method are utilized to analyze crack initiation, propagation, and fracture. Ice is broken when a crack propagates to 14.5% of the ice length. The compressive strength of ice and the contact area are used to calculate the ice load. A collision method was developed based on the Sweep and Prune (SAP) and Gilbert–Johnson–Keerthi (GJK) algorithms. A program for predicting the resistance of ships navigating in pack ice was developed based on MATLAB and the aforementioned theories. The navigation resistance of RV Xuelong at different ice concentrations and speeds was simulated and compared with the model test results from an ice tank. The comparison shows that the simulation results are consistent with the test results, with an average error of 9.05%, indicating the effectiveness and reliability of this numerical method. This study lays a solid foundation for future research on autonomous ship navigation in pack ice.

The effective assessment of risk and speed limitations in ice are critical for the route planning of polar cruise ships. While the Polar Operational Limit Assessment Risk Indexing System (POLARIS) is widely used to evaluate the operational risk in ice, its scope of assessment is limited to the ship’s assigned ice class. For a specific ship with a given ice class, the propulsion performance under varying ice conditions is more essential for assessing the besetting risk and finding an optimal route. To establish a more detailed risk evaluation method for a PC6 class polar cruise ship, propulsion performance under various ice conditions is obtained via model tests in an ice tank. During the tests, the tow force, propeller torque, and thrust are measured under different ship velocities and propeller rotation rates, and the relations between required delivered power (P) and ship velocity (V) under the tested ice thicknesses and concentrations are obtained and extended to other ice conditions by curve fitting. A new risk index outcome, RIO*, is proposed after POLARIS according to the optimized load ranges of the ship’s rated power. Four risk levels, including low, medium, high, and unnavigable, are classified with the required propulsion power in ice being 50%, 85%, 100%, and >100% of the rated power, respectively. The recommended speed limit for each risk level is proposed based on the economical service speed of 11 kn, the operational limit of 3 kn by POLARIS and the minimum speed of 0.5 kn to avoid besetting, respectively. Based on the RIO* and P–V relations, the speed map for varying ice thicknesses and concentrations can be calculated. On these bases, a route planning simulation for the present polar cruise ship is performed. Results of the case study show that 29% of the simulated area is identified as “unnavigable” by the present procedure, while only 8% is detected by POLARIS.