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

Cm Math Rigid

1 · Prerequisites

2 · Summary

The rigid-body layer of the non-relativistic classical mechanics collection, all of it mathematics: the measure-theoretic inertia tensor of a rigid body, the kinematics of a rigid placement, the torque-free Euler equations with their nonlinear stability and instability statements, and the Euler-angle reduction of the heavy symmetric top. Every item carries domain: mathematics and depends only on mathematical suppliers; nothing here asserts that a physical body is rigid, that its motion obeys any balance law, or that the parameters of the heavy-top system are measured quantities. Those readings belong to the consuming classical-mechanics pages, which add their own postulates, interpretations and empirical qualifications.

Coordinates and regularity used on this page and its companion. The mathematics supplier represents R3 as functions on 3={0,1,2}. Write s0,s1,s2 for its standard basis from The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0, and, in the conventional one-based notation of the rigid-body formulas, set e1:=s0, e2:=s1, e3:=s2 and uj:=u(j−1) for j=1,2,3. For a matrix M∈M3(R), write Mjk:=M(j−1,k−1) for j,k=1,2,3; sums indexed from 1 to 3 use these labels. In particular, define tr⁡M:=∑j=13Mjj. Whenever a curve c:I→R3 is called differentiable, this means componentwise differentiability in the sense of The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral, with derivative whose jth component is cj′. Twice differentiable means each component has a first derivative that is differentiable. The Borel sigma-algebra on R3 below is for its usual Euclidean metric topology: d2 is the metric of Rn as the set of functions n→R, and d1, d2, d∞ are metrics on it, and every closed metric ball is closed, hence Borel, by Arbitrary unions and finite intersections of open sets are open, open balls are open and closed balls are closed. Thus the complement of a support ball in the definition of a bounded mass distribution is a measurable set.

Kinematic foundations. Derivative rules for differentiable matrix-valued curves, and skewness of RTR′ for a curve in SO(3) proves the sum, transpose and product rules for matrix-valued curves and the skew-symmetry of RTR′ for a curve in SO(3); The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]× establishes the bijection [ω]×z=ω×z between skew matrices and vectors and its conjugation law R[ω]×RT=[Rω]×; Rigid placements, body and spatial angular velocity, and material velocities then defines a rigid placement x(t,ρ)=c(t)+R(t)ρ, its body and spatial angular velocities, and its material velocities.

Inertia. Mass distributions on R3: total mass, first moment and inertia matrix defines a mass distribution as a bounded finite Borel measure on R3, with total mass, first moment and inertia matrix Ijk=∫(∥ρ∥22δjk−ρjρk) dν; Finiteness of the moments of a bounded finite Borel measure on R3 proves every displayed integral finite; Symmetry, the quadratic identity ωTIω=∫∥ω×ρ∥22 dν, and the positive-definiteness criterion for the inertia matrix proves symmetry, the identity ωTIω=∫∥ω×ρ∥22 dν and the criterion that I is positive definite exactly when every line through the origin has positive mass off it; and Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case supplies the orthonormal principal axes with positive principal moments in the nondegenerate case, using the real spectral theorem.

Dynamics. Kinetic energy and angular momentum of a rigid placement: K=12m∥c′∥22+12ωTIω and Lc=RIω computes K=12m∥c′∥22+12ωTIω and Lc=RIω for a centred distribution, by integrating the pointwise identities of the placement. Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations deduces Euler's equations Iω′+ω×Iω=RTτ from an explicitly adopted torque balance, and in the torque-free case the equation l′=l×(I−1l) with the conserved quantities ∥l∥22=∑jlj2 and ∑jlj2/(2Ij) and the component equations in a principal basis. Lyapunov stability of spin about the largest and smallest principal axes of the torque-free top proves Lyapunov stability of spin about the largest and smallest principal axes by the invariant certificates C/(2I1)−E and E−C/(2I3) together with compact continuation; Instability of spin about the intermediate principal axis of the torque-free top proves instability of the intermediate axis by an explicit invariant branch on which l1′≥c0l1>0, with a logarithmic escape-time bound; and Linearization of the torque-free Euler system at a principal-axis spin, and the oscillatory/exponential dichotomy records the linearization at a principal-axis spin and the oscillatory/exponential dichotomy, with all solution formulas verified by differentiation and uniqueness.

Euler angles and the heavy top. The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives collects the coordinate rotation matrices Rz,Ry, their composition law and their derivatives; The Euler-angle parametrization and the heavy symmetric top system defines the Euler-angle parametrization Θ(φ,θ,ψ)=Rz(φ)Ry(θ)Rz(ψ) and the heavy-top system by its three displayed second-order equations; Angular velocity in Euler angles: ω3=ψ′+φ′cos⁡θ and ω12+ω22=θ′2+φ′2sin⁡2θ computes the body angular velocity in Euler angles, including ω3=ψ′+φ′cos⁡θ and ω12+ω22=θ′2+φ′2sin⁡2θ; Reduction of the heavy symmetric top system: conserved momenta, one-dimensional θ equation and conserved reduced energy derives the two conserved momenta, the explicit rates φ′,ψ′, the conserved reduced energy and the one-dimensional equation Aθ′′=−U′(θ).

The companion page cm-math-rigid-examples records finite arithmetic checks of the inertia computations, the spin-stability certificates and the heavy-top reduction.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Derivative rules for differentiable matrix-valued curves, and skewness of RTR′ for a curve in SO(3)

Statement

Let I⊆R be an interval with at least two points, let t0∈I and let A,B:I→M3(R) be differentiable at t0, meaning (The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations) that every entry function t↦A(t)ij and t↦B(t)ij is differentiable at t0 in the sense of The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set; write A′(t0) for the matrix with entries A′(t0)ij:=(t↦A(t)ij)′(t0), and similarly for B.

  1. A+B and AT are differentiable at t0, with (A+B)′(t0)=A′(t0)+B′(t0) and (AT)′(t0)=(A′(t0))T.
  2. AB is differentiable at t0, with (AB)′(t0)=A′(t0)B(t0)+A(t0)B′(t0).
  3. If A(t)∈SO(3) for every t∈I (The orthogonal group O(3) and the rotation group SO(3)), then the matrix A(t0)TA′(t0) is skew-symmetric, i.e. (ATA′)T=−ATA′ at t0, and likewise A′(t0)A(t0)T is skew-symmetric. In particular the trace of each of these two matrices vanishes.

Facts & Assumptions

Given: An interval I⊆R with at least two points, a point t0∈I, and maps A,B:I→M3(R) differentiable at t0 in the entrywise sense of the statement.

[L1]

If f,g:I→R are differentiable at t0 and α∈R, then f+g, αf and fg are differentiable at t0 with (f+g)′(t0)=f′(t0)+g′(t0), (αf)′(t0)=αf′(t0) and (fg)′(t0)=f′(t0)g(t0)+f(t0)g′(t0) Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0.

[L2]

M3(R) is the function space of entries indexed by {0,1,2}×{0,1,2} with pointwise sum, so (A+B)ij=Aij+Bij; the transpose is (AT)ij=Aji; and the product is (AB)ij=∑k=02AikBkj The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations, The transpose AT of a matrix, Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes.

[L3]

Transposition is involutive and reverses products: (AT)T=A and (AB)T=BTAT Transpose is linear and involutive, and (AB)T=BTAT.

[L4]

For R∈SO(3) one has RTR=I3, det⁡R=1, R invertible with R−1=RT, and also RRT=I3 The orthogonal group O(3) and the rotation group SO(3), SO(3) is a subgroup of the general linear group, with R−1=RT.

[L5]

A constant function I→R is differentiable at every point of I with derivative 0: its difference quotient at any limit point of I is the constant function 0, whose limit is 0 The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set.

[L6]

The identity matrix I3 has entries 1 for equal indices and 0 otherwise, and the matrix all of whose entries are 0 is the zero matrix Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes, The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations.

Proof

technique · direct
1.1givenL1L2

For all i,j the entry (A+B)ij=Aij+Bij is differentiable at t0 with derivative Aij′(t0)+Bij′(t0) by [L1] and [L2]; these are the entries of A′(t0)+B′(t0), so A+B is differentiable with the asserted derivative.

1.2givenL1L2

For all i,j the entry (AT)ij=Aji is differentiable at t0 with derivative Aji′(t0), which is the (i,j) entry of (A′(t0))T; so AT is differentiable with the asserted derivative.

1.3givenL1L2

For all i,j the entry (AB)ij=∑k=02AikBkj is a finite sum of products of functions differentiable at t0, hence differentiable at t0 with derivative ∑k=02(Aik′Bkj+AikBkj′) by two applications of [L1]; this sum is the (i,j) entry of A′(t0)B(t0)+A(t0)B′(t0) by [L2]. Hence AB is differentiable at t0 with (AB)′=A′B+AB′.

2.1step 1.21.3givenL4L5L6

Assume A(t)∈SO(3) for all t∈I. By steps 1.2 and 1.3 the map t↦A(t)TA(t) is differentiable at t0, with derivative (A′)TA+ATA′ at t0. By [L4] this map is the constant map t↦I3, and every entry of this constant entry function is constant, so by [L5] and [L2] its derivative at t0 is the zero matrix [L6]. Hence (A′)TA+ATA′=0 at t0.

2.2step 1.32.1L3L4algebra

By [L4] and [L3] the map t↦A(t)A(t)T is also the constant map t↦I3, so step 1.3 gives A′AT+A(A′)T=0 at t0, and transposing with [L3] gives (A′AT)T=−A′AT.

3.1step 2.1L3algebra

Transposing the identity of step 2.1 and using [L3] gives ATA′+(ATA′)T=0, hence (ATA′)T=−ATA′, so ATA′ is skew-symmetric.

4.1step 3.12.2L2algebra∎

A skew-symmetric matrix has zero diagonal, hence trace 0: the (i,i) entry satisfies Mii=−Mii by step 3.1 applied entrywise, so 2Mii=0 in R; this applies to the two matrices of steps 3.1 and 2.2.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×

Statement

Use the standard basis (e1,e2,e3) of R3 and coordinates ω=ω1e1+ω2e2+ω3e3 (The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0, A finite list v:n→V is an ordered basis if and only if every x∈V equals ∑i<nλivi for exactly one λ:n→F; those scalars are the coordinates of x in that ordered basis). Define the map [⋅]×:R3→M3(R) by [ω]×:=(0−ω3ω2ω30−ω1−ω2ω10). Then:

  1. [ω]×z=ω×z for all ω,z∈R3, the cross product being that of The cross product in R3.
  2. [ω]×T=−[ω]× for every ω; the map ω↦[ω]× is linear and injective, and its image is exactly the set of skew-symmetric matrices of M3(R). Consequently every skew-symmetric Ω∈M3(R) equals [ω]× for exactly one ω, namely ω=(Ω32,Ω13,Ω21).
  3. For every R∈SO(3) (The orthogonal group O(3) and the rotation group SO(3)) and every ω∈R3, R[ω]×RT=[Rω]×.

Facts & Assumptions

Given: Vectors ω,z∈R3, a skew-symmetric matrix Ω∈M3(R) and a rotation R∈SO(3).

[L1]

The cross product has coordinates (u×v)1=u2v3−u3v2, (u×v)2=u3v1−u1v3, (u×v)3=u1v2−u2v1 in the standard basis, and satisfies ⟨u×v,w⟩=det⁡[u v w], where [u v w] is the matrix with these three vectors as columns; moreover u×v is orthogonal to u and to v The cross product in R3, The cross product is bilinear, alternating, and orthogonal to both factors.

[L2]

Use the owning page's one-based labels: uj:=u(j−1), Mij:=M(i−1,j−1), and ej is the standard unit vector at the actual index j−1, for i,j∈{1,2,3}. Relabelling the supplier's product formula gives (AB)ij=∑k=13AikBkj=∑r=02A(i−1,r)B(r,j−1). Sums are entrywise and (AT)ij=Aji; the transpose is involutive and reverses products The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations, Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes, The transpose AT of a matrix, Transpose is linear and involutive, and (AB)T=BTAT, The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0.

[L4]

For R∈SO(3) one has RTR=I3, RRT=I3, det⁡R=1 and R invertible; for all x,y∈R3 therefore ⟨Rx,Ry⟩=⟨x,y⟩, since (Rx)T(Ry)=xTRTRy=xTy The orthogonal group O(3) and the rotation group SO(3), SO(3) is a subgroup of the general linear group, with R−1=RT, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn.

Proof

technique · direct
1.1givenL1L2algebra

Comparing coordinates, ([ω]×z)1=−ω3z2+ω2z3=ω2z3−ω3z2=(ω×z)1, and likewise ([ω]×z)2=ω3z1−ω1z3=(ω×z)2 and ([ω]×z)3=ω1z2−ω2z1=(ω×z)3, all by [L1] and [L2]; this is clause 1.

1.2givenL2algebra

The displayed matrix is manifestly antisymmetric: ([ω]×T)ij=([ω]×)ji=−([ω]×)ij for every i,j by [L2]. Hence [ω]×T=−[ω]×. Entries are linear functions of ω, so ω↦[ω]× is linear. If [ω]×=0, then the (1,3) entry gives ω2=0, the (2,3) entry gives ω1=0 and the (1,2) entry gives ω3=0, so ω=0: the map is injective.

1.3givenL1L2L3L4L5algebra

Let u,v,w∈R3. By [L4], ⟨R(u×v),Rw⟩=⟨u×v,w⟩=det⁡[u v w] by [L1]; by [L2] and [L3], det⁡[u v w]=det⁡(R)det⁡[u v w]=det⁡(R[u v w])=det⁡[Ru Rv Rw]=⟨Ru×Rv,Rw⟩ by [L1] again. Since R is invertible, Rw ranges over all of R3 as w does, so [L5] gives R(u×v)=Ru×Rv.

2.1step 1.2givenL2algebra

Let Ω be skew-symmetric and put ω:=(Ω32,Ω13,Ω21). Comparing entries with the definition of [ω]× and using Ωij=−Ωji and Ωii=0: [ω]× has (1,2) entry −ω3=−Ω21=Ω12, (1,3) entry ω2=Ω13, (2,3) entry −ω1=−Ω32=Ω23, and the lower-triangular entries follow by skew-symmetry of both matrices; so Ω=[ω]×. Together with step 1.2 this shows that the image of ω↦[ω]× is the set of skew-symmetric matrices, with the stated unique preimage; this proves clause 2.

3.1step 1.11.3givenL2L4L6algebra∎

For every z∈R3, R[ω]×RTz=R(ω×(RTz)) by step 1.1, =(Rω)×(RRTz) by step 1.3, =(Rω)×z by [L4], =[Rω]×z by step 1.1. Hence R[ω]×RT=[Rω]× by [L6]; this is clause 3.

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

Mass distributions on R3: total mass, first moment and inertia matrix

Definition

Let B(R3) be the Borel sigma-algebra of R3 (The Borel sigma-algebra of a topological space) and let ν be a measure on the measurable space (R3,B(R3)) (Measures on sigma-algebras, Measure spaces) that is finite, ν(R3)<+∞, and bounded, meaning that there is a real R>0 with ν(R3∖B‾R(0))=0, the closed ball being that of Open ball, closed ball and sphere in a metric space. Throughout, integrals are Lebesgue integrals against ν in the sense of Integral over a measurable subset and Integrable real and complex functions, and their integrals, and for a measurable vector- or matrix-valued function the integral is defined componentwise: ∫F dν:=(∫F1 dν,∫F2 dν,∫F3 dν) and (∫M dν)jk:=∫Mjk dν, whenever every displayed scalar integral is finite.

A mass distribution on R3 is a bounded finite Borel measure ν with m:=ν(R3)∈(0,+∞); the number m is its total mass. Writing points of R3 as ρ and using the Euclidean inner product ⟨⋅,⋅⟩ and norm ∥ρ∥2=⟨ρ,ρ⟩ of The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, define:

  1. the first moment μ(ν):=∫ρ dν∈R3, read componentwise as μ(ν)j=∫ρj dν;
  2. the inertia matrix I(ν)∈M3(R) (The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations) by its entries I(ν)jk:=∫(∥ρ∥22 δjk−ρjρk) dν, where δjk=1 if j=k and δjk=0 otherwise; equivalently I(ν)=∫(∥ρ∥22I3−ρρT) dν with I3 the identity matrix and ρρT the matrix product of Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes.

The mass distribution is centred when μ(ν)=0. Under the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), each of the integrals above is a finite real number: on B‾R(0) one has ∣ρj∣≤R and ∣ρjρk∣≤R2, and ν(R3)<+∞; the precise statement and proof are Finiteness of the moments of a bounded finite Borel measure on R3 ↗, which depends on this definition only for its notation. The assumption is used there exactly once, to obtain the Borel measurability of the polynomial integrands from the published continuity criterion; every later item that invokes this finiteness carries the assumption explicitly.

Remarks

  • Nothing physical is asserted. "Mass distribution", "total mass" and "inertia matrix" are names for the mathematical objects just defined. No postulate that a physical body is described by such a measure, no identification of ν with a measured mass, and no mechanical law is part of this definition; the physical reading is supplied by the classical-mechanics items that consume it, with their own hypotheses.

  • Generality. The measure need not be absolutely continuous with respect to Lebesgue measure and need not be atomless: finitely many point masses, a surface density, and a volume density are all special cases. The bounded-support hypothesis makes every polynomial moment of degree at most two finite; without it the second moment could diverge.

  • Why the first moment matters. Being centred is a condition on the coordinates in which the measure is expressed: it says that the origin of the body coordinate system has been placed at the centre of mass of the distribution. A physical consumer that uses this definition must state that its body coordinates are so placed; the general translation behaviour of μ and I is not asserted here.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Finiteness of the moments of a bounded finite Borel measure on R3

Statement

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)). Let ν be a measure on (R3,B(R3)) (Mass distributions on R3: total mass, first moment and inertia matrix) that is finite and bounded: ν(R3)<+∞ and ν(R3∖B‾R(0))=0 for some R>0 (Open ball, closed ball and sphere in a metric space). Then:

  1. for all j,k∈{1,2,3} the functions ρ↦ρj and ρ↦ρjρk are Borel measurable and ν-integrable, with ∫∣ρj∣ dν≤R ν(R3) and ∫∣ρjρk∣ dν≤R2 ν(R3);
  2. consequently m:=ν(R3) and the scalar integrals ∫ρj dν and ∫(∥ρ∥22δjk−ρjρk) dν are finite, with ∣∫ρj dν∣≤Rm and ∣∫(∥ρ∥22δjk−ρjρk) dν∣≤2R2m. When m>0, these integrals are precisely the components of the first moment μ(ν) and inertia matrix I(ν) of Mass distributions on R3: total mass, first moment and inertia matrix. When m=0, all these integrals vanish; the zero measure is not called a mass distribution.

Facts & Assumptions

Given: A finite Borel measure ν on R3 and a real R>0 with ν(R3∖B‾R(0))=0.

[L1]

Every continuous map R3→R is Borel measurable, under ACω Continuous functions on Euclidean spaces are Borel measurable, Borel measurable and Lebesgue measurable functions on Rn; and sums, scalar multiples and pointwise products of measurable real-valued functions are measurable Arithmetic and lattice operations preserve measurability whenever they are defined.

[L2]

For a nonnegative simple function the nonnegative Lebesgue integral equals the simple integral; in particular the constant function R on R3 has ∫R dν=R ν(R3) The nonnegative Lebesgue integral, The integral of a nonnegative simple function.

[L3]

The nonnegative integral is monotone and homogeneous for nonnegative scalars: if 0≤f≤g then ∫f dν≤∫g dν, and ∫cf dν=c∫f dν for c≥0 Monotonicity and nonnegative homogeneity of the nonnegative integral.

[L4]

A real measurable function f with ∫∣f∣ dν<+∞ is integrable, and then ∣∫f dν∣≤∫∣f∣ dν Integrable real and complex functions, and their integrals, The modulus of an integral is bounded by the integral of the modulus.

[L5]

ρ=ρ1e1+ρ2e2+ρ3e3 satisfies ∣ρj∣≤∥ρ∥2 for each j and ∥ρ∥2≤R exactly on B‾R(0) (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Open ball, closed ball and sphere in a metric space).

Proof

technique · direct
1.1givenL1

Each coordinate function ρ↦ρj is continuous, hence Borel measurable by [L1]; the products ρ↦ρjρk are then measurable by the product clause of [L1].

2.1step 1.1givenL1L2L3L5L6algebra

Let N:=R3∖B‾R(0), a Borel set by [L6], with ν(N)=0. We justify the almost-everywhere comparison from the definition of the nonnegative integral: if a nonnegative measurable f satisfies f≤B on Nc for a finite constant B≥0, every nonnegative simple minorant s≤f has the same simple integral as s1Nc. Indeed, in a finite disjoint level-set representation of s, the contribution of each intersection with N is zero, while ν(A)=ν(A∖N) for every measurable level set A. Pointwise s1Nc≤B, so [L2] and [L3] give ∫s dν≤Bν(R3); taking the supremum over simple minorants gives ∫f dν≤Bν(R3) by The nonnegative Lebesgue integral. Apply this with f=∣ρj∣, B=R, and with f=∣ρjρk∣, B=R2: measurability follows from step 1.1 and [L1], and [L5] gives the bounds on Nc. Thus ∫∣ρj∣ dν≤Rν(R3) and ∫∣ρjρk∣ dν≤R2ν(R3).

3.1step 2.1givenL4

The bounds of step 2.1 are finite because ν is finite, so each of the functions ρj and ρjρk has finite integral of its modulus, hence is ν-integrable by [L4].

4.1step 2.13.1givenL1L4L5algebra∎

Each coordinate integrand is integrable by step 3.1. Each inertia integrand gjk(ρ):=∥ρ∥22δjk−ρjρk is measurable by [L1], and on Nc its modulus is at most 2R2 by [L5]. The simple-minorant comparison proved in step 2.1 gives ∫∣gjk∣ dν≤2R2m<∞, so [L4] makes gjk integrable and gives ∣∫gjk dν∣≤2R2m. Likewise [L4] and step 2.1 give ∣∫ρj dν∣≤Rm. If m>0, ν is a mass distribution and these are its moment and inertia components by Mass distributions on R3: total mass, first moment and inertia matrix. If m=0, the same bounds force every displayed integral to be zero; this case is excluded from the definition's mass-distribution terminology. Thus clause 2 holds for every measure in the hypotheses.

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

Rigid placements, body and spatial angular velocity, and material velocities

Definition

Let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let ν be a mass distribution on R3 with total mass m (Mass distributions on R3: total mass, first moment and inertia matrix), let c:I→R3 be differentiable and let R:I→SO(3) be differentiable (The orthogonal group O(3) and the rotation group SO(3)), differentiability meaning entrywise differentiability of R in the sense of The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set with derivative R′; write c′ for the derivative of c.

The rigid placement with reference shape ν, translation c and rotation R is the map x:I×R3⟶R3,x(t,ρ):=c(t)+R(t)ρ, so that the material point labelled by ρ (a point of the reference shape) is carried to x(t,ρ) at time t; the map x(t,⋅) is the placement at time t, and the pair (c(t),R(t)) is its configuration. The material velocity of the label ρ at time t is the derivative ∂tx(t,ρ):=c′(t)+R′(t)ρ, which exists by differentiability of c and of the entries of R.

Body and spatial angular velocity. By Derivative rules for differentiable matrix-valued curves, and skewness of RTR′ for a curve in SO(3) the matrix R(t)TR′(t) is skew-symmetric at every t∈I. By The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]× there is therefore a unique vector ω(t)∈R3 with R(t)TR′(t)=[ω(t)]×, and the map ω:I→R3 is the body angular velocity of the placement; the spatial angular velocity is ωs(t):=R(t)ω(t). Explicitly ω(t)=((RTR′)32,(RTR′)13,(RTR′)21)(t), a finite sum of products of entries of R and R′; consequently, if R is twice differentiable, then ω is differentiable and its derivative is the body angular acceleration.

Equivalently, by the conjugation clause of The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×, R′(t)R(t)T=[ωs(t)]×.

Remarks

  • The reference shape is a label space. ρ is not the position of the material point at time t; the position is x(t,ρ). The mass distribution ν is a measure on the label space, and no invariance of ν under translations or rotations is assumed.

  • Existence and uniqueness of ω are discharged, not assumed. Skew-symmetry of RTR′ is the differentiated orthogonality relation of the curve R; the vector whose cross-product matrix it is exists and is unique by the explicit bijection of the two cited lemmas. No smoothness beyond differentiability of the entries of R is used to define ω.

  • This is mathematics. No physical law is asserted: not that a physical body is rigid, not that its configuration is described by a curve in SO(3), and not that ω is measured by any apparatus. The physical reading is supplied by the consuming classical-mechanics items.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Symmetry, the quadratic identity ωTIω=∫∥ω×ρ∥22 dν, and the positive-definiteness criterion for the inertia matrix

Statement

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), used only in [L1], through Finiteness of the moments of a bounded finite Borel measure on R3, to make the entries of I finite real numbers. Let ν be a mass distribution on R3 with total mass m=ν(R3) and inertia matrix I=I(ν) (Mass distributions on R3: total mass, first moment and inertia matrix). Write ωTIω:=∑j,k=13ωjIjkωk for the quadratic form of I on R3 and ∥⋅∥2 for the Euclidean norm (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn). Then:

  1. I is symmetric: Ijk=Ikj for all j,k.
  2. For every ω∈R3, ωTIω=∫∥ω×ρ∥22 dν(ρ) ≥0, the cross product being that of The cross product in R3.
  3. The symmetric bilinear form (u,v)↦uTIv on R3 is therefore positive semidefinite (Positive and negative definiteness, the inertia (p,q,r), rank p+q, and signature p−q of a real symmetric bilinear or quadratic form); for each ω≠0, writing Rω={tω:t∈R}, one has ωTIω=0 if and only if ν(R3∖Rω)=0. Thus it is positive definite if and only if ν(R3∖Rω)>0 for every ω≠0, that is, no line through the origin carries the entire mass (Measure-null sets and almost-everywhere statements relative to a measure).
  4. If the form is positive definite, then I is invertible.

Facts & Assumptions

Given: A mass distribution ν with inertia matrix I and a vector ω∈R3.

[L1]

Ijk=∫(∥ρ∥22δjk−ρjρk) dν is a finite real number for every j,k, by the finiteness lemma for bounded finite measures Mass distributions on R3: total mass, first moment and inertia matrix, Finiteness of the moments of a bounded finite Borel measure on R3.

[L2]

The integral is linear on L1(ν): ∫(αf+βg) dν=α∫f dν+β∫g dν for real α,β and integrable f,g The Lebesgue integral is linear on L1(μ).

[L3]

For a measurable h≥0 with ∫h dν=0 one has h=0 ν-almost everywhere; and ∫h dν≥0 for h≥0 A nonnegative measurable function has integral 0 exactly when it vanishes almost everywhere, Monotonicity and nonnegative homogeneity of the nonnegative integral.

[L4]

For all u,v∈R3, ∥u×v∥22=∥u∥22∥v∥22−⟨u,v⟩2, and u×v=0 if and only if u and v are linearly dependent The squared cross-product norm is the Gram determinant of two vectors.

[L6]

The matrix I satisfies N(I)={u:Iu=0}; if N(I)={0} then I is invertible Invertible matrix theorem: invertibility, full pivot rank, RREF I, trivial nullspace and unique solvability are equivalent.

Proof

technique · direct
1.1givenL1algebra

Interchanging j and k leaves the integrand ∥ρ∥22δjk−ρjρk unchanged, because δjk=δkj and ρjρk=ρkρj; the two integrals defining Ijk and Ikj therefore coincide, so I is symmetric.

1.2givenL1L2algebra

Expressing the quadratic form through the entries of [L1] and applying linearity [L2] to the finite sum of integrable functions ωjωk(∥ρ∥22δjk−ρjρk), ωTIω=∑j,kωjωk∫(∥ρ∥22δjk−ρjρk)dν=∫(∑j,kωjωk(∥ρ∥22δjk−ρjρk))dν.

2.1step 1.2givenL4L5algebra

For each ρ, evaluating the δ sum gives ∑j,kωjωk(∥ρ∥22δjk−ρjρk)=∥ω∥22∥ρ∥22−⟨ω,ρ⟩2=∥ω×ρ∥22 by [L4] and bilinearity of the inner product.

3.1step 1.22.1givenL3

Substituting step 2.1 into step 1.2 gives the identity of clause 2, and the integrand ∥ω×ρ∥22 is nonnegative, so its integral is ≥0 by [L3]: the form uTIv is positive semidefinite by definition.

4.1step 3.1givenL1L2L3L4

Let ω≠0 and put h(ρ):=∥ω×ρ∥22. This is a measurable nonnegative integrable polynomial by [L1], [L2] and step 2.1. Its zero set is the closed line Rω, by [L4], so its exceptional set {h≠0}=R3∖Rω is Borel measurable. By step 3.1 and both implications of [L3], ωTIω=0 exactly when h=0 ν-almost everywhere, meaning that its measurable exceptional set has measure zero, or equivalently ν(R3∖Rω)=0. Since the quadratic form is nonnegative and finite, it is strictly positive for this ω exactly when that complement has positive measure. Quantifying over all ω≠0 gives the positive-definiteness criterion in clause 3.

5.1step 4.1givenL6algebra∎

Assume the form is positive definite. If Iω=0 for some ω, then ωTIω=ωT0=0, which for ω≠0 is excluded; hence N(I)={0} and I is invertible by [L6]. This is clause 4.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Kinetic energy and angular momentum of a rigid placement: K=12m∥c′∥22+12ωTIω and Lc=RIω

Statement

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), used only through the bound of Finiteness of the moments of a bounded finite Borel measure on R3 that makes the two integrands of clause 1 ν-integrable. Let J⊆R be an interval with at least two points, let ν be a centred mass distribution on R3 with total mass m and inertia matrix I (Mass distributions on R3: total mass, first moment and inertia matrix), let c:J→R3 and R:J→SO(3) be differentiable, and let x(t,ρ)=c(t)+R(t)ρ be the rigid placement with body angular velocity ω (Rigid placements, body and spatial angular velocity, and material velocities). For t∈J define K(t):=12∫∥∂tx(t,ρ)∥22 dν(ρ),Lc(t):=∫(x(t,ρ)−c(t))×(∂tx(t,ρ)−c′(t)) dν(ρ), the second integral read componentwise (Mass distributions on R3: total mass, first moment and inertia matrix). Then for every t∈J:

  1. ∂tx(t,ρ)=c′(t)+R(t)(ω(t)×ρ) for every ρ∈R3, and R′(t)=R(t)[ω(t)]×=[R(t)ω(t)]×R(t); moreover the two integrands displayed above are ν-integrable, so K(t) and Lc(t) are a well-defined finite real number and vector;
  2. K(t)=12m∥c′(t)∥22+12ω(t)TIω(t);
  3. Lc(t)=R(t)Iω(t).

Facts & Assumptions

Given: A centred mass distribution ν, differentiable c:J→R3 and R:J→SO(3), the rigid placement x(t,ρ)=c(t)+R(t)ρ, its body angular velocity ω with RTR′=[ω]×, and a time t∈J.

[L1]

The body angular velocity satisfies RTR′=[ω]×, hence R′=R[ω]×, and by conjugation R′=[Rω]×R as well; explicitly ω=((RTR′)32,(RTR′)13,(RTR′)21) Rigid placements, body and spatial angular velocity, and material velocities, The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×.

[L3]

The elements of SO(3) preserve the Euclidean inner product and norm: ⟨Ru,Rv⟩=⟨u,v⟩, so ∥Rv∥2=∥v∥2 The orthogonal group O(3) and the rotation group SO(3), The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Real and complex inner product spaces, with the inner product linear in the first argument.

[L4]

The measure ν has finite total mass, ν(R3∖B‾R(0))=0 for some R>0, and vanishing first moment ∫ρj dν=0 for j=1,2,3; every polynomial of degree at most two in the coordinates is ν-integrable, with the bound of Finiteness of the moments of a bounded finite Borel measure on R3 (Mass distributions on R3: total mass, first moment and inertia matrix).

[L5]

The integral is linear on L1(ν) and equals the componentwise integral for vector- or matrix-valued integrands The Lebesgue integral is linear on L1(μ), Integral over a measurable subset, Integrable real and complex functions, and their integrals.

[L7]

For u,v∈R3, the cross product is bilinear, u×v is determined componentwise by (u×v)1=u2v3−u3v2 and its cyclic analogues, and ∥u×v∥22=∥u∥22∥v∥22−⟨u,v⟩2≤∥u∥22∥v∥22 The cross product in R3, The squared cross-product norm is the Gram determinant of two vectors.

Proof

technique · direct
1.1givenL1L2

By [L1], R′=R[ω]×, and [ω]×ρ=ω×ρ by [L2]; hence for every ρ the derivative of s↦x(s,ρ) at t is ∂tx(t,ρ)=c′(t)+R′(t)ρ=c′(t)+R(t)(ω(t)×ρ).

1.2givenL7algebra

For every ρ,ω one has the vector identity ρ×(ω×ρ)=(∥ρ∥22I3−ρρT)ω: for instance the first coordinate of the left side is ρ2(ω1ρ2−ω2ρ1)−ρ3(ω3ρ1−ω1ρ3)=ω1(ρ22+ρ32)−ρ1(ω2ρ2+ω3ρ3), which equals the first coordinate ∥ρ∥22ω1−ρ1⟨ρ,ω⟩ of the right side, and the second and third coordinates are the cyclic analogues.

2.1step 1.1givenL3L4L7algebra

The first integrand satisfies ∥∂tx(t,ρ)∥22≤4∥c′(t)∥22+4∥ω(t)∥22∥ρ∥22 by step 1.1, [L3] and the two elementary bounds ∥a+b∥22≤2∥a∥22+2∥b∥22 and ∥u×v∥2≤∥u∥2∥v∥2; the second integrand satisfies ∥(x−c)×(∂tx−c′)∥2≤∥ρ∥2∥ω(t)×ρ∥2≤∥ω(t)∥2∥ρ∥22 by step 1.1 and [L3]. Both are bounded on the support ball of ν by constants depending on t, so they are ν-integrable by [L4]; this is the integrability clause of the statement.

2.2step 1.1L3L7algebra

Expanding the square with bilinearity and symmetry of the inner product and using [L3], ∥∂tx∥22=∥c′∥22+2⟨c′,R(ω×ρ)⟩+∥ω×ρ∥22.

2.3step 1.1L2algebra

By step 1.1 the relative velocity is (x−c)′=∂tx−c′=R(ω×ρ), so the angular-momentum integrand is (x−c)×(x−c)′=Rρ×R(ω×ρ)=R(ρ×(ω×ρ)) by [L2].

3.1step 2.12.2givenL4L5L6algebra

Integrating step 2.2 with [L5] and pulling the constant R(t) and the linear map ρ↦ω×ρ through the integral, ∫∥∂tx∥22 dν=m∥c′∥22+2⟨c′,R(ω×∫ρ dν)⟩+∫∥ω×ρ∥22 dν. The middle term vanishes because ∫ρ dν=0 by [L4], and the last term equals ωTIω by [L6]; substituting into the definition of K gives clause 2.

4.1step 2.31.2givenL4L5∎

Integrating step 2.3 and substituting step 1.2, Lc(t)=R(t)∫(∥ρ∥22I3−ρρT)dν ω(t)=R(t)Iω(t), the constant matrix R(t) and the finite sum of integrals being interchanged by [L5] componentwise, and the last equality being the definition of I in Mass distributions on R3: total mass, first moment and inertia matrix. This is clause 3.

CorollaryStatement: AI-adaptedProof: AI-adaptedOpen item page →

Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case

Statement

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), used through Symmetry, the quadratic identity ωTIω=∫∥ω×ρ∥22 dν, and the positive-definiteness criterion for the inertia matrix to ensure finiteness of the inertia integrals. Let ν be a mass distribution on R3 with inertia matrix I=I(ν) (Mass distributions on R3: total mass, first moment and inertia matrix), and write Iω for the matrix-vector product and ⟨⋅,⋅⟩ for the Euclidean inner product on R3 (Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn). Then:

  1. there is an orthonormal basis (e1,e2,e3) of R3 (Orthogonal vectors and subspaces, orthogonal and orthonormal sets, and orthonormal bases) and real numbers I1,I2,I3 with Iej=Ijej for j=1,2,3;
  2. for any such basis and any ω=∑jωjej one has Iω=∑jIjωjej, Ij=ejTIej, and ωTIω=∑jIjωj2=∫∥ω×ρ∥22 dν(ρ) ≥0;
  3. each Ij vanishes exactly when ν(R3∖Rej)=0, that is, the whole mass lies on the line Rej; if the inertia form is positive definite then every Ij>0 and I is invertible with I−1ej=Ij−1ej; and if ν is not supported on any line through the origin, meaning that ν(R3∖Rω)>0 for every ω≠0, then every Ij>0.

Facts & Assumptions

Given: A mass distribution ν with inertia matrix I, an orthonormal basis (e1,e2,e3) of eigenvectors of I with eigenvalues I1,I2,I3, and a vector ω=∑jωjej.

[L1]

I is symmetric and ωTIω=∫∥ω×ρ∥22 dν≥0; for each ω≠0, its quadratic form vanishes exactly when ν(R3∖Rω)=0, and the associated bilinear form is positive definite exactly when ν(R3∖Rω)>0 for every ω≠0, and then Iω≠0 for ω≠0 Symmetry, the quadratic identity ωTIω=∫∥ω×ρ∥22 dν, and the positive-definiteness criterion for the inertia matrix.

[L2]

The map ω↦Iω is self-adjoint with respect to the Euclidean inner product, because ⟨Iω,z⟩=(Iω)Tz=ωTITz=ωTIz=⟨ω,Iz⟩ by symmetry of I and the product and transpose rules Symmetry, the quadratic identity ωTIω=∫∥ω×ρ∥22 dν, and the positive-definiteness criterion for the inertia matrix, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn.

Proof

technique · direct
1.1givenL1L2

The map Tω:=Iω is linear and, by [L1] and [L2], self-adjoint on the Euclidean space R3.

2.1step 1.1L3

Applying [L3] to T produces an orthonormal basis (e1,e2,e3) and reals Ij with Iej=Ijej; this is clause 1.

3.1step 2.1givenL1L4algebra

For ω=∑jωjej, linearity of T [L4] gives Iω=∑jIjωjej and hence ωTIω=∑j,kωjωkIj⟨ej,ek⟩=∑jIjωj2 by orthonormality; setting ω=ej gives Ij=ejTIej=∫∥ej×ρ∥22 dν by [L1], and the displayed chain of equalities of clause 2 follows.

4.1step 3.1L1algebra

Each Ij=∫∥ej×ρ∥22 dν is nonnegative, and it is zero exactly when ν(R3∖Rej)=0, by [L1]; this is the vanishing criterion of clause 3.

4.2step 3.14.1givenL4L5algebra

If the inertia form is positive definite, then Ij=ejTIej>0 for every j; moreover Iω≠0 for every ω≠0 by [L1], so the nullspace of T is {0} and I is invertible with unique inverse by [L5]. Define S on the basis by Sej:=Ij−1ej and extend linearly; then SIej=SejIj=Ij⋅Ij−1ej=ej and ISej=I(Ij−1ej)=Ij−1Iej=ej for every j, so SI=IS=Id by [L4], whence S=I−1 by uniqueness and I−1ej=Ij−1ej.

5.1step 4.1givenalgebra∎

If ν is not supported on any line through the origin, then for ej≠0 ν(R3∖Rej)>0, so Ij=∫∥ej×ρ∥22 dν>0 by step 4.1. This completes clause 3.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations

Statement

Assume the Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), used through the inertia suppliers to ensure finiteness of the integrals. Let J⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let ν be a centred mass distribution on R3 that is not supported on any line through the origin, with inertia matrix I (Mass distributions on R3: total mass, first moment and inertia matrix), let c:J→R3 be differentiable and let R:J→SO(3) be twice differentiable, and let x(t,ρ)=c(t)+R(t)ρ be the rigid placement with body angular velocity ω (Rigid placements, body and spatial angular velocity, and material velocities). Write l:=Iω:J→R3; by Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case the inertia form is positive definite and I is invertible.

Adopt the balance hypothesis (B): there is a map τ:J→R3 such that the curve t↦R(t)l(t) is differentiable and (Rl)′(t)=τ(t)(t∈J). Then:

  1. l is differentiable and Euler's equations with torque hold: Iω′+ω×Iω=RTτon J;
  2. if τ=0 on J, then l satisfies the torque-free Euler equations l′=l×(I−1l), and the functions t↦⟨l(t),l(t)⟩ and t↦12⟨l(t),I−1l(t)⟩=12ω(t)TIω(t) are constant on J;
  3. if τ=0 on J and (e1,e2,e3) is a positively oriented orthonormal eigenbasis of I with Iej=Ijej and l=∑jljej, then ⟨l,l⟩=∑jlj2 and 12⟨l,I−1l⟩=∑jlj2/(2Ij) are constant on J, and the torque-free equations read l1′=(1I3−1I2)l2l3,l2′=(1I1−1I3)l1l3,l3′=(1I2−1I1)l1l2.

Facts & Assumptions

Given: The data of the statement, the body angular velocity ω with RTR′=[ω]×, the angular momentum l=Iω, and a map τ satisfying the balance hypothesis (B).

[L1]

The body angular velocity satisfies RTR′=[ω]×, so R′=R[ω]× and R′(t)R(t)T=[Rω]×; it is differentiable because it is a finite sum of products of entries of R and R′ Rigid placements, body and spatial angular velocity, and material velocities, The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×.

[L4]

I is symmetric and positive definite, hence invertible; and for the eigenbasis Iej=Ijej with Ij>0 one has I−1ej=Ij−1ej Symmetry, the quadratic identity ωTIω=∫∥ω×ρ∥22 dν, and the positive-definiteness criterion for the inertia matrix, Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case.

[L5]

The cross product is bilinear and alternating, and ⟨u×v,w⟩=det⁡[u v w], so u×v is orthogonal to u and to v The cross product is bilinear, alternating, and orthogonal to both factors, The cross product in R3.

Proof

technique · direct
1.1givenL1L3

By [L1] the angular velocity ω is differentiable, hence so is l=Iω; the product Rl satisfies (Rl)′=R′l+Rl′ by [L3], and R′l=R[ω]×l=R(ω×l) by [L1] and [L3].

2.1step 1.1givenL2L3L7algebra

By the balance hypothesis (B) and step 1.1, τ=R(ω×l+l′). Multiplying by RT and using [L2] and linearity of matrix-vector multiplication gives RTτ=ω×l+l′. Since l=Iω with I a constant matrix, l′=Iω′ by the finite-sum and scalar rules of [L7], and ω×l=ω×Iω; therefore Iω′+ω×Iω=RTτ, which is clause 1.

3.1step 2.1givenL4L5algebra

Assume τ=0. Then step 2.1 gives Iω′=−ω×Iω; multiplying by I−1 from [L4] gives ω′=−I−1(ω×Iω), hence l′=Iω′=−(ω×Iω)=(Iω)×ω=l×(I−1l), using antisymmetry of the cross product [L5] and I−1l=ω [L4]. This is the first part of clause 2.

4.1step 3.1givenL4L5L6L7algebra

The two functions of clause 2 are constant: the product rule and symmetry of I−1 in [L7] give ddt⟨l,l⟩=2⟨l′,l⟩=2⟨l×(I−1l),l⟩ and ddt12⟨l,I−1l⟩=⟨l′,I−1l⟩=⟨l×(I−1l),I−1l⟩. Both inner products vanish because a cross product is orthogonal to each of its factors [L5]. Also ddt12⟨l,I−1l⟩=ddt12ωTIω because I−1l=ω and I is symmetric [L4], [L6]. These functions are continuous on the solution interval J and differentiable at every interior point, so their zero derivatives make them constant by A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant.

5.1step 4.1givenL4L6algebra

For the eigenbasis, l=∑jljej gives I−1l=∑j(lj/Ij)ej by [L4]; by orthonormality ⟨l,l⟩=∑jlj2 and 12⟨l,I−1l⟩=∑jlj2/(2Ij), and these are constant by step 4.1.

6.1step 5.1givenL2L5algebra∎

Let Q have columns (e1,e2,e3) in the standard basis. The positive orientation and orthonormality give Q∈SO(3); [L2] gives Q(u×v)=(Qu)×(Qv) from the skew/cross correspondence. Thus the standard cross-product coordinate formula applies in this principal basis. Coordinate-wise, l′=l×(I−1l) reads l1′=l2(l3/I3)−l3(l2/I2)=(1/I3−1/I2)l2l3, l2′=l3(l1/I1)−l1(l3/I3)=(1/I1−1/I3)l1l3 and l3′=l1(l2/I2)−l2(l1/I1)=(1/I2−1/I1)l1l2, using the coordinate formula for the cross product [L5] and step 5.1. An arbitrary orthonormal eigenbasis can be made positively oriented by negating one eigenvector; neither its eigenvalues nor the coordinate-free motion/invariant claims change. This is clause 3.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Lyapunov stability of spin about the largest and smallest principal axes of the torque-free top

Statement

Let I1<I2<I3 be positive reals and let I:R3→R3 be the linear map with Iej=Ijej on the standard orthonormal basis (e1,e2,e3) of R3. Consider the autonomous first-order system l′(t)=F(l(t)),F(l):=l×(I−1l), on U:=R3, read as in Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system with the cross product of The cross product in R3 and I−1 the inverse of I (Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case); this is the torque-free Euler system of Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations in angular-momentum variables. Let L≠0. Then:

Le1 and Le3 are equilibria of F and each is Lyapunov stable in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system: for j=1 and j=3 and for every ε>0 there is δ∈(0,ε] such that every maximal solution l with ∥l(0)−Lej∥2<δ is defined for all t≥0 and satisfies ∥l(t)−Lej∥2<ε for every t≥0.

Facts & Assumptions

Given: Positive reals I1<I2<I3, the diagonal map I with Iej=Ijej, the field F(l)=l×(I−1l) on R3, and L≠0.

[L1]

The system z′=F(z) on the open set U=R3 is read as First-order systems, initial value problems, and solutions on intervals with D=R×U; Le1 and Le3 are equilibria exactly when F(Lej)=0, and then Lyapunov stability is the quantified ball statement of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system.

[L2]

The cross product is bilinear and alternating and satisfies ∥u×v∥2≤∥u∥2∥v∥2 The cross product in R3, The squared cross-product norm is the Gram determinant of two vectors, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn.

[L3]

F is locally Lipschitz in the state variable in the sense of Local Lipschitz continuity in the state variable, locally uniform in time and parameters. For ∥x∥2,∥y∥2≤r, bilinearity gives F(x)−F(y)=(x−y)×(I−1x)+y×(I−1(x−y)) [L2]. Since I−1 is diagonal with entries at most 1/I1 in absolute value, ∥I−1z∥2≤∥z∥2/I1; hence ∥F(x)−F(y)∥2≤(2r/I1)∥x−y∥2. A compact time-state cylinder is a compact subset of Euclidean R4, hence bounded by Heine-Borel (Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line, clause 2); choosing r to contain its state projection gives the finite Lipschitz constant required by the definition. The same local estimate gives continuity of F at each state point by the ε-δ definition. Picard-Lindelof therefore applies to the initial value problem of the system Picard-Lindelöf local existence and uniqueness for first-order systems; every datum has a unique maximal solution on an open interval, and every other solution through that datum is a restriction (Extensions, maximal solutions, maximal intervals, and global solutions, Every Picard–Lindelöf initial value problem has one maximal solution on an open interval).

[L4]

Along every solution of l′=l×(I−1l) the quantities C(l)=⟨l,l⟩ and E(l)=12⟨l,I−1l⟩ are constant; for l=∑jljej they equal ∑jlj2 and ∑jlj2/(2Ij) The same algebra as in Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations applies directly to the present diagonal map: I−1 is symmetric, so C′=2⟨l×(I−1l),l⟩=0 and E′=⟨l×(I−1l),I−1l⟩=0 by cross-product orthogonality. Both scalar functions are continuous and differentiable on the solution interval; A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant therefore makes them constant there. These computations require only the given positive diagonal entries, and no realization as a mass-distribution inertia matrix.

[L5]

Any solution of an autonomous system on U=R3 that is trapped near a finite endpoint in a compact subset of U extends past that endpoint: the no-extension hypothesis of A non-extendable solution leaves every compact subset of the domain near its finite endpoint; an autonomous solution trapped in a compact subset of its state domain cannot break down there fails under such trapping.

[L6]

The Euclidean inner product is bilinear and positive definite, so ⟨u,u⟩=∥u∥22 and u≠0 implies ∥u∥2>0 Real and complex inner product spaces, with the inner product linear in the first argument, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn.

[L8]

The Euclidean norm satisfies the triangle inequality, so ∣∥x∥2−∥y∥2∣≤∥x−y∥2; hence the norm is continuous and each sphere {z:∥z∥2=r} is closed (Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation).

Proof

technique · direct
1.1givenL1L2L3

By antisymmetry of the cross product, F(Le1)=Le1×(L/I1)e1=0 and F(Le3)=Le3×(L/I3)e3=0, so Le1 and Le3 are equilibria; and [L3] supplies the existence, uniqueness and continuation theory for the system.

1.2givenL2L4L6L7L8

Let l be the maximal solution with prescribed datum at 0. By [L4] the numbers C(l(t))=∥l(t)∥22=∥l(0)∥22 and E(l(t)) are constant, so the values of l lie in the sphere {z:∥z∥2=∥l(0)∥2}⊆U. This sphere is closed by [L8] and bounded, hence compact by Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line.

1.3givenL4algebra

Define, for l=∑jljej, the quantity D(l):=C(l)/(2I1)−E(l)=12∑j>1(1/I1−1/Ij)lj2. Since 1/I1−1/Ij≥1/I1−1/I2 for j=2,3 by I2<I3 and I1<I2, one has D(l)≥c(l22+l32) with c:=12(1/I1−1/I2)>0.

2.1step 1.2givenL5

Every maximal solution of the system is defined for all t≥0: if a solution were defined on [0,b) with 0<b<+∞ and admitted no extension past b, then by step 1.2 it would be trapped in the compact subset {z:∥z∥2=∥l(0)∥2} of U near b, and [L5] would give an extension past b, a contradiction.

3.1step 2.1step 1.3givenL4L6algebra

Let ε>0 and let c,κ:=1/2I1c be as in step 1.3; choose δ>0 satisfying the finite list of conditions δ≤ε, κ2δ2≤ε2/4, κ2δ2≤∣L∣2/4, δ≤∣L∣/4 and δ(2∣L∣+δ)+κ2δ2≤∣L∣ε/2; such a δ exists because each condition only bounds δ from above by a positive number. Let l be the maximal solution with ∥l(0)−Le1∥2<δ, defined for all t≥0 by step 2.1. Then l2(0)2+l3(0)2≤∥l(0)−Le1∥22<δ2 and, by [L4] and step 1.3 evaluated at t=0, for every t≥0 l2(t)2+l3(t)2≤D(l(t))c=D(l(0))c≤l2(0)2+l3(0)22I1c<κ2δ2≤ε24.

4.1step 3.1givenL4L6algebra

With the data of step 3.1: l1(t)2=C(l(0))−l2(t)2−l3(t)2≥(∣L∣−δ)2−κ2δ2≥∣L∣2/4 by δ≤∣L∣/4 and κ2δ2≤∣L∣2/4, so ∣l1(t)∣≥∣L∣/2 and l1 has the sign of L (it is continuous, ∣l1∣≥∣L∣/2 throughout, and ∣l1(0)−L∣<δ≤∣L∣/4). Hence ∣l1(t)+L∣≥3∣L∣/2 and ∣l1(t)−L∣=∣l1(t)2−L2∣∣l1(t)+L∣≤∣C(l(0))−L2∣+l2(t)2+l3(t)23∣L∣/2≤δ(2∣L∣+δ)+κ2δ23∣L∣/2≤ε3.

5.1step 2.1step 3.1step 4.1givenL1algebra

Combining steps 3.1 and 4.1, ∥l(t)−Le1∥22=(l1(t)−L)2+l2(t)2+l3(t)2<ε29+ε24<ε2 for all t≥0; since for every ε>0 such a δ was produced and every maximal solution starting in Bδ(Le1) is defined for all t≥0 by step 2.1, Le1 is Lyapunov stable in the sense of [L1].

6.1step 1.3step 3.1step 4.1step 5.1givenL4algebra∎

The same argument with the roles of e1 and e3 exchanged uses D3(l):=E(l)−C(l)/(2I3)=12∑j<3(1/Ij−1/I3)lj2≥c3(l12+l22) with c3:=12(1/I2−1/I3)>0, and D3(l)≤(l12+l22)/(2I1). Set κ3:=1/2I1c3; the same positive δ conditions of step 3.1, with κ3 in place of κ, bound l1(t)2+l2(t)2<κ32δ2. The norm invariant then gives ∣l3(t)∣≥∣L∣/2 with the sign of L, and the calculation of step 4.1 gives ∣l3(t)−L∣≤ε/3; the certificate D3 is constant along solutions by [L4], so the argument of steps 1.3, 2.1, 3.1, 5.1 and 4.1 applies verbatim to Le3 and proves its Lyapunov stability.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Instability of spin about the intermediate principal axis of the torque-free top

Statement

Let I1<I2<I3 be positive reals, let I:R3→R3 be the linear map with Iej=Ijej on the standard orthonormal basis, and consider the torque-free Euler system l′=F(l), F(l)=l×(I−1l), on U=R3 (Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations, Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system). Let L≠0. Then l0:=Le2 is an equilibrium of F and is not Lyapunov stable: with ε0:=∣L∣2,a:=1/I1−1/I21/I2−1/I3>0,c0:=a(1I2−1I3)32∣L∣>0, for every δ>0 there is a solution l with ∥l(0)−l0∥2<δ and some t≥0 with ∥l(t)−l0∥2≥ε0; Moreover, with δ0:=∣L∣/(21+a2), the escape can be arranged at a time t≤c0−1log⁡(δ0/η) where η=l1(0)>0 may be chosen as small as desired.

Facts & Assumptions

Given: Positive reals I1<I2<I3, the diagonal map I, the field F(l)=l×(I−1l) on U=R3, and L≠0.

[L1]

The system l′=F(l) is read as First-order systems, initial value problems, and solutions on intervals. For ∥x∥2,∥y∥2≤r, bilinearity gives F(x)−F(y)=(x−y)×(I−1x)+y×(I−1(x−y)); since I−1 is diagonal with entries at most 1/I1 in absolute value and ∥u×v∥2≤∥u∥2∥v∥2, this yields ∥F(x)−F(y)∥2≤(2r/I1)∥x−y∥2. This bound makes F locally Lipschitz in the state variable and also continuous at each state point by the ε-δ definition (Local Lipschitz continuity in the state variable, locally uniform in time and parameters, The cross product in R3, The squared cross-product norm is the Gram determinant of two vectors). Picard-Lindelof gives local existence and uniqueness, hence a unique maximal solution through every datum (Picard-Lindelöf local existence and uniqueness for first-order systems, Extensions, maximal solutions, maximal intervals, and global solutions, Every Picard–Lindelöf initial value problem has one maximal solution on an open interval); a solution trapped in a compact subset of U=R3 near a finite endpoint extends past it (First-order systems, initial value problems, and solutions on intervals, A non-extendable solution leaves every compact subset of the domain near its finite endpoint; an autonomous solution trapped in a compact subset of its state domain cannot break down there).

[L2]

Along every solution the invariants C(l)=∥l∥22 and E(l)=12⟨l,I−1l⟩ are constant; in the eigenbasis Iej=Ijej one has C=∑jlj2, E=∑jlj2/(2Ij), and the component equations l1′=(1/I3−1/I2)l2l3, l2′=(1/I1−1/I3)l1l3, l3′=(1/I2−1/I1)l1l2 The same algebra as in Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations applies directly to the present diagonal map: I−1 is symmetric, so C′=2⟨l×(I−1l),l⟩=0 and E′=⟨l×(I−1l),I−1l⟩=0 by cross-product orthogonality. Both scalar functions are continuous and differentiable on the solution interval; A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant therefore makes them constant there. These computations require only the given positive diagonal entries, and no realization as a mass-distribution inertia matrix.

[L3]

The Euler field is even: F(−l)=(−l)×(I−1(−l))=F(l), by bilinearity and linearity (The cross product in R3). Put σ:=L/∣L∣∈{−1,1}, so σ2=1 and L=σ∣L∣.

[L8]

For r>0, log⁡r is the inverse value of the strictly increasing exponential: exp⁡(log⁡r)=r (The natural logarithm as the inverse of the exponential function, The exponential is a continuous bijection from R onto (0,∞), The exponential function is strictly increasing).

[L7]

The Euclidean norm satisfies the triangle inequality, so ∣∥x∥2−∥y∥2∣≤∥x−y∥2; hence each sphere {z:∥z∥2=r} is closed (Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation).

Proof

technique · direct
1.1givenL1L3

F(Le2)=Le2×(L/I2)e2=0, so l0 is an equilibrium. Put σ=L/∣L∣ as in [L3]. Choose the second and third coordinates of the initial datum with signs σ and −σ, respectively.

2.1step 1.1givenL2L3L4algebra

Let 0<η0<δ0:=∣L∣/(21+a2) and let 0<η<η0. Put s(η):=∣L∣2−(1+a2)η2>0 and lη:=(η, σs(η), −σaη), so that C(lη)=L2 and E(lη)=L2/(2I2): indeed C(lη)=η2+∣L∣2−(1+a2)η2+a2η2=L2 and E(lη)=12(η2/I1+(∣L∣2−(1+a2)η2)/I2+a2η2/I3)=12L2/I2+12η2(1/I1−1/I2−a2(1/I2−1/I3))=L2/(2I2) by the definition of a. Also ∥lη−l0∥2→0 as η↓0, since s(η)→∣L∣.

3.1step 2.1givenL1L2L4algebraL6L7

Let l be the unique maximal solution with l(0)=lη; it is defined for all t≥0, because its invariant sphere C=L2 is closed by [L7] and bounded, hence compact by Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line, and [L1] rules out breakdown inside a compact set. Let τ:=inf⁡{t>0:l1(t)≥δ0 or l1(t)≤0}, so 0<τ≤+∞, and on (0,τ) one has 0<l1<δ0. On that interval l1(t)2+l3(t)2≤l12+a2l12<δ02(1+a2)=L2/4 by [L2] (which gives l3=±al1 on the invariant level set E=L2/(2I2), since E−C/(2I2)=0 and (1/I2−1/I3)>0), so l2(t)2=L2−l12−l32>3L2/4 and l2 has the constant sign σ on [0,τ), its initial sign. Likewise l3=±al1 cannot change sign on [0,τ) without vanishing, which would force l1=l3=0 contrary to l1>0; since l3(0)=−σaη, one has l3(t)=−σal1(t) there.

4.1step 3.1L2algebra

On (0,τ), the first component equation of [L2] and step 3.1 give l1′=(1/I3−1/I2)l2l3=a(1/I2−1/I3)(σl2)l1 ≥c0l1 with c0=a(1/I2−1/I3)32∣L∣>0, because σl2≥32∣L∣.

5.1step 3.14.1givenL5L8algebra

By step 4.1 and the product and chain rules [L5], (e−c0tl1(t))′=e−c0t(l1′(t)−c0l1(t))≥0 on (0,τ); since e−c0tl1(t) is continuous on [0,τ), it is nondecreasing there by [L5], so l1(t)≥ηec0t for t∈[0,τ). If τ=+∞, this would give l1(t)≥ηec0t>δ0 for t>c0−1log⁡(δ0/η), contradicting l1<δ0 on (0,τ); hence τ<+∞.

6.1step 3.15.1givenL4algebra

At the time τ one has l1(τ)=δ0: indeed l1(τ)≠0 by step 5.1 and continuity, so by the definition of τ and continuity it is the alternative l1(τ)≥δ0 that occurs, and on [0,τ] the bound l1<δ0 holds on (0,τ) with l1 continuous, so l1(τ)=δ0; moreover τ≤c0−1log⁡(δ0/η). Therefore ∥l(τ)−l0∥22≥l1(τ)2+l3(τ)2=(1+a2)δ02=L2/4, so ∥l(τ)−l0∥2≥∣L∣/2=ε0.

7.1step 2.16.1givenalgebra∎

It remains to make the initial point as close to l0 as prescribed: for any δ>0 choose η>0 with η<min⁡{δ0,δ} and ∥lη−l0∥2<δ, possible because ∥lη−l0∥2→0 by step 2.1. The corresponding solution then has ∥l(0)−l0∥2<δ and, by step 6.1, escapes the ball of radius ε0 at time τ≤c0−1log⁡(δ0/η). This is instability of l0 in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system, with the explicit time bound of the statement.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Linearization of the torque-free Euler system at a principal-axis spin, and the oscillatory/exponential dichotomy

Statement

Let I1,I2,I3>0 be pairwise distinct, let I:R3→R3 be the linear map with Iej=Ijej on the standard orthonormal basis, and let I−1 denote its inverse, given by I−1ej=Ij−1ej. Define the formal torque-free Euler vector field F(ω):=−I−1(ω×Iω),ω∈R3. When I is an actual rigid-body inertia map written in a principal basis, this is the angular-velocity vector field obtained from Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations. Here the positive diagonal entries are otherwise arbitrary; no realization by a mass distribution is assumed. Let Ω≠0. Then:

  1. F(Ωe3)=0, and for every δ∈R3 F(Ωe3+δ)=Aδ+R(δ),(Aδ)1=Ω(I2−I3)I1δ2,(Aδ)2=Ω(I3−I1)I2δ1,(Aδ)3=0, where ∥R(δ)∥2≤C∥δ∥22 for the explicit constant C:=∥I−1∥F∥I∥F and ∥M∥F:=(∑j,kMjk2)1/2; consequently A is the total derivative of F at Ωe3 (The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder).
  2. Every solution of the linear system δ′=Aδ on R has δ3 constant, and with λ:=Ω2(I2−I3)(I3−I1)I1I2≠0 the pair (δ1,δ2) is determined by its initial value as follows:
    • if λ<0, then with κ:=−λ>0, δ1(t)=δ1(0)cos⁡(κt)+Ω(I2−I3)I1κδ2(0)sin⁡(κt),δ2(t)=δ2(0)cos⁡(κt)+Ω(I3−I1)I2κδ1(0)sin⁡(κt), and every solution is bounded on [0,∞);
    • if λ>0, then the same two formulas hold with cosh⁡ and sinh⁡ in place of cos⁡ and sin⁡ and κ:=λ, the system has unbounded solutions (for instance δ2(0)=0, δ1(0)≠0), and if δ1(0)+Ω(I2−I3)I1κδ2(0)≠0 then ∥δ(t)∥2→+∞ as t→+∞.
  3. Hence the linearization at the spin Ωe3 is oscillatory exactly when (I2−I3)(I3−I1)<0, i.e. when the third principal moment is a largest or a smallest moment, and it has exponentially growing solutions exactly when λ>0, i.e. when I3 lies strictly between I1 and I2; in the latter case the zero solution of the linear system is not stable.

Facts & Assumptions

Given: Pairwise distinct positive reals I1,I2,I3, the diagonal map I, the field F(ω)=−I−1(ω×Iω) and a real Ω≠0.

[L1]

Every x∈R3 has a unique expansion x=∑j=13xjej by the standard-basis result (The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0). Define S(x):=∑j=13(xj/Ij)ej; it is linear and direct substitution gives SI=IS=Id⁡, so S is the inverse linear map of I (Linear map between vector spaces over the same field, Invertible linear maps, linear isomorphisms, and inverse linear maps). The cross product is bilinear with coordinate formula (u×v)1=u2v3−u3v2 and cyclic analogues, and ∥u×v∥2≤∥u∥2∥v∥2 (The cross product in R3, The squared cross-product norm is the Gram determinant of two vectors).

[L2]
[L3]

A is a linear map with Aδ=(Ω(I2−I3)δ2/I1, Ω(I3−I1)δ1/I2, 0); a linear map is totally differentiable at every point with itself as total derivative, and a bound ∥F(ω0+δ)−F(ω0)−Aδ∥2≤C∥δ∥22 exhibits A as the total derivative by the definition The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder.

[L4]

The first-order linear system z′=Az is locally Lipschitz in the state variable and has unique solutions with prescribed initial data; the same holds for z′=Az+b Local Lipschitz continuity in the state variable, locally uniform in time and parameters, Picard-Lindelöf local existence and uniqueness for first-order systems, First-order systems, initial value problems, and solutions on intervals.

[L6]

cosh⁡t=(et+e−t)/2 and sinh⁡t=(et−e−t)/2; the exponential is strictly increasing with image (0,+∞) The six hyperbolic functions and their natural domains, The exponential function is strictly increasing, The exponential is a continuous bijection from R onto (0,∞).

Proof

technique · direct
1.1givenL1algebra

Expanding with [L1], (Ωe3+δ)×I(Ωe3+δ)=Ωe3×Iδ+ΩI3 δ×e3+δ×Iδ, and e3×Iδ=I1δ1e2−I2δ2e1 while δ×e3=δ2e1−δ1e2; collecting components gives (Ωe3+δ)×I(Ωe3+δ)=Ω(I3−I2)δ2e1+Ω(I1−I3)δ1e2+δ×Iδ. Applying −I−1 and [L1] yields F(Ωe3+δ)=Aδ−I−1(δ×Iδ) with A as displayed, and F(Ωe3)=A⋅0=0.

1.2givenL4algebra

Let δ solve δ′=Aδ. The third equation gives δ3′=0, so δ3 is constant; the first two equations read δ1′=cδ2, δ2′=dδ1 with c=Ω(I2−I3)/I1≠0 and d=Ω(I3−I1)/I2≠0 by distinctness and Ω≠0, so δ1′′=cd δ1=λδ1 and δ2=δ1′/c; conversely every solution of x′′=λx produces the solution δ=(x,x′/c,const).

2.1step 1.1givenL2L3

The remainder satisfies ∥I−1(δ×Iδ)∥2≤∥I−1∥F∥δ∥2∥Iδ∥2≤∥I−1∥F∥I∥F∥δ∥22 by [L1] and two applications of [L2]; hence R(δ):=−I−1(δ×Iδ) obeys the bound of clause 1 with C=∥I−1∥F∥I∥F, and A is the total derivative of F at Ωe3 by [L3].

2.2step 1.2givenL4L5algebra

Suppose λ<0 and put κ=−λ. The function x(t):=x(0)cos⁡(κt)+(x′(0)/κ)sin⁡(κt) satisfies x′′=−κ2x=λx by the derivative formulas and chain rule [L5]; the pair (x,x′) solves the first-order system (x,x′)′=(x′,λx), which is locally Lipschitz, so by uniqueness [L4] it is the solution with the given initial data. Taking x=δ1 and using δ2=δ1′/c, δ1′(0)=cδ2(0), gives δ1(t)=δ1(0)cos⁡(κt)+(c δ2(0)/κ)sin⁡(κt) and δ2(t)=δ2(0)cos⁡(κt)−(κ/c)δ1(0)sin⁡(κt); substituting c=Ω(I2−I3)/I1, κ2=−Ω2(I2−I3)(I3−I1)/(I1I2) and κ/c=−Ω(I3−I1)/(I2κ) gives the two displayed formulas. Boundedness follows from ∣sin⁡∣,∣cos⁡∣≤1 [L5].

2.3step 1.2givenL4L5algebra

Suppose λ>0 and put κ=λ. The same argument with cosh⁡,sinh⁡ in place of cos⁡,sin⁡ verifies the two displayed formulas: x(t)=x(0)cosh⁡(κt)+(x′(0)/κ)sinh⁡(κt) satisfies x′′=κ2x by [L5], and uniqueness [L4] identifies it as the solution.

3.1step 2.3givenL6algebra

Let λ>0 and set A0:=δ1(0), B0:=c δ2(0)/κ, so that δ1(t)=A0cosh⁡(κt)+B0sinh⁡(κt)=12(A0+B0)eκt+12(A0−B0)e−κt by [L6]. If A0+B0≠0 then ∣δ1(t)∣≥12∣A0+B0∣eκt−12∣A0−B0∣e−κt→+∞ as t→+∞, because eκt→+∞ and e−κt→0 [L6]; since A0+B0=δ1(0)+cδ2(0)/κ, this is the growth criterion of clause 2, and the choice δ2(0)=0≠δ1(0) satisfies it, so the system has unbounded solutions and clause 3 follows.

4.1step 2.22.33.1givenalgebra∎

Finally the sign of λ is the sign of (I2−I3)(I3−I1) because Ω2/(I1I2)>0; this vanishes exactly when I3 equals I1 or I2, which distinctness excludes, so λ≠0, and λ<0 exactly for extreme, λ>0 exactly for intermediate I3, completing the sign criterion. In the intermediate case, for every δ>0 choose 0<η<min⁡{δ,1} and initial data ηe1. Step 3.1 gives an unbounded solution with initial norm η; by continuity it reaches norm 1 at some t>0. Thus the fixed ε=1 violates the stability condition for every δ, proving that the zero solution is not Lyapunov stable (Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system) and completing clause 3.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives

Statement

For a∈R define, in the standard basis of R3, Rz(a):=(cos⁡a−sin⁡a0sin⁡acos⁡a0001),Ry(a):=(cos⁡a0sin⁡a010−sin⁡a0cos⁡a)∈M3(R). Then:

  1. Rz(a)∈SO(3) and Ry(a)∈SO(3) for every a (The orthogonal group O(3) and the rotation group SO(3)), Rz(a)e3=e3 and Ry(a)e2=e2, and Rz(0)=Ry(0)=I3.
  2. Rz(a)Rz(b)=Rz(a+b) and Ry(a)Ry(b)=Ry(a+b) for all a,b∈R.
  3. For every differentiable a:I→R on an interval I⊆R (The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set, Intervals of R: the nine order-convex forms, nondegeneracy, and length), the matrix-valued curves t↦Rz(a(t)) and t↦Ry(a(t)) are differentiable with (Rz∘a)′=a′ Rz(a) [e3]×=a′ [e3]×Rz(a),(Ry∘a)′=a′ Ry(a) [e2]×=a′ [e2]×Ry(a), where [e3]×,[e2]× are the cross-product matrices of The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×.

Facts & Assumptions

Given: Reals a,b, matrices Rz(a),Ry(a) as displayed, and a differentiable a:I→R.

[L1]

sin⁡ and cos⁡ are differentiable with (sin⁡)′=cos⁡, (cos⁡)′=−sin⁡, and satisfy the addition formulas sin⁡(a+b)=sin⁡acos⁡b+cos⁡asin⁡b, cos⁡(a+b)=cos⁡acos⁡b−sin⁡asin⁡b and the Pythagorean identity sin⁡2a+cos⁡2a=1; also sin⁡(−a)=−sin⁡a and cos⁡(−a)=cos⁡a The derivatives of sine and cosine are cosine and minus sine, The addition formulas for sine and cosine, Parity and the Pythagorean identity for sine and cosine, Sine and cosine defined by their real power series.

[L3]

The determinant is multiplicative and det⁡I3=1; for a block-triangular matrix the determinant is the product of the determinant of the leading block and the remaining diagonal entry For same-sized finite square matrices over a commutative ring, det⁡(AB)=det⁡(A)det⁡(B), For n≥1, the determinant over a commutative ring by the Leibniz formula, and ∣det⁡A∣ for a real matrix.

[L5]

[e3]× and [e2]× are the skew matrices of ω=e3 and ω=e2; in particular [e3]×=(0−10100000) and [e2]×=(001000−100) The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×, The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations.

Proof

technique · direct
1.1givenL1L2L3algebra

Entrywise, Rz(a)TRz(a) has diagonal entries cos⁡2a+sin⁡2a=1, 1, 1 and off-diagonal entries ∓cos⁡asin⁡a±sin⁡acos⁡a=0; hence Rz(a)∈O(3), and expanding the determinant along the third row gives det⁡Rz(a)=det⁡(cos⁡a−sin⁡asin⁡acos⁡a)=(cos⁡a)(cos⁡a)−(−sin⁡a)(sin⁡a)=1, so Rz(a)∈SO(3), with Rz(a)e3=e3 read off the third column.

1.2givenL1L2L3algebra

Similarly Ry(a)TRy(a)=I3 and det⁡Ry(a)=det⁡(cos⁡asin⁡a−sin⁡acos⁡a)=1, so Ry(a)∈SO(3), and Ry(a)e2=e2 is its second column. Rz(0)=Ry(0)=I3 because cos⁡0=1 and sin⁡0=0.

1.3givenL1L2algebra

The (1,1) entry of Rz(a)Rz(b) is cos⁡acos⁡b−sin⁡asin⁡b=cos⁡(a+b), the (1,2) entry is −cos⁡asin⁡b−sin⁡acos⁡b=−sin⁡(a+b), the (2,1) entry is sin⁡acos⁡b+cos⁡asin⁡b=sin⁡(a+b), the (2,2) entry is −sin⁡asin⁡b+cos⁡acos⁡b=cos⁡(a+b), and the remaining entries are 0,0,1; hence Rz(a)Rz(b)=Rz(a+b) by [L1].

1.4givenL1L2algebra

The analogous four planar entries of Ry(a)Ry(b) equal cos⁡(a+b),0,sin⁡(a+b) and 1 in the corresponding positions, so Ry(a)Ry(b)=Ry(a+b) by [L1].

1.5givenL1L4L5algebra

If a:I→R is differentiable, entrywise differentiation using [L1] and [L4] gives (Rz∘a)′=a′(−sin⁡a−cos⁡a0cos⁡a−sin⁡a0000). Direct multiplication of a′Rz(a) with [e3]× from [L5] gives the same matrix: its (1,1),(1,2),(2,1),(2,2) entries are respectively −a′sin⁡a,−a′cos⁡a,a′cos⁡a,−a′sin⁡a, and all remaining entries are zero. Thus (Rz∘a)′=a′Rz(a)[e3]×.

1.6givenL1L4L5algebra

The same comparison for Ry(a(t)): the derivative matrix has entries −a′sin⁡a at (1,1), a′cos⁡a at (1,3), −a′cos⁡a at (3,1), −a′sin⁡a at (3,3) and 0 in the second row and column, which is exactly a′Ry(a)[e2]× by [L5]; hence (Ry∘a)′=a′Ry(a)[e2]×.

2.1step 1.51.6givenL1L2L5algebra∎

Since Rz(a) commutes with [e3]× and Ry(a) commutes with [e2]×, as direct entrywise multiplication using [L1] shows, the two alternative forms a′[e3]×Rz(a) and a′[e2]×Ry(a) of clause 3 hold as well.

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

The Euler-angle parametrization and the heavy symmetric top system

Definition

Euler-angle parametrization. Let Rz,Ry be the coordinate rotation matrices of The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives. The Euler-angle parametrization is the map Θ:R×(0,π)×R⟶SO(3),Θ(φ,θ,ψ):=Rz(φ) Ry(θ) Rz(ψ), which is SO(3)-valued by The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives and SO(3) is a subgroup of the general linear group, with R−1=RT. The numbers (φ,θ,ψ) are the Euler angles and the open set R×(0,π)×R is its domain; θ is restricted to the open interval (0,π) so that sin⁡θ>0, with π defined as twice the smallest positive cosine zero (Pi as twice the smallest positive zero of cosine) and positivity supplied by Pi is the first positive zero of sine. No injectivity and no chart or diffeomorphism property of Θ is asserted or used here, and no periodicity of the angles is imposed — the real parameters stand for themselves.

Parameters. Fix reals A>0, C>0, m>0, d>0 and g>0; they are mathematical parameters, not identified with any physical quantity in this item.

The heavy-top function. Write q=(φ,θ,ψ) for a point of the domain and v=(φ′,θ′,ψ′) for a triple of reals. The heavy-top function is L(q,v):=12A(θ′2+φ′2sin⁡2θ)+12C(ψ′+φ′cos⁡θ)2−mgdcos⁡θ, using the sine and cosine of Sine and cosine defined by their real power series. Repeated application of the sine/cosine derivative identities (The derivatives of sine and cosine are cosine and minus sine) and the finite sum, product and chain rules (Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c)) gives continuous partial derivatives of every order, so it is C∞ on the domain times R3 (Ck Euclidean maps and diffeomorphisms); its relevant partial derivatives are, explicitly, ∂φ′L=Aφ′sin⁡2θ+Ccos⁡θ(ψ′+φ′cos⁡θ),∂φL=0, ∂θ′L=Aθ′,∂θL=Aφ′2sin⁡θcos⁡θ−C(ψ′+φ′cos⁡θ)φ′sin⁡θ+mgdsin⁡θ, ∂ψ′L=C(ψ′+φ′cos⁡θ),∂ψL=0.

The heavy-top system. Let I⊆R be an interval with at least two points. A heavy-top trajectory for these parameters is a twice-differentiable curve q=(φ,θ,ψ):I→R×(0,π)×R such that, with q′=(φ′,θ′,ψ′) its derivative, the three second-order equations ddt(Aφ′sin⁡2θ+Ccos⁡θ(ψ′+φ′cos⁡θ))=0,ddt(Aθ′)=Aφ′2sin⁡θcos⁡θ−C(ψ′+φ′cos⁡θ)φ′sin⁡θ+mgdsin⁡θ, ddt(C(ψ′+φ′cos⁡θ))=0 hold on I, all functions being evaluated at t. These displayed equations are the definition of the system; they are written out so that no separate definition of a variational principle is needed. The first and third are the statements that ∂φ′L and ∂ψ′L are constant along the trajectory.

Remarks

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Reduction of the heavy symmetric top system: conserved momenta, one-dimensional θ equation and conserved reduced energy

Statement

Fix the parameters A>0, C>0, m>0, d>0, g>0 and let q=(φ,θ,ψ):I→R×(0,π)×R be a heavy-top trajectory of The Euler-angle parametrization and the heavy symmetric top system on an interval I with at least two points; write q′=(φ′,θ′,ψ′) for its derivative. Then:

  1. the functions t↦b(t):=C(ψ′(t)+φ′(t)cos⁡θ(t)) and t↦a(t):=Aφ′(t)sin⁡2θ(t)+b(t)cos⁡θ(t) are constant on I; write b,a for their values;
  2. consequently φ′=a−bcos⁡θAsin⁡2θ and ψ′=bC−(a−bcos⁡θ)cos⁡θAsin⁡2θ on I;
  3. the reduced energy E(t):=12Aθ′(t)2+(a−bcos⁡θ(t))22Asin⁡2θ(t)+b22C+mgdcos⁡θ(t) is constant on I;
  4. with the effective potential U(θ):=(a−bcos⁡θ)22Asin⁡2θ+b22C+mgdcos⁡θ on (0,π), one has Aθ′′(t)=−U′(θ(t))(t∈I), so that the angle θ alone obeys the one-dimensional second-order equation Aθ′′=−U′(θ) with the constants a,b determined by the initial data.

Facts & Assumptions

Given: Parameters A,C,m,d,g>0, a heavy-top trajectory q=(φ,θ,ψ) on I with q′=(φ′,θ′,ψ′), and the partial derivatives of L recorded in The Euler-angle parametrization and the heavy symmetric top system.

[L1]

The system defining a heavy-top trajectory is ddt∂φ′L=0, ddt∂θ′L=∂θL, ddt∂ψ′L=0, with ∂φ′L=Aφ′sin⁡2θ+Ccos⁡θ(ψ′+φ′cos⁡θ), ∂θ′L=Aθ′, ∂ψ′L=C(ψ′+φ′cos⁡θ) and ∂θL=Aφ′2sin⁡θcos⁡θ−C(ψ′+φ′cos⁡θ)φ′sin⁡θ+mgdsin⁡θ The Euler-angle parametrization and the heavy symmetric top system.

[L2]

sin⁡θ>0 for θ∈(0,π), and sin⁡,cos⁡ are differentiable with (sin⁡)′=cos⁡, (cos⁡)′=−sin⁡; (sin⁡θ(t))′=cos⁡θ(t) θ′(t) and (cos⁡θ(t))′=−sin⁡θ(t) θ′(t) Sine and cosine defined by their real power series, The derivatives of sine and cosine are cosine and minus sine, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c).

Proof

technique · direct
1.1givenL1

By the first and third equations of [L1], the functions ∂φ′L=Aφ′sin⁡2θ+Ccos⁡θ(ψ′+φ′cos⁡θ) and ∂ψ′L=C(ψ′+φ′cos⁡θ) are constant on I; the second of these is b, and the first equals Aφ′sin⁡2θ+bcos⁡θ=a. Hence a and b are constants by A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant applied to their continuous differentiable restrictions to I, which is clause 1.

1.2givenL1L3L4algebra

Define E(t):=φ′∂φ′L+θ′∂θ′L+ψ′∂ψ′L−L along the trajectory. By [L4] the derivative of t↦L(q(t),q′(t)) is ∑i∂qiL qi′+∑i∂qi′L qi′′, and substituting the three system equations ddt∂qi′L=∂qiL of [L1] turns this into ddt(∑iqi′∂qi′L); hence E′=0. Since E is continuous on the interval I and differentiable at its interior points, A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant makes it constant on I.

1.3givenL2L3algebra

Compute U′: writing N=(a−bcos⁡θ)2 and D=2Asin⁡2θ, the quotient rule [L3] and [L2] give ddθND=N′D−ND′D2=4Absin⁡3θ(a−bcos⁡θ)−4Asin⁡θcos⁡θ(a−bcos⁡θ)24A2sin⁡4θ, which after factoring (a−bcos⁡θ)/(Asin⁡3θ) is a−bcos⁡θAsin⁡3θ(bsin⁡2θ−cos⁡θ(a−bcos⁡θ)); and the remaining terms contribute −mgdsin⁡θ, so U′(θ)=a−bcos⁡θAsin⁡3θ(bsin⁡2θ−cos⁡θ(a−bcos⁡θ))−mgdsin⁡θ.

2.1step 1.1givenL2L3algebra

Solving the two defining identities: b=C(ψ′+φ′cos⁡θ) gives ψ′=b/C−φ′cos⁡θ, and a=Aφ′sin⁡2θ+bcos⁡θ gives φ′=(a−bcos⁡θ)/(Asin⁡2θ) because sin⁡θ>0 by [L2]; substituting the second into the first gives the displayed formula for ψ′. This is clause 2.

2.2step 1.12.11.2givenL1L3algebra

Evaluating E with the partial derivatives of [L1] and L of the definition gives E=12Aθ′2+12Aφ′2sin⁡2θ+12C(ψ′+φ′cos⁡θ)2+mgdcos⁡θ, the cross terms cancelling after collecting φ′(Aφ′sin⁡2θ+Ccos⁡θ(ψ′+φ′cos⁡θ))+ψ′C(ψ′+φ′cos⁡θ)−12C(ψ′+φ′cos⁡θ)2=12C(ψ′+φ′cos⁡θ)2 and the analogous axial term; substituting φ′=(a−bcos⁡θ)/(Asin⁡2θ) and ψ′+φ′cos⁡θ=b/C from steps 1.1 and 2.1 gives the displayed formula for E. Since E is constant by step 1.2, this is clause 3.

3.1step 1.12.11.3givenL1L3algebra∎

Substituting a−bcos⁡θ=Aφ′sin⁡2θ from step 2.1 into the first term of step 1.3 gives U′(θ)=φ′sin⁡θ(b−Aφ′cos⁡θ)−mgdsin⁡θ, hence −U′(θ)=Aφ′2sin⁡θcos⁡θ−bφ′sin⁡θ+mgdsin⁡θ; by step 1.1 and the identity ∂θL=Aφ′2sin⁡θcos⁡θ−C(ψ′+φ′cos⁡θ)φ′sin⁡θ+mgdsin⁡θ of [L1] with C(ψ′+φ′cos⁡θ)=b, the right-hand side equals ∂θL. Since the second system equation reads ddt(Aθ′)=∂θL and ddt(Aθ′)=Aθ′′, we obtain Aθ′′=−U′(θ) on I, which is clause 4.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Angular velocity in Euler angles: ω3=ψ′+φ′cos⁡θ and ω12+ω22=θ′2+φ′2sin⁡2θ

Statement

Let J⊆R be an interval with at least two points, let φ,θ,ψ:J→R be differentiable with θ(J)⊆(0,π), and put R(t):=Θ(φ(t),θ(t),ψ(t))=Rz(φ(t)) Ry(θ(t)) Rz(ψ(t)), the Euler-angle parametrization of The Euler-angle parametrization and the heavy symmetric top system. Then R is a differentiable SO(3)-valued curve, and its body angular velocity ω:J→R3 in the sense of Rigid placements, body and spatial angular velocity, and material velocities (the unique vector with R(t)TR′(t)=[ω(t)]×) is, in the standard basis, ω=(−φ′sin⁡θcos⁡ψ+θ′sin⁡ψ)e1+(φ′sin⁡θsin⁡ψ+θ′cos⁡ψ)e2+(ψ′+φ′cos⁡θ)e3, all functions evaluated at t. Consequently ω3=ψ′+φ′cos⁡θ,ω12+ω22=θ′2+φ′2sin⁡2θ,∥ω∥22=θ′2+φ′2sin⁡2θ+(ψ′+φ′cos⁡θ)2. Moreover, if A,C>0 and I:R3→R3 is the linear map with Ie1=Ae1, Ie2=Ae2, Ie3=Ce3, then ωTIω=A(θ′2+φ′2sin⁡2θ)+C(ψ′+φ′cos⁡θ)2.

Facts & Assumptions

Given: Differentiable φ,θ,ψ:J→R with θ∈(0,π), the curve R=Rz(φ)Ry(θ)Rz(ψ), and the matrices Rz,Ry of The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives.

[L1]

Rz(a),Ry(a)∈SO(3) for all a, Rz(a)Rz(b)=Rz(a+b), Ry(a)Ry(b)=Ry(a+b), and for differentiable a one has (Rz∘a)′=a′Rz(a)[e3]×=a′[e3]×Rz(a) and (Ry∘a)′=a′Ry(a)[e2]×=a′[e2]×Ry(a) The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives.

[L2]

The product rule for matrix-valued curves, the skewness of RTR′ for an SO(3)-valued curve, and the equivalence RTR′=[ω]×⇔R′=R[ω]× hold; moreover [ω]×z=ω×z and RT[u]×R=[RTu]× Derivative rules for differentiable matrix-valued curves, and skewness of RTR′ for a curve in SO(3), The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×, Rigid placements, body and spatial angular velocity, and material velocities.

[L3]

For R∈SO(3), RTR=I3; the matrices e1e1T+e2e2T+e3e3T=I3, [e3]×e3=0 and [e2]×e2=0 The skew-matrix/cross-product correspondence [ω]×z=ω×z and its conjugation law R[ω]×RT=[Rω]×, Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes.

[L5]

Evaluation of Rz(a)e3=e3, Rz(a)e1=cos⁡a e1+sin⁡a e2, Rz(a)e2=−sin⁡a e1+cos⁡a e2 and Ry(a)e1=cos⁡a e1−sin⁡a e3, Ry(a)e3=sin⁡a e1+cos⁡a e3 is the column reading of the displayed matrices The coordinate rotation matrices Rz(a) and Ry(a): values in SO(3), composition law, and derivatives.

Proof

technique · direct
1.1givenL1L2

By [L1] each of the three factors is a differentiable SO(3)-valued curve; by the product rule [L2] and closure of SO(3) under multiplication, R is differentiable with values in SO(3). Hence RTR′ is skew and there is a unique ω with RTR′=[ω]× by [L2].

1.2givenL2L5algebra

By the conjugation identity ST[u]×S=[STu]× of [L2] with S=Ry(θ) and S=Rz(ψ), and the column readings of [L5], Ry(θ)T[e3]×Ry(θ)=[Ry(θ)Te3]×=[−sin⁡θ e1+cos⁡θ e3]×,Rz(ψ)T[e2]×Rz(ψ)=[sin⁡ψ e1+cos⁡ψ e2]×.

2.1step 1.1givenL1L2algebra

Differentiating the product and applying [L1], R˙=φ′ Rz(φ)[e3]×Ry(θ)Rz(ψ)+θ′ Rz(φ)Ry(θ)[e2]×Rz(ψ)+ψ′ Rz(φ)Ry(θ)Rz(ψ)[e3]×.

3.1step 2.1givenL2L3algebra

Multiplying on the left by RT=Rz(ψ)TRy(θ)TRz(φ)T and using Rz(φ)TRz(φ)=I3 [L3] term by term, RTR˙=φ′ Rz(ψ)TRy(θ)T[e3]×Ry(θ)Rz(ψ)+θ′ Rz(ψ)T[e2]×Rz(ψ)+ψ′ [e3]×.

4.1step 3.11.2givenL2L5algebra

Conjugating the first vector of step 1.2 by Rz(ψ)T and using [L5], Rz(ψ)T(−sin⁡θ e1+cos⁡θ e3)=−sin⁡θcos⁡ψ e1+sin⁡θsin⁡ψ e2+cos⁡θ e3; substituting into step 3.1 gives RTR˙=[ω]× with ω exactly as displayed in the statement.

5.1step 4.1givenL4algebra

Reading off the components gives ω3=ψ′+φ′cos⁡θ; expanding the squares of the first two components and using sin⁡2ψ+cos⁡2ψ=1 from [L4] cancels the mixed terms and gives ω12+ω22=φ′2sin⁡2θ+θ′2, whence ∥ω∥22=θ′2+φ′2sin⁡2θ+(ψ′+φ′cos⁡θ)2.

6.1step 4.15.1givenL3algebra∎

For the diagonal map I with Ie1=Ae1, Ie2=Ae2, Ie3=Ce3, the quadratic form is ωTIω=A(ω12+ω22)+Cω32, and substituting step 5.1 gives the displayed identity.

Remarks

5 · Examples, experiments and counterexamples

None yet.

Sources