我编写的一个粉喷桩处理地基，初始始终不平衡，于是我把模型改小了，但是还是不平衡，大家帮我看看怎么回事，谢谢了啊。
;桩体周围网格
;1号桩体
gen zone radcyl p0 (0,0,0) p1 (1,0,0) p2 (0,0,-7) p3 (0,1,0) p4 (1,0,-7) &
p5 (0,1,-7) p6 (1,1,0) p7 (1,1,-7) p8 (0.25,0,0) p9 (0,0.25,0) &
p10 (0.25,0,-7) p11 (0,0.25,-7) size 1 14 4 4 ratio 1 1 1 1 ;7m 粉土
gen zone radcyl p0 (0,0,-7) p1 (1,0,-7) p2 (0,0,-9) p3 (0,1,-7) &
p4 (1,0,-9) p5 (0,1,-9) p6 (1,1,-7) p7 (1,1,-9) p8 (0.25,0,-7) &
p9 (0,0.25,-7) p10 (0.25,0,-9) p11 (0,0.25,-9) size 1 4 4 4 &
ratio 1 1 1 1 ;2m 粉质粘土
gen zone radcyl p0 (0,0,-9) p1 (1,0,-9) p2 (0,0,-15) p3 (0,1,-9) &
p4 (1,0,-15) p5 (0,1,-15) p6 (1,1,-9) p7 (1,1,-15) p8 (0.25,0,-9) &
p9 (0,0.25,-9) p10 (0.25,0,-15) p11 (0,0.25,-15) size 1 12 4 4 &
ratio 1 1 1 1 ;6m 粉砂
;桩体以下网格
gen zone radcyl p0 (0,0,-15) p1 (1,0,-15) p2 (0,0,-21) p3 (0,1,-15) &
p4 (1,0,-21) p5 (0,1,-21) p6 (1,1,-15) p7 (1,1,-21) p8 (0.25,0,-15) &
p9 (0,0.25,-15) p10 (0.25,0,-21) p11 (0,0.25,-21) size 1 12 4 4 &
ratio 1 1 1 1 fill ;6m 粉砂
;2号桩体
gen zone copy 2 0 0
gen zone reflect normal 1 0 0 ori 2 0 0 range x 2 3
;3号桩体
gen zone copy 0 2 0 range x 0 1
gen zone reflect normal 0 1 0 ori 0 2 0 range y 2 3
;4号桩体
gen zone copy 0 2 0 range x 1 3
gen zone reflect normal 0 1 0 ori 0 2 0 range x 1 3 y 2 3
;**********************建立外土体模型********************
;外土体模型的建立
gen zone radtunnel p0 (0,0,0) p1 (5,0,0) p2 (0,0,-21) p3 (0,5,0) &
p4 (5,0,-21) p5 (0,5,-21) p6 (5,5,0) p7 (5,5,-21) &
p8 (3,0,0) p9 (0,3,0) p10 (3,0,-21) p11 (0,3,-21) &
p12 (3,3,0) p13 (3,3,-21) size 6 42 6 10 ratio 1 1 1 1
group ft1 range z -7 0
group fznt2 range z -9 -7
group fs3 range z -15 -9
group fs4 range z -21 -15
;**********************建立接触面模型********************
;1号接触面
interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-7) radius 0.26 &
cylinder end1 (0,0,0) end2 (0,0,-7) radius 0.24 not
interface 2 face range cylinder end1 (0,0,-7) end2 (0,0,-9) radius 0.26 &
cylinder end1 (0,0,-7) end2 (0,0,-9) radius 0.24 not
interface 3 face range cylinder end1 (0,0,-9) end2 (0,0,-15) radius 0.26 &
cylinder end1 (0,0,-9) end2 (0,0,-15) radius 0.24 not
interface 4 face range cylinder end1 (0,0,-14.9) end2 (0,0,-15.1) radius 0.25
;2号接触面
interface 5 face range cylinder end1 (2,0,0) end2 (2,0,-7) radius 0.26 &
cylinder end1 (2,0,0) end2 (2,0,-7) radius 0.24 not
interface 6 face range cylinder end1 (2,0,-7) end2 (2,0,-9) radius 0.26 &
cylinder end1 (2,0,-7) end2 (2,0,-9) radius 0.24 not
interface 7 face range cylinder end1 (2,0,-9) end2 (2,0,-15) radius 0.26 &
cylinder end1 (2,0,-9) end2 (2,0,-15) radius 0.24 not
interface 8 face range cylinder end1 (2,0,-14.9) end2 (2,0,-15.1) radius 0.25
;3号接触面
interface 9 face range cylinder end1 (0,2,0) end2 (0,2,-7) radius 0.26 &
cylinder end1 (0,2,0) end2 (0,2,-7) radius 0.24 not
interface 10 face range cylinder end1 (0,2,-7) end2 (0,2,-9) radius 0.26 &
cylinder end1 (0,2,-7) end2 (0,2,-9) radius 0.24 not
interface 11 face range cylinder end1 (0,2,-9) end2 (0,2,-15) radius 0.26 &
cylinder end1 (0,2,-9) end2 (0,2,-15) radius 0.24 not
interface 12 face range cylinder end1 (0,2,-14.9) end2 (0,2,-15.1) radius 0.25
;4号接触面
interface 13 face range cylinder end1 (2,2,0) end2 (2,2,-7) radius 0.26 &
cylinder end1 (2,2,0) end2 (2,2,-7) radius 0.24 not
interface 14 face range cylinder end1 (2,2,-7) end2 (2,2,-9) radius 0.26 &
cylinder end1 (2,2,-7) end2 (2,2,-9) radius 0.24 not
interface 15 face range cylinder end1 (2,2,-9) end2 (2,2,-15) radius 0.26 &
cylinder end1 (2,2,-9) end2 (2,2,-15) radius 0.24 not
interface 16 face range cylinder end1 (2,2,-14.9) end2 (2,2,-15.1) radius 0.25
;**********************建立桩体模型********************
;1号桩
gen zone cyl p0 (0,0,21) p1 (0.25,0,21) p2 (0,0,6) p3 (0,0.25,21) &
p4 (0.25,0,6) p5 (0,0.25,6) size 1 30 4
;2号桩体
gen zone copy 2 0 0 range z 1 21
gen zone reflect normal 1 0 0 ori 2 0 0 range x 2 3 z 1 21
;3号桩体
gen zone copy 0 2 0 range x 0 1 z 1 21
gen zone reflect normal 0 1 0 ori 0 2 0 range y 2 3 z 1 21
;4号桩体
gen zone copy 0 2 0 range x 1 3 z 1 21
gen zone reflect normal 0 1 0 ori 0 2 0 range x 1 3 y 2 3 z 1 21
group zhuang range z 6 21
ini z add -21.0 range group zhuang
;**************************设置材料属性****************************
model elas range group zhuang
prop bulk 2.9033e6 shear 1.34e6 range group zhuang
model mo range group ft1
prop bulk 6.073e6 shear 2.278e6 coh 1e12 fric 25 ten 1e12 range group ft1
model mo range group fznt2
prop bulk 3.071e6 shear 1.151e6 coh 1e12 fric 21 ten 1e12 range group fznt2
model mo range group fs3
prop bulk 10e6 shear 6.429e6 coh 1e12 fric 31 ten 1e12 range group fs3
model mo range group fs4
prop bulk 10e6 shear 6.429e6 coh 1e12 fric 31 ten 1e12 range group fs4
;**************************分界面性质****************************
interface 1 prop kn 1.6e12 ks 1.6e12 fric 8.3 coh 29000
interface 2 prop kn 1.6e12 ks 1.6e12 fric 5.9 coh 29000
interface 3 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
interface 4 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
interface 5 prop kn 1.6e12 ks 1.6e12 fric 8.3 coh 29000
interface 6 prop kn 1.6e12 ks 1.6e12 fric 5.9 coh 29000
interface 7 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
interface 8 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
interface 9 prop kn 1.6e12 ks 1.6e12 fric 8.3 coh 29000
interface 10 prop kn 1.6e12 ks 1.6e12 fric 5.9 coh 29000
interface 11 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
interface 12 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
interface 13 prop kn 1.6e12 ks 1.6e12 fric 8.3 coh 29000
interface 14 prop kn 1.6e12 ks 1.6e12 fric 5.9 coh 29000
interface 15 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
interface 16 prop kn 1.6e12 ks 1.6e12 fric 8.5 coh 25000
;**************************边界条件****************************
fix x range x 4.9 5.1
fix x range x -0.1 0.1
fix y range y -0.1 0.1
fix y range y 4.9 5.1
fix z range z -21.1 -20.9
set grav 0 0 -10
;**************************初始条件****************************
ini dens 1900 range group ft1
ini dens 1900 range group fznt2
ini dens 1900 range group fs3
ini dens 1900 range group fs4
ini dens 2500 range group zhuang
solve
plot
add surface red
show
