求助！计算不收敛问题

Iris-tectroum

求助！大神帮帮看看这个命令，一直都收敛不了

n

ge zone b size 5 4 1 p0 0 0 18 p1 5 0 18 p2 0 3.95 18 p3 0 0 18.15 gr back-fillingA

ge zone b size 5 4 1 p0 0 0 18.15 p1 5 0 18.15 p2 0 3.95 18.15 p3 0 0 18.35 gr back-fillingA

ge zone b size 5 4 1 p0 0 0 18.35 p1 5 0 18.35 p2 0 3.95 18.35 p3 0 0 18.42 gr back-fillingA

ge zone b size 5 6 1 p0 0 3.95 18 p1 5 3.95 18 p2 0 7.45 18 p3 0 3.95 18.15 gr dianceng

ge zone b size 5 12 1 ratio 1 1.2 1 p0 0 7.45 18 p1 5 7.45 18 p2 0 22.5 18 p3 0 7.45 18.15 gr dianceng

ge zone b size 5 6 1 p0 0 3.95 18.15 p1 5 3.95 18.15 p2 0 7.45 18.15 p3 0 3.95 18.35 gr jiceng

ge zone b size 5 12 1 ratio 1 1.2 1 p0 0 7.45 18.15 p1 5 7.45 18.15 p2 0 22.5 18.15 p3 0 7.45 18.35 gr jiceng

ge zone b size 5 6 1 p0 0 3.95 18.35 p1 5 3.95 18.35 p2 0 7.45 18.35 p3 0 3.95 18.42 gr mianceng

ge zone b size 5 12 1 ratio 1 1.2 1 p0 0 7.45 18.35 p1 5 7.45 18.35 p2 0 22.5 18.35 p3 0 7.45 18.42 gr mianceng

ge zone b size 5 4 2 p0 0 0 15.5 p1 5 0 15.5 p2 0 3.95 15.5 p3 0 0 18 gr cushion

ge zone b size 5 6 2 p0 0 3.95 15.5 p1 5 3.95 15.5 p2 0 7.45 15.5 p3 0 3.95 18 p4 5 7.45 15.5 &

p5 0 4.95 18 p6 5 3.95 18 p7 5 4.95 18 gr back-fillingB

ge zone b size 5 12 2 ratio 1 1.2 1 p0 0 7.45 15.5 p1 5 7.45 15.5 p2 0 22.5 15.5 p3 0 4.95 18 p4 5 22.5 15.5 &

p5 0 22.5 18 p6 5 4.95 18 p7 5 22.5 18 gr plain-fill

attach face ra y 3.94 3.96 z 15.5 18

ge zone b size 5 4 2 p0 0 0 13.5 p1 5 0 13.5 p2 0 3.95 13.5 p3 0 0 15.5 gr back-fillingC

ge zone b size 5 6 2 p0 0 3.95 13.5 p1 5 3.95 13.5 p2 0 4.95 13.5 p3 0 3.95 15.5 gr back-fillingC

ge zone b size 5 12 2 ratio 1 1.2 1 p0 0 4.95 13.5 p1 5 4.95 13.5 p2 0 22.5 13.5 p3 0 4.95 15.5 gr back-fillingC

ge zone b size 5 4 2 p0 0 0 9 p1 5 0 9 p2 0 3.95 9 p3 0 0 13.5 gr clayA

ge zone b size 5 6 2 p0 0 3.95 9 p1 5 3.95 9 p2 0 4.95 9 p3 0 3.95 13.5 gr clayA

ge zone b size 5 12 2 ratio 1 1.2 1 p0 0 4.95 9 p1 5 4.95 9 p2 0 22.5 9 p3 0 4.95 13.5 gr clayA

ge zone b size 5 4 2 p0 0 0 2 p1 5 0 2 p2 0 3.95 2 p3 0 0 9 gr clayB

ge zone b size 5 6 2 p0 0 3.95 2 p1 5 3.95 2 p2 0 4.95 2 p3 0 3.95 9 gr clayB

ge zone b size 5 12 2 ratio 1 1.2 1 p0 0 4.95 2 p1 5 4.95 2 p2 0 22.5 2 p3 0 4.95 9 gr clayB

ge zone b size 5 4 2 p0 0 0 0 p1 5 0 0 p2 0 3.95 0 p3 0 0 2 gr clayC

ge zone b size 5 6 2 p0 0 3.95 0 p1 5 3.95 0 p2 0 4.95 0 p3 0 3.95 2 gr clayC

ge zone b size 5 12 2 ratio 1 1.2 1 p0 0 4.95 0 p1 5 4.95 0 p2 0 22.5 0 p3 0 4.95 2 gr clayC

fix x y z ra z -0.1 0.1

fix y ra y -0.1 0.1

fix y ra y 22.4 22.6

fix x

model mohr

pro shear 7e6 bulk 1.1e7 coh 2.2e4 fric 20 dens 1930 ra gr plain-fill

pro shear 5e6 bulk 1e7 coh 2e4 fric 32 dens 2000 ra gr back-fillingA

pro shear 5e6 bulk 1e7 coh 2e4 fric 32 dens 2000 ra gr back-fillingB

pro shear 7e6 bulk 1.2e7 coh 1.3e4 fric 16.7 dens 1700 ra gr back-fillingC

pro shear 4e6 bulk 6e7 coh 1e4 fric 12 dens 1790 ra gr clayA

pro shear 3e6 bulk 6e6 coh 1.1e4 fric 12 dens 1700 ra gr clayB

pro shear 1.1e7 bulk 1.9e7 coh 5e4 fric 31.9 dens 1890 ra gr clayC

pro shear 2.31e8 bulk 5e8 coh 4.5e5 fric 30 dens 1800 ra gr dianceng

pro shear 1.2e9 bulk 2e9 coh 5e5 fric 35 dens 2200 ra gr jiceng

model e ra gr cushion any gr mianceng

prop shear 1.4e7 bulk 1.6e7 poisson 0.167 ra gr cushion

ini dens 2500 ra gr cushion

prop shear 4.8e8 bulk 8e8 poisson 0.25 dens 2500 ra gr mianceng

set grav 0 0 -9.81

so

ini xdis 0 ydis 0 zdis 0

ini xvel 0 yvel 0 zvel 0

ini state 0

his id 1 gp zdisp 2.5 3.95 18

his id 2 gp zdisp 2.5 4.95 18

his id 3 gp zdisp 2.5 6 18

his id 4 gp zdisp 2.5 7.49 18

his id 5 gp zdisp 2.5 9.5 18

his id 6 gp zdisp 2.5 11.5 18

his id 7 gp zdisp 2.5 13.5 18

his id 8 gp zdisp 2.5 15.45 18

his id 9 gp zdisp 2.5 17.5 18

his id 10 gp zdisp 2.5 19.5 18

his id 11 gp zdisp 2.5 21.5 18

his id 12 gp ydisp 2.5 3.95 18.42

his id 13 gp ydisp 2.5 4.95 18.42

his id 14 gp ydisp 2.5 6 18.42

his id 15 gp ydisp 2.5 7.49 18.42

his id 16 gp ydisp 2.5 9.5 18.42

his id 17 gp ydisp 2.5 11.5 18.42

his id 18 gp ydisp 2.5 13.5 18.42

his id 19 gp ydisp 2.5 15.45 18.42

his id 20 gp ydisp 2.5 17.5 18.42

his id 21 gp ydisp 2.5 19.5 18.42

his id 22 gp ydisp 2.5 21.5 18.42

apply szz -0.7e5 ra x 1 1.5 y 4.97 4.99 z 18.41 18.42

apply szz -0.7e5 ra x 3.5 4 y 4.97 4.99 z 18.41 18.42

apply szz -0.7e5 ra x 1 1.5 y 6.47 6.67 z 18.41 18.42

apply szz -0.7e5 ra x 3.5 4 y 6.47 6.67 z 18.41 18.42

apply szz -0.7e5 ra x 1 1.5 y 8.71 8.91 z 18.41 18.42

apply szz -0.7e5 ra x 3.5 4 y 8.71 8.91 z 18.41 18.42

apply szz -0.7e5 ra x 1 1.5 y 10.21 10.41 z 18.41 18.42

apply szz -0.7e5 ra x 3.5 4 y 10.21 10.41 z 18.41 18.42

apply szz -0.7e5 ra x 1 1.5 y 12.45 12.65 z 18.41 18.42

apply szz -0.7e5 ra x 3.5 4 y 12.45 12.65 z 18.41 18.42

apply szz -0.7e5 ra x 1 1.5 y 13.95 14.15 z 18.41 18.42

apply szz -0.7e5 ra x 3.5 4 y 13.95 14.15 z 18.41 18.42

solve