pile 单元模拟矩形抗滑桩，开挖桩前土体，加shell，计算不了。

图中红色是2个pile抗滑桩，加桩平衡后，开挖桩前土体，加shell防护，再计算出现： *** Overflow error in calculation!
这是怎么回事？

new

gen z b p0 0 0 0 p1 20 0 0 p2 0 12 0 p3 0 0 34 p4 20 12 0 size 20 12 34

gen z b p0 20 0 0 p1 23 0 0 p2 20 12 0 p3 20 0 34 p4 23 12 0 size 3 12 34

gen z b p0 23 0 0 p1 63 0 0 p2 23 12 0 p3 23 0 34 p4 63 12 0 size 40 12 34

gen z b p0 23 0 34 p1 63 0 34 p2 23 12 34 p3 35 0 44 p4 63 12 34 p5 35 12 44 p6 63 0 44 p7 63 12 44 size 40 12 10

pause k

gen merge 1e-6

attach face

model mohr

prop bulk 15e6 shear 5e6 fric 22 coh 1e10 tension 1e10 poisson 0.35

ini dens 1800

fix x y z range z -0.1 0.1

fix y range y -0.1 0.1 

fix y range y 11.9 12.1 

fix x range x -0.1 0.1 

fix x range x 62.9 63.1

set gravity 0 0 -9.8

;pause

solve

ini xdisp 0 ydisp 0 zdisp 0

ini xvel 0 yvel 0 zvel 0

prop coh 30000 tens 0

solve

save tt0-1.sav

ini xdisp 0 ydisp 0 zdisp 0

ini xvel 0 yvel 0 zvel 0

sel pile id 1 begin 20 3 8 end 20 3 34 nseg 26

sel pile id 2 begin 20 9 8 end 20 9 34 nseg 26

sel pile pro d 2500 emod 3e11 nu 0.2 xcarea 6.0 xciy 4.5 xciz 2.0 xcj 6.5 per 10 cs_scoh 21000 cs_sfric 15 cs_sk 2.5e8 cs_ncoh 21000 cs_nfric 15 cs_nk 2.5e8 cs_ng=on slide=on

solve

save tt-st1.sav

ini xdisp 0 ydisp 0 zdisp 0

ini xvel 0 yvel 0 zvel 0

hist unbal

model null range x 0 20 y 0 12 z 26 34

sel shell id 7 range x 20 y 0 12 z 26 34

sel shell id 7 prop iso=(5e10,0.25) thick=0.5

sel join pile id 1 shell id 7

sel join pile id 2 shell id 7

;set large

solve

save tt-st3.sav