参照《FLAC3D基础工程与实例》130页建的桩无法收敛。。。求高手帮忙呀

jasonhulk

将命令流赋在下面了，求解初始平衡的时候一直无法收敛，求解原因哈。。谢谢谢

n
gen zone radcyl p0  (0,0,0)  p1  (3,0,0)  p2  (0,0,-11)  p3  (0,3,0)  &
                p4  (3,0,-11) p5  (0,3,-11) p6  (3,3,0)   p7  (3,3,-11) &
                p8  (.4,0,0)  p9  (0,.4,0)  p10 (.4,0,-11)  p11 (0,.4,-11) &
                size 3 10 4 6  ratio 1 1 1 1.15
gen zone radcyl p0  (0,0,-11)  p1  (3,0,-11)  p2  (0,0,-22)  p3  (0,3,-11)  &
                p4  (3,0,-22)  p5  (0,3,-22)  p6  (3,3,-11)  p7  (3,3,-22)  &
                p8  (.4,0,-11) p9  (0,.4,-11) p10 (.4,0,-22) p11 (0,.4,-22) &
                size 3 5 4 6  ratio 1 1 1 1.15 fill;fill 代表填满洞室
gen zone reflect dd 270 dip 90
gen zone reflect nor 0 1 0
gro 淤泥 ran z -22 -11.9
gro 粉土2 ran z -11.9 -9.5
gro 粘土 ran z -9.5 -7.7
gro 粉土1 ran z -7.7 -2
gro 填土 ran z -2 0
gen zone cyl p0 (0,0,12) p1 (.4,0,12) p2 (0,0,1)  p3 (0,.4,12)  &
             p4 (.4,0,1) p5 (0,.4,1) &
             size 3 10 4
gen zone reflect  normal 0 1 0 ran z 1 12
gen zone reflect  normal 1 0 0 ran z 1 12
group 碎石桩 range z 1 12
interface 1 face range cylinder end1 (0,0,0)  end2 (0,0,-11.1) radius .40 &
                       cylinder end1 (0,0,0)  end2 (0,0,-11.1) radius .39 not
interface 2 face range cylinder end1 (0,0,-10.9)  end2 (0,0,-11.1) radius .40

ini z add -12 range group 碎石桩
save 网格.sav

rest 网格.sav
config dyn fluid
set dyn off fluid off
model mohr ran gro 淤泥 a gro 粉土1 a gro 粉土2 a gro 粘土 a gro 填土 a
prop bulk 0.95E+06 shear 0.83E+06 ten 1e10 c 17.5e10 fri 15.8 ran gro 淤泥
prop bulk 5.00E+06 shear 2.31E+06 ten 1e10 c 3.0e10  fri 29.5 ran gro 粉土2
prop bulk 0.95E+06 shear 0.87E+06 ten 1e10 c 11.9e10 fri 15.5 ran gro 粘土
prop bulk 5.00E+06 shear 2.31E+06 ten 1e10 c 2.6e10 fri 31 ran gro 粉土1
prop bulk 1.11e+06 shear 0.87e+06 ten 1e10 c 10.8e10 fri 12.5 ran gro 填土
model elas ran gro 碎石桩
prop bu 1.67e10 sh 7.69e9 ran gro 碎石桩
interface 1 prop kn 1e8  ks 1e8  fric 20 coh 30000
interface 2 prop kn 1e8  ks 1e8  fric 20 coh 30000
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 碎石桩
fix x ran x 3
fix x ran x -3
fix y ran y 3
fix y ran y -3
fix x y z ran z -22
set grav 9.8
solve
model mohr ran gro 淤泥 a gro 粉土1 a gro 粉土2 a gro 粘土 a gro 填土 a
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 填土
mo e ran gro 碎石桩
prop bu 1.67e10 sh 7.69e9 ran gro 碎石桩
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 碎石桩
solve
sav 初始平衡.sav

jasonhulk

是网格的问题还是土层参数的问题啊？