第一作者占超男，碩士生，1990年生

通信作者李東升男，副教授，1972年生

郵箱：dsli@dlut.edu.cn

基于改進(jìn)交叉模型交叉模態(tài)法的局部損傷識(shí)別方法

占超， 李東升， 任亮， 李宏男

(大連理工大學(xué)海岸及近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧大連116024)

摘要：傳統(tǒng)的交叉模型交叉模態(tài)(CMCM)法由于其核心矩陣的缺秩使其在全局修正時(shí)的解不唯一，為得到唯一解必須人為假定約束。以往利用CMCM法進(jìn)行損傷識(shí)別的研究中，通常將損傷前后質(zhì)量保持不變作為約束求得各單元?jiǎng)偠鹊男拚?，然后將剛度的下降作為判斷損傷的依據(jù)。然而，結(jié)構(gòu)在極端運(yùn)營(yíng)條件下質(zhì)量的變化也是損傷的一種。為此，對(duì)傳統(tǒng)的CMCM法進(jìn)行了改進(jìn)，提出了以下的改進(jìn)方法：首先，由基準(zhǔn)狀態(tài)下的有限元模型和損傷后實(shí)測(cè)的結(jié)構(gòu)低階模態(tài)信息求得核心矩陣，將核心矩陣最小奇異值對(duì)應(yīng)的右奇異向量作為損傷指示向量(DIV)；然后，基于損傷的局部性利用聚類(lèi)分析算法自動(dòng)識(shí)別DIV中的異常元素，將異常元素作為損傷定位的依據(jù)；最后，由DIV估計(jì)各單元質(zhì)量和剛度的損傷程度。改進(jìn)后方法的優(yōu)勢(shì)在于：無(wú)需人為主觀地添加約束即可一次性求出結(jié)構(gòu)各單元質(zhì)量和剛度損傷前后的變化，避免了錯(cuò)誤或不準(zhǔn)確的假設(shè)給損傷識(shí)別結(jié)果帶來(lái)的誤差。通過(guò)數(shù)值實(shí)驗(yàn)對(duì)該改進(jìn)方法的可行性、魯棒性以及損傷敏感性進(jìn)行了研究，并通過(guò)美國(guó)洛斯阿拉莫斯國(guó)家實(shí)驗(yàn)室(LANL)的一個(gè)4自由度板柱結(jié)構(gòu)振動(dòng)臺(tái)實(shí)驗(yàn)對(duì)改進(jìn)方法的有效性做了進(jìn)一步的驗(yàn)證。

關(guān)鍵詞：損傷識(shí)別； 全局修正； 交叉模型交叉模態(tài)法； 結(jié)構(gòu)健康監(jiān)測(cè)

收稿日期：2014-01-02修改稿收到日期：2014-04-03

中圖分類(lèi)號(hào):TP206+.3

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2015.07.021

Abstract:Since the core matrix in the traditional cross-model cross-mode(CMCM) method is rank-deficient, constraints must be applied artificially for getting a unique solution in a complete model-updating situation. In previous damage identification researches based on the CMCM method, structural masses are often supposed to be unchanged for getting the updated coefficients of stiffness, and the reduction of stiffness is regarded as an indicator of damage. However, elemental mass changes under operational conditions are obviously a kind of damages. For detecting damage more effectively, the traditional CMCM method was improved, and an improved method was proposed as follows: evaluate the core matrix from the finite element model under the baseline condition and the measured modal data after damage and take the right singular vector corresponding to the least singular value of the core matrix as a damage indication vector(DIV); identify the abnormal elements in the DIV by using cluster analysis algorithm based on the assumption that damage appears locally and locate damage based on these abnormal elements; derive the extent of damage for each elemental mass and stiffness from the DIV. The advantage of the improved method is its ability to solve the mass and stiffness changes before and after damage for all elements without artificial constraints, thus the error of damage identification caused by imprecise or wrong constraints vanishes. The improved method’s practicability, robustness, and sensitivity to damage were analysed via a numerical simulation, and its effectiveness is further verified by virtue of the experiments on a 4 degrees of freedom test-bed structure done at the Los Alamos National Laboratory.

Local damage identification approach based on improved cross-model cross-mode method

ZHANChao,LIDong-sheng,RENLiang,LIHong-nan(State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024,China)

Key words:damage identification; complete-updating; cross-modal cross-mode; structural health monitoring

隨著我國(guó)經(jīng)濟(jì)的不斷發(fā)展，大量土木工程結(jié)構(gòu)(尤其是高層建筑和大跨橋梁)伴隨著人們的需要如雨后春筍般涌現(xiàn)。這些結(jié)構(gòu)在以后的服役過(guò)程中，由于各種因素(如地震、材料銹蝕等)某些構(gòu)件往往會(huì)出現(xiàn)損傷。損傷的存在勢(shì)必會(huì)降低結(jié)構(gòu)的使用性能，放任損傷的加劇甚至可能引發(fā)重大的土木工程結(jié)構(gòu)事故(如橋梁的突然折斷、房屋的驟然倒塌等)，從而造成重大的人員傷亡和財(cái)產(chǎn)損失。因此，對(duì)重大工程結(jié)構(gòu)的性能進(jìn)行監(jiān)測(cè)和診斷，及時(shí)發(fā)現(xiàn)結(jié)構(gòu)的損傷，對(duì)可能出現(xiàn)的災(zāi)害進(jìn)行預(yù)測(cè)并對(duì)結(jié)構(gòu)的安全性能進(jìn)行正確的評(píng)價(jià)，具有極其重要的意義，同時(shí)也是現(xiàn)代工程的迫切要求和土木工程學(xué)科發(fā)展的一個(gè)重要領(lǐng)域[1]。

結(jié)構(gòu)健康監(jiān)測(cè)和損傷識(shí)別近幾十年來(lái)在硬件和理論上均有很大進(jìn)展[2]，其中基于模型修正的損傷識(shí)別方法受到諸多學(xué)者的關(guān)注[3]。該類(lèi)方法的基本思想為根據(jù)健康狀態(tài)下的實(shí)測(cè)振動(dòng)數(shù)據(jù)對(duì)結(jié)構(gòu)的初始有限元模型進(jìn)行修正，得到基準(zhǔn)模型；結(jié)構(gòu)損傷后，由損傷狀態(tài)下的實(shí)測(cè)振動(dòng)數(shù)據(jù)對(duì)基準(zhǔn)模型進(jìn)行修正，則損傷前后各單元模型參數(shù)的變化即反映出損傷的位置和程度。

傳統(tǒng)的模型修正方法[5]主要有：利用模態(tài)信息的矩陣直接修正法、利用模態(tài)信息的迭代修正法和利用實(shí)測(cè)頻響函數(shù)的修正法等。這些方法或多或少存在各種缺點(diǎn)，比如直接修正法修正后的質(zhì)量、剛度矩陣無(wú)法保證原有的對(duì)稱性、稀疏性、和正定性，缺乏物理意義；迭代修正法存在解的收斂性以及迭代過(guò)程中的誤差累計(jì)問(wèn)題，并且計(jì)算量較大；而頻響函數(shù)修正法的修正參數(shù)估計(jì)為有偏估計(jì)。為了克服傳統(tǒng)方法的不足，近年來(lái)不斷有學(xué)者提出新的模型修正方法。Doebling等[6]提出通過(guò)敏感性分析得到單元連接矩陣，然后由最小秩法求出單元?jiǎng)偠葏?shù)的擾動(dòng)量，并使修正后的剛度矩陣保持了原有的單元連接形式；Kaouk等[7]提出將模型修正分為兩步，即先進(jìn)行誤差定位，再使用最小秩擾動(dòng)法對(duì)結(jié)構(gòu)質(zhì)量、阻尼、剛度矩陣同時(shí)修正；侯吉林[8]提出基于局部模態(tài)的約束子結(jié)構(gòu)模態(tài)修正法，并對(duì)此進(jìn)行了系統(tǒng)地研究。

Hu等[9]提出了交叉模型交叉模態(tài)(CMCM)算法。該方法可以根據(jù)測(cè)得的有限低階模態(tài)信息對(duì)單元質(zhì)量和單元?jiǎng)偠染仃囃瑫r(shí)修正，而且無(wú)需模態(tài)配對(duì)，對(duì)振型沒(méi)有同比例要求，是一種直接的物理屬性修正方法[10]。本文先對(duì)傳統(tǒng)的CMCM法推導(dǎo)過(guò)程進(jìn)行簡(jiǎn)要介紹；然后指出其中的不足并提出改進(jìn)方法；接著通過(guò)一個(gè)數(shù)值實(shí)驗(yàn)對(duì)改進(jìn)方法的可行性、魯棒性和損傷敏感性進(jìn)行研究；最后通過(guò)美國(guó)洛斯阿拉莫斯國(guó)家實(shí)驗(yàn)室(LANL)的一個(gè)4自由度板柱結(jié)構(gòu)實(shí)驗(yàn)對(duì)改進(jìn)方法的可行性做進(jìn)一步驗(yàn)證。

1交叉模型交叉模態(tài)法

記有限元模型的剛度和質(zhì)量矩陣分別為K和M，則

KΦ=MΦΛ

(1)

式中Φ為有限元模型振型矩陣，Λ為與各振型對(duì)應(yīng)的圓頻率的平方值構(gòu)成的Ni階對(duì)角陣。

類(lèi)似式(1)，對(duì)于實(shí)際結(jié)構(gòu)有

K*Φ*=M*Φ*Λ*

(2)

式中，K*和M*分別為實(shí)際結(jié)構(gòu)的剛度和質(zhì)量陣，Φ*實(shí)測(cè)振型矩陣，Λ*為實(shí)測(cè)圓頻率的平方值構(gòu)成的Nj階對(duì)角陣。

式(2)等號(hào)兩邊同時(shí)左乘的Φ轉(zhuǎn)置，即

ΦTK*Φ*=ΦTM*Φ*Λ*

(3)

假設(shè)實(shí)際結(jié)構(gòu)的剛度和質(zhì)量矩陣滿足下列方程，

(4)

式中Kn與Mn分別為在整體坐標(biāo)系下標(biāo)號(hào)為n的單元?jiǎng)偠群蛦卧|(zhì)量矩陣，αn和βn分別為n號(hào)單元?jiǎng)偠染仃囆拚禂?shù)和n號(hào)單元質(zhì)量矩陣修正系數(shù)，NK和NM分別為單元?jiǎng)偠染仃嚭蛦卧|(zhì)量矩陣的個(gè)數(shù)。

將式(4)代入式(3)，可得

(5)

式中

(6)

式(5)可以寫(xiě)成如下矩陣形式，

Gγ=f

(7)

式中

(8)

式(7)中矩陣G的行數(shù)，亦方程數(shù)Nm=NiNj；其列數(shù)，亦未知數(shù)個(gè)數(shù)Nu=NK+NM。由于矩陣G是CMCM算法的核心部分，因此也被稱為核心矩陣[11-12]。需要指出的是，當(dāng)實(shí)測(cè)模態(tài)測(cè)得足夠準(zhǔn)時(shí)，即使式(7)中方程數(shù)Nm大于未知數(shù)個(gè)數(shù)Nu，由于G的缺秩，式(7)仍會(huì)出現(xiàn)解不唯一的情況。要想得到唯一解，必須添加約束。假設(shè)約束方程如下，

(9)

若假設(shè)結(jié)構(gòu)質(zhì)量保持不變，則

(10)

式中,I為維數(shù)為NM的單位方陣。將式(7)與式(10)合并，即得

(11)

式中

(12)

(13)

2改進(jìn)的CMCM局部損傷識(shí)別法

在第1節(jié)介紹的傳統(tǒng)CMCM法中，由于核心矩陣G的缺秩，使其在全局修正時(shí)解不唯一。要想得到唯一解，至少需要人為地添加一個(gè)約束。錯(cuò)誤或者不精確的約束都可能會(huì)引起傳統(tǒng)CMCM法修正結(jié)果出現(xiàn)誤差。以往在對(duì)基于CMCM法的損傷識(shí)別研究中，往往假設(shè)各單元質(zhì)量保持不變，而僅對(duì)各單元?jiǎng)偠冗M(jìn)行修正，并以單元?jiǎng)偠鹊南陆底鳛閾p傷判斷的依據(jù)。然而，結(jié)構(gòu)在極端運(yùn)營(yíng)條件下質(zhì)量的變化(例如暴雪導(dǎo)致房屋頂層質(zhì)量的顯著增加等)也是損傷的一種。本節(jié)對(duì)傳統(tǒng)CMCM法的推導(dǎo)過(guò)程稍加變動(dòng)，提出改進(jìn)的CMCM局部損傷識(shí)別法，如下：

假設(shè)損傷后剛度、質(zhì)量陣K*和M*滿足下列關(guān)系

(14)

式中ηn和εn分別為n號(hào)單元?jiǎng)偠染仃嚭蚽號(hào)單元質(zhì)量矩陣的“損傷比”，與式(4)中的αn和βn具有如下關(guān)系

ηn=1+αn，εn=1+βn

(15)

將式(14)代入式(3)，可得

(16)

將式(16)寫(xiě)成矩陣形式

Gv=0

(17)

其中Cn、Dn和G分別同式(6)和式(8)的定義。v定義如下

(18)

向量v的前半段和后半段分別記錄了各單元?jiǎng)偠取①|(zhì)量損傷前后的比值，求得了v即實(shí)現(xiàn)了對(duì)結(jié)構(gòu)參數(shù)損傷前后變化的全面了解。由式(17)可見(jiàn)，只要v有解，核心矩陣G必缺秩，v必有無(wú)窮多個(gè)解。求解v的問(wèn)題即化為求解G的零空間的問(wèn)題。

數(shù)學(xué)上往往采用奇異值分解的方法求解一個(gè)矩陣的零空間，即

G=USVT

(19)

式中S為Nu×Nu的對(duì)角矩陣，其對(duì)角線上記錄了G的Nu個(gè)從大到小排列的奇異值s1,s2,…,sNu；U(Nm×Nu)和V(Nu×Nu)中各列向量分別稱為G的各奇異值對(duì)應(yīng)的左奇異向量和右奇異向量，并滿足

UTUINu，VTV=INu

(20)

則奇異值s1,s2,…,sNu中0奇異值對(duì)應(yīng)的右奇異向量即為零空間的基。

然而，由于儀器的內(nèi)部噪聲、計(jì)算機(jī)的舍入誤差、環(huán)境噪聲等因素，實(shí)際測(cè)得的模態(tài)中含有噪聲而使得G趨于滿秩，即所有奇異值均大于0。為此，我們?nèi)∑娈愔抵凶钚〉闹祍Nu對(duì)應(yīng)的右奇異向量作為解v的基，記為v*(v*滿足v*Tv*=1)。

v與v*間滿足倍數(shù)關(guān)系，即v=av*，其中a為待定系數(shù)。我們知道，損傷是一種典型的局部現(xiàn)象，兩次動(dòng)力測(cè)試中結(jié)構(gòu)大部分的參數(shù)應(yīng)該基本保持不變，只有損傷處參數(shù)會(huì)有較顯著的變化，即：v中大部分元素都應(yīng)該近似等于1，對(duì)應(yīng)地v*中大部分元素應(yīng)處于同一水平；損傷情況下v*中會(huì)出現(xiàn)少數(shù)異常元素，表現(xiàn)在與其他正常元素不在同一水平線上。v*中的異常元素可以作為損傷定位的依據(jù)。由于向量v*中包含結(jié)構(gòu)的損傷信息，我們又稱v*為損傷指示向量(Damage Indication Vector,DIV)。DIV中正常元素和異常元素的劃分可以通過(guò)聚類(lèi)分析算法[13]自動(dòng)實(shí)現(xiàn)。聚類(lèi)分析能夠從樣本數(shù)據(jù)出發(fā)，計(jì)算樣本之間的空間距離，根據(jù)樣本彼此間的親疏關(guān)系，對(duì)樣本進(jìn)行自動(dòng)分類(lèi)。

(21)

(22)

圖1　四自由度鏈?zhǔn)?結(jié)構(gòu)(數(shù)值實(shí)驗(yàn)) Fig.1 4-DOF chain-like system (for numerical simulation)

則av*即為我們對(duì)解v的估計(jì)，(av*-1)×100%即為我們對(duì)各單元?jiǎng)偠?、質(zhì)量損傷程度的估計(jì)。

3數(shù)值實(shí)驗(yàn)

3.1改進(jìn)方法的可行性

本小節(jié)通過(guò)數(shù)值實(shí)驗(yàn)對(duì)第2節(jié)中改進(jìn)方法的可行性進(jìn)行說(shuō)明。如圖1所示一四自由度鏈?zhǔn)浇Y(jié)構(gòu)，底層柱與小車(chē)相連，結(jié)構(gòu)整體可以在水平地面上自由滑動(dòng)。各集中質(zhì)量?jī)H考慮其在水平方向上的運(yùn)動(dòng)。各層質(zhì)量、小車(chē)質(zhì)量以及各層柱的水平剛度如圖中所示。實(shí)驗(yàn)中，結(jié)構(gòu)的單處損傷通過(guò)將第一層的層間剛度k1折減50%實(shí)現(xiàn)；多處損傷工況下則將第二層的層間剛度k2折減20%，并將第一層的集中質(zhì)量增加20%以模擬極端運(yùn)營(yíng)條件下質(zhì)量的增加。數(shù)值實(shí)驗(yàn)中損傷狀態(tài)下的模態(tài)是通過(guò)損傷后的有限元模型直接計(jì)算得到的。結(jié)構(gòu)基準(zhǔn)狀態(tài)下的模態(tài)和損傷狀態(tài)下的“實(shí)測(cè)”模態(tài)分別如表1和表2所示。結(jié)構(gòu)的第一階模態(tài)為平動(dòng)模態(tài)，在表中不予顯示。

表1　基準(zhǔn)狀態(tài)下有限元模型的模態(tài)

將表1和表2中的模態(tài)數(shù)據(jù)代入式(17)，即可求得損傷指示向量(DIV)，進(jìn)而求得損傷位置和程度。無(wú)噪聲情況下，單損傷工況和多損失工況的損傷識(shí)別結(jié)果如圖2所示。圖中可見(jiàn)，無(wú)論是單損傷工況還是多損傷工況下，DIV中異常元素均正確指示出損傷位置。單損傷工況中，η1被標(biāo)記為異常元素，表明一號(hào)單元?jiǎng)偠瘸霈F(xiàn)異常；多損傷工況中，η2和ε2均被標(biāo)記為異常，表明二號(hào)單元?jiǎng)偠纫约岸?hào)單元質(zhì)量(即圖1中一層質(zhì)量)出現(xiàn)了異常。兩種工況下，損傷程度識(shí)別值均與預(yù)設(shè)值均很好地吻合，從而說(shuō)明改進(jìn)方法可行。

表2　損傷狀態(tài)下的結(jié)構(gòu)模態(tài)

圖2　損傷識(shí)別結(jié)果(數(shù)值實(shí)驗(yàn)) Fig.2 The results of damage identification(numerical simulation)

3.2改進(jìn)方法的魯棒性和損傷敏感性

有關(guān)小損傷的識(shí)別一直是損傷識(shí)別中的難點(diǎn)。小損傷一方面對(duì)結(jié)構(gòu)整體參數(shù)影響不大，難以從結(jié)構(gòu)的動(dòng)力響應(yīng)信號(hào)中辨別；另一方面又對(duì)結(jié)構(gòu)安全造成隱患，可能隨時(shí)發(fā)展成大損傷進(jìn)而導(dǎo)致結(jié)構(gòu)構(gòu)件失效甚至是結(jié)構(gòu)整體的破壞。同時(shí)，實(shí)際工程中來(lái)自方方面面的噪音也對(duì)損傷的正確識(shí)別造成困難。本小節(jié)以多損傷工況(圖1中第二層層間剛度折減y%，第一層集中質(zhì)量增加y%，y為損傷程度)為例，對(duì)改進(jìn)方法的魯棒性和損傷敏感性進(jìn)行研究。

圖3　改進(jìn)方法識(shí)別性能灰度圖 Fig.3 The grayscale of the improved method’s performance

我們以實(shí)測(cè)模態(tài)中的噪聲比例x為橫軸(分辨率0.1%)，以多損傷工況中的損傷程度y為縱軸(分辨率1%)，作出改進(jìn)方法的識(shí)別性能灰度圖，如圖 3所示。圖中各點(diǎn)(共50×20個(gè)點(diǎn))均重復(fù)實(shí)驗(yàn)100次。每次實(shí)驗(yàn)中，首先由該點(diǎn)縱坐標(biāo)對(duì)應(yīng)的損傷程度求得對(duì)應(yīng)的損傷有限元模型和損傷后的模態(tài)數(shù)據(jù)，然后在模擬的振型中加入該點(diǎn)橫坐標(biāo)對(duì)應(yīng)比例的比例噪聲，接著由含噪聲的模態(tài)數(shù)據(jù)利用改進(jìn)方法求得DIV，若DIV正確指示出了損傷位置，則記為一次成功的損傷識(shí)別。100次實(shí)驗(yàn)中，成功識(shí)別次數(shù)z越多，則識(shí)別性能越好，在圖中對(duì)應(yīng)的顏色也就越深。圖中可以看到，識(shí)別性能隨著損傷程度的減小和噪音比例的增大而變差。無(wú)噪聲條件下，任意程度的損傷均可準(zhǔn)確識(shí)別； 1%比例噪聲下，改進(jìn)方法可以對(duì)10%以上的損傷準(zhǔn)確定位(z>95)；而在5%的比例噪聲下，小損傷(y<5)的識(shí)別則完全無(wú)法實(shí)現(xiàn)。

4實(shí)驗(yàn)室實(shí)驗(yàn)

為進(jìn)一步證實(shí)改進(jìn)方法的可行性，本節(jié)對(duì)LANL實(shí)驗(yàn)室的一個(gè)4自由度板柱結(jié)構(gòu)振動(dòng)臺(tái)實(shí)驗(yàn)數(shù)據(jù)[14-15]進(jìn)行分析，并利用改進(jìn)方法對(duì)預(yù)設(shè)的損傷進(jìn)行識(shí)別。實(shí)驗(yàn)裝置如圖4所示。三層“樓板”和“基礎(chǔ)”均為鋁制金屬板(30.5 cm×30.5 cm×2.5 cm)，彼此通過(guò)四根鋁柱使用螺栓連接?！盎A(chǔ)”被限制在軌道上使其僅可以在x方向上滑動(dòng)。激振器采用(20 Hz~150 Hz)的限帶白噪聲對(duì)“基礎(chǔ)”進(jìn)行激勵(lì)。結(jié)構(gòu)上共安裝5個(gè)傳感器(通道1用于測(cè)量激振器對(duì)“基礎(chǔ)”施加的力，通道2、3、4、5分別用于測(cè)量“基礎(chǔ)”和一、二、三層“樓板”的加速度響應(yīng))，各傳感器采樣頻率均為320 Hz，采樣數(shù)目為8 192個(gè)。實(shí)驗(yàn)中將損傷歸為兩類(lèi)，一類(lèi)是線性損傷，另一類(lèi)是非線性損傷。線性損傷通過(guò)在樓板或基礎(chǔ)上添加集中質(zhì)量塊、替換某層的一根或兩根鋁柱(替換后的鋁柱的橫截面面積為原始的一半)以減小層間剛度來(lái)實(shí)現(xiàn)；非線性損傷[16]則通過(guò)調(diào)節(jié)頂層懸掛鋁柱和安裝在第二層樓板上的“碰撞器”之間的間距，在動(dòng)荷載下二者之間不同程度的相互碰撞來(lái)模擬。本文僅對(duì)其中的線性損傷數(shù)據(jù)進(jìn)行分析，如表3[17]所示。結(jié)構(gòu)基準(zhǔn)狀態(tài)下的各傳感器時(shí)程如圖5所示。

圖4　LANL 4自由度振動(dòng)臺(tái)試驗(yàn)裝置圖 Fig.3 Experiment setup for 4-DOF test-bed structure in LANL

表3　實(shí)驗(yàn)分組表 [17]

圖5　基準(zhǔn)狀態(tài)下各傳感器所測(cè)信號(hào)時(shí)程 Fig.5 Time histories of signals in each sensor under baseline condition

我們將激振器施加的力作為系統(tǒng)輸入，將“基礎(chǔ)”和各層“樓板”的加速度信號(hào)作為系統(tǒng)輸出，使用OKID[18](Observer Kalman filter IDentification)算法進(jìn)行模態(tài)識(shí)別。由于實(shí)測(cè)信號(hào)中的噪音和算法階數(shù)的不同取值，測(cè)得的模態(tài)中含有虛假模態(tài)和計(jì)算模態(tài)，即非物理模態(tài)。模態(tài)分析中可以通過(guò)穩(wěn)態(tài)圖[19]的方法鑒別物理模態(tài)和非物理模態(tài)。與非物理模態(tài)不同的是，物理模態(tài)不會(huì)隨算法階數(shù)的增加而變化，因而在穩(wěn)態(tài)圖中呈現(xiàn)近乎一條直線。圖6為基準(zhǔn)狀態(tài)下模態(tài)識(shí)別的穩(wěn)態(tài)圖。圖中每個(gè)實(shí)心點(diǎn)對(duì)應(yīng)一個(gè)識(shí)別的模態(tài)(為了使穩(wěn)態(tài)圖更加清晰，我們?yōu)V去了阻尼比大于5%的模態(tài))；點(diǎn)畫(huà)線為力對(duì)第三層樓板加速度的頻響函數(shù)的幅值。圖中可見(jiàn)，不同算法階數(shù)下物理模態(tài)的頻率具有較好的一致性，如虛線所示。限于篇幅，損傷狀態(tài)下各傳感器的時(shí)程和穩(wěn)態(tài)圖不再列出。不同工況下測(cè)得的模態(tài)具體數(shù)值見(jiàn)表 4。結(jié)構(gòu)的有限元模型參數(shù)與數(shù)值實(shí)驗(yàn)中圖1所示的各參數(shù)相同。將表4中基準(zhǔn)狀態(tài)下的實(shí)測(cè)模態(tài)與表1中的有限元模態(tài)相對(duì)比，可以發(fā)現(xiàn)二者比較吻合，但仍略有差異。為此我們將初始有限元模型和基準(zhǔn)狀態(tài)下的實(shí)測(cè)模態(tài)代入式(17)求得v*，然后由式(22)對(duì)初始有限元模型進(jìn)行全局修正。有限元模型修正參數(shù)及修正效果見(jiàn)表5。表中可見(jiàn)，修正后的有限元模型各階振型與實(shí)測(cè)振型的相似度(即MAC值)相比修正前均有所提高、第二階頻率誤差減小，但第三階和第四階的頻率誤差修正后反而略有增大。這主要是因?yàn)樾拚?，第三階和第四階的頻率誤差已經(jīng)很小(小于0.5%)，而全局修正是對(duì)結(jié)構(gòu)模型的全局優(yōu)化，因而使其誤差反而略有增大。

圖6　結(jié)構(gòu)基準(zhǔn)狀態(tài)下模態(tài)識(shí)別穩(wěn)態(tài)圖 Fig.6 Stability diagram under baseline condition

圖7　損傷識(shí)別結(jié)果(LANL) Fig.7 The results of damage identification(LANL)

將表 4中各損傷狀態(tài)下的實(shí)測(cè)頻率與基準(zhǔn)狀態(tài)下的實(shí)測(cè)頻率相對(duì)比，可以發(fā)現(xiàn)損傷后結(jié)構(gòu)各階頻率均有不同程度的下降。將損傷狀態(tài)下的實(shí)測(cè)模態(tài)和修正后的有限元模型代入式(17)求得DIV，進(jìn)而對(duì)損傷進(jìn)行定位和定量，識(shí)別結(jié)果如圖 7所示。圖中可見(jiàn)，各損傷工況下?lián)p傷指示向量中的異常元素均正確指示出損傷位置，損傷程度的估計(jì)也與表3中給出的預(yù)設(shè)損傷值很好地吻合(兩者間的誤差不超過(guò)5%)。我們還對(duì)工況2利用第1節(jié)中的傳統(tǒng)CMCM法進(jìn)行了計(jì)算， 假設(shè)結(jié)構(gòu)各單元質(zhì)量不發(fā)生變化，結(jié)果得出第一、二層的層間剛度大約下降10%，第三層層間剛度略有增大的錯(cuò)誤結(jié)論。改進(jìn)的CMCM局部損傷識(shí)別法中則無(wú)須人為主觀地給出約束，避免了錯(cuò)誤或者不準(zhǔn)確的約束給損傷識(shí)別結(jié)果帶來(lái)的誤差。

表4　不同工況下實(shí)測(cè)模態(tài)

表5　基準(zhǔn)狀態(tài)下有限元模型全局修正

5結(jié)論

交叉模型交叉模態(tài)法是一種很好的模型修正方法，具有諸多優(yōu)點(diǎn)，然而由于其方法中核心矩陣的缺秩，使得全局修正時(shí)必須人為添加約束。本文對(duì)傳統(tǒng)CMCM法進(jìn)行了改進(jìn)，提出了以下的改進(jìn)方法：首先，由基準(zhǔn)狀態(tài)下的有限元模型和損傷后實(shí)測(cè)的結(jié)構(gòu)低階模態(tài)信息求得核心矩陣，將核心矩陣最小奇異值對(duì)應(yīng)的右奇異向量作為損傷指示向量(DIV)；然后，基于損傷的局部性利用聚類(lèi)分析算法自動(dòng)識(shí)別DIV中的異常元素，將異常元素作為損傷定位的依據(jù)；最后，由DIV估計(jì)各單元質(zhì)量和剛度的損傷程度。改進(jìn)后方法的優(yōu)勢(shì)在于：無(wú)需人為主觀地添加約束即可一次性求出結(jié)構(gòu)各單元質(zhì)量和剛度損傷前后的變化，避免了錯(cuò)誤或不準(zhǔn)確的假設(shè)給損傷識(shí)別結(jié)果帶來(lái)的誤差。數(shù)值實(shí)驗(yàn)和美國(guó)LANL實(shí)驗(yàn)室實(shí)驗(yàn)均證實(shí)該改進(jìn)方法切實(shí)可行，可以對(duì)損傷進(jìn)行正確的定位并對(duì)損傷程度進(jìn)行準(zhǔn)確的估計(jì)(誤差不超過(guò)5%)。數(shù)值實(shí)驗(yàn)還表明該改進(jìn)方法可以在實(shí)測(cè)振型含1%的比例噪聲下對(duì)損傷程度大于10%的多損傷工況進(jìn)行準(zhǔn)確的損傷定位。關(guān)于如何進(jìn)一步提高改進(jìn)方法的魯棒性和對(duì)小損傷的敏感性，還有待進(jìn)一步的研究工作。

參考文獻(xiàn)

[1]李宏男，李東升. 土木工程結(jié)構(gòu)安全性評(píng)估、健康監(jiān)測(cè)及診斷述評(píng)[J]. 地震工程與工程振動(dòng)，2002，22 (3): 82-90.

LI Hong-nan, LI Dong-sheng. Safety assessment, health monitoring and damage diagnosis for structures in civil engineering[J]. Earthquake Engineering and Engineering Vibration, 2002,22(3):82-90.

[2]宗周紅, 任偉新, 阮毅. 土木工程結(jié)構(gòu)損傷診斷研究進(jìn)展[J]. 土木工程學(xué)報(bào), 2003, 36(5): 105-110.

ZONG Zhou-hong, REN Wei-xin, RUAN Yi. Recent advances in research on damage diagnosis for civil engineering structures[J]. China Civil Engineering Journal, 2003, 36(5):105-110.

[3]Mottershead J E, Friswell M I. Model updating in structural dynamics: a survey[J]. Journal of Sound and Vibration, 1993, 167 (2): 347-375.

[4]李輝, 丁樺. 結(jié)構(gòu)動(dòng)力模型修正方法研究進(jìn)展[J]. 力學(xué)進(jìn)展, 2005,35(2): 170-180.

LI Hui, DING Hua. Progress in model updating for structural dynamics[J]. Advances in Mechanics, 2005, 35(2): 170-180.

[5]Friswell M, Mottershead J E. Finite element model updating in structural dynamics[M]. Springer, 1995.

[6]Doebling S W. Minimum-rank optimal update of elemental stiffness parameters for structural damage identification[J]. AIAA journal, 1996, 34 (12): 2615-2621.

[7]Kaouk M, Zimmerman D C, Simmermacher T W. Assessment of damage affecting all structural properties using experimental modal parameters[J]. Journal of Vibration and Acoustics, 2000, 122 (4): 456-463.

[8]侯吉林. 約束子結(jié)構(gòu)模型修正方法[D]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2010.

[9]Hu S L J, Li H, Wang S Q. Cross-model cross-mode method for model updating[J]. Mechanical Systems and Signal Processing, 2007, 21 (4): 1690-1703.

[10]賈輝，李東升. 基于交叉模型交叉模態(tài)法的損傷識(shí)別研究[J]. 土木工程與管理學(xué)報(bào)，2011，28 (3): 133-137.

JIA Hui, LI Dong-sheng. Damage identification inversitigation based on cross-model cross-mode method[J]. Journal of Civil Engineering and Management, 2011, 28(3):133-137.

[11]王俊榮. 海洋平臺(tái)結(jié)構(gòu)物損傷檢測(cè)與模型修正方法研究[D]. 青島：中國(guó)海洋大學(xué), 2009.

[12]王俊榮, 張敏, 李華軍. 一種基于不完備模態(tài)信息的有限元模型修正方法[J]. 工程力學(xué), 2010, 27 (8): 60-65.

WANG Jun-long, ZHANG Min, LI Hua-jun. A finite element model updating approach based on incomplete modal data[J]. Engineering Mechanics, 2010, 27(8): 60-65.

[13]HE Zeng-you, XU Xiao-fei, DENG Sheng-chun. Discovering cluster-based local outliers[J]. Pattern Recognition Letters, 2003, 24 (9-10): 1641-50.

[14]Figueiredo E, Flynn E. Three-story building structure to detect nonlinear effects.http://institute.lanl.gov/ei/damage_id/data/4DOF_Mat_Format.zip. 2009.

[15]Figueiredo E, Park G, Figueiras J, et al. Structural health monitoring algorithm comparisons using standard data sets[R]. LA-14393, Los Alamos, NM (United States): Los Alamos National Laboratory, 2009.

[16]Figueiredo E, Todd M D, Farrar C R, et al. Autoregressive modeling with state-space embedding vectors for damage detection under operational variability[J]. International Journal of Engineering Science, 2010, 48 (10): 822-34.

[17]Miguel R, Hernandez-Garcia, Masri S F, et al. An experimental investigation of change detection in uncertain chain-like systems[J]. Journal of Sound and Vibration, 2010, 329 (12): 2395-2409.

[18]Juang J N. Applied system identification[M]. New Jersey: PTR Prentice Hall Englewood Cliffs Inc, 1994. 175-228.

[19]沃德·海倫，斯蒂芬·拉門(mén)茲，波爾·薩斯. 模態(tài)分析理論與實(shí)驗(yàn)(白化同，郭繼忠譯)[M]. 北京：北京理工大學(xué)出版社，2001.