孟建宇，秦?fù)碥?，謝良甫，于廣明

(1.新疆大學(xué)建筑工程學(xué)院，新疆 烏魯木齊 830047；2.青島理工大學(xué)土木工程學(xué)院，山東 青島 266033；3.青島理工大學(xué)藍(lán)色經(jīng)濟(jì)區(qū)工程建設(shè)與安全協(xié)同創(chuàng)新中心，山東 青島 266033)

0 引 言

隧道工程建設(shè)是通過對地質(zhì)體的改造利用來實(shí)現(xiàn)交通運(yùn)輸功能的地質(zhì)工程活動。隧道結(jié)構(gòu)位于地下獨(dú)特的建筑環(huán)境中，對土體性質(zhì)的評估結(jié)果將直接影響工程的設(shè)計(jì)和建設(shè)，甚至對工程的質(zhì)量安全造成影響[1-4]。在設(shè)計(jì)和施工過程中，通常用給出的土體性質(zhì)參數(shù)均值或設(shè)計(jì)值來代表土性[5-7]。但在實(shí)際工程中，土體是非均質(zhì)材料，空間中每一點(diǎn)的性質(zhì)參數(shù)都不相同，用一個(gè)確定值來概括土體性質(zhì)顯然是不準(zhǔn)確的[4,8,9]。在描述土性在空間變異性方面的理論中[10-13]，Vanmarcke[14]最早提出了土性剖面的隨機(jī)場模型，可以將土層剖面上巖土參數(shù)的不確定性很好地再現(xiàn)出來。國內(nèi)學(xué)者[15-20]對土性參數(shù)不確定性的關(guān)注度也越來越高，對土體參數(shù)不確定性的研究方法進(jìn)行了探討，并且將這些方法應(yīng)用在隧道、基坑和邊坡等工程問題中，得到了一些基于概率統(tǒng)計(jì)的可靠度分析結(jié)論。程紅戰(zhàn)等[15]運(yùn)用隨機(jī)場理論對盾構(gòu)隧道的地表變形進(jìn)行了研究。薛亞東等[16]對邊坡在各向異性隨機(jī)場下的可靠度進(jìn)行了分析。張繼周等[18]對土性參數(shù)的不確定性描述方法進(jìn)行了探討。

本文運(yùn)用隨機(jī)場理論，考慮土體彈性模量的不確定性，對盾構(gòu)隧道開挖引起的隧道變形規(guī)律進(jìn)行研究，并與確定性分析結(jié)果進(jìn)行對比，分析土體彈性模量空間變異性的影響。

1 工程背景

烏魯木齊地鐵1號線修建過程中在多個(gè)區(qū)間采用盾構(gòu)法開挖。盾構(gòu)法具有全自動、掘進(jìn)速度快、不影響地面交通的優(yōu)點(diǎn)。施工過程主要有土體開挖、開挖面支護(hù)、盾構(gòu)推進(jìn)、襯砌拼裝、盾尾脫空。盾構(gòu)機(jī)掘進(jìn)后，地層應(yīng)力得到一定釋放，此時(shí)機(jī)殼可以充當(dāng)初始支護(hù)。在盾尾脫空后，將管片拼裝承受地層作用。本文通過控制應(yīng)力釋放系數(shù)的方法，模擬開挖卸荷后土體與盾構(gòu)機(jī)殼接觸的過程，以及管片拼裝后承受地層作用的效果。盾構(gòu)法施工步驟分析見圖1。

圖1 盾構(gòu)法施工步驟分析

采用FLAC3D有限差分軟件對烏魯木齊1號線某盾構(gòu)區(qū)間進(jìn)行建模。根據(jù)工程地質(zhì)情況，土體材料選擇M-C本構(gòu)，管片及盾構(gòu)機(jī)殼按線彈性考慮。由于此次是對土體彈性模量的空間變異性進(jìn)行研究，將多個(gè)地層細(xì)分計(jì)算意義不大，故按照單一地層進(jìn)行簡化計(jì)算。在收集整理烏魯木齊1號線資料的基礎(chǔ)上，取隧道直徑為10 m，軸線埋深18 m，模型尺寸為50 m×36 m(長×寬)，軸線方向厚度為1 m。模型材料屬性見表1。表1中，null是FLAC3D中用來模擬隧道開挖的空值材料，土體彈性模量為均值，實(shí)際計(jì)算中將參數(shù)隨機(jī)場映射在模型中。

表1 模型材料屬性

2 基于隨機(jī)場理論的盾構(gòu)隧道變形分析方法

2.1 隨機(jī)場理論分析

巖土體性質(zhì)參數(shù)在空間上具有隨機(jī)性和復(fù)雜性的特點(diǎn)，目前基于土體是均一、各向同性材料的假設(shè)，取土體參數(shù)的均值或設(shè)計(jì)值代表土性的方法難以考慮到土體的非均質(zhì)特點(diǎn)。運(yùn)用隨機(jī)場理論生成的土層剖面隨機(jī)場模型，將巖土體參數(shù)按照一定的概率分布離散在隨機(jī)場模型上，用均值、方差及協(xié)方差等來表征其統(tǒng)計(jì)特征，可以較好地概括巖土體參數(shù)的隨機(jī)性。此外，在對土性參數(shù)的試驗(yàn)結(jié)果進(jìn)行統(tǒng)計(jì)分析中可以發(fā)現(xiàn)[8,21-23]，土體性質(zhì)在具有空間上的變異性、不確定性的同時(shí)，也存在空間上的連續(xù)性。傳統(tǒng)意義上將土體參數(shù)看做隨機(jī)變量的方法，沒有考慮到土性參數(shù)在空間位置上的這種特點(diǎn)，而通過建立隨機(jī)場模型可以彌補(bǔ)這一點(diǎn)。

隨機(jī)場模型中利用變異系數(shù)、相關(guān)函數(shù)、相關(guān)距離等特征來反映土體性質(zhì)在空間上的特點(diǎn)。在對土性相關(guān)距離進(jìn)行研究后，目前普遍認(rèn)為相關(guān)距離是具有實(shí)際意義的一項(xiàng)指標(biāo)，用來描述巖土性質(zhì)在空間連續(xù)性上的屬性，在相關(guān)距離內(nèi)巖土體性質(zhì)相近或變化較慢，在相關(guān)距離以外土性變化較大較快。相關(guān)函數(shù)用來描述土體中任意2點(diǎn)之間相關(guān)性，本文采用Markov相關(guān)函數(shù)形式。變異系數(shù)是用來描述分布形狀的一種形式，主要用于數(shù)據(jù)值都是正值的情況。公式如下

(1)

(2)

式中，d是模型中任意2點(diǎn)的距離；λ是相關(guān)距離；ρ(d)是關(guān)于2點(diǎn)距離d的自相關(guān)函數(shù)，本文給出的是二維各向同性隨機(jī)場的相關(guān)函數(shù)；Cv是變異系數(shù)；σ是數(shù)據(jù)樣本的標(biāo)準(zhǔn)方差；μ為樣本均值。

2.2 隨機(jī)場模型建立

由工程經(jīng)驗(yàn)和現(xiàn)有研究可知，土體參數(shù)中彈性模量對隧道變形影響最大。本文考慮土體彈性模量空間變異性對隧道變形的影響，利用協(xié)方差矩陣分解法將MATLAB隨機(jī)場建模與FLAC3D數(shù)值建模，通過fish語言一對一單元映射相結(jié)合，生成二維有限差分各向同性隨機(jī)場。該方法首先需要構(gòu)建協(xié)方差矩陣C，C為隨機(jī)場中n個(gè)離散點(diǎn)之間協(xié)方差所組成的n階矩陣，協(xié)方差通過自相關(guān)函數(shù)計(jì)算，該矩陣正定對稱，進(jìn)行Cholesky分解得到

C=XU=XXT

(3)

式中，X、U分別為下三角矩陣和上三角矩陣，XT為矩陣的轉(zhuǎn)置。生成由n個(gè)相互獨(dú)立且服從標(biāo)準(zhǔn)正態(tài)分布隨機(jī)數(shù)組成的列向量Y，進(jìn)行二階平穩(wěn)標(biāo)準(zhǔn)正態(tài)隨機(jī)場Z的一次實(shí)現(xiàn)，即

Z=XY

(4)

對矩陣Z進(jìn)行相應(yīng)數(shù)學(xué)變換，可以得到具有不同統(tǒng)計(jì)特征(即不同分布形態(tài))的正態(tài)平穩(wěn)隨機(jī)場，多次生成不同的Y向量就可以進(jìn)行隨機(jī)場的多次刷新。考慮土體參數(shù)的非負(fù)特性，此次生成對數(shù)正態(tài)分布隨機(jī)場。最后，結(jié)合Monte-Carlo方法開展隨機(jī)性分析，選取500次作為每組分析的總頻率。基于隨機(jī)場理論的盾構(gòu)隧道變形分析流程見圖2。

圖2 隨機(jī)場下盾構(gòu)隧道變形分析過程

3 確定性分析

根據(jù)建立的數(shù)值模型，首先進(jìn)行不考慮彈性模量空間變異性的確定性分析。在FLAC3D中建立目標(biāo)區(qū)段隧道模型，通過控制應(yīng)力釋放的方式模擬開挖和支護(hù)過程中地層承載能力的效果，以及地應(yīng)力對隧道的階段性影響。本文在模擬中，取開挖后的應(yīng)力釋放系數(shù)為0.2。獲得模型的初始應(yīng)力場后，對隧道進(jìn)行開挖，此時(shí)應(yīng)力釋放率為0.2，之后施加管片支護(hù)，地應(yīng)力全部釋放。計(jì)算結(jié)果見圖3。從圖3可知，隧道變形最顯著的地方在拱頂、拱底和隧道兩側(cè)的中間位置，具體位置對應(yīng)隧道管片襯砌外沿與圍巖的交界面處以及管片襯砌內(nèi)表面上，在相應(yīng)位置布設(shè)測點(diǎn)提取數(shù)據(jù)。確定性分析數(shù)據(jù)見表2。

圖3 位移云圖

表2 隧道變形確定性分析結(jié)果 mm

4 隨機(jī)性分析

按隨機(jī)場的生成方法，建立土體彈性模量的二階平穩(wěn)正態(tài)隨機(jī)場，借助fish語言將隨機(jī)場模型映射到數(shù)值模型中。此時(shí)，模型中各單元的彈性模量都不相同，其他參數(shù)均為確定值。由于FLAC3D中Monte-Carlo本構(gòu)輸入?yún)?shù)是體積模量和剪切模量，需要用彈性模量進(jìn)行轉(zhuǎn)換，此處不作具體說明。圖4是土體彈性模量隨機(jī)場在數(shù)值模型中1次實(shí)現(xiàn)過程中對應(yīng)體積模量的表現(xiàn)，即隨機(jī)場模型示意圖，灰度越深代表該處體積模量越大。

圖4 隨機(jī)場模型示意

多次重復(fù)更新隨機(jī)場并計(jì)算得到結(jié)果，利用Monte-Carlo方法開展隨機(jī)性分析，研究各向同性隨機(jī)場中，彈性模量相關(guān)距離d和變異系數(shù)Cv對隧道變形的影響，隧道變形測點(diǎn)參照確定性分析結(jié)果布設(shè)。此次研究彈性模量均值取50 MPa，按照變異系數(shù)Cv取0.3，相關(guān)距離d分別取1、3、5、10 m和100 m時(shí)，考慮相關(guān)距離d改變對隧道變形的影響。按照相關(guān)距離d取10 m，變異系數(shù)Cv分別取0.1、0.3、0.5時(shí)，考慮變異系數(shù)Cv改變對隧道變形的影響，共計(jì)8組模擬工況，考慮到計(jì)算精度和效率，每組工況計(jì)算500次。圖5為其中1組工況下某個(gè)測點(diǎn)的原始數(shù)據(jù)。

圖5 原始數(shù)據(jù)

4.1 模擬次數(shù)的影響

利用Monte-Carlo方法開展隨機(jī)分析，需要產(chǎn)生具有一定數(shù)量隨機(jī)值的樣本，從中抽樣得到模型的統(tǒng)計(jì)特征。理論上，Monte-Carlo方法需要進(jìn)行大量的試驗(yàn)?zāi)M，模擬次數(shù)越多，所得到結(jié)果越精確。以樣本均值為例，選取其中1組工況，在MATLAB中對1～500次模擬過程中隧道變形測點(diǎn)值進(jìn)行統(tǒng)計(jì)分析，見圖6。從圖6可知，開始時(shí)，模擬次數(shù)較少，均值波動范圍較大，此時(shí)樣本數(shù)字特征誤差很大，難以得到具有實(shí)際意義的結(jié)論。模擬次數(shù)達(dá)到300次時(shí)，隧道側(cè)向變形量、拱頂下沉量、拱底隆起量的數(shù)學(xué)期望開始穩(wěn)定，最終樣本均值在一個(gè)確定值附近小范圍波動，圖6中直線即為該確定值位置，且誤差遠(yuǎn)小于0.05 mm。選取500次計(jì)算可以滿足計(jì)算精度和計(jì)算效率等方面的要求。

圖6 隧道變形平均值隨模擬次數(shù)變化趨勢

4.2 變異系數(shù)的影響

相關(guān)距離d取10 m時(shí)，變異系數(shù)Cv分別取0.1、0.3和0.5，提取測點(diǎn)數(shù)據(jù)進(jìn)行概率統(tǒng)計(jì)分析，對應(yīng)的概率密度分布見圖7。從圖7可知，隧道兩側(cè)位移、拱頂位移、拱底位移的概率密度分布隨著變異系數(shù)的增大越來越分散，同時(shí)分布曲線在向右側(cè)偏移，即概率統(tǒng)計(jì)均值在逐漸變大。變異系數(shù)越大，模型中土體彈性模量變化程度越大，不同單元的剛度差異越明顯，對應(yīng)統(tǒng)計(jì)分析結(jié)果可以看到變形結(jié)果離散程度越來越大，這也說明變異系數(shù)越大，對隧道變形的影響越明顯。隨機(jī)場下計(jì)算得到的隧道變形結(jié)果均值隨著變異系數(shù)的增大而增大，土體彈性模量的空間變異程度對隧道變形有重要影響。

圖7 不同變異系數(shù)下隧道變形值概率密度分布

4.3 相關(guān)距離的影響

在各向同性隨機(jī)場下分析相關(guān)距離對隧道變形的影響。根據(jù)確定性分析結(jié)果，在隧道變形顯著位置布設(shè)測點(diǎn)獲得變形數(shù)據(jù)，即隧道左側(cè)位移值、右側(cè)位移值、拱頂位移值和拱底位移值。對數(shù)據(jù)結(jié)果進(jìn)行統(tǒng)計(jì)分析可以發(fā)現(xiàn)，各組數(shù)據(jù)具有顯著的正態(tài)分布特征，概率密度曲線兩頭低，中間高。圖8為其中1組工況下隧道變形的概率分布特征。對于服從正態(tài)分布的樣本數(shù)據(jù)，樣本均值(期望值)決定分布曲線位置，標(biāo)準(zhǔn)方差決定分布的幅度，即在統(tǒng)計(jì)學(xué)意義上，曲線位置反映隨機(jī)性分析中隧道變形均值大小，曲線分布幅度決定隨機(jī)性分析中變形結(jié)果的離散情況。

圖8 隨機(jī)場下隧道變形值概率密度分布

變異系數(shù)取0.3時(shí)，相關(guān)距離分別取1、3、5、10、50 m和100 m，對應(yīng)分布的特征參數(shù)(均值和標(biāo)準(zhǔn)方差)變化規(guī)律見圖9、10。圖9中，直線是確定性分析中隧道變形值，折線代表隧道變形均值隨相關(guān)距離的變化趨勢。從圖9可知，在隨機(jī)場下，隧道變形量均大于確定性結(jié)果，這也證明了考慮土性參數(shù)時(shí)引入相關(guān)距離的必要性，在傳統(tǒng)分析中，將土性參數(shù)看作確定值往往與實(shí)際情況有所偏差。從圖10可知，隨著相關(guān)距離的增大，隧道變形結(jié)果的標(biāo)準(zhǔn)方差在不斷增大，這表明計(jì)算結(jié)果在趨于離散，變形結(jié)果的概率分布趨于分散。

圖9 隧道變形均值隨相關(guān)距離變化規(guī)律

圖10 隧道變形均方差隨相關(guān)距離變化規(guī)律

對統(tǒng)計(jì)特征進(jìn)行分析，隧道發(fā)生顯著變形時(shí)與圍巖低剛度區(qū)域分布有直接聯(lián)系。將隨機(jī)場參數(shù)映射在數(shù)值模型中時(shí)，模型每1點(diǎn)的土體彈性模量都不完全相同，也就代表模型中某些區(qū)域的剛度將小于平均水平，即在模型中出現(xiàn)了低剛度區(qū)域，這導(dǎo)致隧道發(fā)生顯著變形處的位移值相對確定性分析結(jié)果偏大。相關(guān)距離相對較小時(shí)，低剛度區(qū)域集中范圍較小，對隧道變形的影響在一定程度上并不顯著，但當(dāng)相關(guān)距離超出一個(gè)范圍并不斷增大時(shí)，低剛度區(qū)域范圍也會不斷增大，這樣直接導(dǎo)致隧道變形會隨著相關(guān)距離的增大而增大。相關(guān)距離越大，單次隨機(jī)場模擬計(jì)算中單元差異越小，而多次隨機(jī)場模擬中單元差異會越大，即隨著相關(guān)距離的增大，計(jì)算結(jié)果的統(tǒng)計(jì)特征表現(xiàn)得越來越離散。

5 結(jié) 語

本文將隨機(jī)場計(jì)算機(jī)模擬方法用于盾構(gòu)隧道變形分析，主要考慮土體彈性模量空間變異性對隧道變形的影響，得出以下結(jié)論：

(1)利用Monte-Carlo隨機(jī)模擬方法，結(jié)合隨機(jī)場理論和數(shù)值分析的計(jì)算機(jī)模擬，計(jì)算次數(shù)足夠多時(shí)，可以獲得很好的隨機(jī)分析結(jié)果，進(jìn)行土體彈性模量空間變異性分析是可行的。

(2)隨著土體彈性模量變異系數(shù)的增大，隧道變形結(jié)果越來越離散，同時(shí)變形結(jié)果均值越來越大，說明彈性模量變異程度對隧道變形有重要影響。

(3)考慮土性參數(shù)引入相關(guān)距離是必要的，不同相關(guān)距離下隧道變形結(jié)果的均值大于確定性分析結(jié)果，這與低剛度區(qū)域的分布有關(guān)，同時(shí)相關(guān)距離越大，計(jì)算結(jié)果分布越分散。