Numerical simulation of particle breakage of angular particles using combined DEM and FEM

One of the effective parameters of the behavior of rockfill materials is particle breakage. As a result of particle breakage, both the stress–strain and deformability of materials change significantly. In this article, a novel approach for the two-dimensional numerical simulation of the phenomenon in rockfill (sharp-edge particles) has been developed using combined DEM and FEM. All particles are simulated by the discrete element method (DEM) as an assembly and after each step of DEM analysis, each particle is separately modeled by FEM to determine its possible breakage. If the particle fulfilled the proposed breakage criteria, the breakage path is assumed to be a straight line and is determined by a full finite element stress–strain analysis within that particle and two new particles are generated, replacing the original particle. These procedures are carried out on all particles in each time step of the DEM analysis. Novel approach for the numeric of breakage appears to produce reassuring physically consistent results that improve earlier made unnecessary simplistic assumptions about breakage. To evaluate the effect of particle breakage on rockfill's behavior, two test series with and without breakable particles have been simulated under a biaxial test with different confining pressures. Results indicate that particle breakage reduces the internal friction but increases the deformability of rockfill. Review of the v–p variation of the simulated samples shows that the specific volume has initially been reduced with the increase of mean pressures and then followed by an increase. Also, the increase of stress level reduces the growing length of the v–p path and it means that the dilation is reduced. Generally, any increase of confining stress decreases the internal friction angle of the assembly and the sample fail at higher values of axial stresses and promotes an increase in the deformability. The comparison between the simulations and the reported experimental data shows that the numerical simulation and experimental results are qualitatively in agreement. Overall the presented results show that the proposed model is capable with more accuracy to simulate the particle breakage in rockfill.