- 积分
- 0
- 注册时间
- 2009-12-29
- 仿真币
-
- 最后登录
- 1970-1-1
|
大家好,请求各位好心人,帮帮忙,很急,谢谢!求:李围老师“隧道及地下工程FLAC解析方法”地铁隧道开挖模拟遇到的重要问题,谢谢!问题如下:
1.、三维模型建立后,地铁隧道要分左隧道施工和右隧道施工,在已编好的程序中,为什么左隧道第一步开挖程序运行不了?程序如下:
new
gen zon radcyl
p0 9 0 0 p1 18 0 0 p2 9 36 0 p3 9 0 9 &
size 4 24 8 4 dim 3 3 3 3 rat 1 1 1 1.2
group outr
gen zone cshell p0 9 0 0 p1 12 0 0 p2 9 36 0 p3 9 0 3 &
size 1 24 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr
gen zon reflect dip 90 dd 90 orig 9 0 0
gen zon reflect dip 0
dd 0
ori 9 0 0
;1/2隧道地层部分
gen zon brick
p0 18 0 -9
p1 36 0 -9
p2 18 36 -9
p3 18 0 9
size 10 12 8 group outr1
gen zon brick
p0 18 0 9
p1 36 0 9
p2 18 36 9
p3 18 0 13
size 10 12 2 group outr2
gen zon brick
p0 0 0 9
p1 18 0 9
p2 0 36 9
p3 0 0 13
size 8 12 2 group outr2
gen zon brick
p0 18 0 -23
p1 36 0 -23
p2 18 36 -23
p3 18 0 -9
size 10 12 8 group outr3
gen zon brick
p0 0 0 -23
p1 18 0 -23
p2 0 36 -23
p3 0 0 -9
size 8 12 8 group outr3
gen zon reflect dip 90 dd 270 orig 0 0 0
group linz range group liny x -12.1 -5.99
;绘制模型图
plot add axes red
plot set rotation 20 0 45
plot block group
;0-自重应力场求解
set gravity 0 0 -10
;给定边界条件
fix z range z -23.01,-22.99
fix x range x -36.01,-35.99
fix x range x 35.99,36.01
fix y range y -0.01 0.01
fix y range y 35.99,36.01
model mohr
;围岩的密度
ini density 1800
;体积、剪切、摩擦角、凝聚力、抗拉强度
prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4
step 3000
save Gsol.sav
左隧道施工程序:
;第1步
rest Gsol.sav
initial xdisp=0 ydisp=0 zdisp=0
model null range group inrz y 0 6(此处为什么运行错误,要使其与后面的右隧道开挖顺心进行,如何改进此处,谢谢)
model elas range group linz y 0 6
ini density 2500
range group linz y 0 6
prop bulk=16.67e9,shear=12.5e9 range group linz y 0 6
step 3000
save
step1.sav
;第2步
save
step2.sav
。。。。。。。。
右隧道施工程序:
;第1步
rest step6.sav
model null range group inry y 0 6
model elas range group liny y 0 6(此处也是和上面左隧道一样,出现同样问题怎么解决)
ini density 2500
range group liny y 0 6
prop bulk=16.67e9,shear=12.5e9 range group liny y 0 6
step 3000
save
step7.sav
;第2步
。。。。。。
望各位好心人帮我解决这个问题,使这个程序能顺利运行,本人急用,谢谢大家了,谢谢!祝大家工作顺利,万事如意! |
|