我要模拟的是软岩巷道的流变对底臌作用，用FLAC3D，问题转化为平面应变的形式，软岩分9层，圆拱巷道在同一种岩层内。蠕变时间和步长已有，但在运行时单元体畸形化，导致计算的停止。求高手指导一下，该怎么修改，万分感谢！
new
gen zone radcylinder p0 0 0 0 p1 5 0 0 p2 0 1 0 p3 0 0 5 size 5 1 10 5 dimension 2,2,2 ratio 1 1 1 1 fill group inner;中心圆拱部
gen zone radtunnel p0 0 0 0 p1 0 0 -5 p2 0 1 0 p3 5 0 0 size 5 1 5 5 dimension 2,2,2 ratio 1 1 1 1 fill group inner;中心矩形巷道部
gen zone brick p0 5 0 -5 p1 30 0 -5 p2 5 1 -5 p3 5 0 5 size 25 1 10 ratio 1 1 1;中心右边
gen zone brick p0 0 0 5 p1 30 0 5 p2 0 1 5 p3 0 0 40 size 30 1 30 ratio 1 1 1;中心上部
gen zone brick p0 0 0 -30 p1 30 0 -30 p2 0 1 -30 p3 0 0 -5 size 30 1 30 ratio 1 1 1;中心下部
gen zone reflect normal 1 0 0 origin 0 0 0
model cv
pro bu = 1.34e10 cohesion = 5.68e6 den = 2531 fric = 25 ksh = 3.285e10 kvis = 1.03e12 sh = 1.05e10 mvis = 1.11e11 tension = 1.4e6 range z 30 40
pro bu = 1.47e10 cohesion = 7.33e6 den = 2207 fric = 26.67 ksh = 3.341e10 kvis = 1.03e12 sh = 1.28e10 mvis = 1.11e11 tension = 1.37e6 range z 24 30
pro bu = 9.26e9 cohesion = 0.88e6 den = 1412 fric = 30.91 ksh = 3.178e10 kvis = 1.03e12 sh = 6.1e8 mvis = 1.11e11 tension = 0.33e6 range z 21 24
pro bu = 2.14e10 cohesion = 4.17e6 den = 2567 fric = 30.05 ksh = 3.399e10 kvis = 1.03e12 sh = 1.95e10 mvis = 1.11e11 tension = 1.37e6 range z 14 21
pro bu = 6.25e9 cohesion = 4.36e6 den = 2600 fric = 29.67 ksh = 3.313e10 kvis = 1.03e12 sh = 5.08e9 mvis = 1.11e11 tension = 1.21e6 range z 4 14
pro bu = 9.72e9 cohesion = 5.51e6 den = 2550 fric = 29.88 ksh = 3.102e10 kvis = 1.03e12 sh = 5.56e9 mvis = 1.11e11 tension = 1.24e6 range z -6 4
pro bu = 9.72e9 cohesion = 6.87e6 den = 2578 fric = 28.01 ksh = 3.102e10 kvis = 1.03e12 sh = 5.56e9 mvis = 1.11e11 tension = 1.64e6 range z -10 -6
pro bu = 2.19e10 cohesion = 7.07e6 den = 2400 fric = 26.03 ksh = 3.490e10 kvis = 1.03e12 sh = 2.23e10 mvis = 1.11e11 tension = 1.80e6 range z -20 -10
pro bu = 1.23e10 cohesion = 8.03e6 den = 2200 fric = 28.65 ksh = 3.285e10 kvis = 1.03e12 sh = 9.66e9 mvis = 1.11e11 tension = 1.64e6 range z -30 -20
fix x range x -30.1 -29.9
fix x range x 29.9 30.1
fix y
fix z range z -30.1 -29.9
ini szz -1.6e7 syy -1.6e7 sxx -1.6e7
ini xd 0 yd 0 zd 0
con creep
model null range group inner
set large
hist id 1 unbal
hist id 2 crtime
hist id 3 gp zdisp 0,0,-2
hist id 4 gp xdisp -2,0,-1
hist id 5 gp xdisp 2,0,-1
hist id 6 gp zdisp 0,0,2
plot add bcon szz
plot show
set creep on
set creep mindt = 1 maxdt = 100
set creep dt auto on
solve age 3600
