No Arabic abstract
Rate-effects in sheared disordered solids are studied using molecular dynamics simulations of binary Lennard-Jones glasses in two and three dimensions. In the quasistatic (QS) regime, systems exhibit critical behavior: the magnitudes of avalanches are power-law distributed with a maximum cutoff that diverges with increasing system size $L$. With increasing rate, systems move away from the critical yielding point and the average flow stress rises as a power of the strain rate with exponent $1/beta$, the Herschel-Bulkley exponent. Finite-size scaling collapses of the stress are used to measure $beta$ as well as the exponent $ u$ which characterizes the divergence of the correlation length. The stress and kinetic energy per particle experience fluctuations with strain that scale as $L^{-d/2}$. As the largest avalanche in a system scales as $L^alpha$, this implies $alpha < d/2$. The diffusion rate of particles diverges as a power of decreasing rate before saturating in the QS regime. A scaling theory for the diffusion is derived using the QS avalanche rate distribution and generalized to the finite strain rate regime. This theory is used to collapse curves for different system sizes and confirm $beta/ u$.
Disordered solids respond to quasistatic shear with intermittent avalanches of plastic activity, an example of the crackling noise observed in many nonequilibrium critical systems. The temporal power spectrum of activity within disordered solids consists of three distinct domains: a novel power-law rise with frequency at low frequencies indicating anticorrelation, white-noise at intermediate frequencies, and a power-law decay at high frequencies. As the strain rate increases, the white-noise regime shrinks and ultimately disappears as the finite strain rate restricts the maximum size of an avalanche. A new strain-rate- and system-size-dependent scaling theory is derived for power spectra in both the quasistatic and finite-strain-rate regimes. This theory is validated using data from overdamped two- and three-dimensional molecular dynamics simulations. We identify important exponents in the yielding transition including the dynamic exponent $z$ which relates the size of an avalanche to its duration, the fractal dimension of avalanches, and the exponent characterizing the divergence in correlations with strain rate. Results are related to temporal correlations within a single avalanche and between multiple avalanches.
In soft amorphous materials, shear cessation after large shear deformation leads to structures having residual shear stress. The origin of these states and the distribution of the local shear stresses within the material is not well understood, despite its importance for the change in material properties and consequent applications. In this work, we use molecular dynamics simulations of a model dense non-Brownian soft amorphous material to probe the non-trivial relaxation process towards a residual stress state. We find that, similar to thermal glasses, an increase in shear rate prior to the shear cessation leads to lower residual stress states. We rationalise our findings using a mesoscopic elasto-plastic description that explicitly includes a long range elastic response to local shear transformations. We find that after flow cessation the initial stress relaxation indeed depends on the pre-sheared stress state, but the final residual stress is majorly determined by newly activated plastic events occurring during the relaxation process. Our simplified coarse grained description not only allows to capture the phenomenology of residual stress states but also to rationalise the altered material properties that are probed using small and large deformation protocols applied to the relaxed material.
Granular packings of non-convex or elongated particles can form free-standing structures like walls or arches. For some particle shapes, such as staples, the rigidity arises from interlocking of pairs of particles, but the origins of rigidity for non-interlocking particles remains unclear. We report on experiments and numerical simulations of sheared columns of hexapods, particles consisting of three mutually orthogonal sphero-cylinders whose centers coincide. We vary the length-to-diameter aspect ratio, $alpha$, of the sphero-cylinders and subject the packings to quasistatic direct shear. For small $alpha$, we observe a finite yield stress. For large $alpha$, however, the column becomes rigid when sheared, supporting stresses that increase sharply with increasing strain. Analysis of X-ray micro-computed tomography (Micro-CT) data collected during the shear reveals that the stiffening is associated with a tilted, oblate cluster of hexapods near the nominal shear plane in which particle deformation and average contact number both increase. Simulation results show that the particles are collectively under tension along one direction even though they do not interlock pairwise. These tensions comes from contact forces carrying large torques, and they are perpendicular to the compressive stresses in the packing. They counteract the tendency to dilate, thus stabilize the particle cluster.
When an amorphous solid is deformed cyclically, it may reach a steady state in which the paths of constituent particles trace out closed loops that repeat in each driving cycle. A remarkable variant has been noticed in simulations where the period of particle motions is a multiple of the period of driving, but the reasons for this behavior have remained unclear. Motivated by mesoscopic features of displacement fields in experiments on jammed solids, we propose and analyze a simple model of interacting soft spots -- locations where particles rearrange under stress and that resemble two-level systems with hysteresis. We show that multiperiodic behavior can arise among just three or more soft spots that interact with each other, but in all cases it requires frustrated interactions, illuminating this otherwise elusive type of interaction. We suggest directions for seeking this signature of frustration in experiments and for achieving it in designed systems.
Amorphous solids lack long-range order. Therefore identifying structural defects -- akin to dislocations in crystalline solids -- that carry plastic flow in these systems remains a daunting challenge. By comparing many different structural indicators in computational models of glasses, under a variety of conditions we carefully assess which of these indicators are able to robustly identify the structural defects responsible for plastic flow in amorphous solids. We further demonstrate that the density of defects changes as a function of material preparation and strain in a manner that is highly correlated with the macroscopic material response. Our work represents an important step towards predicting how and when an amorphous solid will fail from its microscopic structure.