In the extreme event such as air-blast or explosion, strong shocks lead to severe damage and fragmentation in structures. Despite the considerable effort made in recent years, reliable numerical prediction of fragmentation processes in materials and solids under blast loading or shock wave remains challenging. The conventional mesh-based methods (eg, finite element method (FEM)) are ineffective due to large deformation-induced mesh distortion issues and exhibit non-convergent solutions in fracture problems. The meshfree method, such as reproducing kernel particle method (RKPM), naturally avoids computational difficulties associated with low-quality meshes, allows efficient adaptive refinement, and provides flexible control of smoothness and locality in numerical approximations.