15.1: Undamped Free Vibrations
- Page ID
- 50618
\( \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}\)Vibrations occur in systems that attempt to return to their resting or equilibrium state when perturbed, or pushed away from their equilibrium state. The simplest vibrations to analyze are undamped, free vibrations with one degree of freedom.
"Undamped" means that there are no energy losses with movement (whether the losses are intentional, from adding dampers, or unintentional, through drag or friction). An undamped system will vibrate forever without any additional applied forces. A simple pendulum has very low damping, and will swing for a long time before stopping. "Damped" means that there are resistive forces and energy losses with movement that cause the system to stop moving eventually.
"Free" means that, after the initial perturbation, the only forces acting on the system are internal to the system (springs, dampers) and/or gravity. A tuning fork continues to vibrate after the one initial perturbation of being struck. In contrast, "forced" means there is an external, typically periodic, force acting on the system. A jackhammer vibrates due to having a supply of compressed air continually forcing the bit up and down, and it stops vibrating very quickly without that external periodic forcing.
"One degree of freedom" means that we will only consider systems with one mass vibrating along one direction (e.g. use variable \(x\)) or about one axis (e.g. use variable \(\theta\)). Systems having more than one mass or vibrating along or about two or more axes have more than one degree of freedom.
We can derive the equation of the system by setting up a free body diagram. Consider a mass sitting on a frictionless surface, attached to a wall via a spring.
The system above is in equilibrium. It is at rest, and will stay at rest unless some other force acts on it.
To start the system vibrating, we need to perturb it. Perturbation is moving the system away from equilibrium by a small amount.
When we perturb this system, we either stretch or compress the spring. This generates a spring force, and the spring force is always in a direction that tends to pull the system back toward equilibrium.
We can generate the equation of motion of the system, and determine the specifics of how it will vibrate, by analyzing this perturbed state. Recall that the spring force or moment is:
\[ \vec{F}_k = k \vec{x} \]
\[ \vec{M}_k = k \vec{\omega} \]
Note that the spring constants in the above equations have different units, depending on whether the spring is linear (Newton/meter) or torsional (Newton-meter/radian), and that \(\theta\) must be given in radians. The magnitude of the spring force depends on \(x\), the distance perturbed from the spring's unstretched length (not necessarily the equilibrium position of the system), and the same is true for the moment of a torsional spring. The spring force or moment is in the direction/orientation opposite that of the displacement. That is, if you pull the mass to the right, the spring force points to the left.
The process for finding the equation of motion of the system is as follows:
- Sketch the system with a small positive perturbation (\(x\) or \(\theta\)).
- Draw the free body diagram of the perturbed system. Ensure that the spring force has a direction opposing the perturbation.
- Find the one equation of motion for the system in the perturbed coordinate using Newton's Second Law. Keep the same positive direction for position, and assign positive acceleration in the same direction.
- Move all terms of the equation to one side, and check that all terms are positive. If all terms are not positive, there is an error in the direction of displacement, acceleration, and/or spring force.
For the example system above, with mass \(m\) and spring constant \(k\), we derive the following:
\[ \sum F_x = m a_x = m \ddot{x} \]
\[ -F_k = m \ddot{x} \]
\[ -kx = m \ddot{x} \]
\[ m \ddot{x} + kx = 0 \]
This gives us a differential equation that describes the motion of the system. We can rewrite it in normal form:
\[ m \ddot{x} + kx = 0 \]
\[ \Rightarrow \ddot{x} + \frac{k}{m} x = 0 \]
\[ \Rightarrow \ddot{x} + \omega_n^2 x = 0 \]
The term \(\omega_n\) is called the angular natural frequency of the system, and has units of radians/second.
\[ \omega_n^2 = \frac{k}{m} \]
\[ \omega_n = \sqrt{\frac{k}{m}} \]
Assuming that the initial perturbation of the system can be described by the position and velocity of the mass at \(t=0\), then:
\[ x(0) = x_0 \]
\[ v(0) = \dot{x}(0) = v_0. \]
The solution to the differential equation, that provides the position \(x(t)\) of the system at time \(t\), is:
\[ x(t) = C \sin \left( \omega_n t + \phi \right), \]
\[ \text{where:} \quad \omega_n = \sqrt{\frac{k}{m}}, \,\, C = \sqrt{ \left( \frac{v_0}{\omega_n} \right)^2 + x_0^2 }, \,\, \phi = \tan ^{-1} \left( \frac{x_0 \omega_n}{v_0} \right). \]
The amplitude \(C\) describes the maximum displacement during the oscillations (i.e. \(x_{max}\)), and the phase \(\phi\) describes how the sine function is shifted in time.
For a system where there is torsional vibration (that is, the oscillation involves a rotation), the equations are similarly:
\[ I \ddot{\theta} + k \theta = 0 \]
\[ \Rightarrow \ddot{\theta} + \frac{k}{I} \theta = 0 \]
\[ \Rightarrow \ddot{\theta} + \omega_n^2 \theta = 0, \]
\[ \text{where} \,\, \omega_n = \sqrt{\frac{k}{I}}. \]
Assuming that the initial perturbation of the system can be described by the position and velocity of the mass at \(t=0\),
\[x(0) = x_0 \]
\[ v(0) = \dot{x}(0) = v_0. \]
The solution to the differential equation, that provides the position \(x(t)\) of the system at time \(t\), is:
\[ x(t) = C \sin (\omega_n t + \phi) \]
\[ \text{where:} \quad \omega_n = \sqrt{\frac{k}{I}}, \,\, C = \sqrt{ \left( \frac{\omega_0}{\omega_n} \right) ^2 + \theta_0^2 }, \,\, \phi = \tan ^{-1} \left( \frac{\theta_0 \omega_n}{\omega_0} \right). \]
Find an expression for the angular natural frequency of the following system, and find the maximum amplitude of vibration of the system with mass \(m\) = 10 kg and spring constant \(k\) = 200 N/m when given an initial displacement of \(x_0\) = 0.1 m and an initial velocity of \(v_0\) = 0.3 m/s.
- Solution
-
Video \(\PageIndex{1}\): Worked solution to example problem \(\PageIndex{1}\). YouTube source: https://youtu.be/J1TVxxVjV_c.
Determine the equation of motion of the system from Newton's Second Law. Assume mass \(m\) = 5 kg and spring constant \(k\) = 500 N/m. Find the initial displacement, \(x_0\), such that the mass oscillates over a total range of 4 meters. Assume the initial perturbation velocity, \(v_0\), is 10 m/s.
- Solution
-
Video \(\PageIndex{2}\): Worked solution to example problem \(\PageIndex{2}\). YouTube source: https://youtu.be/TAy412iVvwE.