VPython - example - 模拟球在封闭箱子内的碰撞

本作品采用知识共享署名-非商业性使用-相同方式共享 3.0 Unported许可协议进行许可。允许非商业转载，但应注明作者及出处。

作者：liuyuan_jq

2011-04-10

From http://maths.anu.edu.au/comptlsci/Tutorial-Gas/bounce2.html

---

 

 

########################################## Import the library(s)#########################################from visual import *########################################### Create Wall(s)##########################################side=4.0thk=0.3s2 = 2*side - thks3 = 2*side + thkwallR = box (pos=( side, 0, 0), length=thk, height=s2, width=s3, color = color.red)wallL = box (pos=(-side, 0, 0), length=thk, height=s2, width=s3, color = color.red)wallB = box (pos=(0, -side, 0), length=s3, height=thk, width=s3, color = color.blue)wallT = box (pos=(0, side, 0), length=s3, height=thk, width=s3, color = color.blue)wallBK = box(pos=(0, 0, -side), length=s2, height=s2, width=thk, color = (0.7,0.7,0.7))########################################### Create Ball(s)##########################################ball_radius=1.0maxpos=side-thk/2-ball_radiusmaxv=2.0ball = sphere (color = color.green, radius = ball_radius)ball.velocity = vector (-1.5, -2.3, +2.7)########################################### Time loop for moving Ball(s)###########################################timestep = 0.05while (1==1): # Set number of times loop is repeated per second rate(100) # Move ball(s) ball.pos = ball.pos + ball.velocity*timestep #check for collisions with the walls #right wall if ball.pos.x > maxpos: ball.velocity.x = -ball.velocity.x #reflect velocity ball.pos.x=2*maxpos-ball.pos.x #reflect position #left wall if ball.pos.x < -maxpos: ball.velocity.x = -ball.velocity.x ball.pos.x=-2*maxpos-ball.pos.x # roof if ball.pos.y > maxpos: ball.velocity.y = -ball.velocity.y ball.pos.y=2*maxpos-ball.pos.y #floor if ball.pos.y < -maxpos: ball.velocity.y = -ball.velocity.y ball.pos.y=-2*maxpos-ball.pos.y #back wall if ball.pos.z > maxpos: ball.velocity.z = -ball.velocity.z ball.pos.z=2*maxpos-ball.pos.z #front wall if ball.pos.z < -maxpos: ball.velocity.z = -ball.velocity.z ball.pos.z=-2*maxpos-ball.pos.z