4.2: Weak Nonlinearity

Most physical oscillators become non-linear with increase in amplitude of the oscillations. Consequences of non-linearity include breakdown of superposition, introduction of additional harmonics, and complicated chaotic motion that has great sensitivity to the initial conditions as illustrated in this chapter. Weak non-linearity is interesting since perturbation theory can be used to solve the non-linear equations of motion.

The potential energy function for a linear oscillator has a pure parabolic shape about the minimum location, that is, \(U= \frac{1}{2}k(x-x_{0})^{2}\) where \(x_{0}\) is the location of the minimum. Weak non-linear systems have small amplitude oscillations \(\Delta x\) about the minimum allowing use of the Taylor expansion

\[U(\Delta x)=U(x_{0})+\Delta x\frac{dU\left( x_{0}\right) }{dx}+\frac{\Delta x^{2}}{2!}\frac{d^{2}U\left( x_{0}\right) }{dx^{2}}+\frac{\Delta x^{3}}{3!} \frac{d^{3}U\left( x_{0}\right) }{dx^{3}}+\frac{\Delta x^{4}}{4!}\frac{ d^{4}U\left( x_{0}\right) }{dx^{4}}+ \dots \label{4.1}\]

By definition, at the minimum \(\frac{dU\left( x_{0}\right) }{dx}=0,\) and thus Equation \ref{4.1} can be written as

\[\Delta U=U(\Delta x)-U(x_{0})=\frac{\Delta x^{2}}{2!}\frac{d^{2}U\left( x_{0}\right) }{dx^{2}}+\frac{\Delta x^{3}}{3!}\frac{d^{3}U\left( x_{0}\right) }{dx^{3}}+\frac{\Delta x^{4}}{4!}\frac{d^{4}U\left( x_{0}\right) }{dx^{4}}+\dots \label{4.2}\]

For small amplitude oscillations the system is linear when only the second-order \(\frac{\Delta x^{2}}{2!}\frac{d^{2}U\left( x_{0}\right) }{dx^{2} }\) term in Equation \ref{4.2} is significant. The linearity for small amplitude oscillations greatly simplifies description of the oscillatory motion in that superposition applies, and complicated chaotic motion is avoided. For slightly larger amplitude motion, where the higher-order terms in the expansion are still much smaller than the second-order term, then perturbation theory can be used as illustrated by the simple plane pendulum which is non linear since the restoring force equals

\[mg\sin \theta \simeq mg(\theta -\frac{\theta ^{3}}{3!}+\frac{\theta ^{5}}{5!} -\frac{\theta ^{7}}{7!}+\dots ) \label{4.3}\]

This is linear only at very small angles where the higher-order terms in the expansion can be neglected. Consider the equation of motion at small amplitudes for the harmonically-driven, linearly-damped plane pendulum

\[\ddot{\theta}+\Gamma \dot{\theta}+\omega _{0}^{2}\sin \theta =\ddot{\theta} +\Gamma \dot{\theta}+\omega _{0}^{2}(\theta -\frac{\theta ^{3}}{6} )=F_{0}\cos \left( \omega t\right) \label{4.4}\]

where only the first two terms in the expansion \ref{4.3} have been included. It was shown in chapter \(3\) that when \(\sin \theta \approx \theta\) then the steady-state solution of Equation \ref{4.4} is of the form

\[\theta \left( t\right) =A\cos \left( \omega t-\delta \right) \label{4.5}\]

Insert this first-order solution into Equation \ref{4.4}, then the cubic term in the expansion gives a term \(\cos^{3}\omega t=\frac{1}{4}(\cos 3\omega t+3\cos \omega t)\). Thus the perturbation expansion to third order involves a solution of the form

\[\theta \left( t\right) =A\cos \left( \omega t-\delta \right) +B\cos 3(\omega t-\delta )\label{4.6}\]

This perturbation solution shows that the non-linear term has distorted the signal by addition of the third harmonic of the driving frequency with an amplitude that depends sensitively on \(\theta\). This illustrates that the superposition principle is not obeyed for this non-linear system, but, if the non-linearity is weak, perturbation theory can be used to derive the solution of a non-linear equation of motion.

Figure \(\PageIndex{1}\) illustrates that for a potential \(U(x)=2x^{2}+x^{4},\) the \( x^{4}\) non-linear term are greatest at the maximum amplitude \(x,\) which makes the total energy contours in state-space more rectangular than the elliptical shape for the harmonic oscillator as shown in figure (3.4.2). The solution is of the form given in Equation \ref{4.6}.

Figure \(\PageIndex{1}\): The left side shows the potential energy for a symmetric potential \(U(x)=2x^2 + x^4\). The right side shows the contours of constant total energy on a state-space diagram.

Example \(\PageIndex{1}\): Non-linear oscillator

Assume that a non-linear oscillator has a potential given by

\[U(x)=\frac{kx^{2}}{2}-\frac{m\lambda x^{3}}{3} \nonumber\]

where \(\lambda\) is small. Find the solution of the equation of motion to first order in \(\lambda\), assuming \(x=0\) at \(t=0\).

Solution

The equation of motion for the nonlinear oscillator is

\[m\ddot{x}=-\frac{dU}{dx}=-kx+m\lambda x^{2} \nonumber\]

If the \(m\lambda x^{2}\) term is neglected, then the second-order equation of motion reduces to a normal linear oscillator with

\[x_{0}=A\sin \left( \omega _{0}t+\varphi \right) \nonumber\]

where

\[\omega _{0}=\sqrt{\frac{k}{m}} \nonumber\]

Assume that the first-order solution has the form

\[x_{1}=x_{0}+\lambda x_{1} \nonumber\]

Substituting this into the equation of motion, and neglecting terms of higher order than \(\lambda ,\) gives

\[\ddot{x}_{1}+\omega _{0}^{2}x_{1}=x_{0}^{2}=\frac{A^{2}}{2}[1-\cos \left( 2\omega _{0}t\right) ] \nonumber\]

To solve this try a particular integral

\[x_{1}=B+C\cos \left( 2\omega _{0}t\right) \nonumber\]

and substitute into the equation of motion gives

\[-3\omega _{0}^{2}C\cos \left( 2\omega _{0}t\right) +\omega _{0}^{2}B=\frac{ A^{2}}{2}-\frac{A^{2}}{2}\cos \left( 2\omega _{0}t\right) \nonumber\]

Comparison of the coefficients gives

\[\begin{aligned} B &=&\frac{A^{2}}{2\omega _{0}^{2}} \\ C &=&\frac{A^{2}}{6\omega _{0}^{2}}\end{aligned} \nonumber\]

The homogeneous equation is

\[\ddot{x}_{1}+\omega _{0}^{2}x_{1}=0 \nonumber\]

which has a solution of the form

\[x_{1}=D_{1}\sin \left( \omega _{0}t\right) +D_{2}\cos \left( \omega _{0}t\right) \nonumber\]

Thus combining the particular and homogeneous solutions gives

\[x_{1}=\left( A+\lambda D_{1}\right) \sin \left( \omega _{0}t\right) +\lambda \left[ \frac{A^{2}}{2\omega _{0}^{2}}+D_{2}\cos \left( \omega _{0}t\right) + \frac{A^{2}}{6\omega _{0}^{2}}\cos \left( 2\omega _{0}t\right) \right] \nonumber\]

The initial condition \(x=0\) at \(t=0\) then gives

\[D_{2}=-\frac{2A^{2}}{3\omega ^{2}} \nonumber\]

and

\[x_{1}=\left( A+\lambda D_{1}\right) \sin \left( \omega _{0}t\right) +\frac{ \lambda A^{2}}{\omega _{0}^{2}}\left[ \frac{1}{2}-\frac{2}{3}\cos \left( \omega _{0}t\right) +\frac{1}{6}\cos \left( 2\omega _{0}t\right) \right] \nonumber\]

The constant \(\left( A+\lambda D_{1}\right)\) is given by the initial amplitude and velocity.

This system is nonlinear in that the output amplitude is not proportional to the input amplitude. Secondly, a large amplitude second harmonic component is introduced in the output waveform; that is, for a non-linear system the gain and frequency decomposition of the output differs from the input. Note that the frequency composition is amplitude dependent. This particular example of a nonlinear system does not exhibit chaos. The Laboratory for Laser Energetics uses nonlinear crystals to double the frequency of laser light.