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
- Absolute and Conditional Convergence; Rearrangement; Products
- Algebraic Closure, Embeddings, and Separability
- Algebraic Extensions, Extension Degree, and Finite Fields
- Areas of Elementary Plane Figures
- Binary Operations, Monoids, Groups and Subgroups
- Cm Math Balances
- Cm Math Geometry
- Compactness in Metric Spaces
- Completeness, Completion, and Uniform Continuity
- Composition Series, the Jordan–Hölder Theorem and Solvable Groups
- Congruences, the Integers Modulo n and the Chinese Remainder Theorem
- Connectedness
- Construction of the Natural Numbers
- Construction of the Real Numbers via Cauchy Sequences
- Construction of the Real Numbers via Dedekind Cuts
- Continuity, IVT, EVT, and Uniform Continuity
- Cosets, Index and Lagrange's Theorem
- Countability and Uncountability
- Cyclic Groups and Direct Products
- Determinants of Matrices over a Commutative Ring
- Diagonalisation and the Minimal Polynomial
- Divisibility, Euclidean Domains, Principal Ideal Domains and Unique Factorisation
- Divisibility, Greatest Common Divisors and Bézout's Identity
- Dual Spaces, Bilinear and Quadratic Forms, and Sylvester's Law of Inertia
- Eigenvalues, Eigenvectors and the Characteristic Polynomial
- Exterior Powers, Orientation and Hodge Duality
- Filters and Ultrafilters
- Finite Counting, Factorials and Binomial Coefficients
- Foundations of the Real Numbers for Analysis
- Gaussian Elimination, Elementary Matrices and Reduced Row Echelon Form
- Group Actions, Orbits, Stabilisers and Cayley's Theorem
- Group Homomorphisms and the Isomorphism Theorems
- Ideals, Quotient Rings and the Isomorphism Theorems for Rings
- Inner Product Spaces, Gram-Schmidt, Projections and Adjoints
- Lebesgue Measure on Euclidean Space
- Limits of Real Functions
- limsup, liminf, and Subsequential Limits
- Linear Independence, Bases and Dimension
- Linear Transformations, Rank-Nullity and Quotient Spaces
- Matrices, the Matrix of a Linear Map, and Change of Basis
- Measurable Functions and Simple Approximation
- Measures and Their Basic Properties
- Metric Spaces
- Mixed Partials, Taylor Formulae, and Extrema
- Monotone Functions, Discontinuities, and Continuity Sets
- Monotone Sequences, Bolzano-Weierstrass, and Cauchy Completeness
- Normal Subgroups and Quotient Groups
- Normed and Banach Spaces
- Order, Zorn's Lemma, and the Axiom of Choice
- Outer Measure and the Caratheodory Extension Theorem
- Picard-Lindelöf and First-Order Ordinary Differential Equations
- Polynomial Rings, the Division Algorithm and Roots
- Power Series and Real-Analytic Functions
- Primes, Euclid's Lemma and the Fundamental Theorem of Arithmetic
- Properties of the Integral and the Working FTC
- Regular Surfaces and Surface Integrals
- Relations, Functions, and Quotients
- Rings, Subrings, Integral Domains and Fields
- Rⁿ as a Normed Space; Vector-Valued Functions
- Roots, Rational Powers, and Classical Inequalities
- Sequences and Limits
- Sequences and Series of Functions; Uniform Convergence
- Series: Convergence and the Nonnegative Tests
- Sigma Algebras and Borel Sets
- Simple Field Extensions and the Construction of the Complex Numbers
- Sine, Cosine, and the Definition of Pi
- Splitting Fields
- Subspaces, Products, and Quotients
- Suprema and Infima
- Sylow's Theorems, p-Groups and Nilpotent Groups
- Symmetric Groups, Cycle Decomposition and the Sign Homomorphism
- The Derivative and the Mean Value Theorems
- The Determinant of a Linear Operator, Cofactors and Cramer's Rule
- The Exponential Function
- The Fundamental Theorem of Algebra
- The Fundamental Theorem of Finite Abelian Groups
- The Galois Correspondence
- The Inverse and Implicit Function Theorems
- The Inverse Function Theorem Completed
- The Lebesgue Integral and the Convergence Theorems
- The Logarithm and General Powers
- The Riemann Integral in Rᵐ and Jordan Content
- The Riemann Integral: Definition and Integrability
- The Spectral Theorem, Positive Operators and Singular Value Decomposition
- The Total Derivative in ℝᵐ → ℝⁿ
- The ZFC Axioms and the Basic Set Constructions
- Topological Spaces and Continuity
- Topology of ℝ
- Triangularisation, Generalised Eigenspaces and Jordan Canonical Form
- Vector Spaces, Linear Subspaces, Span and Direct Sums
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 as functions on . Write for its standard basis from The standard list with and for is an ordered basis of ; hence , and is the zero space with basis and dimension , and, in the conventional one-based notation of the rigid-body formulas, set , , and for . For a matrix , write for ; sums indexed from to use these labels. In particular, define . Whenever a curve 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 th component is . Twice differentiable means each component has a first derivative that is differentiable. The Borel sigma-algebra on below is for its usual Euclidean metric topology: is the metric of as the set of functions , and , , 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 for a curve in proves the sum, transpose and product rules for matrix-valued curves and the skew-symmetry of for a curve in ; The skew-matrix/cross-product correspondence and its conjugation law establishes the bijection between skew matrices and vectors and its conjugation law ; Rigid placements, body and spatial angular velocity, and material velocities then defines a rigid placement , its body and spatial angular velocities, and its material velocities.
Inertia. Mass distributions on : total mass, first moment and inertia matrix defines a mass distribution as a bounded finite Borel measure on , with total mass, first moment and inertia matrix ; Finiteness of the moments of a bounded finite Borel measure on proves every displayed integral finite; Symmetry, the quadratic identity , and the positive-definiteness criterion for the inertia matrix proves symmetry, the identity and the criterion that 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: and computes and for a centred distribution, by integrating the pointwise identities of the placement. Euler's equations under an adopted torque balance, and the torque-free invariants and principal-axis equations deduces Euler's equations from an explicitly adopted torque balance, and in the torque-free case the equation with the conserved quantities and 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 and 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 , 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 and : values in , composition law, and derivatives collects the coordinate rotation matrices , their composition law and their derivatives; The Euler-angle parametrization and the heavy symmetric top system defines the Euler-angle parametrization and the heavy-top system by its three displayed second-order equations; Angular velocity in Euler angles: and computes the body angular velocity in Euler angles, including and ; 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 .
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
Derivative rules for differentiable matrix-valued curves, and skewness of for a curve in
Statement
Let be an interval with at least two points, let and let be differentiable at , meaning (The vector space of by matrices over a field, with entrywise operations) that every entry function and is differentiable at in the sense of The derivative of at a point that is a limit point of , and differentiability on a set; write for the matrix with entries , and similarly for .
- and are differentiable at , with and .
- is differentiable at , with .
- If for every (The orthogonal group and the rotation group ), then the matrix is skew-symmetric, i.e. at , and likewise is skew-symmetric. In particular the trace of each of these two matrices vanishes.
Facts & Assumptions
Given: An interval with at least two points, a point , and maps differentiable at in the entrywise sense of the statement.
If are differentiable at and , then , and are differentiable at with , and Sums, scalar multiples, products and quotients: , , , and when .
is the function space of entries indexed by with pointwise sum, so ; the transpose is ; and the product is The vector space of by matrices over a field, with entrywise operations, The transpose of a matrix, Rectangular matrix multiplication and the identity matrix , including zero-sized shapes.
Transposition is involutive and reverses products: and Transpose is linear and involutive, and .
For one has , , invertible with , and also The orthogonal group and the rotation group , is a subgroup of the general linear group, with .
A constant function is differentiable at every point of with derivative : its difference quotient at any limit point of is the constant function , whose limit is The derivative of at a point that is a limit point of , and differentiability on a set.
The identity matrix has entries for equal indices and otherwise, and the matrix all of whose entries are is the zero matrix Rectangular matrix multiplication and the identity matrix , including zero-sized shapes, The vector space of by matrices over a field, with entrywise operations.
Proof
For all the entry is differentiable at with derivative by [L1] and [L2]; these are the entries of , so is differentiable with the asserted derivative.
For all the entry is differentiable at with derivative , which is the entry of ; so is differentiable with the asserted derivative.
For all the entry is a finite sum of products of functions differentiable at , hence differentiable at with derivative by two applications of [L1]; this sum is the entry of by [L2]. Hence is differentiable at with .
Assume for all . By steps 1.2 and 1.3 the map is differentiable at , with derivative at . By [L4] this map is the constant map , and every entry of this constant entry function is constant, so by [L5] and [L2] its derivative at is the zero matrix [L6]. Hence at .
By [L4] and [L3] the map is also the constant map , so step 1.3 gives at , and transposing with [L3] gives .
Transposing the identity of step 2.1 and using [L3] gives , hence , so is skew-symmetric.
A skew-symmetric matrix has zero diagonal, hence trace : the entry satisfies by step 3.1 applied entrywise, so in ; this applies to the two matrices of steps 3.1 and 2.2.
The skew-matrix/cross-product correspondence and its conjugation law
Statement
Use the standard basis of and coordinates (The standard list with and for is an ordered basis of ; hence , and is the zero space with basis and dimension , A finite list is an ordered basis if and only if every equals for exactly one ; those scalars are the coordinates of in that ordered basis). Define the map by Then:
- for all , the cross product being that of The cross product in .
- for every ; the map is linear and injective, and its image is exactly the set of skew-symmetric matrices of . Consequently every skew-symmetric equals for exactly one , namely .
- For every (The orthogonal group and the rotation group ) and every ,
Facts & Assumptions
Given: Vectors , a skew-symmetric matrix and a rotation .
The cross product has coordinates , , in the standard basis, and satisfies , where is the matrix with these three vectors as columns; moreover is orthogonal to and to The cross product in , The cross product is bilinear, alternating, and orthogonal to both factors.
Use the owning page's one-based labels: , , and is the standard unit vector at the actual index , for . Relabelling the supplier's product formula gives . Sums are entrywise and ; the transpose is involutive and reverses products The vector space of by matrices over a field, with entrywise operations, Rectangular matrix multiplication and the identity matrix , including zero-sized shapes, The transpose of a matrix, Transpose is linear and involutive, and , The standard list with and for is an ordered basis of ; hence , and is the zero space with basis and dimension .
The determinant is multiplicative, , and For same-sized finite square matrices over a commutative ring, , Rectangular matrix multiplication and the identity matrix , including zero-sized shapes.
For one has , , and invertible; for all therefore , since The orthogonal group and the rotation group , is a subgroup of the general linear group, with , The Euclidean inner product on .
The Euclidean inner product is positive definite: if for all then Real and complex inner product spaces, with the inner product linear in the first argument, The Euclidean inner product on .
A matrix is determined by its values on the standard basis: if satisfy for all , then , since the -th column of is by the product formula [L2] and is a basis The standard list with and for is an ordered basis of ; hence , and is the zero space with basis and dimension , A finite list is an ordered basis if and only if every equals for exactly one ; those scalars are the coordinates of in that ordered basis.
Proof
Comparing coordinates, , and likewise and , all by [L1] and [L2]; this is clause 1.
The displayed matrix is manifestly antisymmetric: for every by [L2]. Hence . Entries are linear functions of , so is linear. If , then the entry gives , the entry gives and the entry gives , so : the map is injective.
Let . By [L4], by [L1]; by [L2] and [L3], by [L1] again. Since is invertible, ranges over all of as does, so [L5] gives .
Let be skew-symmetric and put . Comparing entries with the definition of and using and : has entry , entry , entry , 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.
For every , by step 1.1, by step 1.3, by [L4], by step 1.1. Hence by [L6]; this is clause 3.
Mass distributions on : total mass, first moment and inertia matrix
Definition
Let be the Borel sigma-algebra of (The Borel sigma-algebra of a topological space) and let be a measure on the measurable space (Measures on sigma-algebras, Measure spaces) that is finite, , and bounded, meaning that there is a real with 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: and , whenever every displayed scalar integral is finite.
A mass distribution on is a bounded finite Borel measure with ; the number is its total mass. Writing points of as and using the Euclidean inner product and norm of The Euclidean inner product on , define:
- the first moment , read componentwise as ;
- the inertia matrix (The vector space of by matrices over a field, with entrywise operations) by its entries , where if and otherwise; equivalently with the identity matrix and the matrix product of Rectangular matrix multiplication and the identity matrix , including zero-sized shapes.
The mass distribution is centred when . Under the Axiom of Countable Choice (The Axiom of Countable Choice ()), each of the integrals above is a finite real number: on one has and , and ; the precise statement and proof are Finiteness of the moments of a bounded finite Borel measure on ↗, 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 is not asserted here.
Finiteness of the moments of a bounded finite Borel measure on
Statement
Assume the Axiom of Countable Choice (The Axiom of Countable Choice ()). Let be a measure on (Mass distributions on : total mass, first moment and inertia matrix) that is finite and bounded: and for some (Open ball, closed ball and sphere in a metric space). Then:
- for all the functions and are Borel measurable and -integrable, with and ;
- consequently and the scalar integrals and are finite, with and . When , these integrals are precisely the components of the first moment and inertia matrix of Mass distributions on : total mass, first moment and inertia matrix. When , all these integrals vanish; the zero measure is not called a mass distribution.
Facts & Assumptions
Given: A finite Borel measure on and a real with .
Every continuous map is Borel measurable, under Continuous functions on Euclidean spaces are Borel measurable, Borel measurable and Lebesgue measurable functions on ; and sums, scalar multiples and pointwise products of measurable real-valued functions are measurable Arithmetic and lattice operations preserve measurability whenever they are defined.
For a nonnegative simple function the nonnegative Lebesgue integral equals the simple integral; in particular the constant function on has The nonnegative Lebesgue integral, The integral of a nonnegative simple function.
The nonnegative integral is monotone and homogeneous for nonnegative scalars: if then , and for Monotonicity and nonnegative homogeneity of the nonnegative integral.
A real measurable function with is integrable, and then Integrable real and complex functions, and their integrals, The modulus of an integral is bounded by the integral of the modulus.
satisfies for each and exactly on (The Euclidean inner product on , Open ball, closed ball and sphere in a metric space).
The closed ball is closed in the Euclidean metric topology, so is open and belongs to ; thus the given value is defined ( as the set of functions , and , , are metrics on it, Arbitrary unions and finite intersections of open sets are open, open balls are open and closed balls are closed, The metric topology: a set is open when every one of its points has a ball around it inside the set; closed means open complement, The Borel sigma-algebra of a topological space).
Proof
Each coordinate function is continuous, hence Borel measurable by [L1]; the products are then measurable by the product clause of [L1].
Let , a Borel set by [L6], with . We justify the almost-everywhere comparison from the definition of the nonnegative integral: if a nonnegative measurable satisfies on for a finite constant , every nonnegative simple minorant has the same simple integral as . Indeed, in a finite disjoint level-set representation of , the contribution of each intersection with is zero, while for every measurable level set . Pointwise , so [L2] and [L3] give ; taking the supremum over simple minorants gives by The nonnegative Lebesgue integral. Apply this with , , and with , : measurability follows from step 1.1 and [L1], and [L5] gives the bounds on . Thus and .
The bounds of step 2.1 are finite because is finite, so each of the functions and has finite integral of its modulus, hence is -integrable by [L4].
Each coordinate integrand is integrable by step 3.1. Each inertia integrand is measurable by [L1], and on its modulus is at most by [L5]. The simple-minorant comparison proved in step 2.1 gives , so [L4] makes integrable and gives . Likewise [L4] and step 2.1 give . If , is a mass distribution and these are its moment and inertia components by Mass distributions on : total mass, first moment and inertia matrix. If , 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.
Rigid placements, body and spatial angular velocity, and material velocities
Definition
Let be an interval with at least two points (Intervals of : the nine order-convex forms, nondegeneracy, and length), let be a mass distribution on with total mass (Mass distributions on : total mass, first moment and inertia matrix), let be differentiable and let be differentiable (The orthogonal group and the rotation group ), differentiability meaning entrywise differentiability of in the sense of The derivative of at a point that is a limit point of , and differentiability on a set with derivative ; write for the derivative of .
The rigid placement with reference shape , translation and rotation is the map so that the material point labelled by (a point of the reference shape) is carried to at time ; the map is the placement at time , and the pair is its configuration. The material velocity of the label at time is the derivative which exists by differentiability of and of the entries of .
Body and spatial angular velocity. By Derivative rules for differentiable matrix-valued curves, and skewness of for a curve in the matrix is skew-symmetric at every . By The skew-matrix/cross-product correspondence and its conjugation law there is therefore a unique vector with and the map is the body angular velocity of the placement; the spatial angular velocity is . Explicitly , a finite sum of products of entries of and ; consequently, if 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 and its conjugation law , .
Remarks
-
The reference shape is a label space. is not the position of the material point at time ; the position is . 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 is the differentiated orthogonality relation of the curve ; 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 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 , and not that is measured by any apparatus. The physical reading is supplied by the consuming classical-mechanics items.
Symmetry, the quadratic identity , and the positive-definiteness criterion for the inertia matrix
Statement
Assume the Axiom of Countable Choice (The Axiom of Countable Choice ()), used only in [L1], through Finiteness of the moments of a bounded finite Borel measure on , to make the entries of finite real numbers. Let be a mass distribution on with total mass and inertia matrix (Mass distributions on : total mass, first moment and inertia matrix). Write for the quadratic form of on and for the Euclidean norm (The Euclidean inner product on ). Then:
- is symmetric: for all .
- For every , the cross product being that of The cross product in .
- The symmetric bilinear form on is therefore positive semidefinite (Positive and negative definiteness, the inertia , rank , and signature of a real symmetric bilinear or quadratic form); for each , writing , one has if and only if . Thus it is positive definite if and only if for every , that is, no line through the origin carries the entire mass (Measure-null sets and almost-everywhere statements relative to a measure).
- If the form is positive definite, then is invertible.
Facts & Assumptions
Given: A mass distribution with inertia matrix and a vector .
is a finite real number for every , by the finiteness lemma for bounded finite measures Mass distributions on : total mass, first moment and inertia matrix, Finiteness of the moments of a bounded finite Borel measure on .
The integral is linear on : for real and integrable The Lebesgue integral is linear on .
For a measurable with one has -almost everywhere; and for A nonnegative measurable function has integral exactly when it vanishes almost everywhere, Monotonicity and nonnegative homogeneity of the nonnegative integral.
For all , , and if and only if and are linearly dependent The squared cross-product norm is the Gram determinant of two vectors.
A vector is orthogonal to every vector of only if it is , and the Euclidean inner product is bilinear Real and complex inner product spaces, with the inner product linear in the first argument, The Euclidean inner product on .
The matrix satisfies ; if then is invertible Invertible matrix theorem: invertibility, full pivot rank, RREF , trivial nullspace and unique solvability are equivalent.
Proof
Interchanging and leaves the integrand unchanged, because and ; the two integrals defining and therefore coincide, so is symmetric.
Expressing the quadratic form through the entries of [L1] and applying linearity [L2] to the finite sum of integrable functions ,
For each , evaluating the sum gives by [L4] and bilinearity of the inner product.
Substituting step 2.1 into step 1.2 gives the identity of clause 2, and the integrand is nonnegative, so its integral is by [L3]: the form is positive semidefinite by definition.
Let and put . This is a measurable nonnegative integrable polynomial by [L1], [L2] and step 2.1. Its zero set is the closed line , by [L4], so its exceptional set is Borel measurable. By step 3.1 and both implications of [L3], exactly when -almost everywhere, meaning that its measurable exceptional set has measure zero, or equivalently . Since the quadratic form is nonnegative and finite, it is strictly positive for this exactly when that complement has positive measure. Quantifying over all gives the positive-definiteness criterion in clause 3.
Assume the form is positive definite. If for some , then , which for is excluded; hence and is invertible by [L6]. This is clause 4.
Kinetic energy and angular momentum of a rigid placement: and
Statement
Assume the Axiom of Countable Choice (The Axiom of Countable Choice ()), used only through the bound of Finiteness of the moments of a bounded finite Borel measure on that makes the two integrands of clause 1 -integrable. Let be an interval with at least two points, let be a centred mass distribution on with total mass and inertia matrix (Mass distributions on : total mass, first moment and inertia matrix), let and be differentiable, and let be the rigid placement with body angular velocity (Rigid placements, body and spatial angular velocity, and material velocities). For define the second integral read componentwise (Mass distributions on : total mass, first moment and inertia matrix). Then for every :
- for every , and ; moreover the two integrands displayed above are -integrable, so and are a well-defined finite real number and vector;
- ;
- .
Facts & Assumptions
Given: A centred mass distribution , differentiable and , the rigid placement , its body angular velocity with , and a time .
The body angular velocity satisfies , hence , and by conjugation as well; explicitly Rigid placements, body and spatial angular velocity, and material velocities, The skew-matrix/cross-product correspondence and its conjugation law .
The elements of preserve the Euclidean inner product and norm: , so The orthogonal group and the rotation group , The Euclidean inner product on , Real and complex inner product spaces, with the inner product linear in the first argument.
The measure has finite total mass, for some , and vanishing first moment for ; 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 (Mass distributions on : total mass, first moment and inertia matrix).
The integral is linear on and equals the componentwise integral for vector- or matrix-valued integrands The Lebesgue integral is linear on , Integral over a measurable subset, Integrable real and complex functions, and their integrals.
The inertia matrix satisfies Symmetry, the quadratic identity , and the positive-definiteness criterion for the inertia matrix.
For , the cross product is bilinear, is determined componentwise by and its cyclic analogues, and The cross product in , The squared cross-product norm is the Gram determinant of two vectors.
Proof
By [L1], , and by [L2]; hence for every the derivative of at is .
For every one has the vector identity : for instance the first coordinate of the left side is , which equals the first coordinate of the right side, and the second and third coordinates are the cyclic analogues.
The first integrand satisfies by step 1.1, [L3] and the two elementary bounds and ; the second integrand satisfies by step 1.1 and [L3]. Both are bounded on the support ball of by constants depending on , so they are -integrable by [L4]; this is the integrability clause of the statement.
Expanding the square with bilinearity and symmetry of the inner product and using [L3], .
By step 1.1 the relative velocity is , so the angular-momentum integrand is by [L2].
Integrating step 2.2 with [L5] and pulling the constant and the linear map through the integral, The middle term vanishes because by [L4], and the last term equals by [L6]; substituting into the definition of gives clause 2.
Integrating step 2.3 and substituting step 1.2, , the constant matrix and the finite sum of integrals being interchanged by [L5] componentwise, and the last equality being the definition of in Mass distributions on : total mass, first moment and inertia matrix. This is clause 3.
Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case
Statement
Assume the Axiom of Countable Choice (The Axiom of Countable Choice ()), used through Symmetry, the quadratic identity , and the positive-definiteness criterion for the inertia matrix to ensure finiteness of the inertia integrals. Let be a mass distribution on with inertia matrix (Mass distributions on : total mass, first moment and inertia matrix), and write for the matrix-vector product and for the Euclidean inner product on (Rectangular matrix multiplication and the identity matrix , including zero-sized shapes, The Euclidean inner product on ). Then:
- there is an orthonormal basis of (Orthogonal vectors and subspaces, orthogonal and orthonormal sets, and orthonormal bases) and real numbers with for ;
- for any such basis and any one has , , and ;
- each vanishes exactly when , that is, the whole mass lies on the line ; if the inertia form is positive definite then every and is invertible with ; and if is not supported on any line through the origin, meaning that for every , then every .
Facts & Assumptions
Given: A mass distribution with inertia matrix , an orthonormal basis of eigenvectors of with eigenvalues , and a vector .
is symmetric and ; for each , its quadratic form vanishes exactly when , and the associated bilinear form is positive definite exactly when for every , and then for Symmetry, the quadratic identity , and the positive-definiteness criterion for the inertia matrix.
The map is self-adjoint with respect to the Euclidean inner product, because by symmetry of and the product and transpose rules Symmetry, the quadratic identity , and the positive-definiteness criterion for the inertia matrix, The Euclidean inner product on .
Every self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal basis of eigenvectors Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis, Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space; in an orthonormal basis, self-adjointness is symmetry of the matrix In an orthonormal basis, self-adjoint means conjugate-transpose symmetry and normal means commuting with the conjugate transpose.
A linear map is determined by its values on a basis, and for the eigenbasis A finite list is an ordered basis if and only if every equals for exactly one ; those scalars are the coordinates of in that ordered basis.
A square matrix is invertible exactly when its nullspace is , and the inverse is unique Invertible matrix theorem: invertibility, full pivot rank, RREF , trivial nullspace and unique solvability are equivalent, Invertible matrices and the general linear group , Rectangular matrix multiplication and the identity matrix , including zero-sized shapes.
Proof
The map is linear and, by [L1] and [L2], self-adjoint on the Euclidean space .
Applying [L3] to produces an orthonormal basis and reals with ; this is clause 1.
For , linearity of [L4] gives and hence by orthonormality; setting gives by [L1], and the displayed chain of equalities of clause 2 follows.
Each is nonnegative, and it is zero exactly when , by [L1]; this is the vanishing criterion of clause 3.
If the inertia form is positive definite, then for every ; moreover for every by [L1], so the nullspace of is and is invertible with unique inverse by [L5]. Define on the basis by and extend linearly; then and for every , so by [L4], whence by uniqueness and .
If is not supported on any line through the origin, then for , so by step 4.1. This completes clause 3.
Euler's equations under an adopted torque balance, and the torque-free invariants and principal-axis equations
Statement
Assume the Axiom of Countable Choice (The Axiom of Countable Choice ()), used through the inertia suppliers to ensure finiteness of the integrals. Let be an interval with at least two points (Intervals of : the nine order-convex forms, nondegeneracy, and length), let be a centred mass distribution on that is not supported on any line through the origin, with inertia matrix (Mass distributions on : total mass, first moment and inertia matrix), let be differentiable and let be twice differentiable, and let be the rigid placement with body angular velocity (Rigid placements, body and spatial angular velocity, and material velocities). Write ; by Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case the inertia form is positive definite and is invertible.
Adopt the balance hypothesis (B): there is a map such that the curve is differentiable and Then:
- is differentiable and Euler's equations with torque hold:
- if on , then satisfies the torque-free Euler equations , and the functions and are constant on ;
- if on and is a positively oriented orthonormal eigenbasis of with and , then and are constant on , and the torque-free equations read
Facts & Assumptions
Given: The data of the statement, the body angular velocity with , the angular momentum , and a map satisfying the balance hypothesis (B).
The body angular velocity satisfies , so and ; it is differentiable because it is a finite sum of products of entries of and Rigid placements, body and spatial angular velocity, and material velocities, The skew-matrix/cross-product correspondence and its conjugation law .
The product rule for matrix-valued curves gives Derivative rules for differentiable matrix-valued curves, and skewness of for a curve in , and The skew-matrix/cross-product correspondence and its conjugation law .
is symmetric and positive definite, hence invertible; and for the eigenbasis with one has Symmetry, the quadratic identity , and the positive-definiteness criterion for the inertia matrix, Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case.
The cross product is bilinear and alternating, and , so is orthogonal to and to The cross product is bilinear, alternating, and orthogonal to both factors, The cross product in .
The Euclidean inner product is bilinear and symmetric and The Euclidean inner product on , Real and complex inner product spaces, with the inner product linear in the first argument.
The scalar product rule differentiates finite sums of products of differentiable real functions. In an orthonormal eigenbasis of the symmetric positive-definite , has diagonal entries and is therefore symmetric (Sums, scalar multiples, products and quotients: , , , and when , Principal axes and principal moments of the inertia matrix, and positivity in the nondegenerate case).
Proof
By [L1] the angular velocity is differentiable, hence so is ; the product satisfies by [L3], and by [L1] and [L3].
By the balance hypothesis (B) and step 1.1, . Multiplying by and using [L2] and linearity of matrix-vector multiplication gives . Since with a constant matrix, by the finite-sum and scalar rules of [L7], and ; therefore , which is clause 1.
Assume . Then step 2.1 gives ; multiplying by from [L4] gives , hence , using antisymmetry of the cross product [L5] and [L4]. This is the first part of clause 2.
The two functions of clause 2 are constant: the product rule and symmetry of in [L7] give and . Both inner products vanish because a cross product is orthogonal to each of its factors [L5]. Also because and is symmetric [L4], [L6]. These functions are continuous on the solution interval and differentiable at every interior point, so their zero derivatives make them constant by A function continuous on an interval whose derivative vanishes at every interior point of is constant on ; consequently two such functions with the same derivative differ by a constant.
For the eigenbasis, gives by [L4]; by orthonormality and , and these are constant by step 4.1.
Let have columns in the standard basis. The positive orientation and orthonormality give ; [L2] gives from the skew/cross correspondence. Thus the standard cross-product coordinate formula applies in this principal basis. Coordinate-wise, reads , and , 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.
Lyapunov stability of spin about the largest and smallest principal axes of the torque-free top
Statement
Let be positive reals and let be the linear map with on the standard orthonormal basis of . Consider the autonomous first-order system on , read as in Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system with the cross product of The cross product in and the inverse of (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 under an adopted torque balance, and the torque-free invariants and principal-axis equations in angular-momentum variables. Let . Then:
and are equilibria of and each is Lyapunov stable in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system: for and and for every there is such that every maximal solution with is defined for all and satisfies for every .
Facts & Assumptions
Given: Positive reals , the diagonal map with , the field on , and .
The system on the open set is read as First-order systems, initial value problems, and solutions on intervals with ; and are equilibria exactly when , and then Lyapunov stability is the quantified ball statement of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system.
The cross product is bilinear and alternating and satisfies The cross product in , The squared cross-product norm is the Gram determinant of two vectors, The Euclidean inner product on .
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 , bilinearity gives [L2]. Since is diagonal with entries at most in absolute value, ; hence . A compact time-state cylinder is a compact subset of Euclidean , hence bounded by Heine-Borel (Heine-Borel in : with the Euclidean metric a subset of 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 to contain its state projection gives the finite Lipschitz constant required by the definition. The same local estimate gives continuity of 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).
Along every solution of the quantities and are constant; for they equal and The same algebra as in Euler's equations under an adopted torque balance, and the torque-free invariants and principal-axis equations applies directly to the present diagonal map: is symmetric, so and by cross-product orthogonality. Both scalar functions are continuous and differentiable on the solution interval; A function continuous on an interval whose derivative vanishes at every interior point of is constant on ; 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.
Any solution of an autonomous system on that is trapped near a finite endpoint in a compact subset of 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.
The Euclidean inner product is bilinear and positive definite, so and implies Real and complex inner product spaces, with the inner product linear in the first argument, The Euclidean inner product on .
A closed bounded subset of is compact by Heine-Borel in : with the Euclidean metric a subset of 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.
The Euclidean norm satisfies the triangle inequality, so ; hence the norm is continuous and each sphere is closed (Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation).
Proof
By antisymmetry of the cross product, and , so and are equilibria; and [L3] supplies the existence, uniqueness and continuation theory for the system.
Let be the maximal solution with prescribed datum at . By [L4] the numbers and are constant, so the values of lie in the sphere . This sphere is closed by [L8] and bounded, hence compact by Heine-Borel in : with the Euclidean metric a subset of 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.
Define, for , the quantity . Since for by and , one has with .
Every maximal solution of the system is defined for all : if a solution were defined on with and admitted no extension past , then by step 1.2 it would be trapped in the compact subset of near , and [L5] would give an extension past , a contradiction.
Let and let be as in step 1.3; choose satisfying the finite list of conditions , , , and ; such a exists because each condition only bounds from above by a positive number. Let be the maximal solution with , defined for all by step 2.1. Then and, by [L4] and step 1.3 evaluated at , for every
With the data of step 3.1: by and , so and has the sign of (it is continuous, throughout, and ). Hence and
Combining steps 3.1 and 4.1, for all ; since for every such a was produced and every maximal solution starting in is defined for all by step 2.1, is Lyapunov stable in the sense of [L1].
The same argument with the roles of and exchanged uses with , and . Set ; the same positive conditions of step 3.1, with in place of , bound . The norm invariant then gives with the sign of , and the calculation of step 4.1 gives ; the certificate 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 and proves its Lyapunov stability.
Instability of spin about the intermediate principal axis of the torque-free top
Statement
Let be positive reals, let be the linear map with on the standard orthonormal basis, and consider the torque-free Euler system , , on (Euler's equations 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 . Then is an equilibrium of and is not Lyapunov stable: with for every there is a solution with and some with ; Moreover, with , the escape can be arranged at a time where may be chosen as small as desired.
Facts & Assumptions
Given: Positive reals , the diagonal map , the field on , and .
The system is read as First-order systems, initial value problems, and solutions on intervals. For , bilinearity gives ; since is diagonal with entries at most in absolute value and , this yields . This bound makes 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 , 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 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).
Along every solution the invariants and are constant; in the eigenbasis one has , , and the component equations , , The same algebra as in Euler's equations under an adopted torque balance, and the torque-free invariants and principal-axis equations applies directly to the present diagonal map: is symmetric, so and by cross-product orthogonality. Both scalar functions are continuous and differentiable on the solution interval; A function continuous on an interval whose derivative vanishes at every interior point of is constant on ; 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.
The Euler field is even: , by bilinearity and linearity (The cross product in ). Put , so and .
The cross product is bilinear and ; the Euclidean inner product is positive definite The cross product in , The squared cross-product norm is the Gram determinant of two vectors, The Euclidean inner product on , Real and complex inner product spaces, with the inner product linear in the first argument.
If is differentiable on an interval with for all interior , then is nondecreasing; and the product rule holds; for real On an interval , for continuous on and differentiable at every interior point: throughout gives nondecreasing, gives increasing, and give the two decreasing forms; conversely a nondecreasing has and a nonincreasing has wherever it is differentiable, and no strict converse is claimed, Sums, scalar multiples, products and quotients: , , , and when , The chain rule, in one line from Carathéodory: if is differentiable at and is differentiable at , then is differentiable at with , The exponential function is smooth and .
For , is the inverse value of the strictly increasing exponential: (The natural logarithm as the inverse of the exponential function, The exponential is a continuous bijection from onto , The exponential function is strictly increasing).
A closed bounded subset of is compact by Heine-Borel in : with the Euclidean metric a subset of 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.
The Euclidean norm satisfies the triangle inequality, so ; hence each sphere is closed (Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation).
Proof
, so is an equilibrium. Put as in [L3]. Choose the second and third coordinates of the initial datum with signs and , respectively.
Let and let . Put and , so that and : indeed and by the definition of . Also as , since .
Let be the unique maximal solution with ; it is defined for all , because its invariant sphere is closed by [L7] and bounded, hence compact by Heine-Borel in : with the Euclidean metric a subset of 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 , so , and on one has . On that interval by [L2] (which gives on the invariant level set , since and ), so and has the constant sign on , its initial sign. Likewise cannot change sign on without vanishing, which would force contrary to ; since , one has there.
On , the first component equation of [L2] and step 3.1 give with , because .
By step 4.1 and the product and chain rules [L5], on ; since is continuous on , it is nondecreasing there by [L5], so for . If , this would give for , contradicting on ; hence .
At the time one has : indeed by step 5.1 and continuity, so by the definition of and continuity it is the alternative that occurs, and on the bound holds on with continuous, so ; moreover . Therefore , so .
It remains to make the initial point as close to as prescribed: for any choose with and , possible because by step 2.1. The corresponding solution then has and, by step 6.1, escapes the ball of radius at time . This is instability of in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system, with the explicit time bound of the statement.
Linearization of the torque-free Euler system at a principal-axis spin, and the oscillatory/exponential dichotomy
Statement
Let be pairwise distinct, let be the linear map with on the standard orthonormal basis, and let denote its inverse, given by . Define the formal torque-free Euler vector field When is an actual rigid-body inertia map written in a principal basis, this is the angular-velocity vector field obtained from Euler's equations 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 . Then:
- , and for every where for the explicit constant and ; consequently is the total derivative of at (The total (Fréchet) derivative as the linear first-order approximation with remainder).
- Every solution of the linear system on has constant, and with the pair is determined by its initial value as follows:
- if , then with , and every solution is bounded on ;
- if , then the same two formulas hold with and in place of and and , the system has unbounded solutions (for instance , ), and if then as .
- Hence the linearization at the spin is oscillatory exactly when , i.e. when the third principal moment is a largest or a smallest moment, and it has exponentially growing solutions exactly when , i.e. when lies strictly between and ; in the latter case the zero solution of the linear system is not stable.
Facts & Assumptions
Given: Pairwise distinct positive reals , the diagonal map , the field and a real .
Every has a unique expansion by the standard-basis result (The standard list with and for is an ordered basis of ; hence , and is the zero space with basis and dimension ). Define ; it is linear and direct substitution gives , so is the inverse linear map of (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 and cyclic analogues, and (The cross product in , The squared cross-product norm is the Gram determinant of two vectors).
For a matrix , with : this is the row-by-row Cauchy-Schwarz inequality Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation, Rectangular matrix multiplication and the identity matrix , including zero-sized shapes.
is a linear map with ; a linear map is totally differentiable at every point with itself as total derivative, and a bound exhibits as the total derivative by the definition The total (Fréchet) derivative as the linear first-order approximation with remainder.
The first-order linear system is locally Lipschitz in the state variable and has unique solutions with prescribed initial data; the same holds for 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.
Derivatives: , ; , ; the chain and sum rules hold; and The derivatives of sine and cosine are cosine and minus sine, Addition formulas, identities, parity, and derivatives of the hyperbolic functions, Sums, scalar multiples, products and quotients: , , , and when , The chain rule, in one line from Carathéodory: if is differentiable at and is differentiable at , then is differentiable at with , Signs, monotonicity intervals, and ranges of sine and cosine.
and ; the exponential is strictly increasing with image The six hyperbolic functions and their natural domains, The exponential function is strictly increasing, The exponential is a continuous bijection from onto .
Proof
Expanding with [L1], , and while ; collecting components gives . Applying and [L1] yields with as displayed, and .
Let solve . The third equation gives , so is constant; the first two equations read , with and by distinctness and , so and ; conversely every solution of produces the solution .
The remainder satisfies by [L1] and two applications of [L2]; hence obeys the bound of clause 1 with , and is the total derivative of at by [L3].
Suppose and put . The function satisfies by the derivative formulas and chain rule [L5]; the pair solves the first-order system , which is locally Lipschitz, so by uniqueness [L4] it is the solution with the given initial data. Taking and using , , gives and ; substituting , and gives the two displayed formulas. Boundedness follows from [L5].
Suppose and put . The same argument with in place of verifies the two displayed formulas: satisfies by [L5], and uniqueness [L4] identifies it as the solution.
Let and set , , so that by [L6]. If then as , because and [L6]; since , this is the growth criterion of clause 2, and the choice satisfies it, so the system has unbounded solutions and clause 3 follows.
Finally the sign of is the sign of because ; this vanishes exactly when equals or , which distinctness excludes, so , and exactly for extreme, exactly for intermediate , completing the sign criterion. In the intermediate case, for every choose and initial data . Step 3.1 gives an unbounded solution with initial norm ; by continuity it reaches norm at some . Thus the fixed 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.
The coordinate rotation matrices and : values in , composition law, and derivatives
Statement
For define, in the standard basis of , Then:
- and for every (The orthogonal group and the rotation group ), and , and .
- and for all .
- For every differentiable on an interval (The derivative of at a point that is a limit point of , and differentiability on a set, Intervals of : the nine order-convex forms, nondegeneracy, and length), the matrix-valued curves and are differentiable with where are the cross-product matrices of The skew-matrix/cross-product correspondence and its conjugation law .
Facts & Assumptions
Given: Reals , matrices as displayed, and a differentiable .
and are differentiable with , , and satisfy the addition formulas , and the Pythagorean identity ; also and 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.
Matrix product and transpose are entrywise operations: and ; and these elements form a subgroup of The vector space of by matrices over a field, with entrywise operations, Rectangular matrix multiplication and the identity matrix , including zero-sized shapes, The transpose of a matrix, Transpose is linear and involutive, and , The orthogonal group and the rotation group , is a subgroup of the general linear group, with .
The determinant is multiplicative and ; 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, , For , the determinant over a commutative ring by the Leibniz formula, and for a real matrix.
The chain rule and the algebra of derivatives hold for real functions; a matrix-valued curve is differentiable exactly when all entries are, and its derivative is the entrywise derivative The chain rule, in one line from Carathéodory: if is differentiable at and is differentiable at , then is differentiable at with , Sums, scalar multiples, products and quotients: , , , and when , The derivative of at a point that is a limit point of , and differentiability on a set, The vector space of by matrices over a field, with entrywise operations.
and are the skew matrices of and ; in particular and The skew-matrix/cross-product correspondence and its conjugation law , The vector space of by matrices over a field, with entrywise operations.
Proof
Entrywise, has diagonal entries , , and off-diagonal entries ; hence , and expanding the determinant along the third row gives , so , with read off the third column.
Similarly and , so , and is its second column. because and .
The entry of is , the entry is , the entry is , the entry is , and the remaining entries are ; hence by [L1].
The analogous four planar entries of equal and in the corresponding positions, so by [L1].
If is differentiable, entrywise differentiation using [L1] and [L4] gives Direct multiplication of with from [L5] gives the same matrix: its entries are respectively , and all remaining entries are zero. Thus .
The same comparison for : the derivative matrix has entries at , at , at , at and in the second row and column, which is exactly by [L5]; hence .
Since commutes with and commutes with , as direct entrywise multiplication using [L1] shows, the two alternative forms and of clause 3 hold as well.
The Euler-angle parametrization and the heavy symmetric top system
Definition
Euler-angle parametrization. Let be the coordinate rotation matrices of The coordinate rotation matrices and : values in , composition law, and derivatives. The Euler-angle parametrization is the map which is -valued by The coordinate rotation matrices and : values in , composition law, and derivatives and is a subgroup of the general linear group, with . The numbers are the Euler angles and the open set is its domain; is restricted to the open interval so that , 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 , , , and ; they are mathematical parameters, not identified with any physical quantity in this item.
The heavy-top function. Write for a point of the domain and for a triple of reals. The heavy-top function is 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: , , , and when , The chain rule, in one line from Carathéodory: if is differentiable at and is differentiable at , then is differentiable at with ) gives continuous partial derivatives of every order, so it is on the domain times ( Euclidean maps and diffeomorphisms); its relevant partial derivatives are, explicitly,
The heavy-top system. Let be an interval with at least two points. A heavy-top trajectory for these parameters is a twice-differentiable curve such that, with its derivative, the three second-order equations hold on , all functions being evaluated at . 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 and are constant along the trajectory.
Remarks
-
The system as displayed is what later items use. The three equations above are the system studied in Reduction of the heavy symmetric top system: conserved momenta, one-dimensional equation and conserved reduced energy; that item derives the conservation laws and the reduced one-dimensional equation directly from these displayed equations. Naming records the intended reading, but no variational derivation is invoked anywhere on this page, and none of the claims below depends on one.
-
The angular-velocity interpretation is a separate lemma. That is the kinetic-minus-potential reading of a rotating body with principal moments under a uniform field, and that the variables are the same angles that parametrize the rigid placement, is the content of Angular velocity in Euler angles: and and of the physical items that consume this pair, not of this definition.
-
No physical law is asserted here. That a physical symmetric top satisfies such a system, and what the parameters mean physically, are postulates or derived statements of the consuming physics pages.
Reduction of the heavy symmetric top system: conserved momenta, one-dimensional equation and conserved reduced energy
Statement
Fix the parameters , , , , and let be a heavy-top trajectory of The Euler-angle parametrization and the heavy symmetric top system on an interval with at least two points; write for its derivative. Then:
- the functions and are constant on ; write for their values;
- consequently and on ;
- the reduced energy is constant on ;
- with the effective potential on , one has so that the angle alone obeys the one-dimensional second-order equation with the constants determined by the initial data.
Facts & Assumptions
Given: Parameters , a heavy-top trajectory on with , and the partial derivatives of recorded in The Euler-angle parametrization and the heavy symmetric top system.
The system defining a heavy-top trajectory is , , , with , , and The Euler-angle parametrization and the heavy symmetric top system.
Sum, product and quotient rules hold for differentiable real functions; differentiation is linear Sums, scalar multiples, products and quotients: , , , and when .
The derivative of is computed by the chain and product rules, and does not depend on explicitly other than through The chain rule, in one line from Carathéodory: if is differentiable at and is differentiable at , then is differentiable at with , Sums, scalar multiples, products and quotients: , , , and when , The derivative of at a point that is a limit point of , and differentiability on a set, Intervals of : the nine order-convex forms, nondegeneracy, and length.
Proof
By the first and third equations of [L1], the functions and are constant on ; the second of these is , and the first equals . Hence and are constants by A function continuous on an interval whose derivative vanishes at every interior point of is constant on ; consequently two such functions with the same derivative differ by a constant applied to their continuous differentiable restrictions to , which is clause 1.
Define along the trajectory. By [L4] the derivative of is , and substituting the three system equations of [L1] turns this into ; hence . Since is continuous on the interval and differentiable at its interior points, A function continuous on an interval whose derivative vanishes at every interior point of is constant on ; consequently two such functions with the same derivative differ by a constant makes it constant on .
Compute : writing and , the quotient rule [L3] and [L2] give , which after factoring is ; and the remaining terms contribute , so .
Solving the two defining identities: gives , and gives because by [L2]; substituting the second into the first gives the displayed formula for . This is clause 2.
Evaluating with the partial derivatives of [L1] and of the definition gives the cross terms cancelling after collecting and the analogous axial term; substituting and from steps 1.1 and 2.1 gives the displayed formula for . Since is constant by step 1.2, this is clause 3.
Substituting from step 2.1 into the first term of step 1.3 gives , hence ; by step 1.1 and the identity of [L1] with , the right-hand side equals . Since the second system equation reads and , we obtain on , which is clause 4.
Angular velocity in Euler angles: and
Statement
Let be an interval with at least two points, let be differentiable with , and put the Euler-angle parametrization of The Euler-angle parametrization and the heavy symmetric top system. Then is a differentiable -valued curve, and its body angular velocity in the sense of Rigid placements, body and spatial angular velocity, and material velocities (the unique vector with ) is, in the standard basis, all functions evaluated at . Consequently Moreover, if and is the linear map with , , , then .
Facts & Assumptions
Given: Differentiable with , the curve , and the matrices of The coordinate rotation matrices and : values in , composition law, and derivatives.
for all , , , and for differentiable one has and The coordinate rotation matrices and : values in , composition law, and derivatives.
The product rule for matrix-valued curves, the skewness of for an -valued curve, and the equivalence hold; moreover and Derivative rules for differentiable matrix-valued curves, and skewness of for a curve in , The skew-matrix/cross-product correspondence and its conjugation law , Rigid placements, body and spatial angular velocity, and material velocities.
and ; the algebra of derivatives holds The derivatives of sine and cosine are cosine and minus sine, The chain rule, in one line from Carathéodory: if is differentiable at and is differentiable at , then is differentiable at with , Sums, scalar multiples, products and quotients: , , , and when , Sine and cosine defined by their real power series. For the identity used below, [L1] gives , hence by [L3]. Its first column is , so the first diagonal entry of this product is by the matrix-product formula The coordinate rotation matrices and : values in , composition law, and derivatives, Rectangular matrix multiplication and the identity matrix , including zero-sized shapes.
Evaluation of , , and , is the column reading of the displayed matrices The coordinate rotation matrices and : values in , composition law, and derivatives.
Proof
By [L1] each of the three factors is a differentiable -valued curve; by the product rule [L2] and closure of under multiplication, is differentiable with values in . Hence is skew and there is a unique with by [L2].
By the conjugation identity of [L2] with and , and the column readings of [L5],
Differentiating the product and applying [L1],
Multiplying on the left by and using [L3] term by term,
Conjugating the first vector of step 1.2 by and using [L5], ; substituting into step 3.1 gives with exactly as displayed in the statement.
Reading off the components gives ; expanding the squares of the first two components and using from [L4] cancels the mixed terms and gives , whence .
For the diagonal map with , , , the quadratic form is , and substituting step 5.1 gives the displayed identity.
Remarks
- Relation to the kinetic energy. Combined with clause 2 of Kinetic energy and angular momentum of a rigid placement: and , the last identity says that a rigid placement with principal moments has rotational kinetic energy , which is the kinetic part of the heavy-top function of The Euler-angle parametrization and the heavy symmetric top system. No statement about a physical top is made by that comparison; it is the consuming physical items that identify the parameters.
5 · Examples, experiments and counterexamples
None yet.
Sources
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.1.1, printed pp. 47-49: the moving body basis and the antisymmetry computation $R^{\mathsf T}\dot R+\dot R^{\mathsf T}R=0$
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.2 (Kinematics of Rotation), Eq. (2.21): orientation matrix and skew angular-velocity matrix
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.1, Eqs. (4.1)-(4.10): moving frame, antisymmetric angular-velocity matrix and the transport theorem
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.1.1, Eqs. (3.5)-(3.7), printed pp. 47-48: the antisymmetric angular-velocity matrix and its single-index vector
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.2: extracting the angular-velocity vector from an antisymmetric matrix and rotation-covariance of the cross-product identification
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.2, Eqs. (3.19)-(3.22), printed pp. 50-52: inertia-tensor components of a discrete or continuous mass distribution and the symmetry/positivity computation
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.2, Eqs. (4.16)-(4.22): continuous inertia-tensor integral and angular momentum
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.4 (Inertia Tensor), Eq. (2.41) and §2.6 (Vector Angular Momentum)
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.2, Eq. (3.19) and the parallel-axis discussion of §3.2.1, printed pp. 50-53: second moments of a bounded mass distribution are finite
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.2, Eq. (4.16): the inertia-tensor integral of a bounded density
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §§3.1-3.2, Eqs. (3.3)-(3.10) and (3.29)-(3.30), printed pp. 47-54: rigid placement $r(t)=\tilde R(t)+\Delta r(t)$, body angular velocity and centre-of-mass splitting
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §§4.1-4.2, Eqs. (4.1)-(4.10): rigid body as six degrees of freedom, moving frame and angular velocity
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.2 (Kinematics of Rotation): orientation path $M(q(t))$, velocity $\dot M M^{\mathsf T}$ and body components $M^{\mathsf T}\omega$
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.2, Eqs. (3.18)-(3.22), printed pp. 50-52: symmetry of the inertia tensor, positivity and the equality case (all mass on a line)
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.2, Eqs. (4.16)-(4.18): symmetry and the quadratic form of the inertia tensor
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §§2.3-2.5: moments of inertia, inertia tensor and principal moments as a symmetric matrix diagonalized by a rotation
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §§3.2-3.3, Eqs. (3.17)-(3.20) and (3.27)-(3.30), printed pp. 50-54: kinetic energy $T=\tfrac12\omega^{\mathsf T}I\omega$, angular momentum $L=I\omega$, and the centre-of-mass split for a free body
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.2, Eqs. (4.14)-(4.22): rotational kinetic energy and angular momentum through the inertia tensor
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.1 Eq. (2.10), §2.4 Eq. (2.41) and §2.6 Eq. (2.51): kinetic-energy split, inertia tensor and vector angular momentum
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.2, the diagonalization discussion after Eq. (3.20), printed pp. 50-52: principal axes and principal moments via orthogonal diagonalization of the symmetric inertia tensor
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.5 (Principal Moments of Inertia): the eigensystem of the inertia tensor and the diagonal form
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.3: principal axes of the inertia tensor and their positivity
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.3, Eqs. (3.31)-(3.33), printed pp. 53-55: angular-momentum rate, Euler's equations in the principal body frame and the torque extension
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.4, Eq. (4.66): Euler equations $I_j\dot\omega_j+(I_{j''}-I_{j'})\omega_{j'}\omega_{j''}=\tau_j$
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.9 (Euler's Equations), Eqs. (2.80)-(2.83): angular-momentum evolution and Euler equations in principal components
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §§3.4.1-3.4.3, Eqs. (3.34)-(3.45), printed pp. 55-58: free symmetric top, and the asymmetric-top stability discussion with the invariants of the motion
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §§2.8.2-2.9: energy and angular-momentum level curves of the free rigid body and Euler's equations
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.4: free-rotation stability about the extreme and intermediate principal axes
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.4.3, Eqs. (3.40)-(3.45), printed pp. 57-58: linearized perturbation of principal-axis spin and the sign criterion for exponential growth; the nonlinear invariant-branch proof and escape bound are supplied locally
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §§2.8.2-2.9: intermediate-axis behaviour on the angular-momentum level curves
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.4: instability of rotation about the intermediate principal axis
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.4.3, Eqs. (3.41)-(3.45), printed pp. 57-58: perturbation of a principal-axis spin, the linearized equations and the sign criterion for oscillation versus exponential growth
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §§2.8.2-2.9: linear behaviour near principal-axis spins
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.5, Eqs. (3.53)-(3.58), printed pp. 62-65: the three ordered rotations and the Euler-angle parametrization of $SO(3)$
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.5, Eqs. (4.68)-(4.74): three successive rotations about coordinate axes, Euler angles and the angular-velocity components
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.7 (Euler Angles): the rotation matrices about the coordinate axes and their products
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.5 Eq. (3.54) and §3.6 Eqs. (3.62)-(3.70), printed pp. 62-68: the Euler-angle rotation, the heavy-top function and its two cyclic momenta and reduced energy
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.5, Eqs. (4.74)-(4.85): Euler angles for the symmetric top, the exact reduced effective potential and the upright threshold
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.10 (Axisymmetric Tops), Eqs. (2.84)-(2.91): the top Lagrangian, its momenta and reduced energy
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.6, Eqs. (3.63)-(3.70), printed pp. 67-69: the two cyclic momenta $p_\psi,p_\phi$, the reduced energy $E'$ and the effective potential
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.5, Eqs. (4.78)-(4.85): reduced energy and effective potential of the symmetric top
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.10, Eqs. (2.89)-(2.91): the momenta $p_\varphi,p_\psi$ and the energy expressed through the effective potential
- D. Tong, Classical Dynamics (Cambridge lecture notes, 2015), §3.5.2, Eqs. (3.55)-(3.58), printed pp. 64-65: angular velocity in terms of Euler angles, with $\omega_3=\dot\psi+\dot\varphi\cos\theta$ and the squared transverse component
- K. Likharev, Essential Graduate Physics: Classical Mechanics, Ch. 4 §4.5, Eq. (4.74): principal-axis components of the instantaneous angular velocity in Euler angles
- G. J. Sussman and J. Wisdom, Structure and Interpretation of Classical Mechanics, 2nd ed. (MIT Press, 2014), §2.7 (Euler Angles) and §2.10, Eq. (2.86): the kinetic energy $\tfrac12A(\dot\theta^2+\dot\varphi^2\sin^2\theta)+\tfrac12C(\dot\psi+\dot\varphi\cos\theta)^2$