Hybrid-mixed discretization of elasto-dynamic contact problems using consistent singular mass matrices

A. Tkachuk, B. I. Wohlmuth, M. Bischoff

Research output: Contribution to journalArticlepeer-review

14 Scopus citations


An alternative spatial semi-discretization of dynamic contact based on a modified Hamilton's principle is proposed. The modified Hamilton's principle uses displacement, velocity and momentum as variables, which allows their independent spatial discretization. Along with a local static condensation for velocity and momentum, it leads to an approach with a hybrid-mixed consistent mass matrix. An attractive feature of such a formulation is the possibility to construct hybrid singular mass matrices with zero components at those nodes where contact is collocated. This improves numerical stability of the semi-discrete problem: the differential index of the underlying differential-algebraic system is reduced from 3 to 1, and spurious oscillations in the contact pressure, which are commonly reported for formulations with Lagrange multipliers, are significantly reduced. Results of numerical experiments for truss and Timoshenko beam elements are discussed. In addition, the properties of the novel discretization scheme for an unconstrained dynamic problem are assessed by a dispersion analysis.

Original languageEnglish
Pages (from-to)473-493
Number of pages21
JournalInternational Journal for Numerical Methods in Engineering
Issue number5
StatePublished - 4 May 2013


  • Finite element method
  • Hamiltonian
  • Impact
  • Index reduction
  • Solids
  • Spectral
  • Stabilization


Dive into the research topics of 'Hybrid-mixed discretization of elasto-dynamic contact problems using consistent singular mass matrices'. Together they form a unique fingerprint.

Cite this