1、1 工程概况研究此段为中条山隧道 K9+450K10+560 段，此处隧道最大埋深约 540 m，主要由太古界涑水群表壳岩组合解州片麻岩（Hgn）地层组成，构成中条山隧道分水岭北侧的主体；此段组成隧道的围岩岩性主要为变粒岩、花岗片麻岩等。该套地层岩性复杂，组合无规律。岩层产状整体倾向南东，倾角一般在50 70 间变化。在 AK9+900AK10+000 段为区域性断层影响段，此断层为破碎岩石组成，将为基岩裂隙水下渗提供通道，隧道开挖必将引起涌水，同时此段围岩稍差，施工时易引起坍塌。此段同时也是中条山北侧泉水主要涌出段，水文情况复杂。总体评价，本段工程地质条件差。在此处，具体运用 FLAC3D
2、进行模拟的区段均取洞身 YK10+100YK10+180段。该区段为 V 级围岩区域，埋深为 505512 m，为断层，附近太古界涑水群花岗片麻岩、黑云斜长片麻岩、岩石破碎。隧道断面为 SVc 型，如图 2-2 所示。图 2-2 SVc 型隧道断面图隧道衬砌按新奥法原理设计，采用 SVc 型复合式衬砌，该衬砌适用于隧道洞身 V 级断层影响带及软弱破碎围岩段的初期支护及衬砌，超前支护各环采用424 mm 注浆小导管超前预加固围岩，长 4.5 m，环向间距 35 cm，搭接长度1.3 m，斜插角 10 15 ，每环 37 根；初衬以喷、锚、网为主要支护手段：钢拱架为 I20a 型钢，纵向间距 75
3、 cm，每榀钢拱架之间采用 22 钢筋连接，环形间距 1.0 m；锚杆采用 D25 中空注浆锚杆，长 3.5 m，间距 75 cm（纵）100 cm（环） ，与钢拱架交错布置；喷 C25 早强混凝土 26 cm。二次衬砌和仰拱均为C30 钢筋混凝土结构，厚 50 cm。1.2 数值计算模型根据中条山隧道工程的实际状况，为提高计算速度，在保证计算精度的前提下，取桩号 YK10+100 YK10+160 段采用大型有限差分软件 FLAC3D 进行建模分析。对于全断面法、预留核心土法、台阶法，由于整个隧道模型左右对称，为减少计算量，可取隧道模型的一半计算，隧道的计算模型 I 如图 2-3 所示。FL
4、AC3D 3.00Itasca Consulting Group, Inc.Minneapolis, MN USAStep 17407 Model Perspective15:44:44 Mon Dec 05 2011Center: X: 3.000e+001Y: 3.000e+001 Z: 1.000e+001Rotation: X: 0.000Y: 0.000 Z: 50.000Dist: 4.171e+002 Mag.: 0.64Ang.: 22.500Surface Magfac = 0.000e+000Axes LinestyleXYZ图 2-3 计算模型 I 示意图x、y 、 z
5、各方向的长度分别为 60 m、60 m 和 140 m。模型的表面为自由边界，底部为固定边界条件，四周为法向位移约束边界条件。模型共有 10912 个节点，9360 个实体单元。对于单侧壁导坑法，由于整个隧道模型左右不对称，可取隧道模型的一半计算，隧道的计算模型 II 如图 2-4 所示。x、y 、 z 各方向的长度分别为 120 m、60 m 和 140 m。模型的表面为自由边界，底部为固定边界条件，四周为法向位移约束边界条件。模型共有 25637 个节点，23220 个实体单元。FLAC3D 3.00Itasca Consulting Group, Inc.Minneapolis, MN
6、USAStep 20000 Model Perspective15:42:52 Mon Dec 05 2011Center: X: 1.332e+001Y: 3.882e+001 Z: 1.000e+001Rotation: X: 0.000Y: 0.000 Z: 40.000Dist: 4.371e+002 Mag.: 0.64Ang.: 22.500Surface Magfac = 0.000e+000Axes LinestyleXYZ图 2-4 计算模型 II 示意图FLAC3D 提供了丰富的单元库供用户选择，这样可以对各种材料进行模拟。计算中 V 级围岩均采用摩尔 库伦材料，初期支护体
7、系中的喷层混凝土视为线弹性体，二次衬砌混凝土及仰拱均采用 FLAC3D 软件内置的壳单元（shell）进行模拟，系统锚杆和锁脚锚杆均分别采用锚索单元（cable）和桩单元（pile ）进行模拟，超前导管采用桩单元（pile）进行模拟。各材料参数见表 2-1。1.3 数值计算说明需要说明的是，由于研究该段埋深很大，最大处达到 540 m。此处假设隧道埋深为 500 m，假设地应力场仅考虑自重应力场，则按照等效方法在模型顶部竖直方向应施加的垂直向应力为 = 26509.8420 =10.91MPa。地应力平z衡后，模型底部垂直向应力 = 26509.8560 =14.54MPa，隧道中心垂直向z应
8、力 = 26509.8500 =12.99MPa 现按照上述方法进行平衡，得到的竖向应z力云图如图 2-5 及图 2-6 所示。FLAC3D 3.00Itasca Consulting Group, Inc.Mineapolis, MN USAStep 17407 Model Perspective18:51:01 Mon Dec 05 201Center: X: 3.00e+01Y: 3.00e+01 Z: 1.00e+01Rotation: X: 0.00Y: 0.00 Z: 40.00Dist: 4.171e+02 Mag.: 1Ang.: 2.50Contour of SZ Magfa
9、c = 0.00e+00Gradient Calculation-1.3537e+07 to -1.350e+07-1.350e+07 to -1.300e+07-1.300e+07 to -1.250e+07-1.250e+07 to -1.200e+07-1.200e+07 to -1.150e+07-1.150e+07 to -1.100e+07-1.100e+07 to -1.050e+07-1.050e+07 to -1.00e+07-1.00e+07 to -1.00e+07Interval = 5.0e+05FLAC3D 3.00Itasca Consulting Group,
10、Inc.Mineapolis, MN USAStep 17407 Model Perspective18:51:01 Mon Dec 05 201Center: X: 3.00e+01Y: 3.00e+01 Z: 1.00e+01Rotation: X: 0.00Y: 0.00 Z: 40.00Dist: 4.171e+02 Mag.: 1Ang.: 2.50Contour of SZ Magfac = 0.00e+00Gradient Calculation-1.3537e+07 to -1.350e+07-1.350e+07 to -1.300e+07-1.300e+07 to -1.25
11、0e+07-1.250e+07 to -1.200e+07-1.200e+07 to -1.150e+07-1.150e+07 to -1.100e+07-1.100e+07 to -1.050e+07-1.050e+07 to -1.00e+07-1.00e+07 to -1.00e+07Interval = 5.0e+05图 2-5 模型平衡地应力后垂直向应力分布云图（单位：Pa）FLAC3D 3.00Itasca Consulting Group, Inc.Mineapolis, MN USAStep 200 Model Perspective18:57:4 Mon Dec 05 201
12、Center: X: 1.32e+01Y: 3.82e+01 Z: 1.00e+01Rotation: X: 0.00Y: 0.00 Z: 40.00Dist: 4.371e+02 Mag.: 1Ang.: 2.50Contour of SZ Magfac = 0.00e+00Gradient Calculation-1.403e+07 to -1.400e+07-1.400e+07 to -1.350e+07-1.350e+07 to -1.300e+07-1.300e+07 to -1.250e+07-1.250e+07 to -1.200e+07-1.200e+07 to -1.150e
13、+07-1.150e+07 to -1.100e+07-1.100e+07 to -1.050e+07-1.050e+07 to -1.092e+07Interval = 5.0e+05FLAC3D 3.00Itasca Consulting Group, Inc.Mineapolis, MN USAStep 200 Model Perspective18:57:4 Mon Dec 05 201Center: X: 1.32e+01Y: 3.82e+01 Z: 1.00e+01Rotation: X: 0.00Y: 0.0 Z: 4. 0Dist: 4.371e+02 Mag.: 1Ang.:
14、 2.50Contour of SZ Magfac = 0.00e+0Gradient Calculation-1.403e+07 to -1.400e+07-1.400e+07 to -1.350e+07-1.350e+07 to -1.300e+07-1.300e+07 to -1.250e+07-1.250e+07 to -1.200e+07-1.200e+07 to -1.150e+07-1.150e+07 to -1.100e+07-1.100e+07 to -1.050e+07-1.050e+07 to -1.092e+07Interval = 5.0e+05图 2-6 模型平衡地
15、应力后垂直向应力分布云图（单位：Pa）隧道开挖进尺为 2 m。纵向每个网格长为 2 m。纵向有 30 个网格，即不论哪种工法，开挖到 30 步时，隧道模型即贯穿。表 2-1 材料参数汇总表材料 密 度（kg/m 3） 弹性模量（GPa） 泊松比 黏聚力（MPa） 内摩擦角 ()V 级围岩 2650 5 0.4 0.085 38初衬 2438.7 30.11 0.2 二次衬砌 2500 30 0.2 系统锚杆 2479 92.1 锁脚锚杆 2549 68.6 超前导管 2549 68.6 对于初期支护中钢拱架和钢筋网的支护作用采用等效方法计算，即将钢拱架和钢筋网的弹性模量折算给混凝土，其计算方法
16、为： g0cAE（2-1 ）式中：E折算后的混凝土弹性模量（Pa） ；E0原混凝土的弹性模量（Pa ） ；Ag初期支护钢拱架的截面积（m 2） ；Eg初期支护钢拱架的弹性模量（Pa） ；Ac混凝土的截面积（m 2） 。对于系统锚杆及锁脚锚杆的参数选取，亦同样采用折算方法计算获得，其计算方法如下：21AE（2-2）E折算后的锚杆弹性模量（Pa） ；E1水泥砂浆弹性模量（Pa ） ；E2钢管弹性模量（Pa ） ；A锚杆截面积（m 2） ；A1水泥砂浆截面积（m 2） ；A2钢材截面积（m 2） ；对于隧道，实际开挖中，由于卸载打破了原来的平衡，围岩内的各质点在地应力的作用下，将沿最短距离向消除阻力
17、的自由表面方向移动，引起围岩内应力的重新分布，直至达到新的平衡，形成所谓的“二次应力场” 。开挖卸载之前，沿开挖边界上的各点都处于一定的初始应力状态，开挖使这些边界的应力解除，从而引起围岩变形和应力场的变化。从应力路径上看，隧道的开挖过程中应力场的变化是卸载的产物。因此，如何模拟出这种卸载作用显得尤其重要。为更精确的计算出隧道围岩变形及受力变化，Duncan 和 Dunlop（1969） i最早提出了“反转应力释放法” ，之后，国内外学者对其展开了深入研究，孙钧（1988） ii曾建议提出“ 反转应力释放法” ，Mana iii曾提出由挖去单元所产生的释放荷载节点力向量提取等效释放荷载向量的方
18、法，这些都对隧道开挖模拟方面，起到了指导作用。在此处开挖岩体应力释放过程中，采用 Mana 和 Clough57提出的方法，即在每一步开挖之后，开挖体边界的节点将失去开挖体与之相连单元贡献的力，故在这些节点处，不平衡力将不再为 0，其为开挖掉的与之相连单元贡献的力的矢量和，指向开挖面。而其他位置的节点，由于未发生应力增量，故其不平衡力仍然为 0。此时再通过计算得出新的位移场和应力场。再进行下一步开挖计算。在此处的三维模型模拟当中，由于考虑到隧道开挖的时空效应，同时考虑到其与支护的共同作用。对应每一次开挖步，应力并不是一次释放完全的，当刚刚开挖到隧道某个位置时，释放这一部分的 80%的不平衡力，
19、当初衬在这个位置成环时，再释放这一部分的 10%的不平衡力，当仰拱或二衬到达这个位置时，最后释放其 10%的不平衡力。此处运用不同工法进行施工时，具体的模拟过程如下：（1）在模型顶部施加等效应力，平衡地应力；（2）隧道开始开挖，每开挖一步，即计算 1 时步，记录此时开挖边界节点处的不平衡力，即原来开挖体对围岩的支撑力，反向加载这 10%的支撑力，再计算 3000 时步，隧道最大不平衡力很小，即可认为开挖一步后达到平衡。（3）开挖继续向前推进，当施作初衬时，反向加载此位置初始记录时的开挖界面各节点的 10%的不平衡力；同时前方工作面继续开挖，按照（2）释放应力。计算 3000 时步后可认为这一开挖步结束。（4）仰拱或二次衬砌在隧道某一位置施作完后，反向加载此位置初始记录时的开挖界面对应节点的 10%的不平衡力，二衬施作时应力释放与仰拱类似，此时开挖向前推进，掌子面处仍按照（2）释放应力。计算 3000 时步之后可认为这一开挖步结束。
