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.

18 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.