張 濤, 崔長彩, 葉瑞芳, 黃 輝, 黃春棋
(華僑大學機電及自動化學院,福建廈門 361021)
表面紋理的空間參數表征及在磨??煽嘏挪忌拜啽砻娴膽?/p>
張 濤, 崔長彩, 葉瑞芳, 黃 輝, 黃春棋
(華僑大學機電及自動化學院,福建廈門 361021)
磨粒的可控優(yōu)化排布能夠有效改善砂輪的磨削性能。為了表征磨??煽嘏挪嫉男Чu定排布工藝的質量,采用最新的表面三維形貌評定標準ISO 25178-2空間參數(最速衰減自相關長度Sal,表面紋理結構比率Str,表面紋理方向Std)對其表面進行分析。首先研究了空間參數的表征意義及其求解方法,然后用Matlab軟件模擬具有不同磨粒排布類型的砂輪表面,分析了各空間參數在磨??煽嘏挪忌拜啽砻娴膽?。結果表明空間參數可以應用于磨??煽嘏挪夹Ч谋碚骱驮u定,從而為磨??煽嘏挪忌拜喌臋z測和評定提供了一種新方法。
計量學;磨粒可控排布;ISO25178-2標準;空間參數;表征;評定
隨著對表面三維形貌表征和評定需求的日益增強,在有關學者研究成果[1,2]的基礎上,國際標準化組織從2002年開始制訂并陸續(xù)出版了關于表面紋理結構測量和分析的新標準ISO 25178(Geometrical product specifications(GPS)-Surface texture:Areal)[3]。
ISO 25178系列標準不僅涉及有關表面紋理結構評定等操作規(guī)范的定義,同時也包含了測量方法、校準方法以及硬件和軟件校準等規(guī)范。其中,ISO 25178-2定義了用于表面三維形貌表征分析的參數,主要包括幅值參數、空間參數、綜合參數、功能體積參數以及其它用于表面特征表征的特征參數[4]。其中,空間參數主要表示表面紋理在空間上的分布情況,也稱其為紋理參數,主要包括3個參數:Sal(最速衰減自相關長度)可以對表面紋理的頻率組成成分進行分析;Str(表面紋理結構比率)用來表征表面紋理在各方向上的一致性,可以用來對表面形貌各向同性和各向異性的程度進行描述[5];Std(表面紋理方向)用來表征表面紋理的方向。
傳統(tǒng)砂輪磨削面的磨粒是隨機分布的,為了改善砂輪的磨削性能,提高加工質量,有關學者對其表面磨粒的排布和三維形貌進行優(yōu)化設計,即根據實際加工需要使磨粒按照一定的規(guī)律排布[6]。因此這些砂輪的表面必然會存在某種特定的紋理,這些紋理的特性直接對應著磨??煽嘏挪嫉男Ч缒チE挪嫉木鶆蛐?、密度大小及排布方向等,故可以通過對磨??煽嘏挪忌拜啽砻婕y理特性的表征來評定磨粒的實際排布與所設計排布的一致性。
2.1Sal和Str的表征意義與求解
2.1.1Sal和Str的表征意義
ISO 25178-2定義Sal為自相關函數以最快速度衰減到指定值s(0≤s<1,標準里默認值為0.2)時對應的水平距離。Sal是分析表面紋理組成成分的一個重要參數,如果表面低頻成分占主導,則Sal較大;反之,若表面主要由高頻成分組成,Sal則較小。
Str標準定義為自相關函數以最快速度衰減到指定值s時對應的水平距離與自相關函數以最慢速度衰減到指定值s時對應的水平距離之比(0≤s<1)。Str可用于表面各向同性與各向異性的判斷:對于各向異性的表面,Str值趨于0,而對于各向同性的表面,Str值趨于1。并且有學者認為若Str>0.5,表明表面在不同的方向上有較統(tǒng)一的紋理形貌,若Str<0.3,則表明表面是各向異性的[7]。
2.1.2Sal和Str的求解
求解Sal和Str的方法是通過作出表面二維自相關函數圖像來求解,表面二維自相關函數是表面紋理分析的一種有效工具,它描述的是表面某處形貌與另一處形貌間的依賴關系,對于具有周期性和方向性紋理的表面,它的二維自相關函數同樣具有周期性和方向性[8]。Sal和Str求解過程示意圖見圖1。
首先根據自相關函數計算方法作出表面的自相關函數圖,見圖1(a)。得到表面的自相關函數圖后,通過取閾值s并經二值化處理后得到如圖1(b)所示的二值化圖像,然后濾去周邊部分后得到自相關函數的中心區(qū)域輪廓,如圖1(c)所示。圖像的中心(即自相關函數圖像上幅值為1的點)到中心區(qū)域邊界的最短距離Rmin即為自相關函數以最快速度衰減到所取的閾值s時對應的距離,即最速衰減自相關長度Sal=Rmin;圖像的中心到中心區(qū)域邊界的最長距離Rmax即為最慢衰減自相關長度,兩者的比值Rmin/Rmax即為表面紋理結構比率Str=Rmin/Rmax,如圖1(d)所示。
圖1Sal和Str的求解流程圖
2.2Std的意義與求解
表面紋理方向Std表征的是表面紋理相對于所選測量方向的角度。Std可以由表面角功率譜(angular spectrum)密度函數圖求出,所謂的角功率譜密度是由表面功率譜密度函數在各方向上(一般在0°~180°或0°~360°)積分得到:
式中G(r,θ)是表面功率譜密度G(ωp,ωq)的極坐標表達形式,通過坐標變換得到:R(θ)指的是G(r,θ)在θ方向上的半徑的最大值。
表面功率譜密度G(ωp,ωq)可通過表面形貌數據的二維傅里葉變換得到:
其中,F*(ωp,ωq)是傅里葉變換F(ωp,ωq)的復共軛,二維離散傅里葉變換為:
其中,p=0,1,…,M-1;q=0,1,…,N-1;ωp=p/ΔxM,ωq=q/ΔyN;z(xl+1,yl+1)表面形貌高度數據。
在得到角功率譜密度函數A(θ)以后,就可以直觀地看出功率譜密度在各個方向上的分布情況,其中,功率譜最大的方向便對應著表面紋理的方向。
對于一般的砂輪表面,磨粒是隨機分布的,一般不會存在某種特定類型的紋理,可認為是各向同性表面。為了提高砂輪特性和加工質量,可根據需要對磨粒(或磨粒群)按一定的規(guī)律進行排布,因而會在砂輪表面形成特定的紋理,圖2是一種磨??煽嘏挪忌拜喌膶嵨飯D。通過砂輪紋理特性的空間參數表征便可以對砂輪表面磨??煽嘏挪嫉男Чㄈ缒チ7植嫉木鶆蛐?、密度大小、磨粒尺寸大小以及排布的方向等),作出定性或定量的評定。為方便研究磨??煽嘏挪忌拜喌奶匦裕疚睦肕atlab軟件來模擬研究所需的、具有不同紋理特性的砂輪表面。
圖2 一種磨??煽嘏挪嫉纳拜?/p>
3.1Str與磨粒排布均勻性的關系
對于砂輪表面來說,磨粒排布越均勻其表面各向同性程度越高,反之磨粒排布越不均勻,其表面越趨于各向異性。因此可以用Str來分析磨??煽嘏挪嫉木鶆蛐浴D3是用Matlab分別模擬的磨粒均勻排布和非均勻排布的砂輪表面,兩幅圖中的磨粒數量和大小一致,只是排布的均勻性不一樣。表1是閾值s取不同值時,磨粒均勻分布與非均勻分布表面所對應的Str值,可以發(fā)現當s值在0.4附近時,兩個表面Str值的差距比較明顯,文中將閾值設為0.4。圖4、圖5分別是s取0.4時它們的表面自相關函數圖像及Str計算示意圖。
圖3 Matlab模擬的砂輪表面磨粒分布圖
表1 閾值s取不同值時磨粒均勻分布與非均勻分布表面的Str值
圖4 磨粒均勻排布砂輪表面(圖3(a))Str計算示意圖
圖5 磨粒非均勻排布砂輪表面(圖3(b))Str計算示意圖
從圖4與圖5的對比中可以看出:磨粒均勻排布的砂輪表面的自相關函數圖像中心區(qū)域邊界上各點到中心的距離趨于一致,磨粒非均勻排布的砂輪表面的自相關函數對應的區(qū)域邊界上各點到中心的距離差別則較大,由圖4計算得到的Str=0.69也明顯大于由圖5計算得到的Str=0.41。這與前文所述的Str值可以用來表征表面的各向同性和各向異性程度的性質一致,由此可見,Str值可以反映磨粒排布的均勻性程度,因此可以把它作為評定磨粒排布工藝質量的一個重要參數。
3.2 Sal與磨粒大小和分布密度的關系
Sal值大小與表面的頻率組成成分有關,而砂輪表面的磨粒與結合劑分別對應著不同的頻率成分,因此Sal值與磨粒的大小及分布密度也應該有一定的關系。對于圖3所示的磨粒均勻和非均勻分布的砂輪表面,分別保持磨粒的分布特性和磨粒的個數不變,通過增加和減小磨粒的粒徑可以得到如表2所示的Sal隨磨粒大小變化的兩組數據,圖6是對應的關系曲線圖,橫坐標表示磨粒粒徑的大小。同理,保持兩幅圖中磨粒的分布特性及其大小不變,改變相同面積(4 mm2)上磨粒的個數,即改變砂輪表面磨??煽嘏挪嫉拿芏?,可得到表3中兩組Sal隨磨粒分布密度變化的數據,圖7是對應的關系曲線圖,橫坐標表示特定大小的面積內所含磨粒的個數,即磨粒的分布密度。
表2 磨粒大小不同的砂輪表面對應的Sal值μm
表3 磨粒分布密度不同的砂輪表面對應的Sal值
圖6Sal與磨粒大小關系曲線
圖7Sal與磨粒分布密度關系曲線
通過以上的數據分析可以看出,Sal值與磨粒的大小及其在砂輪表面的分布密度有著非常明確的對應關系。從圖6可以看出,在磨粒合理大小的范圍內,磨粒越大Sal值越?。煌瑫r,表示均勻分布與非均勻分布表面Sal值的兩條曲線非常接近,表明Sal值與磨粒的分布特性關系不大。從圖7可以看出:在磨粒大小相同的情況下,磨粒分布越密(即密度越大),Sal值越小。因此,Sal可以用來表征磨粒的大小和密度,進而評定磨粒可控排布工藝的質量。
3.3 用Std表征磨??煽嘏挪嫉姆较?/p>
磨??煽嘏挪忌拜啽砻娴募y理主要是由磨粒按照一定的規(guī)則分布而形成,因而其表面紋理方向指的即是磨粒排布的方向??臻g參數Std用來表征表面紋理的方向,因此可以應用于磨粒可控排布方向的檢測。圖8(a)、(b)是用Matlab分別模擬的磨粒按照與圖示水平方向(測量方向)夾角為arctan3和arctan(1/3)排布的砂輪表面,按照2.2節(jié)中所述的Std的求解過程:首先求出它們的表面功率譜密度函數,然后根據表面功率譜密度函數作出對應的角功率譜密度圖,如圖9所示,通過編程計算可分別得到圖中功率譜密度最大的方向是251.55°和18.42°,即檢測出的磨粒排布方向α1和α2分別為71.55°和18.42°,與理論值arctan 3=71.56°和arctan(1/3)=18.43°均只相差0.01°。
圖8 磨粒按固定方向排布的砂輪表面
模擬計算的結果表明,通過對Std的求解可以準確、方便地檢測出砂輪表面磨??煽嘏挪嫉姆较颍恍枰涍^磨粒及磨粒位置的識別等處理步驟。故Std也可以作為表征磨粒可控排布效果的一個重要參數。
從模擬計算分析可以看出,ISO 25178-2里用于一般表面紋理表征和評定的空間參數可以很好地應用于磨??煽嘏挪忌拜啽砻?,通過對砂輪表面上述3個空間參數的計算分析就可以對磨粒排布的均勻性、密度大小、尺寸大小和排布的方向作出定量的表征,從而為磨??煽嘏挪脊に囐|量的評定提供依據。
[1] Stout K J,Sullivan P J,Dong W P,etal.The development of methods for the characterization of roughness in three dimensions[M].UK:University of Birmingham Press,1993.
[2] Blunt L,Jiang X Q.Advanced techniques for assessment surface topography—development of a basis for the 3D surface texture standards“SURFSTAND”[M].UK:Kogan Page Science,2003.
[3] Blateyron Francois.New 3D parameters and filtration techniques for surface metrology[EB/OL].http://www.qualitymag.com/articles/85037-new-3d-parametersand-filtration-techniques-for-surface-metrology.
[4] ISO/WD 25178-2:2006 Geometrical product specifications(GPS)—Surface texture:Areal—Part 2:Terms,definitions and surface texture parameters[S].
[5] Jiang X Q,Scott P J,Whitehouse D J,etal.Paradigm shifts in surfacemetrology,Part II:The current shift[J].ProceedingsoftheRoyalSocietyA,2007,463(2085):2071-2099.
[6] 陳建毅,黃輝,徐西鵬.釬焊超硬磨料工具的研究進展[J].工具技術,2007,41(2):9-14.
[7] 陳國強.基于GPS標準體系的表面微觀形貌三維表征研究[D].南京:南京農業(yè)大學,2007.
[8] Jiang X Q,Blunt L.Surface analysis techniques to optimise the performance of CNC machine tools[J].TransactionsonEngineeringSciences,2003,44:107-117.
[9] 喬磊.表面三維形貌分析與評定軟件系統(tǒng)研究[D].武漢:華中科技大學,2009.
Characterization of Surface Texture with Spatial Parameters and Its Application on Grinding Wheel Surface with Controllable Grains Arrangement
ZHANG Tao, CUIChang-cai, YE Rui-fang, HUANG Hui, HUANG Chun-qi
(College of Mechanical Engineering and Automation,Huaqiao University,Xiamen,Fujian 361021,China)
The optimized controllable arrangement of grains on a grinding wheel can improve its processing properties effectively.In order to characterize the arrangement effect and so as to evaluate the arrangement quality,spatial parameters(Sal—fastest decay auto-correlation length,Str—texture aspect ratio,Std—surface texture direction)in the latest 3D surface texture standard ISO25178-2 are adopted to analyze the grinding wheel surface.First,themeaning and calculation method of the spatial parameters are studied.Then grinding wheel surfaces with different grain arrangement types are simulated by Matlab and then the application of each spatial parameter on grinding wheelswith particular grain arrangement is analyzed.The application is proved to be reasonable as a result,which provides a new method for the test and evaluation of grinding wheels with controllable grains arrangement.
Metrology;Controllable grains arrangement;ISO 25178-2 standard;Spatial parameters;Characterization;Evaluation
TB92
A
1000-1158(2014)05-0425-05
10.3969/j.issn.1000-1158.2014.05.04
2013-05-29;
2014-02-14
國家自然科學基金(51235004,51075160);教育部新世紀優(yōu)秀人才支持計劃(NCET-10-0116)
張濤(1989-),男,安徽蚌埠人,華僑大學碩士,主要研究方向為智能測控儀器與測控方法。zhangtaoliu1989@163.com崔長彩為通訊作者。cuichc@hqu.edu.cn