Time Integration Algorithm for Flexible Mechanism Dynamics: The DAE alpha-Method

摘要

This paper introduces a new family of second-order methods for solving the index-2 DAE equations of motion for flexible mechanism dynamics. These methods, which extend the alpha-methods for ODEs of structural dynamics to DAEs, possess numerical dissipation that can be controlled by the user. Convergence and stability analysis is given and verifies that the DAE alpha-methods introduce no additional oscillations and preserve the stability of the underlying ODE system. Convergence of the Newton iteration, which can be a source of difficulty in solving nonlinear oscillartoy systems with large stepsizes, is achieved via a coordinate-split modification to the Newton iterations. Numerical results illustrate the effectiveness of the new methods for simulation of flexible mechanisms.