new
SET random ; 随机生成
plo wall
plo add ball
; ----------------------------------------------------
def make_walls ;生成墙
wextend =0.1
hextend =0.1
w_stiff= 1.6e10 ;墙的刚度，取球的1.1倍
_width=-width
_x0 = _width*(1.0 + wextend) ;width是长方体的宽度的一半
_y0 = _width*(1.0 + wextend)
_z0 = 0
_x1 = width*(1.0 + wextend)
_y1 = _width*(1.0 + wextend)
_z1 = 0
_x2 = width*(1.0 + wextend)
_y2 = width*(1.0 +wextend)
_z2 = 0
_x3 = _width*(1.0 + wextend)
_y3 = width*(1.0 + wextend)
_z3 = 0
command
wall id=1 kn=w_stiff face (_x0,_y0,_z0) (_x1,_y1,_z1) (_x2,_y2,_z2) &
(_x3,_y3,_z3) ;要按右手法则建立,而且pfc里面的坐标轴是一定的，要按照里面的来
end_command
_x0 = _width*(1.0 + wextend)
_y0 = _width*(1.0 + wextend)
_z0 = height ;height是长方体的高度
_x1 = _width*(1.0 + wextend)
_y1 = width*(1.0 + wextend)
_z1 = height
_x2 = width*(1.0 + wextend)
_y2 = width*(1.0 + wextend)
_z2 = height
_x3 = width*(1.0 + wextend)
_y3 = _width*(1.0 + wextend)
_z3 = height
command
wall id=2 kn=w_stiff face (_x0,_y0,_z0) (_x1,_y1,_z1) (_x2,_y2,_z2) &
(_x3,_y3,_z3)
end_command
_x0 = width
_y0 = _width*(1.0 + wextend)
_z0 = height*(1.0 + hextend)
_x1 = width
_y1 = width*(1.0 + wextend)
_z1 = height*(1.0 + hextend)
_x2 = width
_y2 = width*(1.0 + wextend)
_z2 = -hextend*height
_x3 = width
_y3 = _width*(1.0 + wextend)
_z3 = -hextend*height
command
wall id=3 kn=w_stiff face (_x0,_y0,_z0) (_x1,_y1,_z1) (_x2,_y2,_z2) &
