FLAC3D计算隧道作业Word格式.docx

FLAC3D计算隧道作业Word格式.docx

《FLAC3D计算隧道作业Word格式.docx》由会员分享，可在线阅读，更多相关《FLAC3D计算隧道作业Word格式.docx（23页珍藏版）》请在冰豆网上搜索。

右上半圆衬砌

genzoneradtunp0000p100-30p2010p33000size1011030dim5555ratio1111groupoutsiderock;

右下矩形部分围岩

genzoneradtunp0000p100-5p2010p3500size101101dim4.94.94.94.9ratio1111groupconcretliner;

右下矩形部分衬砌

genzonebrickp000-4.9p1add4.900p2add010p3add004.9size10110ratio1111groupinsiderock;

隧道内部

genzonreflectdip90dd90orig000;

关于z轴对称

plotblockgroup

plotaddaxesblack

savejianmo.sav

（1）自重应力场作用下隧道开挖后的拱顶下沉和边墙水平收敛大小以及围岩中的塑性区大小。

围岩选用M-C模型，参数的选取：

岩体密度（kg/m3）=2200,重力加速度（kg/N）=10,体积模量（Pa）=5e10，剪切模量（Pa）=1.1e10，摩擦角=30o，粘聚力（N）=1.8e7，抗拉强度（N/m2）=1e7。

在Z=30m平面处施加的应力为γh=22000×

（500-30）=10.34e6N/m2。

命令流如下；

;

边界条件

fixzrangez-30.01-29.99

fixxrangex-30.01-29.99

fixxrangex29.9930.01

fixyrangey-0.010.01

fixyrangey0.991.01

applyszz-10.34e6rangez29.9930.01;

设置初始应力状态

modelmohr

setgravity00-10;

设置重力

inidensity2200;

围岩的密度

propbulk=5e10shear=1.1e10fric=30coh=0.18e7ten=1e7

;

体积模量、剪切模量、摩擦角、粘聚力、抗拉强度

setmechratio=1e-4

plotblockstate

solve;

求解自重应力场

saveGravsol.sav

plotblockgroup

plotaddaxesred

毛洞开挖计算

initialxdisp=0ydisp=0zdisp=0

modelnullrangegroupinsiderockanygroupconcretlinerany

setmechratio=1e-4

histunb

histgpzdisp005

histgpxdisp50-2

plotblockstate

savemaodong0.0.sav

计算结果如下图所示：

未开挖时的隧道开挖后毛洞塑性区

开挖毛洞后z方向上的位移云图开挖毛洞后x方向的位移云图

开挖毛洞后Z方向的位移曲线开挖毛洞后X方向的位移曲线

开挖毛洞后z方向的应力开挖毛洞后x方向的应力

可以看到开挖毛洞后未支护的情况下，拱顶下沉4.775mm，边墙水平收敛2.996mm，围岩塑性区如上图所示。

（2）当侧压力系数为0.5—2.5时的影响

侧压力系数为0.5时，命令流如下：

iniszz-11e6grad0022000;

inisxx-5.5e6grad0011000

inisyy-5.5e6grad0011000

给定边界条件

fixxrangex-30.01-29.99

fixxrangex29.9930.01

fixzrangez-30.01-29.99

fixzrangez29.9930.01

savemaodong0.5.sav

不同的侧压力系数的命令流，只需修改不同的初始应力状态的命令即可，在此对侧压力系数分别为1.0,1.5,2.0,2.5的情况，不再重复添加命令流。

侧压力系数不同开挖后的塑性区如下图所示：

λ=0.5λ=1.0

λ=1.5λ=2.0

λ=2.5

随着侧压力系数的增大，隧道顶部和边墙的塑性区都在逐渐扩大，且拱顶扩大的更多，当侧压力系数为2.5时，塑性区急剧扩大而不能自稳。

不同的侧压力系数开挖后Z方向的位移云图如下：

λ=0.5λ=1.0

不同的侧压力系数拱顶处Z方向的位移曲线如下：

λ=1.5λ=2.0

不同的侧压力系数X方向的位移云图如下:

不同的侧压力系数X方向的位移曲线如下:

不同的侧压力系数X方向的应力云图如下：

不同的侧压力系数Z方向应力云图如下：

侧压力系数λ

0.5

1.0

1.5

2.0

2.5

拱顶竖直位移（mm）

2.861

2.988

4.217

6.222

8.142

边墙水平位移（mm）

2.052

3.253

4.949

7.200

9.963

可以看到随着侧压力系数的增大，拱顶竖直位移和边墙水平位移都相应增加，且水平位移增加的速度比竖直位移增加的速度快。

塑性区随着侧压力系数的增大而增加，且拱顶处的塑性区增加的多。

这说明水平的构造应力对边墙的水平位移和拱顶的塑性区影响较大。

（3）开挖隧道后采用锚喷支护：

命令流如下：

restoreGravsol.sav

step100;

计算开挖100步后进行支护

模筑衬砌计算

modelelasrangegroupconcretlinerany

inidensity2500rangegroupconcretlinerany;

衬砌混凝土的密度

propbulk=26.67e10,shear=22.5e10rangegroupconcretlinerany

衬砌混凝土的体积弹模、剪切弹模

加锚杆

selcableid=1begin-500end-800nseg=4

selcableid=2begin-4.755301.5451end-7.608502.4721nseg=4

selcableid=3begin-4.045102.9389end-6.472104.7023nseg=4

selcableid=4begin-2.938904.0451end-4.702306.4721nseg=4

selcableid=5begin-1.545104.7553end-2.472107.6085nseg=4

selcableid=6begin005end008nseg=4

selcableid=7begin1.545104.7553end2.472107.6085nseg=4

selcableid=8begin2.938904.0451end4.702306.4721nseg=4

selcableid=9begin4.045102.9389end6.472104.7023nseg=4

selcableid=10begin4.755301.5451end7.608502.4721nseg=4

selcableid=11begin500end800nseg=4

selcableid=12begin50-1.25end80-1.25nseg=4

selcableid=13begin50-2.5end80-2.5nseg=4

selcableid=14begin50-3.75end80-3.75nseg=4

selcableid=15begin50-5end80-5nseg=4

selcableid=16begin-50-1.25end-80-1.25nseg=4

selcableid=17begin-50-2.5end-80-2.5nseg=4

selcableid=18begin-50-3.75end-80-3.75nseg=4

selcableid=19begin-50-5end-80-5nseg=4

锚杆力学参数

selcablepropxcarea=49e-4emod=45e10ytens=28e9gr_k=1.5e8gr_coh=16e6

支护后计算结果如下图：

支护前塑性区支护后塑性区

锚杆应力锚杆与塑性区的相对位置

支护后X方向位移云图支护后z方向位移云图

支护后X方向应力云图支护后Z方向应力云图

支护后Z方向位移曲线支护后X方向位移曲线

支护前

拱顶竖直位移（mm）X0

边墙水平位移（mm）Z0

支护后

拱顶竖直位移（mm）X1

边墙水平位移（mm）Z1

3.181

4.775

2.038

3.172

从图中可以看到，支护后塑性区明显比未支护的情况小了很多，支护后围岩的受力情况得到改善，竖直位移和水平位移都减少了。

在位移曲线中可以明显的看到一个转折点，这是在施加初期支护后对围岩变形的限制，阻止了围岩的进一步变形所致。