- 积分
- 0
- 注册时间
- 2010-2-2
- 仿真币
-
- 最后登录
- 1970-1-1
|
我编了一段 就是墙不能动,请高手帮帮忙:
; FNAME: fang.DAT
new
set random
gen id 1 2000 rad 0.02 0.02 x 0 1 y -1 1 z 0 1 tries 40000
wall id 1 face (0,-1,0) (1,-1,0) (1 1 0) (0 1 0)
wall id 2 face (0 -1 0) (0 -1 0.5) (1 -1 0.5) (1 -1 0)
wall id 3 face (0 -1 0) (0 1 0) (0 1 0.5) (0 -1 0.5)
wall id 4 face (0 1 0.5) (0 1 0) (1 1 0) (1 1 0.5)
wall id 5 face (1 1 0) (1 -1 0) (1 -1 0.5) (1 1 0.5)
wall id 6 face (1 -1 0.5) (0 -1 0.5) (0 1 0.5) (1 1 0.5)
wall id 1 kn 1e6 ks 1e6 fric 1.0
wall id 2 kn 1e6 ks 1e6 fric 1.0
wall id 3 kn 1e6 ks 1e6 fric 1.0
wall id 4 kn 1e6 ks 1e6 fric 1.0
wall id 5 kn 1e6 ks 1e6 fric 1.0
wall id 6 kn 1e6 ks 1e6 fric 1.0
prop kn 1e4 ks 1e4 dens 2000 fric 0.3
set grav 0 0 -9.81
def startup
old_time = time
omiga = 100.48
end
def shak_it
while_stepping
real_time = time
if real_time < length
y_vel = ampl*sin(omiga*real_time)
z_vel = ampl*sin(omiga*real_time)
end_if
end
prop yvel y_vel zvel z_vel
set time = 0 ampl=4 length=10
startup
pl set mag 1
pl add ball blue
pl add wall white
pl add axes
plot set rotation (0,0,-120)
plot show
cycle 10
set dt 0.005
cycle 1000 |
|