1
$\begingroup$

I encountered the following differential equation when I tried to derive the equation of motion of a simple pendulum:

$\frac{\mathrm d^2 \theta}{\mathrm dt^2}+g\sin\theta=0$

How can I solve the above equation?

$\endgroup$
  • 3
    $\begingroup$ The solution involves Jacobi's elliptic function "sn". There is a (very sketchy) description of how to integrate the equation at mai.liu.se/~halun/complex/elliptic. $\endgroup$ – Hans Lundmark Feb 1 '12 at 15:22
  • $\begingroup$ you can solve a small angle approximation by assuming $\sin\theta\approx\theta$ $\endgroup$ – yoyo Feb 1 '12 at 15:24
  • $\begingroup$ @yoyo Can't it be done without the small angle approximation? I wanted to derive a general equation for all cases. Anyways I don't know how to solve even with the approximation, so it'll be helpful even if you explain that method. $\endgroup$ – Green Noob Feb 1 '12 at 15:31
2
$\begingroup$

Use substitution : $\theta' =v$ ,therefore we have that :

$$\theta''=\frac{dv}{dt}\cdot \frac{dt}{d\theta}\cdot \frac{d\theta}{dt} \Rightarrow \theta''=\frac{dv}{d\theta}\cdot v \Rightarrow \theta''=v'\cdot v$$

where $v$ is function in terms of variable $\theta$ .So differential equation becomes :

$v' \cdot v +g \cdot \sin \theta=0$

which is separable differential equation .

$\endgroup$
  • $\begingroup$ Can you explain how it is in separable form? I'm a beginner so please bear with me.. $\endgroup$ – Green Noob Feb 1 '12 at 15:45
  • 2
    $\begingroup$ $\int v\,dv=-g\cdot \int \sin \theta \,d\theta$ $\endgroup$ – Peđa Terzić Feb 1 '12 at 15:48
  • 1
    $\begingroup$ then you have to integrate something like $$\frac{\theta'}{\sqrt{2g\cos\theta+C}}?$$ $\endgroup$ – yoyo Feb 1 '12 at 16:07
  • $\begingroup$ @yoyo I simplified a step further & got this equation : $\frac{d\theta}{\sqrt{2C+2gcos\theta}}=dt$ I'm not sure I can integrate this. $\endgroup$ – Green Noob Feb 1 '12 at 16:09
  • 1
    $\begingroup$ @GreenNoob: Writing $\cos\theta=1-2\sin^2(\theta/2)$ reveals that it's an elliptic integral, so it can't be done with elementary functions. When you integrate, $t$ will be an elliptic integral of $\theta$, and inverting that relation gives $\theta$ as an elliptic function of $t$. (Elliptic functions are by definition inverses of elliptic integrals, similar to how one can define the sine function as the inverse of the arcsin integral $\int_0^x \frac{ds}{\sqrt{1-s^2}}$.) $\endgroup$ – Hans Lundmark Feb 2 '12 at 7:19
2
$\begingroup$

Start with $$ \frac{1}{2}\frac{\mathrm{d}\dot{\theta}^{2}}{\mathrm{d}\theta} = \dot{\theta}\frac{\mathrm{d}\dot{\theta}}{\mathrm{d}\theta} = \frac{\mathrm{d}\theta}{\mathrm{d}t}\frac{\mathrm{d}\dot{\theta}}{\mathrm{d}\theta} = \frac{\mathrm{d}\dot{\theta}}{\mathrm{d}t} = \ddot{\theta} $$ Then your equation becomes $$ \frac{1}{2}\frac{\mathrm{d}\dot{\theta}^{2}}{\mathrm{d}\theta} = -\frac{g}{\ell}\sin(\theta) $$ or $$ \mathrm{d}\dot{\theta}^{2} = -\frac{2g}{\ell}\mathrm{d}\sin(\theta) \implies \dot{\theta}^{2} = \frac{2g}{\ell}\cos(\theta) + c_{1} $$ It's a bit easier if we assume initial conditions, say $\dot{\theta}(t_{0}) = \dot{\theta}_{0}$ and $\theta(t_{0}) = \theta_{0}$, so that $$ \dot{\theta}^{2} = \frac{2g}{\ell}\left[\cos(\theta)- \cos(\theta_{0}) + \frac{\ell\dot{\theta}_{0}^{2}}{2g} \right] $$ Then $$ \frac{\mathrm{d}\theta}{\mathrm{d}t} = \sqrt{ \frac{2g}{\ell}}\sqrt{\cos(\theta) - \cos(\theta_{0}) + \frac{\ell\dot{\theta}_{0}^{2}}{2g} } $$ so that $$ \mathrm{d}t = \sqrt{\frac{\ell}{2g}}\frac{\mathrm{d}\theta}{ \sqrt{\cos(\theta) - \cos(\theta_{0}) + \frac{\ell\dot{\theta}_{0}^{2}}{2g}}} $$ or $$ t_{f} - t_{0} = \sqrt{\frac{\ell}{2g}}\int_{\theta_{0}}^{\theta_{f}}\frac{\mathrm{d}\theta}{ \sqrt{\cos(\theta) - \cos(\theta_{0}) + \frac{\ell\dot{\theta}_{0}^{2}}{2g}}} $$ This equation is of the form $t = f(\theta)$. Your solution is given by $\theta = f^{-1}(t)$. That's about as much as you need to know, since it's more efficient to just solve the original equation numerically.

If you really need a closed form for $f$, Mathematica will give you one, in terms of the function EllipticF.

$\endgroup$
1
$\begingroup$

replacing $\sin\theta$ by $\theta$ (physically assuming small angle deflection) gives you a homogeneous second order linear differential equation with constant coefficients, whose general solution can be found in most introductory diff eq texts (or a google search). this new equation represents a simple harmonic oscillator (acceleration proportional to displacement, like a spring force). $$ \theta''+g\theta=0 $$ has solutions $A\cos(\sqrt{g}t)+B\sin(\sqrt{g}t)$. so, for example, if the initial displacement is $\theta_0$ and initial angular velocity is $0$ then the solution is $$ \theta_0\cos(\sqrt{g}t) $$

$\endgroup$
  • 1
    $\begingroup$ This doesn't exactly answer the question since this is just for an approximation to the original differential equation. $\endgroup$ – Garrett Jun 8 '14 at 23:49

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.