隧道襯砌可靠指標的優(yōu)化求解方法

龔彥峰， 陳文博， 房 倩

(1. 中鐵第四勘察設(shè)計院集團有限公司， 湖北 武漢 430063； 2. 水下隧道技術(shù)湖北省工程實驗室， 湖北 武漢 430063；3. 北京交通大學城市地下工程教育部重點實驗室， 北京 100044)

0 引言

JC法是工程結(jié)構(gòu)可靠度計算中最常用的方法，由Rackwitz和Fiessleer于20世紀70年代末提出，是以當量正態(tài)化考慮隨機變量實際分布的可靠指標計算方法，之后被國際安全度聯(lián)合委員會(JCSS)推薦采用[1]。吳世偉于1990年在《結(jié)構(gòu)可靠度分析》[2]中對JC法的簡化和改進算法做了詳細說明； 1999年，徐軍等[3]采用改進的JC法對具有高次非線性和復雜性的三維巖體強度準則進行可靠指標計算； 李典慶等[4]于2002年提出一種改進的JC法，解決了功能函數(shù)為線性時JC法不一定收斂的問題； 李繼祥等[5]在李典慶等的基礎(chǔ)上進一步改進了JC法的求解過程，使得計算精度和收斂速度得到提高。

JC法適用于隨機變量為任意分布下結(jié)構(gòu)可靠指標的求解，計算精度基本能夠滿足工程要求。我國GB 50068—2001《建筑結(jié)構(gòu)可靠度設(shè)計統(tǒng)一標準》和Q/CR 9007—2014《鐵路工程結(jié)構(gòu)可靠度設(shè)計統(tǒng)一標準》均推薦JC法為工程結(jié)構(gòu)可靠指標計算方法[1]。2015年，中國鐵路總公司發(fā)布的Q/CR 9129—2015《鐵路隧道極限狀態(tài)法設(shè)計暫行規(guī)范》[6](簡稱《暫行規(guī)范》)也采用JC法對隧道結(jié)構(gòu)可靠指標進行校準，并優(yōu)化設(shè)計了表達式分項系數(shù)。以往JC法的應(yīng)用和研究均針對單一且連續(xù)光滑的功能函數(shù)，然而，對于隧道素混凝土襯砌或鋼筋混凝土襯砌，其極限狀態(tài)的功能函數(shù)在空間上是分段、不光滑的，分別由抗裂與抗壓極限狀態(tài)、大偏心與小偏心極限狀態(tài)的功能函數(shù)組合而成，不同襯砌的極限狀態(tài)功能函數(shù)在空間上雖然連續(xù)但不光滑，且各段有其可行條件。JC法在迭代過程中不能對驗算點位置有所限制，最終計算結(jié)果可能超出函數(shù)的可行域，因此，選取各段計算的可靠指標中的最小值作為結(jié)構(gòu)可靠指標，也是不合理的。為解決該問題，喻渝等[7]以偏心距均值和受壓區(qū)高度均值分別作為素混凝土和鋼筋混凝土襯砌的極限狀態(tài)功能函數(shù)選用的判據(jù)，但由于圍巖物理力學參數(shù)、襯砌材料性能等的變異性，使得襯砌的受力狀態(tài)具有不確定性，以偏心距或受壓區(qū)統(tǒng)計量作為判據(jù)缺乏理論支撐，其嚴謹性和準確性存疑。

為了能夠考慮功能函數(shù)分段性質(zhì)及各段的約束條件，本文引入基于可靠指標幾何意義的優(yōu)化求解方法。該方法與JC法同源于二階矩方法，其優(yōu)點是能夠通過約束函數(shù)限制解的可行域。近年來，優(yōu)化算法的研究和應(yīng)用已經(jīng)十分成熟， 1991年，Liu等[8]采用梯度投影法、罰函數(shù)法、增廣拉格朗日法、序列二次規(guī)劃法等優(yōu)化方法進行可靠指標的計算； 之后Powell法、遺傳算法、粒子群算法、蟻群優(yōu)化算法等相繼應(yīng)用到可靠指標的優(yōu)化求解中，但其在約束條件上的優(yōu)勢并未受到重視。

本文對高速鐵路隧道二次襯砌(下文中襯砌均指二次襯砌)承載能力極限狀態(tài)各功能函數(shù)的約束條件進行整理，通過抽樣計算得到襯砌力學響應(yīng)的統(tǒng)計特征，分析隧道襯砌受力狀態(tài)的統(tǒng)計特征，建立有約束的可靠指標優(yōu)化求解模型，分別應(yīng)用到素混凝土和鋼筋混凝土襯砌，并與JC法進行比較。

1 功能函數(shù)可行域的判別

1.1 素混凝土襯砌的功能函數(shù)

素混凝土襯砌有抗裂和抗壓2種極限狀態(tài)。對于有抗裂要求的素混凝土矩形截面偏心受壓構(gòu)件，參考地面建筑結(jié)構(gòu)的混凝土結(jié)構(gòu)設(shè)計相關(guān)文獻[9]，考慮裂縫出現(xiàn)瞬間受拉區(qū)混凝土的塑性特征，建立素混凝土矩形構(gòu)件抗裂驗算表達式； 20世紀90年代，原鐵道部建設(shè)司工程建設(shè)科研項目通過大量強度試驗，引入偏心影響系數(shù)α建立了素混凝土襯砌偏心抗壓承載力驗算表達式[10]。素混凝土襯砌的2種極限狀態(tài)功能函數(shù)見式(1)—(2)。

(1)

(2)

式(1)—(2)中：γ為抵抗矩塑性影響系數(shù)；φ為構(gòu)件縱向彎曲系數(shù)；ft為混凝土抗拉強度；b為截面寬度；h為截面高度；e0為軸力偏心距，e0=M/N，其中M為彎矩，N為軸力；α為軸向力偏心影響系數(shù)；fc為混凝土抗壓強度。

根據(jù)素混凝土抗壓、抗裂驗算表達式，得到Ⅲb型素混凝土襯砌抗壓、抗裂承載力隨偏心距變化示意圖(見圖 1)。由圖1可知： 當偏心距較小時，抗裂承載力呈現(xiàn)劇烈波動，但此時拉應(yīng)力并不是混凝土破壞的因素，因此不需要抗裂驗算； 當偏心距在0.2h附近，抗裂與抗壓承載力接近，之后抗裂承載力逐漸下降且低于抗壓承載力。

圖1 素混凝土襯砌抗力變化曲線圖

《暫行規(guī)范》[6]中規(guī)定： 僅當偏心距不小于0.2倍截面高度時，才進行抗裂驗算。選取0.2h附近交點為臨界偏心距e0*，則素混凝土襯砌的抗裂狀態(tài)功能函數(shù)的約束條件為e0≥e0*，抗壓狀態(tài)功能函數(shù)的約束條件為e0≤e0*。

1.2 鋼筋混凝土襯砌的功能函數(shù)

參考《暫行規(guī)范》[6],對鋼筋混凝土襯砌的極限狀態(tài)方程進行換算，消去相對受壓區(qū)高度，以極限軸力、材料性能參數(shù)和截面幾何參數(shù)表示鋼筋混凝土襯砌功能函數(shù)，見式(3)—(4)。

(3)

(4)

σs根據(jù)相對受壓區(qū)高度分別取值，見式(5)。

(5)

式中：fy為鋼筋抗拉強度；ξb為界限相對受壓區(qū)高度；β1為等效矩形應(yīng)力的高度系數(shù)； 3個參數(shù)都只與材料特性有關(guān)。

本文采用界限偏心距e0b來區(qū)分大小偏心極限狀態(tài)，其依據(jù)與計算公式見文獻[11]。鋼筋混凝土襯砌的大、小偏心極限狀態(tài)功能函數(shù)的約束條件分別為e0≥e0b、e0≤e0b。

2 隧道襯砌力學響應(yīng)的統(tǒng)計分析

高速鐵路隧道工程建設(shè)中對不同圍巖條件、設(shè)計行車速度的雙線隧道采取不同的二次襯砌設(shè)計參數(shù)，不同速度、圍巖條件下襯砌結(jié)構(gòu)和材料的類型見表1。

表1不同速度、圍巖條件下的襯砌和材料類型

Table 1 Types of lining and materials under different speeds of train and surrounding rocks

速度/(km/h)圍巖等級襯砌類型材料350、250ⅡⅡa、Ⅱb素混凝土ⅢⅢa、Ⅲb素混凝土ⅣⅣa、Ⅳb鋼筋混凝土ⅤⅤa、Ⅴb鋼筋混凝土

本文選取時速350 km高速鐵路雙線隧道Ⅲa型(二次襯砌為素混凝土結(jié)構(gòu)，拱墻厚40 cm，仰拱厚50 cm)、Ⅳb型(二次襯砌為鋼筋混凝土，拱墻厚45 cm，仰拱厚55 cm)襯砌，采用荷載-結(jié)構(gòu)模型，參考相關(guān)研究對材料性能、圍巖物理力學參數(shù)等的統(tǒng)計數(shù)據(jù)[7]，編制Monte-Carlo隨機有限元計算程序。

從素混凝土和鋼筋混凝土襯砌中，各選取從拱頂?shù)窖龉爸行牡?1個截面，對其力學響應(yīng)特征進行統(tǒng)計，對以偏心距表示的受力狀態(tài)進行分析，截面編號如圖2所示。圖3示出Ⅲa型襯砌拱頂處軸力、彎矩的均值和變異系數(shù)隨抽樣次數(shù)增加的變化趨勢，當抽樣次數(shù)達9 000 時，襯砌軸力和彎矩的均值和變異系數(shù)已經(jīng)趨于穩(wěn)定。經(jīng)抽樣發(fā)現(xiàn)，Ⅲa型和Ⅳb型襯砌各截面的力學響應(yīng)特征均符合此規(guī)律，故確定隨機抽樣的樣本容量為 9 000次。圖 4(a)和4(b)分別示出Ⅲa型襯砌拱頂處軸力和彎矩的分布直方圖、正態(tài)和對數(shù)正態(tài)分布擬合曲線及其擬合優(yōu)度。由圖4可知，正態(tài)和對數(shù)正態(tài)分布的擬合優(yōu)度接近，并逼近于1，說明擬合效果均較好?？紤]到計算的簡便性和2種分布的特點，假定軸力和彎矩服從正態(tài)分布，其抽樣置信水平為95%。

圖2 襯砌單元編號示意圖

圖5示出Ⅲa、Ⅳb型襯砌截面偏心距平均值和變異系數(shù)隨截面位置的變化曲線。由圖5可知： 2種襯砌的偏心距均值及變異系數(shù)的變化趨勢大致相同，從拱頂?shù)焦暗壮尸F(xiàn)波浪式變化，拱頂偏心距最大，其次為拱腰附近，Ⅳb型襯砌的偏心距均值稍大于Ⅲa型襯砌； 2種襯砌各截面偏心距的變異系數(shù)均低于5‰，基本分布在均值附近。各截面偏心距的低變異性表明了極限狀態(tài)出現(xiàn)概率的兩極分化，可以從側(cè)面說明喻渝等[7]判別極限狀態(tài)功能函數(shù)的方法具有一定的合理性，但缺乏對約束條件的明確表達，仍然不夠嚴謹和準確。

3 可靠指標的優(yōu)化求解

3.1 優(yōu)化計算模型的建立

(6)

式中：σX為隨機變量X的標準差；μX為隨機變量X的均值。

(a) 軸力均值

(b) 彎矩均值

(c) 變異系數(shù)

Fig. 3 Variation curves of mean value of axial force and bending moment and coefficient of variation of lining of type Ⅲa

由結(jié)構(gòu)可靠度理論可知，可靠指標的幾何意義是：標準正態(tài)坐標空間中原點到極限狀態(tài)曲面的最短距離[1]。根據(jù)材料強度參數(shù)、幾何參數(shù)和抽樣計算所得彎矩、軸力的分布類型與統(tǒng)計特征(均值、標準差)，通過坐標變換，可以將極限狀態(tài)函數(shù)及其約束條件分別表示為優(yōu)化求解模型中的等式約束與不等式約束，目標函數(shù)則為標準正態(tài)空間中功能函數(shù)曲面到原點的最小距離，如式(7)所示。

(7)

(b) 彎矩

Fig. 4 Distribution fitting of mechanical response of lining of type Ⅲa

圖5 Ⅲa、Ⅳb型襯砌截面偏心距平均值和變異系數(shù)變化曲線

Fig. 5 Variation curves of mean value and coefficient of variation of eccentricity of lining sections of type Ⅲaand Ⅳb

以鋼筋混凝土襯砌為例說明優(yōu)化方法相對于JC法的優(yōu)點。圖6示出鋼筋混凝土襯砌大、小極限狀態(tài)功能函數(shù)在軸力、彎矩空間的投影，其中AOB、COD分別對應(yīng)大偏心和小偏心極限狀態(tài)，僅AOD段代表實際極限狀態(tài)，OB、OC段為無效段，界限偏心距e0b=cotα。當采用JC法計算可靠指標時，由于驗算點在一個固定的功能函數(shù)上持續(xù)迭代，驗算點位置有可能出現(xiàn)在OB、OC段，使得計算結(jié)果有誤； 同時，采用JC法對2種極限狀態(tài)計算得到的2個可靠指標，仍要驗算其是否符合約束條件，過程復雜。優(yōu)化求解方法則可通過約束條件將求解域限定在AO、OD段，所求最小值即為最終求解值，從而優(yōu)化求解值也必將大于等于JC法迭代計算的最小值。

圖6 功能函數(shù)空間示意圖

3.2 素混凝土襯砌可靠指標計算

利用前面抽樣計算得到的襯砌內(nèi)力響應(yīng)的統(tǒng)計特征和喻渝等[7]得出的C30混凝土襯砌厚度和混凝土抗壓、抗拉強度的統(tǒng)計特征(見表2)，以彎矩、軸力、襯砌厚度和混凝土抗壓、抗拉強度為隨機變量，計算Ⅲa型襯砌各截面可靠指標。

表2襯砌性能與幾何參數(shù)統(tǒng)計(C30)[7]

Table 2 Statistics of characteristics and geometric parameters of lining (C30)[7]

參數(shù)軸心抗壓強度/MPa軸心抗拉強度/MPa襯砌厚度/m均值25.982.860.4標準差3.640.400.02分布類型正態(tài)分布正態(tài)分布正態(tài)分布

為測試優(yōu)化算法性能，采用Matlab優(yōu)化工具箱中的Global Search函數(shù)和JC法分別對素混凝土襯砌抗裂和抗壓極限狀態(tài)的可靠指標進行求解計算。Global Search是一種從基于分散搜索生成的試驗點開始，分別進行搜索的全局最優(yōu)化算法，適用于目標函數(shù)和約束函數(shù)均連續(xù)的優(yōu)化問題。在算法測試中，對所有功能函數(shù)均不設(shè)置約束，求解目標均為可靠指標最小值。

素混凝土襯砌2種方法的可靠指標計算結(jié)果見圖7，當偏心距較小時，抗裂承載力較小，采用JC 法計算的抗裂可靠指標為負值； 而優(yōu)化求解計算的是驗算點到原點的空間距離，為正值，因此兩者會出現(xiàn)互為相反數(shù)的情況，本文取其絕對值進行比較。圖7(a)中，Global Search函數(shù)優(yōu)化求解值與JC法求解值完全重合； 圖 7(b)中，JC法和Global Search函數(shù)求解值也完全重合，說明Global Search函數(shù)尋優(yōu)性能良好，能夠滿足素混凝土的可靠指標求解。因此，采用Global Search函數(shù)對可靠指標進行優(yōu)化求解(有約束)，并與 JC法求解值(無約束)作對比，結(jié)果如圖8所示，Global Search優(yōu)化值與JC法求解值分段重合。以上計算結(jié)果說明： 優(yōu)化方法中功能函數(shù)約束條件的引入，解決了采用JC法時功能函數(shù)界限模糊、計算值難以區(qū)分的問題，實現(xiàn)了襯砌各截面可靠指標的準確求解。

(a) 抗裂

(b) 抗壓

圖8 素混凝土襯砌可靠指標優(yōu)化求解計算結(jié)果

3.3 鋼筋混凝土襯砌可靠指標計算

利用前面抽樣計算得到的襯砌內(nèi)力響應(yīng)的統(tǒng)計特征和喻渝等[7]得出的C35混凝土襯砌厚度、混凝土抗壓強度和鋼筋抗拉強度的統(tǒng)計特征(見表3)，以彎矩、軸力、襯砌厚度、混凝土抗壓強度和鋼筋抗拉強度為隨機變量，對Ⅳb型襯砌各截面可靠指標進行計算。

表3襯砌性能與幾何參數(shù)統(tǒng)計(C35)[7]

Table 3 Statistics of characteristics and geometric parameters of lining (C35)[7]

參數(shù)混凝土軸心抗壓強度/MPa鋼筋抗拉強度/MPa襯砌厚度/m均值29.89353.80.45標準差3.8912.030.022 5分布類型正態(tài)分布正態(tài)分布正態(tài)分布

同樣采用Global Search函數(shù)和JC法分別對鋼筋混凝土襯砌的可靠指標進行求解計算，結(jié)果如圖9所示。由圖 9可知： 對于2種功能函數(shù)，JC法和Global Search函數(shù)優(yōu)化求解值完全重合。因此，鋼筋混凝土襯砌的可靠指標優(yōu)化仍選用Global Search函數(shù)。

(a) 小偏心狀態(tài)

(b) 大偏心狀態(tài)

Fig. 9 Comparison of calculation methods forRIof reinforced concrete lining

鋼筋混凝土襯砌可靠指標優(yōu)化求解計算結(jié)果如圖10所示，由圖10可知： 1—39號截面的優(yōu)化求解結(jié)果與不同極限狀態(tài)下的JC法計算結(jié)果分段重合； 40號截面之后的優(yōu)化求解值與JC法計算值不重合，但不小于JC法計算的最小值。分析可得： 求解值重合是因為該處迭代終止點并未落在失效段，且可靠指標最小，優(yōu)化法成功搜尋到最優(yōu)解并與其重合； 不重合的原因是由于缺乏約束，JC法迭代終止點落在了失效段，JC法計算值始終為該功能函數(shù)所求的最小值，而優(yōu)化法僅在可行域之內(nèi)搜尋，故JC法計算最小值結(jié)果小于優(yōu)化值。計算結(jié)果一方面表明考慮極限狀態(tài)約束條件的優(yōu)化求解方法在鋼筋混凝土襯砌可靠指標計算中取得成功應(yīng)用，另一方面也說明采用JC法分別計算2種極限狀態(tài)的可靠指標，再取其中最小值的方法是不可取的。

圖10 鋼筋混凝土襯砌可靠指標優(yōu)化求解計算結(jié)果

Fig. 10 Optimization solution results forRIof reinforced concrete lining

4 結(jié)論與討論

4.1 結(jié)論

本文通過采用隨機有限元方法、Monte-Carlo抽樣計算、優(yōu)化函數(shù)的調(diào)用和對相關(guān)研究的歸納整理等手段，對隧道襯砌受力狀態(tài)和可靠度計算方法進行研究，得到以下結(jié)論。

1)綜合素混凝土襯砌和鋼筋混凝土襯砌極限狀態(tài)功能函數(shù)的相關(guān)研究，給出抗裂與抗壓、大偏心與小偏心極限狀態(tài)的判別條件。

2)對時速350 km雙線隧道Ⅲa、Ⅳb襯砌進行抽樣計算，考慮2類襯砌各截面偏心距的低變異性，分析目前可靠指標計算中極限狀態(tài)功能函數(shù)選用方法的合理性及不足。

3)建立可靠指標的優(yōu)化求解模型，證明了Global Search函數(shù)尋優(yōu)性能較好，能夠考慮各極限狀態(tài)功能函數(shù)的適用條件，消除極限狀態(tài)功能函數(shù)選取的模糊性，實現(xiàn)了對隧道襯砌各截面可靠指標的準確快速求解。

本文研究采用優(yōu)化求解方法，解決了JC法求解隧道襯砌可靠指標中極限狀態(tài)功能函數(shù)選取的問題，與目前采用偏心距或受壓區(qū)高度均值的判別方法相比，更具有嚴謹性和便捷性； 同時，證明了用JC法分別計算2種極限狀態(tài)然后取其最小值的不合理性。

4.2 存在問題及建議

隧道賦存于復雜的地質(zhì)條件中，圍巖參數(shù)及襯砌與圍巖的相互作用模式具有高度的不確定性。本文采用荷載-結(jié)構(gòu)模型和根據(jù)當前隧道規(guī)范提供的地層荷載計算方法，并未考慮計算模型的不確定性和圍巖物理力學參數(shù)的空間變異性。此外，“單失效模式”的可靠性分析方法雖然被廣泛采用，但對于隧道這種高次超靜定結(jié)構(gòu)，其評價指標過于單一，理論基礎(chǔ)也不盡合理； 而對于隧道結(jié)構(gòu)多失效模式的可靠指標計算(如點估計法、界限范圍估計法等)，優(yōu)化方法能夠排除不同極限狀態(tài)的干擾，無需進行極限狀態(tài)判別和計算值的驗算與比選，從而實現(xiàn)快速準確求解，適用性更強。

在當前的研究現(xiàn)狀下，描述和考慮圍巖物理力學參數(shù)及力學計算模型的不確定性，建立基于隧道-圍巖體系的失效模式及可靠性評價體系仍是隧道可靠性研究的重點和難點。本文僅對可靠指標的計算方法進行部分改進，下一步應(yīng)緊緊圍繞以上幾個基本問題，基于目前隧道失效模式、隧道與圍巖相互作用關(guān)系及其不確定性描述的最新研究成果，發(fā)展隧道可靠性評價體系，進一步強化可靠性計算的手段與方法。