An Improved Contact Force Model of Polyhedral Elements for the Discrete Element Method

The discrete element method (DEM) serves as a robust tool for simulating the mechanical behavior of granular materials. The accuracy of the DEM simulation is determined by the calculation of contact forces between elements. Compared to spherical elements, the contact modes of polyhedral elements are more complex. In this work, the contact force between polyhedral elements under different contact modes was investigated by experiment. Based on the experimental results, the normal stiffness coefficient in the Cundall’s contact model was modified. The improved contact force model was then applied in the DEM simulation and validated by means of comparison with the results of packing experiments. The research results demonstrate that the improved contact force model can be effectively applied to the simulation of different contact modes between polyhedral elements. The results of the packing experiment highly coincide with the results of the DEM simulation, which confirms the accuracy and reliability of the improved contact force model.