王洪樂 周青華 熊青春③

(①西安增材制造國家研究院有限公司，陜西 西安 710300；② 四川大學(xué)空天科學(xué)與工程學(xué)院，四川 成都 610065；③成都飛機工業(yè)(集團)有限公司，四川 成都 610092)

精密超精密數(shù)控加工機床精度預(yù)測和誤差補償問題一直是學(xué)術(shù)界研究的熱點[1]。伴隨著國家自然科學(xué)基金委員會《機械工程學(xué)科發(fā)展戰(zhàn)略報告(2011-2020)》的出版[2]和“高檔數(shù)控機床與基礎(chǔ)制造裝備”重大專項的開展，國內(nèi)外研究學(xué)者對多軸數(shù)控機床加工精度預(yù)測和誤差補償?shù)葐栴}進行了大量深入的科學(xué)研究。

為了減小數(shù)控機床的機械加工誤差，提高被加工工件精度，必須對整個機械加工工藝系統(tǒng)誤差影響因素進行測量、評定、建模和補償，因此建立準確可靠的數(shù)控機床精度誤差預(yù)測模型，為后期數(shù)控機床誤差實時動態(tài)補償提供有效補償參量極其重要。國內(nèi)外研究學(xué)者針對三軸、五軸機床以及加工中心，采用了幾何誤差建模法、誤差參數(shù)建模法、二次關(guān)系誤差建模法、動靜態(tài)誤差時變特性建模法、機器人機構(gòu)學(xué)建模法等方法進行了數(shù)控機床精度誤差的建模與誤差補償。Han[3]在1986年采用傅里葉變換方法建立了數(shù)控機床的空間幾何誤差模型。Reshetov和Portman[4]在1988年根據(jù)小角度誤差假設(shè)，采用變分法計算出任意結(jié)構(gòu)的數(shù)控機床誤差精度模型，將數(shù)控機床的誤差參數(shù)值采用函數(shù)形式表達，該精度誤差模型包含了數(shù)控機床成形系統(tǒng)誤差模型、物體小位移誤差模型和空間幾何誤差模型。粟時平[5]在2002年運用多體系統(tǒng)運動學(xué)理論，分析了多軸數(shù)控機床的運動副間的拓撲結(jié)構(gòu)，將三軸、五軸數(shù)控機床各部件的關(guān)聯(lián)關(guān)系用低序體陣列來表示，并將各部件的幾何變換通過四階齊次特征矩陣表示，系統(tǒng)完整地推導(dǎo)出了數(shù)控機床運動誤差精度模型。Hsu[6]在2007年通過將數(shù)控機床各個部件作為剛體，根據(jù)多體系運動學(xué)理論和四階齊次坐標變換矩陣構(gòu)建了RRTTT型的五軸數(shù)控機床幾何誤差模型。2012年范晉偉[7]研究了帶擺角頭的數(shù)控機床幾何誤差的建模與補償方法，主要建立了理想條件與實際條件下刀具路線數(shù)控指令及刀具軌跡三者間的相互映射關(guān)系。2016年李杰[8]根據(jù)多軸數(shù)控機床誤差元素的時變特性，提出了一種數(shù)控機床靜、動態(tài)綜合誤差的數(shù)控機床誤差建模方法，并基于此構(gòu)建了靜、動態(tài)誤差與刀具軌跡誤差的映射關(guān)系模型。

綜合大量學(xué)術(shù)資料發(fā)現(xiàn)，針對數(shù)控機床誤差檢測、精度誤差建模、誤差補償方法等方面仍然存在著理論通用性差、誤差元素表述方式困難、計算參數(shù)易受人為因素影響而出現(xiàn)誤差等方面的問題。本文主要從工程實踐角度出發(fā)，基于天津大學(xué)劉又午[9-10]教授提出的多體系統(tǒng)理論，研究了B-A擺頭五軸龍門數(shù)控機床的誤差產(chǎn)生機理和誤差精度建模方法步驟，并建立了B-A擺頭五軸龍門數(shù)控機床精度誤差預(yù)測模型。通過采用雙頻激光干涉儀測量B-A擺頭五軸數(shù)控機床的X/Y/Z軸的各項誤差，并運用國防科技大學(xué)李圣怡團隊[11]提出的平動軸“十二線法辨識算法”進行計算，得出數(shù)控機床的21項空間幾何誤差。將上述各項誤差值代入到B-A擺頭五軸龍門數(shù)控機床的誤差精度預(yù)測模型中，預(yù)測數(shù)控機床的加工精度，并與三坐標測量儀檢測結(jié)果對比，為后續(xù)開發(fā)數(shù)控機床誤差補償分析計算系統(tǒng)和切削工藝優(yōu)化提供理論指導(dǎo)。

1 B-A擺頭五軸龍門數(shù)控機床精度誤差建模

本文主要研究對象為B-A擺頭的五軸龍門數(shù)控機床，如圖1所示。在結(jié)構(gòu)方面，相對于普通三軸機床，增加了2個擺動軸，使其主軸刀具具有姿態(tài)調(diào)整功能，為復(fù)雜曲面結(jié)構(gòu)件的切削加工提供了便利。然而，在多數(shù)機械加工工廠都普遍存在由于在機床服役期過長、環(huán)境熱影響、運動副磨損、維修保養(yǎng)不到位等因素，使同批大量切削加工工件存在尺寸超差和表面質(zhì)量不達標的誤差現(xiàn)象。從該機床主軸端刀具角度看，其刀具位置和姿態(tài)的誤差包括了機床幾何運動誤差、切削熱影響誤差、切削力誤差、伺服控制誤差等等，并把主軸端刀具誤差，稱之為B-A擺頭五軸龍門數(shù)控機床綜合誤差。因此，從B-A擺頭五軸龍門數(shù)控機床主軸端刀具出發(fā)，建立該機床精度誤差預(yù)測模型，測量并辨識各項空間幾何誤差參數(shù)，可以為后期數(shù)控機床誤差實時動態(tài)補償技術(shù)提供理論依據(jù)。

1.1 多體系統(tǒng)理論

多體系統(tǒng)理論[10]是天津大學(xué)劉又午教授在70年代后期提出的，用于分析機械臂、三坐標測量機、數(shù)控機床等開環(huán)或閉環(huán)多體系統(tǒng)裝備的運動學(xué)和動力學(xué)關(guān)系的理論。其主要采用低序體陣列矩陣方法描述運動部件間的關(guān)系，簡潔方便且適用于計算機自動編程描述。多體系統(tǒng)是對一般復(fù)雜機械系統(tǒng)的完整抽象和有效描述方法，可以用若干柔性或剛性體通過某種方式連接而成，在實際應(yīng)用中的機械系統(tǒng)都可以抽象成多體系統(tǒng)，目前已經(jīng)在數(shù)控機床、機器人坐標測量機等復(fù)雜機械系統(tǒng)中廣泛應(yīng)用。如圖1所示，首先在各典型“體”上建立參考坐標系；以大地為慣性參考系體B0，以此為基礎(chǔ)，在與大地相連的B-A擺頭五軸龍門數(shù)控機床工作臺體1上選定B1體；然后可以按照遠離體1距離的方向，依次劃定工件2上選定B2體、X向?qū)к夁\動軸3上選定體B3、在Y向橫梁運動軸4上選定體B4、在主軸方向的Z向運動軸5上選定體B5、在A偏擺軸6選定體B6、在B偏擺軸7選定體B7；依據(jù)以上分析，建立B-A擺頭五軸龍門數(shù)控機床多體系統(tǒng)拓撲結(jié)構(gòu)關(guān)系圖，如圖2所示。

表1 低序體陣列表

分析該B-A擺頭五軸龍門數(shù)控機床的低序體陣列關(guān)系。從圖2可以清晰發(fā)現(xiàn)，各個典型體之間都是通過單約束關(guān)系聯(lián)系起來的，如工作臺1和X向運動軸3之間是移動副，Z向主軸運動體5與刀具A向偏擺軸之間是轉(zhuǎn)動副。體B2的1階低序體是B1，體B5的2階低序體是B3；將各個體之間的低序體陣列關(guān)系匯總?cè)绫?所示。

在典型體Bj和相鄰低序體Bi上分別在其上建立廣義坐標系Oj-XjYjZj和Oi-XiYiZi，且坐標系Bj和Bi是動態(tài)坐標系，在慣性參考系B0上建立靜坐標系O0-X0Y0Z0。點Oj相對于點Oi的位置變化表征了典型體Bj和相鄰低序體Bi的平移位置變化；坐標系Oj-XjYjZj相對于坐標系Oi-XiYiZi的旋轉(zhuǎn)變化表征了典型體Bj和相鄰低序體Bi的旋轉(zhuǎn)姿態(tài)變化。因此，坐標系Oj-XjYjZj和Oi-XiYiZi的相對變化，就轉(zhuǎn)換為了平移位置變化(Wijp)和旋轉(zhuǎn)姿態(tài)變化(Wijs)的結(jié)合，此變化關(guān)系可以用式(1)表示。

Bj←Wijs·Wijp·Bi

(1)

對于B-A擺頭五軸龍門數(shù)控機床的各個運動部件來說，各個體之間的運動關(guān)系由運動約束和結(jié)構(gòu)約束進行限定，從每個運動部件都有6個自由度角度考慮，結(jié)合圖1給出的B-A擺頭五軸龍門數(shù)控機床各部件關(guān)系，可以分析得到各個運動體之間的約束自由度，匯總?cè)绫?所示。

表2 B-A擺頭五軸龍門數(shù)控機床各運動部件的自由度

在研究了B-A擺頭五軸龍門數(shù)控機床的拓撲結(jié)構(gòu)、各個典型體之間的低序體陣列、各個典型體之間的自由度約束基礎(chǔ)之上，就可以對各個典型體之間的變換進行詳細分析。各體間的變換矩陣采用4階齊次矩陣。B-A擺頭五軸龍門數(shù)控機床的工作臺1到工件2的轉(zhuǎn)換矩陣為：

(2)

由于B-A擺頭五軸龍門數(shù)控機床在切削加工中工件2固定在工作臺1上，式(2)中不包含元素誤差量。

工作臺1到X移動導(dǎo)軌3的轉(zhuǎn)換矩陣為：

(3)

當X軸運動x時，在其6個自由度方向產(chǎn)生誤差，在式(3)中用Δx13、Δy13、Δz13、Δα13、Δβ13和Δγ13表示。以下相同。

由X向運動的立柱到Y(jié)向運動橫梁4的轉(zhuǎn)換矩陣為：

(4)

由Y向運動橫梁到Z向運動軸5的轉(zhuǎn)換矩陣為：

(5)

由Z向運動軸到A偏擺軸6的轉(zhuǎn)換矩陣為：

(6)

由A偏擺軸到B偏擺軸7的轉(zhuǎn)換矩陣為：

(7)

B偏擺軸到刀具8的轉(zhuǎn)換矩陣為：

(8)

在切削加工中，刀具8安裝在主軸端精度極高，在精度誤差模型中暫不考慮其誤差。

1.2 數(shù)控機床幾何誤差精度預(yù)測模型

以B-A擺頭五軸龍門數(shù)控機床為研究對象，根據(jù)前文敘述的方法，結(jié)合主軸端刀具在無誤差理想和有誤差實際條件下的成形運動函數(shù)、成形運動約束方程及實際成形運動關(guān)系，可以建立得出B-A擺頭五軸龍門數(shù)控機床的幾何精度誤差預(yù)測模型。設(shè)定主軸端刀具上任一點P(px,py,pz)，其在刀具坐標系內(nèi)的齊次坐標為P(px,py,pz,1)，則在工作臺坐標系中的表示為：

(9)

工件上被加工點T(tx,ty,tz) ，其在工件自身坐標系內(nèi)的齊次坐標為T(tx,ty,tz,1)，則在工作臺坐標系中的表示為：

(10)

則B-A擺頭五軸龍門數(shù)控機床的精度誤差預(yù)測模型為：

(11)

2 數(shù)控機床誤差測量

影響數(shù)控機床出現(xiàn)切削加工誤差的因素眾多，包括傳動部件磨損誤差、熱環(huán)境影響誤差、切削速度、切削加速度沖擊誤差等等，結(jié)合第2章節(jié)中的精度誤差模型中的誤差參數(shù)，從主軸端刀具誤差可以稱之為B-A擺頭五軸龍門數(shù)控機床的綜合誤差這個角度考慮，采用雙頻激光干涉儀測量B-A擺頭五軸龍門數(shù)控機床空間幾何誤差，并將測量結(jié)果運用粟時平[5，11]提出的“十二線法辨識算法”計算，得出B-A擺頭五軸龍門數(shù)控機床空間21項幾何誤差。

雙頻激光干涉儀測量五軸龍門數(shù)控機床的平動軸誤差。首先在機床工作臺上建立空間虛擬立方體，如圖3所。

在B-A擺頭五軸龍門數(shù)控機床工作臺上設(shè)定沿數(shù)控機床工作空間中的虛擬立方體12條直線運動測量，其包括X、Y、Z軸單獨運動和聯(lián)動，誤差測量結(jié)果如圖4所示。

本文將“十二線誤差辨識算法”[11]應(yīng)用于測量到的B-A擺頭五軸龍門數(shù)控機床的誤差測量數(shù)據(jù)，辨識得到B-A擺頭五軸龍門數(shù)控機床的X/Y/Z運動軸的21項空間幾何誤差，如圖5所示，并將該結(jié)果應(yīng)用于前文的B-A擺頭五軸龍門數(shù)控機床的精度誤差預(yù)測模型式(11)中。

3 數(shù)控機床精度誤差預(yù)測

為了驗證該B-A擺角五軸龍門數(shù)控機床的精度誤差預(yù)測模型的可行性，以及采用平動軸的“十二線辨識算法”[11]的有效性，進行了B-A擺頭五軸數(shù)控機床切削工件精度測量實驗。被切削工件尺寸為：400 mm×400 mm×200 mm，銑削其上表面。經(jīng)計算得到尺寸精度最大誤差為23 μm至-10 μm，如圖6所示。采用三坐標測量機，檢測工件表面上相應(yīng)點的誤差值并繪制在圖6中?？梢钥闯?，精度誤差預(yù)測與測量結(jié)果吻合度很好。表3中給出了在該工件表面上選定的6個測量點的對比結(jié)果，充分表明了該B-A擺角五軸龍門數(shù)控機床精度誤差預(yù)測模型準確度高達34%以上，對于誤差較大點位的準確度更高，為后期該B-A擺頭五軸龍門數(shù)控機床的誤差實時動態(tài)補償提供了非常有效的參考價值和理論指導(dǎo)依據(jù)。

4 結(jié)語

本文采用多體系統(tǒng)運動學(xué)理論，建立了B-A擺頭五軸龍門數(shù)控機床誤差精度預(yù)測模型，并通過對該機床平動軸進行的誤差參數(shù)測量實驗，驗證了本文所建立的預(yù)測模型的準確性，對于后期該B-A擺頭五軸龍門數(shù)控機床精度誤差實時動態(tài)補償提供了極其有效的參考數(shù)據(jù)。采用雙頻激光干涉儀測量該B-A擺頭五軸龍門數(shù)控機床在工作臺空間上標定立方體的12條運動直線產(chǎn)生的誤差，并采用“十二線法誤差辨識算法”計算得到了該B-A擺頭五軸龍門數(shù)控機床的21項空間幾何誤差，將誤差參數(shù)辨識結(jié)果應(yīng)用于該B-A擺頭五軸龍門數(shù)控機床精度誤差預(yù)測模型，預(yù)測結(jié)果和三坐標測量儀測量結(jié)果十分接近。該方法通用性較好，適用于各種類型結(jié)構(gòu)的多軸數(shù)控機床，為開發(fā)B-A擺頭五軸龍門數(shù)控機床誤差實時動態(tài)補償軟件提供了極好的理論基礎(chǔ)，并為后期工程實踐中的B-A五軸龍門數(shù)控機床誤差補償提供了可靠依據(jù)。

表3 誤差精度預(yù)測和測量

[1]周志雄, 周秦源, 任瑩輝. 復(fù)雜曲面加工技術(shù)的研究現(xiàn)狀與發(fā)展趨勢[J]. 機械工程學(xué)報, 2010, 46(17): 105-113.

[2]國家自然科學(xué)基金委員會工程與材料科學(xué)部. 機械工程學(xué)科發(fā)展戰(zhàn)略報告(2011—2020)[M].北京：科學(xué)出版社，2010.

[3]Han Z J, Zhou K. Improvement of positioning accuracy of rotating table microcomputer controlcompensation[J]. Precision Engineering, 1986, 4(8): 115-120.

[4]Reshetov D N, Portman V T. Accuracy of machine tools[M]. ASME Press, New York, 1988.

[5]粟時平. 多軸數(shù)控機床精度建模與誤差補償方法研究[D]. 長沙: 國防科技大學(xué)， 2002.

[6]Hsu YY, Wang S S. A new compensation method for geometry errors of five-axis machine tools[J]. International of Machine Tools and Manufacture, 2007,47(2): 352-360.

[7]范晉偉, 羅建平, 蒙順政, 等. 帶擺角頭五軸數(shù)控機床幾何誤差建模及補償方法研究[J]. 機械設(shè)計與制造, 2012(11): 4-6.

[8]李杰, 劉辛軍, 謝福貴, 等. 基于時變特性的數(shù)控機床綜合誤差建模方法[J]. 科技導(dǎo)報, 2016, 34(2): 65- 70.

[9]劉又午. 多體動力學(xué)的休斯頓方法及其發(fā)展[J]. 中國機械工程, 2000, 11(6): 601-607.

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

[11]李圣抬, 戴一凡, 尹自強，等. 精密和超精密機床精度建模技術(shù)[M]. 長沙: 國防科技大學(xué)出版社, 2007.