Search

Strongly anisotropic geomaterials undergo fracture under compressive loading. This paper applies a phase-field fracture model to study this fracture process. While phase-field fracture models have several advantages, they provide unphysical predictions when the stress state is complex and includes compression that can cause crack faces to contact. Building on a phase-field model that accounts for compressive traction across the crack face, this paper extends the model to anisotropic fracture. The key features include: (1) a homogenized anisotropic elastic response and strongly-anisotropic model for the work to fracture; (2) an effective damage response that accounts consistently for compressive traction across the crack face, that is derived from the anisotropic elastic response; (3) a regularized crack normal field that overcomes the shortcomings of the isotropic setting, and enables the correct crack response, both across and transverse to the crack face. To test the model, we first compare the predictions to phase-field fracture evolution calculations in a fully-resolved layered specimen with spatial inhomogeneity, and show that it captures the overall patterns of crack growth. We then apply the model to previously-reported experimental observations of fracture evolution in laboratory specimens of shales under compression with confinement, and find that it predicts well the observed crack patterns in a broad range of loading conditions. We further apply the model to predict the growth of wing cracks under compression and confinement. The effective crack response model enables us to treat the initial crack simply as a non-singular damaged zone within the computational domain, thereby allowing for easy and general computations., Comment: To appear in International Journal for Numerical and Analytical Methods in Geomechanics

Phase-field models are a leading approach for realistic fracture problems. They treat the crack as a second phase and use gradient terms to smear out the crack faces, enabling the use of standard numerical methods for simulations. This regularization causes cracks to occupy a finite volume in the reference, and leads to the inability to appropriately model the closing or contacting -- without healing -- of crack faces. Specifically, the classical idealized crack face tractions are that the shear component is zero, and that the normal component is zero when the crack opens and identical to the intact material when the crack closes. Phase-field fracture models do not replicate this behavior. This work addresses this shortcoming by introducing an effective crack energy density that endows the regularized (finite volume) phase-field crack with the effective properties of an idealized sharp crack. The approach is based on applying the QR (upper triangular) decomposition of the deformation gradient tensor in the basis of the crack, enabling a transparent identification of the crack deformation modes. By then relaxing over those modes that do not cost energy, an effective energy is obtained that has the intact response when the crack faces close and zero energy when the crack faces are open. A highlight of this approach is that it lies completely in the setting of finite deformation, enabling potential application to soft materials and other settings with large deformation or rotations. The model is applied to numerically study representative complex loadings, including (1) cyclic loading on a cavity in a soft solid that shows the growth and closing of cracks in complex stress states; and (2) cyclic shear that shows a complex pattern of crack branching driven by the closure of cracks.