顶管开挖面支护压力影响研究
0 引言
顶管工程中, 土压平衡式顶管机工作的关键是控制土仓支护压力与所处土层地下水压力及土压力处于平衡状态
在隧道工程中, 关于开挖面极限支护压力和稳定性研究较多, 朱伟等
文章采用数值模拟方法, 研究开挖面水平位移与支护压力的关系, 并对地表不同测点的竖向位移进行分析, 探究了支护压力作用下开挖面及地表位移的影响因素, 为顶管施工控制开挖面稳定和地表变形提供相应参考。
1 开挖面受力分析
Leca和Dormieux
为了模拟因支护压力不足造成的超挖, 须满足:σγ>σT;建立二者之间的关系, 引入支护压力比
土层土压力的选取关系到开挖面初始支护压力的大小。余彬泉等
极限平衡理论计算出的土压力多依赖于土压力系数, 若初始支护压力在pp和pa之间取值, 当埋深较大时, 范围较广, 结果误差偏大, 在计算过程中可能会出现支护压力比还未变化, 开挖面在不平衡力的作用下已发生位移, 使得模拟结果不准确。模型试验中, pc的确定受过程控制的影响较大, 存在一定的误差。
本文研究中, 开挖面的初始支护压力取地层平衡时的原位静止土压力, 为避免理论计算与数值模型求解误差对后期模拟过程的影响, 静止土压力取模型计算结果。
2 顶管开挖面支护压力数值模拟
2.1 有限单元法原理
常规分析方法在处理岩土开挖问题上较为复杂, 基于合理的假设, 通过有限元方法对施工过程进行模拟可反映出开挖过程中土体的变形规律。有限单元法
考虑顶管施工条件的复杂性以及研究侧重点主要在开挖面, 本文数值模拟作如下假定: (1) 不考虑顶管过程的注浆减摩作用, 一次性开挖, 并施加管节支护; (2) 土体是由均质、单一材料构成的各向同性体, 采用Mohr-Coulomb本构模型, 用实体单元模拟; (3) 管节由混凝土材料构成, 采用壳单元模拟。
2.2 模型建立及参数选取
研究表明, 圆形隧道开挖有限元模拟中, 5倍管径以外土体的应力-应变影响在1%以下, 模型建立如下: (1) 土体模型长22m、宽10m、高12m, 采用C3D8R单元类型, 共建立10 032个单元; (2) 管节外径D为2m, 采用S4R单元类型, 共建立270个单元, 埋深管径比R=C/D, 则管节埋深为1m, 埋深为1m时的模型如图3所示, 模型计算参数取值如表1所示。
2.3 数值模拟过程
1) 管节顶进至8m, 在开挖面上施加与地层初始静止土压力大小相等的梯形支护作用力。
2) 逐渐减小开挖面支护压力, 将支护压力比变化量很小, 而开挖面水平位移发生突变时的支护压力作为极限支护压力;将迭代过程不收敛时的支护压力作为开挖面失稳破坏压力。
2.4 模拟结果讨论与分析
支护压力比为0.476时的地表位移云图如图4所示, 此时开挖面水平位移发生突变, 地表竖向位移有向开挖面前方区域聚集的趋势, 范围分别以管节上下表面为边界向外扩散, 大致呈“烟囱”状。
1) 支护压力比与开挖面水平位移关系
为了解支护压力减小过程中开挖面的水平位移变化趋势, 取开挖面中心点为参考点, 结果如图5所示。
从图5中可以看出随着支护压力的减小, 开挖面位移变化过程大致可分为3个阶段: (1) 第1阶段支护压力比满足0.476<λ<1时, 开挖面基本维持稳定, 水平位移几乎为0; (2) 第2阶段支护压力比λ=0.476时, 位移发生突变, 此时的支护压力即为开挖面极限支护压力, 约为9.51kPa; (3) 第3阶段支护压力比λ<0.476, 支护压力以较小的变化量就可引起开挖面较大的水平位移, 随着支护压力的减小, 计算过程继续迭代, 直至最后不收敛, 此时支护压力比约为0.4, 将计算过程不收敛时开挖面的位移作为失稳的状态指标, 约为0.547m。
2) 支护压力比与地表竖向位移关系
开挖面支护压力减小到极限支护压力时, 开挖面水平位移发生突变, 引起地表产生附加竖向位移, 为了解距开挖面不同位置的竖向位移, 分别取距开挖面0, 0.5, 1, 1.5, 2, 2.5m这6个测点, 结果如图6所示。
![图6 沿顶进方向距开挖面不同测点竖向位移Fig.6 Vertical displacement along the jacking direction from different measuring points](/User/GetImg.ashx?f=SGJS/10455//SGJS201901027_07800.jpg&uid=WEEvREcwSlJHSldTTEYzVTFPV2k1VVcwUUppcDkrK3NweHNuUkRVSHAzQT0=$9A4hF_YAuvQ5obgVAqNKPCYcEjKensW4IQMovwHtwkF4VYPoHbKxJw!!)
图6 沿顶进方向距开挖面不同测点竖向位移Fig.6 Vertical displacement along the jacking direction from different measuring points
从图6可以看出, 距开挖面不同测点的地表竖向位移不同, 但趋势大致相同, 先随着支护压力比的减小, 发生微小变化, 后出现拐点, 开挖面达极限支护压力, 此后随着支护压力比的减小, 竖向位移呈现不同程度的增大。地表开始出现较大竖向位移时的极限支护压力比要比图5结果小, 约为0.31。地表沉降量最大位置不在开挖面正上方, 而在其前方1.5m处, 最大沉降量约为58mm。极限支护压力之前, 地表竖向位移几乎为0;极限支护压力之后, 随着极限支护压力的减小, 竖向位移急剧增大。
为了解地表沉降范围及趋势, 在支护压力比为0.31时, 沿顶进方向取一路径, 结果如图7所示。
从图7中可以看出, 埋深为1m时, 地表竖向位移主要集中在开挖面后方2m和开挖面前方4m范围内, 约为3D (D为管节直径) 。沉降曲线呈高斯分布, 其中最大沉降量发生在开挖面前方2m范围内。分析可知, 开挖面前后方分别在原位土层作用和管节支护作用下维持稳定, 位移基本为0;地表位移向开挖面区域聚集, 最终形成稳定的沉降范围。
3 开挖面及地表位移影响因素分析
在本文研究中, 土体遵循Mohr-Coulomb屈服准则, 开挖面滑动变形受材料黏聚力和内摩擦角值的影响;Chambon和Corte研究发现, 不同埋深-管径比对开挖面的变形破坏也有影响。基于此, 分别研究管节埋深C、土体黏聚力c、土体内摩擦角φ对开挖面位移及地表沉降的影响。
3.1 埋深的影响
控制土体参数相同 (c=1.1kPa, φ=10°) , 分别对不同埋深的开挖面支护压力进行模拟, 参数取值如表2所示, 模拟结果如图8所示, 从图8中可以看出, 不同埋深的支护压力比-位移曲线呈现相同的趋势, 开挖面位移变化同样经历上文讨论的3个阶段。但不同埋深的极限支护压力比不同, 达到极限支护压力时开挖面的水平位移也不同。埋深越大, 极限支护压力比越大, 达到极限支护压力时的水平位移也就越大;反之亦然。埋深为1m时, 开挖面最先达到极限支护压力, 支护压力比约为0.476, 水平位移约为3.29mm, 结果如表3所示。
埋深分别为1, 2, 3m时的地表沉降曲线如图9所示, 从图9中可以看出, 不同埋深, 地表沉降趋势大致相同, 但沿纵向呈现不同的沉降范围, 随着埋深的增大, 沉降范围逐渐扩大;地表沉降量随埋深的增大逐渐变小, 埋深为3m时的沉降量约为1m时的1/2。
3.2 内摩擦角影响
控制管节埋深和土体黏聚力不变 (C=3m, c=1.1kPa) , 分别对土体采用不同内摩擦角的开挖面支护压力进行模拟, 参数取值如表4所示, 模拟结果如图10所示。从图10中可以看出, 同一埋深, 不同内摩擦角下, 摩擦角越大, 极限支护压力比越小, 说明开挖面位移发生突变时的支护压力越小, 在支护压力减小过程中, 开挖面的稳定性越好。摩擦角为20°时, 支护压力比约为0.05, 开挖面水平位移约为13.31mm, 当支护压力比减小为0时, 迭代过程收敛, 计算结束, 此时开挖面水平位移约为65mm, 接近于摩擦角为5°时的突变位移 (极限支护压力比时的位移) , 结果如表5所示。
![图1 0 不同内摩擦角开挖面位移曲线Fig.10 Displacement curve of excavation face with different internal friction angle](/User/GetImg.ashx?f=SGJS/10455//SGJS201901027_09500.jpg&uid=WEEvREcwSlJHSldTTEYzVTFPV2k1VVcwUUppcDkrK3NweHNuUkRVSHAzQT0=$9A4hF_YAuvQ5obgVAqNKPCYcEjKensW4IQMovwHtwkF4VYPoHbKxJw!!)
图1 0 不同内摩擦角开挖面位移曲线Fig.10 Displacement curve of excavation face with different internal friction angle
图11为沿顶进方向内摩擦角影响下地表沉降曲线, 从图11可以看出, 不同内摩擦角对地表的纵向影响范围大致相同, 地表竖向位移随内摩擦角的增大而减小。
由以上结果分析可知, 土体的内摩擦角是影响开挖面稳定的重要参数, 内摩擦角较大时, 支护压力的减小对地表位移影响较小, 且计算过程容易收敛。
3.3 黏聚力的影响
控制管节埋深和土体内摩擦角不变 (C=3m, φ=10°) , 分别对土体选取不同黏聚力的开挖面支护压力进行模拟, 参数取值如表6所示, 模拟结果如图12所示。从图12可以看出, 随着土体黏聚力的增大, 支护压力比与位移关系逐渐呈线性变化, 黏聚力越大, 线性变化趋势越明显, 直至计算收敛;黏聚力越小, 较大的支护压力也不能阻挡开挖面的位移变化, 计算过程越不容易收敛。当黏聚力为3.9kPa时, 线性变化趋势最为明显, 此时支护压力比减小为0, 开挖面水平位移约为46mm;当黏聚力为2.5kPa时, 支护压力比减小到接近于0时的位移大约为302mm。
图13为沿顶进方向的地表沉降曲线, 从图13中可以看出, 地表沉降量因黏聚力的不同呈现不同的沉降量, 但沉降槽宽度大致相同。土体黏聚力越大, 地表沉降量越小, 黏聚力为3.9kPa时的地表沉降量最小, 约为3mm;黏聚力为1.1kPa时, 沉降最大位置在开挖面前方1m处, 沉降量约为0.109m。
对以上结果分析可知, 支护压力减小引起开挖面位移变化受到土体黏聚力的影响, 黏聚力越大, 支护压力对开挖面的作用越不明显, 结合本文模型, 黏聚力>3.9kPa时, 开挖面的稳定性较好, 地表沉降可得到有效控制。
4 顶管工程实例分析
某给水工程, 穿越城市西四环, 起点桩号为0+027, 终点桩号为4+763, 管线全长4 736m。属于第四系地层, 第四系上更新统主要以褐黄色黄土状粉土为主, 局部夹有黄土状粉质黏土、粉土和粉质黏土;第四系中更新统顶板埋深在5~30m, 主要土层为褐红色粉质黏土, 呈硬塑~坚硬状, 局部夹有黏土、粉土、中细砂、卵石层和胶结层。经过调研与现场勘察, 在0+303—0+388段要穿越西四环道和道路两侧居民生活区, 综合考虑地质条件、生活和交通等因素, 该段管道采用土压平衡式顶管机进行施工, 施工平面如图14所示。
供水线路管节采用直径400mm的钢管, 顶管机在2台200t液压千斤顶作用下从工作井开始顶进, 接收井回收, 管节紧随其后埋设于土层中, 管节埋深3.7m (地面距顶管机中心的高度) , 随着顶管施工进行, 深层土体因应力释放而发生位移, 开挖面支护压力过小会使得路面和临近建筑物产生不均匀沉降, 须进行实时监测并采取保护措施。在开挖面前方0.5, 1, 1.5, 2m处地表相应位置分别布置监测点, 图中加粗部分表示已埋设的管节, “+”表示顶管机到达时, 地表距开挖面不同位置布置的测点, 采用精密水准仪对地面竖向位移进行监测。土体静止土压力系数取0.55, 土体容重取18kN/m3, 通过理论计算, 开挖面土压力为36.63kPa, 顶进20m时顶管机压力表盘土仓压力为35.94kPa, 支护压力比约为0.98, 顶进前后地表竖向位移几乎无变化。将顶进到5, 10, 15, 20m时的开挖面分别定义为开挖面1、开挖面2、开挖面3和开挖面4, 顶管机到达不同开挖面位置时的土仓压力及测点竖向位移部分数据如表7所示。
对表7中的监测数据进行分析可知, 顶管机到达不同开挖面时的支护压力不同, 监测点的竖向位移也不一样, 支护压力越接近于管节所处土层的土压力, 地表竖向位移越小, 地表沉降控制效果越好;在同一开挖面, 同一支护压力下, 距开挖面越远, 地表沉降越小, 监测结果与上文模拟支护压力影响地表位移趋势吻合。
5 结语
1) 开挖面支护压力减小到极限支护压力时, 开挖面水平位移发生突变, 引起地表附加竖向位移, 极限支护压力的大小受土体内摩擦角和黏聚力影响较大。
2) 地表竖向位移随支护压力的减小而增大, 距开挖面不同位置的竖向位移有相同的变化趋势, 地表竖向位移曲线呈高斯分布, 沉降范围和沉降量大小受管节埋深和土体参数的影响。
3) 通过数值模拟方法得出的支护压力对地表竖向位移影响趋势与施工现场监测数据结果吻合, 在实际工程中, 应合理控制顶管机土仓压力以达到控制地表沉降的目的。
[2]朱伟, 秦建设, 卢廷浩.砂土中盾构开挖面变形与破坏数值模拟研究[J].岩土工程学报, 2005, 27 (8) :897-902.
[3]黄正荣, 朱伟, 梁精华, 等.浅埋砂土中盾构法隧道极限支护压力与稳定研究[J].岩土工程学报, 2006, 28 (11) :2005-2009.
[4]刘宇, 唐春安, 谢欣.砂土地层盾构开挖面极限最小支护压力[J].东北大学学报 (自然科学版) , 2015, 36 (3) :433-436.
[5]乔金丽, 张义同, 许春彦.考虑渗流的盾构隧道开挖面稳定性分析[J].水文地质工程地质, 2009 (1) :80-86.
[6]唐晓武, 李姣阳, 邹金杰, 等.浅埋盾构隧道开挖面失稳发展过程模型试验[J].浙江大学学报学报 (工学版) , 2017, 51 (5) :863-869.
[7]周舒威, 夏才初, 葛金科, 等.黏土中超大直径顶管开挖面主动极限支护压力计算方法[J].岩土工程学报, 2013, 35 (11) :2060-2067.
[8]魏纲, 贺峰.砂性土中顶管开挖面最小支护压力的计算[J].地下空间与工程学报, 2007, 3 (5) :903-908.
[9] E LECA, L DORMIEUX. Upper and lower bound solutions for the face stability of shallow circular tunnels in frictional material[J]. Géotechnique, 1990, 40 (4) :581-606.
[10]王勖成, 邵敏.有限单元法基本原理和数值方法[M].北京:清华大学出版社, 1997.