黃明奎,全 健
(重慶交通大學 土木工程學院,重慶 400074)
考慮區(qū)間影響因素的地下洞室圍巖位移分析
黃明奎,全 健
(重慶交通大學 土木工程學院,重慶 400074)
針對影響圍巖位移的參數(shù)的不確定特點,將傳統(tǒng)圍巖位移計算模型與區(qū)間分析方法相結合,建立了適用于區(qū)間分析的地下洞室圍巖位移計算模型,利用該模型分析了IV級圍巖條件下拱頂位移的區(qū)間分布特點。結果表明:該模型客觀上更能反映圍巖位移的實際變形性態(tài),研究成果可為后續(xù)工程設計、施工以及支護時機的選擇提供參考。
隧道工程;圍巖;位移;區(qū)間分析
地下洞室圍巖周邊位移已成為評判隧道結構穩(wěn)定及安全的重要因素[1]。大量的工程實踐證明,合理的圍巖位移理論分析及可靠的現(xiàn)場位移監(jiān)測,是評判地下洞室圍巖穩(wěn)定性的關鍵。合理有效的圍巖位移計算目前已成為地下洞室穩(wěn)定性分析的重要研究內容,受到國內外學者、工程技術人員的廣泛關注。但是,現(xiàn)有的圍巖位移分析絕大部分是基于確定的參數(shù)、初始條件和邊界條件[2],這必然導致分析結果與工程實際情況存在差異。此外,由于工程巖體的隱蔽性,使得前期現(xiàn)場勘察、室內外試驗得到的圍巖力學參數(shù)與實際圍巖力學參數(shù)不吻合[3]。如果采用確定的圍巖力學參數(shù)、初始條件以及邊界條件對地下洞室圍巖位移進行分析與預測可能無法得到符合實際的位移變化規(guī)律,甚至可能嚴重偏離真實情況,進而導致對圍巖穩(wěn)定性的錯誤判斷[4]。因此,開展地下洞室圍巖位移變化的不確定性分析十分必要。
目前,對于不確定性問題的研究,主要有2種解決辦法[5]:①點估計;②區(qū)間估計。現(xiàn)有的基于確定參數(shù)、初始條件及邊界條件計算圍巖位移的單一數(shù)值解,是一種點估計方法。區(qū)間估計則是從圍巖參數(shù)的上下限來表示其不確定性變化區(qū)間,進而建立圍巖位移不確定性區(qū)間分析模型。巖體是在地質歷史中形成,由巖石單元體和結構面網絡組成的,具有一定結構并賦存于一定的天然應力和地下水等地質環(huán)境中的地質體。其力學參數(shù)的不確定性是由其自身的特點決定的,不僅與結構面、地下水狀況以及地應力有關,而且還與地下洞室的埋深有關,具有明顯的空間特性。因此,區(qū)間估計比點估計更能反映圍巖位移實際變化動態(tài),已逐漸應用于地下工程非確定性分析中。王登剛等[6]考慮現(xiàn)場觀測數(shù)據(jù)的不確定性,采用約束變尺度方法,獲得巷道圍巖地應力和彈性模量的區(qū)間反演模型;蔣沖等[4]根據(jù)巖土參數(shù)的區(qū)間性特點,建立了基于區(qū)間分析的巖土結構穩(wěn)定性非概率可靠性分析模型,并應用于工程實際,驗證了模型的可行性;張永杰等[7]考慮巖溶區(qū)路基穩(wěn)定性影響因素的模糊性、評價指標與計算參數(shù)取值的區(qū)間性,建立了巖溶區(qū)公路路基穩(wěn)定性的區(qū)間模糊評判分析方法。筆者根據(jù)地下洞室結構圍巖的特點,以彈塑性理論和區(qū)間數(shù)學理論為基礎,考慮圍巖力學參數(shù)的隨機性和變異性,建立圍巖位移的區(qū)間分析模型,使位移計算更能反映圍巖位移的實際變形性態(tài)。
1.1 區(qū)間分析
區(qū)間分析方法是在20世紀50年代末提出并發(fā)展起來的[8]。最初主要用于處理計算內的浮點運算,由于它能夠反映工程中參數(shù)的不確定性,很快被引入到工程分析中。主要分析原則如下:
值得注意的是,在區(qū)間運算中,如果同一區(qū)間變量出現(xiàn)多于一次或區(qū)間變量之間存在“相關性”區(qū)間運算,則運算結果可能導致擴張。
1.2 考慮區(qū)間影響的圍巖位移模型的構建
1.2.1 二次應力狀態(tài)下圍巖位移計算模型
根據(jù)彈塑性力學理論,可得地下洞室圍巖在二次應力狀態(tài)下洞壁周邊的最大徑向位移計算公式[8]:
(1)
式(1)可轉化為
式中:E,μ分別為地下洞室圍巖的變形模量和泊松比;λ為側壓力系數(shù)。
由式(1)可以看出,地下洞室周邊圍巖徑向位移不僅與圍巖變形模量、泊松比有關,而且還與隧道所在地的初始應力場有著密切的關系。并且圍巖的變形模量、泊松比隨洞室開挖的推進不斷變化,為不確定量。因此,采用傳統(tǒng)的單一數(shù)值解不能反映實際圍巖的變形性質,應引入不確定性方法進行研究分析。
1.2.2 區(qū)間變形模量E和和區(qū)間泊松比μ確定
地下洞室在開挖過程中,由于圍巖應力重分布而引起圍巖的變形模量和泊松比在整個過程中為一變量。因此,采用區(qū)間函數(shù)能夠反映施工過程圍巖的變形特性。為了分析簡潔,筆者以實際工程中比較典型的IV級圍巖為例進行圍巖變形特性分析。根據(jù)JTGD70—2004《公路隧道設計規(guī)范》(以下簡稱《規(guī)范》),圍巖變形模量區(qū)間E=[1.3,6]GPa,泊松比區(qū)間μ=[0.30,0.35]。
1.2.3 區(qū)間側壓力系數(shù)λ的確定
根據(jù)式(1),地下洞室圍巖洞壁位移與初始應力場有著密切關系。地下洞室在開挖過程中,由于卸荷而引起圍巖內應力場變化,即洞室周邊圍巖水平應力和豎向應力在整個施工過程中是變化的,也就是圍巖側壓力系數(shù)在整個成洞過程中為一變量。由于忽略構造應力場的影響,在洞室埋深一定的情況下采用側壓力系數(shù)來表征初始應力場變化。筆者分析對象為圓形洞室,其高跨比H/B=1<1.7,根據(jù)《規(guī)范》,地下洞室周邊IV級圍巖的側壓力系數(shù)的區(qū)間函數(shù)可確定為λ=[0.15,0.30]。
目前區(qū)間函數(shù)的求解方法較多,但大多數(shù)區(qū)間計算的解區(qū)間容易引起擴張,造成計算區(qū)間寬度大于實際寬度,必須尋求能控制解區(qū)間擴張的區(qū)間分析方法,筆者采用區(qū)間擴張法進行區(qū)間分析[7]。結合以上影響地下洞室圍巖位移的因素——變形模量、泊松比和應力場,通過函數(shù)構造處理來控制區(qū)間擴張,以滿足工程實際精度要求。
2.1 地下洞室圍巖位移區(qū)間函數(shù)的求解
根據(jù)區(qū)間函數(shù)的計算方法,則有
且存在常數(shù)K,使:
計算中可以通過控制N獲得值域的任意精度,即
根據(jù)區(qū)間分析法,在計算過程中誤差η可以按照相對精度要求ε來控制,即
2.2 地下洞室圍巖位移區(qū)間函數(shù)的計算
首先確定地下洞室圍巖位移的計算精度,設定區(qū)間的超寬度小于0.05,即η≥0.95,同時為了分析的典型性,筆者考慮拱頂處圍巖位移的計算,即θ=90°。
地下洞室圍巖位移的基本影響參變量為
式中:X1=E=[1.3,6],X2=μ=[0.30.0.35],X3=λ=[0.15,0.30]
F(E,μ,λ)=
令θ=90°,則有
F(E,μ,λ)=
1)當N=N1=1時,有
式中:ri為地下洞室半徑,mm;σz為地下洞室圍巖豎向應力,MPa。
2)當N=N2=2時,有
3)當N=N3=3時,有
地下洞室施工過程中圍巖位移的不確定性主要體現(xiàn)在施工過程中地質環(huán)境、工程活動和圍巖空間力學特性的不確定性兩個方面。在研究分析上,前者主要體現(xiàn)在圍巖位移計算模型的不確定性;后者則表現(xiàn)為圍巖計算參數(shù)的不確定性。筆者則針對影響圍巖位移的參數(shù)的不確定特點,將傳統(tǒng)圍巖位移計算模型與區(qū)間分析方法相結合,建立了適用于區(qū)間分析的地下洞室圍巖計算模型,利用該模型分析了IV級圍巖條件下拱頂位移的區(qū)間分布特點。結果表明,該模型客觀上更能反映圍巖位移的實際變形性態(tài),研究成果可為后續(xù)工程設計、施工以及支護時機的選擇提供參考。但同時要說明的是筆者對地下洞室圍巖位移不確定性的分析還是初步的,沒有考慮計算模型不確定和圍巖參數(shù)不確定同時存在的情況,因此,研究有待于進一步地深化。
[1] 豐明海,何振寧. 鐵路隧道施工中圍巖變形失穩(wěn)工程地質問題[J]. 工程地質學報,2014,22(4):677-683. FENG Minghai,HE Zhenning. Engineering geological problems of deformation and destabilization of surrounding rocks in railway tunnel construction[J].JournalofEngineeringGeology,2014,22(4):677-683.
[2] 趙新銘,劉寧,張劍. 巖土力學反分析的數(shù)值反演方法[J]. 水利水電科技進展,2003,23(2):55-58. ZHAO Xinming,LIU Ning,ZHANG Jian. Numerical inverse method of the bank analysis in geotechnical mechanics[J].WaterConservancyandHydropowerScienceandTechnologyProgress,2003,23(2):55-58.
[3] 王繼承,蘇永華. 何滿潮.區(qū)間分析方法在深部巖體工程中的應用[J].巖石力學與工程學報,2005,24(21):3835-3840. WANG Jicheng,SU Yonghua,HE Manchao. Application of interval analysis methods to surrounding rock stability of deep roadway engineering[J].ChineseJournalofRockMechanicsandEngineering,2005,24(21):3835-3840.
[4] 蔣沖,趙明華,曹文貴. 基于區(qū)間分析的巖土結構非概率可靠性分析[J]. 湖南大學學報(自然科學版),2008,35(3):1-4. JIANG Chong,ZHAO Minghua,CAO Wengui. Non-probabilistic reliabilityresearch on uncertain parameters in rock mechanics[J].JournalofHunanUniversity(NaturalScience),2008,35(3):1-4.
[5] 李興高,劉維寧.剛性擋墻上土壓力不確定性的計算研究[J].巖土工程學報,2007,29(3):353-359. LI Xinggao,LIU Weining. Study on the indeterminacy of the earth pressure on rigid retaining walls[J].ChineseJournalofGeotechnicalEngineering,2007,29(3):353-359.
[6] 王登剛,劉迎曦,李守巨,等. 巷道圍巖初始應力場和彈性模量的區(qū)間反演方法[J].巖石力學與工程學報,2002,21(3):305-308. WANG Denggang,LIU Yingxi,LI Shouju,et al. Interval back analysis on initial stresses and elastic modulus of surrounding rocks of tunnels[J].ChineseJournalofRockMechanicsandEngineering,2002,21(3):305-308.
[7] 張永杰,曹文貴,趙明華,等. 巖溶區(qū)公路路基穩(wěn)定性的區(qū)間模糊評判分析方法[J].巖土工程學報,2011,33(1):38-45. ZHANG Yongjie,CAO Wengui,ZHAO Minghua,et al. Interval fuzzy judgment method for roadbed stability in Karsts area[J].ChineseJournalofGeotechnicalEngineering,2011,33(1):38-45.
[8] 王德人,張連生,鄧乃揚.非線性方程的區(qū)間算法[M].上海:上海科學技術出版社,1987. WANG Deren,ZHANG Liansheng,DENG Naiyang.IntervalAlgorithmofNonlinearEquations[M].Shanghai: Shanghai Science and Technology Press,1987.
[9] 黃明奎.鋼管混凝土構件力學性質的區(qū)間分析方法[J].四川建筑科學研究,2009,35(2):32-35. HUANG Mingkui. The interval analysis method of the mechanical property for concrete-filled steel tube[J].SichuanBuildingScience,2009,35(2):32-35.
[10] 夏永旭,王永東.隧道結構力學計算[M].北京:人民交通出版社,2004. XIA Yongxu,WANG Yongdong.MechanicalCalculationofTunnelStructure[M].Beijing: China Communications Press,2004.
Analysis on Displacement of Surrounding Rock in Tunnel with Consideration of Interval Influencing Factors
HUANG Mingkui,QUAN Jian
( School of Civil Engineering,Chongqing Jiaotong University,Chongqing 400074,P. R. China)
Surrounding rock displacement was one of the important parameters to assess rock stability in tunnel and the uncertainty of these parameters was seldom considered in the conventional calculation of surrounding rock displacement. Because of the uncertainty of parameters that influence the surrounding rock displacement, a new displacement calculating model suitable for the interval analysis was built by combining the conventional displacement calculating model with the interval analysis method, and then the model was applied to analyze the interval distribution features of arch crown displacement of surrounding rock condition of rating IV. The analysis result shows that the model can objectively reflect the actual deformation behavior of rock due to displacement and the research results can be taken as the reference for determination on subsequent works design, execution and the time for supporting .
tunnel engineering;surrounding rock;displacement;interval analysis
10.3969/j.issn.1674-0696.2016.02.06
2014-10-01;
2015-04-02
國家自然科學基金項目(51308574);重慶市教委基金項目(KJI20407)
黃明奎(1975—),男,四川隆昌人,教授,博士(后),主要從事隧道及巖土工程方面的研究。E-mial:hmksile@163.com。
U453;O319.56
A
1674-0696(2016)02-021-03