Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Similar documents
Lagrangian and Hamiltonian Mechanics

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Lagrangian and Hamiltonian Dynamics

The Principle of Least Action

Chapter 2 Lagrangian Modeling

and from it produce the action integral whose variation we set to zero:

Table of Common Derivatives By David Abraham

Implicit Differentiation

Calculus of Variations

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

221A Lecture Notes Notes on Classica Mechanics I

Lecture 27: Generalized Coordinates and Lagrange s Equations of Motion

Introduction to the Vlasov-Poisson system

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

G j dq i + G j. q i. = a jt. and

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Linear First-Order Equations

Math 342 Partial Differential Equations «Viktor Grigoryan

Calculus of Variations

4. Important theorems in quantum mechanics

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

MATH2231-Differentiation (2)

u t v t v t c a u t b a v t u t v t b a

Euler equations for multiple integrals

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v

The Ehrenfest Theorems

Vectors in two dimensions

Calculus of variations - Lecture 11

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Quantum mechanical approaches to the virial

The Exact Form and General Integrating Factors

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

Short Intro to Coordinate Transformation

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Introduction to variational calculus: Lecture notes 1

arxiv: v1 [physics.class-ph] 20 Dec 2017

Diagonalization of Matrices Dr. E. Jacobs

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere

Topic 2.3: The Geometry of Derivatives of Vector Functions

Schrödinger s equation.

The Principle of Least Action and Designing Fiber Optics

Multivariable Calculus: Chapter 13: Topic Guide and Formulas (pgs ) * line segment notation above a variable indicates vector

Quantum Mechanics in Three Dimensions

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

The continuity equation

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Goldstein Chapter 1 Exercises

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Numerical Integrator. Graphics

Statics. There are four fundamental quantities which occur in mechanics:

Separation of Variables

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

THERE ARE BASICALLY TWO approaches to conservative mechanics, due principally

Lecture XVI: Symmetrical spacetimes

11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

θ x = f ( x,t) could be written as

Math 1271 Solutions for Fall 2005 Final Exam

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

1 Heisenberg Representation

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Math 115 Section 018 Course Note

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Gyroscopic matrices of the right beams and the discs

Tutorial Test 5 2D welding robot

GEODESIC BOUNDARY VALUE PROBLEMS WITH SYMMETRY. Colin J. Cotter. Darryl D. Holm. (Communicated by the associate editor name)

Mathematical Review Problems

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

1 M3-4-5A16 Assessed Problems # 1: Do 4 out of 5 problems

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

2.20 Marine Hydrodynamics Lecture 3

Linear and quadratic approximation

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask

Conservation Laws. Chapter Conservation of Energy

Physics 115C Homework 4

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

Lie symmetry and Mei conservation law of continuum system

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

Partial Differential Equations

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs

Discrete Hamilton Jacobi Theory and Discrete Optimal Control

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

Basic Differentiation Rules and Rates of Change. The Constant Rule

Physics 170 Week 7, Lecture 2

1 Lecture 20: Implicit differentiation

Pure Further Mathematics 1. Revision Notes

Calculus in the AP Physics C Course The Derivative

arxiv: v2 [math.ds] 26 Apr 2017

6 General properties of an autonomous system of two first order ODE

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

6 Wave equation in spherical polar coordinates

Transcription:

Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0, the initial positions an velocities r i t 0 ) an ṙ i t 0 ) for all point masses forming the system. Another formulation for the problem coul be to only give the positions of all point masses in the system, at the initial time t an final time t of the motion uner consieration. To see why this is mathematically equivalent, consier first t such that t = t t 0. Specifying r i t ) an r i t ) is then tantamount to specifying r i t ) an ṙ i t ). By continuity, for t sufficiently small, the motion remains unequivocally etermine if one specifies r i t ) an r i t ). As t grows arbitrarily, however, one can fin several possible motions with the same en points. The principle of stationary action, sometimes also calle less accurately the principle of least action, says that among all possible paths from r i t ) an r i t ), the physically realizable paths are the paths that extremize a functional calle the action S efine as follows Sq,..., q N ] = t Lq t),..., q N t), q t,..., q N, t) ) where L is calle the Lagrangian of the system. Note that the functions q i t) are generalize coorinates; they o not necessarily have to be Cartesian coorinates or other stanar coorinates. q i oes not have to have units of meter, an q i / oes not have to have the units of velocity. What the principle of stationary action is The principle of stationary action can be seen as a reformulation of the question of fining the proper evolution of a system in configuration space from a ifferential formulation to a variational/integral formulation. This has several avantages, which we briefly mentione in the last lecture, an which we will highlight as we start to use the variational formulation in practical examples. What the principle of stationary action is not The principle of stationary action oes not contain any new physics. All the physics is in fining appropriate Lagrangians o escribe a system of interest. We will soon learn a simple metho to construct Lagrangians for a large class of physical systems. However, there is no general metho to construct a Lagrangian for every system. In fact, in new fiels of physics, fining an appropriate Lagrangian to escribe the ynamics of a system can be a major accomplishment in its own right. Examples Before we learn how to construct Lagrangians for a certain class of physical systems, let us first see explicit variational formulations for the first two simple examples we looke at in the previous lecture. Newton s apple We saw that Newton s equations for Newton s apple are z = g ) The Newtonian ynamics for the apple can be written in a variational form as follows. A Lagrangian for the free-fall of an object suject to gravity is Lz, z, t) = ) z gz 3) where the mass of the object has been set to m = without any loss of generality. The action for the motion between times t an t thus is t Sz] = Lz, z t ) t, t) = z gz] t

As we have seen in our review of the calculus of variation, extremizing the action with fixe enpoints zt ) an zt ) leas to the following Euler-Lagrange equation for zt): Lz, z ), t) = Lz, z, t) 4) We have So Equation 4) can be rewritten as as expecte. The simple penulum A Lagrangian for the simple penulum is Minimizing the action Lz, z z, t) =, Lz, z, t) = g z = g Lθ, θ/, t) = ) θ l + gl cos θ 5) Sθ] = t t l ) θ + gl cos θ] leas to the E-L equations Lθ, θ/, t)) = Lθ, θ/, t) 6) We have Lθ, θ/, t) = l θ, Lθ, θ/, t) = gl sin θ so Equation 6) takes the form l g θ = gl sin θ θ + l sin θ = 0 which agrees with the secon-orer ODE we ha previously obtaine for that system. 3 Arbitrariness in the efinition of the Lagrangian Suppose a Lagrangian L escribes the ynamics of a system of point masses. Then, given the form of the E-L equation L) = L it is clear that the Lagrangian L = AL + B, where A an B are two constants, leas to the same equations of motion for the system. More interestingly, the equations of motion are also recovere if one replaces L with the Lagrangian L efine by L = L + ft, qt)) 7) Inee, if S is the action associate with the Lagrangian L an S is the action associate with the Lagrangian L, then Equation 7) implies t S q] = Sq] + t ft, qt)) = Sq] + ft, qt )) ft, qt )) 8) The principle of stationary action states that the physically relevant trajectory in configuration space is obtaine by extremals of the action holing the initial an final times an positions fixe. Since t, t, qt ) an qt ) are fixe, the extremals of S are the same as those of S, so both L an L can be use to escribe the motion of the system.

4 Hamilton s principle 4. Statement of the principle An important class of problems in classical mechanics is the class of problems that involves conservative forces only. Conservative forces are forces that can be expresse as the graient of a scalar function one often calls a potential, that epens only on time an the positions of the particles: F = V t, rt)). For these situations, the construction of a Lagrangian is straightforwar, an given by Hamilton s principle, which we state below: A system of point masses for which the forces are erive from a potential energy that is inepenent of velocity evolves along a path q for which the action Sq] = t t Lt, qt), q ) is stationary with respect to variations of the path q that leave the enpoints fixe, where L = T V is the ifference between the kinetic energy T an the potential energy V. 4. Illustration The examples provie in section are illustrations of Hamilton s principle. Let us stuy one more important example. Newton s law of gravitation tells us that two objects with masses m an M separate by a istance R exert a force with magnitue F G = G mm R on each other, an the force is an attractive force, irecte along the line joining the centers of gravity of the two objects. G 6.67 0 N.m.kg is calle the gravitational constant. We will consier the particular case in which M m, so the object of mass M is essentially at rest, an the lighter object, with mass m, orbits aroun the larger object. This is a very goo approximation for the motion of the Earth aroun the sun, or satellites, natural the Moon) or not, orbiting aroun Earth. It easy to show that the motion of the object of mass m is containe in a plane. If we take the center of gravity of the massive object as the origin of our coorinate system, we can write that the force acting on the light object is F G = GmM r r 3 Let us compute the quantity r r ṙ) = r r = r GM r 3 = 0 where the secon equality was obtaine by applying Newton s law. We conclue that Lt) = r ṙ = Cst = L0) = r0) ṙ0) 9) This means that the motion of the object of mass m is at all times in the plane perpenicular to L0). The analysis of the motion of the object can thus be reuce to stuying a system with two egrees of freeom. Let us consier the polar coorinate system r, θ) in the plane perpenicular to L0), with the origin fixe at the center of gravity of the massive object. In this coorinate system, the force F G can be written as F G = V where V = GmM 0) r is the gravitational potential. We can therefore use Hamilton s principle to write a Lagrangian for the object of mass m. The kinetic energy of the object is T = m r + r θ ) ] 3

so a Lagrangian for the object is Lr, θ, r, θ, t) = T V = m r + r θ ) ] + GmM r ) The equations for the motion of the object are given by the E-L equations. For the r coorinate, the E-L equation is 3Lr, θ, r, θ, t) = Lr, θ, r, θ ), t) r θ = r GM r ) The E-L equation for the θ coorinate is 4Lr, θ, r, θ, t) = Lr, θ, r, θ, t) r θ ) = 0 3) Equation 3) tells us that the quantity r θ is a conserve quantity of the motion. This quantity is calle the angular momentum. It is the magnitue of the vector L introuce in Equation 9). We will soon see that this quantity is conserve for all potentials V such that V = V r ), sometimes calle central potentials the associate forces are calle central forces). Let us write r θ = Γ = Cst Plugging Γ into Equation ), we have r = Γ r 3 GM r Multiplying this equation by r/, one fins ) r + Γ r GM r ] = 0 We conclue that ) r + Γ r GM r = Cst = E 4) The total energy of the system is conserve: the first two terms correspon to the kinetic energy, the last term is the potential energy. Equation 4) can be seen as the energy equation for a system with one egree of freeom, with an effective potential energy V eff r) = Γ r GM r This potential is ominate by the attractive gravitational part at large istances, but ominate by the repulsive angular momentum part at short istances. V eff r) is shown in Figure for the particular case Γ = GM =. Since 4) can be rewritten as ) r = E V eff r) point masses with E < 0 have boune orbits: there exists r in the attractive part of the potential such that E = V eff r ), an this r is a turning point of the trajectory. On the other han, point masses with E > 0 have unboune orbits. We will soon see that there are very efficient ways of fining out from the outset, from the functional epenence of the Lagrangian, that energy an angular momentum are conserve in the system stuie here. We coul have save a bit of algebra. Still, it oes not hurt too much to get practice with the umb ways at first. 4

.5 V eff r) 0.5 0 0.5 0 4 6 8 0 r Figure : V eff r) as a function of r for Γ = GM = A little igression regaring the classical -boy problem It turns out that any system consisting of two boies with masses m an m, at positions r t) an r t) an each exerting the force F = ±G m m r r 3 r r ) on the other boy can be stuie within the framework we presente above, even if the masses m an m are comparable an none of the boies can be consiere static. Here is the reason why. Newton s equations for the boies are Let us efine the position vectors m r = G m m r r 3 r r ) 5) m r = G m m r r 3 r r ) 6) R = m r + m r m + m r = r r R is the position of the center of mass of the system. Aing Equations 5) an 6), one fins the following equation for R: R = 0 7) This is expecte. The system of the two boies is isolate, so in the absence of external forces, the center of mass of the whole system has a constant velocity. Subtracting 6)/m from 5)/m, one fins a secon orer ODE for r: µ r = G m m r 3 r 8) with µ = m m /m + m ), which is often calle the reuce mass. Equation 8) escribes the motion of a particle of mass µ an position vector r with respect to a fixe center of force. This is precisely the situation we have stuie with the Lagrangian formulation above. 5

The separate motion of the two boies is obtaine by solving for r using either Newton s approach or a variational approach), an reconstructing r t) an r t) with the relations m r = R + r m + m m r = R r m + m 6