skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, February 13 until 2:00 AM ET on Friday, February 14 due to maintenance. We apologize for the inconvenience.


Title: Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures
Abstract

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

 
more » « less
Award ID(s):
1753031
PAR ID:
10537558
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Springer
Date Published:
Journal Name:
Applied Mathematics and Mechanics
Volume:
44
Issue:
7
ISSN:
0253-4827
Page Range / eLocation ID:
1125 to 1150
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Many geo‐engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase‐field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase‐field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy‐based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase‐field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase‐field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. The comparison of numerical results for two‐dimensional test cases with existing analytical solutions demonstrates that the proposed phase‐field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three‐dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

     
    more » « less
  2. In this paper, material deformation during ultra-precision machining (UPM) on the C-, R-, and A-planes of sapphire was investigated using the slip/fracture activation model where the likelihood of activation of individual plastic deformation and fracture systems on different crystallographic planes was calculated. The stress data obtained from molecular dynamics (MD) simulations were utilized, and the slip/fracture activation model was developed by incorporating the principal stresses in calculating the plastic deformation and fracture cleavage parameters. The analysis methodology was applied to study material deformation along various cutting orientations in sapphire. The stress field at crack initiation during UPM on C-, R-, and A-planes of sapphire was calculated using molecular dynamics (MD) simulations. An equation describing the relationship between crack initiation and its triggering parameters was formulated considering the systems’ plastic deformation and cleavage fractures. The model can qualitatively predict the crack initiations for various cutting orientations. The proposed model was verified through ultra-precision orthogonal plunge cut experiments along the same cutting orientations as in the MD simulations. 
    more » « less
  3. 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 models 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.

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

    This work proposes a novel approach for coupling non‐isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid‐filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface‐capturing phase‐field fracture method and interface‐tracking thermo‐fluid‐structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase‐field approach to represent fractures and reconstruct the geometry to frame a thermo‐fluid‐structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase‐field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid‐filled fracture with the phase‐field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase‐field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

     
    more » « less