6.S: Lagrangian Dynamics (Summary)
( \newcommand{\kernel}{\mathrm{null}\,}\)
Newtonian plausibility argument for Lagrangian mechanics
A justification for introducing the calculus of variations to classical mechanics becomes apparent when the concept of the Lagrangian L≡T−U is used in the functional and time t is the independent variable. It was shown that Newton’s equation of motion can be rewritten as
ddt∂L∂˙qi−∂L∂qi=FEXqi
where FEXyi are the excluded forces of constraint plus any other conservative or non-conservative forces not included in the potential U. This corresponds to the Euler-Lagrange equation for determining the minimum of the time integral of the Lagrangian.
Equation ??? can be written as
ddt∂L∂˙qi−∂L∂qi=m∑kλk(t)∂gk∂qi+FEXCqi
where the Lagrange multiplier term accounts for holonomic constraint forces, and FEXCqi includes all additional forces not accounted for by the scalar potential U, or the Lagrange multiplier terms FHCqi. The constraint forces can be included explicitly as generalized forces in the excluded term \boldsymbol{F_{q_{i}}^{EXC}\} of Equation ???.
d’Alembert’s Principle
It was shown that d’Alembert’s Principle
N∑i(FAi−˙pi)⋅δri=0
cleverly transforms the principle of virtual work from the realm of statics to dynamics. Application of virtual work to statics primarily leads to algebraic equations between the forces, whereas d’Alembert’s principle applied to dynamics leads to differential equations.
Lagrange equations from d’Alembert’s Principle
After transforming to generalized coordinates, d’Alembert’s Principle leads to
N∑j[{ddt(∂T∂˙qj)−∂T∂qj}−Qj]δqj=0
If all the n coordinates qj are independent, then Equation ??? implies that the term in the square brackets is zero for each individual value of j. That is, this implies the basic Euler-Lagrange equations of motion.
The handling of both conservative and non-conservative generalized forces Qj is best achieved by assuming that the generalized force Qj=∑niFAi⋅∂ˉri∂qj can be partitioned into a conservative velocity-independent term, that can be expressed in terms of the gradient of a scalar potential, −∇Ui, plus an excluded generalized force QEXj which contains the non-conservative, velocity-dependent, and all the constraint forces not explicitly included in the potential Uj. That is,
Qj=−∇Uj+QEXj
Inserting ??? into ???, and assuming that the potential U is velocity independent, allows ??? to be rewritten as
∑j[{ddt(∂(T−U)∂˙qj)−∂(T−U)∂qj}−QEXj]δqj=0
Expressed in terms of the standard Lagrangian L=T−U this gives
N∑j[{ddt(∂L∂˙qj)−∂L∂qj}−QEXj]δqj=0
Note that Equation ??? contains the basic Euler-Lagrange Equation ??? for the special case when U=0. In addition, note that if all the generalized coordinates are independent, then the square bracket terms are zero for each value of j, which leads to the n general Euler-Lagrange equations of motion
{ddt(∂L∂˙qj)−∂L∂qj}=QEXj
where n≥j≥1. Newtonian mechanics has trouble handling constraint forces because they lead to coupling of the degrees of freedom. Lagrangian mechanics is more powerful since it provides the following three ways to handle such correlated motion.
1) Minimal set of generalized coordinates
If the n coordinates qj are independent, then the square bracket equals zero for each value of j in Equation ???, which corresponds to Euler’s equation for each of the n independent coordinates. If the n generalized coordinates are coupled by m constraints, then the coordinates can be transformed to a minimal set of s=n−m independent coordinates which then can be solved by applying Equation ??? to the minimal set of s independent coordinates.
2) Lagrange multipliers approach
The Lagrangian method concentrates solely on active forces, completely ignoring all other internal forces. In Lagrangian mechanics the generalized forces, corresponding to each generalized coordinate, can be partitioned three ways
Qj=−∇U+m∑k=1λk∂gk∂qj(q,t)+QEXCj
where the velocity-independent conservative forces can be absorbed into a scalar potential U, the holonomic constraint forces can be handled using the Lagrange multiplier term ∑mk=1λk∂gk∂qj(q,t), and the remaining part of the active forces can be absorbed into the generalized force QEXCj. The scalar potential energy U is handled by absorbing it into the standard Lagrangian L=T−U. If the constraint forces are holonomic then these forces are easily and elegantly handled by use of Lagrange multipliers. All remaining forces, including dissipative forces, can be handled by including them explicitly in the the generalized force QEXCj.
Combining the above two equations gives
N∑j[{ddt(∂L∂˙qj)−∂L∂qj}−QEXCj−m∑k=1λk∂gk∂qj(q,t)]δqj=0
Use of the Lagrange multipliers to handle the m constraint forces ensures that all n infinitessimals δqj are independent implying that the expression in the square bracket must be zero for each of the n values of j. This leads to n Lagrange equations plus m constraint relations
{ddt(∂L∂˙qj)−∂L∂qj}=QEXCj+m∑k=1λk∂gk∂qj(q,t)
where j=1,2,3,…n.
3) Generalized forces approach
The two right-hand terms in ??? can be understood to be those forces acting on the system that are not absorbed into the scalar potential U component of the Lagrangian L. The Lagrange multiplier terms ∑mk=1λk∂gk∂qj(q,t) account for the holonomic forces of constraint that are not included in the conservative potential or in the generalized forces QEXCj. The generalized force
QEXCj=n∑iFAi⋅∂ri∂pj
is the sum of the components in the qj direction for all external forces that have not been taken into account by the scalar potential or the Lagrange multipliers. Thus the non-conservative generalized force QEXCj contains non-holonomic constraint forces, including dissipative forces such as drag or friction, that are not included in U, or used in the Lagrange multiplier terms to account for the holonomic constraint forces.
Applying the Euler-Lagrange equations in mechanics:
The optimal way to exploit Lagrangian mechanics is as follows:
- Select a set of independent generalized coordinates.
- Partition the active forces into three groups:
- Conservative one-body forces
- Holonomic constraint forces
- Generalized forces
- Minimize the number of generalized coordinates.
- Derive the Lagrangian
- Derive the equations of motion
Velocity-dependent Lorentz force:
Usually velocity-dependent forces are non-holonomic. However, electromagnetism is a special case where the velocity-dependent Lorentz force F=q(E+v×B) can be obtained from a velocity-dependent potential function U(q,.q,t). It was shown that the velocity-dependent potential
U=qΦ−qv⋅A
leads to the Lorentz force where Φ is the scalar electric potential and A the vector potential.
Time-dependent forces:
It was shown that time-dependent forces can lead to complicated motion having both stable regions and unstable regions of motion that can exhibit chaos.
Impulsive forces:
A generalized impulse ˜Qj can be derived for an instantaneous impulsive force from the time integral of the impulsive forces Pi given by equation (3.12.49) using the time integral of equation (7.2.13), that is Δpj=˜Qj=limτ→0∫t+τtQEXCjdτ≡limτ→0∫t+τt∑iF⋅i∂ri∂qjdτ=∑i˜P⋅i∂ri∂qj
Note that the generalized impulse ˜Qj can be a translational impulse ˜Pj with corresponding translational variable qj or an angular impulsive torque ˜Tj with corresponding angular variable ϕj.
Comparison of Newtonian and Lagrangian mechanics:
In contrast to Newtonian mechanics, which is based on knowing all the vector forces acting on a system, Lagrangian mechanics can derive the equations of motion using generalized coordinates without requiring knowledge of the constraint forces acting on the system. Lagrangian mechanics provides a remarkably powerful, and incredibly consistent, approach to solving for the equations of motion in classical mechanics which is especially powerful for handling systems that are subject to holonomic constraints.