Past exam of the physics B course of the University of Cambridge / 2016 / Paper 1 / Section C

For a central potential, angular momentum is conserved and the motion lies in a plane. With $L=m{r}^2\dot{\theta }$, radial energy conservation becomesAllowed radii have $\mathcal{E}\ge V_{\mathrm{eff}}$; equality gives radial turning points. A well with two finite turning points gives regular radially bound motion. Access to arbitrarily large radii gives an unbound branch, while access to the origin can lead to collision in a singular potential. Stationary points are circular-orbit candidates; a strict local minimum gives radial stability. For a potential tending to zero at infinity, energy relative to that asymptote is useful, but a singular inner potential and centrifugal barrier must still be examined.

Integrating the angular-momentum relation with the origin of angle at $t=0$ givesFor $V=-A/r+B/r^2$ at $L=L_0$, the effective potential is $-A/r+(L_0^2+2m{B}_0)/(2m{r}^2)$. It is exactly the Kepler effective potential with $L_{\mathrm{K}}=\sqrt{L_0^2+2m{B}_0}$, provided that expression is real. Equal radial energy and initial radial phase therefore give identical $r(t)$, although the angular motion differs. Writing $\psi$ for the Kepler angle, $\dot{\theta }/\dot{\psi }=L_0/L_{\mathrm{K}}=k$. With matched periapsis origins,Here $r_0=L_{\mathrm{K}}^2/(mA)$ and the eccentricity is fixed by radial energy. For $B_0>0$, $k<1$ and successive periapses are separated by $2\pi k$, so there is retrograde apsidal motion relative to a closed Kepler ellipse. The radial equivalence is the relevant content of the Newton's theorem of revolving orbits.

For the pure inverse-cube central force set $A=0$. The free-particle reference orbit at $B=0$ is a straight line, $r=d/\cos (\psi -{\psi }_0)$, not a finite-$A$ conic with its old $r_0$ held fixed. To cover attractive as well as repulsive cases completely, putwhere primes denote angle derivatives and $L\ne 0$. The inverse-cube central orbit classification gives the following physically positive-radius branches.

For $\lambda >0$, write $u=C\cos [\sqrt{\lambda }(\theta -{\theta }_0)]$ with $C>0$, so $r=C^{-1}\sec [\sqrt{\lambda }(\theta -{\theta }_0)]$. These are scattering orbits with positive energy and an inner turning radius $1/C$; the angular interval between the two infinite-radius ends is $\pi /\sqrt{\lambda }$. Repulsion ($B>0$) gives an interval smaller than $\pi$; weak attraction ($-L^2/(2m)<B<0$) gives a larger interval and can cause winding. The $B=0$ limit is the straight line.

For $\lambda =0$, $u=C+D(\theta -{\theta }_0)$. If $D\ne 0$, reciprocal-linear spirals connect an infinite-radius end to collision at infinite accumulated angle. If $D=0$ and $C>0$, every radius admits a circular solution at zero energy: attraction exactly balances centrifugal acceleration. Each is a marginal circular orbit under an inverse-cube force, with no radial restoring force.

For $\lambda <0$, put $\beta =\sqrt{-\lambda }$ and $u=C_{+}{e}^{\beta \theta }+C_{-}{e}^{-\beta \theta }$. Up to angle shifts and reversal, the three energy classes arewith positive-radius domains understood. The negative-energy branch is radially bounded but reaches the singular centre on either side of its outer turning point; it is not a regular closed orbit. Collision requires additional physics to continue through $r=0$. Zero angular momentum gives radial trajectories and is excluded from the angular parametrization above.

For a body-fixed mass element at position $\mathbf{r}$ relative to the fixed tip, $\mathbf{v}=\mathbf{\omega }\times \mathbf{r}$. ThereforeThe inertia tensor is real symmetric. Its orthogonal eigenvectors define principal axes and its eigenvalues are principal moments of inertia. angular momentum is parallel to angular velocity when rotation is along a principal axis, but generally not for arbitrary rotation.

Write $\mathbf{n}={\hat{\mathbf{e}}}_3$ and $s=I_3{\omega }_3$. Axial symmetry gives $\mathbf{L}=I_1\mathbf{\omega }+(I_3-I_1){\omega }_3\mathbf{n}$. Since $\dot{\mathbf{n}}=\mathbf{\omega }\times \mathbf{n}$, the vector triple-product identity gives $\mathbf{n}\times \dot{\mathbf{n}}=\mathbf{\omega }-{\omega }_3\mathbf{n}$. HenceDifferentiate in the inertial frame, using constant $s$ and $\dot{\mathbf{n}}\times \dot{\mathbf{n}}=0$. The gravitational torque about the fixed tip is $(\ell \mathbf{n})\times (-mg\hat{\mathbf{z}})$, yielding the symmetry-axis equation of a heavy symmetric top:The tip reaction exerts no moment about that tip.

For steady precession, $\mathbf{n}=(\sin \theta \cos {\Omega }_{P}t,\sin \theta \sin {\Omega }_{P}t,\cos \theta )$. At azimuth zero, put $s_{\theta }=\sin \theta$ and $c_{\theta }=\cos \theta$. Then $\dot{\mathbf{n}}=(0,{\Omega }_P{s}_{\theta },0)$, $\ddot{\mathbf{n}}=(-{\Omega }_P^2{s}_{\theta },0,0)$, and the $y$ component gives $s{\Omega }_P{s}_{\theta }-I_1{\Omega }_P^2{s}_{\theta }{c}_{\theta }=mg\ell s_{\theta }$. For a nonvertical axis this becomesReal steady branches require a nonnegative discriminant. In particular, for a top pointing above the horizontal there is a spin threshold for these fixed-angle solutions. At $\cos \theta =0$ use the quadratic's linear limit ${\Omega }_P=mg\ell /s$, rather than dividing by zero. For a vertical axis, azimuth does not describe a physical conical motion and the division by $\sin \theta$ is inappropriate.

For rapid positive spin $s$ with $s^2\gg 4mg\ell I_1\cos \theta$, the slow and fast steady precession roots areThe approximation assumes a nonsingular fixed inclination. For $g=0$, the quadratic instead gives ${\Omega }_P=0$ or ${\Omega }_P=s/(I_1\cos \theta )$. The zero branch has a stationary symmetry axis and pure axial spin. On the other branch angular momentum is constant and parallel or antiparallel to the chosen $z$ axis; $L_z=I_1{\Omega }_P$ and $s=L_z\cos \theta$. This is torque-free precession in space, not the different rate at which transverse components rotate in body coordinates.

Take $\zeta$ as signed distance from the neutral axis within the beam cross-section. Bending into a large radius $R$ changes the length of a fibre proportionally to $\zeta /R$, so linear elasticity gives $\sigma =E\zeta /R$, up to the curvature sign convention. The moment magnitude isHere $I$ is the second moment of area about the neutral axis, with units of length to the fourth power, not the mass moment of inertia. A homogeneous cross-section has its neutral axis through the centroid under pure bending.

For the straight pin-ended beam, the compressive dead load produces bending moment $-Fy$ at transverse offset $y$. In the small-slope approximation, curvature is $y^{\prime \prime }$ with the corresponding sign convention, givingThe stated endpoints fix positions rather than slopes: this interpretation is consistent with the illustrated end forces and the later sinusoidal unstressed shape. Truly clamped endpoints would have a different threshold. A nonzero solution $y=A_n\sin (n\pi x/L)$ requires $F_n=n^2{\pi }^{2}EI/L^2$. Thus the first Euler load is $F_c={\pi }^{2}EI/L^2$. Below it the perfect straight state is stable, at it the first mode is marginal in linear theory, and above it the straight state is unstable. The linear equation selects thresholds but cannot set the amplitude or describe a finite postbuckled equilibrium.

For the initially curved compressed beam, elastic stress responds to the change from its stress-free curvature. With $y=y_0+\eta$, the moment is $EI(y^{\prime \prime }-y_0^{\prime \prime })=EI{\eta }^{\prime \prime }$, so force balance givesPut $q=\pi /L$ and seek the branch continuously connected to the unloaded shape. Since $y_0=Y\sin qx$, substitute $\eta =C\sin qx$ to obtain $-EI{q}^{2}C=-F(Y+C)$, orAway from eigenvalues, the homogeneous solution $A\sin (\sqrt{F/(EI)}x)+B\cos (\sqrt{F/(EI)}x)$ vanishes under the two end conditions, so this particular solution is the unique linear equilibrium. For $0\le F<F_c$ and $Y>0$, the maximum additional displacement is ${\eta }_{\max }=YF/(F_c-F)$ and maximum total displacement is $y_{\max }=Y/(1-F/F_c)$, both at midspan.

This sinusoidal imperfection amplification of a column produces deflection at every positive load and strong amplification near $F_c$. The apparent divergence signals loss of the small-slope model before infinite displacement; nonlinear geometry is needed there. The formal continuation above the first threshold is not a stable physical loading branch. Both total and additional displacement are distinguished in the graph.

irrotational flow has $\nabla \times \mathbf{v}=0$, so locally $\mathbf{v}=\nabla \Phi$; the simply connected fluid domain here permits a globally single-valued velocity potential. Incompressibility gives $\nabla \cdot \mathbf{v}=0$, hence ${\nabla }^2\Phi =0$. In a multiply connected domain, zero local vorticity alone would not rule out circulation and a multivalued potential.

Take $k>0$ so the disturbance decays at depth. Differentiating the given candidate potential yieldsThe $x$ and $z$ second derivatives of $\Phi$ cancel in the Laplacian. At $z=0$, $v_z=h_0\omega \cos (\omega t)\sin (kx)={\partial }_{t}Z$, so the linearized kinematic free-surface condition holds. This uses $k{h}_0\ll 1$; advection of the displaced surface is higher order.

Since ${\cos }^2(kx)+{\sin }^2(kx)=1$, the local squared speed is $h_0^2{\omega }^2{\cos }^2(\omega t)e^{2kz}$. The depth integral is thereforeIt is independent of horizontal position; no further spatial average is needed for this kinetic energy. Replacing the upper boundary $Z$ by zero changes it only beyond quadratic order in amplitude.

Relative to the flat surface, the gravitational energy at each horizontal position changes by $\rho g{\int }_0^{Z}zdz=\rho g{Z}^2/2$. Averaging over one wavelength givesThe energy of a deep-water standing gravity wave, averaged over wavelength, is $T+⟨U{⟩}_x$ and must be constant for a free inviscid standing wave, which requires equal coefficients of the complementary squared sine and cosine:As an independent check, the linearized dynamic boundary condition ${\partial }_t\Phi +gZ=0$ at $z=0$ gives exactly $g-{\omega }^2/k=0$. Infinite depth, negligible surface tension and no dissipative forcing are part of this gravity-wave model.