$\require{cancel}$ $\require{physics}$ $\require{mathtools}$

Mechanics Notes

January 2021 • Notes

Table of Contents

  1. Lagrangians Overview
  2. Hamiltonians
  3. Oscillations
  4. The Central Force Problem
  5. Scattering Theory
  6. Rigid Body Motion
  7. Continuous Media

Lagrangians Overview

Assuming Newton’s second, we derive the euler lagrange equations

\begin{equation} m\ddot{\vb{r}}_\alpha = -\nabla V\eval{}_{\vb{r}_\alpha} + \vb{F}_\alpha^{\text{n-c}} \implies \dv{t}\pdv{L}{\dot{q}}-\pdv{L}{q} = \sum\pdv{\vb{r}_\alpha}{q}\cdot \vb{F}_\alpha^\text{n-c} \end{equation}

where we absorb the conservative forces into $L$. Some nonconservative forces can also be written as a velocity-dependent potential; an example is the Rayleigh dissipation function

\begin{equation} \mathcal{F} = \frac{1}{2}\sum(k_x v_{ix}^2 + k_y v_{iy}^2 + k_z v_{iz}^2) \end{equation}

for which the force is given by

\begin{equation} \pdv{\vb{r}_\alpha}{q}\cdot \vb{F}_\alpha^\text{n-c} = -\pdv{\mathcal{F}}{\dot{q}_j} \end{equation}

if there is a dissipation function $\mathcal{F}$, we have

\begin{equation} \dv{h}{t} = - 2\mathcal{F} - \pdv{L}{t} \end{equation}

as the rate of energy loss.

Constraints and Lagrange Multipliers

Fun vocabulary:

  • A holonomic constraint is one that can be written $f(\vb{r}_1,\, \vb{r}_2,\,\dots;\;t) = 0$
  • A rheonomous constraint is one that contains a dependence on time. Otherwise, the constraint is called scleronomous
  • A semiholonomic constraint is one that can be written in the form $f(q,\, \dot{q}) = 0$ or, often, $a_{i} \dd{q_i} + a_{t} \dd{t} = 0$
    •


    We include these constraints using Lagrange multipliers:

    1. Since the constraint equations $f_\alpha$ are identically zero, we can say
      \begin{equation} \delta\int(L-\lambda^\alpha(t) f_\alpha) \dd{t} = 0 \end{equation}
    2. Solving the variation gives the nonhomogeneous Euler-Lagrange equations for forces of constraint
      \begin{equation} \dv{t}\pdv{L}{\dot{q}}-\pdv{L}{q} = \sum \lambda^\alpha \qty[\pdv{f_\alpha}{q} - \dv{t}\pdv{f_\alpha}{\dot{q}}] - \dv{\lambda^\alpha}{t} \pdv{f_\alpha}{\dot{q}} \end{equation}
      Note that this right hand side simplifies a lot if the constraint is holonomic

    Conserved quantities

    cyclic coords (lagrangian does not explictly depend on them) immediately give conserved quantities, and typically energy is conserved as well.

    Noethers theorem says if

    \begin{equation} q_i\to q_i+\delta q_i \implies \delta L = \dv{f}{t} \end{equation}

    for some $f(q,\,t)$, then $Q=p_i\delta q^i - f$ is constant when the EoM is followed

    Return to Table of Contents

    Hamiltonians

    \begin{equation} H(p,\,q) = p_i\dot{q}^i - L \end{equation}

    with equations of motion

    \begin{equation} \pdv{H}{q}=-\dot{p}\qq{and}\pdv{H}{p}=\dot{q} \end{equation}

    Return to Table of Contents

    Oscillations

    One dimension: driving and damping

    Consider the EoM

    \begin{equation} m\ddot{x} = -kx -\gamma\dot{x}+F_D(t) \end{equation}

    define the undamped frequency and the quality factor

    \begin{equation} \omega_0 = \sqrt{\frac{k}{m}} \qquad Q = \frac{\sqrt{mk}}{\gamma} \implies \ddot{x} + \frac{\omega_0}{Q}\dot{x} + \omega_0^2x = \frac{1}{m}F_D \end{equation}

    if we take $x\sim e^{rt}$ and solve the homogeneous, we get the characteristic equation

    \begin{equation} r^{2}+\frac{\omega_0}{Q}r + \omega_0^2 = 0\implies r_{\pm} = \frac{\omega_0}{2Q}\qty[-1\pm\sqrt{1-4Q^2}] \end{equation}

    which has three cases:

  • $Q> 1/2$ "underdamped" — decaying packet times an oscillatory solution
    • \begin{equation} x(t) = e^{-\omega_0t/2Q} \cos(\omega_0\sqrt{1-1/4Q^2} + \phi) \end{equation}
  • $Q=1/2$ "critically damped"
    • \begin{equation} x(t) = [x_0 + (v_0+x_0\omega_0)t]e^{-\omega_0t} \end{equation}
    you can shove it hard enough to overshoot the equilibrium point but it'll approach it exponentially from there
  • $Q> 1/2$ "overdamped" — real and distinct roots give a lincomb
    • \begin{equation} x(t) = Ae^{r_+t} + Be^{r_-t} \end{equation}


    To get a particular solution in the case of a driving force, generally nice to solve in the complex plane. $x = \Re[z(t)]$ for $z(t) = Ae^{-i\omega t}$ where $\omega$ is the driving frequency and $A$ is in general complex. We can pull this out into a real amplitude and a phase lag. For a periodic driving force, we get resonance: $|A|$ is big when $\omega\sim\omega_0$ and $Q$ large.

    \begin{equation} |A| \text{ maximized when }\omega = \omega_0\sqrt{1-1/2Q^2} \end{equation}

    in $A-\omega$ space this gives a lorentzian with width $\Gamma = \omega/Q$ that looks like

    \begin{equation} A(\omega) \sim \frac{\Gamma^2/4}{(\omega-\omega_0)^2+\Gamma^2/4}E_\text{max} \end{equation}

    Higher dimensions

    If we have a Lagrangian of the form

    \begin{equation} L = \frac{1}{2}T_{ij}\dot{q}_i\dot{q}_j - \frac{1}{2}V_{ij}q_iq_j + \mathcal{O}(q^3) = \frac{1}{2}\dot{\vb{q}}^T\mathbb{T}\dot{\vb{q}} - \frac{1}{2}\vb{q}^T\mathbb{V}\vb{q} \end{equation}

    then the equations of motion read

    \begin{equation} T_{ij}\ddot{q}_j + V_{ij}q_j = 0 \implies q_j = v_je^{i\omega t} \implies (-\omega^2 \mathbb{T} + \mathbb{V})\vb{v} = 0 \end{equation}

    so we find valid $\omega^2$s by $\det(-\omega^2\mathbb{T}-\mathbb{V}) = 0$, and the sign determines the time dependence of $q$:

  • if $\omega^2_\alpha > 0$, $\vb{q} = \vb{v}_\alpha\qty[A_\alpha\cos(\omega t) + B_\alpha\sin(\omega t)]$
  • if $\omega^2_\alpha = 0$, $\vb{q} = \vb{v}_\alpha\qty[A_\alpha t + B_\alpha]$
  • if $\omega^2_\alpha < 0$, $\vb{q} = \vb{v}_\alpha\qty[A_\alpha e^{|\omega_\alpha|t} + B_\alpha e^{-|\omega_\alpha|t}]$
    •

    and the different $\vb{v}_\alpha$ are $\mathbb{T}$-orthornomal, meaning

    \begin{equation} \vb{v}_\alpha^T\mathbb{T}\vb{v}_\beta = \delta_{\alpha\beta} \end{equation}

    so in a 2D system knowing one determines the other

    Return to Table of Contents

    The Central Force Problem

    \begin{equation} L = \frac{1}{2}m\dot{r}^2 + \frac{1}{2}m(r\dot{\phi})^2 - V(r) \end{equation}

    can transform to the hammy or to the routhian ($R= p_\phi \dot{\phi}-L$, but we leave $r$ as in the lagrangian formalism)

    \begin{equation} H = \frac{p_r^2}{2m} + \underbrace{\frac{p_\phi^2}{2mr^2} + V(r)}_{V_\text{eff}} \end{equation}

    we usually care about $r(\phi)$, so we convert the time derivatives to phi derivatives:

    \begin{equation} \dot{r} = \dv{r}{\phi}\dot{\phi}=r'\frac{p_\phi}{mr^2}\implies E = \frac{p_\phi^2}{2m}\frac{r'^2}{r^4} + V_\text{eff} \end{equation}

    note it’s also often nice to swap variables to $u = 1/r$

    The Kepler problem

    \begin{align} E &= \frac{p_\phi^2}{2m}\frac{r'^2}{r^4} +\frac{p_\phi^2}{2mr^2} -\frac{k}{r} \\ &= \frac{p_\phi^2}{2m}u'^2 +\frac{p_\phi^2}{2m}u^2 -ku\\ &= \frac{p_\phi^2}{2m}\qty[u'^2 + \qty(u- \frac{km}{p_\phi^2})^2] - \frac{mk^2}{2p_\phi^2} \end{align}

    which is an offset SHO! customary names:

    \begin{gather} u_0 = \frac{km}{p_\phi^2} = \frac{1}{p} \qquad \varepsilon = Ap \implies r(\phi) = \frac{p}{1+\varepsilon\cos(\phi-\phi_0)}\\ E = \frac{p_\phi^2}{2m}A^2 - \frac{mk^2}{2p_\phi^2} = \frac{k}{2p}(\varepsilon^2-1) \end{gather}

    geometry of an ellipse:

  • $a$: semimajor axis, $2a= r_\text{min} + r_\text{max} = 2p/(1-\epsilon^2)$
  • $b$: semiminor axis, hard to find explicitly so use pythagorean thm
  • $c$: focus to center distance: $2c=r_\text{max}-r_\text{min}=2p\varepsilon/(1-\epsilon^2)=2\varepsilon a$
  • $r_\text{min} = p/(1+\varepsilon)$
  • $r_\text{max} = p/(1-\varepsilon)$
    •

    The Laplace-Runge-Lenz vector

    For a general central force, where $\dot{\vb{p}} = f(r)\frac{\vb{r}}{r}$, we can show

    \begin{equation} \dv{t}\qty(\vb{p}\times\vb{L}) = -mr^2\,f(r) \dv{t}(\frac{\vb{r}}{r}) \end{equation}

    in the Kepler problem, where $f(r) = -\frac{k}{r^2}$, this ensures the conservation of the vector

    \begin{equation} \vb{A} = \vb{p}\times\vb{L} -mk\frac{\vb{r}}{r} \end{equation}

    This vector lies in the plane of orbit, and $A^2 = m^2k^2 + 2mE\ell^2$.

    Return to Table of Contents

    Scattering Theory

    This is kinda a mess of angle names and things, but the crucial thing to remember is

    \begin{equation} \dv{\sigma}{\Omega} = \frac{2\pi b\dd{b}}{2\pi\sin\theta\dd{\theta}} = \frac{b}{\sin\theta}\dv{b}{\theta} \end{equation}

    Return to Table of Contents

    Rigid Body Motion

    tree main types of problems:

  • torque free motion: can precess, rotation about intermediate axis unstable
  • nonzero torque, eg weighted top: precession about vertical, axis of rotation can change (nutation)
  • rolling without slipping
    •

    components of rotational motion:

  • rotation about an axis
  • precession of this axis (change in $\phi$)
  • nutation (change in spherical $\theta$), bobbing as it precesses
    •

    basic strategy: decompose motion into CM + movement about the CM

    \begin{gather} \vb{F}_\text{ext} = M \vb{a}_\text{CM} \qquad \vb{r}_\text{CM}\times\vb{F}_\text{ext} = \dv{t}\vb{L}_\text{CM}=\vb{N}_\text{about CM}\\ \vb{L}_\text{tot} = \underbrace{\vb{L}_\text{CM}}_{\mathclap{\vb{r}_\text{CM}\times M\vb{v}_\text{CM}}} + \vb{L}_\text{about CM} \qquad T = T_\text{CM} + T_\text{rot} \end{gather}

    In general, if a vector $\vb{A}$ of fixed length rotates with an angular velocity $\vb*{\omega}$, then

    \begin{equation} \dv{\vb{A}}{t} = \vb*{\omega}\times\vb{A} \end{equation}

    In particular, when measuring $\vb{r}$ from a point $\vb{R}$ on a rotating body,

    \begin{equation} \qty(\dv{\vb{r}}{t})_f = \qty(\dv{\vb{r}}{t})_r + \vb*{\omega}\times\vb{r} \end{equation}

    and we typically take $p$ to be the CM. Taking the second derivative,

    \begin{equation} \vb{a}_f = \vb{a}_r + \underbrace{\ddot{\vb{R}}_f}_{\mathclap{\text{translation}}} + \underbrace{\dot{\omega}\times\vb{r}}_{\mathclap{\text{rotation}}} + \underbrace{2\omega\times\vb{v}_r}_{\mathclap{\text{coriolis}}} +\underbrace{\omega\times(\omega\times \vb{r})}_{\mathclap{\text{centrifugal}}} \end{equation}

    Note that if we have angular momentum in a fixed frame, we can look at the torque in the fixed frame by considering $\dv{L}{t}$ in the body frame; this gives us Euler’s equations

    Inertia tensor

    We define the inertia tensor by

    \begin{equation} I_{jk} = \int \rho(\vb{r})\,(r^2\delta_{ij} - x_ix_j)\dd{V} \end{equation}

    With $\vb{L} = \mathbb{I}\cdot\vb*{\omega}$, kinetic energy is given by

    \begin{equation} T = \frac{1}{2}\vb*{\omega}\cdot\mathbb{I}\cdot\vb*{\omega} = \frac{1}{2}\vb{L}\cdot\mathbb{I}^{-1}\cdot\vb{L} \end{equation}

    We can diagonalize $\mathbb{I}$ using principle axes, which in general are time-dependent.

    Forget about the parallel axis theorem, but it’s good to know the

    perpendicular axis theorem: For a planar object, the sum of the moments of inertia in the plane equals the moment of inertia perpendicular to the plane

    \begin{equation} I_{1,\,\parallel} + I_{2,\,\parallel} = I_\bot \end{equation}

    Table of $I$s for basic shapes (memorize these):

    Shape Rotation Moment of Inertia
    Thin Ring
    (radius $R$)    		 About center $MR^2$
    Ring of Ribbon
    (radius $R$, width $w$)    		 About central diameter $\frac{1}{2}MR^2 + \frac{1}{12}Mw^2$
    Solid cylinder
    (radius $R$)    		 About center $\frac{1}{2}MR^2$
    Hollow cylinder
    (outer radius $R$, inner radius $r$)    		 About center $\frac{1}{2}M(R^2 + r^2)$
    Solid ball
    (radius $R$)    		 About center $\frac{2}{5}MR^2$
    Solid rod
    (length $\ell$)    		 About center $\frac{1}{12}M\ell^2$
    Solid rod
    (length $\ell$)    		 About end $\frac{1}{3}M\ell^2$
    Thin rectangular plate
    (length $\ell$, width $w$)    		 About center $\frac{1}{12}M(\ell^2 + w^2)$

    Rotation under torque: the massive top

    This is fairly straightforward: construct the lagrangian in terms of three angles: a rotation of the top itself, the precession of the top around in a circle, and a nutation angle that allows the top to bob up and down under the influence of gravity. Introduce the conserved quantities (energy and two conjugate momenta to the rotation and precession angles) and make the substitution $u = -\sin\theta$ for $\theta$ the nutation angle, with $\theta=0$ taken to be horizontal. Then $u=1$ means the top is upright, $u=0$ means it’s horizontal, and $u=-1$ means it’s upside down. We get an equation that looks like $\dot{u} = $ a cubic polynomial, so we have 3 cases for where the roots lie: 3 distinct real roots, a double root and a crossing point, or only one real root.

    Rolling without slipping

    The rolling without slipping constraint is that the velocity of the point of contact has to be equal to the velocity of the surface, where the contact point has motion due to the center of mass and its rotation about the CM, and the surface is usually (but not always!) stationary.

    Use Newton’s 2nd and the torque equation $r\times F = N = \dv{L}{t}$. For the torque equation, we can choose torque wrt center of mass, the contact point, or about some arbitrary point (the most useful is usually CM though). You can also solve this with Lagrangians, so long as you can integrate the semiholonomic constraint into one that doesnt depend on the velocities or include it with lagrange multipliers.

    Note that the angular momentum can often be weird! For the right hand side, you don’t necessarily want to rely on the parallel axis theorem to get $I$ so that you can just use $I\omega$; you should directly find \(\vb{L} = \vb{L}_\text{CM} + \vb{L}_\text{rot}\) and take its derivative. (In particular \(\vb{L}_\text{CM}\) is the weird one; the momentum about the CM is typically fine to think of as $I\omega$)

    Return to Table of Contents

    Continuous Media

    Strings and surfaces

    Can solve these using forces: consider a $\dd{\ell}$ or $\dd{A}$ and look at the forces, including tension

    Or can solve them variationally, by minimizing the potential energy with a fixed boundary and length and stuff. To include the length constraint, do it with a multiplier. You can often find ``energy’’ conserved by legendre transforming.

    It’s good to know that waves on strings follow the wave equation with a velocity

    \begin{equation} v = \sqrt{\frac{T\ell}{m}} \end{equation}

    Fluids

    Generally done by looking at forces and energy in a small part of the fluid.

    Pressure $p$ is such that $-\nabla p$ acts as a force per unit volume: it’s the force on a $\dd{V}$ due to the rest of the fluid. Thus there is a potential energy associated to being surrounded by fluid.

    In static fluids, the potential energy is constant, since there is no energy cost to slow mixing. Typically see things like

    \begin{equation} U = p + \psi = \text{const} \end{equation}

    with $\psi$ the potential energy per unit volume due to external forces (e.g. $\psi = \rho g h$ for gravity, but we can also have centrifugal energy $-\frac{1}{2}\rho\Omega^2r^2$ or something)

    In a flowing fluid with no viscosity, we have a similar energy conservation:

    \begin{equation} p + \psi + \frac{1}{2}\rho u^2 = \text{const} \qq{and}\nabla\cdot(\rho \vb{u}) = -\pdv{\rho}{t} \end{equation}

    for an incompressible fluid, the density never changes, so we have that the divergence of the velocity field is zero.

    If the velocity field is rotationless, we can write it as the gradient of some potential $\vb{u} = \nabla\phi$, and incompressibility makes that potential satisfy laplace’s equation $\nabla^2\phi=0$

    Viscosity makes it so that layers of fluid will try to pull on the surrounding layers; say we look at $u_x(z)$ for some fluid flowing in the $x$ direction where lower layers move slower than higher layers. We have

    \begin{equation} F_x/A = \eta\dv{u_x}{z}\qq{for some constant $\eta$} \end{equation}

    Return to Table of Contents

     Prelims  Notes

    Similarly tagged posts