下边是命令流，也不知道是interface出问题了还是出现零厚度了，各位帮忙看看啊，谢谢了
gen zone radcy size 3 15 5 20 p0 0 0 51 p1 40 0 51 p2 0 0 1 p3 0 40 51 p4 40 0 1 p5 0 40 1 p6 40 40 51 p7 40 40 1 p8 5 0 51 p9 0 5 51 p10 5 0 1 p11 0 5 1
gen zone reflect norm 1 0 0
gen zone reflect norm 0 1 0
interface 1 face range cylinder end1 0 0 0.9 end2 0 0 51.1 radius 5.1 cylinder end1 0 0 0.9 end2 0 0 51.1 radius 4.9 not
gen zone cyl p0 0 0 102 p1 5 0 102 p2 0 0 51 p3 0 5 102 p4 5 0 51 p5 0 5 51 size 1 20 3 group pile
gen zone reflect norm 1 0 0 range group pile
gen zone reflect norm 0 1 0 range group pile
ini z add -50 range group pile
group soil range z 20 51 group pile not
group rock1 range z 1 20 group pile not
gen zone brick size 10 10 4 p0 -40 -40 -6 p1 40 -40 -6 p2 -40 40 -6 p3 -40 -40 0 group rock2
ini z add 1.0 range group rock2
interface 2 face range cylinder end1 0 0 0.9 end2 0 0 1.1 radius 5.1
plot add grid yellow range group soil
plot add grid red range group rock1
plot add grid green range group rock2
plot add grid blue range group pile
plot add interface blue
plot show
model mohr range group pile not
pro bu 0.53e9 sh 0.33e9 coh 4.387e6 fric 38 tension 1.63e6 dens 2200 range group rock1
pro bu 0.54e9 sh 0.35e9 coh 4.154e6 fric 40 tension 2.43e6 dens 2500 range group rock2
pro bu 0.43e9 sh 0.14e9 coh 0.55e6 fric 20 tension 1.15e6 dens 1700 range group soil
model elas range group pile
pro bu 0.77e9 sh 0.44e9 dens 3000 range group pile
fix x range x -39.9 -40.1 39.9 40.1
fix y range y -39.9 -40.1 39.9 40.1
fix x y z range z -5.1 -4.9
set gravity 0 0 -9.81
inter 1 prop kn 300e9 ks 300e9 tens 1e10 sbratio 1 di 30 co 5e6 fric 20
inter 2 prop kn 300e9 ks 300e9 tens 1e10 sbratio 1 di 50 co 5e6 fric 40
hist id 1 unba
hist id 2 gp zdisp 0 0 51
plot add hist 1
step 15000
save 7.sav
ini szz -3.5e6 range z 51.9 52.1
in xdis 0 ydis 0 zdis 0
in xvel 0 yvel 0 zvel 0
set large
step 15000
plot set plane ori 0 0 0 norm 0 1 0
plot contour szz plane
