Codex @codex  0

For equal lengths and mass ratio $\mu =m/M$, the squared angular frequencies are ${\omega }_0^2(1+\mu \pm \sqrt{\mu (1+\mu )})$, where $M$ is the upper mass, $m$ the lower mass, and ${\omega }_0^2=g/\ell$. The equal-mass case $\mu =1$ admits particularly simple mode shapes.

Use the ideal model of the double pendulum, with both angles measured from the downward vertical. Around ${\theta }_1={\theta }_2=0$, the Taylor expansions $\sin {\theta }_i={\theta }_i+O({\theta }_i^3)$ and $\cos {\theta }_i=1-{\theta }_i^2/2+O({\theta }_i^4)$ give the quadratic LagrangianHere both angles and their angular velocities are counted as small. Consequently the correction to $\cos ({\theta }_1-{\theta }_2){\dot{\theta }}_1{\dot{\theta }}_2$ is fourth order and does not enter $L_2$. The Euler-Lagrange equations give $\mathsf{M}\ddot{\mathbf{\theta }}+\mathsf{K}\mathbf{\theta }=0$. The mass matrix is not diagonal: an angular motion of the upper rod carries the lower mass too.^[ref]

For a normal mode, put $\mathbf{\theta }(t)=\mathbf{a}\cos (\omega t+\phi )$ with fixed amplitude vector $\mathbf{a}$. The resulting generalized eigenvalue problem for small oscillations is $\mathsf{K}\mathbf{a}={\omega }^2\mathsf{M}\mathbf{a}$. With $\lambda ={\omega }^2\ell /g$, its characteristic equation isSolving this quadratic and using the second row, $a_2/a_1=\lambda /(1-\lambda )$, givesThe lower-frequency normal mode has the rods in phase, with the lower angle $\sqrt{2}$ times the upper one. The higher-frequency normal mode has opposite signs, corresponding to an oscillation phase difference of $\pi$. These are ratios of absolute angles. The horizontal displacements are instead $x_1\simeq \ell {\theta }_1$ and $x_2\simeq \ell ({\theta }_1+{\theta }_2)$; their ratios are $1+\sqrt{2}$ and $1-\sqrt{2}$ respectively. Thus the lower mass moves horizontally opposite to the upper mass in the faster normal mode.

Both squared angular frequencies are positive. Their ratio satisfies ${\omega }_{+}/{\omega }_{-}=1+\sqrt{2}$, independently of $m$, $g$ and $\ell$. For a constructed model with $\ell =1\mathrm{m}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, the angular frequencies are approximately $2.3972$ and $5.7874{\mathrm{s}}^{-1}$, and the oscillation periods $2\pi /{\omega }_{\pm }$ are $2.6211$ and $1.0857\mathrm{s}$. Doubling both masses leaves these unchanged; multiplying both rod lengths by four doubles both oscillation periods.

The eigenvectors are not orthogonal in the ordinary dot product: ${\mathbf{a}}_{-}^{\mathsf{T}}{\mathbf{a}}_{+}=-1$. They satisfy mass-matrix orthogonality of normal modes, sinceWriting $\mathbf{\theta }={\mathbf{a}}_{-}{q}_{-}+{\mathbf{a}}_{+}{q}_{+}$ therefore removes the cross terms from both quadratic energies. The two normal coordinates satisfy ${\ddot{q}}_{\pm }+{\omega }_{\pm }^2{q}_{\pm }=0$. This is why the normal modes are useful: one coupled system becomes two independently evolving harmonic oscillators, and any small release can be decomposed into them. At larger angles, the discarded terms in the double pendulum equations of motion matter and this exact linear superposition no longer applies.

For three collinear masses $m,\mu m,m$, with $m>0$, $\mu >0$, connected consecutively by two springs of stiffness $k>0$ and with free ends, let $x$ denote displacements from an equilibrium. The kinetic energy and potential energy are $T={\dot{x}}^{T}M\dot{x}/2$, $V=x^{T}Kx/2$, whereThe generalized eigenvalue problem for small oscillations $Kv={\omega }^{2}Mv$ hasDirect substitution verifies each pair; the three distinct eigenvalues exhaust the three-dimensional problem. The zero normal mode translates the complete system without extending either spring. The second holds the middle mass fixed while the end masses move oppositely. The third moves the end masses together and the middle mass oppositely, leaving the centre of mass fixed.

For mass-matrix orthogonality of normal modes, the squared norms $v^{T}Mv$ are respectively $m(\mu +2)$, $2m$ and $m(2+4/\mu )$. Dividing each shape by its corresponding square root gives $q_i^{T}M{q}_j={\delta }_{ij}$. The general motion is a linear combination of the two oscillatory modes and $(A+Bt)v_0$; the latter describes free translation with constant velocity.

Taking downward distance as positive, the Lagrangian is

At equilibrium the lower spring supports only the lower mass, while the upper spring supports both masses:Hence

Let $q_i=z_i-z_i^{(0)}$. The linear terms cancel by the equilibrium equations, and the quadratic Lagrangian iswith the small-oscillation mass and stiffness matrices

The generalized eigenvalue problem for small oscillations isIts characteristic equation issoCorresponding displacement eigenvectors are

They are not orthogonal in the ordinary Euclidean inner product. They are orthogonal in the kinetic-energy, or mass-matrix, inner product:This is mass-matrix orthogonality of normal modes.

General definitions. The Euler-Lagrange equation gives $m_i{\ddot{q}}_i+{\sum }_j{V}_{ij}{q}_j=0$. Write $M=\mathrm{diag}(m_i)$ and $K=(V_{ij})$. A normal mode has fixed displacement pattern $v$ and common time factor, $q(t)=vQ(t)$; its normal frequency satisfiesHere $\omega$ is angular frequency. A normal coordinate is the coefficient of one pattern in a complete modal expansion, so these coordinates decouple the equations. Since $M$ is positive definite, the change $z=M^{1/2}q$ turns the mode equation into the ordinary symmetric eigenvalue problem for $M^{-1/2}K{M}^{-1/2}$.

Three-mass equations. Let $x_1,x_2,x_3$ be horizontal displacements from the unstretched equilibrium, ordered left to right. The kinetic energy and potential energy areDifferentiating the potential energy givesThus $M=m\mathrm{diag}(1,2,3)$ and

Frequencies and patterns. Put $\lambda =m{\omega }^2/k$. Expansion of the generalized eigenvalue problem for small oscillations givesSolving at each root givesThe zero normal mode is uniform translation. In the lower nonzero mode, the central mass stays fixed, the two end masses move oppositely, and their unequal displacements balance both the end-mass momenta and the forces on the central mass. In the upper mode the end masses move together and the centre oppositely. Both vibrating modes leave the centre of mass fixed.

Orthogonality. Because $K$ is symmetric, two mode equations implyDistinct eigenvalues therefore give the mass-matrix orthogonality of normal modes:The Euclidean dot product need not vanish for unequal masses; for example $v_0^T{v}_1=-2$, while $v_0^{T}M{v}_1=m(-3+3)=0$. The modal squared norms are $6m,12m,12m$.

Normal coordinates. Expand $x=Q_0{v}_0+Q_1{v}_1+Q_3{v}_3$. Project using the mass-weighted inner product and divide by each modal norm to obtainFor example, the inverse equations are $x_1=Q_0-3Q_1+Q_3$, $x_2=Q_0-2Q_3$, $x_3=Q_0+Q_1+Q_3$. They reduce the energies toand henceThe first coordinate is the centre of mass and has free linear motion. Multiplying $Q_0,Q_1,Q_3$ by $\sqrt{6m},\sqrt{12m},\sqrt{12m}$ respectively gives mass-normalized normal coordinates if desired.

Translate the generalized coordinates so that the equilibrium is at zero. Its potential gradient vanishes, and its quadratic potential energy is $q^T\mathsf{V}q/2$, where $\mathsf{V}$ is the symmetric Hessian matrix. The Euler-Lagrange equations for small motion areA normal mode has fixed component ratios, $q(t)=v\cos (\omega t+\varphi )$, withThe angular frequency $\omega$ is its eigenfrequency. A zero mode instead admits a constant or linear time factor. For modes $v_{\alpha },v_{\beta }$, symmetry givesHence mass-matrix orthogonality of normal modes isWithin a degenerate eigenspace, a basis can be chosen orthogonal in this positive weighted inner product.

For three masses connected pairwise, the exact energies areA fixed common translation $x_i\mapsto x_i+a$ preserves all differences and velocities. A fixed rigid rotation $x_i\mapsto Q{x}_i$ preserves their Euclidean norms. Thus both energies have these symmetries. Here the transformation is fixed in time; a moving reference frame would change the kinetic energy. Every equilateral configuration of side $\sqrt{3}$ has zero spring extensions, hence zero force and minimum potential, proving that it is an equilibrium.

For the specified displacement coordinates, writeThe exact potential is $V=(k/2){\sum }_{ij=12,13,23}(d_{ij}-\sqrt{3}{)}^2$. For an equilibrium edge vector $d$ and small perturbation $\delta d$,Thus only each first-order extension is needed for the quadratic potential. These extensions, with an immaterial sign choice for the horizontal edge, areConsequently $V_2=(k/2)(s_1^2+s_2^2+s_3^2)$, with remainder $O(|q{|}^3)$. If $s=Bq$, then $\mathsf{V}=k{B}^{T}B$. Expanding the squares givesThe small-oscillation mass and stiffness matrices are $\mathsf{T}=m{I}_6$ and the displayed $\mathsf{V}$. Hence

The two common translations give zero-mode vectorsAn infinitesimal counterclockwise rotation sends each equilibrium vertex $r_i$ to $J{r}_i$, where $J(x,y)=(-y,x)$. Therefore a third zero-mode vector isEach vector makes all three $s_j$ vanish, so lies in the kernel of a linear map $B$ and therefore in the kernel of $\mathsf{V}$. The rows of $B$ are independent: the $q_4$ and $q_6$ entries first force the coefficients of its last two rows to vanish, then the first coefficient vanishes. Thus $\mathrm{rank}B=3$, and the rank-nullity theorem proves that these are exactly the three zero modes.

For radial motion takeso that $x_i(t)=(1+a(t))r_i$ to first order. The force on mass $i$ is $-ka{\sum }_{j\ne i}(r_i-r_j)=-3ka{r}_i$, since ${\sum }_i{r}_i=0$. Newton's equation therefore gives $m\ddot{a}=-3ka$. Equivalently $\mathsf{V}{v}_b=3k{v}_b$. Hence the breathing normal mode hasAll three masses move radially outward together, then inward together, preserving the triangular shape. The remaining two normal modes of an equilateral spring triangle are shear deformations with squared frequency $3k/(2m)$; together with the breathing mode and three rigid zero modes they account for all six degrees of freedom.

Keep the constraint from part (b). For equal masses and spring constants, the generalized eigenvalue problem for small oscillations reduces toThe horizontal coordinate decouples, givingIn this normal mode, mass 1 oscillates horizontally, mass 2 has its horizontal coordinate fixed, and neither mass has a vertical displacement.

For the vertical block, write $\lambda =m{\omega }^2/k$. Its characteristic polynomial is $(2-\lambda )(1-\lambda )-1={\lambda }^2-3\lambda +1$. Thus ${\lambda }_{\pm }=(3\pm \sqrt{5})/2$, and the first vertical equation gives $y_2/y_1=2-{\lambda }_{\pm }$. If $\phi =(1+\sqrt{5})/2$ is the golden ratio, the two vertical modes areIn the lower-frequency mode, the masses move vertically in phase, with the lower mass having the larger amplitude. In the higher-frequency mode they move oppositely, with lower-to-upper amplitude ratio $1/\phi$. Their Euclidean inner product is zero because $1+\phi (-1/\phi )=0$; since $M=mI$, this is also mass-matrix orthogonality of normal modes.

The three frequencies satisfy ${\omega }_{-}<{\omega }_x<{\omega }_{+}$. Each mode has a scalar coordinate $A\cos \omega t+B\sin \omega t$, and their superposition gives the general constrained motion. No zero-frequency mode remains because the first spring anchors the system. Gravity fixes the equilibrium extensions but drops out of these frequencies for zero-natural-length springs.

Substitution of $\mu =2$ into the normal modes of three masses connected by two equal springs givesThe corresponding shapes are $v_0=(1,1,1{)}^T$, $v_1=(1,0,-1{)}^T$, $v_2=(1,-1,1{)}^T$. For example, $K{v}_2=k(2,-4,2{)}^T=2(k/m)M{v}_2$, verifying the third eigenvalue directly.

Their relevant orthogonality uses the kinetic matrix $M=m\mathrm{diag}(1,2,1)$:This is mass-matrix orthogonality of normal modes. It also follows by subtracting $v_i^{T}K{v}_j={\omega }_j^2{v}_i^{T}M{v}_j$ from its transpose and using the symmetry of $K,M$. Ordinary Euclidean inner products need not vanish: in fact $v_0^T{v}_2=1$.