DOCSLIB.ORG

Introduction to Analytical Mechanics

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Introduction to Analytical Mechanics Introduction to Analytical Mechanics Lagrange’s Equations Siamak G. Faal 1/1 07/29/2016 1 Rigid body dynamics Dynamic equations: • Provide a relation between forces and system accelerations • Explain rigid body motion and trajectories Image sources (left to right): https://www.etsy.com/listing/96803935/funnel-shaped-spinning-top-natural-wood https://spaceflightnow.com/2016/06/30/juno-switched-to-autopilot-mode-for-jupiter-final-approach/ https://upload.wikimedia.org/wikipedia/commons/d/d8/NASA_Mars_Rover.jpg 2 http://www.popularmechanics.com/cars/a14665/why-car-suspensions-are-better-than-ever/ Kinematics vs Kinetics Kinematics Kinetics Focuses on pure motion of Relationship between the particles and bodies (without motion of bodies and forces considering masses and and torques forces) Image sources (left to right): https://www.youtube.com/watch?v=PAHjo4DXOjc 3 http://www.timelinecoverbanner.com/facebook-covers/billiard.jpg Dynamic Forces and Reactions Dynamic analysis are essential to determine reaction forces that are cussed by dynamic forces (e.g. Accelerations) Image sources (left to right): http://r.hswstatic.com/w_404/gif/worlds-greatest-roller-coasters-149508746.jpg http://www.murraymitchell.com/wp-content/uploads/2011/04/red_aircraft_performing_aerobatics.jpg 4 Differential Equations of Motion provide a set of differential equations which allow modeling, simulation and control of dynamic mechanical systems Why humanoid robots do not walk or act like humans !? 5 Background Information 6 Notations Throughout this presentation: Parameters Definitions • Vectors are noted with overlines: subscript G The quantity related to (e.g. ) center of mass • The first derivative with respect to� Linear velocity Angular velocity time: = Linear Acceleration ̇ • The second derivative with respect to Resultant force vector Gravitational acceleration time: = Σ� 2 Generalized coordinate 2 ̈ 7 Time Derivative of a Vector Time derivative of a vector in rotating coordinate frame is: = + × ̇ Angular velocity of the Relative rate of change of frame the vector 8 Newtonian Dynamics Newton laws of motion: I. The velocity of a particle can be changed by the application of a force II. The resultant force is proportional to the acceleration of particle with the factor of mass = Σ� � III. The forces acting on a body result from an interaction with another body. The action and reaction forces have the same magnitude, the same line of action but opposite directions 9 Newtonian Equations of Motion For each rigid body: Linear momentum of the center of mass: = = ̅ ̅ Σ� ̅ = /× + Σ� ̅ � �̇ Time derivative of the angular momentum of Vector from arbitrary point the object about point A A to G (center of mass) Acceleration of the = + × arbitrary point A �̇ �̇ � � 10 Newtonian Equations of Motion A simplified form of the equations are achieved by considering the following assumptions: . Mass of the object is constant . Moment equations are calculated at the center of mass = =Σ� +� × Σ� �̇ � � 11 Energy Kinetic energy of a rigid body: 1 1 = , + , 2 2 Potential energy of an object with mass located at altitude due to gravity: ℎ = Potential energy stored in an spring withℎ stiffness : 1 = 2 Δ 12 Degrees of Freedom The number of independent parameters that define system configuration In 3D space In planar motion A free object has A free object has 6 3 DOF DOF 13 Mechanical Joints Revolute joint Planar joint Adds 5 constraints Adds 3 constraints Relative DOF: 1 Relative DOF: 3 Cylindrical joint Prismatic joint Spherical joint Adds 4 constraints Adds 5 constraints Adds 3 constraints Relative DOF: 2 Relative DOF: 1 Relative DOF: 3 14 Generalized Coordinates A set of geometrical parameters which uniquely defines the position / configuration of the system The minimum number of generalized coordinates that required to specify the position of the system is equal to the number of DOF of the system (unconstrained or independent coordinates) Generalized coordinates do not form a unique set! Generalized coordinates are commonly notes with 15 Generalized Coordinates Q) What can be a set of generalized coordinates for a simple rod in on a plane. A) 16 Generalized Coordinates Q) Can you think about any other generalized coordinates set for the same system? A) 17 Generalized Coordinates Q) How about another one? A) This is an ambiguous set of coordinates! 18 Generalized Coordinates Here we have more coordinates than the digress Q) Let’s do one more of freedom of the system, Thus the coordinates can not have any arbitrary value and they need to satisfy two constrain equations: = cos A) = sin − − 19 Constraints Any generalized coordinate which can not have any arbitrary value is a constrained coordinate Relation between constrained generalized coordinates is called constraint equation : Number of constraint equations − Number of constrained Number of degrees generalized coordinates of freedom 20 Virtual Displacement In a virtual movement, the 3 generalized coordinates of the system are considered to be incremented by infinitesimal amount from the values they have at an δ̂ arbitraryδ instant, with time held constant. = 2 δ̂ � δ̂ =1 1 21 Virtual Displacement Example: A B Generalized coordinate: C Q) What are the virtual displacements of pin F and virtual rotation of bar EF resulting from D E δ F 22 Virtual Displacement Solution: A B = cos −1 3 C = +2 4 / / 2 2 2 2 1 2 ̅ ̅ − ̅ D E = cos = 4 / −1 δ δ δ − 2 2 1 2 2 1 − = = / / 2 2 4 / 3 δ̅ ̅ δ ̅ − 2 2 1 2 ̅ − F 23 Virtual Work When a particle is given a virtual displacement , the forces acting on the particle do virtual work δ̅ δ Since is infinitesimal, is infinitesimal as well δ̅ δ Because the change is virtual, time is held constant at an arbitrary value The virtual work done by constraint or internal forces is zero! 24 Generalized Forces = δ � � ⋅ δ̅ = ̅ � → δ � ⋅ � δ =1 = ̅ � δ � � ⋅ δ =1 Generalized force 25 Joseph-Louis Lagrange 25 January 1736 – 10 April 1813, Italian Enlightenment Era mathematician and astronomer. Lagrange made significant contributions to the fields of analysis, number theory, and both classical and celestial mechanics. Lagrange’s Equations 26 Lagrange’s Equations of Motion + = , = 1, 2, … , − ̇ Alternative (common) form: = ℒ − = , = 1, 2, … , ℒ ℒ − ̇ 27 Lagrange’s Equations of Motion . Directly provides the differential equations of motion without any need for solving a system of equations. Provides a much easier method to solve multi-body systems . It is very systematic and can be used to automatically generate differential equations of motion for different systems (As long as the derivatives are computed correctly) 28 Example 1 Find differential equations of motion for the following system using both Newtonian mechanics and Lagrange’s equations 29 Example 1 Newtonian approach: Spring force is equal to: = 0 Assume = 0 − − 0 = Newton’s second law for→ x-axis:− Free-body diagram: = Since the only force appliedΣ to thë system is the spring force: = = So the differential equationΣ of motion− for the system is: + = 0 30 ̈ Example 1 Lagrangian approach: The kinetic energy of the system is: = 1 2 2 ̇ The potential energy of the system is: = 1 2 2 Since there is no force acting on the system, the generalized Calculations: force is zero: = 0 Lagrangian of the system is: = ℒ 1 1 ̇ = ̇ 2 2 = 2 2 ℒ ̇ − ℒ Thus, the differential equation of motion is: ̈ ̇ = = 0 = 0 ℒ − ℒ ℒ − → ̈ − 31 ̇ Example 2 Find differential equations of motion for the following system using both Newtonian mechanics and Lagrange’s equations 32 Example 2 Newtonian approach: Spring force is equal to: = 0 Assume = 0 − − 0 = → − Newton’s second law for x-axis: Free-body diagram: = The summation of forcesΣ along ẍ axis is: = So the differential equationΣ of −motion for the system is: + = ̈ 33 Example 2 Lagrangian approach: The kinetic energy of the system is: = 1 2 2 ̇ The potential energy of the system is: = 1 2 2 The generalized force acting on the system is: Calculations: = = ̅ Σ� ⋅ ⋅ = Lagrangian of the system is: ℒ ̇ 1 1 ̇ = 2 2 = 2 2 ℒ ℒ ̇ − ̈ Thus, the differential equation of motion is: ̇ = = + = ℒ − ℒ ℒ − → ̈ 34 ̇ Example 3 Find differential equations of motion for the following system using both Newtonian mechanics and Lagrange’s equations , 35 Example 3 Newtonian approach: , Moment of inertia of a rod: 1 1 3 12 2 2 Free-body diagram: Moment equation around point O gives us: 1 2 = cos 3 2 2 ̈ − 36 Example 3 Lagrangian approach: Center of mass position in terms of generalized coordinates , = cos + sin 2 2 ̇ =̅ sin ̅+ cos ̅ → ̅ 2 2 ̇ ̇ The kinetic energy− of the system̅ is: ̅ 1 1 1 1 = , + , = + 2 2 8 24 2 2 2 2 1 ̇ ̇ ̇ ̇ = ̅ ̅ � � 6 2 2 The potential ̇ energy of the system is: = = sin 2 ̅ ⋅ ̅ 37 Example 3 Lagrangian approach (cont.): The generalized force is , = = = ̅ � Σ ⋅ Lagrangian of the system is: Calculations: 1 = + sin 1 6 2 = 2 2 3 ℒ ̇ ℒ 2 Thus, the differential equation of motion is: ̇ ̇ 1 = 3 = ℒ 2 ̈ ℒ ℒ ̇ − = cos ̇ 2 1 ℒ cos = 3 2 2 38 → − References [1] Ginsberg, Jerry H. Advanced engineering dynamics. Cambridge University Press, 1998. [2] Greenwood, Donald T. Advanced dynamics. Cambridge University Press, 2006. 39.
Load more

Recommended publications
  • Analytical Mechanics
    A Guided Tour of Analytical Mechanics with animations in MAPLE Rouben Rostamian Department of Mathematics and Statistics UMBC [email protected] December 2, 2018 ii Contents Preface vii 1 An introduction through examples 1 1.1 ThesimplependulumàlaNewton ...................... 1 1.2 ThesimplependulumàlaEuler ....................... 3 1.3 ThesimplependulumàlaLagrange.. .. .. ... .. .. ... .. .. .. 3 1.4 Thedoublependulum .............................. 4 Exercises .......................................... .. 6 2 Work and potential energy 9 Exercises .......................................... .. 12 3 A single particle in a conservative force field 13 3.1 The principle of conservation of energy . ..... 13 3.2 Thescalarcase ................................... 14 3.3 Stability....................................... 16 3.4 Thephaseportraitofasimplependulum . ... 16 Exercises .......................................... .. 17 4 TheKapitsa pendulum 19 4.1 Theinvertedpendulum ............................. 19 4.2 Averaging out the fast oscillations . ...... 19 4.3 Stabilityanalysis ............................... ... 22 Exercises .......................................... .. 23 5 Lagrangian mechanics 25 5.1 Newtonianmechanics .............................. 25 5.2 Holonomicconstraints............................ .. 26 5.3 Generalizedcoordinates .......................... ... 29 5.4 Virtual displacements, virtual work, and generalized force....... 30 5.5 External versus reaction forces . ..... 32 5.6 The equations of motion for a holonomic system . ...
    [Show full text]
  • Lagrangian Mechanics - Wikipedia, the Free Encyclopedia Page 1 of 11
    Lagrangian mechanics - Wikipedia, the free encyclopedia Page 1 of 11 Lagrangian mechanics From Wikipedia, the free encyclopedia Lagrangian mechanics is a re-formulation of classical mechanics that combines Classical mechanics conservation of momentum with conservation of energy. It was introduced by the French mathematician Joseph-Louis Lagrange in 1788. Newton's Second Law In Lagrangian mechanics, the trajectory of a system of particles is derived by solving History of classical mechanics · the Lagrange equations in one of two forms, either the Lagrange equations of the Timeline of classical mechanics [1] first kind , which treat constraints explicitly as extra equations, often using Branches [2][3] Lagrange multipliers; or the Lagrange equations of the second kind , which Statics · Dynamics / Kinetics · Kinematics · [1] incorporate the constraints directly by judicious choice of generalized coordinates. Applied mechanics · Celestial mechanics · [4] The fundamental lemma of the calculus of variations shows that solving the Continuum mechanics · Lagrange equations is equivalent to finding the path for which the action functional is Statistical mechanics stationary, a quantity that is the integral of the Lagrangian over time. Formulations The use of generalized coordinates may considerably simplify a system's analysis. Newtonian mechanics (Vectorial For example, consider a small frictionless bead traveling in a groove. If one is tracking the bead as a particle, calculation of the motion of the bead using Newtonian mechanics) mechanics would require solving for the time-varying constraint force required to Analytical mechanics: keep the bead in the groove. For the same problem using Lagrangian mechanics, one Lagrangian mechanics looks at the path of the groove and chooses a set of independent generalized Hamiltonian mechanics coordinates that completely characterize the possible motion of the bead.
    [Show full text]
  • A Simple Treatment of Constraint Forces and Constraint Moments in the Dynamics of Rigid Bodies
    ASME Applied Mechanics Reviews http://dx.doi.org/10.1115/1.4028099 This paper appeared in Vol.67(1), 014801, 2015. (Online: Aug 25, 2014) (8 pages) A Simple Treatment of Constraint Forces and Constraint Moments in the Dynamics of Rigid Bodies Oliver M. O’Reilly · Arun R. Srinivasa November 6, 2014 Abstract In this expository article, a simple concise a set of 3-1-3 Euler angles ψ,θ, and φ, and the vertical treatment of Lagrange’s prescription for constraint forces coordinate zP of the point of contact XP . The coordi- and constraint moments in the dynamics of rigid bodies nate zP is related to the vertical coordinate z of the is presented. The treatment is suited to both Newton- center of mass by the simple relation Euler and Lagrangian treatments of rigid body dynam- 6 ics and is illuminated with a range of examples from q = zP = z − R sin(θ). (1) classical mechanics and orthopedic biomechanics. After choosing the coordinate system, we can calculate Keywords Lagrange’s equations · Constraints · the kinetic T and potential U energies of the disk. For 6 Constraints forces · Constraint moments the disk in motion on the plane, the coordinate q is constrained to be zero. As a result, the equations of motion have a Lagrange multiplier λ on the right hand 1 Introduction side: d ∂T ∂T ∂U − + = ΦΓ , (Γ =1,..., 6) , (2) Consider formulating the equations of motion for the dt ∂q˙Γ ∂qΓ ∂qΓ classical problem of a thin circular disk shown in Fig- ure 1 that is sliding with a point XP in contact with a where the generalized constraint force is smooth horizontal surface.
    [Show full text]
  • A Contribution on Modeling Methodologies for Multibody Systems
    RENATO MAIA MATARAZZO ORSINO A CONTRIBUTION ON MODELING METHODOLOGIES FOR MULTIBODY SYSTEMS Tese apresentada à Escola Politécnica da Universidade de São Paulo para obtenção do Título de Doutor em Engenharia Mecânica. São Paulo 2016 RENATO MAIA MATARAZZO ORSINO A CONTRIBUTION ON MODELING METHODOLOGIES FOR MULTIBODY SYSTEMS Tese apresentada à Escola Politécnica da Universidade de São Paulo para obtenção do Título de Doutor em Engenharia Mecânica. Área de concentração: Engenharia Mecânica Orientador: Prof. Dr. Tarcisio Antonio Hess Coelho São Paulo 2016 Este exemplar foi revisado e corrigido em relação à versão original, sob responsabilidade única do autor e com a anuência de seu orientador. São Paulo, ______ de ____________________ de __________ Assinatura do autor: ________________________ Assinatura do orientador: ________________________ Catalogação-na-publicação Orsino, Renato Maia Matarazzo A contribution on modeling methodologies for multibody systems / R. M. M. Orsino -- versão corr. -- São Paulo, 2016. 213 p. Tese (Doutorado) - Escola Politécnica da Universidade de São Paulo. Departamento de Engenharia Mecânica. 1.Dinâmica 2.Cinemática 3.Mecanismos I.Universidade de São Paulo. Escola Politécnica. Departamento de Engenharia Mecânica II.t. Acknowledgements I would like to acknowledge everyone who supported, advised and encou- raged me during the last 4 years, since I begun this long journey towards my PhD degree. Professors, family, friends, colleagues, staff (both from University of São Paulo and TU Delft), thank you all! A second acknowledgment to CNPq – Conselho Nacional de Desen- volvimento Científico e Tecnológico, of the Ministry of Science, Technology and Innovation of Brazil for sponsoring this research (processes numbers 163026/2012-2 and 203112/2014-8). I would like to dedicate the remaining of this section for some special ack- nowledgements.
    [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]