[问题] 空气阻力code问题

楼主: hank546288 (hanker8862)   2017-10-29 18:05:24
这是我的code
我该怎么让f会随着ball.v变化呢?
t=0
dt=0.001
m=5
g=-9.8
ball.a=g
p=1.2
c=0.5
A=0.25*pi
f=0.5*p*c*A*ball.v**2
while ball.pos.y>=0.5:
t=t+dt
rate(1000)
f+=0.5*p*c*A*(ball.a*dt)**2
ball.v+=ball.a*dt
ball.pos.y+=ball.v*dt
ball.a=ball.a+f/m
如果照我这样写 f会变成v^2+adt^2 但我要的是f=(v+adt)^2 不知道该怎么写出来
作者: neil987 (R5大小姐-EX人品崩坏)   2017-10-29 20:17:00
把空气阻力简化看 f=kv^2 所以你下面 f+=kadt^2就会变成f = kv^2 +k adt^2v2=(v1+at) f2 = k(v2)^2 所以你只要改变v 再带入f=kv^2你while循环里面 不是有v=v+adt吗= = 就是那个ball.v+=ball.a*dt f = k*ball.v**2你的v有先改吗

Links booklink

Contact Us: admin [ a t ] ucptt.com