harmonic oscillator

1 Introduction

Harmonic oscillator is the simplest model but one of the most important vibrating system. It is often encountered in engineering systems and commonly produced by the unbalance in rotating machinery, isolation, earthquakes, bridges, building, control and atomatization devices, just for naming a few examples. Although pure harmonic oscillation is less likely to occur than general periodic oscillations by the action of arbitrary types of excitation, understanding the behavior of a system undergoing harmonic oscillation is essential in order to comprehend how the system will respond to more general types of excitation. Such harmonic oscillations may be in the form of free vibrations or by the action of an exciter force applied at some point of the system, as we shall see later. Indeed, harmonic oscillator is a single degree of freedom system (it is like a vibration module) which is tipically constituted by an ideal spring in parallel with an ideal (but not conservative) viscous damper, both attached to a mass in one end and the other one grounded.
In this entry we will make an elementary but nonconventional discussion about the fundamental concepts involved in a single vibrating system. Unlike the conventional method, which is usually discussed in inherent Literature about this topic, where, by one side, the solution of the corresponding differential equation of motion is expressed in terms of circular trigonometric functions, having to deal with complex arbitrary constants (obviously depending on the initial conditions) and, on the other hand, studying separately the three cases of damping. Alternatively we shall use hyperbolic functions and their properties to obtain the solution of the mentioned differential equation. This way we gain a little in brevity and hilarity in the exposition, as we may study all the three cases of damping simultaneously, avoiding at the same time calculating complex constants.
Besides preliminars, we will separate the discussion in three paragraphs namely: free damped vibrations, forced damped vibrations and mechanical resonance.

2 Preliminars

The purpose of this section is concept definifitions.
Response: $x=x(t)$, the general solution of the linear differential equation involved in the motion of harmonic oscillator. We will assume $x>0$ downward, like the sense of gravitatory field.
Ideal spring: an unidimensional model, weightless, linear, elastic and that obeys the Hooke’s law $F_s=kx$, where $k$ is an elastic constant so-called stiffness, $F_s$ is the restitutive spring force and $x$ an arbitrary spring elongation.
Ideal damper: an unidimensional model, weightless, where the damping force is proportional to the velocity of oscillator, i.e. $F_d=c\dot{x}$, being $c$ the damping constant.
Natural angular frequency: ${\omega }_n=\sqrt{k/m}$ (rad/sec), a specific property of the system; $m$ is the mass of oscillator.
Damping factor: $\zeta =c/2\sqrt{km}\equiv c/2m{\omega }_n$ (physical dimensionless).
Critical damping: $c_c=2\sqrt{km}\equiv 2m{\omega }_n$. Thus $\zeta =c/c_c$, and it is also called damping ratio.
Static equilibrium configuration: a static position at $t=0^{-}$, reached because the action of gravitatory field over the mass of oscillator, i.e. the weight $mg$ ($g$ is the gravity acceleration), thus deflecting the spring a quantity $\delta$, from its natural length, so-called spring static deflection. As a result, from Newton’s equation $\sum F_x=0$, the spring force $k\delta$ (damper does not work!) due to its static deflection is balanced, at every instant, by the weight of oscillator and therefore both forces do not appear in the dynamical equation.

3 Free Damped Vibrations

With the help of an oscillator’s free body diagram and by applying the Newton’s second law $\sum F_x=m\ddot{x}$, we obtain the differential equation of motion

$$m\ddot{x}+c\dot{x}+kx=F(t),$$ (1)

being $F(t)$ the so-called exciter force. It is convenient to introduce ${\omega }_n$ and $\zeta$. So that (1) becomes

$$\ddot{x}+2\zeta {\omega }_n\dot{x}+{\omega }_n^{2}x=f(t),$$ (2)

where $f(t)=F(t)/m$. But in this section we are interested in free vibrations. So that we set $f(t)\equiv 0$ and thus we get the associate homogeneous differential equation

$$\ddot{x}+2\zeta {\omega }_n\dot{x}+{\omega }_n^{2}x=0.$$ (3)

In this case, about free oscillations, $x$ means an arbitrary displacement from the static equilibrium position. Let us suppose initial conditions $x(0)=x_0,\dot{x}(0)={\dot{x}}_0$. These conditions may be interpreted as follows: from the static equilibrium position we give an initial displacement $x_0$ and an initial velocity ${\dot{x}}_0$ to the mass $m$ and then we let it free at instant $t=0$. Accordingly, the characteristic equation of (3) is

$$s^2+2\zeta {\omega }_{n}s+{\omega }_n^2=0,$$ (4)

whose roots are given by

$$s_1=-\zeta {\omega }_n+{\omega }_n\sqrt{{\zeta }^2-1}\equiv -\zeta {\omega }_n+{\omega }_d,$$ (5)

$$s_2=-\zeta {\omega }_n-{\omega }_n\sqrt{{\zeta }^2-1}\equiv -\zeta {\omega }_n-{\omega }_d,$$ (6)

where we are defining ${\omega }_d:={\omega }_n\sqrt{{\zeta }^2-1}$ which is called frecuency of damped oscillation.
In the nullspace of the linear operator associate to (3), we find the general solution

$$x(t)=A{e}^{s_{1}t}+B{e}^{s_{2}t}⟹\dot{x}(t)=A{s}_1{e}^{s_{1}t}+B{s}_2{e}^{s_{2}t}.$$ (7)

By substituting the initial conditions in (7), we obtain the symmetric linear system of equations

$$x_0=A+B,{\dot{x}}_0=A{s}_1+B{s}_2.$$ (8)

From (5), (6) and solving (8) for $A$ and $B$, leads to

$$A=\frac{{\dot{x}}_0+(\zeta {\omega }_n+{\omega }_d)x_0}{2{\omega }_d},B=-\frac{{\dot{x}}_0+(\zeta {\omega }_n-{\omega }_d)x_0}{2{\omega }_d},$$ (9)

which substituted into the solution (7) and after an elementary algebraic manipulation, leads to

$$x(t)=\frac{e^{-\zeta {\omega }_{n}t}}{{\omega }_d}\left\{({\dot{x}}_0+\zeta {\omega }_n{x}_0)\frac{e^{{\omega }_{d}t}-e^{-{\omega }_{d}t}}{2}+{\omega }_d{x}_0\frac{e^{{\omega }_{d}t}+e^{-{\omega }_{d}t}}{2}\right\},$$ (10)

that, in terms of hyperbolic functions, it is expressed as

$$x(t)=\frac{e^{-\zeta {\omega }_{n}t}}{{\omega }_d}\left\{({\dot{x}}_0+\zeta {\omega }_n{x}_0)\sinh {\omega }_{d}t+{\omega }_d{x}_0\cosh {\omega }_{d}t\right\}.$$ (11)

It is now easily seen that the general solution of the homogeneous equation (3) satifies ${lim}_{t\to \mathrm{\infty }}x(t)=0$, since we are dealing with a transitory solution which agree with the fact that an exciter force $f(t)$ is not present.
Damping factor lets a classification about the kind of damping in concordance with the values that $\zeta$ may take. Thus,

Let us proceed to study all the three cases.

1. $\zeta >1$. Overdamped case. The response is nonperiodic and is given directly by (11), i.e.

   $$x(t)=\frac{e^{-\zeta {\omega }_{n}t}}{{\omega }_d}\left\{({\dot{x}}_0+\zeta {\omega }_n{x}_0)\sinh {\omega }_{d}t+{\omega }_d{x}_0\cosh {\omega }_{d}t\right\}.$$ (12)

   Roots $s_1>s_2\in ℝ$. Hence $e^{s_{1}t},e^{s_{1}t}$ are linearly independent and form a basis in the above mentioned nullspace.

2. $\zeta =1$. Critical damping case. The response is also nonperiodic. In addition, , so that $e^{s_{1}t},e^{s_{2}t}$ are linearly dependent and therefore they do not form a complete basis in the space of solutions of the differential equation (3). Nevertheless, we may find out a second solution linearly independent because in this case ${\omega }_d\equiv 0$, and recalling (11) we note that (by applying L’Hospital),

   $$\underset{{\omega }_d\to 0}{lim}\frac{\sinh {\omega }_{d}t}{{\omega }_d}=t.$$

   Thus, the general solution for critical damping may be expressed as

   $$x(t)\equiv x_c(t)=e^{-{\omega }_{n}t}\left\{x_0+({\dot{x}}_0+{\omega }_n{x}_0)t\right\}.$$ (13)

3. . Underdamped case. Since in this case ${\omega }_d\in ℂ$, we redefine it as ${\omega }_d=i{\omega }_n\sqrt{1-{\zeta }^2}:=i{\omega }_d^{\prime },i=\sqrt{-1}$. Recalling the hyperbolic equations $\sinh {\omega }_{d}t=\sinh i{\omega }_d^{\prime }t=i\sin {\omega }_d^{\prime }t$, $\cosh {\omega }_{d}t=\cosh i{\omega }_d^{\prime }t=\cos {\omega }_d^{\prime }t$ which substituted into (11) gives the general solution on this case. That is,

   $$x(t)=e^{-\zeta {\omega }_{n}t}\left\{\frac{{\dot{x}}_0+\zeta {\omega }_n{x}_0}{{\omega }_d^{\prime }}\sin {\omega }_d^{\prime }t+x_0\cos {\omega }_d^{\prime }t\right\}.$$ (14)

   This solution is clearly harmonics. Alternatively, we may introduce definitions about the amplitude $X$ of the response and the phase angle $\varphi$. They are,

   $$X:=\sqrt{x_0^2+{\left(\frac{{\dot{x}}_0+\zeta {\omega }_n{x}_0}{{\omega }_d^{\prime }}\right)}^2},\tan \varphi :=\frac{{\omega }_d^{\prime }{x}_0}{{\dot{x}}_0+\zeta {\omega }_n{x}_0}\cdot$$ (15)

   With these definitions placed into (14), it becomes

   $$x(t)=X{e}^{-\zeta {\omega }_{n}t}\sin ({\omega }_d^{\prime }t+\varphi ).$$ (16)

4 Forced Damped Vibrations

In this case an exciter force is present. Many useful and important applications, a few mentioned above in the Introduction, agree with an excellent degree of accuracy if we assume an intensive periodic exciter force (exciter force by mass unit) on the form $f(t)=f_0\cos \omega t$, being $f_0:=F_0/m$ the amplitude of the intensive exciter force, and $\omega$ the frecuency of the exciter force. In such case, the differential equation of motion (2) takes the form

$$\ddot{x}+2\zeta {\omega }_n\dot{x}+{\omega }_n^{2}x=f_0\cos \omega t.$$ (17)

It is well-known, from the linear differential equations theory, that the general solution of (17) may be expressed as the superposition of an homegeneous or transitory response and a particular or permanent response, i.e.

$$x(t)=x_h(t)+x_p(t).$$ (18)

The general homogeneous solution $x_h$ already have been stablished, so that in this section we focus our attention in find out a particular solution $x_p$. We search that solution on the form

$$x_p(t)=C\cos \omega t+D\sin \omega t.$$ (19)

Its first and second time derivatives are

$${\dot{x}}_p(t)=-C\omega \sin \omega t+D\omega \cos \omega t⟹{\ddot{x}}_p(t)=-C{\omega }^2\cos \omega t-D{\omega }^2\sin \omega t.$$ (20)

Next we substitute (19) and (20) into (17), setting $x(t)=x_p(t)$ and by collecting terms in $\cos \omega t$ and $\sin \omega t$, we have

$$\{({\omega }_n^2-{\omega }^2)C+(2\zeta {\omega }_n\omega )D\}\cos \omega t+\{({\omega }_n^2-{\omega }^2)D-(2\zeta {\omega }_n\omega )C\}\sin \omega t=f_0\cos \omega t.$$

We are dealing with a subspace of dimension $2$, where the set of functions $\{\cos \omega t,\sin \omega t\}\subset 𝒞[0,\mathrm{\infty })$ form an orthogonal basis in the spaces of solutions of (17). In other words, $\cos \omega t$ and $\sin \omega t$ are there linearly independent. Therefore,

$$({\omega }_n^2-{\omega }^2)D-(2\zeta {\omega }_n\omega )C\equiv 0,\text{and}\mathit{\hspace{1em}}({\omega }_n^2-{\omega }^2)C+(2\zeta {\omega }_n\omega )D\equiv f_0.$$

By solving for $C$ and $D$,

$$C=\frac{({\omega }_n^2-{\omega }^2)f_0}{{({\omega }_n^2-{\omega }^2)}^2+4{\zeta }^2{\omega }_n^2{\omega }^2},D=\frac{2\zeta {\omega }_n\omega f_0}{{({\omega }_n^2-{\omega }^2)}^2+4{\zeta }^2{\omega }_n^2{\omega }^2}\cdot$$

These one are introduced into (19) and we obtain a particular solution of (17), i.e. the permanent response of harmonic oscillator’s forced vibrations. That is,

$$x_p(t)=\frac{({\omega }_n^2-{\omega }^2)f_0}{{({\omega }_n^2-{\omega }^2)}^2+4{\zeta }^2{\omega }_n^2{\omega }^2}\cos \omega t+\frac{2\zeta {\omega }_n\omega f_0}{{({\omega }_n^2-{\omega }^2)}^2+4{\zeta }^2{\omega }_n^2{\omega }^2}\sin \omega t.$$ (21)

One also may define the frecuency ratio $r=\omega /{\omega }_n$ and continue to discuss an important set of physical considerations, but that is not the sight of this entry.

5 Mechanical Resonance

In absence of damping, i.e. $\zeta =0$, (21) reduces to

$$x_p(t)=\frac{f_0}{{\omega }_n^2-{\omega }^2}\cos \omega t.$$ (22)

One realize (theoretically) that as the frequency of exciter force coincides with the natural frequency of the system, i.e. $\omega ={\omega }_n$ (or equivalently $r=1$), then $x_p=\mathrm{\infty }$. This is what is defined as mechanical resonance phenomena. In the practice, however, low values of $\zeta$ ($\zeta \ll 1$) make feasible high values of the permanent response and thus, this fact, is also called resonance. Such phenomena is undesirable in any system, unlike the magnetic resonance that provides beneficial applications in several fields of science.

6 References

P. Fernández, Un método elemental alternativo para el estudio de las vibraciones libres amortiguadas en un oscilador armónico, (Comunicación interna), Departamento de Mecánica, Facultad de Ingeniería, Universidad Central de Venezuela (UCV), Caracas, 1997.

Title	harmonic oscillator
Canonical name	HarmonicOscillator
Date of creation	2013-03-22 17:22:46
Last modified on	2013-03-22 17:22:46
Owner	perucho (2192)
Last modified by	perucho (2192)
Numerical id	10
Author	perucho (2192)
Entry type	Application
Classification	msc 34C05
Classification	msc 34A30
Classification	msc 34-01