请大家看看，怎么平衡不了啊

zxf613

new
gen zone radcyl p0(0,0,0) p1(20,0,0) p2(0,0,-5) p3(0,10,0) p4(20,0,-5) p5(0,10,-5) p6(20,10,0) p7(20,10,-5) p8(1,0,0) p9(0,1,0) p10(1,0,-5) p11(0,1,-5) &
                 size 10 5 5 5  ratio 1 1 1 1
gen zone radcyl p0(0,0,-5) p1(20,0,-5) p2(0,0,-10) p3(0,10,-5) p4(20,0,-10) p5(0,10,-10) p6(20,10,-5) p7(20,10,-10) p8(1,0,-5) p9(0,1,-5) p10(1,0,-10) p11(0,1,-10) &
                size 10 5 5 5  ratio 1 1 1 1 fill
interface 1 face range cylinder end1 (0,0,0)  end2 (0,0,-5.1) radius 1.1 cylinder end1 (0,0,0)  end2 (0,0,-5.1) radius 0.9 not
interface 2 face range cylinder end1 (0,0,-4.9)  end2 (0,0,-5.1) radius 1.1
int 1 maxedge 0.55
int 2 maxedge 0.55
gen zone cyl p0 (0,0,8) p1 (1,0,8) p2 (0,0,1)  p3 (0,1,8)  p4 (1,0,1) p5 (0,1,1) size 3 10 6
gen zone cyl p0 (0,0,8.1) p1 (1,0,8.1) p2 (0,0,8)  p3 (0,1,8.1)  p4 (1,0,8) p5 (0,1,8) size 3 1 6
group pile range z 1 8.1
ini z add -6.0 range group pile
gen zone reflec dip 90 dd 0
gen zone reflec dip 90 dd 90
group 1 range z 0 -1
group 2 range z -1 -9
group 3 range z -9 -10
group pile1 range group pile z -5 2
group pile2 range group pile z 2 2.1
;边界条件
fix x y z ran z -9.9 -10.1
fix x ran x -19.9 -20.1
fix x ran x  19.9  20.1
fix y ran y  -9.9 -10.1
fix y ran y   9.9  10.1
save 1网格.sav
pl bl gr
rest 1网格.sav
config fluid dynamic
set fluid off dynamic off
; Mechanical model and paras
model mohr ran gro 1 any gro 2 any gro 3 any

prop bulk 14.71e6 shear 5.640e6 c 10e3 f 25 ran gro 1
prop bulk 2.941e7 shear 1.128e7 c 0    f 30 ran gro 2
prop bulk 14.71e6 shear 5.640e6 c 15e3 f 25 ran gro 3

model elastic ran gro pile1
prop bulk 1.67e10 shear 7.69e9 ran gro pile1
model null ran gro pile2
ini dens 1040 ran gro 1
ini dens 1440 ran gro 2
ini dens 1100 ran gro 3
ini dens 2400 ran gro pile1
; Interface paras
int 1 prop kn 2e9 ks 2e9 fric 15
int 2 prop kn 2e9 ks 2e9 fric 15
set grav 10
set mech rat 1e-6
solve
save 2初始应力.sav