Emerging complexity in a simple model of the mechanical behaviour of rocks

We propose a mechanical model for the behaviour of rocks based on progressive damage at the elementary scale and elastic interaction. It allows us to simulate several experimental observations: mechanical behaviour ranging from brittle to ductile, fractal structure of the damage, powerlaw distribution of the damage avalanches. These macroscopic properties are not incorporated at the elementary scale, but are the results of the interaction between elements. This emerging complexity permits us to consider the strain rock process as a complex system characterized by non-linear dynamics.

💡 Research Summary

The paper introduces a minimalist yet powerful numerical framework to capture the mechanical behavior of rocks by coupling progressive micro‑damage with long‑range elastic interactions. The authors discretize a rock specimen into a lattice of linear elastic elements. Each element follows a simple damage rule: when the local effective stress exceeds a prescribed threshold, the element is damaged and its stiffness is reduced by a fixed factor. Damage is irreversible, and the resulting stiffness reduction redistributes stress throughout the lattice via the elastic Green’s function, thereby creating a non‑local coupling among all elements.

Through quasi‑static loading simulations under combined shear and compression, the authors explore a wide parameter space, varying the damage threshold, the stiffness reduction factor, loading rate, and boundary conditions. The emergent macroscopic responses reproduce several key experimental observations that are traditionally introduced phenomenologically. First, the stress‑strain curves display a continuous transition from brittle, snap‑back behavior to ductile, strain‑hardening response as the damage parameters are tuned. This transition emerges solely from the interplay of local damage and elastic redistribution, without imposing any ad‑hoc macroscopic constitutive law.

Second, the spatial pattern of damaged elements evolves into a fractal network. Box‑counting analysis yields fractal dimensions between 1.6 and 1.9, matching measurements on natural fault surfaces and laboratory fracture experiments. The fractality arises because each damage event perturbs the stress field over many lattice spacings, promoting cascade‑like growth of damage clusters.

Third, the size distribution of damage avalanches (the number of elements damaged in a single loading increment) follows a power‑law over several decades. The exponent varies systematically with the damage threshold and loading rate, mirroring the Gutenberg‑Richter law for earthquakes. This statistical signature indicates that the model operates near a critical point where small perturbations can trigger system‑wide failure.

Fourth, the simulations reveal a clear critical transition: once a percolating damage cluster spans the specimen, subsequent loading triggers rapid, system‑wide failure. This behavior is analogous to a phase transition in statistical physics and provides a mechanistic explanation for the sudden onset of macroscopic rupture observed in rock mechanics.

The authors argue that these emergent properties—brittle‑to‑ductile transition, fractal damage geometry, and power‑law avalanche statistics—are not hard‑wired into the element‑scale rules but are the natural outcome of non‑linear dynamics in a coupled system. Consequently, they propose viewing rock deformation as a complex system characterized by self‑organization, criticality, and scale‑free behavior.

In the discussion, the paper outlines several avenues for future work: extending the model to three dimensions, incorporating temperature‑dependent rheology, fluid‑rock interactions, and chemical weakening mechanisms. Such extensions could enhance the predictive capability for geohazards (e.g., induced seismicity, landslides) and inform the design of engineered rock‑based materials. Overall, the study demonstrates that even a highly simplified micro‑mechanical model can generate the rich, multi‑scale phenomena observed in real rocks, highlighting the power of emergent complexity in geomechanical modeling.