Abstract
In this work, we present a peridynamic-based simulation method for modeling quasi-static fracture propagation in isotropic and anisotropic rock within the framework of peridynamic least square minimization (PDLSM). The original isotropic elastic PDLSM is further extended to investigate fracture propagation in anisotropic materials in this study. The proposed AN-PDLSM model integrates an anisotropic model in fracture mechanics to analyze the failure process of anisotropic rocks. An important advancement in this research lies in the incorporation of the maximum energy release rate criterion (MERR) into the PDLSM framework for the first time. This enhancement enables accurately determining crack propagation and the associated crack angles. The proposed model utilizes the energy release rate calculated through the J-integral method to assess bond breakage, and it employs a mesh-independent, piecewise linear fracture model to describe crack propagation. The proposed method fully combines the merits of traditional fracture mechanics with the unique capabilities of peridynamics. To demonstrate the effectiveness of the proposed model, a simulation of fracture evolution in isotropic plates subjected to semi-circular bending tests is presented and compared with experimental results. It is shown that the proposed model accurately replicates fracture trajectories in isotropic specimens. In the context of anisotropic rock, the effect of a weak coefficient on crack morphology is discussed in order to obtain a suitable value. Additionally, the impact of bedding angles on fracture paths through our proposed model is also explored, revealing excellent agreement with experimental results. The findings in this research demonstrate that the proposed AN-PDLSM model is exceptionally proficient at capturing the intricate, oscillating crack paths observed in anisotropic rock materials.