Classical MechanicsBSc · Notes

1. Geometry inside the kinetic energy

The term mr2θ˙2/2mr^2\dot\theta^2/2 in polar coordinates is not an extra interaction. It is ordinary translational kinetic energy expressed in variables whose unit displacements correspond to different physical distances. Understanding this structure explains generalized momentum, the inertial coefficients in Lagrange's equations, and the conditions under which the energy function equals T+VT+V.

This chapter uses the chain rule and dot product from Unit 1. Coefficients may depend on coordinates and time, but are held fixed when differentiating with respect to generalized velocities.

2. Expansion into quadratic, linear, and constant parts

For each particle, write

r˙a=∑ieaiq˙i+wa,eai=∂ra∂qi,wa=∂ra∂t.\dot{\mathbf r}_a=\sum_i\mathbf e_{ai}\dot q_i+\mathbf w_a, \qquad \mathbf e_{ai}=\frac{\partial\mathbf r_a}{\partial q_i},\quad \mathbf w_a=\frac{\partial\mathbf r_a}{\partial t}.(10.1)

The vector eai\mathbf e_{ai} describes the positional response to coordinate qiq_i, and wa\mathbf w_a describes transport caused by the prescribed time dependence of the constraints. Squaring gives

r˙a2=∑ij(eai⋅eaj)q˙iq˙j+2∑i(eai⋅wa)q˙i+wa2.\dot{\mathbf r}_a^2=\sum_{ij}(\mathbf e_{ai}\cdot\mathbf e_{aj})\dot q_i\dot q_j +2\sum_i(\mathbf e_{ai}\cdot\mathbf w_a)\dot q_i+\mathbf w_a^2.(10.2)

Multiplying by ma/2m_a/2 and summing yields

T=T2+T1+T0=12∑ijMijq˙iq˙j+∑ibiq˙i+c,\boxed{T=T_2+T_1+T_0 =\frac12\sum_{ij}M_{ij}\dot q_i\dot q_j+\sum_i b_i\dot q_i+c,}(10.3)

where

Mij=∑amaeai⋅eaj,bi=∑amaeai⋅wa,c=12∑amawa2.M_{ij}=\sum_a m_a\mathbf e_{ai}\cdot\mathbf e_{aj},\quad b_i=\sum_a m_a\mathbf e_{ai}\cdot\mathbf w_a,\quad c=\frac12\sum_a m_a\mathbf w_a^2.(10.4)

The subscripts indicate degree in the velocities. T0T_0 does not mean zero kinetic energy: a bead can have q˙=0\dot q=0 relative to a moving support while still moving in the laboratory. If the transformation has no explicit time dependence, all wa\mathbf w_a vanish and T=T2T=T_2.

3. Symmetry and positivity of the inertia coefficients

Because the dot product is symmetric, Mij=MjiM_{ij}=M_{ji}. For any real numbers ziz_i,

∑ijMijzizj=∑ama∣∑ieaizi∣2≥0.\sum_{ij}M_{ij}z_iz_j=\sum_a m_a\left|\sum_i\mathbf e_{ai}z_i\right|^2\geq0.(10.5)

For positive masses and regular independent coordinates, a nonzero vector zz produces a nonzero displacement of at least one particle. The right-hand side is then strictly positive. Thus the matrix MM is positive definite in a regular coordinate patch and invertible there.

Differentiating the quadratic expression uses its symmetry:

∂T∂q˙k=12∑jMkjq˙j+12∑iMikq˙i+bk=∑jMkjq˙j+bk.\frac{\partial T}{\partial\dot q_k} =\frac12\sum_jM_{kj}\dot q_j+\frac12\sum_iM_{ik}\dot q_i+b_k =\sum_jM_{kj}\dot q_j+b_k.(10.6)

Hence ∂2T/∂q˙i∂q˙j=Mij\partial^2T/\partial\dot q_i\partial\dot q_j=M_{ij}. For a velocity-independent VV, this is also the velocity Hessian of LL. Invertibility will permit replacing velocities by momenta in Hamiltonian mechanics. A coordinate singularity, such as polar coordinates at r=0r=0, can make this representation singular even though the Cartesian mechanics is regular.

4. Euler's identity and the energy function

A homogeneous function of degree ss satisfies f(λq˙)=λsf(q˙)f(\lambda\dot q)=\lambda^s f(\dot q). Differentiate both sides with respect to λ\lambda at λ=1\lambda=1:

∑iq˙i∂f∂q˙i=sf.\sum_i\dot q_i\frac{\partial f}{\partial\dot q_i}=sf.(10.7)

Apply this separately to T2,T1,T0T_2,T_1,T_0:

∑iq˙i∂T∂q˙i=2T2+T1.\boxed{\sum_i\dot q_i\frac{\partial T}{\partial\dot q_i}=2T_2+T_1.}(10.8)

For L=T−V(q,t)L=T-V(q,t), pi=∂T/∂q˙ip_i=\partial T/\partial\dot q_i, so

E=∑iq˙ipi−L=2T2+T1−(T2+T1+T0−V)=T2−T0+V.\mathcal E=\sum_i\dot q_i p_i-L =2T_2+T_1-(T_2+T_1+T_0-V) =\boxed{T_2-T_0+V}.(10.9)

For time-independent transformations, T1=T0=0T_1=T_0=0 and E=T+V\mathcal E=T+V. Whether that quantity is conserved is a second question, decided by E˙=∑iQiq˙i−∂tL\dot{\mathcal E}=\sum_iQ_i\dot q_i-\partial_tL. The equality with mechanical energy and the constancy along a motion require different checks.

5. The kinetic-energy derivative identity

One further property explains why differentiating a scalar can recover the projection of acceleration. With constant masses,

∂T∂q˙i=∑amar˙a⋅eai.\frac{\partial T}{\partial\dot q_i}=\sum_a m_a\dot{\mathbf r}_a\cdot\mathbf e_{ai}.(10.10)

Take a total time derivative:

ddt∂T∂q˙i=∑ama(r¨a⋅eai+r˙a⋅deaidt).\frac d{dt}\frac{\partial T}{\partial\dot q_i} =\sum_a m_a\left(\ddot{\mathbf r}_a\cdot\mathbf e_{ai}+\dot{\mathbf r}_a\cdot\frac{d\mathbf e_{ai}}{dt}\right).(10.11)

From differentiating the transformation and commuting smooth mixed partial derivatives, ∂r˙a/∂qi=deai/dt\partial\dot{\mathbf r}_a/\partial q_i=d\mathbf e_{ai}/dt. Therefore

∂T∂qi=∑amar˙a⋅deaidt.\frac{\partial T}{\partial q_i}=\sum_a m_a\dot{\mathbf r}_a\cdot\frac{d\mathbf e_{ai}}{dt}.(10.12)

Subtract the equal velocity terms to obtain

ddt∂T∂q˙i−∂T∂qi=∑amar¨a⋅∂ra∂qi.\boxed{\frac d{dt}\frac{\partial T}{\partial\dot q_i}-\frac{\partial T}{\partial q_i} =\sum_a m_a\ddot{\mathbf r}_a\cdot\frac{\partial\mathbf r_a}{\partial q_i}.}(10.13)

This is a kinetic-energy identity, valid also for smooth explicitly time-dependent transformations. Its use to derive Lagrange's equations from D’Alembert's principle is placed in the optional appendix, respecting the syllabus distinction.

6. Applications

Example 1 · Polar kinetic energy

For x=rcos⁡θx=r\cos\theta, y=rsin⁡θy=r\sin\theta, one finds er=(cos⁡θ,sin⁡θ)\mathbf e_r=(\cos\theta,\sin\theta) and eθ=(−rsin⁡θ,rcos⁡θ)\mathbf e_\theta=(-r\sin\theta,r\cos\theta). Their dot products are 1,0,r21,0,r^2. Thus Mrr=mM_{rr}=m, Mrθ=0M_{r\theta}=0, Mθθ=mr2M_{\theta\theta}=mr^2, and

T=12mr˙2+12mr2θ˙2.T=\frac12m\dot r^2+\frac12mr^2\dot\theta^2.(10.14)

Here pr=mr˙p_r=m\dot r and pθ=mr2θ˙p_\theta=mr^2\dot\theta. Direct calculation gives r˙pr+θ˙pθ=mr˙2+mr2θ˙2=2T\dot r p_r+\dot\theta p_\theta=m\dot r^2+mr^2\dot\theta^2=2T. The angular momentum coefficient is a moment of inertia, not a mass. At r=0r=0, the angular coordinate ceases to distinguish positions.

Example 2 · A pendulum with a moving pivot

Let the pivot move horizontally along X(t)X(t). With x=X(t)+asin⁡θx=X(t)+a\sin\theta and y=−acos⁡θy=-a\cos\theta,

T=12ma2θ˙2+maX˙cos⁡θθ˙+12mX˙2.T=\frac12ma^2\dot\theta^2+ma\dot X\cos\theta\dot\theta+\frac12m\dot X^2.(10.15)

The three terms are T2,T1,T0T_2,T_1,T_0. With V=mga(1−cos⁡θ)V=mga(1-\cos\theta),

pθ=ma2θ˙+maX˙cos⁡θ,E=12ma2θ˙2−12mX˙2+V.p_\theta=ma^2\dot\theta+ma\dot X\cos\theta,\qquad \mathcal E=\frac12ma^2\dot\theta^2-\frac12m\dot X^2+V.(10.16)

This differs from laboratory T+VT+V by the linear transport term and twice the constant term. The momentum is not simply the bob's angular momentum relative to the pivot. Its extra term records the prescribed support motion in the laboratory kinetic energy.

Example 3 · A bead on a rotating radial wire

A straight horizontal wire through the origin rotates at prescribed angular speed Ω\Omega. Choose signed distance qq along it, so r=(qcos⁡Ωt,qsin⁡Ωt)\mathbf r=(q\cos\Omega t,q\sin\Omega t). Then eq⋅w=0\mathbf e_q\cdot\mathbf w=0, ∣w∣2=Ω2q2|\mathbf w|^2=\Omega^2q^2, and

T=12mq˙2+12mΩ2q2,V=0.T=\frac12m\dot q^2+\frac12m\Omega^2q^2,\qquad V=0.(10.17)

Lagrange's equation gives mq¨−mΩ2q=0m\ddot q-m\Omega^2q=0. The energy function is E=mq˙2/2−mΩ2q2/2\mathcal E=m\dot q^2/2-m\Omega^2q^2/2 and is conserved because this reduced LL is autonomous. Laboratory kinetic energy is not conserved: the motor sustaining the rotation supplies work through the wire reaction. Thus autonomy of a reduced Lagrangian need not imply conservation of laboratory mechanical energy.

Example 4 · Mixed kinetic terms

Consider T=12m(q˙12+2αq˙1q˙2+q˙22)T=\tfrac12m(\dot q_1^2+2\alpha\dot q_1\dot q_2+\dot q_2^2), where both coordinates have dimensions of length and ∣α∣<1|\alpha|<1. The coefficient matrix is m(1αα1)m\begin{pmatrix}1&\alpha\\\alpha&1\end{pmatrix}, with determinant m2(1−α2)>0m^2(1-\alpha^2)>0.

The momenta are p1=m(q˙1+αq˙2)p_1=m(\dot q_1+\alpha\dot q_2) and p2=m(αq˙1+q˙2)p_2=m(\alpha\dot q_1+\dot q_2). Multiplying each by its velocity and adding gives 2T2T, including both cross terms. Solving the two equations gives

q˙1=p1−αp2m(1−α2),q˙2=p2−αp1m(1−α2).\dot q_1=\frac{p_1-\alpha p_2}{m(1-\alpha^2)},\qquad \dot q_2=\frac{p_2-\alpha p_1}{m(1-\alpha^2)}.(10.18)

When α=1\alpha=1, only q˙1+q˙2\dot q_1+\dot q_2 affects TT, so the assumed two-velocity description is singular. The restriction is mathematical information about independence, not an arbitrary algebraic inconvenience.

7. Conditions behind the familiar formulas

Before writing ∑ipiq˙i=2T\sum_i p_i\dot q_i=2T, check that the entire kinetic energy is quadratic homogeneous in the generalized velocities. Moving constraints can add lower-degree terms. Before replacing E\mathcal E by T+VT+V, check the same condition and the absence of velocity dependence in VV. These checks will be central to constructing and interpreting the Hamiltonian.

Exercises

Hints and solutions are collected separately.

  1. 10.1 Why can TT be nonzero when every q˙i=0\dot q_i=0?

  2. 10.2 Derive the coefficients Mij,bi,cM_{ij},b_i,c in the kinetic-energy expansion.

  3. 10.3 Prove that the regular inertia matrix is positive definite.

  4. 10.4 For T=Aq˙2+Bq˙+CT=A\dot q^2+B\dot q+C, calculate q˙ ∂T/∂q˙\dot q\,\partial T/\partial\dot q and compare it with 2T2T.

  5. 10.5 Prove E=T2−T0+V\mathcal E=T_2-T_0+V for velocity-independent VV.

  6. 10.6 For a bead on r=(q,Ut)\mathbf r=(q,Ut) with constant UU, find T,p,HT,p,H when V=0V=0.

  7. 10.7 What fails in polar coordinates at r=0r=0?

  8. 10.8 Calculate p1,p2p_1,p_2 for T=12m(2q˙12+2q˙1q˙2+3q˙22)T=\tfrac12m(2\dot q_1^2+2\dot q_1\dot q_2+3\dot q_2^2).

  9. 10.9 Derive the kinetic-energy identity relating its derivatives to acceleration projection.

  10. 10.10 For a pivot moving as X=UtX=Ut, find the term linear in pendulum angular velocity.

  11. 10.11 Does a symmetric inertia matrix have equal diagonal elements?

  12. 10.12 For L=12mq˙2+B(q)q˙−V(q)L=\tfrac12m\dot q^2+B(q)\dot q-V(q), show that B(q)B(q) cancels from the equation of motion.

Reading progress is saved on this device.

Search the notes

Search all four units and the course references.