Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

5.3: Reduced Equations

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

A much more important issue is the stability of the solutions described by Eq. (48). Indeed, Figure 4 shows that within a certain range of parameters, these equations give three different values for the oscillation amplitude (and phase), and it is important to understand which of these solutions are stable. Since these solutions are not the fixed points in the sense discussed in Sec. 3.2 (each point in Figure 4 represents a nearly-sinusoidal oscillation), their stability analysis needs a more general approach that would be valid for oscillations with amplitude and phase slowly evolving in time. This approach will also enable the analysis of non-stationary (especially the initial transient) processes, which are of importance for some dynamic systems.

First of all, let us formalize the way the harmonic balance equations, such as Eqs. (47), are obtained for the general case (38) - rather than for the particular Eq. (43) considered in the last section. After plugging in the 0th  approximation (41) into the right-hand side of equation (38) we have to require the amplitudes of both quadrature components of frequency ω to vanish. From the standard Fourier analysis, we know that these requirements may be represented as ¯f(0)sinΨ=0,¯f(0)cosΨ=0,

where the top bar means the time averaging - in our current case, over the period 2π/ω of the right-hand side of Eq. (52), with the arguments calculated in the 0th  approximation: f(0)f(t,q(0),˙q(0),)f(t,AcosΨ,AωsinΨ,), with Ψ=ωtφ.
Now, for a transient process the contribution of q(0) to the left-hand side of Eq. (38) is not zero any longer, because its amplitude and phase may be both slow functions of time - see Eq. (41). Let us calculate this contribution. The exact result would be ¨q(0)+ω2q(0)(d2dt2+ω2)Acos(ωtφ)=(¨A+2˙φωA˙φ2A)cos(ωtφ)2˙A(ω˙φ)sin(ωtφ).
However, in the first approximation in ε, we may neglect the second derivative of A, and also the squares and products of the first derivatives of A and φ (which are all of the second order in ε ), so that Eq. (54) is reduced to ¨q(0)+ω2q(0)2A˙φωcos(ωtφ)2˙Aωsin(ωtφ).
On the right-hand side of Eq. (53), we can neglect the time derivatives of the amplitude and phase at all, because this part is already proportional to the small parameter. Hence, in the first order in ε, Eq. (38) becomes ¨q(1)+ω2q(1)=f(0)eff(0)(2A˙φωcosΨ2˙AωsinΨ).
Now, applying Eqs. (52) to the function f(0)ef, and taking into account that the time averages of sin2Ψ and cos2Ψ are both equal to 1/2, while the time average of the product sinΨcosΨ vanishes, we get a pair of so-called reduced equations (alternatively called "truncated", or "RWA", or "van der Pol" equations) for the time evolution of the amplitude and phase: ˙A=1ω¯f(0)sinΨ,˙φ=1ωA¯f(0)cosΨ
Extending the definition (4) of the complex amplitude of oscillations to their slow evolution in time, a(t) A(t)exp{iφ(t)}, and differentiating this relation, the two equations (57a) may be also rewritten in the form of either one equation for a : ˙a=iω¯f(0)ei(Ψ+φ)iω¯f(0)eiωt
or two equations for the real and imaginary parts of a(t)=u(t)+iv(t) : ˙u=1ω¯f(0)sinωt,˙v=1ω¯f(0)cosωt.
The first-order harmonic balance equations (52) are evidently just the particular case of the reduced equations (57) for stationary oscillations ( ˙A=˙φ=0).21

Superficially, the system (57a) of two coupled, first-order differential equations may look more complex than the initial, second-order differential equation (38), but actually, it is usually much simpler. For example, let us spell them out for the easy case of free oscillations a linear oscillator with damping. For that, we may reuse the ready Eq. (46) by taking α=f0=0, and thus turning Eqs. (57a) into ˙A=1ω¯f(0)sinΨ1ω¯(2ξωAcosΨ+2δωAsinΨ)sinΨδA,˙φ=1ωA¯f(0)cosΨ1ωA¯(2ξωAcosΨ+2δωAsinΨ)cosΨξ.

The solution of Eq. (58a) gives us the same "envelope" law A(t)=A(0)eδt as the exact solution (10) of the initial differential equation, while the elementary integration of Eq. (58b) yields φ(t)=ξt+ φ(0)ωtω0t+φ(0). This means that our approximate solution, q(0)(t)=A(t)cos[ωtφ(t)]=A(0)eδtcos[ω0tφ(0)],
agrees with the exact Eq. (9), and misses only the correction (8) of the oscillation frequency. (This correction is of the second order in δ, i.e. of the order of ε2, and hence is beyond the accuracy of our first approximation.) It is remarkable how nicely do the reduced equations recover the proper frequency of free oscillations in this autonomous system - in which the very notion of ω is ambiguous.

The result is different at forced oscillations. For example, for the (generally, nonlinear) Duffing oscillator described by Eq. (43) with f00, Eqs. (57a) yield the reduced equations, ˙A=δA+f02ωsinφ,A˙φ=ξ(A)A+f02ωcosφ,

which are valid for an arbitrary function ξ(A), provided that this nonlinear detuning remains much smaller than the oscillation frequency. Here (after a transient), the amplitude and phase tend to the stationary states described by Eqs. (47). This means that φ becomes a constant, so that q(0)Acos(ωt const), i.e. the reduced equations again automatically recover the correct frequency of the solution, in this case equal to the external force frequency.

Note that each stationary oscillation regime, with certain amplitude and phase, corresponds to a fixed point of the reduced equations, so that the stability of those fixed points determines that of the oscillations. In the next three sections, we will carry out such an analysis for several simple systems of key importance for physics and engineering.


21 One may ask why cannot we stick to just one, most compact, complex-amplitude form (57b) of the reduced equations. The main reason is that when the function f(q,˙q,t) is nonlinear, we cannot replace its real arguments, such as q=Acos(ωtφ), with their complex-function representations like aexp{iωt} (as could be done in the linear problems considered in Sec. 5.1), and need to use real variables, such as either {A,φ} or {u,v}, anyway.


This page titled 5.3: Reduced Equations 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?