Аннотация:The development of equipment capable of operating in the Arctic environment could benefit from a mathematical model of interaction with ice. In an attempt to create such a model, a numerical simulation of a collision between a rigid body and an ice field was conducted. Discontinuous Galerkin method was used to solve the equations modeling an elastoplastic body, specifically the Prandtl–Reuss model with the Mises–Schleichert yield condition. The nature of destruction and stress distribution in ice mass were analyzed, and then the momentum change from the collision was used to estimate ice resistance. The results of this analysis show the applicability of the chosen method to the problem of ice destruction.