We consider self-organization and memory formation of a sheared amorphous solid subject to a random shear strain protocol confined to a strain range $\pm \varepsilon_{\rm max}$. We show that, as in the case of applied cyclic shear, the response of the driven system retains a memory of the strain range that can subsequently be retrieved. Our findings suggest that self-organization and memory formation of disordered materials can emerge under more general conditions, such as a system interacting with its fluctuating environment.
The mechanical behavior of disordered materials such as dense suspensions, glasses or granular materials depends on their thermal and mechanical past. Here we report the memory behavior of a quenched mesoscopic elasto-plastic (QMEP) model. After prior oscillatory training, a simple read-out protocol gives access to both the training protocol's amplitude and the last shear direction. The memory of direction emerges from the development of a mechanical polarization during training. The analysis of sample-to-sample fluctuations gives direct access to the irreversibility transition. Despite the quadrupolar nature of the elastic interactions in amorphous solids, a behavior close to Return Point Memory (RPM) is observed. The quasi RPM property is used to build a simple Preisach-like model of directional memory.
The coarse-graining of amorphous plasticity from the atomistic to the mesoscopic scale is studied in the framework of a simple scalar elasto-plastic model. Building on recent results obtained on the atomistic scale, we discuss the interest in a disordered landscape-informed threshold disorder to reproduce the physics of amorphous plasticity. We show that accounting for a rejuvenation scenario allows us to reproduce quasi-quantitatively the evolution of the mean local yield stress and the localization behavior. We emphasize the crucial role of two dimensionless parameters: the relative strength of the yield stress disorder with respect to the typical stress drops associated with a plastic rearrangement, and the age parameter characterizing the relative stability of the initial glass with respect to the rejuvenated glass that emerges upon shear deformation.
While soft quasilocalized vibrational modes are known to populate the low-frequency spectrum of glassy solids, their contribution to thermal properties is still not fully elucidated. We numerically show that, despite their spatially localized nature, these modes are as effective heat carriers as the delocalized ones and can contribute non-negligible to the total thermal conductivity in the low-temperature regime, especially for T < 0.05 Tg where Tg is the glass transition temperature. We further prove that the mutual coherence between the low-frequency quasilocalized modes and other modes explains this high thermal exchange performance. Our finding finally provides a perspective on the thermal transport behaviour of the low-frequency quasilocalized modes in glassy solids.
We develop a mesoscopic model to study the plastic behavior of an amorphous material under cyclic loading. The model is depinning-like and driven by a disordered thresholds dynamics which are coupled by long-range elastic interactions. We propose a simple protocol of "glass preparation" which allows us to mimic thermalisation at high temperature, as well as aging at vanishing temperature. Various levels of glass stabilities (from brittle to ductile) can be achieved by tuning the aging duration. The aged glasses are then immersed into a quenched disorder landscape and serve as initial configurations for various protocols of mechanical loading by shearing. The dependence of the plastic behavior upon monotonous loading is recovered. The behavior under cyclic loading is studied for different ages and system sizes. The size and age dependence of the irreversibility transition is discussed. A thorough characterization of the disorder-landscape is achieved through the analysis of the transition graphs, which describe the plastic deformation pathways under athermal quasi-static shear. In particular, the analysis of the stability ranges of the strongly connected components of the transition graphs reveals the emergence of a phase-separation like process associated with the aging of the glass. Increasing the age and hence stability of the initial glass, results in a gradual break-up of the landscape of dynamically accessible stable states into three distinct regions: one region centered around the initially prepared glass phase, and two additional regions, characterized by well-separated ranges of positive and negative plastic strains, each of which is accessible only from the initial glass phase by passing through the stress peak in the forward, respectively, backward shearing directions.
While deeply supercooled liquids exhibit divergent viscosity and increasingly heterogeneous dynamics as the temperature drops, their structure shows only seemingly marginal changes. Understanding the nature of relaxation processes in this dramatic slowdown is key for understanding the glass transition. Here, we show by atomistic simulations that the heterogeneous dynamics of glass-forming liquids strongly correlate with the local residual plastic strengths along soft directions computed in the initial inherent structures. The correlation increases with decreasing temperature and is maximum in the vicinity of the relaxation time. For the lowest temperature investigated, this maximum is comparable with the best values from the literature dealing with the structure-property relationship. However, the nonlinear probe of the local shear resistance in soft directions provides here a real-space picture of relaxation processes. Our detection method of thermal rearrangements allows us to investigate the first passage time statistics and to study the scaling between the activation energy barriers and the residual plastic strengths. These results shed new light on the nature of relaxations of glassy systems by emphasizing the analogy between the thermal relaxations in viscous liquids and the plastic shear transformation in amorphous solids.
We compare the macroscopic and the local plastic behavior of a model amorphous solid based on two radically different numerical descriptions. On the one hand, we simulate glass samples by atomistic simulations. On the other, we implement a mesoscale elasto-plastic model based on a solid-mechanics description. The latter is extended to consider the anisotropy of the yield surface via statistically distributed local and discrete weak planes on which shear transformations can be activated. To make the comparison as quantitative as possible, we consider the simple case of a quasistatically driven two-dimensional system in the stationary flow state and compare mechanical observables measured on both models over the same length scales. We show that the macroscale response, including its fluctuations, can be quantitatively recovered for a range of elasto-plastic mesoscale parameters. Using a newly developed method that makes it possible to probe the local yield stresses in atomistic simulations, we calibrate the local mechanical response of the elasto-plastic model at different coarse-graining scales. In this case, the calibration shows a qualitative agreement only for an optimized subset of mesoscale parameters and for sufficiently coarse probing length scales. This calibration allows us to establish a length scale for the mesoscopic elements that corresponds to an upper bound of the shear transformation size, a key physical parameter in elasto-plastic models. We find that certain properties naturally emerge from the elasto-plastic model. In particular, we show that the elasto-plastic model reproduces the Bauschinger effect, namely the plasticity-induced anisotropy in the stress-strain response. We discuss the successes and failures of our approach, the impact of different model ingredients and propose future research directions for quantitative multi-scale models of amorphous plasticity.
D. Richard, M. Ozawa, S. Patinet, E. Stanifer, B. Shang, S.A. Ridout, B. Xu, G. Zhang, P.K. Morse, J.-L. Barrat, L. Berthier, M.L. Falk, P. Guan, A.J. Liu, K. Martens, S. Sastry, D. Vandembroucq, E. Lerner, M.L. Manning 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.
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 (AQS) 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_{\rm 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.
We study the structural origin of the Bauschinger effect by accessing numerically the local plastic thresholds in the steady state flow of a two-dimensional model glass under athermal quasistatic deformation. More specifically, we compute the local residual strength, $\Delta\tau^{c}$, for arbitrary loading orientations and find that plastic deformation generically induces material polarization, i.e., a forward-backward asymmetry in the $\Delta\tau^{c}$ distribution. In steady plastic flow, local packings are on average closer to forward (rather than backward) instabilities, due to the stress-induced bias of barriers. However, presumably due to mechanical noise, a significant fraction of zones lie close to reverse (backward) yielding, as the distribution of $\Delta\tau^{c}$ for reverse shearing extends quasilinearly down to zero local residual strength. By constructing an elementary model of the early plastic response, we then show that unloading causes reverse plasticity of a growing amplitude, i.e., reverse softening, while it shifts away forward-yielding barriers. This result in an inversion of polarization in the low-$\Delta\tau^{c}$ region and, consequently, in the Bauschinger effect. This scenario is quite generic, which explains the pervasiveness of the effect.
Structural heterogeneity of amorphous solids present difficult challenges that stymie the prediction of plastic events, which are intimately connected to their mechanical behavior. Based on a perturbation analysis of the potential energy landscape, we derive the atomic nonaffinity as an indicator with intrinsic orientation, which quantifies the contribution of an individual atom to the total nonaffine modulus of the system. We find that the atomic nonaffinity can efficiently characterize the locations of the shear transformation zones, with a predicative capacity comparable to the best indicators. More importantly, the atomic nonaffinity, combining the sign of third order derivative of energy with respect to coordinates, reveals an intrinsic softest shear orientation. By analyzing the angle between this orientation and the shear loading direction, it is possible to predict the protocol-dependent response of plastic events. Employing the new method, the distribution of orientations of shear transformation zones in a model two-dimensional amorphous solids can be measured. The resulting plastic events can be understood from a simple model of independent plastic events occurring at variously oriented shear transformation zones. These results shed light on the characterization and prediction of the mechanical response of amorphous solids.
We develop and extend a method presented in [S. Patinet, D. Vandembroucq, and M. L. Falk, Phys. Rev. Lett., 117, 045501 (2016)] 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 non-perturbative 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.
Granular chain packings exhibit a striking emergent strain-stiffening behavior despite the individual looseness of the constitutive chains. Using indentation experiments on such assemblies, we measure an exponential increase in the collective resistance force $F$ with the indentation depth $z$, and with the square root of the number $\mathcal{N}$ of beads per chain. These two observations are respectively reminiscent of the self-amplification of friction in a capstan or in interleaved books, as well as the physics of polymers. The experimental data are well captured by a novel model based on these two ingredients. Specifically, the resistance force is found to vary according to the universal relation: $\log F \sim \mu \sqrt{\mathcal{N}} \Phi^{11/8}z/ b $, where $\mu$ is the friction coefficient between two elementary beads, $b$ is their size, and $\Phi$ is the volume fraction of chain beads when semi-diluted in a surrounding medium of unconnected beads. Our study suggests that theories normally confined to the realm of polymer physics at a molecular level can be used to explain phenomena at a macroscopic level. This class of systems enables the study of friction in complex assemblies, with practical implications for the design of new materials, the textile industry, and biology.
Statistical models are essential to get a better understanding of the role of disorder in brittle disordered solids. Fiber bundle models play a special role as a paradigm, with a very good balance of simplicity and non-trivial effects. We introduce here a variant of the fiber bundle model where the load is transferred among the fibers through a very compliant membrane. This Soft Membrane fiber bundle mode reduces to the classical Local Load Sharing fiber bundle model in 1D. Highlighting the continuum limit of the model allows to compute an equivalent toughness for the fiber bundle and hence discuss nucleation of a critical defect. The computation of the toughness allows for drawing a simple connection with crack front propagation (depinning) models.
Understanding the role played by the microstructure of materials on their macroscopic failure properties is an important challenge in solid mechanics. Indeed, when a crack propagates at a heterogeneous brittle interface, the front is trapped by tougher regions and deforms. This pinning induces non-linearities in the crack propagation problem, even within Linear Elastic Fracture Mechanics theory, and modifies the overall failure properties of the material. For example crack front pinning by tougher places could increase the fracture resistance of multilayer structures, with interesting technological applications. Analytical perturbation approaches, based on Bueckner-Rice elastic line models, focus on the crack front perturbations, hence allow for a description of these phenomena. Here, they are applied to experiments investigating the propagation of a purely interfacial crack in a simple toughness pattern: a single defect strip surrounded by homogeneous interface. We show that by taking into account the finite size of the body, quantitative agreement with experimental and finite elements results is achieved. In particular this method allows to predict the toughness contrast, i.e. the toughness difference between the single defect strip and its homogeneous surrounding medium. This opens the way to a more accurate use of the perturbation method to study more disordered heterogeneous materials, where the finite elements method is less adequate. From our results, we also propose a simple method to determine the adhesion energy of tough interfaces by measuring the crack front deformation induced by known interface patterns.
The propagation of dislocations in random crystals is evidenced to be governed by atomic-scale avalanches whose the extension in space and the time intermittency characterizingly diverge at the critical threshold. Our work is the very first atomic-scale evidence that the paradigm of second order phase transitions applies to the depinning of elastic interfaces in random media.
A mesoscopic model of amorphous plasticity is discussed in the context of depinning models. After embedding in a d + 1 dimensional space, where the accumulated plastic strain lives along the additional dimension, the gradual plastic deformation of amorphous media can be regarded as the motion of an elastic manifold in a disordered landscape. While the associated depinning transition leads to scaling properties, the quadrupolar Eshelby interactions at play in amorphous plasticity induce specific additional features like shear-banding and weak ergodicity breakdown. The latters are shown to be controlled by the existence of soft modes of the elastic interaction, the consequence of which is discussed in the context of depinning.
The propagation of an adhesive crack through an anisotropic heterogeneous interface is considered. Tuning the local toughness distribution function and spatial correlation is numerically shown to induce a transition between weak to strong pinning conditions. While the macroscopic effective toughness is given by the mean local toughness in case of weak pinning, a systematic toughness enhancement is observed for strong pinning (the critical point of the depinning transition). A self-consistent approximation is shown to account very accurately for this evolution, without any free parameter.
On the basis of the classical dislocation theory, the Solid Solution Hardening (SSH) is commonly ascribed to the pinning of the edge dislocations. At the atomic level, the theoretical study of the dislocation cores contrasts with such a prediction. Using the static molecular simulations with some interatomic effective potentials, we demonstrate numerically that the critical resolved shear stress associated to a screw dislocation in a random Ni(Al) single crystal has same order as the edge one. Such a result is imposed by the details of the dislocation stacking fault and the core dissociation into Shockley partials. The SSH statistical theory is employed to tentatively predict analytically the data acquired through our atomistic simulations at different Al concentration.