命令流如下，小弟做的是个碎石桩单桩的的半个模型，运行后就出现了math error:zero divide invalid 的除数为0的错误
检查了下参数似乎该赋得都赋了吧？麻烦各位老师帮我看看。。。谢谢了
n
gen zone radcyl p0 (0,0,0) p1 (10,0,0) p2 (0,0,-11) p3 (0,10,0) &
p4 (10,0,-11) p5 (0,10,-11) p6 (10,10,0) p7 (10,10,-11) &
p8 (.4,0,0) p9 (0,.4,0) p10 (.4,0,-11) p11 (0,.4,-11) &
size 3 11 6 15 ratio 1 1 1 1.15
gen zone radcyl p0 (0,0,-11) p1 (10,0,-11) p2 (0,0,-22) p3 (0,10,-11) &
p4 (10,0,-22) p5 (0,10,-22) p6 (10,10,-11) p7 (10,10,-22) &
p8 (.4,0,-11) p9 (0,.4,-11) p10 (.4,0,-22) p11 (0,.4,-22) &
size 3 11 6 15 ratio 1 1 1 1.15 fill;fill 代表填满洞室
gen zone reflect dd 270 dip 90
interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-10.9) radius .41 &
cylinder end1 (0,0,0) end2 (0,0,-10.9) radius .39 not
interface 1 face range cylinder end1 (0,0,-10.9) end2 (0,0,-11) radius .41
int 1 prop kn 1e8 ks 1e8 f 20 c 1e3
interface 1 perm = on
interface 1 effe = off
gen zone cyl p0 (0,0,11) p1 (.4,0,11) p2 (0,0,0) p3 (0,.4,11) &
p4 (.4,0,0) p5 (0,.4,0) &
size 3 11 6
group pile range z 0 11
gen zone reflect normal 1 0 0 o 0 0 0 ran gro pile
ini z add -11.0 range group pile
gro 淤泥 ran z -22 -12.9 &
gro pile not
gro 粉土2 ran z -12.9 -9.5 &
gro pile not
gro 粘土 ran z -9.5 -7.7 &
gro pile not
gro 粉土1 ran z -7.7 -2 &
gro pile not
gro 填土 ran z -2 0 &
gro pile not
fix x ran x -.1 .1
fix x ran x 9.9 10.1
fix x y z ran z -21.9 22.1
model mohr
prop bulk 0.95E+06 shear 0.83E+06 c 17.5e3 fri 15.8 ran gro 淤泥
prop bulk 5.00E+06 shear 2.31E+06 c 3.0e3 fri 29.5 ran gro 粉土2
prop bulk 0.95E+06 shear 0.87E+06 c 11.9e3 fri 15.5 ran gro 粘土
prop bulk 5.00E+06 shear 2.31E+06 c 2.6e3 fri 31 ran gro 粉土1
prop bulk 1.11e+06 shear 0.87e+06 c 10.8e3 fri 12.5 ran gro 填土
prop bu 1.67e10 sh 7.69e9 c 2e3 fri 30 ran gro pile
ini dens 1710 ran gro 淤泥
ini dens 1890 ran gro 粉土2
ini dens 1810 ran gro 粘土
ini dens 1910 ran gro 粉土1
ini dens 1870 ran gro 填土
ini dens 2400 ran gro pile
set grav 9.8
solve elas
save solopile.sav
