張 霞, 蔡順燕, 李紅梅, 林 梅
(成都師范學院 物理與工程技術學院, 成都 611130)
磨削是利用砂輪中的磨粒進行微量切削加工的過程,砂輪表面的三維形貌對磨削力、材料去除方式、磨削溫度和工件表面質量等有重要影響[1-2]。準確有效地測量和分析砂輪表面的三維形貌,是提高磨削工件表面質量、優化砂輪修整工藝參數和磨削加工工藝參數的重要依據。特別是一些特殊工件材料的加工,方便、準確地測量加工時砂輪的表面三維形貌參數,將對整個磨削加工過程有重要指導作用[3-4]。
目前,很多國內外學者對此進行了大量研究?;豇P偉等[5-6]使用基于白光干涉原理的三維表面輪廓儀精確測量砂輪的表面形貌,取得了較好效果;XIE等[7-8]利用掃描電子顯微鏡對比了不同砂輪的表面形貌;CAI等[9-11]利用激光共聚焦三維測量儀對砂輪表面形貌進行測量,并計算了其表面形貌參數。
常用測量方法的測量過程往往復雜耗時,還需要將砂輪從機床上拆卸后測量,無法實現砂輪的在線測量。如何方便準確地在線測量機床砂輪的表面三維形貌參數并對其評價,是目前急需解決的磨削問題。
隨著現代超精密測量儀器的發展,采用二維(2D)位移傳感器可實現機床砂輪表面形貌的在線測量。測量中砂輪以一定速度旋轉,設置二維位移傳感器的采樣頻率來獲得砂輪表面的三維(3D)形貌數據[12]。李廈等[12-13]采用激光位移傳感器測量氧化鋁砂輪的三維形貌,獲取了砂輪局部表面形貌模型,評價了多粒度號氧化鋁砂輪的磨粒密度、出刃高度、有效出刃高度等特性參數;但該方法計算過程復雜,需要擬合砂輪出刃高度曲線,其結果的準確性對實驗環境的依賴性較高,難以在實際生產中應用。
為此,提出一種簡單有效分析砂輪表面三維形貌的方法。通過對砂輪表面三維形貌的基準平面進行計算和數字濾波后,獲得砂輪表面三維形貌的評價參數,可對砂輪形貌特征進行描述,最后通過實驗驗證該方法的有效性。
砂輪三維形貌的顯示和評定需要有一個基準面。實際被測砂輪表面不能作為基準面,而砂輪的設計表面是理想表面,在實際使用后其位置變化,也不能作為基準面。我們用最小二乘法來確定被測砂輪的基準平面。圖1為砂輪表面形貌測量示意圖。
如圖1所示:設x正方向為沿砂輪軸線向左方向,y正方向為砂輪表面圓周切線旋轉方向,在x和y方向上的采樣數據點間隔Δx和Δy分別為:
(1)
將2D激光位移傳感器測量的(xi,yj)坐標處的實際高度值定義為f(xi,yj)(xi=iΔx,yj=jΔy;i=1,2,...,M;j=1,2,...,N)
式中:Ls為x方向采樣數據點總長度,υs表示砂輪轉速,fk表示采樣頻率,M和N分別表示在x和y方向的總采樣點數。
設最小二乘法的基準平面函數:
z(x,y)=ax+by+c
(2)
式中:a,b,c分別為實數系數。
設ξ(xi,yj)表示殘差表面函數,是原始數據與基準平面之間的差異,如式(3)所示:
ξ(xi,yj)=f(xi,yj)-z(xi,yj)
=f(xi,yj)-(axi+byj+c)
(3)

(4)
(5)
求解式(5),可求出a,b,c的值;然后代入式(3)和式(4),可得到最小二乘法的基準平面函數z(x,y)和殘差表面數據ξ(xi,yj)。殘差表面數據ξ(xi,yj)將用作濾波和計算三維形貌參數的數據源。
由于實驗測量的砂輪三維形貌數據含有系統干擾和砂輪表面的高頻成分,不能直接用于砂輪表面三維形貌的分析及評價,故采用頻譜分析和數字濾波方法對測得的原始數據進行預處理[13]。在歐式坐標系中,殘差ξ(xi,yj)的M×N數字陣列為ξ(xi,yj),其中xi=iΔx,yj=jΔy;i=1,2,…,M;j=1,2,…,N。
對其進行三維離散傅里葉變換,得到:
(6)
式中:p=0,1,…,M-1;q=0,1,…,N-1。
沿x,y方向的空間頻率up,vq如式(7)所示:
(7)
式中:p和q分別表示x和y方向的采樣點數。
為了去除高頻正弦分量,我們采用理想的低通濾波傳遞函數:
(8)
式中:ω1為低通濾波的截止頻率。
濾波系統的頻域輸出,如式(9)所示:
G(up,vq)=H1p(up,vq)·F(up,vq)
(9)
然后用傅里葉逆變換來重建砂輪表面的三維形貌數據,如式(10)所示:
(10)
式中fh(xi,yj)(xi=iΔx,yj=jΔy;i=1,2,…,M;j=1,2,…,N)是濾波后的表面殘差數據,將被用于砂輪表面三維形貌參數的計算。
砂輪表面的三維形貌參數是對砂輪表面的微觀幾何形狀特性的某些方面做出的精確描述,可有效反映出機床砂輪表面的微觀幾何形狀、不平度高度、寬度、砂輪磨損及摩擦潤滑等特性[3-4]。
1.3.1 與微觀不平度高度有關的參數
(1)表面形貌的算術平均高度ha:
(11)
(2)表面形貌的高度均方根偏差hq:
(12)
式(11)和(12)的ha和hq是對整個表面形貌高度信息的描述,是簡單常用的高度參數。
(3)表面形貌的最大峰值hp:
hp=max(fh(xi,yj))
(13)
即hp是測量區域內最高峰相對于基準平面的高度。
(4)表面形貌的最大谷值hv:
hv=min(fh(xi,yj))
(14)
即hv是測量區域內最低谷相對于基準平面的高度。
(5)表面形貌的最大高度值hz:
hz=hp-hv
(15)
即hz是測量區域內最高峰值和最低谷值的高度差。
式(13)、式(14)和式(15)中:hp可有效反映砂輪表面磨粒的切削特性;hv反映砂輪表面容納冷卻液的體積;hz是砂輪形貌起伏的最大高度,當砂輪發生機制磨損時,hz可快速表達出高度的變化。
1.3.2 與微觀不平度間距有關的參數
(1)表面形貌的均方根斜率dΔq:
(16)

(2)表面形貌的展開界面面積比Sdr:
界面元素為實際表面上點(xi,yj)處的最小采樣四邊形。界面元素的面積Aij為:
(17)
界面面積比Sdr如式(18)所示:
(18)
式(16)和式(18)中dΔq和Sdr值越大,說明表面形貌高度、寬度或兩者都顯著,與砂輪堵塞磨損有一定關系。
1.3.3 與微觀不平度形狀有關的參數
(1)表面高度分布的偏斜度Fsk:
(19)
(2)表面高度分布的峭度Fku
(20)
根據參考文獻[5]可知:式(19)中Fsk為負值時表面呈飽滿形,相應耐磨性好;為正值時表面呈空虛形。故Fsk能表達不同的表面磨損條件類型。式(20)中Fku是與正態分布形狀為標準來比較的參數,當峭度分布符合正態分布時Fku=3,大于3時表示峭度分布較陡,小于3時表示較為平坦。
圖2為砂輪表面形貌測量實驗實物圖。圖2中:實驗機床為德國Schleifring公司K-P48T平面磨床;磨削砂輪為單晶剛玉(SA)砂輪:型號WA60H10V,直徑330 mm;選用KEYENCE公司的超高精度2D激光位移傳感器,型號LJ-G015,波長650 nm,量程5.2 mm,分辨率0.2 μm。
測量時,首先將超高精度2D激光位移傳感器的測頭安置于被測機床砂輪表面的正下方,再調整其激光射線,使其垂直于被測砂輪表面。設定傳感器信號采樣頻率fk=250 Hz,傳感器每次采樣的數據點數為800。實際測量過程中,先啟動機床砂輪作低速轉動,通常設置砂輪轉速υs≤1 m/s,再觸發2D激光位移傳感器進行數據采集,同時將傳感器采集的數據信號通過控制器傳送到計算機進行分析處理。
在砂輪轉速分別為0.2 m/s和0.5 m/s 的條件下測量砂輪表面形貌,經1.1中的砂輪表面基準平面及1.2中的數字濾波計算后,用Matlab軟件處理實驗數據,獲得砂輪表面的三維形貌如圖3。
利用1.3中的砂輪表面形貌參數計算公式得出圖3的三維形貌參數評價如表1所示。從表1可以看出:在不同砂輪轉速下2組對應的各參數基本相同,不受砂輪轉速變化的影響;且表1的評價參數較為全面地描述了砂輪的表面形貌特征,為保證砂輪修整效果及工件磨削質量提供了基本指標和依據。
為驗證上述二維位移傳感器測量評價方法的有效性和實用性,利用美國Zeta自動三維測量系統對砂輪表面進行測量,并將該儀器的測量結果與二維激光位移傳感器的測量計算結果對比,結果如表2。Zeta三維測量系統及實測砂輪表面形貌過程如圖4。

表1 砂輪表面三維形貌評價參數表

表2 測量結果對比表

(a)Zeta自動三維測量系統Zetaautomaticthree-dimensionalmeasurementsystem(b)Zeta測量砂輪表面形貌過程Zetameasurementofsurfacemorphologyofgrindingwheel圖4 Zeta三維測量系統及實測砂輪表面形貌過程Fig.4Zetathree-dimensionalmeasuringsystemandtheprocessofmeasuringgrindingwheelsurfacetopography
表2中2種方法的測量結果基本相同,驗證了本文使用的2D激光位移傳感器測量計算砂輪表面形貌方法的可行性和實用性,對砂輪表面三維形貌測量的后續研究及評價有指導意義。
(1)提出了磨削加工過程中基于二維激光位移傳感器的砂輪表面形貌測量及評價方法,實現其在線測量及評價。
(2)根據二維激光位移傳感器的測量原理,建立三維形貌基準平面和數字信號濾波方法,通過測量數據來計算砂輪表面的三維形貌參數。
(3)用二維激光位移傳感器實驗測量砂輪表面,比較不同砂輪轉速下砂輪表面的形貌參數,并與Zeta自動三維測量系統的測量結果進行對比,驗證了該方法的有效性和實用性;該方法使用簡單,測量時不用拆卸砂輪,具有良好的應用前景。