Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

2nd order non linear ODE

  1. Mar 12, 2010 #1
    So I have never solved a second order NON LINEAR ode before, but my current project has thrown one at me.

    I am going to use matlab to solve it when I have it in the correct state, but at the mo, I cant get it to that...
    Its three coulped equations, in terms of x,theta and x2. They are in a matrix, ma + cv + kx = force

    for each one.

    The way they are coupled though is a trigonometric relationship between x and theta and x2.

    theta = f(x)
    x2 = f(x)
    and x changes with time.

    What should I do to include this relationship?
  2. jcsd
  3. Mar 13, 2010 #2
    The question is too vague. Post the equations
  4. Mar 13, 2010 #3
    Please see attached PDF. It was easier to do it in word than using latex.

    Doing all those differentials correctly is very difficult, and it quicly becomes quite ugly.

    Attached Files:

  5. Mar 15, 2010 #4
    Sorry, but I think I don't understand where the constraint between [tex](\theta,x_{2},x)[/tex] is comming. Could you post the original problem? Anyway, instead of explicitly using the constraints in the equations of motion, you might try using lagrange multipliers, which work better in this cases. Add this to your lagrangian:


    And derive the equations of motion without substituing the constraints. When simulating, the equations for the evolution of the unknowns [tex]\Lambda_{1,2}[/tex] are [tex]\dot{\Lambda_{1,2}}=0[/tex], as they are constants
  6. Mar 15, 2010 #5
    The original problem is two pistons attached to one flywheel, where one of them has an applied force on it, and both are experiencing dampening, x is the displacement of the piston experiencing the external force, and x2 is the displacement on the other piston (attached image).

    then the equations I posted before have come from this diagram.

    I havent heard of Lagrange multipliers, i shall look them up

    Attached Files:

Share this great discussion with others via Reddit, Google+, Twitter, or Facebook