一个公路隧道开挖支护实例的讨论
set log onset logfile q.log
new
;全断面法
;generate primitive components of grid 生成最初网格
;upper tunnel
gen zon radcyl p0 0 0 0 p1 12 0 0 p2 0 20 0 p3 0 0 8 &
dim 8.75 5.75 8.75 5.75 size 8 10 12 5 rat 1 1 1 1 group 围岩1
gen zon cshell p0 0 0 0 p1 8.75 0 0 p2 0 20 0 p3 0 0 5.75 &
dim 8.5 5.5 8.5 5.5 size 3 10 12 6 rat 1 1 1 1 group 初支
gen zon cshell p0 0 0 0 p1 8.5 0 0 p2 0 20 0 p3 0 0 5.5 &
dim 8 5 8 5 size 4 10 12 8 rat 1 1 1 1 group 二衬 fill group 原岩
;lower tunnel
gen zon radcyl p0 0 0 0 p1 0 0 -6 p2 0 20 0 p3 12 0 0 &
dim 3.75 8.75 3.75 8.75 size 6 10 16 5 rat 1 1 1 1 group 围岩2
gen zon cshell p0 0 0 0 p1 0 0 -3.75 p2 0 20 0 p3 8.75 0 0 &
dim 3.5 8.5 3.5 8.5 size 3 10 16 6 rat 1 1 1 1 group 仰拱1
gen zon cshell p0 0 0 0 p1 0 0 -3.5 p2 0 20 0 p3 8.5 0 0 &
dim 3 8 3 8 size 4 10 16 8 rat 1 1 1 1 group 仰拱2 fill group 仰拱原岩
gen zone reflect normal -1 0 0
gen zone reflect dip 90 dd 270 origin -12 0 0
;边界网格
gen zon radtun p0 -12 0 0 p1 -12 0 -25 p2 -12 20 0 p3 40 0 0 &
dim 6 24 6 24 size 8 10 16 10 rat 1 1 1 1.1 group 围岩3
gen zon radtun p0 -12 0 0 p1 -64 0 0 p2 -12 20 0 p3 -12 0 -25 &
dim 24 6 24 6 size 16 10 8 10 rat 1 1 1 1.1 group 围岩3_1
gen zone radtun p0 -12 0 0 p1 40 0 0 p2 -12 20 0 p3 -12 0 25 &
p4 40 20 0 p5 -12 20 25 p6 40 0 15 p7 40 20 15 &
dim 24 8 24 8 size 12 10 6 10 rat 1 1 1 1.1 group 围岩4
gen zone radtun p0 -12 0 0 p1 -12 0 25 p2 -12 20 0 p3 -64 0 0 &
p4 -12 20 25 p5 -64 20 0 p6 -64 0 35 p7 -64 20 35 &
dim 8 24 8 24 size 6 10 12 10 rat 1 1 1 1.1 group 围岩4_1
;合并不贯通网格点
;gen merge 1e-2
;显示块体图
plo block group
;保存
save q0.sav
;paus
;设定摩尔库伦材料参数和密度
model mohr
pro bulk 11.8e6 she 2.5e6 fric 12 coh 46e3 ten 0 dil 0 &
density 1800 range z 25 35
pro bulk 11.1e6 she 3.7e6 fric 18 coh 23e3 ten 0 dil 5 &
density 1900 range z -5 25
pro bulk 2.08e8 she 1.19e8 fric 35 coh 45e3 ten 5e3 dil 15 &
density 2100 range z -25 -5
;给定边界条件
fix z range z -25.1 -24.9
fix x range x -64.1 -63.9
fix x range x 39.9 40.1
fix y range y -.1 .1
fix y range y 19.9 20.1;
;初始应力赋值
;apply szz -5e5 range z 19.9 20.1 深埋或有荷载时
;赋予重力
set grav 0 0 -10
;设置大应变
set large
;记录最大不平衡力
his unbal
hist gp zdis 0 5 8
hist gp zdis -24 5 8
hist gp xdis 8 5 0
hist gp xdis -8 5 0
hist gp xdis -16 5 0
hist gp xdis -32 5 0
;
;绘制位移云图
plo con zdisp outline on
;plo set mod fir pos -8 -200 0
;plo set ma=0.25
plo sho
;求解
solve
save q1.sav
;
;位移初始化(归零）后清除
ini xdisp=0 ydisp=0 zdisp=0
hist purge
;
;全断面法同时开挖两隧道
;exacavate the double tunnel
model null range group 初支
model null range group 二衬
model null range group 原岩
model null range group 仰拱1
model null range group 仰拱2
model null range group 仰拱原岩
;
save q2.sav
;适时支护
;绘制位移云图
plo con zdisp outline off
plo sho
;应力释放
step 10
;初期支护
model el range group 初支
model el range group 仰拱1
prop bulk 15.6e9 shear 11.7e9 dens 2300 range group 初支
prop bulk 15.6e9 shear 11.7e9 dens 2300 range group 仰拱1
set grav 10
set large
;
;绘制位移云图
plo con zdisp outline off
plo sho
;应力释放
step 10
save q4.sav
;
;二衬
model el range group 二衬
model el range group 仰拱2
prop bulk 16.4e9 shear 12.3e9 dens 2500 range group 二衬
prop bulk 16.4e9 shear 12.3e9 dens 2500 range group 仰拱2
;绘制位移云图
plo con zdisp outline off
plo sho
;plo his 1
;求解
set mec ratio=1.0e-4
solve
save q5.sav
cal t 问题1：
应力释放步数只能是10~15，
若step=50以上，隧道断面变形较大，
二衬加上去就出现网格单元错误，
这样下来没有时间效应，各种施工方法没法比较~~啊
请高手指点一下，怎么提高和控制释放应力步数，
而又能加上二衬~~不至于出现网格单元错误...
前提是软弱围岩条件下的...必须考虑大变形啊...
问题2：
加上衬砌后，模型收敛由初始的6000多步变成开挖支护后60000多步啊~~
为什么收敛变的这么缓慢呢?
[ 本帖最后由 wengoo21 于 2008-10-20 12:58 编辑 ] 我也在考虑这样的问题。不知道该怎么办啊
页:
[1]