Title: Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks
Fracture patterns experienced under a dynamic uniaxial compressive load are highly sensitive to rock microstructural defects due to its brittleness and the absence of macroscopic stress concentration points. We propose two different approaches for modeling rock microstructural defects and inhomogeneity. In the explicit realization approach, microcracks with certain statistics are incorporated in the computational domain. In the implicit realization approach, fracture strength values are sampled using a Weibull probability distribution. We use the Mohr-Coulomb failure criterion to define an effective stress in the context of an interfacial damage model. This model predicts crack propagation at angles ±ɸch = ±(45 − ɸ/2) relative to the direction of compressive load, where ɸ is the friction angle. By using appropriate models for fracture strength anisotropy, we demonstrate the interaction of rock weakest plane and ɸch. Numerical results demonstrate the greater effect of strength anisotropy on fracture pattern when an explicit approach is employed. In addition, the density of fractures increases as the angle of the weakest planes approaches ±ɸch. The fracture simulations are performed by an h-adaptive asynchronous spacetime discontinuous Galerkin (aSDG) method that can accommodate crack propagation in any directions. more »« less
Abedi, Reza; Clarke, Philip L.
(, Proceedings of ARMA 2018)
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.
Garrard, J. M.; Abedi, Reza; Clarke, P. L.
(, 52nd U.S. Rock Mechanics/Geomechanics Symposium)
Realistic fracture simulations in rock as a heterogeneous brittle material with significant inherent randomness require the use of models that incorporate its inhomogeneities and statistical variability. The high dependence of their fracture progress on microstructural defects results in wide scatter in their ultimate strength and the so-called size effect. This paper proposes an approach based on statistical volume elements (SVEs) to characterize rock fracture strength at the mesoscale. The use of SVEs ensures that the material randomness is maintained upon averaging of microscale features. Because the fracture strength varies not just spatially, but also by the angle of loading, this work includes angular variability to properly model a heterogeneous rock domain. Two different microcrack distributions, one angularly uniform and one angularly biased towards a specific angle, are used to show that implementing angle into the random field provides the most realistic fracture simulation. An adaptive asynchronous spacetime discontinuous Galerkin (aSDG) finite element method is used to perform the dynamic fracture simulations.
Garrard, J.M.; Abedi, R.; Clarke, P.L.
(, Proceedings of ARMA 2018)
Realistic fracture simulations in rock as a heterogeneous brittle material with significant inherent ran- domness require the use of models that incorporate its inhomogeneities and statistical variability. The high dependence of their fracture progress on microstructural defects results in wide scatter in their ultimate strength and the so-called size effect. This paper proposes an approach based on statistical volume elements (SVEs) to characterize rock fracture strength at the mesoscale. The use of SVEs ensures that the material randomness is maintained upon averaging of microscale features. Because the fracture strength varies not just spatially, but also by the angle of loading, this work includes angular variability to properly model a heterogeneous rock domain. Two different microcrack distributions, one angularly uniform and one angularly biased towards a specific angle, are used to show that implementing angle into the random field provides the most realistic fracture simulation. An adaptive asynchronous spacetime discontinuous Galerkin (aSDG) finite element method is used to perform the dynamic fracture simulations.
Bahmani, Bahador; Clarke, Philip; Abedi, Reza
(, Proceedings of ASME 2018 International Mechanical Engineering Congress and Exposition)
The microstructural design has an essential effect on the fracture response of brittle materials. We present a stochastic bulk damage formulation to model dynamic brittle fracture. This model is compared with a similar interfacial model for homogeneous and heterogeneous materials. The damage models are rate-dependent, and the corresponding damage evolution includes delay effects. The delay effect provides mesh objectivity with much less computational efforts. A stochastic field is defined for material cohesion and fracture strength to involve microstructure effects in the proposed formulations. The statistical fields are constructed through the Karhunen-Loeve (KL) method. An advanced asynchronous Spacetime Discontinuous Galerkin (aSDG) method is used to discretize the final system of coupled equations. Application of the presented formulation is shown through dynamic fracture simulation of rock under a uniaxial compressive load. The final results show that a stochastic bulk damage model produces more realistic results in comparison with a homogenizes model.
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.
Abedi, R., and Clarke, P.L. Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks. Retrieved from https://par.nsf.gov/biblio/10117623. Proceeding 52th U.S. Rock Mechanics/Geomechanics Symposium .
Abedi, R., & Clarke, P.L. Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks. Proceeding 52th U.S. Rock Mechanics/Geomechanics Symposium, (). Retrieved from https://par.nsf.gov/biblio/10117623.
Abedi, R., and Clarke, P.L.
"Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks". Proceeding 52th U.S. Rock Mechanics/Geomechanics Symposium (). Country unknown/Code not available. https://par.nsf.gov/biblio/10117623.
@article{osti_10117623,
place = {Country unknown/Code not available},
title = {Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks},
url = {https://par.nsf.gov/biblio/10117623},
abstractNote = {Fracture patterns experienced under a dynamic uniaxial compressive load are highly sensitive to rock microstructural defects due to its brittleness and the absence of macroscopic stress concentration points. We propose two different approaches for modeling rock microstructural defects and inhomogeneity. In the explicit realization approach, microcracks with certain statistics are incorporated in the computational domain. In the implicit realization approach, fracture strength values are sampled using a Weibull probability distribution. We use the Mohr-Coulomb failure criterion to define an effective stress in the context of an interfacial damage model. This model predicts crack propagation at angles ±ɸch = ±(45 − ɸ/2) relative to the direction of compressive load, where ɸ is the friction angle. By using appropriate models for fracture strength anisotropy, we demonstrate the interaction of rock weakest plane and ɸch. Numerical results demonstrate the greater effect of strength anisotropy on fracture pattern when an explicit approach is employed. In addition, the density of fractures increases as the angle of the weakest planes approaches ±ɸch. The fracture simulations are performed by an h-adaptive asynchronous spacetime discontinuous Galerkin (aSDG) method that can accommodate crack propagation in any directions.},
journal = {Proceeding 52th U.S. Rock Mechanics/Geomechanics Symposium},
author = {Abedi, R. and Clarke, P.L.},
}
Warning: Leaving National Science Foundation Website
You are now leaving the National Science Foundation website to go to a non-government website.
Website:
NSF takes no responsibility for and exercises no control over the views expressed or the accuracy of
the information contained on this site. Also be aware that NSF's privacy policy does not apply to this site.
An official website of the United States government
