李苗，徐霞飛，付貴海

偏壓情況下連拱隧道施工順序的數(shù)值分析

李苗1，徐霞飛2，付貴海1

(1.湖南城市學(xué)院土木工程學(xué)院，湖南益陽(yáng)413000；2.湖南省武靖高速公路建設(shè)開發(fā)有限公司，湖南邵陽(yáng)422400)

采用有限元法數(shù)值模擬云南思小高速公路的雙連拱隧道中導(dǎo)洞施工法分步開挖過程，研究雙連拱隧道在偏壓作用下施工過程的力學(xué)行為．通過動(dòng)態(tài)數(shù)值模擬，得到隧洞在不同開挖順序下圍巖的應(yīng)力變化和塑性區(qū)的發(fā)展范圍以及中隔墻、襯砌的應(yīng)力情況；對(duì)比分析不同施工方案各關(guān)鍵工序的模擬結(jié)果，并從應(yīng)力轉(zhuǎn)移的角度分析洞室開挖的影響，得到偏壓作用下連拱隧道施工順序的優(yōu)選方案．

偏壓連拱隧道；施工順序；應(yīng)力轉(zhuǎn)移；數(shù)值分析

上世紀(jì)30年代開始了連拱隧道的建設(shè)，70年代得到較大發(fā)展．與分離式獨(dú)立隧道相比，其優(yōu)勢(shì)在于[1]：節(jié)省引線占地，改善公路線形，地形受限??；結(jié)構(gòu)跨度大、整體性好；綜合造價(jià)較獨(dú)立隧道低．目前，連拱隧道的設(shè)計(jì)與施工均達(dá)到較高水平．德國(guó)Wuerzburg隧道長(zhǎng)達(dá)5 100 m[2]，1974年日本首次修建伊祖雙連拱隧道[3]．于上世紀(jì)90年代，我國(guó)陸續(xù)采用連拱隧道方案，在云南、湖南、福建和廣東等省份修建公路連拱隧道達(dá)幾十座[4]，已建和在建的總數(shù)已位居世界之首[5]．我國(guó)雙連拱隧道的基本特征是淺埋中短長(zhǎng)度，80%以上埋深在30 m內(nèi)，與歐美的情況較接近[6-8]，較日本略大，一般不超過60 m[7]．連拱隧道相對(duì)分離式隧道結(jié)構(gòu)復(fù)雜，開挖斷面增大，開挖、支護(hù)交錯(cuò)進(jìn)行，施工中體系需多次轉(zhuǎn)換，隧道圍巖應(yīng)力變化復(fù)雜．目前連拱隧道可供借鑒的施工經(jīng)驗(yàn)較少，缺乏完整適用的設(shè)計(jì)施工方法，相應(yīng)的研究少，導(dǎo)致設(shè)計(jì)與施工具有較大的盲目性．我國(guó)連拱隧道多采用三導(dǎo)洞施工法，該方法施工工序復(fù)雜．中導(dǎo)洞施工法施工工序相對(duì)較少，但要求地質(zhì)條件較好，對(duì)隧道圍巖影響也較大，制約了其推廣．因此，有必要對(duì)圍巖及支護(hù)結(jié)構(gòu)在中導(dǎo)洞施工法不同開挖順序下的受力情況進(jìn)行分析，為連拱隧道的設(shè)計(jì)施工提供科學(xué)依據(jù)．

依據(jù)云南思小高速公路連拱隧道設(shè)計(jì)資料，進(jìn)行二維有限元建模．?dāng)?shù)值模擬不同開挖方案下雙連拱隧道洞室的施工過程，分析隧道圍巖應(yīng)力情況、塑性區(qū)發(fā)展以及中隔墻與襯砌的應(yīng)力狀況．

1 計(jì)算模型

1.1模型建立

選取典型斷面進(jìn)行數(shù)值分析，模型處于偏壓狀態(tài)，左洞為埋深較深一側(cè)，左洞最大埋深為40 m，如圖1所示．圍巖、中隔墻與襯砌均使用等參平面四邊形單元，錨桿為二維桿單元，整個(gè)連拱隧道有限元模型共劃分3 466個(gè)單元．?dāng)嗝骓斆鏋榈乇?，即自由面；斷面兩?cè)(左右邊界)采用水平約束，底部計(jì)算邊界采用豎向約束．隧道有限元模型以及支護(hù)結(jié)構(gòu)有限元模型如圖2、3所示．

圖1 計(jì)算模型示意圖

圖2 連拱隧道有限元模型

圖3 支護(hù)結(jié)構(gòu)有限元模型

1.2模型假定

①所有材料均質(zhì)連續(xù)，且各向同性；

②巖體不考慮構(gòu)造應(yīng)力影響，初始應(yīng)力場(chǎng)僅由自重形成；

③按平面應(yīng)變問題分析，不考慮空間效應(yīng)；④不考慮變形時(shí)效性[9]．

1.3圍巖與支護(hù)材料

基于隧道設(shè)計(jì)資料與《公路隧道設(shè)計(jì)規(guī)范》對(duì)材料的物理力學(xué)指標(biāo)加以選取，如表1所示．

1.4材料屈服準(zhǔn)則

選取DP(Drucker-Prager)材料屈服準(zhǔn)則進(jìn)行圍巖應(yīng)力分析，該準(zhǔn)則對(duì)Von Mises屈服準(zhǔn)則進(jìn)行了修正，考慮了靜水壓力因素對(duì)巖土屈服的影響．Drucker-Prager準(zhǔn)則表達(dá)式如下．

式中I1為應(yīng)力張量第一不變量，J2為應(yīng)力偏張量第二不變量，表達(dá)式分別為：

1s、2s、3s分別為第一、第二及第三主應(yīng)力；a、k都是與巖土內(nèi)摩擦角j、粘聚力c相關(guān)的參數(shù)，

表1 材料參數(shù)

1.5連拱隧道施工順序的模擬分析方案

本次計(jì)算根據(jù)隧道左右洞不同的開挖順序(其中左洞為埋深較深一側(cè))，模擬了三種不同的中導(dǎo)洞法施工過程，方案分類如表2所示．針對(duì)不同方案，分別進(jìn)行數(shù)值計(jì)算，詳細(xì)分析關(guān)鍵工序下洞室圍巖應(yīng)力場(chǎng)的變化以及支護(hù)結(jié)構(gòu)的受力情況，然后對(duì)各施工方案進(jìn)行討論，得到較合理的開挖順序．

表2 開挖方案及關(guān)鍵工序說明

2 結(jié)果分析

對(duì)比分析不同開挖方案各關(guān)鍵工序的計(jì)算結(jié)果，分析數(shù)據(jù)的選取應(yīng)遵循：(一)反映洞室斷面開挖后圍巖、中隔墻與襯砌等支護(hù)結(jié)構(gòu)的應(yīng)力突變；(二)反映圍巖、中隔墻與襯砌等在各階段的最不利受力情況．因此，選取各洞室斷面開挖后支護(hù)階段的數(shù)據(jù)進(jìn)行對(duì)比分析．

計(jì)算初始地應(yīng)力，如圖4、5所示，圍巖1s為-1.04 MPa～3.47 KPa，3s為-2.64 MPa～-0.03 MPa，地表為受拉狀態(tài)．

圖4 原巖應(yīng)力1s

圖5 原巖應(yīng)力3s

對(duì)比三種開挖方案的圍巖應(yīng)力值，拉應(yīng)力（圖6）最大是方案二右洞室先開挖，方案三兩側(cè)洞室同時(shí)開挖次之，方案一左洞室先開挖的拉應(yīng)力最?。粐鷰r壓應(yīng)力值(圖7)最大的開挖方案是兩洞室同時(shí)開挖，最小的方案是先開挖左洞室．綜合以上分析結(jié)果，方案一左洞室先開挖的方案有利于控制圍巖拉、壓應(yīng)力最大值；也說明從圍巖的穩(wěn)定角度考慮，采取埋深較大一側(cè)洞室先開挖的方案更有利．各開挖方案施工過程中應(yīng)力最值出現(xiàn)的位置見表3．

表3 施工過程中圍巖應(yīng)力最值的位置

圖6 圍巖max1s

圖7 圍巖min3s

圖8 圍巖最大塑性應(yīng)變

表4 洞室頂拱最終沉降(mm)

各方案的圍巖頂拱最終沉降計(jì)算結(jié)果如表4所示，沉降值由大到小排序依次是：方案三(兩側(cè)洞室同時(shí)開挖)、方案二(右洞室先開挖)、方案一(左洞室先開挖)．因此，采取先開挖左洞室的方案對(duì)控制洞室頂拱沉降更合適．

洞室圍巖的塑性應(yīng)變變化如圖8所示，洞室上臺(tái)階開挖后，方案三的塑性應(yīng)變?cè)龇畲?，方案一的圍巖塑性應(yīng)變較為穩(wěn)定，采用方案二的圍巖塑性應(yīng)變?cè)龇^大．對(duì)比各方案的最終塑性應(yīng)變值，先施工左洞室的方案最優(yōu)，先施工右洞室為最不利方案．各開挖方案施工過程中圍巖最大塑性應(yīng)變區(qū)位置如表5所示．

表5 施工過程中圍巖最大塑性應(yīng)變區(qū)的位置

對(duì)比各方案初次襯砌的應(yīng)力變化如圖9、10所示，先施工左洞室，初襯拉應(yīng)力最小．左右洞室同時(shí)開挖，初襯壓應(yīng)力最?。紤]到施作二襯前初襯對(duì)洞室的穩(wěn)定作用，為防止初襯受拉開裂，宜選擇先開挖左洞室的方案．各開挖方案施工過程中初襯應(yīng)力最值的位置如表6所示．

圖9 初次襯砌1maxs

圖10 初次襯砌3mins

表6 施工過程中初襯應(yīng)力最值的位置

對(duì)比分析各方案二次襯砌的應(yīng)力，如圖11、12所示，方案三的二襯應(yīng)力很小，表明剛施作完的二襯未發(fā)揮承載能力，故不參與比較．采用先施工左洞室方案，二襯的應(yīng)力較先施工右洞室方案大，但仍滿足強(qiáng)度要求．各開挖方案施工過程中二襯應(yīng)力最值的位置見表7．

圖11 二次襯砌1maxs

圖12 二次襯砌3mins

對(duì)比各方案中隔墻的應(yīng)力值如圖13、14所示，先開挖左洞室中隔墻拉應(yīng)力最大，左右洞室同時(shí)開挖次之．采用先開挖右洞室的方案應(yīng)力最小，因此更有利于控制墻體的最大應(yīng)力．中隔墻應(yīng)力最值位置如表8所示．

表7 施工過程中二襯應(yīng)力最值的位置

圖13 中隔墻max1s

圖14 中隔墻min3s

表8 施工過程中中隔墻應(yīng)力最值的位置

3 應(yīng)力轉(zhuǎn)移的影響

對(duì)地下工程進(jìn)行開挖后，為抵抗不均勻變形，巖體會(huì)自我調(diào)節(jié)，形成壓力拱這一現(xiàn)象．若將圍巖視作結(jié)構(gòu)，巖體處于壓力拱中，且承擔(dān)自重及圍巖以上的荷載，并阻止上方巖體的坍塌，具有拱的力學(xué)特征[10]．洞室的開挖使得拱體自重及其上部荷載的傳遞會(huì)發(fā)生路線偏離(向洞室兩側(cè))，即形成應(yīng)力轉(zhuǎn)移[11]，產(chǎn)生應(yīng)力集中的區(qū)域，如圖15所示．由圖16、17所示，較大的壓應(yīng)力區(qū)出現(xiàn)在已開挖洞室埋深較大的一側(cè)(左側(cè))，與應(yīng)力轉(zhuǎn)移的結(jié)果相符．

圖15 應(yīng)力轉(zhuǎn)移示意圖

圖16 開挖左洞后的圍巖3s

圖17 開挖右洞后的圍巖3s

比較兩側(cè)洞室在方案一、二開挖順序下的計(jì)算結(jié)果，先開挖淺埋深一側(cè)洞室會(huì)使得洞室左側(cè)形成較大的地應(yīng)力區(qū)域．下一步開挖深埋深洞室時(shí)，較大的地應(yīng)力對(duì)施工是個(gè)不利因素，甚至?xí)a(chǎn)生“巖爆”的現(xiàn)象．反之，先開挖深埋深洞室，后開挖淺埋深洞室，則避免了施工過程中大地應(yīng)力區(qū)域的形成，如圖17所示．所以，考慮應(yīng)力轉(zhuǎn)移帶來的影響，在偏壓條件下，較好的開挖方案是先開挖埋深較深的洞室．

4 結(jié)論

依托云南省思小高速公路雙連拱隧道工程，對(duì)連拱隧道開挖順序進(jìn)行數(shù)值模擬分析，并結(jié)合應(yīng)力轉(zhuǎn)移分析，得出以下結(jié)論．

(1)先開挖埋深較深一側(cè)洞室后開挖埋深較淺一側(cè)洞室，對(duì)于圍巖穩(wěn)定、拱頂沉降、塑性應(yīng)變、初支結(jié)構(gòu)受力具有較好效果，采用先淺埋深洞室后深埋深洞的開挖順序，有利于二襯結(jié)構(gòu)和中隔墻的受力．

(2)根據(jù)圍巖、中隔墻以及襯砌的主應(yīng)力最值計(jì)算結(jié)果，可知結(jié)構(gòu)易發(fā)生破壞的位置．

圍巖：仰拱底部、頂拱與中導(dǎo)洞交匯處與中導(dǎo)洞墻腳處；中隔墻：與初襯交接處、底座兩端；初次襯砌：上半斷面初襯拱頂、上下半斷面初襯交接內(nèi)、外側(cè)，即交接處存在主應(yīng)力差值最大值；二次襯砌：中隔墻與二襯頂拱、仰拱交接處，二襯頂拱、仰拱內(nèi)側(cè)．

(3)綜合考慮數(shù)值模擬分析結(jié)果，并結(jié)合應(yīng)力轉(zhuǎn)移原理，連拱隧道施工的優(yōu)選方案為先開挖深埋深的洞室，后開挖淺埋深的洞室．

[1]王軍.連拱隧道施工過程的力學(xué)分析[D].同濟(jì)大學(xué),2003.

[2]朱漢華,傅鶴林,潘明軍,范臻輝.公路連體隧道的設(shè)計(jì)與施工[J].中國(guó)鐵道科學(xué),2003(8):96-100.

[3]陳秋南,張永興,劉新榮,包太.連拱隧道的設(shè)計(jì)理論與動(dòng)態(tài)施工力學(xué)研究[J].地下空間與工程學(xué)報(bào),2005(2):156-160.

[4]關(guān)寶樹.隧道工程施工要點(diǎn)集[M].北京:人民交通出版社, 2003.

[5]林剛.連拱公路隧道的合理施工方法研究[D].西南交通大學(xué), 2002.

[6]何川,林剛,張志強(qiáng),等.連拱隧道的主要修建技術(shù)問題和對(duì)策[J].現(xiàn)代隧道技術(shù),2002(39):265-269.

[7]劉洪州,黃倫海.連拱隧道設(shè)計(jì)施工技術(shù)研究現(xiàn)狀[J].西部探礦工程,2001(1):54-55.

[8]馬萬權(quán),程崇國(guó).關(guān)于連拱隧道建設(shè)技術(shù)問題的思考[J].公路, 2003(10):29-32.

[9]謝康和,周鍵.巖土工程有限元分析理論與應(yīng)用[M].北京:科學(xué)出版社,2002.

[10]梁曉丹,劉剛,趙堅(jiān).地下工程壓力拱拱體的確定與成拱分析[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,33(3):314-317.

[11]鄭西貴,李學(xué)華,趙立新,呂衛(wèi)東.巷外掘巷保護(hù)巷道的應(yīng)力轉(zhuǎn)移技術(shù)[J].礦山壓力與頂板管理,2005(4):46-47.

(責(zé)任編校：陳智全)

Numerical Simulation Study on Construction Process of Partial Pressure in JointArch Tunnel

LI Miao1,XU Xia-fei2,FU Gui-hai1
(1.College of Civil Engineering,Hunan City University,Yiyang,Hunan 413000,China;2.Wugang Jingzhou Highway Construction and Development Co.,Ltd.,Shaoyang,Hunan 422400,China)

Advancing middle drift construction process of Si-xiao joint arch tunnel in Yunnan is simulated with FEM,and then the structural mechanics behavior in joint arch tunnel construction under partial pressure is studied.The results under different excavation sequences which are the surrounding rock stress variation,plastic area development,stress of tunnel’s mid-partition and lining are got through dynamic numerical simulation.By means of comparatively analyzing simulation results on key steps of various construction schemes and stress transfer effect on tunnel excavation,optimized construction scheme of partial pressure joint arch tunnel finally is obtained.

partial pressure joint arch tunnel;construction process;stress transfer;numerical simulation

TU455

A

10.3969/j.issn.1672-7304.2015.01.001

1672–7304(2015)01–0001–06

2015-01-10

湖南省教育廳一般科技計(jì)劃資助項(xiàng)目(14C0214)；湖南省教育廳青年科技資助項(xiàng)目(13B010)

李苗(1980-)，男，漢族，湖南益陽(yáng)人，講師，主要從事橋梁健康監(jiān)測(cè)及損傷識(shí)別的研究．