！！！！求助，为什么我初始地应力后的沉降这么大！！！

hidekyo

C:\Users\hidekyo\Desktop\1.png
命令流：
new
gen zone radcyl p0 8.6 0 0 p1 12.2 0 0 p2 8.6 36 0 p3 8.6 0 3.6 &
    size 4 24 8 4 dim 3.1 3.1 3.1 3.1 rat 1 1 1 1.2 group outr5
gen zone cshell p0 8.6 0 0 p1 11.7 0 0 p2 8.6 36 0 p3 8.6 0 3.1 &
    size 1 24 8 4 dim 2.75 2.75 2.75 2.75 rat 1 1 1 1 group liny fill group inry
gen zone reflect dip 90 dd 90 orig 8.6 0 0
gen zone reflect dip 0 dd 0 orig 8.6 0 0
;1/2隧道地层部分
gen zone brick p0 12.2 0 -3.6 p1 36 0 -3.6 p2 12.2 36 -3.6 p3 12.2 0 3.6 size 30 12 8 group outr5
gen zone brick p0 0 0 -3.6 p1 5 0 -3.6 p2 0 36 -3.6 p3 0 0 3.6 size 6 12 8 group outr5
gen zone brick p0 12.2 0 -22 p1 36 0 -22 p2 12.2 36 -22 p3 12.2 0 -3.6 size 30 12 20 group outr5
gen zone brick p0 5 0 -22 p1 12.2 0 -22 p2 5 36 -22 p3 5 0 -3.6 size 8 12 20 group outr5
gen zone brick p0 0 0 -22 p1 5 0 -22 p2 0 36 -22 p3 0 0 -3.6 size 6 12 20 group outr5
gen zone brick p0 12.2 0 3.6 p1 36 0 3.6 p2 12.2 36 3.6 p3 12.2 0 21.6 size 30 12 20 group outr6
gen zone brick p0 5 0 3.6 p1 12.2 0 3.6 p2 5 36 3.6 p3 5 0 21.6 size 8 12 20 group outr6
gen zone brick p0 0 0 3.6 p1 5 0 3.6 p2 0 36 3.6 p3 0 0 21.6 size 6 12 20 group outr6
;完整隧道
gen zone reflect dip 90 dd 270 orig 0 0 0
;绘制模型图
plot add axes red
plot set rotation 20 0 45
plot block group
;===============================================================施工过程模拟
;自重应力场求解
fix z range z -22.1 -21.9
fix x range x -36.1 -35.9
fix x range x 35.9 36.1
fix y range y -0.1 0.1
fix y range y 35.9 36.1
model mohr
;绘制模型图
group linz range group liny x -11.8 -5.49
group inrz range group inry x -11.36 -5.84
group outr4 range group outr6 z 3.59 6.61
group outr3 range group outr6 z 6.59 12.61
group outr2 range group outr6 z 12.59 19.61
group outr1 range group outr6 z 19.59 21.6
plot add axes red
plot set rotation 20 0 45
plot block group
;围岩的密度
ini density 1800 range group outr1
ini density 1920 range group outr2
ini density 1900 range group outr3
ini density 1850 range group outr4
ini density 1800 range group outr5
;体积模量、剪切模量、摩擦角、凝聚力、抗拉强度
prop bulk=1.41e8 shear=5.1e7 fric=20 coh=3.5e4 tension=5.0e4 range group outr1
prop bulk=1.42e8 shear=5.2e7 fric=26 coh=3.0e4 tension=5.0e4 range group outr2
prop bulk=1.43e8 shear=5.3e7 fric=30 coh=3.7e4 tension=5.0e4 range group outr3
prop bulk=1.44e8 shear=5.4e7 fric=22 coh=4.0e4 tension=5.0e4 range group outr4
prop bulk=1.47e8 shear=5.6e7 fric=17 coh=4.5e4 tension=5.0e4 range group outr5
prop bulk=1.47e8 shear=5.6e7 fric=17 coh=4.5e4 tension=5.0e4 range group liny
prop bulk=1.47e8 shear=5.6e7 fric=17 coh=4.5e4 tension=5.0e4 range group linz
prop bulk=1.47e8 shear=5.6e7 fric=17 coh=4.5e4 tension=5.0e4 range group inry
prop bulk=1.47e8 shear=5.6e7 fric=17 coh=4.5e4 tension=5.0e4 range group inrz
ini szz -40e6 grad 0 0 9.2e5 range x -22 21.6
ini syy -20e6 grad 0 0 5.6e5 range y 0 36
ini sxx -20e6 grad 0 0 2.8e5 range z -36 36
set grav 0 0 -10
plot con zdis
step 30000
sav Gsol.sav

hidekyo

自己顶个
我做的是地铁盾构隧道的开挖模拟
的自重应力求解后的沉降就有好几十米
求高人指导下！！
我在赶论文!!!!!!

liqiangre

ini syy -20e6 grad 0 0 5.6e5 range y 0 36
ini sxx -20e6 grad 0 0 2.8e5 range z -36 36 0

这两个range是不是应该一样?

binggehaoyang

参数 有可能是