4.3: Resonance
Resonance is a phenomenon in which an oscillator responds most strongly to a driving force that matches its own natural frequency of vibration. For example, suppose a child is on a playground swing with a natural frequency of 1 Hz. That is, if you pull the child away from equilibrium, release her, and then stop doing anything for a while, she'll oscillate at 1 Hz. If there was no friction, as we assumed in section 2.5 , then the sum of her gravitational and kinetic energy would remain constant, and the amplitude would be exactly the same from one oscillation to the next. However, friction is going to convert these forms of energy into heat, so her oscillations would gradually die out. To keep this from happening, you might give her a push once per cycle, i.e., the frequency of your pushes would be 1 Hz, which is the same as the swing's natural frequency. As long as you stay in rhythm, the swing responds quite well. If you start the swing from rest, and then give pushes at 1 Hz, the swing's amplitude rapidly builds up, as in figure a , until after a while it reaches a steady state in which friction removes just as much energy as you put in over the course of one cycle.
self-check:
In figure a , compare the amplitude of the cycle immediately following the first push to the amplitude after the second. Compare the energies as well. (answer in the back of the PDF version of the book)
What will happen if you try pushing at 2 Hz? Your first push puts in some momentum, \(p\), but your second push happens after only half a cycle, when the swing is coming right back at you, with momentum \(-p\)! The momentum transfer from the second push is exactly enough to stop the swing. The result is a very weak, and not very sinusoidal, motion, b .
Making the math easy
This is a simple and physically transparent example of resonance: the swing responds most strongly if you match its natural rhythm. However, it has some characteristics that are mathematically ugly and possibly unrealistic. The quick, hard pushes are known as impulse forces, c , and they lead to an \(x\)-\(t\) graph that has nondifferentiable kinks.
A second mathematical issue that I glossed over in the swing example was how friction behaves. In section 3.2.4 , about forces between solids, the empirical equation for kinetic friction was independent of velocity. Fluid friction, on the other hand, is velocity-dependent. For a child on a swing, fluid friction is the most important form of friction, and is approximately proportional to \(v^2\). In still other situations, e.g., with a low-density gas or friction between solid surfaces that have been lubricated with a fluid such as oil, we may find that the frictional force has some other dependence on velocity, perhaps being proportional to \(v\), or having some other complicated velocity dependence that can't even be expressed with a simple equation. It would be extremely complicated to have to treat all of these different possibilities in complete generality, so for the rest of this section, we'll assume friction proportional to velocity
\[\begin{equation*} F = -bv , \end{equation*}\]
simply because the resulting equations happen to be the easiest to solve. Even when the friction doesn't behave in exactly this way, many of our results may still be at least qualitatively correct.
3.3.1 Damped, free motion
Numerical treatment
An oscillator that has friction is referred to as damped. Let's use numerical techniques to find the motion of a damped oscillator that is released away from equilibrium, but experiences no driving force after that. We can expect that the motion will consist of oscillations that gradually die out.
\[\begin{equation*} ma+bv+kx = 0 , \end{equation*}\]
which gives symmetric treatment to three terms involving \(x\) and its first and second derivatives, \(v\) and \(a\). Now instead of calculating the time \(\Delta{}t=\Delta{}x/v\) required to move a predetermined distance \(\Delta{}x\), we pick \(\Delta{}t\) and determine the distance traveled in that time, \(\Delta{}x=v\Delta{}t\). Also, we can no longer update \(v\) based on conservation of energy, since we don't have any easy way to keep track of how much mechanical energy has been changed into heat energy. Instead, we recalculate the velocity using \(\Delta{}v=a\Delta{}t\).
import math k=39.4784 # chosen to give a period of 1 second m=1. b=0.211 # chosen to make the results simple x=1. v=0. t=0. dt=.01 n=1000 for j in range(n): x=x+v*dt a=(-k*x-b*v)/m if (v>0) and (v+a*dt<0) : print("turnaround at t=",t,", x=",x) v=v+a*dt t=t+dt
turnaround at t= 0.99 , x= 0.899919262445 turnaround at t= 1.99 , x= 0.809844934046 turnaround at t= 2.99 , x= 0.728777519477 turnaround at t= 3.99 , x= 0.655817260033 turnaround at t= 4.99 , x= 0.590154191135 turnaround at t= 5.99 , x= 0.531059189965 turnaround at t= 6.99 , x= 0.477875914756 turnaround at t= 7.99 , x= 0.430013546991 turnaround at t= 8.99 , x= 0.386940256644 turnaround at t= 9.99 , x= 0.348177318484
The spring constant, \(k=4\pi{}=39.4784\) N/m, is designed so that if the undamped equation \(f=(1/2\pi)\sqrt{k/m}\) was still true, the frequency would be 1 Hz. We start by noting that the addition of a small amount of damping doesn't seem to have changed the period at all, or at least not to within the accuracy of the calculation. 10 You can check for yourself, however, that a large value of \(b\), say 5 \(\text{N}\!\cdot\!\text{s}/\text{m}\), does change the period significantly.
We release the mass from \(x=1\ \text{m}\), and after one cycle, it only comes back to about \(x=0.9\ \text{m}\). I chose \(b=0.211\ \text{N}\!\cdot\!\text{s}/\text{m}\) by fiddling around until I got this result, since a decrease of exactly 10% is easy to discuss. Notice how the amplitude after two cycles is about \(0.81\ \text{m}\), i.e., \(1\ \text{m}\) times \(0.9^2\): the amplitude has again dropped by exactly 10%. This pattern continues for as long as the simulation runs, e.g., for the last two cycles, we have 0.34818/0.38694=0.89982, or almost exactly 0.9 again. It might have seemed capricious when I chose to use the unrealistic equation \(F=-bv\), but this is the payoff. Only with \(-bv\) friction do we get this kind of mathematically simple exponential decay.
Because the decay is exponential, it never dies out completely; this is different from the behavior we would have had with Coulomb friction, which does make objects grind completely to a stop at some point. With friction that acts like \(F=-bv\), \(v\) gets smaller as the oscillations get smaller. The smaller and smaller force then causes them to die out at a rate that is slower and slower.
Contributors and Attributions
Benjamin Crowell (Fullerton College). Conceptual Physics is copyrighted with a CC-BY-SA license.