盾構隧道開挖引起的鄰近群樁豎向位移研究

可文海，管凌霄，薛齊，徐長節,4，耿大新

(1.華東交通大學江西省巖土工程基礎設施安全與控制重點實驗室，江西南昌，330013；2.江西省地下空間技術開發工程研究中心，江西南昌，330013；3.浙江杭海城際鐵路有限公司，浙江嘉興，314499；4.浙江大學濱海和城市巖土工程研究中心，浙江杭州，310058)

眾所周知，盾構隧道在開挖過程中會引起地層損失，當隧道掘進方向存在鄰近樁基時，必然會導致鄰近樁基產生變形甚至破壞，進而威脅到上部結構安全。因此，準確地預測盾構隧道開挖對鄰近群樁的影響，成為城市地鐵工程中重點關注的問題之一[1]。

從20世紀70年代開始，眾多的國內外學者研究了上述問題[2]。在試驗研究方面，MORTON 等[3]利用室內模型試驗發現盾構隧道開挖對鄰近樁基的影響較大，尤其是在軟土中；ONG等[4]采用離心機試驗研究了黏性土中的盾構隧道開挖對鄰近樁基的影響；在數值模擬方面，CHEN 等[5-6]采用數值方法，將盾構隧道開挖引起的土體位移施加到樁基上模擬樁-土相互作用，進而分析盾構隧道開挖引起的樁基內力與位移；在理論分析方面，兩段法是目前研究盾構隧道開挖引起鄰近樁基變形的常用方法[7-8]，即將盾構隧道開挖引起鄰近樁基變形分為2個階段：第1階段分析盾構隧道開挖引起的樁基軸線處土體位移；第2階段將樁基軸線處土體位移直接加到樁基上，再采用合理的模型模擬樁-土相互作用，求解出樁基變形的解析解。兩段法的計算步驟清晰簡潔，物理意義明確[9-10]，能準確計算盾構隧道開挖引起的樁基變形。第1階段中，LOGANATHAN等[11]提出盾構隧道引起周圍土體位移的計算公式，該公式具有良好的計算精度而被廣泛應用[12]；第2 階段中，程康等[13]將樁基視為Pasternak 地基上的鐵木辛柯梁，研究了盾構隧道開挖引起的樁基水平變形；張治國等[14]考慮了樁側土體三維效應對樁基的影響，分析了樁基的水平變形及彎矩；李早等[15-16]基于Winkler 地基梁模型并結合有限差分法，研究了盾構隧道開挖引起的鄰近樁基豎向位移及內力。

綜上所述，眾多學者全面分析了盾構隧道開挖引起鄰近樁基受力變形的問題，但上述的理論分析中，第2階段基于彈性地基梁的求解方法需要結合有限差分法[14-16]，甚至需要擬合盾構隧道開挖引起的土體位移函數，這些求解步驟大幅增加了計算的復雜程度，無法簡便地針對盾構隧道開挖引起的樁基受力變形進行全面地參數分析。

為此，在前人研究的基礎上，采用兩段法在第2階段將盾構隧道開挖引起的土體豎向位移轉化為荷載施加到樁基上，再基于Winkler地基梁模型結合疊加法原理[17]，首先求出集中荷載作用下的樁基位移，其次進行積分得到樁基在分布荷載作用下的樁基位移，最后考慮群樁間的土體遮攔效應[18]求解出盾構隧道開挖引起的鄰近群樁豎向位移。此外，將本文計算結果與有限元模擬結果進行對比驗證，并深入分析了隧道埋深、地層損失比、樁基直徑、土體彈性模量以及群樁間距變化對樁基豎向位移的影響。

1 基本方程的建立與求解

1.1 盾構隧道開挖引起的土體自由場豎向位移

LOGANATHAN等[11]基于線彈性理論提出了盾構隧道開挖引起鄰近土體豎向位移的解析解，其計算結果與實測值吻合較好。故本文根據LOGANATHAN 公式，計算盾構隧道開挖引起的鄰近土體豎向位移U(x,z)[11]

式中：R為隧道半徑；x為距隧道中心線的水平距離；z為距地表的垂直距離；ε0為隧道開挖造成的地層損失比；H為隧道軸線埋深；v為土體泊松比。

1.2 盾構隧道開挖對鄰近單樁的影響

圖1所示為本文采用的簡化計算模型，其中，d為樁基直徑，L為樁基的入土深度，x0為樁軸線到隧道中心的距離，此時盾構隧道開挖引起的樁基軸線處(x=x0)處的土體豎向位移為U(z)=U(x0,z)。在盾構隧道開挖引起的樁-土相互作用方面，基于Winkler理論作如下假定：

1)樁為彈性體；

2)土為彈塑性連續變形體；

3)樁-土之間不發生滑移，且始終保持彈性接觸。

圖1 簡化計算模型Fig.1 Simplified calculation model

基于Winkler 模型，樁基豎向位移的控制方程為[19]

式中：W(z)為樁基的位移函數；Ep為樁基彈性模量；Ap為樁基橫截面面積；P(z)=kU(z)為樁基受到的豎向外荷載；k為樁側土的彈簧剛度，根據RANDOLPH等[18]的方法，彈簧剛度k可?。?/p>

式中：Gs為土體剪切模量；r0=d/2 為樁基半徑；rm為樁的有效影響半徑，rm=χL(1-v)；χ為經驗系數，對于均質土，χ≈2.5[20]。

為求解式(2)中樁基受分布力作用下的位移函數，可先求解受集中力作用的樁基控制方程，圖2所示為集中荷載作用下樁基模型圖。即令P(z)=0，求其通解，得

式中：A1和A2為待定系數；

再假設此時樁基底端反力為0，結合式(4)，可得樁基軸線任意點ξ受集中力dp(見圖2)引起的樁基豎向位移dW*(z)為[21]

圖2 集中荷載作用下的樁基模型Fig.2 Pile foundation under concentrated load

在Winkler 地基模型中，任意點ξ因盾構隧道開挖引起的土體位移產生的集中力為

通過Winkler的疊加原理，對式(5)在樁基入土深度范圍內積分，可得到不考慮樁基底端反力時盾構隧道開挖引起的樁基豎向位移W*(z)

由于實際的樁基底端反力為Pb=kb[W(L)-U(L)]，其中：kb為樁基底端的彈簧剛度；kb=Es為樁基底端土體的彈性模量；hb為樁基底端到基巖的深度；U(L)為樁基底端土體的原位移；W(L)為樁基底端實際豎向位移。根據樁基底部樁-土相互作用的平衡關系[21]，可得：

解得樁基底端反力后，將Pb看作樁基在ξ=L處受到的方向向上的集中力，根據Winkler疊加原理便可求得實際的樁基豎向位移

1.3 盾構隧道開挖對鄰近群樁的影響

樁基對鄰近土體的位移會產生阻礙作用，當盾構隧道開挖引起鄰近土體自由場沉降時，樁位處的土體位移會因鄰近樁基的阻礙作用而減小，這種現象稱為樁基的遮攔效應。因此，土體自由場中的樁基不僅會受到盾構隧道開挖的影響，還會受到鄰近樁基的影響。

圖3所示為群樁簡化計算模型，按照距離隧道軸線的遠近將樁基依次定義為樁1、樁2，…，樁i，…，樁n。若土中不存在樁基，盾構隧道開挖引起的樁1 位置處土體豎向位移為U1(z)，樁1 因盾構隧道開挖引起的豎向位移為W1(z)，則樁1 對土體自由場位移的遮攔效應為ΔW1(z)=U1(z)-W1(z)。根據COOKE 提出的剪切位移法，將樁身周圍土體理想地視作同心圓柱體，設樁間距為s，采用RANDOLPH 等[18]提出的樁間土位移傳遞系數ψ(s)，可得到樁2 位置處因樁1 遮攔效應而產生的土體位移U21(z)[19]

圖3 群樁簡化計算模型Fig.3 Simplified computing model of group pile

由于樁基遮攔效應引起的位移與土體自由場位移方向相反，因此，在盾構隧道與樁1的共同影響下樁2位置處的土體位移為

同理可得，樁1位置處因盾構隧道、樁2共同影響而引起的土體位移為(z)=U1(z)-U12(z)。

基于此，本文計算盾構隧道開挖引起鄰近群樁豎向位移的具體步驟為：

1)采用式(1)計算出盾構隧道開挖引起各樁基位置處的土體豎向位移Ui(z)(i=1，2，…，n)；

2)計算出只考慮單樁時，各樁基在盾構隧道開挖影響下產生的豎向位移Wi(z)；

3)根據式(10)和(11)計算出考慮鄰近樁基的遮攔效應時，樁i軸線處因盾構隧道與鄰近樁基共同影響下的土體豎向位移i，i=1，2，…，n)；

2 算例驗證

為了驗證本文所提方法的準確性，以LOGANATHAN 等[22]應用三維邊界元程序GEPAN模擬隧道開挖對鄰近2×2群樁影響。該計算假設土體為均質且各向同性的彈性體，土體彈性模量Es為24 MPa，泊松比v為0.5，樁基也為各向同性彈性體，樁基彈性模量Ep為3×104MPa，入土深度L為25 m，直徑d為0.8 m，隧道半徑R為3 m，隧道軸線埋深H為20 m，地層損失率ε0為1%。第1 排樁距離隧道軸線為4.5 m，第2 排樁距離隧道為6.9 m，兩樁軸線間距為2.4 m。

圖4所示為本文方法計算結果與GEPAN 程序計算結果的比較。其中圖4(a)所示為只考慮單樁時，盾構隧道開挖引起鄰近單樁沉降的計算結果。由圖4(a)可知：采用本文方法計算的單樁豎向位移曲線與GEPAN 程序計算結果十分吻合，所呈現的規律基本一致，證明了本文方法的準確性。圖4(b)所示為考慮群樁的遮攔效應時，盾構隧道開挖引起鄰近群樁沉降的計算結果。由圖4(b)可知：對于后樁(距隧道軸線更遠)，采用本文方法計算的樁基豎向位移曲線與GEPAN 程序計算結果基本吻合；但是對于前樁(距隧道軸線更近)，本文方法得到的結果整體偏大。這是由于本文方法在分析群樁之間的相互影響時只考慮到鄰近樁基對土體自由場的遮攔效應，而實際工程中應該還存在其他效應。當分析前樁對后樁的影響時，前樁引起的土體遮攔效應較大，為群樁之間相互影響的主要影響因素；但分析后樁對前樁的影響時，后樁引起的土體遮攔效應較小，不是主要影響因素。

3 豎向位移影響因素

本文所提計算方法無需擬合土體位移函數，也無需結合有限差分法，計算過程簡便，更易于研究樁基豎向位移與各物理參量之間關系?，F取如下算例進行分析，隧道的計算參數如下：H為20 m，R為3 m，?0為1%；樁基的計算參數如下：d為1 m，L為25 m，Ep為3×104MPa；取土中存在2根樁，分析盾構隧道對群樁的影響，前樁軸線到隧道的距離x1為5 m，后樁軸線到隧道的距離x2為8 m；土體的計算參數如下：Es為20 MPa，泊松比v為0.3。本文在研究某一參數對樁基豎向位移的影響時，其余參數不變。

圖4 隧道附近2×2承臺群樁計算結果比較Fig.4 Comparison of results of 2×2 pile group with cap adjacent to tunneling

3.1 隧道埋深

為研究隧道埋深變化對群樁豎向位移的影響，分別分析隧道埋深為10，15和20 m的情況。

圖5所示為不同影響因素變化下的樁基豎向位移。由圖5(a)可見：隧道埋深從10 m 增加到20 m時，兩樁的豎向位移均增大，但增幅在減小。這是由于盾構隧道開挖后，隧道中心線上方的土體向下位移，中心線下方的土體向上位移，因此，隧道埋深增加將引起了土體向下位移的范圍增大，此時，土中樁基便受到更多的土體向下位移的荷載，導致其產生了更大的豎向位移。

3.2 地層損失比

為研究地層損失比變化對群樁豎向位移的影響，取地層損失比分別為0.5%，1.0%和2.0%進行分析。從圖5(b)可見：樁基的豎向位移隨著地層損失比增加而等比例增加。這是因為式(1)中，地層損失比與隧道開挖引起的土體位移為正比例關系，當地層損失比增大時，隧道開挖引起的土體位移將隨之等比例增加，進而導致鄰近群樁受盾構隧道開挖影響的程度等比例增大。

3.3 樁基直徑

分別取樁基直徑為0.50，0.75和1.00 m分析樁基直徑變化對群樁豎向位移的影響。從圖5(c)可見：樁基直徑從0.50 m 增加到1.00 m 時，兩樁的豎向位移均有略微減小。這是由于樁基直徑增加會導致其剛度增大，進而引起其抵抗樁周土體位移影響的能力增加。因此，在同樣的土體位移荷載作用下，樁基直徑增加可引起其產生的豎向位移減小。

3.4 土體彈性模量

分別取土體的彈性模量分別為10，20 和30 MPa，研究土體彈性模量變化對群樁豎向位移的影響。從圖5(d)可見：隨著土體彈性模量增加，兩根樁基的頂端豎向位移(最大位移)均增大，但底端豎向位移均減小。結合3.1節分析認為，土體彈性模量增加引起隧道中心線上方的樁體受到土體向下位移的荷載增大，隧道中心線下方的樁體受到土體向上位移的荷載增大，同時樁基底端受到土體的反力也隨之增大，這導致了樁頂位移增大而樁底位移減小的現象。

圖5 不同影響因素變化下的樁基豎向位移Fig.5 Vertical displacement of pile foundation under different influencing factors

3.5 樁基間距

為研究樁基間距變化對群樁豎向位移的影響，取表1所列工況進行分析。

表1 不同樁基間距工況Table 1 Working conditions of different pile foundation spacing m

圖6 樁基間距離變化對樁基豎向位移的影響Fig.6 Influence of different pile spacing on vertical displacement

圖6所示為不同樁基間距下的樁基豎向位移。從圖6(a)可見：隨著后樁的逐漸遠離，前樁的豎向位移增大，且逐漸接近單樁的豎向位移。這是由于后樁與前樁的距離增加時，其對前樁造成的遮攔效應減弱。從圖6(b)可見：隨著后樁到隧道的距離越來越遠，后樁的豎向位移逐漸減少，且與單樁時豎向位移的差值越來越小。這是因為后樁在遠離隧道時，其受到盾構隧道開挖的影響減小，同時其與前樁的間距在增加，因此，受到前樁遮攔效應的影響也逐漸減弱。由此可見：樁間距增加可導致群樁間的遮攔效應減弱，樁基的相對豎向位移增加；樁基與隧道距離增加可導致樁基受到盾構隧道開挖的影響減小、樁基的豎向位移減小。

4 結論

1)采用本文方法計算盾構隧道引起的鄰近單樁豎向位移與有限元結果吻合，但在考慮群樁的遮攔效應時，計算的前樁豎向位移與有限元結果相比整體偏大。

2)控制其余參數不變時，隧道埋深增加及隧道開挖引起的地層損失比增大均會導致鄰近群樁的豎向位移增大。

3)當樁基直徑增大時，鄰近群樁因盾構隧道開挖引起的豎向位移均略微減小。

4)土體彈性模量增加會導致鄰近群樁的頂端豎向位移(最大位移)增大，底端豎向位移減小。

5)樁基與隧道距離增加可減弱盾構隧道開挖對鄰近樁基的影響，減小樁基豎向位移；群樁間距增大導致樁基間的遮攔效應減弱、樁基的相應豎向位移增大。