基于有限元的顶管工作井三维分析.docxVIP

基于有限元的顶管工作井三维分析 1 工作井的变形和应力 隧道法是隧道和大型地下管道的重要非开挖技术，包括重要的地下管道技术，其中，井的施工在施工中起着重要作用，因此，必须承受管道前面膨胀的所有水平推力。水平推力随顶管直径和长度而增大，一般为1000～2500t，最大可达3000t以上。在巨大水平推力作用下，工作井可产生明显的向后倾斜位移，并在井壁特别是靠背后壁中出现高应力分布。工作井的倾斜位移影响顶管的定位精度和纠偏，可使顶管偏离设计轴线；同时过大的倾斜位移使地面隆起，引起邻近管线和路面破坏等一系列问题。井壁中的高应力则常导致井壁开裂漏水，无法正常工作，甚至使施工半途而废，造成重大经济损失。因此，控制工作井的变形和应力是顶管法顺利施工的重要条件之一。 为控制顶管工作井的变形和应力，在设计中进行相应的变形计算和应力分析是非常必要的。然而，承受水平荷载工作井的变形和应力计算比较复杂，我国相关设计规范也未作出相应说明和规定，有关的研究工作也很少见，而实际设计中都未对工作井的变形和应力进行分析。因此，目前工作井的设计主要依据工程经验，存在很大的盲目性和随意性，使得完成的设计可能不满足施工要求而造成失效，也可能过于保守而造成浪费，只有把科学的计算分析方法与实践经验相结合，才能给出既满足施工要求，又经济合理的优良设计。 水平荷载下的工作井变形和应力分布是典型的三维问题，井体的位移和应力与嵌岩情况、周围土层性质、井体结构参数等诸多因素有关。本文采用了MARC有限元软件，对矩形工作井的变形和应力进行了三维计算，给出了工作井后壁位移和应力结果，分析了嵌岩与否、井周土层弹性模量及覆土厚度等参数对变形和应力的影响。计算结果表明，工作井不入岩时，随着井周土层弹性模量增大，工作井位移明显减小，但对应力的影响不大；工作井入岩时，井周土层弹性模量的改变对工作井的位移和应力影响较小；覆土厚度H的改变对井后壁水平位移的影响较小，而工作井的水平压应力随覆土厚度H减小大致成比例减小。文中指出工作井不入岩时，通过土体加固提高井后靠背土层的弹性模量，可有效控制工作井的变形，该计算分析可为工作井设计提供有益的参考，对完善顶管施工设计理论具有指导意义。 2 井周土体单元网格划分 图1为顶管工作井受力的断面简图及其俯视图, 其中L、B分别为工作井平面的长度和宽度;d1、H1分别为连续墙的壁厚和高度;d2、H2分别为后背墙的壁厚和高度, d2=d1;d3、b、H3分别为钢垫块的厚度、宽度和高度;d4为井底板的厚度。井周土体为层状, 各层厚度为hi, 重度为, 弹性模量为Ei, 泊松比为。根据工作井的受力和变形特性, 其变形和应力计算属于三维问题, 为此取图2所示的三维有限元计算模型。在该计算模型中, 对于井周土体在水平面内截取井体尺寸的3～5倍长宽L·B, 而深度方向截至基岩作为计算范围。单元网格划分从井位中心向计算模型外边界由密到疏逐步过渡;对于水平推力作用的井后壁及后壁外土体是变形大和应力集中区域, 因此适当加密该区域的单元网格。计算还考虑了内支撑设置情况。土层-井体的整体计算模型及其单元网格划分如图2所示。计算中土体和岩石选用三维实体单元, 井体、后背墙、垫块及底板选用壳单元, 内支撑选用三维梁单元。模型底部边界 (岩层) 固定;模型侧面边界水平面内无位移, 竖向则可自由变形。 3 计算与分析 3.1 井后壁外土体 设井周土为均匀线弹性介质，土层，厚度为14.0m，计算模型参数分别为：水平顶力1500t，井深10m, L·B=8m×8m，底板d4=0.5m；顶进洞口直径D=4m，洞口中心位置在地面以下7m处，顶管顶面覆土厚度5m；后靠墙d2=0.5m，宽8.0m, H2=8m，为顶进洞口外径的2倍，嵌入底板以下1.0m；千斤顶后的钢垫板b·H3=4.0m×4.0m, d3=0.25m；工作井采用地下连续墙支护，厚600mm，嵌固深度为6m。 利用有限元程序建模截取土体的计算范围时，考虑到直接支承水平顶力的井后壁外土体对井的变形和应力影响更明显，因此截取井后壁外土体30m，而前方土体为20m，井后壁外土体的单元网格也适当加密，计算深度取岩层下2m。岩层参数为：。设有2道内支撑，分别位于-2m和-5m处，采用外径600mm、壁厚10mm的钢管。 3.2 井后壁应力分布 工作井是否入岩及井周土体性质对井体变形和应力分布的影响极为明显，为了分析上述两者对顶管工作井变形和应力的影响，笔者先分为工作井入岩、工作井不入岩（此时设土层厚度为16m）两种情况，其中又令土体弹性模量E依次取6×103, 1×104, 2×104, 3×104, 6×104kPa等5种情况进行分析，给出了矩形工作井变形和应力的三维有限元计算结果。图3～4分别为工作井入岩和不入岩时，后壁中线处的变形与应力计算