DOCSLIB.ORG

Dynamic Modeling and Simulation of Multi-Body Systems Using the Udwadia-Kalaba Theory

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Dynamic Modeling and Simulation of Multi-Body Systems Using the Udwadia-Kalaba Theory CHINESE JOURNAL OF MECHANICAL ENGINEERING Vol. 26, No. 5, 2013 ·839· DOI: 10.3901/CJME.2013.05.839, available online at www.springerlink.com; www.cjmenet.com; www.cjmenet.com.cn Dynamic Modeling and Simulation of Multi-body Systems Using the Udwadia-Kalaba Theory ZHAO Han1, ZHEN Shengchao1, 2, *, and CHEN Ye-Hwa2 1 School of Mechanical Engineering, Hefei University of Technology, Hefei 230009, China 2 School of Mechanical Engineering, Georgia Institute of Technology, Atlanta 30332, USA Received December 20, 2012; revised March 20, 2013; accepted March 27, 2013 Abstract: Laboratory experiments were conducted for falling U-chain, but explicit analytic form of the general equations of motion was not presented. Several modeling methods were developed for fish robots, however they just focused on the whole fish’s locomotion which does little favor to understand the detailed swimming behavior of fish. Udwadia-Kalaba theory is used to model these two multi-body systems and obtain explicit analytic equations of motion. For falling U-chain, the mass matrix is non-singular. Second-order constraints are used to get the constraint force and equations of motion and the numerical simulation is conducted. Simulation results show that the chain tip falls faster than the freely falling body. For fish robot, two-joint Carangiform fish robot is focused on. Quasi-steady wing theory is used to approximately calculate fluid lift force acting on the caudal fin. Based on the obtained explicit analytic equations of motion (the mass matrix is singular), propulsive characteristics of each part of the fish robot are obtained. Through these two cases of U chain and fish robot, how to use Udwadia-Kalaba equation to obtain the dynamical model is shown and the modeling methodology for multi-body systems is presented. It is also shown that Udwadia-Kalaba theory is applicable to systems whether or not their mass matrices are singular. In the whole process of applying Udwadia-Kalaba equation, Lagrangian multipliers and quasi-coordinates are not used. Udwadia-Kalaba theory is creatively applied to dynamical modeling of falling U-chain and fish robot problems and explicit analytic equations of motion are obtained. Key words: Udwadia-Kalaba equation, multi-body systems, falling U-chain, fish robot developed equations of constrained motion when the 1 Introduction∗ constraints satisfy D’ Alember’ s principle. In the treatise [5] on the analytical dynamics, PARS refers to the The general problem of obtaining equations of motion Gibbs-Appell equations as “probably the most for constrained discrete mechanical systems has been an comprehensive equations of motion so far discovered”. But area of considerable interest among scientists and engineers. the Gibbs-Appell equations require a “lucky” choice of It is also one of the central issues in multi-body dynamics. problem-specific quasi-coordinates and they suffer from The problem has been aggressively and continuously similar problems when dealing with systems with a large pursued by many scientists, engineers and mathematicians number of degrees of freedom and many non- integrable [6] since constrained motion was initially described by constraints. DIRAC used Poisson brackets, a recursive LAGRANGE[1]. He invented the special Lagrange scheme for determining the Lagrange multipliers, for multiplier method to deal with constrained motion. singular Hamiltonian systems where the constraints do not However, the Lagrange multiplier method relies on exactly depend on time. [7–8] problem-specific approaches to determine the multipliers; it On the other hand, UDWADIA, et al , obtained a is often very difficult to find the multipliers to obtain the concise, explicit set of equations of motion for constrained explicit equations of motion for systems with large discrete dynamic systems which lead to a simple and new numbers of degrees of freedom and a mass of fundamental view of Lagrangian mechanics. They derived non-integrable constraints. GAUSS[2] introduced a new the fundamental equation of motion that describes the general principle of mechanics for handling constrained dynamics of constrained systems from Gauss’ s principle motion. Gauss’s Principle gives a clear description of the which seems somewhat less popular than the principles of general nature of constrained motion through minimization Lagrange, Hamilton, Gibbs and Appell. The equations can of a function of the accelerations of the particles of a deal with holonomic and also non-holonomic constraints. [9–10] system. GIBBS[3] and APPELL[4] have independently UDWADIA, et al , observed that all the research above has used D’ Alember’ s principle as their starting point. D’ Alember’ s principle assumes that the forces of constraints * Corresponding author. E-mail: [email protected] are considered to be ideal and the total work done by the © Chinese Mechanical Engineering Society and Springer-Verlag Berlin Heidelberg 2013 ZHAO Han, et al: Dynamic Modeling and Simulation of Multi-body Systems ·840· Using the Udwadia-Kalaba Theory forces of constraints under virtual displacement is always Boltzmann, we are not required to develop a system of zero. This assumption works well in many situations and is proper quasi-coordinates. regarded as the core of classical analytical dynamics, but it The significant advantages of Udwadia and Kalaba’ s is not applicable when the constraints are nonideal. Thus, new theory are threefold. First, it allows descriptions of UDWADIA and KALABA generalized their previous mechanical systems that use more than the minimum equations to constrained mechanical systems that may not number of required coordinates which speeds and eases the satisfy D’ Alember’ s principle. Systems with singular setup of the equations of motion of complex systems. mass matrices are not common in classical dynamics when Second, it applies to systems whether or not their mass dealing with unconstrained motion. PARS[5] proposed that matrices are singular, is applicable to systems with when the minimum number of coordinates is employed for holonomic and/or nonholonomic constraints and to systems describing the unconstrained motion of mechanical systems, whose constraint forces may or may not be ideal. Third, the the corresponding Lagrange equations usually yield new theory opens up a new way of modeling complex non-singular, symmetric and positive definite mass multi-body systems. It permits decomposition of such matrices. However, singular mass matrices can arise when systems into sub-systems whose equations of motion are one wants greater flexibility in modeling complex known, and then combine the sub-system equations to get mechanical systems by using more than the minimum the composite system’ s equations of motion in a number of required generalized coordinates. Thus, straightforward and simple manner. UDWADIA, et al[11], developed general and explicit equations of motion to handle systems whether or not their 2 Explicit General Equations of Motion mass matrices are singular. We call these discoveries the for Multi-body System Udwadia-Kalaba theory. About 20 years has passed, we seldom realize its importance only because it is abstract. In Using Udwadia and Kalaba’ s approach, we first fact, the Udwadia-Kalaba theory can be applied to consider an unconstrained discrete dynamical system constrained discrete mechanical systems with unmatched whose configuration is described by the n generalized ease, clarity and elegance, especially to multi-body systems T coordinates q : [q12 ,q , ,qn ]. Its equation of motion such as mechanical manipulators, space vehicles, robots can be obtained, using Newtonian or Lagrangian mechanics, and similar systems. by the relation: As we know, the principles of mechanics[12–16] are so perfect that it is impossible to create a totally new Mq(, tt ) q Qq (, q , ), (1) fundamental principal for the theory of motion and nn equilibrium of discrete, dynamical systems. However, an where Mq( , t ) R is symmetric and positive definite n n additional perspective has been proposed by Udwadia and inertia matrix, q R is the velocity and q R is the n Kalaba which is useful to help us understand nature’ s law acceleration, and Qq(, q , t ) R is the force imposed on from new points of view[7–11]. the system whose constraints are released. The imposed Using Udwadia and Kalaba’ s marvelous theory, we can forces could include centrifugal force, gravitational force formulate the explicit, general equations of motion for and control input. The generalized acceleration of the constrained discrete dynamic systems in three steps：First, unconstrained system, which we denote by aq(, q , t ) is in terms of the generalized coordinates, we consider the thus given by unconstrained discrete dynamic system whose equations of 1 motion can be written using Newtonian or Lagrangian q M(, qtt ) Qqq (, , ) aqq (, , t ). (2) mechanics. Then we form the constraint equations which Second, constraints present in the system should be could include the usual holonomic, non-holonomic, considered. We shall assume that the system is subjected to scleronomic, rheonomic, catastatic and acataetatic variaties h holonomic constraints and mh nonholonomic of constraints or even combinations of such constraints. constraints of the form: Finally, we impose the additional generalized forces of constraint on the system. That is to say, we add the force ϕi (q , ti ) 0, 1, 2, , h , (3) resulting from the presence of constraints to the unconstrained system’ s force. ϕi (qq , , t ) 0,
Load more

Recommended publications
  • Chapter 3 Dynamics of Rigid Body Systems
    Rigid Body Dynamics Algorithms Roy Featherstone Rigid Body Dynamics Algorithms Roy Featherstone The Austrailian National University Canberra, ACT Austrailia Library of Congress Control Number: 2007936980 ISBN 978-0-387-74314-1 ISBN 978-1-4899-7560-7 (eBook) Printed on acid-free paper. @ 2008 Springer Science+Business Media, LLC All rights reserved. This work may not be translated or copied in whole or in part without the written permission of the publisher (Springer Science+Business Media, LLC, 233 Spring Street, New York, NY 10013, USA), except for brief excerpts in connection with reviews or scholarly analysis. Use in connection with any form of information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed is forbidden. The use in this publication of trade names, trademarks, service marks and similar terms, even if they are not identified as such, is not to be taken as an expression of opinion as to whether or not they are subject to proprietary rights. 9 8 7 6 5 4 3 2 1 springer.com Preface The purpose of this book is to present a substantial collection of the most efficient algorithms for calculating rigid-body dynamics, and to explain them in enough detail that the reader can understand how they work, and how to adapt them (or create new algorithms) to suit the reader’s needs. The collection includes the following well-known algorithms: the recursive Newton-Euler algo- rithm, the composite-rigid-body algorithm and the articulated-body algorithm. It also includes algorithms for kinematic loops and floating bases.
    [Show full text]
  • Lecture Notes for PC2132 (Classical Mechanics)
    Version: 27th Sept, 2017 11:27; svn-64 Lecture notes for PC2132 (Classical Mechanics) For AY2017/18, Semester 1 Lecturer: Christian Kurtsiefer, Tutors: Adrian Nugraha Utama, Shi Yicheng, Jiaan Qi, Tan Zong Sheng Disclaimer These notes are by no means a suitable replacement for a proper textbook on classical mechanics, nor should it replace your own notes. It is just a best-effort affair with the intention to be useful. This is a document in process – it hopefully does not contain too many mis- takes, but please contact us if you feel that you spotted one. 1 Version: 27th Sept, 2017 11:27; svn-64 Notation There is an attempt to do use consistent notations through this lecture. Below is a list what symbols typically refer to, unless they are referenced to otherwise. a a vector; often the acceleration e a unit vector, i.e., a vector of length 1 (e e = 1) dA a vector differential representing an oriented· surface area element ds a vector differential representing a line element F a vector representing a force H a scalar representing the Hamilton function of a system I the tensor for the inertia of a rigid body L a scalar representing the Lagrange function L a vector representing an angular momentum N a vector representing a torque m a scalar representing the mass of a particle m mass matrix in a system of coupled masses M the total mass of a system p a vector representing the momentum mv of a particle P a vector representing the total momentum of a system qk a generalized coordinate Q the “quality factor” of a damped harmonic oscillator r (or
    [Show full text]
  • Chapter 2 Lagrange's and Hamilton's Equations
    Chapter 2 Lagrange’s and Hamilton’s Equations In this chapter, we consider two reformulations of Newtonian mechanics, the Lagrangian and the Hamiltonian formalism. The first is naturally associated with configuration space, extended by time, while the latter is the natural description for working in phase space. Lagrange developed his approach in 1764 in a study of the libration of the moon, but it is best thought of as a general method of treating dynamics in terms of generalized coordinates for configuration space. It so transcends its origin that the Lagrangian is considered the fundamental object which describes a quantum field theory. Hamilton’s approach arose in 1835 in his unification of the language of optics and mechanics. It too had a usefulness far beyond its origin, and the Hamiltonian is now most familiar as the operator in quantum mechanics which determines the evolution in time of the wave function. We begin by deriving Lagrange’s equation as a simple change of coordi- nates in an unconstrained system, one which is evolving according to New- ton’s laws with force laws given by some potential. Lagrangian mechanics is also and especially useful in the presence of constraints, so we will then extend the formalism to this more general situation. 35 36 CHAPTER 2. LAGRANGE’S AND HAMILTON’S EQUATIONS 2.1 Lagrangian for unconstrained systems For a collection of particles with conservative forces described by a potential, we have in inertial cartesian coordinates mx¨i = Fi. The left hand side of this equation is determined by the kinetic energy func- tion as the time derivative of the momentum pi = ∂T/∂x˙ i, while the right hand side is a derivative of the potential energy, ∂U/∂x .AsT is indepen- − i dent of xi and U is independent ofx ˙ i in these coordinates, we can write both sides in terms of the Lagrangian L = T U, which is then a function of both the coordinates and their velocities.
    [Show full text]
  • DUALITY SYSTEM in APPLIED MECHANICS and OPTIMAL CONTROL Advances in Mechanics and Mathematics
    DUALITY SYSTEM IN APPLIED MECHANICS AND OPTIMAL CONTROL Advances in Mechanics and Mathematics Volume 5 Series Editors: David Y. Gao Virginia Polytechnic Institute and State University, U.S.A. Ray W. Ogden University of Glasgow, U.K. Advisory Editors: I. Ekeland University of British Columbia, Canada K.R. Rajagopal Texas A&M University, U.S.A. W. Yang Tsinghua University, P.R. China DUALITY SYSTEM IN APPLIED MECHANICS AND OPTIMAL CONTROL by Wan-Xie Zhong Dalian University of Technology KLUWER ACADEMIC PUBLISHERS NEW YORK, BOSTON, DORDRECHT, LONDON, MOSCOW eBook ISBN: 1-4020-7881-1 Print ISBN: 1-4020-7880-3 ©2004 Kluwer Academic Publishers New York, Boston, Dordrecht, London, Moscow Print ©2004 Kluwer Academic Publishers Boston All rights reserved No part of this eBook may be reproduced or transmitted in any form or by any means, electronic, mechanical, recording, or otherwise, without written consent from the Publisher Created in the United States of America Visit Kluwer Online at: http://kluweronline.com and Kluwer's eBookstore at: http://ebooks.kluweronline.com Contents v Preface xi Introduction 1 0.1, Introduction to Precise Integration method 4 0.1.1, Homogeneous equation, algorithm for exponential matrix. 4 0.1.2, Solution of inhomogeneous equation 6 0.1.3, Precision analysis 7 0.1.4, Discussions on time-variant system or non-linear system 9 Chapter 1, Introduction to analytical dynamics 11 1.1, Holonomic and nonholonomic constraints 11 1.2, Generalized displacement, degrees of freedom and virtual displacement 14 1.3, Principle of
    [Show full text]
  • Solid Mechanics at Harvard University
    SOLID MECHANICS James R. Rice School of Engineering and Applied Sciences, and Department of Earth and Planetary Sciences Harvard University, Cambridge, MA 02138 USA Original version: October 1994 This revision: February 2010 Downloadable at: http://esag.harvard.edu/rice/e0_Solid_Mechanics_94_10.pdf TABLE OF CONTENTS provided on last three pages, pp. 87-89 INTRODUCTION The application of the principles of mechanics to bulk matter is conventionally divided into the mechanics of fluids and the mechanics of solids. The entire subject is often called continuum mechanics, particularly when we adopt the useful model of matter as being continuously divisible, making no reference to its discrete structure at microscopic length scales well below those of the application or phenomenon of interest. Solid mechanics is concerned with the stressing, deformation and failure of solid materials and structures. What, then, is a solid? Any material, fluid or solid, can support normal forces. These are forces directed perpendicular, or normal, to a material plane across which they act. The force per unit of area of that plane is called the normal stress. Water at the base of a pond, air in an automobile tire, the stones of a Roman arch, rocks at base of a mountain, the skin of a pressurized airplane cabin, a stretched rubber band and the bones of a runner all support force in that way (some only when the force is compressive). We call a material solid rather than fluid if it can also support a substantial shearing force over the time scale of some natural process or technological application of interest.
    [Show full text]
  • Singular Mass Matrix and Redundant Constraints in Unilaterally Constrained Lagrangian and Hamiltonian Systems Bernard Brogliato, Daniel Goeleven
    Singular mass matrix and redundant constraints in unilaterally constrained Lagrangian and Hamiltonian systems Bernard Brogliato, Daniel Goeleven To cite this version: Bernard Brogliato, Daniel Goeleven. Singular mass matrix and redundant constraints in unilaterally constrained Lagrangian and Hamiltonian systems. Multibody System Dynamics, Springer Verlag, 2015, 35 (1), pp.39-61. 10.1007/s11044-014-9437-4. hal-01088286 HAL Id: hal-01088286 https://hal.inria.fr/hal-01088286 Submitted on 29 Oct 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Singular mass matrix and redundant constraints in unilaterally constrained Lagrangian and Hamiltonian systems Bernard Brogliato · Daniel Goeleven Abstract This article deals with the analysis of the contact complementarity problem for Lagrangian systems subjected to unilateral constraints, and with a singular mass matrix and redundant constraints. Previous results by the authors on existence and uniqueness of solu- tions of some classes of variational inequalities are used to characterize the well-posedness of the contact problem. Criteria involving conditions on the tangent cone and the constraints gradient are given. It is shown that the proposed criteria easily extend to the case where the system is also subjected to a set of bilateral holonomic constraints, in addition to the unilat- eral ones.
    [Show full text]
  • Simple Harmonic Motion Coupled Oscillators Eigenvalue Problem
    Coupled Oscillators Concepts of primary interest: Simple Harmonic Motion Coupled Oscillators Eigenvalue problem Developing the equations of Motion: Intuitive Newton’s Law Method Lagrangian Mechanics Hamiltonian Mechanics Work-Energy for Mechanical Engineers Sample calculations: A matched equal mass pair of coupled oscillator – with variable coupling An oscillator pair with unequal masses Initial condition matching for near degenerate modes Comparison of two choices for generalized coordinates Applications Tools of the Trade Cramer’s Rule Coupled Differential Equation Principles Underlying the Coupled Oscillator Procedures Engineers can use this handout by replacing -2 by . A brief discussion of forced couple oscillations will be added. dd example with damping The Plan: A system of coupled linear oscillators is to be studied. To anchor the discussion, the solution to the single simple harmonic oscillator is reviewed. Next a system of two masses coupled by a linear spring is analyzed by guessing the solution. The result is a pair of two-particle sinusoidal oscillation patterns each with its characteristic frequency. More general systems are analyzed by assuming that sinusoidal solution modes exist and applying a diagonalization approach to identify them. The method can be made more formal by noting that it is a particular example of the matrix solution to systems of coupled linear differential equations. At several points, physical interpretations of results are presented that enhance our confidence in the procedure. Contact: [email protected] Just the single simple harmonic oscillator: The motion of a mass about its equilibrium position while subject to a linear restoring force is well- studied. The model equation is: ma mx kx where x is measured relative to the equilibrium position.
    [Show full text]
  • Rotordynamic Analyses Using Finite Element Method a Thesis
    Rotordynamic Analyses using Finite Element Method A Thesis Presented to the faculty of the School of Engineering and Applied Science University of Virginia In Partial Fulfillment of the requirements for the Degree Master of Science By: Michael Branagan May 2014 APPROVAL SHEET The Thesis is submitted in partial fulfillment of the requirements for the degree of Master of Science AUTHOR The thesis has been read and approved by the examining committee: Houston G. Wood, III Advisor Robert Durgin Rockwell, Jr Advisor Richard Wesley Kent Accepted for the School of Engineering and Applied Science: Dean, School of Engineering and Applied Science May 2014 Abstract The ability to perform critical frequency and forced response analyses is a vital tool for designing and troubleshooting rotordynamic systems. These analyses can be a challenge due to the wide variety of components found in rotordynamic systems. There are a variety of methods that can be used to perform these analyses. The finite element method is one method that can be used to perform these analyses. A finite element approach to rotordynamic analyses is presented in this thesis. First, the finite element model is developed. The equations of mo- tion that are used to analyze the finite element model are developed. Next the development of a MATLAB based rotordynamic tool incorporating this method is presented. The features of this tool, including degree-of-freedom coupling, multiple rotor systems, and the inclusion of tilting pad bearings with full coefficients, aerodynamic cross couplings, thrust bearings, flexible bearing supports, flexible couplings, and the stiffness properties of disks in the model, are presented as well as improvements that have been made to the efficiency of the program.
    [Show full text]
  • Forward and Inverse Dynamics of a Unicycle-Like Mobile Robot
    machines Article Forward and Inverse Dynamics of a Unicycle-Like Mobile Robot Carmine Maria Pappalardo * and Domenico Guida Department of Industrial Engineering, University of Salerno, Via Giovanni Paolo II 132, 84084 Fisciano, Italy; [email protected] * Correspondence: [email protected]; Tel.: +39-089-964-372 Received: 6 November 2018; Accepted: 8 January 2019; Published: 11 January 2019 Abstract: In this research work, a new method for solving forward and inverse dynamic problems of mechanical systems having an underactuated structure and subjected to holonomic and/or nonholonomic constraints is developed. The method devised in this paper is based on the combination of the Udwadia-Kalaba Equations with the Underactuation Equivalence Principle. First, an analytical method based on the Udwadia-Kalaba Equations is employed in the paper for handling dynamic and control problems of nonlinear nonholonomic mechanical systems in the same computational framework. Subsequently, the Underactuation Equivalence Principle is used for extending the capabilities of the Udwadia-Kalaba Equations from fully actuated mechanical systems to underactuated mechanical systems. The Underactuation Equivalence Principle represents an efficient method recently developed in the field of classical mechanics. The Underactuation Equivalence Principle is used in this paper for mathematically formalizing the underactuation property of a mechanical system considering a particular set of nonholonomic algebraic constraints defined at the acceleration level. On the other hand, in this study, the Udwadia-Kalaba Equations are analytically reformulated in a mathematical form suitable for treating inverse dynamic problems. By doing so, the Udwadia-Kalaba Equations are employed in conjunction with the Underactuation Equivalence Principle for developing a nonlinear control method based on an inverse dynamic approach.
    [Show full text]
  • Dynamics in About 100 Pages
    Intermediate Dynamics in about 100 pages Anindya Chatterjee Department of Mechanical Engineering Indian Institute of Technology Kanpur 208016 This version: November 4, 2014 2 Preface These notes started with a course at the Indian Institute of Science (Bangalore) on the dynamics and control of mechanical systems, of which I taught the “dynamics” part several times. I have since moved to IIT Kanpur via IIT Kharagpur. There are many excellent but fat books on dynamics. My notes, at least, are brief. I hope you will find my treatment of rotation matrices easy and useful (“rotate the body, not the axes”). You might enjoy my discussion of nonholonomic constraints in Lagrangian mechanics (“where do the λ’s come from?”). If you like something else, please send me email. Many students initially struggle with simple dynamics problems. It is like riding a bicycle. It looks easy until you try, and later you wonder what the fuss was about; but a journey lies in between. Several excellent people tried to teach me mechanics. Progress was slow. Professor A. C. Gomes of St. Xavier’s College, Kolkata, taught me well in my impressionable teens (1983-85). Professor Andy Ruina of Cornell University, my last formal teacher, helped a lot (1993-96). I took a long time to understand what I describe in these notes. But then dynamics is a great subject developed by great people. I am just a messenger, trying to bring the story to you. Several students, friends, and colleagues have helped with various parts of these notes (typing, figures, general discussion). I am embarrassed to say that I did not write down their names at the time, and so may miss someone now.
    [Show full text]