以武汉某排水通道工程为工程背景，运用有限元分析软件 Ｐｌａｘｉｓ３Ｄ ，对双侧壁导坑法隧道的施工过程开展数值模拟，重点分析不同循环开挖进尺、不同工序对隧道开挖过程中的变形响应及地表变形规律，优化双侧壁导坑法隧  道施工方案。计算结果表明：循环进尺为４ ｍ 时计算得到的地表沉降值及拱顶沉降值显著大于循环进尺为１ ｍ 和２ ｍ时的对应值。施工顺序１（Ⅰ→ Ⅱ→Ⅲ→Ⅳ→Ⅴ→Ⅵ）较施工顺序２（Ⅰ→ Ⅲ→Ⅱ→Ⅳ→Ⅴ→Ⅵ ）对地层扰动影响更小。对比计算结果，综合考虑施工安全及施工效率的因素，建议隧道的循环进尺取２ｍ，施工工序采用顺序１（Ⅰ→ Ⅱ→Ⅲ→Ⅳ→Ⅴ→Ⅵ）。该研究可为武汉地区相似地层环境条件下的暗挖隧道工程提供设计参考。

关键词：双侧壁导坑法；施工方案优化；有限元；变形响应

大断面隧道由于其断面大、地 质条件复杂等特点，往往是重点控制工 程。特别是对于大跨大断面隧道复杂地质条件下的开挖，保证其围岩稳定性具有重大的工程实际意义。双侧壁导坑法也称眼镜工法，是大断面隧道分部开挖的工法之一，在工程 实 践 中 积 累 了 不 少 成 功 的 经 验。但是，在复杂 地质条件下的大断面隧道，双 侧壁导坑法的施工参数选择，需要结合工程实际情况具体研究。

现阶段，针对暗 挖法隧道施工优化研究已经取得了一定的进展［１－９］。 朱正国等［１］以 兰渝铁路仓圆隧道为背景，确定了穿越泥石流堆积体隧道的围岩加固范围、合 理施工方法及台阶参数；王仁杰等［４］通过现场测试及数值模拟的手段，对双侧壁导坑法两种不同开挖顺序的隧道变形及受力进行了对比分析；李讯等［７］采用数值模拟和现场监控相结合的方法，研究浅埋超大断面隧道的破坏形态、开 挖工法、施 工参数及支护体系力学特性。

武汉某排水隧道工程双侧壁导坑法施工段穿越复杂地层，同时下穿既有地下结构或建（构）筑物，施工风险大。本文以武汉某排水隧道工程为背景，采用 Ｐｌａｘｉｓ３Ｄ 有限元分析软件，建立 红黏土地层—隧道结构—周边环境的一体化分析模型，开展三维数值模拟计算。分析在不同开挖循环进尺，不同施工工序条件下，暗挖隧道施工过程中的变形演化规律，探讨武汉红黏土下穿双侧壁导坑法隧道施工最佳参数选择。

选取典型地质剖面，采用 Ｐｌａｘｉｓ３Ｄ 有限元软件，建立数值计算模型，开 展数值模拟计算。模型高３６．９３ ｍ，模型地表以下依次分布着２．９ ｍ 深的杂填土、８．４ｍ 深的红黏土和２８．７ｍ 深的硅质岩。隧道穿越硅质岩层。其中，隧道 采用双侧壁导坑法施工。隧道初期支护采用 ＣＦ３０ 钢纤维喷射混凝土；二次衬砌采用 Ｃ３５Ｐ８ 混凝土。隧道开挖后及时施作相应的初期支护，初期支护稳定后 施作二次衬砌。隧道洞身采用复合式衬砌，以锚杆、钢筋 网喷混凝土、钢拱架为初期支护，模筑混凝土和钢 筋混凝土为二次支护。排水通道工程现场如图１ 所示。

充分考虑开挖影响范围和边界效应，三维计算模型 Ｘ 方向总长取５０ ｍ，Ｙ 方向长度为２０ ｍ，竖直方向视土层深度而定。隧道全长２０ ｍ，隧道掘进方向为Ｙ 轴正方向。计算模型的侧边界采用法向约束，底部边界采用全约束（法向、切向约束）。

计算分析中，土体采用三角形１５节点单元进行

模拟，有限元划分网格如图２ 所示。其中，初衬、二衬及临时衬砌均采用板单元模拟，锚杆等效到周围土体参数中，其余均采用实体单元模拟，衬砌与土 体接触面上设置接触单元。红黏土采用小应变硬化本构模型，其余土体均采用摩尔－库仑本构模型；初衬、临时支撑和二衬采用线弹性模型。土层计算参数见表１，衬砌计算参数见表２。

结合上述实际工况，分别考虑１、２、４ ｍ３ 种循环进尺长度，１（Ⅰ →Ⅱ→Ⅲ→Ⅳ→Ⅴ→Ⅵ），２（Ⅰ → Ⅲ→Ⅱ→Ⅳ→Ⅴ→Ⅵ）两种施工顺序及二衬支护及时性，对双侧壁导坑法隧道下穿天然气管道进行三维数值模拟，对比分析隧道施工扰动引起的高压天然气管道及周边环境变形演化规律，从而得到最优的施工方案。

施工过程严格按照设计资料进行数值模拟。全隧道均采用上下台阶法施工，每部开挖后施作相应的初期支护及相应的临时支护，全断面开挖完成后及时施作仰拱，然后及时进行二次衬砌施工。双侧壁导坑法隧道横截面及施工支护情况如图 ３ 所示。

分别选取１、２、４ ｍ３种循环进尺长度进行三维有限元分析，其余施工参数保持不变。对以上３ 种循环进尺分别分析计算结果如下：  

取ｙ＝１０ ｍ 处，即隧道中点处断面，分析不同循环进尺条件下隧道开挖完成后该断面的横向地表沉降规律。图４、图５ 分别为地层竖向变形位移云图及地表沉降曲线。可以看出，当 循环进尺为４ ｍ时，隧道开挖产生的地表沉降值整体最大，最大值为７．６ ｍｍ；当循环进尺为１ ｍ 时，隧道开挖产生的地表沉降值整体最小，沉降最大值为 ４．５ ｍｍ。选取ｘ＝２５ ｍ 处，即隧道轴线位置处断面分析得到不同循环进尺下隧道开挖完成后纵向地表沉降曲 线，如图６所示。同样地，当循环进尺为４ｍ 时地表纵向沉降曲线整体最大。

取ｙ＝１０ ｍ 处，即隧道中点处断面，分析不同循环进尺条件下隧道开挖完成后衬砌变形规律。图７、图８ 分别为隧道结构竖向变形图及水平变形图。表３汇总了该断面在不同循环进尺条件下隧道结构变形最大值。可以看出，不同循环进尺下隧道结构的变形模式基本相同。其中，循环进尺为１ ｍ 时隧道结构变形最小，循环进尺为４ ｍ 时最大。循环开挖进尺对隧道结构竖向变形影响更大。

 图９和图１０分别展示了隧道开挖完成后，不同循环进尺条件下隧道轴线处衬砌拱顶沉降及拱底隆起曲线。可以看出，循环进尺为１ ｍ 和２ ｍ 的隧道拱顶沉降值差别较小，而循环进尺为４ ｍ 的隧道拱顶沉降值整体明显大于循环进尺为１ ｍ 和２ ｍ的两个工况。３种循环进尺条件下的隧道拱底隆起曲线差别较小。

对于开挖进尺而言，一 次开挖进尺的步距越大，产生的变形和沉降也越大；开挖步距越小，越容 易控制开挖引起的变形。然而，若步距过小，将影 响人员和机械在隧道中施工，影 响施工效率。因此，结合上 述计算结果，建 议开挖进尺不宜过大。综上，考虑到施工安全及施工效率的因素，建议隧 道的循环进尺取２ ｍ。

双侧壁工法对围岩产生多次扰动，引起应力场的重分布与数次演变，洞周位移非线性叠加，产生 损伤性塑性区，并向深部延伸，支护变位也逐渐增大，即所谓群洞效应。为避免群洞影响，一般先安 排一侧导坑先行，另一侧导坑间隔一定距离后进。对双侧壁导坑法采用顺序１（Ⅰ → Ⅱ→ Ⅲ→ Ⅳ → Ⅴ→Ⅵ），顺序２（Ⅰ →Ⅲ→Ⅱ→Ⅳ→Ⅴ→Ⅵ）２ 种导坑开挖工序，隧道循环进尺默认２ ｍ，采用 Ｐｌａｘｉｓ３Ｄ 软件进行三维有限元计算，分析计算结果，确定最佳施工顺序。

取ｙ＝１０ ｍ 处，即隧道中点处断面，分析不同循环进尺条件下隧道开挖完成后该断面的横向地表沉降规律。图１１、图１２ 分别为地层竖向变形位移云图及地表沉降曲线。可以看出，当采用第２ 种施工顺序时，隧道开挖产生 的地表沉降值整体最大，最大值为６．０ ｍｍ；当采用施工顺序１ 时，隧道开挖产生的地表沉降值整体最小，沉 降最大值为５．５ｍｍ。选取ｘ＝２５ｍ 处，即隧道轴线位置处断面分析得到不同循环进尺下隧道开挖完成后纵向地表沉降曲线，如图１３ 所示。同样地，当采用第２ 种施工顺序时地表纵向沉降曲线整体最大。

取ｙ＝１０ ｍ 处，即隧道中点处断面，分析不同施工工序条件下隧道开挖完成后衬砌变形规律。图１４、图１５ 分别为隧道结构竖向变形图及水平变形图。表４汇总了该断面在不同施工工序条件下隧道结构变形最大值。可以看出，整体上，施工工序不同对隧道结构拱顶沉降及拱底隆起影响较小。

图１６和图１７ 分别展示了隧道开挖完成后，不同施工工序条件下隧道轴线处衬砌拱顶沉降及拱底隆起曲线。可以看出，不同施工顺序对于隧道拱底隆起几乎没有影响，而施工顺序２ 造成的隧道拱顶沉降比施工顺序１更大。

从控制周边环境变形以及保证初期支护的安全角度出发，顺序１ 比顺序２ 更有利，所以从施工安全的角度考虑建议采取顺序１ 开挖。

本文基于武汉某排水隧道工程双侧壁导坑法隧道，分别考虑１、２、４ ｍ３种循环进尺长度，１（Ⅰ→ Ⅱ→Ⅲ→ Ⅳ → Ⅴ→ Ⅵ）、２（Ⅰ → Ⅲ→ Ⅱ→ Ⅳ → Ⅴ → Ⅵ）两种施工顺序，采用 Ｐｌａｘｉｓ３Ｄ 软件，建立了红黏土地层－隧道结构－周边环境的一体化分析模型，形成了大跨扁平暗挖隧道全过程施工效应数值模拟方法。通过分析计算结果，确定了最优施工方案， 同时也为武汉地区下穿红黏土暗挖隧道工程提供了有效参考。具体结论如下：

１）当循环进尺为４ ｍ 时，隧道开挖产生的地表横向和纵向沉降值整体最大，循环进尺为 １ ｍ 和２ ｍ的地表沉降值差异较小；隧道开挖完成时，不同循环进尺条件下的隧道拱顶沉降值差异较小，但循环进尺为４ ｍ 的施工方案隧道拱顶沉降值明显大于循环进尺为１ ｍ 和２ ｍ 所对应的值。

２）当循环进尺为２ ｍ 时，施工顺序１ 计算得到的地表沉降值小于施工顺序２ 的地表沉降值，拱顶沉降及拱底隆起值差异较小。

３）对比计算结果，综合考虑施工安全及施工效率的因素，建议隧道的循环进尺取２ ｍ，施工顺序采用顺序１（Ⅰ→ Ⅱ→ Ⅲ→ Ⅳ → Ⅴ→ Ⅵ），并在实际条件允许的情况下及时进行二次衬砌及仰拱的施作。