黃鋒，朱合華，徐前衛(wèi)

(1.重慶交通大學(xué) 土木工程學(xué)院，重慶 400074；

2.重慶交通大學(xué) 山區(qū)橋梁與隧道國家重點實驗室培育基地，重慶 400074；

3.同濟(jì)大學(xué) 土木工程學(xué)院，上海 200092；

4.同濟(jì)大學(xué) 交通運輸工程學(xué)院，上海 200092)

?

不同地應(yīng)力場軟巖隧道漸進(jìn)破壞試驗與分析

黃鋒1,2，朱合華3，徐前衛(wèi)4

(1.重慶交通大學(xué) 土木工程學(xué)院，重慶 400074；

2.重慶交通大學(xué) 山區(qū)橋梁與隧道國家重點實驗室培育基地，重慶 400074；

3.同濟(jì)大學(xué) 土木工程學(xué)院，上海 200092；

4.同濟(jì)大學(xué) 交通運輸工程學(xué)院，上海 200092)

摘要:針對軟弱地層中常常發(fā)生的隧道塌方事故，以Ⅳ級圍巖深埋公路隧道為對象，利用模型試驗方法研究特定應(yīng)力場下無支護(hù)隧道圍巖的漸進(jìn)破壞機(jī)理。采用彈塑性損傷本構(gòu)模型，對模型試驗工況進(jìn)行有限元數(shù)值模擬，并在此基礎(chǔ)上，對不同初始地應(yīng)力場的情況進(jìn)行擴(kuò)展分析。研究結(jié)果表明：1)隧道開挖后圍巖的破壞區(qū)主要集中在隧道拱頂上方，數(shù)值計算中用損傷變量最大值表示的破壞區(qū)與模型試驗吻合較好；2)隨著側(cè)向壓力系數(shù)的增大，松動破壞區(qū)面積隨之?dāng)U大，且主要區(qū)域由隧道邊墻兩側(cè)逐漸轉(zhuǎn)變?yōu)樗淼拦绊敗⒐暗讌^(qū)域；3)隨著側(cè)向壓力系數(shù)的增大，邊墻水平收斂位移逐漸減小，拱頂收斂位移逐漸增大。

關(guān)鍵詞：軟弱圍巖；初始地應(yīng)力；漸進(jìn)破壞；模型試驗；數(shù)值模擬

在軟弱破碎巖體中修建隧道工程，因其賦存介質(zhì)的多變性、不可預(yù)見性而難以駕馭和控制[1]，圍巖失穩(wěn)造成的塌方事故已成為隧道施工中最常見的災(zāi)害現(xiàn)象之一[2]。事實上，隧道塌方是圍巖破壞由量變到質(zhì)變的表現(xiàn)，只有清楚地認(rèn)識其漸進(jìn)形成機(jī)制，才能制定出合理的防治措施。同時，初始地應(yīng)力場是引起隧道開挖后圍巖應(yīng)力重分布的最直接外因，對分析圍巖穩(wěn)定性至關(guān)重要[3]。因此，開展不同初始應(yīng)力場條件下軟弱圍巖隧道塌方破壞的漸進(jìn)形成機(jī)制，對隧道塌方事故的防治具有重要意義。目前，對隧道圍巖漸進(jìn)性破壞問題的研究，主要集中在以下幾方面：1)室內(nèi)試驗方面，以相似模型試驗為主，大多數(shù)試驗是在自重應(yīng)力場下進(jìn)行的漸進(jìn)性破壞分析，李英杰等[4]利用相似模型試驗對深埋隧道圍巖漸進(jìn)性破壞特征進(jìn)行研究；在不同地應(yīng)力條件下的情況則重點研究了圍巖及支護(hù)的開挖力學(xué)效應(yīng)，來弘鵬等[5]研究了開挖對不同應(yīng)力場公路隧道圍巖應(yīng)力分布的影響，分析了不同初始應(yīng)力場公路隧道結(jié)構(gòu)的受力特點。2)數(shù)值模擬方面，以有限元為主的連續(xù)體分析方法為主，Sterpi等[6]采用有限元數(shù)值方法對淺埋隧道圍巖穩(wěn)定性和漸進(jìn)性破壞機(jī)理進(jìn)行了模擬分析。由于難以建立合適的巖體破壞準(zhǔn)則，再加上材料參數(shù)難以確定，僅僅依靠數(shù)值計算往往不易獲得可信結(jié)果[7]。因而，采用模型試驗與數(shù)值計算相結(jié)合的方法成為了一種較好的選擇，王明年等[8]通過大比例尺模型試驗和有限元方法，對三車道公路隧道在不同構(gòu)造應(yīng)力作用下的力學(xué)行為進(jìn)行了研究?；谏鲜稣J(rèn)識，首先，以深埋IV級三車道公路隧道圍巖為對象，采用模型試驗方法模擬了特定初始地應(yīng)力場條件下圍巖的漸進(jìn)破壞過程。然后，采用有限元數(shù)值分析方法，對模型試驗工況進(jìn)行數(shù)值模擬對比分析。在此基礎(chǔ)上，對不同初始地應(yīng)力場的情況進(jìn)行擴(kuò)展分析，探討在不同側(cè)向壓力系數(shù)條件下隧道圍巖的漸進(jìn)破壞機(jī)理和應(yīng)力、應(yīng)變場特征。

1模型試驗方法

1.1　相似比與材料

物理模型試驗選用幾何相似比為Cl=50，根據(jù)相似原理推導(dǎo)可以得到，容重相似比Cγ=1；泊松比、應(yīng)變、摩擦角相似比為Cμ=Cε=Cφ=1；強(qiáng)度、應(yīng)力、黏聚力、彈模的相似比為CRc=CRt=Cσ=Cc=CE=Cc=50。

模型試驗以容易生隧道塌方的Ⅳ級圍巖為對象，其力學(xué)參數(shù)可根據(jù)《公路隧道設(shè)計規(guī)范》選取，見表1中的原型，模型材料的力學(xué)參數(shù)見表1中的模型。經(jīng)過多次反復(fù)試驗，最終研制出滿足要求的IV級圍巖模型材料及配合比為，重晶石粉∶砂∶水∶淀粉=12∶4∶2∶0.92∶0.35。

表1?、艏墖鷰r物理力學(xué)性質(zhì)指標(biāo)

1.2　模型試驗裝置

模型試驗采用了較為簡單的平面應(yīng)變試驗，按照幾何相似比隧道斷面尺寸為0.18 m×0.32 m，隧道兩側(cè)及底部的圍巖選取2倍以上隧道洞徑以減小邊界效應(yīng)，試驗箱尺寸為1.6 m×1.3 m×0.4 m。整個試驗系統(tǒng)由試驗土箱、加載系統(tǒng)、開挖裝置和數(shù)據(jù)采集等子系統(tǒng)構(gòu)成，如圖1所示。試驗?zāi)Ｐ团_用兩組180工字鋼以及25號加勁槽鋼對模型槽前后進(jìn)行約束；箱體前、后2個面為厚2 cm的有機(jī)玻璃板組成，具有足夠剛度的同時還可以觀察土體內(nèi)部的位移變化情況；為了減少模型地層和箱壁的摩擦阻力，在箱體內(nèi)壁貼上一層Teflon薄膜，并涂上一層潤滑油脂；隧道開挖過程，通過預(yù)埋一個可拆卸的特殊裝置來實現(xiàn)；隧道變埋深的情況，可以通過液壓千斤頂在模型表面施加上覆荷載P來實現(xiàn)，二者間關(guān)系如下：

(1)

式中：H是實際隧道埋深，m；γ是巖體容重，kN/m3；P是地表超載，kPa。

圖1　模型試驗系統(tǒng)Fig.1 Model test system

1.3　測試內(nèi)容與步驟

試驗中測量的巖體內(nèi)部的壓力，采用箔式微型壓力盒測量，主要布置在拱頂上方，如圖2所示，主要測試圍巖徑向、周向壓力的改變值，正的表示壓力增加，負(fù)的表示壓力減小。開挖前測量預(yù)埋壓力盒的讀數(shù)，根據(jù)壓力值的大小可初步估計地層初始地應(yīng)力狀態(tài)，作為后面數(shù)值分析的重要依據(jù)。地層表面的位移變形通過位移計讀取，內(nèi)部位移通過圍巖材料中的染色彩沙和試驗裝置上的標(biāo)志線之間的差值來獲得。試驗全過程采用數(shù)字照相機(jī)和人工素描的方式，記錄分析圍巖開挖后的松動破壞范圍。

圖2　測點布置示意圖Fig.2 Layout schematic diagram of measurement items

模型地層采用夯實填筑法制作，其基本流程如下：材料準(zhǔn)備→均勻拌合材料→用落雨法分層攤鋪→按密度控制夯實材料→在測點位置相應(yīng)埋設(shè)壓力盒。試驗加載方式采用“先自重固結(jié)后開挖，再加載至破壞”的試驗方法，即模型材料在自重作用下開挖，至圍巖變形穩(wěn)定后，再施加表面附加荷載，至圍巖破壞形成穩(wěn)定塌落拱。

2數(shù)值模擬方法

2.1　本構(gòu)模型與參數(shù)

為了進(jìn)行模型試驗對比和擴(kuò)展分析，數(shù)值模擬選用有限元軟件Abaqus進(jìn)行。本構(gòu)模型采用了基于Drucker-Prager屈服準(zhǔn)則的損傷模型，該模型的特點是適合于分析材料應(yīng)變軟化和漸進(jìn)性破壞問題[9]。與傳統(tǒng)彈塑性本構(gòu)相比，該損傷本構(gòu)的區(qū)別在計算過程中基于應(yīng)變等效原理進(jìn)行應(yīng)力更新和材料彈性模量弱化，其表達(dá)式如下：

(2)

(3)

根據(jù)Aydan等[11]的研究成果可知，材料破壞過程中各階段的變形間存在以下關(guān)系：

(4)

式中：彈性應(yīng)變極限εe和圍壓σc≈σ3可以在第1個計算子步(常規(guī)彈性塑性增量迭代)中產(chǎn)生。

數(shù)值計算采用的材料力學(xué)參數(shù)，采用模型試驗獲得的相似材料力學(xué)參數(shù)，如表1所示。

2.2　數(shù)值模型與方案

數(shù)值計算采用的有限元模型與物理模型試驗的情況一致，如圖3 (a)所示：模型左右、前后和下邊界設(shè)為均為法向位移約束，上表面為自由邊界；通過施加表面上覆荷載和側(cè)向壓力系數(shù)的方式來實現(xiàn)初始地應(yīng)力場的模擬。為了后面結(jié)果分析的方便，在隧道拱頂、拱肩、邊墻、拱腳和拱底位置沿徑向設(shè)置了測線AA’～ EE’如圖3 (b)所示。

為了研究方便，將初始地應(yīng)力場用側(cè)向壓力系數(shù)k來表示。模型試驗開始前測量拱頂上方各測點位置應(yīng)力值，通過水平應(yīng)力與垂直應(yīng)力之比，確定出地層初始平均側(cè)向壓力系數(shù)k約為1.20，如表2所示。這與《公路隧道設(shè)計規(guī)范(JTG D70—2004)》的建議取值大致相同：側(cè)向壓力系數(shù)k一般為0.5~5.0，大多數(shù)為0.8~1.5[5]。

(a)數(shù)值計算模型;(b) 模型邊界及測線布置圖3　模擬模式試驗工況的有限元模型Fig.3 Finite element model for physical model tests case

表2　壓力盒初始讀數(shù)與平均側(cè)向壓力系數(shù)

數(shù)值模擬方案如表3所示，包括2組工況，1)數(shù)值模擬與模型試驗對比，分析不同附加荷載條件下圍巖的破壞模式，k=1.20；2)在附加荷載不變的情況下，分析不同側(cè)壓力系數(shù)條件下的圍巖漸進(jìn)破壞，共模擬了k=0.54(按金尼克假設(shè)k=μ/(1-μ))，1.00，1.20和1.50共4種情況。

表3　數(shù)值模擬方案

3變荷載下圍巖漸進(jìn)性破壞分析

3.1　漸進(jìn)性破壞過程分析

隧道完成開挖后，圍巖初始階段的變形情況如圖4(a)所示，可見拱頂部分有較大沉降，圍巖侵入隧道輪廓面界線內(nèi)；隨著上覆荷載P的增大，圍巖首先在隧道邊墻、拱腰附近出現(xiàn)局部破壞，但之后破壞范圍向拱頂區(qū)域擴(kuò)展，如圖4(b)和4(c)所示；隨著上覆荷載P的繼續(xù)增大，隧道圍巖破壞區(qū)主要集中在拱頂附近，且呈拱形漸進(jìn)擴(kuò)大如圖4(d)和4(e)所示；最終上覆荷載為P=36 kPa，裂縫發(fā)展至模型表面以及邊墻內(nèi)約2 cm，隧道可見塌落拱高度約11 cm，寬度與隧道跨度相當(dāng)，如圖4(f)所示?？偟膩碚f，隧道破壞區(qū)域主要集中在拱頂上方區(qū)域，邊墻和拱底下方有局部破壞，塌落區(qū)內(nèi)的圍巖自重是隧道松動荷載的主要來源[10]。

(a)P=0 kPa；(b)P=3 kPa；(c)P=12 kPa；(d)P=36 kPa圖4　不同上覆荷載P作用下的圍巖漸進(jìn)性破壞過程Fig.4 Progressive failure of tunnel as overload increasing

不同上覆荷載P作用下，圍巖損傷變量的演化過程如圖5 (a)~5(c)所示?？偟膩碚f，圍巖損傷變量的演化過程與模型試驗類似，隧道拱腳處應(yīng)力集中較大首先出現(xiàn)了破壞，隨著模型表面超載P的增大，破壞區(qū)向拱頂發(fā)展，最后破壞區(qū)域主要集中在拱肩和拱頂附近。上覆荷載P=36 kPa時，數(shù)值計算與模型試驗的破壞區(qū)對比分析，如圖5 (d)所示。可以看出，二者在仰拱以上部分無論是破壞區(qū)的面積大小還是形狀都是相似的，表明基于損傷模型的數(shù)值計算結(jié)果，具有一定的合理性和可靠性；模型試驗中由于側(cè)壁的摩擦作用導(dǎo)致上覆荷載不能順利下傳，以至隧道開挖后的仰拱底部的破壞區(qū)不明顯。

(a)P=3 kPa；(b)P=12 kPa；(c)P=36 kPa；(d)P=36 kPa圖5　不同荷載下圍巖損傷區(qū)演化分析Fig.5 Distribution of damage zone under different overloads

3.2　圍巖應(yīng)力變化分析

圖6是隧道邊墻附近位置測點圍巖應(yīng)力增量的變化情況，其中T1表示的是隧道周向應(yīng)力增量Δσθ，R1表示的是隧道徑向應(yīng)力增量Δσr。隨著上覆荷載P增加，隧道周向應(yīng)力隨之增大，而徑向應(yīng)力隨之減小，這表明距離隧道邊墻100 mm的位置處應(yīng)該在圍巖應(yīng)力松動區(qū)范圍內(nèi)[12]。

圖6　隧道邊墻附近測點應(yīng)力增量隨上覆荷載的變化Fig.6 Stress increments of measure data near tunnel hance as overload increasing

圖7是隧道拱頂上方圍巖豎向(/隧道徑向)應(yīng)力增量Δσr隨上覆荷載的變化情況。隧道開挖后各測點值均處于減小狀態(tài)，離拱頂較近的點(R=100mm)一直是處于減小狀態(tài)，這可能與其處于松動區(qū)范圍之內(nèi)有關(guān)；離拱頂較遠(yuǎn)處的點(R=250mm)在拱腰開始破壞后急劇上升，表明應(yīng)力調(diào)整后可能處于圍巖壓力拱區(qū)域內(nèi)；離拱頂最遠(yuǎn)處的點(R=400mm)在拱頂出現(xiàn)較大塌落拱之后有所增加，但增幅相對較小。

圖8是隧道拱頂上方水平方向(/隧道周向)應(yīng)力增量Δσθ隨上覆荷載的變化情況。各測點處的應(yīng)力同樣受隧道圍巖應(yīng)力重分布的影響，每當(dāng)圍巖產(chǎn)生裂縫或有塌落時，應(yīng)力值便會有突然的小幅增加。距拱頂不同深度處點的應(yīng)力均呈增大趨勢，其中R=250mm處的測點增大最多，表明該位置處應(yīng)該在隧道壓力拱區(qū)(應(yīng)力承載區(qū))范圍內(nèi)[12]。

圖7　拱頂上方各點豎向應(yīng)力增量隨P的變化Fig.7 Vertical stress increments of measure data over tunnel crown as overload increasing

圖8　拱頂上方水平應(yīng)力增量隨P的變化Fig.8 Horizontal stress increments of measure data over tunnel crown as overload increasing

距離隧道拱頂250 mm處(測點T3)，圍巖周向應(yīng)力的數(shù)值模擬與試驗結(jié)果對比分析，如圖9所示。可以看出，總體上二者的結(jié)果較為吻合，上覆荷載較小時，數(shù)值模擬與試驗結(jié)果更為接近；上覆荷載較大時，由于圍巖破壞導(dǎo)致模型試驗中應(yīng)力測點松動，圍巖應(yīng)力傳遞困難所得結(jié)果偏小。

圖9　T3測點應(yīng)力的數(shù)值計算與模型試驗對比Fig.9 Stress comparison of measure point T3 between numerical calculation and model test

數(shù)值計算中圍巖內(nèi)部沿不同測線方向上周向應(yīng)力增量的變化規(guī)律，如圖10所示?？梢钥闯鲋芟驊?yīng)力變化的分布規(guī)律，隧道開挖卸載后松動區(qū)內(nèi)應(yīng)力降低(Δσθ<0)，應(yīng)力不斷地向圍巖深部轉(zhuǎn)移，在壓力拱區(qū)應(yīng)力升高(Δσθ>0)，最后又開始下降逐漸過渡到原巖應(yīng)力區(qū)(Δσθ→0)[12]。

圖10　不同測線周向應(yīng)力增量Fig.10 Distribution of tangential stress increment along different measure lines

3.3　圍巖變形規(guī)律分析

圖11是模型表面上覆荷載P=36 kPa條件下，圍巖位移矢量等值圖(為了研究方便，這里以柱坐標(biāo)顯示)?？梢钥闯觯瑖鷰r開挖后總體上向洞內(nèi)發(fā)生收縮變形，其分布形態(tài)與損傷區(qū)類似在拱頂附近出現(xiàn)拱形梯度帶，位移最大值出現(xiàn)在隧道拱頂上方，約1.4 mm；其余部分的變形量相對較小。

圖11　P=36 kPa條件下圍巖徑向位移等值圖Fig.11 Radical displacement isogram of surrounding rock mass in the case of P=36 kPa

模型表面不同上覆荷載條件下，模型材料表面位移沿橫向分布的對比分析，如圖12所示。可以看出，數(shù)值計算和模型試驗的結(jié)果總體上較為吻合，且隨著上覆荷載的增加，二者間的差異越小；曲線均呈馬鞍形分布，這與經(jīng)典的Peck曲線相似[13]，即隧道開挖后拱頂正上方出現(xiàn)最大沉降，距離隧道軸線越遠(yuǎn)沉降越小。

圖12　模型表面沉降的對比分析Fig.12 Settlements comparison analysis along model surface

4不同地應(yīng)力下圍巖力學(xué)效應(yīng)分析

4.1　圍巖破壞模式分析

圖13是模型表面上覆荷載P=36kPa時，在不同側(cè)壓系數(shù)k情況下的圍巖損傷變量分布圖。隨著側(cè)向壓力系數(shù)k的增大，圍巖損傷、破壞區(qū)由集中在隧道邊墻兩側(cè)區(qū)域(k≤1)，逐漸轉(zhuǎn)變?yōu)橹饕植荚谒淼拦绊?、拱底區(qū)域(k>1)。可見側(cè)向壓力系數(shù)大小對于圍巖的破壞模式有著重要影響。

(a)k=0.54；(b)k=1.00；(c)k=1.20；(d)k=1.50圖13　不同側(cè)向壓力系數(shù)下圍巖開挖損傷區(qū)分布Fig.13 Distribution isogram of damage zones under different lateral pressure coefficients

不同側(cè)向應(yīng)力系數(shù)k條件下，圍巖損傷、破壞區(qū)的對比分析，如表4所示(為了分析方便用面積比值表示)?？梢钥闯?，圍巖損傷區(qū)的面積遠(yuǎn)大于破壞區(qū)，且面積都隨著側(cè)向壓力系數(shù)的增大而增大，但破壞區(qū)的增加速率更快。

表4　不同側(cè)向應(yīng)力系數(shù)下?lián)p傷、破壞區(qū)面積Table 4 Values of damage zones and failure zones under different lateral pressure coefficients

4.2　圍巖應(yīng)力場特征分析

不同側(cè)向應(yīng)力系數(shù)k條件下，隧道開挖后圍巖周向應(yīng)力增量沿邊墻CC’測線的分布規(guī)律，如圖14所示?？梢钥闯?，不同側(cè)向壓力系數(shù)下，圍巖周向應(yīng)力增量幾乎都呈現(xiàn)出先增大后減小的規(guī)律；且側(cè)向應(yīng)力系數(shù)越大，圍巖周向應(yīng)力增量峰值越小。

圖14　隧道開挖后圍巖周向應(yīng)力增量沿CC’測線分布Fig.14 Distribution of surrounding rock mass stress increments after tunneling along measure line CC’

圖15是不同側(cè)壓系數(shù)情況下在模型表面上覆附加荷載P=36 kPa時，隧道圍巖內(nèi)周向應(yīng)力增量Δσθ的分布等值圖，其中，紅色區(qū)域為應(yīng)力松動區(qū)，藍(lán)色區(qū)域為壓力拱承載區(qū)。可以看出，不同側(cè)向壓力系數(shù)情況下，圍巖松動區(qū)、壓力拱的形狀和大小均不同；松動區(qū)集中在隧道輪廓面一定深度范圍內(nèi)，且形態(tài)與損傷區(qū)相似。

4.3　圍巖應(yīng)變場特征分析

不同側(cè)向壓力系數(shù)條件下在模型表面上覆附加荷載P=36 kPa時，圍巖最終徑向位移分布如圖16所示。在側(cè)向壓力系數(shù)k≤1情況下，圍巖徑向位移最大值出現(xiàn)在隧道邊墻附近；在側(cè)向壓力系數(shù)k>1情況下，最大值出現(xiàn)在隧道拱頂。

表5是圍巖變形在不同側(cè)向應(yīng)力系數(shù)k的條件下的計算結(jié)果?？梢钥闯觯S著側(cè)向壓力系數(shù)k的增大，圍巖應(yīng)力環(huán)境發(fā)生改變，邊墻水平收斂位移逐漸減小，拱頂收斂位移逐漸增大；側(cè)向壓力系數(shù)k≤1情況下，圍巖徑向位移最大值出現(xiàn)在隧道邊墻，在側(cè)向壓力系數(shù)k>1情況下，最大值出現(xiàn)在隧道拱頂；側(cè)向壓力系數(shù)k=1.2的情況下，隧道表面變形相對均勻，邊墻和拱頂?shù)淖冃谓咏?/p>

(a)k=0.54；(b)k=1.00；(c)k=1.20；(d)k=1.50圖15　不同側(cè)向應(yīng)力系數(shù)下周向應(yīng)力增量圖(P=36 kPa)Fig.15 Isograms of tangential stress increment under different lateral pressure coefficients (P=36 kPa)

(a)k=0.54；(b)k=1.00；(c)k=1.20；(d)k=1.50圖16　不同側(cè)向應(yīng)力系數(shù)下圍巖徑向應(yīng)力分布圖Fig.16 Radical displacement isograms of surrounding rock mass different lateral pressure coefficients

表5　不同側(cè)向應(yīng)力系數(shù)下豎向、水平收斂變形

5結(jié)論

1)隧道開挖后，圍巖的損傷破壞區(qū)主要集中在隧道拱頂上方，伴隨著邊墻和拱底位置有局部破壞，這是隧道松動荷載的主要來源，隨著模型表面上覆荷載的增加圍巖破壞區(qū)的范圍隨之?dāng)U大。

2)圍巖損傷區(qū)、破壞區(qū)都隨著側(cè)向壓力系數(shù)的增大而增大，且由集中在隧道邊墻兩側(cè)逐漸轉(zhuǎn)變?yōu)樗淼拦绊?、拱底區(qū)域。

3)根據(jù)隧道開挖后的周向應(yīng)力增量，周邊地層可劃分為松動區(qū)→壓力拱區(qū)→原巖應(yīng)力區(qū)，且不同側(cè)向壓力系數(shù)下，圍巖松動區(qū)、壓力拱的形狀和大小均不同，其規(guī)律與損傷變量分布相似。

4)隧道開挖后圍巖向洞內(nèi)發(fā)生收明顯縮變形，隨著側(cè)向壓力系數(shù)的增大，圍巖邊墻水平收斂位移逐漸減小，拱頂收斂位移逐漸增大。

5)模型試驗中受到測量技術(shù)的限制，對圍巖變形的分析較少，下一步考慮利用數(shù)字照相或散斑技術(shù)，探討巖體應(yīng)變與漸進(jìn)破壞之間的關(guān)系。

參考文獻(xiàn)：

[1]Seokwon J, Jongwoo K, Youngho S, et al.Effect of a fault and weak plane on the stability of a tunnel [C]// Regional ISRM Sinorock 2004 Symposium, China: Int J Rock Mech Min Sci: 1-6.

[2] 李釗.隧道坍方突發(fā)事件風(fēng)險原因統(tǒng)計及范例推理[J].鐵道科學(xué)與工程學(xué)報, 2009, 6(4): 54-58.

LI Zhao.Statistics of tunnel collapse accidents and its assessment based on case reasoning [J].Journal of Railway Science and Engineering.2009, 6(4): 54-58.

[3] 羅晶，彭立敏，施成華，等.不同埋深及偏壓角度條件下隧道力學(xué)特性[J].鐵道科學(xué)與工程學(xué)報, 2012, 9(4): 75-78.

LUO Jing， PENG Limini， SHI Chenhua, et al.Stress characteristics of tunnel under different buried depths and bias angle conditions [J].Journal of Railway Science and Engineering.2012, 9(4): 75-78.

[4] 李英杰，張頂立，宋義敏，等.軟弱破碎深埋隧道圍巖漸進(jìn)性破壞試驗研究[J].巖石力學(xué)與工程學(xué)報，2012, 3(16): 1138-1147.

LI Yingjie, ZHANG Dingli, SONG Yimin, et al.Experimental research of progressive damage of surrounding rock for soft fractured deep tunnel [J].Chinese Journal of Rock Mechanics and Engineering, 2014, 2012, 3(16): 1138-1147.

[5] 來弘鵬，林永貴，謝永利，等.不同應(yīng)力場下軟弱圍巖公路隧道的力學(xué)特征試驗[J].中國公路學(xué)報，2008, 2(14): 81-87.

LAI Hongpeng, LIN Yonggui, XIE Yongli, et al.Mechanical characteristic tests of soft-weak surrounding rock highway tunnel under different stress fields [J].China Journal of Highway and Transport, 2008,2(14): 81-87.

[6] Sterpi D, Cividini A.A physical and numerical investigation on the stability of shallow tunnels in strain softening media[J].Rock Mechanics and Rock Engineering, 2004, 37(4): 277-298.

[7] 周太全，華淵，連俊英，等.軟弱圍巖隧道施工全過程非線性有限元分析[J].巖土力學(xué)，2004，25(增2)：338-342.

ZHOU Taiquan, HUAYuan, LIAN Junying， et al.Nonlinear finite element analysis for railway tunnel construction process within soft and weak rock[J].Rock and Soil Mechanics, 2004, 25(Supp 2)：338-342.

[8] 王明年，關(guān)寶樹，何川.三車道公路隧道在不同構(gòu)造應(yīng)力作用下的力學(xué)行為研究[J].巖土工程學(xué)報，1998，20(1)：51-55.

WANG Mingnian, GUAN Baoshu, HE Chuan.Study on mechanical behavior of three lane road tunnel under tectonic stresses [J].Chinese Journal of Geotechnical Engineering, 1998, 20(1): 51-55.

[9] 朱合華，黃 鋒，徐前衛(wèi).變埋深下軟弱破碎隧道圍巖漸進(jìn)性破壞試驗與數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報，2010, 29(6): 1113-1122.

ZHU Hehua, HUANG Feng, XU Qianwei.Model test and numerical simulation for progressive failure of weak and fractured tunnel surrounding rock under different overburden depths [J].Chinese Journal of Rock Mechanics and Engineering, 2010, 29(6): 1113-1122.

[10] Huang F, Zhu H H, Xu Q W, et al.The effect of weak interlayer on the failure pattern of rock mass around tunnel Scaled model tests and numerical analysis [J].Tunnelling and Underground Space Technology, 2013, 35(4): 207-218.

[11] Aydan ?, Akagi T, Kawamoto T.The squeezing potential of rocks around tunnels theory and prediction [J].Rock Mechanics and Rock Engineering, 1993, 26(2): 137-163.

[12] Huang Z P, Broch E, Lu M.Cavern roof stability mechanism of arching and stabilization by rock bolting[J].Tunneling and Underground Space Technology, 2002, 17(3), 249-261.

[13] Peck R B.Deep excavations and tunnelling in soft ground [C]// Proceedings of the 7th International Conference on Soil Mechanics and Foundation Engineering.Mexico: [s.n.], 1969: 225-290.

(編輯陽麗霞)

The test and analysis of progressive failure for soft-weaksurrounding rock mass tunnel in different ground stress fields

HUANG Feng1, 2, ZHU Hehua3, XU Qianwei4

(1.School of Civil Engineering, Chongqing Jiaotong University, Chongqing 400074，China；

2.State Key Laboratory Breeding Base of Mountain Bridge and Tunnel Engineering, Chongqing Jiaotong University, Chongqing 400074，China；

3.College of Civil Engineering,Tongji University，Shanghai 200092，China；

4.College of Transportation Engineering, Tongji University，Shanghai 200092，China)

Abstract:According to the common tunnel collapse inner soft-weak ground, the deep buried road tunnel with surrounding rock mass of grade Ⅳ was chosen as study object.The progressive mechanism of tunnel surrounding rock mass without support system was studied by model test.The case of model test was simulated by FEM with an elasto-plastic damage constitutive mode.On this basis, the cases of different initial ground stress field were carried out.The study shows that: 1) the main area of surrounding rock mass failure after tunneling is centralized above tunnel crown.The failure area denoted by maximum of damage variable in numerical calculation are similar with the results in model test.2) With the increase of cateral pressure coefficients, the area of loose failure becomes more simultaneously, and the main areas gradually change from tunnel hances to crown and invert.3) With the increase of cateral pressure coefficients the convergence displacement of hance becomes less and the displacement of crown becomes larger.

Key words:soft-weak surrounding rock mass; initial ground stress; progressive failure; model test; numerical simulation

通訊作者：黃鋒(1982-)，男，重慶人，副教授，博士，從事隧道及地下工程領(lǐng)域的教學(xué)與研究；E-mail：huangfeng216@126.com

基金項目：國家自然科學(xué)基金資助項目(51308574)；中國博士后科學(xué)基金項目(2014M562286)；重慶市科委基礎(chǔ)與前沿研究項目(cstc2013jcyjA30007)；重慶市教委科學(xué)技術(shù)研究項目(KJ130404)；山區(qū)橋梁與隧道工程國家重點實驗室培育基地開放基金項目(CQSLBF-Y13-4)

收稿日期：2015-09-05

中圖分類號：U45

文獻(xiàn)標(biāo)志碼：A

文章編號：1672-7029(2015)06-1412-08