武建力 高 抗

(中交第二公路勘察設(shè)計(jì)研究院有限公司 武漢 430056)

隨著國(guó)內(nèi)盾構(gòu)技術(shù)的成熟，大直徑盾構(gòu)法在水下隧道工程中得到了廣泛應(yīng)用。武漢三陽(yáng)路長(zhǎng)江隧道、上海長(zhǎng)江隧道、杭州錢江隧道等超大直徑盾構(gòu)隧道越來(lái)越多。盾構(gòu)隧道管片的結(jié)構(gòu)耐久性問(wèn)題關(guān)系到隧道使用壽命、舒適度等多方面，結(jié)構(gòu)使用環(huán)境能明顯地導(dǎo)致結(jié)構(gòu)材料性能隨使用時(shí)間的延長(zhǎng)而劣化，因此水下大直徑盾構(gòu)隧道管片耐久性設(shè)計(jì)是結(jié)構(gòu)設(shè)計(jì)必不可少的重要內(nèi)容。

管片襯砌結(jié)構(gòu)在投入使用后，受到混凝土碳化、氯離子侵蝕等各種環(huán)境因素的作用，最終會(huì)引起鋼筋銹蝕，混凝土截面損傷，以及鋼筋與混凝土之間的黏結(jié)性能退化，從而導(dǎo)致隧道襯砌結(jié)構(gòu)承載能力逐步降低，結(jié)構(gòu)失效。因此，在襯砌結(jié)構(gòu)耐久性損傷過(guò)程中，環(huán)境因素起了決定性的作用，而地下復(fù)雜環(huán)境也給混凝土的耐久性帶來(lái)了極大的挑戰(zhàn)[1]。

隧道襯砌結(jié)構(gòu)耐久性是與時(shí)間相關(guān)的，是由襯砌結(jié)構(gòu)外界環(huán)境和本身因素共同作用下所達(dá)到的使工程結(jié)構(gòu)正常使用的年限[2]。水下隧道襯砌主體結(jié)構(gòu)設(shè)計(jì)使用年限一般為100年，有的隧道甚至提出更高的年為(港珠澳大橋設(shè)計(jì)使用年限為120年)。因此，為了保證水下大直徑盾構(gòu)隧道100年的設(shè)計(jì)使用年限，有必要對(duì)隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)進(jìn)行研究，不僅具有重大的經(jīng)濟(jì)意義，也可為今后類似水下盾構(gòu)隧道耐久性設(shè)計(jì)提供理論依據(jù)。

1 大直徑盾構(gòu)隧道結(jié)構(gòu)耐久性設(shè)計(jì)

李田、劉西拉提出的基于近似概率的耐久性設(shè)計(jì)方法，從計(jì)算和構(gòu)造要求兩部分進(jìn)行耐久性設(shè)計(jì)[3-4]。在計(jì)算和驗(yàn)算過(guò)程中，采用現(xiàn)行結(jié)構(gòu)設(shè)計(jì)規(guī)范GB 50010-2010(2015版)《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》中極限設(shè)計(jì)表達(dá)式中加入耐久性設(shè)計(jì)系數(shù)的方法，與該規(guī)范相協(xié)調(diào)。該方法概念明確，形式簡(jiǎn)單，而且易被工程設(shè)計(jì)人員熟悉和掌握。同時(shí)，考慮了抗力隨時(shí)間衰減的變化規(guī)律，使新建結(jié)構(gòu)有明確的目標(biāo)使用年限，達(dá)到結(jié)構(gòu)安全、經(jīng)濟(jì)、耐久和實(shí)用的建設(shè)目的。基于此，本文建立的水下大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)方法就是基于近似概率法的耐久性設(shè)計(jì)方法，同時(shí)結(jié)合DuraCrete方法，將參數(shù)考慮為隨機(jī)變量。

1.1 基于近似概率法的耐久性設(shè)計(jì)方法

現(xiàn)行的GB 50010-2010(2015版)混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范中關(guān)于計(jì)算與驗(yàn)算部分的公式可以歸納為式(1)所示形式。

S≤R

(1)

式中：S為結(jié)構(gòu)內(nèi)力(作用效應(yīng))設(shè)計(jì)值；R為結(jié)構(gòu)抗力設(shè)計(jì)值。

基于式(1)，可得混凝土結(jié)構(gòu)耐久性設(shè)計(jì)的極限狀態(tài)功能函數(shù)隨機(jī)過(guò)程表達(dá)式如式(2)。

Z(t)=R(t)-S(t)

(2)

式中：Z(t)為混凝土結(jié)構(gòu)極限狀態(tài)功能函數(shù)隨機(jī)過(guò)程；R(t)為混凝土結(jié)構(gòu)抗力衰減隨機(jī)過(guò)程；S(t)為混凝土結(jié)構(gòu)荷載效應(yīng)隨機(jī)過(guò)程。

結(jié)構(gòu)構(gòu)件在內(nèi)外因素的共同作用下，隨著時(shí)間的增長(zhǎng)會(huì)發(fā)生耐久性損傷，結(jié)構(gòu)承載力下降，即結(jié)構(gòu)抗力R(t)是隨時(shí)間增長(zhǎng)而降低的隨機(jī)過(guò)程；而荷載效應(yīng)S(t)在使用期內(nèi)可能會(huì)發(fā)生變化，但由于構(gòu)件受力方式、構(gòu)件結(jié)構(gòu)形式和荷載作用形式基本恒定，所以在正常情況下變化的幅度很小。因此，在考慮耐久性損傷的情況下，引入結(jié)構(gòu)安全系數(shù)和耐久性設(shè)計(jì)系數(shù)，則基于近似概率的耐久性設(shè)計(jì)極限狀態(tài)方程見(jiàn)式(3)。

r0S(t)≤ηR(t)

(3)

式中：r0為結(jié)構(gòu)安全系數(shù)；η為耐久性設(shè)計(jì)系數(shù)，是可靠度指標(biāo)的函數(shù)，可表示為

(4)

式中：β為規(guī)范[5]中規(guī)定的可靠度指標(biāo)；βt為達(dá)到使用年限時(shí)的可靠度指標(biāo)，一般取βt=β0；β(t)為結(jié)構(gòu)使用期限內(nèi)可靠度指標(biāo)下降的過(guò)程，是時(shí)間t的函數(shù)。

綜上所述，進(jìn)行結(jié)構(gòu)耐久性設(shè)計(jì)主要是求得結(jié)構(gòu)抗力R(t)衰減的隨機(jī)過(guò)程、結(jié)構(gòu)荷載效應(yīng)S(t)的隨機(jī)過(guò)程，及可靠度指標(biāo)β(t)。

1.2 大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)流程

大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)流程圖見(jiàn)圖1。

圖1 大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)流程圖

2 大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)耐久性計(jì)算

2.1 襯砌結(jié)構(gòu)抗力衰減統(tǒng)計(jì)特征

盾構(gòu)隧道管片結(jié)構(gòu)，其受力形態(tài)一般為偏心受力。近年來(lái)，國(guó)內(nèi)外學(xué)者對(duì)銹后鋼筋混凝土構(gòu)件的承載力進(jìn)行了大量的研究，研究結(jié)果表明，鋼筋銹蝕對(duì)構(gòu)件承載力的影響主要體現(xiàn)在3個(gè)方面：①鋼筋截面積的減小；②鋼筋屈服強(qiáng)度的降低；③鋼筋與混凝土之間黏結(jié)性能的退化。考慮以上3個(gè)方面，本文在計(jì)算鋼筋銹蝕后隧道襯砌結(jié)構(gòu)抗力時(shí)，分別引入3個(gè)損傷系數(shù)[6]，分別反映鋼筋銹蝕引起的鋼筋截面損失及強(qiáng)度降低、混凝土截面損傷、鋼筋和混凝土間黏結(jié)力下降三方面的損傷效應(yīng)。

因?yàn)榻Y(jié)構(gòu)抗力是隨服務(wù)年限的增加不斷降低的，所以將損傷系數(shù)看成隨機(jī)變量。同時(shí)考慮混凝土強(qiáng)度、鋼筋強(qiáng)度、鋼筋銹蝕率等隨機(jī)變量，從而求得隧道襯砌結(jié)構(gòu)抗力衰減計(jì)算模型。

根據(jù)可靠度原理，隨機(jī)變量的統(tǒng)計(jì)特征為均值和標(biāo)準(zhǔn)差。通過(guò)計(jì)算結(jié)構(gòu)的各隨機(jī)變量參數(shù)的統(tǒng)計(jì)特征，可得到小偏心受壓襯砌結(jié)構(gòu)抗力R(t)統(tǒng)計(jì)特征見(jiàn)式(5)。

μRN(t)=RN(t)|μ

(5)

(6)

式中：μRN(t)為結(jié)構(gòu)抗力平均強(qiáng)度；σRN(t)為結(jié)構(gòu)抗力標(biāo)準(zhǔn)差。

2.2 荷載效應(yīng)S(t)統(tǒng)計(jì)特征

對(duì)大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)而言，由于其受力變形的機(jī)理相當(dāng)復(fù)雜，且缺少基本變量的統(tǒng)計(jì)資料，因此幾乎不能采用解析法求得作用效應(yīng)的顯式解。針對(duì)這些問(wèn)題，國(guó)內(nèi)外已進(jìn)行了大量的工作，目前通用的做法是采用隨機(jī)有限元法，該方法概念明確，計(jì)算量少，但需要對(duì)確定性有限元程序進(jìn)行改造，且對(duì)高度非線性問(wèn)題，存在較大的計(jì)算誤差。為解決以上問(wèn)題，張彌等[7]首先將響應(yīng)面法引進(jìn)到隧道結(jié)構(gòu)可靠度分析中，用以計(jì)算襯砌結(jié)構(gòu)的荷載效應(yīng)。

本文主要采用傳統(tǒng)的有限元響應(yīng)面法，將結(jié)構(gòu)高度的隱式非線性極限狀態(tài)方程用一個(gè)簡(jiǎn)單的響應(yīng)面函數(shù)進(jìn)行替代，即建立荷載效應(yīng)與隨機(jī)變量之間的一個(gè)顯式表達(dá)式，然后以該顯式表達(dá)式作為響應(yīng)面函數(shù)，利用Monte-Carlo方法進(jìn)行抽樣計(jì)算，最后得到荷載效應(yīng)的統(tǒng)計(jì)特征。

計(jì)算襯砌結(jié)構(gòu)荷載效應(yīng)的統(tǒng)計(jì)特征主要是利用ANSYS計(jì)算軟件中的可靠度計(jì)算模塊，首先利用“荷載-結(jié)構(gòu)”法生成可靠度分析文件，然后利用中心點(diǎn)復(fù)合設(shè)計(jì)擬合大直徑盾構(gòu)隧道襯砌管片結(jié)構(gòu)截面上彎矩和軸力的響應(yīng)面方程，最后利用Monte-Carlo方法進(jìn)行抽樣計(jì)算，最后得到該截面彎矩和軸力的均值和方差。

以某大直徑盾構(gòu)隧道為例，隧道的外徑為14.5 m、內(nèi)徑為13.3 mm，選取江中埋深約17.5 m處襯砌管片單位寬度進(jìn)行荷載效應(yīng)的統(tǒng)計(jì)特征計(jì)算，用響應(yīng)面法擬合截面上的彎矩和軸力的表達(dá)式，并進(jìn)行敏感性分析。因荷載效應(yīng)影響因素較為復(fù)雜，所以本文只將幾個(gè)主要參數(shù)看成隨機(jī)變量，其余參數(shù)假定為定值，主要參數(shù)的統(tǒng)計(jì)特征[8]見(jiàn)表1。

表1 影響襯砌結(jié)構(gòu)荷載效應(yīng)隨機(jī)參數(shù)統(tǒng)計(jì)特征表

通過(guò)“荷載-結(jié)構(gòu)”法，荷載取隧道上覆水土壓力，按松散土壓力理論計(jì)算，得到該截面的襯砌結(jié)構(gòu)內(nèi)力圖見(jiàn)圖2。

圖2 襯砌結(jié)構(gòu)內(nèi)力圖

根據(jù)圖2的計(jì)算結(jié)果提取軸力(受壓)的最大值，利用ANSYS有限元的可靠度分析模塊，基于響應(yīng)面法計(jì)算襯砌結(jié)構(gòu)軸力的均值和標(biāo)準(zhǔn)差，其襯砌結(jié)構(gòu)軸力分布圖見(jiàn)圖3。

圖3 襯砌結(jié)構(gòu)軸力分布圖

由圖3可見(jiàn)，襯砌結(jié)構(gòu)軸力的均值為μN(yùn)=6 442.2 kPa，標(biāo)準(zhǔn)差為σN=32.942 kPa。

根據(jù)響應(yīng)面有限元計(jì)算，對(duì)軸力影響因素：管片混凝土彈性模量E、管片混凝土密度DENS、管片厚度H、地基彈性抗力系數(shù)K4個(gè)輸入隨機(jī)變量進(jìn)行敏感性分析，顯著性水平2.5%下，軸力影響因素敏感性分析結(jié)果對(duì)比餅狀圖見(jiàn)圖4。

圖4 軸力影響因素敏感性分析圖

由圖4可見(jiàn)，對(duì)于影響軸力的4個(gè)隨機(jī)因素中，地基彈性抗力系數(shù)對(duì)襯砌結(jié)構(gòu)軸力的影響是最大的，其次是管片密度、彈性模量和管片厚度。

2.3 襯砌結(jié)構(gòu)時(shí)變可靠度指標(biāo)計(jì)算

大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)極限狀態(tài)方程是高度非線性的，直接采用一般的可靠度計(jì)算方法得到的是幾何可靠度，存在較大誤差，需進(jìn)行改進(jìn)?？刹捎没诮Y(jié)構(gòu)風(fēng)險(xiǎn)最小化原理的支持向量機(jī)回歸(SVR)和響應(yīng)面法相結(jié)合的方法，利用支持向量機(jī)回歸對(duì)小樣本數(shù)據(jù)良好的學(xué)習(xí)和泛化能力，用SVR重構(gòu)結(jié)構(gòu)響應(yīng)面方程建立基于SVR和響應(yīng)面結(jié)合的結(jié)構(gòu)時(shí)變可靠度的計(jì)算方法。

通過(guò)結(jié)合室內(nèi)試驗(yàn)和經(jīng)驗(yàn)公式得到結(jié)構(gòu)原始極限狀態(tài)方程及其參數(shù)的統(tǒng)計(jì)特征，通過(guò)二次規(guī)劃，得到利用支持向量回歸機(jī)重構(gòu)的響應(yīng)面極限方程，以響應(yīng)面方程作為近似功能函數(shù)進(jìn)行抽樣，考慮抗力衰減的時(shí)變性，從而計(jì)算結(jié)構(gòu)的時(shí)變可靠度。主要流程如下。

1) 假定初始迭代點(diǎn)x(t)=[x1(t),x2(t),…,xn(t)]T，一般取平均值μX(t)。

2) 通過(guò)結(jié)構(gòu)數(shù)值分析或試驗(yàn)在各個(gè)展開(kāi)點(diǎn)處計(jì)算功能函數(shù)的估計(jì)值gi(t)(i=1,2,…,2n+1)，并形成相應(yīng)的系數(shù)矩陣。

3) 將上述2n+1個(gè)樣本點(diǎn)作為訓(xùn)練樣本集，并對(duì)輸入樣本和輸出樣本進(jìn)行數(shù)據(jù)歸一化處理。

4) 根據(jù)輸入和輸出樣本，利用交叉驗(yàn)證法和試算法選擇最優(yōu)的SVR參數(shù)C和r。

5) 構(gòu)造核函數(shù)K矩陣，將歸一化后的2n+1個(gè)樣本點(diǎn)代入非線性求解問(wèn)題的內(nèi)積核函數(shù)，求解二次規(guī)劃問(wèn)題，求得拉格朗日乘子及常數(shù)項(xiàng)b的數(shù)值，從而求得SVR響應(yīng)面方程式。

6) 通過(guò)SVR響應(yīng)面方程求得設(shè)計(jì)驗(yàn)算點(diǎn)和可靠度指標(biāo)，并判斷其變化幅度。

|βk(t)-βk-1(t)|是否滿足所需精度，若不滿足，則返回第(3)步，重新開(kāi)始迭代計(jì)算。

7) 根據(jù)可靠度指標(biāo)與失效概率、耐久性系數(shù)之間的轉(zhuǎn)換關(guān)系，分別求出結(jié)構(gòu)服務(wù)年限內(nèi)失效概率與耐久性系數(shù)。

2.4 襯砌結(jié)構(gòu)時(shí)變可靠度指標(biāo)計(jì)算

基于上述程序理論，對(duì)某大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)進(jìn)行耐久性設(shè)計(jì)，利用MATLAB軟件，計(jì)算結(jié)構(gòu)的時(shí)變可靠度指標(biāo)及失效概率，從而計(jì)算結(jié)構(gòu)耐久性設(shè)計(jì)系數(shù)η(t)和結(jié)構(gòu)安全系數(shù)，當(dāng)?shù)竭_(dá)結(jié)構(gòu)100年服務(wù)壽命時(shí)，確定其對(duì)應(yīng)的耐久性設(shè)計(jì)系數(shù)，從而據(jù)此建立隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)的極限狀態(tài)方程見(jiàn)式(7)。

γ0S≤η(100)R

(7)

式中：γ0為結(jié)構(gòu)重要性系數(shù)，本文取1.1；S為100年服務(wù)壽命內(nèi)襯砌結(jié)構(gòu)的軸力，本文取初始軸力的均值；η(100)為100年服務(wù)壽命時(shí)襯砌結(jié)構(gòu)耐久性設(shè)計(jì)系數(shù)，根據(jù)程序進(jìn)行計(jì)算；R為100年服務(wù)壽命時(shí)襯砌結(jié)構(gòu)抗力，按程序進(jìn)行計(jì)算。

依據(jù)以上極限狀態(tài)方程可以改進(jìn)隧道襯砌結(jié)構(gòu)設(shè)計(jì)，從而確保隧道襯砌結(jié)構(gòu)主體結(jié)構(gòu)的設(shè)計(jì)使用年限的順利實(shí)現(xiàn)。實(shí)例計(jì)算如下。

隧道盾構(gòu)管片環(huán)寬2 m，管片厚度60 cm，采用雙面楔形形式。管片混凝土標(biāo)號(hào)C60，保護(hù)層厚度為50 mm。

計(jì)算結(jié)構(gòu)失效概率、結(jié)構(gòu)可靠度指標(biāo)、結(jié)構(gòu)耐久性設(shè)計(jì)系數(shù)圖，結(jié)果見(jiàn)圖5～圖7。

圖5 某大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)失效概率計(jì)算圖

圖6 某大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)可靠度指標(biāo)計(jì)算圖

圖7 某大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)系數(shù)計(jì)算圖

從圖5～圖7中可知，某大直徑盾構(gòu)隧道建成投入使用100年后，隧道襯砌結(jié)構(gòu)可靠度指標(biāo)為2.047 7，耐久性設(shè)計(jì)系數(shù)為0.691 3，失效概率為2.03%，滿足GB/T 50476-2019 《混凝土結(jié)構(gòu)耐久性設(shè)計(jì)標(biāo)準(zhǔn)》失效概率不大于5%～10%的要求，說(shuō)明盾構(gòu)隧道100年的設(shè)計(jì)使用年限可以得到保證。同時(shí)由100年后耐久性設(shè)計(jì)系數(shù)，可以建立大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)耐久性設(shè)計(jì)的極限狀態(tài)方程如式(8)。

γ0S≤0.693 1R

(8)

根據(jù)極限狀態(tài)方程的計(jì)算，可知某大直徑盾構(gòu)隧道襯砌主體結(jié)構(gòu)能滿足結(jié)構(gòu)100年服務(wù)壽命的要求。

3 結(jié)論

本章通過(guò)采用支持向量機(jī)(SVR)與響應(yīng)面相結(jié)合的方法作為大直徑盾構(gòu)襯砌結(jié)構(gòu)耐久性設(shè)計(jì)的計(jì)算方法，基于可靠度理論，利用MATLAB計(jì)算軟件，并在某大直徑盾構(gòu)隧道耐久設(shè)計(jì)中得到了應(yīng)用，主要工作及研究成果如下。

1) 在對(duì)比傳統(tǒng)耐久性設(shè)計(jì)方法和相關(guān)規(guī)范規(guī)定的不足之后，利用SVR和響應(yīng)面相結(jié)合的原理，給出了大直徑盾構(gòu)襯砌結(jié)構(gòu)耐久性設(shè)計(jì)的基本步驟，利用MATLAB編程軟件進(jìn)行計(jì)算。

2) 通過(guò)建立結(jié)構(gòu)抗力和結(jié)構(gòu)作用效應(yīng)之間的極限狀態(tài)方程，從而計(jì)算結(jié)構(gòu)時(shí)變可靠度指標(biāo)和失效概率，進(jìn)而計(jì)算結(jié)構(gòu)耐久性設(shè)計(jì)系數(shù)，進(jìn)行耐久性設(shè)計(jì)。

3) 通過(guò)ANSYS計(jì)算軟件中的可靠度計(jì)算模塊，先根據(jù)“荷載-結(jié)構(gòu)”法計(jì)算結(jié)構(gòu)的內(nèi)力，然后進(jìn)行響應(yīng)面模擬，得到結(jié)構(gòu)內(nèi)力的均值和方差，并進(jìn)行影響結(jié)構(gòu)內(nèi)力的參數(shù)敏感性分析。

4) 用SVR重構(gòu)結(jié)構(gòu)響應(yīng)面方程，建立了基于SVR和響應(yīng)面結(jié)合的結(jié)構(gòu)時(shí)變可靠度的計(jì)算方法。利用MATLAB軟件編制了大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)在服務(wù)壽命期限內(nèi)的時(shí)變可靠度指標(biāo)及其失效概率的計(jì)算程序。

5) 根據(jù)耐久性設(shè)計(jì)系數(shù)與時(shí)變可靠度指標(biāo)的相互關(guān)系，計(jì)算隧道襯砌結(jié)構(gòu)達(dá)到100年服務(wù)壽命時(shí)的耐久性設(shè)計(jì)系數(shù)，得到某大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)100年服務(wù)壽命時(shí)的極限狀態(tài)方程：從而得到某大直徑盾構(gòu)隧道襯砌結(jié)構(gòu)滿足100年服務(wù)壽命的結(jié)果。