林　超

(貴州高速公路集團(tuán)有限公司)

?

隧道工程偏壓研究與反分析

林超

(貴州高速公路集團(tuán)有限公司)

摘要：針對(duì)某隧道工程研究了隧道偏壓的成因與判別方法，采用反分析方法確定了隧道偏壓的分布特征，對(duì)預(yù)測(cè)隧道狀態(tài)和反饋設(shè)計(jì)具有重要意義。

關(guān)鍵詞：偏壓隧道；荷載反分析；荷載結(jié)構(gòu)模型；新奧法

1工程概況

羅(源)長(zhǎng)(樂)高速公路馬宅頂隧道屬丘陵地貌，標(biāo)高為+50～110 m，地形呈波狀起伏，進(jìn)口處坡角約為20°，出口處坡角約為37°。隧道進(jìn)出口處有3、8 m殘坡積層出露，其余地段為花崗巖及其風(fēng)化層直接出露，圍巖類別為II類；洞身段以弱風(fēng)化混合花崗巖為主，少部分凝灰?guī)r，局部輝長(zhǎng)(綠)巖脈巖侵入，巖圍巖類別為III、Ⅳ類，其巖性均為硬質(zhì)巖，質(zhì)地較完整，致密而堅(jiān)硬。

隧道結(jié)構(gòu)按新奧法原理進(jìn)行設(shè)計(jì)，采用復(fù)合襯砌，型式為不對(duì)稱雙連拱隧道；由于隧道所處地形起伏，進(jìn)出口段埋深較淺，圍巖類別較低，可能產(chǎn)生偏壓。

2偏壓隧道的判定及荷載計(jì)算

2.1　偏壓隧道的成因

(1)施工原因：因施工方法不當(dāng)引起開挖斷面局部坍塌，從而改變了圍巖壓力的相對(duì)穩(wěn)定性，造成應(yīng)力集中而引起隧道偏壓，如處理得當(dāng)，一般不會(huì)影響正常施工。

(2)地質(zhì)原因：因巖層產(chǎn)狀傾斜，節(jié)理發(fā)育，其間又有軟弱結(jié)構(gòu)面或滑動(dòng)面，自穩(wěn)能力極差，施工中一旦受到擾動(dòng)，巖體就會(huì)沿層理面出現(xiàn)滑動(dòng)。

(3)地形原因：隧道傍山，地面顯著傾斜，側(cè)壓力較大，且隧道理深較淺。

2.2　偏壓隧道的判斷

(1)地形引起的偏壓：圍巖類別、地面坡度和覆蓋層厚度是判別隧道偏壓的3個(gè)重要因素。當(dāng)隧道外側(cè)拱肩至地表面的垂直距離t值等于或小于下表所列數(shù)值時(shí)，應(yīng)視為偏壓隧道。一般在III類以下圍巖中，以地形引起的偏壓為主。

表1　拱肩至地表面垂直距離t值　m

(2)地質(zhì)構(gòu)造引起的偏壓：地質(zhì)構(gòu)造常在多裂隙圍巖(以IV-III類較為突出)中引起隧道偏壓，其壓力分布主要與下列因素有關(guān)。

①圍巖的工程地質(zhì)條件及控制性裂隙、節(jié)理或?qū)永?統(tǒng)稱為弱面)的產(chǎn)狀及其與隧道軸線的組合關(guān)系；

②圍巖擾動(dòng)范圍；

③控制性弱面的強(qiáng)度以及作用在弱面上的法向力大小等。隧道一側(cè)受2個(gè)傾斜的軟弱面(傾角為α)及一組節(jié)理面所切割時(shí)，會(huì)形成不穩(wěn)定塊體，當(dāng)圍巖的內(nèi)摩擦角φ小于弱面傾角α?xí)r，巖層將沿弱面滑動(dòng)并產(chǎn)生偏壓。

(3)施工原因引起的偏壓：由于開挖不當(dāng)或支護(hù)不及時(shí)引起一側(cè)圍巖發(fā)生局部坍塌，或回填不實(shí)造成不穩(wěn)定土體，人為形成了偏壓的地質(zhì)構(gòu)造。

2.3　施工階段偏壓荷載的確定

隧道偏壓的產(chǎn)生，不僅僅是由地形和埋深、圍巖類別決定的，洞型和施工也會(huì)對(duì)它產(chǎn)生影響。要對(duì)這些因素進(jìn)行準(zhǔn)確的評(píng)估是很困難的，最佳的途徑應(yīng)是根據(jù)隧道在施工過程中的各種表現(xiàn)(隧道系統(tǒng)的輸出)，反分析作用在隧道上的各種因素，主要包括荷載和地層參數(shù)(隧道系統(tǒng)的輸入)。

根據(jù)模型試驗(yàn)結(jié)果，可假設(shè)作用于襯砌結(jié)構(gòu)上的壓力荷載呈拋物線分布模式，如圖1所示。圖中，將襯砌左邊緣與上邊緣的交點(diǎn)(左上角點(diǎn))定義為x-y坐標(biāo)系的原點(diǎn)O，且假設(shè)pi，qi(i=1，4) 為待求未知量。則垂直和水平分布力p(x，y)，q(x，y)可寫成以下形式:

圖1　襯砌壓力荷載分布模式

q(x，y)=bo+b1y+b2y2+ b3x

(1)

式中：ai，bi(i=0，3)是與pi，qi相關(guān)的。用pi，qi，代換ai，bi，上式可變換為

p(x,y)=p1+(2dx-1)dxp21+dxp31+4(1-dx)dxp41

q(x,y)=q1+(2dy-1)dyq21+dyq31+4(1-dy)dyq41

(2)

式中：dx=x/XL，dy=y/YL，Pi1=pi-p1，qi1=qi-q1(i=2,4)，XL，YL分別是結(jié)構(gòu)外緣在x，y向的最大尺寸。

這種荷載模式共有8個(gè)參數(shù)，可以模擬各種情況，比如，均勻分布，線性分布，拋物線分布等。在對(duì)反分析初值沒有把握的情況下，可以先進(jìn)行均勻分布荷載的反演，然后，用反演結(jié)果作為拋物線分布荷載反演的初值。

3反分析方法

這里采用直接反分析方法，采用優(yōu)化方法通過不斷調(diào)整需反演參數(shù)的取值使得正算值與實(shí)測(cè)值的差異最小，認(rèn)為這時(shí)的參數(shù)值就是我們所需要的反演值。這種辦法由于反分析過程和正算過程的獨(dú)立性，便于編輯通用程序，適用的范圍更廣。

反分析的目標(biāo)函數(shù)一般選擇如下形式，

(3)

式中：fi為監(jiān)測(cè)量的有限元計(jì)算，ui為實(shí)際量測(cè)值，n為測(cè)點(diǎn)總數(shù)，J與fi為待反分析參數(shù)X的函數(shù)，通過不段優(yōu)化使得目標(biāo)函數(shù)取得最小值。

在本工程中，由于檢測(cè)項(xiàng)目較多，既有量值極小的位移，又有量值很大的襯砌內(nèi)力值和圍巖壓力值，直接采用上述的目標(biāo)函數(shù)，會(huì)出現(xiàn)小量被大量吞噬的現(xiàn)象，即目標(biāo)函數(shù)值由大量控制，而小量不起作用。因此，我們對(duì)目標(biāo)函數(shù)進(jìn)行了改造，采取如下的形式

(4)

fij、uij含義同上，其下標(biāo)的意義為第i個(gè)監(jiān)測(cè)項(xiàng)目的第j個(gè)測(cè)點(diǎn)。

反分析的優(yōu)化方法目前有很多，梯度法總是沿著函數(shù)值下降的方向?qū)ふ易顑?yōu)點(diǎn)，因此開始優(yōu)化時(shí)效率較高，但隨著優(yōu)化的進(jìn)行效率開始降低，同時(shí)它找到的是局部最優(yōu)點(diǎn)，而且對(duì)目標(biāo)函數(shù)的要求很高(可導(dǎo)性、可顯式表示等)；另外，還有很多別的優(yōu)化方法。這里采用單純形法。

正分析和反分析都由程序?qū)崿F(xiàn)，反分析的目標(biāo)參量可以是圍巖壓力，也可以是地層彈簧剛度，或者是二者一起參加反演；正分析可以是普通的山嶺隧道，也可以是盾構(gòu)法拼裝隧道。

4輸入?yún)?shù)的確定

4.1　中墻內(nèi)力計(jì)算

量測(cè)得到的是中墻鋼筋的應(yīng)力或中墻的表面應(yīng)變，而有限元模型中把中墻簡(jiǎn)化為梁，因此，要把應(yīng)力或應(yīng)變轉(zhuǎn)化為軸力、彎矩，才能作為輸入?yún)?shù)。

表2　中墻實(shí)測(cè)內(nèi)力

4.2　襯砌壓力的確定

襯砌壓力是時(shí)間的函數(shù)，隨著襯砌強(qiáng)度的增加，圍巖變形的穩(wěn)定，壓力逐漸達(dá)到穩(wěn)定。反分析時(shí)就采用達(dá)到穩(wěn)定時(shí)的壓力值，這里采用右洞K57+030截面的壓力值，如圖2所示。

圖2　二襯法向壓力圖

5反分析結(jié)果

首先計(jì)算圍巖壓力的理論值，采用規(guī)范提供的公式

q=0.45×26-sγω

ω=1+i(B-5)=1+0.1(28-5)=3.3

q=0.45×6-3×23×3.3=273

e=0.2×273=55

我們將以此為初值進(jìn)行后面的反分析，同時(shí)，地層抗力值的取值也按規(guī)范取值，Kn=5.e5，Ks=1.e5。

首先假定地層抗力取上述值，進(jìn)行荷載反分析，結(jié)果見表3，然后將荷載反演值作為定值進(jìn)行地層參數(shù)反分析，結(jié)果見表4，最后利用反演的地層參數(shù)值再進(jìn)行荷載反分析，其結(jié)果見表5，其對(duì)應(yīng)的預(yù)測(cè)值與實(shí)測(cè)值的對(duì)比如表6所示。

表3　荷載反分析結(jié)果(步長(zhǎng)為1，收斂值0.000 1)

表4　地層參數(shù)反分析結(jié)果

表5　第二次荷載反分析結(jié)果

表6　實(shí)測(cè)值、預(yù)測(cè)值對(duì)比表

由上表可見，預(yù)測(cè)值與實(shí)測(cè)值總體上較為符合，符合較差的是軸力項(xiàng)目和E點(diǎn)的圍巖壓力。我們看到，中墻在模型中只在兩端受力，由隔離體的平衡條件可知，中墻上部和下部的軸力應(yīng)該是相等的，由預(yù)測(cè)值也可以發(fā)現(xiàn)這一點(diǎn)；實(shí)測(cè)值上下部不一樣，可能是有別的原因，比如施工時(shí)將支撐撐在中墻上。而D、E點(diǎn)的圍巖壓力值發(fā)生了一定的突變，顯然預(yù)測(cè)值比實(shí)測(cè)值更為平滑，更為合理。

反演值與前面按規(guī)范的計(jì)算值較為接近，由于圍巖較好，可以看到側(cè)壓力與垂直壓力相比較小；同時(shí)，隧道存在一定的偏壓，大洞所受壓力稍大，不過不是很明顯。

6結(jié)語

通過馬宅頂隧道工程的實(shí)踐，我們進(jìn)行了較為全面、細(xì)致地監(jiān)控量測(cè)，取得了第一手的寶貴資料；在對(duì)量測(cè)資料進(jìn)行了認(rèn)真的整理后，以量測(cè)資料為基礎(chǔ)，進(jìn)行了荷載與地層參數(shù)的反分析；然后以反分析結(jié)果為基礎(chǔ)，進(jìn)行襯砌受力狀態(tài)的預(yù)報(bào)，可以對(duì)反饋設(shè)計(jì)提供指導(dǎo)。最后，我們得到如下結(jié)論。

(1)以監(jiān)測(cè)、反分析為基礎(chǔ)的信息化設(shè)計(jì)方法應(yīng)得到重視。監(jiān)控量測(cè)有雙重意義，它既為反分析提供數(shù)據(jù)，它本身又是隧道及襯砌穩(wěn)定狀態(tài)的判據(jù)。反分析結(jié)果可以直接作為反饋設(shè)計(jì)的輸入?yún)?shù)，其預(yù)測(cè)結(jié)果也可用于判斷隧道的安全狀態(tài)。

(2)反分析結(jié)果與規(guī)范計(jì)算值比較，豎向荷載最大值P=314kN/m2，比規(guī)范計(jì)算值273大15%;最小值P=184，比規(guī)范計(jì)算值小32%。水平向荷載同樣與規(guī)范計(jì)算值差異不很大，但卻反映了隧道偏壓的特征。可見，反分析方法比常規(guī)方法的解能更接近于參數(shù)的真實(shí)狀態(tài)，因而能作出更為準(zhǔn)確的預(yù)測(cè)。

參考文獻(xiàn)：

[1]鐘新樵．土質(zhì)偏壓隧道襯砌模型試驗(yàn)分析[J]．西南交通大學(xué)學(xué)報(bào)，1996，31(6):602-606.

[2]張海亮，董志．黃土質(zhì)偏壓隧道施工技術(shù)[J].鐵道建筑技術(shù)，2000，(4)：25-27.

[3]朱合華，崔茂玉，楊勁松，盾構(gòu)襯砌管片的設(shè)計(jì)模型與荷載分布的研究[J].巖土工程學(xué)報(bào)，2000，22(2):190-194.

[4]夏才初，潘國(guó)榮等，土木工程監(jiān)測(cè)技術(shù)[M].北京：中國(guó)建筑工業(yè)出版社，2001.

作者簡(jiǎn)介：林超(1979-)，男，工程師，研究方向：工程管理。

收稿日期：2014-11-14

中圖分類號(hào)：U442

文獻(xiàn)標(biāo)識(shí)碼：C

文章編號(hào)：1008-3383(2015)08-0115-02