Phase-Field Modeling of Fracture under Compression and Confinement in Anisotropic Geomaterials

TL;DR

This paper develops an advanced phase-field fracture model for anisotropic geomaterials under compression, incorporating anisotropic elasticity and damage responses to accurately predict complex crack patterns observed in experiments.

Contribution

It extends existing phase-field models to anisotropic materials with compressive traction, enabling realistic simulation of fracture processes in layered geomaterials.

Findings

Accurately predicts crack growth patterns in layered specimens.

Matches experimental observations of shale fracture under compression.

Enables efficient simulation of wing crack growth in anisotropic media.

Abstract

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.