Investigation on the Baumgarte Stabilization Method for Dynamic Analysis of Constrained Multibody Systems

2008 ◽  
pp. 305-312 ◽  
Author(s):  
Paulo Flores ◽  
Rui Pereira ◽  
Margarida Machado ◽  
Eurico Seabra
Author(s):  
Paulo Flores ◽  
Margarida Machado ◽  
Eurico Seabra ◽  
Miguel Tavares da Silva

This paper presents and discusses the results obtained from a parametric study on the Baumgarte stabilization method for forward dynamics of constrained multibody systems. The main purpose of this work is to analyze the influence of the variables that affect the violation of constraints, chiefly the values of the Baumgarte parameters, the integration method, the time step, and the quality of the initial conditions for the positions. In the sequel of this process, the formulation of the rigid multibody systems is reviewed. The generalized Cartesian coordinates are selected as the variables to describe the bodies’ degrees of freedom. The formulation of the equations of motion uses the Newton–Euler approach, augmented with the constraint equations that lead to a set of differential algebraic equations. Furthermore, the main issues related to the stabilization of the violation of constraints based on the Baumgarte approach are revised. Special attention is also given to some techniques that help in the selection process of the values of the Baumgarte parameters, namely, those based on the Taylor’s series and the Laplace transform technique. Finally, a slider-crank mechanism with eccentricity is considered as an example of application in order to illustrate how the violation of constraints can be affected by different factors.


Author(s):  
Paulo Flores ◽  
Margarida Machado ◽  
Eurico Seabra ◽  
Miguel Tavares da Silva

This paper presents and discusses the results obtained from a parametric study on the Baumgarte stabilization method for forward dynamics of constrained multibody systems. The main purpose of this work is to analyze the influence of the variables that affect the violation of constraints, chiefly the values of the Baumgarte parameters, the integration method, the time step and the quality of the initial conditions for the positions. In the sequel of this process the formulation of the rigid multibody systems is reviewed. The generalized Cartesian coordinates are selected as the variables to describe the bodies’ degrees of freedom. The formulation of the equations of motion uses the Newton-Euler approach that is augmented with the constraint equations that lead to a set of differential algebraic equations. Furthermore, the main issues related to the stabilization of the violation of constraints based on the Baumgarte approach are revised. Special attention is also given to some techniques that help in the selection process of the values of the Baumgarte parameters, namely those based on the Taylor’s series and Laplace transform technique. Finally, a slider crank mechanism with eccentricity is considered as an example of application in order to illustrated how the violation of constraints can be affected by different factors such as the Baumgarte parameters, integrator, time step and initial guesses.


Author(s):  
Apiwat Reungwetwattana ◽  
Shigeki Toyama

Abstract This paper presents an efficient extension of Rosenthal’s order-n algorithm for multibody systems containing closed loops. Closed topological loops are handled by cut joint technique. Violation of the kinematic constraint equations of cut joints is corrected by Baumgarte’s constraint violation stabilization method. A reliable approach for selecting the parameters used in the constraint stabilization method is proposed. Dynamic analysis of a slider crank mechanism is carried out to demonstrate efficiency of the proposed method.


Author(s):  
E. Bayo ◽  
J. M. Jimenez

Abstract We investigate in this paper the different approaches that can be derived from the use of the Hamiltonian or canonical equations of motion for constrained mechanical systems with the intention of responding to the question of whether the use of these equations leads to more efficient and stable numerical algorithms than those coming from acceleration based formalisms. In this process, we propose a new penalty based canonical description of the equations of motion of constrained mechanical systems. This technique leads to a reduced set of first order ordinary differential equations in terms of the canonical variables with no Lagrange’s multipliers involved in the equations. This method shows a clear advantage over the previously proposed acceleration based formulation, in terms of numerical efficiency. In addition, we examine the use of the canonical equations based on independent coordinates, and conclude that in this second case the use of the acceleration based formulation is more advantageous than the canonical counterpart.


Author(s):  
Andreas Müller ◽  
Zdravko Terze

The stabilization of geometric constraints is vital for an accurate numerical solution of the differential-algebraic equations (DAE) governing the dynamics of constrained multibody systems (MBS). Although this has been a central topic in numerical MBS dynamics using classical vector space formulations, it has not yet been sufficiently addressed when using Lie group formulations. A straightforward approach is to impose constraints directly on the Lie group elements that represent the MBS motion, which requires additional constraints accounting for the invariants of the Lie group. On the other hand, most numerical Lie group integration schemes introduce local coordinates within the integration step, and it is natural to perform the stabilization in terms of these local coordinates. Such a formulation is presented in this paper for index 1 formulation. The stabilization method is applicable to general coordinate mappings (canonical coordinates, Cayley-Rodriguez, Study) on the MBS configuration space Lie group. The stabilization scheme resembles the well-known vectors space projection and pseudo-inverse method consisting in an iterative procedure. A numerical example is presented and it is shown that the Lie group stabilization scheme converges normally within one iteration step, like the scheme in the vector space formulation.


1994 ◽  
Vol 116 (2) ◽  
pp. 423-428 ◽  
Author(s):  
W. Blajer ◽  
D. Bestle ◽  
W. Schiehlen

A method is proposed for the automatic generation of an orthogonal complement matrix to the constraint matrix for the dynamic analysis of constrained multibody systems. The clue for this method lies in the determination of local constraint matrices and their orthogonal complements relative to the local reference frames of particular constrained points. These matrices are then transformed into the system’s configuration space in order to form the final constraint matrix and its orthogonal complement. The avoidance of singularities in the formulation is discussed. The method is especially suited for the dynamic analysis of multibody systems with many constraints and/or closed-loops.


Sign in / Sign up

Export Citation Format

Share Document