基于希爾伯特黃變換的刀具磨損特征提取

第一作者 孫惠斌 男，博士，副教授，1977年生

通信作者 牛偉龍 男，碩士，1989年生

基于希爾伯特黃變換的刀具磨損特征提取

孫惠斌，牛偉龍，王俊陽

(西北工業(yè)大學(xué) 機(jī)電學(xué)院，西安710072)

摘要：概述了希爾伯特黃變換(HHT)的基本理論和算法，對信號經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)后得到的固有模態(tài)函數(shù)(IMF)求取振幅均值，差值篩選出與刀具磨損相關(guān)的IMF分量，并對單分量固有模態(tài)函數(shù)求取邊際譜，獲取邊際譜最大幅值點(diǎn)，建立他們與刀具磨損之間的映射關(guān)系，進(jìn)行特征提取，將其作為神經(jīng)網(wǎng)絡(luò)的輸入特征向量，結(jié)合希爾伯特三維時(shí)頻譜進(jìn)行刀具磨損狀態(tài)的判斷。研究結(jié)果證明，該方法可以作為刀具磨損監(jiān)測中信號特征提取的一種簡單和可靠的方法。

關(guān)鍵詞：希爾伯特黃變換；小波去噪；固有模態(tài)函數(shù)；希爾伯特譜；邊際譜

基金項(xiàng)目：陜西省自然科學(xué)基金(2013JM7001);西北工業(yè)大學(xué)基礎(chǔ)研究基金( JC20110215);西北工業(yè)大學(xué)研究生創(chuàng)業(yè)種子基金

收稿日期：2013-11-01修改稿收到日期：2014-03-03

中圖分類號:TH165；TN911文獻(xiàn)標(biāo)志碼： A

Tool wear feature extraction based on Hilbert-Huang transformation

SUNHui-bin,NIUWei-long，WANGJun-yang(School of Mechanical Engineering, Northwestern Polytechnical University, Xi’an 710072, China)

Abstract:After presenting the basic theory and algorithm of Hilbert-Huang transformation (HHT), a tool signal was decomposed with the empirical mode decomposition (EMD) method and its intrinsic mode functions (IMFs) were gained to obtain their average amplitude. The IMF components related to tool wear were chosen using a difference screen. Meanwhile, the marginal spectrum of a single intrinsic mode function was obtained and its maximum amplitude point was then found. By establishing the mapping relationship between maximum amplitude points and tool wear, the features of tool wear were extracted. Regarding them as input eigen-vectors of a neural network, and combined with Hilbert spectra, the tool wear status judgment was implemented. The study results showed that this approach is a simple and reliable method to detect the level of tool wear.

Key words: Hilbert-Huang transformation; wavelet domain denoising; intrinsic mode function (IMF); Hilbert spectrum; marginal spectrum

刀具是加工過程中最常見也是最重要的加工要素之一，其磨損狀態(tài)直接影響了加工工件的質(zhì)量和精度。研究刀具磨損狀態(tài)的監(jiān)測具有重大的意義，不僅可以大大提高機(jī)械加工效率,而且可以降低加工成本，有著較大的經(jīng)濟(jì)效益。由于直接監(jiān)測刀具磨損難度大，間接監(jiān)測刀具磨損已成為一個(gè)重要研究方向。其主要技術(shù)路線是：

(1) 信號采集：通過力、加速度、聲發(fā)射、功率等一些相應(yīng)的傳感器對刀具加工過程中的一些信號進(jìn)行采集。

(2) 信號處理：通過一些現(xiàn)有的信號處理方法(如：傅里葉變換、小波變換[1]等)對采集的信號進(jìn)行分析。

(3) 特征提?。盒盘柼幚硗戤吅髮⒌毒哒Ｄp和嚴(yán)重磨損狀態(tài)下信號處理結(jié)果進(jìn)行對比，找出它們之間的不同，并與磨損狀態(tài)建立相應(yīng)的映射關(guān)系。

(4) 智能決策：根據(jù)提取的特征對刀具磨損狀態(tài)進(jìn)行識別和判斷，如模式識別、專家系統(tǒng)、神經(jīng)網(wǎng)絡(luò)。

其中，信號處理是一個(gè)關(guān)鍵環(huán)節(jié)，它可以將刀具磨損從數(shù)據(jù)中進(jìn)行定量描述，信號處理方法也決定了可提取到的磨損特征，以及智能決策的可靠性。但傳統(tǒng)的傅里葉變換、短時(shí)傅里葉變換、小波變換等信號處理方法存在以下不足：

(1) 傅里葉變換只能處理線性非平穩(wěn)的信號，小波變換雖然在理論上能處理非線性非平穩(wěn)信號，但在實(shí)際算法實(shí)現(xiàn)中卻只能處理線性非平穩(wěn)信號。

(2) 不具有自適應(yīng)性，傅里葉變換的基是三角函數(shù)，小波變換的基是滿足“可容性條件”的小波基。①一旦選定了基函數(shù)，那么就必須用這種基函數(shù)來分析所有的數(shù)據(jù)；②選定了分解尺度，所得到的結(jié)果是某一固定頻率段的時(shí)域波形，所包含的頻率只與信號的分析頻率有關(guān)，而與信號本身無關(guān)。不能準(zhǔn)確的反映被分析信號的特性。

(3) 受Heisenberg測不準(zhǔn)原理制約，即時(shí)間窗口與頻率窗口的乘積為一個(gè)常數(shù)。這就意味著如果要提

高時(shí)間精度就得犧牲頻率精度，反之亦然，不能在時(shí)間和頻率同時(shí)達(dá)到很高的精度。

(4) 需要預(yù)先選擇基函數(shù)，其計(jì)算方式是通過與基函數(shù)的卷積產(chǎn)生的,這就導(dǎo)致了傅里葉變換的頻率是全局性的，小波變換是區(qū)域性的，不能對信號局部特征進(jìn)行精確的分析。

希爾伯特黃變換則不同于這些傳統(tǒng)方法，它徹底擺脫了線性和平穩(wěn)性束縛[2]，其適用于分析非線性非平穩(wěn)信號，能夠自適應(yīng)產(chǎn)生“基”，即由“篩選”過程產(chǎn)生的IMF為它的基函數(shù)。同時(shí)也不受Heisenberg測不準(zhǔn)原理制約，可以在時(shí)間和頻率同時(shí)達(dá)到很高的精度，并借助Hilbert變換求得相位函數(shù)，再對相位函數(shù)求導(dǎo)產(chǎn)生瞬時(shí)頻率,這樣求出的瞬時(shí)頻率是局部性的，對局部特征反映更為準(zhǔn)確。因此，希爾伯特黃變換從很大程度上彌補(bǔ)了傳統(tǒng)信號分析方法的不足，使得它比傳統(tǒng)信號的時(shí)頻分析更優(yōu)越，因而具有更好的應(yīng)用前景和更廣的應(yīng)用空間[3-4]。這里就希爾伯特黃變換對刀具磨損監(jiān)測過程中振動(dòng)信號的特征提取進(jìn)行了分析和實(shí)例研究。

1希爾伯特黃變換

1.1EMD方法基本原理

Huang等[5]首次提出希爾伯特黃變換(Hilbert-Huang Transform HHT)，是一種全新的信號時(shí)頻分析理論，在這一理論中，引入了固有模態(tài)函數(shù)(Intrinsic Mode Function, IMF)[6]的概念,在經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)的基礎(chǔ)上[7]，對每個(gè)IMF進(jìn)行Hilbert變換得到瞬時(shí)頻率[8]，從而將信號精確表示成頻率-時(shí)間-能量(或者幅度)的分布，稱為Hilbert譜。HHT的特點(diǎn)在于它是基于信號局部特征，能夠自適應(yīng)的高效分解,而且它特別適用于分析非線性非平穩(wěn)信號。

EMD是HHT的一個(gè)關(guān)鍵步驟,首先找到信號x(t)x(t)的極大值和極小值,通過三次樣條擬合,從而得到信號的上包絡(luò)線U(t)和下包絡(luò)線L(t),計(jì)算上、下包絡(luò)線的均值m1，然后再計(jì)算x(t)-m1=h1。

判斷h1是否為一個(gè)IMF的標(biāo)準(zhǔn)是：

(1) 在整個(gè)數(shù)據(jù)段內(nèi)，極值點(diǎn)的個(gè)數(shù)和過零點(diǎn)的個(gè)數(shù)必須相等或最多相差一個(gè)；

(2) 在任意時(shí)刻，由局部極大值點(diǎn)形成的上包絡(luò)線和由局部極小值點(diǎn)形成的下包絡(luò)線的平均值為零，也就是說，上下包絡(luò)線相對于時(shí)間軸局部對稱。理想情況下，如果h1是一個(gè)IMF，那么h1就是x(t)的第1個(gè)IMF分量；如果h1不滿足IMF的條件，把h1作為原始數(shù)據(jù)，重復(fù)上述分解過程，得到上下包絡(luò)線的平均值m11，再判斷h11=h1-m11是否滿足IMF的條件，如不滿足，則重復(fù)循環(huán)k次，得到h1(k-1)-m1k=h1k使得h1k滿足IMF的條件。記c1=h1k，則c1為信號x(t)的第1個(gè)滿足IMF條件的分量。

將c1從x(t)中分離出來，得到：

r1=x(t)-c1

(1)

將r1作為原始數(shù)據(jù)，重復(fù)上述步驟，依次篩選出n個(gè)滿足IMF定義的分量(r1,r2,…,rn)，當(dāng)rn為一個(gè)單調(diào)函數(shù)或只有一個(gè)極值時(shí)停止循環(huán)，于是信號可以被寫成：

(2)

式中:rn為殘余函數(shù)，代表信號的平均趨勢。

1.2Hilbert變換

對EMD 分解得到的每個(gè)IMF 分量應(yīng)用Hilbert 變換，Hilbert變換相當(dāng)于將一個(gè)信號x(t)與1/πt做卷積：

(3)

式中:P為柯西主分量，得到的解析信號為：

Z(t)=X(t)+iY(t)=a(t)eiθ(t)

(4)

其中

(5)

(6)

則瞬時(shí)頻率[9]ω(t)定義為：

(7)

由于EMD分解得到的IMF均為單分量信號，故對其進(jìn)行Hilbert變換可得到信號的瞬時(shí)頻率，瞬時(shí)相位和瞬時(shí)幅值。

1.3希爾伯特譜和邊際譜

把式(2)～式(6)所表示的變換用于每個(gè)固有模態(tài)函數(shù)序列，數(shù)據(jù)便可表示為下式：

(8)

(9)

(10)

把三維時(shí)頻譜經(jīng)過對時(shí)間的積分，便形成了只有頻率和幅值的二維譜圖。這種由HHT得到的邊際譜與Fourier頻譜有相似之處。從統(tǒng)計(jì)觀點(diǎn)上來看，HHT邊際譜表示了該頻率上振幅(能量)在時(shí)間上的累加，能夠反映各頻率上的能量分布，但因?yàn)樗矔r(shí)頻率定義為時(shí)間的函數(shù)，不同以往Fourier等需要完整的振蕩波周期來定義局部的頻率值，而且求取的能量值不是全局定義的，因此HHT邊際譜對信號的局部特征反映更為準(zhǔn)確。尤其是在分析非平穩(wěn)信號時(shí)，這種定義對于頻率隨時(shí)間隨時(shí)變化的信號特征來說，能夠最真實(shí)地反映振動(dòng)特點(diǎn)。而Fourier頻譜的某一點(diǎn)頻率上的幅值則表示的是在整個(gè)信號里有一個(gè)含有此頻率的三角函數(shù)組分。因此，邊際譜和傅氏譜的意義是不同的。

在刀具切削加工中，刀具的磨損是一個(gè)比較復(fù)雜的過程，而對刀具磨損的狀態(tài)監(jiān)測往往需要將刀具在正常磨損狀態(tài)和嚴(yán)重磨損狀態(tài)進(jìn)行一定意義上的區(qū)分，即對信號進(jìn)行特征提取，提取到的特征則要反映刀具磨損的狀態(tài)。而選用希爾伯特黃變換，它能處理非線性非平穩(wěn)信號[11]并具有較強(qiáng)的自適應(yīng)性，而對于它經(jīng)過EMD分解后所得每個(gè)本征模函數(shù)分量IMF都有相應(yīng)的物理意義與其對應(yīng)，這更方便于有針對性的對其進(jìn)行研究與分析。除此之外，希爾伯特譜與希爾伯特邊際譜也都可以在一定程度上反映出刀具磨損狀態(tài)的特征。

2刀具磨損的特征提取方法

2.1IMF的振幅均值法

在實(shí)際應(yīng)用領(lǐng)域中得到的信號總是混雜著一定的噪聲，而噪聲的存在嚴(yán)重干擾了信號的本來面目，不利于進(jìn)一步的信號分析和處理。因此，在信號處理過程中對噪聲加以消除或減少，以便最大程度地提取出有用的信號。為了增加研究分析的精確性和可靠性，先對采集的振動(dòng)信號進(jìn)行小波閾值去噪[12-13]，再對其進(jìn)行希爾伯特分析。

經(jīng)過小波去噪和EMD分解后，任何復(fù)雜的信號都可分解為為數(shù)不多的IMF且相互獨(dú)立，每個(gè)IMF分量也都有相應(yīng)的物理意義與其對應(yīng)。因此，可以對IMF分量求振幅均值，除去那些與刀具磨損不相關(guān)的IMF分量，將加工過程中刀具與工件接觸時(shí)所對應(yīng)的IMF分量提取出，與刀具磨損狀態(tài)建立相應(yīng)的映射關(guān)系，作為刀具磨損的特征，由于離散采樣，則振幅均值為：

(11)

選擇與磨損相關(guān)的IMF分量時(shí)，可以利用差值方法來進(jìn)行篩選，在三種不同磨損狀態(tài)下，加工參數(shù)中只有刀具磨損量改變較大，其它的加工參數(shù)基本不變，因此只有刀具與工件接觸時(shí)所對應(yīng)的IMF分量會(huì)隨著刀具磨損的加劇而產(chǎn)生明顯的變化，即三種狀態(tài)下它們的振幅均值的差值比較大。實(shí)驗(yàn)表明，當(dāng)它們的振幅差值大于1.5時(shí)，都可以認(rèn)為它們是與刀具磨損相關(guān)的IMF分量的振幅均值：

(12)

(13)

式中:IMFki、IMFkm、IMFkn分別表示的是嚴(yán)重磨損、正常磨損、初期磨損下的IMF分量的振幅均值。

邊際譜可以比較準(zhǔn)確的反應(yīng)信號的實(shí)際頻率，由于在加工過程中受到環(huán)境等因素的干擾，經(jīng)過EMD分解后的IMF分量中有些反映的是與刀具加工過程無關(guān)的信號。因此，只要選出能反映加工時(shí)刀具切削工件所對應(yīng)的IMF，該IMF包含的能量和信息是最多的，單獨(dú)對其求邊際譜即可，設(shè)第k個(gè)IMF分量滿足條件，即反應(yīng)刀具切削工件時(shí)所包含最多的能量和信息，用于式(2)～式(6)所表示的變換，則它的邊際譜：

I(t)=ak(t)ejθk(t) =ak(t)ej∫ωk(t)dt

(14)

H(ω,t)=ak(t)ej∫ωk(t)dt

(15)

(16)

對于它的邊際譜，由于磨損狀態(tài)的不同，它們的最大幅值點(diǎn)相差明顯，故可以將其最大幅值點(diǎn)作為磨損特征進(jìn)行提取。

希爾伯特三維時(shí)頻幅值譜，表示在不同頻率和時(shí)間點(diǎn)下的能量振幅，可以根據(jù)它在不同磨損狀態(tài)下表現(xiàn)出來的頻率幅值的變化以及波動(dòng)幅度大小作為判斷刀具磨損狀態(tài)的輔助參考。

刀具磨損監(jiān)測過程其實(shí)也是一個(gè)智能決策[14]的過程，而人工神經(jīng)網(wǎng)絡(luò)是智能決策中一種十分有效的方法。通過神經(jīng)網(wǎng)絡(luò)建立信號特征和刀具磨損狀態(tài)之間的對應(yīng)關(guān)系來識別刀具的磨損狀態(tài),是目前刀具磨損識別研究者廣泛采用的一種方法，對刀具磨損特征提取完畢后，可以將上述特征作為神經(jīng)網(wǎng)絡(luò)的輸入特征向量，對神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練和識別。

圖1為基于希爾伯特黃變換的刀具磨損監(jiān)測系統(tǒng)流程，圖1中IMFk1、IMFk2、IMFk3為加工過程中刀具與工件接觸時(shí)所對應(yīng)的IMF分量的振幅均值，IMFks為包含能量和信息最多的固有模態(tài)函數(shù)，求其邊際譜最大幅值點(diǎn)，將他們作為神經(jīng)網(wǎng)絡(luò)的輸入特征向量，進(jìn)行神經(jīng)網(wǎng)絡(luò)的訓(xùn)練和狀態(tài)識別，隨后反饋給計(jì)算機(jī)，計(jì)算機(jī)再根據(jù)反饋結(jié)果對機(jī)床進(jìn)行自動(dòng)換刀。

圖1　基于希爾伯特黃變換的刀具磨損監(jiān)測系統(tǒng) Fig.1 Tool wear monitoring system based on Hilbert Huang Transform

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

本實(shí)驗(yàn)采用C000031-CNC4240小型數(shù)控銑床，材料采用45#調(diào)制鋼，硬度：HB200，尺寸：70 mm×50 mm× 50mm，刀具采用蘇州AHNO生產(chǎn)的型號為Eco-U13φ×25 mm×75 mm的4刃整體HM平頭硬質(zhì)合金立銑刀。傳感器采用KISTLER公司生產(chǎn)的型號為8702B100M1單軸加速度計(jì)(±100 g，與地絕緣，靈敏度50 mV/g)，對Z軸方向進(jìn)行振動(dòng)信號采集。除此，本實(shí)驗(yàn)還采用了便攜式數(shù)據(jù)采集系統(tǒng)主機(jī)，DEWSOFT公司生產(chǎn)的DEWSOFT-43，并設(shè)置采樣頻率為20 kHz，實(shí)驗(yàn)平臺(見圖2)。

圖2　實(shí)驗(yàn)平臺 Fig.2 Experiment platform

在銑削過程中，刀齒從極薄的切削層下切入工件，有擠壓的滑行并伴隨著高溫高壓，銑刀的前刀面和后刀面經(jīng)常與工件表面發(fā)生強(qiáng)烈的摩擦，因此生產(chǎn)中以后刀面磨損寬度VB作為判斷銑刀磨損程度的一個(gè)標(biāo)準(zhǔn)(見圖3 )，實(shí)驗(yàn)規(guī)定刀具初期磨損量VB為0～0.05 mm，正常磨損量VB為0.05～0.3 mm，當(dāng)大于0.3 mm時(shí)為嚴(yán)重磨損。

圖3　刀具磨損的測量 Fig.3 Measurement of tool wear

為了全面和準(zhǔn)確分析，找出最佳水平組合，使得指標(biāo)最優(yōu)，且要減少試驗(yàn)次數(shù)和試驗(yàn)復(fù)雜程度。實(shí)驗(yàn)采用正交試驗(yàn)方法[15]，根據(jù)切削三要素和刀具磨損量指定如下方案(見表1)：

在1∶5萬航磁ΔT異常等值線平面圖和剖面平面圖上，異常呈近EW向帶狀展布，長度約40km，上桃園以西，ΔT>400nT的區(qū)域NW向長度約18km，NS向?qū)挾燃s2.5km，異常特征為南部磁場變化緩慢、梯度小且無明顯的局部負(fù)異常，而其北部磁場變化較快、梯度較大且有較強(qiáng)的局部負(fù)異常相伴生；上桃園向東，磁場呈現(xiàn)出劇烈跳躍變化的特征，在向城鎮(zhèn)南側(cè)由于鐵礦體埋深較大，磁場變化也變得較平穩(wěn)；而其南支場值變化較復(fù)雜，異常形態(tài)為不規(guī)則狀。經(jīng)對異常向上延拓處理。異常極大值的軸向幾乎沒有改變，只是其水平位置略向南移動(dòng)，由此可知引起磁異常的異常源應(yīng)為向下延深較大，略向南傾的磁性體。

表1　因素水平表(三因素四水平)

以表1為依據(jù)，分別對刀具初期磨損、正常磨損和嚴(yán)重磨損狀態(tài)下的振動(dòng)信號進(jìn)行采集，共采集48組信號，圖4為主軸轉(zhuǎn)速V=800 r/min，進(jìn)給速度f=12.5 mm/min，切削深度ap=1 mm下嚴(yán)重磨損狀態(tài)下的原始信號和經(jīng)過小波去噪后的信號：

圖4　原始信號與去噪后的信號 Fig.4 Original signal and denoised signal

圖5是刀具在嚴(yán)重磨損狀態(tài)下的EMD分解，其切削參數(shù)為：主軸轉(zhuǎn)速V=800 r/min，進(jìn)給速度f=12.5 mm/min，切削深度ap=1 mm。圖6為原始信號的重構(gòu)，因?yàn)?EMD分解在原理上是完備的，也就是說它每次都在原始信號中減去一個(gè)分量，最后剩余一個(gè)趨勢分量，因此將各個(gè)IMF以及余量求和便可得到完全一樣的原始信號。除此之外，可以在重構(gòu)過程中發(fā)現(xiàn)，信號的幾個(gè)高頻分量疊加即可反映原始信號的一系列特征，而不用將信號所有IMF 進(jìn)行疊加，即只要將在原始信號中確實(shí)存在的高頻分量成分疊加即可。從圖6中就可以看出，當(dāng)疊加到高頻分量IMF9時(shí)，信號基本上就與原始信號相同了。

圖5　信號的EMD分解 Fig.5 Signal’s EMD

圖6　信號的重構(gòu) Fig.6 Signal reconstruction

圖7　IMF平均振幅 Fig.7 IMF average amplitude

如圖7為不同磨損狀態(tài)下的各個(gè)IMF分量的振幅均值?？梢园l(fā)現(xiàn)初期磨損、正常磨損與嚴(yán)重磨損狀態(tài)下的IMF1、IMF2、IMF7、IMF8、IMF9、IMF10 、IMF11、IMF12振幅差值都不是很大，而IMF3、IMF4、IMF5、IMF6則有明顯的差值，且大于1.5，因此IMF3、IMF4、IMF5、IMF6為刀具與工件接觸時(shí)所對應(yīng)的IMF分量，同時(shí)隨著刀具磨損的增大,刀具切削工件時(shí)對應(yīng)的IMF分量的振幅均值也在明顯的增大。因此，可以將IMF3、IMF4、IMF5、IMF6振幅均值作為刀具磨損的特征進(jìn)行記錄和提取。盡管這些結(jié)果不是線性的，但足夠訓(xùn)練任何的神經(jīng)網(wǎng)絡(luò)。

經(jīng)過EMD分解后，對每個(gè)IMF分量進(jìn)行希爾伯特變換，可以得到希爾伯特三維時(shí)頻幅值譜，為了便于觀察和分析，從不同狀態(tài)下的信號中隨機(jī)截取1 200個(gè)采樣點(diǎn)并放大信號，如圖8、圖9和圖10，其切削參數(shù)為：主軸轉(zhuǎn)速V=400 r/min，進(jìn)給速度f=10 mm/min，切削深度ap=0.6 mm，比較三種磨損狀態(tài)，可以發(fā)現(xiàn)在初期磨損和正常磨損狀態(tài)下幅值在不同頻率下變化不是很明顯，曲線波動(dòng)變化不是很大，而嚴(yán)重磨損則正好相反，并且在圖中低頻處出現(xiàn)明顯的散亂點(diǎn)，這都可以作為刀具磨損特征提取的參考。

圖8 初期磨損狀態(tài)下三維時(shí)頻幅值譜Fig.8Hilbertspectrumforinitialworn圖9 正常磨損狀態(tài)下三維時(shí)頻幅值譜Fig.9Hilbertspectrumfornormalworn圖10 嚴(yán)重磨損狀態(tài)下三維時(shí)頻幅值譜Fig.10Hilbertspectrumforseverelyworn

表2　磨損特征

圖11、圖12和圖13分別為刀具初期磨損、正常磨損和嚴(yán)重磨損狀態(tài)下的IMF4分量的希爾伯特邊際譜，其切削參數(shù)為：主軸轉(zhuǎn)速V=800 r/min，進(jìn)給速度f=12.5 mm/min，切削深度ap=1 mm，從圖中可以看出初期磨損和正常磨損狀態(tài)下它的邊際譜最大幅值點(diǎn)為0.002 764 2、0.007 451 3，而嚴(yán)重磨損狀態(tài)下邊際譜最大幅值點(diǎn)為0.042 157，即隨著刀具磨損的加劇，它的邊際譜最大幅值點(diǎn)也在隨著增大，且增幅明顯。

經(jīng)過試驗(yàn)得到的不同切削狀態(tài)和不同磨損狀態(tài)下的48組數(shù)據(jù)，從中隨機(jī)選取16組數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的訓(xùn)練樣本，隨機(jī)選取8組數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的目標(biāo)樣本，經(jīng)過希爾伯特黃變換處理提取特征，如表2，將前16組作為訓(xùn)練樣本，后8組作為目標(biāo)樣本：

將上述IMF和邊際譜最大幅值作為神經(jīng)網(wǎng)絡(luò)的輸入特征向量輸入神經(jīng)網(wǎng)絡(luò)，由于BP神經(jīng)網(wǎng)絡(luò)每次訓(xùn)練都有隨機(jī)數(shù)的產(chǎn)生，故將輸入特征向量進(jìn)行5次訓(xùn)練和5次神經(jīng)網(wǎng)絡(luò)判斷，取平均值，結(jié)果(見表3)：

表3　實(shí)驗(yàn)結(jié)果

圖11 初期磨損狀態(tài)下IMF4希爾伯特邊際譜Fig.11MarginalspectrumofIMF4forinitialworn圖12 正常磨損狀態(tài)下IMF4希爾伯特邊際譜Fig.12MarginalspectrumofIMF4fornormalworn圖13 嚴(yán)重磨損狀態(tài)下IMF4希爾伯特邊際譜Fig.13MarginalspectrumofIMF4forseverelyworn

4結(jié)論

針對目前刀具磨損監(jiān)測中信號特征提取方法的不足，提出了一種在小波去噪的基礎(chǔ)上運(yùn)用希爾伯特黃變換來進(jìn)行刀具磨損特征提取的方法，得出以下結(jié)論：

(1) 該方法可以將不同頻率下的信號逐步分離，方便找出與刀具磨損相關(guān)的敏感信號，有助于有目的性的分析。

(2) 彌補(bǔ)了傳統(tǒng)信號分析方法的高分辨率和自適應(yīng)不足，有效的反映和保留了刀具磨損時(shí)的非線性非平穩(wěn)信號的特征。

(3) 所提取的刀具磨損特征IMF3、IMF4、IMF5、IMF6的振幅均值和邊際譜中的最大幅值能準(zhǔn)確可靠地反應(yīng)刀具磨損狀態(tài)。

結(jié)果證明，基于希爾伯特黃變換的刀具磨損特征提取可以作為刀具磨損監(jiān)測的一種簡單和可靠地方法。但由于加工條件的復(fù)雜性和加工過程中的不穩(wěn)定性，周圍環(huán)境的干擾噪聲嚴(yán)重的影響了信號的分析精度，因此也需要通過合理的方式消除或減少噪聲的干擾，僅靠小波去噪的方式是遠(yuǎn)不夠的，除此之外，目前的智能決策方法如神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)等都有它本身的局限性和對樣本的依賴性，對于造價(jià)昂貴且特質(zhì)的刀具，建立刀具磨損樣本數(shù)據(jù)庫是不現(xiàn)實(shí)的，因此，進(jìn)一步在HHT的基礎(chǔ)上研究基于小樣本的刀具磨損智能決策方法也是研究的一大重點(diǎn)。

參 考 文 獻(xiàn)

[1] Zhu K P,Wong Y S, Hong G S. Wavelet analysis of sensor signals for tool condition monitoring:A review and some new results [J]. International Journal of Machine Tools & Manufacture 2009,49:537-553.

[2] 吳琛,周瑞忠. 基于 Hilbert 譜的結(jié)構(gòu)動(dòng)力響應(yīng)非線性特征分析[J]. 振動(dòng)與沖擊,2013, 32(14):71-76.

WU Chen ，ZHOU Rui-zhong. Nonlinear behavior analysis of structural dynamic response based on Hilbert spectrum[J]. Journal of Vibration and Shock,2013,32(14):71-76.

[3] 陳雙喜，林建輝，陳建政. 基于希爾伯特-黃變換提取車輛軌道耦合系統(tǒng)時(shí)頻特性[J].振動(dòng)與沖擊,2013,32(5):43-47.

CHEN Shuang-xi,LIN Jian-hui, CHEN Jian-zheng.Time-frequency characteristics extraction of vehicle-track coupling system based on Hilbert-Huang transform[J]. Journal of Vibration and Shock, 2013,32(5):43-47.

[4] Raja J E,Kiong L C,Soong L W. Hilbert-Huang transform-based emitted sound signal analysis for tool flank wear monitoring[J]. Arabian Journal for Science and Engineering，2013，38:2219-2226,

[5] Huang N E，Shen Zheng，Long S R，et al．The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]．Proc. R. Soc. Lond. 1998，A：903-995．

[6] Bisu C F, Zapciu M,Cahuc O,et al. Envelope dynamic analysis: a new approach for milling process monitoring[J]. The International Journal of Advanced Manufacturing Technology,2012,62(9),471-486,

[7] Du Wei, Liu Quan. A new ball bearing fault diagnosis method based on EMD and SVM[J]. Advanced Electrical and Electronics Engineering,2011, 2: 423-428.

[8] Gupta R,Kumar A,Bahl R. Estimation of instantaneous frequencies using iterative empirical mode decomposition[J]. Signal, Image and Video Processing,2012, 10(5):1863-1703.

[9] 符嬈，張群巖，趙述元. 航空發(fā)動(dòng)機(jī)試飛中轉(zhuǎn)靜子碰摩故障信號處理的希爾伯特-黃變換(HHT)方法[J].噪聲與振動(dòng)控制,2012,2:123-127.

FU Rao, ZHANG Qun-yan,ZHAO Shu-yuan. Signal processing of rubbing trouble between rotor and stator of aeroengines in flight-test using Hilbert-Huang transform [J]. Noise and Vibration Control, 2012,2:123-127.

[10] 鐘佑明,秦樹人,湯寶平. 希爾伯特-黃變換中邊際譜的研究[J].系統(tǒng)工程與電子技術(shù),2004,26(9):1324-1326.

ZHONG You-ming, QIN Shu-ren, TANG Bao-ping. Study on the marginal spectrum in Hilbert-Huang transform[J]. Systems Engineering and Electronics,2004.26(9):1324-1326.

[11] Bassiuny A M,Li X. Flute breakage detection during end milling using Hilbert-Huang transformand smoothed nonlinear energy operator[J].International Journal of Machine Tools and Manufacture 2007,47:1011-1020.

[12] 代海波，單銳，王換鵬，等. 基于改進(jìn)閾值函數(shù)的小波去噪算法研究[J]. 噪聲與振動(dòng)控制,2012,44(6):189-193.

DAI Hai-bo,SHAN Rui,WANG Huan-peng,et al. Study on wavelet denoising algorithm based on iimproved threshold function and maximal overlap discrete wavelet packet transform [J]. Noise and Vibration Control, 2012,44(6):189-193.

[13] 唐進(jìn)元，陳維濤，陳思雨，等. 一種新的小波閾值函數(shù)及其在振動(dòng)信號去噪分析中的應(yīng)用[J]. 振動(dòng)與沖擊，2009,28(7)：118-121.

TANG Jin-yuan, CHEN Wei-tao, CHEN Si-yu, et al.Wavelet-based vibration signal denoising with a new adaptive thresholding function[J]. Journal of Vibration and Shock, 2009,28(7)：118-121.

[14] 何正嘉，陳進(jìn)，王太勇，等.機(jī)械故障診斷理論與應(yīng)用[M].北京：高等教育出版社，2010.

[15] 李志西，杜雙奎. 試驗(yàn)優(yōu)化設(shè)計(jì)與統(tǒng)計(jì)分析[M].北京：科學(xué)出版社，2010.