基于PSO-BP神經(jīng)網(wǎng)絡(luò)的鋼管混凝土拱橋主梁地震易損性分析
2023-04-29 00:00:00龍梅溫干祥
西部交通科技 2023年9期
作者簡介：龍 梅（1974—），高級工程師，講師，主要從事道路橋梁教學、設(shè)計、施工管理等工作。
摘要：為解決大跨度鋼管混凝土拱橋在地震易損性評估中計算效率低、計算成本高的問題，文章提出了一種基于粒子群算法改進BP神經(jīng)網(wǎng)絡(luò)的地震易損性計算方法?；谡龖B(tài)分布函數(shù)建立了鋼管混凝土拱橋在地震激勵下的易損性條件概率表征方式，采用粒子群算法對BP神經(jīng)網(wǎng)絡(luò)的權(quán)值與閾值進行尋優(yōu)，訓(xùn)練了鋼管混凝土拱橋在地震激勵下的響應(yīng)模型。同時，以某大跨度鋼管混凝土拱橋為例，分析了其在考慮不同損傷狀態(tài)下的地震易損性曲線。結(jié)果表明：PSO-BP神經(jīng)網(wǎng)絡(luò)模型可以較為精確地擬合鋼管混凝土拱橋的易損性曲線，預(yù)測精度遠高于未經(jīng)優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)；主梁在輕微、中度和嚴重損傷狀態(tài)下的地震易損性曲線趨勢大致相同，損傷概率增速先增加后減小，當?shù)孛娣逯导铀俣龋?.4 g時，主梁完全破壞概率幾乎為0。
關(guān)鍵詞：鋼管混凝土拱橋；抗震性能；易損性曲線；粒子群算法；BP神經(jīng)網(wǎng)絡(luò)
中圖分類號：U448.22
0 引言
大跨度橋梁在地震下的結(jié)構(gòu)安全一直以來都是專家學者的研究熱點，鋼管混凝土拱橋作為一種典型的大跨度橋梁，其在地震作用下的響應(yīng)規(guī)律和損傷概率備受關(guān)注，眾多專家學者也針對大跨度鋼管混凝土拱橋的地震易損性展開了相關(guān)分析。黃飛鴻等以鋼管混凝土拱橋支座、拱肋和系梁為對象，采用增量動力分析法研究了其在4種損傷狀態(tài)下的易損性曲線［1］；熊程充分考慮了沖刷作用對于鋼管混凝土拱橋結(jié)構(gòu)的影響，分析了沖刷作用下不同地震激勵下的橋梁損傷概率，建立了考慮沖刷深度的構(gòu)建概率地震需求模型［2］；王力等針對異型鋼管混凝土拱橋結(jié)構(gòu)，基于Midas Civil軟件建立了非線性有限元模型，在增量動力分析的基礎(chǔ)上評估了橋梁主要構(gòu)件的損傷情況，提出了相關(guān)建議［3］；馮莉等提出了基于性能的高鐵鋼管混凝土拱橋地震易損性評估方法，并通過有限元模型增量動力分析得到了橋梁在不同地震激勵下的經(jīng)濟風險［4］；劉震等為研究近斷層地震動對鋼管混凝土拱橋結(jié)構(gòu)損傷的影響，在Park-Ang損傷模型的基礎(chǔ)上繪制了鋼管混凝土拱橋的地震易損曲線，分析了近斷層效應(yīng)對結(jié)構(gòu)地震易損性的影響［5］。
從相關(guān)研究可知，對于大跨度鋼管混凝土這類復(fù)雜結(jié)構(gòu)在地震作用下的響應(yīng)分析，傳統(tǒng)的有限元計算方法存在計算效率低、計算成本高等不足［6］，為解決這一問題，已有專家學者開始采用機器學習算法代替原有的有限元模型進行計算分析。楊耀鑫等基于神經(jīng)網(wǎng)絡(luò)算法建立了結(jié)構(gòu)在地震后的框架結(jié)構(gòu)損傷評估模型［7］；許澤坤等采用長短期記憶神經(jīng)網(wǎng)絡(luò)模型計算了非線性框架結(jié)構(gòu)的地震響應(yīng)［8］。從現(xiàn)有的研究中可以發(fā)現(xiàn)，對于鋼管混凝土拱橋地震易損性的評估仍以有限元計算方式為主，基于此，本文提出一種基于粒子群算法改進BP神經(jīng)網(wǎng)絡(luò)的鋼管混凝土拱橋主梁地震易損性計算方法，以期為類似結(jié)構(gòu)的地震易損性分析提供一種新的思路。
1 地震易損性分析原理
1.1 橋梁地震易損性原理
橋梁結(jié)構(gòu)的地震易損性的本質(zhì)是一種損傷破壞概率，指橋梁結(jié)構(gòu)在地震激勵的作用下超越結(jié)構(gòu)極限狀態(tài)而發(fā)生破壞的概率，令其為Pf，則對于該條件概率的描述，一般采用式（1）進行描述：
Pf=P（SD-SC≥0|IM）"" （1）
式中：P（x）——條件概率；
SD——地震激勵下的結(jié)構(gòu)響應(yīng)峰值；
SC——結(jié)構(gòu)的極限承載能力；
IM——地震動強度。
地震激勵下的結(jié)構(gòu)響應(yīng)峰值SD與地震動強度IM的對數(shù)之間存在如下線性相關(guān)關(guān)系：
lnSD=mln（IM）+n""" （2）
式中：m、n——線性相關(guān)系數(shù)。
地震激勵下的橋梁結(jié)構(gòu)響應(yīng)峰值與結(jié)構(gòu)承載能力分別服從對數(shù)正態(tài)分布，如式（3）所示：
式中：μD、μC——結(jié)構(gòu)地震響應(yīng)峰值和構(gòu)建承載力的數(shù)學期望；
βD、βC——結(jié)構(gòu)地震響應(yīng)峰值和構(gòu)建承載力的方差。
根據(jù)對數(shù)運算原則，聯(lián)立式（1），可將橋梁結(jié)構(gòu)地震易損性的條件概率描述為式（4）所示：
Pf=P（lnSD-lnSC≥0|IM）"" （4）
由此，可將橋梁地震易損性條件概率表示為：
式中：?（x）——標準正態(tài)分布函數(shù)；
b——線性相關(guān)系數(shù)。
選取地面峰值加速度PGA為地震動強度IM的表征指標，取β2Cl+β2D值為0.5。
1.2 主梁損傷指標及地震動的選取
合理地定義結(jié)構(gòu)損傷指標是橋梁結(jié)構(gòu)地震易損性分析的關(guān)鍵，通過定義損傷指標，可以在描述橋梁損傷狀態(tài)的前提下計算結(jié)構(gòu)的易損性曲線。對于大跨度鋼管混凝土拱橋的主梁，存在的損傷等級主要分為輕微、中等、嚴重和破壞，參考文獻［9］中對于鋼管混凝土拱橋主梁損傷指標的評定，基于鋼筋曲率定義大跨度鋼管混凝土拱橋的損傷指標如表1所示。
以地面峰值加速度PGA作為篩選指標，地震加速度選取范圍為0.1～2.2 g，從太平洋地震研究中心篩選50條地震波作為結(jié)構(gòu)地震動激勵，由于篩選的地震加速度范圍較大，可以在一定程度上提升大跨度鋼管混凝土拱橋的主梁易損性曲線精度。
2 基于PSO-BP神經(jīng)網(wǎng)絡(luò)的地震響應(yīng)預(yù)測模型
2.1 BP神經(jīng)網(wǎng)絡(luò)原理
采用數(shù)值計算的方式進行結(jié)構(gòu)抗震分析往往存在計算量大，計算時間成本高，計算效率低等不足。為提高結(jié)構(gòu)地震響應(yīng)的計算效率，引入BP神經(jīng)網(wǎng)絡(luò)對結(jié)構(gòu)在地震激勵下的響應(yīng)進行預(yù)測。BP神經(jīng)網(wǎng)絡(luò)是一種由三層神經(jīng)元結(jié)構(gòu)組成的人工預(yù)測網(wǎng)絡(luò)，可以通過隱含層神經(jīng)元之間的信息傳遞學習樣本數(shù)據(jù)與目標數(shù)據(jù)之間的映射關(guān)系，其基本網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示。
從圖1中可以看出，輸入層神經(jīng)元為xn，輸出層神經(jīng)元為ym，定義隱含層神經(jīng)元為bu，則由n個輸入層神經(jīng)元xn到m個輸出層神經(jīng)元ym的傳遞關(guān)系如式（6）所示：
式中：ωnu、vum——輸入層神經(jīng)元n與隱含層神經(jīng)元u，隱含層神經(jīng)元u與輸出層神經(jīng)元m的連接權(quán)重；
ku、pm——隱含層與輸出層的神經(jīng)元閾值。
BP神經(jīng)網(wǎng)絡(luò)的預(yù)測原理為通過將每次訓(xùn)練的實際輸出值與期望輸出值進行比較，對預(yù)測誤差進行修正，并對該修正過程進行反向傳播，傳遞給每一層神經(jīng)元，實現(xiàn)誤差收斂的效果。
基于PSO-BP神經(jīng)網(wǎng)絡(luò)的鋼管混凝土拱橋主梁地震易損性分析/龍 梅，溫干祥
2.2 基于PSO-BP神經(jīng)網(wǎng)絡(luò)的地震響應(yīng)預(yù)測模型
僅依靠BP神經(jīng)網(wǎng)絡(luò)自身對于神經(jīng)元權(quán)值、閾值的搜索往往修正效率較低，故本文引入粒子群算法對訓(xùn)練過程中的最佳權(quán)值、閾值進行尋優(yōu)求解［10］。將BP神經(jīng)網(wǎng)絡(luò)的權(quán)值閾值轉(zhuǎn)化為粒子群在搜索空間中的位置坐標，粒子群在d維搜索空間中的位置與速度更新策略如式（7）、式（8）所示：
基于PSO-BP神經(jīng)網(wǎng)絡(luò)的地震響應(yīng)預(yù)測模型訓(xùn)練過程如圖2所示。
PSO算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的主要流程如下：
步驟1：通過有限元模型形成基礎(chǔ)的地震激勵與主梁響應(yīng)的訓(xùn)練數(shù)據(jù)，形成歸一化的樣本數(shù)據(jù)集。
步驟2：將BP神經(jīng)網(wǎng)絡(luò)權(quán)值閾值轉(zhuǎn)化為粒子群算法的位置坐標，初始化粒子群算法和BP神經(jīng)網(wǎng)絡(luò)基本參數(shù)，對樣本數(shù)據(jù)集進行訓(xùn)練擬合。
步驟3：判斷BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練精度，若達到目標精度則輸出最優(yōu)網(wǎng)絡(luò)權(quán)值閾值；若未達到則更新粒子群速度與位置，尋找新一輪迭代的全局最優(yōu)值。
步驟4：判斷粒子群算法是否達到最大迭代次數(shù)，若達到則輸出最優(yōu)參數(shù)組合返回BP神經(jīng)網(wǎng)絡(luò)重新訓(xùn)練擬合，若未達到則繼續(xù)迭代尋優(yōu)。
通過上述流程，可以高效地實現(xiàn)BP神經(jīng)網(wǎng)絡(luò)對于地震激勵與結(jié)構(gòu)響應(yīng)之間的映射關(guān)系學習，形成大跨度拱橋主梁地震易損性曲線的計算模型。
3 工程算例
3.1 工程概況及有限元模型
某大跨度上承式鋼管混凝土拱橋主跨計算跨徑為450 m，主拱圈采用等寬變高度空間桁架結(jié)構(gòu)，拱軸線為懸鏈線，拱軸線系數(shù)m=1.55，矢高100 m，矢跨比f=1/4.5，鋼管內(nèi)灌注C60自密實微膨脹混凝土。主橋橋面系采用跨徑31.6 m的鋼混疊合梁，梁高2.1 m，橫向間距6.5 m，順橋向每3.95 m設(shè)置一道鋼橫梁，采用焊接工字型截面，端部設(shè)置混凝土橫梁。主橋橋型布置如圖3所示。
采用Midas Civil軟件建立大跨度鋼管混凝土拱橋的有限元計算模型，主梁、主拱圈和立柱均采用梁單元，兩側(cè)拱腳用全固結(jié)約束，全橋共4 360個節(jié)點，6 888個單元。有限元模型如圖4所示。
3.2 PSO-BP神經(jīng)網(wǎng)絡(luò)擬合結(jié)果
在MATLAB軟件中編寫粒子群算法的執(zhí)行程序，并采用MATLAB中BP神經(jīng)網(wǎng)絡(luò)工具箱對主梁在地震激勵下的結(jié)構(gòu)響應(yīng)進行學習。如圖5所示給出了粒子群算法對BP神經(jīng)網(wǎng)絡(luò)權(quán)值與閾值的尋優(yōu)過程。從圖5可以看出，粒子群算法的收斂性良好，在第50次迭代中基本達到全局最優(yōu)值。
如圖6所示給出了10組地震波激勵下的擬合結(jié)果。從圖6可以看出，采用粒子群算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)對結(jié)構(gòu)響應(yīng)的預(yù)測效果明顯優(yōu)于未經(jīng)優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)，10個測試集的預(yù)測結(jié)果表明，粒子群算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)平均相對誤差僅為5.2%，最大相對誤差僅為9.8%，而未經(jīng)粒子群算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)平均相對誤差達到了17.6%，最大相對誤差為23.8%，說明粒子群算法對于BP神經(jīng)網(wǎng)絡(luò)權(quán)值和閾值的優(yōu)化效果十分明顯，該模型可作為結(jié)構(gòu)易損性曲線的預(yù)測模型。
3.3 地震易損性分析
如圖7所示給出了不同損傷等級下的大跨度鋼管混凝土拱橋主梁地震易損性曲線。從圖7可以看出，不同損傷狀態(tài)下PSO-BP神經(jīng)網(wǎng)絡(luò)對于地震易損性曲線的擬合效果較好，與有限元模型計算得到的易損性曲線趨勢基本吻合，在地面峰值加速度不斷增加的情況下主梁損傷概率也隨之增大，說明PSO-BP神經(jīng)網(wǎng)絡(luò)模型可以良好地應(yīng)用于大跨度鋼管混凝土拱橋地震易損性曲線的預(yù)測中。而未經(jīng)粒子群算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)雖然可以大致擬合出主梁易損性的變化趨勢，但其損傷超越概率的預(yù)測結(jié)果與有限元模型的計算結(jié)果相差較大，數(shù)據(jù)存在一定的失真，不具備可參考性。
對比圖7（a）～（c）可以看出，輕微、中度和嚴重損傷下大跨度鋼管混凝土拱橋的主梁損傷概率發(fā)展趨勢大致相同，超越概率的增速體現(xiàn)為先增加后減小的趨勢，當?shù)孛娣逯导铀俣忍幱?～0.2 g時，主梁損傷概率較低，且增速較慢；當?shù)孛娣逯导铀俣忍幱?.2～1.0 g時，主梁損傷概率與地面峰值加速度間大致呈線性相關(guān)關(guān)系；當?shù)孛娣逯导铀俣忍幱?.0～1.2 g時，主梁損傷概率增速放緩，輕微損傷情況下的增速減緩最為明顯。當PGA取1.2 g時，輕微、中度、嚴重損傷概率分別為0.92、0.87和0.86。由此可知當?shù)孛娣逯导铀俣容^大時，結(jié)構(gòu)發(fā)生嚴重損傷的概率較大。從圖7（d）中可以看出，考慮在大跨度鋼管混凝土拱橋完全破壞的情況下，當?shù)孛娣逯导铀俣龋?.4 g時，主梁出現(xiàn)完全破壞的概率幾乎為0，當?shù)孛娣逯导铀俣忍幱?.4～1.2 g時，主梁損傷概率隨PGA的增加呈現(xiàn)出線性增加的趨勢，地面峰值加速度取1.2 g情況下，主梁出現(xiàn)完全破壞的概率為0.43，說明主梁存在一定的失效風險。
4 結(jié)語
本文基于粒子群算法對BP神經(jīng)網(wǎng)絡(luò)進行優(yōu)化，提出了一種大跨度鋼管混凝土拱橋主梁的地震易損性計算方法，對比了不同損傷情況下的主梁地震損傷概率，得到結(jié)論如下：
（1）采用粒子群算法改進BP神經(jīng)網(wǎng)絡(luò)的方法可以有效提升BP神經(jīng)網(wǎng)絡(luò)的預(yù)測精度，降低BP神經(jīng)網(wǎng)絡(luò)對于大跨度鋼管混凝土拱橋結(jié)構(gòu)響應(yīng)的預(yù)測誤差，提高BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練速度。
（2）從不同損傷情況的主梁地震易損性曲線可以看出，未經(jīng)優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)可以大致擬合出主梁在地震激勵下的損傷概率發(fā)展趨勢，但其結(jié)果相對有限元計算結(jié)果存在較大誤差，而經(jīng)過粒子群算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)對主梁地震易損性曲線的預(yù)測效果良好，與有限元模型的計算結(jié)果基本保持一致。
（3）當考慮主梁輕微、中度和嚴重損傷狀態(tài)時，主梁的損傷概率趨勢大致相同，損傷概率增速先增大后減小，當考慮主梁完全破壞情況時，地面峰值加速度＜0.4 g時，主梁破壞概率幾乎為0，地面峰值加速度達到1.2 g時，主梁完全破壞的概率為0.43，結(jié)構(gòu)存在一定失效風險。
參考文獻
［1］黃飛鴻，何沛祥，吳騰飛.下承式鋼管混凝土拱橋地震易損性分析［J］.合肥工業(yè)大學學報（自然科學版），2022，45（6）：801-807.
［2］熊 程.考慮沖刷作用的大跨鋼管混凝土拱橋地震易損性評估［D］.南昌：南昌大學，2021.
［3］王 力，虞廬松，劉世忠，等.橫向地震作用下異型鋼管混凝土拱橋地震易損性分析［J］.蘭州交通大學學報，2020，39（5）：13-19.
［4］馮 莉，樊燕燕，王 力，等.基于性能的高速鐵路鋼管混凝土拱橋地震經(jīng)濟風險分析［J］.鐵道科學與工程學報，2019，16（3）：573-580.
［5］劉 震，張 哲，王英博.基于近斷層效應(yīng)鋼管混凝土拱橋地震易損性分析［J］.武漢理工大學學報（交通科學與工程版），2016，40（5）：815-819.
［6］于京池，金愛云，潘堅文，等.基于GA-BP神經(jīng)網(wǎng)絡(luò)的拱壩地震易損性分析［J］.清華大學學報（自然科學版），2022，62（8）：1 321-1 329.
［7］楊耀鑫，楊永強，楊 游，等.基于神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)震后快速損傷評估［J］.世界地震工程，2023，39（1）：49-58.
［8］許澤坤，陳 雋.非線性結(jié)構(gòu)地震響應(yīng)的神經(jīng)網(wǎng)絡(luò)算法［J］.工程力學，2021，38（9）：133-145.
［9］黃香健，陳雙慶，王 華，等.沖刷作用下鋼管混凝土拱橋地震易損性分析［J］.公路工程，2022，47（5）：40-46，112.
［10］曹檢云，劉國坤，孫劍峰.大跨度懸臂澆筑鋼筋混凝土拱橋多目標索力優(yōu)化研究［J］.公路工程，2022，47（3）：30-34，56.
收稿日期：2023-04-09
