计算2000步保证平衡吗？
new;1/4 上半部分
gen zon radcyl
p0 15.7 0 0 p1 24.5 0 0 p2 15.7 1 0 p3 15.7 0 8.8 &
size 5 2 10 4 dim 5.8 5.8 5.8 5.8 rat 1 1 1 1.2 group outrockys
gen zone cshell p0 15.7 0 0 p1 21.5 0 0 p2 15.7 1 0 p3 15.7 0 5.8 &
size 1 2 10 4 dim 5.4 5.4 5.4 5.4 rat 1 1 1 1 group liningys fill group inrockys
gen zon reflect dip 90 dd 270 orig 15.7 0 0
;1/4 下半部分
gen zon radcyl
p0 15.7 0 0 p1 15.7 0 -6.9
p2 15.7 1 0 p3 24.5 0 0 &
size 5 2 10 4 dim 3.9 5.8 3.9 5.8 rat 1 1 1 1.2 group outrockyx
gen zone cshell p0 15.7 0 0 p1 15.7 0 -3.9 p2 15.7 1 0 p3 21.5 0 0 &
size 1 2 10 4 dim 3.5 5.4 3.5 5.4 rat 1 1 1 1 group liningyx fill group inrockyx
gen zon reflect dip 90 dd 270 orig 15.7 0 0 range group liningyx any &
group inrockyx any group outrockyx any
;右隧道中地层网格
gen zon brick
p0 0 0 -6.9
p1 6.9 0 -6.9
p2 0 1 -6.9
p3 0 0 0
size 6 2 5
group outrock1
gen zon brick
p0 24.5 0 -6.9
p1 65 0 -6.9
p2 24.5 1 -6.9
p3 24.5 0 0
size 20 2 5
group outrock1
gen zon brick
p0 0 0 0
p1 6.9 0 0
p2 0 1 0
p3 0 0 8.8
size 6 2 5
group outrock2
gen zon brick
p0 24.5 0 0
p1 65 0 0
p2 24.5 1 0
p3 24.5 0 8.8
size 20 2 5
group outrock2
;右隧道上侧土体网格
gen zon brick
p0 0 0 8.8
p1 6.9 0 8.8
p2 0 1 8.8
p3 0 0 36
size 6 2 10
group outrock3
gen zon brick
p0 24.5 0 8.8
p1 65 0 8.8
p2 24.5 1 8.8
p3 24.5 0 36
size 20 2 10
group outrock3
gen zon brick
p0 6.9 0 8.8
p1 24.5 0 8.8
p2 6.9 1 8.8
p3 6.9 0 36
size 10 2 10
group outrock3
;右隧道下侧土体网格
gen zon brick
p0 0 0 -24
p1 6.9 0 -24
p2 0 1 -24
p3 0 0 -6.9
size 6 2 8
group outrock4
gen zon brick
p0 24.5 0 -24
p1 65 0 -24
p2 24.5 1 -24
p3 24.5 0 -6.9
size 20 2 8
groupoutrock4
gen zon brick
p0 6.9 0 -24
p1 24.5 0 -24
p2 6.9 1 -24
p3 6.9 0 -6.9
size 10 2 8
group outrock4
gen zon reflect dip 90 dd 270 orig 0 0 0
plot add axes red
plot block group
;设置重力加速度
set gravity 0 0 -10
;给定边界条件
fix z range z -24.01,-23.99
fix x range x -65.01,-64.99
fix x range x 64.99,65.01
fix y range y -0.01 0.01
fix y range y 0.99,1.01
;0-求解自重应力场
model mohr
ini density 2200
;围岩的密度
prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 tension=5e5
;体积、剪切、摩擦角、凝聚力、抗拉强度
step 6000
save Gsol.sav
plot cont zdisp
plot cont szz
;1-左隧道上半断面开挖
rest Gsol.sav
ini xdis 0.0 ydis 0.0 zdis 0.0
ini density 2400 range group outrockzs any
prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 &
range group outrockzs any
model null range group inrockzs any group liningzs any
step 2000
save stepp1.sav
;2-初期支护
计算2000步保证平衡吗？这个有什么依据，欢迎大家讨论
页:
[1]