DEM modeling of the one-dimensional compression of sands incorporating a statistical particle fragmentation scheme

This paper presents a novel framework of modeling crushable granular materials under mechanical loadings based on the discrete element method (DEM). The framework is featured with the construction of the one-to-one model in which every particle in a physical experiment has its own numerical twin and allows the modeling of irregular shaped fragments during the continuous breakage process. First, image processing techniques and spherical harmonic (SH) analysis were adopted, respectively, to segment and label particles and to construct a one-to-one model mathematically in DEM. Then, a particle crushing criterion based on the maximum inter-particle contact force was used to predict the crushing events, showing fitting results that agreed very well with a large number of single particle crushing tests. Next, a statistical approach for the generation of particle fragmentation modes of a given type of sand particles based on the principal component analysis (PCA) was proposed. The aim of the PCA was to analyze the statistical trends of the coefficient matrix, which was composed of the SH coefficients of all the particles involved in the analysis. Finally, a successful modeling of a particle crushing event was achieved by replacing the particle, which was judged by the crushing criterion to undergo crushing, with a few sub-particles chosen randomly from a specific fragment template constructed using the micro-computed tomography (micro-CT) data.