The Lab

The Lecture Room

The Equations of Motion

Oscillations and Resonances

Nonlinear Dynamics

Miscellaneous Topics

Links and Literature

The Harmonic Oscillator

Nonlinear Oscillations

Nonlinear Resonance

Parametric Resonance

The Harmonic Oscillator

The harmonic oscillator is a canonical system discussed in every freshman course of physics. Thus, you might skip this lecture if you are familiar with it.

1. Basic equations of motion and solutions

The linearized equation of motion of an undamped and undriven pendulum is called a harmonic oscillator:


d2phi/dt2 + omega02phi = 0,

Nearly all oscillators and oscillations in physics are modeled by this equation of motion, at least in a first approximation, because it can be solved analytically. It is a linear second-order differential equation with constant coefficients. Such a differential equation can be solved easily by the ansatz phi= exp(lambdat), where lambda is some constant which has to be calculated. Each d/dt puts a factor lambda before the exponential function. Thus, a linear ordinary differential equation like (1) turns into a polynomial because the common exponential function (a nonzero factor) can be dropped. Thus, lambda2 + omega02 = 0. The roots of this so-called characteristic polynomial are lambda = ±iomega0. Hence, we get two solutions. Because the equation of motion is linear, one gets new solutions by multiplying a known solution by an arbitrary factor or by adding a solution to another one. This so-called superposition principle of linear differential equations is very powerful. It says that any solution can be built up by only a few different types of "building blocks" called fundamental solutions. The solutions exp(-iomega0t) and exp(+iomega0t) are a set of two independent fundamental solutions of the harmonic oscillator. Thus, the general solution of (1) reads


phi(t) = C exp(iomega0t) + C* exp(-iomega0t),

where C is a complex number and C* denotes its conjugated complex. Because of the superposition principle, the set of fundamental solutions can be arbitrary, except that a fundamental solution has to be independent from the other ones. That is, each fundamental solution can not be written as a weighted sum of the other fundamental solutions. Instead of the complex valued exponential functions, one often chooses sin(omega0t)=(exp(iomega0t)-exp(-iomega0t))/2i and cos(omega0t)=(exp(iomega0t)+exp(-iomega0t))/2 as another set of fundamental solutions. Therefore,

phi(t) = A sin(omega0t) + B cos(omega0t),

is a general solution, too. The constants A and B can be expressed in term of C and vice versa.

Using the porperties of trigonometric functions, the general solution can also be expressed in a physically more relevant way:

phi(t) = phimax sin(omega0t-psi),

where phimax and psi are the amplitude and the phase of oscillation, respectively. Again they can be expressed in terms of the integration constants of (2) or (3). Note, that these constants are determined by the initial conditions phi(0) and dphi/dt|0. The period of oscillation is T0=2pi/omega0. These three quantities are independent from each other. This isn't a triviality as nonlinear oscillations show. In fact, the independence of the frequency from the amplitude is a more subtile definition of a harmonic oscillator.

2. Damping

What happens if the harmonic oscilliator is damped? In this case, the equation of motions has an additional term which comes from the damping force:


d2phi/dt2 + gammadphi/dt + omega02phi = 0.

Again, with the ansatz phi= exp(lambdat), one gets a characteristic polynomial: lambda2 + gammalambda + omega02 = 0. It has the solutions


lambda± = -gamma/2 ± (gamma2/4 - omega02)1/2.

Depending on wether gamma2/4 is less than, equal, or greater than omega02, one has to distinguish three cases:

The underdamped case: gamma/2 < omega0
The eigenvalues lambda± are a conjugated complex pair. The imaginary part (omega02 - gamma2/4 )1/2 is just the frequency omega of a damped oscillation. The real part determines the rate of decay. The general solution reads

phi(t) = phimax exp(-gammat/2) sin(omegat+psi),

The exponentially decaying amplitude of the oscillation is called the envelope.
The critically damped case: gamma/2 = omega0
Both eigenvalues are identically. In this so-called degenerated case, the ansatz phi = exp(lambdat) is not generally enough because it leads only to one solution. But for a fundamental system two independent solutions are needed. Where do we get a second solution? In such a case the ansatz phi= t exp(lambdat) helps. Thus the general solution reads

phi(t) = (At+B) exp(-gammat/2).

The overdamped case: gamma/2 > omega0
Both eigenvalues are negative and real. The general solution reads

phi(t) = A+ exp(lambda+t) + A- exp(lambda-t).

Any initial condition leads to an exponentially decay. There is only one extrema (minimum or maximum) contrary to the underdamped case with infinitely many maxima and minima. The decay is dominated by the eigenvalue most closest to zero.

3. Resonance

Oscillations in a damped harmonic oscillator decay to zero. When it is driven by a periodic force, one oscillation survives. The equation of motion of a damped and driven harmonic oscillator reads


d2phi/dt2 + gammadphi/dt+ omega02phi = a cos 2pif t.

This is still a linear differential equation, but the sum of two solutions are no longer a solution. The reason for that is the inhomogenous term on the right-hand side. Nevertheless, the superposition principle holds in a modified version. A solution of (9) plus an arbitray solution of (5) (a so-called homogeneous solution) is also a solution of (9). One needs only one solution of (9) (the so-called particular solution) to generate any solution of (9) with the help of the homogeneous solution. Here, a particular solution can be found by the ansatz


phi(t) = A sin(2pif t) + B cos(2pif t).

It turns (9) into a sum with terms proportionally to cos(2pif t) and sin(2pif t). Because the equation should hold for any value of t, the sum of all cosinus terms and the sum of all sinus terms have to be zero independently. Hence, [omega02-(2pif)2] B + 2pifgamma A = a and [omega02-(2pif)2] A - 2pifgamma B = 0, respectively. The solution of these equations reads


A = 2pifgammaa / ([omega02-(2pif)2]2 + (2pifgamma)2),

B = [omega02-(2pif)2] a / ([omega02-(2pif)2]2 + (2pifgamma)2).

In the long-time limit any solution of (9) approaches the particular solution (10) since any solution of (5) decays to zero. The harmonic oscillator, therefore, oscillates not with its eigenfrequency omega0 but with the frequency of the periodic force. The amplitude of oscillations depend on the driving frequency. It has its maximum when the driving frequency matches the eigenfrequency. This phenomenon is called resonance. To see this, rewrite equation (10) in the form phi(t) = phimax cos(2pift-psi), with

 phimax = a / ([omega02-(2pif)2]2 + (2pifgamma)2)1/2


psi = arctan (2pifgamma/[omega02-(2pif)2]).

In the underdamped case, the amplitude as a function of the driving frequency has a maximum near the eigenfrequency f0 = omega0/2pi of the oscillator. The deviation is of quadratic order in gamma. The maximum is a/(omega0gamma) in leading order of gamma. This is by a factor omega0/gamma larger than the amplitude for fto0. That is, the amplitude of oscillation is larger than the amplitude of driving. This can be easily observed in the lab at the horizontally driven pendulum. The width of the so-called resonance line is proportional to gamma. In the critically damped and overdamped case the resonance line disappears.

In the limit fto0 the phase of oscillation is identical with the driving phase. In resonance, the phase is just 90° behind the driving phase. At high frequencies the phase approaches 180°. Thus, the pendulum moves just opposite to the drive. Again, this can be easily observed in the lab at the horizontally driven pendulum.

QUESTIONS worth to think about:
  1. What happens with an undamped harmonic oscillator driven exactly by its eigenfrequency?
  2. Another look at the dynamics of the damped and driven harmonic oscillator is the following one: Instead of discussing the solution as a function of time, discuss it as a function of the initial conditions. For example, calculate the angle and the angular velocity at a given time as a function of the angle and the angluar velocity for the previous period of driving. That is, given are the initial conditions phi(0) and dphi/dt|0. Calculate phi(1/f) and dphi/dt|1/f as a function of the initial conditions. This calculation leads to a Poincaré map or stroboscopic map.







© 1998 Franz-Josef Elmer,  elmer@ubaclu.unibas.ch last modified Sunday, July 19, 1998.