基于概率神經(jīng)網(wǎng)絡(luò)的鋼桁梁橋螺栓脫落損傷定位識別研究

李保龍 房志恒 戰(zhàn)家旺

1.國能朔黃鐵路發(fā)展有限責(zé)任公司，河北肅寧062350；2.北京交通大學(xué)土木建筑工程學(xué)院，北京100044

鋼桁梁橋是鐵路橋梁的常見橋型之一。目前鐵路鋼桁梁橋多使用摩擦型高強螺栓進行連接，螺栓在使用過程中會產(chǎn)生銹蝕、松動、斷裂、脫落等病害。螺栓脫落會導(dǎo)致連接處傳力性能下降，一旦螺栓大量失效將給鋼桁梁橋帶來極大安全隱患，因此對螺栓脫落損傷進行識別具有重要意義。

目前螺栓病害檢測方法有目視檢測法、敲擊回聲法、壓電阻抗法、聲發(fā)射法、超聲波檢測法和基于動力測試的檢測方法［1-4］?，F(xiàn)有的螺栓病害檢測方法大多是局部檢測，而鋼桁梁橋的節(jié)點數(shù)量多，局部檢測難以有效識別螺栓脫落損傷。為減少待檢測節(jié)點的數(shù)量，對螺栓脫落損傷進行定位識別有很必要。因此，提出一種基于概率神經(jīng)網(wǎng)絡(luò)的鋼桁梁橋螺栓脫落損傷定位識別方法。該方法通過將螺栓脫落損傷定位在較小的子結(jié)構(gòu)內(nèi)，縮小損傷檢測范圍，提高檢測效率，為進一步定量評估損傷奠定基礎(chǔ)。

1 工程概況

以一單線下承式簡支鋼桁梁橋（圖1）為工程依托。鋼桁梁橋跨徑64 m，主桁高11 m，由2片主桁組成，主桁中心距5.75 m。該橋采用摩擦型高強螺栓連接，螺紋規(guī)格為M22，預(yù)拉力為200 kN，摩擦因數(shù)為0.45。鋼桁梁橋斷面及節(jié)點編號如圖2所示。

圖1 下承式簡支鋼桁梁橋全貌

圖2 鋼桁梁橋斷面及節(jié)點編號

2 基于概率神經(jīng)網(wǎng)絡(luò)的鋼桁梁橋螺栓脫落損傷定位識別方法

概率神經(jīng)網(wǎng)絡(luò)（Probabilistic Neural Network，PNN）是基于徑向基函數(shù)和經(jīng)典概率密度估計原理的一種前饋型神經(jīng)網(wǎng)絡(luò)。概率神經(jīng)網(wǎng)絡(luò)多用于損傷的定位，根據(jù)貝葉斯最優(yōu)分類決策理論對發(fā)生在區(qū)域內(nèi)的損傷進行模式分類，將損傷的定位識別轉(zhuǎn)化為損傷模式分類問題。使用概率神經(jīng)網(wǎng)絡(luò)進行鋼桁梁橋螺栓脫落損傷定位的流程如圖3所示。

圖3 螺栓脫落損傷定位流程

螺栓脫落損傷定位的具體步驟：①將鋼桁梁橋劃分為若干個子結(jié)構(gòu)，將發(fā)生在各子結(jié)構(gòu)內(nèi)的損傷視為一種損傷模式；②通過參數(shù)敏感性分析確定概率神經(jīng)網(wǎng)絡(luò)的輸入?yún)?shù)；③通過有限元模擬計算得到概率神經(jīng)網(wǎng)絡(luò)的樣本；④選擇部分樣本作為訓(xùn)練樣本，對概率神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練；⑤使用測試樣本對訓(xùn)練好的概率神經(jīng)網(wǎng)絡(luò)進行測試，以檢驗螺栓脫落損傷定位識別的正確率。

3 考慮螺栓脫落損傷的鋼桁梁橋有限元模型

使用ANSYS建立鋼桁梁橋的有限元模型（圖4），單元類型選用空間梁單元beam188。鋼桁梁橋的各個節(jié)點按剛接處理。

為模擬螺栓脫落損傷，將模型中每個主桁桿件劃分為3個單元，靠近節(jié)點的2個單元稱為桿端單元。通過降低桿端單元的抗彎剛度模擬螺栓脫落損傷［5］。由于截面慣性矩I的改變較為復(fù)雜，因此假定I不變，僅通過改變彈性模量E實現(xiàn)截面抗彎剛度EI的變化。

為簡化分析，假定主桁桿件在兩端發(fā)生相同程度的螺栓脫落損傷。損傷程度通過桿端抗彎剛度折減系數(shù)α來表征。α的計算公式［6］為

式中：Eu、Ed分別為損傷前、后桿端單元的彈性模量。

為了通過概率神經(jīng)網(wǎng)絡(luò)對螺栓脫落損傷進行定位識別，須將鋼桁梁橋劃分為若干個子結(jié)構(gòu)。子結(jié)構(gòu)法的產(chǎn)生最初是為了解決艦船、飛行器等復(fù)雜系統(tǒng)的靜動力分析問題，通過將自由度過多的復(fù)雜結(jié)構(gòu)拆分為簡單的子結(jié)構(gòu)，可以減少分析對計算能力的需求。損傷識別中子結(jié)構(gòu)法主要用于減少損傷識別的待識別參數(shù)數(shù)量，實現(xiàn)損傷的初步定位。子結(jié)構(gòu)的劃分方式有按照幾何形態(tài)劃分和按照幾何位置劃分。按照幾何位置劃分是將幾何位置相鄰的單元劃分到同一個子結(jié)構(gòu)中，本文采用該劃分方式。

將每片主桁按照幾何位置劃分為2個對稱的子結(jié)構(gòu)，2片主桁共劃分為4個子結(jié)構(gòu)，如圖5所示。后續(xù)訓(xùn)練概率神經(jīng)網(wǎng)絡(luò)時輸入的參數(shù)均由本節(jié)考慮螺栓脫落損傷的有限元模型計算得到。

圖5 鋼桁梁橋子結(jié)構(gòu)劃分

4 概率神經(jīng)網(wǎng)絡(luò)輸入?yún)?shù)確定

神經(jīng)網(wǎng)絡(luò)常見的輸入?yún)?shù)有標(biāo)準(zhǔn)化的頻率變化率（Normalized Frequency Change Rate，NFCR）、標(biāo)準(zhǔn)化的頻率變化比（Normalized Frequency Change Ratio，NFSR）和歸一化的損傷信號指標(biāo)（Normalized Damage Signal Index，NDSI）。

NFCR和NFSR是與頻率相關(guān)的輸入?yún)?shù)?；陬l率構(gòu)造的輸入?yún)?shù)無法分辨對稱位置發(fā)生的損傷，因此對對稱結(jié)構(gòu)的損傷識別有局限性。NDSI將振型考慮進來，解決了對稱結(jié)構(gòu)的損傷識別問題，但振型的變化對局部損傷并不敏感。

目前已經(jīng)有學(xué)者［7-9］將曲率模態(tài)與神經(jīng)網(wǎng)絡(luò)結(jié)合起來進行損傷識別。曲率模態(tài)對局部損傷敏感，因此本文選擇曲率模態(tài)變化率作為神經(jīng)網(wǎng)絡(luò)的輸入?yún)?shù)。

梁式結(jié)構(gòu)曲率K的計算公式［10］為

式中：ρ為曲率半徑；M為截面彎矩；x為梁長度方向距梁端的距離；y為梁豎向撓度，y'、y″分別為y對x的一階導(dǎo)數(shù)和二階導(dǎo)數(shù)。

考慮到梁的變形為小變形，y'可近似認為是0，則

因為無法對曲率模態(tài)進行測量，對位移模態(tài)進行中心差分計算得到曲率模態(tài)K?i。如果各差分點在x軸上等間距布置，則K?i可表示為

式中：?i為第i個節(jié)點處的模態(tài)位移；h為傳感器間距。

結(jié)構(gòu)損傷會導(dǎo)致曲率模態(tài)在損傷處發(fā)生突變。為反映損傷前后曲率模態(tài)的變化，定義第i階的豎彎曲率模態(tài)變化率KCi為

式中：K?iu、K?id分別為損傷前后第i階曲率模態(tài)。

輸入?yún)?shù)必須將發(fā)生損傷的子結(jié)構(gòu)有效區(qū)分出來，因此先對豎彎曲率模態(tài)變化率KC進行參數(shù)敏感性分析。由于4個子結(jié)構(gòu)具有對稱性，為簡化分析只在子結(jié)構(gòu)1中引入損傷。分析工況見表1。

表1 子結(jié)構(gòu)1參數(shù)KC敏感性分析工況

每個子結(jié)構(gòu)只選取下弦上的節(jié)點計算一階豎彎曲率模態(tài)變化率KC1，對各子結(jié)構(gòu)的KC1平均值進行對比，結(jié)果見圖6?？梢钥闯觯咏Y(jié)構(gòu)1內(nèi)的KC1平均值明顯高于無損傷的子結(jié)構(gòu)（子結(jié)構(gòu)2，3，4），說明曲率模態(tài)變化率可以將損傷子結(jié)構(gòu)有效區(qū)分出來。

圖6 各子結(jié)構(gòu)內(nèi)下弦節(jié)點的KC1平均值

選擇各子結(jié)構(gòu)內(nèi)的下弦節(jié)點（E1、E2、E3、E5、E6、E7、E1'、E2'、E3'、E5'、E6'、E7'）計算一階豎彎曲率模態(tài)變化率KC1，得到12維的神經(jīng)網(wǎng)絡(luò)輸入?yún)?shù)向量。

5 概率神經(jīng)網(wǎng)絡(luò)的訓(xùn)練與識別

使用不同訓(xùn)練樣本進行訓(xùn)練，分析單一子結(jié)構(gòu)發(fā)生單損傷和多損傷時概率神經(jīng)網(wǎng)絡(luò)的定位識別結(jié)果。神經(jīng)網(wǎng)絡(luò)的輸入?yún)?shù)為損傷前后鋼桁梁橋下弦節(jié)點的KC1。曲率模態(tài)采用第3節(jié)有限元模型經(jīng)模態(tài)分析得到。

在對概率神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練時，設(shè)定神經(jīng)網(wǎng)絡(luò)的輸出結(jié)果為4維列向量。輸出的［1，0，0，0］、［0，1，0，0］、［0，0，1，0］、［0，0，0，1］分別代表子結(jié)構(gòu)1、子結(jié)構(gòu)2、子結(jié)構(gòu)3、子結(jié)構(gòu)4發(fā)生損傷。完成訓(xùn)練后，將除訓(xùn)練樣本外的樣本設(shè)為測試樣本，輸入到訓(xùn)練好的概率神經(jīng)網(wǎng)絡(luò)中。通過神經(jīng)網(wǎng)絡(luò)輸出的向量，對發(fā)生損傷的子結(jié)構(gòu)進行識別。

5.1 概率神經(jīng)網(wǎng)絡(luò)樣本設(shè)置

為檢驗神經(jīng)網(wǎng)絡(luò)的泛化能力，設(shè)置α為0.3、0.4、0.5、0.6的4種損傷程度的樣本。根據(jù)同一子結(jié)構(gòu)內(nèi)發(fā)生螺栓脫落的桿件數(shù)量，設(shè)置單損傷、雙損傷及三損傷3種損傷樣本集。本文研究的鋼桁梁橋共58個主桁桿件，單損傷、雙損傷、三損傷工況每種損傷程度分別設(shè)置58、392、1 638個損傷樣本集。采用單損傷樣本對概率神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練，采用單損傷和多損傷樣本對訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)進行測試。工況設(shè)置見表2。

表2 單損傷訓(xùn)練樣本工況設(shè)置

5.2 定位識別結(jié)果

1）單損傷定位識別結(jié)果

按照表2中的單損傷訓(xùn)練樣本對概率神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練，測試樣本為除訓(xùn)練樣本外的單損傷樣本。

概率神經(jīng)網(wǎng)絡(luò)須給出Spread值。Spread值對概率神經(jīng)網(wǎng)絡(luò)的模式分類效果影響很大，可通過試算確定最佳的Spread值。Spread在0.1～1.0取值，對各訓(xùn)練樣本工況進行試算。不同Spread值下單損傷定位識別正確率見表3。

表3 不同Spread值下單損傷定位識別正確率

由表3可見，不同Spread值下單損傷定位識別正確率存在較大差異。以訓(xùn)練樣本工況1為例，Spread值取0.1和1.0時定位識別正確率相差20.69%，可見Spread值對識別結(jié)果有很大影響。

選擇合適的（定位識別正確率最高）Spread值后各訓(xùn)練樣本工況下單損傷定位識別正確率見圖7。

圖7 各訓(xùn)練樣本工況下單損傷定位識別正確率

由表2和圖7可見：①訓(xùn)練樣本數(shù)量不足會導(dǎo)致識別正確率低。工況1的訓(xùn)練樣本數(shù)量僅為8個，識別正確率明顯低于其他工況。②使用同一損傷程度的單損傷樣本進行訓(xùn)練時，樣本數(shù)量達到16個時識別正確率已經(jīng)很高，繼續(xù)增加訓(xùn)練樣本數(shù)量對識別正確率的提升效果很小。③訓(xùn)練樣本數(shù)量相同時其損傷程度（α為0.3、0.4、0.5、0.6）對單損傷樣本的識別結(jié)果沒有明顯影響。

2）多損傷定位結(jié)果

仍按照表2中的單損傷訓(xùn)練樣本對概率神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練。從4種損傷程度（α為0.3、0.4、0.5、0.6）的多損傷樣本集中各選擇160個作為測試樣本。對不同的訓(xùn)練樣本工況進行試算，過程與上述單損傷定位相同。各訓(xùn)練樣本工況下選取合適的Spread值后多損傷定位識別正確率見圖8。

由表2和圖8可見：①訓(xùn)練樣本的損傷程度對多損傷定位識別結(jié)果影響很??；②使用足夠多的單損傷樣本對神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練后，對多損傷定位識別的正確率可高達90%以上。

圖8 各訓(xùn)練樣本工況下多損傷定位識別正確率

6 結(jié)論

1）概率神經(jīng)網(wǎng)絡(luò)的Spread值對識別結(jié)果影響很大，須通過試算合理取值。

2）訓(xùn)練樣本數(shù)量相同時其損傷程度對概率神經(jīng)網(wǎng)絡(luò)識別結(jié)果的影響很小。

3）使用足夠多的單損傷樣本對概率神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練后，對單損傷及多損傷定位識別的正確率均高達90%以上。