A novel discrete element method (DEM) model is proposed to better reproduce the behaviour of porous soft rocks. With the final goal of simulating pile penetration problems efficiency and scalability are two underlining features. The contact model is based on the macro-element theory and employs damage laws to govern the plastic deformations developing at the microscale. To attain (i) high porosity states, (ii) represent irregular shaped grains and (iii) incorporate the physical presence of bond fragments, the model is cast within a far-field interaction framework allowing for non-overlapping particles to transmit forces. After presenting a calibration procedure, the model is used to replicate the behaviour of Maastricht calcarenite. In particular, the mechanical response of this calcarenite is explored within the critical state theory framework. Finally, the efficiency, performance and scalability of the model is tested by simulating physical model experiments of cone-ended penetration tests in Maastricht calcarenite from the literature. To boost efficiency of the 3D numerical simulations, a coupled DEM-FDM (Finite Differential Method) framework is used. The good fit between the experimental and numerical results suggest that the new model can be used to unveil microscopic mechanism controlling the macroscopic response of soft-rock/structure interaction problems.

This study presents a fully coupled thermo-hydro-mechanical (THM) constitutive model for clay rocks. The model is formulated within the elastic-viscoplasticity framework, which considers nonlinearity and softening after peak strength, anisotropy of stiffness and strength, as well as permeability variation due to damage. In addition, the mechanical properties are coupled with thermal phenomena and accumulated plastic strains. The adopted nonlocal and viscoplastic approaches enhance numerical efficiency and provide the possibility to simulate localization phenomena. The model is validated against experimental data from laboratory tests conducted on Callovo-Oxfordian (COx) claystone samples that are initially unsaturated and under suction. The tests include a thermal phase where the COx specimens are subjected to different temperature increases. A good agreement with experimental data is obtained. In addition, parametric analyses are carried out to investigate the influence of the hydraulic boundary conditions (B.C.) and post-failure behavior models on the THM behavior evolution. It is shown that different drainage conditions affect the thermally induced pore pressures that, in turn, influence the onset of softening. The constitutive model presented constitutes a promising approach for simulating the most important features of the THM behavior of clay rocks. It is a tool with a high potential for application to several relevant case studies, such as thermal fracturing analysis of nuclear waste disposal systems. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Micron-scale crack propagation in red-bed soft rocks under hydraulic action is a common cause of engineering disasters due to damage to the hard rock-soft rock-water interface. Previous studies have not provided a theoretical analysis of the length, inclination angle, and propagation angle of micron-scale cracks, nor have they established appropriate criteria to describe the crack propagation process. The propagation mechanism of micron-scale cracks in red-bed soft rocks under hydraulic action is not yet fully understood, which makes it challenging to prevent engineering disasters in these types of rocks. To address this issue, we have used the existing generalized maximum tangential stress (GMTS) and generalized maximum energy release rate (GMERR) criteria as the basis and introduced parameters related to micron-scale crack propagation and water action. The GMTS and GMERR criteria for micron-scale crack propagation in red-bed soft rocks under hydraulic action (abbreviated as the Wmic-GMTS and Wmic-GMERR criteria, respectively) were established to evaluate micron-scale crack propagation in red-bed soft rocks under hydraulic action. The influence of the parameters was also described. The process of micron-scale crack propagation under hydraulic action was monitored using uniaxial compression tests (UCTs) based on digital image correlation (DIC) technology. The study analyzed the length, propagation and inclination angles, and mechanical parameters of micron-scale crack propagation to confirm the reliability of the established criteria. The findings suggest that the Wmic-GMTS and Wmic-GMERR criteria are effective in describing the micron-scale crack propagation in red-bed soft rocks under hydraulic action. This study discusses the mechanism of micron-scale crack propagation and its effect on engineering disasters under hydraulic action. It covers topics such as the internal-external weakening of nano-scale particles, lateral propagation of micron-scale cracks, weakening of the mechanical properties of millimeter-scale soft rocks, and resulting interface damage at the engineering scale. The study provides a theoretical basis for the mechanism of disasters in red-bed soft-rock engineering under hydraulic action. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V.