Discrete element model for general polyhedra

We present a version of the Discrete Element Method considering the particles as rigid polyhedra. The Principle of Virtual Work is employed as basis for a multibody dynamics model. Each particle surface is split into sub-regions, which are tracked for contact with other sub-regions of neighboring particles. Contact interactions are modeled pointwise, considering vertex-face, edge-edge, vertex-edge and vertex-vertex interactions. General polyhedra with triangular faces are considered as particles, permitting multiple pointwise interactions which are automatically detected along the model evolution. We propose a combined interface law composed of a penalty and a barrier approach, to fulfill the contact constraints. Numerical examples demonstrate that the model can handle normal and frictional contact effects in a robust manner. These include simulations of convex and non-convex particles, showing the potential of applicability to materials with complex shaped particles such as sand and railway ballast.