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!

Spring pendulum system fortran program

  1. Sep 30, 2011 #1
    I just want to know if someone has the fortran code for the numerical solution of the pendulum with a spring. And if it is so, can it write it here?
     
  2. jcsd
  3. Oct 1, 2011 #2
    Have you tried to do this on your own already? https://www.physicsforums.com/showthread.php?t=94379" Show what you've done so far so we can help you.
     
    Last edited by a moderator: Apr 26, 2017
  4. Oct 1, 2011 #3
    Ok, this is what I have:

    PROGRAM DIFFERENTIAL_EQUATIONS_SYSTEM

    WRITE (*,*) 'WRITE A, B, N, ALF1 AND ALF2'
    READ (*,*) A,B,N,ALF1,ALF2
    H=(B-A)/N
    T=A
    W1=ALF1
    W2=ALF2

    OPEN(10,FILE='RUNGE_KUTTA_HIGHER_ORDER.TXT')

    DO I=1,N

    WRITE (10,*) T,W2
    WRITE (*,*) T,W2

    W1=W1+H*F1(T,W1,W2)
    W2=W2+H*F2(T,W1,W2)
    T=A+I*H

    END DO
    END

    FUNCTION F1(T,W1,W2)
    REAL W1,W2,T
    F1=W1
    RETURN
    END

    FUNCTION F2(T,W1,W2)
    REAL W1,W2,T
    F2=-3.0*W2-(9.8/3.0)*SIN(W1)
    RETURN
    END

    maybe my way of programming is not the best. What I have are the equations of motion for the pendulum-spring system, they are coupled, and I don't know how to solve them numerically, any help???

    And it would be any advice to improve my programming it is welcomed.
     
  5. Oct 1, 2011 #4
    Could you please explain a bit more about the problem and your code? When you say pendulum-spring system, do you mean a rigid pendulum with a mass on the end is hanging from a spring? Can you show us the equations of motion you're using? And what method are you trying to use, 2nd-order Runge-Kutta?

    As for the variables, I'm guessing A and B are initial and final times, N is number of timesteps, ALF1 and ALF2 are initial angles of the spring and pendulum, respectively. Correct? More descriptive variable names in general are very helpful if anyone besides yourself needs to understand your code.

    A Fortran tip you may not be aware of: declare your variables explicity everywhere (like how you did in your functions: REAL W1,W2,T), and put IMPLICIT NONE as the first line in all of your programs and functions/subroutines, which forces you to declare your varibles. If you do not explicity declare varibles, Fortran implicity does so in the following manner. If the first letter of your variable name is I, J, K, L, M or N, it will be an INTEGER, otherwise it will be a REAL. This is widely considered to be one of the worst features of Fortran, and using IMPLICIT NONE is always recommended. In your code there are no problems with this, but it's easy to make mistakes.
     
  6. Oct 1, 2011 #5
    The theory is on the pages 42-43 of the book.pdf file. I made a mistake with the name of the code, the method it isn't Runge-Kutta it is Euler's Method but for a second order differential equation.

    All you assume on second paragraph it is right, about the variables. And thanks for the tip, I will keep it in mind.
     

    Attached Files:

    • book.pdf
      book.pdf
      File size:
      205.1 KB
      Views:
      166
    Last edited: Oct 1, 2011
  7. Oct 2, 2011 #6
    You'll have to introduce some phase/state variables to convert those two 2nd-order ODEs into four 1st-order ODEs, and then solve those. http://12000.org/my_courses/UC_davis/spring_2011/MAE_121_eng_dynamics/lab/lab_one/report/report.htm" [Broken] (eqn 3 and the ones above it). Hope that helps.
     
    Last edited by a moderator: May 5, 2017
  8. Dec 22, 2011 #7
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Spring pendulum system fortran program
  1. Fortran Programming (Replies: 5)

Loading...