王曉星，張昱輝，魯 彪

(1. 衢州市交通設(shè)計(jì)有限公司，浙江 衢州 324000； 2.浙江數(shù)智交院科技股份有限公司，浙江 杭州 310030)

0 引言

隨著我國高速公路網(wǎng)中車流量的逐年攀升，早年建成的高速公路早已處于超負(fù)荷運(yùn)轉(zhuǎn)狀態(tài)，因此高速公路改擴(kuò)建越來越受到人們的重視。隧道作為高速公路改擴(kuò)建工程的關(guān)鍵節(jié)點(diǎn)，施工過程中經(jīng)常會(huì)遇到臨近施工、擴(kuò)挖施工的情況[1-2]，甚至有時(shí)為了滿足國家規(guī)范中對道路安全與通行量的要求需對隧道進(jìn)行整條挖除，因此研究如何在保證既有臨近隧道結(jié)構(gòu)安全性的前提下圓滿完成改擴(kuò)建任務(wù)就顯得尤為重要。

針對施工活動(dòng)中產(chǎn)生的震動(dòng)對臨近隧道及構(gòu)筑物造成的影響，眾多學(xué)者進(jìn)行了廣泛而深入的研究。方智淳、張暻等[3-4]對采石場爆破震動(dòng)對臨近隧道的影響進(jìn)行分析，評估了采礦隧道結(jié)構(gòu)的安全性。胡剛、唐浩、孫艷軍等[5-7]基于數(shù)值模擬研究了爆破施工對臨近隧道的影響，得到了既有隧道的動(dòng)力響應(yīng)特征。丁祖德等[8]通過數(shù)值模擬分別比較了基坑開挖時(shí)采用不同保護(hù)措施附近既有隧道附加變形和附加內(nèi)力的大小、分布形態(tài)及其變形控制效果。高宏偉、朱蕾、黃茂松等[9-11]對基坑開挖時(shí)臨近及下臥隧道的受影響程度與范圍進(jìn)行分析，得到不同凈距施工時(shí)既有隧道的應(yīng)力與變形特征。還有學(xué)者對臨近隧道爆破施工的控制參數(shù)、爆破措施進(jìn)行了研究[12-16]。針對爆破影響問題，前人已做了廣泛的研究，但目前對路基挖除對另一幅運(yùn)營隧道影響的研究還比較少。本研究以杭金衢高速樊村隧道擴(kuò)建工程為背景，利用Midas GTS NX有限元軟件研究不同方式挖除路基時(shí)對既有隧道的影響，以期為工程的順利實(shí)施提供依據(jù)。

1 工程概況

如圖1所示，既有樊村隧道1號(hào)隧道出口距離五里樞紐距離為109 m，不滿足《公路立體交叉設(shè)計(jì)細(xì)則》(JTG/T D21—2014)單向四車道隧道出口與前方互通式立交最小距離為800 m的要求，故改建工程擬保留既有右線，在既有右線右側(cè)新建兩車道右線，將既有左線隧道整體挖除改建為整體四車道路基。

圖1 杭金衢高速樊村隧道擴(kuò)建工程Fig.1 Fancun tunnel expansion project of Hangzhou-Jinhua-Quzhou expressway

擬擴(kuò)建樊村隧道總長度628 m，洞頂最大埋深約46 m，隧道Ⅲ級(jí)圍巖占47.8%，V級(jí)和IV級(jí)圍巖占52.2%，主要穿越石英砂巖及含礫石英砂巖，局部為泥質(zhì)粉砂巖，圍巖整體性較差。隧道內(nèi)輪廓斷面為單心圓，支護(hù)形式采用錨噴支護(hù)+復(fù)合模注混凝土襯砌，內(nèi)夾防排水層，其中左線隧道ZK419+890處邊坡離既有右線隧道襯砌最近、圍巖最差，隧道間凈距為48.53 m，邊坡擬開挖線與隧道最近距離為23.58 m，隧道埋深31.02 m，是影響隧道路基開挖的最不利斷面，因此選取該斷面作分析模型，研究其在機(jī)械和爆破開挖方式下臨近隧道及其巖土體的變形和受力特性，以保證工程順利實(shí)施。

2 有限元模型及參數(shù)

為消除模型尺寸效應(yīng)引起的誤差，結(jié)合斷面ZK419+890的地質(zhì)狀況，選取斷面ZK419+890的計(jì)算模型尺寸為長×寬×高=190 m×60 m×89 m，隧道尺寸按1∶1選取，計(jì)算網(wǎng)格圖見圖2。

圖2 ZK419+890斷面計(jì)算模型網(wǎng)格圖Fig.2 Calculation model grid of section ZK419+890

選擇莫爾-庫倫理想彈塑性模型作為該計(jì)算的本構(gòu)模型，結(jié)合《公路隧道設(shè)計(jì)規(guī)范》(JTG 3370.1—2018)和樊村隧道地勘資料，隧道圍巖及二襯模型參數(shù)可按表1選取。

表1 本構(gòu)模型的主要參數(shù)Tab.1 Main parameters of constitutive model

3 機(jī)械開挖計(jì)算結(jié)果及分析

機(jī)械開挖采用從上至下分3步挖除路基，結(jié)合有限元分析軟件，按照開挖順序分析每層開挖完成后邊坡巖土體和隧道變形與受力情況。在土壓力作用下，各工況下邊坡巖土體的變形情況如圖3所示?？梢钥闯?，隨著開挖深度不斷加深，邊坡所受卸載效應(yīng)逐漸加重，對既有臨近隧道影響加大。

圖3 不同開挖階段邊坡巖土體位移云圖(單位：mm)Fig.3 Nephograms of displacement of slope rock and soil at different excavation stages(unit：mm)

第1層土方開挖完成后，卸載效應(yīng)主要集中在開挖底面，開挖區(qū)底部土體有明顯的隆起趨勢，同時(shí)受左側(cè)陡坡土壓力作用，巖體沿水平方向向開挖臨空面形成較大位移。該階段巖土體以豎直隆起為主，最大隆起高度8 mm，而既有隧道位于開挖面右下方較遠(yuǎn)距離，所受影響較小。與第1次土方開挖相比，第2次土方開挖完成后左側(cè)陡坡臨空面上的土體位移加大，卸載效應(yīng)在開挖面的水平方向更加明顯。而隨著開挖深度的加深，水平卸載效應(yīng)對既有隧道的影響也在不斷加大，此時(shí)隧道已經(jīng)產(chǎn)生較為明顯的位移。當(dāng)?shù)?層土方開挖完成后，卸載效應(yīng)達(dá)到最大，其中尤以左側(cè)邊坡最為明顯，此時(shí)豎直方向土體隆起位移不大，但邊坡在水平方向上有顯著的向開挖面移動(dòng)的趨勢，帶動(dòng)右側(cè)既有隧道向邊坡移動(dòng)方向產(chǎn)生較大位移。在3次開挖工況中，以第2次巖土體位移變化較快。

隨著土體的開挖，土體的卸載效應(yīng)會(huì)逐漸作用到臨近隧道襯砌結(jié)構(gòu)上，使隧道支護(hù)結(jié)構(gòu)的內(nèi)力發(fā)生重分布，同時(shí)會(huì)向卸載方向產(chǎn)生一定的位移。而在周圍巖體的約束下，隧道襯砌變形趨勢與巖土體變形趨勢基本保持一致，變形位置主要集中在左拱肩部位，其變形量隨著開挖量增大不斷增加。經(jīng)過計(jì)算，3層路基開挖完畢時(shí)襯砌最大位移量約4 mm，隨著后續(xù)施工活動(dòng)的進(jìn)行存在失穩(wěn)的可能，應(yīng)做好襯砌監(jiān)測與邊坡加固工作。

在支護(hù)結(jié)構(gòu)封閉成環(huán)以后，其整體剛度大大提升，在環(huán)向剛度效應(yīng)作用下支護(hù)結(jié)構(gòu)向卸載方向產(chǎn)生的位移會(huì)受到鋼拱架與襯砌的約束，從而使得非卸載方向的支護(hù)結(jié)構(gòu)內(nèi)部產(chǎn)生拉應(yīng)力，但隧道上覆土層的重力作用會(huì)使拉應(yīng)力有所折減，在兩者共同作用下，支護(hù)結(jié)構(gòu)內(nèi)力主要集中在拱頂、拱肩位置，隨著卸載效應(yīng)的不斷完成，應(yīng)力會(huì)在拱頂達(dá)到最大值約34.3 kN·m，使得該位置支護(hù)結(jié)構(gòu)存在破壞的風(fēng)險(xiǎn)，因此在利用機(jī)械開挖法進(jìn)行施工時(shí)，隧道左拱肩位移和拱頂內(nèi)側(cè)的內(nèi)力應(yīng)當(dāng)重點(diǎn)加強(qiáng)監(jiān)測。

4 爆破開挖計(jì)算結(jié)果分析

4.1 爆破荷載

采用面動(dòng)力荷載模擬爆破荷載，所作假定如下:(1)爆破荷載以壓力形式均布作用在坡面上，方向垂直于坡面。(2)根據(jù)計(jì)算和爆破實(shí)際情況，爆炸沖擊波宜采用指數(shù)型荷載衰減模型，見圖3。

在不耦合裝藥情況下，爆轟波會(huì)首先壓縮裝藥與藥室壁之間間隙內(nèi)的空氣形成空氣沖擊波，而后再由空氣沖擊波作用于藥室壁，對藥室壁加載。假定爆轟產(chǎn)物在間隙內(nèi)的膨脹為絕熱膨脹，忽略間隙內(nèi)空氣，則爆轟產(chǎn)物撞擊藥室壁前的炮孔內(nèi)壓力即入射壓力Pi為[17]：

(1)

式中,Pm為爆轟產(chǎn)物開始膨脹時(shí)的平均壓力；Vc為炸藥體積；Vb為藥室體積；ρ0為炸藥密度；D1為爆轟波傳播速度。當(dāng)爆轟產(chǎn)物撞擊藥室壁時(shí)，壓力將明顯增大，增大倍數(shù)n=8～11。取Vc/Vb=1/3，n=10，則可得爆轟峰值壓力Pi=0.427 GPa。邊坡開挖爆破按深孔爆破考慮，孔距為5 m，則可等效深孔爆破的爆破孔壓等效為面壓力Pm=0.083 GPa。

對于爆破動(dòng)力分析中的孔壁壓力衰減，文獻(xiàn)[18]給出了相對較為精確的指數(shù)型荷載衰減模式，即：

Pt=4Pm(e-0.707ωt-e-1.414ωt)，

(2)

(3)

式中,w為波形參數(shù)；vp為介質(zhì)(巖體)縱波波速；r0為藥包半徑；t為爆破時(shí)程。

在該模型作用下，爆破孔壓力時(shí)程曲線見圖4。

圖4 爆破孔壓力時(shí)程曲線Fig.4 Time history curve of blast hole pressure

4.2 邊界條件及阻尼

采用爆破施工，爆破產(chǎn)生的振動(dòng)波會(huì)使固定邊界反復(fù)反彈，故計(jì)算時(shí)采用黏彈性邊界比較合理。在MIDAS/GTS/NX中，通過曲面彈簧定義彈性邊界，彈簧系數(shù)可根據(jù)道路設(shè)計(jì)規(guī)范地基反力系數(shù)計(jì)算，如式(2)、式(3)所示。

豎直彈簧系數(shù)為：

(4)

水平彈簧系數(shù)為：

(5)

材料的阻尼特性會(huì)逐漸消耗施工活動(dòng)震動(dòng)所產(chǎn)生的能量，在體系的運(yùn)動(dòng)方程中，阻尼力與慣性力與抗力相比要小得多，可以近似地計(jì)算阻尼矩陣。在MIDAS/GTX/NX時(shí)域分析中，常采用的是瑞利阻尼。瑞利阻尼與頻率有關(guān)，其特征主要由臨界阻尼比和最小中心頻率決定，在本研究中取臨界阻尼比xmin=3.5%，最小中心頻率可用模型自振頻率f代替，取f=18.26 Hz。

4.3 計(jì)算結(jié)果分析

前面的計(jì)算結(jié)果表明，機(jī)械挖除邊坡土體擴(kuò)建路基過程中既有右線圍巖最大位移量及襯砌的內(nèi)力都處于相對安全范圍內(nèi)。但機(jī)械開挖進(jìn)度較慢，尤其遇到堅(jiān)硬巖石時(shí)很可能會(huì)延誤工期，因此有必要在機(jī)械開挖工法中穿插爆破開挖工法以提高施工進(jìn)度。該研究中開挖邊坡中心距既有隧道49 m，考慮最不利開挖工況，即先用機(jī)械開挖工法依次挖出上2層土體，最后1層進(jìn)行爆破開挖，研究最不利情況下擴(kuò)建路基時(shí)對既有隧道的影響。

圖5 爆破開挖邊坡巖土體各質(zhì)點(diǎn)位移云圖(單位：mm)Fig.5 Nephograms of displacement of each particle of rock and soil during blasting excavation(unit：mm)

與機(jī)械開挖不同，爆破開挖對周圍巖土體的擾動(dòng)更為明顯。從圖5(a)可以明顯看出，在爆破開挖完成后，邊坡兩側(cè)巖土體都向臨空面產(chǎn)生了位移，且最大位移值達(dá)到了2.2 cm，幾乎是機(jī)械開挖作用下巖土體位移值的20倍。與此同時(shí)，也可以發(fā)現(xiàn)，爆炸振動(dòng)產(chǎn)生的橫向沖擊波對邊坡巖土體影響程度更大，如圖5(b)～(c)所示，橫波作用下，邊坡巖土體向臨空面產(chǎn)生的位移最大值達(dá)到了1.8 cm，是縱向沖擊波的5倍，這與已有研究[4]爆破振動(dòng)時(shí)襯砌容易產(chǎn)生水平方向的拉伸破壞的觀點(diǎn)是一致的，縱向沖擊波的邊界效應(yīng)更加明顯，對邊界巖土體產(chǎn)生的影響較大，而對臨空面巖土體的影響有限，因此研究對既有臨近隧道的影響時(shí)應(yīng)重點(diǎn)關(guān)注爆炸產(chǎn)生的橫向沖擊波。

爆破振動(dòng)對支護(hù)結(jié)構(gòu)的影響可從襯砌的動(dòng)應(yīng)力及襯砌結(jié)構(gòu)關(guān)鍵點(diǎn)的震動(dòng)速度來分析[18]，但受施工條件和測量儀器影響，動(dòng)應(yīng)力很難進(jìn)行現(xiàn)場測量，為了便于進(jìn)行現(xiàn)場監(jiān)測，選取襯砌結(jié)構(gòu)關(guān)鍵點(diǎn)的震動(dòng)速度進(jìn)行分析。根據(jù)已有研究[3]，使用爆破開挖工法對既有隧道襯砌迎爆側(cè)的危害程度大于背爆側(cè)，迎爆側(cè)拱腰質(zhì)點(diǎn)的拉應(yīng)力明顯大于其他位置，因此在計(jì)算過程在既有隧道模型迎爆側(cè)襯砌布置4個(gè)測點(diǎn)，如圖6所示。

圖6 既有右線隧道襯砌結(jié)構(gòu)監(jiān)測點(diǎn)布置Fig.6 Layout of monitoring points of lining structure in the existing tunnel

通過計(jì)算可以得到各監(jiān)測點(diǎn)在爆破過程中的振動(dòng)情況，如圖7所示。可以發(fā)現(xiàn)，在襯砌完整的情況下，迎爆側(cè)襯砌結(jié)構(gòu)所受的沖擊力峰值相對于爆破荷載中心峰值略有滯后，但時(shí)間間隔極短，爆炸中心產(chǎn)生的沖擊波在土體材料中幾乎沒有被削弱，這主要是因?yàn)楸ㄖ行呐c隧道的距離很近，短距離內(nèi)土體材料的阻尼效應(yīng)并不明顯，爆炸產(chǎn)生的能量很難被消耗。在迎爆側(cè)各測點(diǎn)中，測點(diǎn)2(側(cè)墻中質(zhì)點(diǎn))處所受影響最大，既有隧道襯砌振動(dòng)速度達(dá)到最大值11.42 cm/s，小于爆破安全規(guī)程的允許振動(dòng)速度15 cm/s，測點(diǎn)3與測點(diǎn)4的弧度較為平緩，爆炸沖擊波在該位置極易發(fā)生反射，繞射反應(yīng)較為顯著，所以質(zhì)點(diǎn)振動(dòng)速度較測點(diǎn)2低。

圖7 各測點(diǎn)振動(dòng)速度歷程曲線Fig.7 Vibration speed history curve of each measuring point

5 結(jié)論

通過現(xiàn)場取樣，結(jié)合室內(nèi)試驗(yàn)測定圍巖及隧道襯砌的力學(xué)參數(shù)，建立三維模型，分析了機(jī)械開挖和爆破開挖過程中巖體變形、既有隧道襯砌變形、既有隧道襯砌內(nèi)力的變化過程，得到以下結(jié)論：

(1)機(jī)械開挖路基的分析結(jié)果表明，采用三步開挖法開挖路基時(shí)，隨著土體的不斷開挖，卸載效應(yīng)不斷增加，對既有隧道的影響也不斷增大，受開挖影響最大的位置在左側(cè)拱肩位置，施工過程中應(yīng)做好既有隧道襯砌的監(jiān)測工作。

(2)在土體開挖范圍內(nèi)，爆破施工振動(dòng)效應(yīng)對隧道襯砌結(jié)構(gòu)的影響大于機(jī)械開挖，在襯砌完整的情況下，襯砌最大應(yīng)力在安全范圍內(nèi)，但當(dāng)襯砌存在裂縫時(shí)，在應(yīng)力疊加效應(yīng)作用下，襯砌內(nèi)最大應(yīng)力會(huì)超出承載力限值，此時(shí)應(yīng)當(dāng)注意提前加固襯砌結(jié)構(gòu)。

(3)受隧道爆破點(diǎn)位置的影響，爆炸沖擊波產(chǎn)生的襯砌應(yīng)力主要集中在側(cè)墻位置，為保證爆破開挖作業(yè)中隧道襯砌結(jié)構(gòu)的安全性，現(xiàn)場監(jiān)測中，應(yīng)以隧道迎爆側(cè)襯砌側(cè)墻處允許的極限抗拉強(qiáng)度對應(yīng)的振速作為控制指標(biāo)。