龔子維，劉希強(qiáng)，張仲寧，楊 京，程建春，劉翔雄

(1. 南京大學(xué)聲學(xué)研究所，江蘇南京210093；2. 人工微結(jié)構(gòu)科學(xué)與技術(shù)協(xié)同創(chuàng)新中心，江蘇南京210093；3. 華辰精密裝備(昆山)股份有限公司，江蘇昆山215337)

0 引 言

聲發(fā)射(Acoustic Emission, AE)是指材料受到外力或者內(nèi)力作用，發(fā)生塑性變形或斷裂，以瞬態(tài)彈性波的形式釋放應(yīng)變能的現(xiàn)象[1-2]。利用聲發(fā)射檢測(cè)磨削加工過程中砂輪的鈍化狀態(tài)，相比于其他間接檢測(cè)方法具有環(huán)境噪聲易分離、可選特征參數(shù)較多，靈敏度高、響應(yīng)速度快等特點(diǎn)[3]，近年來聲發(fā)射檢測(cè)在刀具鈍化狀態(tài)檢測(cè)中得到廣泛應(yīng)用。文獻(xiàn)[4]研究了單磨粒模型中磨粒磨損過程及聲發(fā)射產(chǎn)生機(jī)理；文獻(xiàn)[5]建立了砂輪磨削的簡(jiǎn)化理論模型，推導(dǎo)出塑性功功率與剪切力、接觸面積以及切入深度等加工參數(shù)的關(guān)系，聲發(fā)射能量大小與砂輪鈍化存在相關(guān)性。

與其他刀具類似，砂輪鈍化程度隨磨削時(shí)間推移而增加，文獻(xiàn)[6-7]分別用加工材料單位去除量所需時(shí)間或能量來表示砂輪鈍化程度，加工過程中砂輪從修整結(jié)束到下一次修整之前可以分為初期鈍化、中期鈍化和嚴(yán)重鈍化三個(gè)鈍化階段。初期鈍化階段砂輪剛經(jīng)過修整，磨粒較為鋒利，發(fā)生的變化以磨粒破碎為主，這種鈍化會(huì)產(chǎn)生新的切削刃，即“自銳現(xiàn)象”，對(duì)磨削能力影響較??；在中期鈍化階段，磨粒發(fā)生塑性變形，表面逐漸磨平，與加工材料接觸面積增大，砂輪鈍化程度平穩(wěn)上升；到了嚴(yán)重鈍化階段，更多磨粒開始出現(xiàn)斷裂、脫落，砂輪加速鈍化，需要及時(shí)修整[6]。圖 1為磨削過程中砂輪鈍化程度變化曲線。

圖1 砂輪磨削鈍化曲線Fig.1 Grinding wheel blunting curve

利用聲發(fā)射信號(hào)來檢測(cè)砂輪的鈍化狀態(tài)，首先需要對(duì)采集到的原始信號(hào)進(jìn)行去噪處理，聲發(fā)射信號(hào)的頻率主要分布在100～500 kHz，其他為低頻的摩擦噪聲和較高頻的電氣噪聲。目前主要采用的信號(hào)處理方法包括小波變換(Wavelet Transform, WT)、經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)等。WT方法利用窗函數(shù)平移將原信號(hào)分解成若干小波分量，這種方法需根據(jù)信號(hào)特點(diǎn)選定小波基，易產(chǎn)生異常諧波。EMD方法屬于遞歸分解，根據(jù)信號(hào)本身的特點(diǎn)將一段信號(hào)分解成一組本征模態(tài)函數(shù)(Intrinsic Mode Function, IMF)，適用于非線性、非平穩(wěn)信號(hào)。針對(duì)EMD方法易產(chǎn)生包絡(luò)線誤差累積、模態(tài)混淆及端點(diǎn)效應(yīng)等問題，Dragomiretskiy等[8]提出變分模態(tài)分解(Variational Mode Decomposition, VMD)，該方法具有完備的數(shù)學(xué)理論基礎(chǔ)，克服了小波分解需定義小波基的問題，通過求解約束變分模型將信號(hào)分解為指定數(shù)量的本征模態(tài)分量，能夠避免端點(diǎn)效應(yīng)、抑制模態(tài)混疊，且分解較為徹底。

對(duì)分解后的分量進(jìn)行篩選并重構(gòu)得到聲發(fā)射信號(hào)后，還需要提取聲發(fā)射信號(hào)中能夠表征砂輪鈍化狀態(tài)的特征參數(shù)。為探究砂輪鈍化程度和聲發(fā)射信號(hào)特征參數(shù)的關(guān)系，文獻(xiàn)[9]通過實(shí)驗(yàn)將聲發(fā)射的振鈴計(jì)數(shù)和均方根值(Root Mean Square, RMS)兩個(gè)特征參數(shù)與砂輪鈍化狀態(tài)聯(lián)系起來；文獻(xiàn)[10]提取聲發(fā)射 RMS信號(hào)，利用聲發(fā)射頻域信號(hào)頻譜矩心特征參數(shù)來表征砂輪鈍化狀態(tài)；文獻(xiàn)[11]證明RMS和能量計(jì)數(shù)兩個(gè)統(tǒng)計(jì)數(shù)據(jù)含有砂輪表面狀態(tài)的重要信息；文獻(xiàn)[12]將功率比(Ratio of Power, RoP)引入砂輪磨損評(píng)估中；文獻(xiàn)[13]提取了聲發(fā)射信號(hào)的幅度、振鈴計(jì)數(shù)、能量計(jì)數(shù)等八個(gè)特征參數(shù)構(gòu)建訓(xùn)練集對(duì)砂輪鈍化狀態(tài)進(jìn)行識(shí)別。以上方法基本實(shí)現(xiàn)了對(duì)砂輪鈍化狀態(tài)的識(shí)別，但均未對(duì)不同特征參數(shù)的作用進(jìn)行量化比較，沒有系統(tǒng)性衡量聲發(fā)射信號(hào)多種特征參數(shù)對(duì)砂輪鈍化狀態(tài)的識(shí)別效果。

特征參數(shù)與砂輪鈍化狀態(tài)的具體對(duì)應(yīng)關(guān)系數(shù)學(xué)模型目前較難建立，國(guó)內(nèi)外相關(guān)研究通過訓(xùn)練神經(jīng)網(wǎng)絡(luò)來識(shí)別砂輪鈍化狀態(tài)，而聲發(fā)射信號(hào)本身存在一定的波動(dòng)，會(huì)影響識(shí)別的準(zhǔn)確性。概率神經(jīng)網(wǎng)絡(luò)(Probabilistic Neural Network, PNN)利用密度函數(shù)估計(jì)方法輸出識(shí)別結(jié)果，對(duì)于非線性數(shù)據(jù)也能保持較高精度，適用于處理非線性、非平穩(wěn)的磨削聲發(fā)射信號(hào)。

綜合以上因素，本文提出一種基于VMD-PNN進(jìn)行聲發(fā)射信號(hào)處理和砂輪鈍化狀態(tài)識(shí)別的方法。利用VMD快速精確地得到聲發(fā)射信號(hào)，選取典型特征參數(shù)構(gòu)建多維特征向量數(shù)據(jù)集，通過 PNN網(wǎng)絡(luò)建立磨削聲發(fā)射信號(hào)與砂輪鈍化狀態(tài)的關(guān)系。經(jīng)過測(cè)試，該方法可以較為準(zhǔn)確識(shí)別砂輪鈍化狀態(tài)，并能橫向?qū)Ρ炔煌卣鲄?shù)的檢測(cè)效果。

1 方法理論

1.1 變分模態(tài)分解(VMD)

VMD算法基于維納(Wiener)濾波和希爾伯特(Hilbert)變換方法，能夠?qū)⑿盘?hào)自適應(yīng)分解為 k個(gè)稀疏的本征模態(tài)函數(shù)(IMF)，第k個(gè)本征模態(tài)函數(shù)表示為

K為分解層數(shù)，然后利用交替方向乘法算子算法(Alternate Direction Method of Multipliers, ADMM)求解式(5)的最優(yōu)解，在頻域內(nèi)迭代更新uk、ωk和λ。

最后重復(fù)式(6)～(8)，直到滿足式(9)中的終止條件，并輸出最后結(jié)果。

式中：ε為設(shè)定的終止參數(shù)。

1.2 聲發(fā)射特征參數(shù)選取

聲發(fā)射信號(hào)是有許多聲發(fā)射事件(AE Hit)組成的，從波形首次越過門限值開始，直到落回門限值以下并且一段時(shí)間內(nèi)(該時(shí)間設(shè)置為撞擊定義時(shí)間)不再升到門限值以上，定義為一次AE Hit。一個(gè)典型的AE Hit波形如圖2所示。

圖2 典型聲發(fā)射信號(hào)波形圖Fig.2 Waveform of typical acoustic emission signal

對(duì)于每次AE Hit，有多種常用特征參數(shù)來進(jìn)行描述。幅度(Amplitude)指AE Hit達(dá)到的最大振幅；持續(xù)時(shí)間(Duration)指一次AE Hit從開始到結(jié)束的間隔；振鈴計(jì)數(shù)(ring-down counts, 記為C)，表示一次AE hit越過門限的震蕩次數(shù)；能量計(jì)數(shù)(SE)，表示一個(gè)AE hit的波形包絡(luò)面積；均方根值(RMS)，指的是AE hit信號(hào)幅度的均方根值，表達(dá)式為

式中：u( t)為信號(hào)幅度，以上特征參數(shù)都應(yīng)用于描述單次AE Hit，在一段聲發(fā)射信號(hào)中，需計(jì)算所有AE Hit對(duì)應(yīng)特征參數(shù)并求均值作為該段聲發(fā)射信號(hào)的特征參數(shù)，因此這些特征參數(shù)不能衡量AE Hit的密度和聲發(fā)射信號(hào)整體狀況。

文獻(xiàn)[14]提出一種通過功率比來評(píng)價(jià)聲發(fā)射信號(hào)的特征量，聲發(fā)射功率比：

式(11)是RP頻域和時(shí)域的兩種表示方法，N1和N2為聲發(fā)射的頻率范圍，在本文中RP含義為聲發(fā)射信號(hào)功率與原始信號(hào)功率之比。RP可以排除信號(hào)隨時(shí)間變化引起的功率值整體波動(dòng)影響，適用于衡量聲發(fā)射信號(hào)在時(shí)域的變化規(guī)律。

參考能量計(jì)數(shù)的定義原理，用希爾伯特變換對(duì)信號(hào)求包絡(luò)，在一段信號(hào)中，聲發(fā)射分量包絡(luò)能量和原始信號(hào)包絡(luò)能量的比值定義為聲發(fā)射包絡(luò)能量比REE：

式中： xAE(t)是信號(hào)x( t)經(jīng)過分解后的聲發(fā)射分量,H[ x ( t)]是信號(hào)的希爾伯特變換。

REE含義為聲發(fā)射信號(hào)能量占比，這里的聲發(fā)射能量是一種結(jié)合聲發(fā)射信號(hào)特征定義的包絡(luò)能量，含義為信號(hào)檢波包絡(luò)線下的面積，量綱為mV·ms，與物理意義上的能量不同。

特征參數(shù)選取是聲發(fā)射檢測(cè)的重要環(huán)節(jié)，以上聲發(fā)射特征參數(shù)可以分為兩類，一類是對(duì)單次 AE hit進(jìn)行評(píng)價(jià)，分別從信號(hào)電平、時(shí)間和包絡(luò)能量(電平×?xí)r間)三個(gè)維度進(jìn)行衡量，在計(jì)算時(shí)需對(duì)一段信號(hào)所有AE hit分別計(jì)算并求平均值；另一類為聲發(fā)射信號(hào)與原始信號(hào)的功率或包絡(luò)能量比值，是對(duì)聲發(fā)射信號(hào)進(jìn)行整體衡量，上文所提到的聲發(fā)射特征參數(shù)分類如表1所示。其中幅度和持續(xù)時(shí)間易受異常聲發(fā)射事件影響，波動(dòng)較大，結(jié)合此前相關(guān)研究基礎(chǔ)，本文選取C、SE、RMS、RP和REE作為砂輪鈍化狀態(tài)識(shí)別參考特征參數(shù)。

表1 聲發(fā)射特征參數(shù)分類Table 1 Classification of AE characteristic parameters

1.3 概率神經(jīng)網(wǎng)絡(luò)(PNN)

概率神經(jīng)網(wǎng)絡(luò)是一種建立在徑向基(Radial Basis Function, RBF)網(wǎng)絡(luò)的基礎(chǔ)上，基于貝葉斯決策理論，結(jié)合密度函數(shù)估計(jì)方法的神經(jīng)網(wǎng)絡(luò)[15]。概率神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)如圖3所示。

圖3 PNN神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.3 PNN neural network structure

結(jié)構(gòu)共分為四層，第一層為輸入層，接收訓(xùn)練數(shù)據(jù)集的值，神經(jīng)元個(gè)數(shù)與輸入向量長(zhǎng)度相等，并

傳遞給第二層。第二層是隱含層，該層接收輸入層的數(shù)據(jù)，并計(jì)算輸入向量與中心的距離，輸出一個(gè)標(biāo)量，第i類模式的第j神經(jīng)元的輸入輸出關(guān)系為

2 實(shí)驗(yàn)與結(jié)果分析

2.1 實(shí)驗(yàn)系統(tǒng)搭建

搭建聲發(fā)射檢測(cè)及砂輪鈍化狀態(tài)識(shí)別實(shí)驗(yàn)系統(tǒng)，包括三個(gè)組成模塊，如圖4所示。

圖4 聲發(fā)射檢測(cè)實(shí)驗(yàn)系統(tǒng)Fig.4 Experimental system for acoustic emission detection

砂輪選用 36粒綠色碳化硅砂輪，線速度為42 m·s-1，進(jìn)刀量為0.03 mm；加工工件為合金鑄鋼；傳感器為全接觸式AE傳感器，固介靈敏度較高，固定在磨床砂輪外殼上，響應(yīng)范圍為 0～3 MHz；AE傳感器接收到的信號(hào)由安捷倫u2331數(shù)據(jù)采集設(shè)備采集，存儲(chǔ)為二進(jìn)制文件，再通過計(jì)算機(jī)進(jìn)行信號(hào)分析處理。

2.2 VMD參數(shù)設(shè)定

信號(hào)處理的關(guān)鍵是利用VMD方法對(duì)原始信號(hào)進(jìn)行分解得到聲發(fā)射信號(hào)，VMD有兩個(gè)重要參數(shù)α和K，α是懲罰因子，決定分解出每個(gè)分量的帶寬限制，通常取默認(rèn)值2 000；分解層數(shù)K的確定方法有多種，這里采取峭度判別法。相比于噪聲信號(hào)，聲發(fā)射信號(hào)波形陡峭，峭度值較大。逐步增大K，觀察每個(gè)分量的波形和峭度，直到出現(xiàn)峭度值較小的非聲發(fā)射分量為止，此時(shí)分解層數(shù)為適宜的K值，取其中峭度較大的分量重構(gòu)，即得到聲發(fā)射信號(hào)。VMD各模態(tài)的中心頻率及峭度分布如表2、3所示。VMD 5層分解的IMF波形及頻譜圖如圖5所示。

從表2、表3及圖5觀察到，隨著分解層數(shù)K的增大，各分量的中心頻率基本保持不變，在分解層數(shù)達(dá)到5時(shí)出現(xiàn)中心頻率大于1 000 kHz而峭度較低的電氣噪聲分量，故分解層數(shù)K定為5，并判斷IMF2、IMF3和IMF4是聲發(fā)射信號(hào)分量，將這三個(gè)分量重構(gòu)即得到所需的聲發(fā)射信號(hào)。

表2 VMD分解各模態(tài)中心頻率分布表Table 2 Modal center frequencies of VMD decomposition

表3 VMD分解各模態(tài)峭度分布表Table 3 Modal kurtosis distribution of VMD decomposition

圖5 VMD 5層分解IMF波形圖及頻譜圖Fig.5 IMF waveforms and spectrums of 5 layer VMD decomposition

讀取三段已標(biāo)記砂輪鈍化狀態(tài)的信號(hào)數(shù)據(jù)波形文件并進(jìn)行5層VMD分解，取IMF2、IMF3和IMF4進(jìn)行重構(gòu)，觀察不同砂輪鈍化狀態(tài)下原始信號(hào)和聲發(fā)射信號(hào)波形。圖 6～8分別為初期鈍化、中期鈍化、嚴(yán)重鈍化三個(gè)階段的原始信號(hào)和 VMD分解后重構(gòu)的聲發(fā)射信號(hào)波形圖及對(duì)應(yīng)頻譜。

圖6 初期鈍化原始信號(hào)與聲發(fā)射信號(hào)波形及頻譜圖Fig.6 The waveforms and spectrums of the original and the VMD reconstructed AE signals in initial blunting state

圖7 中期鈍化原始信號(hào)與聲發(fā)射信號(hào)波形及頻譜圖Fig.7 The waveforms and spectrums of the original and the VMD reconstructed AE signals in intermediate blunting state

圖8 嚴(yán)重鈍化原始信號(hào)與聲發(fā)射信號(hào)波形及頻譜圖Fig.8 The waveforms and spectrums of the original and the VMD reconstructed AE signals in severe blunting state

可以看出，磨削加工過程中砂輪處于初期鈍化階段時(shí)，聲發(fā)射事件強(qiáng)度較弱，分布稀疏，20 kHz以上頻段幅值較??；砂輪進(jìn)入中期鈍化階段以及嚴(yán)重鈍化階段，單次聲發(fā)射事件強(qiáng)度增大，聲發(fā)射信號(hào)頻譜幅值增大，故單位時(shí)間內(nèi)聲發(fā)射信號(hào)與原始信號(hào)的能量比也大幅增加。

2.3 數(shù)據(jù)集訓(xùn)練與識(shí)別結(jié)果分析

信號(hào)處理模塊的全過程如圖9所示，取磨削過程中已明確砂輪鈍化狀態(tài)的初期鈍化信號(hào)40 s，中期鈍化信號(hào)40 s，嚴(yán)重鈍化信號(hào)20 s。對(duì)信號(hào)分別進(jìn)行分幀處理，以 0.1 s為一幀，對(duì)每幀信號(hào)進(jìn)行VMD分解、重構(gòu)，得到1 000幀聲發(fā)射信號(hào)。設(shè)定門限值為 300 mV，計(jì)算每幀聲發(fā)射信號(hào)所有正常聲發(fā)射事件的C、SE、RMS，并在一幀內(nèi)求平均值得到平均C、平均SE、平均RMS、對(duì)每幀聲發(fā)射信號(hào)進(jìn)行整體計(jì)算得到RP和REE，共五個(gè)特征參數(shù)，構(gòu)建一組五維特征向量數(shù)據(jù)集作為訓(xùn)練集，將訓(xùn)練集作歸一化處理后輸入PNN網(wǎng)絡(luò)中訓(xùn)練。圖9中將C、SE、RMS、RP、REE5個(gè)特征參數(shù)分別標(biāo)注為①，②，③，④，⑤。

圖9 多維特征向量數(shù)據(jù)集的構(gòu)建過程Fig.9 Construction process of multi-dimensional characteristic vector dataset

改變信號(hào)分幀的起始位置，將訓(xùn)練集每幀信號(hào)后移0.05 s，幀長(zhǎng)依然為0.1 s，得到的數(shù)據(jù)集與訓(xùn)練集是同一段信號(hào)因不同分幀方法所得，因此被認(rèn)為是獨(dú)立同分布的。將該數(shù)據(jù)集作為測(cè)試集，隨機(jī)取三種狀態(tài)各80、80、40幀，共200幀進(jìn)行測(cè)試。

測(cè)試結(jié)果用混淆矩陣表達(dá)，如圖 10所示，初期鈍化狀態(tài)、中期鈍化狀態(tài)和嚴(yán)重鈍化狀態(tài)分別被標(biāo)記為類別1、2和3。同時(shí)利用五種特征進(jìn)行訓(xùn)練和測(cè)試，整體識(shí)別準(zhǔn)確率為94.5%，多次重復(fù)測(cè)試集取樣過程，平均準(zhǔn)確率維持在94%以上，且嚴(yán)重鈍化狀態(tài)對(duì)應(yīng)的數(shù)據(jù)均能被正確判定為嚴(yán)重鈍化?？梢哉J(rèn)為該神經(jīng)網(wǎng)絡(luò)經(jīng)過訓(xùn)練后能夠根據(jù)特征參數(shù)對(duì)砂輪鈍化狀態(tài)進(jìn)行識(shí)別，并對(duì)嚴(yán)重鈍化狀態(tài)作出準(zhǔn)確預(yù)警。

圖10 混淆矩陣，砂輪鈍化識(shí)別結(jié)果Fig.10 Confusion matrix and the identification result of grinding wheel blunting state

為評(píng)估該神經(jīng)網(wǎng)絡(luò)的泛用性，取另一段加工材料為合金鍛鋼支承輥的400 s連續(xù)磨削加工信號(hào)，設(shè)定幀長(zhǎng)為0.1 s，幀移為1 s，即用每1 s信號(hào)中間0.1 s的特征參數(shù)來表征這1 s的砂輪狀態(tài)。得到該段連續(xù)磨削加工信號(hào)的數(shù)據(jù)集作為連續(xù)信號(hào)數(shù)據(jù)集，圖 11為該數(shù)據(jù)集的五維特征參數(shù)變化圖。對(duì)連續(xù)信號(hào)數(shù)據(jù)集作歸一化處理并輸入已通過此前訓(xùn)練集訓(xùn)練好的PNN網(wǎng)絡(luò)，得到識(shí)別結(jié)果如圖12所示。連續(xù)信號(hào)數(shù)據(jù)集被清晰分成了三種鈍化狀態(tài)，不同顏色代表不同鈍化狀態(tài)，且隨時(shí)間推移存在階段性變化。在兩個(gè)鈍化狀態(tài)之間存在短暫的過渡態(tài)，根據(jù)加工要求和材料特點(diǎn)，當(dāng)連續(xù)若干幀信號(hào)有一定幀數(shù)信號(hào)被識(shí)別為嚴(yán)重鈍化信號(hào)，可認(rèn)為砂輪已經(jīng)進(jìn)入嚴(yán)重鈍化狀態(tài)，此時(shí)需要對(duì)砂輪進(jìn)行修整。本文方法應(yīng)用在其他加工工件上依然能夠?qū)ι拜嗏g化狀態(tài)進(jìn)行有效識(shí)別。

圖11 連續(xù)磨削加工信號(hào)的5維特征參數(shù)數(shù)據(jù)集隨時(shí)間變化圖Fig.11 Variations of the 5D characteristic parameter dataset of continuous grinding signal with time

圖12 砂輪鈍化狀態(tài)判別結(jié)果Fig.12 Identifying results of grinding wheel blunting states

2.4 不同特征參數(shù)識(shí)別準(zhǔn)確率比較

探究不同特征參數(shù)對(duì)于識(shí)別準(zhǔn)確性的影響，對(duì)特征參數(shù)進(jìn)行分組，分組采用部分特征參數(shù)進(jìn)行訓(xùn)練。5種特征參數(shù)分別標(biāo)注為①，②，③，④，⑤，按表4分為15種組合方式。

表4 特征參數(shù)分組表Table 4 Grouping table for characteristic parameters

對(duì)不同特征參數(shù)組合的訓(xùn)練集和數(shù)據(jù)集分別進(jìn)行訓(xùn)練與測(cè)試，結(jié)果如圖13和表5所示。其中，組合5因?yàn)橥瑫r(shí)用到了五種特征參數(shù)準(zhǔn)確率最高。單個(gè)特征參數(shù)的識(shí)別準(zhǔn)確率中，組合 10(RMS)的準(zhǔn)確度最高，達(dá)到 91.0%；其次是組合 15(REE)和組合 13(RP)，單個(gè)特征參數(shù)的準(zhǔn)確度達(dá)到 88.0%和85.5%。

圖13 特征參數(shù)組合測(cè)試正確率對(duì)比Fig.13 Comparison of correct rates between different characteristic parameter combination tests

表5 特征參數(shù)單獨(dú)訓(xùn)練測(cè)試結(jié)果表Table 5 Test results of individual characteristic parameter training

由此可得，不同聲發(fā)射特征參數(shù)對(duì)砂輪鈍化狀態(tài)識(shí)別效果存在差異，使用多個(gè)特征參數(shù)的測(cè)試準(zhǔn)確率高于使用其中部分特征參數(shù)，RMS、REE等典型特征參數(shù)與砂輪鈍化狀態(tài)相關(guān)性較強(qiáng)。

進(jìn)一步探究?jī)深惵暟l(fā)射特征參數(shù)在不同鈍化階段的識(shí)別準(zhǔn)確率差異，對(duì)比組合3和組合14，結(jié)果如表6。

表6 兩組特征參數(shù)在不同鈍化階段的識(shí)別準(zhǔn)確率Table 6 Identification accuracy of two groups of characteristic parameters in different blunting stages

組合3特征參數(shù)(C、SE和RMS)衡量單次AE Hit強(qiáng)度變化，組合14特征參數(shù)(RP和REE)整體衡量聲發(fā)射信號(hào)能量大小。從砂輪鈍化原因來看，初期鈍化主要是磨粒的破碎，AE hit密度增加而單次AE hit強(qiáng)度變化較?。恢衅阝g化以磨粒磨損和塑性變形為主，磨削力和接觸面積增大，AE hit強(qiáng)度有明顯增加，使得第一類特征參數(shù)(組合1)在中期鈍化階段識(shí)別準(zhǔn)確率高于初期鈍化階段；嚴(yán)重鈍化階段 AE hit強(qiáng)度和AE hit密度都大幅增加導(dǎo)致區(qū)分明顯，兩組識(shí)別準(zhǔn)確率都很高。由此可知，聲發(fā)射特征參數(shù)變化與砂輪鈍化原因存在對(duì)應(yīng)關(guān)系，通過針對(duì)性選取典型特征參數(shù)可以較準(zhǔn)確地識(shí)別砂輪在特定階段的鈍化情況。

3 結(jié) 論

本文利用VMD方法，將磨削加工過程中產(chǎn)生的聲發(fā)射信號(hào)分離出來，保留了聲發(fā)射信號(hào)的主要特征并避免了噪聲干擾。聲發(fā)射特征參數(shù)變化存在一定的波動(dòng)性，PNN能夠較為準(zhǔn)確地將特征參數(shù)對(duì)應(yīng)到砂輪的不同鈍化狀態(tài)，測(cè)試準(zhǔn)確率較高(94.5%)，且經(jīng)過泛用性評(píng)估，實(shí)現(xiàn)了對(duì)砂輪鈍化狀態(tài)的識(shí)別和嚴(yán)重鈍化狀態(tài)預(yù)警。

聲發(fā)射信號(hào)特征參數(shù)的選取是影響識(shí)別準(zhǔn)確率的關(guān)鍵，本文在聲發(fā)射檢測(cè)相關(guān)研究基礎(chǔ)上對(duì)特征參數(shù)進(jìn)行歸類，參考聲發(fā)射能量的定義提出聲發(fā)射包絡(luò)能量比(REE)，共選取5個(gè)典型特征參數(shù)，通過分組訓(xùn)練并比較測(cè)試準(zhǔn)確率的方式，證實(shí)了所選特征參數(shù)對(duì)砂輪鈍化狀態(tài)有較高的識(shí)別準(zhǔn)確率，并驗(yàn)證了典型特征參數(shù)變化與砂輪不同鈍化原因存在相關(guān)性，可為實(shí)際聲發(fā)射檢測(cè)特征參數(shù)的選取提供參考。