skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Deformation and Frictional Failure of Granular Media in 3D Analog and Numerical Experiments
Abstract Frictional sliding along grain boundaries in brittle shear zones can result in the fragmentation of individual grains, which ultimately can impact slip dynamics. During deformation at small scales, stick–slip motion can occur between grains when existing force chains break due to grain rearrangement or failure, resulting in frictional sliding of granular material. The rearrangement of the grains leads to dilation of the granular package, reducing the shear stress and subsequently leading to slip. Here, we conduct physical experiments employing HydroOrbs, an elasto-plastic material, to investigate grain comminution in granular media under simple shear conditions. Our findings demonstrate that the degree of grain comminution is dependent on both the normal force and the size of the grains. Using the experimental setup, we benchmark Discrete Element Method (DEM) numerical models, which are capable of simulating the movement, rotation, and fracturing of elasto-plastic grains subjected to simple shear. The DEM models successfully replicate both grain comminution patterns and horizontal force fluctuations observed in our physical experiments. They show that increasing normal forces correlate with higher horizontal forces and more fractured grains. The ability of our DEM models to accurately reproduce experimental results opens up new avenues for investigating various parameter spaces that may not be accessible through traditional laboratory experiments, for example, in assessing how internal friction or cohesion affect deformation in granular systems.  more » « less
Award ID(s):
1843676
PAR ID:
10502454
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Pure and Applied Geophysics
Volume:
181
Issue:
7
ISSN:
0033-4553
Format(s):
Medium: X Size: p. 2083-2105
Size(s):
p. 2083-2105
Sponsoring Org:
National Science Foundation
More Like this
  1. Deformation experiments on hematite characterize its slip‐rate dependent frictional properties and deformation mechanisms. These data inform interpretations of slip behavior from exhumed hematite‐coated faults and present‐day deformation at depth. We used a rotary‐shear apparatus to conduct single‐velocity and velocity‐step experiments on polycrystalline specular hematite rock (∼17 μm average plate thickness) at slip rates of 0.85 μm/s to 320 mm/s, displacements of primarily 1–3 cm and up to 45 cm, and normal stresses of 5 and 8.5 MPa. The average coefficient of friction is 0.70; velocity‐step experiments indicate velocity‐strengthening to velocity‐neutral behavior at rates <1 mm/s. Scanning electron microscopy showed experimentally generated faults develop in a semi‐continuous, thin layer of red hematite gouge. Angular gouge particles have an average diameter of ∼0.7 μm, and grain size reduction during slip yields a factor of 10–100 increase in surface area. Hematite is amenable to (U‐Th)/He thermochronometry, which can quantify fault‐related thermal and mechanical processes. Comparison of hematite (U‐Th)/He dates from the undeformed material and experimentally produced gouge indicates He loss occurs during comminution at slow deformation rates without an associated temperature rise required for diffusive loss. Our results imply that, in natural fault rocks, deformation localizes within coarse‐grained hematite by stable sliding, and that hematite (U‐Th)/He dates acquired from ultracataclasite or highly comminuted gouge reflect minor He loss unrelated to thermal processes. Consequently, the magnitude of temperature rise and associated thermal resetting in hematite‐bearing fault rocks based on (U‐Th)/He thermochronometry may be overestimated if only diffusive loss of He is considered. 
    more » « less
  2. Abstract Localized frictional sliding on faults in the continental crust transitions at depth to distributed deformation in viscous shear zones. This brittle‐ductile transition (BDT), and/or the transition from velocity‐weakening (VW) to velocity‐strengthening (VS) friction, are controlled by the lithospheric thermal structure and composition. Here, we investigate these transitions, and their effect on the depth extent of earthquakes, using 2D antiplane shear simulations of a strike‐slip fault with rate‐and‐state friction. The off‐fault material is viscoelastic, with temperature‐dependent dislocation creep. We solve the heat equation for temperature, accounting for frictional and viscous shear heating that creates a thermal anomaly relative to the ambient geotherm which reduces viscosity and facilitates viscous flow. We explore several geotherms and effective normal stress distributions (by changing pore pressure), quantifying the thermal anomaly, seismic and aseismic slip, and the transition from frictional sliding to viscous flow. The thermal anomaly can reach several hundred degrees below the seismogenic zone in models with hydrostatic pressure but is smaller for higher pressure (and these high‐pressure models are most consistent with San Andreas Fault heat flow constraints). Shear heating raises the BDT, sometimes to where it limits rupture depth rather than the frictional VW‐to‐VS transition. Our thermomechanical modeling framework can be used to evaluate lithospheric rheology and thermal models through predictions of earthquake ruptures, postseismic and interseismic crustal deformation, heat flow, and the geological structures that reflect the complex deformation beneath faults. 
    more » « less
  3. Typical granular materials are far from optimal in terms of mechanical performance: Random packing leads to poor load transfer in the form of thin and dispersed force lines within the material, to inhomogeneous jamming, and to strain localization. In addition, localized contacts between individual grains result in low stiffness, strength and brittleness. Here we propose a granular material that simultaneously embodies three approaches to increase strength: geometrical design of individual grains, crystallization, and infiltration by an adhesive. Using mechanical vibrations, we assembled millimeter-scale 3D printed grains with rhombic dodecahedral shapes into fully dense FCC granular crystals. We then infiltrated the granular structure with a tacky, polyacrylic adhesive that is orders of magnitude weaker than the grains, but which provides sustained adhesion over large interfacial displacements. The resulting material is a fully dense, free-standing space filling granular crystal. Compressive tests show that these granular crystals are up to 60 times stronger than randomly packed cohesive spheres and they display a rich set of mechanisms: Nonlinear deformations, crystal plasticity reminiscent of atomistic mechanisms, cross-slip, shear-induced dilatancy, micro-buckling, and tensile strength. To capture some of these mechanisms we developed a multiscale model that incorporates local cohesion between grains, resolved shear and normal stresses on available slip planes, and prediction of compressive strength as function of loading orientation. The predicted strength is highly anisotropic and agrees well with the compression experiments. Once fully understood and harnessed, we envision that these mechanisms will lead to granular engineering materials with unusual combinations of mechanical performances attractive for many applications. 
    more » « less
  4. We report on an open-source, publicly available C++ software module called Chrono::GPU, which uses the Discrete Element Method (DEM) to simulate large granular systems on Graphics Processing Unit (GPU) cards. The solver supports the integration of granular material with geometries defined by triangle meshes, as well as co-simulation with the multi-physics simulation engine Chrono. Chrono::GPU adopts a smooth contact formulation and implements various common contact force models, such as the Hertzian model for normal force and the Mindlin friction force model, which takes into account the history of tangential displacement, rolling frictional torques, and cohesion. We report on the code structure and highlight its use of mixed data types for reducing the memory footprint and increasing simulation speed. We discuss several validation tests (wave propagation, rotating drum, direct shear test, crater test) that compare the simulation results against experimental data or results reported in the literature. In another benchmark test, we demonstrate linear scaling with a problem size up to the GPU memory capacity; specifically, for systems with 130 million DEM elements. The simulation infrastructure is demonstrated in conjunction with simulations of the NASA Curiosity rover, which is currently active on Mars. 
    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