At present, the calculation method of ice load on surface navigation ships has been very mature, but the calculation method of submarine ice load is very few. The reasonable evaluation of submarine ice load has become an urgent problem to be solved. In this paper, the mechanical characteristics of the submarine surfacing ice-breaking process are systematically analyzed. Based on the theory of plate and shell, the theoretical calculation models of ice-breaking resistance of the submarine command tower and hull are established, respectively, and the ice load calculation method of the command tower and the hull is obtained. Then, the submarine model SUBOFF is used to perform the numerical simulation of the submarine’s ice-breaking and surfacing process. The numerical result is compared with the ice-breaking resistance calculation model. The results show that the ice-breaking resistance calculation model proposed in this paper is consistent with the numerical simulation results, and the influence of parameters such as ice mechanical properties, the upper area of the command tower, and initial crack length on ice-breaking resistance is established. The calculation model in this paper can provide a theoretical reference for the optimization design of polar submarine structures.

The increased demand for performing crane ship modeling has led to the necessity for fast and accurate numerical experiments. This paper presents an approach for creating a numerical model of a crane vessel with a suspended load that allows for real-time control of crane parts. The model is developed in the MATLAB/Simulink environment, which makes it possible to extend it further to the user's needs. The authors describe the approach to the calculation of wave-induced ship motions, presents the Simulink block model and describes the features encountered during the simulation process. The possibility of real-time control of the position of crane parts is also shown, keeping the calculation speed of the ship hydrodynamics.

In recent years numerical modelling is widely used for solving ice-structure interaction problems. The most popular numerical method is finite element method. One of the promising approaches for modelling ice destruction is cohesive element method, where ice field is discretized using both bulk and special cohesive elements, which are possible crack paths. The investigation of modelling features and limitations of this method have been performed by several authors but with not appropriate quality. This paper belongs to series of numerical experiments performed to study the cohesive element techniques for ice modelling and includes some considerations and study of influence of finite element shape.