DOCSLIB.ORG

Rigid Body Dynamics Huamin Wang

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Rigid Body Dynamics Huamin Wang Rigid Body Dynamics Huamin Wang The motion of a rigid body contains two parts: linear motion and angular motion. In this docu- ment, we will study both motions. We will also discuss how to detect collisions and estimate collision responses. To simplify our notations, we will follow the standard formats: bold fonts are vectors, italic fonts are scalars, and UPPERCASES are matrices. The physics quantities associated with linear motion are: position, velocity, force, ... The physics quantities associated with angular motion are: orientation, angular velocity, torque, ... 1 Linear Motion Let us first examine the linear motion of a rigid body. Without considering the orientation, this is equivalent to simulating the motion of the body mass center x(t), which changes as a function of time t. Assuming that the initial model is centered at the origin as shown in Figure 1a. We can simply treat x(t) as the position of the rigid body. To model the linear motion, we also define the derivative of x(t) to time t, which is commonly known as the linear velocity v(t). To animate the linear motion, our goal is to update the position and the velocity of the center, each time the Update function gets called. Let the time of the last Update call be t0 and the time of the current call be t1. The time difference ∆t = t1 − t0 is known as the time step. Given the position and the velocity at t0, we calculate the new position and the new velocity as: ( 1 R t1 v(t1) = v(t0) + f(t)dt ≈ v(t0) + ∆tf(t0)=M; M t0 (1) x(t ) = x(t ) + R t1 v(t)dt ≈ x(t ) + ∆tv(t ); 1 0 t0 0 1 in which M is the mass of the object and f(t0) is the total force evaluated at time t0. If we assume that the mass distributed at each vertex i is mi, the total mass is simply their sum: X M = mi: (2) i Note that Equation 6 is just an approximation. Formally, this is called explicit time integration. We will discuss other time integrators later. One question here is how to evaluate f(t0). This vector represents the total force applied on the gravity 2 whole rigid body. It contains the gravity force f (t0) = Mg, where g = [0; −9:8m=s ; 0] is the gravity acceleration. It contains other forces as well, such as the collision force. Many simulators also model the damping effect by some damping force. Let us do not do this for now. 2 Angular Motion To simulate the angular motion of a rigid body, we must define everything in an angular way. 2.1 Quaternion The first and foremost question is how to define the orientation, i.e., \angular position". One idea is to use the matrix. However, only an orthonormal matrix can be a rotational matrix and we must keep enforcing the orthonormality, which is not very convenient. Alternatively, we can define a bunch 1 ܘ௜ ௜ሺݐሻܚ ሺݐሻܠ ۼ ௜ܘ ܘ௜ ௜ሺݐሻܚ ௜ሺݐሻܚ ௜ሺݐሻܚ+ሺݐሻܠ ሺݐሻ ሻ ሻܠ ሺݐሻܠ ௜ሺݐܚ+ ሺݐܠ ௜ሺݐሻܚ+ሺݐሻܠ ௜ሺݐሻܚ ሺݐሻ+࣓ሺݐሻ ൈܞ (a) Rest configuration (b) Moved configuration (c) Collision configuration Figure 1: A rigid body under different configurations. of rotation angles around X, Y, and Z axes, known as Euler angles. Again, that is not easy to derive the formula for their time evolution. Our solution is to use quaternion. A quaternion is a 4D vector q = [x; y; z; w], or q = [v; w], which contains a 3D vector v and a scalar w. A quaternion must be normalized: x2 + y2 + z2 + w2 = 1: (3) A quaternion represents a rotation around axis v with an angle θ, such as w = cos(θ=2) and x2 + y2 + z2 = sin2(θ=2). A quaternion can be converted into a rotation matrix as: 2 w2 + x2 − y2 − z2 2(xy − wz) 2(xz + wy) 3 R = 4 2(xy + wz) w2 − x2 + y2 − z2 2(yz − xw) 5 : (4) 2(xz − wy) 2(yz + xw) w2 − x2 − y2 + z2 Compared with a rotational matrix, a quaternion is very easy to maintain: we just have to make it normalized all the time. Quaternion arithmetic is defined as follows. Let q1 = [v1; w1] and q2 = [v2; w2] be two quaternion vectors and s be a scalar. We have: q1 + q2 = [v1 + v2; w1 + w2]; sq1 = [sv1; sw1]; (5) q1q2 = [w1v2 + w2v1 + v1 × v2; w1w2 − v1 · v2]: Intuitively, q1q2 means rotation by q2, and then rotation by q1. 2.2 Time Evolution The angular motion of a rigid body can be updated as follows: ( R t1 −1 −1 !(t1) = !(t0) + I (t)τ (t)dt ≈ !(t0) + ∆tI (t0)τ (t0); t0 (6) q(t ) = q(t ) + 1 R t1 !^(t)q(t)dt ≈ x(t ) + 1 ∆t!^(t )q(t ); 1 0 2 t0 0 2 1 0 where !^ = [!; 0] is a quaternion extended from ! and !^(t1)q(t0) is the product of two quaternions. The 3D vector !(t) is the angular velocity. Its direction represents the angular velocity axis and its magnitude represents the rotation speed. The matrix I−1(t) is the body inertia. Like the mass, the inertia indicates the resistance of the body to rotational motion. But unlike the mass, the inertia 2 is not constant. Intuitively, when a body rotates, its resistance to the rotation around X, Y, Z axes will also change. So we must update the inertia as a function of time: I(t) = R(t)IbodyRT(t); (7) where R(t) is the rotational matrix calculated using the quaternion q(t) at time t, and Ibody is the inertia in the original configuration: body X T T I = (pi pi)1 − pipi : (8) i Here pi the original position of a vertex i and 1 is the 3-by-3 identity matrix. Now the only question is: how to evaluate the torque τ (t0) at time t0? Let fi(t) be the force applied at an individual vertex i, we calculate the torque as: X X τ (t0) = ri(t) × fi(t) = (R(t)pi) × fi(t); (9) i i where ri(t) is the vector from the center to vertex i, after the body is rotated as shown in Figure 1b. Note that we do not need to consider the gravity force in angular motion. This is because the gravity should not cause the body to spin. Mathematically, this implies the torques caused by the gravity are summed to zero. 3 Collision Handling To make rigid body simulation interesting, we must use interesting forces, specially the force caused by collision. To begin with, let us consider the status a vertex i shown in Figure 1c. Its position at time t is x(t) + ri(t), and its velocity at time t is: vt + !(t) × ri(t), where ri(t) = R(t)pi is the vector from the center to i. We claim that the vertex is in collision with the floor, if two conditions are satisfied: • a. Its position xi(t) is below the ground floor: y = 0. • b. Its velocity vi(t) is still going down. Let N = [0; 1; 0] be the upward vector. We formulate the two conditions as: x (t) · N = (x(t) + r (t)) · N < 0; i i (10) vi(t) · N = (vt + !(t) × ri(t)) · N < 0: By testing these two conditions, we can know if vertex i is in collision at time t. 3.1 Collision Response Once we find the collision, we need to remove it by applying collision response to the rigid body. There are many ways to calculate the collision response. Here we choose to use a collision impulse method, which gives more plausible results than other approaches, such as penalty forces. Our basic idea is to formulate the expected vertex velocity after collision handling. We know that before collision, the vertex is moving toward the ground floor: vi(t) · N < 0. So after collision new handling, we hope the vertex can move away from the ground floor: vi (t) · N = −µvi(t) · N > 0, in which µ 2 [0; 1) is called the restitution coefficient. Assuming that the velocity change happens only in the Y direction (friction free), we have: new new new vi (t) = v (t) + ! (t) × ri(t) = vi(t) − (µ + 1) (vi(t) · N) N: (11) 3 To achieve vnew(t) and !new(t), we assume that there is a sudden but unknown impulse j applied at vertex i and it causes the following changes to the velocities: vnew(t) = v(t) + j=M; new −1 (12) ! (t) = !(t) + I (t)(ri(t) × j) : By combining Equation 11 and Equation 12, we get: −1 j=M + I (t)(ri(t) × j) × ri(t) = −(µ + 1) (vi(t) · N) N: (13) To solve Equation 13, we propose to define the cross product by a matrix product: 2 3 0 −riz riy ∗ ri(t) × j = Ri j = 4 riz 0 −rix 5 j; (14) −riy rix 0 where ri(t) = [rix; riy; riz]. So we get: 1 Kj = 1 − R∗I−1R∗ j = −(µ + 1) (v (t) · N) N; (15) M i i i and we can then solve j using the inverse of the matrix K. Once we get j, we use it to update v(t) and !(t) respectively. That concludes our collision handling process. 3.2 More Details There are more implementation details. What if there are many vertices in collision? Just detect them all and use the average displacement to the center as ri(t) instead. Why oscillation? If you see oscillation when the body sits on the floor, this is because it keeps bouncing back and forth due to the restitution coefficient.
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]
  • 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.
    [Show full text]
  • 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
    [Show full text]
  • Newton Euler Equations of Motion Examples
    Newton Euler Equations Of Motion Examples Alto and onymous Antonino often interloping some obligations excursively or outstrikes sunward. Pasteboard and Sarmatia Kincaid never flits his redwood! Potatory and larboard Leighton never roller-skating otherwhile when Trip notarizes his counterproofs. Velocity thus resulting in the tumbling motion of rigid bodies. Equations of motion Euler-Lagrange Newton-Euler Equations of motion. Motion of examples of experiments that a random walker uses cookies. Forces by each other two examples of example are second kind, we will refer to specify any parameter in. 213 Translational and Rotational Equations of Motion. Robotics Lecture Dynamics. Independence from a thorough description and angular velocity as expected or tofollowa userdefined behaviour does it only be loaded geometry in an appropriate cuts in. An interface to derive a particular instance: divide and author provides a positive moment is to express to output side can be run at all previous step. The analysis of rotational motions which make necessary to decide whether rotations are. For xddot and whatnot in which a very much easier in which together or arena where to use them in two backwards operation complies with respect to rotations. Which influence of examples are true, is due to independent coordinates. On sameor adjacent joints at each moment equation is also be more specific white ellipses represent rotations are unconditionally stable, for motion break down direction. Unit quaternions or Euler parameters are known to be well suited for the. The angular momentum and time and runnable python code. The example will be run physics examples are models can be symbolic generator runs faster rotation kinetic energy.
    [Show full text]
  • Leonhard Euler - Wikipedia, the Free Encyclopedia Page 1 of 14
    Leonhard Euler - Wikipedia, the free encyclopedia Page 1 of 14 Leonhard Euler From Wikipedia, the free encyclopedia Leonhard Euler ( German pronunciation: [l]; English Leonhard Euler approximation, "Oiler" [1] 15 April 1707 – 18 September 1783) was a pioneering Swiss mathematician and physicist. He made important discoveries in fields as diverse as infinitesimal calculus and graph theory. He also introduced much of the modern mathematical terminology and notation, particularly for mathematical analysis, such as the notion of a mathematical function.[2] He is also renowned for his work in mechanics, fluid dynamics, optics, and astronomy. Euler spent most of his adult life in St. Petersburg, Russia, and in Berlin, Prussia. He is considered to be the preeminent mathematician of the 18th century, and one of the greatest of all time. He is also one of the most prolific mathematicians ever; his collected works fill 60–80 quarto volumes. [3] A statement attributed to Pierre-Simon Laplace expresses Euler's influence on mathematics: "Read Euler, read Euler, he is our teacher in all things," which has also been translated as "Read Portrait by Emanuel Handmann 1756(?) Euler, read Euler, he is the master of us all." [4] Born 15 April 1707 Euler was featured on the sixth series of the Swiss 10- Basel, Switzerland franc banknote and on numerous Swiss, German, and Died Russian postage stamps. The asteroid 2002 Euler was 18 September 1783 (aged 76) named in his honor. He is also commemorated by the [OS: 7 September 1783] Lutheran Church on their Calendar of Saints on 24 St. Petersburg, Russia May – he was a devout Christian (and believer in Residence Prussia, Russia biblical inerrancy) who wrote apologetics and argued Switzerland [5] forcefully against the prominent atheists of his time.
    [Show full text]
  • Lecture L27 - 3D Rigid Body Dynamics: Kinetic Energy; Instability; Equations of Motion
    J. Peraire, S. Widnall 16.07 Dynamics Fall 2008 Version 2.0 Lecture L27 - 3D Rigid Body Dynamics: Kinetic Energy; Instability; Equations of Motion 3D Rigid Body Dynamics In Lecture 25 and 26, we laid the foundation for our study of the three-dimensional dynamics of rigid bodies by: 1.) developing the framework for the description of changes in angular velocity due to a general motion of a three-dimensional rotating body; and 2.) developing the framework for the effects of the distribution of mass of a three-dimensional rotating body on its motion, defining the principal axes of a body, the inertia tensor, and how to change from one reference coordinate system to another. We now undertake the description of angular momentum, moments and motion of a general three-dimensional rotating body. We approach this very difficult general problem from two points of view. The first is to prescribe the motion in term of given rotations about fixed axes and solve for the force system required to sustain this motion. The second is to study the ”free” motions of a body in a simple force field such as gravitational force acting through the center of mass or ”free” motion such as occurs in a ”zero-g” environment. The typical problems in this second category involve gyroscopes and spinning tops. The second set of problems is by far the more difficult. Before we begin this general approach, we examine a case where kinetic energy can give us considerable insight into the behavior of a rotating body. This example has considerable practical importance and its neglect has been the cause of several system failures.
    [Show full text]
  • 8.09(F14) Chapter 2: Rigid Body Dynamics
    Chapter 2 Rigid Body Dynamics 2.1 Coordinates of a Rigid Body A set of N particles forms a rigid body if the distance between any 2 particles is fixed: rij ≡ jri − rjj = cij = constant: (2.1) Given these constraints, how many generalized coordinates are there? If we know 3 non-collinear points in the body, the remaing points are fully determined by triangulation. The first point has 3 coordinates for translation in 3 dimensions. The second point has 2 coordinates for spherical rotation about the first point, as r12 is fixed. The third point has one coordinate for circular rotation about the axis of r12, as r13 and r23 are fixed. Hence, there are 6 independent coordinates, as represented in Fig. 2.1. This result is independent of N, so this also applies to a continuous body (in the limit of N ! 1). Figure 2.1: 3 non-collinear points can be fully determined by using only 6 coordinates. Since the distances between any two other points are fixed in the rigid body, any other point of the body is fully determined by the distance to these 3 points. 29 CHAPTER 2. RIGID BODY DYNAMICS The translations of the body require three spatial coordinates. These translations can be taken from any fixed point in the body. Typically the fixed point is the center of mass (CM), defined as: 1 X R = m r ; (2.2) M i i i where mi is the mass of the i-th particle and ri the position of that particle with respect to a fixed origin and set of axes (which will notationally be unprimed) as in Fig.
    [Show full text]
  • 2D Rigid Body Dynamics: Work and Energy
    J. Peraire, S. Widnall 16.07 Dynamics Fall 2008 Version 1.0 Lecture L22 - 2D Rigid Body Dynamics: Work and Energy In this lecture, we will revisit the principle of work and energy introduced in lecture L11-13 for particle dynamics, and extend it to 2D rigid body dynamics. Kinetic Energy for a 2D Rigid Body We start by recalling the kinetic energy expression for a system of particles derived in lecture L11, n 1 2 X 1 02 T = mvG + mir_i ; 2 2 i=1 0 where n is the total number of particles, mi denotes the mass of particle i, and ri is the position vector of Pn particle i with respect to the center of mass, G. Also, m = i=1 mi is the total mass of the system, and vG is the velocity of the center of mass. The above expression states that the kinetic energy of a system of particles equals the kinetic energy of a particle of mass m moving with the velocity of the center of mass, plus the kinetic energy due to the motion of the particles relative to the center of mass, G. For a 2D rigid body, the velocity of all particles relative to the center of mass is a pure rotation. Thus, we can write 0 0 r_ i = ! × ri: Therefore, we have n n n X 1 02 X 1 0 0 X 1 02 2 mir_i = mi(! × ri) · (! × ri) = miri ! ; 2 2 2 i=1 i=1 i=1 0 Pn 02 where we have used the fact that ! and ri are perpendicular.
    [Show full text]
  • Rigid Body Dynamics (Course Slides), M Müller-Fischer 2005, ETHZ Zurich
    From Particles to Rigid Bodies • Particles • Rigid bodies – No rotations – 6 DoFs (translation + rotation) – Linear velocity v only – Linear velocity v – 3N DoFs – Angular velocity ω COMP768- M.Lin Outline • Rigid Body Representation • Kinematics • Dynamics • Simulation Algorithm • Collisions and Contact Response COMP768- M.Lin Coordinate Systems Body Space • Body Space (Local Coordinate System) – Rigid bodies are defined relative to this system – Center of mass is the origin (for convenience) • We will specify body-related physical properties (inertia, …) in this frame COMP768- M.Lin Coordinate Systems World Space • World Space: rigid body transformation to common frame translation COMP768- M.Lin rotation Center of mass • Definition • Motivation: forces (one mass particle:) Image ETHZ 2005 (entire body:) COMP768- M.Lin Rotations • Euler angles: – 3 DoFs: roll, pitch, heading – Dependent on order of application – Not practical Image ETHZ 2005 COMP768- M.Lin Rotations • Rotation matrix – 3x3 matrix: 9 DoFs – Columns: world-space coordinates of body- space base vectors – Rotate a vector: Image ETHZ 2005 COMP768- M.Lin Rotations • Problem with rotation matrices: numerical drift • Fix: use Gram-Schmidt orthogonalization • Drift is easier to fix with quaternions COMP768- M.Lin Unit Quaternion Definition • q = [s,v] : s is a scalar, v is vector • A rotation of θ about a unit axis u can be u represented by the unit quaternion: [cos(θ/2), sin(θ /2) u] θ • Rotate a vector: • Fix drift: – 4-tuple: vector representation of rotation – Normalized quaternion
    [Show full text]
  • Lecture 28 3D Rigid Body Dynamics: Equations of Motion
    J. Peraire, S. Widnall 16.07 Dynamics Fall 2009 Version 3.0 Lecture L28 - 3D Rigid Body Dynamics: Equations of Motion; Euler's Equations 3D Rigid Body Dynamics: Euler's Equations We now turn to the task of deriving the general equations of motion for a three-dimensional rigid body. These equations are referred to as Euler's equations. The governing equations are those of conservation of linear momentum L = MvG and angular momentum, H = [I]!, where we have written the moment of inertia in matrix form to remind us that in general the direction of the angular momentum is not in the direction of the rotation vector !. Conservation of linear momentum requires L_ = F (1) Conservation of angular momentum, about a fixed point O, requires _ H0 = M (2) or about the center of mass G _ HG = M G (3) In our previous application of these equations, we specified the motion and used the equations to specify what moments would be required to produce the prescribed motion. In this more general formulation, we allow the body to execute free motions, possibly under the action of external moments. We consider the general motion of a body about its center of mass, first examining this in a inertial reference frame. 1 At an instant of time, we can calculate the angular momentum of the body as H = [I]!. One possible method to obtain the moments and the motion of the body is to perform our analysis in this inertial coordinate system. We would of course align our coordinate system initially with the principal axes of the body.
    [Show full text]
  • Rigid-Body Dynamics with Friction and Impact∗
    SIAM REVIEW c 2000 Society for Industrial and Applied Mathematics Vol. 42, No. 1, pp. 3–39 Rigid-Body Dynamics with Friction and Impact∗ David E. Stewart† Abstract. Rigid-body dynamics with unilateral contact is a good approximation for a wide range of everyday phenomena, from the operation of car brakes to walking to rock slides. It is also of vital importance for simulating robots, virtual reality, and realistic animation. However, correctly modeling rigid-body dynamics with friction is difficult due to a number of dis- continuities in the behavior of rigid bodies and the discontinuities inherent in the Coulomb friction law. This is particularly crucial for handling situations with large coefficients of friction, which can result in paradoxical results known at least since Painlev´e[C. R. Acad. Sci. Paris, 121 (1895), pp. 112–115]. This single example has been a counterexample and cause of controversy ever since, and only recently have there been rigorous mathematical results that show the existence of solutions to his example. The new mathematical developments in rigid-body dynamics have come from several sources: “sweeping processes” and the measure differential inclusions of Moreau in the 1970s and 1980s, the variational inequality approaches of Duvaut and J.-L. Lions in the 1970s, and the use of complementarity problems to formulate frictional contact problems by L¨otstedt in the early 1980s. However, it wasn’t until much more recently that these tools were finally able to produce rigorous results about rigid-body dynamics with Coulomb friction and impulses. Key words. rigid-body dynamics, Coulomb friction, contact mechanics, measure-differential inclu- sions, complementarity problems AMS subject classifications.
    [Show full text]
  • Euler, Newton, and Foundations for Mechanics
    Euler, Newton, and Foundations for Mechanics Oxford Handbooks Online Euler, Newton, and Foundations for Mechanics Marius Stan The Oxford Handbook of Newton Edited by Eric Schliesser and Chris Smeenk Subject: Philosophy, History of Western Philosophy (Post-Classical) Online Publication Date: Mar 2017 DOI: 10.1093/oxfordhb/9780199930418.013.31 Abstract and Keywords This chapter looks at Euler’s relation to Newton, and at his role in the rise of Newtonian mechanics. It aims to give a sense of Newton’s complicated legacy for Enlightenment science and to point out that key “Newtonian” results are really due to Euler. The chapter begins with a historiographical overview of Euler’s complicated relation to Newton. Then it presents an instance of Euler extending broadly Newtonian notions to a field he created: rigid-body dynamics. Finally, it outlines three open questions about Euler’s mechanical foundations and their proximity to Newtonianism: Is his mechanics a monolithic account? What theory of matter is compatible with his dynamical laws? What were his dynamical laws, really? Keywords: Newton, Euler, Newtonian, Enlightenment science, rigid-body dynamics, classical mechanics, Second Law, matter theory Though the Principia was an immense breakthrough, it was also an unfinished work in many ways. Newton in Book III had outlined a number of programs that were left for posterity to perfect, correct, and complete. And, after his treatise reached the shores of Europe, it was not clear to anyone how its concepts and laws might apply beyond gravity, which Newton had treated with great success. In fact, no one was sure that they extend to all mechanical phenomena.
    [Show full text]