基于HHT的泄流結(jié)構(gòu)損傷在線監(jiān)測(cè)方法研究

練繼建 ，李成業(yè)，劉 昉 ，馬 斌

(天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室，天津 300072)

1 前言

在水利工程中，泄流結(jié)構(gòu)的工作條件極其復(fù)雜，由于水流的強(qiáng)烈紊動(dòng)，產(chǎn)生的脈動(dòng)壓力作用在結(jié)構(gòu)物上，極有可能造成結(jié)構(gòu)物的強(qiáng)烈振動(dòng)，甚至導(dǎo)致結(jié)構(gòu)的破壞。為此，必須尋找有效的監(jiān)測(cè)方法，以快速、準(zhǔn)確地發(fā)現(xiàn)泄流結(jié)構(gòu)中存在的隱患，有效預(yù)報(bào)其安全狀況。現(xiàn)有的損傷檢測(cè)方法大都事先假設(shè)已經(jīng)獲得了比較詳盡、準(zhǔn)確的有限元模型信息，根據(jù)實(shí)測(cè)結(jié)構(gòu)的模態(tài)參數(shù)進(jìn)行損傷識(shí)別，這往往會(huì)給實(shí)際問題的解決帶來困難。基于信號(hào)處理技術(shù)的損傷識(shí)別是當(dāng)前結(jié)構(gòu)健康監(jiān)測(cè)領(lǐng)域的一種新方法，該方法直接對(duì)結(jié)構(gòu)振動(dòng)信號(hào)進(jìn)行分析，從中提取完備信息并以此對(duì)結(jié)構(gòu)損傷進(jìn)行識(shí)別，具有較好的適用性和敏感性，并有助于克服由噪聲、模態(tài)參數(shù)識(shí)別精度、損傷指標(biāo)敏感性等因素對(duì)損傷結(jié)果的影響。泄流結(jié)構(gòu)損傷過程具有強(qiáng)隨機(jī)性和耦合性特點(diǎn)，加之環(huán)境背景噪聲的干擾作用，使得振動(dòng)響應(yīng)中往往含有非線性成分、奇異點(diǎn)及不規(guī)則的突變部分，應(yīng)用傳統(tǒng)的基于傅立葉變換的信號(hào)處理方法存在一定的局限性。

1998年，N.E.Huang等人提出了 Hilbert-Huang 變換 (HHT)[1，2]，該方法由經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition，EMD)和希爾伯特變換(Hilbert Transform，HT)組成。通過經(jīng)驗(yàn)?zāi)B(tài)分解方法對(duì)結(jié)構(gòu)的振動(dòng)響應(yīng)信號(hào)進(jìn)行預(yù)處理，可以把非平穩(wěn)的原始信號(hào)分解成一系列具有不同特征尺度的平穩(wěn)信號(hào)，把原信號(hào)的能量分配到各個(gè)固有模態(tài)函數(shù)分量上，而且這種分解是自適應(yīng)的，因而更能反映損傷的本質(zhì)信息。目前，HHT法已被廣泛應(yīng)用于結(jié)構(gòu)健康監(jiān)測(cè)與損傷識(shí)別領(lǐng)域[3～6]。

筆者等提出了一種基于HHT的泄流結(jié)構(gòu)損傷在線監(jiān)測(cè)方法，利用信號(hào)能量譜構(gòu)造能量分布向量，進(jìn)而定義結(jié)構(gòu)異常指標(biāo)，通過計(jì)算相鄰時(shí)間段內(nèi)的異常指標(biāo)，實(shí)現(xiàn)結(jié)構(gòu)損傷的在線監(jiān)測(cè)和識(shí)別。

2 HHT算法的基本原理

N.E.Huang等人定義滿足下面兩個(gè)條件的函數(shù)為固有模態(tài)函數(shù):a.整個(gè)數(shù)據(jù)序列中，極值點(diǎn)的數(shù)量與過零點(diǎn)數(shù)量相等或至多相差1;b.信號(hào)關(guān)于時(shí)間軸局部對(duì)稱。

實(shí)測(cè)泄流結(jié)構(gòu)振動(dòng)信號(hào)通常不滿足IMF條件，N.E.Huang等人提出經(jīng)驗(yàn)?zāi)B(tài)分解法。對(duì)任意原始信號(hào)x(t)，EMD方法的計(jì)算過程是:找出x(t)所有的極大值點(diǎn)，應(yīng)用三次樣條函數(shù)擬合成原始信號(hào)的上包絡(luò)線;找出x(t)所有的極小值點(diǎn)，應(yīng)用三次樣條函數(shù)擬合成原始信號(hào)的下包絡(luò)線;計(jì)算上下包絡(luò)線的均值m1(t)，將原始信號(hào)x(t)減去m1(t)得到一個(gè)新的信號(hào)h1(t)。通常情況下，h1(t)并不滿足一個(gè)IMF分量的條件，此時(shí)將h1(t)看成原信號(hào)，重復(fù)上述步驟，得到第k次篩選的數(shù)據(jù)h1k(t)=h1(k-1)(t)-m1k(t)。此時(shí)需要有一個(gè)終止準(zhǔn)則來判斷h1k(t)是不是一個(gè)IMF分量，否則過多地重復(fù)處理過程，會(huì)導(dǎo)致IMF分量失去意義。這個(gè)準(zhǔn)則可以用兩個(gè)連續(xù)的處理結(jié)果之間的標(biāo)準(zhǔn)差Sd表示，即:

實(shí)踐證明，當(dāng)Sd的值在0.2～0.3時(shí)，既能保證基本模態(tài)函數(shù)的線性和穩(wěn)定性，又能使其具有相應(yīng)的物理意義。當(dāng)h1k(t)滿足Sd值的要求，則h1k(t)為原信號(hào)的第一階IMF分量，記為c1(t)，將原始信號(hào)x(t)減去c1(t)，重復(fù)上述過程，直到最后一個(gè)殘余信號(hào)rn(t)不可再分解為止。最后，原始信號(hào)x(t)可以表示成:

文獻(xiàn)[7]應(yīng)用EMD算法分析爆破振動(dòng)信號(hào)，驗(yàn)證了該方法的完備性和正交性。應(yīng)用EMD分解可得到信號(hào)x(t)的多個(gè)IMF分量的組合，對(duì)ci(t)作Hilbert變換:

式(3)中，PV代表柯西主值(cauchy principal value)，因此構(gòu)造解析信號(hào)z(t):

在相位函數(shù)θi(t)的基礎(chǔ)上定義瞬時(shí)頻率為:

瞬時(shí)頻率是時(shí)間t的函數(shù)，它表示某一時(shí)刻信號(hào)能量在頻域分布的量度。

3 基于HHT的結(jié)構(gòu)損傷在線監(jiān)測(cè)指標(biāo)

3.1 信號(hào)的能量分布向量

對(duì)原始信號(hào)的每個(gè)IMF分量進(jìn)行Hilbert變換后，可得Hilbert譜的表達(dá)形式:

式(6)中，i=1，2，…，n表示 IMF 分量的個(gè)數(shù)，省略了殘余分量rn(t);Re表示取實(shí)部。如果Hω，()t對(duì)時(shí)間積分，就可以得到Hilbert邊際譜:

邊際譜表達(dá)了每個(gè)頻率在全局上的幅度，它代表了統(tǒng)計(jì)意義上的全部累加幅度。如果振幅的平方對(duì)時(shí)間積分，就可以得到Hilbert能量譜:

Hilbert能量譜表達(dá)了每個(gè)頻率在整個(gè)時(shí)間長度內(nèi)所積累的能量，將信號(hào)的頻率區(qū)間[ωmin，ωmax]分成n個(gè)頻帶 ωi(i=1，2，3，…，n)，則各個(gè)頻帶上的能量可表示為E(ωi)，進(jìn)而定義能量分布向量ηi為:

3.2 結(jié)構(gòu)損傷異常指標(biāo)

結(jié)構(gòu)發(fā)生損傷引起固有頻率的變化，其振動(dòng)響應(yīng)能量在頻域以及各階模態(tài)上的分布也會(huì)發(fā)生改變，且改變量和損傷位置及程度有關(guān)。利用信號(hào)Hilbert能量譜構(gòu)造結(jié)構(gòu)能量分布向量ηi，再進(jìn)一步用它來構(gòu)造結(jié)構(gòu)異常指數(shù)，對(duì)損傷是否發(fā)生進(jìn)行識(shí)別。設(shè)、分別為無損傷和未知情況下的結(jié)構(gòu)能量分布向量，i=1，2，3，…，n。定義和的互相關(guān)系數(shù)為 Φ()，構(gòu)造損傷指標(biāo)N[8]:f

當(dāng)Nf等于0或者非常小，表明響應(yīng)中模態(tài)成分沒有發(fā)生改變或僅發(fā)生微小變化，無損傷發(fā)生;當(dāng)結(jié)構(gòu)出現(xiàn)損傷，能量在各階模態(tài)上的分布會(huì)發(fā)生明顯改變，模態(tài)成分匹配度降低，Nf將遠(yuǎn)大于0。

結(jié)構(gòu)出現(xiàn)損傷，能量在頻域上的分布必然發(fā)生變化，導(dǎo)致能量分布向量ηi峰值所對(duì)應(yīng)的頻率發(fā)生變化，構(gòu)建損傷指標(biāo)Rf:

4 水電站導(dǎo)墻結(jié)構(gòu)在線監(jiān)測(cè)數(shù)值仿真試驗(yàn)

文章以某工程導(dǎo)墻結(jié)構(gòu)為例，導(dǎo)墻截面尺寸如圖1所示。實(shí)際工程導(dǎo)墻順?biāo)鞣较蛎?4 m設(shè)一條結(jié)構(gòu)縫，因此取24 m長的導(dǎo)墻段建立有限元模型，模型中地基模擬深度為54 m，兩側(cè)寬度各取30 m，地基采用全約束，將與水流方向垂直的兩個(gè)導(dǎo)墻面加上法向約束，沿高程自上而下依次布置1#～6#測(cè)點(diǎn)，有限元模型及測(cè)點(diǎn)布置如圖2所示。混凝土材料密度取2 400 kg/m3，動(dòng)彈性模量取35 GPa，泊松比取0.167，在計(jì)算時(shí)將地基材料屬性與導(dǎo)墻取為一致。

圖1 導(dǎo)墻斷面圖(單位:m)Fig.1 The section of the guide wall(unit:m)

為驗(yàn)證指標(biāo)Nf和Rf的穩(wěn)定性，根據(jù)泄流結(jié)構(gòu)激勵(lì)形式的特點(diǎn)，設(shè)計(jì)激勵(lì)位置相同但大小不同的白噪聲點(diǎn)激勵(lì)、激勵(lì)位置和大小均不同的白噪聲點(diǎn)激勵(lì)以及激勵(lì)位置和大小均不同的白噪聲面激勵(lì)三種工況，如表1中工況1～工況3，激勵(lì)值及其改變量均由模型試驗(yàn)[9]確定。由模型試驗(yàn)和有限元計(jì)算結(jié)果可知，導(dǎo)墻在截面A、高程260 m處的動(dòng)應(yīng)力值最大，如圖3所示，該處出現(xiàn)裂縫的可能性也最大，為驗(yàn)證指標(biāo)Nf和Rf對(duì)不同損傷程度的敏感性，文章設(shè)計(jì)了截面A處、裂縫深度分別為導(dǎo)墻厚度1/8、1/4和1/2三種程度的損傷，見表1工況4～工況6。

圖2 導(dǎo)墻有限元模型及測(cè)點(diǎn)布置圖Fig.2 The finite element model and measured points locations

圖3 模型試驗(yàn)最不利工況下導(dǎo)墻動(dòng)應(yīng)力Fig.3 Guide wall dynamic stress for the most unfavorable condition of test

表1 數(shù)值試驗(yàn)工況表Table1 The conditions of numerical simulation test

續(xù)表

4.1 無損傷情況下各指標(biāo)穩(wěn)定性分析

對(duì)無損工況1～工況3條件下各測(cè)點(diǎn)的位移響應(yīng)進(jìn)行EMD分解，應(yīng)用式(6)～(8)計(jì)算信號(hào)的能量譜，應(yīng)用式(9)即可得到能量分布向量。通過計(jì)算分析，各工況下所有測(cè)點(diǎn)激勵(lì)變化(t=25 s)前后ηi的峰值位置以及在各頻率段上的分布規(guī)律基本相同，由于篇幅限制，文中僅列出1工況下3#、4#測(cè)點(diǎn)激勵(lì)變化前后的ηi分布圖，如圖4所示。根據(jù)式(10)可計(jì)算出各測(cè)點(diǎn)不同工況下的Nf值，見表2。從各測(cè)點(diǎn)的能量譜圖可以看出，信號(hào)能量主要在0～25 Hz范圍內(nèi)，找出t時(shí)刻前后能量分布向量ηi各峰值所對(duì)應(yīng)的頻率和，由式(11)計(jì)算出各測(cè)點(diǎn)不同工況下的Rf值，見表2。

圖4 工況1改變激勵(lì)前后3#、4#測(cè)點(diǎn)的能量分布向量圖Fig.4 Energy distribution vector of point 3#，4#before and after excitation changed for condition 1

表2 工況1～工況3不同測(cè)點(diǎn)的Nf和Rf指標(biāo)Table 2 The Nfand Rfof every measured point from condition 1 to condition 3

從表2中可以看出，工況1條件下各測(cè)點(diǎn)的Nf值均接近于0，而工況2、工況3激勵(lì)位置和大小都發(fā)生改變的情況下，各測(cè)點(diǎn)的Nf指標(biāo)均有所增加，但由模型試驗(yàn)確定的最不利荷載條件下(工況3)，各測(cè)點(diǎn)指標(biāo)Nf的最大值為0.089 2，因此本例定義超過此最大值10%的數(shù)值作為損傷預(yù)警指數(shù)，即，當(dāng)Nf＞0.1時(shí)，即認(rèn)為結(jié)構(gòu)出現(xiàn)異常。與Nf值受激勵(lì)位置和激勵(lì)條件變化的影響不同，各工況下各測(cè)點(diǎn)信號(hào)的Rf指標(biāo)均為0，不因外界條件的變化而變化，只與結(jié)構(gòu)自身的特性相關(guān)，當(dāng)Rf＞0時(shí)，即認(rèn)為結(jié)構(gòu)出現(xiàn)異常。在實(shí)際泄流結(jié)構(gòu)損傷監(jiān)測(cè)過程中，可綜合分析Nf和Rf兩項(xiàng)指標(biāo)值，作為判斷結(jié)構(gòu)出現(xiàn)損傷的依據(jù)。

4.2 泄流結(jié)構(gòu)損傷的判定

4.2.1 在線監(jiān)測(cè)信號(hào)IMF分量的奇異點(diǎn)和瞬時(shí)頻率

文獻(xiàn)[10]和文獻(xiàn)[11]指出結(jié)構(gòu)剛度發(fā)生突變，其振動(dòng)響應(yīng)EMD結(jié)果中的高階固有模態(tài)函數(shù)在時(shí)域上將出現(xiàn)奇異點(diǎn)。泄流結(jié)構(gòu)出現(xiàn)損傷必然導(dǎo)致剛度的降低，結(jié)構(gòu)在線監(jiān)測(cè)信號(hào)表現(xiàn)出一個(gè)小的不連續(xù)性，信號(hào)出現(xiàn)了奇異點(diǎn)，直接觀察振動(dòng)信號(hào)很難發(fā)現(xiàn)這一不連續(xù)點(diǎn)，而對(duì)結(jié)構(gòu)響應(yīng)進(jìn)行EMD分解所得的高頻IMF信號(hào)可以檢測(cè)出該信號(hào)的奇異點(diǎn)及其發(fā)生的時(shí)刻。對(duì)4、5、6三個(gè)損傷工況下各測(cè)點(diǎn)響應(yīng)進(jìn)行EMD分解，分解結(jié)果的第一階IMF分量在25 s時(shí)刻均出現(xiàn)了明顯的突變，任取工況4和工況6測(cè)點(diǎn)2的響應(yīng)進(jìn)行分析，如圖5所示，表明t=25 s時(shí)刻結(jié)構(gòu)剛度可能發(fā)生突變。

圖5 工況4和工況6測(cè)點(diǎn)2響應(yīng)EMD分解的第一階IMF分量Fig.5 The first IMF component of point 2 for condition 4 and condition 6 respectively

外界載荷和工況的突然變化也會(huì)使EMD的分解結(jié)果出現(xiàn)奇異點(diǎn)，因此信號(hào)出現(xiàn)奇異點(diǎn)不足以判斷結(jié)構(gòu)出現(xiàn)損傷。信號(hào)IMF分量的瞬時(shí)頻率是時(shí)間的函數(shù)，反映結(jié)構(gòu)固有頻率隨時(shí)間的變化規(guī)律，對(duì)于健康的結(jié)構(gòu)其瞬時(shí)頻率是不變的，瞬時(shí)頻率的衰減反映了結(jié)構(gòu)已受到某種損傷。通過分析，損傷工況下各測(cè)點(diǎn)IMF分量的瞬時(shí)頻率在損傷前后均有一定的下降。任取工況6測(cè)點(diǎn)2分解后第四階IMF分量的瞬時(shí)頻率，如圖6所示。從圖6中可以看出，結(jié)構(gòu)響應(yīng)的瞬時(shí)頻率在第25 s時(shí)突然下降，瞬時(shí)頻率的平均值也由4.12 Hz下降到2.77 Hz，這表明此時(shí)刻結(jié)構(gòu)發(fā)生某種突然的損傷。

圖6 工況6測(cè)點(diǎn)2響應(yīng)EMD分解第四階IMF分量的瞬時(shí)頻率Fig.6 The instantaneous frequency of the fourth IMF component from point 2 for condition 6

4.2.2 異常指標(biāo)對(duì)損傷的敏感性分析

計(jì)算工況4～工況6各測(cè)點(diǎn)的能量分布向量ηi，通過分析可知，各工況下所有測(cè)點(diǎn)損傷前后ηi值的峰值位置及其在各頻率段上的分布均發(fā)生了明顯變化。由于篇幅限制，文中僅列出5工況下3#、4#測(cè)點(diǎn)損傷前后的ηi分布圖，如圖7所示。

按照3.2節(jié)中的方法，分別計(jì)算工況4～工況6各測(cè)點(diǎn)的Nf和Rf指標(biāo)，見表3。從表3可以看出，損傷特征指標(biāo)Nf和Rf對(duì)各種不同損傷模式均較為敏感，可以判斷結(jié)構(gòu)損傷的發(fā)生。當(dāng)結(jié)構(gòu)損傷位置一定時(shí)，隨著損傷程度的增加，各測(cè)點(diǎn)指標(biāo)值相應(yīng)變大，說明根據(jù)損傷指標(biāo)值的大小可以定性地判斷結(jié)構(gòu)的損傷程度。但是相同工況下，各測(cè)點(diǎn)的指標(biāo)值沒有明顯的規(guī)律性，根據(jù)能量定義的特征指標(biāo)Nf和Rf無法判斷損傷發(fā)生的位置，這是因?yàn)榇笮退そY(jié)構(gòu)不同部位的振動(dòng)情況復(fù)雜，對(duì)損傷的敏感程度也不相同。

表3 工況4～工況6不同測(cè)點(diǎn)的Nf和Rf指標(biāo)Table 3 The Nfand Rfof every measured point from condition 4 to condition 6

4.3 指標(biāo)抗噪性研究

信號(hào)采集過程中不可避免地受到噪聲的影響，為了分析噪聲對(duì)損傷指標(biāo)的影響，文章在各工況下1#測(cè)點(diǎn)響應(yīng)信號(hào)中分別加入標(biāo)準(zhǔn)差為10%、20%、30%、40%、50%、60%、70%、80%倍信號(hào)標(biāo)準(zhǔn)差的白噪聲，定義為噪聲水平1～8。EMD分解本身就是一個(gè)由高頻到低頻的濾波過程，由文獻(xiàn)[12]可知，白噪聲EMD分解后能量主要集中在第一階分量(約占50% ～70%)，應(yīng)用信號(hào)能量譜計(jì)算損傷指標(biāo)時(shí)，除去第一階IMF分量以及最后幾階低頻干擾分量，計(jì)算結(jié)果見表4。

圖7 工況5損傷前后3#、4#測(cè)點(diǎn)的能量分布向量圖Fig.7 Energy distribution vector of point 3#，4#before and after damage appearance for condition 5

從表4可以看出，在各種噪聲水平下，無損傷工況的Nf指標(biāo)均小于預(yù)先設(shè)定的預(yù)警指數(shù)0.1，與表2中不含噪聲情況下的Nf值相比沒有明顯的規(guī)律性;而損傷工況下，隨著噪聲水平的不斷提高，信號(hào)Nf指標(biāo)逐漸減小，但減小的最大幅度(25.52%)遠(yuǎn)低于噪聲的增幅。信號(hào)各工況下的Rf指標(biāo)受噪聲的影響較小，基本與無噪情況下的計(jì)算結(jié)果相同。因此，在一定噪聲水平的情況下，文中方法仍可以判斷結(jié)構(gòu)的損傷。

表4 各工況1#測(cè)點(diǎn)含噪信號(hào)的Nf和Rf指標(biāo)Table 4 The Nfand Rfof point 1#responses disturbed by different noise levels for every condition

基于HHT能量分布指標(biāo)的泄流結(jié)構(gòu)在線監(jiān)測(cè)方法無需知道結(jié)構(gòu)的模態(tài)參數(shù)等其他特征信息，也無需激勵(lì)信息，只需對(duì)任一位置處響應(yīng)進(jìn)行監(jiān)測(cè)即可實(shí)現(xiàn)，能夠滿足結(jié)構(gòu)健康監(jiān)測(cè)的動(dòng)態(tài)、實(shí)時(shí)性要求。但對(duì)于損傷位置和損傷程度的進(jìn)一步精確辨識(shí)，還需進(jìn)一步研究。

5 結(jié)語

文章利用HHT法在描述信號(hào)瞬時(shí)特性和信號(hào)處理方面獨(dú)特的優(yōu)勢(shì)，提出了一種基于HHT的泄流結(jié)構(gòu)損傷在線監(jiān)測(cè)方法。首先應(yīng)用HHT法計(jì)算響應(yīng)Hilbert能量譜，構(gòu)造能量分布向量ηi，根據(jù)結(jié)構(gòu)發(fā)生損傷時(shí)其振動(dòng)能量在頻域和各階模態(tài)上的分布也會(huì)隨之改變的特點(diǎn)，定義結(jié)構(gòu)異常指標(biāo)Nf和Rf，通過計(jì)算相鄰時(shí)間段內(nèi)的異常指標(biāo)，實(shí)現(xiàn)結(jié)構(gòu)損傷在線監(jiān)測(cè)。

應(yīng)用文中方法對(duì)一導(dǎo)墻結(jié)構(gòu)進(jìn)行數(shù)值仿真試驗(yàn)，模擬了不同激勵(lì)情況下的三種無損傷工況以及相同損傷位置不同損傷程度下的三種有損工況。分別進(jìn)行了無損傷情況下各指標(biāo)的穩(wěn)定性分析、損傷情況下各指標(biāo)的敏感性分析以及各工況下指標(biāo)的抗噪性分析，結(jié)果表明，文中的異常指標(biāo)能夠較精確地判斷結(jié)構(gòu)的損傷，同時(shí)具有一定的抗噪能力，其中指標(biāo)Rf的穩(wěn)定性、敏感性以及抗噪性更為突出。但是應(yīng)用文中方法對(duì)損傷位置和損傷程度進(jìn)行精確辨識(shí)，還需要進(jìn)一步的深入研究。

[1]Huang N E，Shen Z.The empirical mode decomposition and Hilbert spectrum for nonlinear and non-stationary time series analysis[J].Proceedings of the Royal Society of London，Series A.1998，454:903-995.

[2]Huang N E，Steven R L.A new view of nonlinear water waves:the Hilbert Spectrum[J].Annual Review of Fluid Mesh，1999，31(2):417-457.

[3]Huang N E，Huang K.基于希爾伯特-黃變換的鐵路橋梁結(jié)構(gòu)健康監(jiān)測(cè)[J].中國鐵道科學(xué)，2006，27(l):l-7.

[4]舒忠平.基于希爾波特-黃變換的結(jié)構(gòu)損傷檢測(cè)方法研究[D].西安:西北工業(yè)大學(xué)，2006.

[5]石志曉，李 昕，周 晶.損傷檢測(cè)的經(jīng)驗(yàn)?zāi)B(tài)分解法[J].大連理工大學(xué)學(xué)報(bào)，2005，45(3):401-404.

[6]Qu W Z，Chen M X，Zeng Y L，et al.Structural damage detection using the Hilbert-Huang transforms method[C]//Li A，Sih G C，Nied H F，eds.Proceedings of International Conference on Health Monitoring of Structure，Materials and Environment.Nanjing:Southeast University Press，2007:293-298.

[7]張義平，李夕兵，左宇軍，等.爆破振動(dòng)信號(hào)的HHT分析與應(yīng)用[M].北京:冶金工業(yè)出版社，2008.

[8]李傳習(xí)，陳富強(qiáng).基于HHT的結(jié)構(gòu)損傷特征量與異常診斷[J].長沙理工大學(xué)學(xué)報(bào)(自然科學(xué)版)，2007，4(3):29-33.

[9]練繼建.向家壩水電站消力池底板與導(dǎo)(隔)墻結(jié)構(gòu)水彈性模型試驗(yàn)研究[R].天津:天津大學(xué)水利水電工程系，2010.

[10]Xu Y L，Chen J.Structural damage detection using empirical mode decomposition:experimental investigation[J].Journal of Engineering Mechanics，2004，130(11):1279-1288.

[11]Yang J N，Lei Y，Lin S，et al.Hilbert-Huang based approach for structural damage detection[J].Journal of Engineering Mechanics，2004，130(1):85-95.

[12]Zhao Huawu，Norden E H.A study of the characteristics of white noise using the empirical mode decomposition method[J].The Royal Society，Proceedings Royal Society of London，A(2004)460:1597-1611.