Systematic investigation into the role of particle multi-level morphology in determining the shear behavior of granular materials via DEM simulation

This study aims to numerically investigate the dependence of shear behaviors of granular materials on particle multi-level morphology, encompassing form (mean of elongation and flatness), angularity, and roughness. To this end, four series of particles—ellipsoidal particles with varying elongation indices (EI) and flatness indices (FI), and equiaxed concave particles with varying angularity indices (AI) and roughness (RG)—are adopted to isolate the effects of particle morphology at each level. Through discrete element method (DEM) simulations, granular assemblies with varying values of EI, FI, AI, and RG are consolidated to achieve the maximum packing density under the same confining pressure. Subsequently, they are subjected to triaxial drained compression tests until reaching a critical state. Based on numerical results, the independent effects of EI, FI, AI, and RG on the macro- and micro-mechanical response of granular materials are examined and compared in detail. Furthermore, using the 'Stress-force-fabric' (SFF) relationship, the anisotropic origins of the peak- and critical-state shear strengths induced by EI, FI, AI, and RG are also analyzed.