郭敬濱，曹紅艷，王克新，齊永岳(. 天津大學精密測試技術(shù)及儀器國家重點實驗室，天津 30007；. 溫州市計量技術(shù)研究院，溫州 35000)

轉(zhuǎn)臺分度誤差的檢定及補償模型的建立

郭敬濱1，曹紅艷1，王克新2，齊永岳1
(1. 天津大學精密測試技術(shù)及儀器國家重點實驗室，天津 300072；2. 溫州市計量技術(shù)研究院，溫州 325000)

用誤差補償?shù)姆椒ㄌ岣咿D(zhuǎn)臺分度精度是一種經(jīng)濟、有效的技術(shù)措施，用比較法對轉(zhuǎn)臺的分度誤差進行了檢定，并分析了其測量不確定度；從轉(zhuǎn)臺分度誤差的來源出發(fā)，根據(jù)轉(zhuǎn)臺分度誤差的周期性等特點和諧波分析的結(jié)果，并結(jié)合最小二乘法、三次樣條插值、傅里葉級數(shù)擬合3種模型的特點，通過分析，選擇傅里葉級數(shù)作為轉(zhuǎn)臺分度誤差的補償模型，并對3種模型的擬合效果進行了比較．應用傅里葉級數(shù)模型對精度為17.82″ 的轉(zhuǎn)臺進行補償，補償后轉(zhuǎn)臺的分度誤差最大值為2.36″，從而驗證了分度誤差模型的正確性和可行性．

分度誤差；檢定；數(shù)據(jù)擬合；傅里葉級數(shù)；補償模型

轉(zhuǎn)臺是進行圓分度的常用設備，經(jīng)常作為機床附件進行零件的圓分度加工．為了測量齒輪等回轉(zhuǎn)體零件，也常用到轉(zhuǎn)臺．隨著科學技術(shù)的不斷發(fā)展，對齒輪測量精度的要求也越來越高[1-2]．轉(zhuǎn)臺精度會影響齒輪等回轉(zhuǎn)體零件的測量精度，通過軟件誤差補償?shù)姆椒梢院唵巍⒂行У靥岣咿D(zhuǎn)臺的精度，以較低的成本達到較高的精度要求．

由于在轉(zhuǎn)臺的分度誤差中，系統(tǒng)誤差占主要成分，從而可以用非實時誤差補償?shù)姆椒▽D(zhuǎn)臺分度誤差進行修正，基本思想是以精度更高的儀器作為基準，將轉(zhuǎn)臺分度誤差先檢定出來，選擇合適的擬合或插值函數(shù)作為轉(zhuǎn)臺分度誤差的補償模型，對轉(zhuǎn)臺分度誤差進行補償．由此可知，在轉(zhuǎn)臺分度誤差補償中，誤差的檢定和模型的選擇將直接影響補償效果，是分度誤差補償?shù)年P鍵．

針對采用多面棱體檢定時受多面棱體面數(shù)限制測量間隔不能太小的問題，筆者提出采用多齒分度臺與平面反射鏡組合構(gòu)建可回轉(zhuǎn)的多角度基準的檢定方法，并對該方法的測量不確定度進行了分析．對采用最小二乘多項式擬合、傅里葉級數(shù)擬合和三次樣條插值3種補償模型進行了分析對比，由于分度誤差具有周期性特點，采用傅里葉變換的方法建立補償模型效果更好．

1 分度誤差檢定及測量不確定度分析

1.1 分度誤差來源

轉(zhuǎn)臺的分度誤差主要包括圓光柵本身的制造誤差和機械裝配誤差等．

制造圓光柵的方法有機械刻畫法和光電式刻畫法．采用機械法刻畫圓光柵時，其刻畫誤差與圓刻機的蝸輪副誤差、傳動鏈誤差、調(diào)整誤差等有關[3]．目前，常采用光電式刻畫法，主要有接觸式光刻法、接近式光刻法及投影式光刻法[4-5]，光電式刻畫法精度高，不受齒輪傳動系統(tǒng)誤差、蝸輪副的制造偏差等影響，其制造精度與制造工藝有關．圓光柵的制造誤差是以2π為周期的函數(shù)，因此可以將其展成傅里葉級數(shù)，即多次諧波疊加的形式．

裝配誤差包括圓光柵的安裝偏心、傾斜等，安裝偏心會引起一次諧波誤差，傾斜會引起二次諧波誤

差[6-9]．

以上各誤差分量都具有周期性，空間基周期為2π，因此，轉(zhuǎn)臺的分度誤差函數(shù)為多次諧波的疊加．

1.2 誤差的檢定

為了補償轉(zhuǎn)臺的分度誤差，首先需要對其進行檢定．常用的檢定方法是以正多面棱體作為圓分度標準器，用自準直儀進行讀數(shù)測量．為了提高檢定精度，常采用全組合測量法，但全組合測量法工作量大，而且受多面棱體面數(shù)的限制，測量間隔不能太?。壳俺Ｓ玫亩嗝胬怏w面數(shù)最多為72面，可達到的最小測量間隔為5°，且一個多面棱體只能實現(xiàn)一種或幾種角度間隔的測量．多齒分度臺具有較易獲得高的分度精度、對使用環(huán)境要求不高、操作簡便等優(yōu)點，作為圓分度標準器可以實現(xiàn)較高的精度[10-12]．因此，本文利用多齒分度臺與平面反射鏡組合，形成一種可回轉(zhuǎn)的多角度基準，能夠測量正多面棱體不能實現(xiàn)的角度間隔，且工作量小、精度高．

用多齒分度臺檢定轉(zhuǎn)臺分度誤差的基本思想是：以多齒分度臺為分度基準件、其回轉(zhuǎn)的角度作為標準角度，轉(zhuǎn)臺實際轉(zhuǎn)過的角度與標準角度比較，由光電自準直儀讀取角度偏差．其測量原理如圖1所示．

多齒分度臺安裝在轉(zhuǎn)臺工作面的中心，平面反射鏡固定在多齒分度臺上，利用自準直儀進行瞄準．測量時，先將被檢轉(zhuǎn)臺置零位，讀取自準直儀的初始值f0，然后多齒分度臺及反射鏡隨轉(zhuǎn)臺一起轉(zhuǎn)動φ角度(φ=360°/n，n為整數(shù)，且需滿足φ是多齒分度臺最小分度間隔的整數(shù)倍，實驗中采用360齒精密多齒分度臺，可實現(xiàn)的最小間隔為1°)，當反射鏡隨多齒分度臺向相反方向回轉(zhuǎn)相同的角度后，光電自準直儀再次瞄準反射鏡并讀取測量值f1．兩次讀數(shù)之差即為該角度間隔的分度誤差Δφ=f1-f0=φ′-φ(其中φ′為轉(zhuǎn)臺實際轉(zhuǎn)過的角度)，如此順序?qū)Ω鞅粰z角度進行測量(記第i個測量點讀數(shù)為fi)，直至整周．各測量點讀數(shù)與初始值之差xi(xi=fi-f0)即是轉(zhuǎn)臺相對于零位的分度誤差．

圖1 檢測裝置示意Fig.1 Schematic diagram of detection device

1.3 不確定度分析

當多齒分度臺存在安裝偏心時，通過分析可知，安裝偏心不會引起自準直儀示值變化，僅使光線照射在反射鏡上的位置發(fā)生了改變ΔL．如果反射鏡存在平面度誤差，由位置變化引起的誤差與反射鏡的平面度有關．因此影響測量結(jié)果的因素主要有多齒分度臺的分度誤差和重復性誤差、光電自準直儀的示值誤差和穩(wěn)定性誤差、反射鏡的平面度誤差以及外界環(huán)境的影響等．

本實驗使用的360齒精密多齒分度臺分度誤差u1=0.3′′，重復性誤差u2=0.18′′；雙軸光電自準直儀測量精度為u3=0.5′′，重復性誤差u4=0.3′′；反射鏡反射面寬度H=60,mm，平面度誤差δ=λ/8≈0.079 μm，當光線偏離原來位置的距離ΔL=2,mm時，引入的誤差u5=8δΔLρ/H2=0.07′′，其中ρ=206265′′，表示1弧度對應的秒值．

由于上述各項不確定度分量相互獨立，則合成不確定度(即測量不確定度)為

2 分度誤差補償模型的建立

轉(zhuǎn)臺分度誤差檢定后，給出的是一系列離散點的分度誤差，而在實施誤差修正時，補償點可能不在測量點上，需要對檢定結(jié)果進行插值或擬合處理，找出轉(zhuǎn)臺的分度誤差X與轉(zhuǎn)角θ之間的關系，建立數(shù)學模型，從而求得轉(zhuǎn)臺任意補償點處的分度誤差．數(shù)學上常用的擬合方法有最小二乘法、三次樣條法和分段多項式插值等．

最小二乘多項式擬合可以得到較高的精度，表達式唯一，易于計算，但當擬合的曲線比較復雜時，多項式次數(shù)高，會產(chǎn)生Runge現(xiàn)象，導致擬合精度較低．因此，對于比較復雜的曲線，不適宜采用最小二乘法多項式擬合．

在用插值算法時，高次多項式插值會產(chǎn)生振蕩現(xiàn)象，數(shù)值穩(wěn)定性差，導致插值誤差增大；分段低次多項式插值簡單、收斂性好，但在節(jié)點處光滑性較差．三次樣條插值(記為S(θ))不存在這些問題，它在每個子區(qū)間[θi,θi+1]上都是3次多項式，即S(θ)=aθ3+bθ2+cθ+d．三次樣條插值可以用較低的階數(shù)來對比較復雜的數(shù)據(jù)進行擬合，并且在節(jié)點處具有二階連續(xù)導數(shù)，從而可以保證光滑連接．

根據(jù)檢定得到的各節(jié)點處的分度誤差，可得插值函數(shù)

式中：θ為轉(zhuǎn)臺轉(zhuǎn)角；θi、θi+1分別為第i、i+1個檢定點處的轉(zhuǎn)角；hi=θi+1-θi．以第一類邊界條件為例，mi可根據(jù)線性方程組(2)求得．

根據(jù)檢定得到的數(shù)據(jù)，由式(2)即可求得mi，代入式(1)可得到各子區(qū)間[θi,θi+1]上的3次多項式．為了求得轉(zhuǎn)臺轉(zhuǎn)角θ處的分度誤差，首先要判斷θ所在的區(qū)間，再將θ代入該區(qū)間上的3次多項式即可．

通過以上分析可知，三次樣條插值計算過程復雜，且由于在實測數(shù)據(jù)中總是存在誤差，對數(shù)據(jù)進行三次樣條插值，隨機誤差沒有被均化，影響擬合精度.

傅里葉級數(shù)是以三角函數(shù)表示的無窮級數(shù)，特別適合對周期性函數(shù)進行描述，在工程技術(shù)領域應用廣泛[13-15]．設轉(zhuǎn)臺的分度誤差為()Xθ，由于它是以2,π為周期的函數(shù)，根據(jù)傅里葉級數(shù)展開和收斂定理可知，()Xθ可以展開成傅里葉級數(shù)，即

式中：k為諧波次數(shù)；a0、ak、bk為傅里葉系數(shù)．

式(3)由無限個諧波疊加而成，但由于在實際檢定時測量點數(shù)n有限，因此只能取有限項諧波疊加來逼近X(θ)．設傅里葉級數(shù)擬合中周期項個數(shù)為m(m≤[n/2])，則擬合得到的轉(zhuǎn)臺分度誤差函數(shù)F(θ)為

根據(jù)測量得到的離散數(shù)據(jù)，由式(5)求出傅里葉級數(shù)擬合方程中各系數(shù)，將各次諧波疊加，即可得到連續(xù)函數(shù)F(θ)．

由式(4)和式(5)可知，傅里葉級數(shù)擬合模型計算過程簡單，能很好地反映數(shù)值的周期現(xiàn)象，對復雜曲線的擬合能力強；根據(jù)擬合函數(shù)中各次諧波的幅值大小，可以為加工者尋找、分析誤差源，也可為使用者尋求最佳使用方法提供有用的信息，而且傅里葉級數(shù)擬合可以抑制測量數(shù)據(jù)中隨機誤差的影響．在擬合過程中由于m的取值受測量點數(shù)n的限制(m≤n/2)，因此，擬合精度與n有關，當精度要求較高時，需要通過適當增大n值來獲得較高的擬合精度．

通過以上的分析可知，3種模型各有優(yōu)缺點，實際應用時應根據(jù)擬合對象的特點來選擇合適的擬合模型．

誤差檢定時選定測量間隔為5°，測量得到的離散數(shù)據(jù)及經(jīng)諧波分析得到的各次諧波幅值大小如圖2所示．

圖2 轉(zhuǎn)臺分度誤差檢定結(jié)果Fig.2 Test results of indexing error

表1 3種模型擬合結(jié)果對比Tab.1 Comparison of three kinds of fitting models (″)

由圖2可知，檢定得到的離散數(shù)據(jù)比較復雜，因此，不適宜采用最小二乘多項式擬合．轉(zhuǎn)臺分度誤差為周期函數(shù)，是以2π為基周期的多次諧波的疊加，而三次樣條插值不能反映數(shù)據(jù)的周期性．由圖2(b)諧波分析結(jié)果可以看出，轉(zhuǎn)臺分度誤差以低次諧波為主，由于傅里葉級數(shù)適用于對周期函數(shù)進行描述，對復雜曲線的擬合能力強，針對轉(zhuǎn)臺分度誤差的周期性、誤差來源以及數(shù)據(jù)的復雜程度等特點，選擇傅里葉級數(shù)作為轉(zhuǎn)臺分度誤差的擬合模型．

為了驗證傅里葉級數(shù)擬合方法的有效性和可行性，本文針對檢定數(shù)據(jù)分別建立了轉(zhuǎn)臺分度誤差的最小二乘多項式擬合模型、傅里葉級數(shù)擬合模型和三次樣條插值模型，并對它們的擬合效果進行了比較．

為了使比較具有一般性，從0°起每隔18°或17°選取比較點，從而使選取的點既包括檢定點也包括中間點．分別根據(jù)3種模型計算各點的分度誤差值，與實測得到的數(shù)據(jù)對比，求得各點的殘差，數(shù)據(jù)見表1，殘差曲線如圖3所示．

擬合精度由殘差平方和來評定，認為殘差平方和越小，曲線的擬合效果越好[15]．通過計算求得最小二乘法擬合殘差平方和為17.08，傅里葉級數(shù)擬合為1.35，三次樣條插值為0.83．通過對比可以看出，傅里葉級數(shù)比最小二乘多項式擬合精度高，與三次樣條插值模型的擬合結(jié)果很接近．由于三次樣條插值函數(shù)需要求解線性方程組，且對于給定的θ角，需要判斷θ所在的區(qū)間，計算過程復雜[16]；傅里葉級數(shù)擬合中，根據(jù)式(5)可求得各次諧波的系數(shù)，再將θ代入式(4)求解即可得到轉(zhuǎn)臺轉(zhuǎn)角為θ時的分度誤差，易于計算，且能體現(xiàn)分度誤差的周期性，因此，選擇傅里葉級數(shù)()Fθ作為轉(zhuǎn)臺分度誤差的擬合模型，是一種簡單有效的數(shù)據(jù)處理方法．

圖3 殘差曲線Fig.3 Residual curves

3 實驗驗證

在對分度誤差進行檢定和擬合完成以后，就可以對其進行補償．誤差補償后對轉(zhuǎn)臺的分度誤差進行5次測量，并對轉(zhuǎn)臺反轉(zhuǎn)時的補償效果進行了實驗．補償前后的分度誤差曲線如圖4所示(補償后測量點數(shù)n為144)．可以看出，補償前轉(zhuǎn)臺的分度誤差最大值為17.82''；補償后正轉(zhuǎn)和反轉(zhuǎn)時分度誤差大部分分布在±0.5"之間，且具有隨機性，與補償之前相比，定位精度得到了明顯的改善．5次測量中最大分度誤差的平均值為2.36''，重復性標準差(用實驗標準偏差表示)σ=0.40"．該重復性誤差是由轉(zhuǎn)臺本身的重復性、測角系統(tǒng)的重復性誤差(包括自準直儀讀數(shù)重復性、多齒分度臺測角重復性等)以及環(huán)境等因素引起的．通過以上的分析可知，該補償方法對轉(zhuǎn)臺正、反轉(zhuǎn)時都適用，且重復性較好．

圖4 誤差補償效果Fig.4 Effect of error compensation

4 結(jié) 語

本文對轉(zhuǎn)臺分度誤差的來源進行了分析，利用多齒分度臺和雙軸光電自準直儀對轉(zhuǎn)臺的分度誤差進行檢定，根據(jù)諧波分析的結(jié)果和分度誤差的周期性等特點以及數(shù)據(jù)的復雜程度，選擇傅里葉級數(shù)擬合作為轉(zhuǎn)臺分度誤差的補償模型．與其他兩種模型相比，該模型計算簡單且適用于對周期性函數(shù)的描述．對3種模型的擬合效果進行對比，結(jié)果表明，傅里葉級數(shù)具有較高的擬合精度．最后，對轉(zhuǎn)臺分度誤差進行補償實驗，補償后效果很好，大大提高了轉(zhuǎn)臺的分度精度，從而驗證了分度誤差模型的合理性和正確性．

［1］ 郭敬濱，張大廈，張繼承. 大齒輪齒形測量中補償安裝偏心的理論分析[J]. 天津大學學報，2012，42(5)：435-439. Guo Jingbin，Zhang Dasha，Zhang Jicheng. Compensation theory analysis of installation eccentricity in measurement of large gear tooth[J]. Journal of Tianjin University，2012，42(5)：435-439(in Chinese).

［2］ 郭敬濱，王 嫻，劉海軍，等. 擺線齒輪誤差的測量及修形量計算[J]. 天津大學學報，2011，44(1)：85-89. Guo Jingbin，Wang Xian，Liu Haijun，et al. Measurement of cycloidal gear error and calculation of modification[J]. Journal of Tianjin University，2011，44(1)：85-89(in Chinese).

［3］ 楊永才，陳 林，孫銀春，等. 圓光柵刻劃誤差的自動校正[J]. 華東工業(yè)大學學報，1996，18(3)：85-90. Yang Yongcai，Chen Lin，Sun Yinchun，et al. Auto-matic correction for dividing errors of radial grating [J]. East China University of Technology，1996，18(3)：85-90(in Chinese).

［4］ 曹天寧，周鵬飛. 光學零件制造工藝學[M]. 北京：機械工業(yè)出版社，1987. Cao Tianning，Zhou Pengfei. Optical Parts Manufacturing Technology [M]. Beijing：China Machine Press，1987(in Chinese).

［5］ 曹向群，黃維實，金 彤. 光柵計量技術(shù)[M]. 杭州：浙江大學出版社，1992. Cao Xiangqun，Huang Weishi，Jin Tong. Grating Measurement Technology[M]. Hangzhou：Zhejiang University Press，1992(in Chinese).

［6］ 趙人杰，馬文禮. 利用誤差諧波補償法提高金屬圓光柵測角精度[J]. 儀器儀表用戶，2009，16(3)：69-71. Zhao Renjie，Ma Wenli. Improving the accuracy of newtype encoders using error harmonic compensation [J]. Instrumentation Customer，2009，16(3)：69-71(in Chinese).

［7］ 黃宗升，秦石喬，王省書，等. 光柵角編碼器誤差分析及用激光陀螺標校的研究[J]. 儀器儀表學報，2007，28(10)：1866-1869. Huang Zongsheng，Qin Shiqiao，Wang Xingshu，et al. Error analysis of optical angular encoder and its calibration with ring laser gyro [J]. Chinese Journal of Scientific Instrument，2007，28(10)：1866-1869(in Chinese).

［8］ 張禮松，管炳良. 關節(jié)坐標測量機研制中圓光柵誤差修正技術(shù)[J]. 計測技術(shù)，2007，27(4)：41-43，69. Zhang Lisong，Guan Bingliang. Error correction and its application to multi-joint CMM research [J]. Measuring Technique，2007，27(4)：41-43，69(in Chinese).

［9］ Qin Shiqiao，Huang Zongsheng，Wang Xingshu. Optical angular encoder installation error measurement and calibration by ring laser gyroscope [J]. Instrumentation and Measurement，2010，59(3)：506-511.

［10］ Hatiris E，Orton P A，Poliakoff J F，et al. Eccentricity error compensation models for an incremental motion encoder[C]//Proceedings of the 16th European Simulation Multiconference on Modelling and Simulation 2002. Darmstadt，Germany：ESM，2002：328-331.

［11］ 《角度計量》編寫組. 角度計量[M]. 北京：中國標準出版社，1984. Compilation Group of 《Angle Measurement》. Angle Measurement [M]. Beijing：The Standards Press of China，1984(in Chinese).

［12］ Kaul S K，Tickoo A K，Koul R，et al. Improving the accuracy of low-cost resolver-based encoders using harmonic analysis[J]. Nuclear Instruments and Methods in Physics Research Section A：Accelerators，Spectrometers，Detectors and Associated Equipment，2008，586(A)：345-355.

［13］ 李友峰. N次諧波法擬合供油凸輪型線的探討[J]. 內(nèi)燃機車，2006，386(4)：15-17. Li Youfeng. Study for Nth harmonic method fitting fuel cam profile [J]. Diesel Locomotive，2006，386(4)：15-17(in Chinese).

［14］ Mahmoud F. Harmonic analysis and applications [J]. Kybernetes，2012，41(1)：129-144.

［15］ 魏玉業(yè)，趙鳳陽. 曲線最佳擬合的評價標準[J]. 測繪科學，2010，35(1)：195-196. Wei Yuye，Zhao Fengyang. The evaluation criteria of optimal curve fitting [J]. Science of Surveying and Mapping，2010，35(1)：195-196(in Chinese).

［16］ 高貫斌，王 文，林 鏗，等. 圓光柵角度傳感器的誤差補償及參數(shù)辨識[J]. 光學精密工程，2010，18(8)：1766-1772. Gao Guanbin，Wang Wen，Lin Keng，et al. Error compensation and parameter identification of circular grating angle sensors [J]. Optics and Precision Engineering，2010，18(8)：1766-1772(in Chinese).

(責任編輯：趙艷靜)

Check of Indexing Error for Turntable and Establishment of Its Compensation Model

Guo Jingbin1，Cao Hongyan1，Wang Kexin2，Qi Yongyue1
(1. State Key Laboratory of Precision Measuring Technology and Instruments，Tianjin University，Tianjin 300072，China；2. Wenzhou Institute of Measurement Technology，Wenzhou 325000，China)

Error compensation method is an economical and effective technical measure to improve the indexing accuracy for turntable. The indexing error was checked by comparison method，and the measurement uncertainty was also analyzed in this paper. According to the periodicity of indexing error and the harmonic analysis results，and combined with the characteristics of least square method，cubic spline interpolation and Fourier series，this paper chooses Fourier series as the compensation model for turntable indexing error and compares the fitting effects of the three models. By using the Fourier series，the maximum indexing error is 2.36″ after compensation the turntable with 17.82″，which has verified the validity and feasibility of the compensation model.

indexing error；check；data fitting；Fourier series；compensation model

TB922

A

0493-2137(2014)06-0524-06

10.11784/tdxbz201207027

2012-07-09；

2012-10-12.

天津市自然科學基金資助項目(10JCYBJC06600).

郭敬濱（1959— ），男，副教授.

郭敬濱，guojingbin@tju.edu.cn.