降雨入滲對膨脹土隧道力學特性影響的數(shù)值分析

顧永貴

(中鐵二十四局集團安徽工程有限公司，安徽合肥 230011)

膨脹土具有脹縮性、裂隙性和超固結性等工程性質(zhì)。工程危害性表現(xiàn)在土中含水率變化時，其富含的膨脹性巖土礦物會產(chǎn)生顯著的濕脹干縮變形，對路基、邊坡、隧道支護結構等造成不利影響。隧道穿越膨脹土地區(qū)時，其支護結構受到膨脹土濕脹干縮影響引起的變形損壞已成為目前比較普遍的問題。國內(nèi)外學者針對膨脹土增濕過程中的工程特性進行了很多研究[1-5]。周坤等[6]通過有限元方法模擬隧道周邊膨脹土吸水后對隧道襯砌的影響，認為膨脹率對隧道結構影響最大，其次是膨脹圈厚度，對隧道結構作用最小的是埋深。林剛等[7]基于膨脹接觸壓力的概念，結合成都地鐵2號線膨脹巖土地層結構分析，計算出不同膨脹巖土分布、膨脹力、埋深以及拼裝方式下的結構內(nèi)力值。郭瑞等[8]采用數(shù)值模擬計算，認為膨脹土存在于隧道上、下方時，其膨脹會導致管片結構的變形加劇和內(nèi)力增大，同時提出加強配筋或通縫拼裝等措施，以降低管片結構受力。曾仲毅等[9]結合熱傳導膨脹模擬增濕膨脹的替代方程、膨脹力及滲流參數(shù)，通過有限差分軟件模擬不同膨脹力條件下的增濕過程，得出引發(fā)隧道支護結構破壞的臨界含水率為19.5%，臨界膨脹力為140 kPa。

眾所周知，含水率變化是影響膨脹土力學性質(zhì)的主要因素，而隧道的擾動將導致圍巖裂隙的發(fā)展，滲水將加大圍巖的濕度，從而對襯砌結構的膨脹力產(chǎn)生影響。目前，對不同擾動程度下膨脹土隧道的襯砌結構力學特征缺少全面的定量研究。因此，通過定義隧道施工的松動區(qū)(吸濕范圍)，基于PLAXIS二維數(shù)值模型，以溫度場模擬濕度場[10-11]，仿真模擬隧道施工過程中圍巖膨脹區(qū)滲水增濕的過程，分析不同松動區(qū)的濕度場變化趨勢，研究膨脹土增濕膨脹對隧道襯砌結構的受力變化特征。

1 數(shù)值模型及模擬方案

1.1 研究對象

合肥地鐵五號線某下穿暗挖隧道左線長53.9 m，右線長42.9 m，隧道寬7.6 m，高8 m，埋深15.0 m，如圖1所示。暗挖段施工采用CRD法，以1→3→2→4的順序開挖；3與4和1與3均以5 m間距施工，每步進尺1.5～2.0 m。采用“噴射混凝土(C25，厚25 cm)+鋼筋網(wǎng)+工字鋼(工22a)”進行初期支護。

圖1 新奧法隧道開挖分區(qū)示意

1.2 計算模型及參數(shù)

采用Plaxis 2D建立數(shù)值分析模型，模型“寬×高”為100 m × 40 m，如圖2所示。土層自地表(高程為±0.000)分布如下：填土(厚2 m)、黏土2(厚6 m)、黏土3(厚17 m)和黏土4(厚15 m)。隧道拱頂高程為-15 m，地下水位高程為-30 m。通過15節(jié)點高階三角形單元離散地層，總計分為1 062個單元，8 920個節(jié)點(如圖3所示)。

圖2 土體分層、水位及隧道結構幾何模型

圖3 隧道與土體有限元網(wǎng)格及邊界條件

應用HS模型作為巖土體的本構模型，可模擬關于復雜加卸載開挖類巖土體的力學行為[12]，按照固結快剪指標調(diào)整土體的強度參數(shù)。其中，通過數(shù)值模擬反演聯(lián)合表面浸潤及膨脹力試驗結果確定熱參數(shù)(見表1)；工字鋼(工22a)和初期支護的混凝土噴射則利用板單元進行模擬(見表2)。

表1 土體物理力學參數(shù)

表2 板單元參數(shù)

1.3 松動區(qū)范圍定義

隧道開挖卸載導致應力劇增，一定區(qū)域內(nèi)的圍巖土體將會產(chǎn)生變形松動，在隧道周圍出現(xiàn)“松動區(qū)”。此處定義只在“松動區(qū)”發(fā)生吸濕，表層填土無膨脹性，不計入“松動區(qū)”。在巖石工程領域中，主要從隧道開挖導致的圍巖應力和應變變化兩方面研究“松動區(qū)”，分析和模擬膨脹土的工程特性，從隧道開挖而導致圍巖偏應力及位移變化過程界定“松動區(qū)”的范圍。隧道開挖前，巖土體偏應力大致呈水平分布，隨著開挖作業(yè)展開，在開挖卸載和支護結構阻力共同作用下，圍巖出現(xiàn)偏應力集中現(xiàn)象，偏應力等值線的集中范圍可定義成“松動區(qū)”(見圖4)。初始階段，巖土體尚未發(fā)生移動，隧道開挖后，洞壁周邊土體朝臨空面發(fā)生位移，可定義出現(xiàn)較大位移的區(qū)域為“松動區(qū)”(見圖5)?？紤]增濕膨脹的隧道開挖松動區(qū)界定原則是隧道周圍左、右、下方一倍洞徑處和拱頂上方地表填土層底面處，并認為松動區(qū)域的上層潛水易滲入膨脹性黏土，導致黏土膨脹及強度降低，引起隧道支護結構受力產(chǎn)生變化。

圖4 隧道開挖引起周圍土體偏應力分布等值線

圖5 隧道開挖引起周圍土體位移分布等值線

1.4 模擬方案

溫度場與滲流場具相似性。因此，可將滲流場土體滲透系數(shù)的變化等價為溫度場下土體導熱系數(shù)的變化[6]，利用溫度場的土體升溫膨脹來模擬膨脹土的增濕膨脹特性[9]。隨著隧道開挖后圍巖產(chǎn)生松動區(qū)，松動區(qū)應力場、變形場發(fā)生較大改變，土體出現(xiàn)裂隙，引起滲透系數(shù)的變化。由此可見，松動程度不同，相應的滲透系數(shù)也不同。在設計模擬方案時，通過提高膨脹土導熱系數(shù)至原始導熱系數(shù)的100倍、1 000倍、10 000倍，來考慮隧道開挖松動區(qū)在這三種不同的擾動程度下土體導熱系數(shù)的變化(見表3)。隨開挖斷面的變化，隧道周圍的松動區(qū)范圍也隨之改變。為方便研究，假定松動區(qū)范圍在隧道開挖過程中保持不變。

表3 隧道開挖松動區(qū)不同擾動程度下的土體導熱系數(shù) kW/m/K

(1)邊界條件

為了模擬上層滯水滲入膨脹土松動區(qū)產(chǎn)生增濕膨脹而導致隧道結構出現(xiàn)的附加變形及應力，應設置相應的熱邊界條件。因地下水位在地面以下30 m，所以只需考慮賦存于地表填土層中水體沿松動區(qū)的下滲。將地表填土層底面(即黏土2層頂面)定義為增濕邊界，即于模型高程-2.000 m處設置成303.1 K(與增濕最大飽和度相對應)。上層滯水下滲不作用于模型的左右側面以及下底面，定義為封閉邊界(溫度零變化，熱量零交換，飽和度不變)。

(2)施工階段設置

隧道28 d完成一次全斷面開挖。受降雨等因素影響，賦存于地表填土層中的水體得到充足補給，隨著隧道開挖出現(xiàn)松動區(qū)，上層水體滲入，引起膨脹土松動區(qū)增濕膨脹，從而對隧道結構變形及受力產(chǎn)生影響[13-14]。為模擬隧道開挖后引起膨脹土增濕膨脹的過程和開挖導致支護結構的變形情況，可按照表3所定義松動區(qū)的4種不同松動程度，根據(jù)表4的施工階段設置，分別完成膨脹土隧道的增濕模擬。

表4 膨脹土隧道增濕膨脹模擬施工階段設置

2 計算結果及分析

2.1 土壓力發(fā)展規(guī)律

為觀察隧道周邊松動區(qū)域土體增濕膨脹后膨脹力的發(fā)展情況，在洞周上下左右距離洞壁0.5 m處分別設置一條剖線(如圖6所示)，輸出各施工階段增濕前后剖線上的水平土壓力和豎向土壓力變化情況，如表5～表7所示。

從表5～表7可以看出，隨著擾動程度增大，洞周土體增濕前后的土壓力差異也逐漸增大。

在“輕微松動”條件下，溫度影響深度尚未達到隧道所處的黏土3地層，增濕前后洞周各剖線土壓力基本沒有變化，洞周土體內(nèi)基本不產(chǎn)生膨脹力。

圖6 隧道周邊土壓力剖線布置示意

表5 隧道周邊膨脹土增濕前后應力變化(輕微松動) kPa

表6 隧道周邊膨脹土增濕前后應力變化(顯著松動) kPa

表7 隧道周邊膨脹土增濕前后應力變化(完全松動) kPa

在“顯著松動”條件下，溫度影響深度已經(jīng)達到隧道所處的黏土3地層，但尚未達到增濕飽和，洞周各剖線土壓力有一定提高，表明洞周土體內(nèi)將產(chǎn)生少量的膨脹力。

在“完全松動”條件下，洞周土體達到增濕飽和狀態(tài)，洞周各剖線土壓力顯著提高，表明此時洞周土體增濕后將產(chǎn)生較大的膨脹力，繼而也會使得隧道支護結構所受的附加荷載明顯增大。

2.2 結構內(nèi)力發(fā)展規(guī)律

在不同擾動程度下，各施工階段洞周土體增濕膨脹后隧道初支結構的內(nèi)力和變形響應分布情況如圖7～圖12所示。圖13給出了各施工步序下，土體的不同擾動程度(不同的導熱系數(shù)及產(chǎn)生的膨脹力)對應的彎矩極值變化情況。

由圖7～圖13可以看到，隨著隧道開挖擾動程度的加劇，隧道初期支護的受力和變形與不考慮增濕膨脹條件下的內(nèi)力和變形之間的差異逐步增大。擾動程度越大，洞周松動區(qū)域浸水增濕速度越快，增濕飽和程度越大，地層壓力亦隨之增大，導致隧道初支結構承受較大的附加荷載作用，土體增濕膨脹后隧道初支結構內(nèi)力最大增幅接近20%。

圖7 上層水入滲增濕后隧道結構彎矩分布(輕微松動)

圖8 上層水入滲增濕后隧道結構變形(輕微松動)

圖9 上層水入滲增濕后隧道結構彎矩分布(顯著松動)

圖10 上層水入滲增濕后隧道結構變形(顯著松動)

圖11 上層水入滲增濕后隧道結構彎矩分布(完全松動)

圖13 各施工步序下隧道初支結構彎矩極值隨土體膨脹程度的變化

3 結論

(1)根據(jù)隧道開挖而導致的圍巖偏應力及位移變化過程來界定“松動區(qū)”的范圍。考慮增濕膨脹的隧道開挖松動區(qū)界定為隧道周圍左、右、下方一倍洞徑處和拱頂上方地表填土層底面處，認為松動區(qū)域的上層潛水易滲入膨脹性黏土，導致黏土膨脹及強度降低，引起隧道支護結構受力產(chǎn)生變化。

(2)在“輕微松動”條件下，增濕前后洞周各剖線土壓力基本沒有變化，洞周土體內(nèi)基本不產(chǎn)生膨脹力；在“顯著松動”條件下，洞周各剖線土壓力有一定提高，表明洞周土體內(nèi)將產(chǎn)生少量的膨脹力；在“完全松動”條件下，隧道開挖完成時洞周土體達到增濕飽和狀態(tài)，洞周各剖線土壓力顯著提高，表明此時洞周土體增濕后將產(chǎn)生較大的膨脹力，繼而也會使得隧道支護結構所受的附加荷載明顯增大。

(3)隨著隧道開挖擾動程度的加劇，隧道初期支護的受力和變形與不考慮增濕膨脹條件下的內(nèi)力和變形之間的差異逐步增大，土體增濕膨脹后隧道初支結構內(nèi)力最大增幅接近20%。

(4)在目前的膨脹土隧道數(shù)值模擬研究中，普遍采用指定吸濕膨脹范圍的方法，本文中對指定松動區(qū)內(nèi)考慮吸濕膨脹也屬于這一類。在借鑒松動圈概念的基礎上，根據(jù)變形、應力梯度等指標劃定松動區(qū)范圍，但是相關指標仍然需要人為指定，導致吸濕膨脹范圍的劃定不夠精確，后續(xù)可以對吸濕膨脹范圍的精確確定方法作進一步研究。

(5)膨脹土的膨脹變形與初始含水率、終了含水率和上覆荷載等均有關系。膨脹土含水率變化過程實際是非飽和土滲流過程，具有強非線性。膨脹土吸水膨脹可通過引入膨脹系數(shù)(采用初始應變法)進行有限元分析。但膨脹系數(shù)與多因素相關，如應力狀態(tài)、初始含水率等，對某一種特定膨脹土不是一個常數(shù)。由于試驗資料有限，未考慮增濕過程中膨脹系數(shù)的變化，而是指定了幾種極端的不利工況進行分析，后續(xù)可加強相關試驗研究。

[1] RAO S M， REDDY B V V， MUTTHARAM M. The impact of cyclic wetting and drying on the swelling behaviour of stabilized expansive soils[J]. Engineering Geology， 2001，60(7):223-233

[2] Patrick ChegeKariuki， Freek van der Meer. A unified swelling potential index for expansive soils[J]. Engineering Geology， 2004，72(3):1-8

[3] BanuIkizler S， Mustafa Aytekin， Mustafa Vekli， et al. Prediction of swelling pressures of expansive soils using artificial neural networks[J]. Advances in Engineering Software， 2010，41(4):647-655

[4] 丁振洲，鄭穎人，李利晟.膨脹力變化規(guī)律試驗研究[J].巖土力學，2007，28(7):1328-1332

[5] 朱豪，王柳江，劉斯宏，等.南陽膨脹土膨脹力特性試驗[J].南水北調(diào)與水利科技，2011，9(5):11-14

[6] 周坤.膨脹土隧道襯砌膨脹力數(shù)值模擬研究[D].成都：西南交通大學，2007

[7] 林剛，羅世培，郭俊，等.膨脹巖土地層盾構隧道結構力學行為研究[J].現(xiàn)代隧道技術，2011，48(3):74-79

[8] 郭瑞，何川，方勇.膨脹土地層中盾構隧道管片結構受力分析與對策研究[J].現(xiàn)代隧道技術，2010，47(6):17-22

[9] 曾仲毅，徐幫樹，胡世權，等.增濕條件下膨脹土隧道襯砌破壞數(shù)值分析[J].巖土力學，2014，35(3):871-880

[10] 晁峰，王明年，舒東利，等.膨脹土附加應力場對深基坑土壓力的影響[J].中南大學學報(自然科學版)，2016，47(9):3177-3182

[11] 繆協(xié)興.用濕度應力場理論分析膨脹巖巷道圍巖變形[J].中國礦業(yè)大學學報，1995，24(1):58-63

[12] 徐中華，王衛(wèi)東.敏感環(huán)境下基坑數(shù)值分析中土體本構模型的選擇[J].巖土力學，2010，31(1):258-264

[13] 王清標，張聰，溫小康，等.膨脹土脹縮機理及其在控制隧道沉降中的應用[J].地下空間與工程學報，2014，10(6):1426-1432

[14] 王龍.四川盆地某邊坡膨脹土的力學特性試驗研究[J].水利水電技術，2016，47(5):43-45