一種基于數據間相關性的激光噴丸聲學監測技術

邱辰霖, 程 禮 , 何衛鋒

(1. 空軍工程大學 航空航天工程學院，西安 710038；2.北京航空航天大學 能源與動力工程學院，北京 100191)

一種基于數據間相關性的激光噴丸聲學監測技術

邱辰霖1,2, 程 禮1, 何衛鋒1

(1. 空軍工程大學 航空航天工程學院，西安 710038；2.北京航空航天大學 能源與動力工程學院，北京 100191)

激光噴丸作為一種新興的金屬材料表面強化技術，與普通噴丸相比，具有更顯著的改性作用，使其在航空等領域有著重要的應用價值。隨著該技術的工業化程度不斷加深，需要加強對其工作狀態的監測診斷，保證其良好的運行狀況。由于聲學信號不僅攜帶有豐富的工作特征信息，而且獲取成本較低，并能夠實現無損檢測，所以將聲學信號應用于激光噴丸的狀態監測。通過分析由等離子體沖擊波衰減所產生的聲波信號，提取信號特征，進行過程監測。針對于沖擊信號的非線性特征，從信號內相鄰數據間相關性的角度，提出了一種新的沖擊聲信號特征挖掘方法。一方面對于模擬沖擊信號進行了分析，另一方面，對于激光噴丸實際聲信號進行了處理，表明該方法能夠識別沖擊信號的特征變化，可以用于監測激光噴丸的工作過程，操作簡單且速度較快，具有較好的應用前景。

激光噴丸；聲學信號；狀態監測；相關性

激光噴丸技術，是一種新興的金屬材料表面強化技術，近些年來不斷發展，在航空等領域都具有重要的應用價值[1]。它借助強激光誘導產生的超強沖擊波，產生更多的位錯密度和剩余壓應力，達到比普通噴丸更強的處理效果，顯著延長金屬材料的疲勞壽命、改善材料的抗腐蝕斷裂能力[2]。隨著該技術工業化程度不斷加深，對于其的狀態監測和故障診斷，愈加顯得重要。實驗階段采用儀器對于材料性能參數逐一檢測的方法無法應付繁重的工業生產，需要開發一種簡便、無損的在線監測手段，以滿足工業化生產的要求[3]。

聲學信號獲取成本低，且包含有豐富的特征信息，聲學診斷技術已經在許多方面都有成熟的應用，如旋轉機械的故障診斷等[4-6]。在激光噴丸過程中產生的等離子體沖擊波，及其衰減而得的聲波信號，與強化過程有密切的聯系，國內外學者通過分析沖擊聲波信號的傳播特征[7-10]，將等離子體沖擊波及其衰減產生的聲波信號用于監測激光噴丸的工作狀態[11-13]。由于沖擊波信號的非線性、非平穩特征，給聲波信號分析提出了更高的要求，為了開發一種對于等離子體聲信號處理效果更佳的在線狀態監測方法，本文提出一種基于信號內相鄰數據間相互關聯性的特征分析方法，將聲信號數據列轉換成特征更明顯的相關系數序列，提取特征信息進行狀態識別。通過模擬分析和實驗驗證，表明該方法能夠通過沖擊聲信號對于激光噴丸強化進行監測，效果較好且操作簡便，具有實用潛力。

1 激光噴丸

激光噴丸強化，也稱為激光沖擊強化處理，是一種利用高功率脈沖激光對于材料進行表面改性的新技術[14]，其工作原理圖如圖1所示。

圖1 激光噴丸強化原理圖Fig.1 Schematic diagram of laser peening

圖1是約束模式的激光噴丸強化，也是工業生產中常用的工作模式。高功率的激光脈沖穿過約束層照射到保護層上，保護層瞬間蒸發為等離子體，并且在約束層的束縛下繼續吸收激光能量，產生高溫、高壓的等離子體沖擊波，此沖擊波在約束層的作用下大部分向表層材料傳播，改變材料表層的微觀結構，增加其位錯密度，提高其殘余壓應力，最終達到對于材料表層改性的效果。由于沖擊波的高溫高壓特征，使得其強化效果比一般的噴丸強化有很大的提升。

離子體沖擊波是激光噴丸強化的關鍵，對于最終的材料特性改善起著重要作用。在沖擊波向內傳播作用于材料的同時，也有部分沖擊波穿過約束層，向空氣傳播，它與材料內部的沖擊波具有一定的關系[15]。此沖擊波會迅速衰減，產生聲波，分析此聲波信號，就可以獲取與材料內沖擊波信號相關的特征信息，從而達到對于強化效果的狀態監測。

2 聲學信號分析

2.1 模擬沖擊信號

聲波信號是一組時間序列信號，由于前后時間點具有天然的相關性，使得基于時間順序的聲波信號的數據點間也是相互關聯的。本文通過Pearson相關系數度量相鄰數據間的相關性，依次計算相鄰數據子列間的相關系數，組成相關系數序列，從而將復雜的原始序列轉變為一個特征更明顯的簡單序列，從中提取有用信息。

激光噴丸所產生的聲波信號如圖2所示。

圖2 激光噴丸聲波信號Fig.2 Acoustic signal of laserpeening

此信號具有典型的沖擊波特征，短的上升時間、超高的壓力峰值和快速的衰減過程。通過式(1)對于沖擊波信號進行簡單模擬[16]。

(1)

式中:A0為模擬沖擊波信號的峰值;θ為衰減常數;t為時間。由于在實際應用中不可避免地會混有噪聲信號，所以在模擬沖擊信號中加入白噪聲。

取A0=300，θ=0.01，t∈(0,1)，加入白噪聲信號，其沖擊波形如圖3所示。

圖3 模擬沖擊波信號示意圖(A0=300，θ=0.01)Fig.3 Schematic diagram of simulating shock wave signal(A0=300，θ=0.01)

采用本文的方法，依次計算相鄰數據子列(長度取2 000點)間的相關系數，其變化曲線如圖4所示。

圖4 沖擊波信號相鄰數據相關系數變化曲線(A0=300，θ=0.01)Fig.4 Variation of the neighbor subsequences’ correlation coefficient obtained from the simulating shock wave signal(A0=300，θ=0.01)

由圖4可知，相關系數變化曲線有一個明顯的凸起，對應于時域信號中沖擊波峰值出現的區域，可以將其凸起處的極大值作為特征值，此時該值為0.96。

2.2 沖擊信號特征識別

當θ=0.01時，改變沖擊波幅值，取A0=100，其相鄰數據相關系數變化曲線如圖5所示。

圖5 沖擊波信號相鄰數據相關系數變化曲線(A0=100，θ=0.01)Fig.5 Variation of the neighbor subsequences’ correlation coefficient obtained from the simulating shock wave signal(A0=100，θ=0.01)

由圖5可知，其變化曲線凸起處的極大值為0.74，可見隨著沖擊波幅值的減小，特征值也會明顯減小，所以，可以通過該特征值識別沖擊波模擬信號的幅值變化。

當A0=300時，改變衰減常數，取θ=1/150，其相鄰數據相關系數的變化曲線如圖6所示。

圖6 沖擊波信號相鄰數據相關系數變化曲線(A0=100，θ=1/150)Fig.6 Variation of the neighbor subsequences’ correlation coefficient obtained from the simulating shock wave signal(A0=300，θ=1/150)

由圖6可知，其曲線凸起處的極大值為0.5，可見，隨著衰減常數的減小，特征值明顯減小。所以，可以通過該特征值識別衰減常數的變化。

2.3 分析小結

采用本文提出的基于相鄰數據相關性的信號分析方法，通過對于含噪條件下沖擊波模擬信號進行信號處理，可以提取相關系數變化曲線中凸起處的極大值作為特征值；通過該特征值，能夠識別沖擊波模擬信號的幅值和衰減常數等參量的變化。所以，本文分析方法可以用于分析沖擊信號的特征變化。

3 實驗分析

設計激光噴丸實驗，實測等離子體沖擊聲波數據，采用本文的方法進行信號處理。

實驗在西安天瑞達公司的協助下完成，實驗設備是該公司生產用的激光沖擊強化裝置，實驗時激光器參數設置為波長1 064 nm 、光斑直徑2.2～2.4 mm、脈寬20±2 ns，通過聲傳感器(靈敏度40 mv/Pa)采集聲波信號，通過DEWEsoft進行數據儲存。實驗金屬材料采用鈦合金，用黑漆作為保護層，施加流水約束。

實驗過程中，聲傳感器放置的位置有50 cm、60 cm、70 cm、80 cm四個距離值，每一個位置處，激光能量取3 J、4 J、5 J、6 J四個能量值 ，每一種能量下，分別在有約束狀態和無約束狀態下進行強化處理。

3.1 信號分析

在50 cm處，3 J能量下，有約束層狀態時，截取一段沖擊聲波信號，如圖7所示。

圖7 激光噴丸聲波信號Fig.7 Acoustic signal of laser peening

對于所取的信號序列，依次計算相鄰數據子列(長度取2 000點)的相關系數，其變化曲線如圖8所示。

圖8 激光噴丸聲波信號內相關數據相關系數變化曲線Fig.8 Variation of the neighbor subsequence’ correlation coefficient obtained from the shock acoustic signal of laser peening

由圖8可知，當計算相鄰子序列中開始出現沖擊波信號正壓峰值時，相關系數會出現突變，之后相關系數值為一個較穩定值，當沖擊波信號正壓峰值移出計算子序列時，相關系數會再次出現突變。兩個突變值之間相關系數的穩定值，可以用作表征沖擊波特征。

圖中除了沖擊波的曲線特征，還有周期性出現的曲線凸起，周期大約1 s，對應于激光機器工作過程的計時聲音，這個信號在時域中與噪聲信號混合在一起，無法區分，而在相關系數變化曲線中卻可以比較清晰地分辨?？梢?，相關系數法對于聲波信號的分析，可以將時域分辨不清的信號特征，清楚區分出來。

3.2 能量變化

將相關系數曲線中，相關系數值負值突變和正值突變中的平穩值，作為聲波信號的特征值，以表征相應原始沖擊波的特征信息。不同激光能量下，不同距離處，正常工作狀態時此特征值如表1所示。

表1 不同正常工作狀態時沖擊聲波信號特征值對比

由表1對比可知，隨著距離的增大，沖擊聲波信號的特征值會減??；隨著激光能量的提高，該特征值會增大。

以3 J為基準，其他能量狀態下的特征值變化率如表2所示。

表2 不同能量狀態下的特征值變化率

由表2可知，不同能量狀態之間的特征值變化明顯，可以通過特征值對于能量狀態進行監測。

3.3 約束層狀態

不同采集距離處，在不同的激光能量下，對比有約束層和無約束層狀態，其特征值如表3所示。

表3 不同約束狀態下的沖擊聲波特征值對比

由表3對比可知，采集距離和激光能量相同的情況下，有約束狀態下的聲波特征值要大于無約束狀態下的特征值。

無約束狀態對比有約束狀態，其特征值的減少率如表4所示。

表4 無約束狀態對比有約束狀態的特征值減小率

由表4可知，無約束狀態對比有約束狀態，其特征值變化比較明顯，可以通過此聲波特征值對于約束層狀態進行監測。

3.4 參數選取

采用本文的分析方法，對于激光噴丸沖擊聲信號進行信號處理的過程中，計算相關系數時相鄰子序列的長度值，是一個可變參數，會對分析效果有所影響，需要綜合考慮分析精度和速度，進行合理選擇。

所選取相鄰子序列的長度越長，提取的特征值會趨于穩定，但是運行所消耗的時間會更多；如果選取相鄰子序列的長度越短，運行所用時間更少，但是提取的特征值可能還未達到穩定值。對比總結不同子序列長度下的分析結果，將沖擊波從出現到基本衰減所經歷的大致長度取為相鄰子序列長度比較合適，該子序列長度下提取的特征值已趨于穩定，且不會影響到分析速度。

3.5 分析小結

將基于相鄰數據相關性的信號分析方法，應用于激光噴丸沖擊聲波信號的分析處理，將沖擊聲信號序列轉變為一個特征明顯的相關系數序列進行分析。

通過相關系數曲線的變化情況，可以提取相關系數變化曲線中突變值之間的平穩值作為聲波信號的特征值，還可以將信號中混雜的不明顯信號特征清楚區別出來。

通過聲波特征值的對比分析，該特征值會隨著聲波傳播距離的增加而顯著減小，會隨著激光能量的提高而明顯增大，無約束時的特征值要小于有約束狀態。

所以，通過基于相關性分析的信號處理方法，挖掘聲波信號的有用特征信息，能夠監測激光噴丸過程中的能量及約束層狀態。

4 結 論

針對于激光噴丸強化過程的聲學診斷，本文提出一種基于相鄰數據相關性的聲波信號處理方法，通過模擬信號分析和實驗檢驗，可以得出以下結論:

(1)該方法能夠有效地提取出與沖擊波相關的特征信息，而且能夠將時域信號中不明顯的信號特征顯化出來。

(2)該方法能夠用于監測激光能量和約束層狀態變化，且效果比較明顯。

本文方法在使用過程中，操作簡單，速度較快，具有進一步實用的潛力。

[ 1 ] QIAO Hongchao. Experimental investigation of laser peening on Ti17 titanium alloy for rotor blade applications [J].Applied Surface Science, 2015,351: 524-530.

[ 2 ] MONTROSS C S, TAO Wei, LIN Ye, et al. Laser shock processing and its effects on microstructure and properties of metal alloys: a review [J]. International Journal of Fatigue, 2002,24(10): 1021-1036.

[ 3 ] 曹子文，鄒世坤，車志剛，等. 航空發動機葉片激光沖擊處理過程控制研究[J].航空發動機，2011，37(1)： 60-62. CAO Ziwen, ZOU Shikun, CHE Zhigang, et al. Research on process control of laser shock peening for aeroengine blade[J]. Aeroengine, 2011, 37(1): 60-62.

[ 4 ] CHEE K T, PHIL I, DAVID M. A comparative experimental study on the diagnostic and prognostic capabilities of acoustics emission, vibration and spectrometric oil analysis for spur gears [J].Mechanical System and Signal Processing, 2007,21(1): 208-233.

[ 5 ] JENA D P, PANIGRAHI S N. Automatic gear and bearing fault localization using vibration and acoustic signals [J]. Applied Acoustics, 2015,98: 20-33.

[ 6 ] 魯文波，蔣偉康. 利用聲場空間分布特征診斷滾動軸承故障[J].機械工程學報，2012，48(13)：68-72. LU Wenbo, JIANG Weikang. Diagnosing rolling bearing faults using spatial distribution features of sound field [J]. Journal of Mechanical Engineering, 2012, 48(13): 68-72.

[ 7 ] ZHAO Rui, XU Rongqing, SHEN Zhonghua, et al. Dynamics of laser-induced shock wave by optical probe in air [J]. Optik, 2006,117(7): 299-302.

[ 8 ]AGUILERA J A, ARAGON C. Characteristic of laser-induced plasma during its expansion in air by optical emission spectroscopy: observation of strong explosion self-similar behavior [J]. Spectrochimica Acta Part B Atomic Spectroscopy, 2014,97(9):86-93.

[ 9 ] 卞保民，陳建平，楊玲，等. 空氣中激光等離子體沖擊波的傳輸特性研究[J].物理學報，2000，49(3)：445-448. BIAN Baomin, CHEN Jianping, YANG Ling, et al.The transmission characteristic of airborne laser plasma shock wave [J]. Acta Physica Sinica, 2000, 49(3): 445-448.

[10] 卞保民，楊玲，陳笑，等. 激光等離子體及點爆炸空氣沖擊波前運動方程的研究[J].物理學報，2002，51(4)：809-812. BIAN Baomin, YANG Ling, CHEN Xiao, et al. Study of the laser-induced plasmas and the kinematics of shock waves in air by a way intense explosion [J]. Acta Physica Sinica, 2002, 51(4): 809-812.

[11] BURGER M, NIKOLIC Z. Frequency domain and wavelet analysis of the laser-induced plasma shock waves [J].Spectrochimica Acta Part B Atomic Spectroscopy, 2015,110: 70-78.

[12] 鄒彪，陳建平，倪曉武. 強激光等離子體有關物理量的聲學測量[J]. 激光技術，2000，24(6)：366-369. ZOU Biao, CHEN Jianping, NI Xiaowu. Acoustic measurement of parameters of laser induced plasma on the target[J]. Laser Technology, 2000, 24(6): 366-369.

[13] 喬紅超，趙吉賓. 激光沖擊強化在線檢測系統設計及應用[J].激光與光電子學進展，2013,7：51-55. QIAO Hongchao, ZHAO Jibin. Design and implementation of online laser peening detection system[J]. Laser & Optoelectronics Progress, 2013,7:51-55.

[14] 范勇，王聲波，吳鴻興，等. 航空鋁合金材料激光沖擊強化實驗研究[J]. 激光技術，2003，27(4)：273-275. FAN Yong, WANG Shengbo, WU Hongxing, et al. Experimental research of laser shock processing on aerial aluminum alloy [J]. Laser Technology, 2003,27(4): 273-275.

[15] 鄒彪.關于激光等離子體聲波的數學模型[J].數學的實踐與認識，2007，37(15)：65-69. ZOU Biao. Mathematical model of laser plasma voice wave [J]. Mathematics in Practice and Theory, 2007, 37(15): 65-69.

[16] 杜紅棉，祖靜，馬鐵華，等. 自由場傳感器外形結構對沖擊波測試的影響研究[J].振動與沖擊，2011，30(11)：85-89. DU Hongmian, ZU Jing, MA Tiehua, et al. Effect of mount configuration of free-field transducers on shock wave measurement [J].Journal of Vibration and Shock, 2011, 30(11): 85-89.

A condition monitoring method for laser peening based on the correlation between the adjacent aata

QIUChenlin1,2，CHENGLi1，HEWeifeng1

(1. School of Aeronautics and Astronautics Engineering, Air Force Engineering University, Xi’an 710038， China; 2. School of Jet Propulsion, Beijing University of Aeronautics and Astronautics, Beijing 100191， China)

As an innovative surface treatment of metal material, laser peening makes greater performance in improving the mechanical property comparing with the conventional shot peening, and it can bring significant contributions to various fields like aerocraft industry. In accordance with the developing pace of its industrialization, the condition monitoring, aiming at maintaining a good working state, is arousing more and more attentions. The acoustic signal with lower cost contains abundant information related to the operating process, and it can be used as non-destructive condition detection of the laser peening. The shock acoustic signal, which is propagating from the original plasma shock wave, was analyzed to extract the valuable characteristics to monitor the process. Considering the nonlinearity of the shock acoustic signal, a new signal processing method was proposed, based on the correlation between the neighbor data in the signal array. On one hand, it was used in analyzing the simulating shock-wave signal; on the other hand, it was applied in studying the practical signal obtained from the laser shock processing. It was shown that the novel method could detect the shock wave signal efficiently with easy and quick operation. It thus has potential in applications.

laser peening; acoustic signal; condition monitoring; correlation

國家自然科學基金資助項目(51205406)

2015-10-20 修改稿收到日期：2016-02-23

邱辰霖 男，博士生，1988年12月生

程禮 男，教授，1963年4月生

TP391；TN249

A

10.13465/j.cnki.jvs.2017.04.022