The article provides numerical simulation of ice fracture by compression using of the discrete element method by the software Yade. The Bonded Particle Model with damage was used. The computational calibration of the material parameters was performed on the test model of the sample in the form of a hyperboloid, the difference in strength by tensile and compression was 3.26 times. A comparison of real and computational tests on compression of cylindrical samples with diameter and height 50 mm at a temperature of -10 °С was carried out. The deformation curves are similar in the general character, reactions level, and deformations. The method makes it possible to simulate the ice fracture. According to the estimation of the deformation work, the difference between the typical experimental curve and the calculated is about 2 %.