Abstract: Granular materials exist widely in nature and engineering fields. Particles may break under external loads. The breakage behavior of granular materials will cause changes in their physical and mechanical properties, which will have a great effect on engineering and construction. The study of the breakage behavior of particles at macro and micro scales can not only reveal the mechanical mechanism of particle breakage, but also provide guarantee for the safety and normal operation of the engineering field. Therefore, the analysis of particle breakage process has both practical engineering significance and theoretical research value. In this paper, the numerical analysis methods of particle breakage behavior are reviewed. Among the numerical methods based on discrete element method (DEM), the bonded particle model (BPM) and fragment replacement method (FRM) are introduced. In the method based on discrete element-finite element coupling algorithm, the proportional boundary finite element method, combined finite element-discrete element method (FDEM) and cohesive zone model (CZM) are introduced. In addition, the Peridynamics (PD) method is introduced in detail. The proposal, implementation process, development, key issues, advantages and disadvantages of the above numerical methods are summarized and discussed. In addition, the research results at home and abroad and the main applications in engineering are reviewed for each numerical method, and the key issues of each method are introduced. Finally, the current numerical research on particle breakage is summarized, and the future development direction is briefly prospected.