基于EMD精密球研磨過程振動信號分析研究

趙文宏，謝良江，陳　澤，安　平

(浙江工業(yè)大學(xué) 特種裝備制造與先進(jìn)加工技術(shù)教育部重點(diǎn)實(shí)驗室,浙江 杭州 310014)

基于EMD精密球研磨過程振動信號分析研究

趙文宏，謝良江，陳澤，安平

(浙江工業(yè)大學(xué) 特種裝備制造與先進(jìn)加工技術(shù)教育部重點(diǎn)實(shí)驗室,浙江 杭州 310014)

摘要：精密球在研磨過程中產(chǎn)生的振動，會影響球體的表面質(zhì)量,為了了解精密球研磨過程中的振動，引入經(jīng)驗?zāi)Ｊ椒纸?Empirical mode decomposition,EMD)方法,對精密球研磨過程中產(chǎn)生的振動信號進(jìn)行分析研究.以D-800 Ballpoli球體研磨機(jī)為實(shí)驗對象，系統(tǒng)把采集到的復(fù)雜振動信號轉(zhuǎn)化成簡單電信號，對信號進(jìn)行EMD分析處理，實(shí)驗研究在不同的研磨壓力、研磨速度以及研磨液質(zhì)量分?jǐn)?shù)的情況下，發(fā)現(xiàn)研磨過程中的振幅發(fā)生變化，研磨壓力、研磨速度以及研磨液質(zhì)量分?jǐn)?shù)是研磨振動的影響因素.隨著研磨時間加長，研磨振動幅度越來越小.在相同實(shí)驗條件下，該方法較傅里葉波形分析方法更加直觀方便。

關(guān)鍵詞：精密球加工；振動信號；經(jīng)驗?zāi)Ｊ椒纸?；振動幅?/p>

Study on vibration analysis of ball polish based on

emperical mode decomposition

ZHAO Wenhong, XIE Liangjiang, CHEN Ze, AN Ping

(Key Laboratory of Special Purpose Equipment and Advanced Manufacturing Technology, Ministry

of Education, Zhejiang University of Technology, Hangzhou 310014, China)

Abstract:The vibration caused by precision balls lapping will have a great influence on ball accuracy and surface quality. It is important to introduce empirical mode decomposition (EMD) to analyze the vibration caused during lapping process. The vibration signal of D-800 ballpoli was analyzed by EMD through transforming to electronica signal. The result also shows vibration amplitude is changed in different lapping pressure, lapping speed and lapping fluid. Lapping pressure, lapping speed and lapping fluid are vibration influencing factors. With the increasing of lapping time, the vibration amplitudes will become smaller. EMD is more effective and intuition than Fast Fourier Transform Algorithm with the same experiment conditions。

Keywords:precision ball lapping; vibration signal; empirical mode decomposition; vibration amplitude

研磨球在磨球機(jī)中的運(yùn)動是由研磨球環(huán)繞軸中心的公轉(zhuǎn)和自轉(zhuǎn)產(chǎn)生，同時球與球之間的撞擊，磨盤對球的壓力，設(shè)備的缺陷等都會導(dǎo)致振動[1].振動是超精密球體加工中一個主要的不可避免的阻力.研磨振動是影響球體表面質(zhì)量的一個因數(shù)[2-3]，因此研究研磨過程中振動，對提高被研磨球的表面質(zhì)量有幫助.目前，現(xiàn)有的振動分析方法有很多，主要有FFT、ARMA、Cohen類分布、STFT、小波變換及自適應(yīng)時頻分析等.筆者將經(jīng)驗?zāi)Ｊ椒纸夥椒?Empirical mode decomposition,EMD)引入到超精密球體加工領(lǐng)域.EMD是一種研究非平穩(wěn)數(shù)據(jù)序列的方法，這種信號處理方法首先由美國宇航局的Huang.N.E等于1998年提出[4].EMD方法是一種典型的自適應(yīng)分析方法，該方法把信號本身作為基函數(shù)，比FFT等利用外界函數(shù)作為基函數(shù)的分析方法誤差小，精確度高.目前，自適應(yīng)時頻分析方法，使用領(lǐng)域正在逐漸變多[5-8]。

以D-800 Ballpoli雙自轉(zhuǎn)磨球機(jī)為實(shí)驗對象，利用EMD振動信號處理方法，實(shí)現(xiàn)了對振動信號的監(jiān)測和分析，該方法和FFT分析方法對比，驗證了該方法的優(yōu)點(diǎn).研究了磨球機(jī)振動與磨盤轉(zhuǎn)速、研磨液、加載壓力的關(guān)系以及振動對球體粗糙度的影響。

1磨球機(jī)振動信號采集原理及硬件組成

振動檢測采用電測法，電測法測振動量是通過慣性傳感器得到.把傳感器與振動物體相連接，振動物體振動時，慣性傳感器傳出一個電信號，通過電信號得到振動量的大小.圖1是D-800 Ballpoli結(jié)構(gòu)示意圖.振動檢測傳感器安裝在圖中11上加壓盤，靠近被研磨球體。

1—床身；2—動力部分；3—外軸套；4—立柱；5—橫梁；6—渦輪蝸桿加壓機(jī)構(gòu)；7—絲杠；8—上壓力傳感器；9—加壓彈簧；10—下壓力傳感器；11—上加壓盤；12—內(nèi)盤；13—外盤；14—軸承圖1　D-800 Ballpoli結(jié)構(gòu)示意圖Fig.1　D-800 Ballpoli structure diagram

振動信號被傳感器采集后，通過信號放大電路被單片機(jī)識別，筆者采用了兩級單電源交流信號放大電路.電路前一級是同相放大電路，后一級為反相放大電路.經(jīng)過放大電路處理的信號是一個模擬信號，通過AD轉(zhuǎn)換，信號變?yōu)榭勺R別數(shù)字信號并進(jìn)行分析.本系統(tǒng)選用壓電式加速度傳感器AD100T傳感器，實(shí)驗過程中上研磨盤不動，檢測裝置相對上研磨盤保持靜止，下研磨盤轉(zhuǎn)動，被檢測球體經(jīng)過檢測裝置的頻率不超過1 000 Hz，而加速度傳感器AD100T的檢測頻率最高達(dá)到15 000 Hz,可以在較短時間內(nèi)完成對研磨球的振動檢測，因此該傳感器可以保證振動監(jiān)測的可靠性.安裝位置如圖2所示.加速度傳送來的信號然后通過XK343信號調(diào)理器，阿爾泰USB2812數(shù)據(jù)采集卡.從信號檢測點(diǎn)檢測到信號后，通過信號采集電路，把振動信號轉(zhuǎn)化為簡單電信號，再傳到上位機(jī)上進(jìn)行EMD分析處理，結(jié)果以圖表形式顯示在顯示屏上[9-13]。

圖2　振動檢測頭檢測位置示意圖Fig.2　Vibration detection location diagram

2振動信號的EMD分析

經(jīng)驗?zāi)B(tài)分解的主要內(nèi)容：先根據(jù)信號自身的時間特征分解成多個固有模態(tài)函數(shù)(Intrinsic mode function，IMF)，再對得到的Hilbert變換進(jìn)行處理，得到信號的瞬時頻率和振幅，從而得到時頻譜.最后對每個譜進(jìn)行Hilbert邊際譜的計算,通過邊際譜看出概率意義上幅值在整個數(shù)據(jù)跨度上的積累幅值.在表1的實(shí)驗條件與加工參數(shù)條件下，進(jìn)行實(shí)驗。

表1　實(shí)驗條件與加工參數(shù)

圖3　原始信號及其IMF圖Fig.3　Original signal and IMF chart

檢測到的振動信號原始圖形，并對其進(jìn)行EMD分解，得到IMF如圖3所示，圖形橫坐標(biāo)為時間，傳感器的采樣頻率是1 000 Hz，縱坐標(biāo)是信號幅值，單位是mm.固有模式函數(shù)是為了獲取多分量信號真正的頻率時變特征，將多分量信號分解為單分量信號，方便獲取信號瞬時頻率.圖4(a)是原始信號的邊際譜圖，圖4(b，c)是研磨一段時間后的邊際譜圖。

測量研磨過程中的任意20顆球，得到前面16 h的球徑平均差走勢，如圖5所示。

圖4　實(shí)驗信號邊際譜三個時間階段圖Fig.4　Marginal spectrum of three stages

圖5　球徑平均差走勢圖Fig.5　Average deviation of sphere diameter

瞬時頻率表示信號交變快慢的物理量，任一瞬時頻率都有一定的能量，將所有時刻的某個頻率的能量值加起來就是信號中該頻率的總能量，也就是邊際線的高度.邊際線反應(yīng)了在某個頻率點(diǎn)振動幅度的大小，如圖4(a)所示邊際譜，表示研磨過程中振動最大的是球體在50,150,250,700 Hz的振動幅度最大，排除交流電信號的干擾因素[14]，剩下150，250，700 Hz三個主要振動頻率，在圖4中，以250 Hz振動頻率為例，圖4(a)圖的振動幅度0.1 mm，圖4(b)圖的振動幅度是0.08 mm，圖4(c)圖的振動幅度是0.03 mm.說明隨著研磨時間的增加，研磨振動在不斷減小。

分別對導(dǎo)致振動的其中三個振動影響因素：研磨壓力、研磨速度和研磨液的質(zhì)量分?jǐn)?shù)進(jìn)行實(shí)驗[15]，并對這三個要素進(jìn)行實(shí)驗驗證。

圖6(a)所示是原始實(shí)驗條件邊際譜，原始實(shí)驗條件是平均加載到每球的壓力是10 N，10%的研磨液質(zhì)量分?jǐn)?shù)，研磨轉(zhuǎn)速為10 r/min。

圖6(b)中平均加載到每球的研磨壓力增加到20 N，10%的研磨液質(zhì)量分?jǐn)?shù)，研磨轉(zhuǎn)速為10 r/min。

圖6(c)中平均加載到每球的壓力是10 N，10%的研磨液質(zhì)量分?jǐn)?shù)，研磨盤轉(zhuǎn)速減小到10 r/min。

圖6(d)中平均加載到每球的壓力是10 N，研磨液質(zhì)量分?jǐn)?shù)增加到30%，研磨轉(zhuǎn)速為10 r/min。

圖6　實(shí)驗信號邊際譜不同條件對比圖Fig.6　Marginal spectrum of different conditions

以轉(zhuǎn)速分析為例，隨著轉(zhuǎn)速的減小，邊際譜中250 Hz區(qū)域峰值減小，說明隨著轉(zhuǎn)速的下降，由研磨盤轉(zhuǎn)速帶來的振動在減小.轉(zhuǎn)速對250 Hz的振動頻率影響較大.根據(jù)文獻(xiàn)記載和實(shí)驗驗證，研磨壓力、研磨速度、研磨液的質(zhì)量分?jǐn)?shù)是影響磨球振動的三個影響因素[16]。

以研磨速度為例，與FFT分析方法進(jìn)行比較，每球10 N的研磨壓力，10%的研磨液質(zhì)量分?jǐn)?shù)，研磨轉(zhuǎn)速分別為30，10 r/min，研磨振動邊際譜如圖7(a，b)所示.利用FFT方法分析振動，假設(shè)研磨速度對振動產(chǎn)生影響，圖7(c，d)是不同研磨盤轉(zhuǎn)速的FFT功率譜圖.4個實(shí)驗有相同的每球10 N的研磨壓力，10%的研磨液質(zhì)量分?jǐn)?shù)，不同的研磨轉(zhuǎn)速分別是30，10 r/min。

圖7　EMD邊際譜圖和FFT功率譜圖比較Fig.7　EMD marginal spectrum and FFT power spectrum

隨著轉(zhuǎn)速的減小，邊際譜中250 Hz區(qū)域峰值減小，說明隨著轉(zhuǎn)速的下降，由研磨盤轉(zhuǎn)速帶來的振動在減??；利用FFT方法分析振動，在假設(shè)基礎(chǔ)上，通過單一變量多次實(shí)驗，雖然也能得到振動隨著轉(zhuǎn)速減小而減小，但是沒有邊際譜圖形象直觀。

綜合文獻(xiàn)和實(shí)驗驗證，研磨壓力、研磨速度和研磨液的質(zhì)量分?jǐn)?shù)是磨球機(jī)振動產(chǎn)生的三個影響因素.基于EMD的精密球研磨過程中振動信號監(jiān)測與分析方法能正確反應(yīng)研磨中產(chǎn)生的振動，符合振動監(jiān)測的要求。

3結(jié)論

提出了基于EMD的精密球研磨過程中振動信號監(jiān)測與分析技術(shù)，并對該技術(shù)中的理論分析、硬件設(shè)計與軟件設(shè)計進(jìn)行了研究.在D-800Ballpoli雙自轉(zhuǎn)磨球機(jī)上進(jìn)行實(shí)驗，實(shí)驗證明，該系統(tǒng)具有良好特性，實(shí)驗信號邊際譜很好地表現(xiàn)精密球研磨過程中振動信號的變化.結(jié)合文獻(xiàn)記載，研磨壓力、磨盤的研磨速度以及研磨液的質(zhì)量分?jǐn)?shù)是振動產(chǎn)生的三個影響因素.隨著研磨時間的增加，振動幅值越來越小.相對于傳統(tǒng)的FFT頻譜分析法需要在假設(shè)的基礎(chǔ)上多次實(shí)驗比較得到振動影響因素的方法，這種方法直觀有效，符合振動監(jiān)測的要求，對降低研磨振動，提高精密球的研磨效果有很大幫助.通過對磨球機(jī)的狀態(tài)監(jiān)控可以有效的對球體的加工情況進(jìn)行了解，改進(jìn)精密球的加工工藝。

參考文獻(xiàn)：

[1]夏其表.精密球體研磨加工的在線檢測與優(yōu)化控制[D].杭州：浙江工業(yè)大學(xué),2008。

[2]WON J K, MODARRES M. Improved Bayesian method for diagnosing equipment partialf ailures in process plants[J]. Computers and Chemical Engineering,1998,22(10):1483-1502。

[3]呂冰海,袁巨龍,戴勇.氮化硅陶瓷球研磨過程中磨損形式的研究[J].摩擦學(xué)學(xué)報，2008,28(5):416-421。

[4]曹沖鋒.基于EMD的機(jī)械振動分析與診斷方法研究[D].杭州：浙江大學(xué),2009。

[5]王強(qiáng).基于壓電材料的振動能量采集技術(shù)的研究[D].南京：南京航空航天大學(xué)，2008。

[6]閆鍇.基于EMD的機(jī)械振動分析與診斷方法研究[J].科技資訊,2014(1):47-48。

[7]YU Danjiang, REN Weixin. EMD-based stochastic subspace identification of structures from operational vibration measurements[J].Engineering Structures,2005，27(27):1741-1751。

[8]DU Qinhua, YANG Shunian. Application of the EMD method in the vibration analysis of ball bearings[J]. Mechanical Systems and Signal Processing,2007,21(6):2634-2644。

[9]李穎.基于LabVIEW電動機(jī)振動檢測系統(tǒng)的研究[D].石家莊：河北農(nóng)業(yè)大學(xué),2012。

[10]王明哲.基于Atom和ARM的實(shí)時振動檢測儀[D].哈爾濱：哈爾濱工程大學(xué),2012。

[11]于洋，楊金英.機(jī)械振動檢測儀器的研制[J].儀表技術(shù)與傳感器,2008(9):35-36。

[12]李偉光，張送，朱火美，等.基于虛擬儀器的電主軸振動檢測[J].機(jī)械工程師,2009(2):93-94。

[13]周見行，李存兵，季行建，等.基于虛擬儀器技術(shù)輸油泵遠(yuǎn)程狀態(tài)監(jiān)測系統(tǒng)[J].浙江工業(yè)大學(xué)學(xué)報,2009(4)：469-472。

[14]安平.高精密球加工過程中的狀態(tài)監(jiān)測[D].杭州：浙江工業(yè)大學(xué),2013。

[15]顧盛挺，鮑雨梅，柴國鐘，等.顆粒表面粗糙度對納米壓痕特性影響研究[J].浙江工業(yè)大學(xué)學(xué)報,2013(2)：214-217。

[16]周文委；翁國慶；張有兵.基于ARM的多功能三相電能質(zhì)量在線監(jiān)測器[J].浙江工業(yè)大學(xué)學(xué)報,2012(2)：213-216。

(責(zé)任編輯：陳石平)

中圖分類號：G356.6

文獻(xiàn)標(biāo)志碼：A

文章編號：1006-4303(2015)02-0175-05

作者簡介：趙文宏(1969—)，男，浙江浦江人，教授級高級工程師，主要從事超精密加工設(shè)備的控制系統(tǒng)的設(shè)計，E-mail：916666039@qq.com。

基金項目：國家自然科學(xué)基金資助項目(51375455)

收稿日期：2014-11-10