于大國，寧 磊，孟曉華

(1.中北大學 機械工程與自動化學院，太原 030051； 2.中北大學 山西省深孔加工工程技術研究中心,太原 030051)

?

基于最小二乘法深孔軸線直線度誤差評定*

于大國1，2，寧 磊1，孟曉華1

(1.中北大學 機械工程與自動化學院，太原 030051； 2.中北大學 山西省深孔加工工程技術研究中心,太原 030051)

利用模態(tài)分析中的擬合圓法求得深孔測量截面的最小二乘圓，確定截面孔心坐標，借助MATLAB三維繪圖指令依次連接孔心坐標點得到實際深孔軸線；在得到實際深孔軸線的基礎上，對深孔軸線直線度誤差在任意方向上進行評定，采用最小二乘法原理建立數(shù)學模型求得深孔軸線直線度誤差值；由MATLAB軟件對實驗數(shù)據(jù)進行處理，所得誤差值與要求的誤差值相對比；結果表明，該方法簡單且易于計算機數(shù)據(jù)處理，能滿足深孔軸線直線度的檢測，具有較高的實際應用價值。

測量截面；擬合圓；深孔軸線

0 引言

科學技術迅速發(fā)展和生產(chǎn)水平不斷提高使得機加工工業(yè)的制造精度和檢驗測量器具有了很大程度的提高和改善，因而有關形位誤差的評定方法、評定理論、解算方法等問題的研究已成為計量學領域研究的熱點[1]。

同軸度誤差是指實際被測軸線相對于理想軸線的變動量,理想軸線的位置應符合最小條件[2]。深孔軸線相對本身而言不能獨立存在，只有通過軸孔圓柱面加以體現(xiàn)。因此，深孔軸線的直線度誤差就是實際軸孔圓柱各截面輪廓中心連線的直線度誤差[3]。

深孔零件的直線度誤差對其工作性能有很大影響，準確評定深孔零件的直線度誤差，不但可以作為深孔零件驗收合格的依據(jù)，還可以用來分析誤差產(chǎn)生的原因，為提高深孔零件加工精度以及裝配精度提供可靠的依據(jù)。

目前深孔軸線直線度誤差的評定方法有：最小包容區(qū)域法(最小條件法)、兩端點連線法及最小二乘法[4]。最小包容區(qū)域法符合國家標準精度也最高，但在很多情況下，尋找和判斷符合最小條件的理想要素很困難，所以在實際應用中采用后兩種評定方法更為普遍[5]。兩端點連線法與其它兩種算法相比，計算簡單易求出基準軸線，但精度較低,不適于現(xiàn)代精密測量[6]。最小二乘法理論成熟，算法簡便，誤差值僅略大于最小包容區(qū)域法,能夠滿足一般精度要求，在直線度誤差評定中得到廣泛的應用[7]。

如何評定深孔零件軸線直線度誤差，提高深孔零件加工質(zhì)量，是人們普遍關心的問題，也是該領域研究人員的研究熱點之一[8-9]。因此研究評定深孔零件軸線直線度誤差具有非常重要的意義。

1 最小二乘圓法求解深孔軸線

將被測深孔零件放在三坐標測量機工作臺上，測頭沿內(nèi)孔壁在軸線方向的最高和最低位置打點調(diào)整被測孔軸線同工作臺基本垂直，同時固定被測件，測量模型如圖1。將被測深孔零件沿孔軸向分成若干截面，每一截面布置若干采樣點，采樣點坐標為Pij(xij,yij,zj)(i=1,2,…，m；j=1,2,…，n)，測量模型中擬合圓圓心坐標為Cj(aj,bj,zj)，其中m為每一截面采樣點數(shù)，n為測量的截面數(shù)，zj坐標為測量過程沿z上升的距離。

圖1 測量模型

假設測量過程中取第k個截面(k

圖2 測量截面數(shù)學模型

圖中：

Pik——測量截面上的任意一點；

Ck——測量截面擬合圓的圓心；

R——最小二乘擬合圓的半徑。

如果加工的深孔軸線沒有出現(xiàn)直線度誤差，那么測量點一定落在半徑為R的圓周上。實際深孔加工過程中不可避免出現(xiàn)誤差，因此測量截面上的點不可能剛好落在擬合圓周上，必定會存在誤差eik。

設最小二乘擬合圓的方程為

(x-ak)2+(y-bk)2=R2

那么誤差eik為

eik=(xik-ak)2+(yik-bk)2-R2

當測量點在擬合圓內(nèi)部時eik<0，在外部時eik>0，由隨機誤差的特征可知，正負誤差會相互抵消，為了避免這種情況發(fā)生，我們可以將誤差eik取平方，依據(jù)最小二乘原理必須滿足誤差eik平方和最小。

具體運算過程如下：

分別對u，v，w求偏導有

即

轉(zhuǎn)換成矩陣為

借助于數(shù)學軟件MATLAB求出u，v，w后，得ak=-u/2，bk=-v/2，因此第k個截面擬合圓的坐標為Ck(-u/2，-v/2，zk)。

在上述理論基礎上，根據(jù)每個截面測量的點Pij(xij,yij,zj)分別計算出截面孔心坐標Cj(aj,bj,zj)，依次連接各個截面孔心坐標Cj就可以得到被測深孔零件的實際軸線。

2 最小二乘法評定深孔軸線直線度誤差

與直線度公差帶相對應，直線度誤差分為給定平面內(nèi)、給定方向和任意方向的直線度誤差等三種形式[10]。對于回轉(zhuǎn)精度較高的回轉(zhuǎn)體零件如深孔類零件，其安裝孔心線多用任意方向上的直線度公差控制。因此將深孔軸線直線度誤差的評定歸結為任意方向的直線度誤差評定。

基于最小二乘圓法進行數(shù)學建模，由MATLAB軟件計算我們可以得到各截面孔心坐標為：C1(a1,b1,z1),C2(a2,b2,z2),…,Cj(aj,bj,zj),…,Cn(an,bn,zn), 依次連接這些點可以得到實際深孔軸線ls，用平行于ls的最小二乘中線ll s面包容該實際深孔軸線時，取其中具有最小直徑的圓柱面的直徑φl s作為誤差值，如圖3所示,則該n個孔心坐標點的算術平均中心點Co(ao,bo,zo)的坐標為：

(1)

圖3 深孔軸線直線度誤差坐標系

根據(jù)誤差理論，所有截面孔心坐標值的算術平均中心與其理想直線最為接近，因此以算術平均中心Co為基點的最小二乘擬合直線ll s是一種更為理想的直線。

設擬合直線ll s的方向向量為(l,m,n)，并且直線穿過點Co(ao,bo,zo)，因此可以得到擬合直線ll s的方程為：

(2)

實際深孔軸線上的任意一點Cj(aj,bj,zj)到擬合直線ll s的距離dj為：

dj=

(3)

f=

(4)

于是以上問題轉(zhuǎn)化為求式(4)的最小值問題，令一階偏導數(shù)為零得

(5)

根據(jù)空間解析幾何理論，對空間擬合直線的方向數(shù)l,m,n進行歸一化處理，即

l2+m2+n2=1

(6)

聯(lián)立式(4)和(5)求解得

(7)

將各截面孔心坐標值Cj(aj,bj,zj)(j=1,2,…，n)代入由式(1)得到算術平均中心點Co(ao,bo,zo)，接著將得到的Co點坐標值以及各截面孔心坐標值Cj代入式(7)得到擬合直線方向數(shù)l,m,n，最后將方向數(shù)l、m、n和ao、bo、zo代入式(3)得到以aj、bj、zj為未知數(shù)的方程。

將各截面孔心坐標值分別代入dj，找出其中的最大值dmax，則深孔軸線的直線度誤差值φl s為

φl s= 2(max{dj} )(j=1,2,…,n)

3 實驗測量與結果

實驗采用孔徑D=20mm，長度l=500mm，直線度要求為0.1mm的某一深孔零件，用測量頭對其內(nèi)孔各個截面進行測量,取m=8，n=10。由于測量點比較多，將測量截面采樣點坐標值省略，通過MATLAB計算得到被測截面擬合圓心坐標值，如表1所示。

表1 被測截面擬合圓心坐標計算值(單位：mm)

借助MATLAB三維繪圖命令，輸入表1中的被測截面擬合圓心坐標計算值得到被測深孔零件的實際軸線三維曲線圖，如圖4所示。

圖4 實際深孔軸線三維曲線圖

由MATLAB軟件編程計算得到該深孔零件軸線直線度誤差值為0.085mm，與所要求的直線度誤差數(shù)值相差0.015mm。最小二乘法評定深孔軸線直線度誤差流程圖，如圖5所示。

圖5 深孔軸線直線度誤差評定流程圖

4 結論

綜上所述，本文采用最小二乘擬合圓法確定測量截面孔心坐標點并依次連接孔心坐標點得到實際深孔軸線，以最小二乘法為依據(jù)借助MATLAB編程對實際深孔軸線直線度誤差進行評定。通過實驗論證分析，該方法簡單，數(shù)據(jù)易于處理，適用范圍廣，能滿足深孔軸線直線度的檢測，具有較高的實際應用價值。

[1]ZHANG Si-rong.Research on the integrative modeling theory integrated management method and application of product information and process supporting concurrent design[D].Hangzhou: Zhejiang University,2008.

[2]翟旭軍,張小萍,周圣鏵.基于改進粒子群算法的同軸度誤差評定及其可視化[J].組合機床與自動化加工技術，2012(12):55-59.

[3]粱勝龍.精密配合孔軸線直線度誤差評定[J].機械設計與研究，2009,25(3):94-97.

[4]劉雁蜀，王天琦.管類零件軸線直線度的處理方法研究[J].機械設計與制造,2008(12):144-147.

[5]李小亭，王樹彩．長度計量[M]．北京：中國計量出版社,2002．

[6]李淑娟,劉云霞.基于坐標變換原理的最小區(qū)域法評定空間直線度誤差[J].計測技術,2006,26(1):24-25.

[7]黃富貴，崔長彩.評定直誤差的最小二乘法與最小包容區(qū)域法精度之比較[J]. 光學精密工程,2007,15(6):889-894.

[8]黃富貴,崔長彩.任意方向上直線度誤差的評定新方法[J]. 機械工程學報,2008,44(7):221-224.

[9]戴麗,顧建剛.最小二乘圓法及其在原木定心中的應用[J]. 木材加工機械, 2002,13(2):16-18.

[10]續(xù)永剛，向立明，高國生.零件直線度誤差虛擬檢測系統(tǒng)研究[J].制造業(yè)自動化，2010,32(7):100-103.

(編輯 李秀敏)

Error Evaluation Method for the Axis Straightness of Deep-hole Based on Least-square Principle

YU Da-guo1,2, NING Lei1, MENG Xiao-hua1

(1.College of Mechanical Engineering and Automation, North University of China,Taiyuan 030051, China；2. Engineering Research Center for Deep-hole Drilling of Shanxi Province, North University of China ,Taiyuan 030051, China)

The method of least square principle in error evaluation is deduced from the coincidence circularity of mode analysis. The coordinate of the center of the measuring cross sections can be obtained via the method and the real axis of deep-hole can be gat from connecting the center of the hole coordinate points by using three-dimensional drawing instructions of MATLAB. On the basis of the real axis of deep-hole, carrying on evaluate arbitrary spatial straightness error of the axis of deep-hole. According to the least square principle, we can create mathematical modeling to find the axis straightness error of deep-hole finally. The experimental data can be deal with MATLAB and the resulting error value is compared to the error value of being required. The result indicates that the method is easy for data processing by computer, can satisfy the deep-hole straightness detection and has a high application value.

measuring cross sections; coincidence circularity; axis of deep-hole

1001-2265(2014)01-0039-03

10.13462/j.cnki.mmtamt.2014.01.011

2013-05-03;

2013-06-01

國家自然科學基金項目(51175482)；山西省國際合作項目“中韓合作BTA深孔加工的動態(tài)穩(wěn)定性研究”(2012081030)；中北大學?；?2012111)。

于大國(1964—)，男，江蘇金湖縣人，中北大學副教授，博士，研究方向為深孔加工技術，(E-mail)909886509@qq.com;通訊作者：寧磊(1987—)，男，山東費縣人，中北大學機械工程與自動化學院碩士研究生，主要研究方向為深孔加工技術，(E-mail)nl_158@126.com。

TH122;TG65;

A