Alphabeta Math
ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck pass
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

The equal-arm double pendulum in the linear approximation: the slow same-side mode, the fast opposite-side mode and the bounded modal motions

Physical scope

A worked instance of the coupled quadratic model: the two-angle small-oscillation data M = m l^2 [[2,1],[1,1]], K = m g l diag(2,1) attributed to the linearisation of an equal-mass equal-length planar double pendulum, with its two modal data (the slower same-side mode of frequency sqrt((g/l)(2-sqrt2)) and the faster opposite-side mode of frequency sqrt((g/l)(2+sqrt2))), the general bounded modal motion and the pure-mode data. The example verifies the quadratic model on the declared matrices; it does not derive the pendulum Lagrangians, does not control the error of the linear approximation at finite amplitude and asserts no measured period.

Example

The algebra below uses the fixed-unit numerical convention of The coupled quadratic oscillator model: mass metric, stiffness, generalized force, modal vocabulary and the small-oscillation model at a mechanical equilibrium, with coordinate units radians, time unit 1 s and energy unit 1 J. Physical model parameters are given their stated units; in matrix products, initial-vector identities and coefficient formulas their numerical coefficients in these fixed units are used. Modal amplitudes have common unit a∗=(1 J)1/2(1 s); the j-th component of a mass-normalized modal vector has unit sj/a∗. Modal oscillatory and hyperbolic coefficients have unit a∗, and zero-mode slopes have unit a∗/(1 s). Physical coordinates and energies are reconstructed with the coordinate units and 1 J, respectively. In particular a numerical single-mode datum q(0)=ui means a modal amplitude of one in this common amplitude unit.

Global-time, growth and boundedness claims below refer to the unique maximal extension of the autonomous free equation on R; local motions remain its restrictions. Values at 0 in coefficient formulas are evaluated on that extension when the original interval omits 0.

Let a planar double pendulum consist of two equal particle labels of mass m>0 connected by two equal massless rods of length l>0, pivoted at a fixed point; let θ1,θ2 be the angles of the rods from the downward vertical, read in radians (dimensionless under the angle convention of Dimensions and SI unit conventions of the non-relativistic classical mechanics framework), and let g>0 be the uniform gravitational acceleration of the model. In the linear approximation about the downward equilibrium of the cited treatment, the small-oscillation data are the coupled quadratic model of The coupled quadratic oscillator model: mass metric, stiffness, generalized force, modal vocabulary and the small-oscillation model at a mechanical equilibrium with n=2,

M=ml2(2111),K=mgl(2001),J≡0.

Write γs:=2−2, γf:=2+2. Then:

  1. M is symmetric positive definite (with M1/2 its unique positive definite square root) and K is symmetric positive definite, so the model is oscillatory with both modal stiffnesses positive;
  2. the modal data are λs=glγswithws=(12),λf=glγfwithwf=(1−2), that is Kws=λsMws and Kwf=λfMwf; the M-orthonormal modal vectors are us=ws/2ml2γf and uf=wf/2ml2γs;
  3. the mode frequencies are ωs=gl(2−2)<ωf=gl(2+2), and every motion is bounded on R and has the form q(t)=(Acos⁡(ωst)+Bsin⁡(ωst))us+(Ccos⁡(ωft)+Dsin⁡(ωft))uf;
  4. in the slow mode us the two angles have the same sign (θ2=2 θ1) and in the fast mode uf they have opposite signs (θ2=−2 θ1); the datum q(0)=us, q′(0)=0 excites exactly the slow mode and the datum q(0)=uf, q′(0)=0 exactly the fast mode.

Facts & Assumptions

Given: The parameters m,l,g>0, the matrices M,K of the Example, the vectors ws,wf and the numbers γs,γf.

[F1]

Model dictionary (The coupled quadratic oscillator model: mass metric, stiffness, generalized force, modal vocabulary and the small-oscillation model at a mechanical equilibrium): the data (M,K,0) declare a free coupled quadratic model with two degrees of freedom, kinetic energy 12q′TMq′, potential energy 12qTKq and the modal vocabulary of clause (e).

[F2]

Free-mode theory (The free coupled quadratic model: global existence, the complete modal classification, mode energies, intrinsic degenerate mode subspaces, boundedness and the unstable direction clauses 1, 2, 4 and 5): every free motion is the modal superposition of that theorem; positive definite stiffness gives positive modal stiffnesses and bounded motions; the pure-λ motions are the solutions of q′′+λq=0 with values in the intrinsic subspace Nλ.

[F3]

Direct linear algebra (Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn): M-inner products, matrix products and positive definiteness are checked by the displayed quadratic forms and products. The existence of the positive definite square root is supplied by [F1], not by these arithmetic definitions.

[F4]

Scalar oscillator calculus (The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a), Ck Euclidean maps and diffeomorphisms): for λ>0 the scalar equation z′′+λz=0 has exactly the family Acos⁡(ωt)+Bsin⁡(ωt) with ω=λ.

Verification

Proof technique: direct; check positive definiteness, verify the two eigenpairs by explicit multiplication, and read the motion and boundedness statements from the free-mode theorem.

1.1F1F3algebra

Positive definiteness: M is symmetric and, for w=(x,y)≠0, wTMw=ml2(2x2+2xy+y2)=ml2(x2+(x+y)2)>0, since both squares can vanish only when x=y=0. Likewise wTKw=mgl(2x2+y2)>0, so K is symmetric positive definite; hence both modal stiffnesses are positive and the model is oscillatory.

1.2F1F3algebra

Eigenpairs: multiplying, Kws=mgl(22),Mws=ml2(2+21+2), and 2/(2+2)=2/(1+2)=2−2, so Kws=gl(2−2)Mws; similarly Kwf=mgl(2−2),Mwf=ml2(2−21−2), and 2/(2−2)=−2/(1−2)=2+2, so Kwf=gl(2+2)Mwf. The normalizations give usTMus=ufTMuf=1 and usTMuf=0, because wsTMws=2ml2γf, wfTMwf=2ml2γs and wsTMwf=ml2(2−2)=0; this proves clauses 1 and 2.

2.1F2F4step 1.2

Motions: by [F2] with λs,λf>0 every motion is the displayed two-frequency superposition, with coefficients bijective onto the data, and it is bounded on R because both modal families are bounded; the frequencies are λs and λf by [F4]. This is clause 3.

2.2F2step 1.2algebra

Mode patterns and pure-mode data: reading the coordinates of ws=(1,2) and wf=(1,−2) gives θ2=2 θ1 in the slow mode and θ2=−2 θ1 in the fast mode; the data q(0)=us, q′(0)=0 correspond to A=1, B=C=D=0, hence to q(t)=cos⁡(ωst)us, and analogously for uf. This is clause 4.

2.3step 1.2algebra

Frequency comparison: γf−γs=22>0 gives ωf>ωs: the mode in which the rods swing in opposite directions is the faster one, consistent with the cited treatment's conclusion.

3.1givenstep 1.1step 2.1step 2.2∎

Clauses 1-4 are verified on the declared matrices, which are the small-oscillation data attributed to the double pendulum in the cited linearization. The verification does not derive those matrices from the pendulum Lagrangian, makes no claim about the accuracy of the linear approximation at finite angle, and asserts no measured period; the gravitational acceleration g is a model parameter.

Remarks

  • Why the slow mode is the same-side mode. The potential energy 12qTKq=mgl(θ12+12θ22) penalizes θ1 more strongly than θ2, while the kinetic metric couples the two angles; the resulting slow direction is the one in which both rods move together and the fast direction the one in which they move oppositely. The computation above verifies the frequencies and patterns; the physical reading is the interpretation of the model data.

  • The linear model is not the pendulum. The matrices are the coefficients of the second-order Taylor expansion of the full Lagrangian in displacement and velocity at the downward rest point: M is the zeroth-order value of the kinetic metric at that equilibrium, whereas K is the Hessian of the potential. Derivatives of the kinetic metric multiply terms of total degree at least three in displacement and velocity, so they do not enter this quadratic Lagrangian. The exact pendulum is not a quadratic model, and the finite-time accuracy of the linear approximation is controlled only by Finite-time error control for the small-oscillation approximation: the quadratic remainder of the mechanical field and the variation-of-constants bound on small data, not by this example.

Depends on

Dependency tree · 2 levels

Sources