new
set random 10001
set orientation on
def basic_data
height = 0.1;io.in('请输入试样高度') ;0.1
width = 0.05;io.in('请输入试样宽度') ;0.05
cylinder_axis_vec = vector(0,0,1)
cylinder_base_vec = vector(0.0,0.0,[-0.7*height])
cylinder_top_vec = vector(0.0,0.0,[0.7*height])
cylinder_height = 2 * 0.7 * height
cylinder_rad = 0.5 * width
bottom_disk_position_vec = vector(0.0,0.0,[-0.5*height])
top_disk_position_vec = vector(0.0,0.0,[0.5*height])
w_resolution = 0.5
poros =0.16;io.in('请输入试样孔隙率');ball:0.35,clump:0.35
rd_min=1.5e-3;io.in('请输入颗粒最小半径')
rd_max=3.0e-3;io.in('请输入颗粒最大半径')
end
@basic_data
domain extent ([-width],[width]) ([-width],[width]) ([-height],[height])
;cmat default model linearcbond
cmat default model linear
wall generate id 1 cylinder radius [0.5*width] height [1.4*height] base 0 0 [-0.7*height] axis 0 0 1 ...
cap false false onewall ...
resolution @w_resolution
wall generate id 2 plane position 0 0 [-0.5*height] dip 0 ddir 0
wall generate id 3 plane position 0 0 [0.5*height] dip 0 ddir 0
geometry set geo_cylinder
geometry generate cylinder axis (0,0,1) &
base (0.0,0.0,[-0.5*height]) height [2*0.5*height] &
cap true true &
radius [cylinder_rad] &
resolution 0.1
;-----------------------------------------------------------------------------------------------------
;-----------------------------------------------------------------------------------------------------
ball distribute porosity [poros] radius [rd_min] range geometry geo_cylinder count odd direction (0.0 0.0 1.0)
ball attribute xposition multiply 0.85 yposition multiply 0.85 zposition multiply 0.95
ball attribute density 2.66e3 damp 0.7
ball property kn 1.0e5
wall property kn 1e7 ks 1e7

可以看一下手册ball delete range 命令