This paper presents a dynamic formulation for the simulation of nearly incompressible structures using a mixed finite element method with equal-order interpolation pairs. Specifically, the nodal unknowns are the displacement and the volumetric strain component, something that makes possible the reconstruction of the complete stain at the integration point level and thus enables the use of strain-driven constitutive laws. Furthermore, we also discuss the resulting eigenvalue problem and how it can be applied for the modal analysis of linear elastic solids. The article puts special emphasis on the stabilisation technique used, which becomes crucial in the resolution of the generalised eigenvalue problem. In particular, we prove that using a variational multiscale method assuming the sub-grid scales to lie in the finite element space orthogonal to that of the approximation, namely the Orthogonal Sub-Grid Scales (OSGS), results in a convenient linear and symmetric generalised eigenvalue problem. The correctness, convergence and performance of the method are proven by solving a series of two- and three-dimensional examples.
Read full abstract