Crossover behavior in burst avalanches: signature of imminent failure

Critical Scaling through Gini Index

In the systems showing critical behavior, various response functions have a singularity at the critical point. Therefore, as the driving field is tuned toward its critical value, the response functions change drastically, typically diverging with universal critical exponents. In this Letter, we quantify the inequality of response functions with measures traditionally used in economics, namely by constructing a Lorenz curve and calculating the corresponding Gini index. The scaling of such a response function, when written in terms of the Gini index, shows singularity at a point that is at least as universal as the corresponding critical exponent. The critical scaling, therefore, becomes a single parameter fit, which is a considerable simplification from the usual form where the critical point and critical exponents are independent. We also show that another measure of inequality, the Kolkata index, crosses the Gini index at a point just prior to the critical point. Therefore, monitoring these two inequality indices for a system where the critical point is not known can produce a precursory signal for the imminent criticality. This could be useful in many systems, including that in condensed matter, bio- and geophysics to atmospheric physics. The generality and numerical validity of the calculations are shown with the Monte Carlo simulations of the two dimensional Ising model, site percolation on square lattice, and the fiber bundle model of fracture.

Inequality of avalanche sizes in models of fracture

Prediction of an imminent catastrophic event in a driven disordered system is of paramount importance-from the laboratory scale controlled fracture experiment to the largest scale of mechanical failure, i.e., earthquakes. It has long been conjectured that the statistical regularities in the energy emission time series mirror the "health" of such driven systems and hence have the potential for forecasting imminent catastrophe. Among other statistical regularities, a measure of how unequal avalanche sizes are is potentially a crucial indicator of imminent failure. The inequalities of avalanche sizes are quantified using inequality indices traditionally used in socioeconomic systems: the Gini index g, the Hirsch index h, and the Kolkata index k. It is shown analytically (for the mean-field case) and numerically (for the non-mean-field case) with models of quasi-brittle materials that the indices show universal behavior near the breaking points in such models and hence could serve as indicators of imminent breakdown of stressed disordered systems.

Analysis of Acoustic Emission Energy Distribution and Avalanche Dynamics of Sandstone with Different Particle Sizes

The mechanical properties of sandstone have an important impact on the stability of the coal mine roof and floor, sandstone gas mining, and underground engineering safety. In order to study the critical characteristics on the failure process of sandstone with different particle sizes under uniaxial compression conditions, avalanche dynamics theory and a critical model are used to analyze the distribution of acoustic emission (AE) parameters, and the maximum likelihood estimation is used to accurately estimate the critical parameters. The results showed that the AE phenomenon of sandstone can be divided into four stages: initial compaction period, quiet period, crack stable growth period and outbreak period. During the process of compression failure, the larger the particle size is, the more seriously the sandstone is damaged. The AE energy probability density distribution follows single power-law distribution, and the AE energy critical exponent is 1.20 and follows the characteristics of scale-free regarding the power-law distribution on the particle sizes. When the stress runs up to 90% of peak stress, the bifurcation ratio increases sharply and shows the characteristics of the critical state. The waiting time and the avalanche size distribution follow double power-law distribution, and the inflection points are 0.03 and 37. Before and after the inflection point, the waiting time critical exponent and the avalanche size critical exponent are 1.90, 0.40 and 2.40, 1.60. This shows that the dynamic evolution process of sandstone under uniaxial compression condition can be characterized well by the fiber bundle model.

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.

Energy exponents of avalanches and Hausdorff dimensions of collapse patterns

A simple numerical model to simulate athermal avalanches is presented. The model is inspired by the "porous collapse" process where the compression of porous materials generates collapse cascades, leading to power law distributed avalanches. The energy (E), amplitude (A_{max}), and size (S) exponents are derived by computer simulation in two approximations. Time-dependent "jerk" spectra are calculated in a single avalanche model where each avalanche is simulated separately from other avalanches. The average avalanche profile is parabolic, the scaling between energy and amplitude follows E∼A_{max}^{2}, and the energy exponent is ε = 1.33. Adding a general noise term in a continuous event model generates infinite avalanche sequences which allow the evaluation of waiting time distributions and pattern formation. We find the validity of the Omori law and the same exponents as in the single avalanche model. We then add spatial correlations by stipulating the ratio G/N between growth processes G (linked to a previous event location) and nucleation processes N (with new, randomly chosen nucleation sites). We found, in good approximation, a power law correlation between the energy exponent ε and the Hausdorff dimension H_{D} of the resulting collapse pattern H_{D}-1∼ɛ^{-3}. The evolving patterns depend strongly on G/N with the distribution of collapse sites equally power law distributed. Its exponent ɛ_{topo} would be linked to the dynamical exponent ε if each collapse carried an energy equivalent to the size of the collapse. A complex correlation between ɛ,ɛ_{topo}, and H_{D} emerges, depending strongly on the relative occupancy of the collapse sites in the simulation box.

Role of spatial patterns in fracture of disordered multiphase materials

Multiphase materials, such as composite materials, exhibit multiple competing failure mechanisms during the growth of a macroscopic defect. For the simulation of the overall fracture process in such materials, we develop a two-phase spring network model that accounts for the architecture between the different components as well as the respective disorders in their failure characteristics. In the specific case of a plain weave architecture, we show that any offset between the layers reduces the delocalization of the stresses at the crack tip and thereby substantially lowers the strength and fracture toughness of the overall laminate. The avalanche statistics of the broken springs do not show a distinguishable dependence on the offsets between layers. The power-law exponents are found to be much smaller than that of disordered spring network models in the absence of a crack. A discussion is developed on the possibility of the avalanche statistics being those near breakdown.

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.

Brittle to quasibrittle transition in a compound fiber bundle

The brittle to quasibrittle transition has been studied for a compound of two different kinds of fibrous materials, having distinct difference in their breaking strengths under the framework of the fiber bundle model. A random fiber bundle model has been devised with a bimodal distribution of the breaking strengths of the individual fibers. The bimodal distribution is assumed to consist of two symmetrically placed rectangular probability distributions of strengths p and 1-p, each of width d, and separated by a gap 2s. Different properties of the transition have been studied varying these three parameters and using the well-known equal load-sharing dynamics. Our study exhibits a brittle to quasibrittle transition at the critical width d_{c}(s,p)=p(1/2-s)/(1+p) confirmed by our numerical results.

Fiber-bundle model with time-dependent healing mechanisms to simulate progressive failure of snow

Snow is a heterogeneous material with strain- and/or load-rate-dependent strength. In particular, a transition from ductile-to-brittle failure behavior with increasing load rate is observed. The rate-dependent behavior can partly be explained with the existence of a unique healing mechanism in snow that stems from its high homologous temperature (temperature close to melting point). As soon as broken elements in the ice matrix get in contact, they start sintering and the structure may regain strength. Moreover, the ice matrix is subjected to viscous deformation, inducing a relaxation of local load concentrations and, therefore, further counteracting the damage process. Ideal tools for studying the failure process of heterogeneous materials are the fiber-bundle models (FBMs), which allow investigating the effects of basic microstructural characteristics on the general macroscopic failure behavior. We present an FBM with two concurrent time-dependent healing mechanisms: sintering of broken fibers and relaxation of load inhomogeneities. Sintering compensates damage by creating additional intact, load-supporting fibers which lead to an increase of the bundle strength. However, the character of the failure is not changed by sintering alone. With combined sintering and load relaxation, load is distributed from old stronger fibers to new fibers that carry fewer load. So as we additionally incorporated load redistribution to the FBM, the failure occurred suddenly without decrease of the order parameter-describing the amount of damage in the bundle-and without divergence of the fiber failure rate. Moreover, the b value, i.e., the power-law exponent of frequency-magnitude statistics of fibers breaking in load redistribution steps, at failure converged to b≈2, a value higher than that of a classical FBM without healing (b=3/2). These results indicate that healing, as the combined effect of sintering and load relaxation, changes the type of the phase transition at failure. This change of the phase transition is important for quantifying or predicting the failure (e.g., by monitoring acoustic emissions) of snow or other materials for which healing plays an important role.

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.

Precursory changes in seismic velocity for the spectrum of earthquake failure modes

Temporal changes in seismic velocity during the earthquake cycle have the potential to illuminate physical processes associated with fault weakening and connections between the range of fault slip behaviors including slow earthquakes, tremor and low frequency earthquakes1. Laboratory and theoretical studies predict changes in seismic velocity prior to earthquake failure2, however tectonic faults fail in a spectrum of modes and little is known about precursors for those modes3. Here we show that precursory changes of wave speed occur in laboratory faults for the complete spectrum of failure modes observed for tectonic faults. We systematically altered the stiffness of the loading system to reproduce the transition from slow to fast stick-slip and monitored ultrasonic wave speed during frictional sliding. We find systematic variations of elastic properties during the seismic cycle for both slow and fast earthquakes indicating similar physical mechanisms during rupture nucleation. Our data show that accelerated fault creep causes reduction of seismic velocity and elastic moduli during the preparatory phase preceding failure, which suggests that real time monitoring of active faults may be a means to detect earthquake precursors.

Blending stiffness and strength disorder can stabilize fracture

Quasibrittle behavior, where macroscopic failure is preceded by stable damaging and intensive cracking activity, is a desired feature of materials because it makes fracture predictable. Based on a fiber-bundle model with global load sharing we show that blending strength and stiffness disorder of material elements leads to the stabilization of fracture, i.e., samples that are brittle when one source of disorder is present become quasibrittle as a consequence of blending. We derive a condition of quasibrittle behavior in terms of the joint distribution of the two sources of disorder. Breaking bursts have a power-law size distribution of exponent 5/2 without any crossover to a lower exponent when the amount of disorder is gradually decreased. The results have practical relevance for the design of materials to increase the safety of constructions.

Failure criterion for materials with spatially correlated mechanical properties

The role of spatially correlated mechanical elements in the failure behavior of heterogeneous materials represented by fiber bundle models (FBMs) was evaluated systematically for different load redistribution rules. Increasing the range of spatial correlation for FBMs with local load sharing is marked by a transition from ductilelike failure characteristics into brittlelike failure. The study identified a global failure criterion based on macroscopic properties (external load and cumulative damage) that is independent of spatial correlation or load redistribution rules. This general metric could be applied to assess the mechanical stability of complex and heterogeneous systems and thus provide an important component for early warning of a class of geophysical ruptures.

Fiber bundle model with highly disordered breaking thresholds

We present a study of the fiber bundle model using equal load-sharing dynamics where the breaking thresholds of the fibers are drawn randomly from a power-law distribution of the form p(b)∼b-1 in the range 10-β to 10β. Tuning the value of β continuously over a wide range, the critical behavior of the fiber bundle has been studied both analytically as well as numerically. Our results are: (i) The critical load σc(β,N) for the bundle of size N approaches its asymptotic value σc(β) as σc(β,N)=σc(β)+AN-1/ν(β), where σc(β) has been obtained analytically as σc(β)=10β/(2βeln10) for β≥βu=1/(2ln10), and for β<βu the weakest fiber failure leads to the catastrophic breakdown of the entire fiber bundle, similar to brittle materials, leading to σ_{c}(β)=10-β; (ii) the fraction of broken fibers right before the complete breakdown of the bundle has the form 1-1/(2βln10); (iii) the distribution D(Δ) of the avalanches of size Δ follows a power-law D(Δ)∼Δ-ξ with ξ=5/2 for Δ≫Δc(β) and ξ=3/2 for Δ≪Δc(β), where the crossover avalanche size Δc(β)=2/(1-e10-2β)2.

Approach to failure in porous granular materials under compression

We investigate the approach to catastrophic failure in a model porous granular material undergoing uniaxial compression. A discrete element computational model is used to simulate both the microstructure of the material and the complex dynamics and feedbacks involved in local fracturing and the production of crackling noise. Under strain-controlled loading, microcracks initially nucleate in an uncorrelated way all over the sample. As loading proceeds the damage localizes into a narrow damage band inclined at 30°-45° to the load direction. Inside the damage band the material is crushed into a poorly sorted mixture of mainly fine powder hosting some larger fragments. The mass probability density distribution of particles in the damage zone is a power law of exponent 2.1, similar to a value of 1.87 inferred from observations of the length distribution of wear products (gouge) in natural and laboratory faults. Dynamic bursts of radiated energy, analogous to acoustic emissions observed in laboratory experiments on porous sedimentary rocks, are identified as correlated trails or cascades of local ruptures that emerge from the stress redistribution process. As the system approaches macroscopic failure consecutive bursts become progressively more correlated. Their size distribution is also a power law, with an equivalent Gutenberg-Richter b value of 1.22 averaged over the whole test, ranging from 3 to 0.5 at the time of failure, all similar to those observed in laboratory tests on granular sandstone samples. The formation of the damage band itself is marked by a decrease in the average distance between consecutive bursts and an emergent power-law correlation integral of event locations with a correlation dimension of 2.55, also similar to those observed in the laboratory (between 2.75 and 2.25).

Self-organized dynamics in local load-sharing fiber bundle models

We study the dynamics of a local load-sharing fiber bundle model in two dimensions under an external load (which increases with time at a fixed slow rate) applied at a single point. Due to the local load-sharing nature, the redistributed load remains localized along the boundary of the broken patch. The system then goes to a self-organized state with a stationary average value of load per fiber along the (increasing) boundary of the broken patch (damaged region) and a scale-free distribution of avalanche sizes and other related quantities are observed. In particular, when the load redistribution is only among nearest surviving fiber(s), the numerical estimates of the exponent values are comparable with those of the Manna model. When the load redistribution is uniform along the patch boundary, the model shows a simple mean-field limit of this self-organizing critical behavior, for which we give analytical estimates of the saturation load per fiber values and avalanche size distribution exponent. These are in good agreement with numerical simulation results.

Creep rupture as a non-homogeneous Poissonian process

Creep rupture of heterogeneous materials occurring under constant sub-critical external loads is responsible for the collapse of engineering constructions and for natural catastrophes. Acoustic monitoring of crackling bursts provides microscopic insight into the failure process. Based on a fiber bundle model, we show that the accelerating bursting activity when approaching failure can be described by the Omori law. For long range load redistribution the time series of bursts proved to be a non-homogeneous Poissonian process with power law distributed burst sizes and waiting times. We demonstrate that limitations of experiments such as finite detection threshold and time resolution have striking effects on the characteristic exponents, which have to be taken into account when comparing model calculations with experiments. Recording events solely within the Omori time to failure the size distribution of bursts has a crossover to a lower exponent which is promising for forecasting the imminent catastrophic failure.

Temporal and spacial evolution of bursts in creep rupture

We investigate the temporal and spacial evolution of single bursts and their statistics emerging in heterogeneous materials under a constant external load. Based on a fiber bundle model we demonstrate that when the load redistribution is localized along a propagating crack front, the average temporal shape of pulses has a right-handed asymmetry; however, for long range interaction a symmetric shape with parabolic functional form is obtained. The pulse shape and spatial evolution of bursts proved to be correlated, which can be exploited in materials' testing. The probability distribution of the size and duration of bursts have power law behavior with a crossover to higher exponents as the load is lowered. The crossover emerges due to the competition of the slow and fast modes of local breaking being dominant at low and high loads, respectively.

Noise-induced rupture process: phase boundary and scaling of waiting time distribution

A bundle of fibers has been considered here as a model for composite materials, where breaking of the fibers occur due to a combined influence of applied load (stress) and external noise. Through numerical simulation and a mean-field calculation we show that there exists a robust phase boundary between continuous (no waiting time) and intermittent fracturing regimes. In the intermittent regime, throughout the entire rupture process avalanches of different sizes are produced and there is a waiting time between two consecutive avalanches. The statistics of waiting times follows a Γ distribution and the avalanche distribution shows power-law scaling, similar to what has been observed in the case of earthquake events and bursts in fracture experiments. We propose a prediction scheme that can tell when the system is expected to reach the continuous fracturing point from the intermittent phase.

Avalanche process of the fiber-bundle model with stick-slip dynamics and a variable Young modulus

In order to more accurately describe the fracture process of extensive biological fibers, a fiber-bundle model with stick-slip dynamics and a variable Young modulus is constructed. In this model, the Young modulus of a fiber is assumed to increase or decrease by multiplying with a changing ratio after local sliding events. So, the maximum number of stick-slip events of a single fiber and the changing ratio of the Young modulus are the two key parameters of the model. By means of analytical theory and numerical simulation, the constitutive law, the critical stress, the average size of the largest avalanche, and the avalanche size distribution are shown against the two parameters of the model. From a macroscopic viewpoint, the constitutive curves show different morphologies varying from a local plastic state to a unimodal parabola, while from a microscopic viewpoint, the avalanche size distributions can be well fitted into a power law relationship, which is in accord with the classical fiber-bundle model.

Prediction of the collapse point of overloaded materials by monitoring energy emissions

A bundle of many fibers with stochastically distributed breaking thresholds is considered as a model of composite materials. The fibers are assumed to share the load equally and to obey Hookean elasticity up to the breaking point. The bundle is slightly overloaded which leads to complete failure. We study the properties of emission bursts in which an amount of energy E is released. The analysis shows that the size of the energy bursts has a minimum when the system is halfway from the collapse point.

Statistical model with two order parameters for ductile and soft fiber bundles in nanoscience and biomaterials

Traditional fiber bundles models (FBMs) have been an effective tool to understand brittle heterogeneous systems. However, fiber bundles in modern nano- and bioapplications demand a new generation of FBM capturing more complex deformation processes in addition to damage. In the context of loose bundle systems and with reference to time-independent plasticity and soft biomaterials, we formulate a generalized statistical model for ductile fracture and nonlinear elastic problems capable of handling more simultaneous deformation mechanisms by means of two order parameters (as opposed to one). As the first rational FBM for coupled damage problems, it may be the cornerstone for advanced statistical models of heterogeneous systems in nanoscience and materials design, especially to explore hierarchical and bio-inspired concepts in the arena of nanobiotechnology. Applicative examples are provided for illustrative purposes at last, discussing issues in inverse analysis (i.e., nonlinear elastic polymer fiber and ductile Cu submicron bars arrays) and direct design (i.e., strength prediction).

Kertész line of thermally activated breakdown phenomena

Based on a fiber bundle model we substantially extend the phase-transition analogy of thermally activated breakdown of homogeneous materials. We show that the competition of breaking due to stress enhancement and due to thermal fluctuations leads to an astonishing complexity of the phase space of the system: varying the load and the temperature a phase boundary emerges, separating a Griffith-type regime of abrupt failure analogous to first-order phase transitions from disorder dominated fracture where a spanning cluster of cracks emerges. We demonstrate that the phase boundary is the Kertész line of the system along which thermally activated fracture appears as a continuous phase transition analogous to percolation. The Kertész line has technological relevance setting the boundary of safe operation for construction components under high thermal loads.

Hierarchical simulations for the design of supertough nanofibers inspired by spider silk

Biological materials such as spider silk display hierarchical structures, from nano to macro, effectively linking nanoscale constituents to larger-scale functional material properties. Here, we develop a model that is capable of determining the strength and toughness of elastic-plastic composites from the properties, percentages, and arrangement of its constituents, and of estimating the corresponding dissipated energy during damage progression, in crack-opening control. Specifically, we adopt a fiber bundle model approach with a hierarchical multiscale self-similar procedure which enables to span various orders of magnitude in size and to explicitly take into account the hierarchical topology of natural materials. Hierarchical architectures and self-consistent energy dissipation mechanisms (including plasticity), both omitted in common fiber bundle models, are fully considered in our model. By considering one of the toughest known materials today as an example application, a synthetic fiber composed of single-walled carbon nanotubes and polyvinyl alcohol gel, we compute strength and specific energy absorption values that are consistent with those experimentally observed. Our calculations are capable of predicting these values solely based on the properties of the constituent materials and knowledge of the structural multiscale topology. Due to the crack-opening control nature of the simulations, it is also possible to derive a critical minimal percentage of plastic component needed to avoid catastrophic behavior of the material. These results suggest that the model is capable of helping in the design of new supertough materials.

Avalanche dynamics of fiber bundle models

We present a detailed analytical and numerical study of the avalanche distributions of the continuous damage fiber bundle model (CDFBM). Linearly elastic fibers undergo a series of partial failure events which give rise to a gradual degradation of their stiffness. We show that the model reproduces a wide range of mechanical behaviors. We find that macroscopic hardening and plastic responses are characterized by avalanche distributions, which exhibit an algebraic decay with exponents between 5/2 and 2 different from those observed in mean-field fiber bundle models. We also derive analytically the phase diagram of a family of CDFBM which covers a large variety of potential avalanche size distributions. Our results provide a unified view of the statistics of breaking avalanches in fiber bundle models.

Fiber bundle model with stick-slip dynamics

We propose a generic model to describe the mechanical response and failure of systems which undergo a series of stick-slip events when subjected to an external load. We model the system as a bundle of fibers, where single fibers can gradually increase their relaxed length with a stick-slip mechanism activated by the increasing load. We determine the constitutive equation of the system and show by analytical calculations that on the macroscale a plastic response emerges followed by a hardening or softening regime. Releasing the load, an irreversible permanent deformation occurs which depends on the properties of sliding events. For quenched and annealed disorder of the failure thresholds the same qualitative behavior is found, however, in the annealed case the plastic regime is more pronounced.

Breaking-rate minimum predicts the collapse point of overloaded materials

As a model of composite materials, we choose a bundle of fibers with stochastically distributed breaking thresholds for the individual fibers. The fibers are assumed to share the load equally, and to obey Hookean elasticity right up to the breaking point. We study the evolution of the fiber breaking rate at a constant load in excess of the critical load. The analysis shows that the breaking rate reaches a minimum when the system is half-way from its complete collapse.

Size scaling and bursting activity in thermally activated breakdown of fiber bundles

We study subcritical fracture driven by thermally activated damage accumulation in the framework of fiber bundle models. We show that in the presence of stress inhomogeneities, thermally activated cracking results in an anomalous size effect; i.e., the average lifetime t{f} decreases as a power law of the system size t{f} approximately L{-z}, where the exponent z depends on the external load sigma and on the temperature T in the form z approximately f(sigma/T{3/2}). We propose a modified form of the Arrhenius law which provides a comprehensive description of thermally activated breakdown. Thermal fluctuations trigger bursts of breakings which have a power law size distribution.

Random fiber bundle with many discontinuities in the threshold distribution

We study the breakdown of a random fiber bundle model (RFBM) with n discontinuities in the threshold distribution using the global load sharing scheme. In other words, n+1 different classes of fibers identified on the basis of their threshold strengths are mixed such that the strengths of the fibers in the ith class are uniformly distributed between the values sigma2i-2 and sigma2i-1, where 1< or =i< or =n+1 . Moreover, there is a gap in the threshold distribution between ith and (i+1)-th class. We show that although the critical stress depends on the parameter values of the system, the critical exponents are identical to that obtained in the recursive dynamics of a RFBM with a uniform distribution and global load sharing. The avalanche size distribution, on the other hand, shows a nonuniversal, non-power-law behavior for smaller values of avalanche sizes which becomes prominent only when a critical distribution is approached. We establish that the behavior of the avalanche size distribution for an arbitrary n is qualitatively similar to a RFBM with a single discontinuity in the threshold distribution (n=1) , especially when the density and the range of threshold values of fibers belonging to strongest (n+1)-th class is kept identical in all the cases.

Subcritical crack growth: the microscopic origin of Paris' law

We investigate the origin of Paris' law, which states that the velocity of a crack at subcritical load grows like a power law, da/dt ~ (DeltaK)(m), where DeltaK is the stress-intensity-factor amplitude. Starting from a damage-accumulation function proportional to (Deltasigma)(gamma), Deltasigma being the stress amplitude, we show analytically that the asymptotic exponent m can be expressed as a piecewise-linear function of the exponent gamma, namely, m=6-2gamma for gamma<gamma(c), and m=gamma for gamma> or =gamma(c), reflecting the existence of a critical value gamma(c)=2. We perform numerical simulations to confirm this result for finite sizes. Finally, we introduce bounded disorder in the breaking thresholds and find that below gamma(c) disorder is relevant, i.e., the exponent m is changed, while above gamma(c) disorder is irrelevant.

Continuous damage fiber bundle model for strongly disordered materials

We present an extension of the continuous damage fiber bundle model to describe the gradual degradation of highly heterogeneous materials under an increasing external load. The breaking of a fiber in the model is preceded by a sequence of partial failure events occurring at random threshold values. In order to capture the subsequent propagation and arrest of cracks, furthermore, the disorder of the number of degradation steps of material constituents, the failure thresholds of single fibers, are sorted into ascending order and their total number is a Poissonian distributed random variable over the fibers. Analytical and numerical calculations showed that the failure process of the system is governed by extreme value statistics, which has a substantial effect on the macroscopic constitutive behavior and on the microscopic bursting activity as well.

Energy dissipation statistics in the random fuse model

We study the statistics of the dissipated energy in the two-dimensional random fuse model for fracture under different imposed strain conditions. By means of extensive numerical simulations we compare different ways to compute the dissipated energy. In the case of an infinitely slow driving rate (quasistatic model), we find that the probability distribution of the released energy shows two different scaling regions separated by a sharp energy crossover. At low energies, the probability of having an event of energy E decays as approximately E(-1/2), which is robust and independent of the energy quantifier used (or lattice type). At high energies, fluctuations dominate the energy distribution, leading to a crossover to a different scaling regime, approximately E(-2.75), whenever the released energy is computed over the whole system. On the contrary, strong finite-size effects are observed if we consider only the energy dissipated at microfractures. In a different numerical experiment, the quasistatic dynamics condition is relaxed, so that the system is driven at finite strain load rates, and we find that the energy distribution decays as P(E) approximately E(-1) for all the energy range.

Universality behind Basquin's Law of Fatigue

Basquin's law of fatigue states that the lifetime of the system has a power-law dependence on the external load amplitude, tf approximately sigma 0- alpha, where the exponent alpha has a strong material dependence. We show that in spite of the broad scatter of the exponent alpha, the fatigue fracture of heterogeneous materials exhibits universal features. We propose a generic scaling form for the macroscopic deformation and show that at the fatigue limit the system undergoes a continuous phase transition. On the microlevel, the fatigue fracture proceeds in bursts characterized by universal power-law distributions. We demonstrate that the system dependent details are contained in Basquin's exponent for time to failure, and once this is taken into account, remaining features of failure are universal.

Critical ruptures in a bundle of slowly relaxing fibers

We study the damage enhanced creep rupture of disordered materials by means of a fiber bundle model. Broken fibers undergo a slow stress relaxation modeled by a Maxwell element whose stress exponent m can vary in a broad range. Under global load sharing we show that due to the strength disorder of fibers, the lifetime t(f) of the bundle has sample-to-sample fluctuations characterized by a log-normal distribution independent of the type of disorder. We determine the Monkman-Grant relation of the model and establish a relation between the rupture life t(f) and the characteristic time t(m) of the intermediate creep regime of the bundle where the minimum strain rate is reached, making possible reliable estimates of t(f) from short term measurements. Approaching macroscopic failure, the deformation rate has a finite time power law singularity whose exponent is a decreasing function of m. On the microlevel the distribution of waiting times is found to have a power law behavior with m-dependent exponents different below and above the critical load of the bundle. Approaching the critical load from above, the cutoff value of the distributions has a power law divergence whose exponent coincides with the stress exponent of Maxwell elements.

Energy bursts in fiber bundle models of composite materials

A bundle of many fibers with stochastically distributed breaking thresholds for the individual fibers is considered as a model of composite materials. The bundle is loaded until complete failure, to capture the failure scenario of composite materials under external load. The fibers are assumed to share the load equally, and to obey Hookean elasticity right up to the breaking point. We determine the distribution of bursts in which an amount of energy E is released. The energy distribution follows asymptotically a universal power law E(-5/2) , for any statistical distribution of fiber strengths. A similar power law dependence is found in some experimental acoustic emission studies of loaded composite materials.

Damage process of a fiber bundle with a strain gradient

We study the damage process of fiber bundles in a wedge-shape geometry which ensures a constant strain gradient. To obtain the wedge geometry we consider the three-point bending of a bar, which is modeled as two rigid blocks glued together by a thin elastic interface. The interface is discretized by parallel fibers with random failure thresholds, which become elongated when the bar is bent. Analyzing the progressive damage of the system we show that the strain gradient results in a rich spectrum of novel behavior of fiber bundles. We find that for weak disorder an interface crack is formed as a continuous region of failed fibers. Ahead of the crack a process zone develops which proved to shrink with increasing deformation, making the crack tip sharper as the crack advances. For strong disorder, failure of the system occurs as a spatially random sequence of breakings. Damage of the fiber bundle proceeds in bursts whose size distribution shows a power law behavior with a crossover from an exponent 2.5 to 2.0 as the disorder is weakened. The size of the largest burst increases as a power law of the strength of disorder with an exponent 23 and saturates for strongly disordered bundles.

Critical behavior of random fibers with mixed Weibull distribution

A random fiber bundle model with a mixed Weibull distribution is studied under the global load sharing scheme. The mixed model consists of two sets of fibers. The threshold strength of one set of fibers is randomly chosen from a Weibull distribution with a particular Weibull index, and another set of fibers with a different index. The mixing tunes the critical stress of the bundle and the variation of critical stress with the amount of mixing is determined using a probabilistic method where the external load is increased quasistatically. In a special case which we illustrate, the critical stress is found to vary linearly with the mixing parameter. The critical exponents and power-law behavior of burst avalanche size distribution is found to remain unaltered due to mixing.

Effect of discontinuity in the threshold distribution on the critical behavior of a random fiber bundle

The critical behavior of a random fiber bundle model with mixed uniform distribution of threshold strengths and global load sharing rule is studied with a special emphasis on the nature of distribution of avalanches for different parameters of the distribution. The discontinuity in the threshold strength distribution of fibers nontrivially modifies the critical stress as well as puts a restriction on the allowed values of parameters for which the recursive dynamics approach holds good. The discontinuity leads to a nonuniversal behavior in the avalanche size distribution for smaller values of avalanche size. We observe that apart from the mean field behavior for larger avalanches, a new behavior for smaller avalanche size is observed as a critical threshold distribution is approached. The phenomenological understanding of the above result is provided using the exact analytical result for the avalanche size distribution. Most interestingly, the prominence of nonuniversal behavior in avalanche size distribution depends on the system parameters.

Local load sharing fiber bundles with a lower cutoff of strength disorder

We study the failure properties of fiber bundles with a finite lower cutoff of the strength disorder varying the range of interaction between the limiting cases of completely global and completely local load sharing. Computer simulations revealed that at any range of load redistribution there exists a critical cutoff strength where the macroscopic response of the bundle becomes perfectly brittle, i.e., linearly elastic behavior is obtained up to global failure, which occurs catastrophically after the breaking of a small number of fibers. As an extension of recent mean field studies [Phys. Rev. Lett. 95, 125501 (2005)], we demonstrate that approaching the critical cutoff, the size distribution of bursts of breaking fibers shows a crossover to a universal power law form with an exponent 3/2 independent of the range of interaction.

Crossover behavior in failure avalanches

Composite materials, with statistically distributed thresholds for breakdown of individual elements, are considered. During the failure process of such materials under external stress (load or voltage), avalanches consisting of simultaneous rupture of several elements occur, with a distribution D(Delta) of the magnitude Delta of such avalanches. The distribution is typically a power law D(Delta) proportional to Delta (-xi). For the systems we study here, a crossover behavior is seen between two power laws, with a small exponent xi in the vicinity of complete breakdown and a larger exponent xi for failures away from the breakdown point. We demonstrate this analytically for bundles of many fibers where the load is uniformly distributed among the surviving fibers. In this case xi=3/2 near the breakdown point and xi=5/2 away from it. The latter is known to be the generic behavior. This crossover is a signal of imminent catastrophic failure of the material. Near the breakdown point, avalanche statistics show nontrivial finite size scaling. We observe similar crossover behavior in a network of electric fuses, and find xi=2 near the catastrophic failure and xi=3 away from it. For this fuse model power dissipation avalanches show a similar crossover near breakdown.

Failure process of a bundle of plastic fibers

We present an extension of fiber bundle models considering that failed fibers still carry a fraction 0 < or = alpha < or = 1 of their failure load. The value of alpha interpolates between the perfectly brittle failure (alpha = 0) and perfectly plastic behavior (alpha = 1) of fibers. We show that the finite load bearing capacity of broken fibers has a substantial effect on the failure process of the bundle. In the case of global load sharing it is found that for alpha --> 1 the macroscopic response of the bundle becomes perfectly plastic with a yield stress equal to the average fiber strength. On the microlevel, the size distribution of avalanches has a crossover from a power law of exponent approximately 2.5 to a faster exponential decay. For localized load sharing, computer simulations revealed a sharp transition at a well-defined value alpha(c) from a phase where macroscopic failure occurs due to localization as a consequence of local stress enhancements, to another one where the disordered fiber strength dominates the damage process. Analyzing the microstructure of damage, the transition proved to be analogous to percolation. At the critical point alpha(c), the spanning cluster of damage is found to be compact with a fractal boundary. The distribution of bursts of fiber breakings shows a power-law behavior with a universal exponent approximately 1.5 equal to the mean-field exponent of fiber bundles of critical strength distributions. The model can be relevant to understand the shear failure of glued interfaces where failed regions can still transmit load by remaining in contact.