我建模，施加应力之后，显示洞室顶部的竖向位移过大，请问是怎么回事，是应力施加的问题，还是未加支护？还是什么？
PS：我把y轴和z轴互换了，为了和另一个例子作对比，不影响结果……
new
gen z radtunnel size 4 5 4 14 ratio 1 1 1 1.2 p0 50 50 0 p1 50 0 0 p2 50 50 50 p3 0 50 0 &
dim 10 10 10
gen z reflect normal 0 0 1 origin 0 0 50
gen z reflect normal 0 1 0 origin 50 50 50
gen z reflect normal 1 0 0 origin 50 0 0
plot bl gr
model mohr
prop shear 1e10 bulk 1.67e10 fric 37
ini den 2500
fix x y z range y -0.01 0.01
plot set rotation 20 0 30
hist gp ydisp (40 60 0)
hist gp ydisp (40 60 50)
hist gp ydisp (50 60 0)
hist gp ydisp (50 60 50)
hist gp ydisp (40 40 0)
hist gp ydisp (40 40 50)
hist gp ydisp (50 40 0)
hist gp ydisp (50 40 50)
hist unbal
;设置
ini syy -7.5e6 grad 0 7.5e4 0
ini szz -3.75e6 grad 0 3.75e4 0
ini sxx -3.75e6 grad 0 3.75e4 0
apply syy -7.5e6 range y -0.1 0.1 x 0 100 z 0 100
apply sxx -3.75e6 grad 0 3.75e4 0 range x -0.1 0.1 y 0 100 z 0 100
apply sxx -3.75e6 grad 0 3.75e4 0 range x 99.9 100.1 y 0 100 z 0 100
apply szz -3.75e6 grad 0 3.75e4 0 range z -0.1 0.1 x 0 100 y 0 100
apply szz -3.75e6 grad 0 3.75e4 0 range z 99.9 100.1 x 0 100 y 0 100
set grav 0 -10 0
;计算
step 3000
plot contour ydisp
plot add axes black
plot show
