A Voronoi strain-based method for granular materials and continua

Abstract

In a recent article (Frenning in Comp Part Mech 24:1–4, 2021), we demonstrated that a Delaunay-based strain estimate could be used as a starting point for the development of a particle-based method for continua. In this article, we argue that the Voronoi diagram, dual to the previously used Delaunay tetrahedralization, provides a more natural description of the underlying particulate system. For this reason, a Voronoi-based estimate of the deformation gradient is derived and used to the same effect. Although the gradient vectors cease to be antisymmetric, sums over nearest neighbors vanish, which results in a formulation that not only is linearly complete but also satisfies the patch test irrespective of initial particle placement. Pairwise forces, inferred from the local (nonaffine) deformation of each bond or contact, impart a physical stabilization