超大断面扁平结构隧道开挖方法数值模拟及优化研究.docxVIP

? ? 超大断面扁平结构隧道开挖方法数值模拟及优化研究 ? ? 李 盼 (湖南省交通科学研究院有限公司, 湖南 长沙 410015) 当隧道的开挖断面面积超过140 m2、隧道的扁平率(高宽比)在含仰拱时0.65或在不含仰拱时0.55且隧道断面的形状近似为椭圆形扁平状时，一般称为超大断面扁平结构隧道。随着我国交通流量和对交通条件需求的日益增长，超大断面扁平结构公路隧道的建设呈现出逐渐增多趋势。超大断面扁平结构公路隧道能有效缓解高密度的交通流量，并提高断面利用率，逐渐成为未来隧道结构形式的主流[1]。然而，超大断面扁平结构隧道的受力特性和围岩应力分布情况较为复杂，隧道在开挖过程中的应力重分布趋于不利，这导致隧道底脚处的应力集中过大，支护结构的承载力较小[2]，从而使得隧道拱顶的稳定性降低，这些问题的产生必然会对超大断面扁平结构隧道的施工带来更为严峻的考验。在现阶段研究中，黄朱林[3]通过二维弹塑性有限元数值模拟就扁平率对大跨隧道围岩及支护的稳定性影响进行了分析，得出了扁平率的大小对隧道安全系数的影响机理；范君黎等[4]利用扩展有限元法(XFEM)对不同扁平率隧道的稳定性进行了研究，得出了隧道扁平率的大小与隧道强度、损伤破坏以及裂纹产生的相互关系；熊根贵[5]应用ABAQUS对超大断面隧道施工开挖方法和分步开挖过程进行了有限元数值分析，得出了隧道的最优开挖工序；瞿东明等[6]应用数值有限元法，对比分析了不同隧道开挖方法下隧道受力以及发生变形大小的差异性，得出超大断面隧道的安全施工方法。隧道的开挖是一个复杂的系统工程，因此，有必要针对超大断面扁平结构隧道的开挖施工问题进行数值模拟分析，得出隧道在开挖过程中的受力特性及分布情况，从而提出一种最优的超大断面扁平结构隧道的开挖方法，为隧道开挖施工提供一定的理论依据。 本文结合某超大断面扁平结构公路隧道工程，采用MIDAS/GTS有限元分析软件，建立该隧道的三维开挖施工模型，并详细分析该隧道在3种不同开挖方法下围岩位移场和应力场的变化，揭示该超大断面扁平结构隧道的力学特性及位移变化规律，最后通过对比分析，得到该超大断面扁平结构隧道的最优开挖方法，进而保障隧道的施工安全。 1 工程概况 该公路隧道为温州市某隧道工程，隧道洞长1360m，开挖高度13.8m，宽度16.5m，洞净高约6.9m，开挖断面面积为172m2，扁平率为0.5(不含仰拱)，隧道为上下双层四心卧式椭圆结构。根据地质测绘、钻探、槽探资料，该隧道口自地面以下以残坡积物和全风化基岩为主，风化物为粉质黏土，风化土呈可塑状；遇水后，强度明显降低，属中软土地基土层，隧道围岩等级为Ⅲ～Ⅴ级。 2 不同开挖方法及数值模拟分析 2.1 隧道开挖方案 随着对隧道工程领域研究不断深入，我国逐步发展形成了一套较为完善的隧道开挖方法理论体系[7]。其中针对超大断面公路隧道常用的开挖方法主要有台阶法、CD工法、双侧壁导坑法[8-10]等。台阶法工期短、施工方便，但容易产生大变形；双侧壁导坑法对围岩扰动较小，但工序复杂、成本高。因此，根据该隧道工程的实际特点，提出了双侧壁导坑法结合台阶法的施工方案，并应用MIDAS/GTS软件对三台阶七步法、双侧壁导坑法和双侧壁导坑法结合台阶法进行对比分析，以寻求最优的开挖施工方案。3种开挖方案的开挖工序如图1所示。 图1 隧道不同开挖方案的开挖工序 2.2 计算模型及参数设置 由于该公路隧道规模较大，本文仅选取隧道进口段围岩条件较差的典型断面ZK1+375～ZK1+405进行数值模拟分析，该断面所处的围岩级别为Ⅴ级，其物理力学参数主要来源工程地质勘察报告和设计报告，见表1。 表1 围岩及隧道结构支护参数项目弹性模量/GPa泊松比重度/(kN·m-3)黏聚力/kPa内摩擦角/(°)粉质黏土0.080.3014.324.315.1全风化凝灰岩1.100.3018.325.314.8喷射混凝土250.2522——锚杆2000.3076—— 考虑边界效应，根据现场工程的实际情况和圣维南原理，隧道边界取3倍洞径宽度作为有限元分析范围，计算模型的左右边界施加水平方向约束，底部施加竖向约束，上边界为自由边界[11]。分析中围岩体选用理想弹塑性本构关系，其屈服准则采用莫尔-库伦(Mohr-Coulomb)强度破坏准则，其它支护结构如锚杆、喷混凝土等均使用弹性本构关系，其中钢拱架采用I25b拱架，纵向间距0.5 m×0.5 m，钢拱架的作用采用等效方法就行分析，即将工字钢的弹性模量折算到喷射混凝土上，计算公式为： E=E0+Ag×Eg/Ac (1) 式中：E为折算后混凝土弹模；E0为原混凝土弹性模量；Ag为钢拱架截面积；Eg为钢材弹模；Ac为喷射混凝土截面积。 由于该隧道埋深较浅，计算初始应力场仅考虑自重应力场。最终，由MIDAS