Oscillations III: Damped Driven Oscillator
A Driven Damped Oscillator: Equation of Motion
Here we take the damped oscillator analyzed in the previous lecture and add a periodic external driving force. If this driving force has the same period as the oscillator, the amplitude can increase, perhaps to disastrous proportions, as in the famous case of the Tacoma Narrows Bridge.
The equation of motion for a driven damped oscillator is:
We shall be using for the driving frequency, and for the natural frequency of the oscillator (meaning that ignoring damping, so )
The Driven Steady State Solution and Initial Transient Behavior
To get a general idea of how a damped driven oscillator behaves under a wide variety of conditions, start by exploring with our applet.
(Another approach we have used is a spreadsheet—it gives instant precise results, and can be animated, but is perhaps less intuitive.)
The solution to the differential equation above is not unique: as with any second order differential equation, there are two constants of integration, which are determined by specifying the initial position and velocity.
However, as we shall prove below using complex numbers, the equation does have a unique steady state solution with oscillating at the same frequency as the external drive. How can that be fitted to arbitrary initial conditions? The key is that we can add to the steady state solution any solution of the undriven equation and we’ll clearly still have a solution of the full damped driven equation.
We know what those undriven solutions look like: they all die away as time goes on. So, we can add such a solution to fit the specified initial conditions, and after a while the system will lose memory of those conditions and settle into the steady driven solution. The initial deviations from the steady solution needed to satisfy initial conditions are termed transients.
Here’s a pair of examples: the same driven damped oscillator, started with zero velocity, once from the origin and once from 0.5:
Notice that after about 70 seconds, the two curves are the same, both in amplitude and phase.
A much clearer illustration is the initial conditions applet: two quite different complicated curves are drawn in real time, then a click reveals that the difference between the two behaves very smoothly. Try it—then put in different initial conditions.
Using Complex Numbers to Solve the Driven Steady State Equation Easily
We begin by writing:
external driving force =
with real, so the actual physical driving force is just the real part of this, that is,
So now we’re trying to solve the equation
We’ll try the complex function, , with a real number, cycling at the same frequency as the driving force. We can always take the amplitude to be real: that is not a restriction, since we’ve added the adjustable phase factor Physically, this factor allows the solution to differ from the driver in phase, which it usually does. If we succeed in finding an that satisfies the equation, the real parts of the two sides of the equation must be equal:
If is a solution to the equation with the complex driving force, its real part, will be a solution to the equation with the real physical driving force, .
It’s very easy to check that is a solution to the equation, with the right and ! Just plug it in and see what happens. The differentiations are simple, giving
To nail down and we begin by cancelling out the common factor then shift the to the other side to find
Now is a real number, and the right hand side of this equation looks alarmingly complex, so what’s going on?
Let’s begin to untangle this by diagramming that complex number in the denominator,
It has real part and imaginary part .
Like any complex number, it can be expressed in terms of its amplitude and its phase
Putting it in the equation in this notation gives
Now, and are real, so will be real (as it must be) if is real: so and
recalling , and here (see diagram).
So we’ve already solved the differential equation: the amplitude is proportional to the strength of the driving force, and that ratio is determined by the parameters of the undriven oscillator, and the driving oscillation frequency.
The important thing to note about the amplitude is that if the damping is small, gets very large when the frequency of
the driver approaches the natural frequency of the oscillator! This is called resonance, and is what
happened to the
The phase lag of the oscillations behind the driver, , is completely determined by the frequency together with the physical constants of the undriven oscillator: the mass, spring constant, and damping strength. So, when the driving force generates the motion , the lag angle is independent of the strength of the driving force: a stronger force doesn’t get the oscillator more in sync, it just increases the amplitude of the oscillations.
Note that at low frequencies, the oscillator lags behind by a small angle, but at resonance and for driving frequencies above
Back to Reality
To summarize: we’ve just established that with and is a solution to the driven damped oscillator equation with the complex driving force
So, equating the real parts of the two sides of the equation, since are all real,
is a solution of the equation with the real driving force .
We could have found this out without complex numbers, by using a trial solution However, it’s not that easythe left hand side becomes a mix of sines and cosines, and one needs to use trig identities to sort it all out. With a little practice, the complex method is easier and is certainly more direct.
Now the total energy of the oscillator is
Note that this is not constant through the cycle unless the oscillator is at resonance,
We can see from the above that at the resonant frequency , and from the previous section
so the energy in the oscillator at the resonant frequency is
So Q, the quality factor, the measure of how long an oscillator keeps ringing, also measures the strength of response of the oscillator to an external driver at the resonant frequency.
But what happens on going away from the resonant frequency? Let’s assume that is large, and the driving force is kept constant. It won’t take much change in from for the denominator in the expression for to double in size. In fact, for large it’s a good approximation to replace by over that variation, and it is then straightforward to check that the energy in the oscillator drops to one-half its resonant value for
Exercise: prove this.
The bottom line is that for increasing the response at the resonant frequency gets larger, but this large response takes place over a narrower and narrower range in driving frequencies.
How Hard is the Driver Working to Keep this Thing Going?
It’s simplest to work with the real solution. Suppose the oscillator moves through in a time the driving force does work , so
The important thing is the average rate of working of the driving force, the mean power input, found by averaging over a complete cycle:
From , ,
averaging the power input (the bar above means average over a complete cycle) and denoting average power by
since over one cycle the average and .
(Remembering at all times, and sine is just cosine moved over, so they must have the same average over a complete cycle.)
This can be expressed entirely in terms of the driving force and frequency. Since
Exercise 1: Prove that for a lightly damped oscillator, at resonance the oscillator extracts the most work from the driving force.
Exercise 2: Prove that any solution of the damped oscillator equation (with ) can be added to the driven oscillator solution, and gives another solution to the driven oscillator. How do you pick the “right solution”?