Seismiclike organization of avalanches in a driven long-range elastic string as a paradigm of brittle cracks

Loss of memory of an elastic line on its way to limit cycles

Oscillatory-driven amorphous materials forget their initial configuration and converge to limit cycles. Here we investigate this memory loss under a nonquasistatic drive in a minimal model system, with quenched disorder and memory encoded in a spatial pattern, where oscillating protocols are formally replaced by a positive-velocity drive. We consider an elastic line driven athermally in a quenched disorder with biperiodic boundary conditions and tunable system size, thus controlling the area swept by the line per cycle as would the oscillation amplitude. The convergence to disorder-dependent limit cycle is strongly coupled to the nature of its velocity dynamics depending on system size. Based on the corresponding phase diagram, we propose a generic scenario for memory formation in disordered systems under finite driving rate.

Critical Scaling of Solid Fragmentation at Quasistatic and Finite Strain Rates

Using two-dimensional simulations of sheared, brittle solids, we characterize the resulting fragmentation and explore its underlying critical nature. Under quasistatic loading, a power-law distribution of fragment masses emerges after fracture which grows with increasing strain. With increasing strain rate, the maximum size of a grain decreases and a shallower distribution is produced. We propose a scaling theory for distributions based on a fractal scaling of the largest mass with system size in the quasistatic limit or with a correlation length that diverges as a power of rate in the finite-rate limit. Critical exponents are measured using finite-size scaling techniques.

Effect of thresholding on avalanches and their clustering for interfaces with long-range elasticity

Avalanches are often defined as signals higher than some detection level in bursty systems. The choice of the detection threshold affects the number of avalanches, but it can also affect their temporal correlations. We simulated the depinning of a long-range elastic interface and applied different thresholds including a zero one on the data to see how the sizes and durations of events change and how this affects temporal avalanche clustering. Higher thresholds result in steeper size and duration distributions and cause the avalanches to cluster temporally. Using methods from seismology, the frequency of the events in the clusters was found to decrease as a power-law of time, and the size of an event in a cluster was found to help predict how many events it is followed by. The results bring closer theoretical studies of this class of models to real experiments, but also highlight how different phenomena can be obtained from the same set of data.

Statistical properties of avalanches via the c-record process

We study the statistics of avalanches, as a response to an applied force, undergone by a particle hopping on a one-dimensional lattice where the pinning forces at each site are independent and identically distributed (i.i.d.), each drawn from a continuous f(x). The avalanches in this model correspond to the interrecord intervals in a modified record process of i.i.d. variables, defined by a single parameter c>0. This parameter characterizes the record formation via the recursive process R_{k}>R_{k-1}-c, where R_{k} denotes the value of the kth record. We show that for c>0, if f(x) decays slower than an exponential for large x, the record process is nonstationary as in the standard c=0 case. In contrast, if f(x) has a faster than exponential tail, the record process becomes stationary and the avalanche size distribution π(n) has a decay faster than 1/n^{2} for large n. The marginal case where f(x) decays exponentially for large x exhibits a phase transition from a nonstationary phase to a stationary phase as c increases through a critical value c_{crit}. Focusing on f(x)=e^{-x} (with x≥0), we show that c_{crit}=1 and for c<1, the record statistics is nonstationary. However, for c>1, the record statistics is stationary with avalanche size distribution π(n)∼n^{-1-λ(c)} for large n. Consequently, for c>1, the mean number of records up to N steps grows algebraically ∼N^{λ(c)} for large N. Remarkably, the exponent λ(c) depends continuously on c for c>1 and is given by the unique positive root of c=-ln(1-λ)/λ. We also unveil the presence of nontrivial correlations between avalanches in the stationary phase that resemble earthquake sequences.

Effects of external noise on threshold-induced correlations in ferromagnetic systems

In the present paper we investigate the impact of the external noise and detection threshold level on the simulation data for the systems that evolve through metastable states. As a representative model of such systems we chose the nonequilibrium athermal random-field Ising model with two types of the external noise, uniform white noise and Gaussian white noise with various different standard deviations, imposed on the original response signal obtained in model simulations. We applied a wide range of detection threshold levels in analysis of the signal and show how these quantities affect the values of exponent γ_{S/T} (describing the scaling of the average avalanche size with duration), the shift of waiting time between the avalanches, and finally the collapses of the waiting time distributions. The results are obtained via extensive numerical simulations on the equilateral three-dimensional cubic lattices of various sizes and disorders.

Controlling crackling dynamics by triggering low-intensity avalanches

We examine the effect of small, spatially localized excitations applied periodically in different manners, on the crackling dynamics of a brittle crack driven slowly in a heterogeneous solid. When properly adjusted, these excitations are observed to radically modify avalanche statistics and considerably limit the magnitude of the largest events. Surprisingly, this does not require information on the front loading state at the time of excitation; applying it either at a random location or at the most loaded point gives the same results. Subsequently, we unravel how the excitation amplitude, spatial extent, and frequency govern the effect. We find that the excitation efficiency is ruled by a single reduced parameter, namely the injected power per unit front length; the suppression of extreme avalanches is maximum at a well-defined optimal value of this control parameter. analysis opens another way to control the largest events in crackling dynamics. Beyond fracture problems, it may be relevant for crackling systems described by models of the same universality class, such as the wetting of heterogeneous substrates or magnetic walls in amorphous magnets.

Dynamic crack growth along heterogeneous planar interfaces: Interaction with unidimensional strips

We examine theoretically and numerically fast propagation of a tensile crack along unidimensional strips with periodically evolving toughness. In such dynamic fracture regimes, crack front waves form and transport front disturbances along the crack edge at speed less than the Rayleigh wave speed and depending on the crack speed. In this configuration, standing front waves dictate the spatiotemporal evolution of the local crack front speed, which takes a specific scaling form. Analytical examination of both the short-time and long-time limits of the problem reveals the parameter dependency with strip wavelength, toughness contrast and overall fracture speed. Implications and generalization to unidimensional strips of arbitrary shape are lastly discussed.

Understanding interfacial fracture behavior between microinterlocked soft layers using physics-based cohesive zone modeling

We examine the underlying fracture mechanics of the human skin dermal-epidermal layer's microinterlocks using a physics-based cohesive zone finite-element model. Using microfabrication techniques, we fabricated highly dense arrays of spherical microstructures of radius ≈50μm without and with undercuts, which occur in an open spherical cavity whose centroid lies below the microstructure surface to create microinterlocks in polydimethylsiloxane layers. From experimental peel tests, we find that the maximum density microinterlocks without and with undercuts enable the respective ≈4-fold and ≈5-fold increase in adhesion strength as compared to the plain layers. Critical visualization of the single microinterlock fracture from the cohesive zone model reveals a contact interaction-based phenomena where the primary propagating crack is arrested and the secondary crack is initiated in the microinterlocked area. Strain energy energetics confirmed significantly lower strain energy dissipation for the microinterlock with the undercut as compared to its nonundercut counterpart. These phenomena are completely absent in a plain interface fracture where the fracture propagates catastrophically without any arrests. These events confirm the difference in the experimental results corroborated by the Cook-Gordon mechanism. The findings from the cohesive zone simulation provide deeper insights into soft microinterlock fracture mechanics that could prominently help in the rational designing of sutureless skin grafts and electronic skin.

Universal Scaling of the Velocity Field in Crack Front Propagation

The propagation of a crack front in disordered materials is jerky and characterized by bursts of activity, called avalanches. These phenomena are the manifestation of an out-of-equilibrium phase transition originated by the disorder. As a result avalanches display universal scalings which are, however, difficult to characterize in experiments at a finite drive. Here, we show that the correlation functions of the velocity field along the front allow us to extract the critical exponents of the transition and to identify the universality class of the system. We employ these correlations to characterize the universal behavior of the transition in simulations and in an experiment of crack propagation. This analysis is robust, efficient, and can be extended to all systems displaying avalanche dynamics.