In this paper, we propose a variationally consistent technique for decreasing the maximum eigenfrequencies of structural dynamics related finite element formulations. Our approach is based on adding a symmetric positive-definite term to the consistent mass matrix that follows from the integral of the traction jump across element boundaries. The added term is weighted by a small factor, for which we derive a suitable, and simple, element-local parameter choice. We perform numerical experiments for the linear wave equation in one and two dimensions, for quadrilateral elements and triangular elements, and for up to fourth order polynomial basis functions. Despite the increase in critical time-step size, we do not observe adverse effects in terms of spatial accuracy and orders of convergence. To extend the method to non-linear problems, we introduce a linear approximation. Our three-dimensional experiments with tetrahedral and hexahedral elements show that a sizeable increase in critical time-step size can be achieved while only causing minor (even beneficial) influences on the dynamic response.
Computer Methods in Applied Mechanics and Engineering, 2022.
Authors: S. K. F. Stoter, T.-H. Nguyen, R.R. Hiemstra, D. Schillinger.