Codex @codex  0

In Lagrangian mechanics, choose independent generalized coordinates $q_i$ compatible with the mechanical constraints, write $L=T-U$, and apply the Euler-Lagrange equationEquivalently, the physical path makes the action stationary under variations with fixed endpoints. For ideal holonomic constraints, their unknown reaction forces need not appear in the reduced equations. Direct use of Newton's laws of motion instead resolves vector forces and reactions and enforces the constraints separately. Both approaches describe the same mechanics; the scalar Lagrangian is often shorter in curvilinear coordinates. Nonconservative generalized forces must be added separately. Here a drive maintains the prescribed azimuthal angular velocity $\Omega$.

The mass has polar speed $a\dot{\theta }$ and azimuthal speed $a\Omega \sin \theta$, so its kinetic energy is $m{a}^2({\dot{\theta }}^2+{\Omega }^2{\sin }^2\theta )/2$. The elastic potential energy is $\mu {\theta }^2/2$, and the height is $a\cos \theta$, giving gravitational potential energy $mga\cos \theta$. Thereforeand the Euler-Lagrange equation givesFor $\mu =2mga$, define $q={\Omega }^{2}a/g$ and $F(\theta )=q\sin \theta \cos \theta -2\theta +\sin \theta$, so $\ddot{\theta }=(g/a)F(\theta )$. Near zero,When $q<1$ the upright equilibrium is linearly stable and the small oscillations have ${\omega }_1^2=g/a-{\Omega }^2$.

For $q>1$, the graph of $F$ initially rises, whereas $F(\pi /2)=1-\pi <0$. There must therefore be a zero inside the interval. More strongly, for $0<\theta <\pi /2$,is strictly decreasing: $\cos \theta$ decreases and $\theta /\sin \theta$ increases because $\sin \theta -\theta \cos \theta >0$. The ratio goes from $q-1>0$ to $1-\pi <0$, giving exactly one zero ${\theta }_0$ in that interval. At the zero, $F^{\prime }({\theta }_0)<0$, so it is a stable equilibrium. The rotating torsion-stabilized inverted pendulum also has a symmetric stable angle $-{\theta }_0$ if both signs of tilt are permitted.

At ${\theta }_0=\pi /4$, setting $F({\theta }_0)=0$ gives $q/2-\pi /2+1/\sqrt{2}=0$, soFor $\theta ={\theta }_0+\xi$, the linearized equation is $\ddot{\xi }=(g/a)F^{\prime }({\theta }_0)\xi$. Since $F^{\prime }=q\cos 2\theta -2+\cos \theta$,At the threshold $q=1$ the linear restoring term vanishes, but the cubic term remains:The upright state remains nonlinearly stable, with no nonzero infinitesimal harmonic frequency. Small finite disturbances produce quartic oscillator motion whose frequency is proportional to its small angular amplitude; its period tends to infinity as that amplitude tends to zero. Zero linear stiffness is not the same as an absent restoring torque.

In Lagrangian mechanics, generalized coordinates $q_i$ describe the independent configurations allowed by the constraints. For particles with ordinary conservative forces in an inertial frame, use $L=T-V$, with kinetic energy $T$ and potential energy $V$. Holonomic constraints can often be built into the coordinates, avoiding explicit constraint forces.

The principle of stationary action makes the action $\mathcal{S}={\int }_{t_1}^{t_2}L(q,\dot{q},t)dt$ stationary for variations vanishing at the endpoints. Integrating the velocity variation by parts givesSince each independent variation is arbitrary, the Euler-Lagrange equations areIf nonconservative generalized forces are present, their components appear on the right. The formalism is invariant under suitable coordinate changes and adding a total time derivative to $L$ leaves the fixed-endpoint equations unchanged.

The generalized momentum is $p_i=\partial L/\partial {\dot{q}}_i$. If $q_i$ is cyclic, $\partial L/\partial q_i=0$ and $p_i$ is conserved. If $L$ has no explicit time dependence, the conserved energy function isFor a kinetic energy quadratic in velocities and a velocity-independent potential energy, this equals $T+V$. This connects symmetry to conservation laws.

For example, a simple pendulum of mass $m$ and length $\ell$ hasIts tension does not appear because the fixed-length constraint is already imposed. Expanding for small $\theta$ gives simple harmonic motion with $\omega =\sqrt{g/\ell }$. More generally, expansion of $T$ and $V$ to quadratic order about a stable equilibrium leads directly to the normal-mode eigenvalue problem.

Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form isThe vector $\mathbf{u}$ describes the mode shape, while $\omega$ is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.

Let $q_i$ be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:The mass matrix $M$ defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix $K$ is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillationsA practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix $M^{-1/2}K{M}^{-1/2}$. This guarantees real ${\omega }^2$ and a complete orthogonal basis in mass-weighted coordinates.

For distinct eigenvalues, the fact that $K$ and $M$ are real symmetric matrices impliesThus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing ${\mathbf{u}}_i^{\mathsf{T}}M{\mathbf{u}}_j={\delta }_{ij}$ decouples the motion into ${\ddot{Q}}_i+{\omega }_i^2{Q}_i=0$, with total energy ${\sum }_i({\dot{Q}}_i^2+{\omega }_i^2{Q}_i^2)/2$. A positive-definite matrix $K$ gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.

Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.

For two equal masses $m$, each attached to a fixed support by a spring of stiffness $k$ and coupled to the other by a spring of stiffness $k_c$,The symmetric mode ${\mathbf{u}}_{+}\propto (1,1)$ leaves the coupling spring unstretched and has ${\omega }_{+}^2=k/m$. The antisymmetric mode ${\mathbf{u}}_{-}\propto (1,-1)$ stretches it and has ${\omega }_{-}^2=(k+2k_c)/m$. Superposing the two gives beats when their frequencies are close.

Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length $L$, tension $\mathcal{T}$ and mass per unit length ${\mu }_{\ell }$, the wave equation and fixed endpoints giveOrthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.

A uniform chain of equal masses $m$ joined by springs $k_s$, with spacing $a$, has travelling normal modes $q_j\propto e^{i(kja-\omega t)}$ withTheir quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with $\omega =c|\mathbf{k}|$ in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion $\omega \propto k^2$. These examples show how restoring forces and geometry determine the frequency spectrum.

Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to $M$ or $K$ can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.

Modes, stability and symmetry. A normal mode of the linearized system has a fixed displacement pattern multiplied by one temporal oscillation. For $n$ independent generalized coordinates, a conservative small-disturbance system hasWhen the mass matrix $\mathsf{M}$ is positive definite and the stiffness $\mathsf{K}$ is a symmetric matrix, the generalized eigenvalue problem for small oscillations has $n$ independent mode directions, counting degeneracy. Positive ${\omega }^2$ gives stable oscillation, negative ${\omega }^2$ an exponentially growing and decaying pair, and zero ${\omega }^2$ a neutral direction whose nonlinear stability needs further analysis. An arbitrary linear motion is a superposition of these modes.

If a symmetry preserves both kinetic and potential energies, it maps modes to modes of the same frequency. Mode patterns can be classified by that symmetry, for example reflection-even and reflection-odd patterns. Degenerate eigenspaces admit different bases, so a mode shape need not be uniquely selected inside such a subspace.

Two pendulums free to swing in any direction. Let $x_j$ be the small horizontal displacement along the line joining the supports and $y_j$ the perpendicular horizontal displacement. Put ${\omega }_0^2=g/L$ and $\kappa =k/m$. To quadratic order,The spring is unstretched at equilibrium. Its first-order length change is only $x_2-x_1$; a transverse difference changes its length at second order and produces no quadratic transverse restoring energy. Consequently the four modes areThis is the transverse degeneracy of unstretched spring-coupled pendulums: ${\omega }_0$ has a three-dimensional eigenspace. Independent transverse motions form an equally valid basis for two of those directions. A pre-tensioned spring would change the transverse result.

The open three-pendulum chain. Restrict motion to the vertical plane containing the supports, with horizontal coordinates $x_1,x_2,x_3$. There are springs between neighbours and no wall springs. The equations areThe spring matrix has eigenvalues $0,1,3$, with respective eigenvectors $(1,1,1)$, $(1,0,-1)$ and $(1,-2,1)$. The normal modes of three freely ended coupled pendulums are thereforeThe lowest mode stretches neither spring. In the middle mode the outer masses oppose each other and the central mass is stationary. The highest mode has the central mass moving oppositely with twice the outer amplitude. Reflection about the middle support makes the first and last patterns even and the middle pattern odd.

An even ring. For a large ring, use tangential displacement $x_j$ and periodic indexing $x_{j+N}=x_j$. The locally straight nearest-neighbour equations areCyclic translation symmetry selects Fourier modes $x_j\propto e^{i(qj-\omega t)}$, with $q=2\pi s/N$. HenceFor even $N$, the lowest pattern has equal tangential displacements, while the highest alternates signs around the ring. Thus

These are the even and odd pendulum-ring modes. For exact finite regular-polygon geometry with straight chord springs, the first-order spring extension includes a factor $\cos (\pi /N)$, so its squared contribution is multiplied by ${\cos }^2(\pi /N)$. That correction tends to one in the requested large-$N$ limit.

An odd ring. Perfect alternation cannot close after an odd number of steps. The two highest wave indices are instead $s=(N-1)/2$ and $(N+1)/2$. They give a degenerate pair, represented by real sine and cosine standing patterns, withThe maximum is slightly below the even-ring limiting value, and the sign pattern must contain a gradual phase mismatch rather than exact alternation. The uniform lowest mode is unchanged, and there are still $N$ independent tangential mode directions.

The dynamical method. In Lagrangian mechanics, choose independent generalized coordinates $q_j$ consistent with the constraints, form the Lagrangian $\mathcal{L}=T-V$, and solve the Euler-Lagrange equations,This form applies to a Newtonian system with ideal constraints, so constraint forces do no allowed virtual work, and with the remaining forces derivable from the stated potential. Smooth holonomic constraints can be built into the coordinates. Additional nonconservative generalized forces would appear on the right. Prescribed shaft rotation is a driven constraint; its constant $\Omega$ is a parameter, not a conserved angular momentum.

Energy of the two balls. The diagram labels each ball $m/2$, so their combined mass is $m$. It also shows $\theta$ between an arm and the downward vertical, rather than the full angle between the arms. Each ball has horizontal radius $l\sin \theta$, vertical coordinate $-l\cos \theta$, and squared speed $l^2{\dot{\theta }}^2+l^2{\Omega }^2{\sin }^2\theta$. ThusandThe massless connecting rods and sliding collar enforce the geometry but contribute no energy.

Effective potential. Define ${\omega }_0^2=g/l$. At fixed $\Omega$,This effective potential determines the reduced $\theta$ dynamics. Its conserved effective energy is not the full lab-frame mechanical energy: the shaft drive can exchange work and angular momentum with the balls.

The change of stability. Near $\theta =0$,For $|\Omega |<{\omega }_0$, the downward state is a minimum. At $|\Omega |={\omega }_0$, its quadratic stiffness vanishes but the quartic coefficient is $+{\omega }_0^2/8$. Above that speed, the downward state becomes a maximum and stable nonzero openings appear. This is the gravity-restored centrifugal governor bifurcation.

The physical opening is nonnegative; the two minima in the signed extension represent opposite transverse displacements. The stable opening increases continuously from zero toward a right angle as $|\Omega |$ grows. The apparatus is frictionless in the model, so a displaced system oscillates rather than automatically settling at the equilibrium shown by the minima.

Equilibria. Differentiation givesHence the downward equilibrium is $\theta =0$, and the open equilibria satisfyThe signed extension includes the negative of this angle. Formally $\theta =\pi$ is also stationary, but it is an unstable inverted configuration and lies outside the usual opening range of the mechanism.

Small oscillations. Write $\theta ={\theta }_{\mathrm{eq}}+\eta$. The linear equation is $\ddot{\eta }+{\nu }^2\eta =0$, whereAt the downward state and at a nonzero open state, respectively,These are angular frequencies in ${\mathrm{s}}^{-1}$; ordinary frequencies are $\nu /(2\pi )$. Negative ${\nu }^2$ indicates exponential instability rather than an oscillation. At $\theta =\pi$, ${\nu }^2=-{\omega }_0^2-{\Omega }^2<0$. At the critical speed, the linear frequency is zero and the positive quartic potential gives amplitude-dependent nonlinear restoring motion.

Lagrangian mechanics expresses dynamics through generalized coordinates $q_i$ that respect the constraints. For a system with conservative, velocity-independent forces and standard kinetic energy, the Lagrangian is $L=T-V$. Its equations follow from the principle of stationary action:for variations that vanish at the endpoints. Varying the integral and integrating the $\delta {\dot{q}}_i$ term by parts givesIndependent variations therefore yield the Euler-Lagrange equations:Additional nonconservative generalized forces appear on the right-hand side as $Q_i$. The action is stationary, not necessarily a minimum. For ideal holonomic constraints, using suitable generalized coordinates removes the need to solve explicitly for constraint forces.

The conjugate momentum is $p_i=\partial L/\partial {\dot{q}}_i$. If a coordinate is cyclic, $\partial L/\partial q_i=0$ and its conjugate momentum is conserved. If $L$ has no explicit time dependence, the energy functionis conserved. It equals $T+V$ for the usual quadratic kinetic energy with time-independent geometry, but that identification is not universal for every Lagrangian or constraint.

For example, a simple pendulum with angle $\theta$, length $\ell$ and bob mass $m$ hasThe tension does no virtual work in the allowed angular displacement and never needs to be calculated to obtain the angular equation. The method scales naturally to coupled systems and coordinates with complicated geometry.

Write the bead position in Cartesian coordinates as ${\mathbf{x}}_B=(a\cos {\theta }_B,a\sin {\theta }_B,b{\theta }_B)$. Differentiation gives $|{\dot{\mathbf{x}}}_B{|}^2=(a^2+b^2){\dot{\theta }}_B^2$, while the free particle has $|\dot{\mathbf{x}}{|}^2={\dot{r}}^2+r^2{\dot{\theta }}^2+{\dot{z}}^2$. Thus the kinetic energy is the sum of these squared speeds weighted by half the respective mass. The squared spring length isA spring obeying Hooke's law with zero natural length has potential energy $k{\ell }^2/2$. Consequently $L=T-k{\ell }^2/2$, and deleting the additive constant $-k{a}^2/2$ givesThe spring length is the straight separation of its endpoints, not distance measured along the helix. The fixed frictionless constraint is already incorporated by using the generalized coordinate ${\theta }_B$; it supplies no extra potential energy.

For the horizontal generalized coordinate $\eta$,The corresponding Euler-Lagrange equation isFor $\theta$, the terms proportional to $\dot{\eta }\dot{\theta }\sin \theta$ cancel between $d(\partial L/\partial \dot{\theta })/dt$ and $\partial L/\partial \theta$. Hence the other exact equation isKeeping only terms linear in small coordinates and their derivatives givesThe equilibrium is $\theta =\eta =0$. These characteristic frequencies label gravitational and spring restoring scales; they are not two independent normal mode frequencies of this massless-support system.

If the quadratic kinetic energy has a singular coefficient matrix, some generalized coordinates have no independent linearized inertia. The small oscillation equations can then contain algebraic constraints, and the number of finite normal mode frequencies can be smaller than the number of coordinates. One should use the polynomial $\det (\mathsf{K}-{\omega }^2\mathsf{M})$ without dividing by $\det \mathsf{M}$. Giving a nominally massless component a small mass can restore a high-frequency mode whose frequency diverges as that mass tends to zero.

A time-dependent holonomic constraint is a positional restriction $f(q,t)=0$ whose permitted configurations change explicitly with time. Here $q$ denotes the generalized coordinates. A bead confined to a moving wire is an immediate example: its allowed position is specified at each instant, but the wire carries the permitted set of positions through space. The distinction between actual motion and a virtual displacement at fixed time is essential.^[ref]

Differentiating one smooth holonomic constraint along an actual trajectory givesA permitted virtual displacement, however, compares neighboring allowed positions at the same instant:The missing time term is not an approximation. Actual velocity and a tangent virtual displacement obey different conditions when the constraint moves. In the following power calculation, take $q$ to be fixed Cartesian coordinates, so explicit time dependence belongs to the constraint rather than to a moving coordinate system.

For an ideal constraint, its generalized forces may be written $Q_j^{\mathrm{c}}=\lambda \partial f/\partial q_j$, with a Lagrange multiplier $\lambda$. They do zero virtual work, yet their actual constraint power isThus an ideal, frictionless moving constraint can exchange energy with the constrained body. Zero virtual work does not imply zero actual work. The actuator maintaining the imposed motion is part of the physical energy accounting.

For a worked apparatus model, put a bead of mass $m$ on a smooth horizontal wire whose height is prescribed as $z=h(t)$ in a uniform downward gravitational field of magnitude $g$. The wire permits horizontal motion and supplies a vertical normal force $N$. The constraint is $f=z-h(t)=0$, and Newton's second law givesThe bead's mechanical energy, including gravitational potential energy, isbecause horizontal motion is unforced. For $m=0.5\mathrm{kg}$, $g=9.81m/{\mathrm{s}}^2$ and $h(t)=\frac{1}{2}a{t}^2$ with $a=2m/{\mathrm{s}}^2$, the wire exerts $5.905\mathrm{N}$; at $t=1\mathrm{s}$ it delivers $11.81\mathrm{W}$ to the bead. A vertical normal force is perpendicular to every allowed horizontal virtual displacement, but parallel to part of the actual velocity. Measuring wire motion and bead load would expose precisely the energy transfer that is lost by treating the moving wire as a stationary, workless support.

A virtual displacement is an infinitesimal change of configuration compatible with the instantaneous mechanical constraints, evaluated at fixed time. It is a variation among admissible configurations, rather than necessarily the displacement made during a time step. In generalized coordinates, writeSubstitution into the applied virtual work identifies the generalized forces conjugate to the coordinates.^[ref]

An original moving-constraint example makes the fixed-time condition visible. Consider a particle constrained to the circleFor a virtual displacement, freeze $t$ and therefore $R$; first-order compatibility gives $x\delta x+y\delta y=0$. For an actual motion, differentiation with respect to time instead givesAt $(x,y)=(R,0)$ with $R>0$, every admissible virtual displacement has $\delta x=0$, while an actual velocity must have $\dot{x}=\dot{R}$. If the circle is expanding at $0.020m/s$, the actual radial velocity is nonzero even though every virtual displacement is tangent to the instantaneous circle.

Consequently, an ideal radial constraint force can have zero virtual work but nonzero instantaneous power during this driven motion. Zero work on admissible variations is not the statement that a moving support never transfers energy. The example is a time-dependent holonomic constraint; the variations satisfy its instantaneous linearized condition. In a beam, the same distinction lets one impose zero endpoint virtual displacement while leaving the virtual endpoint slope free, as in the moment-free natural boundary condition for a beam.