吳飛潔, 邵生俊,2, 佘芳濤,2

(1.西安理工大學 巖土工程研究所， 陜西 西安 710048；2.陜西省黃土力學與工程重點實驗室， 陜西 西安 710048)

?

黃土隧道圍巖壓力的一種極限平衡理論計算方法研究

吳飛潔1, 邵生俊1,2, 佘芳濤1,2

(1.西安理工大學 巖土工程研究所， 陜西 西安 710048；2.陜西省黃土力學與工程重點實驗室， 陜西 西安 710048)

隧道圍巖壓力的計算方法研究是隧道工程支護結構設計中的關鍵問題。本文依據(jù)黃土隧道圍巖的破壞模式調查和圍巖壓力的現(xiàn)場實測結果，提出了黃土隧道破壞模式，進而按照隧道上方土體與側面楔形體的極限平衡條件，解析得到一種適用于黃土隧道的圍巖壓力計算方法，并將實際黃土隧道的圍巖壓力計算結果，與以往土質隧道的5種常用方法計算結果進行對比。研究結果表明：該方法的計算結果介于以往圍巖壓力計算結果的最大值與最小值之間，且隨著埋深增大存在峰值，該峰值可用于判斷深淺埋臨界深度，并且隨黃土強度指標變化，圍巖壓力計算結果仍處于前述變化范圍內(nèi)。

極限平衡理論； 黃土隧道； 破壞模式； 圍巖壓力

長期以來，圍巖壓力理論是隧道與其它地下工程研究領域中的一項重要研究課題。按荷載結構模型進行洞室襯砌設計時，如何合理地確定圍巖壓力是該問題的關鍵所在?，F(xiàn)有的圍巖壓力計算方法主要有全土柱理論、普氏公式、太沙基公式、比爾鮑曼公式、規(guī)范法以及計算形變圍巖壓力的卡柯公式與芬納公式[1-3]。其中，全土柱理論適用于明挖回填工程或埋深較淺、地質條件差的隧道;普氏公式的建立是基于自然平衡拱理論，即假定巖土體開挖后洞頂可以自然成拱，因此僅適用于具有一定自承載能力的深埋隧道;太沙基公式和比爾鮑曼公式主要適用于淺埋隧道?！惰F路隧道設計規(guī)范》[4]的圍巖壓力計算公式是根據(jù)不同的隧道埋深而分為超淺埋、淺埋和深埋三種情況，其中超淺埋時采用全土柱理論，淺埋時采用謝家烋[5]公式，深埋時采用基于樣本統(tǒng)計和圍巖分級的經(jīng)驗公式。對于工程技術人員，該方法便于實際應用，適用范圍相對較廣?？鹿脚c芬納公式是基于對圓形洞室圍巖形變壓力分析所提出的兩種圍巖壓力計算方法。這兩個公式中，塑性區(qū)域半徑的確定大多憑借經(jīng)驗公式，缺乏一定的理論基礎。

綜上分析可見，現(xiàn)有的圍巖壓力計算方法均具有其自身的適用性和局限性[6]，目前尚沒有一種完全適用于黃土隧道圍巖壓力的計算方法。對此，本文通過調查分析黃土隧道圍巖的破壞模式和圍巖壓力的現(xiàn)場實測結果，提出一種基于極限平衡理論的黃土隧道圍巖壓力計算方法，并通過實例計算，比較不同圍巖壓力計算方法對黃土隧道的適用性。該研究對于隧道圍巖壓力計算的理論研究和黃土隧道工程建設都具有一定的實際意義。

1　黃土隧道圍巖的破壞模式和圍巖壓力監(jiān)測結果分析

1.1黃土隧道圍巖的破壞模式

在依托鄭西客運專線大量實地調查資料[7-8]的基礎上，發(fā)現(xiàn)隧道上方的地表裂縫基本對稱分布在隧道中線兩側，且地表裂縫離隧道中線的距離隨隧道埋深增加而增大，現(xiàn)場調查的典型裂縫情況如圖1所示。這一現(xiàn)象與淺埋隧道圍巖理論破壞模式的發(fā)展趨勢一致。與此同時，也對破裂角進行了大量的現(xiàn)場調查，最終確定實際破裂角是從地表開始到隧道邊墻即最大寬度處為止，其典型斷面實際破裂角如圖2所示。

圖1　黃土隧道典型裂縫分布圖[9]Fig.1　Distribution of typical cracks of loess tunnel

圖2　典型斷面實際破裂角[9]Fig.2　Angle between typical failure and horizontal planes

1.2圍巖壓力監(jiān)測結果分析

通過分析大量的鄭西客運專線初期支護圍巖壓力的現(xiàn)場實測資料發(fā)現(xiàn)，圍巖壓力基本都表現(xiàn)為壓力，即隧道開挖后圍巖有向洞內(nèi)擠壓的趨勢，這在理論上是符合實際的。

圍巖壓力沿洞室周邊分布并不均勻，其最大值可能出現(xiàn)在洞周的任一部位，這是由洞室周圍的地層條件和圍巖強度的不確定性引起的。但總體而言，黃土[10]隧道圍巖壓力“貓耳朵”分布性狀較為明顯，圍巖壓力最大值出現(xiàn)在拱腰而非拱頂，如圖3所示(圖中Yi為各監(jiān)測點的位置，相應圍巖接觸壓力單位均為kPa)。這與前面通過實地調查確定的實際破裂角是相吻合的，即從地表開始到隧道邊墻即最大寬度處(拱腰)為止。

圖3　函谷關隧道圍巖壓力監(jiān)測分布圖[9]Fig.3　Measured surrounding rock pressure of Hanguguan tunnel

2　基于極限平衡理論的黃土隧道圍巖壓力計算方法

2.1基本原理

依據(jù)黃土隧道實地調查的破壞模式和圍巖壓力的現(xiàn)場實測結果，提出了基于極限平衡理論的黃土隧道圍巖破壞模型。由于洞室開挖,圍巖自身的承載力不足以負荷逐漸增大的圍巖壓力，引起上覆巖體向下移動，同時帶動兩側三棱體下沉，最終形成達到極限平衡條件的破壞區(qū)，如圖4所示。

圖4　新模型示意圖Fig.4　Diagram of the new model

圖中假定ECDF沿側墻豎直面FD、EC呈剛體運動，當洞室開挖時ECDF下沉，并同時帶動兩側三棱體ACE與FDK下沉，最終導致兩側楔形土體發(fā)生塑性滑移破壞。此時洞室兩側形成的破裂面與水平面的夾角為β。假定拱頂上方剛體和兩側塑性楔形體服從極限平衡條件，拱頂上方剛體對稱，豎直面上只有正應力(ML平面上剪應力等于零)，且服從靜止土壓力變化規(guī)律。

2.2計算公式推導

截取破壞區(qū)的一半作為脫離體，其作用力示于圖4中的(b)、(c)(圖中T′、N′與T″、N″為兩個塊體之間的作用力與反作用力)。分別建立塊體MLFD與塊體FDK的力平衡方程，并聯(lián)立求解得：

N=Gsinα-Psinα+Ecosα

(1)

T=Gcosα-Pcosα-Esinα

(2)

式中，E為中軸線上的靜止土壓力；T、N分別為破裂面上的切向力與法向力；G1、G2分別為塊體MLFD與塊體FDK的重力，且G=G1+G2；P為圍巖壓力；β為破裂面與水平向夾角(即破裂角， 通常β=45°+φ/2， 其中φ為黃土的摩擦角)， 且定義α=90°-β。

在洞室上覆巖體MLDK中取一厚度為dh的微元條，其沿斜邊DK方向長度為dl，洞室埋深為h，微元條的寬度為b,由圖示幾何關系可知dl=dh/cosα；則作用在微元條兩端的靜止土壓力dE、破裂面上的摩擦力dT及正壓力dN分別為：

dE=σxdh

(3)

dT=τfdl

(4)

dN=σdl

(5)

式中，σx為作用在MLDK上的水平向應力；τf、σ分別為破裂面上的抗剪強度與正應力。

因此，從地面到洞室頂部巖體中，E、T、N分別為：

(6)

(7)

(8)

式中，K0為靜止土壓力系數(shù)；σz為作用在MLDK上的豎向應力；z為沿深度方向的任意埋深；c為黃土的粘聚力(kPa)；γ為黃土的容重(kN·m-3)；H為黃土隧道的埋深(m)。

在洞室頂部土體下沉過程中，作用于洞室上的總重力G為塊體MKDL的自重，即圖4中塊體MLFD、塊體FDK的重力之和：

(9)

將E、N、T、G代入式(1)、(2)，求得圍巖壓力P：

(10)

則單位寬度作用的圍巖荷載：

(11)

若令：

(12)

pa=KγH

(13)

式中，B為黃土隧道圍巖破裂面拱頂上方的有效寬度(m)。

水平圍巖壓力采用朗肯主動土壓力理論進行計算。

3　新模型圍巖壓力合理性驗證

參考寶蘭客運專線的勘察設計資料，黃土隧道土層的容重取為18 kN/m3,內(nèi)摩擦角取為25°，粘聚力取為30 kPa，隧道跨度取為14.72 m，隧道高度取為12.58m，靜止土壓力系數(shù)K0取為0.58(按《建筑基礎工程技術規(guī)范》中的經(jīng)驗公式計算)。設隧道埋深取值為0～120 m，計算得到的新模型圍巖壓力與隧道埋深的關系如圖5所示。由圖5可以看出，新模型圍巖壓力計算結果介于比爾鮑曼公式與太沙基公式之間，且隨埋深的增加呈拋物線形式變化。圍巖壓力峰值出現(xiàn)在埋深60～70 m之間，超過峰值埋深之后，圍巖壓力迅速衰減，當隧道埋深超過一定的范圍后，圍巖壓力甚至會出現(xiàn)負值，這與工程實際不符。因此，此方法適用于淺埋隧道圍巖壓力計算，其圍巖壓力峰值對應的埋深是深淺埋臨界埋深。

圖5　新模型計算圍巖壓力與隧道埋深的關系Fig.5　Relationship between the surrounding rock pressure via the new model and the bedded depth of tunnel

為了進一步驗證新模型理論的合理性，計算不同c值、φ值下的圍巖壓力，其變化趨勢分別如圖6、圖7所示。

圖6　不同黏聚力條件下的新模型圍巖壓力Fig.6　Surrounding rock pressure of the new model under different conditions of cohesion

圖7　不同內(nèi)摩擦角條件下的新模型圍巖壓力Fig.7　Surrounding rock pressure of the new model under different conditions of internal frictional angle

由圖6可知，當c值增大時，圍巖壓力逐漸減小，與圍巖越穩(wěn)定圍巖壓力越小這一客觀規(guī)律吻合，表明了新模型的合理性。同時，不管c值怎么變化，峰值點之前的圍巖壓力總是介于全土柱理論與太沙基公式之間，即現(xiàn)有圍巖壓力計算方法結果的最大值與最小值之間，符合黃土隧道圍巖壓力量值的大小及其變化規(guī)律。

圖7對比了不同內(nèi)摩擦角的圍巖壓力變化趨勢，可以看出，隨著內(nèi)摩擦角的增大，圍巖壓力逐漸減小，這一現(xiàn)象符合圍巖壓力隨內(nèi)摩擦角變化的客觀規(guī)律?？傮w上看，不論φ值變化范圍多大，峰值點之前的圍巖壓力總是介于全土柱理論與太沙基公式之間，符合黃土隧道圍巖壓力量值的大小及其變化規(guī)律。因此，新模型對淺埋黃土隧道圍巖壓力的計算具有明顯的合理性。

與此同時，新模型提供了一種新的界定黃土隧道深淺埋的方法。本文認為曲線的拐點，即為界定隧道深淺埋的分界點，與以往深淺埋的分界點相比，該分界點處隧道埋深稍微偏大，故在計算圍巖壓力時，偏于保守。但考慮到黃土隧道地區(qū)普遍存在的不良地質條件及頻繁發(fā)生的隧道工程問題，建議采用。

圖8給出了粘聚力與深淺埋臨界深度的關系。針對同一種黃土，不同初始濕度的粘聚力變化較大，故含水率的變化將引起深淺埋臨界值的改變。隨著含水率的減小，粘聚力逐漸增大，則深淺埋臨界深度減小，這與實際工程一致。

圖8　粘聚力與深淺埋臨界埋深的關系Fig.8　Relationship between the critical bedded depth of shallow and deep tunnel and the cohesion

4　結　論

1) 本模型依據(jù)黃土隧道實地調查的破壞模式和圍巖壓力的現(xiàn)場實測結果，在推導公式時，綜合考慮了隧道有效寬度、隧道埋深、黃土力學特性參數(shù)、破裂角、靜止土壓力系數(shù)以及容重對圍巖壓力的影響，較為全面地反映了各相關因素對黃土隧道圍巖壓力的影響。因此，本式可以更好地反映黃土隧道的破壞機理，從而使計算結果更為合理。

2) 新模型圍巖壓力計算結果介于比爾鮑曼公式與太沙基公式之間，且隨埋深的增加呈拋物線形式變化。當隧道埋深超過某一深度時，圍巖壓力為負值，這與工程實踐相違背。因此，此方法僅適用于淺埋隧道圍巖壓力計算。

3) 本式計算得出的圍巖壓力峰值可用于判斷深淺埋臨界深度，隨著黃土粘聚強度的減小，深淺埋臨界深度減小。同時，隨著黃土內(nèi)摩擦角的減小，深淺埋臨界深度減小。

[1]孫鈞,侯學淵.地下結構[M].北京:科學出版社,1987.

[2]趙明階.巖石力學[M].北京：人民交通出版社，2011：8.

[3]沈明榮，陳建峰.巖體力學[M].上海:同濟大學出版社，2006.

[4]鐵道第二勘察設計院.鐵路隧道設計規(guī)范:TB10003—2005[S].北京：中國鐵道出版社，2005．

[5]謝家烋．淺埋隧道的地層壓力[J]．土木工程學報，1964，(6)：58-70．

XIE Jiaxiao. Earth pressure on shallow burial tunnel[J]. China Civil Engineering Journal, 1964，(6): 58-70.

[6]李鵬飛,周燁，伍冬.隧道圍巖壓力計算方法及其適用范圍[J].中國鐵道科學,2013,34(6):55-60.

LI Pengfei, ZHOU Ye,WU Dong. Calculation method for surrounding rock pressure and application scopes[J].China Railway Science , 2013,34(6): 55-60.

[7]王明年,郭軍，羅祿森,等.高速鐵路大斷面深埋黃土隧道圍巖壓力計算方法[J].中國鐵道科學,2009,30(5):53-58.

WANG Mingnian, GUO Jun, LUO Lusen,et al. Calculation method for the surrounding rock pressure of deep buried large sectional loess tunnel of high-speed railway[J]. China Railway Science, 2009,30(5): 53-58.

[8]王明年,郭軍，羅祿森,等.高速鐵路大斷面黃土隧道深淺埋分界深度研究[J].巖土力學,2010,31(4)：1157-1162.

WANG Mingnian, GUO Jun, LUO Lusen,et al. Study of critical buried depth of large cross-section loess tunnel for high speed railway[J]. Rock and Soil Mechanics, 2010,31(4)：1157-1162.

[9]趙勇，李國良，喻渝.黃土隧道工程[M].北京:中國鐵道出版社,2011.

[10]葛苗苗,李寧，鄭建國,等.基于蠕變試驗的黃土高填方工后沉降規(guī)律數(shù)值研究[J].西安理工大學學報,2015,31(3):295-300.

GE Miaomiao, LI Ning,ZHENG Jianguo,et al. Numerical study on the settlement regularity of loess high fill engineering based on creep tests[J]. Journal of Xi’an University of Technology, 2015,31(3):295-300.

(責任編輯周蓓)

Study of the calculation method for surrounding rock pressure in loess tunnel based on the limit equilibrium theory

WU Feijie1， SHAO Shengjun1,2， SHE Fangtao1,2

(1.Institute of Geotechnical Engineering, Xi’an University of Technology, Xi’an 710048, China; 2.Shaanxi Province Key Laboratory of Loess Mechanics and Engineering, Xi’an 710048, China)

The calculation method for surrounding rock pressure is a key problem in designing the supporting structure of tunnel engineering. Based on the failure mode of loess tunnel and the field measuring results of surrounding rock pressure, a new failure mode is proposed in this paper. According to the limit equilibrium condition of the soil mass upward loess tunnel and the soil mass on both sides, the calculation method is derived. A comparison between the new method and five common methods is conducted by calculating the surrounding rock pressure of the same loess tunnel, revealing that result by the new method is between the maximum value and the minimum value. And with the increase of the bedded depth, the peak depth of the loess tunnel is indicated, it shows the critical bedded depth of shallow and deep tunnel with the same strength characteristics. At the same time, the variation of surrounding rock pressure with the change of soil strength index is also studied, which is still between the maximum value and the minimum value of common surrounding rock pressure.

limit equilibrium theory; loess tunnel; failure mode; surrounding rock pressure

10.19322/j.cnki.issn.1006-4710.2016.03.015

2015-12-11

國家自然科學基金資助項目(41272320)；陜西省黃土力學與工程重點實驗室重點科研計劃資助項目(2010JS084)；陜西省科技廳基金資助項目(2014JQ5173)；陜西省教育廳基金資助項目(14JK1519)

吳飛潔，女，碩士生，研究方向為黃土動力學與隧道工程。E-mail：1550659940@qq.com

邵生俊，男，教授，博導，研究方向為黃土力學、黃土動力學等。E-mail：sjshao@xaut.edu.cn

TU43

A

1006-4710(2016)03-0338-05