Discontinuous variational time integrators for complex multibody collisions. (3rd October 2014)

• Record Type: Journal Article
• Authors: Johnson, G.; Leyendecker, S.; Ortiz, M.
• Is Part Of: International journal for numerical methods in engineering. Volume 100:Number 12(2014)
• Page Start: 871
• Page End: 913
• Publication Date: 2014-10-03
• Subjects: Numerical analysis -- Periodicals; Engineering mathematics -- Periodicals; 620.001518
• DOI: 10.1002/nme.4764
• Languages: English
• ISSNs: 0029-5981

Abstract

<abstract abstract-type="main" id="nme4764-abs-0001"> <title>SUMMARY</title> <p id="nme4764-para-0001">The objective of the present work is to formulate a new class of <italic>discontinuous variational time integrators</italic> that allow the system to adopt two possibly different configurations at each sampling time <italic>t</italic><sub><italic>k</italic></sub>, representing predictor and corrector configurations of the system. The resulting sequence of configuration pairs then represents a <italic>discontinuous</italic>—or <italic>non‐classical</italic>—trajectory. Continuous or classical trajectories are recovered simply by enforcing a continuity constraint at all times. In particular, in systems subject to one‐sided contact constraints simulated via discontinuous variational time integrators, the predictor configuration is not required to satisfy the one‐sided constraints, whereas the corrector configuration is obtained by a closest‐point projection (CPP) onto the admissible set. The resulting trajectories are generally discontinuous, or <italic>non‐classical</italic>, but are expected to converge to classical or continuous solutions for decreasing time steps. We account for dissipation, including friction, by means of a discrete Lagrange–d'Alembert principle, and make extensive use of the <italic>spacetime formalism</italic> in order to ensure exact energy conservation in conservative systems, and the right rate of energy decay in dissipative systems. The structure, range and scope of the discontinuous variational time integrators, and their accuracy characteristics are illustrated by means of examples of application concerned with rigid multibody dynamics. Copyright © 2014 John Wiley & Sons, Ltd.</p> </abstract>