skip to main content


Title: Microstructural Controls on Mixed Mode Dynamic Fracture Propagation in Crystalline and Porous Granular Rocks
Abstract

Brittle fracture propagation in rocks is a complex process due to significant grain‐scale heterogeneity and evolving stress states under dynamic loading conditions. In this work, we use digital image correlation and linear elastic fracture mechanics to make instantaneous measurements of the opening (mode I) and in plane shear (mode II) components of the stress intensity field during dynamic mixed mode crack initiation and propagation in crystalline and granular rocks. Both rock types display some similar fracture behaviors as observed in engineered materials, including rate dependent fracture initiation toughness and a direct relationship between propagation toughness and crack velocity; however, measured propagation toughness is higher than quasi‐static values at crack velocities well below the branching velocity in both rocks. Additionally, due to grain scale controls on the fracture process, mixed mode crack propagation is fundamentally different between these two rock types. Mixed mode propagation is energetically more favorable than pure opening mode propagation in sandstone, while the opposite is true in granite. Furthermore, following initiation, propagation in granite occurs so as to minimize the mode II contribution, irrespective of the initiation conditions, while fractures in sandstone maintain a non‐negligible mode II contribution during propagation across the sample.

 
more » « less
NSF-PAR ID:
10368388
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
127
Issue:
3
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT:

    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 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.

    1. INTRODUCTION

    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).

     
    more » « less
  2. SUMMARY

    A better understanding of damage accumulation before dynamic failure events in geological material is essential to improve seismic hazard assessment. Previous research has demonstrated the sensitivity of seismic velocities to variations in crack geometry, with established evidence indicating that initial crack closure induces rapid changes in velocity. Our study extends these findings by investigating velocity changes by applying coda wave interferometry (CWI). We use an array of 16 piezoceramic transducers to send and record ultrasonic pulses and to determine changes in seismic velocity on intact and faulted Westerly granite samples. Velocity changes are determined from CWI and direct phase arrivals. This study consists of three sets of experiments designed to characterize variations in seismic velocity under various initial and boundary conditions. The first set of experiments tracks velocity changes during hydrostatic compression from 2 and 191 MPa in intact Westerly granite samples. The second set of experiments focuses on saw-cut samples with different roughness and examines the effects of confining pressure increase from 2 to 120 MPa. The dynamic formation of a fracture and the preceding damage accumulation is the focus of the third type of experiment, during which we fractured an initially intact rock sample by increasing the differential stress up to 780 MPa while keeping the sample confined at 75 MPa. The tests show that: (i) The velocity change for rough saw cut samples suggests that the changes in bulk material properties have a more pronounced influence than fault surface apertures or roughness. (ii) Seismic velocities demonstrate higher sensitivity to damage accumulation under increasing differential stress than macroscopic measurements. Axial stress measured by an external load cell deviates from linearity around two-third through the experiment at a stress level of 290 MPa higher than during the initial drop in seismic velocities. (iii) Direct waves exhibit strong anisotropy with increasing differential stress and accumulating damage before rock fracture. Coda waves, on the other hand, effectively average over elastic wave propagation for both fast and slow directions, and the resulting velocity estimates show little evidence for anisotropy. The results demonstrate the sensitivity of seismic velocity to damage evolution at various boundary conditions and progressive microcrack generation with long lead times before dynamic fracture.

     
    more » « less
  3. null (Ed.)
    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. 
    more » « less
  4. null (Ed.)
    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. 
    more » « less
  5. Abstract

    The operation of fracture, diffusion, and intracrystalline‐plastic micromechanisms during semibrittle deformation of rock is directly relevant to understanding mechanical behavior across the brittle‐plastic transition in the crust. An outstanding question is whether (1) the micromechanisms of semibrittle flow can be considered to operate independently, as represented in typical crustal strength profiles across the brittle to plastic transition, or (2) the micromechanisms are coupled such that the transition is represented by a distinct rheology with dependency on effective pressure, temperature, and strain rate. We employ triaxial stress‐cycling experiments to investigate elastic‐plastic and viscoelastic behaviors during semibrittle flow in two distinctly different monomineralic, polycrystalline, synthetic salt‐rocks. During semibrittle flow at high differential stress, granular, low‐porosity, work‐hardened salt‐rocks deform predominantly by grain‐boundary sliding and wing‐crack opening accompanied by minor intragranular dislocation glide. In contrast, fully annealed, near‐zero porosity salt‐rocks flow at lower differential stress by intragranular dislocation glide accompanied by grain‐boundary sliding and opening. Grain‐boundary sliding is frictional during semibrittle flow at higher strain rates, but the associated dispersal of water from fluid inclusions along boundaries can activate fluid‐assisted diffusional sliding at lower strain rates. Changes in elastic properties with semibrittle flow largely reflect activation of sliding along closed grain boundaries. Observed microstructures, pronounced hysteresis and anelasticity during cyclic stressing after semibrittle flow, and stress relaxation behaviors indicate coupled operation of micromechanisms leading to a distinct rheology (hypothesis 2 above).

     
    more » « less