劉小兵，陳有亮，孟偉波

（上海理工大學(xué) 環(huán)境與建筑學(xué)院，上海 200093）

隨著城市的發(fā)展，高層建筑、超高層建筑、地下建筑等工程大量興起，而與之對(duì)應(yīng)的是越來越多、越來越深的深基坑工程的出現(xiàn)。在城市復(fù)雜環(huán)境條件下基坑開挖引起的地下管線變形直接影響到人民的生命財(cái)產(chǎn)安全，因此受到極大的重視。所謂基坑開挖的環(huán)境效應(yīng)問題就是基坑開挖使土體由初始應(yīng)力場(chǎng)狀態(tài)變?yōu)榈诙?yīng)力狀態(tài)。土體內(nèi)應(yīng)力重新分布，致使圍護(hù)結(jié)構(gòu)產(chǎn)生變形、位移，引起基坑周圍地表沉陷，從而對(duì)臨近建筑物和地下設(shè)施帶來不利影響，諸如臨近建筑物的開裂、傾斜、道路開裂、地下管線的變形和開裂等[1]。在基坑開挖對(duì)周圍地下管線影響的評(píng)估中，分析基坑引起的管線變形是非常重要的工作。

在管線變形的眾多影響因素中，深基坑圍護(hù)結(jié)構(gòu)變形和其基坑周圍地表沉陷，是引起管線位移、變形，甚至破壞的根本原因。目前國(guó)內(nèi)外很多學(xué)者對(duì)基坑開挖引起周圍地下管線的位移進(jìn)行了研究[2-9]，并得出了不同的結(jié)論。

為此，本文針對(duì)深基坑開挖引起的周圍地下管線位移，結(jié)合Winkler 彈性地基梁模型，采用兩階段法，建立了受土體卸載附加變形影響的地下管線豎向和水平方向位移方程，通過有限單元法解出地下管線位移，經(jīng)過與工程實(shí)例對(duì)比分析，表明理論計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)吻合。

1 基坑開挖引起周圍地下管線位移的兩階段法分析

兩階段法將基坑開挖對(duì)管線的影響分成兩個(gè)階段。第一階段，分析開挖引起的管線位置處的自由土體位移，包括豎向位移和水平位移；第二階段，基于Winkler 地基模型，將上述自由土體位移視為作用在管線上的附加應(yīng)力，分析計(jì)算管線的附加位移。

1.1 作用在管線上的豎向附加應(yīng)力

根據(jù)李大勇等[1]的研究，地下管線受到自重、上覆土壓力、地面超載、地基反力等作用力。但由于自重和上覆土壓力在基坑開挖前就已經(jīng)達(dá)到平衡，因此在計(jì)算基坑開挖引起的管線變形時(shí)，應(yīng)將其忽略。李大勇等[1]將基坑周圍地表分為沉陷區(qū)和非沉陷區(qū)分別計(jì)算，并將沉陷區(qū)的沉降曲線設(shè)為拋物線，并假定沉陷區(qū)和非沉陷區(qū)交界處土體表面不發(fā)生沉降，這與實(shí)際情況差距較大，易產(chǎn)生較大誤差。

張陳蓉等[10]在前人研究的基礎(chǔ)上提出了坑外地表土體沉降預(yù)測(cè)公式

式中：ωv,max為地面最大沉降值；H 為基坑開挖深度；A 為變形影響半徑，與基坑開挖深度H 和基坑開挖深度長(zhǎng)度L 有關(guān)。

基于位移控制理論，張陳蓉等[10]由擬合結(jié)果得到上海軟土地區(qū)板式維護(hù)基坑體系基坑外的任一點(diǎn)土體沉降曲線公式

因此，地表沉降引起的管線處附加應(yīng)力q1為

式中： Kv為地基豎向基床換算系數(shù)（N/m2）； K0為地基豎向基床系數(shù)（N/m2）；D 為地下管線外徑。

目前很多學(xué)者均采用Vesic[11]提出的地基基床系數(shù)公式

式中： Es為地基土彈性模量； EI為管線的等效抗彎剛度； νs為地基土泊松比。

Vesic 公式是基于置于地表的彈性地基梁提出的，沒有考慮埋深對(duì)地基豎向機(jī)床系數(shù)的影響。俞劍等[12]改進(jìn)了Vesic 的地基基床系數(shù)，可以考慮地基土埋深的影響，本文將地下管線看做埋于地下的彈性地基梁，因此俞劍等[12]提出的模型公式更為符合本文實(shí)際情況，公式如下：

式中：η 為與管線埋深h 和管線外徑D 相關(guān)的系數(shù)。當(dāng)h/D ≤0.5時(shí)， η=2.18；當(dāng)h/D ＞0.5， η=1+。

1.2 地下管線的豎向變形理論

在Winkler 彈性地基梁模型中，假定土體與管線之間保持彈性接觸，用連續(xù)分布的彈簧來模擬管-土之間的相互作用，土體與管線之間滿足變形協(xié)調(diào)，不可發(fā)生分離。計(jì)算模型見圖1。

圖 1 Winkler 彈性地基梁模型Fig.1 Winkler elastic foundation beam model

根據(jù)Winkler 彈性地基梁理論，地下管線豎向位移的微分方程為

式中， q0為地下管線的自重、管內(nèi)水的重量以及上附土壓力的合力。計(jì)算時(shí)，取q0=0。式（6）為高階微分方程，直接求解較為復(fù)雜，本文采用有限單 元法求解該方程。

圖 2 管線離散示意圖Fig.2 Schematic diagram of pipeline discrete

圖2 為對(duì)地下管線做了有限差分離散化的示意圖。管線離散為(n+5)個(gè)節(jié)點(diǎn)單元（其中管線兩端各設(shè)置2 個(gè)虛擬節(jié)點(diǎn)單元），每個(gè)單元的長(zhǎng)度設(shè)為l。根據(jù)有限差分法原理，式（6）中的微分項(xiàng)可表示為

式中： yi，yi+1，yi-1，yi+2，yi-2分別為節(jié)點(diǎn)i，i+1，i-1，i +2，i -2處的位移。

由此可得式（6）的有限差分表達(dá)式為

邊界條件：假定管線兩端自由，則有兩端的彎矩M、剪力Q 均為0，即

由式（9）可得4 個(gè)虛節(jié)點(diǎn)的位移

式（6）可寫成矩陣的形式

并可簡(jiǎn)化為

式中：Kt為管線變形剛度矩陣；Ks為地基剛度矩陣； y為管線豎向位移列向量；Q 為管線處土體附加變形列向量。其中各變量的具體表達(dá)形式如下：

令K=Kt+Ks， 則式（12）可改寫為

式中，K-1為 K的逆矩陣。因此，在已求得管線上附加應(yīng)力的情況下，可通過式（13）求得基坑開挖引起周圍地下管線的變形。

1.3 作用在管線上的水平附加應(yīng)力

計(jì)算管線水平方向位移時(shí)，可將其看做一個(gè)水平放置的彈性地基梁?；娱_挖，支護(hù)結(jié)構(gòu)會(huì)向基坑內(nèi)側(cè)方向隆起從而達(dá)到新的平衡，埋在土中的管線也會(huì)隨之發(fā)生變形。設(shè)在管線埋深處基坑支護(hù)結(jié)構(gòu)變形曲線方程為

式中，a，b，c 均為常數(shù)，其值可由以下邊界條件確定：

因此，基坑支護(hù)變形引起的管線處水平附加應(yīng)力為

式中， Kh為地基水平基床換算系數(shù)。根據(jù)劉建航等[4]的研究， Kh的 取值為 Kv的1.5～2.0 倍。其余參數(shù)與上文相同。

因此，可以得到管線水平方向附加應(yīng)力為

1.4 地下管線的水平方向變形理論

將管線看做水平放置的彈性地基梁，可得其位移微分方程為

求解該方程，可參照上文管線豎向位移微分方程的解法。

2 工程實(shí)例

即墨變電站位于浦東新區(qū)陸家嘴地區(qū)銀城路南側(cè)、昌邑路北側(cè)、日照路東側(cè)三角形地塊內(nèi)。地塊用地面積7 100 m2，其站體平面形狀布置為L(zhǎng) 型，如圖3 所示。具體平面尺寸為：矩形101.5 m×56.4 m 切去東北角35.3 m×21.4 m，站體占地面積為4 969.18 m2。地下變電站為剪力墻結(jié)構(gòu)體系，在地下一層頂板設(shè)轉(zhuǎn)換層。建筑物為地下三層、地上五層（局部四層）布置，基坑面積4 569 m2，基坑開挖深度22.0 m?；幽蟼?cè)沿昌邑路分布有一條上水管線，材質(zhì)為鋼材，管徑300 mm，管厚4 mm，彈性模量20.1 GPa，頂埋深1.0 m，距基坑邊約5.87 m。管線上設(shè)10 個(gè)監(jiān)測(cè)點(diǎn)，分別為M1，M2，M3，M4，M5，M6，M7，M8，M9，M10。

圖 3 管線監(jiān)測(cè)點(diǎn)示意圖Fig.3 Schematic diagram of pipeline monitoring points

為了計(jì)算簡(jiǎn)便而準(zhǔn)確，根據(jù)地勘報(bào)告表，土體假定為均質(zhì)土層，采用線彈性模型，重度為17.6 kN/m3，壓縮模量為6.1 MPa，泊松比為0.3。圖4 為采取本文提出的解析方法計(jì)算得到的基坑開挖引起周圍地下管線豎向位移的計(jì)算值與工程實(shí)測(cè)值的對(duì)比。可以看出，本文方法計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)的地下管線縱向變形曲線分布規(guī)律基本吻合，但數(shù)值上還存在一定差別。理論計(jì)算得到的管線沉降最大值為35.42 mm，實(shí)測(cè)管線最大沉降值為36.8 mm，相對(duì)誤差為3.8%。該誤差產(chǎn)生的主要原因?yàn)椋篴. 土層性質(zhì)的不均勻性；b. 基坑形狀不是規(guī)則的矩形，沒有對(duì)稱性；c. 地基基床系數(shù)是根據(jù)前人公式算得，沒有經(jīng)過現(xiàn)場(chǎng)實(shí)測(cè)，可能存在差異。

圖 4 基坑周圍管線豎向位移理論值與實(shí)測(cè)值Fig. 4 Theoretical and measurement values for the vertical displacement of pipelines surrounding the foundation pit

3 基坑開挖引起地下管線變形的因素分析

3.1 管線埋深

本文將管線埋深取8 種情況，分別為：1，2，3，4，5，10，15，20 m。通過圖5、圖6 可以看出，管線豎向位移隨著其埋深的增大而減小，減小速度與管線埋深增大速度不呈正比關(guān)系。管線埋深越大，其最大豎向位移越小，且其隨埋深的增大而減小的速率也越小。當(dāng)管線埋深達(dá)到20 m時(shí)，管線的豎向位移幾乎為0。

圖 5 不同埋深地下管線的豎向位移情況Fig. 5 Vertical displacements of pipelines with different depths

圖 6 埋深對(duì)地下管線的最大豎向位移的影響Fig. 6 Influence of depth on the vertical displacement of pipelines

3.2 地下管線到基坑邊緣的水平距離

圖7 為管線最大豎向位移與管線到基坑的水平距離之間的變化規(guī)律，兩者的變化規(guī)律呈V 字形，即管線與基坑邊緣的水平距離隨管線最大豎向位移先增大后減小，當(dāng)水平距離增加到一定程度后（本文為50 m），管線豎向位移趨近于0。最大豎向位移出現(xiàn)在水平距離10～15 m 之間，即0.5H 左右，這與前人[13-14]研究的坑外土體沉降規(guī)律基本一致。

圖 7 水平距離對(duì)管線最大豎向位移的影響Fig.7 Influence of horizontal distance on the maximum vertical distance of pipeline

3.3 地基土性質(zhì)（基床系數(shù)）

圖8 為地基土基床系數(shù)對(duì)管線豎向位移的影響情況?？梢钥闯觯叵鹿芫€的最大豎向位移隨地基土基床系數(shù) Kv增 加而逐漸減少，但當(dāng) Kv增大到某一值后(本例 Kv=1.2×107kN/m2)，地下管線的豎向位移增加幅度趨于平穩(wěn)，此時(shí)地基基床系數(shù)Kv的增加對(duì)豎向位移的影響已逐漸減小直至可忽略不計(jì)。

圖 8 地基土基床系數(shù)對(duì)管線豎向位移的影響Fig. 8 Influence of foundation soil coefficient of subgrade reaction on the vertical distance of pipeline

4 結(jié) 論

a. 基于Winkler 彈性地基梁理論，利用兩階段法建立了受基坑開挖影響的地下管線豎向位移及水平位移方程，并通過工程實(shí)例驗(yàn)算，證明計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)的地下管線變形值基本吻合。

b. 管線豎向位移隨著其埋深的增大而減小，且其減小的速度與管線埋深的增大速度不是線性關(guān)系。管線埋深越大，其最大豎向位移越小，且其隨埋深的增大而減小的速率也越小。

c. 管線最大豎向位移與管線到基坑的水平距離之間的變化規(guī)律為凹槽形，即先增大后減小，當(dāng)水平距離增加到一定程度后（本文為50 m），管線豎向位移趨近于零。最大豎向位移出現(xiàn)在水平距離0.5H（H 為基坑開挖深度）左右。

d. 地下管線的最大沉降值隨地基土基床系數(shù)Kv增加而減小的，但當(dāng) Kv增大到一定程度后(本例 Kv=1.2×107N/m2)，地下管線的最大沉降增加趨勢(shì)趨于平穩(wěn)，即 Kv的增加對(duì)管線沉降的影響已逐漸減小。