Abstract
This paper presents a new method for the integration of the equations of motion of constrained multibody systems in descriptor form. The method is based on the penalty-Augmented Lagrangian formulation and uses massorthogonal projections for the solution to satisfy the kinematic constraint conditions. The number of equations being solved is equal to the number of states, and does not depend on the number of constraint conditions. Therefore, the method is particularly suitable for systems with redundant constraints, singular configurations or topology changes. The major advantage of the new method relies on the fact that for a low computational cost, the constraints in positions, velocities and accelerations are satisfied to machine precision during the numerical integration. This process is efficiently done by means of a mass-orthogonal projection without the need for coordinate partitioning or reduction to a minimum set of coordinates. The projection scheme allows for a more accurate and robust integration of the equations of motion since constraint violations constitute one of the primary sources of numerical errors and instabilities during the integration process. The proposed projection is also applied to the classical Lagrangian approach, thus eliminating the need for further stabilization as well as the selection of parameters in Baumgarte's method.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Baumgarte, J., ‘Stabilization of constraints and integrals of motion in dynamical systems’,Computer Methods in Applied Mechanics and Engineering 1, 1972, 1–16.
Park, K. C. and Chiou, J. C., ‘Stabilization of computational procedures for constrained dynamical systems’,Journal of Guidance, Control, and Dynamics 11, 1988, 365–370.
Bayo, E., García de Jalón, J., and Serna, M. A., ‘A modified Lagrangian formulation for the dynamic analysis of constrained mechanical systems’,Computer Methods in Applied Mechanics and Engineering 71, 1988, 183–195.
Kurdila, A. J., Junkins, J. L., and Hsu, S., ‘Lyapunov stable penalty methods for imposing nonholonomic constraints in multibody system dynamics’,Nonlinear Dynamics (to appear).
Bayo, E. and Avello, A., ‘Singularity free augmented Lagrangian algorithms for constraint multibody dynamics’,Nonlinear Dynamics 5, 1994, 209–231.
Wehage, R. A. and Haug, E. J., ‘Generalized coordinate partitioning for dimension reduction in analysis of constrained dynamic systems’,Journal of Mechanical Design 104, 1982, 247–255.
Kane, T. R. and Levinson, D. A.,Dynamics: Theory and Applications, McGraw-Hill, New York, 1985.
Huston, R. L.,Multibody Dynamics, Butterworth-Heinemann, New York, 1990.
Serna, M. A., Avilés, R., and García de Jalón, J., ‘Dynamic analysis of planar mechanisms with lower-pairs in basic coordinates’,Mechanism and Machine Theory 17, 1982, 397–403.
Jerkovsky, W., ‘The structure of multibody dynamic equations’,Journal of Guidance and Control 1, 1978, 173–182.
Kim, S. S. and Vanderploeg, M. J., ‘A general and efficient method for dynamic analysis of mechanical systems using velocity transformations’,Journal of Mechanisms, Transmissions and Automation in Design 108, 1986, 176–182.
Nikravesh, P. E. and Gim, G., ‘Systematic construction of the equations of motion for multibody systems containing closed kinematic loops’,Advances in Design Automation 3, 1989, 27–33.
García de Jalón, J., Avello, A., Jiménez, J. M., Martín, F., and Cuadrado, J., ‘Real-time simulation of complex 3-D multibody systems with realistic graphics’, InReal-Time Integration Methods for Mechanical System Simulation, NATO ASI Series, Vol. 69, Springer-Verlag, Heidelberg, 1990, pp. 265–292.
Bae, D.-S. and Won, Y. S., ‘A Hamiltonian equation of motion for real time vehicle simulation’,Advances in Design Automation 2, 1990, 151–157.
Avello, A., Jiménez, J. M., Bayo, E., and García de Jalón, J., ‘A simple and highly parallelizable method for real-time dynamic simulation based on velocity transformations’,Computer Methods in Applied Mechanics and Engineering 107, 1993, 313–339.
Brenan, K. E., Campbell, S. L., and Petzold, L. R.,The Numerical Solution of Initial Value Problems in Differential-Algebraic Equations, Elsevier, New York, 1989.
Griepentrog, E., Hanke, M., and Marz, R., ‘Berlin seminar on algebraic equations’, Fachbereich Mathematik der Humboldt-Universität z Berlin, 1992.
Fuhrer, C. and Schwertassek, R., ‘Generation and solution of multibody systems equations’,International Journal of Non-Linear Mechanics 25, 1990, 127–141.
Fuhrer, C. and Leimkuhler, B. J., ‘Numerical solution of differential algebraic equations for constrained mechanical motion’,Numerische Mathematik 59, 1991, 55–69.
Matsson, S. E. and Soderlind, G., ‘Index reduction in differential-algebraic equations using dummy derivatives’,SIAM Journal of Scientific Computation 14, 1993, 667–692.
Shampine, L. F. and Gordon, M. K.,Computer Solution of Ordinary Differential Equations: The Initial Value Problem, W. J. Freeman, San Francisco, CA, 1975.
García de Jalón, J. and Bayo, E.,Kinematic and Dynamic Simulation of Multibody Systems, Springer-Verlag, New York, 1994.
Sanz-Serna, J. M., ‘Symplectic integrators for Hamiltonian problems: An overview’, inActa Numerica 1992, Cambridge University Press, Cambridge, England, 1992, pp. 243–286.
Schiehlen, W. O.,Multibody System Handbook, Springer-Verlag, Heidelberg, 1990.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Bayo, E., Ledesma, R. Augmented lagrangian and mass-orthogonal projection methods for constrained multibody dynamics. Nonlinear Dyn 9, 113–130 (1996). https://doi.org/10.1007/BF01833296
Received:
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/BF01833296