矩形頂管頂推力對上部地表位移影響規律研究

張建安

（上海市政工程設計研究總院（集團）有限公司，上海市 200092）

矩形頂管頂推力對上部地表位移影響規律研究

張建安

（上海市政工程設計研究總院（集團）有限公司，上海市 200092）

矩形頂管憑借空間利用率高、覆土淺等特點在市政工程項目中應用得越來越多，但類似的工程案例與研究卻有限。以某過街通道項目為研究對象，采用Midas/GTS軟件進行三維數值模擬，分析矩形頂管施工過程中對上部地表位移的影響。研究發現頂推力為1.0P（P為開挖面中心處土體豎向壓力）時矩形頂管對掘進面前方土體影響范圍相對最小，為1.0D（D為管片長邊尺寸），開挖斷面橫向影響范圍為±3.5D。同時探討了頂管上部地表沿頂進方向水平位移的變化機理。

矩形頂管；數值模擬；影響范圍

0 引言

近年來隨著城市地下空間的開發，頂管施工越來越多的應用到各類市政基礎設施建設中，其中矩形頂管對比傳統的圓形頂管有空間利用率高、覆土淺等優點，因此在過街通道、地鐵出入口等項目中得到越來越多的應用[3]。

目前對于矩形頂管施工引起的土體變形計算方法主要有：經驗法、解析法、隨機介質理論、數值模擬計算等。其中郭亮[1]結合鄭州中州大道工程探討了頂管施工中的土體沉降規律。溫鎖林[2]根據上海軌道交通2號線東延伸段工程金科路站4號出入口地下聯通道工程分析了矩形頂管施工對環境的影響。龐臣軍等[4]結合包頭市阿爾丁大街(內蒙古科技大學)地下通道建設工程，采用數值模擬研究分析矩形頂管施工影響范圍。王日東[5]利用理論推導對矩形頂管推進引起的三維土體變形提供了簡單的理論計算方法。林曉慶[6]采用FLAC-3D進行數值模擬分析了矩形頂管施工引起的土體及鄰近地下管線變形和應力分布特點。

本文以某過街通道項目（以下簡稱本工程）為例，采用Midas/GTS軟件進行三維數值模擬方法，結合其他研究人員的結論進行驗證，分析矩形頂管施工頂推力等對上部地表位移的影響，為今后類似工程設計施工提供參考。

1 工程概況

本工程于下穿道路南側設置始發井，平面投影尺寸為9.4 m×9.2 m；北側設置接收井，平面尺寸為5.5 m×8.6 m。頂進段全長44.9 m，矩形管片每節長1.5 m，共30節，管片斷面凈尺寸為5.5 m× 3.3 m（寬×高），壁厚500 mm，頂面覆土厚度為4.85 m。通道及管片示意圖見圖1、圖2。

工程范圍內土層分布自上而下依次為：1-1層雜填土，土質松軟；2-2層粉土，中壓縮性，工程地質特性較差；2-4層粉砂，中等壓縮性，工程地質特性一般；3層粉質黏土，工程地質特性較好，5層中風化巖，工程地質特性好。各土層相應物理參數詳見表1，新建頂管通道主要在粉土及粉砂層內頂進施工。

圖1 通道平面位置示意圖

圖2 頂推管節斷面圖（單位：mm)

表1 土層及注漿材料物理參數

2 三維數值模擬

2.1 計算假設

參考其他研究人員對頂管施工的數值模擬研究工作[5-7],文中對計算模型做如下的假定:（1）土體為均質、各向同性、理想彈塑性體；（2）不考慮管道接頭的影響，將其為各向同性的線彈性體；

（3）頂管正面推進力為矩形均布荷載，地層損失沿管道軸向均勻分布；

（4）頂管推進過程中不考慮土體時間效應，只考慮頂進空間位置的變化。

2.2 計算模型

分析采用Midas/GTS有限元分析軟件，綜合考慮本項目的范圍，考慮三維模型尺寸為122 m（沿通道長度）×70 m（寬度）×29.7 m（深度），如圖3所示。始發井、接收井結構、頂管管片材料參數等詳見表2。

圖3 分析模型示意圖

表2 結構材料物理力學參數

2.3 計算工況

本次分析僅模擬頂管頂進施工過程，始發井及接收井在頂進前已施工完成。模擬分析時，通過激活、鈍化不同區域單元模擬矩形頂管頂進過程，每頂進一節管片為一步施工工況。通過修改管片外一層單元參數來模擬管片外注漿。由于共頂進30節管節故本工程可分為30個推進步，結合初始階段和最后機頭進入接收井的施工，可將本工程分為32個施工工況。施工工況示意如圖4所示。

圖4 施工工況示意圖

3 結果分析

3.1 地面沉降分析

頂管施工引起土體位移的主要因素有以下幾種：（1）土體損失；（2）正面附加頂推力；（3）頂管機頭、管片與土體間的摩擦力。其中正面頂推力對穩定頂進面，減小地面沉降變形尤為重要。根據土力學基本理論，開挖面土體的側壓力可取K0γh（K0為土體側壓力系數，一般可取0.5～0.7）。

選取正面頂推力為頂管中心深度土體側向力0.7P（kPa），頂進22.5 m時（工況16，機頭行至地面道路中心線處）沿頂進軸線方向地表沉降曲線如圖5所示，根據曲線分析可知，地表沉降最大值約11.2 mm。頂進面前方約10 m，即約1.5D（D為管片長邊尺寸）時地表土體受到擾動開始沉降，頂進面后方7.5 m（約1.2D）之后土體沉降基本達到穩定的最大值，可見土體在頂進面前后2.7D范圍之內完成大部分沉降變形，后續隨著頂進機頭的遠離，土體的沉降幾乎不再有所變化，該影響范圍較其他研究人員的結論[4，6]（約3D）相似。

圖5 頂推力為0.7P時頂進至22.5 m處地表沉降曲線

為分析機頭正面頂推力對沉降的影響，分別計算頂推力為0.5P、0.7P、1.0P、1.5P、2.0P時機頭頂進至22.5 m處的頂管中心線上方地表沉降，如圖6所示，頂進面后方土體沉降變形及沉降范圍基本相似，但機頭前方土體變形則有一定差異，放大該區域如圖7所示。

圖6 不同頂推力時頂進至22.5 m處地面沉降曲線

圖7 不同頂推力時頂進至22.5 m處機頭前方地面沉降曲線

由圖7可知，在本工程中，當機頭頂推力大于1.0P時，頂進面前方地表由于土體被擠壓發生隆起，2.0P時隆起最大值約1.0 mm；當機頭頂推力小于1.0P時，頂進面前方土體發生較大范圍沉降，當選取頂推力為1.0P時，前方土體雖然也會發生沉降，但其沉降的范圍僅為1.0D，較其他頂推力情況（1.5D～3.5D）為最小。因此對于本工程，為盡量減小頂推施工對前方土體的影響范圍，選取頂推力為1.0P是合理的。

不同頂推力情況下，頂進面地表橫向沉降曲線如圖8所示。由圖可知，不同的頂推力作用下，地面橫向曲線的沉降最大值有所不同，但橫向影響范圍幾乎一致，約為3.5D，這與其他研究人員的研究結果[6]相近似。

圖8 不同頂推力時頂進至22.5 m處地面道路橫向沉降曲線

3.2 地表水平變形分析

現階段的研究更多地關注頂管施工對地表沉降的影響，對于地表的水平位移研究有限。當選用適當的頂推力時，地表會隨著頂管施工的推進而發生沿頂進方向的水平位移，如圖9所示。由圖可知，頂管上部地表會由于頂管施工發生沿頂進方向的位移，變形曲線也近似于正態分布曲線，影響范圍與地面沉降槽范圍幾乎相同，約為3.5D。

圖9 頂進不同距離時22.5 m斷面處地表沿頂進方向水平位移

提取22.5 m斷面頂管中心軸上方地表點的位移隨頂進過程的變化曲線如圖10所示。可發現當頂管機頭靠近該點正下方時，地表位移會發生與頂進相反的回縮。究其原因，當機頭還未到達地表該點正下方時，該點的水平位移主要因逐漸靠近的機頭頂進面的頂推力造成，而當機頭到達并通過該點正下方時，由于切削土體以及土層損失的作用，地表發生向機頭開挖面沉陷的位移，造成地表沉降的同時，發生與頂進方向相反的水平位移，待機頭通過后管片開始在后座千斤頂的作用下推進，雖然有注漿等減阻措施，但地表土體仍然會受到管片的摩擦作用而產生沿頂進方向的水平位移。因此減小管片與土體之間的摩擦阻力是減小地表水平位移的重要措施。

圖10 22.5 m斷面處地表中心點沿頂進方向水平位移

4 結論

本文通過數值模擬分析，研究矩形頂管施工對上部地表位移的影響，可得結論如下：

（1）矩形頂管施工時不同的頂推力作用下導致地表的沉降范圍大大不同。針對本工程而言當選取頂推力為頂進面中心處土體自重壓力的一倍時，對機頭前方土體的影響范圍最小，約1.0D（D為管片長邊尺寸）。

（2）就本工程而言矩形頂管上部地表的沉降槽及地表水平位移影響范圍約為±3.5D。

（3）造成頂管上部地表沿頂進方向的水平位移的主要因素可能有機頭頂進面頂推力的作用以及后續管片推進過程中對土體的摩擦帶動作用，其中后者的影響更大。因此選用適當的注漿措施減小管片與土體間的摩擦是控制地表水平位移變形的重要措施。

[1]郭亮.大斷面矩形頂管施工中的土體沉降規律分析[J].建筑施工,2014,36(6):731-732.

[2]溫鎖林.大斷面矩形頂管施工對環境影響研究[J].中國市政工程,2011（5）：37-39.

[3]彭立敏，王哲，葉藝超，等.矩形頂管技術發展與研究現狀[J].隧道建設,2015，35(1):1-8.

[4]龐臣軍，鮑先凱.矩形頂管施工的數值模擬研究[J].施工技術, 2013(42):410-412.

[5]王日東.矩形頂管施工引起的土體變形計算方法研究[J].佳木斯大學學報(自然科學版),2014，32(5):711-722.

[6]林曉慶.矩形頂管施工對鄰近地下管線的影響分析[D].廣州:廣州大學,2012.

[7]楊朝帥.大斷面矩形頂管隧道施工土體穩定性分析[J].低溫建筑技術,2016(4)：131-134.

U412.37+3.1

A

1009-7716（2017）08-0266-03

10.16799/j.cnki.csdqyfh.2017.08.084

2017-04-10

張建安（1987-），男，江蘇揚州人，碩士，工程師，主要從事地下工程的設計和研究。