Alphabeta Math
Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck pass
How statement and proof provenance work

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

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

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

Strictly damped oscillator: explicit decay bounds, convergence to rest and asymptotic stability of the origin

Physical scope

The strictly damped oscillator with m>0, k>0, c>0: every motion has explicit position and velocity bounds in normalized SI coordinates, tends to zero together with its velocity, and the origin of the phase plane is a globally asymptotically stable rest point in the phase norm defined by the numerical SI coordinates. The result is a conditional deduction from the declared model data; it asserts no value of the parameters, no measured decay rate and no behaviour of any real damper, and its constants are the explicit ones produced by the classification.

Statement

Adopt the damped oscillator model of The damped oscillator model: a harmonic restoring force with a phenomenological linear drag −cx˙, c≥0 with m>0, k>0 and c>0 (strict damping), force law F(x,v)=−kx−cv and energy E(x,v)=12mv2+12kx2. Position and velocity coordinates in every phase norm below are their normalized numerical SI values, as defined in clause (a) of The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points. For the datum (0,x0,v0), put r0:=∥(x0,v0)∥2 in that norm. In displayed time-dependent formulas, t is read numerically in seconds, and rates in exponents are read in inverse seconds. Write Δ=c2−4mk, γ=c/(2m)>0, λ±=(−c±Δ)/(2m) for Δ≥0, λ=−c/(2m) for Δ=0 and ω=4mk−c2/(2m) for Δ<0. Then:

  1. Motion bounds. For every motion x:R→R of the model with datum (0,x0,v0) and phase norm r0 the following bounds hold for all t≥0, with the coefficients A,B of the classification of Classification of the motions of the damped oscillator by the discriminant, and nonincreasing energy E′=−c x′2:
    • Δ>0: ∣x(t)∣≤(∣A∣+∣B∣)eλ+t and ∣x′(t)∣≤(∣λ+∣∣A∣+∣λ−∣∣B∣)eλ+t;
    • Δ=0: ∣x(t)∣≤(∣A∣+∣B∣t)eλt and ∣x′(t)∣≤(∣λA+B∣+∣λB∣t)eλt;
    • Δ<0: ∣x(t)∣≤(∣A∣+∣B∣)e−γt and ∣x′(t)∣≤(∣ωB−γA∣+∣γB+ωA∣)e−γt; and there are constants C1,C2>0, depending only on (m,k,c), with ∣x(t)∣≤C1r0,∣x′(t)∣≤C2r0(t≥0), for the numerical SI coordinate values; equivalently the corresponding physical position and velocity are bounded by (1 m)C1r0 and (1 m/s)C2r0. Here C1,C2>0 are dimensionless constants in the normalized SI coordinates and depend only on the numerical values of (m,k,c) in the units declared by Dimensions and SI unit conventions of the non-relativistic classical mechanics framework.
  2. Convergence. For every maximal motion, x(t)→0 and x′(t)→0 as t→+∞ (Limits at +∞ and −∞, and infinite limits at a point).
  3. Asymptotic stability. The rest point (0,0) is Lyapunov stable and globally asymptotically stable in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system for the model's phase system: for every ε>0 there is δ>0 such that every maximal motion with ∥(x(t0),x′(t0))∥2<δ for some t0 is defined for all t≥t0 and satisfies ∥(x(t),x′(t))∥2<ε for all t≥t0; and every maximal motion satisfies ∥(x(t),x′(t))∥2→0 as t→+∞. The norm is the Euclidean norm on the normalized numerical SI coordinates of The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points.

The constants are explicit in the proof and no value of m,k,c is asserted.

Facts & Assumptions

Given: The strictly damped model with m>0, k>0, c>0, a motion x:R→R with datum (0,x0,v0), the phase norm r0=∥(x0,v0)∥2 of clause (a) of The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points, and the dimensionless parameters M:=m/(1 kg), K:=k/(1 N/m) and C:=c/(1 kg/s). Numerical rate values are denoted by Λ±:=(1 s)λ±, Λ:=(1 s)λ, Γ:=(1 s)γ and W:=(1 s)ω where applicable.

[L1]

Model vocabulary (The damped oscillator model: a harmonic restoring force with a phenomenological linear drag −cx˙, c≥0, The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria, Classification of the motions of the damped oscillator by the discriminant, and nonincreasing energy E′=−c x′2): all motions extend to R and are classified by the sign of Δ into the three families, with the coefficients A,B of clause 4 of The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data applied to the datum; the energy satisfies E(x(t),x′(t))≤E(x(0),x′(0)) for every t≥0 because E is nonincreasing, and E(0)=E(x0,v0)=12mv02+12kx02; writing X:=x and V:=x′ for the numerical SI coordinates, its numerical energy in joules is E:=E/(1 J)=12MV2+12KX2; the equation is m x′′+c x′+k x=0; the phase system is the one of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data.

[L2]

Strictly damped homogeneous oscillator with k>0: every solution tends to zero: in the three cases, with the coefficients A,B of [L1], ∣x(t)∣≤(∣A∣+∣B∣)eλ+t for Δ>0, ∣x(t)∣≤(∣A∣+∣B∣t)eλt for Δ=0, and ∣x(t)∣≤(∣A∣+∣B∣)e−γt for Δ<0, for t≥0; in each case the decay exponent is negative: λ+<0 for Δ>0 (because Δ<c when k>0), λ=−c/(2m)<0 for Δ=0, and γ>0 for Δ<0.

[L4]

A series of nonnegative terms converges iff its partial sums are bounded, and then the sum is their supremum with The real exponential function and the number e by a power series: for u≥0 the exponential series has nonnegative terms, and the partial sum 1+u is at most the sum, so u≤eu(u≥0).

[L5]

Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}, Squaring is monotone on the nonnegatives: square roots of nonnegative reals exist and are unique, and the square root is strictly increasing on [0,∞); in particular ∣x′(t)∣≤2E(x(t),x′(t))/m for every motion.

Proof

technique · cases on the sign of $\Delta$; use the classification and the decay corollary for the position in normalized SI coordinates, the energy bound for the normalized velocity, and then assemble the uniform bound and stability statement
1.1givenL1L5algebra

Preliminary. Because the phase system is autonomous, all statements are invariant under translating time, and it suffices to consider the datum at t0=0; the functions x and x′ of the statement are then defined for all t≥0 by [L1]. Let X0:=x0 and V0:=v0 be the normalized numerical SI coordinates of the datum, and let M,K be the dimensionless numerical values of m,k defined in Given; then r02=X02+V02. With E:=E/(1 J), the normalized initial energy is E(0)=12MV02+12KX02≤12max⁡(M,K) r02, since X02,V02≤r02. By [L1] the energy is nonincreasing and E(t)≥12MV(t)2, where V(t):=x′(t) is the normalized numerical SI velocity coordinate; hence ∣V(t)∣≤2E(t)M≤max⁡(M,K)M  r0(t≥0). This is the uniform normalized velocity bound with C2E:=max⁡(M,K)/M; and every maximal motion is defined on all of R, while arbitrary local motions are its restrictions by [L1].

2.1assume-case 1givenstep 1.1L1L2L3algebra

Case Δ>0. Let Λ±:=(1 s)λ± and d:=Λ+−Λ−>0 be the numerical SI rates and their difference. Write Aˉ,Bˉ for the numerical SI values of the position coefficients A,B in [L1]. Their formulas give Aˉ=V0−Λ−X0d,Bˉ=Λ+X0−V0d. Hence ∣Aˉ∣≤(1+∣Λ−∣)r0/d and ∣Bˉ∣≤(1+∣Λ+∣)r0/d, since ∣X0∣,∣V0∣≤r0. By [L2], ∣x(t)∣≤(∣Aˉ∣+∣Bˉ∣)eΛ+t, so ∣x(t)∣≤C1r0 with C1:=(2+∣Λ−∣+∣Λ+∣)/d, using eΛ+t≤1 for t≥0 by [L3]. Differentiating x(t)=AˉeΛ+t+BˉeΛ−t gives the numerical velocity x′(t)=Λ+AˉeΛ+t+Λ−BˉeΛ−t, so ∣x′(t)∣≤(∣Λ+∣∣Aˉ∣+∣Λ−∣∣Bˉ∣)eΛ+t, the normalized form of the displayed case bound; it tends to 0 because eΛ+t→0 by [L3].

2.2assume-case 2givenstep 1.1L1L2L3L4algebra

Case Δ=0. Write Λ:=(1 s)λ<0, and let Aˉ,Bˉ be the numerical SI values of the critical-case coefficients in [L1]. The coefficient formulas give Aˉ=X0 and Bˉ=V0−ΛX0, so ∣Aˉ∣≤r0 and ∣Bˉ∣≤(1+∣Λ∣)r0. For t≥0 put u=∣Λ∣t≥0; by [L4], u≤eu, so teΛt≤1/∣Λ∣. By [L2], the normalized position satisfies ∣x(t)∣≤(∣Aˉ∣+∣Bˉ∣t)eΛt≤∣Aˉ∣+∣Bˉ∣/∣Λ∣≤(2+1/∣Λ∣)r0; take C1:=2+1/∣Λ∣. For convergence use te−∣Λ∣t≤(2/∣Λ∣)e−∣Λ∣t/2, obtained from [L4] with u=∣Λ∣t/2. Thus ∣x(t)∣≤∣Aˉ∣e−∣Λ∣t+(2∣Bˉ∣/∣Λ∣)e−∣Λ∣t/2→0 by [L3]. Differentiating x(t)=(Aˉ+Bˉt)eΛt gives the numerical velocity x′(t)=(ΛAˉ+Bˉ+ΛBˉt)eΛt; its t-term tends to 0 by the same inequality and its other term tends to 0 by [L3]. This proves the displayed case bounds and x′(t)→0.

2.3assume-case 3givenstep 1.1L1L2L3algebra

Case Δ<0. Write Γ:=(1 s)γ>0 and W:=(1 s)ω>0, and let Aˉ,Bˉ be the numerical SI values of the oscillatory-case coefficients in [L1]. Then Aˉ=X0 and Bˉ=(V0+ΓX0)/W, so ∣Aˉ∣≤r0 and ∣Bˉ∣≤(1+Γ)r0/W. By [L2], the numerical position satisfies ∣x(t)∣≤(∣Aˉ∣+∣Bˉ∣)e−Γt, which tends to 0 by [L3]; take C1:=1+(1+Γ)/W. Differentiating x(t)=e−Γt(Aˉcos⁡(Wt)+Bˉsin⁡(Wt)) gives the numerical velocity x′(t)=e−Γt[(WBˉ−ΓAˉ)cos⁡(Wt)−(ΓBˉ+WAˉ)sin⁡(Wt)], so its absolute value is bounded by the normalized form of the displayed case bound and tends to 0 by [L3].

3.1cases-exhaustivegivenstep 1.1step 2.1step 2.2step 2.3algebra∎

Uniform bound and stability. By the three cases of steps 2.1-2.3 the numerical position bound ∣x(t)∣≤C1r0 holds for all t≥0 with the constant C1 of the case at hand, and by step 1.1 the numerical velocity bound ∣x′(t)∣≤C2r0 holds with the universal constant C2:=C2E from energy monotonicity; both dimensionless constants depend only on the numerical SI values of (m,k,c). Let ε>0 and set δ:=ε/(C1+C2)>0. If a maximal motion satisfies r0=∥(x(t0),x′(t0))∥2<δ for some t0, then the shifted curve s↦x(t0+s) is again a motion of the autonomous model, with datum (0,x(t0),x′(t0)) and the same initial phase norm r0. Applying the uniform bounds of steps 1.1 and 2.1-2.3 gives ∥(x(t),x′(t))∥2=x(t)2+x′(t)2≤∣x(t)∣+∣x′(t)∣≤(C1+C2)r0<ε(t≥t0), where the norm is on the normalized SI phase coordinates of The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points. Thus (0,0) is Lyapunov stable in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system; and clauses 1 and 2 give x(t)→0, x′(t)→0 for every maximal motion, so every maximal motion tends to the rest point in the phase-plane norm and (0,0) is globally asymptotically stable.

Remarks

  • The three regimes at a glance. For Δ>0 the motion is a combination of two decaying exponentials (overdamped), for Δ=0 of eλt and teλt (critically damped) and for Δ<0 of a decaying oscillation (underdamped); the overdamped case has rate ∣λ+∣, the underdamped case rate γ, and the critical case has a polynomial factor multiplying eλt. In that critical case the proof bounds this factor by an exponential with rate ∣λ∣/2; these rates are positive and the bounds are uniform for data of bounded phase norm.

  • Stability is a statement about the declared model. The asymptotic stability of (0,0) follows from the explicit families of the linear model; it is not a claim that any real oscillator inevitably comes to rest, nor does it carry any empirical qualification beyond the model's own.

Depends on

Dependency tree · 2 levels

Sources