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 Variations — Examples

1 · Prerequisites

2 · Summary

The companion page to cm-math-variations records explicit worked checks for its variational, Legendre and constraint results. Everything here is mathematics: the examples carry domain: mathematics, depend only on the mathematical items of the parent page and on mathematical suppliers, and make no physical or empirical claim. The models are stated with their adopted data (carrier, Lagrangian, constraint, mass matrix, force); each conclusion is a conditional deduction from those data.

A stationary path that is not a local minimum: the quadratic model L=12(v2−4q2) on [0,π] verifies the scope of the stationarity theorem with the indefinite quadratic Lagrangian L=12(v2−4q2) on [0,π]: the curve q∗(t)=sin⁡2t is stationary and the exact identity S(q∗+ϵsin⁡)=S(q∗)−3π4ϵ2 shows that it is not a local minimum and that the model has no minimiser. Regular versus hyperregular: the bounded-gradient Lagrangian 1+v2 and the degenerate zero Lagrangian separates regularity from hyperregularity: L1(v)=1+v2 has a strictly positive velocity Hessian but its velocity gradient has range (−1,1), the zero Lagrangian is neither regular nor hyperregular, and 12v2 is both.

The remaining three examples exercise the constraint machinery. The ideal planar pendulum: equation of motion, conserved energy, Lyapunov stability and the frequency g/l derives ml2θ′′+mglsin⁡θ=0 from the reduced d'Alembert equation, proves conservation of the reduced energy, Lyapunov stability of the hanging equilibrium through the energy-stability theorem, and the linearised frequency g/l. The rotating hoop: reduced Lagrangian, the equation θ′′=Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ, and the reduced energy versus the inertial mechanical energy treats a circle rotating about the vertical axis with prescribed constant angular speed: differentiating the carrier gives the velocity and acceleration identities, the reduced Euler–Lagrange equation θ′′=Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ, conservation of the reduced energy, the comparison with the inertial mechanical energy K+Ug (whose rate of change is the reaction power under an ideal reaction), and the subcritical stability threshold aΩ2≤g. The planar rolling disk: the holonomic constraint x−aθ=const and the acceleration f/(m+I/a2) models the planar rolling disk as a regular holonomic ideal-constraint system: the constraint x−aθ=const integrates, the projected equation gives (m+I/a2)x′′=f with multiplier λ=−fI/(I+ma2), so that x′′=2f/(3m) for a uniform thin disk and x′′=f/(2m) for a hoop, and an explicit quadratic motion realises the acceleration. The scope remarks on each example separate the ideal model from its physical realisation.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

None yet.

5 · Examples, experiments and counterexamples

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

A stationary path that is not a local minimum: the quadratic model L=12(v2−4q2) on [0,π]

Example

Let n=1, [a,b]=[0,π], U=R, and L(q,v)=12(v2−4q2), a C∞ (hence C2) Lagrangian on R×R×R (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations, Ck Euclidean maps and diffeomorphisms). Let

q∗(t):=sin⁡2t,η(t):=sin⁡t.

Use fixed-endpoint admissible C2 curves with q(0)=q(π)=0 and the norm ∥q∥C2:=max⁡0≤j≤2sup⁡t∈[0,π]∣q(j)(t)∣. A local minimum means S(q∗)≤S(q) for all admissible q in some ball about q∗ in this norm. Then q∗ is a stationary curve of the action S(q)=∫0π12(q′2−4q2)dt, and for the one-parameter family q∗+ϵη one has the exact identity

S(q∗+ϵη)=S(q∗)−3π4ϵ2(ϵ∈R).

In particular S(q∗+ϵη)<S(q∗) for every ϵ≠0: the stationary curve is not a local minimum of the action. The model has no minimiser at all, since S(λη)=−3π4λ2→−∞ as ∣λ∣→∞.

Facts & Assumptions

Given: The Lagrangian L(q,v)=12(v2−4q2) on [0,π], the curves q∗(t)=sin⁡2t and η(t)=sin⁡t, and the action S of A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations.

[F1]

The functions sin⁡ and cos⁡ are differentiable on R with sin⁡′=cos⁡, cos⁡′=−sin⁡, sin⁡0=0, cos⁡0=1 (The derivatives of sine and cosine are cosine and minus sine, Sine and cosine defined by their real power series); the chain rule gives q∗′(t)=2cos⁡2t and q∗′′(t)=−4sin⁡2t (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The addition formulas cos⁡(2t)=cos⁡2t−sin⁡2t, sin⁡(2t)=2sin⁡tcos⁡t and the Pythagorean identity cos⁡2t+sin⁡2t=1 hold for all real t (The addition formulas for sine and cosine, Parity and the Pythagorean identity for sine and cosine), and sin⁡(2π)=0 because 2π is a period of sin⁡ (The zero sets of sine and cosine and the least positive common period 2 pi).

[F2]

The Euler–Lagrange expression of L along a C2 curve q is EL(q)=−q′′−4q, because Lq=−4q and Lv=q′; stationarity of a C2 curve for S is equivalent to EL(q)=0 (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations, Stationarity of the time-dependent action is equivalent to the vanishing of the Euler-Lagrange expression, The first variation of a time-dependent action on a finite interval: admissibility, the integral formula and continuity of the Euler-Lagrange coefficient).

Verification

technique · direct; check the Euler–Lagrange equation and expand the quadratic action exactly
1.1F1F2given

The curve q∗ is C2 on an open interval containing [0,π] with q∗(0)=q∗(π)=0, and EL(q∗)=−q∗′′−4q∗=4sin⁡2t−4sin⁡2t=0; hence q∗ is stationary for S by [F2].

2.1F1F2F3step 1.1algebra

For every ϵ the curve q∗+ϵη is C2 with zero endpoint values, and the integrand of its action is a quadratic polynomial in ϵ: 12((q∗′+ϵη′)2−4(q∗+ϵη)2)=12(q∗′2−4q∗2)+ϵ(q∗′η′−4q∗η)+ϵ22(η′2−4η2). By [F3] the action is the sum of the integrals of the three summands, and the linear term is ϵ δS(q∗)[η]=0 by step 1.1 and [F2].

3.1F1F3step 2.1algebra

The quadratic term: with η=sin⁡, η′=cos⁡ and η′2−4η2=cos⁡2t−4sin⁡2t. By [F1], cos⁡2t=12(1+cos⁡2t) and sin⁡2t=12(1−cos⁡2t), so cos⁡2t−4sin⁡2t=12(1+cos⁡2t)−2(1−cos⁡2t)=−32+52cos⁡2t. Integrating over [0,π] and using ∫0πcos⁡2t dt=0 from [F3] and ∫0π1 dt=π gives ∫0π(η′2−4η2)dt=−3π2. Substituting into step 2.1 gives S(q∗+ϵη)=S(q∗)−3π4ϵ2 for every real ϵ.

4.1F1F2F3step 2.1step 3.1algebra

A positive direction is ζ(t):=sin⁡3t, which vanishes at both endpoints. The same expansion as step 2.1 has zero linear term by stationarity and gives S(q∗+ϵζ)−S(q∗)=12ϵ2∫0π(9cos⁡23t−4sin⁡23t)dt=5π4ϵ2. Use the double-angle formulas and ∫0πcos⁡6t dt=[sin⁡6t/6]0π=0, by periodicity from [F1]. Thus arbitrarily small nonzero variations in this direction increase the action, whereas step 3.1 decreases it in direction η.

5.1step 3.1algebra∎

Since ∥ϵη∥C2=∣ϵ∣∥η∥C2→0 and −3π4ϵ2<0 for every ϵ≠0, arbitrarily small nonzero variations of the stationary curve strictly decrease the action; thus q∗ is not a local minimum. Moreover S(λη)=−3π4λ2 tends to −∞, so the model has no minimising curve at all.

Remarks

  • The example is a mathematical model with an indefinite quadratic Lagrangian; it shows that the Euler–Lagrange equation is a necessary condition only, and that "principle of least action" is a misnomer in general. The stationary curve here is a saddle in the infinite-dimensional sense realised by the decreasing family q∗+ϵsin⁡t and the increasing family q∗+ϵsin⁡3t of steps 3.1 and 4.1.

  • The zero curve q=0 is also stationary for this model, and it is likewise not a local minimum, by the same computation with q∗ replaced by 0. Nothing here depends on the interpretation of L as a physical energy; no positivity of the kinetic term is claimed.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

Regular versus hyperregular: the bounded-gradient Lagrangian 1+v2 and the degenerate zero Lagrangian

Example

Take n=1, J=U=R, and consider the three fibre-full C2 Lagrangians L1(v)=1+v2, L2(v)=0 and L3(v)=12v2 on J×U×R (The Legendre map of a time-dependent Lagrangian: regularity and hyperregularity). Then:

  1. L1 is regular on the chart but not hyperregular: its velocity Hessian is (1+v2)−3/2>0 at every v, while the fibre Legendre map v↦L1′(v)=v/1+v2 has range (−1,1), so it is not surjective onto R;
  2. L2 is neither regular nor hyperregular: L2′≡0 and L2′′≡0, so the velocity Hessian is singular everywhere and the fibre map is constant, hence not injective;
  3. L3 is both regular and hyperregular: L3′′≡1 and the fibre map v↦v is a C∞ bijection of R onto itself.

Thus regularity does not imply hyperregularity, and neither condition is automatic.

Facts & Assumptions

Given: The three C2 Lagrangians L1(v)=1+v2, L2(v)=0, L3(v)=12v2, their fibre maps v↦Lj′(v), and their second derivatives.

[F3]

Elementary estimates: 0<1+v2 and ∣v(1+v2)−1/2∣=∣v∣/1+v2<1 for every real v, because ∣v∣<1+v2 follows from squaring; in particular p=1 is not in the range of φ1, since v/1+v2=1 would give v≥0 and v2=1+v2, i.e. 0=1.

Verification

technique · direct; compute the derivatives, compare the ranges with $\mathbb R$, and check regularity and hyperregularity against the definition
1.1F1algebra

By [F1] the three Lagrangians are C2, and the velocity Hessians at v are (1+v2)−3/2, 0 and 1 respectively. Hence L1 and L3 are regular on the chart, while L2 is regular at no point, because regularity means a nonsingular one-by-one velocity Hessian (The Legendre map of a time-dependent Lagrangian: regularity and hyperregularity, Nonsingular velocity Hessian is equivalent to invertibility of the Legendre derivative, hence gives a local C^1 inverse).

2.1F2F3step 1.1

By [F2], φ1 is injective and φ3 is bijective R→R, while φ2 is constant. Moreover, by [F3] the range of φ1 is contained in (−1,1) and omits p=1; for every p∈(−1,1), setting v=p/1−p2 gives 1+v2=1/(1−p2) and v/1+v2=p. Thus the range is exactly (−1,1), and φ1 is not surjective onto R. Therefore L1 is regular but not hyperregular, L3 is hyperregular, and L2 fails hyperregularity by non-injectivity (and fails regularity as well).

3.1step 1.1step 2.1∎

The three clauses of the Example are proved: clause 1 by the regularity of L1 from step 1.1 and the failure of surjectivity from step 2.1; clause 2 by the vanishing Hessian and constancy from steps 1.1 and 2.1; clause 3 by the invertibility of the identity map and the nonvanishing Hessian.

Remarks

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

The ideal planar pendulum: equation of motion, conserved energy, Lyapunov stability and the frequency g/l

Example

Let m,l,g>0 and let Ω:=(−π,π)⊆R. Consider the embedding

X:Ω→R2,X(θ):=(lsin⁡θ, −lcos⁡θ),

the ambient mass M:=mI2 and the constant force F:=(0,−mg). Let I⊆R be a nonempty open interval (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let θ:I→Ω be a C2 curve and put x(t):=X(θ(t)). Then:

  1. X is a C3 immersion with ∥X′(θ)∥2=l for every θ; if the ambient motion satisfies Mx′′=F+R with R(t)⋅X′(θ(t))δ=0 for all δ∈R, then

ml2θ′′(t)+mglsin⁡θ(t)=0(t∈I);

  1. the reduced energy E(θ,θ′):=12ml2θ′2−mglcos⁡θ is constant along every solution of the equation;
  2. θ=0 is a strict local minimum of the reduced potential Ured(θ):=−mglcos⁡θ on Ω, and the equilibrium (0,0) of the pendulum system θ′=v, v′=−(g/l)sin⁡θ is Lyapunov stable;
  3. the linearisation at θ=0 is ml2ξ′′+mglξ=0, a quadratic mode with angular frequency

ω=g/l.

Facts & Assumptions

Given: The parameters m,l,g>0, the embedding X, the matrix M=mI2, the force F=(0,−mg), the nonempty open interval I, a C2 curve θ:I→Ω with x=X∘θ, and a reaction R as in clause 1.

[F1]

Derivatives of X: since sin⁡′=cos⁡, cos⁡′=−sin⁡ (The derivatives of sine and cosine are cosine and minus sine), X′(θ)=(lcos⁡θ,lsin⁡θ) and X′′(θ)=(−lsin⁡θ,lcos⁡θ); hence ∥X′(θ)∥22=l2, X′(θ)⋅X′′(θ)=0 and X′(θ)⋅F=−mglsin⁡θ (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). Since ∥X′(θ)∥2=l≠0, the derivative X′(θ) is injective and X is an immersion (Submersions and immersions between Euclidean open sets).

[F2]

The reduced d'Alembert corollary (Reduced d'Alembert equation for an immersion: generalized forces and the quadratic correction): under the hypothesis on R in clause 1, X′(θ)TMX′(θ) θ′′=X′(θ)TF−X′(θ)TMX′′(θ) θ′2.

[F3]

The batch-3 energy-stability theorem (A strict local minimum of the potential gives a Lyapunov stable equilibrium of the mechanical system, by energy conservation, The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart): for an open Ω⊆R, a constant positive mass M0, a C2 potential U and a strict local minimum θ0 of U with U′(θ0)=0, the equilibrium (θ0,0) of θ′=v, v′=−U′(θ)/M0 is Lyapunov stable.

[F4]

The quadratic-mode theorem (Decoupling of the constant-coefficient quadratic system Mq′′+Kq=0 into normal modes, with the zero-mode and energy identities): for a symmetric positive definite mass M and a symmetric non-negative stiffness K, the linear system Mξ′′+Kξ=0 has modes of frequency λ with λ the eigenvalues of M−1/2KM−1/2.

Verification

technique · direct; contract the ambient equation with $X'$ and identify the resulting one-dimensional mechanical system, then invoke the stability and mode theorems
1.1F1F2algebra

Clause 1: by [F2] and [F1], X′TMX′ θ′′=l2m θ′′=X′TF−X′TMX′′θ′2=−mglsin⁡θ−0⋅θ′2, because M=mI2 and X′⋅X′′=0. Dividing by ml2>0 gives ml2θ′′+mglsin⁡θ=0.

2.1F5step 1.1algebra

Clause 2: along a solution, by the product and chain rules [F5], ddtE(θ,θ′)=ml2θ′θ′′+mglsin⁡θ θ′=θ′(ml2θ′′+mglsin⁡θ)=0, so by [F5] the energy is constant on I.

2.2F5step 1.1algebra

Clause 3: the reduced potential Ured(θ)=−mglcos⁡θ is C2 with Ured′(θ)=mglsin⁡θ, so Ured′(0)=0; and for 0<∣θ∣<π one has θ∉2πZ, hence cos⁡θ<1 by [F5], so Ured(θ)−Ured(0)=mgl(1−cos⁡θ)>0. Thus 0 is a strict local minimum of Ured on Ω. With M0:=ml2>0, U:=Ured and θ0:=0, the mechanical vector field of [F3] is (θ,v)↦(v, (ml2)−1(−mglsin⁡θ))=(v,−(g/l)sin⁡θ), which is the pendulum system; [F3] makes (0,0) Lyapunov stable.

2.3F4step 1.1algebra

Clause 4: the right-hand side of the pendulum equation is −(g/l)sin⁡θ; since sin⁡0=0 and sin⁡′(0)=cos⁡0=1 (The derivatives of sine and cosine are cosine and minus sine), the linearisation at θ=0 is ml2ξ′′+mglξ=0. This is the constant-coefficient quadratic system of [F4] with M=ml2 and K=mgl; the transformed stiffness is λ=K/M=g/l>0, so [F4] gives a mode of angular frequency ω=λ=g/l.

3.1step 1.1step 2.1step 2.2step 2.3∎

Clauses 1-4 are proved.

Remarks

  • The example is a mathematical model of an ideal frictionless planar pendulum in a uniform downward force field; statements about a real pendulum require apparatus, friction and measurement assumptions that are not part of this item. The stability conclusion is Lyapunov stability, not asymptotic attraction, and it concerns the exact pendulum equation in the chart Ω=(−π,π).

  • The frequency g/l is the frequency of the linearised mode by [F4]; no claim is made that the nonlinear pendulum oscillates periodically with that frequency for finite amplitude.

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

The rotating hoop: reduced Lagrangian, the equation θ′′=Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ, and the reduced energy versus the inertial mechanical energy

Example

Let a,m,g>0 and Ω∈R, and let

X:R×R→R3,X(θ,t):=(asin⁡θcos⁡(Ωt), asin⁡θsin⁡(Ωt), a(1−cos⁡θ)),

the position of a point of a circle of radius a rotating about the vertical axis with constant angular speed Ω. For an open interval I and a C2 curve θ:I→R put x(t):=X(θ(t),t). Stationarity on I means stationarity of the restricted curve on every [t1,t2]⊆I with t1<t2, using fixed-endpoint variations. Then:

  1. x˙=Xθθ′+Xt and 12m∣x˙∣2=12ma2(θ′2+Ω2sin⁡2θ), while x¨⋅Xθ=a2θ′′−a2Ω2sin⁡θcos⁡θ;
  2. the curves θ that are stationary for the action of the time-dependent Lagrangian chart L(θ,v):=12ma2(v2+Ω2sin⁡2θ)−mga(1−cos⁡θ) on O:=R×R×R are exactly the C2 solutions of

θ′′=Ω2sin⁡θcos⁡θ−gasin⁡θ;

  1. the reduced energy E(θ,v):=12ma2v2−12ma2Ω2sin⁡2θ−mgacos⁡θ is constant along every C2 solution, whereas the inertial kinetic-plus-potential energy K+Ug=12m∣x˙∣2+mga(1−cos⁡θ) satisfies

ddt(K+Ug)=2ma2Ω2sin⁡θcos⁡θ θ′=ma2Ω2ddtsin⁡2θ

along every solution, and equals E+ma2Ω2sin⁡2θ+mga; it is therefore in general not conserved, and it differs from E by ma2Ω2sin⁡2θ+mga; 4. if the ambient motion is additionally modelled by mx¨=(0,0,−mg)+R with an ideal reaction satisfying R(t)⋅Xθ(θ(t),t)=0, then the same equation of part 2 holds, and the reaction power is

R⋅x˙=ddt(K+Ug)=2ma2Ω2sin⁡θcos⁡θ θ′,

generally nonzero: the prescribed rotation can do work on the point; 5. if aΩ2≤g, then 0 is a strict local minimum of Ueff(θ):=−12ma2Ω2sin⁡2θ−mgacos⁡θ with Ueff′(0)=0, and the equilibrium (0,0) of θ′=v, v′=Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ is Lyapunov stable; if aΩ2>g, then 0 is not a local minimum of Ueff and the stability criterion does not apply.

Facts & Assumptions

Given: The constants a,m,g>0 and Ω∈R, the map X displayed in the Example, an open interval I (Intervals of R: the nine order-convex forms, nondegeneracy, and length) and a C2 curve θ:I→R with x:=X(θ(⋅),⋅); the Lagrangian L and the reduced potential Ueff of the Example.

[F1]

Derivatives and inner products of the carrier (The derivatives of sine and cosine are cosine and minus sine, Sine and cosine defined by their real power series, Parity and the Pythagorean identity for sine and cosine, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a), 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 Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral, Ck Euclidean maps and diffeomorphisms): X is C∞, and with c:=cos⁡(Ωt), s:=sin⁡(Ωt), Xθ=(acos⁡θ c, acos⁡θ s, asin⁡θ),Xt=(−aΩsin⁡θ s, aΩsin⁡θ c, 0), Xθθ=(−asin⁡θ c, −asin⁡θ s, acos⁡θ),Xθt=(−aΩcos⁡θ s, aΩcos⁡θ c, 0),Xtt=(−aΩ2sin⁡θ c, −aΩ2sin⁡θ s, 0), and the six identities ∣Xθ∣2=a2,∣Xt∣2=a2Ω2sin⁡2θ,Xθ⋅Xt=0, Xθ⋅Xθθ=0,Xθ⋅Xθt=0,Xθ⋅Xtt=−a2Ω2sin⁡θcos⁡θ hold; each is expansion of the displayed components together with cos⁡2+sin⁡2=1.

[F2]

Chain rule along the curve (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)): if x(t)=X(θ(t),t) with θ of class C2, then x˙=Xθθ′+Xt and x¨=Xθθθ′2+2Xθtθ′+Xtt+Xθθ′′, all derivatives being evaluated at (θ(t),t).

[F3]

The chart and its Euler–Lagrange expression (A time-dependent Lagrangian chart, its action functional and its fixed-endpoint first variations): L is C∞ on O=R×R×R, so (O,L) is a time-dependent Lagrangian chart over every [t1,t2]×R; its gradients are Lθ=ma2Ω2sin⁡θcos⁡θ−mgasin⁡θ, Lv=ma2v, and its Euler–Lagrange expression along a C2 curve is EL(θ)=ma2Ω2sin⁡θcos⁡θ−mgasin⁡θ−ma2θ′′. By Stationarity of the time-dependent action is equivalent to the vanishing of the Euler-Lagrange expression a C2 curve is stationary for the action if and only if EL(θ)=0.

[F5]

Ideal reactions (A regular holonomic constraint, its virtual displacements and the ideal-reaction law): a reaction that performs no virtual work on the moving constraint satisfies R(t)⋅δx=0 for every fixed-time virtual displacement; for the circle at time t the virtual displacements are the multiples of Xθ(θ(t),t), so the hypothesis used below is R⋅Xθ=0.

[F6]

The batch-3 energy-stability theorem (A strict local minimum of the potential gives a Lyapunov stable equilibrium of the mechanical system, by energy conservation, The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart): for a C2 potential U:R→R and q0 with U′(q0)=0 a strict local minimum of U, the equilibrium (q0,0) of q′=v, v′=−U′(q)/M0 with the constant M0=ma2>0 is Lyapunov stable.

Verification

technique · direct; compute the kinematics of the carrier, identify the reduced Lagrangian and its Euler–Lagrange equation, differentiate the two energies along solutions, and contract the ambient equation with the tangent vector for the reaction power
1.1F1F2algebra

By [F2] and [F1], x˙=Xθθ′+Xt and ∣x˙∣2=∣Xθ∣2θ′2+2θ′ Xθ⋅Xt+∣Xt∣2=a2θ′2+a2Ω2sin⁡2θ, which gives the velocity formula of part 1.

1.2F1F2algebra

By [F2] and [F1], x¨⋅Xθ=θ′′∣Xθ∣2+2θ′ Xθt⋅Xθ+θ′2 Xθθ⋅Xθ+Xtt⋅Xθ=a2θ′′−a2Ω2sin⁡θcos⁡θ, which is the acceleration formula of part 1.

1.3F3algebra

By [F3], EL(θ)=ma2Ω2sin⁡θcos⁡θ−mgasin⁡θ−ma2θ′′, so EL(θ)=0 if and only if ma2θ′′=ma2Ω2sin⁡θcos⁡θ−mgasin⁡θ, that is, if and only if θ′′=Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ; by the equivalence of [F3] the stationary curves are exactly the C2 solutions, which is part 2.

2.1F4step 1.3algebra

Let θ be a C2 solution. By [F4] the function t↦E(θ(t),θ′(t)) is differentiable on I with ddtE(θ,θ′)=ma2θ′θ′′−ma2Ω2sin⁡θcos⁡θ θ′+mgasin⁡θ θ′=θ′(ma2θ′′−ma2Ω2sin⁡θcos⁡θ+mgasin⁡θ)=0 by step 1.3; hence E(θ,θ′) is constant on I by [F4]. This is the first assertion of part 3.

2.2F1F5step 1.2algebra

Assume now the ambient model mx¨=(0,0,−mg)+R with R⋅Xθ=0 as in [F5]. Contracting with Xθ kills R, and (0,0,−mg)⋅Xθ=−mgasin⁡θ by [F1]; by step 1.2, m(a2θ′′−a2Ω2sin⁡θcos⁡θ)=−mgasin⁡θ, i.e. θ′′=Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ after division by ma2>0. This is the first assertion of part 4.

2.3F1F6step 1.3algebra

Part 5: Ueff′(θ)=−ma2Ω2sin⁡θcos⁡θ+mgasin⁡θ, so Ueff′(0)=0, and the mechanical system of [F6] with M0=ma2 and this potential is exactly θ′=v, v′=Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ; moreover the exact factorization Ueff(θ)−Ueff(0)=mga(1−cos⁡θ)−12ma2Ω2sin⁡2θ=ma(1−cos⁡θ)[g−12aΩ2(1+cos⁡θ)] together with cos⁡θ≤1 and 1−cos⁡θ=2sin⁡2(θ/2) (The addition formulas for sine and cosine), which with The zero sets of sine and cosine and the least positive common period 2 pi gives cos⁡θ=1  ⟺  θ∈2πZ shows that for aΩ2≤g one has Ueff(θ)>Ueff(0) for 0<∣θ∣<2π: the product is nonnegative and the bracket is at least g−aΩ2≥0, vanishing only if g=aΩ2 and cos⁡θ=1. Hence 0 is a strict local minimum of Ueff, and [F6] makes (0,0) Lyapunov stable. If aΩ2>g, the bracket is negative for ∣θ∣ small (continuity at θ=0, where it equals g−aΩ2), while 1−cos⁡θ>0 for small nonzero θ, so Ueff(θ)<Ueff(0) there and 0 is not a local minimum; [F6] gives no conclusion in that case.

3.1F1F4step 1.1step 2.1algebra

By steps 1.1 and 2.1, K+Ug=12ma2(θ′2+Ω2sin⁡2θ)+mga(1−cos⁡θ)=E(θ,θ′)+ma2Ω2sin⁡2θ+mga; along a solution, by [F4] and step 1.3, ddt(K+Ug)=ddt(ma2Ω2sin⁡2θ)=2ma2Ω2sin⁡θcos⁡θ θ′. The phase field (θ,v)↦(v,Ω2sin⁡θcos⁡θ−(g/a)sin⁡θ) is C1 on R2, hence locally Lipschitz (A C1 map is locally Lipschitz on compact coordinate subsets); Picard-Lindelöf local existence and uniqueness for first-order systems supplies a solution through (π/4,1). If Ω≠0 the displayed derivative is nonzero at the initial time of this solution, where it equals ma2Ω2, so K+Ug is not conserved along all solutions; as functions on the full phase space, K+Ug agrees with E up to the constant mga precisely when the term ma2Ω2sin⁡2θ is absent, i.e. when Ω=0. This is the second assertion of part 3.

4.1F1F2F4F5step 3.1step 2.2algebra

Under the hypotheses of step 2.2, K+Ug has derivative mx¨⋅x˙+mgasin⁡θ θ′ by [F4], and z˙=asin⁡θ θ′ for the third coordinate z=a(1−cos⁡θ) of x; substituting mx¨=(0,0,−mg)+R gives mx¨⋅x˙=−mgz˙+R⋅x˙, so ddt(K+Ug)=R⋅x˙−mgz˙+mgasin⁡θ θ′=R⋅x˙. Since x˙=Xθθ′+Xt and R⋅Xθ=0, also R⋅x˙=R⋅Xt; by step 3.1 the common value is 2ma2Ω2sin⁡θcos⁡θ θ′, generally nonzero. This is the second assertion of part 4.

5.1step 1.1step 1.2step 1.3step 2.1step 3.1step 2.2step 4.1step 2.3∎

Parts 1-5 are proved.

Remarks

  • The reduced energy is Jacobi's generalized energy of the chart: E=vLv−L with the additive constant mga dropped. Since L does not depend on t, E is conserved; the inertial energy K+Ug instead contains the term ma2Ω2sin⁡2θ of the prescribed rotation and is not conserved. This is the mechanism of the rheonomic constraint of The ideal holonomic system in a possibly time-dependent chart: explicit acceleration, constraint preservation, local existence and reaction power, whose reaction power R⋅q′=−Λ⋅Gt need not vanish; part 4 computes that power directly for this carrier.

  • The model is the ideal frictionless threading of a point on a circle rotating with prescribed constant angular speed; the reaction R and the power identity are properties of that ideal model. Statements about a real hoop (friction, forcing torque, deformations, measured frequencies) require apparatus and empirical assumptions that are not part of this item.

  • The stability clause is Lyapunov stability of the equilibrium of the exact reduced equation on the chart R; no asymptotic stability, and no analysis of the supercritical equilibria cos⁡θ=g/(aΩ2) or of the inverted point, is asserted.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passOpen item page →

The planar rolling disk: the holonomic constraint x−aθ=const and the acceleration f/(m+I/a2)

Example

Let a,m,I>0 and f∈R, and consider the two-coordinate model

G(x,θ):=x−aθ,M:=(m00I),Φ:=(f,0),

on O:=R×R2, where x is a centre coordinate and θ an oriented spin. The constraint G=0 is the ideal planar no-slip law of a disk of radius a rolling in a straight line; M is the adopted kinetic matrix (mass m, spin inertia I), and Φ the adopted applied force (force f along x, zero spin torque). Then:

  1. G is a regular holonomic constraint: Gq=(1,−a) has rank one and its kernel is the line R(a,1) of virtual displacements. The associated differential constraint x′−aθ′=0 is exact: on an interval it holds exactly when x−aθ is constant, so in this planar straight-line case the rolling law is holonomic, and after translating x the constant may be taken to be 0;
  2. if q=(x,θ):J→R2 is a C2 ideal-constraint motion of (G,M,Φ) (A regular holonomic constraint, its virtual displacements and the ideal-reaction law), then the d'Alembert condition gives (mx′′−f)a+Iθ′′=0, and differentiating the constraint twice gives x′′=aθ′′; hence

(m+Ia2)x′′=f,x′′=fa2I+ma2,θ′′=afI+ma2;

  1. the multiplier of the motion is λ=−fII+ma2 (uniquely determined), the reaction is R=GqTλ=(λ,−aλ) with ∣λ∣=∣f∣II+ma2≤∣f∣, with strict inequality exactly when f≠0, and the proposition The ideal holonomic system in a possibly time-dependent chart: explicit acceleration, constraint preservation, local existence and reaction power reproduces these numbers from C=GqM−1GqT=1m+a2I and Λ=−C−1GqM−1Φ;
  2. for a uniform thin disk, I=12ma2 and x′′=2f3m; for a hoop (all mass at the rim), I=ma2 and x′′=f2m;
  3. the quadratic motion x(t)=x0+v0t+fa2t22(I+ma2), θ(t)=θ0+ω0t+aft22(I+ma2) with x0=aθ0 and v0=aω0 is an ideal-constraint motion with constant multiplier λ, so the computed acceleration is realised by an explicit solution.

Facts & Assumptions

Given: The constants a,m,I>0 and f∈R, the constraint G(x,θ)=x−aθ, the matrix M=diag(m,I), the force Φ=(f,0), and an open interval J (Intervals of R: the nine order-convex forms, nondegeneracy, and length) carrying a C2 curve q=(x,θ).

[F1]

Constraint data (Ck Euclidean maps and diffeomorphisms, The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case, Submersions and immersions between Euclidean open sets): G is C∞ on R2 with gradient Gq=(1,−a); since 1≠0, Gq has rank 1 and is surjective onto R; and M is symmetric positive definite with inverse diag(m−1,I−1), while Φ is a constant C1 field (Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes).

[F2]

The ideal-constraint model (A regular holonomic constraint, its virtual displacements and the ideal-reaction law): a C2 curve q with G(q(t))=0 and permitted velocity Gq(q(t))q′(t)=0 is an ideal-constraint motion of (G,M,Φ) when there is a continuous λ:J→R with Mq′′=Φ+GqTλ; for every such motion the reaction R:=Mq′′−Φ=GqTλ satisfies the d'Alembert condition R⋅δ=0 for every δ∈ker⁡Gq.

[F3]

Kernel and virtual displacements (Linear map between vector spaces over the same field, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn): δ∈ker⁡Gq means δx−aδθ=0, i.e. δ=(aδθ,δθ)=δθ(a,1); thus the kernel is the line spanned by (a,1) and is one-dimensional.

[F6]

The multiplier lemma (Reactions annihilating the virtual displacements are transpose multiples, with a positive definite multiplier matrix): the matrix C:=GqM−1GqT is symmetric positive definite and invertible (The transpose AT of a matrix, Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn), and a vector R orthogonal to ker⁡Gq is R=GqTλ for a unique λ, namely λ=C−1GqM−1R.

[F7]

The chart proposition (The ideal holonomic system in a possibly time-dependent chart: explicit acceleration, constraint preservation, local existence and reaction power): with the present time-independent data it defines Λ=−C−1[GqM−1Φ+Gtt+2Gtq[v]+Gqq[v,v]] and states that the ideal-constraint motion with permitted velocity satisfies Mq′′=Φ+GqTΛ; here all time derivatives of G vanish.

Verification

technique · direct; compute the kernel, project the ideal equation onto the virtual displacement, differentiate the constraint twice and solve the resulting linear system, then check the multiplier against the proposition and exhibit the motion
1.1F1F3algebra

By [F1], G is C∞ with Gq=(1,−a) of rank one; the virtual displacements are the kernel of [F3], the line R(a,1); the velocity constraint x′−aθ′=0 is exactly Gqq′=0. This is the first assertion of part 1.

1.2F5algebra

Let J⊆R be an interval and q∈C1. If x′−aθ′=0 on J, then the derivative of t↦x(t)−aθ(t) vanishes identically, so the function is constant on J by [F5]; conversely a constant difference has derivative 0. Hence the differential constraint is equivalent to x−aθ≡c, and replacing x by x−c (a translation of the coordinate origin) turns it into G(x,θ)=x−aθ=0. This is the second assertion of part 1.

1.3F2F3algebra

Let q be an ideal-constraint motion. Then Mq′′−Φ=GqTλ and, by [F2], (Mq′′−Φ)⋅δ=0 for every δ∈ker⁡Gq. Taking the virtual displacement δ=(aδ0,δ0), δ0∈R, and writing out Mq′′−Φ=(mx′′−f, Iθ′′), gives δ0[(mx′′−f)a+Iθ′′]=0; choosing δ0=1 yields (mx′′−f)a+Iθ′′=0. This is the first equation of part 2.

1.4F4algebra

Since G(q(t))=x(t)−aθ(t) is identically zero and both x and θ are C2, its second derivative vanishes; by [F4] this second derivative is x′′−aθ′′, so x′′=aθ′′. This is the second relation of part 2.

2.1step 1.3step 1.4algebra

Substituting x′′=aθ′′ of step 1.4 into the identity (mx′′−f)a+Iθ′′=0 of step 1.3 gives (ma2+I)θ′′=af, hence θ′′=afI+ma2 and x′′=aθ′′=fa2I+ma2; substituting θ′′=x′′/a into (ma2+I)θ′′=af and dividing by a gives (m+Ia2)x′′=f. This is the chain of formulas of part 2.

3.1F2F6step 2.1algebra

Multiplier: from Mq′′=Φ+GqTλ in components, mx′′=f+λ and Iθ′′=−aλ; by step 2.1, λ=mx′′−f=mfa2I+ma2−f=−fII+ma2, and then −aλ=afII+ma2=Iθ′′ reproduces the second component. The uniqueness of λ is [F6], applied to R=Mq′′−Φ, which is orthogonal to the kernel by [F2]; thus λ is the unique multiplier of the motion and ∣λ∣=∣f∣I/(I+ma2)≤∣f∣, with equality when f=0 and strict inequality when f≠0, since 0<I/(I+ma2)<1. This is the first assertion of part 3.

3.2step 2.1algebra

Constants: at I=12ma2, x′′=fa212ma2+ma2=fa232ma2=2f3m; at I=ma2, x′′=fa22ma2=f2m. This is part 4.

4.1F1F6F7step 2.1step 3.1algebra

Chart cross-check: by [F1], GqM−1GqT=1m+a2I=I+ma2mI=C and GqM−1Φ=fm, so [F7] gives Λ=−C−1fm=−mII+ma2⋅fm=−fII+ma2, and Mq′′=Φ+GqTΛ gives the same acceleration as step 2.1 and the same multiplier as step 3.1. This is the second assertion of part 3.

4.2F2step 2.1step 3.1algebra

Realisation: let x0,θ0,v0,ω0∈R with x0=aθ0 and v0=aω0, put c:=afI+ma2 and x(t):=x0+v0t+12ac t2,θ(t):=θ0+ω0t+12c t2. Then x(t)−aθ(t)=x0−aθ0+(v0−aω0)t+0=0, so the constraint holds; q is C∞ with x′′=ac, θ′′=c, so x′′=aθ′′; the velocity is permitted because Gt=0 and x′−aθ′=0; and with λ:=−fII+ma2 one has f+λ=mx′′ and −aλ=Iθ′′, so Mq′′=Φ+GqTλ. Hence q is an ideal-constraint motion of (G,M,Φ) with constant multiplier λ, realising the acceleration of part 2. This is part 5.

5.1step 1.1step 1.2step 1.3step 1.4step 2.1step 3.1step 4.1step 3.2step 4.2∎

Parts 1-5 are proved.

Remarks

  • The model is the ideal planar no-slip law with an adopted kinetic matrix; the multiplier λ computed in part 3, i.e. the component of the reaction GqTλ=(λ,−aλ), is the force the model requires along the line of rolling. Whether a real contact can supply that force (sufficient static friction in the contact, no deformation or rolling resistance) is an empirical condition on the apparatus, not a consequence of this mathematical development; when it fails, the no-slip constraint itself is not a valid model and the present equations do not apply.

  • The reduction to a single holonomic constraint uses that the disk rolls in a straight line with fixed vertical plane; the general three-dimensional rolling problem has velocity-dependent constraints that are not integrable in the same way, and no acceleration formula of the present form is asserted for it. The precise planar case is exactly the one identified in Section 5.7 of the cited textbook.

  • The multiplier is the model's constraint force and is not a measured quantity; only the reaction GqTλ=(λ,−aλ) enters the equations, in agreement with the general scope note of A regular holonomic constraint, its virtual displacements and the ideal-reaction law.

Sources