This paper covers the development of prediction techniques for ice load on ice-breakers operating in continuous ice-breaking under level ice conditions using particle-based continuum mechanics. Ice is assumed to be a linear elastic material until the fracture occurs. The maximum normal stress theory is used for the criterion of fracture. The location of the crack can be expressed using a local scalar function consisting of the gradient of the first principal stress and the corresponding eigen-vector. This expression is used to determine the relative position of particle pair to the new crack. The Hertz contact model is introduced to consider the collisions between ice fragments and the collisions between hull and ice fragments. In order to verify the developed technique, the simulation results for the three-point bending problems of ice-specimen and the continuous ice-breaking problem around a wedge-type model ship with bow angle of 20°are compared with the experimental results carrying out at Korea Research Institute of Ships and Ocean Engineering (KRISO).