基坑三维渗流对紧邻区间隧道影响的数值分析.docVIP

基坑三维渗流对紧邻区间隧道影响的数值分析 胡国新刘庭金陈俊生邓飞皇莫海鸿 摘要：基坑工程施工改变地下水渗流场，导致紧邻区间 隧道受力和变形发生改变。为此，应用数值分析手段，采用 流一固耦合计算模式，对某基坑施工过程的渗流场和变形场 进行三维模拟，分析了场地的水位下降规律和区间隧道的变 形规律。并通过建立能反映止水条和螺栓作用的三维管片 环实体计算模型，计算了管片变形和接头纵缝张开量。研究 结果表明：该基坑施工诱发的最大水位下降约为2m;水位下 降2m时管片的纵缝张开增量为，为纵缝张开控制总量的％。 因此，认为该基坑施工诱发的三维渗流对紧邻地铁区间隧道 的影响较小。 关键词:基坑工程三维渗流数值模拟地铁隧道水位下降 管片环接头张开量 0前言 某深基坑工程紧邻广州地铁一号线某地铁车站和区间 隧道，由于该基坑工程采用地下连续墙逆作法施工，基坑的 侧向位移对紧邻区间隧道的影响相对较小。因此，需重点分 析该基坑工程施工诱发的三维渗流对该区间隧道的影响。 为此，首先采用流一固耦合计算模式，对该基坑施工引 起的地下水渗流场和区间隧道变形场进行三维数值模拟;然 后基于三维渗流计算结果，通过建立能反映止水条和螺栓作 用的三维管片环实体计算模型，分析管片环变形和接头纵缝 张开量，为评估基坑施工对紧邻区间隧道的结构安全提供依 据。 1三维渗流计算 有限元模型及其参数 基坑附近区域场地的主要地层资料为:上部4 m杂填土， 中间10m左右的砂土层，下部1 8m左右的强风化砂岩，再往 下则为中风化砂岩。场地地下水埋深约为_2m。地下连续墙 深为21m。 根据该基坑工程及紧邻区间隧道的相关资料，对实际场 地的几何形状和地层分布进行适当简化，建立三维渗流计算 模型。地层分布如图1所示，基坑与地铁车站、区间隧道的 三维模型如图2所示。整个计算模型沿水平方向均为250 m， 沿垂直方向为40 m。计算模型中地层、地下连续墙及盾构 隧道均采用实体单元，共6 6155个单元，72 125个节点。 一般情况下，中风化砂岩可视为不透水层和不可压缩层， 因此，在计算时将下部约10m深的中风化砂岩视为强风化砂 岩，这对渗流和变形计算结果都是偏安全的。 表1为场地地层的渗透系数和力学计算参数。由于杂 填土和砂土的渗透系数远大于强风化砂岩的渗透系数，因此 当地下水绕地下连续墙发生渗流时，相对不透水层强风化砂 岩成为渗流计算的关键控制土层。根据现场分段常水头注 水试验，获得该场地强风化砂岩的渗透系数为x 10_7~X 10- 8m/s。考虑到强风化砂岩中裂隙发育和分布的不确定性，在 三维渗流计算中，将强风化砂岩的渗透系数放大约10倍，即 2 X10-6m/so 在进行三维渗流场和应力场耦合计算时，假定： 1）地下连续墙墙体本身不发生位移； 2 ）地下连续墙墙体和接头部位不发生渗漏； 3）不考虑基坑开挖诱发荷载释放对地层和地铁区间隧 道的影响，而仅考虑地下水渗流的影响。 渗流计算结果及其分析 为分析地下水三维渗流及隧道变形计算结果，建立如图 3所示的分析点、线示意图。 图4为横向分析线B水位下降曲线。从渗流场计算结 果可知，该基坑施工在区间隧道附近诱发的最大水位下降约 为2m。 另外，对强风化砂岩实测渗透系数放大100倍，即2X 10-5m/s进行了三维渗流计算，此时基坑施工在区间隧道附 近诱发的最大水位下降约为4m。 隧道变形计算结果及其分析 根据区间隧道变形计算结果可知：区间隧道最大变形部 位发生距离地铁车站北端约20 m处，且下行线隧道变形量大 于上行线，其中下行线隧道最大竖向位移为，上行线隧道最 大竖向位移为。 2管片环受力与纵缝张开量 计算模型及参数取值 水位下降导致作用在管片环外壁的有效压力增大，水压 力减小。图5示水位下降导致管片环外壁的压力增量。其 中Y为地下水重度，H为水位下降深度,K0为静止土压力侧 向系数，可根据砂土的有效内摩擦角确定，取值为。 为计算水位下降与管片环变形和接头张开量之间的关 系，建立管片环实体计算模型，如图6所示。在模型中充分 考虑了纵、环缝止水条及纵、环向螺栓的影响，如图7所示 并在管片环外侧设置lm厚的土体，采用接触分析考虑土体 对管片环的地层抗力作用，计算时通过将土体的弹性模量进 行折减以反映土体的扰动。 计算结果及分析 图8为水位下降2m时管片环的三维变形，变形放大120 倍。图9为水位下降与管片环接头纵缝最大张开量关系曲 线。根据管片环三维变形计算，管片环两侧腰部位外侧张开， 且水平方向变形略大于垂直方向变形。因此，水平方向收敛 位移成为控制管片环纵缝张开量的主导因素。水位下降2m 时管片环接头纵缝最大张开增量可控制在，为接头纵缝张开 总控制量3 mm（广州地铁总公司提供）的％。 3结论及建议 采用流一固耦合模式，计算了基坑施工