Predicting sample lifetimes in creep fracture of heterogeneous materials

Ductile-to-brittle transition and yielding in soft amorphous materials: perspectives and open questions

Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.

Prediction of depinning transitions in interface models using Gini and Kolkata indices

The intermittent dynamics of driven interfaces through disordered media and its subsequent depinning for large enough driving force is a common feature for a myriad of diverse systems, starting from mode-I fracture, vortex lines in superconductors, and magnetic domain walls to invading fluid in a porous medium, to name a few. In this work, we outline a framework that can give a precursory signal of the imminent depinning transition by monitoring the variations in sizes or the inequality of the intermittent responses of a system that are seen prior to the depinning point. In particular, we use measures traditionally used to quantify economic inequality, i.e., the Gini index and the Kolkata index, for the case of the unequal responses of precritical systems. The crossing point of these two indices serves as a precursor to imminent depinning. Given a scale-free size distribution of the responses, we calculate the expressions for these indices, evaluate their crossing points, and give a recipe for forecasting depinning transitions. We apply this method to the Edwards-Wilkinson, Kardar-Parisi-Zhang, and fiber bundle model interface with variable interaction strengths and quenched disorder. The results are applicable for any interface dynamics undergoing a depinning transition. The results also explain previously observed near-universal values of Gini and Kolkata indices in self-organized critical systems.

Prediction of imminent failure using supervised learning in a fiber bundle model

Prediction of a breakdown in disordered solids under external loading is a question of paramount importance. Here we use a fiber bundle model for disordered solids and record the time series of the avalanche sizes and energy bursts. The time series contain statistical regularities that not only signify universality in the critical behavior of the process of fracture, but also reflect signals of proximity to a catastrophic failure. A systematic analysis of these series using supervised machine learning can predict the time to failure. Different features of the time series become important in different variants of training samples. We explain the reasons for such a switch over of importance among different features. We show that inequality measures for avalanche time series play a crucial role in imminent failure predictions, especially for imperfect training sets, i.e., when simulation parameters of training samples differ considerably from those of the testing samples. We also show the variation of predictability of the system as the interaction range and strengths of disorders are varied in the samples, varying the failure mode from brittle to quasibrittle (with interaction range) and from nucleation to percolation (with disorder strength). The effectiveness of the supervised learning is best when the samples just enter the quasibrittle mode of failure showing scale-free avalanche size distributions.

Prediction of creep failure time using machine learning

A subcritical load on a disordered material can induce creep damage. The creep rate in this case exhibits three temporal regimes viz. an initial decelerating regime followed by a steady-state regime and a stage of accelerating creep that ultimately leads to catastrophic breakdown. Due to the statistical regularities in the creep rate, the time evolution of creep rate has often been used to predict residual lifetime until catastrophic breakdown. However, in disordered samples, these efforts met with limited success. Nevertheless, it is clear that as the failure is approached, the damage become increasingly spatially correlated, and the spatio-temporal patterns of acoustic emission, which serve as a proxy for damage accumulation activity, are likely to mirror such correlations. However, due to the high dimensionality of the data and the complex nature of the correlations it is not straightforward to identify the said correlations and thereby the precursory signals of failure. Here we use supervised machine learning to estimate the remaining time to failure of samples of disordered materials. The machine learning algorithm uses as input the temporal signal provided by a mesoscale elastoplastic model for the evolution of creep damage in disordered solids. Machine learning algorithms are well-suited for assessing the proximity to failure from the time series of the acoustic emissions of sheared samples. We show that materials are relatively more predictable for higher disorder while are relatively less predictable for larger system sizes. We find that machine learning predictions, in the vast majority of cases, perform substantially better than other prediction approaches proposed in the literature.

Record statistics of bursts signals the onset of acceleration towards failure

Forecasting the imminent catastrophic failure has a high importance for a large variety of systems from the collapse of engineering constructions, through the emergence of landslides and earthquakes, to volcanic eruptions. Failure forecast methods predict the lifetime of the system based on the time-to-failure power law of observables describing the final acceleration towards failure. We show that the statistics of records of the event series of breaking bursts, accompanying the failure process, provides a powerful tool to detect the onset of acceleration, as an early warning of the impending catastrophe. We focus on the fracture of heterogeneous materials using a fiber bundle model, which exhibits transitions between perfectly brittle, quasi-brittle, and ductile behaviors as the amount of disorder is increased. Analyzing the lifetime of record size bursts, we demonstrate that the acceleration starts at a characteristic record rank, below which record breaking slows down due to the dominance of disorder in fracturing, while above it stress redistribution gives rise to an enhanced triggering of bursts and acceleration of the dynamics. The emergence of this signal depends on the degree of disorder making both highly brittle fracture of low disorder materials, and ductile fracture of strongly disordered ones, unpredictable.

System-size-dependent avalanche statistics in the limit of high disorder

We investigate the effect of the amount of disorder on the statistics of breaking bursts during the quasistatic fracture of heterogeneous materials. We consider a fiber bundle model where the strength of single fibers is sampled from a power-law distribution over a finite range, so that the amount of materials' disorder can be controlled by varying the power-law exponent and the upper cutoff of fibers' strength. Analytical calculations and computer simulations, performed in the limit of equal load sharing, revealed that depending on the disorder parameters the mechanical response of the bundle is either perfectly brittle where the first fiber breaking triggers a catastrophic avalanche, or it is quasibrittle where macroscopic failure is preceded by a sequence of bursts. In the quasibrittle phase, the statistics of avalanche sizes is found to show a high degree of complexity. In particular, we demonstrate that the functional form of the size distribution of bursts depends on the system size: for large upper cutoffs of fibers' strength, in small systems the sequence of bursts has a high degree of stationarity characterized by a power-law size distribution with a universal exponent. However, for sufficiently large bundles the breaking process accelerates towards the critical point of failure, which gives rise to a crossover between two power laws. The transition between the two regimes occurs at a characteristic system size which depends on the disorder parameters.

Machine learning plastic deformation of crystals

Plastic deformation of micron-scale crystalline solids exhibits stress-strain curves with significant sample-to-sample variations. It is a pertinent question if this variability is purely random or to some extent predictable. Here we show, by employing machine learning techniques such as regression neural networks and support vector machines that deformation predictability evolves with strain and crystal size. Using data from discrete dislocations dynamics simulations, the machine learning models are trained to infer the mapping from features of the pre-existing dislocation configuration to the stress-strain curves. The predictability vs strain relation is non-monotonic and exhibits a system size effect: larger systems are more predictable. Stochastic deformation avalanches give rise to fundamental limits of deformation predictability for intermediate strains. However, the large-strain deformation dynamics of the samples can be predicted surprisingly well.

Crack growth and energy dissipation in paper

Here, we follow the stable propagation of a roughening crack using simultaneously Digital Image Correlation and Infra-Red imaging. In a quasi-two-dimensional paper sample, the crack tip and ahead of that the fracture process zone follow the slowly, diffusively moving “hot spot” ahead of the tip. This also holds when the crack starts to roughen during propagation. The well-established intermittency of the crack advancement and the roughening of the crack in paper are thus subject to the dissipation and decohesion in the hot spot zone. They are therefore not only a result of the depinning of the crack in a heterogeneous material.

Avalanche Behavior in Creep Failure of Disordered Materials

We present a mesoscale elastoplastic model of creep in disordered materials, which considers temperature-dependent stochastic activation of localized deformation events that are coupled by internal stresses, leading to collective avalanche dynamics. We generalize this stochastic plasticity model by introducing damage in terms of a local strength that decreases, on statistical average, with increasing local plastic strain. The model captures failure in terms of strain localization in a catastrophic shear band concomitant with a finite-time singularity of the creep rate. The statistics of avalanches in the run-up to failure is characterized by a decreasing avalanche exponent τ that, at failure, approaches the value τ=1.5 typical of a critical branching process. The average avalanche rate exhibits an inverse Omori law as a function of time to failure. The distribution of interavalanche times turns out to be consistent with the epidemic-type aftershock sequences (ETAS) model of earthquake statistics.

Time-dependent fracture under unloading in a fiber bundle model

We investigate the fracture of heterogeneous materials occurring under unloading from an initial load. Based on a fiber bundle model of time-dependent fracture, we show that depending on the unloading rate the system has two phases: for rapid unloading the system suffers only partial failure and it has an infinite lifetime, while at slow unloading macroscopic failure occurs in a finite time. The transition between the two phases proved to be analogous to continuous phase transitions. Computer simulations revealed that during unloading the fracture proceeds in bursts of local breakings triggered by slowly accumulating damage. In both phases the time evolution starts with a relaxation of the bursting activity characterized by a universal power-law decay of the burst rate. In the phase of finite lifetime the initial slowdown is followed by an acceleration towards macroscopic failure where the increasing rate of bursts obeys the (inverse) Omori law of earthquakes. We pointed out a strong correlation between the time where the event rate reaches a minimum value and of the lifetime of the system which allows for forecasting of the imminent catastrophic failure.

Microscopic dynamics and failure precursors of a gel under mechanical load

Material failure is ubiquitous, with implications from geology to everyday life and material science. It often involves sudden, unpredictable events, with little or no macroscopically detectable precursors. A deeper understanding of the microscopic mechanisms eventually leading to failure is clearly required, but experiments remain scarce. Here, we show that the microscopic dynamics of a colloidal gel, a model network-forming system, exhibit dramatic changes that precede its macroscopic failure by thousands of seconds. Using an original setup coupling light scattering and rheology, we simultaneously measure the macroscopic deformation and the microscopic dynamics of the gel, while applying a constant shear stress. We show that the network failure is preceded by qualitative and quantitative changes of the dynamics, from reversible particle displacements to a burst of irreversible plastic rearrangements.

Catastrophic Failure and Critical Scaling Laws of Fiber Bundle Material

This paper presents a spring-fiber bundle model used to describe the failure process induced by energy release in heterogeneous materials. The conditions that induce catastrophic failure are determined by geometric conditions and energy equilibrium. It is revealed that the relative rates of deformation of, and damage to the fiber bundle with respect to the boundary controlling displacement ε₀ exhibit universal power law behavior near the catastrophic point, with a critical exponent of −1/2. The proportion of the rate of response with respect to acceleration exhibits a linear relationship with increasing displacement in the vicinity of the catastrophic point. This allows for the prediction of catastrophic failure immediately prior to failure by extrapolating the trajectory of this relationship as it asymptotes to zero. Monte Carlo simulations are completed and these two critical scaling laws are confirmed.

Time-dependent breakdown of fiber networks: Uncertainty of lifetime

Materials often fail when subjected to stresses over a prolonged period. The time to failure, also called the lifetime, is known to exhibit large variability of many materials, particularly brittle and quasibrittle materials. For example, a coefficient of variation reaches 100% or even more. Its distribution shape is highly skewed toward zero lifetime, implying a large number of premature failures. This behavior contrasts with that of normal strength, which shows a variation of only 4%-10% and a nearly bell-shaped distribution. The fundamental cause of this large and unique variability of lifetime is not well understood because of the complex interplay between stochastic processes taking place on the molecular level and the hierarchical and disordered structure of the material. We have constructed fiber network models, both regular and random, as a paradigm for general material structures. With such networks, we have performed Monte Carlo simulations of creep failure to establish explicit relationships among fiber characteristics, network structures, system size, and lifetime distribution. We found that fiber characteristics have large, sometimes dominating, influences on the lifetime variability of a network. Among the factors investigated, geometrical disorders of the network were found to be essential to explain the large variability and highly skewed shape of the lifetime distribution. With increasing network size, the distribution asymptotically approaches a double-exponential form. The implication of this result is that, so-called "infant mortality," which is often predicted by the Weibull approximation of the lifetime distribution, may not exist for a large system.

Predicting and assessing rupture in protein gels under oscillatory shear

Soft materials may break irreversibly upon applying sufficiently large shear oscillations, a process whose physical mechanism remains largely elusive. In this work, the rupture of protein gels made of sodium caseinate under an oscillatory stress is shown to occur in an abrupt, brittle-like manner. Upon increasing the stress amplitude, the build-up of harmonic modes in the strain response can be rescaled for all gel concentrations. This rescaling yields an empirical criterion to predict the rupture point way before the samples are significantly damaged. "Fatigue" experiments under stress oscillations of constant amplitude can be mapped onto the former results, which indicates that rupture is independent of the temporal pathway in which strain and damage accumulate. Finally, using ultrasonic imaging, we measure the local mechanical properties of the gels before, during and after breakdown, showing that the strain field remains perfectly homogeneous up to rupture but suddenly gives way to a solid-fluid phase separation upon breakdown.