An efficient damage-plasticity DEM contact model for highly porous rocks

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.