曾志江，高貫斌，馬文金

（650500 云南省 昆明市 昆明理工大學 機電工程學院）

0 引言

隨著現(xiàn)代制造水平的不斷提高，對精密儀器的精度要求也越來越高，溫度變化對精密儀器的測量精度影響十分顯著[1]。關(guān)節(jié)臂式坐標測量機是一種新型的多自由度非正交式坐標系測量機，由于其靈活性高和便攜性好，現(xiàn)已逐步應用于航空航天、機械制造、汽車船舶等領(lǐng)域。測量精度是制約關(guān)節(jié)臂式坐標測量機發(fā)展的一個主要因素，溫度變化帶來的熱變形誤差是影響關(guān)節(jié)臂式坐標測量機精度的重要因素[2-3]，采取相應措施減小其影響對提高關(guān)節(jié)臂式坐標測量機的精度具有重要意義。

由于熱誤差補償研究對硬件條件要求高，且需要長時間、大數(shù)據(jù)樣本的測試分析，目前關(guān)節(jié)臂式坐標測量機熱誤差研究的文獻較少。Santolaria[4]等在不同的溫度下對關(guān)節(jié)臂式坐標測量機的運動學參數(shù)進行辨識，以20 ℃下的運動學參數(shù)作為參考基準值，采用四次多項式擬合預測其它溫度下運動學參數(shù)相對于參考值的偏差，并補償?shù)疥P(guān)節(jié)臂式坐標測量機的控制器中，實現(xiàn)熱誤差補償。胡毅[5]通過對一標準桿件的長度測量，分析了在環(huán)境溫度改變及內(nèi)部熱源對測量精度的影響，并建立了基于神經(jīng)網(wǎng)絡(luò)的補償模型對關(guān)節(jié)臂式坐標測量機熱變形誤差進行修正。于連棟[6]等提出了一種新型圓光柵測角誤差補償方法，基于諧波方法建立了含有環(huán)境溫度影響因子的圓光柵測角誤差補償模型，以提高關(guān)節(jié)臂式坐標測量機的測量精度。

綜上所述，目前對關(guān)節(jié)臂式坐標測量機熱誤差研究均是通過實驗的方法獲得溫度變化后的運動學參數(shù)，并研究運動學參數(shù)變化量與溫度的關(guān)系，尚缺乏數(shù)值模擬與理論上的支撐，溫度變化對關(guān)節(jié)臂式坐標測量機測量精度影響的機理還不夠明確。

為探明關(guān)節(jié)臂式坐標測量機熱變形規(guī)律及其對運動學參數(shù)的影響，本文建立了關(guān)節(jié)臂式坐標測量機的有限元分析模型，得到了熱變形分布圖，提出一種虛擬坐標測量方法，在有限元分析后測得關(guān)節(jié)臂式坐標測量機的運動學參數(shù)，并以此建立熱誤差模型，驗證了模型的正確性。

1 關(guān)節(jié)臂式坐標測量機建模與有限元分析

1.1 運動學建模

關(guān)節(jié)臂式坐標測量機是由連桿串聯(lián)多個關(guān)節(jié)組成的一種空間開鏈式機構(gòu)。通常采用D-H（Denavit-Hartenberg）法來建立關(guān)節(jié)臂式坐標測量機的運動學模型[7]。運動學模型實現(xiàn)了關(guān)節(jié)臂式坐標測量機關(guān)節(jié)空間到坐標空間的變換，是關(guān)節(jié)臂式坐標測量機測量的數(shù)學基礎(chǔ)。

本文實驗對象如圖1 所示。關(guān)節(jié)臂的型號為ROMER RA-7125。根據(jù)D-H 法在各關(guān)節(jié)處建立了關(guān)節(jié)臂式坐標測量機的坐標系統(tǒng)，如圖2 所示。通過計算關(guān)節(jié)臂式坐標測量機兩個相鄰關(guān)節(jié)間的齊次變換矩陣Ti-1,i，見式（1），逐個相鄰關(guān)節(jié)變換矩陣連乘后再乘以測頭偏置（lx，ly，lz，1）T，見式（2），就可得到末端測頭中心在基坐標系下的坐標（x,y,z,1）T。

圖1 關(guān)節(jié)臂式坐標測量機Fig.1 Articulated arm coordinate measuring machine

圖2 關(guān)節(jié)臂式坐標測量機的坐標系統(tǒng)圖Fig.2 Coordinate systems of articulated arm coordinate measuring machine

式中：i=1～6——關(guān)節(jié)序號；θi——關(guān)節(jié)轉(zhuǎn)角，為變量，其值可由安裝在軸端的編碼器測得；αi——關(guān)節(jié)扭轉(zhuǎn)角；di——連桿偏距；ai——桿長。

式中：lx，ly，lz——測頭中心在第6 關(guān)節(jié)坐標系中的位置。

1.2 關(guān)節(jié)臂式坐標測量機熱誤差有限元分析

關(guān)節(jié)臂式坐標測量機的工作溫度范圍通常為10～40 ℃，本文以25 ℃為參考基準，在10～40 ℃范圍內(nèi)分析關(guān)節(jié)臂式坐標測量機各關(guān)節(jié)運動學參數(shù)隨環(huán)境溫度變化的規(guī)律。

相比實際關(guān)節(jié)臂式坐標測量機，對外置平衡機構(gòu)、外殼、防護罩部分進行簡化，建立簡化的三維模型。這些零件對運動學參數(shù)的熱變形沒有影響，簡化的模型不影響熱誤差分析結(jié)果。關(guān)節(jié)臂式坐標測量機主要由2 種材料組成，關(guān)節(jié)處的零件主要為7075 鋁合金，連桿為碳纖維，相關(guān)屬性如表1 所示。關(guān)節(jié)臂式坐標測量機的運動學參數(shù)如表2 所示。

表1 關(guān)節(jié)臂式坐標測量機主要材料屬性Tab.1 Properties of main materials of articulated arm coordinate measuring machine

表2 運動學參數(shù)值Tab.2 Values of kinematic parameters

將關(guān)節(jié)臂式坐標測量機簡化三維模型導入有限元分析軟件，對模型的各部分材料進行設(shè)置。ANSYS 自動網(wǎng)格劃分的網(wǎng)格自適應能力較強，本文采用其自動網(wǎng)格劃分，網(wǎng)格大小也采用智能設(shè)置。固定關(guān)節(jié)臂式坐標測量機基座，并施加溫度載荷，分析施加溫度載荷后關(guān)節(jié)臂式坐標測量機的整體形變。

本文在3 種位姿下分析關(guān)節(jié)臂式坐標測量機熱變形情況。3 種位姿分別采取較近、適中、較遠位置，以此覆蓋關(guān)節(jié)臂式坐標測量機的主要工作空間。在關(guān)節(jié)臂式坐標測量機的工作溫度范圍10～40 ℃下每隔3 ℃采集一次數(shù)據(jù)做分析，由于40 ℃下關(guān)節(jié)臂式坐標測量機的熱變形誤差最大，所以選用40 ℃下3 個位姿的熱變形圖作為例子展示其熱變形情況，分別如圖3—圖5 所示。

圖3 40 ℃下位姿一熱變形Fig.3 Thermal deformation of pose 1 at 40 ℃

圖4 40 ℃下位姿二熱變形Fig.4 Thermal deformation of pose 2 at 40 ℃

圖5 40 ℃下位姿三熱變形Fig.5 Thermal deformation of pose 3 at 40 ℃

由圖3—圖5 可以看到，隨著末端位置遠離基座，熱變形越大，熱變形隨著機構(gòu)呈現(xiàn)逐級放大的情況[8-9]，在較遠位置的最大變形量為0.213 mm。

2 熱誤差補償建模與仿真驗證

2.1 關(guān)節(jié)空間參數(shù)虛擬測量及熱誤差補償模型建立

為量化熱誤差對關(guān)節(jié)臂式坐標測量機參數(shù)的影響，提出一種虛擬坐標測量方法?；谧鴺藴y量原理提取有限元分析后的節(jié)點位置數(shù)據(jù)，利用點擬合線和面等幾何元素，再通過對擬合的幾何元素之間距離和夾角的計算，得到關(guān)節(jié)臂式坐標測量機的運動學參數(shù)。如圖6 所示，通過虛擬測量，發(fā)現(xiàn)關(guān)節(jié)臂式坐標測量機的角度參數(shù)θi，αi不變，而長度參數(shù)di，ai變化較明顯。

圖6 有限元節(jié)點提取及幾何元素擬合Fig.6 Extraction of finite element nodes and geometric element fitting

由于溫度變化對長度為0 的幾何參數(shù)不會產(chǎn)生影響，所以本文只對關(guān)節(jié)臂式坐標測量機中不為0 的長度量d1，d3，d5，d6，a2，a4，a6進行測量，分析溫度變化對它們的影響規(guī)律。經(jīng)分析發(fā)現(xiàn)，關(guān)節(jié)臂式坐標測量機3 個位姿在同一溫度下的各參數(shù)的測量值基本一致，取3 個位姿測量的平均值。在不同溫度下測量得到的結(jié)果如表3 所示。可以看出，隨著溫度的升高，長度參數(shù)總體呈逐漸增加的趨勢。

表3 各溫度下長度參數(shù)Tab.3 Length parameters at different temperatures

根據(jù)表3 中長度參數(shù)值的變化可發(fā)現(xiàn)，長度參數(shù)隨溫度升高呈近似線性增加的現(xiàn)象。由材料平均線膨脹系數(shù)計算式[10]對關(guān)節(jié)臂式坐標測量機6 個關(guān)節(jié)各用一個對應的溫度系數(shù)Ci來描述長度參數(shù)隨溫度變化的關(guān)系，如式（3）所示。根據(jù)虛擬測量的結(jié)果計算得到各關(guān)節(jié)的溫度系數(shù)如表4 所示。通過各關(guān)節(jié)溫度系數(shù)可以計算各關(guān)節(jié)在工作溫度范圍內(nèi)任意溫度下的長度參數(shù)。

表4 各關(guān)節(jié)溫度系數(shù)Tab.4 Temperature coefficients of the joints

式中：L0——基準溫度下的長度；ΔL ——長度的變化量；ΔT——溫度的變化量；i=1～6——關(guān)節(jié)1～6。

由有限元分析后得到溫度系數(shù)Ci可計算關(guān)節(jié)臂式坐標測量機各關(guān)節(jié)在工作溫度范圍內(nèi)任意溫度下的長度參數(shù)變化量，從而建立基于運動學參數(shù)的關(guān)節(jié)臂式坐標測量機熱誤差補償模型，如式（4）和式（5）。

式中：i=1～6 ——關(guān)節(jié)1～6。

2.2 熱誤差仿真補償

使用式（5）中給出的誤差補償模型，對關(guān)節(jié)臂式坐標測量機熱誤差進行補償仿真。為了避免使用建立誤差補償模型時用過的溫度數(shù)據(jù)，分別在15，20，30，35 ℃下通過虛擬測量的方法采集數(shù)據(jù)并進行驗證。將25 ℃作為參考溫度，在此溫度下位姿一、二、三的測頭中心坐標分別為P1（-44.026，-705.601，448.854），P2（105.438，-1 186.060，230.614），P3（252.507，-812.426，356.097）。六自由度關(guān)節(jié)臂式坐標測量機是冗余機構(gòu)，測頭中心保持在同一個位置時，可以通過改變各關(guān)節(jié)的轉(zhuǎn)角來改變其姿態(tài)。為了將姿態(tài)對測量結(jié)果的影響納入進來，在上述P1，P2，P3位置處，通過運動學逆解計算出每個位置的任意5 組關(guān)節(jié)轉(zhuǎn)角，將這些關(guān)節(jié)角代入式（5），可以分別計算出熱誤差補償前后15，20，30，35℃下的末端位置。對每個位置的5 組末端位置數(shù)據(jù)分別求其平均值作為該溫度下補償后的位置，3 個補償后的位置再計算相互之間的空間距離，通過在15，20，30，35 ℃下的仿真距離，分析與補償后的距離差值驗證空間距離精度提高的情況。圖7—圖9 分別為距離1、距離2、距離3在上述驗證組溫度下補償前后空間距離精度對比圖，熱誤差補償后3個空間距離精度分別從0.025，0.024，0.035 mm 提高到了0.008，0.012，0.014 mm，空間距離精度得到顯著提高。

圖7 距離1 補償前后精度對比Fig.7 Accuracy comparison before and after compensation for distance 1

圖8 距離2 補償前后精度對比Fig.8 Accuracy comparison before and after compensation for distance 2

3 結(jié)論

對關(guān)節(jié)臂式坐標測量機的熱誤差進行有限元分析，得到了熱變形分布圖；提出一種虛擬坐標測量方法測得熱變形后的運動學參數(shù)隨溫度的變化規(guī)律，對每個關(guān)節(jié)提出一個溫度系數(shù)并建立相應的熱誤差補償模型，使關(guān)節(jié)臂式坐標測量機的精度及穩(wěn)定性得到了極大的提高。

采用有限元分析手段及數(shù)值仿真技術(shù)對關(guān)節(jié)臂式坐標測量機熱誤差進行分析，彌補了目前缺乏數(shù)值模擬與理論支撐的缺陷；通過運動學參數(shù)的虛擬測量進一步明確了熱誤差對關(guān)節(jié)臂式坐標測量機精度的影響機理，對關(guān)節(jié)臂式坐標測量機后續(xù)的熱誤差分析提供了有效指導。