富水隧道襯砌結(jié)構(gòu)受力特性與合理參數(shù)分析

鄧來， 李化云*， 雷中成， 曹蘇亞

(1.西華大學(xué)土木建筑與環(huán)境學(xué)院， 成都 610039； 2.中鐵隆工程集團(tuán)有限公司， 成都 610041)

近年來，在中國(guó)綜合交通網(wǎng)絡(luò)建設(shè)高速發(fā)展的背景下，中國(guó)西南地區(qū)交通網(wǎng)絡(luò)建設(shè)也在加速發(fā)展。西南地區(qū)具有復(fù)雜的地形地貌，以及多變的地質(zhì)條件，在這種環(huán)境下修建山嶺隧道，富水高水壓?jiǎn)栴}成為一大難題。近些年，中國(guó)大力倡導(dǎo)綠色建筑，由于大量排水會(huì)影響到地區(qū)的水利條件，所以隧道防排水設(shè)計(jì)理念開始從“以排為主”向“以堵為主”轉(zhuǎn)化[1]。這就向隧道襯砌結(jié)構(gòu)的抗水壓能力提出了更高的要求。因此，研究抗水壓襯砌結(jié)構(gòu)顯得尤為關(guān)鍵。

一些學(xué)者針對(duì)此問題進(jìn)行了相關(guān)研究，丁燕平等[2]采用荷載結(jié)構(gòu)模型對(duì)抗水壓襯砌結(jié)構(gòu)進(jìn)行分析，基于三種隧道襯砌結(jié)構(gòu)的受力特點(diǎn)，設(shè)計(jì)了高水壓下富水隧道的襯砌結(jié)構(gòu)形式；陳耀華[3]根據(jù)石板山富水隧道工程，研究了隧道襯砌背后水壓力的變化規(guī)律，并對(duì)原設(shè)計(jì)的二襯厚度進(jìn)行了優(yōu)化；梁巍[4]在同時(shí)考慮水壓和圍壓的情況下，使用ANSYS有限元軟件分析和計(jì)算不同的襯砌結(jié)構(gòu)斷面，從而獲得更為經(jīng)濟(jì)合理的襯砌結(jié)構(gòu)形式；陳五二[5]依據(jù)實(shí)際工程的襯砌結(jié)構(gòu)參數(shù)，采用荷載-結(jié)構(gòu)模型，模擬三種襯砌斷面的受力情況，得到了抗水壓襯砌結(jié)構(gòu)的內(nèi)力分布特征，選出了合理的襯砌結(jié)構(gòu)形式，以滿足不同的圍巖級(jí)別要求；丁浩等[6]在龍?zhí)端淼篮桶自扑淼拦こ痰幕A(chǔ)上，通過數(shù)值仿真，研究了富水公路隧道襯砌結(jié)構(gòu)的改進(jìn)方案。一些學(xué)者采用“限制排放”的隧道防排水設(shè)計(jì)理念，對(duì)富水隧道的開裂機(jī)理進(jìn)行模型試驗(yàn)研究[7-10]，得出了一些更為直觀的結(jié)論，為相關(guān)隧道結(jié)構(gòu)設(shè)計(jì)提供了依據(jù)。

綜上所述，針對(duì)抗水壓隧道襯砌結(jié)構(gòu)開展了諸多研究，也取得了一定的成果。但既有研究大多是采用模型試驗(yàn)和數(shù)值模擬的方法確定水荷載，其結(jié)果與工程實(shí)際有一定出入。相較于模型試驗(yàn)和數(shù)值模擬，現(xiàn)場(chǎng)試驗(yàn)具有試驗(yàn)結(jié)果更直觀，試驗(yàn)數(shù)據(jù)更準(zhǔn)確的特點(diǎn)。基于此，現(xiàn)以峨漢高速豹貍崗隧道工程為依托，在試驗(yàn)斷面安裝水壓監(jiān)測(cè)裝置，獲得水壓力分布規(guī)律，為水荷載的確定提供依據(jù)。基于現(xiàn)場(chǎng)試驗(yàn)結(jié)果，采用ANSYS有限元軟件，建立荷載-結(jié)構(gòu)模型，研究不同水壓下襯砌結(jié)構(gòu)的受力特性及合理的支護(hù)結(jié)構(gòu)參數(shù)，為類似隧道的結(jié)構(gòu)設(shè)計(jì)提供參考。

1 現(xiàn)場(chǎng)監(jiān)測(cè)方案及襯砌水壓力確定

1.1 工程概況

豹貍崗隧道進(jìn)口位于峨眉市金口河區(qū)和平鄉(xiāng)解放村園木溝左岸中下部，出口位于枕頭壩水電站江溝料場(chǎng)尾端的岔河右岸陡坡中下部。隧道全長(zhǎng)3 650 m，屬于特長(zhǎng)高速公路隧道。豹貍崗隧道監(jiān)測(cè)斷最大開挖跨度14 m，高度11 m，最大埋深696 m。隧道穿越的地層主要為炭質(zhì)板巖和砂質(zhì)板巖，圍巖級(jí)別為Ⅴ級(jí)。巖層及支護(hù)結(jié)構(gòu)的主要物理力學(xué)指標(biāo)如表1所示。部分施工段巖體受復(fù)雜地質(zhì)構(gòu)造影響嚴(yán)重，隧址區(qū)域地下水豐富，且常年有穩(wěn)定的補(bǔ)給水源，因此隧道存在承受高水壓的問題?；谝陨瞎こ趟牡刭|(zhì)特點(diǎn)，依據(jù)相關(guān)規(guī)范要求，選擇C25噴射砼作為初期支護(hù)，厚240 mm，二襯選用C35鋼筋混凝土結(jié)構(gòu)，厚400 mm。隧道開挖輪廓如圖1所示。

表1 物理力學(xué)參數(shù)Table 1 Physical and mechanical parameter

圖1 隧道開挖輪廓示意圖Fig.1 Schematic diagram of tunnel excavation outline

1.2 斷面位置及測(cè)點(diǎn)布置

選取隧道左洞ZK73+586、ZK73+486和ZK73+386三個(gè)典型斷面作為試驗(yàn)斷面。由工程地質(zhì)勘查報(bào)告可知，三個(gè)試驗(yàn)斷面的水頭高度依次為29、25、20 m。試驗(yàn)儀器選用振弦式滲壓計(jì)，將其安裝在初期支護(hù)與防水板之間，用以測(cè)量二襯背后的水壓力。為了避免圍巖收斂時(shí)將滲壓計(jì)壓壞，滲壓計(jì)埋設(shè)時(shí)與初期支護(hù)間預(yù)留一定空隙。測(cè)點(diǎn)分別布置于拱頂、拱腰、邊墻、墻腳及仰拱，共10個(gè)點(diǎn)，具體布置方案如圖2所示，現(xiàn)場(chǎng)滲壓計(jì)的布設(shè)如圖3所示。

圖2 測(cè)點(diǎn)布置圖Fig.2 Layout of measuring points

圖3 現(xiàn)場(chǎng)滲壓計(jì)布設(shè)Fig.3 Layout of site water pressure gauge

1.3 試驗(yàn)斷面監(jiān)測(cè)數(shù)據(jù)分析

本次現(xiàn)場(chǎng)監(jiān)測(cè)共耗時(shí)75 d，通過對(duì)數(shù)據(jù)的采集和整理，繪制出水壓時(shí)程變化曲線如圖4所示。在滲壓計(jì)安裝1～20 d，三個(gè)斷面的仰拱部位水壓值變化突出，并在后期不斷增長(zhǎng)，在50～60 d趨于穩(wěn)定，其穩(wěn)定值分別為0.28、0.25、0.20 MPa；二襯結(jié)構(gòu)的水壓值均增長(zhǎng)緩慢，一般在40～50 d趨于穩(wěn)定，其水壓力均明顯小于仰拱部位?？赡苁嵌r背后的排水盲管起到了較好的卸壓作用。ZK73+586斷面的邊墻部位水壓值一直在0 MPa上下波動(dòng)，可能是滲壓計(jì)已經(jīng)損壞。綜上所述，試驗(yàn)斷面滲水明顯，襯砌承受較大的水壓力。三個(gè)試驗(yàn)斷面的襯砌水壓力呈典型的“燈泡型”分布，其水壓力分布特點(diǎn)如圖5所示。

圖4 監(jiān)測(cè)斷面水壓時(shí)程曲線Fig.4 Time history curve of water pressure monitoring section

圖5 監(jiān)測(cè)斷面水壓力分布Fig.5 Water pressure distribution of monitoring section

1.4 外水壓力計(jì)算方法的確定

襯砌水壓力的確定是一個(gè)復(fù)雜的問題，目前工程界對(duì)水壓的具體計(jì)算方法還沒有定論，大多是以既有工程經(jīng)驗(yàn)為依據(jù)對(duì)抗水壓襯砌結(jié)構(gòu)進(jìn)行設(shè)計(jì)，其設(shè)計(jì)結(jié)果普遍偏于保守。過于保守的二次襯砌厚度雖然有利于提高二襯剛度，但會(huì)造成材料的浪費(fèi)。因此，采用合適的外水壓計(jì)算方法對(duì)于抗水壓襯砌結(jié)構(gòu)的設(shè)計(jì)具有重要意義。

由于圍巖裂隙對(duì)地下水滲透過程的消散作用，以及排水系統(tǒng)的卸壓作用[11]，導(dǎo)致襯砌背后的水壓力并不完全等于靜水壓力，在設(shè)計(jì)過程中應(yīng)對(duì)其進(jìn)行折減[12]。對(duì)三個(gè)試驗(yàn)斷面的現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)進(jìn)行線性擬合，得出了二次襯砌結(jié)構(gòu)及仰拱的水壓折減系數(shù)，擬合結(jié)果如圖6所示。

圖6 線性擬合曲線Fig.6 Linear fitting curve

基于擬合結(jié)果，得出水壓力計(jì)算公式為：ps=ρghα；其中：ps為襯砌水壓力；ρ為水密度；h為水頭高度；g為重力加速度；α為水壓折減系數(shù)，二襯結(jié)構(gòu)取0.31，仰拱取0.86。此結(jié)論可為豹貍崗隧道二次襯砌結(jié)構(gòu)的設(shè)計(jì)提供參考。

2 模型的建立及計(jì)算工況設(shè)置

采用全斷面承受外水壓設(shè)計(jì)形式，利用ANSYS有限元軟件分別對(duì)水頭高度為10、20、30、40、50 m時(shí)的5種工況進(jìn)行計(jì)算。

利用ANSYS軟件建立荷載-結(jié)構(gòu)二維模型的基本過程為：首先選取模型輪廓線，再對(duì)輪廓線進(jìn)行網(wǎng)格劃分，隨后引入邊界條件對(duì)模型的變形及位移進(jìn)行約束，最后對(duì)模型施加荷載進(jìn)行計(jì)算。本次建模選取水壓監(jiān)測(cè)斷面的襯砌中線斷面為隧道輪廓線；根據(jù)計(jì)算精度要求，將襯砌模型劃分為72段，利用Beam3單元進(jìn)行模擬；由于圍巖對(duì)襯砌的變形具有約束作用，因此引入邊界條件——地層彈簧對(duì)模型的變形進(jìn)行約束。根據(jù)規(guī)范要求及現(xiàn)場(chǎng)圍巖情況，取彈簧的彈性系數(shù)為150 MPa/m，此過程采用Link10單元進(jìn)行模擬；最后將圍巖壓力和水壓力以等效節(jié)點(diǎn)荷載的形式施加在襯砌模型上。計(jì)算采用彈性本構(gòu)模型，假定襯砌為各向同性且只承受壓力。根據(jù)《公路隧道設(shè)計(jì)規(guī)范》[13]中的相關(guān)公式計(jì)算圍巖荷載；基于襯砌水壓力分布特點(diǎn)，將水荷載等效為梯形荷載，并對(duì)其進(jìn)行折減?；跀?shù)值計(jì)算結(jié)果，確定出不同水頭高度作用下結(jié)構(gòu)的安全系數(shù)，并以此判斷襯砌結(jié)構(gòu)設(shè)計(jì)參數(shù)的合理性。對(duì)應(yīng)的工況及相關(guān)計(jì)算參數(shù)如表2所示。

表2 工況及計(jì)算參數(shù)Table 2 Working conditions and calculation parameters

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

利用ANSYS有限元軟件對(duì)二襯結(jié)構(gòu)進(jìn)行數(shù)值模擬，同時(shí)考慮圍巖荷載和水荷載，獲得不同工況下襯砌結(jié)構(gòu)的內(nèi)力分布特征。

3.1 軸力分析

由于在各種工況下，襯砌結(jié)構(gòu)的內(nèi)力分布特征相似，因此以工況三(水頭高度為30 m)為例，對(duì)結(jié)構(gòu)受力特征進(jìn)行分析。當(dāng)水頭高度為30 m時(shí)，數(shù)值計(jì)算得到的二襯軸力分布特點(diǎn)如圖7所示。分析圖7可知，在結(jié)構(gòu)全斷面承壓的情況下，隧道襯砌結(jié)構(gòu)的軸力呈對(duì)稱分布，這與丁燕平等[14]的研究結(jié)論相一致。云圖顯示，襯砌結(jié)構(gòu)在水壓作用下表現(xiàn)為全斷面受壓。這一狀態(tài)有利于混凝土結(jié)構(gòu)受力。

圖7 水頭高度30 m時(shí)的軸力云圖Fig.7 Axial force nephogram at water head height of 30 m

圖8 結(jié)構(gòu)各部位軸力變化曲線Fig.8 Change curve of axial force at each part of the structure

結(jié)構(gòu)各部位軸力絕對(duì)值與水頭關(guān)系曲線如圖8所示。在各種工況下，墻腳所承受的軸力最大；且隨著水頭高度的增加，墻腳，仰拱和邊墻的軸力增長(zhǎng)最多，增長(zhǎng)率分別為76%、75%和75%；拱腰的軸力增加最少，增長(zhǎng)率為64%；當(dāng)水頭高度在30 m以下時(shí)，隨著水頭高度的增加，結(jié)構(gòu)各部位軸力增長(zhǎng)緩慢。當(dāng)水頭高度超過30 m時(shí)，結(jié)構(gòu)各部位軸力增長(zhǎng)速度明顯加快，其中墻腳和仰拱的軸力增加最快。因此在結(jié)構(gòu)設(shè)計(jì)時(shí)，應(yīng)注意墻腳和仰拱的承載能力。

3.2 彎矩分析

彎矩分析同樣以工況三為例，二次襯砌結(jié)構(gòu)的彎矩圖如圖9所示。襯砌結(jié)構(gòu)的彎矩同樣呈對(duì)稱分布，拱腰和墻腳部位靠近圍巖一側(cè)承受拉應(yīng)力。拱頂、拱腰和仰拱部位靠近隧道凈空一側(cè)受拉。在不同工況下，墻腳部位靠近圍巖一側(cè)的拉應(yīng)力最大，這一受力狀態(tài)有利于發(fā)揮出圍巖的緊箍效應(yīng)。邊墻的彎矩最小，在同一工況下，僅為墻腳的1/4。

圖9 水頭高度30 m時(shí)的彎矩云圖Fig.9 Bending moment nephogram when the head height is 30 m

結(jié)構(gòu)各部位彎矩絕對(duì)值與水頭關(guān)系曲線如圖10所示。各部位彎矩值的變化趨勢(shì)與軸力變化趨勢(shì)類似：當(dāng)水頭高度小于30 m時(shí)，各部位彎矩值增長(zhǎng)緩慢；當(dāng)水頭高度超過30 m時(shí)，彎矩值的增長(zhǎng)速度加快。其中拱頂和墻腳的彎矩增長(zhǎng)速度最快。

3.3 安全系數(shù)分析

在結(jié)構(gòu)設(shè)計(jì)過程中，安全系數(shù)是判斷結(jié)構(gòu)設(shè)計(jì)可靠度的一個(gè)重要參數(shù)，《公路隧道設(shè)計(jì)規(guī)范》規(guī)定，當(dāng)荷載組合為永久荷載+基本可變荷載時(shí)，鋼筋混凝土結(jié)構(gòu)的安全系數(shù)不應(yīng)小于2.0[13]。

本設(shè)計(jì)結(jié)構(gòu)的安全系數(shù)云圖如圖11所示。在全截面承受水壓的情況下，結(jié)構(gòu)安全系數(shù)呈對(duì)稱分布，且在不同工況下，最小值均位于墻腳處。安全系數(shù)變化曲線如圖12所示。當(dāng)水壓增大時(shí)，結(jié)構(gòu)各部位安全系數(shù)均出現(xiàn)明顯下降；當(dāng)水頭高度超過50 m時(shí)，墻腳部位的安全系數(shù)已不符合相關(guān)要求，可能會(huì)出現(xiàn)開裂現(xiàn)象[15]，此時(shí)應(yīng)該增加二襯厚度或提高二襯配筋率。

圖10 結(jié)構(gòu)各部位彎矩變化曲線Fig.10 Bending moment change curve of each part of the structure

圖11 水頭高度為30 m時(shí)的安全系數(shù)圖Fig.11 Safety factor diagram when the head height is 30 m

圖12 結(jié)構(gòu)各部位安全系數(shù)變化曲線Fig.12 Change curve of safety factor of each part of the structure

4 結(jié)論與討論

(1)三個(gè)試驗(yàn)斷面的襯砌水壓力呈典型的“燈泡”型分布，且水壓力最大位置均出現(xiàn)在仰拱處；二次襯砌結(jié)構(gòu)的水壓力折減系數(shù)為0.31，仰拱的水壓力折減系數(shù)為0.86。

(2)當(dāng)結(jié)構(gòu)全斷面承受水荷載時(shí)，軸力、彎矩及安全系數(shù)均呈對(duì)稱分布；墻腳為結(jié)構(gòu)的薄弱部位，在設(shè)計(jì)時(shí)應(yīng)加強(qiáng)墻腳部位的承載能力。

(3)在水頭高度從10 m增加到50 m的過程中，墻腳、仰拱和邊墻的軸力值增量最大，增長(zhǎng)率分別為76%、75%和75%。當(dāng)水壓增大時(shí)，應(yīng)著重控制這三個(gè)部位的安全性。

(4)當(dāng)水頭高度小于30 m時(shí)，結(jié)構(gòu)各部位軸力值和彎矩值隨水壓的增大，變化并不明顯；當(dāng)水頭高度大于30 m時(shí)，結(jié)構(gòu)各部位軸力值和彎矩值的增長(zhǎng)速度明顯加快，軸力增加最快的部位為墻腳，彎矩增加最快的部位為拱頂。

(5)隨著水頭高度的增加，結(jié)構(gòu)各部位安全系數(shù)均出現(xiàn)明顯減小的現(xiàn)象；當(dāng)水頭高度為50 m時(shí)，墻腳部位的安全系數(shù)小于2.0，說明原設(shè)計(jì)結(jié)構(gòu)不能承受超過50 m水頭高度的水壓。