A nonlocal kernel-based continuum damage model for compaction band formation in porous sedimentary rock

Enrique M. del Castillo,Jun Geng,Ronaldo I. Borja
DOI: https://doi.org/10.1007/s00466-024-02540-x
IF: 4.391
2024-10-05
Computational Mechanics
Abstract:We propose a novel continuum damage model for describing the propagation of compaction bands formed by grain crushing in porous sedimentary rock using the kernel-based approximation and discretization of the meshfree Lagrangian smoothed particle hydrodynamics (SPH) method. In the model, damage is assumed to be caused by grain crushing, and the effects of plasticity are incorporated through a critical state type model which depends on a degradation function designed to capture the abrupt onset of pore collapse following grain crushing. In doing so, we account for the two main forces understood to drive compaction band formation, brittle failure at the grain scale, and plastic dissipation. Through the smoothing length parameter, the SPH method possesses nonlocal properties intrinsic to the method, allowing for the development of a nonlocal integral form of the equivalent strain variable used to compute the damage in addition to a distinct gradient enhanced approach that when discretized using SPH is endowed with two characteristic length scales. We compare and contrast both formulations, as well as the roles of the different length scales. We evaluate our model by performing numerical experiments on Bentheim and Berea sandstone as well as on Tuffeau de Maastricht, in both notched and unnotched specimens, matching experimentally observed results such as a transition from high-angled bands to horizontal bands with increasing confining pressure, and similar compaction band styles to those visible in the field. We lastly consider the effects of material heterogeneity on samples of Tuffeau de Maastricht noting that our model produces precursory low-magnitude localization events which may develop into persistent compaction bands involving localized grain crushing, which is consistent with our understanding of compaction band formation.
mechanics,mathematics, interdisciplinary applications
What problem does this paper attempt to address?