蒙曉新 張凱

作者簡介：蒙曉新（1983—），工程師，主要從事交通基礎(chǔ)設(shè)施施工建設(shè)管理、項目投資運(yùn)營管理等工作。

為研究膨脹巖地區(qū)隧道圍巖力學(xué)性能及破壞模式，文章依托某高速公路隧道實(shí)際工程，通過室內(nèi)巖石力學(xué)試驗分析了膨脹巖的力學(xué)特性和膨脹特性，然后通過FLAC3D軟件對該隧道進(jìn)行建模與計算，分析膨脹性圍巖對隧道變形的影響規(guī)律，得到如下結(jié)論：（1）膨脹巖試樣在三軸壓縮條件下的應(yīng)力-應(yīng)變曲線主要經(jīng)歷了彈性、屈服和破壞三個階段，三軸抗壓強(qiáng)度隨著圍壓的增大而增大；（2）試樣初始含水率較低時，其膨脹系數(shù)更大，前期變形速率更快，進(jìn)入變形減速階段的時間更晚，達(dá)到平衡狀態(tài)所需的時間更長，最終的膨脹系數(shù)數(shù)值更大；（3）膨脹性圍巖的膨脹作用對隧道拱底隆起變形的影響最大，圍巖變形程度隨圍巖初始含水率增大而逐漸減小，圍巖的初始含水率＜15%時，圍巖的變形程度較大，在進(jìn)行支護(hù)時應(yīng)重點(diǎn)考慮膨脹圍巖膨脹力的影響。

膨脹巖隧道；力學(xué)性能；圍巖變形規(guī)律；數(shù)值計算；FLAC3D軟件

U456.3+1A381344

0?引言

隨著我國經(jīng)濟(jì)的快速發(fā)展和交通強(qiáng)國戰(zhàn)略的不斷推進(jìn)，我國現(xiàn)已成為世界上隧道數(shù)量最多、規(guī)模最大、范圍最廣、工程條件最復(fù)雜、發(fā)展最迅速的國家［1］。由于我國西南地區(qū)山區(qū)面積大，地質(zhì)情況復(fù)雜，在建設(shè)隧道的過程中往往需要穿越諸多不良地質(zhì)區(qū)域，如軟弱圍巖、膨脹圍巖、富水地層等，隧道的施工安全面臨巨大挑戰(zhàn)。在膨脹巖地區(qū)開挖隧道時，圍巖卸荷及膨脹作用導(dǎo)致的應(yīng)力重分布會使隧道結(jié)構(gòu)產(chǎn)生劇烈擠壓，隧道圍巖的膨脹和收縮會加劇隧道內(nèi)部的滲水和襯砌剝落等現(xiàn)象，對于隧道的建設(shè)產(chǎn)生極大危害。對于膨脹巖的力學(xué)特性以及穿越膨脹巖地層的隧道圍巖穩(wěn)定性的研究已有了較多成果。王軍等［2］對不同含水率的紅砂巖進(jìn)行了室內(nèi)直剪試驗，得到了膨脹巖的粘聚力和內(nèi)摩擦角與含水率之間的對數(shù)關(guān)系。趙二平等［3］對南水北調(diào)中線工程的膨脹巖開展三軸循環(huán)加卸荷試驗，試驗結(jié)果發(fā)現(xiàn)可用負(fù)指數(shù)函數(shù)來表征膨脹巖的應(yīng)變增加量與循環(huán)次數(shù)之間的關(guān)系。曾志雄等［4］對不同初始含水率的延吉膨脹巖壓實(shí)試樣開展了3種水力路徑的一維壓縮試驗，研究了不同水力路徑下的膨脹巖變形特性。郭永春等［5］研制了膨脹巖三軸膨脹力的測試裝置，通過先放試樣后施加約束的試驗方式克服了巖石與容器壁有釋放變形空隙的局限性，解決了現(xiàn)有規(guī)范中對膨脹巖的膨脹力測試結(jié)果偏小的問題。楊軍平等［6］利用南寧地區(qū)的膨脹土進(jìn)行了膨脹圍巖的隧道物理模型，并對模型進(jìn)行了五次干濕循環(huán)，研究了干濕循環(huán)條件下膨脹圍巖的脹縮演變規(guī)律，分析了干濕循環(huán)作用對膨脹圍巖的劣化作用。趙凡等［7］依托某穿越石膏質(zhì)巖地層的隧道實(shí)際工程，基于隧道典型斷面病害的監(jiān)測數(shù)據(jù)，利用FINAL軟件對圍巖的膨脹參數(shù)反演，發(fā)現(xiàn)隧道圍巖的最大膨脹半徑與膨脹系數(shù)為正相關(guān)的關(guān)系。曾仲毅等［8］推導(dǎo)出熱傳導(dǎo)膨脹模擬膨脹巖增濕膨脹的替代方程，利用FLAC 3D軟件的熱-力耦合模塊對某膨脹圍巖隧道進(jìn)行數(shù)值仿真分析，得出支護(hù)結(jié)構(gòu)的受力變形特性隨含水率及膨脹力的變化規(guī)律。

本文依托某膨脹性圍巖地區(qū)的高速公路隧道實(shí)際工程，通過室內(nèi)巖石力學(xué)試驗分析了膨脹巖的力學(xué)特性和膨脹特性，然后通過FLAC 3D軟件對該隧道進(jìn)行建模與計算，分析膨脹性圍巖對隧道變形的影響規(guī)律。

1?膨脹巖力學(xué)特性試驗研究

1.1?工程概況

本文依托某高速公路隧道實(shí)際工程，對膨脹巖的力學(xué)特性及膨脹巖地區(qū)隧道圍巖的變形規(guī)律進(jìn)行研究。隧道全長2 505 m，所處的地貌類型為中低山區(qū)，最大埋深約為94.3 m，隧道圍巖的巖性主要是塊狀構(gòu)造的弱風(fēng)化安山巖，巖體完整，節(jié)理裂隙不發(fā)育，因此在勘察階段判定其為Ⅱ級圍巖，通過全斷面法進(jìn)行隧道開挖。但在施工過程中發(fā)現(xiàn)部分開挖段出現(xiàn)局部初期支護(hù)噴射混凝土開裂和滲水現(xiàn)象，進(jìn)一步調(diào)查發(fā)現(xiàn)隧道穿越安山巖和凝灰質(zhì)礫巖的接觸帶，接觸帶位置處的圍巖較為破碎且?guī)r性混雜，通過現(xiàn)場膨脹性試驗結(jié)果判定接觸帶位置處的圍巖具有弱膨脹性，屬于膨脹性巖體。

1.2?膨脹巖三軸壓縮試驗

由于已有研究發(fā)現(xiàn)，原狀膨脹巖的物理力學(xué)性質(zhì)與實(shí)驗室制得的重塑膨脹巖間存在較大區(qū)別，因此本試驗通過對原狀膨脹巖巖塊進(jìn)行鉆孔取芯和切割機(jī)切割的方法得到高度100 mm，直徑為50 mm的標(biāo)準(zhǔn)試樣。試驗設(shè)備采用高溫高壓巖石三軸儀，分別在圍壓σ3=2.0 MPa、3.0 MPa、4.0 MPa條件下進(jìn)行三軸壓縮試驗。

對膨脹巖的三軸壓縮試驗結(jié)果進(jìn)行分析，得到不同圍壓條件下三軸壓縮試驗的應(yīng)力-應(yīng)變曲線，如圖1所示。試樣在三軸壓縮條件下的應(yīng)力-應(yīng)變曲線主要經(jīng)歷了3個階段：第一個階段為彈性階段，應(yīng)力-應(yīng)變曲線基本為一條直線，即在該階段應(yīng)力和應(yīng)變之間為正比關(guān)系；第二個階段為屈服階段，試樣內(nèi)部的裂隙數(shù)量逐漸增加，隨著應(yīng)力的增大，應(yīng)變增大的速率開始增大；第三個階段為破壞階段，試樣內(nèi)部的裂隙互相貫通，發(fā)生剪切破壞，此時應(yīng)力達(dá)到峰值并出現(xiàn)明顯下降，但仍具有一定的殘余強(qiáng)度。膨脹巖的三軸抗壓強(qiáng)度隨著圍壓的增大而增大，當(dāng)圍壓分別為2.0 MPa、3.0 MPa、4.0 MPa時，試樣三軸抗壓強(qiáng)度值分別為25.8 MPa、29.1 MPa、35.5 MPa，即圍壓由2.0 MPa增大至4.0 MPa，試樣三軸壓縮抗壓強(qiáng)度相對于圍壓為2 MPa時分別增大了12.8%（3 MPa）和37.6%（4 MPa）。

根據(jù)《工程巖體試驗方法標(biāo)準(zhǔn)》（GB/T50266-201）中提供的計算巖石各力學(xué)參數(shù)的方法計算巖石的彈性模量、泊松比、粘聚力和內(nèi)摩擦角，如表1所示。

1.3?側(cè)向約束無荷膨脹率試驗

巖石的側(cè)向約束膨脹率是指巖石標(biāo)準(zhǔn)試樣在有側(cè)限條件下浸水后軸向變形與試樣原高度的比值。鑒于膨脹巖具有遇水泥化的特性，取芯具有很高的難度，且重塑后的巖樣內(nèi)部的親水性礦物的成分與含量均與原狀巖樣一致，重塑巖樣可以很大程度表征原狀巖石的膨脹特性，故本次采用61.8 mm×20 mm的重塑巖樣進(jìn)行膨脹巖的側(cè)向約束無荷膨脹率試驗，試驗在固結(jié)儀中進(jìn)行，設(shè)置初始含水率分別為5%、10%、15%、20%、25%。經(jīng)過分級加水后，發(fā)現(xiàn)當(dāng)含水率達(dá)到35%后繼續(xù)加水，水無法進(jìn)入試樣內(nèi)部，即該膨脹巖的飽和含水率為35%。

對試驗結(jié)果進(jìn)行分析，得到膨脹巖的側(cè)向約束無荷膨脹率試驗結(jié)果如表2所示，繪制不同初始含水量的膨脹巖試樣的膨脹曲線如圖2所示。

如圖2所示，隨著試樣初始含水率升高，其膨脹系數(shù)減小，前期變形速率減小，進(jìn)入減速變形階段的時間提前，達(dá)到平衡狀態(tài)所需的時間縮短，最終膨脹系數(shù)的數(shù)值減小，膨脹率增大速率出現(xiàn)由快變慢的拐點(diǎn)時對應(yīng)的相對含水率增大。不同初始含水率的試樣吸水膨脹后具有不同的膨脹系數(shù)，即試樣的膨脹系數(shù)與其初始含水率有關(guān)，二者間大致呈線性關(guān)系，因此采用線性函數(shù)對試樣的膨脹系數(shù)和初始含水率之間的關(guān)系進(jìn)行擬合，得到擬合函數(shù)式（1）。

可根據(jù)式（1）得到該類膨脹巖在不同初始含水率條件下吸水膨脹達(dá)到飽和的膨脹系數(shù)。

2?考慮圍巖膨脹性的隧道數(shù)值模擬研究

2.1?數(shù)值模型的建立

本文通過FLAC 3D軟件進(jìn)行數(shù)值模型的建立與計算，利用FLAC 3D軟件的熱-力耦合模塊模擬膨脹圍巖的膨脹效應(yīng)。依據(jù)實(shí)際隧道斷面建立數(shù)值模型，其中隧道洞徑r=14.0 m。為減小由邊界效應(yīng)引起的計算誤差，選取的數(shù)值模型的尺寸為140 m×70 m×140 m（長×寬×高），其中x軸為水平方向，y軸為隧道軸線方向，z軸為豎直方向。選取圍巖松動范圍作為計算的膨脹范圍，根據(jù)董方庭等［9］提供的經(jīng)驗公式計算得本文所研究的隧道的松動圈厚度約為1.22 m，故對開挖后1.22 m厚度范圍內(nèi)的圍巖施加溫度場模擬圍巖膨脹作用，建立數(shù)值模型如圖3所示。為減少邊界效應(yīng)的影響，取隧道縱向的中間斷面（Y=35.0 m）作為監(jiān)測面，監(jiān)測點(diǎn)設(shè)置如圖4所示，其中拱頂、拱底監(jiān)測圍巖的豎向位移，其余測點(diǎn)監(jiān)測圍巖水平位移。

邊界條件設(shè)置為：約束模型四周的法向位移，約束模型底面的水平和法向位移，頂部邊界考慮埋深的作用，施加方向向下的等效埋深94 m的巖土體重力的應(yīng)力邊界。模型的初始地應(yīng)力場僅考慮為重力場，在模型完成初始地應(yīng)力平衡后，執(zhí)行隧道開挖的命令。采用全斷面進(jìn)行隧道開挖，循環(huán)進(jìn)尺為2.0 m，初期支護(hù)滯后掌子面2.0 m。

本文選用FLAC 3D軟件中的應(yīng)變軟化摩爾-庫倫模型作為巖土體的本構(gòu)模型，選用彈塑性本構(gòu)模型作為初期支護(hù)的本構(gòu)模型，數(shù)值模型的計算參數(shù)如表3所示。

2.2?計算結(jié)果分析

根據(jù)數(shù)值計算結(jié)果，得到僅考慮開挖條件和考慮開挖和圍巖膨脹共同作用兩種工況下監(jiān)測斷面各監(jiān)測點(diǎn)的最終變形值，如圖5所示。

如圖5所示，膨脹性圍巖的膨脹作用對隧道變形具有不利影響，其中對拱底隆起變形的影響最大，隧道斷面產(chǎn)生向內(nèi)擠壓的變形，尤其對未封底的隧道，圍巖膨脹導(dǎo)致的隧道拱底隆起變形更為嚴(yán)重，隧道拱頂及兩側(cè)圍巖的膨脹壓力均向拱底傳遞，使隧道拱底位置在強(qiáng)烈的膨脹作用下首先失穩(wěn)，拱底位置處的圍巖產(chǎn)生流動現(xiàn)象，向隧道內(nèi)部移動，隧道拱頂及兩側(cè)的圍巖向拱底流動，進(jìn)一步加劇拱底的隆起變形。

2.3?初始含水率對膨脹巖隧道影響性分析

根據(jù)膨脹巖的側(cè)向約束無荷膨脹率試驗結(jié)果可知，初始含水率對膨脹巖的膨脹性具有很大影響，上一節(jié)中分析了初始含水率為5%的膨脹巖在隧道開挖過程中的變形，發(fā)現(xiàn)膨脹巖的膨脹作用對隧道的穩(wěn)定性會產(chǎn)生很大的影響，故本節(jié)對初始含水率分別為5%、10%、15%、20%、25%的膨脹圍巖5種工況進(jìn)行數(shù)值計算，分析初始含水率對膨脹圍巖隧道開挖的影響規(guī)律。根據(jù)數(shù)值計算結(jié)果，得到不同膨脹圍巖初始含水率下的隧道開挖面圍巖洞周位移如圖6所示。

由圖6可知，隧道圍巖的變形程度大致呈現(xiàn)隨初始含水率增大逐漸減小的規(guī)律，分析原因為：當(dāng)圍巖初始含水率較低時，圍巖中存在較多孔隙，遇水后的吸水能力更強(qiáng)，吸收的水分更多，產(chǎn)生更大的膨脹壓力。當(dāng)圍巖的初始含水率為5%和10%時，隧道圍巖的變形程度較大且二者數(shù)值較為接近；當(dāng)圍巖的初始含水率為15%時，隧道圍巖的變形程度出現(xiàn)較大程度的降低；當(dāng)圍巖的初始含水率為20%和25%時，隧道圍巖的變形程度再次出現(xiàn)一定程度的降低且二者數(shù)值較為接近。即當(dāng)圍巖的初始含水率＜15%時，圍巖的變形程度較大，在進(jìn)行支護(hù)時應(yīng)重點(diǎn)考慮膨脹圍巖膨脹力的影響，避免后期產(chǎn)生較大的膨脹變形，導(dǎo)致隧道結(jié)構(gòu)發(fā)生破壞。

3?結(jié)語

（1）膨脹巖試樣在三軸壓縮條件下的應(yīng)力-應(yīng)變曲線主要經(jīng)歷了彈性、屈服和破壞3個階段。三軸抗壓強(qiáng)度隨著圍壓的增大而增大，圍壓由2.0 MPa增大至4.0 MPa，試樣三軸壓縮抗壓強(qiáng)度相對于圍壓為2 MPa時分別增大了12.8%（3 MPa）和37.6%（4 MPa）。

（2）試樣初始含水率較低時，其膨脹系數(shù)更大，前期變形速率更快，進(jìn)入變形減速階段的時間更晚，達(dá)到平衡狀態(tài)所需的時間更長，最終的膨脹系數(shù)數(shù)值更大，膨脹率增大速率出現(xiàn)由快變慢的拐點(diǎn)時對應(yīng)的相對含水率越小。

（3）膨脹性圍巖的膨脹作用對隧道變形具有不利影響，其中對拱底隆起變形的影響最大。隧道圍巖的變形程度大致呈現(xiàn)隨初始含水率增大逐漸減小的規(guī)律，其原因為初始含水率較低的圍巖中存在較多孔隙，遇水后的吸水能力更強(qiáng)，吸收的水分更多，產(chǎn)生更大的膨脹壓力。當(dāng)圍巖的初始含水率＜15%時，圍巖的變形程度較大，在進(jìn)行支護(hù)時應(yīng)重點(diǎn)考慮膨脹圍巖膨脹力的影響。

參考文獻(xiàn)

［1］王夢恕.中國盾構(gòu)和掘進(jìn)機(jī)隧道技術(shù)現(xiàn)狀、存在的問題及發(fā)展思路［J］.隧道建設(shè)，2014，34（3）：179-187.

［2］王?軍，何?淼，汪中衛(wèi).膨脹砂巖的抗剪強(qiáng)度與含水量的關(guān)系［J］.土木工程學(xué)報，2006（1）：98-102.

［3］趙二平，李建林，王瑞紅.不同應(yīng)力路徑下膨脹巖力學(xué)特性試驗研究［J］.人民長江，2014，45（3）：83-86.

［4］曾志雄，孔令偉，黃玨皓，等.延吉膨脹巖變形特性的水-力路徑效應(yīng)與機(jī)制淺析［J］.巖土工程學(xué)報，2016，38（S2）：222-227.

［5］郭永春，趙峰先，閆圣龍，等.紅層泥巖三軸膨脹力的試驗研究［J］.水文地質(zhì)工程地質(zhì)，2022，49（3）：87-93.

［6］楊軍平，王沾義，周立新，等.干濕交替條件下膨脹圍巖脹縮規(guī)律的試驗研究［J］.鐵道科學(xué)與工程學(xué)報，2017，14（1）：117-125.

［7］趙凡，高曉靜，楊香英，等.膨脹巖隧道圍巖膨脹參數(shù)反演分析［J］.現(xiàn)代隧道技術(shù)，2021，58（S1）：20-28.

［8］曾仲毅，徐幫樹，胡世權(quán)，等.增濕條件下膨脹土隧道襯砌破壞數(shù)值分析［J］.巖土力學(xué)，2014，35（3）：871-880.

［9］董方庭，宋宏偉，郭志宏，等.巷道圍巖松動圈支護(hù)理論［J］.煤炭學(xué)報，1994（1）：21-32.