Research Papers

Determination of Holonomic and Nonholonomic Constraint Reactions in an Index-3 Augmented Lagrangian Formulation With Velocity and Acceleration Projections

[+] Author and Article Information
Daniel Dopico

Department of Computer Science and
Department of Mechanical Engineering,
Virginia Tech,
2200 Kraft Drive,
Blacksburg, VA 24060
e-mail: ddopico@udc.es

Francisco González

Department of Mechanical Engineering
and Centre for Intelligent Machines,
McGill University,
817 Sherbrooke Street West,
Montréal, QC H3A 0C3, Canada
e-mail: franglez@cim.mcgill.ca

Javier Cuadrado

Laboratorio de Ingeniería Mecánica,
University of La Coruña,
Mendizábal s/n,
Ferrol 15403, Spain
e-mail: javicuad@cdf.udc.es

József Kövecses

Department of Mechanical Engineering and Centre for Intelligent Machines,
McGill University,
817 Sherbrooke Street West,
Montréal, QC H3A 0C3, Canada
e-mail: jozsef.kovecses@mcgill.ca

1Corresponding author.

Manuscript received July 20, 2013; final manuscript received May 9, 2014; published online July 11, 2014. Assoc. Editor: Dan Negrut.

J. Comput. Nonlinear Dynam 9(4), 041006 (Jul 11, 2014) (9 pages) Paper No: CND-13-1185; doi: 10.1115/1.4027671 History: Received July 20, 2013; Revised May 09, 2014

Index-3 augmented Lagrangian formulations with projections of velocities and accelerations represent an efficient and robust method to carry out the forward-dynamics simulation of multibody systems modeled in dependent coordinates. Existing formalisms, however, were only established for holonomic systems, for which the expression of the constraints at the position-level is known. In this work, an extension of the original algorithms for nonholonomic systems is introduced. Moreover, projections of velocities and accelerations have two side effects: they modify the kinetic energy of the system and they contribute to the constraint reaction forces. Although the effects of the projections on the energy have been studied by several authors, their role in the calculation of the reaction forces has not been described so far. In this work, expressions to determine the constraint reactions from the Lagrange multipliers of the dynamic equations and the Lagrange multipliers of the velocity and acceleration projections are introduced. Simulation results show that the proposed strategy can be used to expand the capabilities of index-3 augmented Lagrangian algorithms, making them able to deal with nonholonomic constraints and provide correct reaction efforts.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Haug, E. J., 1989, Computer Aided Kinematics and Dynamics of Mechanical Systems, Allyn and Bacon, Inc., Boston, MA.
García de Jalón, J., and Bayo, E., 1994, Kinematic and Dynamic Simulation of Multibody Systems. The Real-Time Challenge, Springer, New York.
Paul, B., and Krajcinovic, D., 1970, “Computer Analysis of Machines With Planar Motion,” ASME J. Appl. Mech., 37(3), pp. 697–712. [CrossRef]
Sheth, P., and Uicker, J., 1972, “IMP (Integrated Mechanism Program): A Computer-Aided Design Analysis System for Mechanisms and Linkages,” ASME J. Eng. Ind., 94(2), pp. 454–464. [CrossRef]
Smith, D., Chace, M., and Rubens, A., 1973, “The Automatic Generation of a Mathematical Model for Machinery Systems,” ASME J. Eng. Ind., 95(2), pp. 629–635. [CrossRef]
Brenan, K., Campbell, S., and Petzold, L., 1989, Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations, North-Holland, New York.
Baumgarte, J., 1972, “Stabilization of Constraints and Integrals of Motion in Dynamical Systems,” Comput. Methods Appl. Mech. Eng., 1(1), pp. 1–16. [CrossRef]
Serna, M. A., Avilés, R., and García de Jalón, J., 1982, “Dynamic Analysis of Plane Mechanisms With Lower Pairs in Basic Coordinates,” Mech. Mach. Theor., 17(6), pp. 397–403. [CrossRef]
Wehage, R., and Haug, E., 1982, “Generalized Coordinate Partitioning for Dimension Reduction in Analysis of Constrained Mechanical Systems,” ASME J. Mech. Des., 104, pp. 247–255. [CrossRef]
Bayo, E., García de Jalón, J., and Serna, M. A., 1988, “A Modified Lagrangian Formulation for the Dynamic Analysis of Constrained Mechanical Systems,” Comput. Methods Appl. Mech. Eng., 71(2), pp. 183–195. [CrossRef]
Bayo, E., and Ledesma, R., 1996, “Augmented Lagrangian and Mass-Orthogonal Projection Methods for Constrained Multibody Dynamics,” Nonlinear Dyn., 9(1–2), pp. 113–130. [CrossRef]
Cuadrado, J., Cardenal, J., Morer, P., and Bayo, E., 2000, “Intelligent Simulation of Multibody Dynamics: Space-State and Descriptor Methods in Sequential and Parallel Computing Environments,” Multibody Syst. Dyn., 4(1), pp. 55–73. [CrossRef]
Cuadrado, J., Gutierrez, R., Naya, M., and Morer, P., 2001, “A Comparison in Terms of Accuracy and Efficiency Between a MBS Dynamic Formulation With Stress Analysis and a Non-Linear FEA Code,” Int. J. Numer. Methods Eng., 51(9), pp. 1033–1052. [CrossRef]
García Orden, J., and Goicolea, J., 2000, “Conserving Properties in Constrained Dynamics of Flexible Multibody Systems,” Multibody Syst. Dyn., 4(2–3), pp. 225–244. [CrossRef]
Goicolea, J., and García Orden, J., 2000, “Dynamic Analysis of Rigid and Deformable Multibody Systems With Penalty Methods and Energy-Momentum Schemes,” Comput. Methods Appl. Mech. Eng., 188(4), pp. 789–804. [CrossRef]
Goicolea, J., and García Orden, J., 2002, “Quadratic and Higher-Order Constraints in Energy-Conserving Formulations of Flexible Multibody Systems,” Multibody Syst. Dyn., 7(1), pp. 3–29. [CrossRef]
García Orden, J., and Ortega, R., 2006, “A Conservative Augmented Lagrangian Algorithm for the Dynamics of Constrained Mechanical Systems,” Mech. Based Des. Struct. Mach., 34(4), pp. 449–468. [CrossRef]
García Orden, J. C., and Dopico, D., 2007, “Multibody Dynamics. Computational Methods and Applications,” On the Stabilizing Properties of Energy-Momentum Integrators and Coordinate Projections for Constrained Mechanical Systems, Springer, Dordrecht, The Netherlands, pp. 49–67.
García Orden, J., 2010, “Energy Considerations for the Stabilization of Constrained Mechanical Systems With Velocity Projection,” Nonlinear Dyn., 60(1–2), pp. 49–62. [CrossRef]
García Orden, J., and Conde, S., 2012, “Controllable Velocity Projection for Constraint Stabilization in Multibody Dynamics,” Nonlinear Dyn., 68(1–2), pp. 245–257. [CrossRef]
Cuadrado, J., Dopico, D., Naya, M. A., and González, M., 2004, “Penalty, Semi-Recursive and Hybrid Methods for MBS Real-Time Dynamics in the Context of Structural Integrators,” Multibody Syst. Dyn., 12(2), pp. 117–132. [CrossRef]
Géradin, M., and Cardona, A., 2001, Flexible Multibody Dynamics—A Finite Element Approach, 2nd ed., Wiley, Chichester, England.
Newmark, N. M., 1959, “A Method of Computation for Structural Dynamics,” ASCE J. Eng. Mech. Div., 85(EM3), pp. 67–94.


Grahic Jump Location
Fig. 3

3D motion of a disk rolling on a plane

Grahic Jump Location
Fig. 2

Reaction forces acting on mass A in the x and y directions due to the nonholonomic constraint in the first example, computed with direct solution (54) and with the ALI3-P formulation

Grahic Jump Location

Reaction fy at contact point C during 2D motion of the disk, obtained with the projection of velocities (fy vel) and accelerations (fy acc)

Grahic Jump Location
Fig. 5

Reaction forces fx and fy at contact point C during motion of the disk, compared to mx··P and my··P

Grahic Jump Location
Fig. 1

Nonholonomic constraint restraining the velocity of mass A to be collinear with segment A–B

Grahic Jump Location
Fig. 7

Actuation forces: (a) Bayo's (top) and (b) evolved Cuadrado's (bottom)



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In