渗流—应力场耦合作用下苏州工业园区某地下车库基坑的变形分析

摘 要：地下水渗流对基坑变形的影响成为当前研究的热点，以苏州工业园区某地下车库基坑为例，采用ABAQUS模拟基坑开挖及支护过程，分析基坑开挖过程中的变形及渗流场规律。结果表明：在开挖间歇期的坑外地表沉降量均比同期开挖结束后的沉降量要小，而基底隆起量比同期开挖结束后的隆起量要大。每步开挖间歇结束时，围护墙的水平位移有所减小。随着开挖的进行，围护墙周围的水头等势线越来越密，地面沉降形状为下凹的盆地形状。

关键词：基坑 渗流-应力耦合分析 有限元模拟

中图分类号：TU43 文献标识码：A 文章编号：1672-3791（2013）05（b）-0039-04

基坑开挖时，坑内外通常存在着水头差，地下水将在坑内外水头差作用下发生渗流。基坑开挖过程是地下水渗流与岩土变形动态耦合的过程。利用渗流-应力耦合理论研究开挖过程中地下水的渗流形态和孔隙水压力场的分布，分析地下水渗流对基坑稳定性的影响具有重要意义。

近几年来，许多研究者在分析渗流稳定问题时，引入了渗流场与应力场的关系，即渗流-应力耦合关系，并在岩土工程的各个领域取得了一定的成果和进展，渗流-应力耦合问题已成了研究的热点问题。谢兼量[1]进行了渗流应力耦合条件下的海堤边坡稳定性研究；贾善坡等[2]进行了泥岩隧道施工过程中，渗流场与应力场完全耦合的损伤模型研究；张巍等[3]对大型地下洞室群围岩进行了应力-损伤-渗流的耦合分析；张媛媛[4]，苗丽等，周建国等[6]在土坝的渗流场与应力场的耦合应用方面的研究获得了一些进展；王强等[7]，杨永恒[8]，郭娟[9]，周舒威等[10]基于渗流-应力耦合对尾矿坝的稳定性进行了研究；李筱艳[11]、纪佑军等[12]采用渗流-应力耦合分析，求解基坑的渗流场以及位移场。

本文结合苏州工业园区星海街站南北两侧公共地块地下车库项目，利用ABAQU

S有限元软件进行了基坑工程在渗流-应力耦合作用下的变形分析，可为基坑工程的设计和施工提供参考。

1 工程概况和地质条件

星海街站南北两侧公共地块地下车库场地，位于苏州工业园区星海街及其以西、苏华路南北两侧的公共地块内。项目主要包括绿地、下沉式广场、地下两、三层停车场及局部商业、预留地下通道，部分地面建（构）筑物。地下车库北基坑北侧5 m为苏雅路，东侧基坑边线位于星海街慢车道上，距离建园大厦4～7 m，南侧基坑边线位于苏华路慢车道上，西侧较为空旷，现为绿化草坪。南基坑北侧基坑边线位于苏华路慢车道上，东侧基坑边线位于星海街慢车道上，距离星海大厦3～7 m，南侧基坑边线以南7 m有一近东西走向河道（相门塘），C25孔以南现为在建工地，西侧较为空旷，为绿化草坪。地下室结构采用现浇钢筋混凝土框架体系，地下室结构底板（无论地下两层或三层）均处于同一标高，地下室结构顶板上部覆土（至自然地面）地下三层处约1.3 m，地下两层处约3.0 m，基坑深度为自然地面以下约14.5 m。

3 有限元分析

3.1 模型尺寸与计算参数

2.2 边界条件

模型左边界为轴对称边界条件，右边设置水平方向位移约束，底部边界设置水平、竖直方向位移约束。模型右边界假定孔隙水压力不发生变化为定水头边界，即水源源不断地补充，模型底部为不透水边界，模型左边界孔隙水压力随着降水水头的变化而变化。

2.3 有限元模拟的实施步骤

3 计算结果与分析

3.1 降水开挖引起的基坑变形

地下水渗流作用贯穿于基坑工程整个施工过程中，无论是降水阶段，还是开挖加撑阶段。而基坑的变形主要就是桩体的变形、桩后地面沉降和基坑底部隆起。图2为二维基坑模型降水开挖过程中的土体沉降变形图。可以发现：（1）随着开挖过程的不断进行，桩后地面的沉降量和基坑底部的隆起量都在不断增加。（2）在考虑渗流作用的情况下，基坑在开挖间歇期的沉降变形值均比同期开挖结束后的沉降变形值要小，而基坑隆起量却比同期开挖结束后的隆起量要大，但两者的变化量不大。

3.2 支护结构水平位移

图3为不同开挖阶段围护墙的水平位移变化曲线，每步开挖和间歇期的趋势一致，最大水平位移发生的位置随着开挖深度的增大而逐渐下移，与各步开挖面基本保持一致，最终在基坑底面处达到最大。从图中还可以看出：（1）随着开挖深度的增加，墙身下部的位移随之增大，而顶部位移有减小的趋势，其最大位移发生在第三次开挖结束之后。（2）在每步开挖间歇结束时围护墙的水平位移均有所减小，这是由于超静孔隙水压力的消散和渗流造成坑外水压力的减小，使得墙体的位移有所回落，有利于基坑的稳定性。

3.3 基坑底隆起变形

3.4 孔隙水压力变化

基坑降水开挖引起的渗流场总静孔隙水压力如图5所示，a为地应力平衡时的初始孔隙水压力分布图，此时孔压不存在超静孔压，孔压分布和重力场平衡，所以不发生渗流；b～d为每步开挖结束时和开挖间歇期结束时土体中的总静孔压的分布图，由于降水开挖后基坑内外产生水头差变化，开挖卸载也同时在土中产生负的超静孔隙水压力，土中水在重力势（水头差）和压力势（超静孔压）的共同作用下发生渗流。等势线从入水边界到出水边界逐渐变密集，在围护墙底部附近分布最密，随着开挖的进行，围护墙周围的等势线也是越来越密，水力梯度也就越来越大，故此时地下水的流速也就越大。

3.5 地面沉降

图6是支护墙外地面沉降的变化曲线，可以看出，三层基坑开挖引起的地面沉降，从整个变形趋势来看，它们的形状都是一个类似于向下凹的盆地形状，最大沉降量并不是出现在坑壁，而是在离基坑一定距离的地方，再随着距支护结构边缘距离的增加，沉降变形值逐渐减小。所以，在坑外附近有建筑物的情况下，要严格注意地表的沉降，防止产生允许值之外的不均匀沉降而造成严重损失。

4 结论

在考虑渗流作用的情况下，基坑在开挖间歇期的坑外地表沉降量均比同期开挖结束后的沉降量要小，而基底隆起量却比同期开挖结束后的隆起量要大。围护墙在开挖过程中的最大水平位移发生的位置随着开挖深度的增大而逐渐下移，与各步开挖面基本保持一致，最终在基坑底面处达到最大。而在每步开挖间歇结束时围护墙的水平位移均有所减小。随着开挖的进行，围护墙周围的水头等势线越来越密，流速也越来越大。地面沉降形状为下凹的盆地形状，最大沉降量不是出现在坑壁，而是离基坑一定距离的地方。

参考文献

[1] 谢兼量.考虑渗流应力耦合的边坡稳定性分析[D].南京：河海大学，2007.

[2] 贾善坡，陈卫忠，于洪丹.泥岩隧道施工过程中渗流场与应力场全耦合损伤模型研究[J].岩土力学，2009，30（1）：19-26.

[3] 张巍，肖明，范国邦.大型地下洞室群围岩应力-损伤-渗流应力耦合分析[J].岩土力学，2008，29（7）：1813-1818.

[4] 张媛媛.ANSYS在土坝渗流场和应力场及其耦合分析中的应用研究[D].南京：河海大学，2006.

[5] 苗丽，郭雪莽，王复明.基于应力场与渗流场耦合的土坝稳定性分析[J].人民黄河，2007，29（9）：75-77.

[6] 周建国，李淼，郭雪莽.东湖书库非均质土坝渗流场与应用场耦合分析[J].水利电力机械，2007，29（3）：28-30.

[7] 王强，鲁柄强，王水平.尾矿坝渗流-应力耦合场的有限元分析[J].现代矿业，2010（3）：74-77.

[8] 杨永恒.渗流场与应力场的耦合分析及其在尾矿坝工程中的应用[D].西安：西安理工大学，2006.

[9] 郭娟.木梓沟尾矿堆积坝的物理力学特性及坝体渗流场与应力耦合分析[D].西安：西安理工大学，2009.

[10] 周舒威，李庶林，李青石，等.基于渗流-应力耦合分析的野鸡尾尾矿坝稳定性研究[J].防灾减灾工程学报，2012，32（4）：494-501.

[11] 李筱艳.基于位移反分析的深基坑渗流场与应力场的完全耦合分析[J].岩石力学与工程学报，2004，23（8）：1269-1274.

[12] 纪佑军，刘建军，薛强.基坑开挖中渗流-应力耦合模拟[J].岩土力学：增刊，2007，28：630-634.