A PD-FEM Coupling Approach for Modeling Cracks Propagation in Brittle Rock under Compressive Load

Considering the low computational efficiency and accuracy of peridynamics (PD), a concurrent multiscale method coupling PD and finite element method (FEM) is proposed for modeling crack propagation of brittle rock under compressive load. In the coupling method, the fracture behavior is solved by the ordinary state-based peridynamics (OSBPD), while the elastic deformation of the rock mass is simulated by FEM. The implementation of the approach is as the following steps: first, a hybrid region is introduced in the framework to realize the strain energy equivalence between PD and FEM and eliminate the boundary effect. Then, the short-range force is utilized in PD to model the contact of the crack surface and prevent the particles from penetrating each other. In addition, the tangential force of the short-range force is introduced to simulate the friction sliding effect of the crack surface. Finally, the dynamic relaxation method is used to solve the displacement in the PD-FEM coupling model. The crack propagation of rock samples with a single pre-existing closed fracture under uniaxial compression is simulated by the PD-FEM coupling approach, and the numerical calculation results are in good agreement with the rule of the experimental results. The proposed coupling approach can capture the failure process of rocks under uniaxial compression and reduce the computational cost, simultaneously.