郭昱葵， 馬啟和， 郭震山

(1. 中交第二公路勘察設(shè)計(jì)研究院有限公司， 湖北 武漢 430056; 2. 華中科技大學(xué) 土木工程與力學(xué)學(xué)院， 湖北 武漢 430074; 3. 山西省交通科學(xué)研究院 巖土與地下工程山西省重點(diǎn)實(shí)驗(yàn)室， 山西 太原 030006)

隧道施工中，由降雨引起的雨水入滲不僅會(huì)導(dǎo)致土體含水率和孔隙水壓力上升，使土體自重增加，而且會(huì)使土中基質(zhì)吸力和土體強(qiáng)度參數(shù)隨含水率的增大而明顯減小，進(jìn)而降低土體的抗剪強(qiáng)度，破壞隧道圍巖原有的穩(wěn)定狀態(tài)，可能使隧道支護(hù)體系產(chǎn)生過大的變形和內(nèi)力，嚴(yán)重時(shí)使支護(hù)結(jié)構(gòu)開裂破壞，工程安全性大大降低。因此，探討降雨入滲對(duì)隧道變形及支護(hù)結(jié)構(gòu)受力特性的影響規(guī)律具有非常重要的意義。

降雨過程中的雨水滲透大都為非飽和滲流，而非飽和土中的滲流問題一直是巖土工程研究中的熱點(diǎn)和難點(diǎn)問題[1,2]。目前，國(guó)內(nèi)外學(xué)者在非飽和土的性狀及非飽和滲流方面已進(jìn)行了相關(guān)研究。Fredlund等[3]針對(duì)非飽和土強(qiáng)度計(jì)算公式，在飽和土的摩爾庫(kù)倫抗剪強(qiáng)度公式基礎(chǔ)上考慮基質(zhì)吸力的影響，得到了非飽和土的計(jì)算公式。包承綱等[4,5]以基質(zhì)吸力為中心，結(jié)合實(shí)際工程問題討論了非飽和土的性狀。朱偉等[6]通過歸納總結(jié)，從力學(xué)角度探討了幾個(gè)非飽和土體滲流的基本問題。吳宏偉等[7,8]指出降雨入滲對(duì)工程安全有較大威脅，研究降雨入滲問題具有很重要的理論和工程應(yīng)用意義。

然而，現(xiàn)有文獻(xiàn)主要集中于對(duì)邊坡在降雨入滲過程中變形及破壞形態(tài)的研究[9～11]，而鮮有文獻(xiàn)研究降雨對(duì)隧道變形及支護(hù)受力的影響。另外，在數(shù)值模擬過程中，未能充分考慮降雨入滲過程中基質(zhì)吸力與土體強(qiáng)度參數(shù)隨含水率的變化。鑒于此，本文通過采用FLAC3D軟件，編制相應(yīng)FISH(Friendly Interactive Shell)語(yǔ)言將降雨入滲過程中基質(zhì)吸力及土體強(qiáng)度參數(shù)隨含水率的變化關(guān)系嵌入到計(jì)算過程中，對(duì)某淺埋膨脹性圍巖隧道在降雨入滲過程發(fā)生的非飽和滲流進(jìn)行數(shù)值模擬，并對(duì)降雨過程中隧道的變形及支護(hù)受力進(jìn)行深入研究。

1 非飽和滲流的數(shù)值實(shí)現(xiàn)

FLAC3D軟件中的流固耦合模塊可以完全解決飽和情況下的滲流問題[12]，該計(jì)算中非飽和區(qū)的孔壓被強(qiáng)制設(shè)為零，不能考慮基質(zhì)吸力的作用。然而，在非飽和土體中，基質(zhì)吸力隨含水率的不同而發(fā)生明顯改變，進(jìn)而會(huì)對(duì)土體抗剪強(qiáng)度產(chǎn)生不可忽略的影響[13]。因此，本文使用軟件內(nèi)置的FISH語(yǔ)言編制相應(yīng)的程序，利用孔隙介質(zhì)中基質(zhì)吸力與含水率之間的函數(shù)關(guān)系(即土-水特征曲線)改變基質(zhì)吸力。白福青等[14]以南陽(yáng)膨脹土的實(shí)驗(yàn)研究為基礎(chǔ)，通過使用濾紙法對(duì)土-水特征曲線進(jìn)行分析，發(fā)現(xiàn)了基質(zhì)吸力與含水率的經(jīng)驗(yàn)關(guān)系為：

lg(hm)=61.473w-0.357

(1)

式中：hm為基質(zhì)吸力；w為含水率。

同時(shí)，考慮降雨增濕過程中隧道圍巖的遇水軟化現(xiàn)象，將圍巖的強(qiáng)度參數(shù)(即內(nèi)摩擦角和黏聚力)與含水率的變化關(guān)系嵌入到滲流計(jì)算中，實(shí)現(xiàn)降雨過程中圍巖強(qiáng)度參數(shù)隨計(jì)算時(shí)步的實(shí)時(shí)更新。李丹等[15]以大量的直剪試驗(yàn)數(shù)據(jù)為基礎(chǔ)進(jìn)行擬合分析，得到膨脹土的強(qiáng)度參數(shù)與含水率關(guān)系表達(dá)式為：

(2)

式中：μc及μf分別為黏聚力和內(nèi)摩擦角的均值。

本文將基于式(1)，(2)，建立某淺埋膨脹性圍巖隧道的數(shù)值計(jì)算模型，模擬其降雨入滲過程。數(shù)值計(jì)算流程(圖1)為：

圖1 數(shù)值計(jì)算流程

(1)建立計(jì)算模型，賦予材料相應(yīng)的模型參數(shù)，設(shè)置初始條件及相應(yīng)的邊界條件，并平衡初始地應(yīng)力；

(2)將地應(yīng)力平衡后的位移場(chǎng)及速度場(chǎng)清零，并進(jìn)行隧道的開挖與支護(hù)；

(3)打開滲流模塊并設(shè)置滲流參數(shù)及降雨邊界條件；

(4)在每一計(jì)算時(shí)步中，提取節(jié)點(diǎn)飽和度值Sr，若飽和度不等于1，則將節(jié)點(diǎn)飽和度轉(zhuǎn)化為節(jié)點(diǎn)含水率，按照基質(zhì)吸力與含水率之間的經(jīng)驗(yàn)公式，求得節(jié)點(diǎn)的吸力值，并將其賦值給節(jié)點(diǎn)的孔隙水壓(負(fù)孔隙水壓)；若飽和度等于1，則直接進(jìn)入下一步計(jì)算；

(5)通過反距離加權(quán)差值法，將節(jié)點(diǎn)飽和度轉(zhuǎn)化為單元飽和度，并通過單元含水率求得相應(yīng)的摩擦角和黏聚力，并將其賦值給該單元；

(6)檢測(cè)地表節(jié)點(diǎn)孔隙水壓力p′，若p′>0(達(dá)到飽和)，則可通過固定該節(jié)點(diǎn)的孔隙水壓力(取p′為0)來施加壓力邊界。

(7)檢查是否達(dá)到計(jì)算終止條件，若未達(dá)到則重復(fù)步驟(4)～(7)，否則結(jié)束計(jì)算。

2 隧道非飽和滲流數(shù)值模擬

2.1 工程概況

某山嶺隧道全長(zhǎng)430 m，洞口段埋深8.4 m，是雙洞雙向行駛分離式淺埋隧道。本文取洞口段進(jìn)行研究，該區(qū)段圍巖為V級(jí)膨脹性中風(fēng)化玄武巖。隧道襯砌內(nèi)輪廓采用半徑為5.5 m的單心圓，單洞凈寬10.5 m，凈高5.0 m。

2.2 數(shù)值計(jì)算模型及參數(shù)

如圖2所示，為簡(jiǎn)化計(jì)算，將問題轉(zhuǎn)化為平面應(yīng)變問題來分析，取模型厚度1 m，寬度100 m，埋深8 m。對(duì)模型四周施加法向位移約束，并將底部固定；模型四周、底部均為不透水邊界，頂部為降雨邊界。模擬過程中，圍巖為摩爾庫(kù)倫模型，其摩擦角和黏聚力在降雨過程中隨含水率的變化取值采用文獻(xiàn)[15,16]中的試驗(yàn)結(jié)果。初期支護(hù)厚度為26 cm，采用Shell單元模擬。錨桿沿隧道環(huán)向布置，錨桿直徑32 mm，長(zhǎng)為4 m，采用Cable單元模擬。圍巖及支護(hù)結(jié)構(gòu)的力學(xué)參數(shù)匯總于表1。為了研究不同強(qiáng)度降雨對(duì)隧道變形及支護(hù)結(jié)構(gòu)受力的影響，分別取降雨強(qiáng)度為200，400，600，800，1000 mm/d，并模擬持續(xù)降雨50 h。為簡(jiǎn)化計(jì)算，圍巖滲透系數(shù)均取其飽和滲透系數(shù)5.78×10-4cm/s。

圖2 隧道數(shù)值計(jì)算模型

表1 圍巖及支護(hù)結(jié)構(gòu)力學(xué)參數(shù)

3 模擬結(jié)果分析

3.1 基質(zhì)吸力及含水率分析

文獻(xiàn)[3]基于試驗(yàn)發(fā)現(xiàn)非飽和土的抗剪強(qiáng)度可用凈法向應(yīng)力和基質(zhì)吸力兩個(gè)應(yīng)力狀態(tài)變量來描述，并提出了基于庫(kù)倫公式的非飽和土強(qiáng)度計(jì)算公式：

τf=c′+(σ-μa)tanφ′+(μa-μw)tanφb

(3)

式中：c′為有效凝聚力；φ′為有效內(nèi)摩擦角；σ為正應(yīng)力；μa為孔隙氣壓力；μw為孔隙水壓力；σ-μa為凈法向應(yīng)力；μa-μw為基質(zhì)吸力；tanφb為非飽和土的抗剪強(qiáng)度隨基質(zhì)吸力增加的速率。

從式中可以看出，在非飽和土中，基質(zhì)吸力對(duì)抗剪強(qiáng)度的影響不可忽視。在持續(xù)降雨時(shí)，雨水入滲使非飽和區(qū)土體的含水率顯著增加，導(dǎo)致基質(zhì)吸力明顯減小，進(jìn)而使得土體抗剪強(qiáng)度下降，導(dǎo)致隧道圍巖原有的穩(wěn)定狀態(tài)受到擾動(dòng)；同時(shí)隧道上覆巖體含水率增加，開挖面上覆壓力增大，支護(hù)結(jié)構(gòu)受力變形增大，工程安全性降低。因此，對(duì)降雨條件下土體含水率及基質(zhì)吸力隨降雨歷時(shí)變化規(guī)律的分析有著重要工程意義。

圖3為降雨強(qiáng)度600 mm/d工況下地表處土體含水率與基質(zhì)吸力隨降雨歷時(shí)的變化曲線。從圖中可看出，地表含水率變化與降雨歷時(shí)基本成線性關(guān)系；降雨前期土體入滲率較大，雨水全部自地表入滲，此時(shí)地表含水率增長(zhǎng)速度較快并與降雨歷時(shí)呈線性正相關(guān)；當(dāng)降雨歷時(shí)進(jìn)一步增加時(shí)，雨水入滲率隨之降低，地表由于雨水入滲率降低到小于降雨強(qiáng)度而出現(xiàn)積水，此時(shí)地表土體達(dá)到飽和，降雨50 h后地表處含水率達(dá)到33.87%，接近其飽和含水率(34.19%)。同時(shí)，圖中還顯示：基質(zhì)吸力在降雨初期減小速度較快；當(dāng)降雨歷時(shí)進(jìn)一步增加時(shí)，基質(zhì)吸力減小速度逐漸緩和，最終趨于平穩(wěn)，持續(xù)降雨50 h后，地表處基質(zhì)吸力由初始的151.90 kPa減小到9.69 kPa。

圖3 地表基質(zhì)吸力及含水率隨降雨歷時(shí)變化曲線

圖4為降雨前土體初始含水率及孔隙水壓力分布的等值線云圖。初始土體含水率及孔隙水壓力沿深度方向呈線性分布，地表處土體較為干燥，含水率較低，基質(zhì)吸力最大。圖5～7分別為降雨強(qiáng)度600 mm/d工況下持續(xù)降雨10， 30， 50 h后土體含水率及孔隙水壓力分布的等值線云圖。圖中顯示：當(dāng)降雨持續(xù)時(shí)間增加時(shí)，浸潤(rùn)鋒線不斷向土體深處移動(dòng)，降雨50 h后浸潤(rùn)鋒線推進(jìn)到約15 m深位置處，在浸潤(rùn)鋒線的深度范圍內(nèi)，土體的含水率和孔隙水壓力都有所增大。相同降雨時(shí)間下，地表處含水率增加量最大，浸潤(rùn)鋒線以下土體含水率基本不變。

圖4 初始含水率與孔隙水壓力分布云圖(降雨強(qiáng)度=600 mm/d)

圖5 降雨10 h后含水率與孔隙水壓力分布云圖(降雨強(qiáng)度=600 mm/d)

圖6 降雨30 h后含水率與孔隙水壓力分布云圖(降雨強(qiáng)度=600 mm/d)

圖7 降雨50 h后含水率與孔隙水壓力分布云圖(降雨強(qiáng)度=600 mm/d)

3.2 拱頂沉降分析

圖8為不同降雨強(qiáng)度下，隧道拱頂沉降隨降雨持續(xù)時(shí)間的變化曲線。從圖8可以看出，降雨初期隧道拱頂沉降并無明顯變化；當(dāng)降雨持續(xù)一定時(shí)間后，拱頂沉降開始顯著增加。分析其原因：雨水入滲的影響深度在降雨初期較小，隧道圍巖強(qiáng)度及土中基質(zhì)吸力的減小并不明顯，圍巖仍具有較好的自支護(hù)能力。當(dāng)降雨歷時(shí)進(jìn)一步增加時(shí)，浸潤(rùn)鋒線不斷向土體深處移動(dòng)，即向拱頂發(fā)展。當(dāng)降雨持續(xù)時(shí)間達(dá)到一定程度時(shí)，浸潤(rùn)鋒線接近拱頂，在浸潤(rùn)鋒線范圍內(nèi)，圍巖中基質(zhì)吸力隨含水率的增大而顯著減小，進(jìn)而使土體強(qiáng)度顯著降低，極大地削弱了圍巖承載能力，導(dǎo)致拱頂沉降顯著增加。此外，從圖8可以看出：拱頂沉降隨降雨強(qiáng)度的增加而明顯增大，但是當(dāng)降雨強(qiáng)度增大到一定程度(600 mm/d)后，拱頂沉降對(duì)降雨強(qiáng)度變化的敏感性減弱。分析其原因：當(dāng)降雨強(qiáng)度較小時(shí)，雨水入滲率小于土體的飽和滲透系數(shù)，因此降雨過程中的雨水全部入滲，拱頂沉降基本隨降雨強(qiáng)度線性增加；當(dāng)降雨強(qiáng)度增大到某一值后，部分雨水來不及入滲，不久將會(huì)在地表形成積水，這部分積水將以地面徑流的方式排走，此時(shí)雨水入滲率與土體的飽和滲透系數(shù)相等，因此此時(shí)拱頂沉降與降雨強(qiáng)度關(guān)系不大。降雨前隧道拱頂沉降為14.59 mm，各降雨強(qiáng)度下，降雨結(jié)束后拱頂沉降依次增大至15.79，17.33，18.44，19.10，19.56 mm，相比于降雨前分別增加了8.22%，18.77%，26.39%，30.91%，34.06%。顯然，長(zhǎng)時(shí)間的強(qiáng)降雨對(duì)隧道圍巖變形和穩(wěn)定的影響較大。

圖8 拱頂沉降隨降雨歷時(shí)變化曲線

3.3 初期支護(hù)彎矩分析

圖9為不同降雨強(qiáng)度下隧道初期支護(hù)最大彎矩隨降雨歷時(shí)的變化曲線。其中，各工況下初期支護(hù)的最大彎矩均出現(xiàn)在邊墻與仰拱交界處。與拱頂沉降的變化規(guī)律相似，即在降雨初期，初期支護(hù)彎矩基本不變，降雨持續(xù)一定時(shí)間后，初期支護(hù)彎矩突然增大，最終趨于穩(wěn)定。初期支護(hù)彎矩隨降雨強(qiáng)度的增加而明顯增大，但是當(dāng)降雨強(qiáng)度增大到一定程度(600 mm/d)后，初期支護(hù)彎矩隨降雨強(qiáng)度變化不再明顯并逐漸趨于穩(wěn)定。降雨前初期支護(hù)最大彎矩為24.83 kN·m，各降雨強(qiáng)度下，降雨結(jié)束后，初期支護(hù)最大彎矩依次增大至25.83，27.89，29.37，29.95，30.24 kN·m，相比于降雨前分別增加了4.03%，10.31%，18.28%，20.62%，21.79%?？梢姡L(zhǎng)時(shí)間的強(qiáng)降雨對(duì)初期支護(hù)的彎矩影響很大。

圖9 初期支護(hù)最大彎矩隨降雨歷時(shí)變化曲線

3.4 錨桿軸力分析

圖10為不同降雨強(qiáng)度下錨桿最大軸力隨降雨歷時(shí)的變化曲線。同樣地，各工況下錨桿最大軸力值均出現(xiàn)在邊墻與仰拱交界處。從圖中可看出，在降雨初期，錨桿軸力基本保持不變，隨著降雨時(shí)間的增加，錨桿的最大軸力反而變小。此外，隨著降雨強(qiáng)度的增大，錨桿軸力同樣減小。分析其原因：錨桿主要是通過與周圍巖土體之間的黏結(jié)摩擦作用發(fā)揮功效，降雨導(dǎo)致錨桿錨固段土體軟化，對(duì)錨桿的粘結(jié)錨固作用減弱，錨桿軸力減小。

圖10 錨桿最大軸力隨降雨歷時(shí)變化曲線

4 結(jié) 論

本文通過采用FLAC3D軟件，編制相應(yīng)的FISH語(yǔ)言程序?qū)⒔涤耆霛B過程中基質(zhì)吸力及土體強(qiáng)度參數(shù)隨含水率的變化關(guān)系嵌入到計(jì)算過程中，對(duì)圍巖隧道在降雨入滲過程中發(fā)生的非飽和滲流進(jìn)行數(shù)值模擬，研究了不同降雨強(qiáng)度下隧道變形及支護(hù)受力的變化規(guī)律，得到以下主要結(jié)論：

(1)隨降雨歷時(shí)的增加，浸潤(rùn)鋒線不斷向土體深處移動(dòng)，在浸潤(rùn)鋒線的深度范圍內(nèi)土體的含水率逐漸增加，基質(zhì)吸力降低；降雨持續(xù)一定時(shí)間后，地表土體含水率接近飽和，基質(zhì)吸力趨于零。

(2)降雨初期，隧道拱頂沉降以及初期支護(hù)彎矩隨降雨歷時(shí)的變化并不明顯；當(dāng)降雨持續(xù)一定時(shí)間，浸潤(rùn)鋒線接近拱頂時(shí)，拱頂沉降及初期支護(hù)彎矩均顯著增加。降雨對(duì)錨桿軸力影響不大，當(dāng)降雨歷時(shí)增加時(shí)，圍巖有所軟化，對(duì)錨桿的粘結(jié)錨固效果減弱，錨桿軸力減小。

(3)當(dāng)降雨強(qiáng)度較小時(shí)，隧道拱頂沉降和支護(hù)受力隨降雨強(qiáng)度增加變化明顯；當(dāng)降雨強(qiáng)度增大到一定值后，此時(shí)雨水入滲率與土體飽和滲透系數(shù)相等，拱頂沉降和支護(hù)受力隨降雨強(qiáng)度的增加變化較小并逐漸趨于穩(wěn)定。