Simulation of Shear and Tensile Fractures Using Ductile Phase Field Modelling with the Calibration of P Wave Velocity Measurement and Moment Tensor Inversion

Xu Li,Guangyao Si,Joung Oh,Ismet Canbulat
DOI: https://doi.org/10.1007/s00603-024-03880-z
IF: 6.518
2024-04-29
Rock Mechanics and Rock Engineering
Abstract:The appropriate understanding and formulation of rock discontinuities via FEM is still challenging for rock engineering, as continuous algorithms cannot handle the discontinuities in rock mass. Also, different failure modes of rock samples, containing tensile and shear failure, need to be computed separately. In this study, a novel double-phase field damage model was introduced with two independent phase field damage variables. The construction of the proposed model follows the thermodynamics framework from the overall Helmholtz free energy, with elastic, plastic and surface damage components. The proposed model is calibrated via traditional damage variables, based on ultrasonic wave velocity measurement and acoustic emission monitoring, and both show great consistency between simulation results and laboratory observations. Then the double-phase field damage model is applied to COMSOL software to simulate microcrack propagation in a pre-fractured rock sample. Both lateral and wing cracks are observed in this study, manifested as shear- and tensile-dominated cracks. We also observed different microcracking mechanisms in the proposed numerical models, such as tensile and shear cracking, the influence of plastic strain and the percolation between tensile and shear microcracks. Overall, this study provides valuable insights into the mechanics of microcracking in rocks, and the proposed model shows promising results in simulating crack propagation.
geosciences, multidisciplinary,engineering, geological
What problem does this paper attempt to address?
This paper attempts to address the problem of accurately simulating shear and tensile fractures in rocks using the Finite Element Method (FEM) in rock mass engineering. Traditional methods face difficulties in handling discontinuities in rocks, especially in simultaneously simulating different failure modes (such as tensile and shear failures). Therefore, this paper proposes a new dual-phase field damage model that includes two independent phase field damage variables to simulate tensile and shear microcracks, respectively. The model is calibrated through ultrasonic velocity measurements and acoustic emission monitoring, and numerical simulations are conducted in COMSOL software to verify its effectiveness in simulating the propagation of microcracks in pre-cracked rock samples. ### Main Issues 1. **Accurate simulation of rock discontinuities**: Traditional FEM methods cannot effectively handle discontinuities in rocks, especially when simulating tensile and shear failures. 2. **Model calibration**: Current phase field damage models lack effective calibration methods, particularly in the reliable selection of length scale parameters. 3. **Distinguishing tensile and shear cracks**: Existing phase field damage models fail to adequately distinguish between tensile and shear cracks, while rock materials exhibit significant differences in resistance to these two types of cracks. ### Solutions 1. **Dual-phase field damage model**: A new dual-phase field damage model with two independent phase field damage variables is proposed to simulate tensile and shear microcracks, respectively. 2. **Model calibration**: The model is calibrated through ultrasonic velocity measurements and acoustic emission monitoring to ensure consistency between model parameters and experimental data. 3. **Numerical simulation**: Numerical simulations are conducted in COMSOL software to verify the model's effectiveness in simulating the propagation of microcracks in pre-cracked rock samples. ### Main Contributions 1. **Model innovation**: A new dual-phase field damage model with two independent phase field damage variables is proposed, which can more accurately simulate tensile and shear microcracks in rocks. 2. **Calibration method**: Length scale parameters are rigorously selected based on laboratory observation data, including complete stress-strain curves, ultrasonic velocity measurements, and acoustic emission data. 3. **Numerical validation**: The model's effectiveness is validated through numerical simulations, observing shear and tensile dominated microcracks and different microcrack mechanisms, such as the penetration phenomena of tensile and shear cracks. ### Keywords - Elasto-plastic model - Plastic fracture propagation - Phase field damage model - Ultrasonic - Acoustic emission This study provides valuable insights into the accurate simulation of the mechanical behavior of rock microcracks in rock mass engineering and demonstrates the potential of the proposed model in simulating crack propagation.