No Arabic abstract
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.
Athermal systems across a large range of length scales, ranging from foams and granular bead packings to crumpled metallic sheets, exhibit slow stress relaxation when compressed. Experimentally they show a non-monotonic stress response when decompressed somewhat after an initial compression, i.e. under a two-step, Kovacs-like protocol. It turns out that from this response one can tell the age of the system, suggesting an interpretation as a memory effect. In this work we use a model of an athermal jammed solid, specifically a binary mixture of soft harmonic spheres, to explore this phenomenon through in-silico experiments. Using extensive simulations under conditions analogous to those in experiment, we observe identical phenomenology in the stress response under a two--step protocol. Our model system also recovers the behaviour under a more recently studied three-step protocol, which consists of a compression followed by a decompression and then a final compression. We show that the observed response in both two-step and three-step protocols can be understood using Linear Response Theory. In particular, a linear scaling with age for the two-step protocol arises generically for slow linear responses with power law or logarithmic decay and does not in itself point to any underlying aging dynamics.
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$.
We analyze the fluctuations in particle positions and inter-particle forces in disordered jammed crystals in the limit of weak disorder. We demonstrate that such athermal systems are fundamentally different from their thermal counterparts, characterized by constrained fluctuations of forces perpendicular to the lattice directions. We develop a disorder perturbation expansion in polydispersity about the crystalline state, which we use to derive exact results to linear order. We show that constrained fluctuations result as a consequence of local force balance conditions, and are characterized by non-Gaussian distributions which we derive exactly. We analytically predict several properties of such systems, including the scaling of the average coordination with polydispersity and packing fraction, which we verify with numerical simulations using soft disks with one-sided harmonic interactions.
It is well known that jammed soft materials will flow if sheared above their yield stress - think mayonnaise spread on bread - but a complete microscopic description of this seemingly sim- ple process has yet to emerge. What remains elusive is a microscopic framework that explains the macroscopic flow, derived from a 3-D spatially resolved analysis of the dynamics of the droplets or particles that compose the soft material. By combining confocal-rheology experiments on compressed emulsions and numerical simulations, we unravel that the primary microscopic mechanisms for flow are strongly influenced by the rate of the imposed deformation. When shearing fast, small coordinated clusters of droplets move collectively as in a conga line, while at low rates the flow emerges from bursts of droplet rearrangements, correlated over large domains. These regions exhibit complex spatio-temporal correlation patterns that reflect the long range elasticity embedded in the jammed material. These results identify the three-dimensional structure of microscopic rearrangements within sheared soft solids, revealing that the characteristic shape and dynamics of these structures are strongly determined by the rate of the external shear.
Surface stress and surface energy are two fundamental parameters that determine the surface properties of any material. While it is commonly believed that the surface stress and surface energy of liquids are identical, the relationship between the two parameters in soft polymeric gels remains debatable. In this work, we measured the surface stress and surface energy of soft silicone gels with varying crosslinking densities in soft wetting experiments. Above a critical crosslink density, $k_0$, the surface stress is found to increase significantly with crosslinking density while the surface energy, by contrast, remains unchanged. In this regime, we can estimate a non-zero surface elastic modulus that also increases with the ratio of crosslinkers. By comparing the surface mechanics of the soft gels to their bulk rheology, the surface properties near the critical density $k_0$ are found to be closely related to the underlying percolation transition of the polymer networks.