根据abaqus无渗流偏压隧道开挖分析(包含地应力平衡分析).pdfVIP

基于 ABAQUS 的无渗流偏压隧道开挖分析（包含地应力平 衡分析） 1 隧道建模及地应力平衡 1.1 工程概况 本模型截取的一段隧道通过山体坡度 30°左右一侧，隧道开挖及初衬断面为五心圆各项 参数如表 2.1 所示。隧道跨度 13.36m，高 11.71m，偏压部分覆盖层厚度约为 20m 。隧道区 地下水主要为第四系松散层孔隙水和基岩裂隙水，垂直隧道中轴线向东 120 米勘察孔水位为 3 65 米，向西 120 米勘察孔水位为 10 米。施工过程中隧道单位正常涌水量 3.43m /d ，中等富 水。具体围岩初衬参数将在后文建模时给出。 表 2.1 隧道断面参数(单位：cm) r[1] R[1] r[4] R[4] h[1] h[2] h[3] a b h H B 603 653 1315 1370 212 191 912 0 220 815 1171 1366 1.2 隧道建模及偏压判定 结合隧道工程勘察报告和规范取围岩和初衬参数值如表2.2 所示，隧道初衬厚度为30cm 。 表 2.2 围岩和初衬参数取值 弹性模量/GPa 泊松比 密度/Kg/m3 内摩擦角/° 内聚力/MPa V 级围岩 1.8 0.38 2200 25 0.3 初衬 25 0.2 2400 / / 文献指出，取模型边界为隧道开挖直径 8 倍时，其地震和静力计算对隧道的影响可以忽 略不计。本模型取隧道开挖直径的 10 倍，力求将边界影响降低到最小，同时兼顾计算工作 量。建模时先简化山体坡脚 30°一侧为一顶角 120°的圆锥，按照隧道实际位置确定隧道后， 在山体上切下符合隧道开挖直径 10 倍的计算模型（如图 2.2 ），这样会使模型上表面近似山 体的曲面，较符合实际地形条件，模型自前向后偏压角度略有减小。与直接将上表面简化为 平面的模型相比，其计算结果和实际工程的符合性有较大提升，对隧道内部应力的分布也反 映的更为准确。本章将对无渗流条件下的静力偏压情况进行模拟，以便和后文中渗流作用下 偏压隧道的应力场变化进行比较。模型上隧道圈线和穿越隧道的斜线是为方便后续渗流计算 所做分割面，不影响计算。 图2.2 模型切取示意图 模型尺寸如上文所述，首先不考虑渗流场，仅考虑重力作用。岩体采用摩尔-库伦模型， 初衬选用弾性模型。其中锚杆注浆加固区影响范围约为 4m ，加固效果在软件模拟中以提高 影响范围内岩土参数的方式进行处理。将加固区围岩的弹性模量、密度、粘聚力和内摩擦角 提高 10%~20%，泊松比降低 10%左右。建立模型如图 2.3 所示，对模型前后左右四个面和底 面进行位移约束，上表面自由，模型整体施加重力作用。采用六面体单元形状，C3D8R 单元， 其中山体和隧道岩体共划分 18237 个单元（如图2.3-a ）；初衬为图2.3-b 中红色区域，划分4100 个单元。隧道开挖计算前利用 MODEL CHANGE，REMOVE 功能移除初衬单元，开挖后重新激 活，不影响前置计算。模型网格划分较密，虽加大了工作量，但可以增加计算精度，使结果 更加准确。 （a ）模型整体示意图 （b ）隧道部位放大图 图2.3 模型建成图 结合《铁路隧道设计规范》、《公路隧道设计规范》可知影响偏压隧道结构应力分布的主 要因素有横坡坡度和最大埋深。其中最大埋深可由普式理论进行确定：普式理论提出当隧道 埋深至一定深度后，隧道开挖后在其上方将形成自然拱，仅自然拱内的围岩松散体对支护结 构施加压力，自然拱外的围岩压力由自然拱承担。由此可知因为