中国资源正逐步集中于密集的城市地区,在繁华地带土体资源出现了紧缺的局面,使得基坑工程周围的施工工况更加复杂,因此, 对基坑开挖引起的周围环境的影响的研究更加紧迫。土体开挖对周围环境的影响主要体现在基坑外土体的移动以及土体的不均匀沉降。Migliazza等[1]利用分析和经验公式,并通过建立三维有限元模型,分析隧道在不同开挖阶段及支护阶段由开挖引起的地表沉降。Gordon等[2]提出了一种预测支护墙体最大挠度及最大地表沉降的半经验模型,研究表明,该模型可以准确地预测在软黏土地区由于支护开挖引起的支护墙体最大挠度及地表沉降。Hsieh等[3]提出了2种沉降剖面模型,并通过实际案例对比验证,研究表明可以通过沉降剖面曲线预测沉降的近似值。李琳等[4]利用分层总和法求出由降水引起的地面沉降,同时在经验假设的基础上,利用剪切位移关系求出由于围护结构对地面沉降的约束作用产生的沉降,把这两者叠加作为总的地面沉降。李琳等[5]研究了成层土中基坑开挖降水引起的土中应力变化及周围地表沉降的求解方法。王卫东等[6]根据上海地区大量工程实测数据提出了地表沉降曲线的预估方法并验证了其合理性,提出了预估基坑开挖引起的建筑物附加变形的计算方法。目前,大部分学者对基坑施工引起周围土体不均匀沉降主要集中在考虑单一因素,但工程实际中周围土体的沉降主要由降水、支护结构变形等其他因素共同引起的,只有把几个重要因素结合起来考虑才能更能有效地反映出实际。施成华等[7]推导出了一种考虑渗流作用的基坑降水地表沉降计算公式,利用叠加原理最终得出由开挖和降水引起的地表沉降分布计算公式。本文综合考虑了降水、支护结构变形与基坑底部隆起这3个因素对周围土体沉降的影响,利用降水漏斗曲线来对土体进行分层沉降计算,然后叠加算出由于降水引起的土体沉降,根据支护结构的刚度和类型,推导出由于支护结构变形和坑底隆起两者共同引起的土体沉降,最终把由于降水和支护结构变形及其坑底隆起引起的土体沉降叠加修正后,得出基坑周围土体沉降的理论解析解。
降水前,地下水位以下土体是一种由固体土颗粒和液体孔隙水组成的两相介质,土体所受到的力由固体土颗粒和孔隙水共同承担,土体所受的应力为土体有效自重γ′在土体中引起的应力。当土体内的水由于降水被疏干后,孔隙水压力便发生转移、消散,土体所受的应力为土体天然重度γ在土体中引起的应力,土体所承担的应力将增大,即有效应力增加。在浅部土层中进行疏干降水,基坑周围含水层的水不断流向基坑,经过一段时间之后,在基坑周围形成漏斗状的弯曲水面,即所谓的“降水漏斗”曲线[7]。如图 1,假设地基降水前地下水位于地表下h处,降水达到稳定状态后,地下水降至ho处。设在ho处一微单元dxdz,微单元所受应力为
式中:γ为土体天然重度,γ′为土体有效重度。
降水后,ho处的水被疏干,空隙水压力分散给土体颗粒承担,土体受到的应力为p′=γho,则土体在降水后有效应力增量为
式中:γw为水的重度。
f(x)为降水漏斗曲线。
式中:R为图 2中所示漏斗曲线与原水位线的交点,即抽水影响半径;h为原地下水位至地表的距离;hw为降水井中水面距离基坑底部的距离;rw为为降水井的半径;Hw为基坑外原地面至降水井水面的距离。
上述降水漏斗曲线方程中的x是距离降水井中轴线的距离,在计算降水引起的土体沉降时,坐标轴一般建立在支护结构处或便于计算的点处,且坑外地下水位一般是控制性降水,坑外地下水不一定降低在基坑底部以下0.5 m以下,故在计算中需要对f(x)中的x进行修正。
式中:xo为根据已知的实际坑外地下水位f(x),通过参数Hw、hw、h反推算出的一个定值, 其他符号意义同上。
基坑周围土体降水后,土体有效应力增加,由有效应力增加引起的沉降可以用式(5)计算
式中:Δhi为第i土层的厚度,Esi为第i土层土体的压缩模量。
如图 1中微单元dxdz产生的微元压缩可表示为
采用分层总和法计算由降水引起的土体沉降,则在原地下水位与降水漏斗曲线之间的土体沉降为
式中:m为原地下水位与降水后地下水位之间的土层数;Zi为第i层土顶部距离地表的深度;αi为第i层土范围内平均附加应力系数,可查《建筑地基基础设计规范》。hi为第i层土的厚度。
降水漏斗曲线上部因降水而增加的平均竖向附加应力,本应为三角形分布形式,但为计算方便取其平均值并假设为矩形分布,则由于附加应力及其渗流作用引起降水曲线下的土体沉降量为
式中:n为原地下水位面至支护桩底面之间的总土层数。
综上所述,采用分层总和法计算由降水引起的基坑外部土体沉降,则距离坑边任意距离x处的由降水引起的基坑外土体沉降为
式中:ψc为沉降计算经验系数,根据地区沉降观测资料及经验确定,无地区经验时可采用《建筑地基基础设计规法》表 5.3.3的数值;αi为矩形均布荷载作用下地基中竖向应力分量的平均附加应力系数。
基坑开挖引起坑外土体沉降主要考虑2个因素:支护结构变形及基坑底部隆起。根据分析大量学者[8-12]对基坑外地表沉降的研究,对基坑开挖引起的坑外土体沉降计算公式采用
并作以下假设:
1)基坑在沿支护结构纵向方向上,地表土体在理想状态下沉降曲线应按照支护结构中部截面对称。有学者曾假设纵向沉降曲线按照高斯曲线分布,但高斯分布曲线在距支护结构中部对称面一定距离后,收敛非常迅速,而实际沉降的收敛速度并没有那么迅速,故假设纵向沉降曲线按照抛物线分布更符合实际沉降曲线,如图 3。
2)研究表明,基坑支护结构后土体横向沉降在距支护结构一定距离后,沉降曲线迅速收敛,基坑开挖对此处的土体沉降影响表现不再明显,故采用高斯函数来近似土体横向沉降曲线更加合乎土体实际沉降。
3)在大部分对基坑外地表沉降的研究中,均未考虑基坑底部隆起的因素,得到的沉降曲线大部分是三角形,抛物线或分段直线函数。考虑到基坑底部隆起对坑外地表沉降也有影响,故采用横向高斯分布函数与纵向抛物线联合来计算由于基坑开挖和基坑底部隆起2个因素引起的坑外地表沉降,如图 3、4。
4)考虑到由于基坑底部隆起,导致基坑开挖引起的沉降最大值出现在距支护结构一定距离的xo处(δv(xo, 0)=δv max),故横向沉降曲线的假设运用
5)大量学者的研究表明[13-17],在基坑支护边缘处的土体沉降并不为0,而是一个与最大沉降值δv max具有一定比例关系的一个范围值;
6)基坑具有很明显的空间效应,实际工程中支护结构的变形会出现,在支护结构中部变形最大,在基坑角部则变形比较微小,这就是边角效应。考虑到边角效应,在计算上述式子的未知系数a、b时,假设一个边界条件,在基坑角部位于支护结构处土体沉降值为零,即
根据以上边界条件即可求得上述沉降式(10)的未知系数a、b
把以上系数代入式(10)即得围护结构背后地表沉降曲线的方程。
基坑降水、支护结构变形及坑底隆起引起的沉降之间密不可分,基坑内外的水位压力差也会引起支护结构变形、基坑底部隆起,如果单独考虑各因素,把最后的沉降叠加则可能会引起沉降量叠加,故应对沉降和进行修正。因降水对支护结构变形和基坑隆起的影响主要是通过水位压力差引起的,故采用修正系数η
则计算基坑开挖引起土体沉降的理论公式为
某工程开挖深度10.35~12.70 m,基坑为长方形,开挖面积约55 m×35 m,根据“安全、经济、方便施工”的原则,基坑采用排桩(钻孔灌注桩)+两道钢筋砼对撑及角撑(外围三轴水泥搅拌桩止水)坑内外简易管井降水(坑外控制性降水)的围护形式。场地地下水类型上部属孔隙潜水,下部碎石混粘土中含孔隙承压水,实测潜水位埋深0.70~1.40 m,坑内地下水位控制在坑底以下0.5 m,基坑周边实施控制性降水,坑外水位控制在4 m以下。基坑土质条件一般,开挖深度影响范围内,基坑上部含有5~6 m的砂质粉土。
基坑坑外地下水位控制在4 m以下,根据式(13)可计算出修正系数η=0.747,第i层土范围内平均附加应力系数αi根据《建筑地基基础设计规范》附录K矩形面积上均布荷载作用下角点的平均附加应力系数α取用。根据本基坑工程的支护条件可取最大沉降值δvmax为0.1%H,则根据以上工程条件可以计算出最终的理论沉降曲线并与实测值进行对比分析,见表 2。
如图 5所示,根据计算结果可以分析得出以下结论:
1)基坑降水在支护结构边缘引起最大沉降,由基坑开挖引起的最大沉降发生在距支护结构边缘一定距离处。
2)两者沉降和修正后产生的最大沉降在本工程中发生在基坑边缘,且产生的坑后沉降曲线近似三角。
3)工程的坑外降水深度比较浅,故坑外土体沉降主要受基坑开挖影响比较大。
4)实测数据与理论计算值非常接近,可以有效预测坑外土体沉降量。
分析了由降水引起的土体沉降的机理,采用修正的分层总和法计算由降水引起的土体沉降,合理地推算出由降水引起的坑外土体沉降计算公式。把支护结构变形和坑底隆起二者结合起来考虑,合理假设由二者引起的坑外土体沉降形式,推导得出基坑开挖引起的坑外土体沉降理论计算公式。把降水沉降及基坑施工沉降二者计算结果叠加并加以修正,最终得到坑外土体沉降的理论解析解。经过实际工程验证,实测值分布在理论值范围内,可以有效地预估坑外土体沉降量。在预估基坑外土体沉降的研究中,尽可能多地考虑引起沉降的重要因素以及引起沉降的机理,把这些因素综合起来考虑合理假设,再通过实际工程进行验证、修正、验证。如此反复,最终可以得到合理有效的理论结果,推动基坑方面的理论发展,指导工程实践,使得理论研究具有重大现实意义。