范志超，張啟龍

（中國(guó)電建市政建設(shè)集團(tuán)有限公司，天津 300384）

1 引言

城市水利管線建設(shè)中， 常使用頂管掘進(jìn)法開(kāi)展地下管線工程建設(shè)，尤其是遇到穿越城市湖泊、江河等工程，泥水平衡頂管掘進(jìn)機(jī)得到廣泛應(yīng)用。 當(dāng)泥水平衡掘進(jìn)機(jī)在穿越江河湖泊施工工程時(shí)， 需要在臨近湖泊附近修建沉井結(jié)構(gòu)作為頂管掘進(jìn)機(jī)始發(fā)井。 修建沉井結(jié)構(gòu)時(shí)，通常會(huì)降低湖泊水位；修建完成后，恢復(fù)湖泊原有水位。 這就容易造成沉井結(jié)構(gòu)不同外表面承受不同水壓力作用，而且這些水壓力是動(dòng)水壓力。 在沉井結(jié)構(gòu)受力研究方面也有了一些初步探討。李雨[1]鑒于當(dāng)前部分市政工程沉井結(jié)構(gòu)施工不合理的現(xiàn)狀， 提出了與之相適應(yīng)的市政結(jié)構(gòu)頂管工程中沉井結(jié)構(gòu)施工方案。肖剛[2]結(jié)合某工程實(shí)例，開(kāi)展了市政結(jié)構(gòu)頂管工程沉井結(jié)構(gòu)設(shè)計(jì)相關(guān)探討研究。 鄭一葦[3]以某市政結(jié)構(gòu)頂管工程為主要研究對(duì)象，分析市政工程頂管施工中的沉井結(jié)構(gòu)優(yōu)化設(shè)計(jì)。梁峰[4]提出了通過(guò)對(duì)沉井設(shè)計(jì)工作的優(yōu)化和完善，提高沉井結(jié)構(gòu)的完整性和安全性，提高整個(gè)頂管工程施工工作的質(zhì)量水平。夏國(guó)燕[5]以沉井結(jié)構(gòu)實(shí)際案例為例，分析了市政頂管工程沉井結(jié)構(gòu)的設(shè)計(jì)要點(diǎn)，討論了沉井結(jié)構(gòu)設(shè)計(jì)的相關(guān)應(yīng)用。

本文結(jié)合典型頂管始發(fā)沉井結(jié)構(gòu)工程， 采用數(shù)值計(jì)算分析方法，建立了動(dòng)水壓作用下沉井結(jié)構(gòu)數(shù)值計(jì)算模型，研究了動(dòng)水壓作用下沉井結(jié)構(gòu)整體受力性能， 為臨近湖泊地下工程結(jié)構(gòu)設(shè)計(jì)與施工提供重要的參考。

2 工程概況

珠三角水資源配置工程?hào)|莞配套松山湖水廠一期管線工程，位為松山湖、大朗鎮(zhèn)、大嶺山鎮(zhèn)。 工作內(nèi)容包括5 條配水管線和1 條原水管線。 原水管線管徑為DN3 200， 管材采用鋼管，管道長(zhǎng)度約為420 m，配管線總長(zhǎng)度約12 864 m，管徑范圍為1 000～2 400 mm， 采用球墨鑄鐵管和鋼管， 分為新城路段（不含涉鐵段）、新城路段（涉鐵段）、西線環(huán)湖路段、北線環(huán)湖路段、青田路段、楊朗路段、美景西路段。

頂管掘進(jìn)機(jī)始發(fā)與接收井均采用沉井結(jié)構(gòu)， 深井位置如圖1 所示，其地層主要包括素填土、強(qiáng)風(fēng)化花崗巖（土）、中風(fēng)化花崗巖。 風(fēng)化花崗巖數(shù)弱透水層。 頂管工作井為圓形，直徑（內(nèi)徑）14.0 m，上部區(qū)域井壁厚1.0 m，底部韌腳壁厚1.2 m，井壁及底板采用C30 鋼筋混凝土抗?jié)B等級(jí)P6，封底采用C25混凝土等。 其施工基本方法：沉井下沉根據(jù)地下水位情況，采用不排水挖土下沉方法。 采用挖機(jī)挖掘井底中央部分的土，使沉井形成鍋底，挖土須分層、對(duì)稱、均勻地進(jìn)行，由沉井中間開(kāi)始逐漸挖向四周，每層高0.4～0.5 m。 沉井每節(jié)的混凝土達(dá)到100%強(qiáng)度后方可下沉，下沉?xí)r首先需將沉井井壁上全部預(yù)留孔加以封堵，對(duì)于較大的孔用M10 水泥砂漿砌機(jī)磚封堵并設(shè)加強(qiáng)筋，在內(nèi)外面用1∶2 水泥砂漿抹面厚20 mm。 沉井下沉至設(shè)計(jì)標(biāo)高后觀察其穩(wěn)定性， 在8 h 內(nèi)沉井自沉累計(jì)量不大于10 mm 時(shí)，才能進(jìn)行封底。

圖1 沉井結(jié)構(gòu)位置

3 數(shù)值計(jì)算模型與分析方法

3.1 數(shù)值計(jì)算模型

采用ABAQUS 有限元軟件對(duì)工程進(jìn)行模擬研究。 對(duì)深井的模擬分析，需要先建立深井及周?chē)临|(zhì)的本構(gòu)模型。 根據(jù)工程項(xiàng)目背景， 模擬不同水位對(duì)深井的受力影響以及不同工況對(duì)深井沉降的影響。 各土層均采用摩爾-庫(kù)倫本構(gòu)模型，深井則采用彈性本構(gòu)模型進(jìn)行模擬。 模型示意圖細(xì)節(jié)如下：下部總長(zhǎng)度80 m，左側(cè)深30 m，右側(cè)深18 m，深井依據(jù)結(jié)構(gòu)尺寸，高22 m，寬16 m，壁厚1 m；以深井為分界，深井左側(cè)中風(fēng)化花崗巖深度27 m，粉質(zhì)黏土深3 m，右側(cè)依次為16 m 和2 m。 土層和深井均采用2D 單元。 網(wǎng)格密度采用1 m 進(jìn)行劃分，部分進(jìn)行網(wǎng)格加密。

為建立有效的數(shù)值模擬， 需要對(duì)模型賦予基本物理力學(xué)參數(shù)。 根據(jù)地質(zhì)勘察報(bào)告， 對(duì)物理模型中各計(jì)算參數(shù)賦值如下，深井：彈性模量3.35×104MPa，泊松比0.3，密度23 kN/m3。各層土體的力學(xué)參數(shù)見(jiàn)表1。

表1 模型計(jì)算參數(shù)

3.2 數(shù)值計(jì)算分析方法

深井施工前和投入使用的過(guò)程中， 周?chē)w產(chǎn)生的水頭壓力對(duì)深井的使用有顯著的影響。 基于大型數(shù)值計(jì)算程序ABAQUS 建立深井在不同水位的力學(xué)計(jì)算模型，有序改變水頭壓力，根據(jù)工程實(shí)際分為3 個(gè)工況，分析3 個(gè)工況下深井的變形與滲透規(guī)律，以及內(nèi)力分布的變化規(guī)律。

工況一：考慮深井施工時(shí)的工作水位，湖水深度為8 m；

工況二：深井使用期間的正常水位，①面為水位線，深度24 m，③面深度為18 m，在③面上施加3×104Pa 的應(yīng)力邊界和3×104Pa 的孔隙壓；

工況三：深井使用期間的最高水位，①面為水位線，深度24 m， 在③面上施加7×104Pa 的應(yīng)力邊界和7×104Pa 的孔隙壓。 模型水位線示意圖如圖2 所示（②面為枯水期的水位線）。

圖2 模型水位線示意圖

4 沉井結(jié)構(gòu)整體受力性能分析

4.1 深井水壓力隨水位的變化特征

深井結(jié)構(gòu)在不同工況下， 深井的側(cè)壁孔隙水壓力分布出現(xiàn)明顯變化。 深井前后壁的孔隙水壓力隨著水位的升高而增大，其中，最大值位于深井的底部右側(cè)。 工況一由于水位在深井下面，所以，深井的孔隙水壓力在軟件云圖中顯示是負(fù)值。在最大水頭壓力 （工況三） 下深井底部右側(cè)的孔隙水壓力為78.09 kPa， 工況二深井底部右側(cè)的孔隙水壓力為62.4 kPa，增大了約25.14%。 在深井上布置孔隙水壓力測(cè)點(diǎn)，選擇底部以及靠近下部的點(diǎn)作為測(cè)點(diǎn)， 測(cè)點(diǎn)1、2、3 和測(cè)點(diǎn)5、6、7 分別位于深井左側(cè)和右側(cè)，相鄰側(cè)點(diǎn)相隔1 m，測(cè)點(diǎn)4 位于底部中點(diǎn)。

在深井布置7 處孔隙水壓力測(cè)點(diǎn)， 表2 列出了7 個(gè)測(cè)點(diǎn)在兩種工況下的孔隙水壓力值。 通過(guò)比較可以得出：深井的整體孔隙水壓力與水頭壓力是呈正比例關(guān)系，其中，測(cè)點(diǎn)7 的增長(zhǎng)比例達(dá)到32.5%，在工況二的水頭壓力影響下，深井右側(cè)所承受的孔隙水壓力比左側(cè)大了10%～20%。 在相同工況下，深井底部的孔隙水壓力最大，右側(cè)次之，左側(cè)最小。

表2 不同工況下測(cè)點(diǎn)的孔隙水壓力

4.2 深井隨水位的變化特征

根據(jù)3 種工況下流固耦合作用，在不同水頭壓力作用下，選取深井結(jié)構(gòu)為分析對(duì)象， 可以得到3 種工況下的內(nèi)力分布云圖。 從深井應(yīng)力和剪力云圖分析可知：深井的應(yīng)力和剪力在不同工況下的分布規(guī)律大致相同， 最大應(yīng)力均位于深井左壁底部，最大剪力在底部靠左壁的位置。 由深井剪力云圖可知正剪應(yīng)力主要分布在深井側(cè)壁和底部左半部分， 負(fù)剪應(yīng)力主要分布在底部右半部分。

4.3 深井隨水位的沉降變形

深井結(jié)構(gòu)在不同工況下， 深井的豎向位移分布出現(xiàn)明顯變化。在流固耦合作用下，3 種工況不同動(dòng)水壓力下，工況一的位移沉降整體最小，因?yàn)樗辉谏罹?，由軟件云圖可知右壁的位移比左臂的位移大。 工況二與工況三的位移變化分布大致相似，但由于工況三的水位高3 m，整體位移比工況二大了5 mm 左右。

根據(jù)深井布置的7 個(gè)測(cè)點(diǎn)， 列出了7 個(gè)測(cè)點(diǎn)在3 種工況下的位移。 通過(guò)比較可以得出：深井的右壁位移與測(cè)點(diǎn)位置是呈正比例關(guān)系，隨著測(cè)點(diǎn)逐漸向上移動(dòng)，位移逐漸增加；工況三測(cè)點(diǎn)7 位移2.93 cm 比工況二2.49 cm 增大了4.4 mm，對(duì)比測(cè)點(diǎn)1 和測(cè)點(diǎn)7 在工況二和工況三變化， 測(cè)點(diǎn)1 增大了16.8%，測(cè)點(diǎn)7 增大了17.7%，在動(dòng)水壓的增大的影響下，深井右壁的位移變化比左壁更大，深井的不均勻變形，會(huì)導(dǎo)致深井產(chǎn)生開(kāi)裂，產(chǎn)生的裂縫會(huì)加大周?chē)臐B透程度，影響深井的使用。

5 結(jié)語(yǔ)

針對(duì)動(dòng)水壓力作用下頂管始發(fā)沉井結(jié)構(gòu)整體受力特征，采用數(shù)值計(jì)算分析方法， 建立了動(dòng)水壓作用下沉井結(jié)構(gòu)數(shù)值計(jì)算模型，研究了動(dòng)水壓作用下沉井結(jié)構(gòu)整體受力性能，為臨近湖泊地下工程結(jié)構(gòu)設(shè)計(jì)與施工提供重要的參考。