Skip to content

Lagrangian Mechanics

sources:

  • text: “Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison Wesley.”
  • text: “Landau, L. D., & Lifshitz, E. M. (1976). Mechanics (3rd ed.). Butterworth-Heinemann.”

Attempt these from memory before reading.

  • Write down the Euler-Lagrange equation. Which term is the “generalised momentum”?
  • What does L=T−VL = T - V assume about the forces and the constraints?
  • State Hamilton’s principle. Is the action minimised or only stationary?
  • What is a cyclic coordinate, and what is conserved when you have one?
  • For a bead on a rotating hoop, when does a new stable equilibrium appear?

Motivation: why not just use F=maF = ma

Section titled “Motivation: why not just use F=maF = maF=ma”

Newton’s second law is a vector equation. For a system of NN particles in three dimensions that is 3N3N coupled second-order equations, and the forces of constraint — normal forces, tension, the normal reaction of a rail — all appear in them. Most of the effort in a Newtonian solution goes into computing forces you do not care about, only to eliminate them at the end.

The Lagrangian formulation removes them at the outset by working in generalised coordinates: the smallest set of numbers that specifies the configuration once the constraints are accounted for. A bead constrained to a hoop needs one number, not three. A double pendulum needs two, not six.

The price is that the equations become second-order in the generalised coordinates with no explicit force bookkeeping. The reward is threefold:

  • The constraint forces never appear.
  • The choice of coordinates is free, so you pick the ones that make the problem easy.
  • Symmetries become visible, and by Noether’s theorem each one hands you a conserved quantity for nothing.

The structure of the whole subject is a single scalar function — the Lagrangian — plus one equation.

The Lagrangian of a system is defined as

L(q1,…,qn,q˙1,…,q˙n,t)=T−VL(q_1, \ldots, q_n, \dot{q}_1, \ldots, \dot{q}_n, t) = T - V

Where TT is the kinetic energy and VV is the potential energy.

3.2 Derivation of the Euler-Lagrange Equation from D”Alembert’s Principle

Section titled “3.2 Derivation of the Euler-Lagrange Equation from D”Alembert’s Principle”

Theorem 3.1 (Euler-Lagrange from D’Alembert). The equations of motion for a holonomic system with ideal constraints are:

ddt(∂L∂q˙j)−∂L∂qj=0,j=1,…,n\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_j}\right) - \frac{\partial L}{\partial q_j} = 0, \quad j = 1, \ldots, n

Proof. Start from D’Alembert’s principle with only applied forces (ideal constraints):

∑i(Fi(app)−mir¨i)⋅δri=0\sum_i (\mathbf{F}_i^{(\mathrm{app})} - m_i\ddot{\mathbf{r}}_i) \cdot \delta\mathbf{r}_i = 0

Express the virtual displacement in terms of generalised coordinates:

δri=∑j∂ri∂qjδqj\delta\mathbf{r}_i = \sum_j \frac{\partial \mathbf{r}_i}{\partial q_j}\delta q_j

First term (applied forces). For a conservative system, Fi(app)=−∇iV\mathbf{F}_i^{(\mathrm{app})} = -\nabla_i V So:

∑iFi(app)⋅δri=−∑i∇iV⋅∑j∂ri∂qjδqj=−∑j∂V∂qjδqj\sum_i \mathbf{F}_i^{(\mathrm{app})} \cdot \delta\mathbf{r}_i = -\sum_i \nabla_i V \cdot \sum_j \frac{\partial \mathbf{r}_i}{\partial q_j}\delta q_j = -\sum_j \frac{\partial V}{\partial q_j}\delta q_j

Defining the generalised force Qj=∑iFi⋅∂ri∂qjQ_j = \sum_i \mathbf{F}_i \cdot \frac{\partial \mathbf{r}_i}{\partial q_j} For conservative forces Qj=−∂V/∂qjQ_j = -\partial V/\partial q_j.

Second term (inertia). Using ∂r˙i∂q˙j=∂ri∂qj\frac{\partial \dot{\mathbf{r}}_i}{\partial \dot{q}_j} = \frac{\partial \mathbf{r}_i}{\partial q_j} (which holds when ri=ri(q,t)\mathbf{r}_i = \mathbf{r}_i(q, t)):

∑imir¨i⋅δri=∑imir¨i⋅∑j∂ri∂qjδqj=∑j[∑imir¨i⋅∂ri∂qj]δqj\sum_i m_i\ddot{\mathbf{r}}_i \cdot \delta\mathbf{r}_i = \sum_i m_i\ddot{\mathbf{r}}_i \cdot \sum_j \frac{\partial \mathbf{r}_i}{\partial q_j}\delta q_j = \sum_j \left[\sum_i m_i\ddot{\mathbf{r}}_i \cdot \frac{\partial \mathbf{r}_i}{\partial q_j}\right]\delta q_j

Now:

∑imir¨i⋅∂ri∂qj=ddt(∑imir˙i⋅∂ri∂qj)−∑imir˙i⋅ddt∂ri∂qj\sum_i m_i\ddot{\mathbf{r}}_i \cdot \frac{\partial \mathbf{r}_i}{\partial q_j} = \frac{d}{dt}\left(\sum_i m_i\dot{\mathbf{r}}_i \cdot \frac{\partial \mathbf{r}_i}{\partial q_j}\right) - \sum_i m_i\dot{\mathbf{r}}_i \cdot \frac{d}{dt}\frac{\partial \mathbf{r}_i}{\partial q_j}

Using ddt∂ri∂qj=∂r˙i∂qj\frac{d}{dt}\frac{\partial \mathbf{r}_i}{\partial q_j} = \frac{\partial \dot{\mathbf{r}}_i}{\partial q_j} and ∂r˙i∂q˙j=∂ri∂qj\frac{\partial \dot{\mathbf{r}}_i}{\partial \dot{q}_j} = \frac{\partial \mathbf{r}_i}{\partial q_j}:

∑imir¨i⋅∂ri∂qj=ddt(∑imir˙i⋅∂r˙i∂q˙j)−∑imir˙i⋅∂r˙i∂qj\sum_i m_i\ddot{\mathbf{r}}_i \cdot \frac{\partial \mathbf{r}_i}{\partial q_j} = \frac{d}{dt}\left(\sum_i m_i\dot{\mathbf{r}}_i \cdot \frac{\partial \dot{\mathbf{r}}_i}{\partial \dot{q}_j}\right) - \sum_i m_i\dot{\mathbf{r}}_i \cdot \frac{\partial \dot{\mathbf{r}}_i}{\partial q_j} =ddt∂T∂q˙j−∂T∂qj= \frac{d}{dt}\frac{\partial T}{\partial \dot{q}_j} - \frac{\partial T}{\partial q_j}

Combining both terms in D’Alembert’s principle:

∑j[Qj−ddt∂T∂q˙j+∂T∂qj]δqj=0\sum_j \left[Q_j - \frac{d}{dt}\frac{\partial T}{\partial \dot{q}_j} + \frac{\partial T}{\partial q_j}\right]\delta q_j = 0

For conservative forces, Qj=−∂V/∂qjQ_j = -\partial V/\partial q_j. Since L=T−VL = T - V and VV is independent of q˙j\dot{q}_j:

∑j[∂L∂qj−ddt∂L∂q˙j]δqj=0\sum_j \left[\frac{\partial L}{\partial q_j} - \frac{d}{dt}\frac{\partial L}{\partial \dot{q}_j}\right]\delta q_j = 0

Since the δqj\delta q_j are independent (and we have nn degrees of freedom):

ddt(∂L∂q˙j)−∂L∂qj=0\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_j}\right) - \frac{\partial L}{\partial q_j} = 0

■\blacksquare

3.3 Hamilton’s Principle (Action Formulation)

Section titled “3.3 Hamilton’s Principle (Action Formulation)”

Theorem 3.2 (Hamilton’s Principle). The actual path of a system between times t1t_1 and t2t_2 is The one that makes the action

S=∫t1t2L(q,q˙,t) dtS = \int_{t_1}^{t_2} L(q, \dot{q}, t)\, dt

Stationary.

Theorem 3.3 (Euler-Lagrange Equation from Hamilton’s Principle). The path q(t)q(t) that makes SS stationary satisfies

ddt(∂L∂q˙j)−∂L∂qj=0,j=1,…,n\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_j}\right) - \frac{\partial L}{\partial q_j} = 0, \quad j = 1, \ldots, n

Proof (for one degree of freedom). Consider a variation q(t)+ϵη(t)q(t) + \epsilon \eta(t) where η(t1)=η(t2)=0\eta(t_1) = \eta(t_2) = 0. The variation of the action:

δS=∫t1t2(∂L∂qη+∂L∂q˙η˙)dt\delta S = \int_{t_1}^{t_2} \left(\frac{\partial L}{\partial q}\eta + \frac{\partial L}{\partial \dot{q}}\dot{\eta}\right) dt

Integrating the second term by parts:

δS=∫t1t2(∂L∂q−ddt∂L∂q˙)η dt+[∂L∂q˙η]t1t2\delta S = \int_{t_1}^{t_2} \left(\frac{\partial L}{\partial q} - \frac{d}{dt}\frac{\partial L}{\partial \dot{q}}\right) \eta\, dt + \left[\frac{\partial L}{\partial \dot{q}}\eta\right]_{t_1}^{t_2}

The boundary term vanishes since η(t1)=η(t2)=0\eta(t_1) = \eta(t_2) = 0. For δS=0\delta S = 0 for all η\eta By the fundamental lemma of the calculus of variations:

∂L∂q−ddt∂L∂q˙=0\frac{\partial L}{\partial q} - \frac{d}{dt}\frac{\partial L}{\partial \dot{q}} = 0

■\blacksquare

Intuition. Hamilton’s principle says nature is “lazy”: out of all possible paths connecting two configurations, the actual path taken is the one that makes the action stationary. This is a profound generalisation of Fermat’s principle of least time in optics.

Problem. Derive the equation of motion for a simple pendulum of length ll and mass mm.

Solution. Take θ\theta as the generalised coordinate. The position of the bob is (lsin⁡θ,−lcos⁡θ)(l\sin\theta, -l\cos\theta).

T=12m(x˙2+y˙2)=12ml2θ˙2T = \frac{1}{2}m(\dot{x}^2 + \dot{y}^2) = \frac{1}{2}ml^2\dot{\theta}^2 V=−mglcos⁡θV = -mgl\cos\theta L=T−V=12ml2θ˙2+mglcos⁡θL = T - V = \frac{1}{2}ml^2\dot{\theta}^2 + mgl\cos\theta

Euler-Lagrange equation:

∂L∂θ=−mglsin⁡θ,∂L∂θ˙=ml2θ˙\frac{\partial L}{\partial \theta} = -mgl\sin\theta, \quad \frac{\partial L}{\partial \dot{\theta}} = ml^2\dot{\theta} ddt(ml2θ˙)+mglsin⁡θ=0  ⟹  θ¨+glsin⁡θ=0\frac{d}{dt}(ml^2\dot{\theta}) + mgl\sin\theta = 0 \implies \ddot{\theta} + \frac{g}{l}\sin\theta = 0

For small angles (sin⁡θ≈θ\sin\theta \approx \theta): θ¨+glθ=0\ddot{\theta} + \frac{g}{l}\theta = 0 Giving simple Harmonic motion with ω=g/l\omega = \sqrt{g/l}. ■\blacksquare

Problem. Derive the equations of motion for a double pendulum: mass m1m_1 on rod l1l_1 Mass m2m_2 on rod l2l_2 attached to m1m_1.

Solution. Generalised coordinates: angles θ1,θ2\theta_1, \theta_2 from the vertical. Position of m1m_1:

x1=l1sin⁡θ1,y1=−l1cos⁡θ1x_1 = l_1\sin\theta_1, \quad y_1 = -l_1\cos\theta_1

Position of m2m_2:

x2=l1sin⁡θ1+l2sin⁡θ2,y2=−l1cos⁡θ1−l2cos⁡θ2x_2 = l_1\sin\theta_1 + l_2\sin\theta_2, \quad y_2 = -l_1\cos\theta_1 - l_2\cos\theta_2

Velocities:

x˙1=l1θ˙1cos⁡θ1,y˙1=l1θ˙1sin⁡θ1\dot{x}_1 = l_1\dot{\theta}_1\cos\theta_1, \quad \dot{y}_1 = l_1\dot{\theta}_1\sin\theta_1 x˙2=l1θ˙1cos⁡θ1+l2θ˙2cos⁡θ2,y˙2=l1θ˙1sin⁡θ1+l2θ˙2sin⁡θ2\dot{x}_2 = l_1\dot{\theta}_1\cos\theta_1 + l_2\dot{\theta}_2\cos\theta_2, \quad \dot{y}_2 = l_1\dot{\theta}_1\sin\theta_1 + l_2\dot{\theta}_2\sin\theta_2

Kinetic energy:

T=12m1(x˙12+y˙12)+12m2(x˙22+y˙22)T = \frac{1}{2}m_1(\dot{x}_1^2 + \dot{y}_1^2) + \frac{1}{2}m_2(\dot{x}_2^2 + \dot{y}_2^2) =12(m1+m2)l12θ˙12+12m2l22θ˙22+m2l1l2θ˙1θ˙2cos⁡(θ1−θ2)= \frac{1}{2}(m_1 + m_2)l_1^2\dot{\theta}_1^2 + \frac{1}{2}m_2 l_2^2\dot{\theta}_2^2 + m_2 l_1 l_2\dot{\theta}_1\dot{\theta}_2\cos(\theta_1 - \theta_2)

Potential energy:

V=−m1gl1cos⁡θ1−m2g(l1cos⁡θ1+l2cos⁡θ2)V = -m_1 g l_1\cos\theta_1 - m_2 g(l_1\cos\theta_1 + l_2\cos\theta_2)

The Euler-Lagrange equations for θ1\theta_1 and θ2\theta_2 yield two coupled second-order ODEs. For equal masses and lengths (m1=m2=mm_1 = m_2 = m, l1=l2=ll_1 = l_2 = l):

(m+m)l2θ¨1+ml2θ¨2cos⁡(θ1−θ2)+ml2θ˙22sin⁡(θ1−θ2)+2mglsin⁡θ1=0(m + m)l^2\ddot{\theta}_1 + ml^2\ddot{\theta}_2\cos(\theta_1 - \theta_2) + ml^2\dot{\theta}_2^2\sin(\theta_1 - \theta_2) + 2mgl\sin\theta_1 = 0 ml2θ¨2+ml2θ¨1cos⁡(θ1−θ2)−ml2θ˙12sin⁡(θ1−θ2)+mglsin⁡θ2=0ml^2\ddot{\theta}_2 + ml^2\ddot{\theta}_1\cos(\theta_1 - \theta_2) - ml^2\dot{\theta}_1^2\sin(\theta_1 - \theta_2) + mgl\sin\theta_2 = 0

■\blacksquare

Problem. Two masses m1m_1 and m2m_2 (m1>m2m_1 > m_2) are connected by a massless inextensible string over a frictionless pulley. Find the acceleration using the Lagrangian.

Solution

Choose the vertical displacement xx of m1m_1 (downward positive) as the generalised coordinate. Since the string is inextensible, m2m_2 moves up by xx.

T=12m1x˙2+12m2x˙2=12(m1+m2)x˙2T = \frac{1}{2}m_1\dot{x}^2 + \frac{1}{2}m_2\dot{x}^2 = \frac{1}{2}(m_1 + m_2)\dot{x}^2 V=−m1gx+m2gx=−(m1−m2)gxV = -m_1 g x + m_2 g x = -(m_1 - m_2)gx L=12(m1+m2)x˙2+(m1−m2)gxL = \frac{1}{2}(m_1 + m_2)\dot{x}^2 + (m_1 - m_2)gx

Euler-Lagrange equation:

∂L∂x=(m1−m2)g,∂L∂x˙=(m1+m2)x˙\frac{\partial L}{\partial x} = (m_1 - m_2)g, \quad \frac{\partial L}{\partial \dot{x}} = (m_1 + m_2)\dot{x} (m1+m2)x¨=(m1−m2)g(m_1 + m_2)\ddot{x} = (m_1 - m_2)g a=x¨=m1−m2m1+m2ga = \ddot{x} = \frac{m_1 - m_2}{m_1 + m_2}g

■\blacksquare

3.7 Worked Example: Bead on a Rotating Hoop

Section titled “3.7 Worked Example: Bead on a Rotating Hoop”

Problem. A bead of mass mm slides without friction on a circular hoop of radius RR. The hoop rotates about a vertical diameter with constant angular velocity ω\omega. Find the equilibrium positions and their stability.

Solution

Use the angle θ\theta from the bottom of the hoop as the generalised coordinate. The position of the bead in cylindrical coordinates (ρ,ϕ,z)(\rho, \phi, z):

ρ=Rsin⁡θ,ϕ=ωt,z=−Rcos⁡θ\rho = R\sin\theta, \quad \phi = \omega t, \quad z = -R\cos\theta

Velocity:

ρ˙=Rθ˙cos⁡θ,ϕ˙=ω,z˙=Rθ˙sin⁡θ\dot{\rho} = R\dot{\theta}\cos\theta, \quad \dot{\phi} = \omega, \quad \dot{z} = R\dot{\theta}\sin\theta

Kinetic energy:

T=12m(ρ˙2+ρ2ϕ˙2+z˙2)=12mR2θ˙2+12mR2ω2sin⁡2θT = \frac{1}{2}m(\dot{\rho}^2 + \rho^2\dot{\phi}^2 + \dot{z}^2) = \frac{1}{2}m R^2\dot{\theta}^2 + \frac{1}{2}mR^2\omega^2\sin^2\theta

Potential energy:

V=−mgRcos⁡θV = -mgR\cos\theta

Lagrangian:

L=12mR2θ˙2+12mR2ω2sin⁡2θ+mgRcos⁡θL = \frac{1}{2}mR^2\dot{\theta}^2 + \frac{1}{2}mR^2\omega^2\sin^2\theta + mgR\cos\theta

Euler-Lagrange equation:

mR2θ¨=mR2ω2sin⁡θcos⁡θ−mgRsin⁡θmR^2\ddot{\theta} = mR^2\omega^2\sin\theta\cos\theta - mgR\sin\theta θ¨=sin⁡θ(ω2cos⁡θ−gR)\ddot{\theta} = \sin\theta\left(\omega^2\cos\theta - \frac{g}{R}\right)

Equilibrium (θ¨=0\ddot{\theta} = 0, θ˙=0\dot{\theta} = 0): sin⁡θ=0\sin\theta = 0 giving θ=0\theta = 0 (bottom), or cos⁡θ=g/(Rω2)\cos\theta = g/(R\omega^2) which exists only when ω2>g/R\omega^2 \gt g/R.

For ω2<g/R\omega^2 \lt g/R: only θ=0\theta = 0 is stable. For ω2>g/R\omega^2 \gt g/R: the bottom becomes unstable and the new equilibria at cos⁡θ=g/(Rω2)\cos\theta = g/(R\omega^2) are stable.

■\blacksquare

3.8 Cyclic Coordinates and Conserved Quantities

Section titled “3.8 Cyclic Coordinates and Conserved Quantities”

Definition. A coordinate qjq_j is cyclic (or ignorable) if it does not appear explicitly in the Lagrangian: ∂L/∂qj=0\partial L / \partial q_j = 0.

Theorem 3.4. If qjq_j is cyclic, the conjugate generalised momentum pj=∂L/∂q˙jp_j = \partial L / \partial \dot{q}_j is a constant of motion.

Proof. The Euler-Lagrange equation for a cyclic coordinate is:

ddt(∂L∂q˙j)=0  ⟹  pj=const\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_j}\right) = 0 \implies p_j = \mathrm{const}

■\blacksquare

Intuition. Cyclic coordinates correspond to symmetries of the system. Each symmetry gives a conserved quantity --- this is the essence of Noether’s theorem (Section 5).

When holonomic constraints cannot be eliminated by coordinate choice, introduce Lagrange multipliers λa\lambda_a:

ddt∂L∂q˙j−∂L∂qj=∑aλa∂fa∂qj\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_j} - \frac{\partial L}{\partial q_j} = \sum_a \lambda_a \frac{\partial f_a}{\partial q_j}

The multipliers λa\lambda_a are proportional to the constraint forces.

3.10 The Energy Function and Its Relation to the Hamiltonian

Section titled “3.10 The Energy Function and Its Relation to the Hamiltonian”

Definition. The energy function (also called the Jacobi integral) is:

h=∑jq˙j∂L∂q˙j−Lh = \sum_j \dot{q}_j \frac{\partial L}{\partial \dot{q}_j} - L

Theorem 3.5. If LL does not depend explicitly on time, then hh is conserved. Furthermore, if the transformation ri=ri(q)\mathbf{r}_i = \mathbf{r}_i(q) does not depend explicitly on time and VV is velocity-independent, then h=T+Vh = T + V (the total energy).

Proof. Taking the total time derivative:

dhdt=∑jq¨j∂L∂q˙j+∑jq˙jddt∂L∂q˙j−∑j∂L∂qjq˙j−∑j∂L∂q˙jq¨j−∂L∂t\frac{dh}{dt} = \sum_j \ddot{q}_j \frac{\partial L}{\partial \dot{q}_j} + \sum_j \dot{q}_j \frac{d}{dt}\frac{\partial L}{\partial \dot{q}_j} - \sum_j \frac{\partial L}{\partial q_j}\dot{q}_j - \sum_j \frac{\partial L}{\partial \dot{q}_j}\ddot{q}_j - \frac{\partial L}{\partial t}

The q¨j\ddot{q}_j terms cancel. Using the Euler-Lagrange equation ddt∂L∂q˙j=∂L∂qj\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_j} = \frac{\partial L}{\partial q_j}:

dhdt=−∂L∂t\frac{dh}{dt} = -\frac{\partial L}{\partial t}

If ∂L/∂t=0\partial L/\partial t = 0 Then dh/dt=0dh/dt = 0.

For the second part, when ri=ri(q)\mathbf{r}_i = \mathbf{r}_i(q) (scleronomic) and V=V(q)V = V(q):

T=12∑i,j,kmi∂ri∂qj∂ri∂qkq˙jq˙kT = \frac{1}{2}\sum_{i,j,k} m_i \frac{\partial \mathbf{r}_i}{\partial q_j}\frac{\partial \mathbf{r}_i}{\partial q_k}\dot{q}_j\dot{q}_k

Is a homogeneous quadratic form in q˙j\dot{q}_j. By Euler’s theorem for homogeneous functions:

∑jq˙j∂T∂q˙j=2T\sum_j \dot{q}_j \frac{\partial T}{\partial \dot{q}_j} = 2T

Since ∂L/∂q˙j=∂T/∂q˙j\partial L/\partial \dot{q}_j = \partial T/\partial \dot{q}_j (as VV is velocity-independent):

h=∑jq˙j∂T∂q˙j−T+V=2T−T+V=T+Vh = \sum_j \dot{q}_j \frac{\partial T}{\partial \dot{q}_j} - T + V = 2T - T + V = T + V

■\blacksquare

The clearest single picture of what the formalism buys is the rotating hoop of worked example 3.7, whose equilibria bifurcate as the hoop speeds up.

stable equilibria of the bead
versus rotation rate omega
theta_eq
pi | * * *
| *
pi/2 *
| *
0 *-------* <- bottom stays stable below
+-------------------------------- omega
g/sqrt(R) omega^2 R > g

Below ω2R=g\omega^2 R = g the only equilibrium is the bottom of the hoop. Above it, the bottom becomes unstable and two new stable equilibria appear at cos⁡θ=g/(Rω2)\cos\theta = g/(R\omega^2) — a pitchfork bifurcation, and the same mechanism that keeps a rider upright on a rotating fairground cylinder.

Multiple representations of the same fact, for the same system:

  • Words. Spin the hoop fast enough and the bead prefers a position off the bottom.
  • Equation. θ¨=sin⁡θ(ω2cos⁡θ−g/R)\ddot\theta = \sin\theta(\omega^2\cos\theta - g/R).
  • Equilibria. sin⁡θ=0\sin\theta = 0 always; cos⁡θ=g/(Rω2)\cos\theta = g/(R\omega^2) when ω2R>g\omega^2 R > g.
  • Stability. The second derivative of the effective potential changes sign at ω2R=g\omega^2 R = g, which is why the bottom turns over.

The Lagrangian formulation is like giving nature a GPS route planner. Instead of tracking every force individually, you describe the landscape of kinetic and potential energy, and nature automatically picks the path that makes the action stationary. Think of a ball rolling down a hill: you do not need to compute each force component if you know the shape of the terrain. The Lagrangian L = T - V encodes this terrain, and the Euler-Lagrange equations are nature’s optimization algorithm. Cyclic coordinates are like flat stretches of the landscape where nothing changes in a particular direction, so the corresponding momentum stays constant. The whole framework is coordinate-independent, meaning you can describe the same physical system using any convenient set of variables without changing the physics.

ClaimTrue?Counterexample
L=T−VL = T - V always gives correct equationsnoa charged particle in an electromagnetic field needs L=T−qϕ+q v⋅AL = T - q\phi + q\,\mathbf{v}\cdot\mathbf{A}
hh always equals T+VT + Vnothe rotating hoop: rheonomic constraints make h≠T+Vh \neq T + V
The action is minimisednoit is stationary; a saddle point is enough
Constraint forces appear in the equationsnothey are eliminated by generalised coordinates
Cyclic coordinates always existnothe double pendulum has none

The electromagnetic Lagrangian is the one to remember. It shows that T−VT - V is not the definition of a Lagrangian but a special case, and that the general object is whatever scalar makes the action stationary.

Mix these with material from the generalised coordinates page, Noether’s theorem, and the ODE module — selecting the method is the part blocked practice never trains.

  1. (Constraints.) A particle slides on a sphere of radius RR under gravity. Write the Lagrangian, show the constraint is holonomic, and eliminate it. At what angle does the particle leave the sphere?
  2. (Noether.) Show that translational symmetry of LL implies conservation of total momentum, and rotational symmetry implies conservation of angular momentum, using $\frac{d}{dt}\left(\sum_j \dot q_j \frac{\partial L}{\partial \dot q_j}
    • L\right)$ for the generator of the symmetry.
  3. (ODEs.) Linearise the double pendulum equations about the equilibrium θ1=θ2=0\theta_1 = \theta_2 = 0, find the normal mode frequencies, and show they reduce to g/l\sqrt{g/l} and 2g/l\sqrt{2g/l} when the two bobs are identical.
  4. (Calculus of variations.) Show that the shortest path between two points on a cylinder is a helix, by writing the arc length as an action and applying the Euler-Lagrange equation.
  5. (Energy.) For the rotating hoop, compute hh, show it is conserved, and show explicitly that h≠T+Vh \neq T + V.

Mistake 1: Including constraint forces in the Lagrangian The Lagrangian formulation eliminates constraint forces (like normal forces or tension in a string) by choosing appropriate generalised coordinates. A common error is adding the constraint force back into the Euler-Lagrange equations. Constraint forces do no work in holonomic systems and should not appear in L=T−VL = T - V.

Mistake 2: Assuming L=T−VL = T - V always gives the correct equations The Lagrangian L=T−VL = T - V is valid when the potential VV is velocity-independent and the coordinate transformation is scleronomic (time-independent). For rheonomic constraints or velocity-dependent potentials, the correct Lagrangian may differ. Always verify the conditions under which L=T−VL = T - V applies.

Mistake 3: Forgetting that the action is stationary, not necessarily minimised Hamilton’s principle states the action is stationary (δS=0\delta S = 0), not necessarily minimised. The actual path can be a saddle point or even a local maximum of the action. The stationarity condition is what yields the Euler-Lagrange equations, and the second variation determines whether the path is a minimum.