1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Break Second order ODE into a system of first order ODE's

  1. Sep 1, 2011 #1
    1. The problem statement, all variables and given/known data

    I havent done this for several years and have forgotten. Kicking myself now over it since it looks like something so simple but i cannot figure it out.... I need to break this second order ODE into a system of first order ODE's in matrix form to use within a crank nicolson method.

    [itex] \frac{d\Theta^{2}}{dt^{2}} + c\frac{d\Theta}{dt} + \frac{g}{L}sin(\Theta) = 0[/itex]



    3. The attempt at a solution

    let
    [itex] \phi_{1} = \Theta[/itex]
    [itex] \frac{\phi_{1}}{dt} = \phi_{2} [/itex]
    [itex] \frac{\phi_{2}}{dt} = -c\phi_{2} - \frac{g}{L}sin{\phi_{1}} [/itex]


    now problem being the [itex] \sin{\phi} [/itex], how do i take the phi out! K is meant to be the coefficients of the terms infront of phi, but in this case its within the sin :S



    [itex] \left[ {\begin{array}{cc}
    \frac{\phi_{1}}{dt} \\
    \frac{\phi_{2}}{dt} \\
    \end{array} } \right]
    =
    \left[ {\begin{array}{cc}
    0 & 1 \\
    unknown & -c \\
    \end{array} } \right]


    \left[ {\begin{array}{cc}
    \phi_{1} \\
    \phi_{2} \\
    \end{array} } \right]
    [/itex]

    Cheers Trent
     
    Last edited: Sep 1, 2011
  2. jcsd
  3. Sep 1, 2011 #2

    lanedance

    User Avatar
    Homework Helper

    This is a non-linear DE, hence the difficulties

    If theta was very small you could use the small angle approximation to linearise the equation
    [tex] sin(\theta(t))\approx \theta(t) [/tex]
     
  4. Sep 1, 2011 #3

    lanedance

    User Avatar
    Homework Helper

    now looking at crank-nicholson which is finite difference method, seems to be set up for partial DEs
    http://en.wikipedia.org/wiki/Crank–Nicolson_method

    as this is a 2nd order nonlinear ordinary DE, why not something like runge kutta?
     
  5. Sep 1, 2011 #4
    unfortunately its not, we are modelling a pendulum which is lubricated well** :(
     
  6. Sep 2, 2011 #5
    we have been instructed to use Crank-Nicholson for some stupid reason in our design brief.

    Im looking up other methods iterative techniques such as newtons method which may be applied within the method....

    Any suggestions or thoughts....?????
     
  7. Sep 2, 2011 #6

    lanedance

    User Avatar
    Homework Helper

    i haven't used it but eve4rything i see on crank-nicholson is for 2d (x,t) differntial equations, so not really sure how it applies here
     
  8. Sep 2, 2011 #7
    All good

    Using the Crank-Nicholson approach, and than applying newtons law, which involves taking the jacobian etc and solving for the residue to be zero works.

    Thanks for your time and help lanedance
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Break Second order ODE into a system of first order ODE's
Loading...