/Python/Http://techforcurious.website/simulation-of-pendulum-vpython-tutorial-visual-python/

Http://techforcurious.website/simulation-of-pendulum-vpython-tutorial-visual-python/

from visual import*
display(width=600,height=600,center=vector(0,12,0),background=color.white)
g=9.8 # acceleration due to gravity
bob=sphere(pos=vector(5,2,0),radius=0.5,color=color.blue)
pivot=vector(0,20,0)
roof=box(pos=pivot,size=vector(10,0.5,10),color=color.green)
rod=cylinder(pos=pivot,axis=bob.pos-pivot,radius=0.1,color=color.red)
t=0 # time 
dt=0.01 # time interval 
l=mag(bob.pos-pivot) # length of pendulum
cs=(pivot.y-bob.pos.y)/l # calculation of cos(theta) 
theta=acos(cs) # angle with vertical direction
vel=0.0 # angular velocity
while (t<100):
  rate(100) # maximum 100 calculations per second
  acc=-g/l*sin(theta) # updating of angular acceleration
  theta=theta+vel*dt # updating of angular position
  vel=vel+acc*dt # updating of angular velocity
  bob.pos=vector(l*sin(theta),pivot.y-l*cos(theta),0) # cal. position
  rod.axis=bob.pos-rod.pos # updating other end of rod of pendulum
  t=t+dt # updating time

Similar codes

35 votes294 views6 answers
Create pandas dataframe from dictionary
Python
1 vote285 views1 answer
Best live score api
Python
1 vote308 views1 answer
For idx, col_name in enumerate(X_train.columns): print("The coefficient for {} is {}".format(file_name, regression_model.coef_[0][idx]))
Python
0 vote339 views1 answer
'datetime.date' object has no attribute 'split'
Python
0 vote277 views1 answer
How to record pyttsx3 file using python
Python
90 votes309 views6 answers
Python get current directory
Python
0 vote296 views1 answer
Take substring of every element in dataframe
Python
34 votes302 views2 answers
Python shuffle list with seed
Python
0 vote282 views1 answer
Python - Comment lire une ligne de fichier par ligne
Python
1 vote281 views2 answers
Setup VS Code for python
Python