Mesh topology-based spurious pressure stabilization in 3D finite elasticity using Voronoi tessellations

Abstract In this paper, we present a mesh topology-based stabilization approach to suppress spurious pressure modes in 3D nearly-incompressible finite elasticity. The focus lies on a mixed formulation with lowest-order approximation for the displacement and pressure fields. Motivated by the fact that the popular H1/P0 element does not fulfill the inf-sup condition, all possible local spurious pressure modes are derived on a patch of elements. The nullspace method is used to determine all spurious pressure solutions. From this, the topological requirements of the finite element mesh are established. We conclude that no more than four elements are allowed to intersect in the same vertex to overcome local checkerboarding. To fulfill this requirement, we employ non-degenerate 3D Voronoi diagrams with several different site distributions. These result in random, centroidal, and honeycomb Voronoi meshes. The resulting convex polyhedral elements are discretized by a polyhedral mixed finite element based on the lowest possible interpolation pair. The numerical examples illustrate that spurious pressure modes do not occur for any degree of mesh refinement as long as the topological mesh requirements are met. Furthermore, it is shown that the numerical inf-sup test is passed. By violating the topological requirements, it is shown that a stable pressure field cannot be guaranteed and the checkerboard phenomenon is provoked.