Analytic Solution

Before attempting a numerical solution of the equations of motion of any dynamical system, it is a good idea to, first, investigate the equations as thoroughly as possible via standard analytic techniques. Unfortunately, Equations (1237) and (1238) constitute a non-linear dynamical system--because of the presence of the $\sin\theta$ term on the right-hand side of Equation (1238). This system, like most non-linear systems, does not possess a simple analytic solution. Fortunately, however, if we restrict our attention to small amplitude oscillations, such that the approximation

$$\sin\theta \simeq \theta$$ (1239)

is valid, then the system becomes linear, and can easily be solved analytically.

The linearized equations of motion of the pendulum take the form:

$$\frac{d\theta}{dt} = v,$$ (1240)

$$\frac{dv}{dt} = -\frac{v}{Q} -\theta + A\,\cos(\omega \,t).$$ (1241)

Suppose that the pendulum's position, $\theta (0)$, and velocity, $v(0)$, are specified at time $t=0$. As is well-known, in this case, the above equations of motion can be solved analytically to give:

$$\theta(t) = \left\{\theta(0) - \frac{A\,(1-\omega^2)}{[(1-\omega^2)^2+\omega^2/Q^2]}\right\} {\rm e}^{-t/2Q}\,\cos(\omega_\ast\, t)$$

$$+\frac{1}{\omega_\ast}\left\{ v(0) + \frac{\theta(0)}{2Q} -\frac{... ...} {[(1-\omega^2)^2+ \omega^2/Q^2]}\right\}{\rm e}^{-t/2Q}\,\sin(\omega_\ast\,t)$$

$$+\frac{A\left[(1-\omega^2)\,\cos(\omega\, t)+ (\omega/Q)\,\sin(\omega\, t)\right]} {\left[(1-\omega^2)^2+\omega^2/Q^2\right]},$$ (1242)

$$v(t) = \left\{v(0) - \frac{A\,\omega^2/Q}{[(1-\omega^2)^2+\omega^2/Q^2]}\right\} {\rm e}^{-t/2Q}\,\cos\omega_\ast t$$

$$-\frac{1}{\omega_\ast}\left\{ \theta(0) + \frac{v(0)}{2Q} -\frac{... ... {[(1-\omega^2)^2+ \omega^2/Q^2]}\right\}{\rm e}^{-t/2Q}\,\sin(\omega_\ast\, t)$$

$$+\frac{\omega\,A\left[-(1-\omega^2)\,\sin(\omega\, t)+ (\omega/Q)\,\cos(\omega\, t)\right]} {\left[(1-\omega^2)^2+\omega^2/Q^2\right]}.$$ (1243)

Here,

$$\omega_\ast = \sqrt{1-\frac{1}{4Q^2}},$$ (1244)

and it is assumed that $Q>1/2$. It can be seen that the above expressions for $\theta$ and $v$ both consist of three terms. The first two terms clearly represent transients, since they depend on the initial conditions, and decay exponentially in time--see Section 3.9. In fact, the e-folding time for the decay of these terms is $2\,Q$ (in normalized time units). The final term represents the time-asymptotic motion of the pendulum, and is manifestly independent of the initial conditions--see Section 3.5.

Figure 58: A phase-space plot of the periodic attractor for a linear damped periodically driven pendulum. Data calculated analytically for $Q=4$ and $\omega =2$.

It is often convenient to visualize the motion of a dynamical system as an orbit, or trajectory, in phase-space, which is defined as the space of all of the dynamical variables required to specify the instantaneous state of the system. For the case in hand, there are two dynamical variables, $v$ and $\theta$, and so phase-space corresponds to the $\theta$-$v$ plane. Note that each different point in this plane corresponds to a unique instantaneous state of the pendulum. [Strictly speaking, we should also consider $t$ to be a dynamical variable, since it appears explicitly on the right-hand side of Equation (1238).]

It is clear, from Equations (1242) and (1243), that if we wait long enough for all of the transients to decay away then the motion of the pendulum settles down to the following simple orbit in phase-space:

$$\theta(t) = \frac{A\left[(1-\omega^2)\,\cos(\omega\, t)+ (\omega/Q)\,\sin(\omega\, t)\right]} {\left[(1-\omega^2)^2+\omega^2/Q^2\right]},$$ (1245)

$$v(t) = \frac{\omega\, A\left[-(1-\omega^2)\,\sin(\omega\, t)+ (\omega/Q)\,\cos(\omega\, t)\right]} {\left[(1-\omega^2)^2+\omega^2/Q^2\right]}.$$ (1246)

This orbit traces out the closed curve

$$\left(\frac{\theta}{\tilde{A}}\right)^2 + \left(\frac{v}{\omega\,\tilde{A}}\right)^2 =1,$$ (1247)

in phase-space, where

$$\tilde{A} = \frac{A}{\sqrt{(1-\omega^2)^2+\omega^2/Q^2}}.$$ (1248)

As illustrated in Figure 58, this curve is an ellipse whose principal axes are aligned with the $v$ and $\theta$ coordinate axes. Observe that the curve is closed, which suggests that the associated motion is periodic in time. In fact, the motion repeats itself exactly every

$$\tau = \frac{2\pi}{\omega}$$ (1249)

normalized time units. The maximum angular displacement of the pendulum from its undriven rest position ($\theta=0$) is $\tilde{A}$. As illustrated in Figure 59, the variation of $\tilde{A}$ with driving frequency $\omega$ [see Equation (1248)] displays all of the features of a classic resonance curve. The maximum amplitude of the driven oscillation is proportional to the quality-factor, $Q$, and is achieved when the driving frequency matches the natural frequency of the pendulum (i.e., when $\vert\omega\vert=1$). Moreover, the width of the resonance in $\omega$-space is proportional to $1/Q$--see Section 3.7.

Figure 59: The maximum angular displacement of a linear damped periodically driven pendulum as a function of driving frequency. The solid curve corresponds to $Q=1$. The short-dashed curve corresponds to $Q=5$. The long-dashed curve corresponds to $Q=10$. Analytic data.

The phase-space curve shown in Figure 58 is called a periodic attractor. It is termed an ``attractor'' because, irrespective of the initial conditions, the trajectory of the system in phase-space tends asymptotically to--in other words, is attracted to--this curve as $t\rightarrow \infty$. This gravitation of phase-space trajectories towards the attractor is illustrated in Figures 60 and 61. Of course, the attractor is termed ``periodic'' because it corresponds to motion which is periodic in time.

Figure 60: The phase-space trajectory of a linear damped periodically driven pendulum. Data calculated analytically for $Q=1$ and $\omega =2$. Here, $v(0)/A=0$ and $\theta (0)/A=0$.

Figure 61: The phase-space trajectory of a linear damped periodically driven pendulum. Data calculated analytically for $Q=1$ and $\omega =2$. Here, $v(0)/A=0.5$ and $\theta (0)/A=0.5$.

Let us summarize our findings so far. We have discovered that if a damped pendulum is subject to a low amplitude periodic drive then its time-asymptotic response (i.e., its response after any transients have died away) is periodic, with the same period as the driving torque. Moreover, the response exhibits resonant behaviour as the driving frequency approaches the natural frequency of oscillation of the pendulum. The amplitude of the resonant response, as well as the width of the resonant window, is governed by the amount of damping in the system. After a little reflection, we can easily appreciate that all of these results are a direct consequence of the linearity of the pendulum's equations of motion in the low amplitude limit. In fact, it is easily demonstrated that the time-asymptotic response of any intrinsically stable linear system (with a discrete spectrum of normal modes) to a periodic drive is periodic, with the same period as the drive. Moreover, if the driving frequency approaches one of the natural frequencies of oscillation of the system then the response exhibits resonant behaviour. But, is this the only allowable time-asymptotic response of a dynamical system to a periodic drive? It turns out that it is not. Indeed, the response of a non-linear system to a periodic drive is generally much richer and far more diverse than simple periodic motion. Since the majority of naturally occurring dynamical systems are non-linear, it is clearly important that we gain a basic understanding of this phenomenon. Unfortunately, we cannot achieve this goal via a standard analytic approach, since non-linear equations of motion generally do not possess simple analytic solutions. Instead, we must employ numerical methods. As an example, let us investigate the dynamics of a damped pendulum, subject to a periodic drive, with no restrictions on the amplitude of the pendulum's motion.