抗滑桩计算不收敛

陈夏雨

命令流如下，计算了一万多步都不收敛。。不晓得是怎么回事
我觉得是接触面参数设置问题
哪位大侠帮我看下，这个接触面参数怎么设置啊，有什么要求啊
new
ge zo b p0  -24 0 0   p1  0 0 0   p2  -24 32 0   p3  -24 0 32  p4  0 32 0 &
p5  -24 32 32  p6  0 0 11 p7  0 32 11  size 8  16  11  gro 1
ge zo b p0 -24 0 32  p1 0 0 11 p2 -24 32 32  p3 -24 0 48   p4 0 32 11 p5 -24 32 48  p6 0 0 27  p7 0 32 27 size 8  16 16  gro 2
;del  ran  x  -3  0  y  10.5  12.5
;del  ran  x  -3  0  y  22.5  24.5
gro 1_1 ran gro 1 x  -3  0  y  10  12
gro 1_2 ran gro 1 x  -3  0  y  22  24
gro 2_1 ran gro 2 x  -3  0  y  10  12
gro 2_2 ran gro 2 x  -3  0  y  22 24
ge zo b p0 -3 10 29.43 p1 0 10 27 p2 -3 12 29.43 p3 -3 10 32 p4 0 12 27 p5 -3 12 32 p6 0 10 32 p7 0 12 32 size 1 1 6 gro 00
ge zo b p0 -3 22 29.43 p1 0 22 27 p2 -3 24 29.43 p3 -3 22 32 p4 0 24  27 p5 -3 24 32 p6 0 22 32 p7 0 24 32 size 1 1 6 gro 00
gro 3 ran gro 00 any gro 2_1 any gro 2_2  any
gro 4 ran  gro 1_1  any gro 1_2  any
ge merge 0.2
ge separate 1
int 1 wrap 1 2
Int  1  maxedge 0.55
ge separate 3
int 2 wrap 2 3
Int  2  maxedge 0.55
ge separate 4
int 3 wrap 1 4
Int  3  maxedge 0.55

pl bl gro
pl add int yel

model mohr  range  group  1
model mohr  range    group  2
prop  bulk 1.515e6   shear  7.813e5   coh 400   fric  38   range   group  1
ini dens 2600 range     group  1
prop  bulk 1.818e4   shear  9.375e3   coh 14   fric 18   range   group  2
ini dens 1900  range     group  2

model elas  range group  3
model elas  range group  4
prop   bulk   4.012e6   shear  1.082e6   range  group  3
ini  dens  2500  range    group  3
prop   bulk   4.012e6   shear  1.082e6   range  group  4
ini  dens  2500  range    group  4
int 1  prop  ks 2.57e7  kn  2.57  coh 7   fric 9
int 2  prop  ks 5.46e7  kn  5.46e7  coh  7   fric 9
int 3  prop  ks 5.46e7  kn  5.46e7  coh  200   fric  19

set  grav  0,0,-9.81
fix  z  range  z  -0.1  0.1  
fix  y  range  y  -0.1  0.1
fix  y  range  y  31.9    32.1   

hist  unbal
solve
save  pile1.sav