Physics:Kapitza's pendulum

Drawing showing how a Kapitza pendulum can be constructed: a motor rotates a crank at a high speed, the crank vibrates a lever arm up and down, which the pendulum is attached to with a pivot.

Kapitza's pendulum or Kapitza pendulum is a rigid pendulum in which the pivot point vibrates in a vertical direction, up and down. It is named after Russian Nobel laureate physicist Pyotr Kapitza, who in 1951 developed a theory which successfully explains some of its unusual properties. The unique feature of the Kapitza pendulum is that the vibrating suspension can cause it to balance stably in an inverted position, with the bob above the suspension point. In the usual pendulum with a fixed suspension, the only stable equilibrium position is with the bob hanging below the suspension point; the inverted position is a point of unstable equilibrium, and the smallest perturbation moves the pendulum out of equilibrium. In nonlinear control theory the Kapitza pendulum is used as an example of a parametric oscillator that demonstrates the concept of "dynamic stabilization".

The pendulum was first described by A. Stephenson in 1908, who found that the upper vertical position of the pendulum might be stable when the driving frequency is fast. Yet until the 1950s there was no explanation for this highly unusual and counterintuitive phenomenon. Pyotr Kapitza was the first to analyze it in 1951. He carried out a number of experimental studies and as well provided an analytical insight into the reasons of stability by splitting the motion into "fast" and "slow" variables and by introducing an effective potential. This innovative work created a new subject in physics – vibrational mechanics. Kapitza's method is used for description of periodic processes in atomic physics, plasma physics and cybernetical physics. The effective potential which describes the "slow" component of motion is described in "Mechanics" volume (§30) of Landau's Course of Theoretical Physics.

Another interesting feature of the Kapitza pendulum system is that the bottom equilibrium position, with the pendulum hanging down below the pivot, is no longer stable. Any tiny deviation from the vertical increases in amplitude with time. Parametric resonance can also occur in this position, and chaotic regimes can be realized in the system when strange attractors are present in the Poincaré section.

Notation

Kapitza's pendulum scheme

Denote the vertical axis as [math]\displaystyle{ y }[/math] and the horizontal axis as [math]\displaystyle{ x }[/math] so that the motion of pendulum happens in the ([math]\displaystyle{ x }[/math]-[math]\displaystyle{ y }[/math]) plane. The following notation will be used

• [math]\displaystyle{ \nu }[/math]—frequency of the vertical oscillations of the suspension,
• [math]\displaystyle{ a }[/math] — amplitude of the oscillations of the suspension,
• [math]\displaystyle{ \omega_0 = \sqrt{g/l} }[/math] — proper frequency of the mathematical pendulum,
• [math]\displaystyle{ g }[/math] — free fall acceleration,
• [math]\displaystyle{ l }[/math] — length of rigid and light pendulum,
• [math]\displaystyle{ m }[/math] — mass.

Denoting the angle between pendulum and downward direction as [math]\displaystyle{ \varphi }[/math] the time dependence of the position of pendulum gets written as

[math]\displaystyle{ \begin{cases} x &= l \sin \varphi\\ y &= - l \cos \varphi - a \cos \nu t \end{cases} }[/math]

Energy

The potential energy of the pendulum is due to gravity and is defined by, in terms of the vertical position, as

[math]\displaystyle{ E_\mathrm{POT} = - m g (l \cos \varphi + a \cos \nu t). \, }[/math]

The kinetic energy in addition to the standard term [math]\displaystyle{ E_\mathrm{KIN}=m l^2 \dot \varphi^2 /2 }[/math], describing velocity of a mathematical pendulum, there is a contribution due to vibrations of the suspension

[math]\displaystyle{ E_\mathrm{KIN} = \frac{m l^2 }{2} \dot \varphi^2 + m a l \nu ~\sin(\nu t) \sin(\varphi)~\dot\varphi + \frac{m a^2 \nu^2}{2} \sin^2(\nu t)\;. }[/math]

The total energy is given by the sum of the kinetic and potential energies [math]\displaystyle{ E = E_\mathrm{KIN} + E_\mathrm{POT} }[/math] and the Lagrangian by their difference [math]\displaystyle{ L = E_\mathrm{KIN} - E_\mathrm{POT} }[/math].

The total energy is conserved in a mathematical pendulum, so time [math]\displaystyle{ t }[/math] dependence of the potential [math]\displaystyle{ E_\mathrm{POT} }[/math] and kinetic [math]\displaystyle{ E_\mathrm{KIN} }[/math] energies is symmetric with respect to the horizontal line. According to the virial theorem the mean kinetic and potential energies in harmonic oscillator are equal. This means that the line of symmetry corresponds to half of the total energy.

In the case of vibrating suspension the system is no longer a closed one and the total energy is no longer conserved. The kinetic energy is more sensitive to vibration compared to the potential one. The potential energy [math]\displaystyle{ E_\mathrm{POT} = mgy }[/math] is bound from below and above [math]\displaystyle{ -mg(l+a)\lt E_\mathrm{POT}\lt mg(l+a) }[/math] while the kinetic energy is bound only from below [math]\displaystyle{ E_\mathrm{KIN}\ge 0 }[/math]. For high frequency of vibrations [math]\displaystyle{ \nu }[/math] the kinetic energy can be large compared to the potential energy.

Equations of motion

Motion of pendulum satisfies Euler–Lagrange equations. The dependence of the phase [math]\displaystyle{ \varphi }[/math] of the pendulum on its position satisfies the equation:

[math]\displaystyle{ \frac{d}{dt} \frac{\partial L}{\partial \dot \varphi} = \frac{\partial L}{\partial \varphi}, }[/math]

where the Lagrangian [math]\displaystyle{ L }[/math] reads

[math]\displaystyle{ L = \frac{m l^2 }{2} \dot \varphi^2 + ml( g + a~\nu^2\cos\nu t) \cos \varphi, }[/math]

up to irrelevant total time derivative terms. The differential equation

[math]\displaystyle{ \ddot \varphi = - (g+a~\nu^2\cos\nu t) \frac{\sin \varphi}{l}, }[/math]

which describes the movement of the pendulum is nonlinear due to the [math]\displaystyle{ \sin\varphi }[/math] factor.

Equilibrium positions

Kapitza's pendulum model is more general than the simple pendulum. The Kapitza model reduces to the latter in the limit [math]\displaystyle{ a = 0 }[/math]. In that limit, the tip of the pendulum describes a circle: [math]\displaystyle{ x^2+y^2 = l^2 = \text{constant} }[/math]. If the energy in the initial moment is larger than the maximum of the potential energy [math]\displaystyle{ E \gt mgl }[/math] then the trajectory will be closed and cyclic. If the initial energy is smaller [math]\displaystyle{ E \lt mgl }[/math] then the pendulum will oscillate close to the only stable point [math]\displaystyle{ (x,y) = (0,-l) }[/math].

When the suspension is vibrating with a small amplitude [math]\displaystyle{ a \ll l }[/math] and with a frequency [math]\displaystyle{ \nu\gg \omega_0 }[/math] much higher than the proper frequency [math]\displaystyle{ \omega_0 }[/math], the angle [math]\displaystyle{ \varphi }[/math] may be viewed as a superposition [math]\displaystyle{ \varphi=\varphi_0+\xi }[/math] of a "slow" component [math]\displaystyle{ \varphi_0 }[/math] and a rapid oscillation [math]\displaystyle{ \xi }[/math] with small amplitude due to the small but rapid vibrations of the suspension. Technically, we perform a perturbative expansion in the "coupling constants" [math]\displaystyle{ (a/l),(\omega_0/\nu) \ll 1 }[/math] while treating the ratio [math]\displaystyle{ (a/l)(\nu/\omega_0) }[/math] as fixed. The perturbative treatment becomes exact in the double scaling limit [math]\displaystyle{ a \to 0 ,\nu\to \infty }[/math]. More precisely, the rapid oscillation [math]\displaystyle{ \xi }[/math] is defined as

[math]\displaystyle{ \xi = \frac{a}{l}\sin\varphi_0 ~\cos\nu t. \, }[/math]

The equation of motion for the "slow" component [math]\displaystyle{ \varphi_0 }[/math] becomes

[math]\displaystyle{ \begin{align} \ddot\varphi_0 = \ddot\varphi - \ddot\xi &= -(g+a~\nu^2\cos\nu t)\frac{\sin\varphi}{l} \\ & {}\quad{} - \frac{a}{l}\left(\ddot\varphi_0 \cos \varphi_0 ~\cos\nu t -\dot\varphi_0^2\sin\varphi_0 ~\cos\nu t - 2\nu\dot\varphi_0\cos\varphi_0~\sin\nu t - \nu^2\sin \varphi_0 ~\cos\nu t \right) \\[8pt] &= -\frac{g}{l}\sin\varphi_0 -(g+a~\nu^2\cos\nu t) \frac{1}{l}\left(\xi\cos\varphi_0 + O(\xi^2)\right) \\ & {}\quad{} - \frac{a}{l}\left( \ddot\varphi_0 \cos\varphi_0 ~\cos\nu t -\dot\varphi_0^2\sin\varphi_0 ~\cos\nu t - 2\nu\dot\varphi_0\cos\varphi_0~\sin\nu t \right). \end{align} }[/math]

Time-averaging over the rapid [math]\displaystyle{ \nu }[/math]-oscillation yields to leading order

[math]\displaystyle{ \ddot \varphi_0 = - \frac{g}{l}\sin \varphi_0 - \frac{1}{2}\left(\frac{a\nu}{l}\right)^2\sin \varphi_0 \cos \varphi_0. }[/math]

The "slow" equation of motion becomes

[math]\displaystyle{ m l^2\ddot \varphi_0 = -\frac{\partial V_{\mathrm{eff}}}{\partial \varphi_0} \;, }[/math]

by introducing an effective potential

[math]\displaystyle{ V_{\mathrm{eff}} = - mgl \cos \varphi_0 + m \left(\frac{a\nu}{2}\sin \varphi_0\right)^2. }[/math]

It turns out that the effective potential [math]\displaystyle{ V_{\mathrm{eff}} }[/math] has two minima if [math]\displaystyle{ (a\nu)^2 \gt 2gl }[/math], or equivalently, [math]\displaystyle{ (a/l)(\nu/\omega_0) \gt \sqrt{2} }[/math]. The first minimum is in the same position [math]\displaystyle{ (x,y)=(0,-l) }[/math] as the mathematical pendulum and the other minimum is in the upper vertical position [math]\displaystyle{ (x,y)=(0,l) }[/math]. As a result the upper vertical position, which is unstable in a mathematical pendulum, can become stable in Kapitza's pendulum.

Bifurcation and chaos

For fixed forcing frequency [math]\displaystyle{ \nu }[/math], and increasing forcing amplitude [math]\displaystyle{ a }[/math], the inverted point [math]\displaystyle{ \varphi = \pi }[/math] becomes a stable equilibrium, and then it would become unstable again. The point undergoes a Hopf bifurcation: the stable equilibrium expands into a stable oscillation around [math]\displaystyle{ \varphi = \pi }[/math], at frequency [math]\displaystyle{ \nu /2 }[/math], and amplitude [math]\displaystyle{ \propto \sqrt{a - a_c} }[/math] where [math]\displaystyle{ a_c }[/math] is the critical amplitude.

For even larger amplitude, the stable oscillation becomes unstable, and the pendulum would start rotating. At even larger amplitudes, the rotating mode is also destroyed, and a strange attractor appears by period-doubling cascade.

Rotating solutions

The rotating solutions of the Kapitza's pendulum occur when the pendulum rotates around the pivot point at the same frequency that the pivot point is driven. There are two rotating solutions, one for a rotation in each direction. We shift to the rotating reference frame using [math]\displaystyle{ \varphi \rightarrow \varphi^{\prime} \pm \nu t }[/math] and the equation for [math]\displaystyle{ \varphi }[/math] becomes:

[math]\displaystyle{ \ddot \varphi^{\prime} = -\frac{1}{l} \left[ \frac{1}{2}a \nu^{2} \sin(\varphi^{\prime}) + g\sin(\varphi^{\prime} \pm \nu t) + \frac{1}{2}a \nu^{2} \sin(\varphi^{\prime} \pm 2\nu t) \right] \;. }[/math]

Again considering the limit in which [math]\displaystyle{ \nu }[/math] is much higher than the proper frequency [math]\displaystyle{ \omega_0 }[/math], we find that the rapid-[math]\displaystyle{ \nu }[/math] slow-[math]\displaystyle{ \varphi_0^{\prime} }[/math] limit leads to the equation:

[math]\displaystyle{ \ddot \varphi_0^{\prime} = -\frac{1}{2l} a \nu^2 \sin\varphi_0^{\prime} \;. }[/math]

The effective potential is just that of a simple pendulum equation. There is a stable equilibrium at [math]\displaystyle{ \varphi_0^{\prime} = 0 }[/math] and an unstable equilibrium at [math]\displaystyle{ \varphi_0^{\prime} = \pi }[/math].

Phase portrait

Interesting phase portraits might be obtained in regimes which are not accessible within analytic descriptions, for example in the case of large amplitude of the suspension [math]\displaystyle{ a \approx l }[/math]. Increasing the amplitude of driving oscillations to half of the pendulum length [math]\displaystyle{ a = l/2 }[/math] leads to the phase portrait shown in the figure.

Further increase of the amplitude to [math]\displaystyle{ a\approx l }[/math] leads to full filling of the internal points of the phase space: if before some points of the phase space were not accessible, now system can reach any of the internal points. This situation holds also for larger values of [math]\displaystyle{ a }[/math].

Interesting facts

• Kapitza noted that a pendulum clock with a vibrating pendulum suspension always goes faster than a clock with a fixed suspension.
• Walking is defined by an 'inverted pendulum' gait in which the body vaults over the stiff limb or limbs with each step. Increased stability during walking might be related to stability of Kapitza's pendulum. This applies regardless of the number of limbs - even arthropods with six, eight or more limbs.

External links

• Demonstration video at Kapitza's Pendulum - YouTube
• Interactive demonstration at Wolfram Demonstrations Project