Crack propagation through disordered materials as a depinning transition: A critical test of the theory

Percolation versus depinning transition: The inherent role of damage hardening during quasibrittle failure

The intermittent damage evolution preceding the failure of heterogeneous brittle solids is well described by scaling laws. In deciphering its origins, failure is routinely interpreted as a critical transition. However at odds with expectations of universality, a large scatter in the value of the scaling exponents is reported during acoustic emission experiments. Here we numerically examine the precursory damage activity to reconcile the experimental observations with critical phenomena framework. Along with the strength of disorder, we consider an additional parameter that describes the progressive damageability of material elements at mesoscopic scale. This hardening behavior encapsulates the microfracturing processes taking place at lower length scales. We find that damage hardening can not only delay the final failure but also affect the preceding damage accumulation. When hardening is low, the precursory activity is strongly influenced by the strength of the disorder and is reminiscent of damage percolation. On the contrary, for large hardening, long-range elastic interactions prevail over disorder, ensuring a rather homogeneous evolution of the damage field in the material. The power-law statistics of the damage bursts is robust to the strength of the disorder and is reminiscent of the collective avalanche dynamics of elastic interfaces near the depinning transition. The existence of these two distinct universality classes also manifests as different values of the scaling exponent characterizing the divergence of the precursor size on approaching failure. Our finding sheds new light on the connection between the level of quasibrittleness of materials and the statistical features of the failure precursors. Finally, it also provides a more complete description of the acoustic precursors and thus paves the way for quantitative techniques of damage monitoring of structures-in-service.

Edwards-Wilkinson depinning transition in fractional Brownian motion background

There are various reports about the critical exponents associated with the depinning transition. In this study, we investigate how the disorder strength present in the support can account for this diversity. Specifically, we examine the depinning transition in the quenched Edwards-Wilkinson (QEW) model on a correlated square lattice, where the correlations are modeled using fractional Brownian motion (FBM) with a Hurst exponent of H.We identify a crossover time [Formula: see text] that separates the dynamics into two distinct regimes: for [Formula: see text], we observe the typical behavior of pinned surfaces, while for [Formula: see text], the behavior differs. We introduce a novel three-variable scaling function that governs the depinning transition for all considered H values. The associated critical exponents exhibit a continuous variation with H, displaying distinct behaviors for anti-correlated ([Formula: see text]) and correlated ([Formula: see text]) cases. The critical driving force decreases with increasing H, as the host medium becomes smoother for higher H values, facilitating fluid mobility. This fact causes the asymptotic velocity exponent [Formula: see text] to increase monotonically with H.

Thermal effects on fracture and the brittle-to-ductile transition

The fracture behavior of brittle and ductile materials can be strongly influenced by thermal fluctuations, especially in micro- and nanodevices as well as in rubberlike and biological materials. However, temperature effects, in particular on the brittle-to-ductile transition, still require a deeper theoretical investigation. As a step in this direction we propose a theory, based on equilibrium statistical mechanics, able to describe the temperature-dependent brittle fracture and brittle-to-ductile transition in prototypical discrete systems consisting in a lattice with breakable elements. Concerning the brittle behavior, we obtain closed form expressions for the temperature-dependent fracture stress and strain, representing a generalized Griffith criterion, ultimately describing the fracture as a genuine phase transition. With regard to the brittle-to-ductile transition, we obtain a complex critical scenario characterized by a threshold temperature between the two fracture regimes (brittle and ductile), an upper and a lower yield strength, and a critical temperature corresponding to the complete breakdown. To show the effectiveness of the proposed models in describing thermal fracture behaviors at small scales, we successfully compare our theoretical results with molecular dynamics simulations of Si and GaN nanowires.

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.

Roughening of the anharmonic elastic interface in correlated random media

We study the roughening properties of the anharmonic elastic interface in the presence of temporally correlated noise. The model can be seen as a generalization of the anharmonic Larkin model, recently introduced by Purrello, Iguain, and Kolton [Phys. Rev. E 99, 032105 (2019)2470-004510.1103/PhysRevE.99.032105], to investigate the effect of higher-order corrections to linear elasticity in the fate of interfaces. We find analytical expressions for the critical exponents as a function of the anharmonicity index n, the noise correlator range θ∈[0,1/2], and dimension d. In d=1 we find that the interface becomes faceted and exhibits anomalous scaling for θ>1/4 for any degree of anharmonicity n>1. Analytical expressions for the anomalous exponents α_{loc} and κ are obtained and compared with a numerical integration of the model. Our theoretical results show that anomalous roughening cannot exist for this model in dimensions d>1.

Effective Toughness of Heterogeneous Materials with Rate-Dependent Fracture Energy

We investigate the dynamic fracture of heterogeneous materials experimentally by measuring displacement fields as a rupture propagates through a periodic array of obstacles of controlled fracture energy. Our measurements demonstrate the applicability of the classical equation of motion of cracks at a discontinuity of fracture energy: the crack speed jumps at the entrance and exit of an obstacle, as predicted by the crack-tip energy balance within the brittle fracture framework. The speed jump amplitude is governed by the fracture energy contrast and by the combination of the rate dependency of the fracture energy and the inertia of the medium, which allows the crack to cross a fracture energy discontinuity at a constant energy release rate. This discontinuous dynamics and the rate dependence cause higher effective toughness, which governs the coarse-grained behavior of these cracks.

Non-Gaussian tail in the force distribution: a hallmark of correlated disorder in the host media of elastic objects

Inferring the nature of disorder in the media where elastic objects are nucleated is of crucial importance for many applications but remains a challenging basic-science problem. Here we propose a method to discern whether weak-point or strong-correlated disorder dominates based on characterizing the distribution of the interaction forces between objects mapped in large fields-of-view. We illustrate our proposal with the case-study system of vortex structures nucleated in type-II superconductors with different pinning landscapes. Interaction force distributions are computed from individual vortex positions imaged in thousands-vortices fields-of-view in a two-orders-of-magnitude-wide vortex-density range. Vortex structures nucleated in point-disordered media present Gaussian distributions of the interaction force components. In contrast, if the media have dilute and randomly-distributed correlated disorder, these distributions present non-Gaussian algebraically-decaying tails for large force magnitudes. We propose that detecting this deviation from the Gaussian behavior is a fingerprint of strong disorder, in our case originated from a dilute distribution of correlated pinning centers.

Depinning Dynamics of Crack Fronts

We investigate experimentally and theoretically the dynamics of a crack front during the microinstabilities taking place in heterogeneous materials between two successive equilibrium positions. We focus specifically on the spatiotemporal evolution of the front, as it relaxes to a straight configuration, after depinning from a single obstacle of controlled strength and size. We show that this depinning dynamics is not controlled by inertia, but instead by the rate dependency of the dissipative mechanisms taking place within the fracture process zone. This implies that the crack speed fluctuations around its average value v_{m} can be predicted from an overdamped equation of motion (v-v_{m})/v_{0}=[G-G_{c}(v_{m})]/G_{c}(v_{m}) involving the characteristic material speed v_{0}=G_{c}(v_{m})/G_{c}^{'}(v_{m}) that emerges from the variation of fracture energy with crack speed. Our findings pave the way to a quantitative description of the critical depinning dynamics of cracks in disordered solids and open up new perspectives for the prediction of the effective failure properties of heterogeneous materials.

Crack growth in heterogeneous brittle solids: intermittency, crackling and induced seismicity

Crack growth is the basic mechanism leading to the failure of brittle materials. Engineering addresses this problem within the framework of continuum mechanics, which links deterministically the crack motion to the applied loading. Such an idealization, however, fails in several situations and in particular cannot capture the highly erratic (earthquake-like) dynamics sometimes observed in slowly fracturing heterogeneous solids. Here, we examine this problem by means of innovative experiments of crack growth in artificial rocks of controlled microstructure. The dynamical events are analysed at both global and local scales, from the time fluctuation of the spatially averaged crack speed and the induced acoustic emission, respectively. Their statistics are characterized and compared with the predictions of a recent approach mapping fracture onset to the depinning of an elastic interface. Finally, the overall time-size organization of the events is characterized to shed light on the mechanisms underlying the scaling laws observed in seismology.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Diffusion in Mesoscopic Lattice Models of Amorphous Plasticity

We present results on tagged particle diffusion in a mesoscale lattice model for sheared amorphous material in athermal quasistatic conditions. We find a short time diffusive regime and a long time diffusive regime whose diffusion coefficients depend on system size in dramatically different ways. At short time, we find that the diffusion coefficient, D, scales roughly linearly with system length, D∼L^{1.05}. This short time behavior is consistent with particle-based simulations. The long-time diffusion coefficient scales like D∼L^{1.6}, close to previous studies which found D∼L^{1.5}. Furthermore, we show that the near-field details of the interaction kernel do not affect the short time behavior but qualitatively and dramatically affect the long time behavior, potentially causing a saturation of the mean-squared displacement at long times. Our finding of a D∼L^{1.05} short time scaling resolves a long standing puzzle about the disagreement between the diffusion coefficient measured in particle-based models and mesoscale lattice models of amorphous plasticity.