Classical MechanicsBSc · Notes

Coordinates, constraints, and virtual work

Symbols follow the course notation. A formula's conditions are part of the result. In particular, distinguish fixed-time virtual variations from actual time evolution, and distinguish remaining forces from forces already included in a potential.

Formula Meaning Symbols Conditions
fα(r1,…,rN,t)=0f_\alpha(\mathbf r_1,\ldots,\mathbf r_N,t)=0 Holonomic equality constraint α\alpha labels constraints A regular independent set is needed for simple counting
n=3N−kn=3N-k Configuration degrees of freedom NN particles, kk independent equalities Three-dimensional point particles; local regular constraints
ra=ra(q,t)\mathbf r_a=\mathbf r_a(q,t) Coordinate transformation aa particle, qiq_i independent coordinates A valid local coordinate patch
va=∑i∂qiraq˙i+∂tra\mathbf v_a=\sum_i\partial_{q_i}\mathbf r_a\dot q_i+\partial_t\mathbf r_a Actual velocity Last term is support transport Smooth transformation
δra=∑i∂qiraδqi\delta\mathbf r_a=\sum_i\partial_{q_i}\mathbf r_a\delta q_i Virtual displacement δt=0\delta t=0 Holonomic coordinates; admissible variation
Qi=∑aFa⋅∂qiraQ_i=\sum_a\mathbf F_a\cdot\partial_{q_i}\mathbf r_a Generalized force Forces appropriate to the stated work In Units 2–3, Fa\mathbf F_a here denotes remaining forces
δW=∑iQiδqi\delta W=\sum_iQ_i\delta q_i Virtual work QiδqiQ_i\delta q_i has energy units Same chosen force set
∑aRa⋅δra=0\sum_a\mathbf R_a\cdot\delta\mathbf r_a=0 Ideality of constraints Ra\mathbf R_a reactions Zero total virtual work for all admissible variations
∑aFa⋅δra=0\sum_a\mathbf F_a\cdot\delta\mathbf r_a=0 Virtual-work equilibrium Applied forces only Static equilibrium, ideal constraints; unilateral qualifications apply
∑a(Fa−maaa)⋅δra=0\sum_a(\mathbf F_a-m_a\mathbf a_a)\cdot\delta\mathbf r_a=0 D’Alembert principle aa\mathbf a_a physical acceleration Inertial frame, constant masses, ideal constraints

Lagrangian mechanics

Formula Meaning Symbols Conditions
L=T−VL=T-V Ordinary mechanical Lagrangian TT kinetic, VV potential energy Velocity-independent position potential in the main treatment
dLq˙i/dt−Lqi=QidL_{\dot q_i}/dt-L_{q_i}=Q_i Lagrange equation Subscripts denote partial derivatives Independent holonomic coordinates, ideal eliminated reactions
pi=Lq˙ip_i=L_{\dot q_i} Conjugate momentum pip_i belongs to qiq_i Defined by the chosen LL
Lqi=0, Qi=0⇒p˙i=0L_{q_i}=0,\ Q_i=0\Rightarrow\dot p_i=0 Cyclic-coordinate conservation qiq_i may still vary The remaining force must vanish
T=12∑ijMijq˙iq˙j+∑ibiq˙i+cT=\tfrac12\sum_{ij}M_{ij}\dot q_i\dot q_j+\sum_ib_i\dot q_i+c General kinetic-energy structure MijM_{ij} inertia matrix, bi,cb_i,c transport terms Ordinary particle kinetic energy and smooth transformations
Mij=∑ama∂qira⋅∂qjraM_{ij}=\sum_am_a\partial_{q_i}\mathbf r_a\cdot\partial_{q_j}\mathbf r_a Symmetric inertia coefficients Mij=MjiM_{ij}=M_{ji} Positive definite in regular independent coordinates with positive masses
∑iq˙iTq˙i=2T2+T1\sum_i\dot q_iT_{\dot q_i}=2T_2+T_1 Euler homogeneity identity TjT_j has velocity degree jj Coefficients held fixed in velocity derivatives
dTq˙i/dt−Tqi=∑amaaa⋅∂qiradT_{\dot q_i}/dt-T_{q_i}=\sum_am_a\mathbf a_a\cdot\partial_{q_i}\mathbf r_a Inertial projection identity Same coordinates as the transformation Constant masses, smooth mixed derivatives
E=∑iq˙ipi−L=T2−T0+V\mathcal E=\sum_i\dot q_ip_i-L=T_2-T_0+V Lagrangian energy function E\mathcal E need not be laboratory energy Last equality assumes L=T−V(q,t)L=T-V(q,t)
E˙=∑iQiq˙i−Lt\dot{\mathcal E}=\sum_iQ_i\dot q_i-L_t Energy-function balance LtL_t is explicit-time derivative Along solutions
R=12∑ijcijq˙iq˙j\mathcal R=\tfrac12\sum_{ij}c_{ij}\dot q_i\dot q_j Rayleigh function Symmetric damping matrix Linear viscous damping; nonnegative matrix for dissipation
Qi(d)=−Rq˙i,Pd=−2RQ_i^{(d)}=-\mathcal R_{\dot q_i},\quad P_d=-2\mathcal R Damping force and power PdP_d dissipative power Quadratic R\mathcal R; stationary coordinates for full physical power

Hamiltonian mechanics

Formula Meaning Symbols Conditions
H(q,p,t)=∑ipiq˙i−LH(q,p,t)=\sum_ip_i\dot q_i-L Legendre transformation Eliminate all dynamical q˙i\dot q_i Nonsingular velocity Hessian
q˙i=Hpi,p˙i=−Hqi\dot q_i=H_{p_i},\quad\dot p_i=-H_{q_i} Canonical Hamilton equations Other phase variables held fixed No remaining generalized forces
q˙i=Hpi,p˙i=−Hqi+Qi\dot q_i=H_{p_i},\quad\dot p_i=-H_{q_i}+Q_i Forced Hamilton equations QiQ_i expressed in q,p,tq,p,t Same regular underlying LL
H˙=Ht+∑iQiq˙i\dot H=H_t+\sum_iQ_i\dot q_i Hamiltonian balance Explicit-time plus force terms Along forced motion
H=T+VH=T+V Mechanical-energy identification Velocity replaced by momentum Sufficient: stationary transformations, T=T2T=T_2, VV velocity-independent
Hqi=0, Qi=0⇒p˙i=0H_{q_i}=0,\ Q_i=0\Rightarrow\dot p_i=0 Hamiltonian cyclicity Conserved conjugate momentum Unforced in that coordinate
H=pr2/(2μ)+pθ2/(2μr2)+V(r)H=p_r^2/(2\mu)+p_\theta^2/(2\mu r^2)+V(r) Polar central-force Hamiltonian pr=μr˙p_r=\mu\dot r, pθ=μr2θ˙p_\theta=\mu r^2\dot\theta Regular region r>0r>0

Central forces and radial motion

Formula Meaning Symbols Conditions
F=F(r)r^=−V′(r)r^\mathbf F=F(r)\hat{\mathbf r}=-V'(r)\hat{\mathbf r} Central interaction FF signed radial magnitude Time-independent radial position force
μ=m1m2/(m1+m2)\mu=m_1m_2/(m_1+m_2) Reduced mass m1,m2m_1,m_2 body masses Isolated two-body relative motion
L=r×p,L˙=0\mathbf L=\mathbf r\times\mathbf p,\quad\dot{\mathbf L}=0 Angular momentum p=μr˙\mathbf p=\mu\dot{\mathbf r} Radial force; constant mass
ℓ=μr2θ˙,S˙=ℓ/(2μ)\ell=\mu r^2\dot\theta,\quad\dot S=\ell/(2\mu) Angular momentum and area rate Choose positive orbital orientation Nonradial motion; area magnitude
ar=r¨−rθ˙2,aθ=rθ¨+2r˙θ˙a_r=\ddot r-r\dot\theta^2,\quad a_\theta=r\ddot\theta+2\dot r\dot\theta Polar acceleration components Radial and tangential basis Plane polar coordinates
Veff=V+ℓ2/(2μr2)V_{\rm eff}=V+\ell^2/(2\mu r^2) Effective radial potential Angular momentum fixed Central conservative problem
E=μr˙2/2+VeffE=\mu\dot r^2/2+V_{\rm eff} Radial energy equation EE relative energy Allowed radii satisfy E≥VeffE\geq V_{\rm eff}
Veff′(rc)=0V'_{\rm eff}(r_c)=0 Circular radius rcr_c constant radius Also E=Veff(rc)E=V_{\rm eff}(r_c)
ωr2=Veff′′(rc)/μ\omega_r^2=V''_{\rm eff}(r_c)/\mu Small radial oscillation frequency squared ωr\omega_r local frequency Positive: linear stability; negative: instability; zero: inconclusive
u′′+u=−μF(1/u)/(ℓ2u2)u''+u=-\mu F(1/u)/(\ell^2u^2) Binet orbit equation u=1/ru=1/r, prime d/dθd/d\theta ℓ≠0\ell\ne0, r>0r>0
E=ℓ2(u′2+u2)/(2μ)+V(1/u)E=\ell^2(u'^2+u^2)/(2\mu)+V(1/u) Energy in orbit variables Same uu convention Same nonradial assumptions
ωr2/Ω2=3+rcV′′/V′\omega_r^2/\Omega^2=3+r_cV''/V' Near-circular frequency ratio Ω=ℓ/(μrc2)\Omega=\ell/(\mu r_c^2) Stable circle; derivatives at rcr_c

Inverse-square orbits

In this table k>0k>0, V=−k/rV=-k/r, V(∞)=0V(\infty)=0, and ℓ>0\ell>0. Italic pp denotes semi-latus rectum; bold p\mathbf p denotes momentum.

Formula Meaning Symbols Conditions
r=p/(1+ecos⁡θ)r=p/(1+e\cos\theta) Conic orbit θ=0\theta=0 toward periapsis Positive denominator on the physical branch
p=ℓ2/(μk)p=\ell^2/(\mu k) Semi-latus rectum Length pp Attractive inverse-square orbit
e2=1+2Eℓ2/(μk2)e^2=1+2E\ell^2/(\mu k^2) Energy–eccentricity relation e≥0e\geq0 Excludes degenerate radial classification
rp=p/(1+e)r_p=p/(1+e) Periapsis radius rpr_p minimum radius Nonradial conic
ra=p/(1−e)r_a=p/(1-e) Apoapsis radius rar_a maximum radius Ellipse 0≤e<10\leq e<1
p=a(1−e2),b=a1−e2p=a(1-e^2),\quad b=a\sqrt{1-e^2} Ellipse geometry a,ba,b semiaxes Ellipse
E=−k/(2a)E=-k/(2a) Elliptic relative energy aa relative semimajor axis Bound orbit
v2=(k/μ)(2/r−1/a)v^2=(k/\mu)(2/r-1/a) Elliptic vis-viva vv relative speed For hyperbola with positive aha_h, replace −1/a-1/a by +1/ah+1/a_h
vc=k/(μr)v_c=\sqrt{k/(\mu r)} Circular relative speed Radius rr Exact circle
vesc=2k/(μr)v_{\rm esc}=\sqrt{2k/(\mu r)} Escape threshold speed Zero energy at infinity Ideal inverse-square model; collision restrictions separate
τ2=4π2μa3/k\tau^2=4\pi^2\mu a^3/k Kepler third law τ\tau period Relative ellipse
τ2=4π2a3/[G(m1+m2)]\tau^2=4\pi^2a^3/[G(m_1+m_2)] Gravitational form k=Gm1m2k=Gm_1m_2 Isolated Newtonian pair
A=p×L−μkr^\mathbf A=\mathbf p\times\mathbf L-\mu k\hat{\mathbf r} Runge–Lenz vector Fixed periapsis direction when nonzero Exact attractive inverse-square force
A⋅L=0,A2=μ2k2+2μEℓ2\mathbf A\cdot\mathbf L=0,\quad A^2=\mu^2k^2+2\mu E\ell^2 LRL constraints A=∣A∣A=|\mathbf A| Same model
e=A/(μk)\mathbf e=\mathbf A/(\mu k) Eccentricity vector Magnitude ee, periapsis direction Direction undefined for a circle
tan⁡(χ/2)=k/(μbimpv∞2)\tan(\chi/2)=k/(\mu b_{\rm imp}v_\infty^2) Attractive deflection angle Impact parameter and incoming speed Hyperbolic encounter, bimp>0b_{\rm imp}>0

The formula sheet is a reference after studying the chapters. The derivation index links each important argument to its assumptions and worked development.

Search the notes

Search all four units and the course references.