基于有限元法和一次二階矩法的橋梁結(jié)構(gòu)可靠度研究

王虎軍 賀小衛(wèi) 高宇 扈軒誠

西安市政設(shè)計研究院有限公司 710068

引言

以概率論和數(shù)理統(tǒng)計為基礎(chǔ)的可靠度計算方法，如中心點法、驗算點法、JC 法、MC 法、響應(yīng)面法等可為橋梁結(jié)構(gòu)的可靠性評估提供便利。但由于實際橋梁結(jié)構(gòu)體系復(fù)雜，基于抗力—荷載效應(yīng)的結(jié)構(gòu)功能函數(shù)可能高度非線性，甚至不能顯式表達。這就使得以功能函數(shù)的解析表達式為基礎(chǔ)的中心點法、驗算點法、JC法和MC法計算結(jié)構(gòu)可靠度變得困難，甚至無法進行下去。響應(yīng)面法雖然能夠用一個簡單的函數(shù)擬合隱含、復(fù)雜的真實功能函數(shù)，但當(dāng)隨機變量較多時，目前采用的多項式形式的響應(yīng)面函數(shù)的系數(shù)計算工作量大且擬合誤差大而變得不適用。

本文結(jié)合確定性有限元法和可靠度計算中改進的一次二階矩法，在不必求得結(jié)構(gòu)功能函數(shù)明確表達式的情況下，進行結(jié)構(gòu)可靠指標(biāo)的分析計算，通過MATLAB編程，計算了一鋼管混凝土拱橋縱梁跨中截面正常使用極限狀態(tài)下的可靠度指標(biāo)，并對各隨機變量靈敏度進行分析，為結(jié)構(gòu)設(shè)計中提高可靠度給出了方向。

1 改進的一次二階矩法

改進的一次二階矩法是將功能函數(shù)的線性化展開點選在失效的邊界Z＝0 上，而且選在結(jié)構(gòu)最大可能失效概率的點上，故又稱設(shè)計驗算點法。其計算原理在文獻［11］中有詳細(xì)論述，在此不再累述。設(shè)結(jié)構(gòu)的功能函數(shù)為Z＝Z（X1，X2，X3，…Xn），其中Xi（i＝1，2，3…n）為隨機變量，如結(jié)構(gòu)幾何尺寸、材料物理參數(shù)及外荷載等，且各隨機變量間相互獨立，則改進的一次二階矩法的迭代計算步驟如下：

2 有限元法與改進的一次二階矩法的結(jié)合

基于抗力—荷載效應(yīng)的結(jié)構(gòu)功能函數(shù)中，荷載效應(yīng)通常為幾何尺寸、荷載、彈性模量、泊松比等隨機變量的隱式函數(shù)，即S＝S（X1，X2，X3，…Xn），其中Xi（i＝1，2，3…n）為各隨機變量。而抗力R則通常為規(guī)定的一限值或具有明確表達式，則結(jié)構(gòu)的功能函數(shù)可表示為Z＝R－S。設(shè)各隨機變量的分布類型已知，則可根據(jù)改進的一次二階矩法迭代求解可靠指標(biāo)和驗算點，由于功能函數(shù)中荷載效應(yīng)不能顯式表達，導(dǎo)致無法直接計算功能函數(shù)在驗算點處的偏導(dǎo)數(shù)，本文利用有限元法的數(shù)值模擬功能及數(shù)值分析中的向前差分法［2］很好地解決了偏導(dǎo)數(shù)的求解問題，即：

3 工程實例

3.1 工程概況及有限元建模

某中承式鋼管混凝土拱橋，計算跨徑為260m，計算矢跨比1／4，拱軸線為懸鏈線，拱軸系數(shù)m＝1.5。拱肋為四肢鋼管混凝土桁架結(jié)構(gòu)，上下弦桿為平放的啞鈴型斷面。拱肋由外徑950mm的Q345Dq螺旋焊管及其內(nèi)填充的C50 混凝土組成，截面總高5.35m，總寬2.75m；吊桿采用OVM.GJ15-15 拉索，橫向雙吊桿體系，縱向間距10m，橫向間距0.4m ＋18.5m ＋0.4m；橫梁全橋共設(shè)置29 道（19 道吊桿橫梁＋2 道拱上橫梁＋8 道立柱橫梁），吊桿橫梁為等高度預(yù)制“Ⅰ”字型梁，高1.8m；橋面系采用吊桿橫梁上設(shè)置縱向“T”型行車道梁，預(yù)制施工，梁高0.7m，腹板寬0.2m，翼緣寬1m，縱向長10m；橋面鋪裝采用10cm 厚聚酯纖維瀝青混凝土，設(shè)置1.5%的雙向橫坡和3%的雙向縱坡。

采用大型通用有限元分析軟件Midas Civil建立全橋模型。鋼管混凝土拱肋、拱上橫撐、橫梁、拱上立柱、行車道板均為梁單元；吊桿采用桁架單元；拱腳位置處用一般支撐中約束所有自由度模擬。全橋共離散為4654個節(jié)點和5460個單元。實橋有限元模型見圖2。

圖2 實橋有限元模型Fig.2 Finite element model of real bridge

3.2 結(jié)構(gòu)可靠度分析

采用規(guī)范規(guī)定的撓度限值作為鋼管混凝土拱橋正常使用極限狀態(tài)下的失效準(zhǔn)則，以鋼管混凝土拱橋縱梁跨中撓度超過限值uL為主要失效模式，根據(jù)《公路鋼管混凝土拱橋設(shè)計規(guī)范》（JTG／T D65—2015）計算得uL＝0.325m，建立結(jié)構(gòu)正常使用極限狀態(tài)功能函數(shù)Z為：

式中：A1為拱肋截面面積；A2為吊桿截面面積；A3為橫梁截面面積；A4為縱梁截面面積；E1為Q345Dq 鋼彈性模量；E2為C50 混凝土彈性模量；E3為拉索鋼絞線彈性模量；I1為拱肋截面慣性矩；I2為橫梁截面慣性矩；I3為縱梁截面慣性矩；F為活載。

各參數(shù)的分布類型、均值、變異系數(shù)主要參照實橋檢測結(jié)果，并結(jié)合《公路工程結(jié)構(gòu)可靠度設(shè)計統(tǒng)一標(biāo)準(zhǔn)》及相關(guān)文獻確定，見表1。

表1 各隨機變量的統(tǒng)計參數(shù)Tab.1 Statistical parameters of random variables

按圖1 迭代步驟進行迭代計算，共迭代16次滿足規(guī)定的誤差要求，有限元計算171 次。部分迭代結(jié)果見表2。

圖1 有限元與改進的一次二階矩法結(jié)合解算流程Fig.1 Flow chart of the combination of finite element and improved first order second moment method

表2 迭代計算過程Tab.2 Iterative calculation process

鋼管混凝土拱橋在設(shè)計過程中，材料性能和幾何參數(shù)等不確定性因素勢必對結(jié)構(gòu)可靠性產(chǎn)生 影響，為了準(zhǔn)確掌握各隨機變量對結(jié)構(gòu)可靠度指標(biāo)的影響程度，對鋼管混凝土拱橋各隨機變量進行靈敏性分析，給出了各隨機變量重要性程度的比較值，如表3 所示。

表3 各隨機變量的靈敏度系數(shù)Tab.3 Sensitivity coefficient of each random variable

為了更直觀地分析各隨機變量對結(jié)構(gòu)可靠指標(biāo)的影響程度，將各隨機變量的靈敏度值的絕對值進行比較，如圖3 所示。

圖3 各隨機變量靈敏度絕對值比較Fig.3 Comparison of sensitivity absolute value of each random variable

由表3 及圖3 可見，拱肋面積A1的靈敏度最大為0.925，對結(jié)構(gòu)可靠指標(biāo)的影響最大，縱梁截面面積A4的靈敏度絕對值最小為0.014，對可靠指標(biāo)影響較小，實際計算時可不考慮其影響。

4 結(jié)語

本文以實橋工程為依托，采用有限元法和改進的一次二階矩法相結(jié)合的方法對其可靠度進行研究，驗證了該方法的可行性，并通過各隨機變量對結(jié)構(gòu)可靠指標(biāo)的靈敏度分析，得到各隨機變量對結(jié)構(gòu)可靠性的影響程度，為結(jié)構(gòu)的優(yōu)化設(shè)計提供方向。