Hello and welcome to our community! Is this your first visit?
Register
Enjoy an ad free experience by logging in. Not a member yet? Register.
Results 1 to 3 of 3
  1. #1
    New to the CF scene
    Join Date
    May 2009
    Posts
    2
    Thanks
    0
    Thanked 0 Times in 0 Posts

    Python/C/C++ or similar program question

    I am programming in python and know some C++ as well so I can interpret something similar to those two. I have written a programm that models a mass on a spring in harmonic motion with a damping force. It displays a graph in phase space and on a velocity and position vs. time graph. I would like to add in an extra force that will only activate if the amplitude of the velocity and position curve drops below some minimum. Essentially like giving the spring a little push when it's velocity drops too low. Any help on how to do so would be greatly appreciated. The program is below in python programming language which is similar to C++, comments follow #



    from scipy import *
    from pylab import *

    k=0.05 #Define constants (k=spring constant,b=damping force,m=mass on spring)
    b=0.1
    m=5
    npts= 100000
    v=zeros(npts) #Define arrays for velocity, position, and time
    x=zeros(npts)
    t=zeros(npts)
    v[0]=30.00 #Initialize arrays
    x[0]=5.00
    t[0]=0.00
    dt=0.01
    for i in arange(npts-1):
    v[i+1]=v[i]+dt*(-k*x[i]-b*v[i]+sin(0.01*t[i]))/m #Difference equation for velocity
    x[i+1]=x[i]+ v[i]*dt #Difference equation for position
    t[i+1]=t[i]+dt #time euation
    figure(1) #the phase plot
    plot(v,x)
    grid(True) #shows grid on the graph
    ylabel('velocity')
    xlabel('position')
    figure(2) #time plots of v and x
    plot(t,v,'-',label='velocity')
    plot(t,x,'--',label='position')
    legend() #shows the labels on graph
    xlabel('Time')
    show()

  • #2
    Supreme Overlord Spookster's Avatar
    Join Date
    May 2002
    Location
    Marion, IA USA
    Posts
    6,280
    Thanks
    4
    Thanked 83 Times in 82 Posts
    Being that Python is whitespace sensitive you should really put code or PHP tags around your code so it retains the code formatting.
    Spookster
    CodingForums Supreme Overlord
    All Hail Spookster

  • #3
    New to the CF scene
    Join Date
    May 2009
    Posts
    2
    Thanks
    0
    Thanked 0 Times in 0 Posts
    But can you help me with programming in a force that only activates when the amplitude of the function drops below a critical value?


  •  

    Tags for this Thread

    Posting Permissions

    • You may not post new threads
    • You may not post replies
    • You may not post attachments
    • You may not edit your posts
    •