skip to main content


This content will become publicly available on July 21, 2024

Title: Brittle to ductile transition during compression of glassy nanoparticles studied in molecular dynamics simulations
Understanding how nanoparticles deform under compression is not only of scientific importance, but also has practical significance in various applications such as tribology, nanoparticle-based probes, and dry grinding of raw materials. In this study, we conducted compression tests on model brittle glassy nanoparticles using molecular dynamics simulations. We found that during the early stages of plastic deformation, shear bands formed in a similar pattern regardless of nanoparticle size. However, as the deformation continued, dominant cracks emerged in large nanoparticles while being suppressed in smaller ones. This size-dependent brittle to ductile transition can be explained by a simple model based on Griffith's theory. We also investigated the effect of surface stress state on fracture using thermally tempered nanoparticles. We observed that the presence of compressive surface stress strengthened the nanoparticle by suppressing crack formation, even when a pre-notch was present. On the other hand, tensile surface stress had the opposite effect. Interestingly, nanoparticles with both tensile and compressive surface stress promoted shear deformation, which could potentially compromise the mechanical performance of tempered glass despite delayed crack formation.  more » « less
Award ID(s):
1936368
NSF-PAR ID:
10492394
Author(s) / Creator(s):
; ;
Publisher / Repository:
AIP Publishing
Date Published:
Journal Name:
Journal of Applied Physics
ISSN:
0021-8979
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The discrete damage model presented in this paper accounts for 42 non-interacting crack microplanes directions. At the scale of the representative volume element, the free enthalpy is the sum of the elastic energy stored in the non-damaged bulk material and in the displacement jumps at crack faces. Closed cracks propagate in the pure mode II, whereas open cracks propagate in the mixed mode (I/II). The elastic domain is at the intersection of the yield surfaces of the activated crack families, and thus describes a non-smooth surface. In order to solve for the 42 crack densities, a Closest Point Projection algorithm is adopted locally. The representative volume element inelastic strain is calculated iteratively using the Newton–Raphson method. The proposed damage model was rigorously calibrated for both compressive and tensile stress paths. Finite element method simulations of triaxial compression tests showed that the transition between brittle and ductile behavior at increasing confining pressure can be captured. The cracks’ density, orientation, and location predicted in the simulations are in agreement with experimental observations made during compression and tension tests, and accurately show the difference between tensile and compressive strength. Plane stress tension tests simulated for a fiber-reinforced brittle material also demonstrated that the model can be used to interpret crack patterns, design composite structures and recommend reparation techniques for structural elements subjected to multiple damage mechanisms. 
    more » « less
  2. ABSTRACT:

    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 of 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.

    1. INTRODUCTION

    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).

     
    more » « less
  3. Abstract Carbonate sediments play a prominent role on the global geological stage as they store more than $$60\%$$ 60 % of world’s oil and $$40\%$$ 40 % of world’s gas reserves. Prediction of the deformation and failure of porous carbonates is, therefore, essential to minimise reservoir compaction, fault reactivation, or wellbore instability. This relies on our understanding of the mechanisms underlying the observed inelastic response to fluid injection or deviatoric stress perturbations. Understanding the impact of deformation/failure on the hydraulic properties of the rock is also essential as injection/production rates will be affected. In this work, we present new experimental results from triaxial deformation experiments carried out to elucidate the behaviour of a porous limestone reservoir analogue (Savonnières limestone). Drained triaxial and isotropic compression tests were conducted at five different confining pressures in dry and water-saturated conditions. Stress–strain data and X-ray tomography images of the rock indicate two distinct types of deformation and failure regimes: at low confinement (10 MPa) brittle failure in the form of dilatant shear banding was dominant; whereas at higher confinement compaction bands orthogonal to the maximum principal stress formed. In addition to the pore pressure effect, the presence of water in the pore space significantly weakened the rock, thereby shrinking the yield envelope compared to the dry conditions, and shifted the brittle–ductile transition to lower effective confining pressures (from 35 MPa to 29 MPa). Finally, permeability measurements during deformation show a reduction of an order of magnitude in the ductile regime due to the formation of the compaction bands. These results highlight the importance of considering the role of the saturating fluid in the brittle–ductile response of porous rocks and elucidate some of the microstructural processes taking place during this transition. 
    more » « less
  4. Abstract

    Crack interactions leading to shear localization were quantified using microstructural analysis for brittle faults and high‐temperature ductile faults formed during experiments on quartz sandstone. In both faulting regimes, the nucleation of macroscopic faults results from the interactions of microfractures at two length scalesin ensemble. Brittle faults nucleate when the longest mesoscale shear fractures and transgranular tensile cracks critically interact. In contrast, ductile faults nucleate when the longest mesoscale shear fractures and multi‐grain scale intergranular shear cracks critically interact. For both faulting regimes, we conclude the interaction and coalescence of the longest mesoscale shear fractures is the fundamental process responsible for fault nucleation. Hence, mesoscale shear fractures, which accommodate the majority of axial strain prior to shear localization in both faulting regimes, also serve as the nucleus of macroscopic faults. Locally, the growth of the mesoscale shear fractures is promoted by the interaction and coalescence of the multi‐grain scale cracks in both faulting regimes. We hypothesize that attainment of a critical microstructure for shear localization (i.e., local clustering of the longest microfractures) requires a characteristic amount of plastic axial strain, which depends on deformation conditions. In brittle faulting, distributed microfracturing is confined within limited regions of the rock volume, which expedites crack clustering and fault nucleation at low characteristic strains. In ductile faulting, distributed microfracturing occurs more uniformly throughout the rock volume, delaying shear localization to high characteristic strains. Accurate prediction of shear localization requires models that describe crack interactions of thelargestflaws that account for crackclustering.

     
    more » « less
  5. Abstract

    Experimentally quantifying the viscoplastic rheology of olivine at the high stresses and low temperatures of the shallow lithosphere is challenging due to olivine's propensity to deform by brittle mechanisms at these conditions. In this study, we use microscale uniaxial compression tests to investigate the rheology of an olivine single crystal at room pressure and temperature. Pillars with nominal diameters of 1.25 μm were prepared using a focused ion beam milling technique and were subjected to sustained axial stresses of several gigapascal. The majority of the pillars failed after dwell times ranging from several seconds to a few hours. However, several pillars exhibited clear evidence of plastic deformation without failure after 4–8 hr under load. The corresponding creep strain rates are estimated to be on the order of 10−6to 10−7 s−1. The uniaxial stresses required to achieve this deformation (4.1–4.4 GPa) are in excellent agreement with complementary data obtained using nanoindentation techniques. Scanning transmission electron microscopy observations indicate that deformation occurred along amorphous shear bands within the deformed pillars. Electron energy loss spectroscopy measurements revealed that the bands are enriched in Fe and depleted in Mg. We propose that inhomogeneities in the cation distribution in olivine concentrate stress and promote the amorphization of the Fe‐rich regions. The time dependence of catastrophic failure events suggests that the amorphous bands must grow to some critical length scale to generate an unstable defect, such as a shear crack.

     
    more » « less