PHY5200 F06

Chapter 5: Oscillations

Reading

Taylor 5.4-5.5 (today) and 5.5-5.6 (Monday).

Recall:

Damped Harmonic Oscillator

[D + γ - q][D + γ + q]x = 0
x = A± exp[-(γ±q)t]

where q² = γ² - ω0². Then we have 4 cases to consider:

Overdamped

This is mathematically the easiest case. The solution is:

x(t) = A1e-(γ-q)t + A2e-(γ+q)t
v(t) = -(γ-q)A1e-(γ-q)t - (γ+q)A2e-(γ+q)t

Critically Damped

x(t) = Ate-γt + Be-γt = (At + B)e-γt
v(t) = (A - γB - γAt)e-γt

This solution is the transition from overdamped to underdamped. It is of interest because for some devices (shock absorbers on a vehicle) this is the optimal regime.

Underdamped

When q² < 0, the solution to the DE is an exponential with a complex exponent. Recall that e = cosθ + isinθ. The solution is then:

x(t) = e-γt [Csin(ωdt) + Dcos(ωdt)] = e-γt Acos(ωdt + δ)
v(t) = -Ae-γt [γsin(ωdt + δ) - ωdcos(ωdt + δ)]

These two forms of the solution are equivalent. In the first case, C and D are the unknown constants while in the second case they are A and δ. The frequency ωd = iq = (ω0² - γ²)½ = (k/m - c²/4m²)½ and is called the damped frequency (ω0 is the undamped or natural frequency).

This solution is almost the same as the solution of the undamped oscillator. There are two differences:

  1. The amplitude decays exponentially to zero, and

  2. the frequency of oscillation is ωd rather than ω0.

Energy of the Underdamped Oscillator and Quality Factor

As when there is no damping, the energy of the oscillator at any instant in time is given by:

E = ½mv² + ½kx².

When there is no damping the energy is constant. With damping, we expect that the energy will tend to decrease with time and we find that:

dE/dt = mvv[dot] + kxx[dot] = (ma + kx)v = -cv²

since ma + cv + kx = 0 is the equation of motion. The frictional constant c is positive as well as the square of the velocity, so the change of energy with time is always negative, that is, the energy is always decreasing. The instantaneous rate of energy loss varies with the velocity, but the fractional energy loss averaged over one period of oscillation is approximately constant.

ΔE = ∫0Td dE/dt dt ≈ -γmω0² A²e-2γt Td

where Td = 2π/ωd. The approximation is valid when the exponential changes little during one cycle so that it can be treated as a constant in the integral. This result is valid when the damping is small, the so-called weak damping limit. The energy of the oscillator at time t is E(t) = ½mA²ω0²e-2γt, so we find that the fraction of energy lost is

|ΔE/E| = 2γTd

We identify 2γ to be an inverse time such that the product is unitless. Calling the time &tau, then τ=(2γ)-1, and

|ΔE/E| = Td / τ .

The quality factor, Q, is defined as 2π times the energy stored in the oscillator divided by the energy lost in a single period of oscillation. Therefore we get immediately that, in the case of weak damping,

Q = 2π τ/ Td

Damped Harmonic Oscillator, Summary of Relations

mx[ddot] + cx[dot] + kx = 0

The following table summarizes relations between the constants used in the damped oscillator results.

quantity relations comments
ω0 ω02 = k/m positive
ωd ωd2 = ω02 - γ2 positive
γ γ = c/2m positive

Four cases to consider: no damping, underdamped, critically damped, and overdamped. These differ in the amount of damping as determined by the relative values of γ and ω0.

damping γ decay parameter
none γ=0 0
under γ<ω0 γ
critical γ=ω0 γ
over γ>ω0 γ - sqrt(γ² - ω0²)

For an oscillator with a natural (undamped) frequency of ω0, let's look at how quickly the system returns to its equilibrium state as a function of the damping characterized by γ. This is determined by how quickly the exponential factor decreases to zero, and is shown in Fig. 5.13 (β = γ). To construct this graph, we imagine a system with a fixed natural frequency ω. Now we arrange things so that we can continuously vary the damping as denoted by the value of γ = 2c/m. When γ = 0 there is no damping. When γ < ω0, the damping goes like e-γt. When γ = ω0, the damping factor is e-γt = e0t. When γ > ω0, there are two damping factors, and, to consider the most general situation where the initial conditions are random, the time to return to equilibrium is determined by the slowest damping factor. The slowest damping factor is e-[γ - (γ² - ω0²)½]t, and for γ sufficiently large, this is approximately e-(ω0/γ)t. This is the result shown in Fig. 5.13. The system returns to equilibrium most quickly at critical damping, making this the optimal regime for shock absorbers on vehicles.


© Robert Harr 2006