Coarse graining atomistic simulations of plastically deforming amorphous solids

Stress Relaxation and Creep Response of Glassy Hydrogels with Dense Physical Associations

Various glassy hydrogels are developed by forming dense physical associations within the matrices, which exhibit forced elastic deformation and possess high stiffness, strength, and toughness. Here, the viscoplastic behaviors of the glassy hydrogel of poly(methacrylamide-co-methacrylic acid) are investigated by stress relaxation and creep measurements. We found that the characteristic time of stress relaxation of the glassy gel is much smaller than that of amorphous polymers. The varying hydrogen bond strength leads to a broad distribution of structural activation energies, which in turn affects the range of characteristic time. In the presence of water, the weak hydrogen bond associations are easily disrupted under applied strain, enhancing segmental mobility and reducing relaxation time in the preyield regime, while in the postyield regime, the relaxation time increases slightly since the chain stretching increases the energy barrier. In creep tests, the creep strain rate accelerates at the initial stage due to stress-activated segments and then decelerates as chains are extensively stretched. The stress required for structural activation during creep is much lower than the Young's modulus of the gel, reflecting the poor structural stability. To further analyze the underlying mechanism of the glassy gel, a micromechanical model is established based on an extension on shear transformation zone theory. By incorporating a state variable for hydrogen bond density, this model can capture the intricate mechanical responses of glassy gels. Our findings reveal that glassy hydrogels are far from the thermodynamic equilibrium state, exhibiting rapid segment activation under external loading. This work provides insights to the dynamics and structural stability of glassy materials and can promote the design and applications of tough hydrogels.

Avalanches in Cu-Zr-Al metallic glasses

Metallic glasses have mechanical properties, which exhibit avalanches in the disguise of stress drops. We study these phenomena in a classical metallic glass system Cu-Zr-Al by athermal quasistatic shear and varying the element concentrations and for pure Cu-Zr 50/50 case the cooling rate. The resulting mechanical properties are close to the behavior found experimentally. At small strains, the pristine systems are akin to other glassy systems with a so-called gap distribution with a small positive exponent. Critical avalanching behavior is found only approaching the yield point. The post-yield stress drops are universal, and the gap distribution becomes flat.

Topological characterization of rearrangements in amorphous solids

In amorphous materials, plasticity is localized and occurs as shear transformations. It was recently shown by Wu et al. that these shear transformations can be predicted by applying topological defect concepts developed for liquid crystals to an analysis of vibrational eigenmodes [Z. W. Wu et al., Nat. Commun. 14, 2955 (2023)10.1038/s41467-023-38547-w]. This study relates the -1 topological defects to the displacement fields expected of an Eshelby inclusion, which are characterized by an orientation and the magnitude of the eigenstrain. A corresponding orientation and magnitude can be defined for each defect using the local displacement field around each defect. These parameters characterize the plastic stress relaxation associated with the local structural rearrangement and can be extracted using the fit to either the global displacement field or the local field. Both methods provide a reasonable estimation of the molecular-dynamics-measured stress drop, confirming the localized nature of the displacements that control both long-range deformation and stress relaxation.

Permanent shear localization in dense disordered materials due to microscopic inertia

In this work using computer simulations of 3D model of dense disordered solids we show, for the first time, the appearance of shear localization in the stationary flow under homogeneous driving conditions. To rationalize our simulation results we develop a continuum model, that couples the dynamics of the local flow to the evolution of a kinetic temperature field related to the local inertial dynamics. Our model predicts that the coupling of the flow field to this additional destabilizing field appears only as a necessary condition for shear localization, a minimum system size is necessary to accommodate the flow instability. Moreover we show that this size criterion resulting from our continuum description is in quantitative agreement with our particle-based simulation results.

Identifying microscopic factors that influence ductility in disordered solids

There are empirical strategies for tuning the degree of strain localization in disordered solids, but they are system-specific and no theoretical framework explains their effectiveness or limitations. Here, we study three model disordered solids: a simulated atomic glass, an experimental granular packing, and a simulated polymer glass. We tune each system using a different strategy to exhibit two different degrees of strain localization. In tandem, we construct structuro-elastoplastic (StEP) models, which reduce descriptions of the systems to a few microscopic features that control strain localization, using a machine learning-based descriptor, softness, to represent the stability of the disordered local structure. The models are based on calculated correlations of softness and rearrangements. Without additional parameters, the models exhibit semiquantitative agreement with observed stress-strain curves and softness statistics for all systems studied. Moreover, the StEP models reveal that initial structure, the near-field effect of rearrangements on local structure, and rearrangement size, respectively, are responsible for the changes in ductility observed in the three systems. Thus, StEP models provide microscopic understanding of how strain localization depends on the interplay of structure, plasticity, and elasticity.

Necking and failure of a particulate gel strand: signatures of yielding on different length scales

"Sticky" spheres with a short-ranged attraction are a basic model of a wide range of materials from the atomic to the granular length scale. Among the complex phenomena exhibited by sticky spheres is the formation of far-from-equilibrium dynamically arrested networks which comprise "strands" of densely packed particles. The aging and failure of such gels under load is a remarkably challenging problem, given the simplicity of the model, as it involves multiple length- and time-scales, making a single approach ineffective. Here we tackle this challenge by addressing the failure of a single strand with a combination of methods. We study the mechanical response of a single strand of a model gel-former to deformation, both numerically and analytically. Under elongation, the strand breaks by a necking instability. We analyse this behaviour at three different length scales: a rheological continuum model of the whole strand; a microscopic analysis of the particle structure and dynamics; and the local stress tensor. Combining these different approaches gives a coherent picture of the necking and failure. The strand has an amorphous local structure and has large residual stresses from its initialisation. We find that neck formation is associated with increased plastic flow, a reduction in the stability of the local structure, and a reduction in the residual stresses; this indicates that the system loses its solid character and starts to behave more like a viscous fluid. These results will inform the development of more detailed models that incorporate the heterogeneous network structure of particulate gels.

First-contact-breaking distributions in strained disordered crystals

We derive exact probability distributions for the strain (ε) at which the first stress drop event occurs in uniformly strained disordered crystals, with quenched disorder introduced through polydispersity in particle sizes. We characterize these first stress drop events numerically as well as theoretically and identify them with the first-contact-breaking event in the system. Our theoretical results are corroborated with numerical simulations of quasistatic volumetric strain applied to disordered near-crystalline configurations of athermal soft particles. We develop a general technique to determine the distribution of strains at which the first stress drop events occur, through an exact mapping between the cumulative distribution of first-contact-breaking events and the volume of a convex polytope whose dimension is determined by the number of defects N_{d} in the system. An exact numerical computation of this polytope volume for systems with small numbers of defects displays a remarkable match with the distribution of strains generated through direct numerical simulations. Finally, we derive the distribution of strains at which the first stress drop occurs, assuming that individual contact-breaking events are uncorrelated, which accurately reproduces distributions obtained from direct numerical simulations.

Elastoplastic Approach Based on Microscopic Insights for the Steady State and Transient Dynamics of Sheared Disordered Solids

We develop a framework to study the mechanical response of athermal amorphous solids via a coupling of mesoscale and microscopic models. Using measurements of coarse-grained quantities from simulations of dense disordered particulate systems, we present a coherent elastoplastic model approach for deformation and flow of yield stress materials. For a given set of parameters, this model allows us to match consistently transient and steady state features of driven disordered systems with diverse preparation histories under both applied shear-rate and creep protocols.

Coordinate transformation methodology for simulating quasistatic elastoplastic solids

Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypoelastoplastic mechanical model, and develop a projection method to enforce quasistatic equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and which imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.

Rejuvenation and shear banding in model amorphous solids

We measure the local yield stress, at the scale of small atomic regions, in a deeply quenched two-dimensional glass model undergoing shear banding in response to athermal quasistatic deformation. We find that the occurrence of essentially a single plastic event suffices to bring the local yield stress distribution to a well-defined value for all strain orientations, thus essentially erasing the memory of the initial structure. It follows that in a well-relaxed sample, plastic events cause the abrupt (nucleation-like) emergence of a local softness contrast and thus precipitate the formation of a band, which, in its early stages, is measurably softer than the steady-state flow. Moreover, this postevent yield stress ensemble presents a mean value comparable to that of the inherent states of a supercooled liquid around the mode-coupling temperature T_{MCT}. This, we argue, explains that the transition between brittle and ductile yielding in amorphous materials occurs around a comparable parent temperature. Our data also permit to capture quantitatively the contributions of pressure and density changes and demonstrate unambiguously that they are negligible compared with the changes of softness caused by structural rejuvenation.

Strain localization in dry sheared granular materials: A compactivity-based approach

Shear banding is widely observed in natural fault zones as well as in laboratory experiments on granular materials. Understanding the dynamics of strain localization under different loading conditions is essential for quantifying strength evolution of fault gouge and energy partitioning during earthquakes and characterizing rheological transitions and fault zone structure changes. To that end, we develop a physics-based continuum model for strain localization in sheared granular materials. The grain-scale dynamics is described by the shear transformation zone (STZ) theory, a nonequilibrium statistical thermodynamic framework for viscoplastic deformation in amorphous materials. Using a finite strain computational framework, we investigate the initiation and growth of complex shear bands under a variety of loading conditions and identify implications for strength evolution and the ductile to brittle transition. Our numerical results show similar localization patterns to field and laboratory observations and suggest that shear zones show more ductile response at higher confining pressures, lower dilatancy, and loose initial conditions. Lower pressures, higher dilatancy, and dense initial conditions favor a brittle response and larger strength drops. These findings shed light on a range of mechanisms for strength evolution in dry sheared granular materials and provide a critical input to physics-based multiscale models of fault zone instabilities.

Local yield stress statistics in model amorphous solids

We develop and extend a method presented by Patinet, Vandembroucq, and Falk [Phys. Rev. Lett. 117, 045501 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.045501] to compute the local yield stresses at the atomic scale in model two-dimensional Lennard-Jones glasses produced via differing quench protocols. This technique allows us to sample the plastic rearrangements in a nonperturbative manner for different loading directions on a well-controlled length scale. Plastic activity upon shearing correlates strongly with the locations of low yield stresses in the quenched states. This correlation is higher in more structurally relaxed systems. The distribution of local yield stresses is also shown to strongly depend on the quench protocol: the more relaxed the glass, the higher the local plastic thresholds. Analysis of the magnitude of local plastic relaxations reveals that stress drops follow exponential distributions, justifying the hypothesis of an average characteristic amplitude often conjectured in mesoscopic or continuum models. The amplitude of the local plastic rearrangements increases on average with the yield stress, regardless of the system preparation. The local yield stress varies with the shear orientation tested and strongly correlates with the plastic rearrangement locations when the system is sheared correspondingly. It is thus argued that plastic rearrangements are the consequence of shear transformation zones encoded in the glass structure that possess weak slip planes along different orientations. Finally, we justify the length scale employed in this work and extract the yield threshold statistics as a function of the size of the probing zones. This method makes it possible to derive physically grounded models of plasticity for amorphous materials by directly revealing the relevant details of the shear transformation zones that mediate this process.