Processing math: 100%
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

5.1: Free and Forced Oscillations

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

In Sec. 3.2 we briefly discussed oscillations in a keystone Hamiltonian system - a 1D harmonic oscillator described by a very simple Lagrangian 1 LT(˙q)U(q)=m2˙q2κ2q2,

whose Lagrange equation of motion, 2

 Harmonic  oscillator:  equation m¨q+κq=0, i.e. ¨q+ω20q=0, with ω20κm0,

is a linear homogeneous differential equation. Its general solution is given by (3.16), which is frequently recast into another, amplitude-phase form: q(t)=ucosω0t+vsinω0t=Acos(ω0tφ),

where A is the amplitude and φ the phase of the oscillations, which are determined by the initial conditions. Mathematically, it is frequently easier to work with sinusoidal functions as complex exponents, by rewriting the last form of Eq. (3a) in one more form: 3 q(t)=Re[Aei(ω0tφ)]=Re[aeiω0t],
aAeiφ,|a|=A,Rea=Acosφ=u,Ima=Asinφ=v.
For an autonomous, Hamiltonian oscillator, Eq. (3) gives the full classical description of its dynamics. However, it is important to understand that this free-oscillation solution, with a constant amplitude A, is due to the conservation of the energy ET+U=κA2/2 of the oscillator. If its energy changes for any reason, the description needs to be generalized.

First of all, if the energy leaks out of the oscillator to its environment (the effect usually called the energy dissipation), the free oscillations decay with time. The simplest model of this effect is represented by an additional linear drag (or "kinematic friction") force, proportional to the generalized velocity and directed opposite to it: Fv=η˙q,

where constant η is called the drag coefficient. 4 The inclusion of this force modifies the equation of motion (2) to become m¨q+η˙q+κq=0.
This equation is frequently rewritten in the form ¨q+2δ¨q+ω20q=0, with δη2m,
where the parameter δ is called the damping coefficient (or just "damping"). Note that Eq. (6) is still a linear homogeneous second-order differential equation, and its general solution still has the form of the sum (3.13) of two exponents of the type exp{λt}, with arbitrary pre-exponential coefficients. Plugging such an exponent into Eq. (6), we get the following algebraic characteristic equation for λ : λ2+2δλ+ω20=0.
Solving this quadratic equation, we get λ±=δ±iω0, where ω0(ω20δ2)1/2,
so that for not very high damping (δ<ω0)5 we get the following generalization of Eq. (3): qfree (t)=c+eλ+t+ceλt=(u0cosω0t+v0sinω0t)eδt=A0eδtcos(ω0tφ0).
The result shows that, besides a certain correction to the free oscillation frequency (which is very small in the most interesting low damping limit, δ<<ω0 ), the energy dissipation leads to an exponential decay of oscillation amplitude with the time constant τ=1/δ: A=A0et/τ, where τ1δ=2mη
A very popular dimensionless measure of damping is the so-called quality factor Q (or just the Q-factor ) that is defined as ω0/2δ, and may be rewritten in several other useful forms: Qω02δ=mω0η=(mκ)1/2η=πττ=ω0τ2,
where τ=2π/ω0 is the oscillation period in the absence of damping - see Eq. (3.29). Since the oscillation energy E is proportional to A2, i.e. decays as exp{2t/τ}, with the time constant τ/2, the last form of Eq. (11) may be used to rewrite the Q-factor in one more form: Q=ω0E(˙E)ω0EP,
where P is the dissipation power. (Two other practical ways to measure Q will be discussed below.) The range of Q-factors of important oscillators is very broad, all the way from Q10 for a human leg (with relaxed muscles), to Q104 of the quartz crystals used in electronic clocks and watches, all the way up to Q1012 for carefully designed microwave cavities with superconducting walls.

In contrast to the decaying free oscillations, the forced oscillations, induced by an external force F(t), may maintain their amplitude (and hence energy) infinitely, even at non-zero damping. This process may be described using a still linear but now inhomogeneous differential equation m¨q+η˙q+κq=F(t),

or, more conveniently for analysis, the following generalization of Eq. (6b):

 Forced  oscillator  with  damping ¨q+2δ˙q+ω20q=f(t), where f(t)F(t)/m.

For a mechanical linear, dissipative 1D oscillator (6), under the effect of an additional external force F(t), Eq. (13a) is just an expression of the 2nd  Newton law. However, according to Eq. (1.41), Eq. (13) is valid for any dissipative, linear 6 1D system whose Gibbs potential energy (1.39) has the form UG(q,t)= κq2/2F(t)q.

The forced-oscillation solutions may be analyzed by two mathematically equivalent methods whose relative convenience depends on the character of function f(t).

(i) Frequency domain. Representing the function f(t) as a Fourier sum of sinusoidal harmonics: 7 f(t)=ωfωeiωt,

and using the linearity of Eq. (13), we may represent its general solution as a sum of the decaying free oscillations (9) with the frequency ω0, independent of the function f(t), and forced oscillations due to each of the Fourier components of the force: 8 q(t)=qfree (t)+qforced (t),qforced (t)=ωaωeiωt
Plugging Eq. (15) into Eq. (13), and requiring the factors before each eiωt on both sides to be equal, we get aω=fωχ(ω),
where the complex function χ(ω), in our particular case equal to χ(ω)=1(ω20ω2)2iωδ,
is called either the response function or (especially for non-mechanical oscillators) the generalized susceptibility. From here, and Eq. (4), the amplitude of the oscillations under the effect of a sinusoidal force is Aω|aω|=|fωχ(ω)|, with |χ(ω)|=1[(ω20ω2)2+(2ωδ)2]1/2
This formula describes, in particular, an increase of the oscillation amplitude Aω at ωω0 see the left panel in Figure 1. In particular, at the exact equality of these two frequencies, |χ(ω)|ω=ω0=12ω0δ,
so that, according to Eq. (11), the ratio of the response magnitudes at ω=ω0 and ω=0(|χ(ω)|ω=0= 1/ω02 ) is exactly equal to the Q-factor of the oscillator. Thus, the response increase is especially strong in the low damping limit (δ<<ω0, i.e. Q>>1); moreover at Q and ωω0 the response diverges. (This fact is very useful for the methods to be discussed later in this section.) This is the classical description of the famous phenomenon of resonance, so ubiquitous in physics.

Screen Shot 2022-01-26 at 9.11.55 PM.png

Figure 5.1. Resonance in the linear oscillator, for several values of Q.

Due to the increase of the resonance peak height, its width is inversely proportional to Q. Quantitatively, in the most interesting low-damping limit, i.e. at Q>>1, the reciprocal Q-factor gives the normalized value of the so-called full-width at half-maximum (FWHM) of the resonance curve: 9 Δωω0=1Q.

Indeed, this Δω is defined as the difference (ω+ω)between the two values of ω at that the square of the oscillator response function, |χ(ω)|2 (which is proportional to the oscillation energy), equals a half of its resonance value (19). In the low damping limit, both these points are very close to ω0, so that in the linear approximation in |ωω0|<<ω0, we may write (ω02ω2)(ω+ω0)(ωω0)2ωξ2ω0ξ, where ξωω0
is a very convenient parameter called detuning, which will be repeatedly used later in this chapter. In this approximation, the second of Eqs. (18) is reduced to 10 |χ(ω)|2=14ω20(δ2+ξ2).
As a result, the points ω±correspond to ξ2=δ2, i.e. ω±=ω0±δ=ω0(1±1/2Q), so that Δωω+ω= ω0/Q, thus proving Eq. (20).

(ii) Time domain. Returning to arbitrary external force f(t), one may argue that Eqs. (9), (15)-(17) provide a full solution of the forced oscillation problem even in this general case. This is formally correct, but this solution may be very inconvenient if the external force is far from a sinusoidal function of time, especially if it is not periodic at all. In this case, we should first calculate the complex amplitudes fω participating in the Fourier sum (14). In the general case of a non-periodic f(t), this is actually the Fourier integral, 11 f(t)=+fωeiωtdt,

so that fω should be calculated using the reciprocal Fourier transform, fω=12π+f(t)eiωtdt
Now we may use Eq. (16) for each Fourier component of the resulting forced oscillations, and rewrite the last of Eqs. (15) as qforced (t)=+aωeiωtdω=+χ(ω)fωeiωtdω=+dωχ(ω)12π+dtf(t)eiω(tt)=+dtf(t)[12π+dωχ(ω)eiω(tt)]
with the response function χ(ω) given, in our case, by Eq. (17). Besides requiring two integrations, Eq. (25) is conceptually uncomforting: it seems to indicate that the oscillator’s coordinate at time t depends not only on the external force exerted at earlier times t<t, but also at future times. This would contradict one of the most fundamental principles of physics (and indeed, science as a whole), the causality: no effect may precede its cause.

Fortunately, a straightforward calculation (left for the reader’s exercise) shows that the response function (17) satisfies the following rule: 12 +χ(ω)eiωτdω=0, for τ<0.

This fact allows the last form of Eq. (25) to be rewritten in either of the following equivalent forms: qforced (t)=tf(t)G(tt)dt0f(tτ)G(τ)dτ,
where G(τ), defined as the Fourier transform of the response function, G(τ)12π+χ(ω)eiωτdω,
is called the (temporal) Green’s function of the system. According to Eq. (26), G(τ)=0 for all τ<0.

While the second form of Eq. (27) is frequently more convenient for calculations, its first form is more suitable for physical interpretation of the Green’s function. Indeed, let us consider the particular case when the force is a delta function f(t)=δ(tt), with t<t, i.e. τtt>0

representing an ultimately short pulse at the moment t ’, with unit "area" f(t)dt. Substituting Eq. (29a) into Eq. (27), 13 we get q(t)=G(tt).
Thus the Green’s function G(tt) is just the oscillator’s response, as measured at time t, to a short force pulse of unit "area", exerted at time t ’. Hence Eq. (27) expresses the linear superposition principle in the time domain: the full effect of the force f(t) on a linear system is a sum of effects of short pulses of duration dt ’ and magnitude f(t), each with its own "weight" G(tt) - see Figure 2 .

Screen Shot 2022-01-26 at 9.13.33 PM.png

Figure 5.2. A schematic, finite-interval representation of a force f(t) as a sum of short pulses at all times t<t, and their contributions to the linear system’s response q(t), as given by Eq. (27).

This picture may be used for the calculation of Green’s function for our particular system. Indeed, Eqs. (29)-(30) mean that G(τ) is just the solution of the differential equation of motion of the system, in our case, Eq. (13), with the replacement tτ, and a δ-functional right-hand side: d2G(τ)dτ2+2δdG(τ)dτ+ω20G(τ)=δ(τ).

Since Eqs. (27) describes only the second term in Eq. (15), i.e. only the forced, rather than free oscillations, we have to exclude the latter by solving Eq. (31) with zero initial conditions: G(0)=dGdτ(0)=0,
where τ=0 means the instant immediately preceding τ=0.

This calculation may be simplified even further. Let us integrate both sides of Eq. (31) over an infinitesimal interval including the origin, e.g. [- dτ/2,+dτ/2], and then follow the limit dτ0. Since the Green’s function has to be continuous because of its physical sense as the (generalized) coordinate, all terms on the left-hand side but the first one vanish, while the first term yields dG/dτ|+0dG/dτ|0. Due to the second of Eqs. (32), the last of these two derivatives equals zero, while the right-hand side of Eq. (31) yields 1 upon the integration. Thus, the function G(τ) may be calculated for τ>0 (i.e. for all times when it is different from zero) by solving the homogeneous version of the system’s equation of motion for τ>0, with the following special initial conditions: G(0)=0,dGdτ(0)=1.

This approach gives us a convenient way for the calculation of Green’s functions of linear systems. In particular for the oscillator with not very high damping (δ<ω0, i.e. Q>1/2), imposing the boundary conditions (33) on the homogeneous equation’s solution (9), we immediately get G(τ)=1ω0eδτsinω0τ
(The same result may be obtained directly from Eq. (28) with the response function χ(ω) given by Eq. (19). This way is, however, a little bit more cumbersome, and is left for the reader’s exercise.)

Relations (27) and (34) provide a very convenient recipe for solving many forced oscillations problems. As a very simple example, let us calculate the transient process in an oscillator under the effect of a constant force being turned on at t=0, i.e. proportional to the theta-function of time: f(t)=f0θ(t){0, for t<0,f0, for t>0,

provided that at t<0 the oscillator was at rest, so that in Eq. (15), qfree (t)0. Then the second form of Eq. (27), and Eq. (34), yield q(t)=0f(tτ)G(τ)dτ=f0t01ω0eδτsinω0τdτ.
The simplest way to work out such integrals is to represent the sine function under it as the imaginary part of exp{iω0t}, and merge the two exponents, getting q(t)=f01ω0Im[1δ+iω0eδτ+iω0τ]t0=F0k[1eδt(cosω0t+δω0sinω0t)]

This result, plotted in Figure 3, is rather natural: it describes nothing more than the transient from the initial position q=0 to the new equilibrium position q0=f0/ω02=F0/κ, accompanied by decaying oscillations. For this particular simple function f(t), the same result might be also obtained by introducing a new variable ˜q(t)q(t)q0 and solving the resulting homogeneous equation for ˜q (with appropriate initial condition  ˜q(0)=q0. However, for more complicated functions  f(t) the Green’s function approach is irreplaceable.

Screen Shot 2022-01-26 at 9.16.52 PM.png

Figure 5.3. The transient process in a linear oscillator, induced by a step-like force f(t), for the particular case δ/ω0=0.1 (i.e., Q=5 ).

Note that for any particular linear system, its Green’s function should be calculated only once, and then may be repeatedly used in Eq. (27) to calculate the system response to various external forces either analytically or numerically. This property makes the Green’s function approach very popular in many other fields of physics with the corresponding generalization or re-definition of the function. 14


1 For the notation brevity, in this chapter I will drop indices "ef" in the energy components T and U, and parameters like m,κ, etc. However, the reader should still remember that T and U do not necessarily coincide with the actual kinetic and potential energies (even if those energies may be uniquely identified) - see Sec. 3.1.

2ω0 is usually called the own frequency of the oscillator. In quantum mechanics, the Germanized version of the same term, eigenfrequency, is used more. In this series, I will use either of the terms, depending on the context.

3 Note that this is the so-called physics convention. Most engineering texts use the opposite sign in the imaginary exponent, exp{iωt}exp{iωt}, with the corresponding sign implications for intermediate formulas, but (of course) similar final results for real variables.

4 Here Eq. (5) is treated as a phenomenological model, but in statistical mechanics, such dissipative term may be derived as an average force exerted upon a system by its environment, at very general assumptions. As discussed in detail elsewhere in this series (SM Chapter 5 and QM Chapter 7), due to the numerous degrees of freedom of a typical environment (think about the molecules of air surrounding the usual mechanical pendulum), its force also has a random component; as a result, the dissipation is fundamentally related to fluctuations. The latter effects may be neglected (as they are in this course) only if E is much higher than the energy scale of the random fluctuations of the oscillator - in the thermal equilibrium at temperature T, the larger of kBT and ω0/2.

5 Systems with high damping (δ>ω0) can hardly be called oscillators, and though they are used in engineering and physics experiment (e.g., for the shock, vibration, and sound isolation), for their detailed discussion I have to refer the interested reader to special literature - see, e.g., C. Harris and A. Piersol, Shock and Vibration Handbook, 5th  ed., McGraw Hill, 2002. Let me only note that dynamics of systems with very high damping (δ≫> ω0 ) has two very different time scales: a relatively short "momentum relaxation time" 1/λ1/2δ=m/η, and a much longer "coordinate relaxation time" 1/λ+2δ/ω02=η/κ.

6 This is a very unfortunate, but common jargon, meaning "the system described by linear equations of motion".

7 Here, in contrast to Eq. (3b), we may drop the operator Re, assuming that fω=fω, so that the imaginary components of the sum compensate each other.

8 In physics, this mathematical property of linear equations is frequently called the linear superposition principle.

9 Note that the phase shift φarg[χ(ω)] between the oscillations and the external force (see the right panel in Figure 1) makes its steepest change, by π/2, within the same frequency interval Δω.

10 Such function of frequency is met in many branches of science, frequently under special names, including the "Cauchy distribution", "the Lorentz function" (or "Lorentzian line", or "Lorentzian distribution"), "the BreitWigner function" (or "the Breit-Wigner distribution"), etc.

11 Let me hope that the reader knows that Eq. (23) may be used for periodic functions as well; in such a case, fω is a set of equidistant delta functions. (A reminder of the basic properties of the Dirac δ-function may be found, for example, in MA Sec. 14.)

12 Eq. (26) remains true for any linear physical systems in which f(t) represents a cause, and q(t) its effect. Following tradition, I discuss the frequency-domain expression of this causality relation (called the KramersKronig relations) in the Classical Electrodynamics part of this lecture series - see EM Sec. 7.2.

13 Technically, for this integration, t ’ in Eq. (27) should be temporarily replaced with another letter, say t ".

14 See, e.g., EM Sec. 2.7, and QM Sec. 2.2.


This page titled 5.1: Free and Forced Oscillations is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Konstantin K. Likharev via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?