Elastic–Plastic Fracture Propagation Modeling in Rock Fracturing via Punch Through Shear Test

Fracture initiation and propagation from a wellbore within a rock formation exhibit nonlinear and inelastic behaviors. When the rock material undergoes plastic deformation prior to failure, the classical Griffith theory is no longer valid. In this study, a variational phase-field approach is applied to model the inelastic behavior of granite rock in a punch through shear test. The rock failure and the fracture initiation and propagation during the loading was simulated and compared to the corresponding experimental investigations. In this numerical approach, the total local free energy is fully coupled with solid deformation and computes the plastic strain rate. The code is scripted in Multiphysics Object Oriented Simulation Environment (MOOSE). The model is shown capable of reproducing the evidenced phenomena from Punch Through Shear (PTS) test encompassing mixed mode fracture pattern Mode I, and Mode II. The numerical results show a good agreement in the stress–displacement curve with experimental data for the critical energy release rate of Gc=600N/m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{c}=600\mathrm{N}/\mathrm{m}$$\end{document}. Therefore, the granite sample’s fracture toughness for Mode II is calculated to be 4.85 MPam\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{MPa}\sqrt{\mathrm{m}}$$\end{document} at no confining pressure.