朱天煦 臧朝平

高溫結(jié)構(gòu)激光掃描模態(tài)測(cè)試的降噪與參數(shù)識(shí)別

朱天煦 臧朝平

(南京航空航天大學(xué)能源與動(dòng)力學(xué)院，南京，210016)

針對(duì)高溫結(jié)構(gòu)激光掃描模態(tài)測(cè)試中存在的噪聲問題，提出了一種基于模態(tài)峰值漢克爾奇異值分解的降噪處理與模態(tài)參數(shù)識(shí)別方法。首先，將測(cè)試得到的頻率響應(yīng)函數(shù)(Frequency Response Function(FRF))經(jīng)過FFT逆變換得到對(duì)應(yīng)的時(shí)間域脈沖響應(yīng)函數(shù)(Impulse Response Function(IRF))，并通過漢克爾奇異值分解(Hankel Singular Value Decomposition (HSVD))，進(jìn)一步將脈沖響應(yīng)函數(shù)轉(zhuǎn)化為按能量從大到小順序排列的一系列分量信號(hào)組合;其次，以恢復(fù)所有關(guān)心的模態(tài)峰值為基準(zhǔn)，將分量信號(hào)從前到后累加，并在所關(guān)心的模態(tài)峰值完全恢復(fù)后，將剩余分量信號(hào)當(dāng)作噪聲舍棄掉;該步驟會(huì)去除掉信號(hào)中包含的大部分噪聲，但仍會(huì)有一些殘余噪聲不可避免地被恢復(fù);再次，對(duì)步驟二中提取得到的分量信號(hào)，采用基于模態(tài)峰值頻率通帶的迭代選取進(jìn)行二次濾波，以分離出屬于模態(tài)峰值的分量信號(hào)，進(jìn)而將它們累加為降噪后的IRF信號(hào)，并轉(zhuǎn)換至頻域以獲取降噪后FRF信號(hào)；最后，對(duì)降噪后的頻響函數(shù)進(jìn)行模態(tài)辨識(shí)以提取模態(tài)參數(shù)。本方法應(yīng)用于600度高溫環(huán)境下一個(gè)直板葉片的激光掃描模態(tài)測(cè)試數(shù)據(jù)的處理，結(jié)果表明頻響函數(shù)中的噪聲被有效濾去，模態(tài)參數(shù)可準(zhǔn)確地提取，顯示了方法的有效性和優(yōu)越性。

高溫；非接觸測(cè)量；漢克爾奇異值分解；模態(tài)峰值；降噪處理

0 引言

在航空、航天工程中，如飛行器等結(jié)構(gòu)受到高溫的作用，其剛度特性會(huì)發(fā)生變化，從而會(huì)改變結(jié)構(gòu)的動(dòng)力學(xué)特性。因此，開展高溫結(jié)構(gòu)的模態(tài)測(cè)試與參數(shù)識(shí)別，以獲取高溫條件下的結(jié)構(gòu)動(dòng)力學(xué)特性，具有重要的意義。在高溫環(huán)境下，接觸式傳感器的使用受到限制，采用非接觸式的激光傳感器進(jìn)行高溫環(huán)境的模態(tài)測(cè)試，具有很高的優(yōu)越性。但是，高溫環(huán)境會(huì)改變結(jié)構(gòu)的表面狀態(tài)，會(huì)不可避免地使結(jié)構(gòu)表面反射而來的激光帶有一定的噪聲，進(jìn)而使測(cè)試得到的頻響函數(shù)包含噪聲，并干擾模態(tài)參數(shù)的準(zhǔn)確提取。因此，需要研究一種相應(yīng)的去噪方法，實(shí)現(xiàn)強(qiáng)噪聲影響的測(cè)試頻響函數(shù)的模態(tài)參數(shù)識(shí)別。

在過去的幾十年內(nèi)，信號(hào)降噪手段得到了很大的發(fā)展，常用的信號(hào)處理方法有奇異值分解(Singular Value Decomposition(SVD))降噪[1]，小波降噪[2]，經(jīng)驗(yàn)?zāi)B(tài)分解降噪(Empirical Modal Decomposition(EMD))[3]等。在這些方法中，由于奇異值分解具有簡單，非參數(shù)的特性，其應(yīng)用最為廣泛。近年來，基于脈沖響應(yīng)信號(hào)(IRF)的漢克爾奇異值分解(HSVD)方法在動(dòng)力學(xué)測(cè)試數(shù)據(jù)降噪中得到了廣泛應(yīng)用[4]。其工作原理是將含有噪聲的IRF信號(hào)分解為一系列的奇異值和對(duì)應(yīng)的分量信號(hào)，通過選取真實(shí)信號(hào)的奇異值和分量信號(hào)進(jìn)行重構(gòu)，舍棄由噪聲信號(hào)產(chǎn)生的奇異值和分量信號(hào)以實(shí)現(xiàn)降噪處理。然而，如何選取有用的分量信號(hào)或奇異值至關(guān)重要。Zhao[5]提出了一種基于奇異值差分譜(Difference Spectrum of Singular Values(DSSV))的降噪方法，該方法提取了DSSV最大值及之前所對(duì)應(yīng)的奇異值以完成降噪處理。Bao[6]提出了一種基于模型階數(shù)指標(biāo)(Model Order Indicator(MOI))的奇異值選取方法，該方法選取了MOI最大值所對(duì)應(yīng)的奇異值及之前的數(shù)據(jù)進(jìn)行重構(gòu)以完成降噪。以上兩種方法主要利用了真實(shí)信號(hào)和噪聲所對(duì)應(yīng)的奇異值之間會(huì)有一個(gè)巨大的差異這一特點(diǎn)來選取有效的奇異值。但是，當(dāng)信號(hào)中噪聲能量較大而待提取的信號(hào)特征能量較小時(shí)，上述方法就不再適用。近年來，學(xué)者們進(jìn)一步提出了一些基于分量信號(hào)特性的選取方法來提升HSVD的降噪效果。文獻(xiàn)[7]提出了一種基于周期調(diào)制強(qiáng)度(Periodic Modulation Intensity(PMI))的分量信號(hào)選取方法，有效地提升了旋轉(zhuǎn)機(jī)械故障信號(hào)的故障特征。文獻(xiàn)[8]提出了一種基于信號(hào)頻率分量的奇異值選取方法，通過分析分量信號(hào)、奇異值和頻率分量之間的關(guān)系，有效地從多頻率分量信號(hào)中提取出了單一頻率分量。上述方法均在旋轉(zhuǎn)機(jī)械故障信號(hào)處理中取得了較好的效果?；谀B(tài)峰值的漢克爾奇異值分解的降噪方法可用于單個(gè)FRF信號(hào)的降噪處理，其采用累加迭代和二次選取迭代來確定IRF信號(hào)中與模態(tài)信息相關(guān)的分量信號(hào)，進(jìn)而將它們與噪聲產(chǎn)生的分量信號(hào)分離出來，以合成降噪后IRF信號(hào)，再將其轉(zhuǎn)至頻域的FRF信號(hào)。該方法無需提取與頻率分量相關(guān)的奇異值，并以其作為中間參數(shù)提取相應(yīng)的分量信號(hào)來完成頻率分量的提取和降噪，而是采用待提取頻率分量的頻率和振幅，直接選取對(duì)應(yīng)的分量信號(hào)來完成頻率成分的提取和降噪處理，因此具有更高的可操縱性。本文將該方法應(yīng)用于600度高溫環(huán)境下直板葉片的激光掃描模態(tài)測(cè)試數(shù)據(jù)的降噪處理和模態(tài)參數(shù)提取，以驗(yàn)證該方法對(duì)高溫結(jié)構(gòu)激光掃描模態(tài)測(cè)試數(shù)據(jù)的可適用性。結(jié)果表明測(cè)試數(shù)據(jù)中的噪聲被有效地消除，同時(shí)降噪后的頻響函數(shù)可用于模態(tài)參數(shù)的準(zhǔn)確提取。

1 基于模態(tài)峰值漢克爾奇異值分解的降噪方法與參數(shù)識(shí)別

基于模態(tài)峰值漢克爾奇異值分解(Modal Peak-based Hankel Singular Value Decomposition (MPHSVD))的降噪方法是采用了模態(tài)峰值分量信號(hào)選取方法的漢克爾奇異值分解濾波器，可用于單通道脈沖響應(yīng)函數(shù)或頻率響應(yīng)函數(shù)信號(hào)的降噪處理，降噪處理后的頻率響應(yīng)函數(shù)可以直接通過模態(tài)參數(shù)識(shí)別方法，獲得準(zhǔn)確的模態(tài)參數(shù)。

1.1 脈沖響應(yīng)函數(shù)的漢克爾奇異值分解

它所對(duì)應(yīng)的漢克爾(Hankel)矩陣為

式（3）中，為Hankel矩陣，為Hankel矩陣的行數(shù)，為Hankel矩陣的列數(shù)，同時(shí)+-1=。為了使Hankel矩陣近似為一個(gè)方陣，的大小設(shè)置為(/2)，‘’的意思為向下取整。設(shè)置Hankel矩陣接近方陣的目的是為了達(dá)到HSVD最大的降噪效果。在形成Hankel矩陣后，矩陣將被SVD分解為分量矩陣

1.2 分量信號(hào)累加迭代

對(duì)包含了一系列模態(tài)信息及部分噪聲的IRF信號(hào)來說，漢克爾奇異值分解可將其分解為一系列按奇異值大小從大到小排列的分量信號(hào)，IRF信號(hào)中的某個(gè)模態(tài)信息可分解為一對(duì)或多對(duì)相鄰的單頻率成分的分量信號(hào)。通常，某階模態(tài)信息分解得到的分量信號(hào)所對(duì)應(yīng)的奇異值的大小要大于大部分噪聲信號(hào)所對(duì)應(yīng)的奇異值，因此，若按照奇異值的大小順序迭代累加分量信號(hào)，即可在一個(gè)未知的迭代次數(shù)恢復(fù)IRF信號(hào)中所有關(guān)心的模態(tài)信息，并且通過IRF的傅里葉變換在FRF中恢復(fù)所關(guān)心的各階模態(tài)。迭代累加過程中舍棄掉的分量信號(hào)，為大部分低能量的噪聲。

式中，E代表了所對(duì)應(yīng)的能量譜。之后，計(jì)算模態(tài)峰值幅度和能量譜E在峰值頻率處的幅度差異

累加迭代將會(huì)停止，否則，迭代將繼續(xù)進(jìn)行。在累加結(jié)束后，所關(guān)心的模態(tài)特性被恢復(fù)，并記錄選取的分量信號(hào)個(gè)數(shù)為。

該步驟會(huì)去除大部分低能量的高斯噪聲和少部分諧波噪聲，然而，一些高能量噪聲所對(duì)應(yīng)的分量信號(hào)在分量信號(hào)序列中的位置要高于某些模態(tài)峰值頻率成分所對(duì)應(yīng)的分量信號(hào)，因此，在累加迭代中，這些高能量噪聲的恢復(fù)要早于某些模態(tài)峰值頻率成分，進(jìn)而會(huì)產(chǎn)生一定的殘余噪聲，因此，需要進(jìn)行二次濾波去除殘余的高能量噪聲。

1.3 基于模態(tài)峰值頻率通帶的分量信號(hào)迭代選取

對(duì)累加迭代所恢復(fù)的殘余高能量噪聲，其存在形式大多形如諧波分量，所對(duì)應(yīng)的分量信號(hào)也近似為單頻率成分信號(hào)。進(jìn)而，由于高能量噪聲和模態(tài)峰值頻率成分所對(duì)應(yīng)的分量信號(hào)為單頻率成分信號(hào)，故均可通過其在頻域能量譜中峰值振幅所對(duì)應(yīng)的峰值頻率指代該分量信號(hào)。同時(shí)，模態(tài)峰值頻率成分所對(duì)應(yīng)的分量信號(hào)的峰值頻率會(huì)在模態(tài)峰值頻率附近，而高能量噪聲對(duì)應(yīng)的分量信號(hào)的峰值頻率則會(huì)在一定程度上遠(yuǎn)離模態(tài)峰值頻率，因此，若在模態(tài)峰值頻率附近設(shè)置一系列的頻率通帶，即可通過選取峰值頻率在通帶內(nèi)部的分量信號(hào)用于重構(gòu)，并舍棄不在通帶內(nèi)部的分量信號(hào)以去除殘余高能量噪聲。

由于累加迭代共選取了個(gè)分量信號(hào)，迭代次數(shù)為的迭代選取被用于去除殘余噪聲，迭代中的第次迭代如下。帶寬的大小靠經(jīng)驗(yàn)選取，同時(shí)通過后續(xù)的試驗(yàn)數(shù)據(jù)分析表明帶寬不能選取太大，也不能選取太小，5Hz或10Hz的帶寬便足夠?yàn)V波器保留模態(tài)特性并去除噪聲。

1.4 基于有理分式多項(xiàng)式擬合的模態(tài)辨識(shí)

2 高溫結(jié)構(gòu)激光掃描測(cè)量數(shù)據(jù)降噪和模態(tài)參數(shù)辨識(shí)

2.1 試驗(yàn)裝置和測(cè)試數(shù)據(jù)分析

本文所采用的試驗(yàn)裝置為一端固支的直板葉片，材料為45鋼。直板葉片尺寸為200mm×60mm×3mm。左下角帶有一個(gè)3mm的通孔，用于與激振器相連，孔中心距自由端底部6mm，左邊12mm。固支端是一個(gè)60×60×7mm的鋼板，共有四個(gè)10mm的通孔分布在固支端的四周，用于固定，其孔中心距相鄰兩邊的距離均為12mm，固支端與自由端之間采用2mm的圓弧接觸連接在一起。試驗(yàn)時(shí)，固支端被四個(gè)螺栓固定在一鋼板上，鋼板與試驗(yàn)臺(tái)采用螺栓鎖緊以模擬剛性支撐。直板葉片后方固定有一連有電磁加熱器的加熱板，以用于加熱直板來模擬高溫環(huán)境。激振器采用懸掛支撐，激振器的頂針和直板葉片之間通過一碳桿透過加熱板相連，以保證在高溫環(huán)境下激振器可正常工作，并將高溫環(huán)境與激振器隔離開來。采用非接觸掃描多普勒激光測(cè)振儀(Scanning Laser Doppler Vibrometer(SLDV))進(jìn)行測(cè)量，獲取激勵(lì)的響應(yīng)信號(hào)。試驗(yàn)的測(cè)點(diǎn)為通過激光對(duì)焦得到的9×5個(gè)測(cè)點(diǎn)，如圖1(a)所示，這些測(cè)點(diǎn)均勻在直板葉片的自由端四周。試驗(yàn)激勵(lì)為電腦控制的偽隨機(jī)激勵(lì)，采樣頻率設(shè)置為6400Hz，并通過激勵(lì)點(diǎn)和響應(yīng)點(diǎn)的測(cè)試得到測(cè)點(diǎn)的頻率響應(yīng)函數(shù)。

圖1 測(cè)點(diǎn)及結(jié)構(gòu)表面情況

測(cè)量時(shí)用電磁加熱器加熱，熱電偶測(cè)溫。當(dāng)直板葉片加熱到600度時(shí)，鋼板出現(xiàn)了圖1(b)所示的明顯的紅熱現(xiàn)象，同時(shí)可發(fā)現(xiàn)由于加熱不均勻，鋼板左右兩側(cè)的溫度明顯高于中間及上下兩端。進(jìn)而通過檢查600度下的頻響函數(shù)數(shù)據(jù)可發(fā)現(xiàn)，溫度較高的左右兩端及其邊緣的部分測(cè)點(diǎn)出現(xiàn)了嚴(yán)重的噪聲，其余部分測(cè)點(diǎn)則噪聲較小。推測(cè)噪聲產(chǎn)生的原因?yàn)檩^高的溫度改變了結(jié)構(gòu)的表面反射情況，使測(cè)試得到的頻響函數(shù)帶有一定的散斑噪聲，因?yàn)樯咴肼暯?jīng)常被認(rèn)為是激光測(cè)振中出現(xiàn)的主要噪聲[10]，而左右兩側(cè)極高的溫度同測(cè)試用到的紅色激光發(fā)生了干涉，進(jìn)一步增強(qiáng)了測(cè)試中存在的散斑噪聲，使得該區(qū)域部分測(cè)點(diǎn)的噪聲嚴(yán)重。本文選取了圖1(a)所示包含極大噪聲的43號(hào)測(cè)點(diǎn)和包含較小噪聲的38號(hào)測(cè)點(diǎn)的600度頻響函數(shù)進(jìn)行分析，頻響函數(shù)的對(duì)比見圖2。

圖2 測(cè)試頻響函數(shù)的對(duì)比

圖2 (a)和(b)給出了頻響函數(shù)的對(duì)比圖像。圖2 (a)中，淡藍(lán)色曲線是測(cè)試得到的43號(hào)測(cè)點(diǎn)頻響函數(shù)，黑色曲線為38號(hào)測(cè)點(diǎn)的頻響函數(shù)，可以明顯地觀察到43號(hào)測(cè)點(diǎn)的測(cè)試數(shù)據(jù)包含極大的噪聲，并且多數(shù)噪聲的振幅要高于模態(tài)頻率的振幅。圖2 (a)中紅色和藍(lán)色圓環(huán)所包含的峰值是所關(guān)心的模態(tài)峰值，紅色圓環(huán)內(nèi)的模態(tài)峰值是在38號(hào)測(cè)點(diǎn)和43號(hào)測(cè)點(diǎn)均能觀察到的模態(tài)峰值，而藍(lán)色圓環(huán)是只能在38號(hào)測(cè)點(diǎn)中觀察到的模態(tài)峰值，43號(hào)測(cè)點(diǎn)相應(yīng)的模態(tài)峰值已經(jīng)完全湮沒在噪聲當(dāng)中，故無法恢復(fù)。從圖2 (b)中可以發(fā)現(xiàn)，38號(hào)測(cè)點(diǎn)頻響函數(shù)相位譜包含一定程度得噪聲，同時(shí)，相比于38號(hào)測(cè)點(diǎn)的數(shù)據(jù)，可發(fā)現(xiàn)43號(hào)測(cè)點(diǎn)的大部分相位譜已湮沒在噪聲當(dāng)中，無法分清各階模態(tài)的真實(shí)相位。

2.2 測(cè)量數(shù)據(jù)降噪處理

在累加迭代結(jié)束后，將進(jìn)行基于模態(tài)峰值頻率通帶的迭代選取，選取峰值頻率在設(shè)置的頻率通帶內(nèi)的分量信號(hào)做進(jìn)一步的降噪處理。考慮到38號(hào)測(cè)點(diǎn)共有8個(gè)關(guān)心的模態(tài)峰值頻率而43號(hào)測(cè)點(diǎn)共有4個(gè)關(guān)心的模態(tài)峰值頻率。8個(gè)在模態(tài)峰值頻率附近帶寬為10Hz的頻率帶被設(shè)置為38號(hào)測(cè)點(diǎn)的頻率通帶。相應(yīng)的，4個(gè)帶寬為10Hz的頻率帶被設(shè)置為43號(hào)測(cè)點(diǎn)的頻率通帶。下面，進(jìn)行迭代選取以分離出屬于模態(tài)峰值的分量信號(hào)，以模態(tài)5為例，屬于38號(hào)測(cè)點(diǎn)頻響函數(shù)模態(tài)5的分量信號(hào)選取方法如圖4所示。

圖4 模態(tài)分量信號(hào)選取方法

如圖4所示，模態(tài)5的模態(tài)峰值頻率為1044Hz，根據(jù)該頻率設(shè)置的頻率通帶為[1039Hz，1049Hz]，如圖4中黑色虛線所示。該模態(tài)主要分解為了圖4中紫色和綠色曲線所代表的分量信號(hào)57和59，同時(shí)，該模態(tài)還有小于百分之1的能量混疊在了一峰值頻率為531.5的高能量噪聲中。分量信號(hào)57和59的峰值頻率為1044Hz，在黑色虛線所圍成的通帶區(qū)間內(nèi)部，進(jìn)而被選為屬于模態(tài)5的分量信號(hào)以用于恢復(fù)降噪后IRF信號(hào)，而由于分量信號(hào)55所代表的高能量噪聲峰值頻率為531.5Hz，不在頻率通帶內(nèi)，分量信號(hào)55即被舍棄。

對(duì)38號(hào)測(cè)點(diǎn)在累加迭代中選取得到的59個(gè)分量信號(hào)，按照?qǐng)D4所示的方法共選取出了36個(gè)屬于模態(tài)特性的分量信號(hào)，對(duì)43號(hào)測(cè)點(diǎn)在累加迭代中選取得到的181個(gè)分量信號(hào)，則選取出了10個(gè)屬于模態(tài)特性的分量信號(hào)。這些選取得到的屬于模態(tài)特性的分量信號(hào)被進(jìn)一步加合在一起以構(gòu)成降噪后的IRF信號(hào)，進(jìn)而被轉(zhuǎn)換至頻域以獲取相應(yīng)的降噪后FRF信號(hào)。降噪后FRF信號(hào)與測(cè)量FRF信號(hào)的對(duì)比分析見圖5。圖5 (a)和(b)給出了38號(hào)測(cè)點(diǎn)降噪后頻響函數(shù)的對(duì)比圖像，可以發(fā)現(xiàn)，圖像上所能觀察到的峰值只有所關(guān)心的模態(tài)峰值，F(xiàn)RF中存在的噪聲和不關(guān)心的峰值被很好的濾除掉了。圖5 (c)給出了43號(hào)測(cè)點(diǎn)的降噪前后FRF的對(duì)比，可以發(fā)現(xiàn)，F(xiàn)RF信號(hào)中大部分噪聲都被去除掉了，但仍剩余一些能量遠(yuǎn)低于模態(tài)特性的噪聲，這是由于噪聲較大，部分噪聲混疊在了模態(tài)特性所對(duì)應(yīng)的分量信號(hào)中，進(jìn)而隨著模態(tài)特性的恢復(fù)而恢復(fù)所造成的，但這些未能去除的噪聲相比于模態(tài)特性能量較小，已不影響進(jìn)一步的模態(tài)分析。圖5 (d)給出了降噪前后FRF相位譜的對(duì)比。相比于降噪前的FRF相位譜，相位譜中大部分噪聲被去除掉了，同時(shí)屬于四個(gè)模態(tài)峰值的相位突變被很好的恢復(fù)，同理，由于存在小能量峰值，F(xiàn)RF信號(hào)的能量譜并不光滑，并存在部分從180至-180的相位突變。從降噪前后的FRF對(duì)比來看，降噪保留了所有關(guān)心的模態(tài)特性，并去除了大部分噪聲，這證明了方法的有效性和優(yōu)越性。

圖5 降噪后信號(hào)對(duì)比分析

2.3 降噪后數(shù)據(jù)的模態(tài)辨識(shí)

在完成了降噪處理后，要對(duì)降噪后的數(shù)據(jù)進(jìn)行模態(tài)辨識(shí)，本文選取的模態(tài)辨識(shí)方法為RFP法，該方法對(duì)降噪后得到的FRF信號(hào)進(jìn)行擬合以獲取相應(yīng)的模態(tài)參數(shù)，降噪后FRF信號(hào)的擬合是分段處理的。擬合所選取的頻段和擬合結(jié)果見圖6所示。圖6 (a)和(b)給出了38號(hào)測(cè)點(diǎn)降噪后FRF擬合結(jié)果，可以觀察到曲線與降噪后曲線吻合度很高。圖6 (c)和(d)給出了43號(hào)測(cè)點(diǎn)降噪后FRF信號(hào)的擬合結(jié)果，同樣可以觀察到降噪后曲線和擬合曲線有較高的吻合度，同時(shí)可以發(fā)現(xiàn)，擬合頻段內(nèi)能量較小的噪聲并未對(duì)擬合結(jié)果造成影響。

在完成了曲線擬合后，將從擬合得到的曲線中提取模態(tài)參數(shù)，提取的參數(shù)為模態(tài)頻率和模態(tài)阻尼，并計(jì)算了兩測(cè)點(diǎn)模態(tài)參數(shù)之間的差異，提取和計(jì)算結(jié)果見表1所示。表中，模態(tài)頻率差異計(jì)算的是提取前后模態(tài)頻率的絕對(duì)差，模態(tài)阻尼比的差計(jì)算如下

式中，為模態(tài)阻尼比的差異，為由38號(hào)測(cè)點(diǎn)提取得到的降噪后FRF信號(hào)模態(tài)阻尼比，則是由43號(hào)測(cè)點(diǎn)提取得到的降噪后FRF信號(hào)的模態(tài)阻尼比。由表1知，由38號(hào)測(cè)點(diǎn)和43號(hào)測(cè)點(diǎn)降噪后頻響函數(shù)提取出來的模態(tài)頻率最大差異為1.6Hz，模態(tài)阻尼比的最大差異約為16%，表明模態(tài)頻率和模態(tài)阻尼比的提取是準(zhǔn)確的，驗(yàn)證了降噪手段的有效性和優(yōu)越性。

表1 提取得到的模態(tài)參數(shù)

3 結(jié)論

針對(duì)高溫結(jié)構(gòu)激光掃描模態(tài)測(cè)試中出現(xiàn)的噪聲問題，提出了一種基于模態(tài)峰值漢克爾奇異值分解的降噪和模態(tài)識(shí)別方法。該方法首先通過分量信號(hào)累加迭代恢復(fù)所關(guān)心的模態(tài)信息，同時(shí)會(huì)一并恢復(fù)一定的殘余噪聲，進(jìn)而通過基于模態(tài)頻率通帶的迭代選取分離出屬于模態(tài)特性的分量信號(hào)，并將其累加實(shí)現(xiàn)2次降噪處理。本文采用600度高溫環(huán)境下一直板葉片結(jié)構(gòu)的激光掃描模態(tài)測(cè)試數(shù)據(jù)進(jìn)行降噪和模態(tài)參數(shù)識(shí)別，結(jié)果表明，本方法能有效地去除頻響函數(shù)中的強(qiáng)噪聲，并保留結(jié)構(gòu)的模態(tài)特性，準(zhǔn)確識(shí)別結(jié)構(gòu)的模態(tài)參數(shù)。該方法對(duì)航空、航天工程中的高溫模態(tài)試驗(yàn)數(shù)據(jù)的分析有指導(dǎo)和借鑒作用。

[1] Golafshan R, Yuce Sanliturk K. SVD and Hankel matrix based de-noising approach for ball bearing fault detection and its assessment using artificial faults[J]. Mechanical Systems and Signal Processing, 2016(70-71): 36-50.

[2] Fang H-T, Huang D-S. Noise reduction in lidar signal based on discrete wavelet transform [J]. Optics Communications, 2004, 233(1-3): 67-76.

[3] Rong H, Gao Y , Guan L, et al. GAM-Based Mooring Alignment for SINS Based on An Improved CEEMD Denoising Method[J]. Sensors, 2019, 19(16): 3564.

[4] Sanliturk K Y, Cakar O. Noise elimination from measured frequency response functions[J]. Mechanical systems and signal processing, 2005, 19(3): 615-631.

[5] Zhao X, Ye B. Selection of effective singular values using difference spectrum and its application to fault diagnosis of headstock[J]. mechanical systems & signal processing, 2011, 25(5): 1617-1631.

[6] Bao X X, Li C L, Xiong C B. Noise elimination algorithm for modal analysis [J]. Applied Physics Letters, 2015, 107(4).

[7] Zhao M, Jia X . A novel strategy for signal denoising using reweighted SVD and its applications to weak fault feature enhancement of rotating machinery[J]. Mechanical Systems & Signal Processing, 2017, 94: 129-147.

[8] Zhao X, Ye B. Separation of Single Frequency Component Using Singular Value Decomposition[J]. Circuits, Systems, and Signal Processing, 2018, 38(1): 191-217.

[9] Formenti, DaveRichardson, et.al. Parameter estimation from frequency response measurements using rational fraction polynomials (twenty years of progress)[C]. Proceedings of SPIE - The International Society for Optical Engineering, 1982.

[10] MartinP, RothbergS. Introducing speckle noise maps for Laser Vibrometry[J]. Optics and Lasers in Engineering, 2009, 47 (3-4): 431-442.

[11] 王智勇, 王則力, 宮文然, 等. 熱結(jié)構(gòu)高溫應(yīng)變光學(xué)測(cè)量技術(shù)發(fā)展探討[J]. 強(qiáng)度與環(huán)境,2019,46(6): 1-8. [WANG Zhi-yong, WANG Ze-li, GONG Wen-ran, et al. Discussion on the Development of High Temperature Strain Optical Measurement Technology for Thermal Structure[J]. Structure & Environment Engineering, 2019,46(6): 1-8.]

[12] 白志富, 張忠, 朱儀凡, 等.高溫動(dòng)力學(xué)模型誤差定位與修正方法研究[J]. 強(qiáng)度與環(huán)境,2019,46(2):21-26.[BAI Zhi-fu, ZHANG Zhong, ZHU Yi-fan, et al. Discussion on the Development of High Temperature Strain Optical Measurement Technology for Thermal Structure[J]. Structure & Environment Engineering, 2019,46(2): 21-26.]

Noise Reduction and Parameter Identification of Modal Test on a High Temperature Structure Measured by Scanning Laser Doppler Vibrometer

ZHU Tian-xu ZANG Chao-ping

(College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

The noise reduction and modal identification method of the modal peak based Hanke lSingular Value Decomposition(HSVD) is proposed in this paper to deal with the noise contaminated in modal testing of structures in the high temperature environment using a non-contact Scanning Laser Doppler Vibrometer (SLDV). Firstly, the measured frequency response function (FRF) signal is transformed to the impulse response function (IRF) by the inverse Fast Fourier Transform (iFFT), and the IRF signal is subsequently decomposed to a set of component signals using the HSVD method. These component signals are ordered according to the singular values from large to small. Secondly, the component signals are sequentially accumulated together with the ordered amplitudes of singular values until the all concerned modal modes are recovered and the other rest abandoned component signals are treated to be the noise. This step can eliminate most parts of the noise contaminated in the signal, but some residue noise is still inevitably recovered. Thirdly, the further iterative noise reduction based on the narrow frequency pass bands of the modal mode is conducted to separate the component signals of the modal features from the noise and the final deposing IRF signal is transformed to the FRF signal. Fourthly, the modal analysis is conducted to the noise reduced FRF signal to extract the modal parameters. The method is applied to measured data from modal testing of a straight plate at 600 centigrade degrees. Results shows that the strong noise contaminated in FRFs is effectively removed and modal parameters are precisely extracted. It demonstrates the validity and superiority of the method.

high temperature; noncontact measurement; the modal peak; HankelSingularValueDecomposition; noise reduction

V414.1

A

1006-3919(2020)02-0046-10

10.19447/j.cnki.11-1773/v.2020.02.008

2019-11-20；

2020-01-17

國家自然科學(xué)基金委員會(huì)與中國工程物理研究院聯(lián)合基金(U1730129)

朱天煦（1998—），男，本科生，研究方向：結(jié)構(gòu)動(dòng)力學(xué)；(210016)南京航空航天大學(xué)能源與動(dòng)力學(xué)院.