陳紹軍，張強，楊鈞杰

（1.四川省公路規(guī)劃勘察設(shè)計研究院有限公司,成都610041；2.北京建筑大學(xué)土木與交通工程學(xué)院,北京100044）

1 引言

近年來，大跨懸索橋不斷興建，其跨徑也不斷增加。目前，世界上跨度超過1 000 m的已建成的橋梁均為懸索橋，例如，土耳其的博斯普魯斯一橋和二橋（主跨分別為1 074 m和1 090 m）、中國香港的青馬大橋（主跨為1 377 m）、中國的潤揚長江大橋南汊懸索橋（主跨為1 490 m）等。大跨懸索橋的分析與設(shè)計的理論研究也在不斷地深入，這些研究大部分都只是采用傳統(tǒng)的方法并基于結(jié)構(gòu)的確定性參數(shù)進行確定性的靜力分析。但是，顯而易見，對于懸索橋這樣一種大型、復(fù)雜的結(jié)構(gòu)，包含了紛繁的不確定性，比如荷載、材料、幾何尺寸的不確定性，這些不確定性必然會對懸索橋的可靠性產(chǎn)生很大的影響，同時，由于懸索橋是一種具有強幾何非線性的結(jié)構(gòu)形式，這也會大大增加可靠度分析的計算量。因此，需要發(fā)展一種準(zhǔn)確而有效的懸索橋可靠度計算方法，這既是大跨懸索橋結(jié)構(gòu)可靠度設(shè)計中需要解決的關(guān)鍵問題，也是大跨懸索橋結(jié)構(gòu)健康監(jiān)測與安全評估系統(tǒng)研究的基礎(chǔ)理論問題[1，2]。

本文以某大跨度懸索橋的扁平鋼箱梁結(jié)構(gòu)為分析對象，研究扁平鋼箱梁結(jié)構(gòu)的可靠度分析方法。為了解決傳統(tǒng)響應(yīng)面方法在分析大型復(fù)雜結(jié)構(gòu)可靠度時可能遇到的不收斂或誤差較大的問題，通過對迭代步長取值的改進，提出了基于改進響應(yīng)面的結(jié)構(gòu)可靠度分析方法。在此基礎(chǔ)上分別采用改進響應(yīng)面法和基于響應(yīng)面的蒙特卡羅法分析了鋼箱梁的靜力可靠度。最后，討論了2種計算方法的有效性以及懸索橋結(jié)構(gòu)的幾何非線性對可靠度指標(biāo)的影響。

2 改進響應(yīng)面法

本文采用響應(yīng)面法分析大跨懸索橋扁平鋼箱梁結(jié)構(gòu)的靜力可靠度問題。取不含交叉項的二次多項式作為響應(yīng)面函數(shù)，即：

式中，X為相互獨立的基本隨機變量向量；b0、bi和bii為待定系數(shù)；xi為隨機變量。

可靠度分析時，響應(yīng)面法采用了向驗算點逼近的迭代策略，即首先從變量均值點附近擬合響應(yīng)面函數(shù)，經(jīng)過多次迭代使得響應(yīng)面函數(shù)逐漸逼近和吻合驗算點附近的真實功能函數(shù)。

傳統(tǒng)響應(yīng)面法迭代的具體步驟為[3-6]：

1）假定初始點X(1)=(x1(1)，x2(1)，...，xn(1))，一般取均值點。

2）近似功能函數(shù)g(x1(1)，x2(1)，...，xn(1))以及g(x1(1)，x2(1)，...，xi(1)±fσi，...，xn(1))得到2n+1個點估計值。其中，σi為隨機變量xi的標(biāo)準(zhǔn)差；f為迭代步長。通常第一步迭代時f可取2～3，此后的迭代過程中f可取為1。

3）利用2n+1個點估計值求解式（1）中的待定系數(shù)，從而得到當(dāng)前迭代點處功能函數(shù)的近似極限狀態(tài)方程。

4）求解近似極限狀態(tài)方程的驗算點X(k)以及可靠指標(biāo)β(k)。其中，上標(biāo)k表示第k步迭代。

傳統(tǒng)響應(yīng)面法存在收斂失敗和計算誤差大的缺點，起主要原因在于迭代步長f的取值。當(dāng)?shù)介Lf較大時，由于多項式對實際功能函數(shù)的擬合度較差，導(dǎo)致響應(yīng)面法無法收斂到真實驗算點，而是在其附近陷入局部收斂解；當(dāng)?shù)介Lf很小時，可以使響應(yīng)面函數(shù)與實際功能函數(shù)在真實驗算點吻合良好，計算精度較高。然而，對于非線性程度較大的極限狀態(tài)方程，迭代步長f如果從迭代過程一開始就很小，將造成迭代過程發(fā)生不收斂的現(xiàn)象。因此，響應(yīng)面函數(shù)的擬合空間應(yīng)由大到小逐漸收縮到驗算點附近較小的范圍內(nèi)。根據(jù)上述分析，本文對傳統(tǒng)響應(yīng)面法的迭代格式加以改進，即采用漸近迭代步長序列f的方法，稱為改進的響應(yīng)面法。也就是在傳統(tǒng)響應(yīng)面法的迭代步驟中，第一步迭代時f取為2，此后的迭代過程中f隨迭代步數(shù)逐漸減小直至取為0.01。漸近迭代步長序列f的定義如下：

式中，n為迭代步數(shù)。

3 實橋可靠度分析

3.1 懸索橋有限元分析模型

該橋是由懸索橋和斜拉橋組合而成的特大型纜索支承型橋梁。其中，南汊懸索橋為主跨1 490 m的單跨雙鉸簡支鋼箱梁橋，為中國第一，世界第三。懸索橋橋塔高約210 m，橋主梁采用全焊扁平流線形封閉鋼箱梁斷面，主纜由平行鋼絲索股組成，吊索采用預(yù)制平行鋼絲束股（PWS），外包PE材料進行防護。

懸索橋結(jié)構(gòu)分析常用脊骨梁模型，通常是將復(fù)雜的主梁簡化為一根梁單元。但是，脊骨梁模型只考慮了截面的幾何性質(zhì)，沒有考慮橫向連接系對主梁剛度的影響，即脊骨梁模型不能考慮橫隔板對主梁剛度的作用，對于擁有密集橫隔板的扁平鋼箱梁，模擬將嚴(yán)重失真。同時，懸索橋是一種具有強幾何非線性的橋梁結(jié)構(gòu)形式。這種幾何非線性來自4個方面：（1）重力剛度，即恒載在主纜中產(chǎn)生的內(nèi)力會抵抗后續(xù)荷載所產(chǎn)生的變形；（2）結(jié)構(gòu)的大位移，即整個結(jié)構(gòu)體系的平衡狀態(tài)應(yīng)建立在變形后的狀態(tài)上；（3）主纜自重垂度，即主纜受力后的變形由彈性變形和自重引起的垂度變化組成；（4）主梁和橋塔的P-Δ效應(yīng)（重力二階效應(yīng)）。在懸索橋結(jié)構(gòu)的可靠性分析中這些非線性必須得到考慮。

根據(jù)上述分析，本文對該懸索橋的主纜和吊索采用多直桿法進行模擬，橋塔的模擬則采用三維梁單元。對于扁平鋼箱梁，采取簡化扁平鋼箱梁縱向U形加勁肋的方法，將橋面板及其加勁肋采用網(wǎng)格劃分尺寸較大的正交異性殼單元進行模擬，圖1給出了該懸索橋的有限元分析模型。對于主梁中的橫隔板，也采用板殼單元模擬，這樣形成了橫隔板、底板和腹板等組成空間箱形結(jié)構(gòu)分析模型，采用空間箱形結(jié)構(gòu)分析模型能準(zhǔn)確地模擬扁平鋼箱梁的截面特性和質(zhì)量分布特性，同時可以真實地反映橫隔板、橫截面畸變和剪力滯后等對全橋響應(yīng)分析的影響[7]。

圖1 懸索橋有限元分析模型

3.2 可靠度分析

本文考慮的隨機變量有懸索橋結(jié)構(gòu)的纜、索、主梁和塔的彈性模量Ei、材料重度γi，纜、索的截面面積Ai，簡化的正交異性板的厚度Ti及活荷載q，隨機變量總數(shù)目為20。變量定義如表1所示。

表1 懸索橋隨機變量定義

基于前述懸索橋有限元分析模型及定義的隨機變量，進行了懸索橋結(jié)構(gòu)靜力可靠度分析。根據(jù)JTG/T D65-05—2015《公路懸索橋設(shè)計規(guī)范》，在正常使用極限狀態(tài)，懸索橋加勁梁由汽車荷載（不計沖擊力）引起的最大豎向撓度值不宜大于跨徑的1/250～1/300。本文取為1/300，[δ]=L/300=4 967 mm（L為主跨跨度）。以此建立極限狀態(tài)方程：

式中，dmid為主跨跨中豎向位移；[δ]為容許豎向位移。

表2 分別給出了采用改進響應(yīng)面法和基于響應(yīng)面的蒙特卡羅法（RSM-MCIS法）[8]計算得到的鋼箱梁跨中位移的可靠度指標(biāo)?；陧憫?yīng)面的蒙特卡羅法的基本原理是：采用響應(yīng)面法構(gòu)造近似極限狀態(tài)方程，在已獲得近似的顯式結(jié)構(gòu)極限狀態(tài)方程上，采用蒙特卡羅法得到可靠度指標(biāo)。計算過程表明，由于改進響應(yīng)面法改進了迭代步長的取值，從而提高了收斂精度。同時，若設(shè)定可靠度指標(biāo)的迭代精度為0.001，則一般迭代2～3次即可收斂，計算時間能夠滿足實際工程的要求。從表2可以看出，采用改進響應(yīng)面法可以較準(zhǔn)確地計算懸索橋在活荷載作用下的靜力可靠指標(biāo)，而基于響應(yīng)面的蒙特卡羅法則計算誤差較大。另外，懸索橋的幾何非線性對扁平鋼箱梁的靜力可靠度會產(chǎn)生較大的影響。注：RSM-MCIS法結(jié)果為抽樣10萬次分析得到。

表2 鋼箱梁跨中豎向位移可靠度指標(biāo)

下面討論2種可靠度分析方法在原理上的區(qū)別。響應(yīng)面法的基本思想是通過一系列有限元數(shù)值計算來擬合一個響應(yīng)面以代替未知的、真實的極限狀態(tài)曲面。當(dāng)極限狀態(tài)曲面比較復(fù)雜的時候，這樣一個簡單的二次多項表達式是不能很好地描述一個復(fù)雜的多維曲面的，每一次迭代所得到的二次表達式僅僅只能反映迭代點附近的功能函數(shù)的情況，因此，要得到足夠精確的驗算點的坐標(biāo)，必然是一個迭代的過程。響應(yīng)面方法的實質(zhì)是通過一系列的迭代得到具有一定精度的可靠指標(biāo)，并不是為了得到結(jié)構(gòu)的極限狀態(tài)方程的實際表達式，實際上也是不可能的。本文所提出的改進響應(yīng)面法并沒有改變響應(yīng)面法的基本思路，但通過對迭代步長f的改進，使得響應(yīng)面法不收斂和精度差的問題得到了很好的改善，所以，采用改進響應(yīng)面法可以適用于類似懸索橋這種大型復(fù)雜結(jié)構(gòu)隱式極限狀態(tài)方程的可靠度分析。

基于響應(yīng)面上的蒙特卡羅法(RSM-MCIS)對采用響應(yīng)面法擬合出的近似極限狀態(tài)方程進行抽樣，當(dāng)實際的極限狀態(tài)方程比較簡單時，特別是當(dāng)實際極限狀態(tài)方程不含交叉項的二次多項式，RSM-MCIS法是精確且有效的。但是，當(dāng)近似極限狀態(tài)方程與真實極限狀態(tài)方程誤差較大時，RSM-MCIS法必然會產(chǎn)生較大誤差。如表2所示，RSM-MCIS法的分析結(jié)果甚至出現(xiàn)了負可靠指標(biāo)的情況，沒有實際的物理意義。因此，對于復(fù)雜的懸索橋結(jié)構(gòu)，采用RSM-MCIS法分析得到的可靠度結(jié)果是不合理的。

4 結(jié)語

本文以某大跨度懸索橋的扁平鋼箱梁結(jié)構(gòu)為分析對象，采用改進的響應(yīng)面法開展了扁平鋼箱梁的可靠度分析方法研究。分析結(jié)果表明：（1）采用改進響應(yīng)面法能有效地解決傳統(tǒng)響應(yīng)面法存在的收斂失敗和計算誤差大的缺點，適用于大型復(fù)雜結(jié)構(gòu)隱式極限狀態(tài)方程的可靠度分析；（2）采用改進響應(yīng)面法可以較準(zhǔn)確地計算懸索橋在活荷載作用下的靜力可靠指標(biāo)，而基于響應(yīng)面的蒙特卡羅法則計算誤差較大；（3）懸索橋的幾何非線性對扁平鋼箱梁的靜力可靠度會產(chǎn)生較大的影響。