Solving a Second-Order Differential Equation

  • Thread starter Saketh
  • Start date
  • #1
261
2
I've been trying to derive the period of a physical pendulum, and I ended up with a differential equation that boils down to this:

[tex]
\frac{d^2 \theta}{dt^2} = C \sin{\theta}
[/tex]

where C is some constant.

With no experience in differential equations, I have no idea how to solve this. I can't find anything about how to solve this type of differential equation, so if someone could point me in the right direction, that would be helpful.

Thanks!
 

Answers and Replies

  • #2
StatusX
Homework Helper
2,571
1
It can't be solved nicely, which is why we usually assume the angle is small and make the approximation [itex] \sin \theta \approx \theta[/itex]. If you want, there's a way to reduce it to finding the integral of something like [itex]1/\sqrt{a-\cos(\theta)}[/tex], but that's not a nice integral to do (ie, you need special functions).
 
  • #3
HallsofIvy
Science Advisor
Homework Helper
41,847
965
That is an extremly non-linear equation and cannot be solved in any simple form. As StatusX said, for smal values of the angle, you can make the approximation [itex]sin(\theta)= \theta[/itex] so the equation becomes [itex]d^2\theta/dt^2= C\theta[/itex].

Another technique is called "quadrature". Since t does not appear explictly in the equation, let [itex]\omega= d\theta/dt[/itex]. Then apply the chain rule: [itex]d^2\theta/dt^2= d\omega/dt= d\omega/d\theta d\theta/dt= \omega d\omega/dt= Csin \theta[/itex]. That's a separable first order differential equation which can be written as [itex]\omega d\omega= Csin(\theta)d\theta[/itex] and can be integrated directly:
[itex](1/2) \omega^2= -Ccos(\theta)+ D[/itex]. The rub is when you replace [itex]\omega[/itex] by [itex]d/t\eta/dt[/itex] and solve for [itex]d\theta/dt[/itex]: [itex] d\theta/dt= \sqrt{D -2C cos\theta}[/itex] so
[itex]\frac{d\theta}{\sqrt{D- 2C cos(\theta}}= dt[/itex] and the left side cannot be integrated in terms of any elementary function. (That's called an "elliptic" integral. I've seen whole book cases of tables of values.)
 
  • #4
261
2
Does this mean that it is impossible to generate a theoretical formulation for pendulums that includes large-angle values?

I suspected that I would have to make a small-angle assumption. If I did have [itex]d^2\theta/dt^2= C\theta[/itex], how would I solve that? (Sorry for my lack of experience - once again, I just need to be pointed in the right diretcion.)
 
  • #5
HallsofIvy
Science Advisor
Homework Helper
41,847
965
For small angles, yes. Of course, in order to get periodic solutions, C must be negative.
 
  • #6
daniel_i_l
Gold Member
867
0
what function equals itself times a constant when differentiating twice?
(HINT: whats the derivative of e^x ?)
 
  • #8
arildno
Science Advisor
Homework Helper
Gold Member
Dearly Missed
10,025
134
Saketh said:
Does this mean that it is impossible to generate a theoretical formulation for pendulums that includes large-angle values?

I suspected that I would have to make a small-angle assumption. If I did have [itex]d^2\theta/dt^2= C\theta[/itex], how would I solve that? (Sorry for my lack of experience - once again, I just need to be pointed in the right diretcion.)
You can get an approximate solution to this, to whichever degree of accuracy you desire, by perturbation theory.

Alternatively, you may use numerics
 
  • #9
261
2
daniel_i_l said:
what function equals itself times a constant when differentiating twice?
(HINT: whats the derivative of e^x ?)
I don't understand - how does that help me solve the differential equation?

Thanks, everyone! http://scienceworld.wolfram.com/physics/Pendulum.html" [Broken] answers all of my questions.
 
Last edited by a moderator:
  • #10
HallsofIvy
Science Advisor
Homework Helper
41,847
965
daniel_i_l said:
what function equals itself times a constant when differentiating twice?
(HINT: whats the derivative of e^x ?)
I can think of 4 such (independent) functions. That's why I said " Of course, in order to get periodic solutions, C must be negative."
 
  • #11
682
1
just a curious question, has there ever been a series solution to the exact motion of undamped pendulum?
 

Related Threads on Solving a Second-Order Differential Equation

Replies
33
Views
1K
Replies
12
Views
3K
Replies
1
Views
1K
Replies
5
Views
2K
Replies
6
Views
2K
Replies
14
Views
3K
Replies
3
Views
8K
Replies
2
Views
4K
Replies
7
Views
487
Replies
1
Views
2K
Top