- 积分
- 0
- 注册时间
- 2011-9-30
- 仿真币
-
- 最后登录
- 1970-1-1
|
new
set random
set disk on
wall id 1 ks 1e19 kn 1e19 fric 0.2 nodes -0.2,-0.1 50.2,-0.1
wall id 2 ks 1e19 kn 1e19 fric 0.2 nodes 50.2,-0.1 50.2,50
wall id 3 ks 1e19 kn 1e19 fric 0.2 nodes 50.2,50 -0.2,50
wall id 4 ks 1e19 kn 1e19 fric 0.2 nodes -0.2,50 -0.2,-0.1
def meiceng
xc = x0
yc = y0
rc = radius
idc = id_start
r2 = 2.0 * radius
yinc = radius * sqrt(3.0)
loop row (1,n_row)
loop col (1,n_col)
command
ball id=idc x=xc y=yc rad=rc
end_command
idc = idc + 1
xc = xc + r2
end_loop
yc = yc + yinc
xc = x0 + radius * (row - (row/2) * 2)
end_loop
end
set echo off
set x0=0.25 y0=0.25 radius=0.25
set id_start=1 n_col=100 n_row=100
meiceng
set echo on
plot set cap size 20
plot add wall brown
plot add axes red
plot add ball black
plot show
用group分组分不了,那位大神给分析下 非常感谢!
|
|