Convergence of the generalized-alpha scheme for constrained mechanical systems

A variant of the generalized-alpha scheme is proposed for constrained mechanical systems represented by index-3 DAEs. Based on the analogy with linear multistep methods, an elegant convergence analysis is developed for this algorithm. Second-order convergence is demonstrated both for the generalized coordinates and the Lagrange multipliers, and those theoretical results are illustrated by numerical tests.