Application of the Bezier integration technique with enhanced stability in forward dynamics of constrained multibody systems with Baumgarte stabilization method

This paper deals with the study and application of a numerical integration scheme based on the Bezier curves to obtain stable and accurate solutions of the governing differential-algebraic equations (DAEs) of constrained multibody systems. For this purpose, the standard Lagrange multiplier method is utilized to derive the equations of motion for constrained mechanical systems. It is well known that constraints violation and instability are amongst the main computational difficulties of various numerical integration algorithms to provide accurate solutions for DAEs. In the present study, the Baumgarte stabilization technique is employed, as well as a criterion to select the Baumgarte parameters based on the stability domain analysis. Stability regions based on dimensionless Baumgarte parameters are obtained for different orders of Bezier curves together with other classic numerical integration algorithms, namely the Adams-Bashforth and Runge-Kutta techniques. It is shown that the proper identification and selection of Baumgarte parameters is highly dependent on the integration scheme utilized. A comparative analysis performed within the present work reveals that the Bezier method provides substantially wider region of stability for the same level of computational effort. A planar slider-crank mechanism is considered as an example of application to demonstrate and implement the Bezier technique, which permits to examine the benefits of the presented stability analysis. It is demonstrated that for certain values of Baumgarte parameters, the Bezier method remains stable, while other methods become unstable. Due to the good performance of the applied method based on the Bezier approach, both in terms of efficiency and stability, it is expected that the introduced integration method to be used for various benchmark studies performed under the umbrella of multibody dynamics.