Kapitsa Pendulum: Stability Criterion Homework

  • Thread starter Thread starter nmatthies
  • Start date Start date
  • Tags Tags
    Stability
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 4K views
nmatthies
Messages
1
Reaction score
0
Hi! I'm a maths student taking a maths by computer course and this week's homework deals with the kapitsa pendulum (reversed pendulum).

Homework Statement



Kapitsa Pendulum is a rigid-rod pendulum the suspension point of which is vibrated. It turns out that vibration can stabilize the upper equilibrium point which is unstable in absence of vibration. Here, we will study this effect using a simplified linear non-autonomous equation [tex]\ddot{x} = x(1+Asin(wt))[/tex] (1), where constants A >> 1 and w >> 1 model the amplitude and cyclic frequence of the suspension point vibrations. [tex]x_{0}=0[/tex] represents the upper equilibrium of the pendulum.

Homework Equations



Theory for Kapitsa pendulum predicts that the motion consists of 2 parts, [tex]x(t)=X(t)+ \tilde{x}(t)[/tex] (2), with fast oscillations [tex]\tilde{x}(t)= \frac {-AX(t)} {w^2} sin(wt)[/tex] (3), and a slowly varying motion X(t) which satisfies the following equation, [tex]\ddot{X}=(1-\frac {A^2} {2w^2}) X[/tex] (4).


The Attempt at a Solution



Now one of the questions is to find the analytical solution of the equation (4) and derive the stability of the point X=0 in terms of the parameters A and w.

The analytical solution I have found is [tex]X(t)=X(0)cos(\sqrt{\frac {A^2} {2w^2} -1}t)[/tex] (the t is outside of the square root). but I honestly do not know where to start for the stability criterion. What defines stability at a point? and how to express it in terms of A and w? Help will be greatly appreciated.

edit:

does it have something to do with [tex]\frac {A^2} {2w^2} -1 < 0 or \frac {A^2} {2w^2} -1 > 0[/tex]?
 
Last edited:
Physics news on Phys.org
nmatthies said:
[tex]\frac {A^2} {2w^2} -1 < 0 or \frac {A^2} {2w^2} -1 > 0[/tex]?
Bingo. You were given that the slow solution follows

[tex]\ddot X = \left(1-\frac{A^2}{2\omega^2}\right)X[/tex]

To make it easier to see what is going on, denote [itex]b\equiv 1-A^2/(2\omega^2)[/itex]. With this,

[tex]\ddot X = bX[/tex]

The solutions to this equation behave very differently if b is positive or negative.