王 慧，王 樂，田潤澤

(1. 長安大學(xué)理學(xué)院，西安 710064；2. 西北工業(yè)大學(xué)航空學(xué)院，西安 710072)

基于振動的結(jié)構(gòu)損傷檢測因其測試成本較低且容易實現(xiàn)在線監(jiān)測，近些年來一直受到研究者及工程技術(shù)人員的廣泛關(guān)注。根據(jù)損傷檢測過程中是否需要結(jié)構(gòu)的理論模型，基于振動的結(jié)構(gòu)損傷檢測方法一般可分為基于模型的方法及不基于模型的方法兩類[1]。基于模型的方法通常利用了修正后的理論模型，檢測精度較好，但其涉及的模型修正過程一般需要較大的計算工作量，目前國內(nèi)外學(xué)者在提升損傷檢測中模型修正的計算效率方面作了大量的工作，如智能算法[2]、云計算[3]、隔離損傷[4]等技術(shù)。相比于基于模型的方法而言，不基于模型的方法具有兩大優(yōu)點：1)不需要完好結(jié)構(gòu)準(zhǔn)確的理論模型；2)檢測流程一般比較簡單。因此，在結(jié)構(gòu)健康監(jiān)測領(lǐng)域，不基于模型的方法具有一定的優(yōu)越性。不基于模型的方法常用的損傷特征有：1)模態(tài)參數(shù)(包括模態(tài)振型及其導(dǎo)出參數(shù))[5?7]；2)頻域響應(yīng)(包括頻譜[8?10]、頻響函數(shù)[11?12]、傳遞率函數(shù)[13?14]、相干函數(shù)[15]等)；3)時域響應(yīng)(包括時域歷程、相關(guān)函數(shù)等)[1,16?18]?？紤]到結(jié)構(gòu)健康監(jiān)測一般需要在結(jié)構(gòu)工作狀態(tài)下開展，而在結(jié)構(gòu)工作狀態(tài)下一般不容易對結(jié)構(gòu)人為施加激勵或進行激勵測量，因此僅利用結(jié)構(gòu)工作狀態(tài)下振動響應(yīng)的結(jié)構(gòu)健康監(jiān)測方法顯得尤為重要，例如：顏王吉等[13?14]首先介紹了傳遞率函數(shù)的特點及分類，進而在闡明傳遞率函數(shù)與頻響函數(shù)的內(nèi)在聯(lián)系的基礎(chǔ)之上，綜述了基于傳遞率函數(shù)的結(jié)構(gòu)損傷檢測方法的研究現(xiàn)狀及其發(fā)展方向；陳闖等[16]以測量的時域加速度響應(yīng)或應(yīng)變響應(yīng)為原始信號，基于馬氏距離累積量和經(jīng)驗?zāi)７纸馓岢隽艘环N結(jié)構(gòu)損傷識別的兩步法，并通過簡支梁數(shù)值模擬和工字鋼模型試驗驗證了方法的可行性及有效性。

利用帶通白噪聲激勵下的時域響應(yīng)互/自相關(guān)分析，本文作者課題組也曾提出了兩種不基于模型的結(jié)構(gòu)損傷檢測方法，即基于內(nèi)積向量(inner product vector，IPV)的結(jié)構(gòu)損傷檢測方法[1,19?20]及基于互相關(guān)函數(shù)幅值向量(cross correlation function amplitude vector，CorV)的結(jié)構(gòu)損傷檢測方法[21]。針對框架結(jié)構(gòu)、復(fù)合材料梁結(jié)構(gòu)以及航空加筋壁板結(jié)構(gòu)的損傷檢測仿真算例及試驗驗證研究表明，這兩種方法僅利用結(jié)構(gòu)在帶通白噪聲激勵下的振動時域響應(yīng)，并經(jīng)過簡單的數(shù)學(xué)運算，就可以檢測出相關(guān)結(jié)構(gòu)的損傷位置，非常適用于結(jié)構(gòu)在線損傷檢測。然而，上述所提出的兩種方法均只采用了互/自相關(guān)函數(shù)中的部分信息。為了從結(jié)構(gòu)時域響應(yīng)互/自相關(guān)函數(shù)中提取更多的結(jié)構(gòu)健康狀態(tài)信息，通過利用不同的結(jié)構(gòu)響應(yīng)組合，本文將IPV 及CorV 擴展到了多種結(jié)構(gòu)特征向量(structural characteristic vector，SCV)，進而建立多種結(jié)構(gòu)損傷指標(biāo)，并采用數(shù)據(jù)融合算法，解決單一結(jié)構(gòu)損傷指標(biāo)在檢測準(zhǔn)確度上的不足，以實現(xiàn)對結(jié)構(gòu)損傷位置的準(zhǔn)確定位。

1 IPV 及CorV

大多結(jié)構(gòu)工作狀態(tài)下所受的環(huán)境激勵都可以假設(shè)為典型隨機激勵(即帶通白噪聲激勵)，IPV及CorV 就是利用帶通白噪聲激勵下結(jié)構(gòu)時域響應(yīng)的互/自相關(guān)函數(shù)來定義的。

1.1 IPV

1.2 CorV

利用結(jié)構(gòu)響應(yīng)互相關(guān)函數(shù)時間延遲為零的值可以定義IPV，而CorV 是利用結(jié)構(gòu)響應(yīng)互相關(guān)函數(shù)的幅值(即最大絕對值)來定義的，如下：

同樣可以證明，在帶通白噪聲激勵下，CorV僅取決于結(jié)構(gòu)頻響函數(shù)矩陣[21]。類似地，結(jié)構(gòu)局部物理參數(shù)的變化會體現(xiàn)在結(jié)構(gòu)頻響函數(shù)矩陣的變化上，進而CorV 也會出現(xiàn)相應(yīng)的變化，也就是說，CorV 也可以當(dāng)作SCV 來構(gòu)建結(jié)構(gòu)損傷指標(biāo)。

1.3 基于IPV 或CorV 的結(jié)構(gòu)損傷檢測方法

根據(jù)IPV 及CorV 的定義，在其計算過程中，可以自動剔除具有如下特性的測量噪聲[1]：1)測量噪聲均值為零；2)某測點的測量噪聲與該測點的真實響應(yīng)無關(guān)；3)不同測點之間的測量噪聲相互獨立。然而，上述對測量噪聲特性描述的1)及3)在工程實際中很難滿足，即測量噪聲的均值可能非零且不同測點的測量噪聲可能具有一定的關(guān)系。因此，我們可以給出一種更加符合工程實際的測量噪聲假設(shè)，即不同測點的測量噪聲具有相同的統(tǒng)計特性(包括均值及方差)。在這種更加符合工程實際的測量噪聲假設(shè)下，可以證明，IPV或CorV 中每一個元素都額外增加了一個與噪聲統(tǒng)計特性相關(guān)的且各個測點都相等的常數(shù)?？紤]到是利用IPV 或CorV 的突變來檢測結(jié)構(gòu)損傷位置，因此給IPV 或CorV 每一個元素都增加一個相等的常數(shù)并不會影響IPV 或CorV 突變的位置。也就是說，滿足上述測量噪聲假設(shè)的情況下，仍可以直接利用IPV 或CorV 的突變來檢測結(jié)構(gòu)損傷位置。

然而，實際振動測試過程中，測量噪聲并不能準(zhǔn)確滿足上述對測量噪聲的假設(shè)，即IPV 或CorV 不可避免地會受到測量噪聲的影響。同時，結(jié)構(gòu)的微小損傷導(dǎo)致的IPV 或CorV 的突變可能會很平緩，難以直接從IPV 或CorV 的突變中判別。因此，為了降低測量噪聲對IPV 或CorV 的影響并有效區(qū)分IPV 或CorV 因局部損傷而導(dǎo)致的突變，可以利用結(jié)構(gòu)損傷前后IPV 或CorV 對應(yīng)位置元素的差值來定義損傷指標(biāo)：

2 基于相關(guān)性分析的SCV

為了從結(jié)構(gòu)時域響應(yīng)互/自相關(guān)函數(shù)中提取更多的結(jié)構(gòu)健康狀態(tài)信息，通過利用不同的結(jié)構(gòu)響應(yīng)組合，可將IPV 及CorV 擴展到多種SCV。

2.1 相同時刻采集的振動響應(yīng)構(gòu)建的互/自相關(guān)函數(shù)

2.2 不同時刻采集的振動響應(yīng)構(gòu)建的互相關(guān)函數(shù)

3 基于數(shù)據(jù)融合的結(jié)構(gòu)損傷檢測

4 試驗驗證

圖1 基于數(shù)據(jù)融合及時域響應(yīng)相關(guān)性分析的損傷檢測流程Fig. 1 Damage detection procedure based on correlation functions of time domain vibration responses and data fusion

試驗驗證選取我們之前損傷檢測研究所采用的8 層剪切框架結(jié)構(gòu)模型[19]，如圖2 所示，每層用一厚度為13 mm 的口字型隔框來模擬樓板(隔框外尺寸為250 mm×250 mm，隔框內(nèi)尺寸為225 mm×225 mm，層間4 個立柱分別由3 根鋼片(尺寸為139 mm×27 mm×1 mm)疊加而成，立柱通過耳片與隔框用螺栓連接。根據(jù)模型特征，其隔框剛度遠大于層間支柱，則可假定其每層的質(zhì)量都集中在層間隔框上，而剛度則集中在層間支柱上，同時考慮到層間支柱在其寬度方向的彎曲剛度遠大于厚度方向的彎曲剛度，所以只需考慮模型在垂直于鋼片寬度方向的運動，即激勵及響應(yīng)測量的方向均為垂直于鋼片寬度方向。對于該實驗?zāi)Ｐ?，可以通過拆卸層間鋼片的方法模擬結(jié)構(gòu)的局部損傷，即去掉兩層之間的1 個鋼片，該局部損傷的程度大致為1/12。文獻[19]模擬了三個不同位置三種不同程度的損傷，分別為D31、D32、D33、D51、D52、D53、D71、D72 及D73(Dij表示卸掉第i?1 層與第i層之間的j個鋼片)。文獻[19]中的方法可以準(zhǔn)確地檢測出除D71 之外的所有結(jié)構(gòu)損傷狀態(tài)，然而在D71 上出現(xiàn)了漏檢，究其原因，對于此類懸臂類結(jié)構(gòu)，相比于其他損傷狀態(tài)，D71 所模擬的結(jié)構(gòu)局部損傷(損傷程度較小且損傷位置靠近結(jié)構(gòu)的自由端)對整體結(jié)構(gòu)的動力學(xué)特征影響很小，因此本文重點研究對D71 損傷狀態(tài)的檢測，即卸掉第6 層與第7 層之間的一個鋼片。

圖2 剪切框架模型的損傷檢測試驗布置[19]Fig. 2 Shear frame structure and test setup for damage detection[19]

在試驗過程中，使用Agilent 33220A 信號發(fā)生器產(chǎn)生0 Hz～5 Hz 的帶通白噪聲激勵，并經(jīng)MB Dynamics SS250VCF 功率放大器放大后，輸入給MB Dynamics Modal 50A 激振器，在框架的第一層進行激勵，并在框架的8 個隔框上分別布置8 個PCB Piezotronics 333B30 加速度計，最后利用LMS SCADAS 以及LMS Test.lab 信號采集系統(tǒng)以512 Hz 的采樣頻率采集框架結(jié)構(gòu)8 個隔框處的加速度響應(yīng)，響應(yīng)采集時長為32 s。

與文獻[19]相同，分別對完好結(jié)構(gòu)及損傷結(jié)構(gòu)進行激勵，并利用原始損傷指標(biāo)的一階差分D′構(gòu)建損傷指標(biāo)，圖3～圖6 分別給出了利用不同損傷指標(biāo)的檢測結(jié)果，圖7 給出了最終的數(shù)據(jù)融合結(jié)果?？梢钥闯觯悍謩e利用幾組SCV 并進行數(shù)據(jù)融合(見圖3～圖6)，雖然也可以確定損傷狀態(tài)D71 的損傷位置，但出現(xiàn)了明顯的誤檢(即其他未損傷位置也出現(xiàn)了損傷概率)；同時利用所有SCV并進行數(shù)據(jù)融合(見圖7)，可以準(zhǔn)確定位D71 的損傷位置，且沒有任何誤檢出現(xiàn)(即第6 層與第7 層之間的損傷概率非常接近1，其余位置均幾乎為0)。針對框架結(jié)構(gòu)微小損傷的檢測結(jié)果表明，利用所有SCV 作為原始信息進行數(shù)據(jù)融合，可以顯著提高檢測的準(zhǔn)確程度。

圖3 利用 RIPV 和 RIPV,Diag的損傷檢測結(jié)果Fig. 3 Damage detection results using RIPV and RIPV,Diag

圖4 利用 RCorV 和 RCorV,Diag的損傷檢測結(jié)果Fig. 4 Damage detection results using RCorV and RCorV,Diag

圖5 利用 IPV 和 IPV,Diag的損傷檢測結(jié)果Fig. 5 Damage detection results using IPV and IPV,Diag

圖6 利用CorV 和 CorV,Diag的損傷檢測結(jié)果Fig. 6 Damage detection results using CorV and CorV,Diag

圖7 利用所有SCV 進行數(shù)據(jù)融合的損傷檢測結(jié)果Fig. 7 Damage detection results by data fusion using all the SCVs

5 結(jié)論

利用結(jié)構(gòu)時域響應(yīng)的相關(guān)性分析，提出了基于時域響應(yīng)相關(guān)性分析及數(shù)據(jù)融合的結(jié)構(gòu)損傷檢測方法，通過框架結(jié)構(gòu)損傷檢測的試驗研究驗證了方法的有效性。本文的主要工作及結(jié)論如下：

(1) 根據(jù)IPV 及CorV 的定義及特性，利用時域響應(yīng)定義了兩個互/自相關(guān)函數(shù)矩陣，即：利用相同時刻各個測點時域響應(yīng)的互/自相關(guān)函數(shù)矩陣M(τ)，以及利用不同時刻各個測點時域響應(yīng)的互相關(guān)函數(shù)矩陣(τ)。

(2) 通過設(shè)置互/自相關(guān)函數(shù)矩陣的時間延遲為零獲得了兩個IPV 矩陣MIPV及IPV，通過對互/自相關(guān)函數(shù)矩陣的元素求最大絕對值獲得了兩個CorV 矩陣MCorV及CorV，并根據(jù)IPV 矩陣及CorV矩陣的特性，定義了八種SCV，即RIPV、RIPV,Diag、RCorV、RCorV,Diag、IPV、IPV,Diag、CorV及CorV,Diag。

(3) 利用八種SCV 并結(jié)合D-S 證據(jù)理論，提出了基于數(shù)據(jù)融合的結(jié)構(gòu)損傷檢測方法，并通過框架結(jié)構(gòu)損傷檢測試驗研究表明，本文方法可以對框架結(jié)構(gòu)上的微小損傷進行定位。

考慮到本文方法進行損傷定位的基礎(chǔ)是利用損傷指標(biāo)的局部最大值，其不足之處就是需要各個可能出現(xiàn)損傷位置的測試數(shù)據(jù)，因此當(dāng)缺失損傷位置附近測點的測量數(shù)據(jù)時，本文方法容易出現(xiàn)漏檢現(xiàn)象，這也是本文方法面向具體工程應(yīng)用時需要重點考慮的研究內(nèi)容。