齊立哲

(北京航空航天大學(xué)機(jī)器人研究所，北京100191)

王 偉

(北京航空航天大學(xué)機(jī)器人研究所，北京100191)

陳 磊

(上海飛機(jī)制造有限公司，上海200436)

余蕾斌

(上海飛機(jī)制造有限公司，上海200436)

工業(yè)機(jī)器人重復(fù)定位精度很高，但絕對定位精度很差.為了提高機(jī)器人的性能及拓展工業(yè)機(jī)器人的應(yīng)用范圍，需要對機(jī)器人運(yùn)動學(xué)模型的參數(shù)進(jìn)行標(biāo)定來降低它的絕對定位誤差.標(biāo)定是建模、測量、參數(shù)識別和誤差補(bǔ)償4個步驟的集成過程［1］.目前已有很多國內(nèi)外學(xué)者基于矩陣變換、矢量等表達(dá)形式給出了多種機(jī)器人的定位誤差分析與補(bǔ)償方法［2-6］，但由于影響機(jī)器人定位誤差的因素很多，基于理論分析公式對機(jī)器人定位誤差進(jìn)行實際標(biāo)定的結(jié)果都不太理想.為了有效實際標(biāo)定出機(jī)器人的定位誤差，Ting K等人提出了從參數(shù)誤差到末端誤差映射的靈敏性雅可比概念，證明了要提高機(jī)構(gòu)的末端誤差，只需要提高關(guān)鍵因素的精度，可以適當(dāng)放寬次要因素誤差［7］;郭劍鷹等人提出將影響轉(zhuǎn)動誤差和位置誤差的因素分離開發(fā)，采用兩步線性識別方法對參數(shù)進(jìn)行標(biāo)定的方法，在一定程度上簡化了標(biāo)定過程［8］;王偉等人采用正交試驗法綜合了3R機(jī)器人機(jī)構(gòu)精度，繪制出了某一復(fù)雜曲面任務(wù)下機(jī)器人末端的誤差分布圖和影響末端誤差因素的顯著程度順序圖，計算出了當(dāng)前目標(biāo)點下誤差因素對末端誤差的影響顯著程度［9］.

本文針對因影響定位誤差因素較多而不易對機(jī)器人定位誤差進(jìn)行準(zhǔn)確標(biāo)定的問題，充分考慮了定位誤差測量環(huán)節(jié)中對定位誤差有影響的坐標(biāo)系參數(shù)因素及因機(jī)器人末端點分布不同而導(dǎo)致定位誤差也不一樣的情況，首先基于雅克比矩陣建立了包含33個運(yùn)動學(xué)參數(shù)的6自由度串聯(lián)機(jī)器人全閉環(huán)定位誤差數(shù)學(xué)模型;然后給出了用于分析各因素單點敏感度的機(jī)器人的單因素微分定位誤差(SFDPE，Single Factor Differential Positioning Error)及相對單因素微分定位誤差(RSFDPE，Relative Single Factor Differential Positioning Error)公式;接著采用正交實驗法設(shè)計了能全面反映機(jī)器人定位誤差分布的數(shù)據(jù)采樣空間并基于區(qū)間概率密度概念給出了單因素綜合微分定位誤差(SFCDPE，Single Factor Comprehensive Differential Positioning Error)及相對單因素綜合微分定位誤差(RSFCDPE，Relative Single Factor Comprehensive Differential Positioning Error)公式;最后基于前面的公式，采用MATLAB分別仿真了影響由FARO ARM和ABB 1410工業(yè)機(jī)器人形成的測量系統(tǒng)定位誤差的各因素的單點敏感度(SPS，Single Point Sensitivity)及多點綜合敏感度(MPCS，More Points Comprehensive Sensitivity)，為進(jìn)一步進(jìn)行機(jī)器人定位誤差的標(biāo)定打下了基礎(chǔ).

1 機(jī)器人全閉環(huán)定位誤差數(shù)學(xué)模型

1.1 全閉環(huán)運(yùn)動學(xué)模型

為了測量出機(jī)器人定位誤差，需要建立機(jī)器人定位誤差測量系統(tǒng).如圖1所示，一般機(jī)器人定位誤差測量系統(tǒng)由測量設(shè)備、測量靶標(biāo)及機(jī)器人組成.

圖1 機(jī)器人定位誤差測量系統(tǒng)組成原理圖

在機(jī)器人定位誤差測量系統(tǒng)中，測量設(shè)備可以直接測量出靶標(biāo)中心P在Measure坐標(biāo)系下的坐標(biāo)，同時機(jī)器人本身也可以作為測量裝置獲得被機(jī)器人抓著的靶標(biāo)中心P在機(jī)器人Base坐標(biāo)系下的坐標(biāo)，兩者進(jìn)行坐標(biāo)變換便可以統(tǒng)一到同一個坐標(biāo)系下，進(jìn)而進(jìn)行比較便可獲得機(jī)器人在其工作空間內(nèi)各個測量點的定位誤差.基于機(jī)器人的定位誤差測量過程，便形成了一個包含2個轉(zhuǎn)換矩陣的全閉環(huán)運(yùn)動學(xué)模型，其表示如下:

其中，PM為測量點P在Measure坐標(biāo)系下的坐標(biāo);BTM為機(jī)器人Base坐標(biāo)系與Measure坐標(biāo)系下的轉(zhuǎn)換矩陣;T0TB為工業(yè)機(jī)器人Tool0坐標(biāo)系在機(jī)器人Base坐標(biāo)系下的轉(zhuǎn)換矩陣;PT0為P點在Tool0坐標(biāo)系下的坐標(biāo).

采用工業(yè)機(jī)器人D-H方法建立的連桿坐標(biāo)系.在運(yùn)動學(xué)分析時，建立連桿坐標(biāo)系.每個連桿坐標(biāo)系轉(zhuǎn)換矩陣包含4個運(yùn)動學(xué)參數(shù)，即θi，αi，ai，di;6自由度工業(yè)機(jī)器人T0TB有24個運(yùn)動學(xué)參數(shù);BTM為一4×4的轉(zhuǎn)換矩陣，可用3個平移變量PBx，PBy，PBz及3個用歐拉角表示的旋轉(zhuǎn)變量RBx，RBy，RBz表示;PT0為一 4×1 的矩陣，根據(jù)其含義可用3個平移變量PTx，PTy，PTz表示.則機(jī)器人定位誤差測量系統(tǒng)的全閉環(huán)運(yùn)動學(xué)模型共包含33個運(yùn)動學(xué)參數(shù)，這33個參數(shù)誤差也是影響機(jī)器人靜態(tài)定位誤差的主要因素.

1.2 全閉環(huán)定位誤差數(shù)學(xué)模型

根據(jù)公式(1)，機(jī)器人定位誤差測量系統(tǒng)全閉環(huán)運(yùn)動學(xué)模型可變形為

其中，(x，y，z)為被測量點在測量坐標(biāo)系下的絕對坐標(biāo)值;gi(i=1～33)為運(yùn)動學(xué)模型中的33個運(yùn)動學(xué)參數(shù);f為以x，y，z為變量，以系統(tǒng)中33個運(yùn)動學(xué)參數(shù)為參數(shù)的三元方程組.當(dāng)運(yùn)動學(xué)參數(shù)無誤差時，也不存在機(jī)器人定位誤差;但實際情況運(yùn)動學(xué)參數(shù)總會存在一定的誤差，從而帶來機(jī)器人的定位誤差.

為簡化分析過程，對式(2)進(jìn)行一階微分，令

從而便獲得基于雅克比矩陣的機(jī)器人全閉環(huán)定位誤差數(shù)學(xué)模型為

2 定位誤差影響因素敏感度分析

影響機(jī)器人全閉環(huán)定位誤差的因素有33個，通過建立的基于雅可比矩陣的機(jī)器人全閉環(huán)定位誤差數(shù)學(xué)模型可將這些因素對定位誤差的影響分離開來，這種分離避免了因這些因素的交叉影響而帶來的對敏感度分析的復(fù)雜性.

為了分析出機(jī)器人在某一位姿時各個因素對定位誤差的SPS，令

其中，Θgk被稱為第 k個因素的機(jī)器人SFDPE;R(gi:gj)被稱為第i及第j個因素的RSFDPE;Θgi為機(jī)器人基準(zhǔn)單因素微分定位誤差，對應(yīng)的因素為基準(zhǔn)因素.

可用各因素間的RSFDPE來分析機(jī)器人在某個位姿時影響定位誤差的各因素的SPS，即機(jī)器人某個因素相對于基準(zhǔn)因素的RSFDPE越大，則表明此因素的SPS較高;反之，較低.

為了分析出機(jī)器人整個工作空間內(nèi)的各個因素的 MPCS，令

其中，ΣΘgk為第 k個因素在機(jī)器人N個位姿的SFCDPE;ΣR(gi:gj)為第i及第j個因素在機(jī)器人N個位姿的RSFCDPE;ΣΘgi為機(jī)器人基準(zhǔn)單因素綜合微分定位誤差，對應(yīng)的因素為基準(zhǔn)因素;為機(jī)器人第 m個位姿的第 i個因素的SFDPE;npim，npjm分別表示機(jī)器人第m個位姿的第i及第j個因素的SFDPE所在區(qū)間n的區(qū)間概率密度，可表示為

可用各因素間的RSFCDPE來分析各因素在機(jī)器人多個位姿對機(jī)器人定位誤差影響的MPCS，即機(jī)器人某個因素相對于基準(zhǔn)因素的RSFCDPE越大，則表明此因素的MPCS較高;反之，較低.

此外，為充分體現(xiàn)機(jī)器人實際定位誤差分布情況，可采用正交試驗法設(shè)計原理設(shè)計機(jī)器人N個位姿，即設(shè)計機(jī)器人位姿樣本空間.在此應(yīng)用中可采用6因素、5水平及25實驗次數(shù)的正交表，即L52(56).從而可獲得機(jī)器人位姿樣本空間:

其中，θmn表示樣本空間中第m個樣本點對應(yīng)的機(jī)器人第n個關(guān)節(jié)角值;vmn表示正交表中第m行第 n列對應(yīng)的值;minθn，maxθn分別表示機(jī)器人第n個關(guān)節(jié)角的最小值和最大值.

從定位誤差傳遞數(shù)學(xué)模型可知，影響機(jī)器人全閉環(huán)定位誤差的運(yùn)動學(xué)因素可以分為3類，下面分別分析各類因素的敏感度.

1)關(guān)節(jié)角因素分析:關(guān)節(jié)角因素敏感度反映了機(jī)器人6個關(guān)節(jié)角誤差對機(jī)器人全閉環(huán)定位誤差的影響程度，通過分析可以得到在機(jī)器人位姿樣本空間內(nèi)關(guān)節(jié)角敏感度分布情況.

2)各結(jié)構(gòu)角度量因素分析:各結(jié)構(gòu)角度量因素敏感度反映了機(jī)器人6個連桿坐標(biāo)系角度參數(shù)、3個機(jī)器人Base坐標(biāo)系與Measure坐標(biāo)系下的轉(zhuǎn)換矩陣的姿態(tài)分量誤差對機(jī)器人全閉環(huán)定位誤差的影響程度，通過分析可以得到在機(jī)器人位姿樣本空間內(nèi)各結(jié)構(gòu)角度量敏感度分布情況.

3)結(jié)構(gòu)平移量因素分析:各結(jié)構(gòu)平移量因素敏感度反映了機(jī)器人12個連桿坐標(biāo)系平移參數(shù)、3個機(jī)器人Base坐標(biāo)系與Measure坐標(biāo)系下的轉(zhuǎn)換矩陣的位置分量及3個測量點P在Measure坐標(biāo)系下的坐標(biāo)分量誤差對機(jī)器人全閉環(huán)定位誤差的影響程度，通過分析可以得到在機(jī)器人位姿樣本空間內(nèi)各結(jié)構(gòu)平移量敏感度分布情況.

3 定位誤差影響因素敏感度仿真

建立的工業(yè)機(jī)器人定位誤差測量實驗系統(tǒng)見圖2，它包含被標(biāo)定機(jī)器人、FARO ARM及測量靶標(biāo).工業(yè)機(jī)器人采用ABB IRB 1410機(jī)器人，此機(jī)器人是一種機(jī)身緊湊的機(jī)器人，承受載荷最高可達(dá)5 kg，具有較高的重復(fù)定位精度，通過對其絕對定位誤差進(jìn)行標(biāo)定后，可以大大拓展其應(yīng)用領(lǐng)域;誤差測量設(shè)備采用USB型FARO ARM，該測量設(shè)備具有測量精度高、操作簡單等優(yōu)點，可以直接測得機(jī)器人手臂上靶標(biāo)在其坐標(biāo)系下的坐標(biāo).

基于正交實驗法設(shè)計的ABB IRB1410機(jī)器人的位姿樣本空間的所有子空間中心點分布情況如圖3所示.

圖2 機(jī)器人定位誤差測量實驗系統(tǒng)圖

圖3 樣本子空間中心點分布圖

各因素的初始值都是已知量.根據(jù)這些信息并結(jié)合第3節(jié)的定量計算公式，便可計算出各個因素對機(jī)器人定位誤差的敏感度分布情況.經(jīng)由仿真結(jié)果看，機(jī)器人各結(jié)構(gòu)平移量因素的SPS和MPCS都一樣.因此，這里只給出機(jī)器人各角度因素SPS及MPCS的仿真結(jié)果，見圖4～圖7.

圖4 關(guān)節(jié)角因素敏感度仿真結(jié)果

圖5 連桿偏角因素敏感度仿真結(jié)果

圖4a為機(jī)器人關(guān)節(jié)角因素在機(jī)器人位姿樣本空間中的SPS仿真結(jié)果.其中，橫坐標(biāo)表示6個關(guān)節(jié)角，第1個關(guān)節(jié)角為基準(zhǔn)因素;縱坐標(biāo)為各因素相對于基準(zhǔn)因素的RSFDPE;25條線給出了不同樣本點對應(yīng)的各因素RSFDPE.圖4b為機(jī)器人關(guān)節(jié)角因素在機(jī)器人位姿樣本空間中的MPCS仿真結(jié)果.其中，橫坐標(biāo)表示6個關(guān)節(jié)角，第一個關(guān)節(jié)角為基準(zhǔn)因素;縱坐標(biāo)為各因素相對于基準(zhǔn)因素的RSFCDPE.

圖5a為機(jī)器人連桿偏角因素在機(jī)器人位姿樣本空間中的SPS仿真結(jié)果.其中，第1個連桿偏角為基準(zhǔn)因素.圖5b為機(jī)器人連桿偏角因素在機(jī)器人位姿樣本空間中的MPCS仿真結(jié)果.

圖6a為機(jī)器人基坐標(biāo)系姿態(tài)因素在機(jī)器人位姿樣本空間中的SPS仿真結(jié)果.其中，第1個基坐標(biāo)系姿態(tài)分量為基準(zhǔn)因素.圖6b為機(jī)器人基坐標(biāo)系姿態(tài)因素在機(jī)器人位姿樣本空間中的MPCS仿真結(jié)果.

圖7a為機(jī)器人3類角度基準(zhǔn)因素在機(jī)器人位姿樣本空間中的SPS仿真結(jié)果.其中，橫坐標(biāo)表示3個角度分量(依次為關(guān)節(jié)角-連桿偏角-基坐標(biāo)系姿態(tài))，第1個角度分量為基準(zhǔn)因素.圖7b為機(jī)器人3類角度因素在機(jī)器人位姿樣本空間中的MPCS仿真結(jié)果.

圖6 基坐標(biāo)系姿態(tài)因素敏感度仿真結(jié)果

圖7 3類因素敏感度對比仿真結(jié)果

4 結(jié)論

1)各個平移因素對機(jī)器人定位誤差影響的敏感程度是一樣的.

2)從各因素的SPS仿真結(jié)果看，在同一個機(jī)器人位姿樣本點，各個角度因素對機(jī)器人定位誤差影響的敏感程度是不一樣的;同一個角度因素，在機(jī)器人不同的位姿樣本點對機(jī)器人定位誤差影響的敏感程度也是不一樣的.

3)從各因素的MPCS仿真結(jié)果看，各關(guān)節(jié)角、各連桿扭角及基坐標(biāo)系姿態(tài)等因素對機(jī)器人定位誤差影響的敏感程度是不一樣的，且與所選的姿態(tài)樣本空間有關(guān).

這些分析與仿真結(jié)果為進(jìn)一步進(jìn)行機(jī)器人定位誤差的有效標(biāo)定打下了基礎(chǔ).

References)

［1］ Bai Y，Zhuang H Q，Roth Z S.Experiment study of PUMA calibration using a laser tracking system［C］//Proceedings of the IEEE International Workshop on Soft Computing in Industrial Application.Piscataway，NJ:IEEE，2003:139-144

［2］ Alici G，Shirinzadeh B.A systematic technique to estimate positioning errors for robot accuracy improvement using laser interferometry based sensing［J］.Mechanism and Machine Theory，2005，40:879-906

［3］ Yu Alexander，Bonev Ilian A，Zsombor-Murray Paul.Geometric approach to the accuracy analysis of a class of 3-DOF planar parallel robots ［J］.Mechanism and Machine Theory，2008，43:364-375

［4］ Fabricio Tadeu Paziani，Benedito Di Giacomo，Roberto Hideaki Tsunaki.Robot measuring form errors［J］.Robotics and Computer-Integrated Manufacturing，2009，25:168-177

［5］焦國太.機(jī)器人位姿誤差的分析與綜合［D］.北京:北京工業(yè)大學(xué)機(jī)械工程與應(yīng)用電子技術(shù)學(xué)院，2002 Jiao Guotai.Analysis and synthesis of robot pose errors［D］.Beijing:College of Mechanical Engineering and Applied Electronics Technology，Beijing University of Technology，2002(in Chinese)

［6］龔星如，沈建新，田威，等.工業(yè)機(jī)器人的絕對定位誤差模型及其補(bǔ)償算法［J］.南京航空航天大學(xué)學(xué)報，2012，44(4):60-64 Gong Xingru，Shen Jianxin，Tian Wei，et al.Absolute positioning error model of industrial robot and compensation algorithm［J］.Journal of Nanjing University of Aeronautics ＆ Astronautics，2012，44(4):60-64(in Chinese)

［7］ Ting K，Long Y.Performance quality and tolerance sensitivity of mechanisms［J］.Transactions of the ASME，1996，118:144-150

［8］郭劍鷹，呂恬生，王樂天.工業(yè)機(jī)器人運(yùn)動學(xué)參數(shù)誤差兩步識別法［J］.上海交通大學(xué)學(xué)報，2003，37(11):1670-1674 Guo Jianying，Lü Tiansheng，Wang Letian.Kinematic parameter error identification for industrial robot with two-step procedure［J］.Journal of Shanghai Jiaotong University，2003，37(11):1670-1674(in Chinese)

［9］王偉，贠超.機(jī)器人機(jī)構(gòu)精度綜合的正交實驗法［J］.機(jī)械工程學(xué)報，2009，45(11):18-24 Wang Wei，Yun Chao.Orthogonal experimental design to synthesizethe accuracy ofrobotic mechanism ［J］.Journalof Mechanical Engineering，2009，45(11):18-24(in Chinese)