何定橋,楊 軍
(清華大學(xué) 土木工程安全與耐久教育部重點(diǎn)實(shí)驗(yàn)室,北京 100084)
近年來(lái),隨著計(jì)算機(jī)技術(shù)、傳感器技術(shù)的快速發(fā)展,結(jié)構(gòu)健康監(jiān)測(cè)逐漸成為研究熱點(diǎn),結(jié)構(gòu)健康監(jiān)測(cè)的一個(gè)重要目的是結(jié)構(gòu)損傷識(shí)別與定位[1-2]。結(jié)構(gòu)損傷識(shí)別與定位由以下部分構(gòu)成:(1) 判別結(jié)構(gòu)是否存在損傷;(2) 判別結(jié)構(gòu)損傷發(fā)生的位置;(3) 判定結(jié)構(gòu)的破壞程度;(4) 對(duì)結(jié)構(gòu)健康狀況進(jìn)行評(píng)估。其中,對(duì)結(jié)構(gòu)損傷位置及損傷程度的判斷是難點(diǎn)?;谏鲜?2)、(3)部分可以實(shí)現(xiàn)結(jié)構(gòu)狀況的整體評(píng)價(jià)。通過(guò)無(wú)線智能傳感器配合云平臺(tái)可以實(shí)現(xiàn)建筑結(jié)構(gòu)的全天候不間斷監(jiān)測(cè)。在地震荷載、風(fēng)荷載等作用下結(jié)構(gòu)出現(xiàn)損傷時(shí),監(jiān)測(cè)數(shù)據(jù)可以有效反映結(jié)構(gòu)的損傷信息[3-4]。
基于監(jiān)測(cè)數(shù)據(jù)的結(jié)構(gòu)損傷識(shí)別與定位可分為模態(tài)參數(shù)方法和非模態(tài)參數(shù)方法。模態(tài)參數(shù)法是根據(jù)結(jié)構(gòu)模態(tài)參數(shù)在損傷前后的變化反映結(jié)構(gòu)損傷的信息。常見(jiàn)的模態(tài)參數(shù)包括固有頻率[5]、模態(tài)振型[6]、振型曲率[7]及應(yīng)變曲率模態(tài)曲率[8]等?;谀B(tài)參數(shù)的方法因?qū)p傷不敏感、模態(tài)振型測(cè)量復(fù)雜、容易受外界干擾等問(wèn)題,在實(shí)際工程應(yīng)用中受到很大的限制,實(shí)用性還需進(jìn)一步的研究[9]。非模態(tài)參數(shù)法原理是對(duì)結(jié)構(gòu)振動(dòng)響應(yīng)進(jìn)行變換后得到可以反映結(jié)構(gòu)損傷的參數(shù)來(lái)進(jìn)行損傷識(shí)別與定位。非模態(tài)參數(shù)法包括基于小波變換[10-12]、希爾伯特黃變換[13]、時(shí)間序列等模型方法。時(shí)間序列模型最初應(yīng)用于經(jīng)濟(jì)學(xué)和電氣工程領(lǐng)域,后在結(jié)構(gòu)監(jiān)測(cè)領(lǐng)域被大量應(yīng)用。Fugate 等[14]建立了結(jié)構(gòu)的自回歸模型,提取模型殘差作為損傷特征量,采用質(zhì)量控制圖法監(jiān)測(cè)結(jié)構(gòu)狀態(tài)。H.Sohn等[15]建立了自回歸-外加輸入自回歸模型(AR-ARX),以模型殘差作為特征量,采用統(tǒng)計(jì)模式識(shí)別方法進(jìn)行損傷識(shí)別與定位。劉毅等[12]建立了自回歸滑動(dòng)平均模型(ARMA),基于自回歸部分參數(shù),采用主成分分析法提出結(jié)構(gòu)損傷特征指標(biāo)。然后采用t-檢驗(yàn)考察該指標(biāo)在損傷前后的變化進(jìn)行損傷定位。王真等[16]建立了自回歸系數(shù)的損傷靈敏度矩陣,通過(guò)該矩陣反映的自回歸系數(shù)變化與損傷系數(shù)變化之間的關(guān)系進(jìn)行損傷識(shí)別與定位。
盧宏彬[17]基于ARMA模型構(gòu)造了殘差指標(biāo)、前三階指標(biāo)和馬氏距離指標(biāo)在內(nèi)的損傷指標(biāo)體系對(duì)結(jié)構(gòu)損傷進(jìn)行了識(shí)別,并將該方法應(yīng)用于主跨428 m的廣州新光大橋健康監(jiān)測(cè)系統(tǒng)上。杜永峰等[18]建立了AR模型,將待識(shí)別工況的殘差與AR預(yù)測(cè)參考模型的殘差的方差之比作為損傷指標(biāo)。張凱瑋等[19]基于ARMA模型和馬氏距離定義了損傷指標(biāo),根據(jù)各工況下隧道結(jié)構(gòu)的加速度響應(yīng)進(jìn)行了損傷識(shí)別分析,實(shí)現(xiàn)對(duì)盾構(gòu)隧道結(jié)構(gòu)損傷的檢測(cè)和定位。
可見(jiàn),當(dāng)前基于ARMA模型的結(jié)構(gòu)損傷識(shí)別研究普遍利用自回歸系數(shù)的統(tǒng)計(jì)學(xué)特征,理論意義不明確,主觀性較強(qiáng)。而核嶺回歸算法具有非線性、擬合性能強(qiáng)、泛化能力強(qiáng)等特點(diǎn)。本文結(jié)合ARMA時(shí)間序列模型與核嶺回歸提出了一種結(jié)構(gòu)損傷識(shí)別與定位的新方法,并用算例驗(yàn)證了其有效性。
結(jié)構(gòu)損傷會(huì)改變結(jié)構(gòu)的動(dòng)力特性,也會(huì)改變結(jié)構(gòu)在時(shí)域響應(yīng)的統(tǒng)計(jì)特性[20]。時(shí)間序列模型可以將大量結(jié)構(gòu)響應(yīng)數(shù)據(jù)所蘊(yùn)含的信息凝聚為少量的模型參數(shù),例如自回歸系數(shù)、滑動(dòng)平均系數(shù)和模型殘差方差等,其中的自回歸系數(shù)內(nèi)含了結(jié)構(gòu)的固有動(dòng)力特性,可以通過(guò)機(jī)器學(xué)習(xí)算法建立結(jié)構(gòu)自回歸系數(shù)變化與損傷的聯(lián)系,進(jìn)行損傷識(shí)別與定位。
對(duì)于一個(gè)平穩(wěn)、正態(tài)、零均值的時(shí)間序列{Xt},若xt與前p步取值xt-1,xt-2,…xt-p,前q步輸入激勵(lì)at-1,at-2,…at-q有關(guān),則時(shí)間序列模型可以表示為:
(1)
式中:p是自回歸模型的階數(shù);q是滑動(dòng)平均模型的階數(shù),模型記為ARMA(p,q),φ和θ分別是時(shí)間序列和輸入激勵(lì)的各階特定系數(shù)。對(duì)式(1)引入后移算子B,使得Bixt=xt-i,模型可以表示為:
(2)
(3)
上式反映了系統(tǒng)的脈沖響應(yīng)函數(shù),分子θ(B)包含了結(jié)構(gòu)與外部之間的交換關(guān)系,分母φ(B)包含了結(jié)構(gòu)的固有特性信息。利用結(jié)構(gòu)在外界激勵(lì)下的位移/加速度時(shí)程建立ARMA模型,模型的自回歸系數(shù)包含了結(jié)構(gòu)固有特性信息。
對(duì)于發(fā)生損傷的結(jié)構(gòu),其總體剛度矩陣為:
(4)
式中:ki為第i個(gè)構(gòu)件未發(fā)生損傷時(shí)的剛度矩陣;M為構(gòu)件的數(shù)量;定義αi為第i個(gè)構(gòu)件的損傷系數(shù),αi∈[0,1],未發(fā)生損傷時(shí)αi=0,完全損傷時(shí)αi=1。
結(jié)構(gòu)發(fā)生損傷后,模態(tài)參數(shù)會(huì)發(fā)生變化,反映結(jié)構(gòu)固有特性信息的自回歸系數(shù)會(huì)發(fā)生變化,但是自回歸系數(shù)的變化與結(jié)構(gòu)損傷系數(shù)αi不一定是線性變化的。假設(shè)結(jié)構(gòu)有p個(gè)測(cè)點(diǎn),對(duì)每個(gè)測(cè)點(diǎn)信號(hào)可以提取前q階自回歸系數(shù),共可獲得s=p×q個(gè)自回歸系數(shù)。定義自回歸系數(shù)向量為:
{φ}=[φ1,φ2…φs]T
(5)
對(duì)自回歸系數(shù)φi(i∈[1,s])做一階變分:
(6)
自回歸系數(shù)向量的一階變分可以表示為:
{δφ}=[δφ1,δφ2,…δφs]T=
(7)
式中:P為s行、M列矩陣;Pij反映了損傷系數(shù)αj的單位變化引起第i個(gè)自回歸系數(shù)的變化量(i∈[1,s],j∈[1,M])。
自回歸系數(shù)的一階變分反映了損傷結(jié)構(gòu)自回歸系數(shù)與未損傷結(jié)構(gòu)的自回歸系數(shù)之差:
(8)
(9)
可以得到:
{δφ}={φ}-{φ0}=P{δα}=P{α}
(10)
進(jìn)而求出結(jié)構(gòu)損傷系數(shù)向量{α}:
{α}=P+[{φ}-{φ0}]
(11)
式中:P+為P的偽逆矩陣,定義為損傷識(shí)別矩陣。該矩陣只與結(jié)構(gòu)自身的動(dòng)力特性有關(guān),若能求出P+矩陣,便可以乘上損傷結(jié)構(gòu)與未損傷結(jié)構(gòu)的自回歸系數(shù)向量之差求出結(jié)構(gòu)的損傷系數(shù)向量,從而識(shí)別損傷的位置與大小。
有大量的損傷樣本與對(duì)應(yīng)的自回歸系數(shù)時(shí),可以通過(guò)機(jī)器學(xué)習(xí)(Machine Learning)對(duì)損傷識(shí)別矩陣進(jìn)行估計(jì)。自1950年圖靈提出機(jī)器學(xué)習(xí)的概念以來(lái),其理論和方法發(fā)展迅速,在土木工程的多個(gè)細(xì)分領(lǐng)域已被廣泛應(yīng)用[21-22],對(duì)損傷程度的預(yù)測(cè)適用于連續(xù)型有監(jiān)督的回歸算法(Regression),首先考慮使用普通的線性回歸對(duì)結(jié)構(gòu)損傷系數(shù)進(jìn)行回歸預(yù)測(cè):
(12)
(13)
式中:N是訓(xùn)練集中的樣本數(shù)量。其使用最小二乘法求得的解為:
(14)
求解線性回歸的前提條件是對(duì)給定數(shù)據(jù)集X,XTX可逆。結(jié)構(gòu)監(jiān)測(cè)系統(tǒng)存在多個(gè)測(cè)點(diǎn),每個(gè)測(cè)點(diǎn)獲得的自回歸系數(shù)階數(shù)較高,樣本的特征數(shù)量s較高,部分特征之間可能高度相關(guān),則XTX接近于奇異的病態(tài)矩陣。上述問(wèn)題變?yōu)椴贿m定問(wèn)題,計(jì)算誤差很大,在最小化目標(biāo)中加一個(gè)懲罰項(xiàng),使得目標(biāo)函數(shù)變?yōu)?
(15)
這種回歸稱(chēng)為嶺回歸(Ridge Regression)。嶺回歸可以限制模型的復(fù)雜度,防止過(guò)擬合(Overfitting),使得模型在復(fù)雜度和性能間達(dá)到平衡,可以有效應(yīng)用于結(jié)構(gòu)損傷定位。λ是待定的參數(shù),嶺回歸的關(guān)鍵是找到合理的λ來(lái)平衡模型的損失函數(shù)和偏差。使用嶺回歸求得系數(shù)的解為:
(16)
(17)
這種回歸稱(chēng)為核嶺回歸,目標(biāo)函數(shù)可以表示為:
(18)
K(x(i),x(j))=Φ(x(i))·Φ(x(j))=
(19)
利用時(shí)間序列與核嶺回歸進(jìn)行結(jié)構(gòu)損傷識(shí)別與定位,需建立有限元模型,在模型中模擬結(jié)構(gòu)的損傷。假設(shè)結(jié)構(gòu)有p層,識(shí)別精度僅限于損傷發(fā)生的樓層,損傷模式記為{α}=[α1,α2,…,αp],共隨機(jī)生成m種破壞模式。
對(duì)m個(gè)損傷模型分別輸入環(huán)境激勵(lì),每個(gè)模型輸出p條響應(yīng)。對(duì)每條輸出響應(yīng)均建立ARMA模型,每個(gè)ARMA模型提取前q個(gè)自回歸系數(shù),得到s=p×q個(gè)自回歸系數(shù)作為后續(xù)核嶺回歸輸入特征,自回歸系數(shù)記為{φ}=[φ1,φ2…φs]T。
將每種損傷模式的s個(gè)自回歸系數(shù)作為核嶺回歸的輸入特征,損傷系數(shù)作為輸出變量,對(duì)模型進(jìn)行訓(xùn)練。完成訓(xùn)練的模型可以接受結(jié)構(gòu)發(fā)生未知損傷后提取的自回歸系數(shù),輸出預(yù)測(cè)的損傷系數(shù)。
本文采用4層鋼筋混凝土框架結(jié)構(gòu)作為數(shù)值算例,該結(jié)構(gòu)長(zhǎng)3.0 m、寬2.0 m、層高2.6 m,梁均為高0.66 m、寬0.25 m的矩形梁,柱為0.4 m×0.6 m矩形截面柱,采用開(kāi)源有限元軟件OpenSees建立結(jié)構(gòu)有限元模型。該模型采用線彈性模型,樓面荷載3 kN/m2,混凝土采用C30,彈性模量為3×104MPa。模型如圖1所示。
圖1 混凝土框架OpenSees有限元模型Fig.1 OpenSees FEM model of the concrete frame
沿結(jié)構(gòu)短邊方向(Y向)在底部輸入白噪聲激勵(lì),激勵(lì)時(shí)長(zhǎng)30 s,采樣間隔為0.01 s,采樣頻率為100 Hz。未損傷結(jié)構(gòu)在白噪聲激勵(lì)下一層(圖1中A點(diǎn))及4層(圖1中D點(diǎn))加速度時(shí)程如圖2所示。
圖2 混凝土框架白噪聲激勵(lì)下加速度時(shí)程Fig.2 Acceleration time history of concrete frame under white noise excitation
對(duì)受損傷的結(jié)構(gòu)輸入白噪聲后得到4層加速度時(shí)程(圖1中A,B,C,D點(diǎn)),分別建立ARMA模型,該模型要求輸入的時(shí)間序列是平穩(wěn)的,且首先對(duì)時(shí)間序列的穩(wěn)定性進(jìn)行單位根檢驗(yàn)(ADF),然后對(duì)模型進(jìn)行定階。
選擇合適的模型階數(shù)非常重要,模型的階數(shù)越大擬合性越強(qiáng),但階數(shù)過(guò)高誤差會(huì)加大,造成過(guò)擬合[23]。常用的定階方法包括自相關(guān)系數(shù)法、偏相關(guān)系數(shù)法、赤池信息量準(zhǔn)則法(AIC)及貝葉斯信息準(zhǔn)則(BIC)等。本文對(duì)模型采用AIC熱力圖定階,根據(jù)經(jīng)驗(yàn)取AR階數(shù)1~30,MA階數(shù)1~5,分別計(jì)算對(duì)應(yīng)的ARMA模型的AIC,選取使得AIC最小的AR與MA階數(shù)組合。
如圖3所示為未損傷結(jié)構(gòu)4層加速度時(shí)程ARMA模型的AIC熱力圖。隨著AR與MA階數(shù)的增大,AIC基本穩(wěn)定于2.6。模型的AR與MA階數(shù)可選擇(19,1)、(14,2)、(12,3)、(12,4)、(9,5)等組合,實(shí)際計(jì)算過(guò)程中優(yōu)先選擇MA階數(shù)較小的模型(19,1)。
圖3 ARMA模型AIC定階熱力圖Fig.3 AIC heatmap of ARMA model
為保證輸入信息的完整性,核嶺回歸模型取每條信號(hào)ARMA模型的前10階自回歸系數(shù),損傷系數(shù)為柱剛度折減系數(shù)。對(duì)于第三層柱剛度折減30%的損傷模型(α1= 0,α2=0,α3=30%,α4=0)自回歸系數(shù)分布如圖4所示。圖5中可以看到結(jié)構(gòu)受到損傷后自回歸系數(shù)系數(shù)發(fā)生了顯著改變,且不同層自回歸系數(shù)變化不同。
圖4 損傷結(jié)構(gòu)前10階自回歸系數(shù)Fig.4 Top 10 autoregressive coefficients of damaged structure
圖5 損傷結(jié)構(gòu)與未損傷結(jié)構(gòu)前10階自回歸系數(shù)之差Fig.5 Difference between the top 10 autoregressive coeffici- ents of damaged and undamaged structures
本文通過(guò)對(duì)上述有限元模型進(jìn)行隨機(jī)構(gòu)造損傷,共生成損傷模型1 200個(gè),分別輸入白噪聲激勵(lì)后得到4 800條加速度時(shí)程數(shù)據(jù),對(duì)每條加速度時(shí)程分別建立ARMA模型,每個(gè)ARMA模型提取前10階自回歸系數(shù),每個(gè)損傷模型對(duì)應(yīng)40個(gè)自回歸系數(shù)以及4個(gè)輸出變量(1層、2層、3層、4層的結(jié)構(gòu)損傷系數(shù))。將數(shù)據(jù)集進(jìn)行分割,隨機(jī)選取1 000個(gè)損傷模型作為訓(xùn)練集,其余200個(gè)模型作為驗(yàn)證集。
機(jī)器學(xué)習(xí)中超參數(shù)選擇對(duì)模型的預(yù)測(cè)精度至關(guān)重要。核嶺回歸中需要確定的超參數(shù)包括:(1)正則化系數(shù)λ,表征正則化的強(qiáng)弱,λ過(guò)小可能出現(xiàn)過(guò)擬合,λ過(guò)大可能使得模型回歸系數(shù)太小,出現(xiàn)欠擬合;(2)核函數(shù)(Kernel),不同核函數(shù)泛化能力不同,平滑程度不同,適用的樣本也不同,本文依據(jù)經(jīng)驗(yàn)選取最為常見(jiàn)的核函數(shù)徑向基函數(shù);(3)徑向基函數(shù)參數(shù)γ代表RBF的幅寬,會(huì)影響模型的泛化性能。需要確定的超參數(shù)只有λ和γ,且超參數(shù)范圍已知,可以利用窮舉搜索找出使得模型表現(xiàn)最好的超參數(shù)。
本文中對(duì)4個(gè)輸出變量的預(yù)測(cè)分為4個(gè)模型,確定超參數(shù)后對(duì)模型進(jìn)行訓(xùn)練,如圖6所示,單個(gè)模型訓(xùn)練時(shí)長(zhǎng)約為4~5 s,模型整體訓(xùn)練總時(shí)長(zhǎng)約16~20 s。
圖6 模型訓(xùn)練樣本量與訓(xùn)練時(shí)間Fig.6 Model training sample size and training time
如圖7所示,單層的核嶺回歸模型R2在訓(xùn)練集上可以達(dá)到一層95%、二層84%、三層81%、四層88%,在驗(yàn)證集上可以達(dá)到一層81%、二層67%、三層68%、四層66%。訓(xùn)練樣本數(shù)量超過(guò)400后模型表現(xiàn)趨于穩(wěn)定。
圖7 模型R2學(xué)習(xí)曲線(R2小于0時(shí),圖中體現(xiàn)為0)Fig.7 Model R2 learning curves
為了驗(yàn)證核嶺回歸的優(yōu)勢(shì),使用無(wú)正則化無(wú)核函數(shù)的線性回歸(Liner Regression)、有正則化無(wú)核函數(shù)的嶺回歸(Ridge Regression)、無(wú)正則化有核函數(shù)的支持向量回歸(Support Vector Regression)三種算法與其進(jìn)行對(duì)比。
如圖8、圖9所示,在訓(xùn)練集上,線性回歸表現(xiàn)最好,R2在訓(xùn)練過(guò)程中維持在95%以上,均方誤差(MSE)維持在0.001以上,但在驗(yàn)證集上線性回歸表現(xiàn)很差,R2維持在0及0以下,模型出現(xiàn)過(guò)擬合。
圖8 不同回歸算法模型R2學(xué)習(xí)曲線 (R2小于0時(shí),圖中體現(xiàn)為0)Fig.8 Model R2 learning curves with different regression methods
圖9 不同回歸算法模型MSE曲線Fig.9 Model MSE curves with different regression methods
表1為最終幾種回歸模型的表現(xiàn),嶺回歸、支持向量回歸的表現(xiàn)介于線性回歸與核嶺回歸之間,結(jié)果說(shuō)明添加正則化與核函數(shù)后可以顯著提升模型的擬合性能與泛化性能。
表1 不同回歸模型訓(xùn)練時(shí)間及表現(xiàn)對(duì)比Table 1 Training time and performance of different regression models
為驗(yàn)證模型損傷識(shí)別與定位的準(zhǔn)確性,選取5種工況進(jìn)行驗(yàn)證,分別為未受損工況、兩種單損傷工況(結(jié)構(gòu)僅一處發(fā)生損傷)及兩種多損傷工況(結(jié)構(gòu)多處發(fā)生損傷)。損傷情況如表2所列。
表2 模型損傷工況Table 2 Damage conditions for model
對(duì)未受損的原結(jié)構(gòu),核嶺回歸模型識(shí)別一至四層的損傷系數(shù)分別為4.6%、0.4%、3.2%、4.2%,平均誤差為3.1%;對(duì)于單損傷工況1,一層損傷10%,模型識(shí)別一層的損傷為8.9%,一至四層平均識(shí)別誤差為2.7%;對(duì)于單損傷工況2,一層損傷30%,模型識(shí)別一層的損傷為29.0%,一至四層平均識(shí)別誤差2.7%;對(duì)于多損傷工況1,模型一層、二層、三層發(fā)生10%、30%、30%損傷,模型識(shí)別結(jié)果為10.1%、28.9%、24.3%,一至四層平均識(shí)別誤差2.7%;對(duì)于多損傷模型2,模型二層、三層、四層發(fā)生10%、30%、50%損傷,模型識(shí)別結(jié)果為9.3%、30.6%、55.2%,一至四層平均識(shí)別誤差3.3%。由表2可看出,模型在多種工況下均表現(xiàn)出較高的識(shí)別準(zhǔn)確率。損傷識(shí)別完整結(jié)果列于表3。
表3 損傷識(shí)別結(jié)果Table 3 Damage identification results
本文通過(guò)理論分析并結(jié)合數(shù)值模擬,建立了核嶺回歸模型,利用結(jié)構(gòu)時(shí)間序列模型中自回歸系數(shù)實(shí)現(xiàn)了結(jié)構(gòu)損傷識(shí)別與定位:
(1) 提出結(jié)構(gòu)的損傷識(shí)別矩陣P+,該矩陣只與結(jié)構(gòu)自身的動(dòng)力特性有關(guān),若能求出該矩陣便可以通過(guò)損傷結(jié)構(gòu)與未損傷結(jié)構(gòu)的自回歸系數(shù)向量之差識(shí)別結(jié)構(gòu)的損傷位置與大小。
(2) 通過(guò)4層框架結(jié)構(gòu)數(shù)值模型對(duì)理論進(jìn)行驗(yàn)證,模型R2在訓(xùn)練集上可以達(dá)到80%,在驗(yàn)證集上可以達(dá)到65%。對(duì)多個(gè)單損傷和多損傷工況進(jìn)行驗(yàn)證,模型識(shí)別平均誤差不超過(guò)3.7%。
(3) 相比無(wú)正則化無(wú)核函數(shù)的線性回歸、有正則化無(wú)核函數(shù)的嶺回歸、無(wú)正則化有核函數(shù)的支持向量回歸三種算法,核嶺回歸準(zhǔn)確率最高,說(shuō)明核函數(shù)和正則化可以有效提高模型的擬合性能與泛化性能,更好地應(yīng)用于結(jié)構(gòu)損傷識(shí)別。
本文不足之處在于需要建立結(jié)構(gòu)較為精確的有限元模型,且需對(duì)結(jié)構(gòu)各種損傷下的情況進(jìn)行數(shù)值模擬,建立ARMA模型,需要較高的計(jì)算成本。