Determination of ice reaction to dynamic loading is a crucial step for the solution of many practical problems in the Arctic region. However, this task is complicated by the time-dependent structural changes in ice that lead to localization of the deformations near the area of load application. This work is focused on the commonly occurring low-speed collisions, which are studied based on the laboratory experiment consisting of a direct impact by a spherical indentor. In order to reconstruct ice complex behavior, it is geometrically divided into three zones: a semi-spherical hydrostatic core beneath the contact patch included in a spherical elastoplastic zone, and the surrounding area of elasticity. Moreover, the development of fractures is taken into account. The main goal is the investigation of three (viscoelastic, elastoviscoplastic and sand) rheological models for the core, where such elaborate processes as microfracturing and recrystallization can take place. As a result of the calculations, wave and stress patterns are analyzed. Quantitative comparison is conducted based on the striker’s velocity and coordinate. The simulations reconstructed non-linear waves and different types of destruction, such as side and median cracks, cone fractures, inner and surface spalling. The obtained results can help to get a better understanding of ice rheology.