• Support PF! Buy your school textbooks, materials and every day products Here!

Pendulum differential equation question

  • Thread starter okalakacheekee
  • Start date
  • #1
okalakacheekee
Alright, for the second order differential d^2 theta/dt^2 = -g/l sin theta, where l is length of pendulum, g is gravity, etc...how do you solve that exactly for theta as a function of time?

I substituted (w=omega) dw/dt for d2theta/dt2...and eventually got 1/2 w^2 = g/l cos theta + C...but that doesn't give me theta as a function of time?

Any help is appreciated...thanks
 

Answers and Replies

  • #2
ahrkron
Staff Emeritus
Gold Member
736
1
This problem usually asks for an approximation for small theta.

When that is the case, you can use the fact that

[tex]\sin \theta \rightarrow \theta[/tex],

[tex]\cos \theta \rightarrow 1 - \frac{\theta^2}{2}[/tex]
 
Last edited:
  • #3
okalakacheekee
but is there a way to solve it exactly?
 
  • #4
okalakacheekee
alright here's what i did
so let's say O is theta, w is omega, for ease of writing

first i have
d2O/dt2 = -g/l sin O
I said w = dO/dt and then dw/dt=d2O/dt2
eventually i got to a point where i had
w dw = -g/l sin O dO
so i integrated and got
1/2 w^2 = g/l cos O + C

problem being, i had no time in there...
so then i put back in dO/dt for w and got
dO/dt = sqrt(2g/l * cosO)
so when you separate everything you get
dO/sqrt (cosO) = sqrt(2g/l) dO
then you integrate..here's where i ran into trouble yet again...how do you integrate the left hand side?
 
  • #5
HallsofIvy
Science Advisor
Homework Helper
41,770
911
No, there is no way to solve the "pendulum problem" exactly.

A standard attack is "linearization"- for small values of θ, replace sin(θ) by its linear approximation θ to get the linear equation d2θ/dt2= -(g/l)θ.

Another method is "quadrature" which is basically what you are doing. Let ω= dθ/dt so that d2θ/dt2= dω/dt= (dω/dθ)(dθ/dt)= &omega dω/dt= -(g/l)sinθ.

That can be integrated to get (1/2)ω2= (g/l)cos(θ)+ C, a "first integral" (which physicists would associate with the "total energy" of the pendulum).

You could, of course, rewrite that as ω= dθ/dt= √((2g/l)cos(θ)+ C); but the resulting integral is an "elliptic integral" which cannot be integrated exactly.
 
  • #6
krab
Science Advisor
896
1
Originally posted by okalakacheekee
alright here's what i did
so let's say O is theta, w is omega, for ease of writing

first i have
d2O/dt2 = -g/l sin O
I said w = dO/dt and then dw/dt=d2O/dt2
eventually i got to a point where i had
w dw = -g/l sin O dO
so i integrated and got
1/2 w^2 = g/l cos O + C

problem being, i had no time in there...
so then i put back in dO/dt for w and got
dO/dt = sqrt(2g/l * cosO)
so when you separate everything you get
dO/sqrt (cosO) = sqrt(2g/l) dO
then you integrate..here's where i ran into trouble yet again...how do you integrate the left hand side?
What you've done turns up the very useful relation of pendulum speed versus amplitude. Another way to derive it is to write KE + PE = [tex]{1\over 2}mv^2-mgl\cos\theta=\mbox{constant}[/tex].

You can plot curves in v vs. θ space (this space is called phase space); these will look like ellipses but become distorted into eye-shaped as the pendulum amplitude reaches large enough angles.

If you know enough "special functions", then yes, the problem is solvable in closed form. Look up info on "Jacobian Elliptic Functions". They're not as common as sines and cosines, but they are just as legitimate in the "finding a closed form solution" sense.
 
Last edited:

Related Threads for: Pendulum differential equation question

  • Last Post
Replies
5
Views
810
Replies
1
Views
222
  • Last Post
Replies
19
Views
6K
Replies
12
Views
2K
  • Last Post
Replies
3
Views
9K
  • Last Post
Replies
9
Views
1K
  • Last Post
Replies
8
Views
11K
  • Last Post
Replies
3
Views
8K
Top