FLA数值模拟分析.docVIP

1工程概况 研宄此段为中条山隧道K9+450?K10+560段，此处隧道最大埋深约540 m, 主要由太古界涑水群表壳岩组合解州片麻岩（Hgn）地层组成，构成中条山隧道 分水岭北侧的主体；此段组成隧道的围岩岩性主要为变粒岩、花岗片麻岩等。该 套地层岩性复杂，组合无规律。岩层产状整体倾向南东，倾角一般在50°?70° 间变化。在AK9+900?AK10+000段为区域性断层影响段，此断层为破碎岩石组 成，将为基岩裂隙水下渗提供通道，隧道开挖必将引起涌水，同时此段围岩稍差, 施工时易引起坍塌。此段同时也是中条山北侧泉水主要涌出段，水文情况复杂。 总体评价，木段工程地质条件差。 在此处，具体运用FLAC3D进行模拟的区段均取洞身YK10+100?YK10+180段。 该区段为V级围岩区域，埋深为505?512 m,为断层，附近太古界涑水群花岗 片麻岩、黑云斜讼片麻岩、岩石破碎。隧道断面为SVc型，如图2-2所示。 图2-2 SVc型隧道断面图 隧道衬砌按新奥法原理设计，采用SVc型复合式衬砌，该衬砌适用于隧道洞 身V级断层影响带及软弱破碎围岩段的初期支护及衬砌，超前支护各环釆用42 X4 mm注漿小导管超前预加固围岩，松4.5 m，环向间距35 cm,搭接长度1.3 m， 斜插角10°?15°,每环37根；初衬以喷、锚、网为主耍支护手段：钢拱架为120a 型钢，纵向间距75 cm，每榀钢拱架之间采用*22钢筋连接，环形间距1.0 m; 锚杆采用D25中空注浆锚杆，长3.5m,间距75 cm （纵）X 100 cm （环），与钢 拱架交错布置；喷C25早强混凝土 26 cm。二次衬砌和仰拱均为C30钢筋混凝土 结构，厚50 cm。 1.2数值计算模型 根据中条山隧道工程的实际状况，为提高计算速度，在保证计算精度的前提 下，取桩号YK10+100?YK10+160段采用大型有限差分软件FLAC3D进行建模分析。 对于全断面法、预留核心土法、台阶法，由于整个隧道模型左右对称，为减 少计算量，可取隧道模型的一半计算，隧道的计算模型I如图2-3所示。 x、y、z各方向的长度分别为60 m、60 m和140 m。模型的表面为自由边界， 底部为固定边界条件，四周为法向位移约束边界条件。模型共有10912个节点, 9360个实体单元。 对于单侧壁导坑法，由于整个隧道模型左右不对称，可取隧道模型的一半计 算，隧道的计算模型II如图2-4所示。 x、y、z各方向的长度分别为120 m、60 m和140 m。模型的表面为自由边 界，底部为固定边界条件，四周为法向位移约束边界条件。模型井有25637个节 点，23220个实体单元。 图2-4计算模型II示意图 FLAC3D提供了丰富的单元库供用户选择，这样可以对各种材料进行模拟。计 算中V级围岩均采用摩尔一库伦材料,初期支护体系中的喷层混凝土视为线弹性 体，二次衬砌混凝土及仰拱均采用FLAC3D软件内置的壳单元(shell)进行模拟， 系统锚杆和锁脚锚杆均分别采用锚索单元(cable)和桩单元(pile)进行模拟， 超前导管采用桩单元(pile)进行模拟。各材料参数见表2-1。 1.3数值计算说明 需要说明的是，由于研宂该段埋深很大，最大处达到540 m。此处假设隧道 埋深为500 m,假设地应力场仅考虑自重应力场，则按照等效方法在模型顶部竖 直方向应施加的垂直向应力为0； = 2650X9.8X420 =10.91MPa。地应力平衡后， 模型底部垂直向应力f = 2650X9.8X560 =14.54MPa，隧道中心垂直向应力= 2650X9.8X500 =12.99MPa现按照上述方法进行平衡，得到的竖向应力云阁如阁 2-5及阁2-6所示。 Contour of SZZMagfac = 0.000e+000 Gradient CalculationI-1.3537e+007to-1.3500e+00 -1.3500e+007to-1.3000e+00 -1.3000e+007to-1.2500e+00 -1.2500e+007to-1.2000e+00 _l-1.2000e+007to-1.1500e+00 Z-1.1500e+007to-1.1000e+00 -1.1000e+007to-1.0500e+00 Contour of SZZ Magfac = 0.000e+000 Gradient Calculation ?1.0500e+007to-1.0000e+00 ■ -1.0000e+007 to -1.0000e+00 Interval = 5.0e+005 图2-5模型I平衡地应力后垂直向应力分布云图（单位：Pa） Contour of SZZMagfac = 0.000e+000 G