錢文喜，耿大新，梁國卿

(1. 江西省公路工程有限責(zé)任公司，江西 南昌 330006；2. 華東交通大學(xué) 土木建筑學(xué)院，江西 南昌 330013；3. 江西省天馳高速科技發(fā)展有限公司，江西 南昌 330103)

0 引言

公路隧道圍巖極其復(fù)雜的非線性、各向異性和隨時(shí)間變化的力學(xué)屬性，決定了圍巖與襯砌結(jié)構(gòu)相互作用的復(fù)雜性。在軟弱圍巖隧道運(yùn)營期間，二次襯砌與圍巖可以同時(shí)作為承載主體，而圍巖蠕變對隧道穩(wěn)定性及襯砌結(jié)構(gòu)安全性具有重要的影響。因此，許多學(xué)者開展了相關(guān)研究：劉甲榮等[1]考慮圍巖蠕變效應(yīng)，分析了不同埋深及偏壓條件下襯砌結(jié)構(gòu)表面拉、壓應(yīng)變與損傷值之間的關(guān)系；李建軍等[2]基于Burgers模型分析了圍巖蠕變作用對襯砌結(jié)構(gòu)內(nèi)力的影響；王迎超等[3]基于廣義Kelvin模型分析了不同時(shí)機(jī)支護(hù)下襯砌結(jié)構(gòu)的受力規(guī)律；徐國文等[4]建立錨桿-圍巖復(fù)合體流變模型分析了圍巖蠕變作用對裂損襯砌長期安全性的影響；Guan等[5-6]通過對Ureshino隧道進(jìn)行跟蹤監(jiān)測，分析了圍巖蠕變與隧道持續(xù)變形的關(guān)系；王中文等[7]考慮隧道圍巖的蠕變特性，利用初襯變形理論公式分析了隧道初襯最終變形及二次襯砌合理的支護(hù)時(shí)機(jī)；唐葭等[8]基于Burgers模型推導(dǎo)出了考慮圍巖蠕變效應(yīng)的圍巖抗力系數(shù)公式，進(jìn)而分析了紅砂巖隧道圍巖的抗力系數(shù)；Fan等[9]采用室內(nèi)模型試驗(yàn)對泥質(zhì)軟巖隧道圍巖的蠕變力學(xué)行為進(jìn)行了研究。此外，左昌群等[10]考慮軟弱圍巖蠕變特性，分析了不同支護(hù)形式下圍巖與初期支護(hù)結(jié)構(gòu)的相互作用；文獻(xiàn)[11-14]對軟巖隧道圍巖的蠕變特性、蠕變模型及蠕變參數(shù)進(jìn)行了大量研究。

目前，常用的模擬隧道圍巖蠕變的元件組合模型有Kelvin模型、西原模型、Burgers模型、Bingham模型等[15]。但不同地區(qū)、不同類別的巖石力學(xué)性狀差異較大，需要選擇適用的元件組合模型來描述巖石的蠕變特性。本研究以江西九景高速公路隧道為依托，基于該隧道Ⅳ級(jí)圍巖區(qū)段泥質(zhì)粉砂巖蠕變試驗(yàn)結(jié)果，利用Cvisc黏彈塑性模型對其蠕變模型及參數(shù)進(jìn)行辨識(shí)，進(jìn)而分析圍巖蠕變效應(yīng)對運(yùn)營隧道襯砌結(jié)構(gòu)安全性的影響，以期為公路隧道運(yùn)營期的管理和養(yǎng)護(hù)提供參考。

1 工程概況

九景高速公路隧道為雙向兩車道，隧道凈寬10.25 m，凈高7.425 m，洞身襯砌采用三心圓曲墻式，運(yùn)營時(shí)間已超過10 a。洞身穿越地層巖性主要為泥盆系細(xì)砂巖和志留系粉砂巖，巖體完整性以較完整～較破碎為主。其中，Ⅳ級(jí)圍巖隧道斷面如圖1所示。

圖1 Ⅳ級(jí)圍巖復(fù)合式襯砌Fig.1 Composite lining of Grade IV surrounding rock

2 砂巖蠕變模型辨識(shí)

2.1 蠕變試驗(yàn)

通過單軸抗壓強(qiáng)度試驗(yàn)，獲得了5組φ50×100 mm 圓柱體巖樣破壞時(shí)最大軸向壓力Pi和抗壓強(qiáng)度Rci(見表1)，其中，i=1,2,…,5。巖樣取自九景高速公路隧道，為泥質(zhì)粉砂巖，無明顯層理，且節(jié)理裂隙不發(fā)育。

表1 泥質(zhì)粉砂巖單軸抗壓強(qiáng)度Tab.1 Uniaxial compressive strength of argillaceous siltstone

通過YSR-300型巖石三軸蠕變試驗(yàn)系統(tǒng)對泥質(zhì)粉砂巖進(jìn)行單軸壓縮蠕變試驗(yàn)，獲得了巖樣在0.7Rc，0.8Rc，0.9Rc應(yīng)力水平下較能反映其蠕變特性的試驗(yàn)數(shù)據(jù)曲線(見圖2)。

圖2 不同應(yīng)力水平下的非線性擬合結(jié)果Fig.2 Nonlinear fitting results at different stress levels

2.2 蠕變模型與參數(shù)確定

Cvisc模型由Kelvin體、Maxwell體和塑性體共同構(gòu)成(見圖3)。當(dāng)應(yīng)力水平較低時(shí)，變形速率逐漸減小進(jìn)入衰減蠕變階段，之后逐漸趨于某一恒定值進(jìn)入穩(wěn)定蠕變階段；當(dāng)應(yīng)力水平等于或超過某一臨界應(yīng)力值σs后，逐漸轉(zhuǎn)化為加速蠕變過程。

在一維應(yīng)力狀態(tài)下，Cvisc模型蠕變本構(gòu)方程為：

(1)當(dāng)σ<σs時(shí)，

(1)

(2)當(dāng)σ≥σs時(shí)，

(2)

式中，EK、EM分別為Kelvin體、Maxwell體的彈性模量；ηK、ηM分別為Kelvin體、Maxwell體的黏性系數(shù)；σ0為常應(yīng)力；εP為塑性體應(yīng)變。

圖3 Cvisc模型Fig.3 Cvisc model

蠕變方程(1)，(2)屬于非線性函數(shù)形式，為求其中的蠕變參數(shù)，可采用最小二乘法中Boltzmann法處理試驗(yàn)數(shù)據(jù)，即采用線性疊加的間接方法求解[16]。結(jié)果表明，圖2中的擬合曲線與試驗(yàn)數(shù)據(jù)曲線的相關(guān)系數(shù)R2在0.92～0.96之間，因此，可以認(rèn)為Cvisc模型能夠描述泥質(zhì)粉砂巖的衰減蠕變和穩(wěn)定蠕變關(guān)系。根據(jù)室內(nèi)蠕變試驗(yàn)結(jié)果，得到Cvisc模型在0.7Rc，0.8Rc，0.9Rc應(yīng)力水平下的蠕變參數(shù)(見表2)。

3 數(shù)值模擬與結(jié)果分析

3.1 模型建立

以Ⅳ級(jí)圍巖公路隧道襯砌結(jié)構(gòu)為研究對象，利用FLAC3D建立兩車道公路隧道三維數(shù)值模型，如圖4所示。隧道模型左、右側(cè)均距離隧道中心線50 m(約4倍隧道開挖跨度)，洞周拱頂以上取24 m，仰拱以下取36 m，縱向取60 m。圍巖采用實(shí)體單元模擬，本構(gòu)模型采用Mohr-Coulomb模型，初始地應(yīng)力場僅考慮重力作用。此時(shí)，邊界條件采用位移約束，即左右兩側(cè)采用水平位移約束，頂部自由無約束，底部采用豎向位移約束。

表2 Cvisc模型參數(shù)Tab.2 Parameters of Cvisc model

圖4 三維隧道數(shù)值模型Fig.4 3D tunnel numerical model

在隧道施工期間，假定圍巖在開挖條件下的變形完全為彈性或同時(shí)含有塑性變形。因此，Cvisc模型中用以描述圍巖時(shí)效變形的黏性元件應(yīng)不發(fā)生作用，圍巖變形僅由模型中的彈性或塑性元件來體現(xiàn)，故此時(shí)蠕變方程(1)、(2)中的t=0。于是，通過命令(SET CREEP OFF)來屏蔽蠕變模型的時(shí)效作用，即僅對圍巖進(jìn)行靜力計(jì)算分析。其中，Ⅳ級(jí)圍巖的物理力學(xué)參數(shù)見表3。

表3 圍巖物理力學(xué)參數(shù)Tab.3 Physical and mechanical parameters of surrounding rock

通過提高洞周圍巖變形模量、黏聚力等指標(biāo)的方法來模擬錨桿加固區(qū)的力學(xué)行為，加固區(qū)厚度取4 m，Ⅳ級(jí)圍巖錨桿加固區(qū)域的物理力學(xué)參數(shù)見表4。

假定圍巖與初期支護(hù)分擔(dān)60%釋放荷載，二次襯砌分擔(dān)40%釋放荷載。于是，在施加襯砌前，對隧道周邊圍巖進(jìn)行某一比例卸荷，即選定應(yīng)力釋放

表4 錨桿加固區(qū)圍巖物理力學(xué)參數(shù)Tab.4 Physical and mechanical parameters of surrounding rock in anchor reinforcement area

率λ=0.6。采用實(shí)體單元模擬噴混凝土層，本構(gòu)模型采用線彈性模型；二次襯砌采用襯砌結(jié)構(gòu)單元(Liner)模擬，初期支護(hù)的物理力學(xué)參數(shù)見表5。

表5 初期支護(hù)物理力學(xué)參數(shù)Tab.5 Physical and mechanical parameters of initial support

考慮到噴混凝土層與二次襯砌之間防水層的影響，單元界面接觸參數(shù)選取見表6。

表6 單元界面接觸參數(shù)Tab.6 Contact parameters of unit interface

在隧道運(yùn)營期間，假定圍巖總變形由開挖導(dǎo)致的彈塑性變形和蠕變變形共同構(gòu)成。因此，Cvisc模型中用以描述圍巖時(shí)效變形的黏性元件開始發(fā)生作用。于是，基于靜力計(jì)算分析結(jié)果，通過命令(SET CREEP ON)來開啟蠕變模型的時(shí)效作用，即對圍巖進(jìn)行蠕變計(jì)算分析。其中，Cvisc模型所需的Mohr-Coulomb強(qiáng)度參數(shù)見表7。對泥質(zhì)粉砂巖在0.7Rc,0.8Rc，0.9Rc應(yīng)力水平下的蠕變參數(shù)進(jìn)行平均值計(jì)算，得到圍巖蠕變參數(shù)。

表7 圍巖蠕變參數(shù)Tab.7 Creep parameters of surrounding rock

Cvisc模型的最大蠕變時(shí)間步長為：

(3)

為了使計(jì)算模型更好地收斂，設(shè)置蠕變計(jì)算時(shí)間步長為2.5×10-3，最大蠕變時(shí)間步長為5.0×10-3，最小蠕變時(shí)間步長為5.0×10-4。將蠕變計(jì)算總歷時(shí)定為3 650 d(即蠕變計(jì)算的總歷時(shí)為10 a)[17]。

3.2 結(jié)果分析

我國《公路隧道設(shè)計(jì)規(guī)范》[18](JTG3370.1—2018)規(guī)定，公路隧道支護(hù)結(jié)構(gòu)需滿足驗(yàn)算的安全系數(shù)要求。規(guī)范規(guī)定混凝土偏心受壓構(gòu)件的抗壓和抗拉強(qiáng)度驗(yàn)算需滿足式(4)和式(5)。

KcN≤φαRαbh,

(4)

(5)

在隧道運(yùn)營10 a過程中，單考慮圍巖蠕變作用，將數(shù)值模擬的計(jì)算結(jié)果代入式(4)和式(5)得到二次襯砌拱頂、拱肩和拱腰位置安全系數(shù)隨時(shí)間的變化曲線(見圖5)。從圖中可以看出，隨著二次襯砌服役時(shí)間的增加，拱頂、拱肩和拱腰位置安全系數(shù)總體呈逐漸減小的趨勢，其中拱頂處安全系數(shù)降低最為顯著，表明圍巖蠕變變形壓力隨著時(shí)間逐漸增加。

圖5 襯砌結(jié)構(gòu)安全系數(shù)隨時(shí)間變化曲線Fig.5 Curves of lining structural safety factor varying with time

隧道運(yùn)營10 a后，二次襯砌軸力、彎矩、安全系數(shù)的分布，見圖6～圖8。從圖中可以看出，襯砌結(jié)構(gòu)最大軸力出現(xiàn)在邊墻位置，為2 700 kN，并沿拱頂和拱底位置環(huán)向減小，最小軸力出現(xiàn)在拱底位置，為757 kN；最大正彎矩出現(xiàn)在拱腳位置附近，為312 kN·m，最大負(fù)彎矩出現(xiàn)在拱底位置，為40.9 kN·m；安全系數(shù)最小值出現(xiàn)在邊墻位置，為2.9，最大安全系數(shù)出現(xiàn)在仰拱底位置，為72.3。這表明，在同一支護(hù)時(shí)間，二次襯砌支護(hù)結(jié)構(gòu)承載較小的位置(即安全系數(shù)較高的位置)，圍巖的蠕變變形量較大，因而支護(hù)結(jié)構(gòu)承受的圍巖壓力相對較小，但過小的支護(hù)承載又會(huì)導(dǎo)致圍巖蠕變變形而增加圍巖壓力，進(jìn)而不利于襯砌結(jié)構(gòu)的長期安全。

圖6 襯砌結(jié)構(gòu)軸力(單位: kN)Fig.6 Axial force of lining structure(unit：kN)

圖7 襯砌結(jié)構(gòu)彎矩(單位: kN·m)Fig.7 Moment of lining structure(unit：kN·m)

圖8 襯砌結(jié)構(gòu)安全系數(shù)Fig.8 Safety factor of lining structure

3.3 隧道埋深的影響分析

在隧道運(yùn)營10 a過程中，考慮圍巖蠕變作用，在11.7 m(hq)，23.4 m(2hq)，35.1 m(3hq)，46.8 m(4hq)的隧道埋深下，將數(shù)值模擬的計(jì)算結(jié)果代入式(4)和式(5)得到二次襯砌拱頂、拱肩、拱腰位置的安全系數(shù)隨其服役時(shí)間的變化曲線(見圖9)。從圖中可以看出，在圍巖蠕變變形壓力作用下，隨著隧道埋深的增加，襯砌結(jié)構(gòu)拱頂、拱肩、拱腰位置在同一服役年限的安全系數(shù)逐漸增大，其中拱頂位置的增長幅度最大。這表明，在同一支護(hù)時(shí)間，二次襯砌支護(hù)結(jié)構(gòu)的承載隨著隧道埋深的增加而減小，即圍巖自身承擔(dān)了較大部分的因蠕變變形而增加的圍巖壓力，從而對二次襯砌的長期安全有利。

圖9 不同埋深下襯砌安全系數(shù)隨時(shí)間變化曲線Fig.9 Curves of lining safety factor varying with time under different depths

4 結(jié)論

以九景高速公路隧道為背景，利用有限差分軟件FLAC3D建立了兩車道公路隧道三維數(shù)值模型，考慮了圍巖蠕變作用，對Ⅳ級(jí)圍巖公路隧道在10 a運(yùn)營期間襯砌結(jié)構(gòu)的安全性進(jìn)行了分析，得出以下結(jié)論：

(1)通過對泥質(zhì)粉砂巖進(jìn)行單軸壓縮蠕變試驗(yàn)，并對試驗(yàn)數(shù)據(jù)曲線進(jìn)行非線性擬合。結(jié)果表明，其相關(guān)系數(shù)R2在0.92～0.96之間，可認(rèn)為Cvisc模型能夠很好地描述泥質(zhì)粉砂巖的衰減蠕變和穩(wěn)定蠕變關(guān)系。

(2)單考慮圍巖蠕變作用，在同一支護(hù)時(shí)間，二次襯砌安全系數(shù)較高的位置，支護(hù)結(jié)構(gòu)承受的圍巖壓力相對較小，圍巖的蠕變變形量較大，但過小的支護(hù)承載又會(huì)導(dǎo)致圍巖蠕變變形而增加圍巖壓力，進(jìn)而不利于運(yùn)營隧道襯砌結(jié)構(gòu)的長期安全。

(3)同時(shí)考慮隧道埋深的影響，在同一支護(hù)時(shí)間，二次襯砌支護(hù)結(jié)構(gòu)的承載隨著隧道埋深的增加而減小，即圍巖能夠承擔(dān)較大部分的因蠕變變形而增加的圍巖壓力，從而對運(yùn)營隧道襯砌結(jié)構(gòu)的長期安全有利。