Crack propagation through phase-separated glasses: effect of the characteristic size of disorder

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Unified scenario for the morphology of crack paths in two-dimensional disordered solids

A combined experimental and numerical investigation of the roughness of intergranular cracks in two-dimensional disordered solids is presented. We focus on brittle materials for which the characteristic length scale of damage is much smaller than the grain size. Surprisingly, brittle cracks do not follow a persistent path with a roughness exponent ζ≈0.6-0.7 as reported for a large range of materials. Instead, we show that they exhibit monoaffine scaling properties characterized by a roughness exponent ζ=0.50±0.05, which we explain theoretically from linear elastic fracture mechanics. Our findings support the description of the roughening process in two-dimensional brittle disordered solids by a random walk. Furthermore, they shed light on the failure mechanism at the origin of the persistent behavior with ζ≈0.6-0.7 observed for fractures in other materials, suggesting a unified scenario for the geometry of crack paths in two-dimensional disordered solids.

Higher Toughness of Metal-nanoparticle-implanted Sodalime Silicate Glass with Increased Ductility

In this report, we propose a novel framework for toughening brittle oxide glass originated from enhanced ductility by implanting a secondary material comprising different mechanical properties. To do so, copper-metal nanoparticles are implanted into the subsurface layer of commercial soda-lime silica glass by using the electrofloat method. The crack initiation load of the implanted glass is found to be comparable to the glass chemically strengthened in ordinary tempering conditions. By observing crack propagation and stress distribution from cross-section, it is found that the crack propagation stops within the metal nanoparticle implanted layer, due to the stress dissipation or relaxation. The copper-implanted glass shows improved toughness with decreased hardness. The toughening mechanism of the composite glass is theoretically studied using molecular dynamics calculations on an amorphous silica model with copper nanoparticles embedded, and Peridynamics fracture simulations for indentation on a glass sheet model whose surface was implicitly modeled as the copper-implanted oxide glass. The experimentally observed phenomena of intrinsic toughening were well explained by the series of the conducted simulations.

Modes of failure in disordered solids

The two principal ingredients determining the failure modes of disordered solids are the strength of heterogeneity and the length scale of the region affected in the solid following a local failure. While the latter facilitates damage nucleation, the former leads to diffused damage-the two extreme natures of the failure modes. In this study, using the random fiber bundle model as a prototype for disordered solids, we classify all failure modes that are the results of interplay between these two effects. We obtain scaling criteria for the different modes and propose a general phase diagram that provides a framework for understanding previous theoretical and experimental attempts of interpolation between these modes. As the fiber bundle model is a long-standing model for interpreting various features of stressed disordered solids, the general phase diagram can serve as a guiding principle in anticipating the responses of disordered solids in general.

Effect of the porosity on the fracture surface roughness of sintered materials: from anisotropic to isotropic self-affine scaling

To unravel how the microstructure affects the fracture surface roughness in heterogeneous brittle solids like rocks or ceramics, we characterized the roughness statistics of postmortem fracture surfaces in homemade materials of adjustable microstructure length scale and porosity, obtained by sintering monodisperse polystyrene beads. Beyond the characteristic size of disorder, the roughness profiles are found to exhibit self-affine scaling features evolving with porosity. Starting from a null value and increasing the porosity, we quantitatively modify the self-affine scaling properties from anisotropic (at low porosity) to isotropic (for porosity >10%).

(Finite) statistical size effects on compressive strength

The larger structures are, the lower their mechanical strength. Already discussed by Leonardo da Vinci and Edmé Mariotte several centuries ago, size effects on strength remain of crucial importance in modern engineering for the elaboration of safety regulations in structural design or the extrapolation of laboratory results to geophysical field scales. Under tensile loading, statistical size effects are traditionally modeled with a weakest-link approach. One of its prominent results is a prediction of vanishing strength at large scales that can be quantified in the framework of extreme value statistics. Despite a frequent use outside its range of validity, this approach remains the dominant tool in the field of statistical size effects. Here we focus on compressive failure, which concerns a wide range of geophysical and geotechnical situations. We show on historical and recent experimental data that weakest-link predictions are not obeyed. In particular, the mechanical strength saturates at a nonzero value toward large scales. Accounting explicitly for the elastic interactions between defects during the damage process, we build a formal analogy of compressive failure with the depinning transition of an elastic manifold. This critical transition interpretation naturally entails finite-size scaling laws for the mean strength and its associated variability. Theoretical predictions are in remarkable agreement with measurements reported for various materials such as rocks, ice, coal, or concrete. This formalism, which can also be extended to the flowing instability of granular media under multiaxial compression, has important practical consequences for future design rules.

Experimental study of the effect of disorder on subcritical crack growth dynamics

The growth dynamics of a single crack in a heterogeneous material under subcritical loading is an intermittent process, and many features of this dynamics have been shown to agree with simple models of thermally activated rupture. In order to better understand the role of material heterogeneities in this process, we study the subcritical propagation of a crack in a sheet of paper in the presence of a distribution of small defects such as holes. The experimental data obtained for two different distributions of holes are discussed in the light of models that predict the slowing down of crack growth when the disorder in the material is increased; however, in contradiction with these theoretical predictions, the experiments result in longer lasting cracks in a more ordered scenario. We argue that this effect is specific to subcritical crack dynamics and that the weakest zones between holes at close distance to each other are responsible for both the acceleration of the crack dynamics and the slightly different roughness of the crack path.

Quantitative prediction of effective toughness at random heterogeneous interfaces

The propagation of an adhesive crack through an anisotropic heterogeneous interface is considered. Tuning the local toughness distribution function and spatial correlation is numerically shown to induce a transition between weak to strong pinning conditions. While the macroscopic effective toughness is given by the mean local toughness in the case of weak pinning, a systematic toughness enhancement is observed for strong pinning (the critical point of the depinning transition). A self-consistent approximation is shown to account very accurately for this evolution, without any free parameter.

Fracture roughness in three-dimensional beam lattice systems

We study the scaling of three-dimensional crack roughness using large-scale beam lattice systems. Our results for prenotched samples indicate that the crack surface is statistically isotropic, with the implication that experimental findings of anisotropy of fracture surface roughness in directions parallel and perpendicular to crack propagation is not due to the scalar or vectorial elasticity of the model. In contrast to scalar fuse lattices, beam lattice systems do not exhibit anomalous scaling or an extra dependence of roughness on system size. The local and global roughness exponents (ζ(loc) and ζ, respectively) are equal to each other, and the three-dimensional crack roughness exponent is estimated to be ζ(loc)=ζ=0.48±0.03 . This closely matches the roughness exponent observed outside the fracture process zone. The probability density distribution p[Δh(ℓ)] of the height differences Δh(ℓ)=[h(x+ℓ)-h(x)] of the crack profile follows a Gaussian distribution, in agreement with experimental results.

Effects of finite probe size on self-affine roughness measurements

The roughness of fracture surfaces exhibits self-affinity for a wide variety of materials and loading conditions. The universality and the range of scales over which this regime extends are still debated. The topography of these surfaces is however often investigated with a finite contact probe. In this case, we show that the correlation function of the roughness can only be measured down to a length scale Deltax{c} which depends on the probe size R, the Hurst exponent zeta of the surface and its topothesy l, and exhibits spurious behavior at smaller scales. First, we derive the dependence of Deltax{c} on these parameters from a simple scaling argument. Then, we verify this dependence numerically. Finally, we establish the relevance of this analysis from AFM measurements on an experimental glass fracture surface and provide a metrological procedure for roughness measurements.

Fracture surface statistics of filled elastomers

Roughness profiles of fracture surfaces formed as a result of the fast crack propagation through a filled rubber were analyzed by means of the height-height correlation functions. The fracture surface was found to be anisotropic in a certain domain of values of length scales; i.e., different values of roughness exponents are observed across and along the crack propagation direction. A two-dimensional analysis reveals a Family-Vicsek scaling in this domain characterized as well by two exponents. These characteristic values of the roughness exponents are found to be close to those observed for fracture surfaces of certain nonrubber materials at length scales smaller than the size of the fracture process zone. Hence, a ductile fracture process can be surmised to occur within the domain of the corresponding length scales.