各位大侠，看看我建模之后的运算，为什么硐室周围的点的位移如此之大，是施加的力出问题了，还是什么？
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
PS:我弄得y轴和z轴互换了，对结果没影响，只是和另一个例子作对比
