# 3.1.3: Harmonic Motion and Resonance Simulation

- Page ID
- 26146

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)The following simulation shows a driven, damped harmonic oscillator; a \(1\text{ kg}\) mass on a spring with spring constant \(2\text{ N/m}\). The displacement of the motion is graphed versus time. The driving frequency, \(f\), can be adjusted so we expect that for one particular frequency we will see the amplitude (maximum displacement) of the motion to be very large; in other words, resonance will occur.

Several other parameters can also be adjusted. \(b\) is the amount of friction in \(\text{Ns/m}\) (this could be air resistance or sliding friction or friction in the spring itself); \(v_{o}\) in \(\text{m/s}\) is the initial velocity of the mass, and \(F_{o}\) is the magnitude of the driving force in Newtons.

In this (and many other simulations we will use) it is easier to write \(\omega =2\pi f\) where \(\omega\) is the **angular frequency**** **in radians per second instead of having to write \(2\pi f\) everywhere.

## Simulation Questions:

- Start with the default parameters, drag the ball to the maximum starting position (or enter \(10.0\text{ m}\) in the box for \(x_{o}\)) and hit play. How does this motion compare with simple harmonic motion (in the last chapter)?
- Reset the simulation, change the friction parameter, \(b\) to \(0.5\text{ Ns/m}\), drag the ball to the maximum starting position (\(10\text{ m}\)) and hit play. (You can also click the button at the top for under damped motion.) What happens? This motion is called
.**damped harmonic motion** - Try different values for \(b\). How is the behavior of the mass for small values of \(b\) (less than one) different than for values larger than \(2\text{ Ns/m}\) (be sure to use the same starting position each time)?
- If the mass oscillates at least once before stopping the damping is called
motion. If the mass never quite gets back to the equilibrium position the motion is called*underdamped*. The case where there is just enough damping so that an oscillation does not occur (the mass just barely makes it back to equilibrium) is called*overdamped***critically damped** - Reset the simulation, set the friction parameter, \(b\) to \(1.0\text{ Ns/m}\) and the magnitude of the driving force, \(F_{o}\) to \(1.0\text{ N}\). The angular frequency should be \(1.0\text{ rad/s}\) also. Drag the ball to the maximum starting position and hit play. In this and the next few questions you will need to wait \(10\text{ s}\) or so for the motion to stabilize. Describe the stable motion after the initial oscillations for this case. This motion is called
.*driven, damped harmonic motion* - Experiment with different amounts of force, keeping friction and angular frequency equal to one. Also start from the same position each time. What is the effect of larger values of force amplitude, \(F_{o}\) on the final, stable motion of the mass?
- With damping, \(b\) set to \(0.2\text{ Ns/m}\) and \(F_{o}\) set to \(1.0\text{ N}\) try several values of the angular driving frequency, \(\omega\). Start with a value of \(1.1\text{ rad/s}\) and increase by \(0.1\) each time until you get to \(1.8\text{ rad/s}\). In each case, wait until the animation ends and measure the amplitude by clicking on the top of the curve near the end. The second number in the yellow box is the amplitude in \(\text{m}\). Write down the amplitude for each driving frequency. Which driving frequency ended up giving the largest amplitude? (This is the resonance frequency of the system.)
- The
, written as \(f_{o}\) is given by the stiffness of the spring, \(κ\), and the mass; \(f_{o} = (κ/m)^{1/2}/2π\) and is measured in \(\text{Hz}\). In this simulation the mass is \(1\text{ kg}\) and the spring constant is \(2\text{ N/m}\) so \(f_{o} = 0.225\text{ Hz}\) and the natural angular frequency, \(\omega _{o} = 2\pi f\) equals \(1.41\text{ rad/s}\). In the previous question you should have seen the maximum amplitude for a driving frequency of \(1.41\text{ rad/s}\). In other words a driving frequency, \(\omega\) of \(1.41\text{ rad/s}\) leads to resonance (maximum amplitude) because it equals the natural frequency, \(\omega _{o}\). In the previous question, did you get a maximum amplitude for a driving frequency of \(1.41\text{ rad/s}\)?**natural frequency** - Use a calculator to find the natural frequency for a spring and mass system with \(m=2.0\text{ kg}\) and \(κ = 5.0\text{ N/m}\). What do you expect the resonance frequency to be for this case?

## Advanced Questions:

- The mathematical formula that describes damped harmonic motion is \(Ae^{-\gamma t}\cos (\omega t+\phi )\) where \(\gamma =b/2m\). Notice that this is the same cosine function for simple harmonic motion but the amplitude, \(A\), is multiplied by an exponentially decreasing function of time, \(e^{-\gamma t}\). So we expect the oscillation of a damped harmonic oscillator to be an up and down cosine function with an amplitude that decreases over time. Check to see if the formula \(Ae^{-\gamma t}\cos (\omega t+\phi )\) really does describe the behavior of damped harmonic motion in the simulation. To do this, use a graphing calculator (or go to meta calculator and chose graphing calculator or use desmos calculator) and plot \(y = 10*exp(-.2 * x) * cos(1.*x)\). This is the case \(A = 10\text{ m}; b = 0.4\text{ Ns/m}; m = 1.0\text{ kg}\); and \(\omega = 1.0\text{ rad/s}\). (You can cut and paste the equation into the online calculator). How does this graph compare with the simulation for these same parameters (note: you are only interested in the positive
*x*values)? - The mathematical formula that describes driven, damped harmonic motion is \(A_{o}\cos (\omega t+\phi )\) where \(A_{o} = (F_{o} /m) / ((\omega ^{2} -\omega _{o}^{2})^{2} + 4 \gamma ^{2}\omega ^{2})^{1/2}\). In this case the amplitude, \(A_{o}\) does not change over time but it is dependent on the driving frequency, \(\omega\). What happens to the amplitude, \(A_{o}\), when \(\omega =\omega _{o}\) (assume all the other factors are constant numbers)? Is there any other combination of \(\omega\) and \(\omega _{o}\) that gives a larger amplitude?
- Verify your previous answer by making a plot of \(A_{o} = (F_{o}/m) / ((\omega ^{2} -\omega _{o}^{2})^{2} + 4 \gamma ^{2}\omega ^{2})^{1/2}\) for a range of driving frequency. To do this, use a graphing calculator to graph \(y = 1.0 / ((1.4^2 - x^2)^2 + 4*.04*1.4^2)^.5\). Where does this graph have a maximum for positive values of \(x\)?
- Based on the previous two questions, what is the resonance frequency of this system and how does this compare with the natural frequency?