Past exam of the mathematics course of the year 2 of the natural sciences course of the University of Cambridge / 2011 / Paper 2 / 7

The Lagrangian is $\mathcal{L}=T-\mathcal{V}$, and the Euler-Lagrange equations areAt a rest equilibrium $q=0$, ${\partial }_i\mathcal{V}(0)=0$ for all three coordinates. The Taylor expansion therefore starts asThe symmetric Hessian matrix is the matrix $\mathsf{V}$ of the quadratic potential energy. If a quadratic coefficient vanishes along a direction, higher-order terms can be the first nonconstant terms there; equilibrium by itself does not imply strict stability.

Write ${\mathsf{T}}_0=\mathsf{T}(0)$. Terms involving the variation of $\mathsf{T}(q)$ are cubic in small displacements and velocities, so the quadratic small oscillation Lagrangian isA normal mode has fixed displacement ratios, $q(t)=v\cos (\omega t+\varphi )$. Its angular frequency $\omega$ is an eigenfrequency, or normal frequency, and satisfies the generalized eigenvalue problem for small oscillationsHere ${\mathsf{T}}_0$ is positive definite for independent physical coordinates and positive masses. A stable quadratic minimum makes ${\omega }^2\ge 0$; a negative value would describe exponential instability rather than oscillation.

For the suspended bar, let ${\kappa }_0$ and ${\kappa }_1$ denote the actual incremental spring constants of the upper and each lower spring. The physical displacement of the bar is $q_1$, while the lower spring extensions change by $q_2-q_1$ and $q_3-q_1$. ConsequentlyThe linear gravitational and elastic terms cancel at the loaded equilibrium. Expanding the remaining squares gives the small-oscillation mass and stiffness matricesThusSource convention: the printed matrix is obtained by setting ${\kappa }_0=\lambda /L$ and ${\kappa }_1=\mu /l$. Although the PDF prose calls $\lambda ,\mu$ spring constants, that matrix treats them as elastic parameters whose ratios to the stated lengths are the incremental stiffnesses. If they are read literally as Hooke-law spring constants, the divisions by $L,l$ should instead be omitted. The numerical case has $L=l=1$, so this discrepancy does not affect its results.

With the numerical parameters, ${\mathsf{T}}_0=\mathrm{diag}(2,1,1)$ andUse the exchange symmetry of the lower masses to find the vertical modes of a horizontally constrained suspended bar. The antisymmetric vector $v_a=(0,1,-1{)}^T$ obeys $\mathsf{V}{v}_a={\mathsf{T}}_0{v}_a$, giving ${\omega }_a=1$. In the symmetric sector put $v=(A,B,B{)}^T$ and $s={\omega }^2$. The two independent equations areTheir determinant gives $s^2-4s+2=0$, so $s=2\pm \sqrt{2}$. Taking $A=1$, the full set isIn the lowest normal mode, all three masses move in phase, with the lower masses moving about $2.414$ times as far as the bar. In the antisymmetric normal mode, the bar stays still and the lower masses move equally in opposite directions; their changing forces on the bar cancel. In the highest normal mode, the lower masses move together against the bar, with magnitude about $0.414$ of its displacement. Any small motion is a superposition of these three modes.