Alphabeta Math
How statement and proof provenance work

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

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

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

One-dimensional conservative and dissipative models

1 · Prerequisites

2 · Summary

This page develops the one-dimensional conservative and dissipative models of the framework as declared restricted model classes and proves their motion theory. The class is fixed by The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria: a single scalar coordinate with an open configuration interval, a fixed positive mass and a scalar force law with the standing continuity and state-Lipschitz regularity, together with the vocabulary of motions, the phase system, the potential presentation with residual forces, the mechanical energy, the equilibria of an autonomous force law and the oscillator-type data (m,c,k,f), with their dimensions. The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points supplies the geometric reading: the phase plane and phase curves, the energy function with its level and sublevel sets, turning positions and points, and the equilibrium and stability vocabulary carried over from the mathematical definitions.

The conservative theory is Existence, uniqueness and conservation of energy for a one-dimensional conservative model: phase curves lie in energy level sets and the velocity is determined by the position and energy up to sign: for a time- and velocity-independent force law F=−U′ with a C2 potential, every datum has exactly one maximal motion, the mechanical energy 12mx′2+U(x) is constant along every motion, the position stays in the sublevel set {U≤E0}, the velocity is determined by the position up to sign wherever U<E0, and the constant motions are exactly the rest motions at equilibria. At a strict local minimum of U, A strict local minimum of the potential is a Lyapunov stable equilibrium of the conservative one-dimensional model, and its linearisation is an oscillator-type model gives Lyapunov stability of the rest point and identifies the linearised model as the oscillator-type model with k=U′′(x∗), with no instability claim in any other case.

The three linear models are adopted postulates with their declared regimes and empirical qualifications: The harmonic (linear restoring) model: a declared linear force law with parameter k>0 (force −kx, potential kx2/2, the quadratic approximation at a stable minimum), The damped oscillator model: a harmonic restoring force with a phenomenological linear drag −cx˙, c≥0 (the harmonic restoring force with the phenomenological drag −cv, fluctuations neglected in its stated regime) and The sinusoidally driven oscillator model: a prescribed external force f0cos⁡(Ωt) added to the damped oscillator (the damped oscillator with the prescribed external force f0cos⁡(Ωt)). Their motion theory is exact and conditional: The motions of the harmonic model are the sinusoids Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)): amplitude and phase, the least period, the conserved energy kR2/2 and the elliptical phase curves classifies the motions as the sinusoids Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)) with amplitude and phase, least period 2π/ω, conserved energy 12kR2 and the elliptical phase curves of the energy level sets; Classification of the motions of the damped oscillator by the discriminant, and nonincreasing energy E′=−c x′2 classifies the damped motions by the sign of c2−4mk and proves the dissipation identity E′=−c x′2 with the characterisation of the rest motion; Strictly damped oscillator: explicit decay bounds, convergence to rest and asymptotic stability of the origin bounds every strictly damped motion by its initial phase norm, proves the convergence to rest together with the velocity and the global asymptotic stability of the origin; and The driven oscillator: steady response and amplitude profile, bounded damped and off-resonant motions, and unbounded undamped resonance gives the steady sinusoidal response with its amplitude ∣f0∣/(k−mΩ2)2+c2Ω2 and, when f0≠0 and c2<2mk, its unique maximum at a positive drive frequency, the reduction of every motion to the steady response plus a homogeneous motion, boundedness for c>0 with exponential approach to the steady response, boundedness at c=0 away from Ω=k/m, and the unbounded resonant motions at c=0, Ω=k/m with nonzero drive. The closing audit Scope and limits of the one-dimensional conservative and dissipative models: declared class, approximation status, empirical parameters and open extensions separates the declared assumptions, the conditional deductions, the approximation status of the quadratic models and the empirical handling of measured parameters, and records the open extensions. Worked checks of the models on explicit data are collected on cm04-one-dimensional-conservative-and-dissipative-models-examples.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria

Physical scope

The declared one-dimensional model class of the framework: a single scalar coordinate with an open configuration interval, a fixed positive mass, and a scalar force law with the standing continuity and state-Lipschitz regularity; the motions, the phase system, the potential presentation with a scalar potential and residual forces, the mechanical energy, the equilibria of a time- and velocity-independent force law, and the oscillator-type data (m,c,k,f). The item fixes vocabulary, domains, regularity and dimensions for the conservative and dissipative models of this page; it adopts no force law of nature, asserts no existence, uniqueness, conservation or stability statement, and records no measured value.

Definition

(a) Data of a one-dimensional model. A one-dimensional model is declared by the following data: a fixed positive mass m>0, of dimension μ read in kilograms (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework, Particle labels, worldlines, the frame-time parameter condition and fixed inertial masses clause (v)); a time window I0⊆R, a nondegenerate open interval read in seconds (Intervals of R: the nine order-convex forms, nondegeneracy, and length); an open configuration domain Q⊆R, a nonempty open interval of the real line read in metres; and a force law F:I0×Q×R⟶R,(t,x,v)⟼F(t,x,v), of dimension μℓτ−2, read in newtons, subject to the standing regularity of the framework: F is continuous on I0×Q×R and locally Lipschitz in the state variable (x,v) locally uniformly in time (Local Lipschitz continuity in the state variable, locally uniform in time and parameters). The framework's three-dimensional force branch is Adopted force model: force maps, the equation of motion, additive composition of force contributions and frame covariance; the class declared here is a restricted model class adopted at the same primitive level, in which the physical position of a single label along one coordinate is the whole state. A model whose force law fails the standing regularity is not a one-dimensional model of this class. When a statement differentiates the force law or the potential, the additional smoothness it needs is stated in that statement.

(b) Motions. Let J⊆I0 be an interval with nonempty interior (endpoints read with one-sided derivatives, as in First-order systems, initial value problems, and solutions on intervals). A motion of the model on J is a curve x:J→Q of class C2 (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral) with m x′′(t)=F(t,x(t),x′(t))(t∈J). In the notation of the three-dimensional branch this is the scalar restriction of its equation of motion clause (iv) to a single coordinate; the motions of the present class are declared to be these curves. A datum at time t0∈I0 is a triple (t0,x0,v0) with x0∈Q and v0∈R; a motion x takes the datum when x(t0)=x0 and x′(t0)=v0.

(c) Phase system. The phase system of the model is the first-order system on the open set I0×Q×R⊆R3 with right-hand side GF(t,x,v):=(v, F(t,x,v)m), read as a system of First-order systems, initial value problems, and solutions on intervals; its state variable is z=(x,v). For an interval J⊆I0 with nonempty interior, a phase solution on J means a differentiable curve z:J→Q×R satisfying z′(t)=GF(t,z(t)) for every t∈J, with domain-relative derivatives at included endpoints. Equivalently, for any t0∈J it is a solution in the cited IVP sense with datum (t0,z(t0)). For a C2 curve x:J→Q, being a motion is equivalent to t↦(x(t),x′(t)) being a phase solution: the first component is x′=v and the second is x′′=F(t,x,x′)/m. Conversely a phase solution (x,v) has x′=v and v′=F(t,x,v)/m; continuity of F makes v′ continuous, so x is C2 and is a motion. This is the argument of clause (c) of The associated first-order system, motions, maximal motions and the Newtonian flow of a force model read in one coordinate. Here x is the position or displacement coordinate and v the velocity coordinate; both are read in the fixed inertial frame in which the model is declared. The phase system inherits the standing regularity: GF is continuous and locally Lipschitz in the state variable locally uniformly in time, so the Picard–Lindelöf hypotheses of Picard-Lindelöf local existence and uniqueness for first-order systems hold at every point of its domain.

(d) Potential presentation, residual forces and energy. The model presents its force law through a potential when it declares a function U:Q→R of class C2 in the sense of clause (a) of Potential presentations of a force model, potential and mechanical energy, closed versus exact force fields, local potentials and equipotential level sets of dimension μℓ2τ−2, read in joules, and a residual force d:I0×Q×R→R of the standing regularity of (a), with F(t,x,v)=−U′(x)+d(t,x,v)((t,x,v)∈I0×Q×R); here U′ is the ordinary derivative. This is the one-coordinate case of the potential presentation of Potential presentations of a force model, potential and mechanical energy, closed versus exact force fields, local potentials and equipotential level sets clause (a) for a time- and velocity-independent potential, and the scalar U is a potential energy of the model's potential part. The model is a global potential model (equivalently, in the framework's usage, a conservative model) when d≡0; the term "conservative" is used only in this sense, as recorded in that definition. For a conservative model the mechanical energy is the function E:Q×R⟶R,E(x,v):=12mv2+U(x), of dimension μℓ2τ−2. For a model with declared residual forces, the mechanical energy along a motion x:J→Q is E(t):=12m x′(t)2+U(x(t)); it is defined along curve families even though the force law is not presented through U alone.

(e) Equilibria. A time- and velocity-independent force law on Q is one of the form F(t,x,v)=f(x) for a function f:Q→R; the model is then autonomous. For such a law, a point x∗∈Q is an equilibrium (or rest position) when f(x∗)=0, and the phase-plane point (x∗,0) is the corresponding rest point of the phase system; the constant curve x≡x∗ is then a motion. When the force law is presented as f=−U′ with U of class C2, the equilibria are exactly the critical points of U: f(x∗)=0  ⟺  U′(x∗)=0. A point x0∈Q is a strict local minimum of U when U(x)>U(x0) for every x∈Q with 0<∣x−x0∣<r0 for some r0>0. The Lyapunov stability and asymptotic stability of a rest point of the autonomous phase system are the notions of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system applied to the field GF; they are not defined here anew.

(f) Oscillator-type data. Let m>0, c≥0, k∈R and let f:I0→R be continuous. A model is of oscillator type with data (m,c,k,f) when its force law has the form F(t,x,v)=f(t)−c v−k x((t,x,v)∈I0×Q×R) with Q a nonempty open interval containing 0 when the equation is read about the origin, and the parameters are read with dimensions [c]=μτ−1 (kilograms per second), [k]=μτ−2 (kilograms per second squared) and [f]=μℓτ−2 (newtons). For such a model the equation of motion is the second-order equation m x′′+c x′+k x=f(t), and the motions are exactly the restrictions of the solutions of the mathematical model of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data with the same data, by the identification of that definition's phase system with the phase system (c). The parameters c and k are model data: they enter the model's declaration and are not derived by this item.

(g) Status. This item introduces vocabulary, objects, regularity and dimensions; it contains no existence, uniqueness, conservation, stability or boundedness assertion, no force law of nature and no measured value. In particular, clause (a) declares the one-dimensional class rather than deriving it from the three-dimensional branch, and clauses (d)-(f) only fix which data a model of the class carries. Every later statement of this page names the model data, the potential representative and the initial datum it uses.

Remarks

DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points

Physical scope

The geometric vocabulary of the one-dimensional models of this page: the phase plane with position and velocity expressed as numerical SI coordinates, the Euclidean norm on that normalized pair, the phase curve of a motion, the mechanical energy function of a conservative model, its energy level sets and sublevel sets, and turning positions and turning points. The item also fixes how the equilibrium, Lyapunov-stability and linear-stability vocabulary of the mathematical definitions is read for the scalar phase system. It asserts no conservation, no containment of phase curves in level sets, no boundedness and no measured value; those are statements of the theorems of this page.

Definition

Adopt a one-dimensional model (m,I0,Q,F) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria, with phase system GF(t,x,v)=(v,F(t,x,v)/m).

(a) Phase plane and phase curves. The phase plane of the model is R2 with coordinates (x,v), where x is the position and v the velocity coordinate of clause (b) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria. They are real numerical SI coordinates: if Xphys and Vphys are the corresponding physical position and velocity, then x=Xphys/(1 m) and v=Vphys/(1 m/s) (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework). The phase domain of the model is the open set I0×Q×R of (time, position, velocity) triples. The phase plane is read with the Euclidean structure and norm ∥(x,v)∥2=x2+v2 on these dimensionless numerical coordinates (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn); physical displacement and velocity are not added as dimensionful quantities. Distances used below are those of this norm. For a curve family x:J→Q of class C1 its phase curve is t↦(x(t),x′(t)); when x is a motion of the model the phase curve is the corresponding solution of the phase system.

(b) The energy function of a conservative model. Let the model be a global potential model of clause (d) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with potential U of class C2, so that F=−U′. The energy function of the model is E:Q×R⟶R,E(x,v):=12mv2+U(x), of dimension μℓ2τ−2 (joules). For E0∈R the energy level set of E0 and the energy sublevel set of E0 are ΣE0:={(x,v)∈Q×R:E(x,v)=E0},SE0:={(x,v)∈Q×R:E(x,v)≤E0}.

(c) Turning positions and turning points. Let the model be as in (b). A point x∈Q is a turning position of the energy E0 when U(x)=E0, and the phase-plane point (x,0) is the corresponding turning point of E0. A turning position at which U′ does not vanish is a simple turning position. In the one-coordinate reading of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart the potential difference U(x)−E0 is the sign-controlling quantity for the velocity of a motion of energy E0; the statements of that kind are the content of the conservation theorem of this page and are not made here.

(d) Rest points and stability vocabulary. For a model with a time- and velocity-independent force law f, the equilibria, rest points, strict local minima and the Lyapunov and asymptotic stability of a rest point are those of clause (e) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria, read for the phase system GF; Lyapunov stability is the notion of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system and linear stability is that of Invariant sets, orbital stability of an invariant set, and linear stability of an equilibrium, applied to GF at its rest point. For an oscillator-type model with data (m,c,k,f) of clause (f) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria the phase system is the phase-plane system of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data, and the phase-plane reading of its equilibria is the one recorded there.

(e) Status. This item fixes geometric objects and reads mathematical stability vocabulary for the scalar phase system; it asserts no conservation of energy, no containment of a phase curve in a level set, no boundedness of level sets and no measured value.

Remarks

  • The phase plane is two-dimensional, and the model is not the pendulum. The coordinates (x,v) are the position and velocity of the declared scalar model; a model with an angle coordinate would carry a quotient of the line by the periodicity of the physical configuration, which is not part of this definition. The pendulum model of the companion page declares its own coordinate interval.

  • Level sets are not orbits. ΣE0 is a subset of the phase plane defined by an equation; that the phase curve of a motion with initial energy E0 lies in it is a conservation statement, proved separately. A level set may be empty, a single rest point, a closed curve or a union of curves, and this item asserts nothing about which case occurs.

PostulateStatement: AI-adaptedProof: Not applicableOpen item page →

The harmonic (linear restoring) model: a declared linear force law with parameter k>0

Physical scope

Adoption of the harmonic oscillator as a restricted one-dimensional model: configuration domain the whole line, fixed mass m>0, force law F(x)=-k x with declared spring constant k>0 of dimension mass per time squared, and potential U(x)=k x^2/2. The model is a linear idealization of a restoring force, adopted for a stated regime (small displacements from a stable equilibrium, in which the quadratic term of the potential is the adequate approximation, as in the sources' linearization procedure); the spring constant is model data whose value, if taken from measurement, is an empirical premise with its own conditions and uncertainty. The item asserts no value of k, no claim that the model applies to a given spring or oscillation, and no measured result.

Postulate

(i) Setting. Work in the one-dimensional model class of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with time window I0=R, configuration domain Q=R, a fixed mass m>0 and a spring constant k>0 of dimension μτ−2, read in newtons per metre (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework). The harmonic model is the model whose force law is the linear restoring force F(x,v):=−k x. Its equation of motion is m x′′=−k x, equivalently m x′′+k x=0; the model is autonomous and velocity-independent, and its potential is U(x):=12k x2, of dimension μℓ2τ−2 (joules), which presents the force law in the sense of clause (d) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria: U is of class C2 in that clause's sense, and −U′(x)=−k x=F(x) (For a natural n≥1 the function x↦xn is differentiable everywhere with derivative ι(n) x n−1; for n=0 it is the constant 1, with derivative 0; for a natural n≥1 the function x↦x−n is differentiable at every x≠0 with derivative −ι(n) x−n−1; consequently every polynomial function is differentiable at every real, with the derivative computed term by term, The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set). The harmonic model is therefore a conservative model of the class, with the potential representative U declared above and no residual force.

(ii) Regularity and framework status. Let X:=xphys/(1 m) and V:=vphys/(1 m/s) be the dimensionless numerical SI state coordinates (the real state coordinates in the displayed force-map signature), and let K:=kphys/(1 N/m) be the numerical spring constant. In this representation the numerical force map F^:=Fphys/(1 N) is F^(X,V)=−KX, so ∣F^(X,V)−F^(X′,V′)∣≤K ∥(X,V)−(X′,V′)∥2. It is therefore locally Lipschitz in the state variable with global numerical-coordinate constant K, uniformly in time (Local Lipschitz continuity in the state variable, locally uniform in time and parameters). The harmonic model is thus an admissible model of the class of clause (a) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria, and the existence, uniqueness and conservation statements of this page apply to it. It is adopted as a restricted model class of the framework, at the same primitive level as the force maps of Adopted force model: force maps, the equation of motion, additive composition of force contributions and frame covariance: the force law is model input.

(iii) Physical reading and regime. The model idealizes a restoring force linear in the displacement: an ideal spring, or the quadratic approximation of a potential U near a strict local minimum x∗ with U′′(x∗)>0, where the Taylor expansion of the effective potential begins with the quadratic term 12U′′(x∗)(x−x∗)2 and the linearized equation carries the effective spring constant U′′(x∗) (the sources' linearization procedure). The model is adopted for a stated regime — displacements so small that the neglected remainder of the potential is adequate for the purpose at hand — and it is not claimed to describe large-amplitude oscillations, nonlinear springs, plasticity, or any given real system. The coefficient k is a model parameter, not a constant of the framework; when its value is taken from a calibration experiment, that value is an empirical premise carrying the apparatus conditions and the measurement uncertainty, and the framework asserts none of it here. The linear form is the simplest restoring law satisfying the standing regularity; that is a modelling convenience, not a derivation.

(iv) Status. Adopted assumption of the framework when a model invokes it; it selects no law of nature on its own. Every consequence below is a conditional statement about the declared equation of motion.

Remarks

  • Zero, negative and nonlinear alternatives are different models. For k=0 the model is the free model m x′′=0; for k<0 the force is repulsive near the origin; a force law with a cubic term is a different model. For the ideal planar point pendulum, declare a bob mass M>0 (kilograms), length l>0 (metres), fixed pivot, and constant downward gravitational acceleration g>0 (metres per second squared), with no friction. Its angular coordinate is a C2 curve θ:J→R on an open time interval J read in seconds; θ is the angle from the downward vertical, read in dimensionless radians. The angular model adopts inertia I:=Ml2 (kilogram metres squared) and the scalar torque map T:R→R, T(θ):=−Mglsin⁡θ (newton metres), with equation Iθ′′=T(θ), equivalently θ′′=−(g/l)sin⁡θ. Here sine is the everywhere-convergent real power series of Sine and cosine defined by their real power series. This is the recorded source's point-pendulum equation (Likharev, Section 1.2, Eq. (1.18), multiplied by l), adopted as an angular model, rather than a translational force law on the length coordinate x. The companion page's small-angle check quantifies its difference from the linearized angular equation; the linear theory does not describe its exact motions.

  • The potential is declared, not measured. Only the value of k enters the physics; the additive constant of U is fixed by the convention U(0)=0 recorded above, which is part of the model's declaration, exactly as permitted by Potential presentations of a force model, potential and mechanical energy, closed versus exact force fields, local potentials and equipotential level sets clause (c); that clause requires a declared representative, not this particular zero.

Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

Existence, uniqueness and conservation of energy for a one-dimensional conservative model: phase curves lie in energy level sets and the velocity is determined by the position and energy up to sign

Physical scope

The basic dynamical content of a conservative one-dimensional model: every datum has exactly one maximal motion; the mechanical energy m x-dot^2/2 + U(x) is constant along every motion; a motion of initial energy E_0 keeps its position in the sublevel set {U <= E_0}, and where U is strictly below E_0 its velocity is one of the two square-root values determined by the position and E_0; every non-turning phase point of the motion lies on the energy level set of E_0; and the constant motions are exactly the rest motions at equilibria. The result is a conditional deduction from the declared potential presentation; it asserts no periodicity, no boundedness of the motion, no asymptotic behaviour and no measured value.

Statement

Adopt a conservative one-dimensional model: the class of clause (a) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with a time- and velocity-independent force law presented through a potential, that is F(x)=−U′(x)(x∈Q) for a function U:Q→R of class C2 (so the residual force of clause (d) of that definition vanishes identically and the model is a global potential model). Write E(x,v)=12mv2+U(x) for the energy function of The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points clause (b). Then:

  1. Existence and uniqueness. Every datum (t0,x0,v0) with t0∈I0, x0∈Q, v0∈R has exactly one maximal motion: a motion x:Jmax⁡→Q with x(t0)=x0 and x′(t0)=v0 that has no extension to a motion on a strictly larger interval. Its interval Jmax⁡ is open, contains t0 and is contained in I0, and every motion of the model with the same datum is the restriction of x to a subinterval.
  2. Energy conservation. For every motion x:J→Q the function t↦E(x(t),x′(t)) is constant on J.
  3. Level-set containment and velocity determination. Let x:J→Q be a motion with datum (t0,x0,v0) and let E0:=E(x0,v0). Then U(x(t))≤E0 for every t∈J; if U(x(t))<E0 then x′(t)=±2(E0−U(x(t)))m, the two signs giving the only two possible values; and x′(t)=0 holds exactly when U(x(t))=E0. In particular the phase curve of x is contained in the energy level set ΣE0 of clause (b) of The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points.
  4. Rest motions. A point x∗∈Q is an equilibrium of the model if and only if U′(x∗)=0; in that case the constant curve x≡x∗ is the maximal motion with datum (t0,x∗,0), and conversely the maximal motion with datum (t0,x∗,0) is constant exactly when U′(x∗)=0.

The statement is conditional on the declared potential U: replacing U by U+c with a constant c changes E on the whole phase plane by c and leaves clauses 1, 3 and 4 unchanged. No claim of periodicity, boundedness of the position or of the velocity, or approach to an equilibrium is made.

Facts & Assumptions

Given: The conservative one-dimensional model (m,I0,Q,F) with F=−U′, its phase system GF(x,v)=(v,−U′(x)/m), its energy E, and a datum (t0,x0,v0).

[L1]

Model vocabulary (The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria clauses (a)-(e)): Q is nonempty open, F is continuous and locally Lipschitz in the state variable; motions are the C2 curves x:J→Q with mx′′=F(x) on an interval J⊆I0 with nonempty interior; x is a motion if and only if t↦(x(t),x′(t)) is a solution of z′=GF(z); and x∗ is an equilibrium if and only if F(x∗)=0, equivalently U′(x∗)=0 because U is C1 and F=−U′.

[L2]

Phase-field regularity: GF(x,v)=(v,−U′(x)/m) is continuous on the open set Q×R, and it is locally Lipschitz in the state variable: the first component is linear, and U′ is of class C1 because U is C2, hence locally Lipschitz on the compact coordinate subsets of the open set Q by A C1 map is locally Lipschitz on compact coordinate subsets with the identity charts; a locally Lipschitz map is locally Lipschitz in the sense of Local Lipschitz continuity in the state variable, locally uniform in time and parameters.

[L3]

Maximal solutions (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, First-order systems, initial value problems, and solutions on intervals): for a continuous field on an open domain that is locally Lipschitz in the state variable, every initial datum has exactly one maximal solution; its maximal interval is open; and every solution with the same datum is a restriction of it.

[L6]

Square roots: for a≥0 all solutions of y2=a are y=±a, and those two values are distinct exactly when a>0 (Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}).

Proof

Proof technique: direct; apply the maximal-solution theorem to the phase system, differentiate the energy along a motion, and read the consequences off the conservation identity.

1.1givenL1L2L3

The phase field GF satisfies the hypotheses of [L3] on the open time-state domain I0×Q×R by [L2]. Hence every initial state (x0,v0) has exactly one maximal solution z:Jmax⁡→Q×R of z′=GF(z) with z(t0)=(x0,v0), with Jmax⁡ open and containing t0, and every solution of the phase system with the same initial state is a restriction of z; the maximal interval is contained in the time projection I0 of this domain. By the biconditional of [L1] the first component x of z is a motion and the second is x′; and if y is any motion with the same datum, then t↦(y(t),y′(t)) is a solution of the phase system with the same initial state and hence a restriction of z, so y is a restriction of x. This is clause 1.

1.2givenL1L4L5algebra

Let x:J→Q be a motion and put e(t):=E(x(t),x′(t))=12m x′(t)2+U(x(t)) for t∈J. At every interior point t of J, x′ is differentiable with derivative x′′ and x is differentiable, so by [L4] the function e is differentiable there with e′(t)=m x′(t)x′′(t)+U′(x(t))x′(t)=x′(t)(m x′′(t)+U′(x(t)))=x′(t)(F(x(t))+U′(x(t)))=0, where the equation of motion m x′′(t)=F(x(t)) and F=−U′ were used. The curve e is continuous on the interval J, so [L5] gives that e is constant on J. This is clause 2.

2.1givenstep 1.2L6algebra

Let x:J→Q be a motion with datum (t0,x0,v0) and E0=E(x0,v0). By step 1.2, E(x(t),x′(t))=E0 for every t∈J; subtracting U(x(t)) gives 12m x′(t)2=E0−U(x(t)). The left side is nonnegative, so U(x(t))≤E0; and x′(t)2=2(E0−U(x(t)))/m. If U(x(t))<E0 then the right side is strictly positive, and [L6] applied to a=2(E0−U(x(t)))/m>0 gives x′(t)=±a; while x′(t)=0 holds exactly when U(x(t))=E0. Finally (x(t),x′(t))∈ΣE0 for every t∈J because E(x(t),x′(t))=E0, by the definition of the level set. This is clause 3.

3.1givenstep 1.1L1∎

Clause 4: by [L1], x∗ is an equilibrium exactly when F(x∗)=0, that is U′(x∗)=0. If so, the constant curve x≡x∗ has x′′≡0 and m⋅0=F(x∗)=0, so it is a motion with datum (t0,x∗,0); by the uniqueness clause of step 1.1 it is the maximal motion with that datum. Conversely, if the maximal motion with datum (t0,x∗,0) is the constant curve x≡x∗, then substituting into the equation of motion gives 0=m⋅0=F(x∗), so U′(x∗)=0. This is clause 4.

Remarks

  • What conservation does and does not give. Clause 3 confines the motion to the sublevel set {U≤E0} and fixes the speed at each position where the level is crossed, up to sign; it does not by itself show that the motion oscillates between turning positions, that a period exists, or that the motion stays bounded for all forward times. Those require additional hypotheses on U and are not claimed here.

  • The sample-independent sign. On an interval on which the velocity does not vanish, its sign is constant, because a continuous function with no zero has constant sign; this is why the phase curve travels along one branch of the level set between consecutive turning points. No turning points are asserted to exist for a general U.

  • Normalization. Adding a constant to U shifts E by that constant on the whole phase plane; ΣE0 moves to ΣE0+c and the velocity formula of clause 3 is unchanged, as it must be, since the force is unchanged. The energy conservation statement is a statement about the declared representative U.

Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

A strict local minimum of the potential is a Lyapunov stable equilibrium of the conservative one-dimensional model, and its linearisation is an oscillator-type model

Physical scope

The equilibrium theory of a conservative one-dimensional model near a strict local minimum of its potential: the rest point is Lyapunov stable, with the Euclidean phase-plane norm, and its linearisation is the oscillator-type model with data (m,0,U''(x_),0), which is the harmonic model with spring constant U''(x_) when that number is positive. The statement is conditional on the declared potential and its strict local minimum; it asserts no asymptotic stability, no attraction and no instability in any case, and no measured value.

Statement

Adopt a conservative one-dimensional model as in Existence, uniqueness and conservation of energy for a one-dimensional conservative model: phase curves lie in energy level sets and the velocity is determined by the position and energy up to sign: the class of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with force law F(x)=−U′(x) for a potential U:Q→R of class C2, phase system GF(x,v)=(v,−U′(x)/m) on Q×R, and energy E=12mv2+U. Let x∗∈Q be a strict local minimum of U: there is r0>0 with B‾(x∗,r0)⊆Q and U(x)>U(x∗) for 0<∣x−x∗∣<r0. Then:

  1. Lyapunov stability. (x∗,0) is a Lyapunov stable rest point of the phase system in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system: for every ε>0 with B‾((x∗,0),ε)⊆Q×R there is δ>0 such that the unique maximal solution of the autonomous phase equation through the datum of every motion x of the model with ∥(x(t0),x′(t0))−(x∗,0)∥2<δ for some t0∈I0 is defined for all t≥t0 and satisfies ∥(x(t)−x∗, x′(t))∥2<ε(t≥t0). Here x in the displayed all-time bound denotes that maximal extension; on the original motion interval it agrees with the given motion. The autonomous equation is read on R×Q×R, even when the declared time window I0 is smaller; no enlargement of that declared window is asserted. The norm is the Euclidean norm on R2 of The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn; the balls are those of Open ball, closed ball and sphere in a metric space. No asymptotic stability and no attraction are asserted.

  2. Linearisation. U′′(x∗) exists (Ck Euclidean maps and diffeomorphisms), the phase field GF is totally differentiable at (x∗,0) with DGF(x∗,0)(δx,δv)=(δv, −U′′(x∗)m δx), and the linearised equation is m ξ′′+U′′(x∗) ξ=0, which is the equation of the oscillator-type model of clause (f) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with data (m,0,U′′(x∗),0) (The damped driven oscillator model: second-order equation, phase-plane system and characteristic data). If in addition U′′(x∗)>0, then that linearised model is the harmonic model of The harmonic (linear restoring) model: a declared linear force law with parameter k>0 with spring constant U′′(x∗), and its undamped frequency in the sense of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data is U′′(x∗)/m.

The statement is conditional on the potential U and on the strict local minimum; no conclusion is drawn for a critical point that is only a weak minimum, a maximum, or a point where U′′ vanishes identically on one side.

Facts & Assumptions

Given: The conservative one-dimensional model with potential U, the strict local minimum x∗, and the quantity D2U:=U′′(x∗).

[L1]

The phase field GF(x,v)=(v,−U′(x)/m) is exactly the mechanical vector field of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart with n=1, Ω:=Q, the constant kinetic metric field M(x)≡m (symmetric positive definite, of class C2 because it is constant), and the potential U: for a constant field M, DqM=0 and ∇q(12vTMv)=0, so the second component of G is M(x)−1(0−0−U′(x))=−U′(x)/m. Its energy E(x,v)=12vTMv+U(x) is the energy of The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points clause (b). [given, The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart]

[L2]

A strict local minimum of the potential gives a Lyapunov stable equilibrium of the mechanical system, by energy conservation: for n≥1, an open Ω⊆Rn, a C2 field M of symmetric positive definite matrices, a C2 potential U and a strict local minimum q0 of U, the point (q0,0) is a Lyapunov stable equilibrium of the mechanical system of [L1]'s definition: for every ε>0 with B‾((q0,0),ε)⊆Ω×Rn there is δ>0 such that every maximal solution (q,v) with ∥(q(0),v(0))−(q0,0)∥2<δ is defined for all t≥0 and satisfies ∥(q(t),v(t))−(q0,0)∥2<ε for every t≥0. The statement uses no asymptotic stability.

[L3]

Linearisation of the mechanical system at an equilibrium: M(q_0) xi'' + D²U(q_0) xi = 0, and the small-oscillation quadratic model: at a mechanical equilibrium q0 of the data of [L1]'s definition, the mechanical vector field is totally differentiable at (q0,0) with DG(q0,0)(δq,δv)=(δv,−M(q0)−1D2U(q0)δq); the linearised equation is M(q0)ξ′′+D2U(q0)ξ=0.

[L4]

Model definitions (The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria clauses (b), (e), (f), The damped driven oscillator model: second-order equation, phase-plane system and characteristic data, The harmonic (linear restoring) model: a declared linear force law with parameter k>0): motions of the model are exactly the curves whose phase curves solve z′=GF(z) on intervals J⊆I0; x∗ is an equilibrium exactly when U′(x∗)=0, which holds at a strict local minimum; the oscillator-type data (m,0,k,0) with k=U′′(x∗) have equation mξ′′+kξ=0; and the harmonic model of the postulate has force law −kx and potential 12kx2.

[L5]

Since U is C2, the unrestricted autonomous phase field (t,x,v)↦(v,−U′(x)/m) is C1 on R×Q×R. It is continuous and locally Lipschitz in the state, locally uniformly in time, by A C1 map is locally Lipschitz on compact coordinate subsets. Thus Every Picard–Lindelöf initial value problem has one maximal solution on an open interval gives a unique maximal solution through each datum; every solution through that datum is its restriction.

Proof

Proof technique: direct; identify the scalar phase system with the one-dimensional instance of the mathematical mechanical system, apply the energy-stability theorem to it, and read the linearisation from the linearisation proposition.

1.1givenL1L4

By [L1] the phase field GF of the model is the mechanical vector field of the mathematical data (n,Ω,M,U)=(1,Q,m,U) with the constant positive definite metric field M≡m. On intervals J⊆I0, the phase curves of model motions are exactly the solutions (q,v) of that mechanical system; unrestricted solutions are model motions only on their restrictions to the declared time window. Its energy is E. At the strict minimum, the difference quotient (U(x∗+h)−U(x∗))/h is nonnegative for small h>0 and nonpositive for small h<0. Its common two-sided limit U′(x∗) is therefore zero, as required by [L2] and [L3].

2.1givenstep 1.1L2L4L5

Clause 1. Let ε>0 with B‾((x∗,0),ε)⊆Q×R. Apply [L2] with Ω=Q, M≡m, the potential U and q0=x∗: it produces δ>0 such that every maximal solution (q,v) of the mechanical system with ∥(q(0),v(0))−(x∗,0)∥2<δ is defined for all t≥0 and satisfies ∥(q(t),v(t))−(x∗,0)∥2<ε for all t≥0. Now let x:J→Q be a motion of the model and t0∈J with ∥(x(t0),x′(t0))−(x∗,0)∥2<δ. The phase system is autonomous, so the translate s↦(x(t0+s),x′(t0+s)), defined for s with t0+s∈J, is a solution of the same system with initial value (x(t0),x′(t0)) at s=0; it has a unique maximal extension for this unrestricted autonomous equation by [L5], and by step 1.1 and the defining property of δ that extension is defined for all s≥0 and stays at distance <ε from (x∗,0). Translating the maximal extension back yields a maximal curve whose domain contains [t0,∞) and which agrees with x on J, and ∥(x(t)−x∗,x′(t))∥2<ε for all t≥t0. Since ε was arbitrary, (x∗,0) is Lyapunov stable in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system.

3.1givenstep 1.1L3L4∎

Clause 2, differentiability and linearised equation. U is C2, so D2U(x∗)=U′′(x∗) exists and x∗ is a mechanical equilibrium because U′(x∗)=0 at a strict local minimum. Applying [L3] with q0=x∗, M(q0)=m and D2U(q0)=U′′(x∗) gives the displayed total derivative and the linearised equation mξ′′+U′′(x∗)ξ=0. By [L4] this is exactly the oscillator-type model with data (m,0,U′′(x∗),0). If U′′(x∗)>0, then k:=U′′(x∗)>0 and [L4] identifies the equation with that of the harmonic model with spring constant k; the frequency k/m is the undamped frequency ω0 attached to the data (m,0,k,0) by The damped driven oscillator model: second-order equation, phase-plane system and characteristic data, and it is quoted here only as the linearised frequency, not as a property of the nonlinear motion.

Remarks

  • Stability is not attraction. Clause 1 asserts only Lyapunov stability: nearby motions stay nearby, with no claim that they approach the rest point. For the conservative model the energy is conserved, and the companion results do not make the rest point attracting; the dissipative model is where attraction is proved.

  • The weak-minimum and saddle cases are not decided here. If a critical point with vanishing U′′ is still a strict minimum, clause 1 applies; at other degenerate critical points or a strict local maximum this item says nothing about stability or instability of the rest point; the linearised model with U′′(x∗)<0 is the repulsive oscillator, whose explicit growing solution is recorded in the companion pages of the mathematics track, but transferring that growth to the nonlinear model requires a separate nonlinear argument that this page does not supply.

  • Why the identification is legitimate. The mechanical system of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart is a mathematical chart model with an arbitrary positive definite C2 metric field; the model of this page is its one-dimensional constant-metric instance, and the identification in step 1.1 is a computation, not a modelling assumption.

Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

The motions of the harmonic model are the sinusoids Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)): amplitude and phase, the least period, the conserved energy kR2/2 and the elliptical phase curves

Physical scope

The complete motion analysis of the harmonic model with m>0 and k>0: every motion is a sinusoid about any chosen base time, with the initial data as coefficients; the amplitude and phase form; the least period 2 pi/omega and the frequency in cycles per second; the conserved mechanical energy k R^2/2; the phase curve is exactly the ellipse of the energy level set, with the two semiaxes; and the origin is the unique rest point and is Lyapunov stable. The result is a conditional deduction from the adopted linear force law; it asserts no value of k or m, no application to a real spring, and no asymptotic stability.

Statement

Adopt the harmonic model of The harmonic (linear restoring) model: a declared linear force law with parameter k>0: the one-dimensional model of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with m>0, k>0, force law F(x)=−kx and potential U(x)=12kx2. Write ω:=km>0 for its undamped frequency in the sense of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data. Then:

  1. Motions. Let J⊆R be an interval with nonempty interior and let t0∈J. A curve x:J→R is a motion of the model on J if and only if there are reals A,B with x(t)=Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0))(t∈J). For every datum (t0,x0,v0) there is exactly one maximal motion, it is defined on all of R, and its coefficients are A=x0, B=v0/ω; every motion with that datum is its restriction. All claims involving every real time, periods or the whole phase ellipse below refer to the unique maximal extension of clause 1; on a local motion interval the same formulas hold by restriction.

  2. Amplitude and phase. With A,B as in clause 1 put R:=A2+B2. Then ∣x(t)∣≤R for every t; if R>0 there is a real φ with cos⁡φ=A/R and sin⁡φ=B/R, and then x(t)=Rcos⁡(ω(t−t0)−φ)(t∈R). If R=0 the motion is the rest motion x≡0.

  3. Period and frequency. For every motion x(t+T)=x(t) holds for all t whenever ωT∈2πZ. If R>0, conversely every period T>0 of x satisfies ωT∈2πZ, so the least positive period of the motion is Tper=2π/ω, and the frequency in cycles per second is f=1/Tper=ω/(2π). If R=0 every T is a period.

  4. Energy and phase curves. Along every motion with coefficients A,B the mechanical energy is E(x(t),x′(t))=12k(A2+B2)=12kR2 for every t, the phase curve is contained in the ellipse ΣR:={(x,v)∈R2:x2R2+v2ω2R2=1} when R>0, and it is exactly that ellipse: every point of ΣR lies on the phase curve of the motion and is attained at exactly the times t=t0+(θ1+φ+2πn)/ω, n∈Z, for the appropriate θ1. For R=0 the phase curve is the single point (0,0), which is also the level set of energy 0.

  5. Rest point. The origin (0,0) is the only rest point of the model, and it is Lyapunov stable, by A strict local minimum of the potential is a Lyapunov stable equilibrium of the conservative one-dimensional model, and its linearisation is an oscillator-type model; no asymptotic stability is asserted.

The statement uses the declared data m,k only through ω and the energy; no value of m or k is asserted, and no real oscillator is claimed to obey the model.

Facts & Assumptions

Given: The harmonic model with m>0, k>0, ω=k/m, a base time t0, and a motion x of the model.

[L1]

Model vocabulary (The harmonic (linear restoring) model: a declared linear force law with parameter k>0, The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria clauses (b), (d), (f)): motions are the C2 curves x:J→R with m x′′+k x=0; the model is conservative with potential U(x)=12kx2, mechanical energy E(x,v)=12mv2+12kx2; and ω=k/m satisfies ω2=k/m, hence mω2=k (Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}).

[L2]

The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data with c=0: for m>0, k>0 the discriminant is Δ=−4mk<0, γ=0 and ω=4mk−02/(2m)=k/m, and a function on an interval J with nonempty interior solves m x′′+k x=0 if and only if x(t)=Acos⁡(ωt)+Bsin⁡(ωt) on J for some reals A,B; for every datum at 0 there is exactly one solution R→R, with A=x(0) and B=x′(0)/ω, and every solution extends to R.

[L3]

Trigonometric values and identities (Quarter-turn values and shifts by pi/2 and pi, Parity and the Pythagorean identity for sine and cosine, The addition formulas for complex trigonometric and hyperbolic functions, Signs, monotonicity intervals, and ranges of sine and cosine): cos⁡0=1, sin⁡0=0, cos⁡π=−1, sin⁡(π/2)=1, cos⁡(x+π)=−cos⁡x, sin⁡(x+π)=−sin⁡x, cos⁡(−x)=cos⁡x, sin⁡(−x)=−sin⁡x, sin⁡2+cos⁡2=1, cos⁡(z+w)=cos⁡zcos⁡w−sin⁡zsin⁡w for real z,w, cosine is strictly decreasing on [0,π] with range [−1,1], sine is nonnegative on [0,π], and ∣sin⁡∣≤1, ∣cos⁡∣≤1.

[L4]

The zero sets of sine and cosine and the least positive common period 2 pi: cos⁡u=0 iff u=(m+1/2)π; both sine and cosine have period 2π and no smaller positive common period. Intermediate value theorem, by bisection with a canonical left-half rule: a continuous function on [a,b] takes every value between f(a) and f(b): a continuous real function on [a,b] attains every value between its endpoint values.

[L5]

Conserved energy and level sets (Existence, uniqueness and conservation of energy for a one-dimensional conservative model: phase curves lie in energy level sets and the velocity is determined by the position and energy up to sign clauses 2 and 3, The phase plane of a one-dimensional model: phase curves, the energy function and its level sets, and turning points): E(x(t),x′(t)) is independent of t; the phase curve lies in the level set of its initial energy; the level set of energy E0≥0 is {(x,v):12kx2+12mv2=E0}.

Proof

Proof technique: direct; transport the mathematical classification of the damped oscillator to the base point t0, prove the unit-circle parametrisation used for the phase and the ellipse, and compute the energy along the explicit family.

1.1L3L4algebra

Unit-circle parametrisation. For reals c,s with c2+s2=1 there is a real θ with cos⁡θ=c and sin⁡θ=s. Indeed, c∈[−1,1] because c2≤c2+s2=1; cosine is continuous on [0,π] with cos⁡0=1 and cos⁡π=−1, so by [L4] there is θ1∈[0,π] with cos⁡θ1=c. Then sin⁡2θ1=1−c2=s2 by [L3], and sin⁡θ1≥0 on [0,π] by [L3]; hence sin⁡θ1=∣s∣. If s≥0 take θ=θ1; if s<0 take θ=−θ1, so that cos⁡θ=cos⁡θ1=c and sin⁡θ=−sin⁡θ1=−∣s∣=s by the parity identities of [L3].

1.2L3L4algebra

The periods of cosine. If T is a real number with cos⁡(u+T)=cos⁡u for all u, then T∈2πZ. Indeed, taking u=π/2 gives cos⁡(π/2+T)=0, so π/2+T=(m+1/2)π for an integer m by [L4], that is T=mπ; taking u=0 gives cos⁡T=1. From cos⁡(x+π)=−cos⁡x and cos⁡0=1 of [L3], induction on ∣m∣ gives cos⁡(mπ)=(−1)m for every integer m; hence (−1)m=1 and m is even, so T∈2πZ. Conversely T∈2πZ is a period of cosine by [L4].

1.3givenL1L2algebra

Clause 1. Let x:J→R be a motion and t0∈J. The shifted curve y(s):=x(t0+s) is defined on the interval J−t0 with nonempty interior and satisfies m y′′+k y=0 there (the equation has no explicit time dependence). By [L2] there are A,B with y(s)=Acos⁡(ωs)+Bsin⁡(ωs), that is x(t)=Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)) on J. Conversely each function of this form is obtained by shifting a solution of [L2] and is therefore a motion. For the datum (t0,x0,v0) apply [L2] to y with y(0)=x0 and y′(0)=v0: there is exactly one solution on R, namely A=x0, B=v0/ω; translating back gives the unique maximal motion, defined on R, and every motion with the datum is a restriction by the uniqueness in clause 1 of Existence, uniqueness and conservation of energy for a one-dimensional conservative model: phase curves lie in energy level sets and the velocity is determined by the position and energy up to sign. This is clause 1.

2.1givenstep 1.1L1L3algebra

Clause 2. Write θ(t):=ω(t−t0). For every real θ, R2−(Acos⁡θ+Bsin⁡θ)2=A2sin⁡2θ+B2cos⁡2θ−2ABsin⁡θcos⁡θ=(Asin⁡θ−Bcos⁡θ)2≥0, where only the expansion of the squares and sin⁡2+cos⁡2=1 of [L3] were used; hence ∣x(t)∣≤R. If R>0, apply step 1.1 to (c,s)=(A/R,B/R), which satisfies c2+s2=1: there is φ with cos⁡φ=A/R and sin⁡φ=B/R. The addition formula of [L3] gives cos⁡(θ−φ)=cos⁡θcos⁡φ+sin⁡θsin⁡φ=Acos⁡θ+Bsin⁡θR, so x(t)=Rcos⁡(ω(t−t0)−φ). If R=0 then A=B=0 and x≡0. This is clause 2.

2.2step 1.3L1L3L5algebra

Clause 4, the energy. With θ=ω(t−t0) the derivative of x is x′(t)=ω(−Asin⁡θ+Bcos⁡θ) by The derivatives of sine and cosine are cosine and minus sine, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c) and 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. Using mω2=k from [L1], E=12m x′2+12k x2=12k[(−Asin⁡θ+Bcos⁡θ)2+(Acos⁡θ+Bsin⁡θ)2]=12k(A2+B2)=12kR2, expanding the squares and using sin⁡2+cos⁡2=1 of [L3]; in particular the energy is the constant 12kR2, in agreement with the conservation clause of [L5]. This is the first assertion of clause 4.

3.1step 1.2step 2.1L4algebra

Clause 3. If ωT∈2πZ, then by [L4] the cosine is 2π-periodic, so cos⁡(θ+ωT)=cos⁡θ and x(t+T)=x(t) for all t; this covers the case R=0 as well. Suppose R>0 and T>0 is a period of x. Dividing clause 2 by R, cos⁡(ω(t+T−t0)−φ)=cos⁡(ω(t−t0)−φ) for all t; as t ranges over R, u:=ω(t−t0)−φ ranges over all of R, so ωT is a period of cosine; by step 1.2, ωT∈2πZ. Hence the positive periods of a nonconstant harmonic motion are exactly the numbers 2πn/ω, n≥1, and the least one is 2π/ω. This is clause 3.

3.2step 1.1step 1.2step 2.1step 2.2L3L5algebra

Clause 4, the ellipse. Let R>0 and use the phase form of clause 2, θ=ω(t−t0): x(t)=Rcos⁡(θ−φ) and x′(t)=−Rωsin⁡(θ−φ), so x(t)2R2+x′(t)2ω2R2=cos⁡2(θ−φ)+sin⁡2(θ−φ)=1 by [L3]: the phase curve is contained in ΣR. Conversely let (x1,v1)∈ΣR, so that (x1/R)2+(v1/(ωR))2=1; by step 1.1 applied to (c,s)=(x1/R, −v1/(ωR)) there is θ1 with cos⁡θ1=x1/R and sin⁡θ1=−v1/(ωR). The map t↦θ(t)=ω(t−t0) is a bijection of R onto R because ω>0, so the times t with θ(t)−φ=θ1+2πn are exactly t=t0+(θ1+φ+2πn)/ω, n∈Z; at those times x(t)=Rcos⁡θ1=x1 and x′(t)=−Rωsin⁡θ1=v1. Hence every point of ΣR lies on the phase curve, at the displayed times. Conversely, attainment of (x1,v1) gives both cos⁡(θ−φ)=cos⁡θ1 and sin⁡(θ−φ)=sin⁡θ1. The addition identities yield cos⁡d=1 and sin⁡d=0 for d=θ−φ−θ1, so cos⁡(u+d)=cos⁡u for every real u; step 1.2 gives d∈2πZ. Hence these are exactly the attainment times. If R=0 then x≡0 by clause 2 and the phase curve is {(0,0)}, which is the level set of energy 0 because 12kx2+12mv2=0 forces x=v=0. This is clause 4.

4.1givenL1L5L6∎

Clause 5. The rest points are the equilibria, and by [L6] the equation U′(x)=kx=0 has the unique solution x=0 because k>0; the origin is a strict global minimum of U, so [L6] gives Lyapunov stability. No asymptotic stability is asserted, and none follows from the conservation clause of [L5]. This is clause 5.

Remarks

  • Amplitude and energy are equivalent data. 12kR2 is the energy, so the amplitude and the energy determine each other; this is the source's E=κA2/2. The phase φ is determined by the datum only up to 2π, corresponding to the least period of clause 3.

  • The ellipse is the whole level set. By clause 4 the phase curve fills its ellipse; the semiaxes are R=2E0/k in the position direction and ωR=2E0/m in the velocity direction. The equilibrium (0,0) is the degenerate level set of energy 0 and is a rest point, not a periodic motion.

  • What the harmonic model does not say. Nothing here asserts that a physical oscillator is harmonic; the model's adequacy is a separate modelling question, and the companion page's pendulum check exhibits a nonlinear model together with a finite-time bound on its difference from the harmonic approximation, and shows that the harmonic energy is not conserved by the nonlinear model.

PostulateStatement: AI-adaptedProof: Not applicableOpen item page →

The damped oscillator model: a harmonic restoring force with a phenomenological linear drag −cx˙, c≥0

Physical scope

Adoption of the damped oscillator as a restricted one-dimensional model: the harmonic restoring force with spring constant k>0 combined, through additive composition of force contributions, with a linear drag force -c x-dot of declared drag coefficient c>=0 (kilograms per second). The drag term is phenomenological: it represents the net effect of a medium that is not modelled, is adopted for a stated regime (viscous, velocity-proportional resistance at small speeds), neglects the random component of the environment's force, and is a different model from quadratic drag or Coulomb friction outside that regime. The item asserts no value of c, no derivation of the drag law, no fluctuation theory and no measured result.

Postulate

(i) Setting. Work in the one-dimensional model class of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with time window I0=R, configuration domain Q=R, a fixed mass m>0, the spring constant k>0 of The harmonic (linear restoring) model: a declared linear force law with parameter k>0, and a drag coefficient c≥0 of dimension μτ−1, read in kilograms per second (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework). The model's force law is declared through the additive composition of force contributions, in the sense of clause (v) of Adopted force model: force maps, the equation of motion, additive composition of force contributions and frame covariance, read in one coordinate: F(x,v)=−k x−c v.

The damped oscillator model is the one-dimensional model with this force law; its equation of motion is m x′′+c x′+k x=0, and its data (m,c,k,0) are oscillator-type data in the sense of clause (f) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria. The model is autonomous and depends on the velocity; it is undamped when c=0, in which case it is exactly the harmonic model of The harmonic (linear restoring) model: a declared linear force law with parameter k>0, and strictly damped when c>0.

(ii) Potential presentation with a residual force. The model presents its force law through the harmonic potential U(x)=12kx2 of The harmonic (linear restoring) model: a declared linear force law with parameter k>0 together with the residual force d(t,x,v)=−c v, in the sense of clause (d) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria: F(x,v)=−U′(x)+d(t,x,v). The mechanical energy of the model is therefore the function E(x,v)=12mv2+12kx2 of clause (d), defined along every motion; the drag force is not presented through a potential and the model is not a global potential model when c>0.

(iii) Regularity. Let X:=xphys/(1 m) and V:=vphys/(1 m/s) be the dimensionless numerical SI state coordinates (the real state coordinates in the displayed force-map signature), and let K:=kphys/(1 N/m), C:=cphys/(1 kg/s) be the numerical coefficients. The numerical force map F^:=Fphys/(1 N) is F^(X,V)=−KX−CV, so ∣F^(X,V)−F^(X′,V′)∣≤K∣X−X′∣+C∣V−V′∣≤(K+C) ∥(X,V)−(X′,V′)∥2. The real number K+C is a Lipschitz constant for this normalized numerical coordinate map; the dimensionful coefficients k and c are not added. The force map is continuous and locally Lipschitz in the state variable, uniformly in time (Local Lipschitz continuity in the state variable, locally uniform in time and parameters). It is an admissible force law of the class of clause (a) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria, and its equation is exactly the homogeneous model m x′′+cx′+kx=0 of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data.

(iv) Physical status and regime. The drag term is a phenomenological constitutive model, adopted rather than derived: it represents the net force of a medium whose microscopic state is not part of the model. The linear velocity-proportional form is a restricted regime — slow motion through a viscous medium, of which the Stokes force c=6πηa is the recorded example in the framework's drag postulate The linear drag model: a phenomenological velocity-proportional resistance with a declared coefficient — and the coefficient c is a model parameter depending on body and medium. Other regimes are different models: drag growing quadratically with speed, or Coulomb sliding friction, whose force contains the factor v/∣v∣ and is undefined at v=0 and is not of the standing regularity. The model also neglects the random component of the environment's force; the recorded source states that this neglect is admissible only in a stated energy regime, in which the oscillator's energy greatly exceeds the energy scale of the environment's fluctuations. No value of c, m or k is asserted, and a measured value used in an application is an empirical premise with its own conditions and uncertainty.

(v) Status. Adopted assumption of the framework when a model invokes it; the damped oscillator selects no law of nature on its own, and no fluctuation–dissipation relation is asserted. Its conditional consequences are the subject of the damped-motion theorems of this page.

Remarks

  • Why a postulate and not a theorem. The equation of motion determines no force law; the drag contribution is additional model input, exactly as in The linear drag model: a phenomenological velocity-proportional resistance with a declared coefficient. The linear form is the simplest velocity-proportional law of the standing regularity, and it is the form of the recorded viscous-regime treatments; neither reason makes it a derivation.

  • c=0 is a parameter value, not a change of model class. The model is declared for all c≥0; at c=0 the residual force vanishes and the model is a global potential model, while for c>0 the mechanical energy is not conserved and the page's dissipation theorem applies.

Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

Classification of the motions of the damped oscillator by the discriminant, and nonincreasing energy E′=−c x′2

Physical scope

The motion classification of the damped oscillator model with m>0, k>0, c>=0: the three solution families determined by the sign of the discriminant c^2-4mk (distinct real roots, double root, oscillatory roots), uniqueness and global existence for every datum, the identity E'=-c x-dot^2 for the mechanical energy E=m x-dot^2/2+k x^2/2, nonincrease of the energy, strict negativity of its derivative when c>0 and x-dot is nonzero, and the characterisation of the constant-energy motions as the rest motion when c>0. The result is a conditional deduction from the declared model data; it asserts no value of the parameters, no measured decay and no behaviour of any real damper.

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, c≥0, force law F(x,v)=−kx−cv and mechanical energy E(x,v)=12mv2+12kx2. Write Δ:=c2−4mk,γ:=c2m,λ±:=−c±Δ2m (Δ≥0),ω:=4mk−c22m (Δ<0), and note that Δ=−4mk<0 and ω=k/m when c=0. Then:

  1. Classification and uniqueness. Let J⊆R be an interval with nonempty interior. A curve x:J→R is a motion of the model if and only if it is one of the following:
  2. Energy dissipation. For every motion x:J→R the energy e(t):=E(x(t),x′(t)) is differentiable at every interior point of J with e′(t)=−c x′(t)2≤0, and for t1<t2 in J one has e(t2)−e(t1)=−c∫t1t2x′(s)2 ds≤0: the energy is nonincreasing on J. If c>0, then the following are equivalent: (i) e is constant on J; (ii) x′(t)=0 at every interior point of J; (iii) x is constant on J; (iv) x≡0 on J, the rest motion of the model. Moreover, at every interior point with x′(t)≠0 one has e′(t)<0, so the energy strictly decreases there.
  3. Rest points. The rest points of the model are exactly the pairs (x,0) with kx=0, that is (0,0), and the constant curve x≡0 is its maximal motion with datum (t0,0,0).

The statement is read for the declared data; it asserts no value of m,k,c and no empirical decay law.

Facts & Assumptions

Given: The damped model with data (m,c,k,0), m>0, k>0, c≥0, its motions x:J→R, its energy E, and the constants Δ,γ,λ±,ω.

[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 clauses (b), (d), (f), The damped driven oscillator model: second-order equation, phase-plane system and characteristic data): the motion equation is m x′′+c x′+k x=0, which is the homogeneous damped oscillator model of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data with the same data; the force law is presented through U(x)=12kx2 with residual force d=−cv; the mechanical energy is E=12mv2+12kx2.

[L2]

The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data with m>0, c≥0, k>0: the model's solutions on an interval with nonempty interior are exactly the three displayed families by the sign of Δ; for every datum at 0 there is exactly one solution R→R, with the displayed coefficients; and every solution extends to R. For c=0 one has Δ=−4mk<0 and ω=4mk/(2m)=k/m because m>0 (Square roots exist: a unique a≥0 with (a)2=a; the positives are {x2:x≠0}).

[L4]

Integration (Newton–Leibniz needs only continuity on [a,b], differentiability on (a,b), and a Riemann-integrable extension of the interior derivative, A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a)): if G is continuous on [t1,t2] and differentiable on (t1,t2) with integrable derivative, then ∫t1t2G′=G(t2)−G(t1); a continuous real function on [t1,t2] is Riemann integrable; and if h1≤h2 are integrable on [t1,t2] with h1≤h2 pointwise, then ∫h1≤∫h2. In particular a nonpositive continuous integrand has nonpositive integral.

[L5]

Proof

technique · direct; transport the mathematical classification of the cited supplier (whose case split by the discriminant is performed there) to the model, then differentiate and integrate the energy along a motion
1.1givenL1L2

Clause 1: by [L1] the motions of the damped model are the solutions of the homogeneous damped oscillator equation with data (m,c,k) under the identification clause (f) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria; [L2] gives the three unshifted families. For a datum at t0, the curve y(s)=x(t0+s) solves the same autonomous equation with y(0)=x0, y′(0)=v0; clause 4 of its supplier fixes the coefficients of y in s. Translating that unique all-real solution back gives existence, uniqueness and the restriction property for the datum at t0, and the special case c=0 with Δ=−4mk<0, γ=0 and ω=k/m. In that case the family of the classification is x(t)=Acos⁡(ωt)+Bsin⁡(ωt), which is exactly the motion family of the harmonic model in clause 1 of The motions of the harmonic model are the sinusoids Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)): amplitude and phase, the least period, the conserved energy kR2/2 and the elliptical phase curves; hence the motions of the model with the parameter value c=0 are exactly those of the harmonic model with the same m and k.

1.2givenL1L3algebra

Clause 2: let x:J→R be a motion and e(t):=E(x(t),x′(t)). At every interior point t of J, x is twice differentiable and x′ is differentiable, so by [L3] e′(t)=m x′(t)x′′(t)+k x(t)x′(t)=x′(t)(m x′′(t)+k x(t))=−c x′(t)2, where the equation of motion was used in the last step. This is nonpositive, and it is strictly negative whenever c>0 and x′(t)≠0.

2.1step 1.2L4

For t1<t2 in J the function e is continuous on [t1,t2] and differentiable on (t1,t2) with continuous derivative by step 1.2. The derivative is Riemann integrable on this closed interval by [L4], so Newton–Leibniz gives e(t2)−e(t1)=∫t1t2e′(s) ds=−c∫t1t2x′(s)2 ds≤0, the inequality by the pointwise nonnegativity of the integrand and the monotonicity of the integral in [L4]. Hence e is nonincreasing on J.

2.2step 1.2L1L5algebra

Let c>0. If e is constant on J, then e′(t)=0 at every interior point, so x′(t)2=0 there by step 1.2, giving (ii); if (ii) holds then x is constant on J by [L5], giving (iii); if x is constant on J, then mx′′=0 and the equation reads kx=0 on J, so x≡0 because k>0 and J is nonempty, giving (iv); and (iv) trivially gives (i) because E(0,0)=0. This proves the equivalence of (i)-(iv).

3.1step 1.1L1algebra∎

Clause 3: substituting the constant curve x≡x∗ into the equation gives kx∗=0, that is x∗=0 because k>0; conversely x≡0 solves the equation, so (0,0) is a rest point and, by uniqueness in clause 1, the maximal motion with datum (t0,0,0).

Remarks

  • Two time scales in the heavily damped case. For c>0 the roots λ± are real and negative when Δ>0, with λ−<−c/(2m)<λ+<0; the source reads the reciprocal rates as a short momentum-relaxation time and a long coordinate-relaxation time. Those readings are remarks about the explicit families, not additional claims of this item.

  • Energy is not the only decreasing quantity. Clause 2 identifies the exact rate of energy loss at every instant; the monotone decrease on J follows by integration. Whether the motion tends to rest as t→∞, and how the decay compares with the explicit exponential bounds, is the subject of the following item for c>0.

Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

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.

PostulateStatement: AI-adaptedProof: Not applicableOpen item page →

The sinusoidally driven oscillator model: a prescribed external force f0cos⁡(Ωt) added to the damped oscillator

Physical scope

Adoption of the driven oscillator as a restricted one-dimensional model: the damped oscillator with parameters m>0, k>0, c>=0 together with an external driving force f(t)=f_0 cos(Omega t) of declared amplitude f_0 and angular frequency Omega>0, prescribed as a function of time and independent of the motion. The drive is an idealization of an apparatus or background that maintains a fixed sinusoidal force; it is adopted for a stated regime, its amplitude and frequency are model parameters, and no measured value, no back-reaction model, no non-sinusoidal waveform and no claim that a given apparatus realizes the model is asserted.

Postulate

(i) Setting. Work in the one-dimensional model class of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria with time window I0=R, configuration domain Q=R, mass m>0, spring constant k>0 (The harmonic (linear restoring) model: a declared linear force law with parameter k>0) and drag coefficient c≥0 (The damped oscillator model: a harmonic restoring force with a phenomenological linear drag −cx˙, c≥0). Let a real drive amplitude f0∈R of dimension μℓτ−2 (newtons) and an angular drive frequency Ω>0 of dimension τ−1 (radians per second; the radian is dimensionless in the framework's angle convention, Dimensions and SI unit conventions of the non-relativistic classical mechanics framework) be declared. Adopt, through the additive composition of force contributions of clause (v) of Adopted force model: force maps, the equation of motion, additive composition of force contributions and frame covariance and the external-contribution reading of Restricted interaction models: pair decompositions, antisymmetry and centrality options, potential forces and the gravitational branch, the force law F(t,x,v):=f0cos⁡(Ωt)−k x−c v. The sinusoidally driven oscillator model is the one-dimensional model with this force law; its equation of motion is m x′′+c x′+k x=f0cos⁡(Ωt), and its data (m,c,k,f) with f(t)=f0cos⁡(Ωt) are oscillator-type data in the sense of clause (f) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria. The term f0cos⁡(Ωt) is the drive or driving force; the model is undamped when c=0 and strictly damped when c>0.

(ii) The drive is prescribed external data. The driving force is a prescribed function of time alone: it does not depend on x or on x′, and the object producing it is not one of the model's labels and is not modelled dynamically. This is the same modelling status as the medium of The linear drag model: a phenomenological velocity-proportional resistance with a declared coefficient: the drive is an external contribution declared as part of the model. A real apparatus, background or force generator that reacts to the motion is not represented by this model; representing such a back-reaction requires a larger model with additional states.

(iii) Potential presentation with residual force. The model presents its force law through the harmonic potential U(x)=12kx2 with residual force d(t,x,v)=f0cos⁡(Ωt)−c v, in the sense of clause (d) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria: F(t,x,v)=−U′(x)+d(t,x,v). The mechanical energy of the model is E(x,v)=12mv2+12kx2 of that clause, defined along every motion; the model is driven and not a global potential model whenever f0≠0.

(iv) Regularity. Let X:=xphys/(1 m), V:=vphys/(1 m/s), and T:=tphys/(1 s) be the dimensionless numerical SI coordinates (the real inputs of the force map). Write K:=kphys/(1 N/m), C:=cphys/(1 kg/s), f^0:=f0/(1 N), and Ω^:=Ω(1 s). The numerical force map is F^(T,X,V)=f^0cos⁡(Ω^T)−KX−CV, so it is continuous and globally Lipschitz in the state variable with numerical-coordinate constant K+C, uniformly in time. This is a bound on dimensionless SI coordinates and force values, not a sum of the dimensionful k and c (Local Lipschitz continuity in the state variable, locally uniform in time and parameters). It is an admissible force law of clause (a) of The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria, and its equation is the forced model m x′′+cx′+kx=f of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data with the continuous function f(t)=f0cos⁡(Ωt).

(v) Physical status and regime. The model is a restricted idealization adopted rather than derived: exactly sinusoidal, with constant amplitude and frequency and no dependence on the motion. The drive amplitude and frequency are model parameters; when their values come from an apparatus calibration, that use is an empirical premise carrying the apparatus conditions and the measurement uncertainty, and this item asserts no such value. The sinusoidal form is the elementary component used by the recorded treatments; the response to a non-sinusoidal periodic drive is a superposition statement belonging to the mathematical theory, not an additional claim here. The model's undamped case c=0 is included as a parameter value, and it is where the resonant unboundedness of the following theorem occurs.

(vi) Status. Adopted assumption of the framework when a model invokes it; it selects no law of nature on its own, and no measured response, no resonance curve of any apparatus and no empirical resonance phenomenon is asserted by the postulate.

Remarks

  • Why a postulate and not a theorem. The equation of motion determines no force law; the drive is additional model input, exactly as the drag of The damped oscillator model: a harmonic restoring force with a phenomenological linear drag −cx˙, c≥0. The linearity of the resulting equation, which makes the explicit response formulas available, is a property of the model's declaration, not a derivation from mechanics.

  • Two readings of Ω. Ω>0 is an angular frequency in radians per second; the frequency in cycles per second is Ω/(2π), and the period of the drive is 2π/Ω. The framework's angle convention makes the radian dimensionless, so [Ω]=τ−1.

Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

The driven oscillator: steady response and amplitude profile, bounded damped and off-resonant motions, and unbounded undamped resonance

Physical scope

The response of the driven oscillator model with m>0, k>0, c>=0, drive amplitude f_0 and angular frequency Omega>0: the explicit steady sinusoidal response and its amplitude |f_0|/sqrt((k-m Omega^2)^2+c^2 Omega^2) whenever the denominator does not vanish; the reduction of every motion to the steady response plus a homogeneous motion; boundedness of every motion in the strictly damped case, with exponential convergence of the transient to the steady response; boundedness for c=0 away from Omega=sqrt(k/m); unbounded motions at c=0 and Omega=sqrt(k/m) with nonzero amplitude, so boundedness fails exactly in the off-singular-case sense; and the positive-frequency amplitude peak when f_0 is nonzero and c^2<2mk. For c^2>=2mk and f_0!=0 the amplitude strictly decreases on Omega>0 and has no maximum on that domain; for f_0=0 it is constant. The result is a conditional deduction from the declared drive and damping models; it asserts no measured resonance curve, no value of any parameter and no behaviour of an apparatus.

Statement

Adopt the sinusoidally driven oscillator model of The sinusoidally driven oscillator model: a prescribed external force f0cos⁡(Ωt) added to the damped oscillator with data m>0, k>0, c≥0, f0∈R and Ω>0, force law F(t,x,v)=f0cos⁡(Ωt)−kx−cv and equation m x′′+cx′+kx=f0cos⁡(Ωt). Put ω0:=km,D:=k−mΩ2+i c Ω∈C, and define the steady response xp(t):=Re⁡(f0eiΩtD)=f0[(k−mΩ2)cos⁡(Ωt)+c Ωsin⁡(Ωt)](k−mΩ2)2+c2Ω2(D≠0), with amplitude A(Ω):=∣f0∣/(k−mΩ2)2+c2Ω2, and in the exceptional case c=0, Ω=ω0 define xp(t):=f0tsin⁡(ω0t)/(2mω0). All all-time boundedness, convergence and unboundedness claims below refer to maximal motions on R; local motions are their restrictions as in clause 2. Then:

  1. Steady response. If D≠0 — that is c>0, or c=0 and Ω≠ω0 — then xp is a motion of the model on R, of period 2π/Ω, with amplitude A(Ω) in the sense ∣xp(t)∣≤A(Ω) for all t and equality at suitable times.
  2. Reduction to the homogeneous model. Every datum has exactly one maximal motion, it is defined on all of R, and every motion x of the model on an interval with nonempty interior can be written x=xp+y, where xp is the particular motion of clause 1, or the resonant motion in the exceptional case, and y is a motion of the damped homogeneous model with the same m, c, k (and is a harmonic motion of The harmonic (linear restoring) model: a declared linear force law with parameter k>0 when c=0).
  3. Strictly damped case c>0. For every Ω>0 the following hold: every motion is bounded on [0,∞); and x(t)−xp(t)→0 and x′(t)−xp′(t)→0 as t→+∞. Moreover xp is the unique motion with datum (t0,xp(t0),xp′(t0)) for any t0∈R. The amplitude function Ω↦A(Ω) has the following shape when f0≠0: if c2<2mk it is strictly increasing on (0,Ω∗) and strictly decreasing on (Ω∗,∞), so has its unique maximum there, with Ω∗=km−c22m2,A(Ω∗)=2m∣f0∣c4mk−c2; if c2≥2mk it is strictly decreasing on (0,∞) and does not attain a maximum on the declared positive-frequency domain. Its continuous extension to Ω=0 has its maximum there. For f0=0 the amplitude is identically zero on the domain, so every frequency has the same value.
  4. Undamped case away from resonance. For c=0 and Ω≠ω0 every motion is bounded on R.
  5. Undamped resonance. For c=0, Ω=ω0 and f0≠0 the motion xp(t)=f0tsin⁡(ω0t)/(2mω0) is unbounded on [0,∞), and every motion of the model is unbounded on [0,∞), although the drive t↦f0cos⁡(ω0t) is bounded.
  6. Off the singular cases. For c>0 the denominator never vanishes and every Ω>0 is admissible; for c=0 the exceptional frequency is exactly Ω=ω0. Consequently the boundedness assertions of clauses 3 and 4 hold for c>0, and for c=0 away from Ω=ω0, and also at undamped resonance when f0=0; for a nonzero drive clause 5 gives unboundedness at undamped resonance.

The statement is read for the declared model data; no value of any parameter and no measured response is asserted.

Facts & Assumptions

Given: The driven model with data (m,c,k,f0,Ω), its motions on intervals with nonempty interior, the constants ω0,D, the functions xp of clause 1, and the homogeneous damped model with the same m,c,k.

[L1]

Model vocabulary (The sinusoidally driven oscillator model: a prescribed external force f0cos⁡(Ωt) added to the damped oscillator, The one-dimensional force models: data, motions, phase system, potential presentation, energy and equilibria clauses (b), (f), The damped driven oscillator model: second-order equation, phase-plane system and characteristic data): motions of the driven model are the solutions of m x′′+cx′+kx=f0cos⁡(Ωt), which is the forced oscillator model of The damped driven oscillator model: second-order equation, phase-plane system and characteristic data with the continuous function f(t)=f0cos⁡(Ωt) and k>0; motions of the homogeneous damped model are the solutions with f≡0, and for c=0 these are the motions of the harmonic model by clause 1 of Classification of the motions of the damped oscillator by the discriminant, and nonincreasing energy E′=−c x′2.

[L2]

Sinusoidal forcing: the steady response, its amplitude and phase, and resonance at Ω=k/m: if D≠0 then xp=Re⁡(f0eiΩt/D) is a solution of the forced equation on R, with the displayed real form and amplitude ∣f0∣/(k−mΩ2)2+c2Ω2; if c=0 and Ω=ω0 then xp(t)=f0tsin⁡(ω0t)/(2mω0) is a solution on R and is unbounded on [0,∞) when f0≠0; if c>0 then every solution x:[0,∞)→R equals xp+xh with xh homogeneous, and x(t)−xp(t)→0, with x bounded.

[L3]

The forced oscillator with k≥0: unique global solution and the convolution (Green's function) formula: for m>0, c≥0, k>0 and continuous f on an interval I, the forced model has exactly one solution on I for every datum; and The homogeneous damped oscillator: solution classification by the discriminant and uniqueness for initial data classifies the homogeneous solutions.

[L4]

Strictly damped oscillator: explicit decay bounds, convergence to rest and asymptotic stability of the origin clauses 2 and 3: a strictly damped homogeneous motion and its derivative tend to 0 and are bounded; Classification of the motions of the damped oscillator by the discriminant, and nonincreasing energy E′=−c x′2 clause 1: the maximal motion of the homogeneous model with a given datum is unique.

[L5]

The motions of the harmonic model are the sinusoids Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)): amplitude and phase, the least period, the conserved energy kR2/2 and the elliptical phase curves clauses 2 and 3: a harmonic motion is bounded by its amplitude R, and ∣y′(t)∣≤ω0R for the motion of amplitude R because y′=−Rω0sin⁡(θ−φ).

[L6]

The amplitude of the sinusoidally forced damped oscillator: location and value of its maximum: for c>0 and f0≠0 the function A(Ω)=∣f0∣/(k−mΩ2)2+c2Ω2 is strictly increasing then strictly decreasing with the maximum Ω∗=k/m−c2/(2m2) and value 2m∣f0∣/(c4mk−c2) when c2<2mk, and strictly decreasing on (0,∞) when c2≥2mk.

[L7]

Derivatives (The derivatives of sine and cosine are cosine and minus sine, Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c)): the real form of xp is a sum of cos⁡(Ωt) and sin⁡(Ωt) with constant coefficients, so xp′(t)=f0Ω[−(k−mΩ2)sin⁡(Ωt)+c Ωcos⁡(Ωt)](k−mΩ2)2+c2Ω2, which is bounded in absolute value by Ω∣f0∣(∣k−mΩ2∣+cΩ)/((k−mΩ2)2+c2Ω2).

Proof

technique · cases; verify the explicit particular motions, reduce an arbitrary motion to the particular one plus a homogeneous motion by linearity, and use the decay, boundedness and amplitude results of the suppliers
1.1assume-case 1givenL1L2L7

Clause 1. Assume D≠0. By [L2] the function xp is a solution of the forced equation on R, hence a motion of the model by [L1], with the displayed real form and amplitude ∣xp(t)∣≤A(Ω). Its periodicity: the real form is a constant linear combination of cos⁡(Ω⋅) and sin⁡(Ω⋅), and for every real u one has cos⁡(u+2π)=cos⁡u and sin⁡(u+2π)=sin⁡u by The zero sets of sine and cosine and the least positive common period 2 pi; taking u=Ωt and u=Ω(t+2π/Ω) gives xp(t+2π/Ω)=xp(t).

1.2givenL1L3algebra

Clause 2. The forced model has k>0≥0 and the continuous drive f0cos⁡(Ωt) on I=R, so [L3] gives exactly one solution on R for every datum; hence every maximal motion is defined on R. Let x be a motion on an interval J with nonempty interior and let xp be the particular motion of clause 1 or the resonant motion of [L2]. The function y:=x−xp is of class C2 on J and m y′′+c y′+k y=(m x′′+cx′+kx)−(m xp′′+cxp′+kxp)=f0cos⁡(Ωt)−f0cos⁡(Ωt)=0 by the linearity of the differential expression (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 two equations hold on J). Hence y is a motion of the homogeneous damped model on J, and for c=0 it is a harmonic motion by [L1]. This is clause 2.

1.3givenL6

Clause 3, amplitude shape. For f0≠0 use [L6] on the amplitude function A of clause 1; for f0=0, the strictly positive denominator for c>0 gives A(Ω)=0 for every Ω directly.

2.1givenstep 1.2L2L4L7

Clause 3, transient and boundedness. Let c>0; then D≠0 because D=0 would force cΩ=0 and hence Ω=0, contrary to Ω>0. Let x be a motion; by clause 2, y=x−xp is a strictly damped homogeneous motion, so by [L4] y(t)→0 and y′(t)→0, and xp, xp′ are bounded by clause 1 and [L7]; hence x(t)→xp(t) in the sense x(t)−xp(t)→0 and x′(t)−xp′(t)→0, and x is bounded on [0,∞). If x has datum (t0,xp(t0),xp′(t0)), then y is a homogeneous motion with y(t0)=y′(t0)=0; by the uniqueness of [L4] (with the datum (t0,0,0)) y≡0, so x=xp.

2.2assume-case 2givenstep 1.1step 1.2L1L5

Clause 4. Let c=0 and Ω≠ω0. Then D=k−mΩ2=m(ω02−Ω2)≠0, so xp is a motion with ∣xp(t)∣≤A(Ω) by clause 1. Let x be a motion on an interval with nonempty interior; by clause 2, y=x−xp is a harmonic motion, so by [L5] there is Ry≥0 with ∣y(t)∣≤Ry for all t; hence ∣x(t)∣≤A(Ω)+Ry on the interval, and since clause 2 extends the motion to R, the motion is bounded on R.

2.3assume-case 3givenstep 1.2L2L5

Clause 5. Let c=0, Ω=ω0 and f0≠0. By [L2] the function xp(t)=f0tsin⁡(ω0t)/(2mω0) is a motion on R and is unbounded on [0,∞). Let x be any motion; by clause 2, y=x−xp is a harmonic motion, bounded by some Ry by [L5]. Choose tk→+∞ with ∣xp(tk)∣→+∞; then ∣x(tk)∣≥∣xp(tk)∣−∣y(tk)∣≥∣xp(tk)∣−Ry→+∞, so x is unbounded on [0,∞) as well, although ∣f0cos⁡(ω0t)∣≤∣f0∣ is bounded by Parity and the Pythagorean identity for sine and cosine. This is clause 5.

3.1cases-exhaustivegivenstep 2.1step 2.2step 2.3algebra∎

Clause 6. For c>0 the computation of step 2.1 gives D≠0 for every Ω>0, and clause 3 supplies the boundedness. For c=0, D=0  ⟺  k=mΩ2  ⟺  Ω=ω0 since Ω>0 and ω0>0; clauses 4 and 5 give boundedness away from Ω=ω0 and unboundedness at it when f0≠0. If f0=0 the equation is homogeneous, so the harmonic classification [L5] gives boundedness also at that frequency.

Remarks

  • What resonance is here. The exceptional case is a property of the declared model: with zero damping and a nonzero drive at the undamped frequency, the particular solution grows linearly and every maximal motion is unbounded. Clause 5 states no physical resonance phenomenon; an apparatus comparison would be a separate, empirically qualified step.

  • The quality factor. For c>0 and f0≠0, use the continuous extension of A to 0 and define Q:=A(ω0)/A(0)=mω0/c=ω0/(2γ), where γ=c/(2m). Likharev §5.1, Eqs. (5.20)-(5.22), PDF page 91, concerns the normalized full width of the squared response, not a half width. In the source's low-damping, small-detuning approximation the squared amplitude is replaced by f02/[4m2ω02((Ω−ω0)2+γ2)]. Its half-maximum points are Ω=ω0±γ, so their separation ΔΩ=2γ satisfies ΔΩ/ω0=1/Q. This is an identity for that approximate profile; no exact width formula for the original amplitude curve or apparatus measurement is asserted.

RemarkRemark: AI-adaptedProof: Not applicableOpen item page →

Scope and limits of the one-dimensional conservative and dissipative models: declared class, approximation status, empirical parameters and open extensions

Physical scope

An audit of the page's conditional content: which statements are declared model assumptions, which are exact conditional deductions of the model, and which are approximation or empirical interfaces; the status of the linear models as idealizations, the finite-time character of the small-angle bound, the phenomenological character of the drag and the source-reported restriction that oscillator energy greatly exceed the fluctuation energy scale, without a universal quantitative tolerance, the prescribed-drive idealization, and the handling of measured parameters as empirical premises. The remark asserts no new mathematical or physical claim and no measured value.

Scope and limits

5 · Examples, experiments and counterexamples

None yet.

Sources