Alphabeta Math
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.

Rigid-body dynamics and stability — Examples

1 · Prerequisites

2 · Summary

All numerical formulas on this page use the SI convention of Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top and A steady precession of the heavy symmetric top checked on explicit data: the pivot balance, the two conserved momenta and the reduced energy: t=T/(1 s) is dimensionless for physical elapsed time T; inertia entries are numerical values in kg m2, angular velocities in s−1, angular momenta in kg m2 s−1, and energies in joules. The angle paths φ(t)=t and ψ(t)=2t therefore have physical rates 1 s−1 and 2 s−1, and every trigonometric argument is dimensionless. The escape bound means T≤(1 s)4log⁡100; no measured rates or times are reported.

These examples and counterexamples exercise the definitions and theorems of cm11-rigid-body-dynamics-and-stability on explicit model data. The first example (Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top) declares two inertia matrices, the symmetric diag(1,1,2) and the asymmetric diag(1,2,3), and checks principal-axis spins and their invariants, the explicit non-uniform symmetric-top motion ω(t)=(12cos⁡t,12sin⁡t,1) with conserved angular momentum and energy, and the linearized frequencies of the asymmetric model at its three principal-axis spins. The counterexample (A torque-free motion of a symmetric rigid-material model with constant angular momentum and energy whose body angular velocity is not constant: the invariants do not force a uniform rotation) uses that same explicit motion to refute the claim that constant angular momentum and constant rotational kinetic energy force a constant body angular velocity: the invariants are C=174 and E=98 while ω(π2)≠ω(0).

The escape example (Explicit escape from the intermediate-axis spin of the model with inertia diag(1,2,3): the constants, the perturbed invariant level, and the escape-time bound) instantiates the instability of the intermediate-axis spin of the model with diag(1,2,3): the constants are a=3, c0=14, δ0=14, ε0=12; the perturbed datum with q=1/400 lies on the same invariant level at distance below 10−2 from the spin, and the motion escapes the ball of radius 12 within 4log⁡100≈18.42 time units. The last example (A steady precession of the heavy symmetric top checked on explicit data: the pivot balance, the two conserved momenta and the reduced energy) verifies a steady precession of the heavy-top branch on declared parameters Ac=2, Cc=1, m=d=g=1 (so A=3) and angles φ(t)=t, θ=π/3, ψ(t)=2t: the pivot balance holds identically in time, the conserved momenta are a=72, b=52, the reduced energy is 194, and the heavy-top system holds with constant nutation angle π/3.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

None yet.

5 · Examples, experiments and counterexamples

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top

Physical scope

Worked checks inside the torque-free branch of the rigid-material model of Adoption of the rigid-material model branch: reference mass distribution, rigid placements, orientation and scope on the declared time window I0=R for both models and two declared inertia matrices: a symmetric model with body inertia matrix diag(1,1,2) and an asymmetric model with diag(1,2,3). The checks verify principal-axis spins and their invariants, an explicit non-uniform symmetric-top motion with conserved angular momentum and energy, and the linearized frequencies at the principal-axis spins of the asymmetric model. No measured value, no concrete body and no experiment is asserted; the numerical data are declared model parameters under the explicit numerical SI convention dimensionless t = T / (one second) in the Example.

Example

Numerical SI convention. Let T denote physical elapsed time and fix T∗=1 s, M∗=1 kg and D∗=1 m. Below t:=T/T∗ is dimensionless and all quantities are numerical SI values: inertia is divided by M∗D∗2, angular velocity and its perturbations by T∗−1, angular momentum by M∗D∗2/T∗, and energy by M∗D∗2/T∗2. Primes differentiate with respect to t. Restore physical angular velocity by ωphys(T)=ω(T/T∗)/T∗ and angular momentum by lphys(T)=(M∗D∗2/T∗)l(T/T∗). The displayed numerical rates 1, 2, 1/3 mean physical rates 1/T∗, 2/T∗, 1/(3T∗); λj is restored by λj,phys=λj/T∗2. Thus cos⁡t, sin⁡t, cos⁡(t/3) and the hyperbolic functions have dimensionless arguments; t=π/2 means T=πT∗/2. For the reconstructed motion use xphys(T,ρphys)=D∗x(T/T∗,ρphys/D∗); the chain rule restores RphysTdRphys/dT=[ωphys]× and the dimensioned torque-free equation. This convention uses Dimensions and SI unit conventions of the non-relativistic classical mechanics framework and reports no measurement.

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), inherited for finite mass moments and inertial integrals through the rigid-material and torque-balance suppliers, and for the double-cover surjectivity used in the quaternion constructions below. Consider two declared models of the rigid-material branch Adoption of the rigid-material model branch: reference mass distribution, rigid placements, orientation and scope, each with body inertia matrix diagonal in the declared positively oriented orthonormal basis (e1,e2,e3) and no other inertia data, each with time window I0=R (equivalently all physical times T∈R), each used in the torque-free branch of Adoption of the external-torque branch of the rigid-material model: declared torque data and the balance of angular momentum about the centre of mass clause (iv):

Verify the following.

  1. Principal-axis spins. For each model, each j and each Ω∈R, the constant path ω≡Ωej is the body angular velocity of a torque-free motion; along it l=ΩIjej is constant, C=Ω2Ij2 and E=12Ω2Ij. In Model S with j=3 and Ω=2 this gives l=4e3, C=16, E=4; in Model A with j=1 and Ω=1 it gives l=e1, C=1, E=12.
  2. A non-uniform symmetric-top motion. For Model S the path ω(t)=(12cos⁡t, 12sin⁡t, 1)(t∈R) satisfies the torque-free Euler equations l′=l×(I−1l) and is the body angular velocity of a torque-free motion of Model S; along it l(t)=(12cos⁡t, 12sin⁡t, 2),C=174,E=98,∥ω(t)∥22=54, and ω is not constant: ω(π2)−ω(0)=(−12,12,0).
  3. Linearization of Model A at its principal-axis spins. At the spin Ωe1 with Ω=1 the linearized system of Linearization of the torque-free spin of the rigid-material model at a principal-axis spin: the oscillatory/exponential dichotomy and linear stability has λ1=−13, κ1=1/3, and the pair (δ2,δ3)=(cos⁡(t/3), −13sin⁡(t/3)) solves the transverse equations; at the intermediate spin e2 with Ω=1 the linearized system has λ2=13>0 and the solution with δ1(0)=1, δ3(0)=0 has δ1(t)=cosh⁡(t/3), which is unbounded on [0,∞); at the extreme spin e3 with Ω=1 one has κ3=1.

Facts & Assumptions

Given: The two models S and A of the Example with time window I0=R and their declared diagonal inertia matrices, the torque-free branch, and the data j, Ω of each claim.

[L1]

Model data and nondegeneracy: for each model the body inertia matrix I=∫(∥ρ∥22I3−ρρT) dν is diagonal with the displayed positive entries in the declared basis, and the model is nondegenerate (Inertia matrix, spatial inertia tensor, inertia tensor about a point and moment of inertia about an axis of a rigid-material model clauses (a) and (d)-(e), Principal axes and principal moments of the inertia data of a rigid-material model, axis moments and extremal values clauses (a) and (d), The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

[L2]

Principal-axis spins: for a nondegenerate torque-free model, every axis j, rate Ω and initial orientation there is a torque-free motion with constant body angular velocity Ωej, and along the motions l=ΩIjej, C=Ω2Ij2, E=12Ω2Ij (Principal-axis uniform rotations of the torque-free rigid-material model: existence, a fixed spatial direction, the invariants, and the converse classification clauses 1-3).

[L3]

Torque-free equations and invariants: Lc=Rl, and l=Iω satisfies l′=l×(I−1l), and ⟨l,l⟩ and 12⟨l,I−1l⟩ are constant along solutions; the latter is E=12ωTIω (Torque-free motion of the rigid-material model: conservation of the angular momentum about the centre of mass, the energy-momentum invariants, the principal-axis equations and global existence of the spin dynamics clauses 2-3, Euler's equations with torque for the rigid-material model, the rotational power identity, the torque-free reduction and frame covariance clause 4).

[L4]

Realisability: for every continuous ω:R→R3 there is a C2 rotation path with body angular velocity ω if ω is C1; and a C2 rotation path with any prescribed initial orientation exists by the surjectivity of the quaternion cover (The quaternion lift of a prescribed continuous angular velocity: q′=12q φ(ω(t)) keeps S3 invariant, its rotation path has body angular velocity ω, and the initial value problem has a unique solution on the whole interval clauses 1-2, The quaternion rotation map ρ: conjugation on the imaginary quaternions, the unit sphere S3, and quaternion lifts and its second remark).

[L5]

Linearization: the derivative of the torque-free spin field at Ωej and the dichotomy with λj=Ω2(Ia−Ij)(Ij−Ib)/(IaIb) and κj=∣λj∣, with the explicit cos⁡,sin⁡ and cosh⁡,sinh⁡ solutions (Linearization of the torque-free spin of the rigid-material model at a principal-axis spin: the oscillatory/exponential dichotomy and linear stability clauses 1-2).

[L7]

Quarter-turn values: π=2γ>0 for the smallest positive cosine zero γ, and sin⁡(π/2)=1, cos⁡(π/2)=0 (Pi as twice the smallest positive zero of cosine, Cosine has a smallest positive zero, lying strictly between zero and two, Quarter-turn values and shifts by pi/2 and pi); square roots denote unique nonnegative roots (Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}).

Proof

technique · direct; instantiate the general statements on the declared diagonal data and check the displayed formulas by differentiation and expansion
1.1givenL1L2algebra

By [L2] applied to each model, each j and each Ω there is a torque-free motion with constant body angular velocity Ωej, and l=ΩIjej, C=Ω2Ij2, E=12Ω2Ij along it; for Model S with I3=2, j=3, Ω=2 this gives l=4e3, C=(2⋅2)2=16, E=12⋅4⋅2=4, and for Model A with I1=1, j=1, Ω=1 it gives l=e1, C=1, E=12, which is claim 1.

1.2givenL1L3L6algebra

For the path ω of claim 2 in Model S, l=Iω=(12cos⁡t,12sin⁡t,2) by [L1], so l′=(−12sin⁡t,12cos⁡t,0) by [L6], while l×(I−1l)=l×(12cos⁡t,12sin⁡t,1) has components (12sin⁡t⋅1−2⋅12sin⁡t, −(12cos⁡t⋅1−2⋅12cos⁡t), 12cos⁡t⋅12sin⁡t−12sin⁡t⋅12cos⁡t)=(−12sin⁡t,12cos⁡t,0) by the coordinate formula for the cross product in [L6]; hence l′=l×(I−1l), so l solves the torque-free equations and ω satisfies the Euler equations by [L3].

1.3givenL5L6algebra

For claim 3, [L5] gives λj=Ω2(Ia−Ij)(Ij−Ib)/(IaIb); for Model A and j=1, with a=2, b=3: λ1=(2−1)(1−3)2⋅3=−13 and κ1=1/3; for j=2, with the cyclic choice a=3, b=1: λ2=(3−2)(2−1)3⋅1=13 and κ2=1/3; for j=3, with a=1, b=2: λ3=(1−3)(3−2)1⋅2=−1 and κ3=1. In the first case the linearized field of [L5] satisfies Ae2=−13e3 and Ae3=e2, and the pair δ2(t)=cos⁡(t/3), δ3(t)=−13sin⁡(t/3) satisfies δ2′=δ3 and δ3′=−13δ2 by [L6], so it solves the transverse linearized equations.

2.1step 1.2givenL1L6L7algebra

The quantities of claim 2: ∥ω(t)∥22=14cos⁡2t+14sin⁡2t+1=14+1=54 and C=⟨l,l⟩=14cos⁡2t+14sin⁡2t+4=174 by the Pythagorean identity in [L6]; and E=12ωTIω=12(14cos⁡2t+14sin⁡2t+2⋅1)=12(14+2)=98. The difference ω(π2)−ω(0)=(−12,12,0)≠0 has norm 22, so ω is not constant.

3.1step 1.2step 2.1givenL3L4L6

Since ω is C∞, hence C1, [L4] supplies a C1 rotation path R:R→SO(3) with body angular velocity ω and any prescribed initial orientation; then R′=R[ω]× is C1, so R is of class C2, and x:I0×R3→R3, x(t,ρ):=R(t)ρ (that is, c≡0), is a motion of Model S on its declared window I0=R by Adoption of the rigid-material model branch: reference mass distribution, rigid placements, orientation and scope clause (ii) with body angular velocity ω and l as in step 1.2. Its balance holds: Lc=Rl by [L3], so (Rl)′=R′l+Rl′=R(ω×l+l′)=0 because l′=l×(I−1l)=l×ω=−ω×l from step 1.2 and the alternation of the cross product in [L6]; hence the motion is torque-free and realizes claim 2.

4.1step 1.3givenL5L6algebra∎

In the intermediate case j=2, the same formulas give Ae1=−13e3 and Ae3=−e1, so the transverse system is δ1′=−δ3, δ3′=−13δ1; with δ1(0)=1, δ3(0)=0 the cosh⁡,sinh⁡ solution of [L5] is δ1(t)=cosh⁡(t/3), δ3(t)=−13sinh⁡(t/3), and δ1 is unbounded on [0,∞) by [L6]; hence the linearized system at the intermediate spin of Model A has an unbounded solution, which is claim 3.

Remarks

CounterexampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

A torque-free motion of a symmetric rigid-material model with constant angular momentum and energy whose body angular velocity is not constant: the invariants do not force a uniform rotation

Physical scope

A refutation inside the torque-free branch of the rigid-material model of Adoption of the rigid-material model branch: reference mass distribution, rigid placements, orientation and scope: on the declared symmetric model with body inertia matrix diag(1,1,2) there is a torque-free motion whose angular momentum about the centre of mass and rotational kinetic energy are both constant while its body angular velocity varies with time. The item refutes only the displayed false statement; it asserts no measurement, no concrete body and no failure of the torque-free equations, of which the exhibited motion is an exact solution.

Statement refuted

Let a model of the rigid-material branch Adoption of the rigid-material model branch: reference mass distribution, rigid placements, orientation and scope be used in the torque-free branch of Adoption of the external-torque branch of the rigid-material model: declared torque data and the balance of angular momentum about the centre of mass clause (iv), and let x(t,ρ)=c(t)+R(t)ρ be a motion of the resulting model on a time window containing [0,∞). Suppose that along the motion the angular momentum about the centre of mass Lc and the rotational kinetic energy 12ωTIω are constant in time. Then the body angular velocity ω of the motion is constant in time.

Facts & Assumptions

Use the numerical SI convention of Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top: t=T/T∗ with T∗=1 s, inertia entries in M∗D∗2 with M∗=1 kg and D∗=1 m, ωphys(T)=ω(T/T∗)/T∗ and energy in M∗D∗2/T∗2. Thus t=π/2 means T=πT∗/2; constancy or nonconstancy is preserved by these fixed nonzero scales. Assume ACω (The Axiom of Countable Choice (ACω)), carried from the supplier for its integral and quaternion constructions. No measured data enter this counterexample (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework).

Given: The model S of Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top, namely the rigid-material model with body inertia matrix diag(1,1,2) in the declared basis, used torque-free, and the path ω(t)=(12cos⁡t,12sin⁡t,1).

[F2]

The path ω is the body angular velocity of a torque-free motion of model S; along that motion Lc is constant with C=∥Lc∥22=174 by Torque-free motion of the rigid-material model: conservation of the angular momentum about the centre of mass, the energy-momentum invariants, the principal-axis equations and global existence of the spin dynamics clause 1, and the rotational kinetic energy is the constant E=98; and ω satisfies the torque-free Euler equations l′=l×(I−1l) with l=Iω=(12cos⁡t,12sin⁡t,2) (Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top claim 2, Torque-free motion of the rigid-material model: conservation of the angular momentum about the centre of mass, the energy-momentum invariants, the principal-axis equations and global existence of the spin dynamics clauses 2-3).

[F3]

Values of the path: ω(0)=(12,0,1) and ω(π2)=(0,12,1), using cos⁡0=sin⁡π2=1 and sin⁡0=cos⁡π2=0 with the sine and cosine of Sine and cosine defined by their real power series, the definition of π in Pi as twice the smallest positive zero of cosine, and the exact quarter-turn values of Quarter-turn values and shifts by pi/2 and pi together with the zero values of The derivatives of sine and cosine are cosine and minus sine; the Euclidean norm of the difference is computed with The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn and Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}.

Counterexample

technique · direct; instance the hypotheses of the refuted statement on the explicit model and motion of [Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top](/physics/item/ex-cm11-free-top-checks) and exhibit two times at which the body angular velocity differs
1.1givenF1F2

The motion of [F2] is a torque-free motion of model S on R, its angular momentum about the centre of mass is constant (with ∥Lc∥22=174), and its rotational kinetic energy is constantly 98; hence all hypotheses of the refuted statement are satisfied.

2.1step 1.1givenF3algebra

By [F3], ω(π2)−ω(0)=(−12,12,0), whose Euclidean norm is 22≠0; hence ω(π2)≠ω(0), so the body angular velocity of the motion is not constant in time.

3.1step 1.1step 2.1given∎

The hypotheses of the refuted statement hold along the motion by step 1.1, while its asserted conclusion fails by step 2.1; therefore the statement is false.

Remarks

  • What the counterexample shows and does not show. It shows that conservation of the angular momentum about the centre of mass and of the rotational kinetic energy does not force a constant body angular velocity, so no proof of uniform rotation may rely on those invariants alone. It does not show that the torque-free equations are inconsistent: the exhibited path is an exact solution of them, and the two invariants are precisely the reason the motion is bounded in l and ω.

  • Where uniform rotation does occur. Constant body angular velocity characterises the principal-axis spins, as proved in Principal-axis uniform rotations of the torque-free rigid-material model: existence, a fixed spatial direction, the invariants, and the converse classification clause 4; the present motion is not one of them because ω(t) is not parallel to a fixed principal direction.

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

Explicit escape from the intermediate-axis spin of the model with inertia diag(1,2,3): the constants, the perturbed invariant level, and the escape-time bound

Physical scope

A worked instance of the intermediate-axis instability of Nonlinear stability of the torque-free spin of the rigid-material model about an extreme principal axis, and instability about the intermediate axis, with an explicit escape time and the orientation caveat on the declared model with body inertia matrix diag(1,2,3) and spin datum l0=e2: the instability constants, a perturbed datum on the same invariant level with the checked proximity bound, and the resulting escape time and escape radius. The numerical data are declared model parameters; no measurement, no concrete body and no experiment is asserted.

Example

Use the numerical SI convention of Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top: t=T/(1 s) is dimensionless, inertia entries are in 1 kg m2, and l,L,q,δ0,ε0 and all distances in angular-momentum space are numerical values in 1 kg m2 s−1. Restore the growth rate by c0,phys=c0/(1 s) and the escape time by T=t(1 s); δ0/q is dimensionless. The datum below is model data, not an observation.

Adopt the model A of Two explicit free-rotation models checked: principal-axis spins, a non-uniform symmetric-top motion with its invariants, and the linearized frequencies of an asymmetric top: the rigid-material model of Adoption of the rigid-material model branch: reference mass distribution, rigid placements, orientation and scope with body inertia matrix I=diag(1,2,3) in the declared orthonormal basis, used in the torque-free branch of Adoption of the external-torque branch of the rigid-material model: declared torque data and the balance of angular momentum about the centre of mass clause (iv). Its principal moments are I1=1<I2=2<I3=3 (Principal axes and principal moments of the inertia data of a rigid-material model, axis moments and extremal values clauses (a) and (d)). Take the intermediate-axis spin datum l0:=Le2 with L=1, and let q:=1400 and lq:=(q, 1−4q2, −3 q), so that q=1/400 and 1−4q2=39999/40000. Verify:

  1. The instability constants of Nonlinear stability of the torque-free spin of the rigid-material model about an extreme principal axis, and instability about the intermediate axis, with an explicit escape time and the orientation caveat clause 3.

a=3,c0=14,δ0=14,ε0=12.

  1. The perturbed datum lies on the invariant level of the spin. C(lq)=C(l0)=1 and E(lq)=E(l0)=14, where C(z)=⟨z,z⟩ and E(z)=12⟨z,I−1z⟩ are the invariants of Torque-free motion of the rigid-material model: conservation of the angular momentum about the centre of mass, the energy-momentum invariants, the principal-axis equations and global existence of the spin dynamics clause 2.
  2. The perturbed datum is close to the spin. ∥lq−l0∥2<1100 and q<min⁡{δ0,1100}.
  3. The escape. There is a torque-free motion of model A whose l-path has l(0)=lq and satisfies ∥l(0)−e2∥2<1100,∥l(t)−e2∥2≥12 at some time 0≤t≤4log⁡100≈18.42; in particular the intermediate-axis spin e2 of model A is not Lyapunov stable, with an explicit effective escape at the displayed constants.

Facts & Assumptions

Given: The model A with I=diag(1,2,3) and principal moments I1=1<I2=2<I3=3, the spin datum l0=e2, the number q=1/400 and the datum lq of the statement.

[L2]

Instability statement: for a nondegenerate torque-free model with I1<I2<I3 and L≠0, the spin Le2 is not Lyapunov stable, and with a, c0, δ0, ε0 as in Nonlinear stability of the torque-free spin of the rigid-material model about an extreme principal axis, and instability about the intermediate axis, with an explicit escape time and the orientation caveat clause 3, for every δ>0 and every q∈(0,min⁡{δ0,δ}) whose prescribed datum additionally satisfies ∥lq−Le2∥2<δ, there is a torque-free motion with ∥l(0)−Le2∥2<δ and a time t≤c0−1log⁡(δ0/q) with ∥l(t)−Le2∥2≥ε0; the motion's datum is the l-path with l(0)=lq of the statement (Nonlinear stability of the torque-free spin of the rigid-material model about an extreme principal axis, and instability about the intermediate axis, with an explicit escape time and the orientation caveat clause 3).

[L3]

Invariants and global existence: C and E are constant along every solution of l′=l×(I−1l), which has a unique global maximal solution through every datum (Torque-free motion of the rigid-material model: conservation of the angular momentum about the centre of mass, the energy-momentum invariants, the principal-axis equations and global existence of the spin dynamics clauses 2 and 4).

[L4]

Arithmetic: 3⋅3=3 and the nonnegative square root is unique, 4=2, and for 0≤s≤1 the elementary ordered-field inequalities used below hold; squares of nonnegative reals are ordered as their bases are (Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Intervals of R: the nine order-convex forms, nondegeneracy, and length, Dimensions and SI unit conventions of the non-relativistic classical mechanics framework).

Proof

technique · direct; compute the instability constants from the principal moments, check the invariants and the proximity of the perturbed datum by explicit arithmetic, and instantiate the escape statement
1.1givenL1L4algebra

With I1=1, I2=2, I3=3 and L=1: 1/I1−1/I2=12, 1/I2−1/I3=16, so a2=1/21/6=3 and a=3>0 by [L4]; then c0=a (1/I2−1/I3) 32∣L∣=3⋅16⋅32=312=14 by [L4]. Further 1+a2=1+3=4, so δ0=∣L∣21+a2=12⋅2=14 by [L4] and ε0=∣L∣2=12, which is claim 1.

2.1step 1.1givenL1L3L4algebra

Writing s=1−4q2, one has s2=1−4q2; hence C(lq)=q2+s2+3q2=4q2+1−4q2=1 and E(lq)=12(q2+s22+3q23)=12(q2+1−4q22+q2)=12⋅12=14 by [L4]. For l0=e2: C(l0)=1 and E(l0)=12⟨e2,I−1e2⟩=12⋅12=14. This is claim 2, and by [L3] the perturbed datum and the spin datum lie on the same invariant level.

2.2step 1.1givenL4algebra

For the proximity bound: ∥lq−l0∥22=q2+(s−1)2+3q2=4q2+(1−s)2, and 1−s=1−s21+s=4q21+s≤4q2 because s≥0; with q=1400, 4q2=4160000=140000 and (1−s)2≤(4q2)2=164004=116⋅108<1109, so ∥lq−l0∥22≤140000+1109<110000, whence ∥lq−l0∥2<1100 by [L4] applied to the nonnegative squares. Moreover q=1400<1100=min⁡{δ0,1100} by step 1.1. This is claim 3.

3.1step 1.1step 2.2givenL2L3algebra∎

Apply clause 3 of Nonlinear stability of the torque-free spin of the rigid-material model about an extreme principal axis, and instability about the intermediate axis, with an explicit escape time and the orientation caveat with δ=1100 and q=1400, which is admissible by step 2.2 and [L2]: there is a torque-free motion of model A with ∥l(0)−e2∥2<1100 and a time 0≤t≤c0−1log⁡(δ0/q) with ∥l(t)−e2∥2≥ε0. With the constants of step 1.1, c0−1log⁡(δ0/q)=4log⁡(1/41/400)=4log⁡100 and ε0=12, so the escape-time bound and escape radius of claim 4 follow; by [L3] the motion's l-path is defined for all times, so the escape time is finite as displayed.

Remarks

  • The escape is effective. The initial perturbation has ∥l(0)−e2∥2<10−2 while the escape radius is 12 at a time at most 4log⁡100≈18.421 in the declared units; this is a quantitative instance of the statement that the intermediate-axis spin is not Lyapunov stable, not merely a linearized indication.

  • What is declared. The inertia matrix, the spin datum and the perturbation parameter q are declared model data in the framework's units; no measurement is reported and no concrete body is asserted to have inertia diag(1,2,3).

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

A steady precession of the heavy symmetric top checked on explicit data: the pivot balance, the two conserved momenta and the reduced energy

Physical scope

A worked steady precession of the heavy symmetric top branch of Adoption of the heavy symmetric top branch: axisymmetric rigid material, a fixed pivot material point, a uniform gravitational field and the balance of angular momentum about the pivot on declared parameters Ac=2, Cc=1, m=1, d=1, g=1 (so the pivot transverse moment is A=3) and the explicit angles φ(t)=t, θ=π/3, ψ(t)=2t: the pivot balance is verified identically in time, the two conserved momenta are a=7/2 and b=5/2, the reduced energy is 19/4, and the heavy-top system holds, so the reduction conclusions apply. All formulas use numerical SI values with dimensionless dimensionless t = T / (one second) as defined in the Example; no measurement, no concrete top and no experiment is asserted.

Example

Numerical SI convention. Let T denote physical elapsed time and fix T∗=1 s, M∗=1 kg and D∗=1 m. All formulas below use numerical SI values: t:=T/T∗ is dimensionless, angles are in radians, positions and d are divided by D∗, masses by M∗, inertia entries by M∗D∗2, and g by D∗/T∗2. Primes denote differentiation with respect to t. Restore physical angular velocity as ωphys(T)=ω(T/T∗)/T∗, angular momentum as Jphys(T)=(M∗D∗2/T∗)J(T/T∗), energy as (M∗D∗2/T∗2)E, and torque as (M∗D∗2/T∗2)τ(T/T∗). Thus the displayed rates φ′=1 and ψ′=2 mean physical rates 1 s−1 and 2 s−1, and every trigonometric argument is dimensionless. The physical placement is xphys(T,ρphys)=D∗x(T/T∗,ρphys/D∗); its derivatives restore the pivot balance and heavy-top equations with the same SI data. These are unit conventions of Dimensions and SI unit conventions of the non-relativistic classical mechanics framework, not measured inputs.

Adopt the heavy symmetric top branch Adoption of the heavy symmetric top branch: axisymmetric rigid material, a fixed pivot material point, a uniform gravitational field and the balance of angular momentum about the pivot with declared parameters Ac=2,Cc=1,m=1,d=1,g=1 in the framework's units (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework), so that the body inertia matrix is I=diag(2,2,1) in the declared basis and the pivot transverse moment is A=Ac+md2=2+1=3,C=Cc=1. Let θ0:=π/3 and consider the angles φ(t)=t,θ(t)=θ0,ψ(t)=2t(t∈R), the orientation path R(t):=Θ(φ(t),θ(t),ψ(t)) of The Euler-angle parametrization and the heavy symmetric top system, the centre of mass c(t):=d R(t)e3, and the motion x(t,ρ):=c(t)+R(t)ρ. Verify:

  1. It is a motion of the branch satisfying the pivot balance. The pivot constraint c(t)+R(t)(−de3)=0 holds, and with ω(t)=(−32cos⁡(2t), 32sin⁡(2t), 52) and γ(t)=R(t)Te3, I~ ω′+ω×I~ ω=−mgd (e3×γ)(t∈R), where I~=diag(3,3,1); hence the motion satisfies the adopted balance of Adoption of the heavy symmetric top branch: axisymmetric rigid material, a fixed pivot material point, a uniform gravitational field and the balance of angular momentum about the pivot clause (v) and the conclusions of Reduction of the heavy symmetric top branch: pivotal angular momentum, Euler-angle equations, the two conserved momenta and the one-dimensional θ equation clauses 1-2 hold along it.
  2. The conserved momenta and the reduced energy. With φ′=1, θ′=0, ψ′=2, a=Aφ′sin⁡2θ0+bcos⁡θ0=72,b=C(ψ′+φ′cos⁡θ0)=52,E=U(θ0)=194.
  3. The heavy-top system holds, and θ is constant. The three equations of The Euler-angle parametrization and the heavy symmetric top system hold for the triple, and Aθ′′=−U′(θ0)=0, so the motion is a steady precession: φ′=1 and ψ′=2 are constant and the nutation angle is constantly π/3.

Facts & Assumptions

Given: The branch parameters Ac=2, Cc=1, m=1, d=1, g=1, the derived constants A=3, C=1, the angle data θ0=π/3, φ(t)=t, ψ(t)=2t, and the path R=Θ(φ,θ0,ψ) with c=dRe3.

[L1]

Branch kinematics and dynamics: the motions of the branch satisfy the pivot constraint c=dRe3 and the balance Jp′=τp with Jp=RI~ω, τp=−mgd(Re3)×e3, so the reduction conclusions of Reduction of the heavy symmetric top branch: pivotal angular momentum, Euler-angle equations, the two conserved momenta and the one-dimensional θ equation clauses 1-4 hold on any interval carrying an Euler-angle description (Adoption of the heavy symmetric top branch: axisymmetric rigid material, a fixed pivot material point, a uniform gravitational field and the balance of angular momentum about the pivot clauses (ii)-(v), Reduction of the heavy symmetric top branch: pivotal angular momentum, Euler-angle equations, the two conserved momenta and the one-dimensional θ equation clauses 1-4).

[L2]

Euler-angle kinematics: for R=Θ(φ,θ,ψ) one has ω=(−φ′sin⁡θcos⁡ψ+θ′sin⁡ψ, φ′sin⁡θsin⁡ψ+θ′cos⁡ψ, ψ′+φ′cos⁡θ) and γ=RTe3=(−sin⁡θcos⁡ψ, sin⁡θsin⁡ψ, cos⁡θ) (Angular velocity in Euler angles: ω3=ψ′+φ′cos⁡θ and ω12+ω22=θ′2+φ′2sin⁡2θ, Reduction of the heavy symmetric top branch: pivotal angular momentum, Euler-angle equations, the two conserved momenta and the one-dimensional θ equation clause 2, The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives).

[L3]

The heavy-top system and its reduction: the three displayed equations of The Euler-angle parametrization and the heavy symmetric top system, and the conclusions Aθ′′=−U′(θ), φ′=(a−bcos⁡θ)/(Asin⁡2θ), ψ′=b/C−(a−bcos⁡θ)cos⁡θ/(Asin⁡2θ) of Reduction of the heavy symmetric top system: conserved momenta, one-dimensional θ equation and conserved reduced energy clauses 2-4.

[L5]

The constant π=2γ, where γ is the smallest positive zero of cosine; cos⁡π=−1, cos⁡0=1, cos⁡(π/2)=0, sin⁡0=0; cosine is strictly decreasing on [0,π] and sine strictly increasing on [0,π/2]; sine and cosine obey their addition formulas; and every nonnegative real has a unique nonnegative square root (Pi as twice the smallest positive zero of cosine, Cosine has a smallest positive zero, lying strictly between zero and two, Quarter-turn values and shifts by pi/2 and pi, Signs, monotonicity intervals, and ranges of sine and cosine, The addition formulas for sine and cosine, Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}).

Verification

technique · direct; compute the body angular velocity, the vertical direction in body coordinates and the two sides of the pivot balance on the declared data, then read off the conserved momenta and the reduced energy
1.1L4L5algebra

Put θ=π/3, c=cos⁡θ and s=sin⁡θ. Since 0<θ<π/2, [L5] gives 0<c<1 and s>0. The addition formulas [L5] and Pythagorean identity [L4] give cos⁡π=cos⁡(3θ)=4c3−3c; as cos⁡π=−1, (c+1)(2c−1)2=0. Because c>0, c=1/2. Then s2=1−c2=3/4 and s>0, so the uniqueness of the nonnegative square root [L5] gives s=3/2.

2.1givenstep 1.1L2L4algebra

By [L2] with φ′=1, θ′=0, ψ′=2, θ0=π/3 and ψ(t)=2t: ω(t)=(−sin⁡θ0cos⁡(2t), sin⁡θ0sin⁡(2t), 2+cos⁡θ0)=(−32cos⁡(2t), 32sin⁡(2t), 52) and γ(t)=(−32cos⁡(2t), 32sin⁡(2t), 12) by step 1.1 and [L2].

2.2step 1.1givenL1L4algebra

With φ′=1, θ′=0, ψ′=2, sin⁡2θ0=34 and cos⁡θ0=12 by [L4]: b=C(ψ′+φ′cos⁡θ0)=1⋅(2+12)=52 and a=Aφ′sin⁡2θ0+bcos⁡θ0=3⋅34+52⋅12=94+54=72.

3.1step 2.1givenL1L4algebra

Compute both sides of the balance with I~=diag(3,3,1): by [L4], ω′=(3sin⁡(2t), 3cos⁡(2t), 0), so I~ω′=(33sin⁡(2t), 33cos⁡(2t), 0), and with I~ω=(−332cos⁡(2t), 332sin⁡(2t), 52), the coordinate formula for the cross product gives (ω×I~ω)1=534sin⁡(2t)−1534sin⁡(2t)=−532sin⁡(2t), (ω×I~ω)2=−1534cos⁡(2t)+534cos⁡(2t)=−532cos⁡(2t) and (ω×I~ω)3=0. Hence the left side is (32sin⁡(2t), 32cos⁡(2t), 0). On the right, e3×γ=(−γ2,γ1,0)=(−32sin⁡(2t), −32cos⁡(2t), 0) and mgd=1, so the right side is the same vector; the balance holds for every t.

3.2step 2.2givenL3L4algebra

The reduced energy: E=12Aθ′2+(a−bcos⁡θ0)22Asin⁡2θ0+b22C+mgdcos⁡θ0=(9/4)22⋅3⋅(3/4)+25/42+12=98+258+48=388=194, where a−bcos⁡θ0=72−54=94 by step 2.2.

4.1step 2.1step 3.1givenL1

The constraint c(t)+R(t)(−de3)=dR(t)e3−dR(t)e3=0 holds by the definition of c; hence with step 3.1 the data constitute a motion of the branch satisfying the adopted pivot balance, and clauses 1-2 of Reduction of the heavy symmetric top branch: pivotal angular momentum, Euler-angle equations, the two conserved momenta and the one-dimensional θ equation apply.

5.1step 4.1step 2.2step 3.2givenL3L4algebra∎

The three equations of the heavy-top system of [L3] hold for the triple: the first and third are the statements that a and b of step 2.2 are constant, and the second is 0=Aφ′2sin⁡θ0cos⁡θ0−C(ψ′+φ′cos⁡θ0)φ′sin⁡θ0+mgdsin⁡θ0=3⋅34−534+234 with φ′=1, which vanishes by [L4]. Hence (φ,θ,ψ) is a heavy-top trajectory with parameters (A,C,m,d,g)=(3,1,1,1,1), and Aθ′′=−U′(θ0)=0; the conservation of a and b and the formulas φ′=(a−bcos⁡θ0)/(Asin⁡2θ0)=9/49/4=1 and ψ′=b/C−(a−bcos⁡θ0)cos⁡θ0/(Asin⁡2θ0)=52−9/4⋅1/29/4=52−12=2 of [L3] are consistent with the data, so the motion is a steady precession with constant nutation angle π/3.

Remarks

Sources