李瑞敏，侯延星，侯春華

(沈陽機床(集團)有限責(zé)任公司設(shè)計研究院，沈陽 110142)

平面度的水平儀檢測及其數(shù)據(jù)的MATLAB處理*

李瑞敏，侯延星，侯春華

(沈陽機床(集團)有限責(zé)任公司設(shè)計研究院，沈陽 110142)

針對機床裝配工程實踐中經(jīng)常遇到的平面度檢測及修合問題，結(jié)合一具體實例詳細介紹了平面度水平儀檢測、最小二乘法評定及用MATLAB進行可視化數(shù)據(jù)處理的整套方法。該方法不僅可以快速方便地獲得平面度值，而且可以進一步為人工刮研平面等操作提供直觀的圖形依據(jù)。

平面度;水平儀;最小二乘法;MATLAB

0 引言

機床裝配工程實踐中，經(jīng)常遇到需要對平面進行平面度檢測及修合的問題。根據(jù)被測平面精度等級、尺寸大小和生產(chǎn)批量的不同，平面度的檢測方法和檢測設(shè)備也不同。常用的檢測方法大體可以分為直接測量法和間接測量法。直接測量法可直接獲得被測平面多點坐標或直接評定其平面度誤差，一般是將被測平面與理想平面進行比較。間接測量法是測量設(shè)備的讀數(shù)值需經(jīng)數(shù)據(jù)處理才能獲得平面多點坐標值的測量方法。直接或間接獲得被測平面多點坐標后，其平面度的評定方法也有多種，包括最小包容區(qū)域法、對角線平面法、三遠點平面法及最小二乘法等［1-3］。

本文將結(jié)合一具體實例介紹用水平儀間接測量、最小二乘法評定、MATLAB進行可視化數(shù)據(jù)處理的平面度檢測方法。

1 平面度的水平儀檢測

首先將平臺調(diào)整到大致水平位置，以保證水平儀在每一測量位置都有具體示數(shù)，然后按照一定的布點形式在平臺上首尾銜接地拖動固定有水平儀的橋板并記錄水平儀的讀數(shù)，最后根據(jù)布點形式的不同，采用一定的方法將讀數(shù)值轉(zhuǎn)換成坐標值進而進行平面度的評價［4-5］。

現(xiàn)有2m×2m的工作臺需要進行平面度檢測，將其放置好并大致調(diào)平，將精度P=0.02mm/1m的水平儀固定在一長L=0.2m的橋板上對其進行平面度檢測，如圖1所示。

圖1 平面度誤差的水平儀檢測

采用如圖2所示的布點形式，長寬方向各9m，測量順序為:

圖2 布點形式

表1 水平儀示數(shù)值g(單位：格數(shù))

表2 各點高度值h及DC向水平儀示數(shù)(單位：格數(shù))

將表2數(shù)據(jù)乘以4×10－6m(P×L)，即為圖2中各點相對于基面的真實z坐標，然后用MATLAB中的scatter3函數(shù)將各點(xi，yi，zij)描繪于空間坐標系中，并用meshgrid和surf函數(shù)［6-7］粗略擬合出實測平面如圖3所示。

2 平面度的最小二乘法評定

最小二乘評定法即以最小二乘中心平面為基面進行平面度評定的方法，所謂最小二乘中心平面是指空間點距其距離的平方和最小的平面［1］。

平面方程的一般表達式為:

圖3 實測平面

可算出所有點到最小二乘中心平面的距離，其中一側(cè)的點距其距離為正，另一側(cè)的點距其距離為負。實際平面的平面度為:max(d) － min(d)［1，10］。針對上述例子，用MATLAB編程將表2中的數(shù)據(jù)代入式(2)，然后用左除函數(shù)“”或線性回歸函數(shù)“regress”即可解出:所以表2中各點的最小二乘中心平面方程為:

最后用“meshgrid”及“mesh”函數(shù)將此平面繪制到圖3所示的空間坐標系中，結(jié)果如圖4所示。

圖4 實測平面及最小二乘中心平面對比

根據(jù)式(3)用MATLAB編程計算表2中100個點到最小二乘中心平面的距離dij，可得max(dij)=1.0562 ×10－5，min(dij)= －8.3618 ×10－6，所以平面度值max(dij)－min(dij)=1.9×10－5m。圖4中高點及低點位置一目了然，如若平面度值不符合要求，對高點進行刮研即可。

3 結(jié)束語

尺寸較大工作臺平面的平面度檢測一般采用間接測量法，測得的數(shù)據(jù)需要按照一定的方法經(jīng)過一系列變換和計算才能得到平面度值。MATLAB是一個功能強大的數(shù)學(xué)計算工具，帶有一個極大的預(yù)定義函數(shù)庫，其中還包含許多繪圖和圖象處理命令。用MATLAB對平面度間接測量所得數(shù)據(jù)進行處理，可以非常簡單地獲得平面度值并繪制出形象直觀的圖形，大大地方便了工作臺平面的平面度鑒定和修合工作。

［1］中國標準出版社.GB/T 11337－2004平面度誤差檢測［S］.北京:中國標準出版社，2004.

［2］廖念釗.互換性與測量技術(shù)基礎(chǔ)［M］.北京:中國計量出版社，2002.

［3］佘勃強，崔亞輝，劉西安.中大平面的平面度檢測［J］.裝備制造技術(shù)，2007(2):44－45.

［4］現(xiàn)代實用機床設(shè)計手冊編委會.現(xiàn)代實用機床設(shè)計手冊(上冊)［M］.北京:機械工業(yè)出版社，2006.

［5］杜明芳.平臺平面度檢測及數(shù)據(jù)的計算機處理［J］.北京印刷學(xué)院學(xué)報，2001，9(1):22 －24.

［6］陳杰.MATLAB寶典［M］.北京:電子工業(yè)出版社，2007.

［7］［美］Gerald Recktenwald，著，伍衛(wèi)國，等，譯.數(shù)值方法和MATLAB實現(xiàn)與應(yīng)用［M］.北京:機械工業(yè)出版社，2004.

［8］田樹耀，黃富貴，等.一種基于MATLAB的形位誤差評定方法［J］.工具技術(shù)，2008(4):96－98.

［9］周景亮，林志熙.大型平面平面度誤差的可視化數(shù)據(jù)處理［J］.機床與液壓，2009，37(10):130 －132，136.

［10］陳永鵬.基于MATLAB的平面度評定方法［J］.工具技，2004(2):52－54.

(編輯 趙蓉)

Flatness Measure w ith Gradienter and The Data Processing w ith MATLAB

LIRui-min，HOU Yan-xing，HOU Chun-hua
(Design and Research Department，Shenyang Machine Tool(Group)Co.，Ltd.，Shenyang 110142，China)

Aimed at the problems of flatnessmeasure and plane repair encountered in the machine tool assembly and the like engineering practice，onemethod of flatnessmeasure w ith grandienter，least square evaluation and visual data processing w ith MATLAB is introduced.W ith themethod，not only the flatness value can be acquired fast and easily，but also the graphics reference can be provided for scraping and the like flatness improvementwork.

flatness;gradienter;least squaremethod;MATLAB

TG83;TP391.7

A

1001－2265(2013)03－0039－03

2012－08－02;

2012－10－08

國家高檔數(shù)控機床與基礎(chǔ)制造裝備科技重大專項資助項目(2010ZX04140－014)

李瑞敏(1984—)，男，內(nèi)蒙古烏蘭察布人，沈陽機床(集團)有限責(zé)任公司設(shè)計研究院工程師，碩士，主要從事數(shù)控機床等機電一體化系統(tǒng)設(shè)計與研究工作，(E－mail)lrm2002@163.com。