skip to main content

This content will become publicly available on June 26, 2023

Title: Characterization of Fracture Process Zone in Short Beam Compression Tests on Barre Granite

Due to rock mass being commonly subjected to compressive or shear loading, the mode II fracture toughness is an important material parameter for rocks. Fracturing in rocks is governed by the behavior of a nonlinear region surrounding the crack tip called the fracture process zone (FPZ). However, the characteristics of mode II fracture are still determined based on the linear elastic fracture mechanics (LEFM), which assumes that a pure mode II loading results in a pure mode II fracture. In this study, the FPZ development in Barre granite specimens under mode II loading was investigated using the short beam compression (SBC) test. Additionally, the influence of lateral confinement on various characteristics of mode II fracture was studied. The experimental setup included the simultaneous monitoring of surface deformation using the two-dimensional digital image correlation technique (2D-DIC) to identify fracture mode and characterize the FPZ evolution in Barre granite specimens. The 2D-DIC analysis showed a dominant mixed-mode I/II fracture in the ligament between two notches, irrespective of confinement level on the SBC specimens. The influence of confinement on the SBC specimens was assessed by analyzing the evolution of crack displacement and changes in value of mode II fracture toughness. Larger levels of more » damage in confined specimens were observed prior to the failure than the unconfined specimens, indicating an increase in the fracture resistance and therefore mode II fracture toughness with the confining stress.


The fracturing in laboratory-scale rock specimens is often characterized by the deformation of the inelastic region surrounding the crack tips, also known as the fracture process zone (FPZ) (Backers et al., 2005; Ghamgosar and Erarslan, 2016). While the influence of the FPZ on mode I fracture in rocks has been extensively investigated, there are limited studies on FPZ development in rocks under pure mode II loading (Ji et al., 2016; Lin et al., 2020; Garg et al., 2021; Li et al., 2021).

« less
; ;
Award ID(s):
Publication Date:
Journal Name:
56th U.S. Rock Mechanics/Geomechanics Symposium
Sponsoring Org:
National Science Foundation
More Like this
  1. Fracturing in brittle rocks exhibits a significant nonlinear region surrounding the crack tip called the fracture process zone (FPZ). In this study, the evolution of the FPZ under pure mode II loading using notched deep beam under three-point loading was investigated. The experimental setup included the simultaneous monitoring of surface deformation using the two-dimensional digital image correlation technique to characterize various crack characteristics such as its type and FPZ evolution in Barre granite specimens. Both displacement and strain approaches of the two-dimensional digital image correlation were used to identify the mode of fracture under pure mode II loading. Both approaches showed that the crack initiation occur under mode I despite the pure mode II loading at the notch tip. The displacement approach was used for characterizing the evolution of the FPZ which analyzed the crack tip opening displacement and crack tip sliding displacement to identify the transition between the three stages of FPZ evolution, namely, (a) elastic stage, (b) formation of the FPZ, and (c) the macro-crack initiation. The results showed that the evolution of the FPZ of mode I fracture under pure mode II loading is similar to cases of pure mode I loading of the same rock.
  2. Fracturing in brittle rocks with an existing crack results in the development of a significant nonlinear region surrounding the crack tip called the fracture process zone. Various experimental and numerical studies have shown that the crack tip parameters such as the crack tip opening displacement (CTOD) and the fracture energy are critically important in characterizing the fracture process zone. In this study, numerical simulations of rock specimens with a center notch subjected to three-point bending were conducted using the extended finite element method (XFEM) along with the cohesive zone model (CZM) to account for fracture process zone. The input parameters of CZM such as the elastic and critical crack opening displacements were first estimated based on the results of three-point bending tests on the center notched Barre granite specimens. Displacements were measured using the two dimensional digital image correlation technique and used to characterize the evolution of the fracture process zone and estimate the parameters of the cohesive zone model. The results from the numerical simulations showed that CZM provided a good agreement with experimental data as it predicted all three stages of cracking from fracture process initiation to macro-crack growth.

    The chemo-mechanical loading of rocks causes the dissolution and precipitation of multiple phases in the rock. This dissolution and precipitation of load-bearing mineral phases lead to the stress redistribution in neighboring phases, which in turn results in deformational changes of the sample composite. The aim of this study is to investigate the link between microstructural evolution and creep behavior of shale rocks subjected to chemo-mechanical loading through modeling time-dependent deformation induced by the dissolution-precipitation process. The model couples the microstructural evolution of the shale rocks with the stress/strain fields inside the material as a function of time. The modeling effort is supplemented with an experimental study where shale rocks were exposed to CO2-rich brine under high temperature and pressure conditions. 3D snapshots of the sample microstructure were generated using segmented micro-CT images of the shale sample. The time-evolving microstructures were then integrated with the Finite element-based mechanical model to simulate the creep induced by dissolution and precipitation processes independent of the intrinsic viscoelasticity/viscoplasticity of the mineral phases. After computation of the time-dependent viscoelastic properties of the shale composite, the combined microstructure model and finite element model were utilized to predict the time-dependent stress and strain fields in different zones ofmore »reacted shale.


    Determination of viscous behavior of shale rocks is key in wide range of applications such as stability of reservoirs, stability of geo-structures subjected to environmental forcing, underground storage of hazardous materials and hydraulic fracturing. Short-term creep strains in hydraulic fracturing can change stress fields and in turn can impact the hydraulic fracturing procedures(H. Sone & Zoback, 2010; Hiroki Sone & Zoback, 2013). While long-term creep strains can hamper the reservoir performance due to the reduction in permeability of the reservoir by closing of fractures and fissures(Du, Hu, Meegoda, & Zhang, 2018; Rybacki, Meier, & Dresen, 2016; Sharma, Prakash, & Abedi, 2019; Hiroki Sone & Zoback, 2014). Owing to these significance of creep strain, it is important to understand the viscoelastic/viscoplastic behavior of shales.

    « less
  4. We observed and modeled the elastic, inelastic and time-dependent viscous properties of damaged Berea Sandstone samples to investigate the impact of damage on the rheological properties of rocks. Cylindrical Berea Sandstone plugs were prepared both parallel and perpendicular to bedding. We impacted the samples with Split Hopkinson Pressure Bar to pervasively fracture the specimens at different strain rates. Longitudinal mode-I fractures are dominant in specimens impacted at relatively low strain rates (about 130 /s), whereas shear fractures also form in specimens deformed at high strain rates (up to 250 /s). The damaged rocks were subjected to multiple steps of differential stress loading and hold stages under 15 MPa confining pressure. A key observation is that higher damaged specimens showed greater axial and volumetric creep strain deformation during loading and hold stages. Poisson ratio also increase with increasing damage. We modeled the volumetric strain of the sandstone specimens using a Perzyna viscoplasticity law that employs the Modified Cam Clay model as the yield criterion (Haghighat et al. 2020). We deduced that fractured rocks undergo substantial bulk time-dependent deformation due to volumetric compaction and fracture closure. Damage increase results in decrease of the effective viscosity of the material.
  5. The line crack models, including linear elastic fracture mechanics (LEFM), cohesive crack model (CCM), and extended finite element method (XFEM), rest on the century-old hypothesis of constancy of materials’ fracture energy. However, the type of fracture test presented here, named the gap test, reveals that, in concrete and probably all quasibrittle materials, including coarse-grained ceramics, rocks, stiff foams, fiber composites, wood, and sea ice, the effective mode I fracture energy depends strongly on the crack-parallel normal stress, in-plane or out-of-plane. This stress can double the fracture energy or reduce it to zero. Why hasn’t this been detected earlier? Because the crack-parallel stress in all standard fracture specimens is negligible, and is, anyway, unaccountable by line crack models. To simulate this phenomenon by finite elements (FE), the fracture process zone must have a finite width, and must be characterized by a realistic tensorial softening damage model whose vectorial constitutive law captures oriented mesoscale frictional slip, microcrack opening, and splitting with microbuckling. This is best accomplished by the FE crack band model which, when coupled with microplane model M7, fits the test results satisfactorily. The lattice discrete particle model also works. However, the scalar stress–displacement softening law of CCM and tensorial modelsmore »with a single-parameter damage law are inadequate. The experiment is proposed as a standard. It represents a simple modification of the three-point-bend test in which both the bending and crack-parallel compression are statically determinate. Finally, a perspective of various far-reaching consequences and limitations of CCM, LEFM, and XFEM is discussed.

    « less