Physlib

Physlib.ClassicalMechanics.DampedHarmonicOscillator.Basic

The Damped Harmonic Oscillator

i. Overview

The damped harmonic oscillator is a classical mechanical system corresponding to a mass `m` under a restoring force `- k x` and a damping force `- γ ẋ`, where `k` is the spring constant, `γ` is the damping coefficient, `x` is the position, and `ẋ` is the velocity.

The equation of motion for the damped harmonic oscillator is: ``` m ẍ + γ ẋ + k x = 0 ```

Depending on the relationship between the damping coefficient and the natural frequency, the system exhibits three different behaviors: - **Underdamped** (γ² < 4mk) : Oscillatory motion with exponentially decaying amplitude - **Critically damped** (γ² = 4mk) : Fastest return to equilibrium without oscillation - **Overdamped** (γ² > 4mk) : Slow return to equilibrium without oscillation

ii. Key results

This module is currently a placeholder for future implementation. The following results are planned to be formalized:

- `DampedHarmonicOscillator`: Structure containing the input data (mass, spring constant, damping coefficient) - `EquationOfMotion`: The equation of motion for the damped harmonic oscillator - Solutions for underdamped, critically damped, and overdamped cases - Energy dissipation properties - Quality factor and relaxation time

iii. Table of contents

- A. The input data (to be implemented) - B. The damped angular frequency (to be implemented) - C. The energies and energy dissipation (to be implemented) - D. The equation of motion (to be implemented) - E. Solutions (to be implemented) - E.1. Underdamped case - E.2. Critically damped case - E.3. Overdamped case - F. Quality factor and decay time (to be implemented)

iv. References

References for the damped harmonic oscillator include: - Landau & Lifshitz, Mechanics, page 76, section 25. - Goldstein, Classical Mechanics, Chapter 2.

A. The input data (placeholder)

The input data for the damped harmonic oscillator will consist of: - Mass `m > 0` - Spring constant `k > 0` - Damping coefficient `γ ≥ 0`

B. The natural angular frequency (placeholder)

The natural angular frequency ω₀ = √(k/m) will be defined here.

C. Equation of motion (Tag: DHO03)

The damped harmonic oscillator with mass `m`, spring constant `k`, and damping coefficient `γ` satisfies

m ẍ + γ ẋ + k x = 0,

where `x : Time → ℝ` is the position as a function of time.

D. The energies and energy dissipation (Tag: DHO04)

For the damped harmonic oscillator, the mechanical energy is

E(t) = ½ S.m (ẋ(t))^2 + ½ S.k (x(t))^2,

where `x : Time → ℝ` is the position as a function of time.

If `x` satisfies the equation of motion

S.m * x¨ + S.γ * ẋ + S.k * x = 0,

then differentiating `E` with respect to time and substituting the equation of motion yields

dE/dt = - S.γ * (ẋ(t))^2 ≤ 0

Thus the energy is non-increasing in time, and it is strictly decreasing whenever `S.γ > 0` and `ẋ(t) ≠ 0`. In particular, for `S.γ > 0` the energy is not conserved, and the energy dissipation rate is proportional to the squared velocity.

E. Damping regimes (placeholder)

The three damping regimes will be defined based on the discriminant γ² - 4mk.

43 declarations

theorem

The spring constant k0k \neq 0

For a damped harmonic oscillator with spring constant kk, the spring constant is non-zero, i.e., k0k \neq 0.

theorem

Mass m0m \neq 0

For a damped harmonic oscillator with mass mm, the mass is non-zero, i.e., m0m \neq 0.

definition

Natural angular frequency ω0=k/m\omega_0 = \sqrt{k/m}

For a damped harmonic oscillator with mass mm and spring constant kk, the natural (undamped) angular frequency ω0\omega_0 is defined as the square root of the ratio of the spring constant to the mass: ω0=km\omega_0 = \sqrt{\frac{k}{m}}

theorem

ω0>0\omega_0 > 0

For a damped harmonic oscillator, the natural angular frequency ω0\omega_0 is strictly positive, i.e., ω0>0\omega_0 > 0.

theorem

ω02=k/m\omega_0^2 = k/m

For a damped harmonic oscillator with mass mm and spring constant kk, the square of the natural angular frequency ω0\omega_0 is equal to the ratio of the spring constant to the mass: ω02=km\omega_0^2 = \frac{k}{m}

definition

Equation of motion mx¨+γx˙+kx=0m \ddot{x} + \gamma \dot{x} + kx = 0

For a damped harmonic oscillator with mass mm, damping coefficient γ\gamma, and spring constant kk, a position function x:TimeRx : \text{Time} \to \mathbb{R} satisfies the equation of motion if for all times tTimet \in \text{Time}, the following differential equation holds: mx¨(t)+γx˙(t)+kx(t)=0m \ddot{x}(t) + \gamma \dot{x}(t) + k x(t) = 0 where x˙(t)\dot{x}(t) and x¨(t)\ddot{x}(t) represent the first and second time derivatives of xx at time tt, respectively.

definition

Kinetic energy T(t)=12mx˙(t)2T(t) = \frac{1}{2} m \dot{x}(t)^2

For a damped harmonic oscillator with mass mm, given a position function x:TimeRx: \text{Time} \to \mathbb{R}, its kinetic energy is the function mapping each time tt to 12m(x˙(t))2 \frac{1}{2} m (\dot{x}(t))^2 where x˙(t)\dot{x}(t) is the time derivative of the position xx at time tt.

definition

Potential energy V(t)=12kx(t)2V(t) = \frac{1}{2} k x(t)^2

For a damped harmonic oscillator with spring constant kk, given a position function x:TimeRx: \text{Time} \to \mathbb{R}, its potential energy is the function t12k(x(t))2t \mapsto \frac{1}{2} k (x(t))^2.

definition

Mechanical energy E(t)=12mx˙(t)2+12kx(t)2E(t) = \frac{1}{2} m \dot{x}(t)^2 + \frac{1}{2} k x(t)^2

For a damped harmonic oscillator with mass mm and spring constant kk, given a position function x:TimeRx: \text{Time} \to \mathbb{R}, its mechanical energy is the function tE(t)t \mapsto E(t) defined as the sum of its kinetic and potential energies: E(t)=12m(x˙(t))2+12k(x(t))2 E(t) = \frac{1}{2} m (\dot{x}(t))^2 + \frac{1}{2} k (x(t))^2 where x˙(t)\dot{x}(t) is the time derivative of the position xx at time tt.

theorem

Kinetic energy T(t)=12mx˙(t)2T(t) = \frac{1}{2} m \dot{x}(t)^2

For a damped harmonic oscillator with mass mm, and for any position function x:TimeRx: \text{Time} \to \mathbb{R}, the kinetic energy is given by the function t12m(x˙(t))2t \mapsto \frac{1}{2} m (\dot{x}(t))^2, where x˙(t)\dot{x}(t) is the time derivative of xx at time tt.

theorem

Potential Energy V(t)=12kx(t)2V(t) = \frac{1}{2} k x(t)^2

For a damped harmonic oscillator SS with spring constant kk, and for any position function x:TimeRx: \text{Time} \to \mathbb{R}, the potential energy is given by the function t12k(x(t))2t \mapsto \frac{1}{2} k (x(t))^2.

definition

Energy dissipation rate γx˙2-\gamma \dot{x}^2

For a damped harmonic oscillator with damping coefficient γ\gamma, and for any position function x:TimeRx : \text{Time} \to \mathbb{R}, the energy dissipation rate is defined as the function tγ(x˙(t))2t \mapsto -\gamma (\dot{x}(t))^2, where x˙(t)\dot{x}(t) is the time derivative of xx at time tt.

theorem

Energy Dissipation Rate dEdt=γx˙2\frac{dE}{dt} = -\gamma \dot{x}^2 for the Damped Harmonic Oscillator

For a damped harmonic oscillator with damping coefficient γ\gamma, let x:TimeRx : \text{Time} \to \mathbb{R} be a position function of class CC^\infty that satisfies the equation of motion mx¨+γx˙+kx=0m \ddot{x} + \gamma \dot{x} + kx = 0. Then the time derivative of the mechanical energy E(t)E(t) at any time tt is given by dEdt(t)=γ(x˙(t))2,\frac{dE}{dt}(t) = -\gamma (\dot{x}(t))^2, where x˙(t)\dot{x}(t) denotes the time derivative of xx at time tt.

theorem

dEdt<0\frac{dE}{dt} < 0 for Non-zero Velocity and Positive Damping in a Damped Harmonic Oscillator

Consider a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ>0\gamma > 0. Let x:TimeRx : \text{Time} \to \mathbb{R} be a position function of class CC^\infty that satisfies the equation of motion mx¨+γx˙+kx=0m \ddot{x} + \gamma \dot{x} + kx = 0. If the velocity at time tt is non-zero, x˙(t)0\dot{x}(t) \neq 0, then the time derivative of the mechanical energy E(t)=12m(x˙(t))2+12k(x(t))2E(t) = \frac{1}{2} m (\dot{x}(t))^2 + \frac{1}{2} k (x(t))^2 is strictly negative: dEdt(t)<0.\frac{dE}{dt}(t) < 0.

definition

Discriminant γ24mk\gamma^2 - 4mk of a damped harmonic oscillator

The discriminant of a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma is defined as γ24mk\gamma^2 - 4mk. This value determines the damping regime (underdamped, critically damped, or overdamped) of the system.

definition

Underdamped condition γ2<4mk\gamma^2 < 4mk for a harmonic oscillator

A damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma is said to be **underdamped** if its discriminant γ24mk\gamma^2 - 4mk is strictly less than zero.

definition

Critically damped condition γ2=4mk\gamma^2 = 4mk for a harmonic oscillator

A damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma is said to be **critically damped** if its discriminant γ24mk\gamma^2 - 4mk is equal to zero.

definition

Overdamped condition γ2>4mk\gamma^2 > 4mk for a harmonic oscillator

A damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma is said to be **overdamped** if its discriminant γ24mk\gamma^2 - 4mk is strictly greater than zero.

theorem

Spring Constant k=mω2k = m\omega^2

For a damped harmonic oscillator with mass mm, spring constant kk, and natural angular frequency ω\omega, the spring constant is equal to the mass times the square of the angular frequency, i.e., k=mω2k = m\omega^2.

theorem

In a Critically Damped Oscillator, k=mβ2k = m \beta^2

For a damped harmonic oscillator with mass mm, spring constant kk, and decay rate β\beta, if the system is critically damped (i.e., γ2=4mk\gamma^2 = 4mk), then the spring constant satisfies the relation k=mβ2k = m \beta^2.

theorem

Natural Frequency is Less than Decay Rate (ω<ζ\omega < \zeta) for Overdamped Harmonic Oscillators

For a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma, if the system is overdamped (i.e., γ2>4mk\gamma^2 > 4mk), then its natural angular frequency ω=k/m\omega = \sqrt{k/m} is strictly less than its decay rate ζ=γ2m\zeta = \frac{\gamma}{2m}. That is, ω<ζ\omega < \zeta.

theorem

The Angular Frequency of an Underdamped Oscillator is Δ2m\frac{\sqrt{-\Delta}}{2m}

For a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma, if the system is underdamped (i.e., its discriminant Δ=γ24mk<0\Delta = \gamma^2 - 4mk < 0), then its damped angular frequency ω\omega is given by ω=Δ2m=4mkγ22m.\omega = \frac{\sqrt{-\Delta}}{2m} = \frac{\sqrt{4mk - \gamma^2}}{2m}.

theorem

Angular Frequency of Critically Damped Oscillator is Zero

For a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma, if the system is critically damped (i.e., γ2=4mk\gamma^2 = 4mk), then its angular frequency is zero.

theorem

Angular Frequency of an Overdamped Oscillator is Δ2m\frac{\sqrt{\Delta}}{2m}

For a damped harmonic oscillator SS that is overdamped, its angular frequency is given by ω=Δ2m,\omega = \frac{\sqrt{\Delta}}{2m}, where Δ=γ24mk\Delta = \gamma^2 - 4mk is the discriminant and mm is the mass of the oscillator.

theorem

ωd2=ω2β2\omega_d^2 = \omega^2 - \beta^2 for underdamped oscillators

For an underdamped harmonic oscillator, the square of the damped angular frequency ωd\omega_d is equal to the square of the natural angular frequency ω\omega minus the square of the decay rate β\beta: ωd2=ω2β2\omega_d^2 = \omega^2 - \beta^2 where ω=k/m\omega = \sqrt{k/m} and β\beta is the decay rate associated with the damping coefficient.

theorem

ωd>0\omega_d > 0 for an Underdamped Oscillator

For a damped harmonic oscillator SS in the underdamped regime (where the mass mm, spring constant kk, and damping coefficient γ\gamma satisfy the condition γ2<4mk\gamma^2 < 4mk), the damped angular frequency ωd\omega_d is strictly positive, i.e., ωd>0\omega_d > 0.

theorem

The damped angular frequency Ω0\Omega \neq 0 in the underdamped regime

Consider a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma. If the oscillator is in the underdamped regime (i.e., its parameters satisfy γ2<4mk\gamma^2 < 4mk), then its damped angular frequency Ω\Omega is non-zero, Ω0\Omega \neq 0.

theorem

Ω2=β2ω2\Omega^2 = \beta^2 - \omega^2 for an Overdamped Harmonic Oscillator

For a damped harmonic oscillator SS in the overdamped regime (where γ2>4mk\gamma^2 > 4mk), the square of the angular frequency Ω\Omega is equal to the square of the decay rate β\beta minus the square of the natural angular frequency ω\omega. That is: Ω2=β2ω2\Omega^2 = \beta^2 - \omega^2 where Ω\Omega is the angular frequency of the oscillator, β\beta is its decay rate, and ω=k/m\omega = \sqrt{k/m} is its natural angular frequency.

theorem

The angular frequency of an overdamped harmonic oscillator is positive (0<ωd0 < \omega_d)

For a damped harmonic oscillator SS with mass mm, spring constant kk, and damping coefficient γ\gamma, if the oscillator is overdamped (i.e., γ2>4mk\gamma^2 > 4mk), then its associated angular frequency is strictly positive (0<angularFrequency0 < \text{angularFrequency}).

theorem

The angular frequency Ω\Omega of an overdamped oscillator is non-zero (Ω0\Omega \neq 0)

For a damped harmonic oscillator SS with mass mm, spring constant kk, and damping coefficient γ\gamma, if the system is overdamped (meaning γ2>4mk\gamma^2 > 4mk), then its associated angular frequency Ω\Omega is non-zero (Ω0\Omega \neq 0).

definition

`HarmonicOscillator` associated with an undamped system SS

Given a damped harmonic oscillator SS characterized by mass mm, spring constant kk, and damping coefficient γ\gamma, this function constructs the corresponding `HarmonicOscillator` (an undamped system) provided a proof hShS that the system is undamped (i.e., γ=0\gamma = 0). The resulting harmonic oscillator retains the same mass mm and spring constant kk as the original system.

theorem

Equivalence of Damped and Undamped Equations of Motion when γ=0\gamma = 0

Consider a damped harmonic oscillator SS with mass mm, spring constant kk, and damping coefficient γ\gamma. If SS is undamped (i.e., γ=0\gamma = 0), then for any smooth trajectory x:TimeR1x: \text{Time} \to \mathbb{R}^1, xx satisfies the damped equation of motion mx¨(t)+γx˙(t)+kx(t)=0m \ddot{x}(t) + \gamma \dot{x}(t) + k x(t) = 0 if and only if it satisfies the equation of motion for the corresponding undamped harmonic oscillator.

definition

Caldirola–Kanai Lagrangian L(t,x,v)L(t, x, v) of the damped harmonic oscillator

The Caldirola–Kanai Lagrangian for a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma is a real-valued function of time tt, position xR1x \in \mathbb{R}^1, and velocity vR1v \in \mathbb{R}^1. It is defined as the Lagrangian of the corresponding undamped harmonic oscillator multiplied by a time-dependent exponential factor: L(t,x,v)=eγmt(12mv,v12kx,x) L(t, x, v) = e^{\frac{\gamma}{m} t} \left( \frac{1}{2} m \langle v, v \rangle - \frac{1}{2} k \langle x, x \rangle \right) where ,\langle \cdot, \cdot \rangle denotes the standard inner product on the 1-dimensional Euclidean space.

theorem

Explicit formula for the Caldirola–Kanai Lagrangian L(t,x,v)L(t, x, v)

For a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma, the Lagrangian LL at time tt, position xR1x \in \mathbb{R}^1, and velocity vR1v \in \mathbb{R}^1 is given by: L(t,x,v)=exp(γmt)(12mv,v12kx,x) L(t, x, v) = \exp\left(\frac{\gamma}{m} t\right) \left( \frac{1}{2} m \langle v, v \rangle - \frac{1}{2} k \langle x, x \rangle \right) where ,\langle \cdot, \cdot \rangle denotes the standard inner product on the 1-dimensional Euclidean space.

theorem

L(t,x,x˙)=eγt/m(TV)L(t, x, \dot{x}) = e^{\gamma t / m} (T - V) for the Damped Harmonic Oscillator

For a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma, the Caldirola–Kanai Lagrangian LL evaluated at time tt, position x(t)R1x(t) \in \mathbb{R}^1, and velocity x˙(t)R1\dot{x}(t) \in \mathbb{R}^1 is equal to the product of an exponential factor eγmte^{\frac{\gamma}{m} t} and the difference between the kinetic energy T(t)T(t) and potential energy V(x(t))V(x(t)): L(t,x(t),x˙(t))=eγmt(T(t)V(x(t))) L(t, x(t), \dot{x}(t)) = e^{\frac{\gamma}{m} t} \left( T(t) - V(x(t)) \right) where T(t)=12mx˙(t),x˙(t)T(t) = \frac{1}{2} m \langle \dot{x}(t), \dot{x}(t) \rangle and V(x(t))=12kx(t),x(t)V(x(t)) = \frac{1}{2} k \langle x(t), x(t) \rangle.

theorem

The Caldirola–Kanai Lagrangian of an undamped oscillator equals the harmonic oscillator Lagrangian

Let SS be a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma. If SS is undamped (i.e., γ=0\gamma = 0), then its Caldirola–Kanai Lagrangian L(t,x,v)L(t, x, v) is equal to the Lagrangian of the corresponding simple harmonic oscillator: L(t,x,v)=12mv,v12kx,x, L(t, x, v) = \frac{1}{2} m \langle v, v \rangle - \frac{1}{2} k \langle x, x \rangle, where tt is time, xR1x \in \mathbb{R}^1 is position, vR1v \in \mathbb{R}^1 is velocity, and ,\langle \cdot, \cdot \rangle denotes the standard inner product on 1-dimensional Euclidean space.

theorem

The Lagrangian of the Damped Harmonic Oscillator is CC^\infty

For a damped harmonic oscillator SS with mass mm, spring constant kk, and damping coefficient γ\gamma, the Lagrangian function L(t,x,v)L(t, x, v) (defined as the Caldirola–Kanai Lagrangian L(t,x,v)=eγmt(12mv,v12kx,x)L(t, x, v) = e^{\frac{\gamma}{m} t} \left( \frac{1}{2} m \langle v, v \rangle - \frac{1}{2} k \langle x, x \rangle \right)) is CnC^n (continuously differentiable to order nn) for any nN{}n \in \mathbb{N} \cup \{\infty\} with respect to time tt, position xx, and velocity vv.

theorem

Gradient of the Damped Harmonic Oscillator Lagrangian with respect to position xx

For a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma, the gradient of the Caldirola–Kanai Lagrangian L(t,x,v)L(t, x, v) with respect to the position xR1x \in \mathbb{R}^1 at time tt and velocity vv is given by xL(t,x,v)=(eγmtk)x. \nabla_x L(t, x, v) = -\left( e^{\frac{\gamma}{m} t} k \right) x.

theorem

vL=meγmtv\nabla_v L = m e^{\frac{\gamma}{m} t} v for the damped harmonic oscillator

For a damped harmonic oscillator with mass mm and damping coefficient γ\gamma, the gradient of the Caldirola–Kanai Lagrangian L(t,x,v)L(t, x, v) with respect to the velocity vv at time tt and position xx is: vL(t,x,v)=(meγmt)v \nabla_v L(t, x, v) = \left( m e^{\frac{\gamma}{m} t} \right) v where xx and vv are vectors in a 1-dimensional Euclidean space R1\mathbb{R}^1.

definition

Variational gradient of the action δSδx\frac{\delta S}{\delta x} for the damped harmonic oscillator

Given a trajectory x(t)x(t) mapping from time to a 1-dimensional Euclidean space R1\mathbb{R}^1, this function computes the variational gradient of the action functional S[x]=L(t,x(t),x˙(t))dtS[x] = \int L(t, x(t), \dot{x}(t)) \, dt for the damped harmonic oscillator. The action is defined using the Caldirola–Kanai Lagrangian: L(t,x,v)=eγmt(12mv,v12kx,x) L(t, x, v) = e^{\frac{\gamma}{m} t} \left( \frac{1}{2} m \langle v, v \rangle - \frac{1}{2} k \langle x, x \rangle \right) where mm is the mass, kk is the spring constant, and γ\gamma is the damping coefficient. The result is a function of time corresponding to the functional derivative δSδx\frac{\delta S}{\delta x} (the Euler–Lagrange expression).

theorem

δSδx\frac{\delta S}{\delta x} equals the Euler-Lagrange Operator for the Damped Harmonic Oscillator

For a damped harmonic oscillator with mass mm, spring constant kk, and damping coefficient γ\gamma, let x:TimeR1x : \text{Time} \to \mathbb{R}^1 be an infinitely differentiable (CC^\infty) trajectory. The variational gradient of the action, denoted as δSδx\frac{\delta S}{\delta x}, is equal to the Euler-Lagrange operator applied to the Caldirola–Kanai Lagrangian L(t,x,v)=eγmt(12mv,v12kx,x)L(t, x, v) = e^{\frac{\gamma}{m} t} \left( \frac{1}{2} m \langle v, v \rangle - \frac{1}{2} k \langle x, x \rangle \right). That is, δSδx=Lxddt(Lx˙) \frac{\delta S}{\delta x} = \frac{\partial L}{\partial x} - \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{x}} \right)

theorem

The variational gradient of the Caldirola–Kanai action equals eγmt(Fmx¨)e^{\frac{\gamma}{m} t} (F - m \ddot{x})

For a damped harmonic oscillator with mass mm and damping coefficient γ\gamma, let x(t)x(t) be a smooth (CC^\infty) trajectory in 1-dimensional Euclidean space. The variational gradient of the action functional δSδx\frac{\delta S}{\delta x} based on the Caldirola–Kanai Lagrangian is given at each time tt by: δSδx(t)=exp(γmt)(F(t)mx¨(t)) \frac{\delta S}{\delta x}(t) = \exp\left(\frac{\gamma}{m} t\right) \left( F(t) - m \ddot{x}(t) \right) where F(t)F(t) is the total physical force (the sum of the restoring and damping forces) and x¨(t)\ddot{x}(t) is the acceleration of the oscillator at time tt.

theorem

Equation of Motion iff δSδx=0\frac{\delta S}{\delta x} = 0

For a damped harmonic oscillator with mass mm, damping coefficient γ\gamma, and spring constant kk, let x:TimeRx: \text{Time} \to \mathbb{R} be an infinitely differentiable trajectory. Then xx satisfies the equation of motion mx¨(t)+γx˙(t)+kx(t)=0m \ddot{x}(t) + \gamma \dot{x}(t) + k x(t) = 0 for all tt if and only if the variational gradient of the action δSδx\frac{\delta S}{\delta x} (defined using the Caldirola–Kanai Lagrangian) is equal to zero.