Premium
Unconditionally energy stable implicit time integration: application to multibody system analysis and design
Author(s) -
Chen Shanshin,
Hansen John M.,
Tortorelli Daniel A.
Publication year - 2000
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/(sici)1097-0207(20000630)48:6<791::aid-nme859>3.0.co;2-z
Subject(s) - multibody system , stability (learning theory) , newmark beta method , constraint (computer aided design) , acceleration , spurious relationship , energy (signal processing) , equations of motion , kinematics , dissipative system , control theory (sociology) , direct integration of a beam , computer science , mathematics , mathematical optimization , finite element method , engineering , structural engineering , classical mechanics , physics , statistics , geometry , control (management) , quantum mechanics , machine learning , artificial intelligence
This paper focuses on the development of an unconditionally stable time‐integration algorithm for multibody dynamics that does not artificially dissipate energy. Unconditional stability is sought to alleviate any stability restrictions on the integration step size, while energy conservation is important for the accuracy of long‐term simulations. In multibody system analysis, the time‐integration scheme is complemented by a choice of co‐ordinates that define the kinematics of the system. As such, the current approach uses a non‐dissipative implicit Newmark method to integrate the equations of motion defined in terms of the independent joint co‐ordinates of the system. In order to extend the unconditional stability of the implicit Newmark method to non‐linear dynamic systems, a discrete energy balance is enforced. This constraint, however, yields spurious oscillations in the computed accelerations and therefore, a new acceleration corrector is developed to eliminate these instabilities and hence retain unconditional stability in an energy sense. An additional benefit of employing the non‐linearly implicit time‐integration method is that it allows for an efficient design sensitivity analysis. In this paper, design sensitivities computed via the direct differentiation method are used for mechanism performance optimization. Copyright © 2000 John Wiley & Sons, Ltd.