Filtering of high modal frequencies for stable real-time explicit integration of deformable objects using the Finite Element Method

Abstract

The behavior, performance, and run-time of mechanical simulations in interactive virtual surgery depend heavily on the type of numerical differential equation solver used to integrate in time the dynamic equations obtained from simulation methods, such as the Finite Element Method. Explicit solvers are fast but only conditionally stable. The condition number of the stiffness matrix limits the highest possible time step. This limit is related to the geometrical properties of the underlying mesh, such as element shape and size. In fact, it can be governed by a small set of ill-shaped elements. For many applications this issue can be solved a priori by a careful meshing. However, when meshes are cut during interactive surgery simulation, it is difficult and computationally expensive to control the quality of the resulting elements. As an alternative, we propose to modify the elemental stiffness matrices directly in order to ensure stability. In this context, we first investigate the behavior of the eigenmodes of the elemental stiffness matrix in a Finite Element Method. We then propose a simple filter to reduce high model frequencies and thus allow larger time steps, while maintaining the general mechanical behavior.