skip to main content

Title: Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks
Fracture in rock as a heterogeneous brittle material, having significant inherent randomness, requires including probabilistic considerations at different scales. Crack growth in rocks is generally associated with complex features such as crack path oscillations, microcrack and crack branching events. Two methods will be presented to address rock inhomogeneity and anisotropy. First, microcracks are explicitly realized in a domain based on specific statistics of crack length and location. Second, a statistical model is used to implicitly represent an inhomogeneous field for fracture strength. Both approaches can be used for rocks in which the natural fractures are oriented in a specific angle, i.e. an aspect for modeling bedding planes in sedimentary rocks.
Award ID(s):
Publication Date:
Journal Name:
Proceedings of ARMA 2018
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.
  2. 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.

    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 ofmore »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
  4. Natural arches are culturally valued rock landforms common in sedimentary rocks of the Colorado Plateau and additionally occur broadly around the world. Recent notable collapses of some of these landforms have highlighted the need to better understand the mechanics of their failure. While environmentally driven weathering has been the focus of most previous studies of arch collapse, comparably little attention has been given to anthropogenic vibration sources and how these often slight- to moderate-magnitude shaking events might steadily weaken arches over time. We collected 12–15 months of continuous ambient vibration data from arches and nearby bedrock in both anthropogenically ‘noisy’ and ‘quiet’ locations and used these datasets to develop an annual model of arch peak ground velocity based on magnitude-cumulative frequency distributions. Working from these models, we added vibration events of varying magnitude or frequency of occurrence, informed by field data, imitating arch vibration in response to different anthropogenic activities such as helicopter flights or induced earthquakes. We then applied subcritical fracture mechanics principles to predict annual crack growth rates in an idealized arch under these different vibration conditions. Our results demonstrate that in a single year, cracks grow minimally longer (∼1%) in ‘noisy’ environments than in areas not experiencing anthropogenicmore »vibration energy. Few (1+) 30-s moderate-magnitude events (∼15 mm/s) or many (>37,000) 30-s low-magnitude events (∼2 mm/s) cause markedly increased crack growth. Our approach provides a valuable new framework for assessing the range and frequency of occurrence of vibrations experienced by an arch, and for predicting arch damage. Our results, in turn, yield important new outputs applicable in support of conservation management of these and similar landforms world-wide under exposure to a range of human-induced vibration activity.« less
  5. International Ocean Discovery Program (IODP) Expedition 357 successfully cored an east–west transect across the southern wall of Atlantis Massif on the western flank of the Mid-Atlantic Ridge to study the links between serpentinization processes and microbial activity in the shallow subsurface of highly altered ultramafic and mafic sequences that have been uplifted to the seafloor along a major detachment fault zone. The primary goals of this expedition were to (1) examine the role of serpentinization in driving hydrothermal systems, sustaining microbial communities, and sequestering carbon; (2) characterize the tectonomagmatic processes that lead to lithospheric heterogeneities and detachment faulting; and (3) assess how abiotic and biotic processes change with variations in rock type and progressive exposure on the seafloor. To accomplish these objectives, we developed a coring and sampling strategy based around the use of seabed rock drills—the first time that such systems have been used in the scientific ocean drilling programs. This technology was chosen in hopes of achieving high recovery of the carbonate cap sequences and intact contact and deformation relationships. The expedition plans also included several engineering developments to assess geochemical parameters during drilling; sample bottom water before and after drilling; supply synthetic tracers during drilling for contaminationmore »assessment; gather downhole electrical resistivity and magnetic susceptibility logs for assessing fractures, fluid flow, and extent of serpentinization; and seal boreholes to provide opportunities for future experiments. Seventeen holes were drilled at nine sites across Atlantis Massif, with two sites on the eastern end of the southern wall (Sites M0068 and M0075), three sites in the central section of the southern wall north of the Lost City hydrothermal field (Sites M0069, M0072, and M0076), two sites on the western end (Sites M0071 and M0073), and two sites north of the southern wall in the direction of the central dome of the massif and Integrated Ocean Drilling Program Site U1309 (Sites M0070 and M0074). Use of seabed rock drills enabled collection of more than 57 m of core, with borehole penetration ranging from 1.3 to 16.44 meters below seafloor and core recoveries as high as 75% of total penetration. This high level of recovery of shallow mantle sequences is unprecedented in the history of ocean drilling. The cores recovered along the southern wall of Atlantis Massif have highly heterogeneous lithologies, types of alteration, and degrees of deformation. The ultramafic rocks are dominated by harzburgites with intervals of dunite and minor pyroxenite veins, as well as gabbroic rocks occurring as melt impregnations and veins, all of which provide information about early magmatic processes and the magmatic evolution in the southernmost portion of Atlantis Massif. Dolerite dikes and basaltic rocks represent the latest stage of magmatic activity. Overall, the ultramafic rocks recovered during Expedition 357 revealed a high degree of serpentinization, as well as metasomatic talc-amphibole-chlorite overprinting and local rodingitization. Metasomatism postdates an early phase of serpentinization but predates late-stage intrusion and alteration of dolerite dikes and the extrusion of basalt. The intensity of alteration is generally lower in the gabbroic and doleritic rocks. Chilled margins in dolerite intruded into talc-amphibole-chlorite schists are observed at the most eastern Site M0075. Deformation in Expedition 357 cores is variable and dominated by brecciation and formation of localized shear zones; the degree of carbonate veining was lower than anticipated. All types of variably altered and deformed ultramafic and mafic rocks occur as components in sedimentary breccias and as fault scarp rubble. The sedimentary cap rocks include basaltic breccias with a carbonate sand matrix and/or fossiliferous carbonate. Fresh glass on basaltic components was observed in some of the breccias. The expedition also successfully applied new technologies, namely (1) extensively using an in situ sensor package and water sampling system on the seabed drills for evaluating real-time dissolved oxygen and methane, pH, oxidation-reduction potential, temperature, and conductivity during drilling; (2) deploying a borehole plug system for sealing seabed drill boreholes at four sites to allow access for future sampling; and (3) proving that tracers can be delivered into drilling fluids when using seabed drills. The rock drill sensor packages and water sampling enabled detection of elevated dissolved methane and hydrogen concentrations during and/or after drilling, with “hot spots” of hydrogen observed over Sites M0068–M0072 and methane over Sites M0070–M0072. Shipboard determination of contamination tracer delivery confirmed appropriate sample handling procedures for microbiological and geochemical analyses, which will aid all subsequent microbiological investigations that are part of the science party sampling plans, as well as verify this new tracer delivery technology for seabed drill rigs. Shipboard investigation of biomass density in select samples revealed relatively low and variable cell densities, and enrichment experiments set up shipboard reveal growth. Thus, we anticipate achieving many of the deep biosphere–related objectives of the expedition through continued scientific investigation in the coming years. Finally, although not an objective of the expedition, we were serendipitously able to generate a high-resolution (20 m per pixel) multibeam bathymetry map across the entire Atlantis Massif and the nearby fracture zone, Mid-Atlantic Ridge, and eastern conjugate, taking advantage of weather and operational downtime. This will assist science party members in evaluating and interpreting tectonic and mass-wasting processes at Atlantis Massif.« less