Material point method to simulate the evolution characteristics of loading damage in fractured sandstone

Cracks in rocks can lead to different load responses and mechanical properties. Primary and secondary cracks are altered during loading, with varying degrees of impact on the failure of underground engineering rock masses. However, the precise modeling of fracture initiation and propagation resulting in failure remains challenging, and the macroscopic evolution characteristics of cracks and stress fields have not been accurately described. In this study, a numerical simulation tool based on the material point method (MPM) and a strain softening model were used to establish a plane loading model with different cutting angles. A simulation of the damage evolution of fractured sandstone under constant-rate loading was conducted for different confining pressures and cutting angles. The results show the following. 1) A comparison of the simulation and experimental results indicates that the proposed numerical method accurately simulates the macroscopic crack generation and the complex physical processes of fractured sandstone during loading. 2) The stress concentration area originates from both ends of the crack and extends toward the edge in a diagonal direction until the crack has penetrated the rock. 3) The ranking of the cutting angle based on the sandstone strength is 0(degrees) > 90(degrees) > 60(degrees) > 30(degrees). The cumulative dissipated energy exhibits an inflection point near the peak stress, and the increase after this point is primarily attributed to friction at the crack surface. The damage degree and dissipated energy of the rock in different stages correspond to the rock's strength.