How to find an equilibrium of the second order system?

I have following exercise to prepare to my exam : #x^{..}=1-3x^2# (i mean second derivative of x by #x^{..}#). And i need to "find equilibria of the system" I know that for first order it is enough to check points where x'=0, but how to extend it to second order if we dont have x' at all? Would be very grateful if someone could clarify it and help to solve the exercise

1 Answer
Jun 23, 2018

See below

Explanation:

#ddot x=1-3x^2# is non-linear in #t# so you either change variable or use a different analysis.

You might know that every conservative system has a potential function #V(x)#, where:

  • # ddot x = - (dV)/(dx)#

If you don't know what a potential function is, I've put a simple doodle at the bottom. They are ubiquitous in physical applications.

So here:

#V(x) =- int 1-3x^2 \ dx= - x + x^3 + V_o #

We don't really care what #V_o# is, we are interested in the shape of the potential function, so we plot it, with #V_o = 0#:

graph{-x + x^3 [-2.738, 2.738, -1.37, 1.366]}

This potential function has 2 equilibrium points corresponding to the stationary points where #(dV)/(dx) = 0#, ie #x = +- 1/sqrt3#

The one on the left is unstable . A particle placed at #x = - 1/sqrt3#, if it underwent a slight displacement, would "roll off the hill" in the direction of the displacement. You can look at it like a picture.

A particle placed at #x = 1/sqrt3#, if it underwent a slight displacement, would oscillate back and forward about #x = 1/sqrt3#. That is stable equilibrium.

There is more, and you can actually solve that DE too, but it does not give the same information.

Potential Function rough idea

In terms of the potential function, we have a conservative sytem so energy #E# is conserved, though it may switch between kinetic energy #T# and potential energy #V#.

For a unit mass:

#E = T + V = const#

#(dE)/(dx) = 0#

#d/(dx) (1/2 v^2) + (dV)/(dx) = 0#

# v(dv)/(dx) + (dV)/(dx) = 0#

  • #v (dv)/(dx) = v (dv )/(dt)(dt)/(dx) = v (dv )/(dt)1/v = (dv)/(dt) = ddot x#

#implies ddot x = - (dV)/(dx)#