張?zhí)媒?/p>

(1.蘇交科集團(tuán)股份有限公司,江蘇 南京 210017; 2.江蘇省地下空間探測(cè)技術(shù)工程研究中心,江蘇 南京 210017)

0 引言

對(duì)于黃土隧道工程,在圍巖處于較低含水量或非飽和黃土?xí)r,隧道圍巖和隧道支護(hù)結(jié)構(gòu)穩(wěn)定性一般較好,一旦巖土體含水率增大,圍巖結(jié)構(gòu)自穩(wěn)能力則急劇下降,造成隧道支護(hù)結(jié)構(gòu)變形,更有甚者導(dǎo)致支護(hù)結(jié)構(gòu)破壞。近十幾年來,高含水率或飽和(可簡(jiǎn)稱富水)黃土造成的施工技術(shù)難題和運(yùn)營(yíng)階段的隧道結(jié)構(gòu)病害越來越多,給隧道施工設(shè)計(jì)和后期運(yùn)營(yíng)養(yǎng)護(hù)帶來很多挑戰(zhàn)。

目前國(guó)內(nèi)一些專家學(xué)者對(duì)高含水率黃土隧道進(jìn)行了理論研究和工程實(shí)踐,劉曉杰等[1]通過現(xiàn)場(chǎng)監(jiān)測(cè)對(duì)黃土隧道圍巖變形規(guī)律進(jìn)行了分析,并提出了采用拱腳鋼管樁與帷幕注漿來控制滲漏水及圍巖大變形控制的措施。來弘鵬等[2]通過數(shù)值模擬了富水黃土隧道上覆軟塑黃土層距隧道拱頂不同距離下隧道開挖卸載過程中圍巖水分遷移規(guī)律。梁慶國(guó)等[3]結(jié)合該富水黃土隧道圍巖變形特征現(xiàn)場(chǎng)監(jiān)測(cè)與對(duì)比分析,討論了經(jīng)帷幕注漿加固效果。張民慶等[4]對(duì)深埋富水黃土隧道地表深孔剛性袖閥管注漿試驗(yàn)研究。王釬[5]以某大斷面黃土隧道為工程背景,通過數(shù)值模擬和現(xiàn)場(chǎng)監(jiān)測(cè)的方式,研究了高含水率黃土隧道初期支護(hù)的力學(xué)特性。劉保成等[6]以賈家灣隧道富水段為研究背景,分析了圍巖變形的原因,提出了富水黃土隧道變形控制關(guān)鍵措施。劉志強(qiáng)等[7]針對(duì)寶蘭線高含水率黃土隧道,通過對(duì)不同施工情況下的黃土隧道進(jìn)行數(shù)值計(jì)算得出:早封閉對(duì)控制沉降變形極為有利。蘇春暉等[8]通過PLAXIS有限元軟件,對(duì)富水黃土隧道圍巖含水率的變化對(duì)隧道穩(wěn)定性影響進(jìn)行了研究。宋維龍等[9]針對(duì)某黃土隧道襯砌開裂病害,通過有限元軟件進(jìn)行數(shù)值模擬,分析了多種圍巖浸水工況下隧道結(jié)構(gòu)的變形規(guī)律和受力特征。陳鶴等[10]通過對(duì)富水黃土隧道的現(xiàn)場(chǎng)監(jiān)測(cè),分析了隧道應(yīng)力、位移的變化規(guī)律,發(fā)現(xiàn)初期支護(hù)產(chǎn)生一定的安全變形反而有利于隧道的施工安全。王新東[11]進(jìn)行了實(shí)驗(yàn)研究,得到了高含水率黃土隧道在不同含水率條件下的變形特征。周小賓[12]結(jié)合富水黃土張茅隧道,通過數(shù)值模擬和現(xiàn)場(chǎng)實(shí)測(cè)的方法,分析了黃土隧道的圍巖受力機(jī)理及變形規(guī)律。王學(xué)敬[13]基于進(jìn)口段淺埋偏壓及高含水率黃土隧道,對(duì)其施工變形機(jī)理進(jìn)行了研究。劉濤等[14]采用現(xiàn)場(chǎng)監(jiān)測(cè)的方式,對(duì)黃土塬區(qū)早勝三號(hào)隧道圍巖含水率進(jìn)行研究,結(jié)果表明:黃土隧道圍巖含水率呈現(xiàn)“增大-波動(dòng)-平穩(wěn)”的變化趨勢(shì)。劉俊平[15]通過現(xiàn)場(chǎng)調(diào)研,總結(jié)出了董志塬區(qū)地下水分布特征,并分析了地下水對(duì)黃土隧道工程的影響。

綜上所述,目前行業(yè)內(nèi)專家學(xué)者雖然對(duì)富水黃土隧道作了一定研究,但是針對(duì)黃土隧道圍巖軟化過程的研究成果較少,尤其是黃土隧道圍巖在水作用下軟化過程、軟化范圍及襯砌變化規(guī)律缺乏相應(yīng)的研究成果。因此,對(duì)黃土隧道圍巖軟化過程進(jìn)行精細(xì)化分析還有待近一步研究。本文以榆林子隧道為工程背景,通過數(shù)值模擬和現(xiàn)場(chǎng)調(diào)研,重點(diǎn)分析圍巖的不同軟化階段和軟化范圍對(duì)黃土隧道圍巖及襯砌的變形規(guī)律,以期為富水黃土隧道工程施工設(shè)計(jì)提供一定的參考和借鑒。

1 工程概況

2 數(shù)值模型建立

2.1 模型建立

針對(duì)黃土的工程地質(zhì)特性,采用有限元數(shù)值分析軟件ABAQUS,對(duì)于圍巖遇水軟化后圍巖強(qiáng)度降低的工況進(jìn)行數(shù)值建模分析,模擬隧道開挖施工過程,利用其溫度場(chǎng)實(shí)現(xiàn)施工完成后的軟化過程模擬;可以得到隧道開挖完成后產(chǎn)生的塑性區(qū),如圖2所示,塑性區(qū)最大影響范圍為3 m。因此,我們軟化范圍以隧道輪廓為初始邊界,向外擴(kuò)展取3 m,如圖3所示。隧道開挖后會(huì)產(chǎn)生應(yīng)力重分布,其影響范圍一般為隧道周圍3倍～5倍洞徑,因此,為避免邊界效應(yīng)的影響,隧道左右邊界和下邊界取3倍的洞徑,所選斷面隧道埋深按照實(shí)際工程取90 m,模型范圍取值為110 m(長(zhǎng))×150 m(寬),圍巖、初期支護(hù)、二次襯砌采用平面應(yīng)變單元CPE4模擬,計(jì)算中不考慮錨桿影響,將其認(rèn)為是對(duì)隧道開挖擾動(dòng)圍巖的補(bǔ)償。

隧道施工采用三臺(tái)階七步開挖,具體施工步驟為:1)上臺(tái)階弧形導(dǎo)坑開挖并支護(hù);2)左中臺(tái)階并支護(hù);3)右中臺(tái)階開挖并支護(hù);4)左下臺(tái)階開挖并支護(hù);5)右下臺(tái)階開挖并支護(hù);6)核心土開挖;7)仰拱開挖并支護(hù);8)仰拱填充層施作;9)二次襯砌施作。模型網(wǎng)格劃分和支護(hù)見圖4。

2.2 模型假定

1)假定土體為各向同性的均質(zhì)土,彈塑性體,本構(gòu)模型采用Mohr-Coulomb;2)初始應(yīng)力只考慮自重應(yīng)力;3)忽略時(shí)間效應(yīng),認(rèn)為每一步開挖和支護(hù)一步到位;4)模型不考慮滲流場(chǎng)的作用,將地下水的影響簡(jiǎn)化為對(duì)圍巖力學(xué)參數(shù)的折減作用。

2.3 計(jì)算參數(shù)

將隧道開挖后圍巖在水的浸泡作用下軟化過程分為三個(gè)階段,分別為初始軟化階段(K1)、中軟化階段(K2)及強(qiáng)軟化階段(K3),通過數(shù)值模擬著重分析圍巖軟化過程對(duì)隧道結(jié)構(gòu)的力學(xué)效應(yīng)的影響。根據(jù)圍巖的結(jié)構(gòu)特征,參照現(xiàn)場(chǎng)地勘資料和JTG 3370.1—2018公路隧道設(shè)計(jì)規(guī)范,給出隧道結(jié)構(gòu)和圍巖的各類物理力學(xué)參數(shù)指標(biāo)如表1,表2所示。

表1 隧道結(jié)構(gòu)力學(xué)參數(shù)

表2 圍巖力學(xué)參數(shù)

2.4 計(jì)算工況

分別考慮軟化階段和軟化范圍兩個(gè)因素,計(jì)算工況共有9種,具體的工況計(jì)算表如表3所示。

表3 計(jì)算工況表

3 結(jié)果分析

3.1 位移分析

3.1.1 圍巖變形

圍巖變形量影響著隧道開挖和支護(hù)是否及時(shí)可靠,為分析隧道開挖和軟化后的圍巖變形規(guī)律,提取拱頂和隧底位置的變形,圍巖豎向位移云圖如圖5所示。圍巖變形隨開挖分析步的變化曲線如圖6所示。圍巖軟化后不同軟化階段隨軟化范圍變化曲線如圖7所示。

由圖5可以看出,對(duì)比圍巖軟化前后的位移云圖,軟化后圍巖隆起和沉降數(shù)值均在變大,圍巖豎向位移變化范圍逐漸擴(kuò)大,其中圍巖隆起最大值由隧底處下移,說明軟化后圍巖的豎向變形及其影響范圍變化顯著。

由圖6可以看出,拱頂沉降和隧底隆起值隨開挖分析步呈現(xiàn)先快速增大后逐漸趨于穩(wěn)定,從具體開挖步來看,隨著上臺(tái)階的開挖拱頂沉降和隧底隆起變化較快,沉降值約為隆起值的2倍,此后隨著中臺(tái)階和下臺(tái)階的開挖和支護(hù),圍巖變形逐漸增大;在仰拱開挖時(shí),拱頂沉降已基本穩(wěn)定在107.39 mm,而隧底隆起有較小幅度的增長(zhǎng),在仰拱封閉后變形趨于穩(wěn)定,最終隆起值為104.05 mm。

從圖7可以看出,在隧道開挖完成后,圍巖軟化過程中,總體來看,隨著軟化范圍的變化,拱頂沉降逐漸增大,而隧底隆起值呈現(xiàn)減小的趨勢(shì)。在軟化階段K1圍巖變形變化小,基本與初始值相接近;軟化階段K2和軟化階段K3,在軟化范圍1 m時(shí)變化較小,在軟化范圍2 m后變化幅度較大;以軟化階段K3為例,拱頂沉降由開挖完成后的數(shù)值107 mm增大為152 mm,變化幅度達(dá)到42%,隧底隆起值由103 mm減小為58 mm,變化幅度為44%?？梢妵鷰r軟化對(duì)于隧道的變形影響很大,因此,在黃土隧道施工過程中應(yīng)及時(shí)排除地下水,避免圍巖進(jìn)一步劣化,對(duì)隧道的穩(wěn)定產(chǎn)生不良影響。

3.1.2 仰拱填充層隆起變形

仰拱填充層是隧道重要的組成部分,仰拱填充層是否穩(wěn)定直接影響隧道施工質(zhì)量安全和運(yùn)營(yíng)期的行車安全性及舒適性。為研究隧道圍巖軟化對(duì)仰拱填充層的影響,提取仰拱填充層中心單元E3,1/4位置單元E2以及拱腳單元E1的變形值,繪制其與軟化階段和軟化范圍的關(guān)系曲線如圖8所示。

從圖8可以看出,仰拱填充層中心、1/4位置以及拱腳處的隆起變形規(guī)律基本一致,都隨軟化階段和軟化范圍的增大不斷增大,在軟化階段K1,仰拱填充層單元E1,E2,E3的隆起變形的增大速率較快,在軟化范圍為3 m時(shí),單元E1,E2,E3的隆起變形都達(dá)到最大值,其中仰拱填充層中心隆起量最大;在軟化階段K2和K3,仰拱填充層單元E1,E2,E3的隆起變形的增大速率在軟化范圍為2 m后放緩,說明軟化范圍在2 m之后對(duì)仰拱填充層的影響較小,并逐漸趨于穩(wěn)定,說明仰拱填充層的隆起變形主要在塑性區(qū)內(nèi)變化。在三個(gè)軟化階段中的仰拱填充層隆起最大值均出現(xiàn)在仰拱填充層中心,隨著隆起增大,可能會(huì)導(dǎo)致仰拱填充層內(nèi)部脫空現(xiàn)象發(fā)生。因此,在高含水率黃土隧道設(shè)計(jì)中應(yīng)重點(diǎn)關(guān)注仰拱填充層中心位置變形,施工過程中密切監(jiān)測(cè)仰拱填充層中心的隆起變形值。

3.2 應(yīng)力分析

以軟化階段K3為例,提取不同軟化范圍下的襯砌最大主應(yīng)力云圖,如圖9所示。

從圖9可以看出,在軟化范圍0即未軟化時(shí),隧道結(jié)構(gòu)最大主應(yīng)力很小,結(jié)構(gòu)處于安全狀態(tài),隨著軟化范圍的增大應(yīng)力逐漸增大,最大應(yīng)力值主要集中在拱腳和仰拱填充層中部,在軟化深度3 m時(shí)達(dá)到1.55 MPa,可以看出圍巖遇水后對(duì)結(jié)構(gòu)的影響極為不利。

3.2.1 仰拱填充層應(yīng)力分析

由應(yīng)力云圖(如圖9所示)可以看出仰拱填充層應(yīng)力呈對(duì)稱分布,因此,為分析三個(gè)軟化階段下仰拱填充層單元隨軟化范圍的變化規(guī)律,可取結(jié)構(gòu)的一半進(jìn)行分析。通過提取仰拱填充層中心單元E3,1/4位置單元E2以及拱腳單元E1的最大和最小主應(yīng)力,繪制其與軟化階段和軟化范圍的關(guān)系曲線見圖10。

由圖10可以看出,軟化階段和軟化范圍的改變,對(duì)仰拱填充層的應(yīng)力變化有一定影響。分別對(duì)仰拱填充層單元E1,E2,E3位置的受力進(jìn)行分析,拱腳單元E1在和軟化階段K1和K2的最大主應(yīng)力均表現(xiàn)為受拉,隨著軟化范圍的增大而增大,但是增大的趨勢(shì)較小,在軟化階段K3單元E1的最大主應(yīng)力隨著軟化范圍的增大反而減小;單元E1在三個(gè)軟化階段的最小主應(yīng)力均表現(xiàn)為受壓,且隨軟化階段和軟化范圍的增大而不斷增大,總的來看,仰拱填充層拱腳位置的應(yīng)力變化受到圍巖軟化的影響較小。

仰拱1/4和隧底位置單元E2,E3的最大主應(yīng)力隨軟化階段和軟化范圍的變化較為一致,都隨軟化階段和軟化范圍的增大而增大,但其增大趨勢(shì)隨軟化階段逐漸減小,最大值出現(xiàn)在軟化階段K3軟化范圍為3 m時(shí),最大值位置為仰拱填充層中心,大小為1 555 kPa;單元E2,E3的最小主應(yīng)力隨軟化階段和軟化范圍的增大也不斷增大;總的來看,軟化范圍增大到3 m時(shí),圍巖軟化對(duì)仰拱填充層的受力影響逐漸減小,圍巖軟化對(duì)仰拱填充層中心影響最大,在設(shè)計(jì)施工中必須引起重視。

3.2.2 初期支護(hù)應(yīng)力分析

初期支護(hù)作為隧道的主要承載結(jié)構(gòu),通過研究不同軟化階段下隨軟化范圍變化的初期支護(hù)應(yīng)力來分析隧道結(jié)構(gòu)的受力變化。從上述云圖(如圖9所示)可看出襯砌應(yīng)力分布較為對(duì)稱,因此,取一半結(jié)構(gòu)提取拱頂單元E1、邊墻單元E2和拱腳單元E3的最大和最小主應(yīng)力進(jìn)行分析,如圖11所示。

由圖11可以看出,初期支護(hù)最大主應(yīng)力均為正,最大值出現(xiàn)在拱腳位置,數(shù)值為1 339 kPa。拱頂、拱墻、拱腳的最大主應(yīng)力隨軟化階段和軟化范圍的增大而增大,拱腳處應(yīng)力的增大較為明顯,而拱墻和拱頂?shù)膽?yīng)力變化較小,說明圍巖軟化對(duì)拱墻和拱頂受力的影響較小,而對(duì)拱腳受力的影響較大,外圈軟化后受到圍巖的擠壓作用,并且隨軟化范圍擴(kuò)大擠壓作用愈發(fā)明顯,導(dǎo)致拱腳的應(yīng)力明顯增大,受拉狀態(tài)不斷增強(qiáng)。拱墻處最小主應(yīng)力在軟化范圍增大到2 m后趨于穩(wěn)定,拱頂、拱腳的最小主應(yīng)力隨軟化范圍的擴(kuò)大逐漸增大?？偟膩砜?圍巖軟化會(huì)導(dǎo)致初期支護(hù)拱腳位置產(chǎn)生較大的拉應(yīng)力,對(duì)于拱腳結(jié)構(gòu)穩(wěn)定產(chǎn)生不利影響。

4 結(jié)論

本文以實(shí)際工程為背景,通過現(xiàn)場(chǎng)調(diào)研和數(shù)值摸擬,設(shè)置不同軟化階段和軟化范圍,分析了隧道施工過程中圍巖軟化導(dǎo)致的結(jié)構(gòu)變形和應(yīng)力變化,得出了如下幾個(gè)結(jié)論:

1)拱頂沉降隨軟化范圍的擴(kuò)大呈現(xiàn)增大趨勢(shì),而隧底隆起逐漸減小。

2)仰拱填充層隆起變形隨軟化范圍擴(kuò)大逐漸增大,當(dāng)軟化范圍超過隧道圍巖塑性區(qū)后逐漸趨于穩(wěn)定;由于仰拱填充層的隆起可能導(dǎo)致隧道仰拱內(nèi)部產(chǎn)生脫空現(xiàn)象。

3)圍巖軟化對(duì)拱墻和拱頂受力的影響較小,而在初期支護(hù)拱腳位置產(chǎn)生較大的拉應(yīng)力,導(dǎo)致拱腳穩(wěn)定性受到較大影響。

4)由于地下水滲流后在拱腳處匯集,使得拱腳圍巖快速軟化,導(dǎo)致拱腳部位承載能力下降,且在施工擾動(dòng)作用下,引起附加沉降,結(jié)構(gòu)應(yīng)力增大,如果不采取措施會(huì)造成支護(hù)結(jié)構(gòu)開裂、侵限和塌方等災(zāi)害發(fā)生,因此施工現(xiàn)場(chǎng)應(yīng)加強(qiáng)隧道防排水體系的設(shè)計(jì),重視隧道基底的加固處理。

5)通過研究發(fā)現(xiàn)圍巖軟化對(duì)于隧道的影響主要在塑性區(qū)范圍內(nèi),在富水黃土隧道施工和后期運(yùn)營(yíng)過程中,應(yīng)密切注意塑性區(qū)內(nèi)的地下水賦存狀態(tài),避免圍巖進(jìn)入強(qiáng)軟化階段。