下穿高速公路浅埋大跨度黄土隧道施工措施研究

1、第32卷第9期 岩 土 力 学 V ol.32 No. 9 2011年9月 Rock and Soil Mechanics Sep. 2011收稿日期:2010-05-31基金项目:铁道部重大科技研发计划资助项目(No. 2005K001-D (G-2。文章编号:1000-7598 (2011 09-2803-07下穿高速公路浅埋大跨度黄土隧道施工措施研究李 健1,谭忠盛1,喻 渝2,倪鲁肃1(1. 北京交通大学 土建学院,北京 100044;2. 铁道第二勘测设计院,成都 610031摘 要:依托郑西客运专线阌乡隧道,开展浅埋大跨黄土隧道优化施工研究。对隧道进入下穿高速公路前的地表沉降,拱顶

2、下沉及初支钢架应力进行现场监测,结果表明,土体沉降及初支应力不能得到有效控制,现有的施工方案无法满足下穿段施工安全要求,需要及时对其调整;经综合分析,提出了减小开挖面积、侧导多台阶开挖、预留核心土、快速封闭及加强初期支护等一系列优化措施,并采用数值计算软件FLAC 3D 验证其安全性和可靠性。通过现场实施及监测结果的反馈,证明优化措施效果明显,保障了阌乡隧道顺利下穿连霍高速公路。研究成果为类似施工中采取相应措施控制土体沉降、保护周边建筑设施的安全使用提供了依据。关 键 词:浅埋暗挖;下穿施工;黄土隧道;现场监测;施工优化 中图分类号:U 456.3 文献标识码:AResearch on con

3、struction procedure for shallow large-spantunnel undercrossing highwayLI Jian 1,TAN Zhong-sheng 1,YU Yu 2,NI Lu-su 1(1. School of Civil Engineering and Architecture, Beijing Jiaotong University, Beijing 100044, China;2. The Second Survey and Design Institute of China Railway, Chengdu 610031, China A

4、bstract : Wenxiang tunnel for Zheng-Xi Special Passenger Railway has been done a series of research on optimizing the construction method for shallow large-span loess tunnel. According to the field monitoring of ground settlement, vault settlement and primary support stress of the tunnel section bef

5、ore tunnel undercrossing the Lian-Huo expressway. The result indicates that soil settlement and primary support stress can not be controlled very well; and the current construction scheme needs to be reformed since it can not satisfy construction safety requirement. Through synthetical analysis, som

6、e optimized methods such as reducing excavation area, increasing excavation steps of pilot tunnel, reserving core soil, quick closure, reinforcing primary support, etc. have been used in construction; and their security and reliability are verified by numerical analysis software fast Lagrangian anal

7、ysis of continua in 3 dimensions (FLAC 3D . From feedback of field implement and monitoring result, the effect of optimized methods is significant; and these methods ensure the undercrossing project implements successfully. These research results have high referential value for analogous project to

8、control soil settlement and protect the safety of neighbouring building facilities.Key words: shallow burying mining; undercrossing construction; loess tunnel; field monitoring; construction optimization1 引 言大跨隧道施工方法,国内外从20世纪80年代后主要采用台阶分部、中壁法以及双侧壁法及环形开挖留核心土法开挖,而对于大跨黄土隧道施工,国内成功采用预留核心土短台阶法解决开挖面积近140 m

9、 2大跨老黄土隧道施工(黄延高速公路道南隧道紧急停车带,并在浅埋新黄土双线铁路隧道中采用CRD 法成功下穿既有铁路(宝兰铁路新曲儿岔隧道1-2。阌乡隧道集大跨(跨度15.6 m 、浅埋(埋深10 m 、软岩(Q 3砂质黄土地层、下穿建构筑物(连霍高速公路等各种情况于一身,在国内外尚无先例,因而,工程施工难度极大,工程施工安全的控制极其重要。按照工程要求实现安全和快速施工,需要对目前修建大跨隧道的主要施工方法进行对比分析和综合评价,在此基础上提出合理的施工方法。目前,国内外学者虽然对隧道施工引起的地层岩 土 力 学 2011年 变形和支护内力的计算及其安全性评价进行了大量的研究3-9,但针对黄土

10、隧道施工不同施工工序及间距对地层变形影响,以及采用现场监测数据结合数值模拟的相关研究,还未见相关文献报道。本文将结合现场试验和数值分析两种手段,对试验段地表沉降、拱顶下沉及初支内力等试验数据进行分析,得出各个施工步开挖沉降所占比重。通过FLAC 3D 软件在已有监测数据的基础上进行验证计算,推导出符合理论计算的土体参数,并应用于优化施工方案中,验证其控制地层变形效果。2 工程概况郑西客运专线阌乡隧道位于郑西客运专线西段的灵宝市,隧道全长为770 m ,为双线大跨浅埋黄土隧道。隧道在以15°的小角度下穿连霍高速公路,下穿段长达270 m ,并分为两部分:正式下穿段为160 m ,从公路

11、正下方穿过,此段的埋深仅有10 m 左右;两侧公路路堑边坡下穿段,共计110 m ,埋深为1024 m 。试验段在公路边坡一侧,长40 m ,隧道与高速公路具体位置关系见图1。隧道所处的地层为Q 3砂质黄土地层,土质松散,土体稳定性极差,而隧道开挖断面特大(阌乡隧道的最大开挖跨度为15.6 m 、高为13.6 m ,最大开挖面积达175 m 2,施工难度极大。图1 隧道与高速公路位置关系Fig.1 Position relationship betweentunnel and highway与岩体相比,黄土土体强度低、变形大、自承能力小、工程性质差,受水的影响十分强烈,一旦被水浸泡达到饱和状态

12、,其强度会明显降低,工程性质发生很大变化。在黄土隧道中,因其土质比较松散,开挖后易形成分层坍塌,层厚一般为20 60 cm 之间,故需合适的施工工法。在试验段将隧道的施工工法进行转换,改CRD 法为双侧壁导坑法,具体施工流程简化如下:(1先开挖右导坑上台阶,开挖高度为7 m ,每次进尺0.8 m ,开挖后立即进行支护。仰拱开挖后立即支护闭合,并及时回填。(2开挖左导坑,开挖方法与右侧相同,开挖面与右侧维持相差10 m 左右的间距。(3开挖中间部分土体,采用台阶法开挖,每次进尺0.8 m ,台阶长度为4 m 左右,开挖约2 m 后,下台阶与仰拱一次开挖,并及时闭合仰拱及回填。图2为试验段双侧壁导

13、坑法施工工序图。试验段具体的支护参数如下:(1超前支护:采用壁厚为8 mm 的159 mm 单层大管棚,管棚的间距为40 cm ,在拱顶110°范围施作,长度为70 m 。掌子面采用25 mm 纤维锚杆,长为12 m ,间距为1.2 m ,梅花形布置。(2初期支护:采用I25a 型钢钢架、挂网喷混凝土,钢架间距为0.8 m 。在边墙采用22 mm 药包锚杆加强支护,药包锚杆长为4 m ,间距为1.0 m ,梅花形布置。图2 隧道施工步序图(单位:mm Fig.2 The excavation procedures of tunnel(unit :mm3 测点布置3.1 地表与拱顶测点

14、布置地表与拱顶沉降布置,在试验段共设置了8个地表沉降观测断面,每个断面间隔为5 m ,采用精密水准仪测量,测量频率1次/d ,具体布置图见图3、4。图3 地表沉降监测点布置图Fig.3 Layout of monitoring points for groundsurface settlementsI18横撑(必要时设置I18横撑(必要时设置I18临时钢架 I18临时钢架I18横撑I18横撑 22水平锚杆 22水平锚杆(L =3 500 (L =3 500 约3000正式下穿前的40 m 作为试验段阌乡隧道郑州西安15°2804第9期 李 健等:下穿高速公路浅埋大跨度黄土隧道施工措施

15、研究(a 测点布置(DK298+775、+785(b 测点布置(DK298+760、+765、+770、+790、+795、+800图4 地表横断面测点布置图Fig.4 Layout of lateral ground surface monitoring points3.2 初支内力测点布置在DK289+785断面的钢架共布设了12处测点,每个测点的钢架内外侧各布设1个应变计测试钢架应力,共24个应变计。在断面内侧布设4个混凝土表面应变计,测试喷混凝土的应力,具体布置如图5所示。 图5 钢拱架应力测点布置图Fig.5 Layout of stress measuring points ons

16、teel framework4 试验结果及分析4.1 拱顶及地表沉降分析经半年多时间的监测,得到试验段典型断面拱顶下沉及地表沉降曲线如图6所示。由图可以看出:(1隧道拱顶总沉降比地表总沉降小,拱顶沉降在310460 mm 之间,地表沉降在480540 mm 之间,拱顶总沉降为地表总沉降的65%94%。(2左导洞超前右导洞5 m 左右、超前中间土体12 m 左右。左导洞拱顶超前沉降占总沉降11%,左导洞地表超前沉降占总沉降18%;右导洞拱顶超前沉降占总沉降15%,右导洞地表超前沉降占总沉降22%;中间部分拱顶超前沉降占总沉降的72%,中间部分地表超前沉降占总沉降的58%,主要受两侧导洞开挖的影响

17、。12-41-32-23-3 4-2 5-2 6-1沉降量/m m左洞拱顶 中洞拱顶 右洞拱顶日期/月-日 左导封闭 右导封闭(a DK298+780断面拱顶下沉全位移历时曲线010020030040050060011-2912-191-81-282-173-8 3-284-175-7沉降量/m m1#2#3#4#5#6#7#8#日期/月-日左导封闭 右导封闭(b DK298+775断面地表沉降全位移历时曲线图6 隧道施工过程拱顶及地表沉降全位移历时曲线Fig.6 Vault and ground settlements time history curvesduring tunnel con

18、struction progress(3左、右导洞上下台阶法开挖,拱顶沉降的超前影响距离均为4 m 左右,地表沉降的超前影响距离均为8 m 左右;中间土体3台阶法开挖,拱顶沉降的超前影响距离为8 m 左右,地表沉降的超前影响距离为12 m 左右。(4地表沉降与拱顶沉降曲线规律一致,在侧导洞开挖到DK298+780断面之后到中间土体也开挖到该断面时,这段时间内发生较大沉降,占总沉降的65%80%,中间部分土体开挖引起的沉降占15%20%,全断面封闭后的沉降仅占3%。 4.2 初支内力分析钢架应力的时间历时曲线如图7所示,应力分布如图8所示,混凝土应力如图9所示。由图79可以看出:(1初支钢架的拱

19、顶内外侧均受拉,其他部分钢架内外侧均为受压。拱顶应力大小在20 MPa 左右;右导坑拱顶压应力最大,达-182.1 MPa ;边墙应力大小达-81 MPa ;仰拱中间受力较小、两侧受力较大达98 MPa 。(2先开挖左右导洞的钢架受中间部开挖影响明显,中间拱部土体的开挖能使两侧导洞钢架的应力改变1040 MPa 。(3全断面封闭能使钢架应力稳定,在拆除临时支撑的双侧壁时,断面钢架应力会有明显变大,直到二衬施工完成后才最终稳定(9月底施作二衬。隧道中线9点8点 7点 6点5点4点 3点2点地面1点5.0 m 3.5 m 3.0 m 2.5 m 3.0 m 3.5 m 5.0 m隧道中线地面5点4

20、点3点2点1点5.0 m 5.5 m 5.5 m 5.0 m2805岩 土 力 学 2011年 -140-100-60 -20 204-2 5-2 6-1 7-1 7-318-30 9-2910-2911-28日期/月-日 应力/M P aWE1WE2WE4WE5WE6WE7WE8WE9WE10拱部开挖 全段面封闭临时支撑拆除左导封闭 右导封闭(a 外环应力-200-160-120 -80 -40 040 4-2 5-2 6-1 7-1 7-318-309-29 10-29 11-28日期/月-日 应力/MP aNE1NE2NE3NE4NE5NE6NE7NE8NE9NE10拱部开挖全段面封闭

21、临时支撑拆除 右导封闭左导封闭(b 内环应力图7 DK298+785断面钢架应力时态曲线图 Fig.7 Time history curves of stresses in steelframe at section DK298+785(a 钢架外侧应力 (b 钢架内侧应力图8 DK298+785断面钢架应力分布图(单位:MPa Fig.8 Stress distributions of steel frame arch atsection DK298+785 (unit: MPa-8-6-4-20 2 4 6 8 5-76-67-68-59-410-411-3 12-3日期/月-日混凝土应力

22、/M P aNC1NC2NC3NC4图9 DK298+785断面混凝土应力时态曲线Fig.9 Time history curves of stress in concreteat section DK298+785(4拱顶混凝土受拉,拉应力大小在1.8 MPa 左右;拱腰处(接近侧导洞拱顶混凝土受压,压应力为-7 MPa 左右;仰拱处的喷混凝土受压,应力在4 MPa 左右。(5由钢架应力和喷混凝土应力的测试结果可知,钢架与喷混凝土的应力分布规律一致,但钢架应力远远大于喷混凝土应力,钢架应力大约是喷混凝土应力的10倍。 4.3 试验段施工分析为了得到隧道施工过程拱顶的沉降规律以及拱顶沉降与地表

23、沉降的关系,根据图6中数据计算出各开挖步产生的沉降及比例,见表1。表1 各开挖步产生的沉降及比例(单位:mm Table 1 Amount and proportion of settlement in everyexcavation step (unit: mm左导洞隧道中线右导洞施工步观测点沉降值占总 沉降 /%沉降值 占总 沉降 /%沉降值占总沉降/%拱顶沉降51.511 0.87 0 5.5 2 左导坑开挖到DK298+780地表沉降85.518 73.7 14 57.912 拱顶沉降115.226 24.7 5 45.815 右导坑开挖到DK298+780地表沉降141.829 12

24、8.5 23 105.822 拱顶沉降369.882 361.4 77 260.884 中间土开挖到DK298+780地表沉降393.482 443.7 81 392.482 拱顶沉降438.397 451.9 97 300.197 全断面封闭后地表沉降465.997 527.4 97 463.197 拱顶总沉降 451.394 467.5 86 309.565 地表总沉降481.2544.1 478.0由表1可知,试验段地表沉降并没有得到很好的控制,按照试验段施工方案很难达到下穿段要求的50 mm 的控制标准10,从而通过分析数据提出以下优化措施:(1由于侧导洞开挖产生的沉降所占的比例最大,

25、因此,为了控制沉降,两侧导洞开挖时分为4步(即3个台阶和仰拱,减小了每步的开挖面积,并预留核心土,减小上中台阶长度,上台阶长度为2.43 m ,中台阶长度为2.43 m ,下台阶长度为5 m 左右,这样能有效地控制掌子面的稳定,减小施工对地表的影响。(2中间土体开挖基本上与试验段一样,也是分上中下3台阶开挖,改进的是在开挖上台阶时预留核心土,以及中下台阶同时开挖。这样能有效地控制掌子面的稳定,并快速地封闭仰供,减小地表的沉降。(3为了能快速封闭开挖面,减少开挖面暴露的时间,减少地应力的释放,从而减小地表沉降,2806第9期李健等:下穿高速公路浅埋大跨度黄土隧道施工措施研究开挖进尺由原来的0.8

26、 m减小为0.6 m。(4侧导洞下台阶及仰供、中间土体上中下台阶都采用机械开挖,这样有利于快速施工,快速封闭仰供,能更有效地控制地表沉降。(5在阌乡隧道施工中,初期支护的刚度对地层变形有较大影响,一方面由于采用双侧壁导坑法施工,当施作两侧弧形导坑内的钢架时,很难将其控制在同一平面内,使得顶部弧形导坑开挖后,钢架无法连成整体,支护承载能力大大降低;另一方面由于隧道上方的高速公路经常有重载车辆通行,而埋深只有10 m。故采用加密钢架、双层初支的措施加强初期支护,其中第2层初支采用格栅喷混凝土有利于使两层支护紧密黏结,能比模筑混凝土更快地提供刚度;同时,由于试验段没有进行初支背后填充注浆而使地表沉降

27、很大,因此,在下穿段必须增加初支背后的填充注浆。5 数值分析为了论证优化施工方案对地层变形的控制效果,进一步深入研究优化施工工序下各施工步的开挖对土体变形计结构内力的影响,利用三维有限差分软件FLAC3D,根据实际开挖步骤,对隧道施工全过程采用Mohr-Coulomb准则进行动态模拟11-12。5.1 模型的建立在有限元计算中,边界约束条件对计算结果影响较大,根据潘昌实等13-14的研究,本次计算范围选取为:试验段上部取至地表24 m,下穿段取至地表10 m;左右部及下部边界取至离隧道外缘周边2倍洞径,约30 m,模型长为30 m;左右为水平约束边界,下部为垂直约束边界,地表为自由边界;试验段

28、模型每施工步台阶长度为4 m,每次进尺为0.8 m,下穿段模型每施工步台阶长度为2.4 m,每次进尺0.6 m,研究断面取在距前边界10 m处。计算中用8节点六面体实体单元模拟围岩、初期支护及临时支护,管棚的作用采用等效方法予以考虑,将管棚的弹性模量折算给地层15。下穿段计算模型总单元数为68 160个,总节点数为71 126个。5.2 围岩与支护参数计算模型中,对锚杆加固的土体围岩,按经验采取提高其c,值来加以仿真模拟16-17,围岩及支护材料的物理力学参数取值如表2所示。在试验段模型计算时,根据现场监测数据对围岩物性指标弹性模量E进行调整,当模型计算结果中地层变形及初支应力与监测数据相同时

29、,将该物性指标用于下穿段优化施工的模型计算中,以保证数值计算的准确性。表2 围岩物性指标及支护结构参数Table 2 Mechanical parameters of ground and support 围岩及结构容重/(kN/m3弹性模量E/GPa泊松比基床系数K/(MPa/m凝聚力c/MPa内摩擦角/(° C30钢筋混凝土C20钢筋混凝土5.3 结果分析(1拱顶沉降分析分析计算结果,得到每步开挖造成的沉降量及占总沉降量的比例,如图10和表3所示。由图9及表3可以发现:随着开挖推进,在20 m距离以前,沉降逐渐增长;之后拱顶沉降趋于收敛;初期支护拱顶沉降不大,最终拱顶沉降值约为4

30、7 mm;由于开挖跨度大,开挖跨度D =16 m,故收敛距离约在1.5D左右; 1020304050600102030 40 50推进距离/m拱顶沉降/mm图10 初期支护拱顶沉降变化曲线图Fig.10 Vault settlement curve of primary support表3 各工序开挖阶段造成的沉降及其比例(单位:mmTable 3 Amount and proportion of settlement caused by every excavation procedure (unit: mm 计算值工序1 工序2 工序3 工序4 前4步工序5 工序6 工序7 工序8 封闭后

31、阶段值/mm 8.79 4.39 1.75 2.5217.454.405.637.876.264.80 占百分比/% 18.94 9.45 3.775.4337.599.4812.1316.9613.4910.3428072808 岩 土 力 学 2011 年 前 4 步工序造成的拱顶沉降为 17 mm 左右， 占总沉降的并不大，说明减小开挖面积、小步距开 后 4 步工序中要注意控制整体沉降，可以采用对钢 架拱脚进行加固的方式，加固方法为增设锁脚锚管 （4 根） 、大拱脚、牛腿和槽钢垫板。 （2）初支内力分析 分析计算结果，得到最终初期支护应力及初期 支护结构安全性，如表 4 所示。 表 4

32、初期支护最终应力（单位：MPa） Table 4 Final stresses of primary support (unit: MPa 部位 拱顶 左拱脚 右拱脚 左墙脚 右墙脚 仰拱 最大主应力 0.4 11.9 6.0 9.2 7.4 1.6 安全性 安全 安全 安全 安全 安全 安全 最小主应力 -1.7 -13.4 -3.8 -7.3 -8.1 -6.5 安全性 安全 安全 安全 安全 安全 安全 累积沉降量/mm 10-31 0 10 20 30 40 50 11-30 日期/月-日 12-30 1-29 2-28 3-30 875-3 880-2 885-4 900-2 915

33、-2 920-2 890-2 895-2 挖、早封闭的优化施工方式取得了预期的效果。在 图 12 高速公路隔离带测点沉降历时曲线 Fig.12 Settlement time history curves of spots on sides of expressway 7 结 论 （1）为了稳定浅埋大跨黄土隧道掌子面，防止 塌方，必须减小每步的开挖面积，将侧导洞分上、 中、下台阶及仰拱开挖，中间土体分上、中、下台 阶开挖，且上台阶预留核心土；同时，采用双层初 期支护，提高初期支护刚度以控制地层变形；采用 人工和机械相结合的开挖方法，实现快速施工、快 速封闭仰供，能更有效地控制地表沉降。阌乡隧道

34、 下穿段采用优化设计方案后，地表沉降基本被控制 在 50 mm 以内。 （2）加快开挖面封闭、导洞断面封闭、隧道全 断面封闭，能非常有效地控制地表沉降。从大量黄 土隧道的施工监测数据可以看出，断面封闭后的沉 降占总沉降 10%以下，阌乡隧道试验段仰供封闭后 的地表沉降占总沉降不到 5%，而且断面封闭后很 快稳定。因此，断面的及时封闭对控制地表沉降十 分重要。 参 考 文 献 1 谭忠盛, 喻渝, 王明年, 等. 大断面浅埋黄土隧道锚杆 作用效果的试验研究J. 岩土力学, 2008, 29(2: 491 496. TAN Zhong-sheng, YU Yu, WANG Ming-nian, e

35、t al. Experimental study on bolt effect on large section shallow 由表可以发现： 初期支护在施工过程中出现了应力集中，拉 压应力集中均发生在供脚、墙脚。 施工过程中，最大主应力为 11.9 MPa，最小 主应力为-13.4 MPa，小于钢筋的屈服强度，远远小 于试验段钢架实测内力，符合安全性要求。 6 工程应用 采用施工优化方案后，下穿段隧道成功穿越连 霍高速公路，并且地表沉降得到了有效控制。 由图 11、 可以得到： 12 截止阌乡隧道顺利贯通， 测点的地表沉降最大值在为 56 mm（因为下穿段施 工施做双层管棚，所以停工 1

36、个月，处于正上方 DK298+885 处地表沉降累积值较大） ，沉降最大的 测点距隧道全断面封闭距离为 20 m，并且收敛较 快，在开挖后 2D 范围距离处收敛稳定，不会再发 生较大的变化，满足行车要求。 日期/月-日 1-29 3-2 11-26 0 累计沉降量/mm 16 32 48 64 80 12-28 4-3 左导 865-2 左导 870-2 左导 875-2 左导 880-1 左导 885-3 中导 890-1 中导 895-1 中导 900-1 中导 905-2 右导 905-1 右导 910-1 右导 915-1 右导 920-1 右导 925-1 depth loess t

37、unnelsJ. Rock and Soil Mechanics, 2008, 29(2: 491496. 2 赵学动, 倪玉山. 黄土隧道围岩变形的长期观测及其 稳定性的分析J. 现代隧道技术, 2002, (增刊: 371 374. ZHAO Xue-dong, NI Yu-shan. Long-term observation of deformation and analysis of stability in loess tunnelJ. Modern Tunnelling Technology, 2002, (Supp.: 371 374. 图 11 高速公路路面测点沉降历时曲线

38、Fig.11 Settlement time history curves of spots on expressway 第9期 李 健等：下穿高速公路浅埋大跨度黄土隧道施工措施研究 2809 3 刘惠敏, 赵月, 吴波. CRD 法不同施工工序的比较与分 析J. 隧道建设, 2005, 25(5: 610. LIU Hui-min, ZHAO Yue, WU Bo. Comparison and analysis of different working procedures of CRD methodJ. Tunnel Construction, 2005, 25(5: 610. Univ

39、ersity, 2008, 32(4: 137140. 11 张志强, 何川. 偏压连拱隧道优化施工的研究J. 岩土 力学, 2007, 28(4: 723727. ZHANG Zhi-qiang, HE Chuan. Research on optimized construction procedure for double arched tunnel in unsymmetrical pressureJ. Rock and Soil Mechanics, 2007, 28(4: 723727. 12 杨小礼, 眭志荣. 浅埋小净距偏压隧道施工工序的数 值分析J. 中南大学学报（自然科学版

40、）, 2007, 38(4: 764770. YANG Xiao-li, SUI Zhi-rong. Numerical simulation of construction sequence for shallow embedded bias tunnels with small clear distanceJ. J. Cent. South Univ. (Science and Technology, 2007, 38(4: 764770. 13 徐干成, 白洪才, 郑颖人, 等. 地下工程支护结构M. 北京: 中国水利水电出版社, 2002. 14 潘昌实. 隧道力学数值方法M. 北京:

41、 中国铁道出版 社, 1995. 15 李术才, 朱维申. 弹塑性大位移有限元方法在软岩隧 道变形预估系统研究中的应用J. 岩石力学与工程学 报, 2002, 21(4: 466470. LI Shu-cai, ZHU Wei-shen. Application of elasto-plastic large displacement finite element method to the study of deformation prediction of soft rock tunnelJ. Chinese Journal of Rock Mechanics and Engineerin

42、g, 2002, 21(4: 466470. 16 李强, 王明年, 李玉文. 双跨连拱隧道 2 种中墙的空间 力学效应分析J. 岩土力学, 2006, 27(4: 667672. LI Qiang, WANG Ming-nian, LI Yu-wen. Threedimensional mechanical effect analysis of two types of midboard for double span-arch-connecting tunnelJ. Rock and Soil Mechanics, 2006, 27(4: 667672. 17 朱合华, 丁文其. 地下结构施工过程的动态仿真模拟 分析J. 岩石力学与工程学报, 1999, 18(5: 558562. ZHU He-hua, DING Wen-qi. The