褚 寧 張為民,2

(1.同濟大學(xué)機械與能源工程學(xué)院，上海 201804； 2.同濟大學(xué)中德學(xué)院，上海 201804)

基于張量理論的數(shù)控機床誤差補償模型

褚 寧1張為民1,2

(1.同濟大學(xué)機械與能源工程學(xué)院，上海 201804； 2.同濟大學(xué)中德學(xué)院，上海 201804)

為了進一步提高數(shù)控機床多體運動學(xué)誤差模型的精度，提出了基于張量理論的機床誤差補償模型。首先面向數(shù)控系統(tǒng)列表插值補償方式，提出了應(yīng)用張量理論設(shè)定基礎(chǔ)階、擴展階和誤差階來建立機床誤差張量模型。然后，提出了采用設(shè)定機床信息、建立張量模型及設(shè)定檢測參數(shù)3個流程實施誤差補償，并通過多元回歸分析提出了采用4階多項式擬合機床空間誤差。最后，以某重型車銑復(fù)合機床實驗為例，提出了分軸步進對角線法下的張量誤差模型，采用球桿儀進行了圓度誤差檢測實驗，結(jié)果表明該方法能夠有效地提高機床精度，比多體運動學(xué)誤差模型的補償效果更優(yōu)。

誤差補償； 張量理論； 分軸步進對角線法； 車銑復(fù)合機床

引言

數(shù)控機床的精度模型和誤差補償模型是數(shù)控機床研究的熱點之一，國內(nèi)外學(xué)者對誤差建模方法進行了廣泛而深入的研究。DONMEZ等[1]于1986年提出建立包含幾何誤差和熱誤差的機床綜合誤差模型。CHANA等[2]采用多體運動學(xué)理論建立包含受力變形誤差的三軸機床誤差模型，LIN等[3]以及TUTUNEA-FATAN等[4]將多體運動學(xué)擴展到五軸數(shù)控機床的應(yīng)用。SCHWENKE等[5]在CIRP上系統(tǒng)地回顧了誤差數(shù)值補償?shù)幕A(chǔ)以及測量幾何誤差的可用方法。劉又午等[6]將多體理論運用到數(shù)控機床的全誤差建模。文獻[7-8]通過熱誤差分析結(jié)合幾何誤差建立綜合誤差模型及實施誤差補償，并擴展到五軸機床的應(yīng)用[9]。粟時平等[10]提出12線法的誤差辨識方法。李小力[11]`對數(shù)控機床綜合誤差的建模及補償研究進行了系統(tǒng)論述。郭然等[12]采用切比雪夫多項式分解幾何誤差。此外，還有許多學(xué)者分別采用粒子群算法[13]、線性插值算法[14]、人工神經(jīng)網(wǎng)絡(luò)[15]、流變理論[16]、微分構(gòu)造[17]和最優(yōu)權(quán)系數(shù)[18]等建模方法來對機床誤差進行建模。多體運動學(xué)誤差模型[19-20]由于良好的通用性和有效性，仍是最為常用的機床誤差建模方式，隨著現(xiàn)代數(shù)控機床各部件的制造精度越來越高，螺距誤差已不是最大的誤差項，如運動部件的直線度誤差及變形誤差都成為重要的誤差來源，以重型龍門機床為例，懸垂誤差可達400 μm/2 000 mm[21]，基于多體運動學(xué)建模多種誤差來源在公式表達和計算上形式過于復(fù)雜，已逐漸不適應(yīng)誤差補償?shù)男枰?/p>

本文在多體運動學(xué)機床模型的基礎(chǔ)上，參考數(shù)控系統(tǒng)列表插值補償方式，提出應(yīng)用張量理論設(shè)定基礎(chǔ)階、擴展階和誤差階來建立機床誤差張量模型。然后基于激光干涉儀制定機床誤差張量模型檢測與辨識方法。最后以重型車銑復(fù)合數(shù)控機床為實例，驗證數(shù)控機床應(yīng)用張量理論建模與補償?shù)目尚行浴?/p>

1 基于張量理論的數(shù)控機床誤差模型

1.1 數(shù)控系統(tǒng)的誤差補償模型

數(shù)控機床誤差補償功能按照參數(shù)可以劃分為機床定義、機床軸設(shè)定和插值補償表3部分。其中插值補償表采用插值列表形式來實現(xiàn)具體的誤差補償，是數(shù)控系統(tǒng)誤差補償功能底層實施的核心?，F(xiàn)有機床誤差模型多面向擬合精度或運動匹配，很少關(guān)注數(shù)控系統(tǒng)插值補償表，以多體運動學(xué)誤差模型為例，主要和機床結(jié)構(gòu)定義相關(guān)，在機床軸設(shè)定和插值補償表上都只有部分相關(guān)，如圖1所示。

圖1 多體運動學(xué)誤差模型對比數(shù)控系統(tǒng)誤差補償功能Fig.1 Multi-body kinematics model of machine-tool and CNC system error compensation function

由圖1可知，多體運動學(xué)誤差模型能較好地匹配數(shù)控機床誤差補償功能的機床定義，但在機床軸設(shè)定和插值補償表功能方面有所欠缺，無法全面地指導(dǎo)誤差補償實施。

1.2 面向數(shù)控系統(tǒng)的誤差張量建模

張量理論是數(shù)學(xué)的一個分支學(xué)科，是力學(xué)和物理學(xué)一個有力的數(shù)學(xué)工具，一般理解上，張量實為標量、向量、線性算子的母概念，1 階張量即為向量，2 階張量即為矩陣，3 階張量就是三維數(shù)組[22]。

圖2 數(shù)控機床誤差張量模型建模流程圖Fig.2 NC machine-tool error tensor model diagram

以多體運動學(xué)誤差模型為基礎(chǔ)，綜合機床信息和機床插值點誤差兩部分，建立數(shù)控機床空間誤差的誤差張量模型。數(shù)控機床誤差張量模型的建模流程包括設(shè)定機床信息、建立張量模型及設(shè)定檢測參數(shù)3部分，具體建模流程如圖2所示。

由圖2可知，建立張量模型主要包括設(shè)定張量階數(shù)和設(shè)定張量維數(shù)，因此數(shù)控機床誤差張量模型最重要的參數(shù)為張量的階數(shù)、維數(shù)以及檢測參數(shù)，可將張量誤差模型定義為基礎(chǔ)階、擴展階和誤差階3部分。基礎(chǔ)階為定義誤差來源項，前3階為X、Y、Z軸，維數(shù)即為直線軸的插值點數(shù)，如果是五軸機床則包含2個擴展張量階；擴展階為定義變形誤差，主要包括負載變形、切削力變形、熱變形等，維數(shù)為變形的狀態(tài)數(shù)；誤差階為機床TCP點的誤差，誤差階維數(shù)為誤差項數(shù)量，如X、Y和Z軸向定位誤差、階數(shù)及維數(shù)設(shè)定如表1所示。

表1 數(shù)控機床誤差張量模型維數(shù)及階數(shù)定義Tab.1 Dimension and order definition of error tensor model

實際應(yīng)用中數(shù)控機床的誤差張量模型不限于表1所列的11階，一般根據(jù)機床結(jié)構(gòu)和誤差檢測設(shè)備條件進行簡化，以AC五軸機床誤差張量模型為例，全誤差張量模型為十階，即M∈RIiIjIkIoIpIqIuIvIw，其中i、j、k為基礎(chǔ)階的維數(shù)，表征X、Y、Z三軸的插補點數(shù)，o、p為五軸旋轉(zhuǎn)階的維數(shù)，表征A、C軸插補點數(shù)，q、u、v為變形誤差階的維數(shù)，表征負載變形、切削力變形和熱變形的變化狀態(tài)，w為誤差階的維數(shù)，表征機床的定位誤差、正反向誤差等。令各軸插補點數(shù)為n，在不考慮變形誤差，僅計算三軸向空間誤差情況下，可建立簡化的誤差張量模型M0∈RInInInInInI3，如僅計算直線軸則可以建立簡化直線軸誤差張量模型ML∈RInInInI3。

2 面向誤差張量模型的檢測流程及檢測關(guān)鍵參數(shù)

2.1 機床誤差張量模型的檢測流程

數(shù)控機床的誤差檢測是誤差建模與補償?shù)闹匾h(huán)節(jié)，對機床的幾何誤差檢測已有國家標準[23]、國際標準[24]等，但幾何誤差檢測標準的目的在于機床定標，不適于指導(dǎo)精確的誤差補償，無法直接用于誤差補償。激光干涉儀是檢測機床誤差精度最高的設(shè)備之一，本文基于ISO 230檢測標準，采用激光干涉儀來建立面向數(shù)控系統(tǒng)誤差補償功能的機床誤差張量模型，步驟及待測參數(shù)如下：

(1)設(shè)定誤差張量模型階數(shù)：按照圖2的建模方法初步建立機床誤差綜合張量模型，關(guān)鍵參數(shù)為運動鏈和誤差張量模型階數(shù)等。

(2)設(shè)定誤差張量模型維數(shù)及插值點間距：根據(jù)直線軸誤差分布規(guī)律和機床運動行程，設(shè)定各直線軸的插值點間距和測量點數(shù)，關(guān)鍵參數(shù)為誤差張量模型維數(shù)和插值點間距等。

(3)設(shè)定采樣路徑：根據(jù)誤差張量模型維數(shù)和直線軸誤差分布規(guī)律設(shè)定激光干涉儀檢測的定位路徑，關(guān)鍵參數(shù)為采樣路徑和采樣次數(shù)。

(4)檢測定位誤差：參考ISO 230-2的標準檢驗循環(huán)，根據(jù)采樣路徑、測量點數(shù)、插值點間距和采樣次數(shù)，實施定位誤差檢測。

(5)檢測直線度誤差：同定位誤差檢測方法流程相似，但需采用不同的反射鏡組來分別檢測直線度和垂直度誤差。

(6)數(shù)據(jù)處理：基于ISO 230標準，通過計算正向定位誤差，反向定位誤差，正向直線度誤差，反向直線度誤差，正向垂直度誤差和反向垂直度誤差等關(guān)鍵參數(shù)，建立機床誤差綜合張量模型。

2.2 檢測關(guān)鍵參數(shù)的實驗

理論上誤差采樣數(shù)據(jù)量越多張量模型精度越高，但過多采樣數(shù)據(jù)量會使得檢測過程過于繁復(fù)以致無法進行，因此需要通過分析機床插值點誤差的分布特點對插值點間距和采樣次數(shù)進行合理優(yōu)化。

插值點的主要測量誤差項為插值點的正向平均位置偏差，反向平均位置偏差，插值點的平均位置偏差，插值點的反向誤差，單向定位標準不確定度和單向重復(fù)定位精度等，本文以某型車銑復(fù)合機床為實驗對象，如圖3所示，為簡化模型，僅計算直線軸的綜合誤差，直線軸的全工作空間為：x∈[-540,-80] mm，y∈[-500,-40] mm，z∈[-420,-20] mm。

圖3 車銑復(fù)合機床結(jié)構(gòu)簡圖Fig.3 Topological structure of mill-turn machine

實驗采用ML10型激光干涉儀，按照ISO 230-2的檢測標準，每隔20 mm采樣一個點，根據(jù)標準檢測循環(huán)實行檢測，重復(fù)實驗5次，實驗結(jié)果如表2所示，篇幅所限，僅列出部分X軸定位誤差的測量記錄。

插值點間距的設(shè)定和相鄰插值點的X、Y、Z三向定位誤差有關(guān)。一般認為機床誤差為連續(xù)函數(shù)，需要將誤差檢測值通過函數(shù)擬合來定義定位誤差函數(shù)。在數(shù)控系統(tǒng)內(nèi)通常采用多項式擬合來計算誤差值，本文也采用多項式擬合來定義誤差函數(shù)，取10次測量的均值為誤差補償值，以X軸定位誤差為例，為了計算方便，統(tǒng)一以正向檢測誤差為準。

表2 X軸定位誤差測量值Tab.2 X axis positioning error measurement results μm

圖4 X軸定位誤差測量均值的多項式擬合分析Fig.4 Polynomial fitting analysis of measuring mean of X-axis positioning error

圖5 多項式擬合的殘差分析Fig.5 Residual analysis of polynomial fitting

由圖4、圖5可知常用的1階擬合和殘差為-4.20～5.28 μm，超過±4.0 μm的出廠標準；采用2～5階多項式擬合的殘差絕對值的最大值分別為4.93、4.64、3.37、2.64 μm，采用4階多項式擬合即可達到標準；采用2～5階多項式擬合的殘差絕對值的均值分別為1.98、1.96、1.60、1.47 μm，采用5階多項式擬合計算時間增加近50%，但精度提高不足10%，綜合分析采用4階多項式擬合已經(jīng)可以較好地對誤差進行擬合。

插值點間距的設(shè)定由誤差函數(shù)和誤差補償精度和數(shù)控系統(tǒng)功能決定，一般數(shù)控機床插值點之間的誤差為線性插值。以X軸為例，先通過三階多項式擬合實際檢測誤差建立X軸的X、Y、Z三向誤差函數(shù)Px(x)、Py(x)和Pz(x)，即

(1)

式中amn(x)——X軸n向誤差擬合函數(shù)m階項

由式(1)及插值點間距l(xiāng)x和插補段數(shù)i可得數(shù)控系統(tǒng)直線插補所需的插補函數(shù)

(2)

式中iNm(x)——X軸第i段m向誤差插值函數(shù)

計算所得的X軸誤差可以換算為X軸的三向幾何誤差Ex(x)、Ey(x)和Ez(x)，即

(3)

插值點間距越小，誤差補償效果越好，但檢測時間變長且數(shù)控系統(tǒng)負擔更重。取最大的插值點間距的條件為

(4)

式中Lx——X軸插值全長

Nxmax——X軸最大插值點數(shù)

Em(x,Pn(x))——X軸第n點的m向定位誤差

Rm(x,Pn(x))——X軸第n點的m向重復(fù)定位誤差

Epx——機床X軸定位精度要求

Rpx——機床X軸重復(fù)定位精度要求

3 車銑復(fù)合機床實驗驗證

3.1 誤差檢測

實驗采用沈陽機床的某型五軸車銑復(fù)合加工中心[25]，其主要特點在于X、Y、Z構(gòu)成三維空間，C軸為帶分度的車主軸，B軸為繞Y軸旋轉(zhuǎn)的軸，形成刀具與主軸回轉(zhuǎn)中心的夾角，可以換裝銑刀和車刀，總體結(jié)構(gòu)如圖6所示。

圖6 五軸車銑中心結(jié)構(gòu)示意圖Fig.6 Structure diagram of five-axis mill-turn machine

實驗以直線軸為例，誤差檢測實驗選用光動公司的MCV-5005套件，檢測過程采用分軸步進對角線檢測技術(shù)[26]，特點是所測得的位移誤差是平行于運動軸線方向的誤差和垂直于運動軸線方向的誤差的矢量和。實驗裝置及光路設(shè)定如圖7所示。

圖7 分軸步進對角線法(PNP/NPN路徑)Fig.7 Step-diagonal measurement method (PNP/NPN)

設(shè)定分軸步進對角線檢測的采樣路徑為PPP和NNN，然后PNP和NPN，然后NNP和PPN，最后NPP和PNN，重復(fù)檢測3次，測量區(qū)間分別為x=[-400,-10] mm，y=[-500,-20] mm，z=[100,1 000] mm，插值點每個軸為31，每條對角線檢測3遍取平均值，如圖8所示。

圖8 分軸步進對角線檢測結(jié)果Fig.8 Test results of step-diagonal measurement method

3.2 基于分軸步進對角線法的張量誤差模型

數(shù)控機床誤差綜合張量模型包括正向定位誤差，反向定位誤差，正向直線度誤差，反向直線度誤差，正向垂直度誤差和反向垂直度誤差。通過檢測8條對角線后統(tǒng)一計算可得，令PPP路徑上對角線矢量的X、Y、Z分量為a、b、c，任意測量點i上的正向定位誤差計算公式為

(5)

(6)

(7)

式中Em(ni)——n軸第i點的m向正向誤差

εl(ni)k——n軸第i點的k路徑上的激光長度

普通機床以正向誤差加反向間隙補償即可，如數(shù)控系統(tǒng)包含獨立的正反向誤差補償程序(SINUMERIK 840D)，則任意測量點i上的反向定位誤差計算公式為

(8)

(9)

(10)

式中E′m(ni)——n軸第i點的m向反向誤差

將測量數(shù)據(jù)代入式(5)～(10)可得直線軸的誤差簡化張量模型，即建立EL∈RI31×I31×I31×I6。以X軸為例，X軸的誤差簡化張量模型為E(x)=E:11:，參照插值補償如表3所示。

表3 X軸誤差簡化張量模型Tab.3 Error tensor model of X axis mm

840D型數(shù)控系統(tǒng)的誤差補償函數(shù)為高階多項式，理論上選用較高的多項式擬合精度較高，但會造成數(shù)控系統(tǒng)負載過重，魯棒性較差，需要對多項式的階數(shù)進行優(yōu)化。以EL所得的誤差簡化張量模型為測試數(shù)據(jù)，分析多項式擬合的擬合效果，以殘差平方和Esse為對比標準，Esse定義為

(11)

考慮到數(shù)控系統(tǒng)的應(yīng)用實際，對比采用1～6階多項式擬合的擬合效果，分為螺距誤差和直線度誤差，結(jié)果如圖9所示。

圖9 多項式擬合效果Fig.9 Polynomial fitting effect

圖9顯示了多項式擬合機床誤差的擬合效果，螺距誤差的多項式擬合效果較好，3階以上的殘差平方和均小于0.000 03，而直線度的誤差非線性的特征更為明顯，如X軸Z向直線度需要5階以上的擬合才達到0.000 3，說明誤差補償應(yīng)當從只關(guān)注螺距誤差和垂直度誤差轉(zhuǎn)換為更關(guān)注于直線度誤差。

以線性插補的補償效果為基準，高階多項式的擬合效果如圖10所示。

圖10顯示了多項式擬合階數(shù)對機床誤差擬合效果的影響，螺距誤差的多項式擬合效果較好，采用3階即可提高擬合精度59.4%～92.5%，5階可提高64.0%～93.2%；各軸直線度的誤差各不相同，X軸3階提高26.4%～54.3%，而Y軸可提高75.5%～95.5%，Z軸為37.4%～57.3%。普通數(shù)控系統(tǒng)可選用3階多項式擬合即可達到滿意的補償效果，高檔數(shù)控機床則可選擇更高階的補償函數(shù)，進一步提高誤差補償?shù)男Ч?/p>

3.3 基于機床誤差張量模型的誤差補償

令插值點間距l(xiāng)x=20 mm，以3階多項式擬合正向誤差為例，建立X、Y和Z軸的三向誤差函數(shù)

圖10 多項式階數(shù)擬合效果對比Fig.10 Comparisons of polynomial fitting effect by polynomial order

(12)

(13)

(14)

該型號車銑復(fù)合機床采用SINUMERIK 840D型數(shù)控系統(tǒng)，通過擴展840D型數(shù)控系統(tǒng)SAG懸垂誤差補償功能[27]來實現(xiàn)數(shù)控機床誤差張量補償。在840D型系統(tǒng)內(nèi)SAG懸垂誤差補償屬插值補償，提供基準軸和補償軸選擇，如定義基準軸$AN_CEC_INPUT_AXIS[6]=(X)，定義補償軸$AN_CEC_OUTPUT_AXIS[6]=(Y)，則代表第6項SAG誤差補償為相對于X軸Y方向上的誤差，將懸垂誤差以直線度誤差代替，則可實現(xiàn)X軸Y向直線度誤差的補償，其他各軸以此類推。

以球桿儀的圓度誤差檢測實驗驗證機床張量誤差模型，即將誤差張量模型轉(zhuǎn)換為表3所列誤差表，通過西門子840D的SAG補償功能實現(xiàn)誤差補償，以XY軸插補運動做圓度誤差檢測為例，結(jié)果如圖11所示。

分別檢測XY平面、XZ平面和YZ平面的圓度誤差，為避免動態(tài)誤差影響，優(yōu)化動態(tài)參數(shù)，并將進給率設(shè)定為1 000 mm/min，檢測3次取均值，檢測結(jié)果如表4所示。

XY、XZ和YZ平面的原始圓度誤差在7.4～10.3之間，表明機床原始精度較高，在單獨補償各軸螺距誤差后，3個平面的圓度誤差增加12.2%～27.7%，表明空間誤差的補償是三軸誤差的綜合，單獨的螺距誤差無法優(yōu)化多軸的聯(lián)動精度。通過采用多項式擬合的機床誤差張量模型實施補償后，XY平面圓度誤差降低22.3%，XZ平面圓度誤差降低18.9%，XZ平面圓度誤差增加7.2%，由于XY為主要加工面，檢測數(shù)據(jù)較多，考慮到誤差張量模型的補償為綜合優(yōu)化，對XY平面的優(yōu)化效果較好。另一方面，實際誤差的誤差函數(shù)的擬合階數(shù)越高，誤差補償效果越好。圓度誤差的檢測表明，在較高的機床原始精度下單純補償激光干涉儀測量的定位誤差已被證明不僅對圓度誤差沒有改善反而會增加圓度誤差，但采用張量模型實施誤差補償確實提高了機床的插補精度。

4 結(jié)論

(1)提出了數(shù)控機床的誤差張量模型建模方法。通過設(shè)定機床信息、建立張量模型及設(shè)定檢測參數(shù)3個流程建立數(shù)控機床誤差張量模型。張量誤差模型可以定義為基礎(chǔ)階、擴展階和誤差階3部分。

(2)建立了機床誤差張量模型的誤差檢測流程和主要檢測參數(shù)。通過設(shè)定誤差張量模型階數(shù)、設(shè)定誤差張量模型維數(shù)及插值點間距、設(shè)定采樣路徑、檢測定位誤差、檢測直線度誤差和數(shù)據(jù)處理6個步驟建立機床空間誤差張量模型。由實驗數(shù)據(jù)分析可得，定位誤差采用低階多項式擬合的精度不足，采用5階多項式擬合計算時間增加近50%，但精度提高不足10%，采用4階多項式擬合可以較好地對誤差進行擬合。

圖11 誤差補償效果(圓度誤差檢測)Fig.11 Compensation effects by roundness error test

(3)實驗驗證誤差張量模型實施誤差補償?shù)挠行?。建立基于分軸步進對角線法的機床誤差張量模型。螺距誤差的多項式擬合效果較好，3階以上的殘差平方和小于0.000 03，而直線度的誤差需要5階以上的擬合。采用球桿儀進行圓度誤差檢測實驗，在機床原始精度較高(圓度誤差10 μm)條件下，通過誤差張量模型仍能進一步提高20%左右，而通用的螺距誤差補償則會降低多軸聯(lián)動的精度12.2%～27.7%。

表4 球桿儀圓度誤差Tab.4 Roundness error by ballbar test

1 DONMEZ M A, BLOMQUIST D S, HOCKEN R J, et al.A general methodology for machine tool accuracy enhancement by error compensation[J].Precision Engineering, 1986, 8(4): 187-196.

2 CHANA R, MANUKID P.Geometric and force errors compensation in a 3-axis CNC milling machine[J].International Journal of Machine Tools and Manufacture, 2004, 44(12-13):1283-1291.

3 LIN Y,SHEN Y.Modelling of five-axis machine tool metrology models using the matrix summation approach [J].International Journal of Advanced Manufacturing Technology, 2003, 21(4): 243-248.

4 TUTUNEA-FATAN O R,FENG H Y.Configuration analysis of five-axis machine tools using a generic kinematic model [J].International Journal of Machine Tools & Manufacture, 2004, 44(11): 1235-1243.

5 SCHWENKE H,KNAPP W,HAITJEMA H,et al.Geometric error measurement and compensation of machines—an update[J].CIRP Annals-Manufacturing Technology,2008,57(2):660-675.

6 劉又午,章青,趙小松,等.數(shù)控機床全誤差模型和誤差補償技術(shù)的研究[J].制造技術(shù)與機床,2003, 7(4): 46-50.

LIU Youwu,ZHANG Qing,ZHAO Xiaosong,et al.Research on total error model and error compensation of CNC machine tool[J].Manufacturing Technology & Machine Tool,2003, 7(4): 46-50.(in Chinese)

7 沈金華,李永祥,魯志政,等.數(shù)控車床幾何和熱誤差綜合實時補償方法應(yīng)用[J].四川大學(xué)學(xué)報:工程科學(xué)版,2008,40(1):163-166.

SHEN Jinhua,LI Yongxiang,LU Zhizheng,et al.A real-time compensation method for the geometric and thermal error on CNC turning machine[J].Journal of Sichuan University:Engineering Science Edition,2008,40(1):163-166.(in Chinese)

8 王維,楊建國,姚曉棟，等.數(shù)控機床幾何誤差與熱誤差綜合建模及其實時補償[J].機械工程學(xué)報,2012,48(7):165-170,179.

WANG Wei,YANG Jianguo,YAO Xiaodong,et al.Synthesis modeling and real-time compensation of geometric error and thermal error for CNC machine tools[J].Journal of Mechanical Engineering,2012,48(7):165-170,179.(in Chinese)

9 王秀山,楊建國,閆嘉鈺.基于多體系統(tǒng)理論的五軸機床綜合誤差建模技術(shù)[J].上海交通大學(xué)學(xué)報,2008,42(5):761-764,769.

WANG Xiushan,YANG Jianguo,YAN Jiayu.Synthesis error modeling of five-axis machine tools based on multi-body system theory[J].Journal of Shanghai Jiaotong University,2008,42(5):761-764,769.(in Chinese)

10 粟時平,李圣怡,王貴林.基于空間誤差模型的加工中心幾何誤差辨識方法[J].機械工程學(xué)報,2002,38(7):121-125.

SU Shiping,LI Shengyi,WANG Guilin.Identification method for errors of machining center based on volumetric error model[J].Journal of Mechanical Engineering,2002,38(7):121-125.(in Chinese)

11 李小力.數(shù)控機床綜合幾何誤差的建模及補償研究[D].武漢:華中科技大學(xué),2006.

LI Xiaoli.Modelling and compensating for comprehensive geometric error of NC machine tools[D].Wuhan:Huazhong University of Science and Technology,2006.(in Chinese)

12 郭然,付國強,孫磊,等.基于切比雪夫多項式的數(shù)控機床幾何誤差參數(shù)化建模[J/OL].農(nóng)業(yè)機械學(xué)報,2015,46(5):336-343.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20150547&journal_id=jcsam.DOI：10.6041/j.issn.1000-1298.2015.05.047.

GUO Ran,FU Guoqiang,SUN Lei,et al.Chebyshev polynomials based parametric modeling of geometric errors for CNC machine tools[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2015,46(5):336-343.(in Chinese)

13 張成新,高峰,李艷,等.基于分段擬合的機床大尺寸工作臺熱誤差補償模型[J].機械工程學(xué)報,2015,51(3):146-152.

ZHANG Chengxin,GAO Feng,LI Yan,et al.Model of thermal error compensation of large size worktable for machine tools based on piece wise fitting [J].Journal of Mechanical Engineering,2015,51(3):146-152.(in Chinese)

14 王維,楊建國.基于插值算法的數(shù)控機床復(fù)合誤差補償技術(shù)[J].上海交通大學(xué)學(xué)報,2014,48(1):12-15.

WANG Wei,YANG Jianguo.Compound error compensation technique for CNC machine tools based on interpolation method[J].Journal of Shanghai Jiao Tong University,2014,48(1):12-15.(in Chinese)

15 陽紅,方輝,劉立新,等.基于熱誤差神經(jīng)網(wǎng)絡(luò)預(yù)測模型的機床重點熱剛度辨識方法[J].機械工程學(xué)報,2011，47(11):117-124.

YANG Hong，F(xiàn)ANG Hui，LIU Lixin，et al.Method of key thermal stiffness identification on a machine tool based on the thermal error neural network prediction model[J].Journal of Mechanical Engineering,2011,47(11):117-124.(in Chinese)

16 TANG H, DUAN J,LAN S, et al.A new geometric error modeling approach for multi-axis system based on stream of variation theory[J].International Journal of Machine Tools and Manufacture, 2015, 92: 41-51.

17 齊繼寶,楊偉民.基于微分變化構(gòu)造法的數(shù)控機床幾何誤差補償方法[J/OL].農(nóng)業(yè)機械學(xué)報,2016,47(9):398-405.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20160953&journal_id=jcsam.DOI：10.6041/j.issn.1000-1298.2016.09.053.

QI Jibao,YANG Weimin.Differential change construction based geometric error compensation for machine tools[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2016,47(9):398-405.(in Chinese)

18 陽紅,向勝華,劉立新,等.基于最優(yōu)權(quán)系數(shù)組合建模的數(shù)控機床熱誤差在線補償[J/OL].農(nóng)業(yè)機械學(xué)報,2012,43(5):216-221.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20120537&journal_id=jcsam.DOI：10.6041/j.issn.1000-1298.2012.05.037.

YANG Hong,XIANG Shenghua,LIU Lixin,et al.Online compensation for CNC machine thermal error based on optimal weights-based combined modeling[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2012,43(5):216-221.(in Chinese)

19 邢元,張連洪,何柏巖,等.基于多體理論的數(shù)控機床精度逆設(shè)計方法[J/OL].農(nóng)業(yè)機械學(xué)報,2014,45(3):282-287.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20140346&journal_id=jcsam.DOI：10.6041/j.issn.1000-1298.2014.03.046.

XING Yuan,ZHANG Lianhong,HE Baiyan,et al.Precision reverse design of numerical controlled (NC) machine on the basis of multibody theory[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2014,45(3):282-287.(in Chinese)

20 MCKEOWN P A, LOXHAM J.Some aspects of the design of high precision measuring machines[J].Annals of CIRP, 1973, 22(1): 139.

21 WU Fenghe, QIAO Lijun,YAO Lingpan.Deformation compensation of ram components of super-heavy-duty CNC floor type boring and milling machine[J].Chinese Journal of Aeronautics,2012(25):269-275.

22 KOLDA T G, BADER B W.Tensor decompositions and applications[J].SIAM Review,2009, 51(3): 455-500.

23 GB/T 17421.2—2016 機床檢驗通則 第2部分：數(shù)控軸線的定位精度和重復(fù)定位精度的確定[S].

GB/T 17421.2—2016 Test code for machine tools—Part 2:determination of accuracy and repeatability of positioning numerically controlled axes[S].(in Chinese)

24 ISO 230-2-2014 Test code for machine tools—Part 2:determination of accuracy and repeatability of positioning numerically controlled axes[S].

25 李德珍,李憲凱.五軸車銑復(fù)合加工技術(shù)的現(xiàn)狀與發(fā)展趨勢[J].航空制造技術(shù),2009(12):47-50.

LI Dezhen,LI Xiankai.Current status and development trend of five-axis turning/milling complex maching technology[J].Aeronautical Manufacturing Technology,2009(12):47-50.(in Chinese)

26 WANG C.A definition of volumetric error[J].Manufacturing Engineering,2005,134(1): 16-17.

27 SINUMERIK 840D sl Manufacturer/Service[DB/OL].http:∥www.automation.siemens.com.

ErrorCompensationModelofCNCMachine-toolBasedonTensorTheory

CHU Ning1ZHANG Weimin1,2

(1.SchoolofMechanicalEngineering,TongjiUniversity,Shanghai201804,China2.ChineseGermanSchoolforPostgraduateStudies,TongjiUniversity,Shanghai201804,China)

In order to improve the geometric accuracy of high-end CNC machine tools, a machine tool error compensation model based on tensor theory was proposed.Firstly, the error tensor model was established by setting the basic order, the extended order and the error order, which was applicable for error compensation table.Secondly, the error compensation was carried out in three processes, which were setting machine tool information, establishing tensor model and setting detection parameter.Referring to the ISO 230-2 standard, the modeling flow and detection parameters of the error tensor model were presented.Using multiple regression analysis, the linear fitting and residual errors were 2.5 μm, while the four polynomial fit residuals were 1 μm.Finally, taking CNC turn-milling machine experiment as an example, the error tensor model by step-diagonal measurement was proposed.The analysis of experimental data showed that the pitch error by three order fitting can improve the fitting precision by 59.4%～92.5%, the sum of squared residuals was less than 0.000 03; the straightness error needed five order fitting for nonlinear characteristics.The ballbar experiment of roundness error showed that roundness error was increased by 12.2%～27.7%, while only compensated the pitch error.For example, the compensation using error tensor model can be reduced by 22.3% of roundness error atXYplane by three order fitting, which was optimized by 26.4%～95.5% compared with that of the linear interpolation.The experiment showed that the error compensation should be paid more attention to the straightness error.The method of error tensor model can effectively improve the accuracy of machine tools for both pitch error and straightness error.

error compensation; tensor theory; step-diagonal measurement method; CNC turn-milling machine

10.6041/j.issn.1000-1298.2017.10.052

TH161+.2； TH161+.5

A

1000-1298(2017)10-0408-09

2017-06-08

2017-07-28

國家科技重大專項(2012ZX04005-031)

褚寧(1984—)，男，博士生，主要從事機床誤差建模、誤差檢測與誤差優(yōu)化研究，E-mail：chuning0126@163.com

張為民(1965—)，男，教授，博士生導(dǎo)師，主要從事金屬切削機床與加工工藝和制造系統(tǒng)等研究，E-mail：iamt@#edu.cn