rt，建模如下，为什么solve不收敛呢
new
;衬砌管片
gen zone cshell p0 0,0,0 p1 3,0,0 p2 0,48,0 p3 0,0,3 dim 2.7,2.7,2.7,2.7 &
size 2,48,10,8 rat 1,1,1,1
gen zone cshell p0 0,0,0 p1 0,0,-3 p2 0,48,0 p3 3,0,0 dim 2.7,2.7,2.7,2.7 &
size 2,48,10,8 rat 1,1,1,1
group liner
;浆体
gen zone cshell p0 0,0,0 p1 3.15,0,0 p2 0,48,0 p3 0,0,3.15 &
dim 3,3,3,3 size 1,48,10,9 rat 1,1,1,1
gen zone cshell p0 0,0,0 p1 0,0,-3.15 p2 0,48,0 p3 3.15,0,0 &
dim 3,3,3,3 size 1,48,10,9 rat 1,1,1,1
group cement
;内部土体
gen zone cyl p0 0,0,0 p1 2.7,0,0 p2 0,48,0 p3 0,0,2.7 size 6,48,10 &
rat 1.2,1,1
gen zone cyl p0 0,0,0 p1 0,0,-2.7 p2 0,48,0 p3 2.7,0,0 size 6,48,10 &
rat 1.2,1,1
group 'inner'
;隧道外层上部（L3）
gen zone radcyl p0 0,0,0 p1 10,0,0 p2 0,48,0 p3 0,0,8.49 &
dim 3.15,3.15,3.15,3.15 size 9,48,10,10 rat 1,1,1,1.2
;隧道外层下部（L3）
gen zone radcyl p0 0,0,0 p1 0,0,-8.24 p2 0,48,0 p3 10,0,0 &
dim 3.15,3.15,3.15,3.15 size 9,48,10,10 rat 1,1,1,1.2
gen zone brick p0 10,0,8.49 p1 10,0,0 p2 10,48,8.49 p3 32,0,8.49 &
size 5,48,10
gen zone brick p0 10,0,0 p1 10,0,-8.24 p2 0,48,0 p3 32,0,0 &
size 5,48,10
group layer3
;表层土（L1）
gen zone brick p0 0,0,18.8 p1 0,0,16.39 p2 0,48,18.8 p3 10,0,18.8 &
size 1,48,5
gen zone brick p0 10,0,18.8 p1 10,0,16.39 p2 10,48,18.8 p3 32,0,18.8 &
size 1,48,10
group layer1
;第二层（L2）
gen zone brick p0 0,0,16.39 p1 0,0,8.49 p2 0,48,16.39 p3 10,0,16.39 &
size 4,48,5 rat 0.8,1,1
gen zone brick p0 10,0,16.39 p1 10,0,8.49 p2 10,48,16.39 p3 32,0,16.39 &
size 4,48,10 rat 0.8,1,1
group layer2
;第四层（L4）
gen zone brick p0 0,0,-8.24 p1 0,0,-17.2 p2 0,48,-8.24 p3 10,0,-8.24 &
size 4,48,5 rat 1.2,1,1
gen zone brick p0 10,0,-8.24 p1 10,0,-17.2 p2 10,48,-8.24 p3 32,0,-8.24 &
size 4,48,10 rat 1.2,1,1
group layer4
;材料属性
mo mohr
prop bulk 22.2e6 shear 7.41e6 coh 10e3 fric 12.5 range z 16.39 18.8
prop bulk 21.92e6 shear 10.11e6 coh 20.1e3 fric 26.5 range z 8.49 16.39
prop bulk 36.75e6 shear 16.96e6 coh 36.3e3 fric 32.2 range z -8.24 8.49
prop bulk 22e6 shear 13.2e6 coh 19.5e3 fric 25.1 range z -17.2 -8.24
;边界条件
fix x range x -.1 .1
fix x range x 31.9 32.1
fix y range y -.1 .1
fix y range y 47.9 48.1
fix z range z -17.3 -17.1
;初始应力
set grav 0,0,-10
ini density 1620 range z 16.39 18.8
ini density 1940 range z 8.49 16.39
ini density 2010 range z -8.24 8.49
ini density 1910 range z -17.2 -8.24
ini szz -304560 grad 0,0,16200 range z 16.39 18.8
ini sxx -163993.8 grad 0,0,8723.08 range z 16.39 18.8
ini syy -163993.8 grad 0,0,8723.08 range z 16.39 18.8
ini szz -357008 grad 0,0,19400 range z 8.49 16.39
ini sxx -153003.4 grad 0,0,8314.29 range z 8.49 16.39
ini syy -153003.4 grad 0,0,8314.29 range z 8.49 16.39
ini szz -362951 grad 0,0,20100 range z -8.24 8.49
ini sxx -155550.4 grad 0,0,8614.29 range z -8.24 8.49
ini syy -155550.4 grad 0,0,8614.29 range z -8.24 8.49
ini szz -371191 grad 0,0,19100 range z -17.2 -8.24
ini sxx -123730.3 grad 0,0,6366.67 range z -17.2 -8.24
ini syy -123730.3 grad 0,0,6366.67 range z -17.2 -8.24
hist unbal
hist gp zdisp 0,0,18.8
hist gp zdisp 0,0,3.15
set large
step 3000
plot add surface green
plot show
