More Like this

1. Random field realization and fracture simulation of rocks with angular bias for fracture strength

   Garrard, J. M. ; Abedi, Reza ; Clarke, P. L. ( January 2018 , 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. 

   more » « less
2. Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks

   Abedi, Reza ; Clarke, Philip L. ( January 2018 , 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. 

   more » « less
3. Comparison of Interfacial and Continuum Models for Dynamic Fragmentation Analysis

   https://doi.org/10.1115/IMECE2018-88294  

   Bahmani, Bahador ; Clarke, Philip ; Abedi, Reza ( November 2018 , 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.

    

   more » « less
4. Modeling of rock inhomogeneity and anisotropy by explicit and implicit representation of microcracks

   Abedi, R. ; Clarke, P.L. ( January 2018 , Proceeding 52th U.S. Rock Mechanics/Geomechanics Symposium)

   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
5. Controlling the fracture response of structures via topology optimization: From delaying fracture nucleation to maximizing toughness

   https://doi.org/10.1016/j.jmps.2023.105227  

   Jia, Yingqi ; Lopez-Pamies, Oscar ; Zhang, Xiaojia Shelly ( April 2023 , Journal of the Mechanics and Physics of Solids)

   It is now a well-established fact that even simple topology variations can drastically change the fracture response of structures. With the objective of gaining quantitative insight into this phenomenon, this paper puts forth a density-based topology optimization framework for the fracture response of structures subjected to quasistatic mechanical loads. One of the two key features of the proposed framework is that it makes use of a complete phase-field fracture theory that has been recently shown capable of accurately describing the nucleation and propagation of brittle fracture in a wide range of nominally elastic materials under a wide range of loading conditions. The other key feature is that the framework is based on a multi-objective function that allows optimizing in a weighted manner: ( ) the initial stiffness of the structure, ( ) the first instance at which fracture nucleates, and ( ) the energy dissipated by fracture propagation once fracture nucleation has occurred. The focus is on the basic case of structures made of a single homogeneous material featuring an isotropic linear elastic behavior alongside an isotropic strength surface and toughness. Novel interpolation rules are proposed for each of these three types of material properties. As a first effort to gain quantitative insight, the framework is deployed to optimize the fracture response of 2D structures wherein the fracture is bound to nucleate in three different types of regions: within the bulk, from geometric singularities (pre-existing cracks and sharp corners), and from smooth parts of the boundary. The obtained optimized structures are shown to exhibit significantly enhanced fracture behaviors compared to those of structures that are optimized according to conventional stiffness maximization. Furthermore, the results serve to reveal a variety of strengthening and toughening mechanisms. These include the promotion of highly porous structures, the formation of tension-compression asymmetric regions, and the removal of cracks and sharp corners. The particular mechanism that is preferred by a given structure, not surprisingly, correlates directly to the elastic, strength, and toughness properties of the material that is made of. 

   more » « less