孔系鉆削振動(dòng)信號(hào)特征波動(dòng)可視化研究及其應(yīng)用

周友行　楊文佳　謝賽元　張　俏　章

本毅湘潭大學(xué)，湘潭,411105

?

孔系鉆削振動(dòng)信號(hào)特征波動(dòng)可視化研究及其應(yīng)用

周友行楊文佳謝賽元張俏章

本毅湘潭大學(xué)，湘潭,411105

摘要：為解決工程應(yīng)用中切削參數(shù)一致的孔系加工質(zhì)量一致性評(píng)估的難題，提出了一種基于振動(dòng)信號(hào)特征波動(dòng)可視化的聚類分析方法。首先采用振動(dòng)傳感器監(jiān)控孔系鉆削過程，提取各孔振動(dòng)信號(hào)小波包能量譜和高階統(tǒng)計(jì)量特征；然后利用雷達(dá)圖得到各孔振動(dòng)信號(hào)特征矩陣分布圖，提取信號(hào)特征雷達(dá)圖多邊形重心特征；最后采用模糊C-均值(FCM)算法對(duì)雷達(dá)圖平面重心點(diǎn)集進(jìn)行聚類分析。理論分析結(jié)果與人工檢測結(jié)果對(duì)比表明：該方法可直觀呈現(xiàn)孔系鉆削質(zhì)量分布情況，簡便、可靠地實(shí)現(xiàn)孔系鉆削質(zhì)量的一致性評(píng)估。

關(guān)鍵詞：鉆削質(zhì)量；振動(dòng)信號(hào)；小波包能量譜；高階統(tǒng)計(jì)量；雷達(dá)圖；聚類分析

0引言

航空航天產(chǎn)品制造過程中，存在大量切削參數(shù)一致的孔系零件加工，但往往由于零部件材料成分不均勻和切削加工過程中的不可控因素(如機(jī)床振動(dòng)、刀具磨損破裂等)影響，導(dǎo)致孔系零件加工質(zhì)量差異較大，帶來嚴(yán)重的安全隱患。同時(shí)，孔系鉆削加工常處于封閉或半封閉的環(huán)境，很難直接檢測每一個(gè)鉆孔尺寸偏差、幾何精度、表面粗糙度和微觀物化性能等質(zhì)量表征參數(shù)。因此，研究可快速有效監(jiān)測和分析孔系鉆削加工質(zhì)量一致性評(píng)估的方法十分必要。

近年來，國內(nèi)外學(xué)者通過研究和分析鉆削過程監(jiān)測信號(hào)，在有效預(yù)測和控制單一鉆削質(zhì)量方面取得了一定成就。其中許多學(xué)者對(duì)鉆削過程中的振動(dòng)監(jiān)測信號(hào)進(jìn)行了廣泛研究：Ramirez等[1]基于振動(dòng)信號(hào)分析了鉆孔質(zhì)量和刀具磨損；Upadhyay等[2]基于振動(dòng)信號(hào)分析了加工表面粗糙度；Abu-Mahfouz[3]利用振動(dòng)信號(hào)對(duì)鉆頭磨損進(jìn)行了評(píng)估與分類；Messaoud等[4]研究了深孔鉆削顫動(dòng)振動(dòng)信號(hào)對(duì)鉆孔表面質(zhì)量的影響；文獻(xiàn)[5-7]融合振動(dòng)、聲發(fā)射、切削力、切削溫度多種傳感器信號(hào)，分析了其與鉆削質(zhì)量間的關(guān)系。這些研究表明：鉆削過程中的加速度振動(dòng)信號(hào)時(shí)頻域特征與鉆孔質(zhì)量之間存在耦合關(guān)系。

孔系鉆削過程中監(jiān)測信號(hào)時(shí)頻域特征的方法多種多樣，采用何種方法能有效并簡單地分析信號(hào)特征變化與孔系質(zhì)量波動(dòng)之間的關(guān)系，就成為了基于信號(hào)時(shí)頻域特征孔系質(zhì)量一致性評(píng)估的關(guān)鍵步驟。另外，由于孔系鉆削過程中獲得的信號(hào)特征本質(zhì)上是一堆龐大的多元數(shù)據(jù)，故可應(yīng)用數(shù)據(jù)可視化技術(shù)[8]，將信號(hào)多元特征數(shù)據(jù)置于視覺空間中，直觀了解信號(hào)時(shí)頻特征變化與孔系加工質(zhì)量波動(dòng)間的關(guān)聯(lián)。李惠君等[9]證明，雷達(dá)圖能直觀地表示多元數(shù)據(jù)特征間的變化關(guān)系。

本文以孔系鉆削過程振動(dòng)監(jiān)測信號(hào)為研究對(duì)象，提取監(jiān)測信號(hào)的小波包能量譜和偏度-峰度特征，然后采用基于雷達(dá)圖和模糊C-均值(FCM)算法的可視化聚類分析法，獲得孔系鉆削振動(dòng)信號(hào)波動(dòng)特征分類結(jié)果，最后通過與人工檢測結(jié)果對(duì)比分析，實(shí)現(xiàn)孔系鉆削加工質(zhì)量的一致性評(píng)估。

1孔系鉆削過程振動(dòng)信號(hào)描述

本文實(shí)驗(yàn)所采用設(shè)備如下：JOINT-5KI數(shù)控銑床、三向振動(dòng)加速度計(jì)Kistler8793A、Kennametal鉆頭004D-0800C(鉆頭直徑8mm)。工件材料為鈦合金。具體工藝參數(shù)如下：主軸轉(zhuǎn)速n=750 r/min，進(jìn)給速度f=30 mm/min，通孔鉆孔深度h=15 mm；采樣頻率20 kHz；干切削。實(shí)驗(yàn)中，將三軸加速度振動(dòng)傳感器固定在主軸箱壁上，如圖1所示。

圖1　傳感器安裝位置

在眾多鉆孔監(jiān)測數(shù)據(jù)中隨機(jī)選取24組鉆孔實(shí)驗(yàn)數(shù)據(jù)，其中由于x軸向和y軸向加速度振動(dòng)信號(hào)源自工件橫向的微小振動(dòng)，兩者信號(hào)變化基本一致，而z軸向加速度振動(dòng)信號(hào)源于工件縱向振動(dòng)，故選擇x軸向和z軸向加速度振動(dòng)信號(hào)作為本文的研究對(duì)象，其時(shí)域波形如圖2、圖3所示。

圖2　x軸向振動(dòng)信號(hào)

圖3　z軸向振動(dòng)信號(hào)

對(duì)24個(gè)鉆削加工質(zhì)量如孔的尺寸、幾何精度、表面粗糙度等依次進(jìn)行人工檢測，具體結(jié)果如表1所示。其中，A表示加工質(zhì)量較好，如鉆孔表面平整光潔，Ra≤3.2 μm；B表示質(zhì)量合格，如鉆孔表面較平整光潔，3.26.3 μm。表1中，第23號(hào)鉆孔內(nèi)表面嚴(yán)重刮擦，第24號(hào)鉆孔鉆頭崩刃，鉆削過程中斷，產(chǎn)生內(nèi)表面嚴(yán)重擦刮的錐形孔。

表1　孔系鉆削質(zhì)量人工檢測結(jié)果

圖2、圖3所示的孔系鉆削過程振動(dòng)信號(hào)可用矩陣S描述如下：

(1)

Xa、Za表示第a(a=1,2,…,24)號(hào)鉆孔x軸向、z軸向振動(dòng)信號(hào)向量：

(2)

式中Xab、Zab為第a號(hào)鉆孔第b(b=1,2,…,c)個(gè)x軸向、z軸向振動(dòng)信號(hào)采樣點(diǎn)數(shù)值；c為第a號(hào)鉆孔振動(dòng)信號(hào)采樣總點(diǎn)數(shù)。

2孔系鉆削振動(dòng)信號(hào)波動(dòng)分析

由于鉆削過程復(fù)雜，振動(dòng)監(jiān)測信號(hào)波動(dòng)影響因素多，故很難直接定量、全面地得出信號(hào)與單個(gè)鉆孔加工質(zhì)量表征參數(shù)間的內(nèi)在聯(lián)系。對(duì)于孔系振動(dòng)監(jiān)測信號(hào)，若在時(shí)頻域?qū)用嫔咸崛⌒盘?hào)特征，通過統(tǒng)計(jì)分析和數(shù)據(jù)挖掘，基于各鉆孔振動(dòng)信號(hào)特征變化與加工質(zhì)量表征參數(shù)波動(dòng)間的耦合關(guān)系，即可開展孔系加工質(zhì)量一致性評(píng)估。

孔系鉆削加工異常引起的信號(hào)突變一定會(huì)反映在信號(hào)中某些頻帶能量的變化上。采用小波包分解技術(shù)將振動(dòng)信號(hào)在全局范圍內(nèi)正交地分解在獨(dú)立的相應(yīng)頻帶內(nèi)，提取頻帶能量特征，分析各孔振動(dòng)信號(hào)能量在不同頻帶的變化規(guī)律，從而間接反映孔系鉆削質(zhì)量的波動(dòng)。因此，信號(hào)小波包能量譜特征變化可作為評(píng)估孔系鉆削質(zhì)量波動(dòng)的一個(gè)指標(biāo)。

此外，孔系鉆削過程中，若工件材料均勻，鉆削狀態(tài)平穩(wěn)，各孔材料去除過程中信號(hào)表現(xiàn)形式將會(huì)基本一致，如圖2、圖3中的1～22組信號(hào)，信號(hào)總體上呈現(xiàn)較強(qiáng)的規(guī)則性和平穩(wěn)性。此時(shí)，若視各孔x軸向加速度信號(hào)為一組非平穩(wěn)、非高斯隨機(jī)過程，則各孔x軸向加速度信號(hào)偏離高斯分布的程度將會(huì)一致；若其中一個(gè)孔鉆削過程異常，則該孔x軸向加速度振動(dòng)信號(hào)會(huì)產(chǎn)生突變，如圖2、圖3中的23組、24組信號(hào)，信號(hào)偏離高斯分布狀況會(huì)明顯與1～22組信號(hào)不一致，所以，描述信號(hào)偏離高斯分布的特征也可以作為評(píng)估孔系鉆削質(zhì)量波動(dòng)的另一個(gè)指標(biāo)。

因此，通過分析各孔鉆削監(jiān)測信號(hào)的小波包能量譜變化和信號(hào)的偏離高斯分布狀況，可比較全面地實(shí)現(xiàn)孔系鉆削質(zhì)量一致性檢測控制。

2.1孔系鉆削振動(dòng)信號(hào)小波包能量特征提取

信號(hào)小波包分解遵循能量守恒原理。小波包頻帶能量譜(wavelet packet energy spectrum，WPES)就是按照能量形式表示小波包分解結(jié)果。以第a號(hào)鉆削過程振動(dòng)信號(hào)序列Za為例，對(duì)其進(jìn)行j層小波包分解，分為2j個(gè)頻帶，則其小波包頻帶能量譜分析原理[10-11]如下：

(3)

k=1,2,…，j

式中，k為小波包分解層數(shù)；m為分解頻帶的位置序號(hào)；E()為數(shù)學(xué)期望。

(4)

則第a號(hào)鉆孔小波包分解的第k層m頻帶的相對(duì)能量為

(5)

第a號(hào)鉆孔小波包分解能量譜(即以各頻帶能量比為元素的特征向量)表示如下：

(6)

據(jù)此，本文實(shí)驗(yàn)采集的z軸振動(dòng)信號(hào)的小波包分解能量譜可以表示為

(7)

采用此方法，可以得到采集信號(hào)在不同頻帶內(nèi)能量的分布特征，同時(shí)獲得不同頻帶信號(hào)能量占總能量的百分比。同理，可獲得x軸振動(dòng)信號(hào)的小波包能量譜。

本文采用db5小波進(jìn)行4層小波包分解，將其分解在16個(gè)頻帶內(nèi)。據(jù)式(7)所示的小波包頻帶能量譜特征矩陣，可繪制x軸和z軸振動(dòng)信號(hào)的能量譜柱狀圖，見圖4。

(a)x軸振動(dòng)信號(hào)頻帶能量圖

(b)z軸振動(dòng)信號(hào)頻帶能量圖圖4　小波包分解頻帶能量柱狀圖

從圖4中可看出，各孔鉆削振動(dòng)信號(hào)小波包分解的各頻帶能量都有不同程度的變化，其中1～22號(hào)鉆孔x軸振動(dòng)信號(hào)的小波包頻帶能量主要集中在第1頻帶，23號(hào)和24號(hào)鉆孔x軸振動(dòng)信號(hào)在低頻部分能量大幅度下降，能量分散轉(zhuǎn)移到第6、第8頻帶。而14號(hào)、15號(hào)以及19～24號(hào)鉆孔z軸振動(dòng)信號(hào)的能量分散在第1、第5、第13頻帶。

2.2孔系鉆削振動(dòng)信號(hào)高階統(tǒng)計(jì)量特征提取

在信號(hào)處理分析中，常采用概率密度函數(shù)描述信號(hào)偏離高斯分布的狀況，而對(duì)于非高斯的鉆削加速度振動(dòng)監(jiān)測信號(hào)，其概率密度函數(shù)的獲取是非常困難的，一般采用信號(hào)三階和四階統(tǒng)計(jì)量來描述信號(hào)的概率密度特征。三階統(tǒng)計(jì)量稱為信號(hào)偏度，偏度體現(xiàn)信號(hào)概率分布曲線的非對(duì)稱性；四階統(tǒng)計(jì)量稱為信號(hào)峰度，峰度值可用來描述概率分布曲線較之高斯信號(hào)表現(xiàn)的尖削或平坦的程度，峰度指標(biāo)對(duì)振動(dòng)時(shí)域信號(hào)的瞬時(shí)脈沖非常敏感。以第a號(hào)鉆孔軸向振動(dòng)信號(hào)Za為例，計(jì)算Za的偏度SZα、峰度KZα的表達(dá)式如下：

SZa=E(Za-μ)3/σ3

(8)

KZa=E(Za-μ)4/σ4

(9)

式中，μ為樣本均值;σ為樣本方差。

由式(8)、式(9)分別計(jì)算孔系鉆削過程中的x軸、z軸振動(dòng)信號(hào)的偏度值、峰度值，如圖5所示。

(a)x軸振動(dòng)信號(hào)的偏度

(b)x軸振動(dòng)信號(hào)峰度

(c)z軸振動(dòng)信號(hào)偏度

(d)z軸振動(dòng)信號(hào)峰度圖5　鉆削振動(dòng)信號(hào)偏度值、峰度值

從圖5a、圖5b中可以看出，23號(hào)、24號(hào)鉆孔x軸振動(dòng)信號(hào)明顯異于其他鉆孔的偏度、峰度值；孔系鉆削z軸振動(dòng)信號(hào)的偏度、峰度值波動(dòng)較大，如圖5c、圖5d所示，14號(hào)、15號(hào)、19～24號(hào)鉆孔的偏度絕對(duì)值以及1號(hào)、14號(hào)、15號(hào)、20～24號(hào)鉆孔的峰度值較其他鉆孔的值偏大。

3基于雷達(dá)圖的鉆削振動(dòng)信號(hào)特征數(shù)據(jù)可視化

從圖4、圖5中可以看出，23號(hào)、24號(hào)鉆孔振動(dòng)監(jiān)測信號(hào)小波包能量譜和偏度、峰度特征變化都異于其他鉆孔，可初步判斷23號(hào)、24號(hào)鉆孔加工過程中出現(xiàn)了異常；而14號(hào)、15號(hào)、20號(hào)、21號(hào)、22號(hào)鉆孔z軸振動(dòng)信號(hào)小波包能量和偏度、峰度卻表現(xiàn)出不一致，如圖4b、圖5c、圖5d所示。為進(jìn)一步討論研究孔系鉆削質(zhì)量的一致性，可視孔系鉆削過程各信號(hào)的上述時(shí)頻特征為一個(gè)數(shù)據(jù)集合樣本，采用聚類算法綜合考慮所有信號(hào)特征，對(duì)特征數(shù)據(jù)進(jìn)行分類識(shí)別。

對(duì)于孔系鉆削監(jiān)測信號(hào)頻帶能量特征，從概率統(tǒng)計(jì)的角度看，可僅選擇頻帶能量變化最為敏感的若干頻帶作為特征向量。從圖4中可看出，x軸振動(dòng)信號(hào)第1、第6、第8頻帶，z軸振動(dòng)信號(hào)第1、第5、第13頻帶信號(hào)能量變化比較明顯， 而且這些特征頻段的能量占據(jù)了總能量的60%左右，以這些頻帶構(gòu)建特征向量，可覆蓋原始信號(hào)的絕大部分信息。在此基礎(chǔ)上，再考慮加上x軸和z軸振動(dòng)信號(hào)偏度、峰度特征共10項(xiàng)特征，可構(gòu)造第a號(hào)鉆孔振動(dòng)信號(hào)的綜合特征向量Ta：

SZa,KXa,KZa)

(10)

將Ta進(jìn)行歸一化處理后，由于綜合特征空間的維數(shù)較多，其空間點(diǎn)集規(guī)律的表達(dá)較困難，本文引入雷達(dá)圖聚類可視化方法，不僅可綜合多維特征分析孔系鉆削質(zhì)量，同時(shí)可直觀地表示多維綜合特征矩陣數(shù)據(jù)。

3.1孔系鉆削振動(dòng)信號(hào)特征矩陣?yán)走_(dá)圖表示

雷達(dá)圖又稱星點(diǎn)圖或蜘蛛圖，是一種將多維空間中的點(diǎn)映射到二維空間的一種可視化方法，即將多維數(shù)據(jù)的各個(gè)屬性值用二維平面圖形表示。本文孔系鉆削振動(dòng)信號(hào)多特征數(shù)據(jù)可視化的雷達(dá)圖繪制方法如下[12]：

(1)作一個(gè)圓，并把圓周分為10等分；

(2)連接圓心和各分點(diǎn)，將這10條半徑依次定義為各指標(biāo)的坐標(biāo)軸，并標(biāo)以適當(dāng)?shù)目潭龋?/p>

(3)對(duì)給定的每一次鉆削過程信號(hào)綜合特征向量Ta，將10個(gè)指標(biāo)值分別取在相應(yīng)的坐標(biāo)軸上，然后將它們連接成一個(gè)10邊形；

(4) 24組數(shù)據(jù)即可形成24個(gè)10邊形，如圖6所示。

3.2雷達(dá)圖重心點(diǎn)集聚類分析

(a)1～4號(hào)鉆孔　　　　 (b)5～8號(hào)鉆孔

(c)9～12號(hào)鉆孔　　　　(d)13～16號(hào)鉆孔

(e)17～20號(hào)鉆孔　　　(f)21～24號(hào)鉆孔 圖6　孔系鉆削振動(dòng)監(jiān)測信號(hào)時(shí)頻特征雷達(dá)圖

雷達(dá)圖表示法直觀地表示了信號(hào)時(shí)頻特征的數(shù)據(jù)結(jié)構(gòu)，圖6中，不同鉆孔信號(hào)時(shí)頻特征矩陣數(shù)據(jù)對(duì)應(yīng)不同的雷達(dá)圖圖形特征。此處所指的圖形特征是指能夠較準(zhǔn)確地描述一個(gè)圖形的形狀，并且可以用來區(qū)別另一個(gè)類別的圖形屬性。顯然，由圖6中雷達(dá)圖的對(duì)比發(fā)現(xiàn)，當(dāng)質(zhì)量出現(xiàn)異常時(shí)，監(jiān)測信號(hào)特征的雷達(dá)圖多邊形重心點(diǎn)位置分布存在著很大的差異。因此，本文選取雷達(dá)圖的重心坐標(biāo)作為聚類分析的特征量。

將上述24個(gè)10邊形雷達(dá)圖的重心坐標(biāo)繪制于一個(gè)極坐標(biāo)系中，利用雷達(dá)圖將高維數(shù)據(jù)的聚類轉(zhuǎn)化成平面上點(diǎn)集的聚類。本文采用FCM算法對(duì)重心點(diǎn)進(jìn)行聚類，每個(gè)聚類簇的數(shù)據(jù)用不同標(biāo)注符號(hào)和顏色表示，并在旁邊標(biāo)記鉆削序列。聚類結(jié)果如圖7所示，其中23號(hào)、24號(hào)鉆孔的雷達(dá)圖重心點(diǎn)明顯偏離了第1～22號(hào)雷達(dá)圖重心點(diǎn)的聚類中心。

圖7　24組數(shù)據(jù)雷達(dá)圖重心點(diǎn)集聚類分析結(jié)果

由表1、表2可知，人工檢測結(jié)果與圖7中雷達(dá)圖重心點(diǎn)集聚類結(jié)果基本一致。在此次實(shí)驗(yàn)中，聚類分析各類別準(zhǔn)確率分別為100%、83.3%和100%，總的準(zhǔn)確率為94.4%。

表2　FCM聚類結(jié)果

4實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證此方法的可行性，改變鉆削參數(shù)和實(shí)驗(yàn)材料再次進(jìn)行了多組鉆削實(shí)驗(yàn)。實(shí)驗(yàn)中具體參數(shù)如下：通孔鉆削深度h=14.44 mm，刀具為1534SU03-0600；工件材料為45鋼；進(jìn)給量f=30 mm/min；主軸轉(zhuǎn)速n=500 r/min；采樣頻率為20 kHz；干切削。

采用本文所述方法對(duì)隨機(jī)抽取的38組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析，并將分類結(jié)果(圖8)與人工檢測結(jié)果進(jìn)行對(duì)比，對(duì)比結(jié)果如表3所示。本次驗(yàn)證實(shí)驗(yàn)聚類分析各類別準(zhǔn)確率分別為100%、80%和100%，總的準(zhǔn)確率為93.3%。

圖8　38組數(shù)據(jù)雷達(dá)圖重心點(diǎn)集聚類分析結(jié)果

表3　聚類分析結(jié)果與人工檢測結(jié)果

因此，由兩次實(shí)驗(yàn)結(jié)果得知，若采用文中所述方法對(duì)鉆削過程振動(dòng)信號(hào)進(jìn)行分類后，只需從不同類別中各找一個(gè)鉆孔做質(zhì)量檢測，就能有效判斷孔系鉆削質(zhì)量。

5結(jié)論

(1)孔系鉆削過程加速度振動(dòng)監(jiān)測信號(hào)的小波包能量譜及高階統(tǒng)計(jì)量特征變化與孔系加工質(zhì)量波動(dòng)間存在聯(lián)系。

(2)將監(jiān)測信號(hào)的小波包能量譜和高階統(tǒng)計(jì)量特征組成綜合特征矩陣，利用雷達(dá)圖實(shí)現(xiàn)了多維綜合特征數(shù)據(jù)的可視化，不同鉆孔監(jiān)測信號(hào)的時(shí)頻特征矩陣數(shù)據(jù)可對(duì)應(yīng)不同的雷達(dá)圖重心坐標(biāo)特征。

(3)采用FCM算法聚類方法對(duì)信號(hào)時(shí)頻域特征進(jìn)行聚類融合，可克服傳感器單一特征可能會(huì)對(duì)孔系質(zhì)量分布產(chǎn)生誤判的不足，提高分類準(zhǔn)確率。

(4)采用雷達(dá)圖可視化技術(shù)處理信號(hào)時(shí)頻特征，可直觀地觀察孔系加工質(zhì)量的分布狀況，加工質(zhì)量特別異常的鉆孔從重心點(diǎn)聚類圖中直觀顯示。

(5)采用基于雷達(dá)圖的可視化聚類方法進(jìn)行孔系零件鉆孔質(zhì)量一致性評(píng)估檢測，聚類結(jié)果顯示直觀，為人工檢測提供參考，從而準(zhǔn)確、快速地實(shí)現(xiàn)孔系質(zhì)量分類。

(6)實(shí)驗(yàn)結(jié)果表明，本文所提出的方法適用于切削工藝參數(shù)、傳感器檢測環(huán)境一致的鉆削加工過程，其他切削過程是否同樣適用則尚需進(jìn)一步研究。

參考文獻(xiàn)：

[1]RamirezC,PoulachonG,RossiF,etal.ToolWearMonitoringandHoleSurfaceQualityduringCFRPDrilling[J].ProcediaCIRP,2014(13): 163-168.

[2]UpadhyayV,JainPK,MehtaNK.In-processPredictionofSurfaceRoughnessinTurningofTi-6Al-4VAlloyUsingCuttingParametersandVibrationSignals[J].Measurement,2013,46(1):154-160.

[3]Abu-MahfouzI.DrillingWearDetectionandClassificationUsingVibrationSignalsandArtificialNeuralNetwork[J].InternationalJournalofMachineToolsandManufacture, 2003,43(7):707-720.

[4]MessaoudA,WeihsC,HeringF.DetectionofChatterVibrationinaDrillingProcessUsingMultivariateControlCharts[J].ComputationalStatistics&DataAnalysis, 2008, 52(6): 3208-3219.

[5]AxinteDA,GindyN,FoxK,etal.ProcessMonitoringtoAssisttheWorkpieceSurfaceQualityinMachining[J].InternationalJournalofMachineToolsandManufacture, 2004,44(10):1091-1098.

[6]CostesJP.APredictiveSurfaceProfileModelforTurningBasedonSpectralAnalysis[J].JournalofMaterialsProcessingTechnology, 2013,213(1):94-100.

[7]RawatS,AttiaH.CharacterizationoftheDryHighSpeedDrillingProcessofWovenCompositesUsingMachinabilityMapsApproach[J].CIRPAnnals, 2009, 58(1):105-108.

[8]YauN.VisualizeThis:theFlowingDataGuidetoDesign,VisualizationandStatistics[M].LosAngeles:Wiley, 2011.

[9]李惠君,李志全.基于改進(jìn)雷達(dá)圖的可視化聚類方法研究[J]. 燕山大學(xué)學(xué)報(bào), 2013,37(1): 58-62.

LiHuijun,LiZhiquan.ResearchonVisualizationClusteringMethodBasedonImprovedRadarChart[J].JournalofYanshanUniversity,2013,37(1):58-62.

[10]葛繼平,李胡生.基于小波包能量累積變異的梁橋損傷識(shí)別試驗(yàn)[J]. 振動(dòng)、測試與診斷,2011,31(1):85-88，131.

GeJiping,LiHusheng.ExperimentonDamageDetectionofGirderBridgeUsingWaveletPacketEnergyAccumulatedVariation[J]. 2011,31(1):85-88，131.

[11]王冬云,張文志. 基于小波包變換的滾動(dòng)軸承故障診斷[J]. 中國機(jī)械工程,2012,23(3):295-298.

WangDongyun,ZhangWenzhi.FaultDiagnosisStudyofBallBearingBasedonWaveletPacketTransform[J].ChinaMechanicalEngineering2012,23(3): 295-298.

[12]李昕,楊麗娟.基于圖表達(dá)原理的多維數(shù)據(jù)特征的選取[J].國外電子測量技術(shù),2007(10):17-20.

LiXin,YangLijuan.Multi-dimensionalDataFeatureSelectionBasedonChartTheory[J].ForeignElectronicMeasurementTechnology, 2007(10):17-20.

(編輯袁興玲)

收稿日期：2015-08-18

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目(51375419，51375418)；湖南省自然科學(xué)基金資助項(xiàng)目(2016JJ2084)

中圖分類號(hào)：TH166

DOI：10.3969/j.issn.1004-132X.2016.12.010

作者簡介：周友行，男，1971年生。湘潭大學(xué)機(jī)械工程學(xué)院教授、博士研究生導(dǎo)師。主要研究方向?yàn)閿?shù)字化設(shè)計(jì)與制造。發(fā)表論文50余篇。楊文佳，男，1990年生。湘潭大學(xué)機(jī)械工程學(xué)院碩士研究生。謝賽元，女，1989年生。湘潭大學(xué)機(jī)械工程學(xué)院碩士研究生。張俏，女，1990年生。湘潭大學(xué)機(jī)械工程學(xué)院碩士研究生。章本毅，男，1991年生。湘潭大學(xué)機(jī)械工程學(xué)院碩士研究生。

VisualizationofHoleSeriesDrillingVibrationSignalFeatureFluctuationandItsApplications

ZhouYouhangYangWenjiaXieSaiyuanZhangQiaoZhangBenyi

XiangtanUniversity,Xiangtan,Hunan,411105

Abstract：To analyse the consistency evaluation of hole series part drilling quality, a clustering analysis method was presented based on fluctuation of vibration signal features. A vibration sensor was used to monitor the drilling process, and the vibration signal wavelet packet energy spectrum and higher-order statistics features of each hole in drilling process were extracted to construct the monitoring signal characteristic matrix. Then radar chart was used to reconstruct the distribution graphs of wavelet packet energy spectrum features and high-order statistics features, and the center of gravity features were extracted from these radar charts. Finally, fuzzy C-means(FCM) algorithm was used to complete the clustering analyses of these center of gravity features points. By comparison, it shows that hole drilling quality distribution can be visually presented and hole drilling quality consistency evaluation can be realized simply and reliably by clustering analysis based on radar chart.

Key words:drilling quality；vibration signal；wavelet packet energy spectrum；high-order statistics；radar chart；clustering analysis