Leonhard Euler (1707-1783) and Rigid Body Dynamics
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Glossary Physics (I-Introduction)
1 Glossary Physics (I-introduction) - Efficiency: The percent of the work put into a machine that is converted into useful work output; = work done / energy used [-]. = eta In machines: The work output of any machine cannot exceed the work input (<=100%); in an ideal machine, where no energy is transformed into heat: work(input) = work(output), =100%. Energy: The property of a system that enables it to do work. Conservation o. E.: Energy cannot be created or destroyed; it may be transformed from one form into another, but the total amount of energy never changes. Equilibrium: The state of an object when not acted upon by a net force or net torque; an object in equilibrium may be at rest or moving at uniform velocity - not accelerating. Mechanical E.: The state of an object or system of objects for which any impressed forces cancels to zero and no acceleration occurs. Dynamic E.: Object is moving without experiencing acceleration. Static E.: Object is at rest.F Force: The influence that can cause an object to be accelerated or retarded; is always in the direction of the net force, hence a vector quantity; the four elementary forces are: Electromagnetic F.: Is an attraction or repulsion G, gravit. const.6.672E-11[Nm2/kg2] between electric charges: d, distance [m] 2 2 2 2 F = 1/(40) (q1q2/d ) [(CC/m )(Nm /C )] = [N] m,M, mass [kg] Gravitational F.: Is a mutual attraction between all masses: q, charge [As] [C] 2 2 2 2 F = GmM/d [Nm /kg kg 1/m ] = [N] 0, dielectric constant Strong F.: (nuclear force) Acts within the nuclei of atoms: 8.854E-12 [C2/Nm2] [F/m] 2 2 2 2 2 F = 1/(40) (e /d ) [(CC/m )(Nm /C )] = [N] , 3.14 [-] Weak F.: Manifests itself in special reactions among elementary e, 1.60210 E-19 [As] [C] particles, such as the reaction that occur in radioactive decay. -
Non-Local Momentum Transport Parameterizations
Non-local Momentum Transport Parameterizations Joe Tribbia NCAR.ESSL.CGD.AMP Outline • Historical view: gravity wave drag (GWD) and convective momentum transport (CMT) • GWD development -semi-linear theory -impact • CMT development -theory -impact Both parameterizations of recent vintage compared to radiation or PBL GWD CMT • 1960’s discussion by Philips, • 1972 cumulus vorticity Blumen and Bretherton damping ‘observed’ Holton • 1970’s quantification Lilly • 1976 Schneider and and momentum budget by Lindzen -Cumulus Friction Swinbank • 1980’s NASA GLAS model- • 1980’s incorporation into Helfand NWP and climate models- • 1990’s pressure term- Miller and Palmer and Gregory McFarlane Atmospheric Gravity Waves Simple gravity wave model Topographic Gravity Waves and Drag • Flow over topography generates gravity (i.e. buoyancy) waves • <u’w’> is positive in example • Power spectrum of Earth’s topography α k-2 so there is a lot of subgrid orography • Subgrid orography generating unresolved gravity waves can transport momentum vertically • Let’s parameterize this mechanism! Begin with linear wave theory Simplest model for gravity waves: with Assume w’ α ei(kx+mz-σt) gives the dispersion relation or Linear theory (cont.) Sinusoidal topography ; set σ=0. Gives linear lower BC Small scale waves k>N/U0 decay Larger scale waves k<N/U0 propagate Semi-linear Parameterization Propagating solution with upward group velocity In the hydrostatic limit The surface drag can be related to the momentum transport δh=isentropic Momentum transport invariant by displacement Eliassen-Palm. Deposited when η=U linear theory is invalid (CL, breaking) z φ=phase Gravity Wave Drag Parameterization Convective or shear instabilty begins to dissipate wave- momentum flux no longer constant Waves propagate vertically, amplitude grows as r-1/2 (energy force cons.). -
10. Collisions • Use Conservation of Momentum and Energy and The
10. Collisions • Use conservation of momentum and energy and the center of mass to understand collisions between two objects. • During a collision, two or more objects exert a force on one another for a short time: -F(t) F(t) Before During After • It is not necessary for the objects to touch during a collision, e.g. an asteroid flied by the earth is considered a collision because its path is changed due to the gravitational attraction of the earth. One can still use conservation of momentum and energy to analyze the collision. Impulse: During a collision, the objects exert a force on one another. This force may be complicated and change with time. However, from Newton's 3rd Law, the two objects must exert an equal and opposite force on one another. F(t) t ti tf Dt From Newton'sr 2nd Law: dp r = F (t) dt r r dp = F (t)dt r r r r tf p f - pi = Dp = ò F (t)dt ti The change in the momentum is defined as the impulse of the collision. • Impulse is a vector quantity. Impulse-Linear Momentum Theorem: In a collision, the impulse on an object is equal to the change in momentum: r r J = Dp Conservation of Linear Momentum: In a system of two or more particles that are colliding, the forces that these objects exert on one another are internal forces. These internal forces cannot change the momentum of the system. Only an external force can change the momentum. The linear momentum of a closed isolated system is conserved during a collision of objects within the system. -
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. -
Fluids – Lecture 7 Notes 1
Fluids – Lecture 7 Notes 1. Momentum Flow 2. Momentum Conservation Reading: Anderson 2.5 Momentum Flow Before we can apply the principle of momentum conservation to a fixed permeable control volume, we must first examine the effect of flow through its surface. When material flows through the surface, it carries not only mass, but momentum as well. The momentum flow can be described as −→ −→ momentum flow = (mass flow) × (momentum /mass) where the mass flow was defined earlier, and the momentum/mass is simply the velocity vector V~ . Therefore −→ momentum flow =m ˙ V~ = ρ V~ ·nˆ A V~ = ρVnA V~ where Vn = V~ ·nˆ as before. Note that while mass flow is a scalar, the momentum flow is a vector, and points in the same direction as V~ . The momentum flux vector is defined simply as the momentum flow per area. −→ momentum flux = ρVn V~ V n^ . mV . A m ρ Momentum Conservation Newton’s second law states that during a short time interval dt, the impulse of a force F~ applied to some affected mass, will produce a momentum change dP~a in that affected mass. When applied to a fixed control volume, this principle becomes F dP~ a = F~ (1) dt V . dP~ ˙ ˙ P(t) P + P~ − P~ = F~ (2) . out dt out in Pin 1 In the second equation (2), P~ is defined as the instantaneous momentum inside the control volume. P~ (t) ≡ ρ V~ dV ZZZ ˙ The P~ out is added because mass leaving the control volume carries away momentum provided ˙ by F~ , which P~ alone doesn’t account for. -
The Bernoulli Edition the Collected Scientific Papers of the Mathematicians and Physicists of the Bernoulli Family
Bernoulli2005.qxd 24.01.2006 16:34 Seite 1 The Bernoulli Edition The Collected Scientific Papers of the Mathematicians and Physicists of the Bernoulli Family Edited on behalf of the Naturforschende Gesellschaft in Basel and the Otto Spiess-Stiftung, with support of the Schweizerischer Nationalfonds and the Verein zur Förderung der Bernoulli-Edition Bernoulli2005.qxd 24.01.2006 16:34 Seite 2 The Scientific Legacy Èthe Bernoullis' contributions to the theory of oscillations, especially Daniel's discovery of of the Bernoullis the main theorems on stationary modes. Johann II considered, but rejected, a theory of Modern science is predominantly based on the transversal wave optics; Jacob II came discoveries in the fields of mathematics and the tantalizingly close to formulating the natural sciences in the 17th and 18th centuries. equations for the vibrating plate – an Eight members of the Bernoulli family as well as important topic of the time the Bernoulli disciple Jacob Hermann made Èthe important steps Daniel Bernoulli took significant contributions to this development in toward a theory of errors. His efforts to the areas of mathematics, physics, engineering improve the apparatus for measuring the and medicine. Some of their most influential inclination of the Earth's magnetic field led achievements may be listed as follows: him to the first systematic evaluation of ÈJacob Bernoulli's pioneering work in proba- experimental errors bility theory, which included the discovery of ÈDaniel's achievements in medicine, including the Law of Large Numbers, the basic theorem the first computation of the work done by the underlying all statistical analysis human heart. -
The Bernoulli Family
Mathematical Discoveries of the Bernoulli Brothers Caroline Ellis Union University MAT 498 November 30, 2001 Bernoulli Family Tree Nikolaus (1623-1708) Jakob I Nikolaus I Johann I (1654-1705) (1662-1716) (1667-1748) Nikolaus II Nikolaus III Daniel I Johann II (1687-1759) (1695-1726) (1700-1782) (1710-1790) This Swiss family produced eight mathematicians in three generations. We will focus on some of the mathematical discoveries of Jakob l and his brother Johann l. Some History Nikolaus Bernoulli wanted Jakob to be a Protestant pastor and Johann to be a doctor. They obeyed their father and earned degrees in theology and medicine, respectively. But… Some History, cont. Jakob and Johann taught themselves the “new math” – calculus – from Leibniz‟s notes and papers. They started to have contact with Leibniz, and are now known as his most important students. http://www-history.mcs.st-andrews.ac.uk/history/PictDisplay/Leibniz.html Jakob Bernoulli (1654-1705) learned about mathematics and astronomy studied Descarte‟s La Géometrie, John Wallis‟s Arithmetica Infinitorum, and Isaac Barrow‟s Lectiones Geometricae convinced Leibniz to change the name of the new math from calculus sunmatorius to calculus integralis http://www-history.mcs.st-andrews.ac.uk/history/PictDisplay/Bernoulli_Jakob.html Johann Bernoulli (1667-1748) studied mathematics and physics gave calculus lessons to Marquis de L‟Hôpital Johann‟s greatest student was Euler won the Paris Academy‟s biennial prize competition three times – 1727, 1730, and 1734 http://www-history.mcs.st-andrews.ac.uk/history/PictDisplay/Bernoulli_Johann.html Jakob vs. Johann Johann Bernoulli had greater intuitive power and descriptive ability Jakob had a deeper intellect but took longer to arrive at a solution Famous Problems the catenary (hanging chain) the brachistocrone (shortest time) the divergence of the harmonic series (1/n) The Catenary: Hanging Chain Jakob Bernoulli Galileo guessed that proposed this problem this curve was a in the May 1690 edition parabola, but he never of Acta Eruditorum. -
Lie Group Formulation of Articulated Rigid Body Dynamics
Lie Group Formulation of Articulated Rigid Body Dynamics Junggon Kim 12/10/2012, Ver 2.01 Abstract It has been usual in most old-style text books for dynamics to treat the formulas describing linear(or translational) and angular(or rotational) motion of a rigid body separately. For example, the famous Newton's 2nd law, f = ma, for the translational motion of a rigid body has its partner, so-called the Euler's equation which describes the rotational motion of the body. Separating translation and rotation, however, causes a huge complexity in deriving the equations of motion of articulated rigid body systems such as robots. In Section1, an elegant single equation of motion of a rigid body moving in 3D space is derived using a Lie group formulation. In Section2, the recursive Newton-Euler algorithm (inverse dynamics), Articulated-Body algorithm (forward dynamics) and a generalized recursive algorithm (hybrid dynamics) for open chains or tree-structured articulated body systems are rewritten with the geometric formulation for rigid body. In Section3, dynamics of constrained systems such as a closed loop mechanism will be described. Finally, in Section4, analytic derivatives of the dynamics algorithms, which would be useful for optimization and sensitivity analysis, are presented.1 1 Dynamics of a Rigid Body This section describes the equations of motion of a single rigid body in a geometric manner. 1.1 Rigid Body Motion To describe the motion of a rigid body, we need to represent both the position and orien- tation of the body. Let fBg be a coordinate frame attached to the rigid body and fAg be an arbitrary coordinate frame, and all coordinate frames will be right-handed Cartesian from now on. -
Lecture 1: Introduction
Lecture 1: Introduction E. J. Hinch Non-Newtonian fluids occur commonly in our world. These fluids, such as toothpaste, saliva, oils, mud and lava, exhibit a number of behaviors that are different from Newtonian fluids and have a number of additional material properties. In general, these differences arise because the fluid has a microstructure that influences the flow. In section 2, we will present a collection of some of the interesting phenomena arising from flow nonlinearities, the inhibition of stretching, elastic effects and normal stresses. In section 3 we will discuss a variety of devices for measuring material properties, a process known as rheometry. 1 Fluid Mechanical Preliminaries The equations of motion for an incompressible fluid of unit density are (for details and derivation see any text on fluid mechanics, e.g. [1]) @u + (u · r) u = r · S + F (1) @t r · u = 0 (2) where u is the velocity, S is the total stress tensor and F are the body forces. It is customary to divide the total stress into an isotropic part and a deviatoric part as in S = −pI + σ (3) where tr σ = 0. These equations are closed only if we can relate the deviatoric stress to the velocity field (the pressure field satisfies the incompressibility condition). It is common to look for local models where the stress depends only on the local gradients of the flow: σ = σ (E) where E is the rate of strain tensor 1 E = ru + ruT ; (4) 2 the symmetric part of the the velocity gradient tensor. The trace-free requirement on σ and the physical requirement of symmetry σ = σT means that there are only 5 independent components of the deviatoric stress: 3 shear stresses (the off-diagonal elements) and 2 normal stress differences (the diagonal elements constrained to sum to 0). -
Rigid Body Dynamics 2
Rigid Body Dynamics 2 CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2017 Cross Product & Hat Operator Derivative of a Rotating Vector Let’s say that vector r is rotating around the origin, maintaining a fixed distance At any instant, it has an angular velocity of ω ω dr r ω r dt ω r Product Rule The product rule of differential calculus can be extended to vector and matrix products as well da b da db b a dt dt dt dab da db b a dt dt dt dA B dA dB B A dt dt dt Rigid Bodies We treat a rigid body as a system of particles, where the distance between any two particles is fixed We will assume that internal forces are generated to hold the relative positions fixed. These internal forces are all balanced out with Newton’s third law, so that they all cancel out and have no effect on the total momentum or angular momentum The rigid body can actually have an infinite number of particles, spread out over a finite volume Instead of mass being concentrated at discrete points, we will consider the density as being variable over the volume Rigid Body Mass With a system of particles, we defined the total mass as: n m m i i1 For a rigid body, we will define it as the integral of the density ρ over some volumetric domain Ω m d Angular Momentum The linear momentum of a particle is 퐩 = 푚퐯 We define the moment of momentum (or angular momentum) of a particle at some offset r as the vector 퐋 = 퐫 × 퐩 Like linear momentum, angular momentum is conserved in a mechanical system If the particle is constrained only -
Leonhard Euler: His Life, the Man, and His Works∗
SIAM REVIEW c 2008 Walter Gautschi Vol. 50, No. 1, pp. 3–33 Leonhard Euler: His Life, the Man, and His Works∗ Walter Gautschi† Abstract. On the occasion of the 300th anniversary (on April 15, 2007) of Euler’s birth, an attempt is made to bring Euler’s genius to the attention of a broad segment of the educated public. The three stations of his life—Basel, St. Petersburg, andBerlin—are sketchedandthe principal works identified in more or less chronological order. To convey a flavor of his work andits impact on modernscience, a few of Euler’s memorable contributions are selected anddiscussedinmore detail. Remarks on Euler’s personality, intellect, andcraftsmanship roundout the presentation. Key words. LeonhardEuler, sketch of Euler’s life, works, andpersonality AMS subject classification. 01A50 DOI. 10.1137/070702710 Seh ich die Werke der Meister an, So sehe ich, was sie getan; Betracht ich meine Siebensachen, Seh ich, was ich h¨att sollen machen. –Goethe, Weimar 1814/1815 1. Introduction. It is a virtually impossible task to do justice, in a short span of time and space, to the great genius of Leonhard Euler. All we can do, in this lecture, is to bring across some glimpses of Euler’s incredibly voluminous and diverse work, which today fills 74 massive volumes of the Opera omnia (with two more to come). Nine additional volumes of correspondence are planned and have already appeared in part, and about seven volumes of notebooks and diaries still await editing! We begin in section 2 with a brief outline of Euler’s life, going through the three stations of his life: Basel, St. -
Navier-Stokes-Equation
Math 613 * Fall 2018 * Victor Matveev Derivation of the Navier-Stokes Equation 1. Relationship between force (stress), stress tensor, and strain: Consider any sub-volume inside the fluid, with variable unit normal n to the surface of this sub-volume. Definition: Force per area at each point along the surface of this sub-volume is called the stress vector T. When fluid is not in motion, T is pointing parallel to the outward normal n, and its magnitude equals pressure p: T = p n. However, if there is shear flow, the two are not parallel to each other, so we need a marix (a tensor), called the stress-tensor , to express the force direction relative to the normal direction, defined as follows: T Tn or Tnkjjk As we will see below, σ is a symmetric matrix, so we can also write Tn or Tnkkjj The difference in directions of T and n is due to the non-diagonal “deviatoric” part of the stress tensor, jk, which makes the force deviate from the normal: jkp jk jk where p is the usual (scalar) pressure From general considerations, it is clear that the only source of such “skew” / ”deviatoric” force in fluid is the shear component of the flow, described by the shear (non-diagonal) part of the “strain rate” tensor e kj: 2 1 jk2ee jk mm jk where euujk j k k j (strain rate tensro) 3 2 Note: the funny construct 2/3 guarantees that the part of proportional to has a zero trace. The two terms above represent the most general (and the only possible) mathematical expression that depends on first-order velocity derivatives and is invariant under coordinate transformations like rotations.