Abstract: Abstract：[Objectives] In order to study the distribution of ship ice pressure during the interaction between ship and ice breaking, a numerical simulation analysis was carried out for polar ships in the presence of the ice-breaking channel. [Methods] The discrete element method was used to model the ship and the ice-breaking unit. It was assumed that the ice-breaking was composed of an ideal two-dimensional discrete, and the effects of ocean currents on the buoyancy, drag and additional mass of the ice-breaking unit were taken into account. The numerical model was validated by the experimental results of model of MT Uikku to compare and analyze the effects of different speeds and ice conditions on the ice pressure in the ship area. [Results] The results showed that the ice loads were mainly concentrated on the bow of the ship when the ship moved in the ice-broken channel. Under the comparative analysis of different ice conditions and speed, the ice pressure in the bow area increased with the increase of ice thickness, speed and ice broken unit density, among which ice thickness was the most important factor affecting the amplitude of ice pressure. The comparative analysis under different ice conditions and speeds indicated that the ice pressure in the ice-breaking area had the greatest influence on the ice pressure in the bow area, and the ice pressure in the bow transition area has a significant effect. [Conclusions] The numerical analysis method can provide some certain references for the safe navigation and structural design of polar ships.