丁士龍,孫志浩,張恒志,竇炳珺,鄭凱達,徐長節,4,5
(1.浙江省大成建設集團有限公司,浙江 杭州 310000;2.浙江大學建筑工程學院,浙江 杭州 310058;3.浙江大學平衡建筑研究中心,浙江 杭州 310058;4.華東交通大學江西省巖土工程基礎設施安全與控制重點實驗室,江西 南昌 330013;5.華東交通大學江西省地下空間技術開發工程研究中心,江西 南昌 330013)
基坑開挖會引起周圍土層沉降[1-3],進一步引發周圍建筑沉降[4-5]。 基坑開挖引起的地表沉降曲線形式一般為凹槽型或三角形,這兩種形式的地表沉降均會造成周圍建筑的不均勻沉降進而造成建筑物墻體產生裂縫甚至破壞。 目前已有不少學者對基坑開挖引起周圍建筑物的影響進行了研究。
在理論解析方面,木林隆等[6]利用考慮土體小應變特性的簡化計算方法計算了基坑開挖引起的土體位移場,再結合被動樁兩階段分析法,提出了一種計算基坑開挖對鄰近樁基礎影響的簡化分析方法;趙延林[7]應用彈塑性大變形理論研究了樁錨支護基坑開挖引起的鄰近建筑物不均勻沉降,并提出了基坑開挖存在一個臨界深度,當開挖超過臨界深度時, 建筑物的不均勻沉降變化率會急劇增加。在監測數據分析方面,劉念武等[8-9]根據監測數據詳細分析了基坑開挖對淺基礎建筑及樁基礎建筑沉降的影響;孟憲國等[10]利用Verhulst 分析方法預測了地鐵端頭井基坑開挖引起的地表沉降,并與實測數據進行了對比,證明了該方法預測地表沉降的有效性;王池[11]結合上海某風井基坑地表沉降監測值,對不同地表沉降預測模型進行了對比分析,并比選出了精度最高的地表預測模型;Cheng 等[12]根據實測數據分析了某逆作法基坑開挖對周邊多個既有建筑物的影響,并認為建筑物基礎形式、尺寸及與基坑的距離都會影響鄰近基坑建筑物的變形。 在數值模擬方面,高丙麗等[13]和崔鐵軍等[14]利用FLAC有限差分軟件對盾構工作井基坑開挖引起的周邊土體及管線變形規律進行了研究;章潤紅等[15-16]利用PLAXIS 2D 軟件對基坑開挖卸載作用下臨近地鐵結構附加彎矩、 位移響應和臨近邊坡安全性、位移響應進行了研究;蘆友明[17]和劉睿[18]利用ABAQUS有限元分析軟件對基坑開挖引起鄰近建筑物的變形規律進行了詳細的分析,并綜合考慮了基坑與建筑物距離及基坑支護形式等因素的影響;為精細化分析建筑物與基坑間空間位置差異對建筑物變形產生的影響,鄭剛等[19]利用有限元軟件PLAXIS 分析了建筑物與基坑邊成不同距離及角度時建筑物的變形性狀。
以往的研究主要基于平面應變狀態下基坑開挖對臨近建筑物影響,或者僅考慮三維條件下基坑開挖對非坑角區域建筑物的影響。 因存在坑角效應,即坑角處支護結構的受力變形較小,學者較少對該區域建筑物的響應進行研究。 但實際工程監測中發現,當基坑開挖較深時,對坑角外建筑物(圖1)也會產生較大影響。 為了研究坑角外建筑物對深基坑開挖的響應問題,以杭州市某地鐵風井深基坑為研究對象,采用有限元軟件PLAXIS 3D 對風井基坑及周邊淺基礎建筑進行了三維數值模擬。 通過數值模擬結果與現場監測數據對比驗證模型的有效性。在此模型基礎上, 對淺基礎建筑物形狀進行了簡化,并通過改變淺基礎建筑與坑角距離,研究了不同距離下坑角外建筑物墻體沉降、傾斜、扭轉和基礎水平位移的變化規律。 得出的結論可為類似工程監測方案的制定提供參考。

圖1 坑角外區域示意圖Fig.1 Schematic diagram of the area outside the corner of the foundation pit
某城際鐵路中間風井基坑位于浙江省杭州市。風井基坑開挖深度為27 m,采用咬合樁加3 道內支撐的支護形式,咬合樁樁徑為1.0 m,樁間距為0.8 m,樁長為30 m。 第1、3 道支撐為鋼筋混凝土支撐,第2道支撐為鋼支撐,基坑支護結構剖面示意圖如圖2 所示。 4 道支撐分別設置于-1.05,-8.05,-15.05,23.00 m處,基坑支撐平面布置圖如圖3 所示,其中鋼支撐安裝時施加990 kN 的預加軸力。 風井基坑周圍有多幢磚混結構建筑(均為淺基礎),南側緊鄰河道,風井基坑平面位置示意圖如圖4 所示。 該場地土層上軟下硬,上部為素填土層,下部為卵石層和巖石層。 現場巖土層的基本物理力學參數見表1。

圖2 基坑支護結構示意圖Fig.2 Schematic diagram of foundation pit support structure

圖3 基坑支撐布置圖(單位:mm)Fig.3 Layout of foundation pit support(Unit:mm)

圖4 風井基坑平面位置圖Fig.4 Plan view of foundation pit of subway wind tunnel

表1 現場巖土層的基本物理力學參數Tab.1 Basic physical and mechanical parameters of field rock and soil layers
使用有限元軟件PLAXIS 3D 對上述基坑及周圍建筑進行模擬。 土體和全(強)風化巖采用小應變土體硬化模型 (hardening soil model with smallstrain stiffness ,HSS 模型),該模型能適用于多種土類(軟土和較硬土層)的破壞和變形行為的描述。 中外已有許多學者[15,20-22]使用HSS 模型來模擬基坑開挖時的土體變形,其模擬結果與實測結果具有很好的一致性。
HSS 模型是土體硬化模型(hardening soil model,HS 模型)的一個改進,兩者均可以描述土體的剪切硬化和體積硬化行為,其區別在于小應變條件下剛度值的選用。 HS 模型假設土體在卸載和重加載時的剪切剛度是定值,但實際上隨著土體應變的增大,土體的剪切剛度G 會呈非線性衰減(圖5)。HSS模型通過增加兩個參數G0,γ0.7來模擬土體的這種特性,其中G0為土體在應變很小時的剪切模量,γ0.7為割線剪切模量退化為0.7G0時的剪切應變。

圖5 土體的剪切剛度-應變曲線Fig.5 Shear modulus-strain curve of soil
中風化巖層與土體在性質上有較大差異,用一般土體的本構模型進行模擬較難得到合理結果,故本文采用可以描述巖石強度與變形行為的霍克-布朗模型進行模擬。 霍克-布朗模型所采用的破壞準則是一種非線性強度近似準則,在其連續方程中不僅包含剪切強度,也包含拉伸強度[23]。該破壞準則可以用最大主應力和最小主應力關系來描述

式中:σ1為最大主應力;σ3為最小主應力;mb是對完整巖石參數mι的折減;σci為完整巖石的單軸抗壓強度;s 和α 均為巖塊的輔助材料參數。
各巖土層的參數取值見表1 和表2, 巖土層參數主要由地勘報告獲得, 地勘報告中未給出的參數根據相關文獻[23-24]和參數反分析獲得。 根據文獻[24], 關于小應變本構模型各參數的經驗取值方法如表3 所示。 使用板單元模擬咬合樁,界面單元模擬樁土相互作用(界面單元強度Rinter 取0.9,Rinter 定義見文獻[23]),梁單元模擬內支撐和圍檁,板單元模擬周邊建筑物的墻體、樓板和基礎,使用面荷載模擬坑邊施工荷載,面荷載根據實際工況取30 kPa。 各結構計算參數見表4,有限元模型如圖6所示,有限元網格劃分如圖7 所示。

圖6 有限元模型Fig.6 Finite element model

圖7 數值分析模型網格圖Fig.7 Mesh of model for numerical simulation

表2 模型中巖土層的物理力學參數Tab.2 Physical and mechanical parameters of the rock and soil layer in the model

表3 HSS 模型參數Tab.3 HSS model parameters

表4 結構計算參數Tab.4 Calculation parameters of structure
為考慮建模的尺寸效應, 進行了模型尺寸試算,最終確定X 和Y 方向計算寬度自開挖邊界向外取開挖深度的5 倍,Z 方向計算寬度自開挖邊界向外取開挖深度的4 倍。 模型邊界設置為標準邊界,即側邊約束水平向位移,豎向位移自由,底邊同時約束水平和豎向位移,頂面無約束。
模型嚴格模擬實際施工工況, 基坑開挖的具體實現步驟為:①建立整個場地土體、建筑物及支護結構模型;②初始地應力的平衡,建立初始應力場,激活建筑物,同時支護結構的剛度消失,即支護結構單元失去活性,使土體自重沉降過程中支護結構對土體自重沉降無影響;③初始應力場引起的位移值清零, 激活支護結構板單元,激活坑邊荷載; ④降水至基坑開挖面以下1 m,分層挖土并激活相應支撐(若為鋼支撐則同時激活預應力),降水是通過使土體變為干實現的,土體開挖是通過逐層使土體單元失去活性實現的[22]。主要開挖工況如表5 所示。

表5 開挖階段工況Tab.5 Working conditions at the excavation stage
選取剖面A-A(圖7)與基坑ab 邊交點處圍護樁為研究對象,分別繪制出開挖至坑底時該處圍護樁水平位移計算值與實測值隨深度變化的曲線。 如圖8 所示,圍護樁水平位移計算值與實測值變化趨勢一致,隨樁深增加均呈現出先增大后減小,然后緩慢增大再減小的趨勢, 且均在樁深為10 m 左右時達到最大值。 由于開挖采用自上而下的方式,支撐架設及時,支護樁最大水平位移并沒有出現在最終開挖標高附近[25]。

圖8 支護樁水平位移計算值與實測值對比Fig.8 Comparison of calculated horizontal displacement of supporting pile with measured value
圍護樁水平位移實測最大值為16.79 mm,模擬最大值為14.85 mm。 深度0 m 處模擬值大于實測值,主要原因是實際工程中對基坑周邊地表進行了硬化,硬化路面與第一道支撐的冠梁部分剛接從而抑制了樁頂的位移。 深度10 m 處模擬值小于實測值, 主要原因在于實際工程中的坑邊荷載較為復雜, 而模擬過程中僅采用了均布面荷載來模擬,導致不能完全反映坑邊荷載對支護結構位移的影響。除此之外, 土的各向異性對開挖分析有一定影響,但HSS 模型不能考慮這種特性,這也是產生誤差的原因之一[26]。 總體來說模擬曲線與實測曲線較為吻合,證明本文的有限元模型參數選取合理,可用于后續研究。
為研究坑角外建筑物與坑角距離對建筑物墻體和基礎的影響,在上述模型的基礎上,改變建筑物與基坑的距離,設置了建筑物與基坑角部凈距分別為0.2h,0.5h,1.0h,2.0h(h 為基坑開挖深度)共4組模型,如圖9 所示。 為了使分析結果更具有普遍意義,對坑邊建筑物進行了適當的簡化,將長條形的淺基礎建筑調整為等面積的正方形建筑,并設置建筑物位于坑角45°方向,如圖9 所示。

圖9 建筑物距基坑不同距離示意圖Fig.9 Diagram of different distance between building and foundation pit
提取不同距離下臨近建筑物靠近基坑側墻體(圖9)沉降值,如圖10 所示。

圖10 不同距離下臨近建筑物墻體(靠近基坑側)沉降Fig.10 Settlement of wall near building ( near foundation pit ) at different distance
由圖10 可知, 當臨近建筑物與基坑距離逐漸增大時, 其墻體最大沉降值呈先增大后減小的趨勢。 主要原因在于基坑開挖時,產生的坑外地表沉降曲線為凹槽形,其上建筑物也會產生與地表沉降相同的沉降趨勢。 當建筑物與坑邊較近時,建筑物墻體左側沉降小,右側沉降大,而當建筑物與坑邊較遠時,建筑物墻體左側沉降大,右側沉降小。其與基坑距離的增大呈先增大后減小再反向增大的趨勢。 當距離為1.5h 時,最大差異沉降值為-0.13 mm,表明此時墻體靠近基坑側沉降大,遠離基坑側沉降小。

表6 不同距離下臨近建筑物墻體最大差異沉降Tab.6 Maximum differential settlement of walls of adjacent buildings at different distance
提取不同距離下坑角外建筑物靠近基坑側墻體位置A 處水平位移值(圖9),如圖11 所示。 建筑物距離基坑較近時, 墻體向遠離基坑方向傾斜,而當建筑物離坑邊超過一定距離后,墻體向基坑方向傾斜。 該結果表明,受基坑開挖形成的沉降槽影響,與坑角距離較近處建筑物有“后仰”變形趨勢,而與坑角距離較遠處建筑物有“前傾”變形趨勢。

圖11 不同距離下臨近建筑物墻體(A 點)傾斜Fig.11 Inclination of wall ( point A ) of adjacent buildings at different distance
定義墻體整體傾斜為墻體傾斜方向兩端點的水平方向位移差與墻體高度的比值,表7 為不同距離下臨近建筑物墻體整體傾斜對比。 由表7 可知,隨著建筑物離基坑距離的增大,建筑物的傾斜呈現先增大后減小的趨勢,最后會出現反向傾斜,建筑物墻體傾斜與墻體的差異沉降變化趨勢是一致的。

表7 不同距離下臨近建筑物墻體整體傾斜對比Tab.7 Overall inclination comparison of walls of adjacent buildings at different distance
提取不同距離下坑角外建筑物靠近基坑側墻體位置B 處水平位移值,如圖12 所示。 由圖12 可知,受Y 方向地表沉降變化影響,位置B 處的墻體水平位移值與A 處存在一定差異,這將造成臨近基坑側墻體發生扭轉變形。 定義扭轉程度=|A 處墻體整體傾斜-B 處墻體整體傾斜|/A 處墻體整體傾斜,絕對值越大則墻體的扭轉程度越大,不同距離下臨近建筑物的墻體扭轉程度如表8 所示。 從表8 中可以看出,隨著建筑物與坑邊距離的增大,建筑物墻體的扭轉程度呈先減小后增大再減小的趨勢。 三維空間條件下,建筑物受基坑開挖產生的墻體扭轉程度變化趨勢較為復雜,主要原因在于建筑物墻體傾斜同時受到X,Y 兩個方向開挖卸荷引起的地表沉降影響。

表8 不同距離下臨近建筑物墻體扭轉程度對比Tab.8 Torsion degree comparison of walls of adjacent buildings at different distance

圖12 不同距離下臨近建筑物墻體(B 點)傾斜Fig.12 Inclination of wall (point B) of adjacent buildings at different distance
提取不同距離下建筑物靠近基坑側基礎水平方向位移值(XY 平面內位移),如圖13 所示。

圖13 不同距離下臨近建筑物基礎水平方向位移Fig.13 Horizontal displacement of foundation of adjacent buildings at different distance
由圖13 可知, 基坑開挖會引起臨近建筑物基礎向基坑反方向位移,但隨著建筑物與基坑距離的增大, 建筑物基礎的水平位移呈逐漸減小趨勢,但當距離增大到1.5h 時,建筑物基礎開始朝向基坑側位移。 如3.2 節分析,當建筑物與坑角距離較近時,建筑物將發生“后仰”的變形趨勢,此時建筑物基礎也將發生遠離基坑側的位移;而當建筑物與坑角距離較遠時,建筑物將發生“前傾”的變形趨勢,此時建筑物基礎會發生向基坑側靠近的位移。
由圖13 可知, 墻體的水平位移曲線沿建筑物基礎長度方向基本沒有變化,即沿X 方向建筑物基礎的水平位移變化較小,主要原因在于基礎的剛度較大導致其呈整體位移形式。
采用三維數值模擬方法,研究了風井深基坑開挖對坑角外臨近淺基礎建筑物的影響,主要得出以下結論:
1) 建筑物與基坑距離較近時,墻體沉降呈臨近基坑側小,遠離基坑側大的趨勢,但當建筑物與基坑超過一定距離(本算例約為1.5h),墻體沉降會呈相反趨勢。 墻體不均勻沉降呈先增大后減小再反向增大趨勢。
2) 建筑物與基坑距離較近時, 建筑物墻體向基坑外傾斜,但當距離增大到一定值時(本算例為1.5h),建筑物墻體向基坑內傾斜,墻體傾斜值亦呈先增大后減小再反向增大的趨勢。
3) 因基坑開挖導致的坑角外建筑物墻體兩端傾斜程度不一致,墻體會產生扭轉變形。 隨著建筑物與坑邊距離的增大,建筑物墻體的扭轉程度呈先減小后增大再減小的趨勢。
4) 建筑物與基坑距離較近時,建筑物基礎向遠離基坑側水平位移, 隨著建筑物與基坑距離的增大,建筑物基礎的水平位移將逐漸減小,最終會發生向基坑方向靠近的位移。