黃昔光

(北方工業(yè)大學(xué) 機電工程學(xué)院,北京 100041)

廖啟征

(北京郵電大學(xué) 自動化學(xué)院,北京 100876)

空間 6R串聯(lián)機器人機構(gòu)位置逆解新算法

黃昔光

(北方工業(yè)大學(xué) 機電工程學(xué)院,北京 100041)

廖啟征

(北京郵電大學(xué) 自動化學(xué)院,北京 100876)

將倍四元數(shù)的復(fù)指數(shù)形式應(yīng)用于串聯(lián)機構(gòu)位置逆解分析中,提出了空間 6R(R代表轉(zhuǎn)動副)串聯(lián)機構(gòu)位置逆解新算法.基于倍四元數(shù)建立了空間 6R串聯(lián)機構(gòu)位置逆解的數(shù)學(xué)模型;然后,使用線性消元和 Dixon結(jié)式消元法,得到了 6×6的結(jié)式;由于采用未知轉(zhuǎn)角的復(fù)指數(shù)形式,不需要提取任何公因式,可直接獲得該機構(gòu)位置逆解的一元 16次輸入輸出方程和全部 16組封閉解.最后通過數(shù)字實例證明了該方法無增根無漏根.算例表明算法簡潔,易于程序?qū)崿F(xiàn),為串聯(lián)機構(gòu)位置逆解分析提供了新的理論基礎(chǔ).

串聯(lián)機構(gòu);倍四元數(shù);位置逆解

由轉(zhuǎn)動副、移動副以及圓柱副組成的空間串聯(lián)機構(gòu),結(jié)構(gòu)類型有數(shù)十種之多,其運動學(xué)位置逆解能否用統(tǒng)一的算法求解,一直是國內(nèi)外機構(gòu)學(xué)和機器人領(lǐng)域的研究目標(biāo).1980年,文獻[1]最早嘗試運用球面三角法進行統(tǒng)一求解,但最后剩下了 15種機構(gòu)無法解決.1986年,文獻[2]解決了喻為機構(gòu)運動分析中的珠穆朗瑪峰難題——7R機構(gòu)位移分析以后,理論上完成了全部串聯(lián)機構(gòu)的位置逆解,但針對不同類型的串聯(lián)機構(gòu)進行數(shù)學(xué)建模時,將產(chǎn)生不同的數(shù)學(xué)模型,需要不同的消元技巧,難以實現(xiàn)計算機程序自動求解.2006年,文獻[3]提出把串聯(lián)運動鏈拆成幾個簡單部分的組合,但該方法只適于某些解耦的特殊情況.以往用于串聯(lián)機構(gòu)位置逆解數(shù)學(xué)建模的方法主要有D-H矩陣法、球面三角法、實矩陣法、對偶數(shù)法等,得到了各不相同的逆解算法,不具有通用性.究其原因,主要是缺乏有效的數(shù)學(xué)建模方法,不能將D-H矩陣中的旋轉(zhuǎn)和平移統(tǒng)一表達(dá),難以建立串聯(lián)機構(gòu)逆解統(tǒng)一數(shù)學(xué)模型.隨著串聯(lián)機器人模塊化、可重構(gòu)技術(shù)的發(fā)展,研究 6R串聯(lián)機構(gòu)位置逆解新方法,探索串聯(lián)機構(gòu)統(tǒng)一逆解建模及通用算法,具有重要的理論價值與實際意義.

倍四元數(shù)把剛體三維空間的旋轉(zhuǎn)和平移統(tǒng)一為四維空間的純轉(zhuǎn)動,從而減少了機構(gòu)的種類.文獻[4]將三維空間中位移表示為四維空間中倍四元數(shù)表示的雙旋轉(zhuǎn);文獻[5]將倍四元數(shù)應(yīng)用于運動插值中;文獻[6]應(yīng)用倍四元數(shù)對空間 RR機器人進行綜合;文獻[7]運用倍四元數(shù)完成了空間 6自由度串聯(lián)機構(gòu)位置逆解,但該方法在構(gòu)造Dixon結(jié)式時,需要通過人工參與查找公因式來消除增根,未能實現(xiàn)該機構(gòu)位置逆解的自動求解.

本文采用未知轉(zhuǎn)角的復(fù)指數(shù)形式,改進了一種空間 6R串聯(lián)機構(gòu)位置逆解新算法.基于倍四元數(shù)建立了空間 6R串聯(lián)機構(gòu)位置逆解數(shù)學(xué)模型,構(gòu)造了 6×6的 Dixon結(jié)式,新算法不需要提取任何公因式,可直接獲得該機構(gòu)位置逆解的一元 16次方程及其全部封閉解,方法簡單實用,易于程序?qū)崿F(xiàn).

1 倍四元數(shù)形式的齊次變換矩陣

串聯(lián)機構(gòu)位置逆解數(shù)學(xué)建模通常采用 D-H矩陣法,相鄰關(guān)節(jié)兩坐標(biāo)系 i-1和 i之間的關(guān)系,可通過圖 1中各個參數(shù)來描述.

圖1 D-H相鄰坐標(biāo)變換

圖1表示坐標(biāo)系 i-1繞自身 z軸旋轉(zhuǎn) θi且沿z軸平移 si后得到一新的坐標(biāo)系;新坐標(biāo)系再繞其自身 x軸旋轉(zhuǎn) αi且沿 x軸平移 ai后得到 i坐標(biāo)系.其中,坐標(biāo)系 i-1繞自身 z軸旋轉(zhuǎn) θi且沿 z軸平移 si的齊次變換矩陣,運用倍四元數(shù)可表示為[5]

其中

由式(1)可知,坐標(biāo)系 i-1繞自身 z軸的旋轉(zhuǎn)和平移變換,運用倍四元數(shù)可以近似表達(dá)為四維空間中旋轉(zhuǎn)(θ+γ)角和(θ-γ)角的雙旋轉(zhuǎn),從而將 D-H矩陣中的旋轉(zhuǎn)和平移統(tǒng)一為四維空間的雙旋轉(zhuǎn).

為便于運算,式(1)可表示為 Clifford代數(shù)形式:

同理,坐標(biāo)系繞其 x軸旋轉(zhuǎn)角 αi且沿 x軸平移 ai的齊次變換矩陣用倍四元數(shù)形式表示為

其中

式中 ρi=ai/R.參考文獻[5],R=L/δ1/2,其中 L為機器人手臂所能達(dá)到的空間尺寸的最大值,δ為指定的精度.

2 運動學(xué)位置逆解

2.1 數(shù)學(xué)建模

6R串聯(lián)機構(gòu)位置逆解就是已知操作臂的末端空間位姿和機構(gòu)參數(shù) si,ai,αi,求各關(guān)節(jié)的輸入角 θi(i=1,2,…,6).6R串聯(lián)機構(gòu)有 6個關(guān)節(jié),需要經(jīng)過 6次圖 1中的空間三維運動,運動學(xué)方程可表示為

G為機器人操作末端位姿的倍四元數(shù)形式:

根據(jù) Clifford代數(shù)運算法則[8],式(6)可以表示為

由式(7)和式 (8)得

將式(2)、式 (5)代入式 (9)、式 (10),根據(jù)四元數(shù)各元素對應(yīng)相等,得

其中,Ei為 4×16的矩陣;Di為 4×4的矩陣,其元素均為已知結(jié)構(gòu)參數(shù)及輸入?yún)?shù)確定;

2.2 消元過程

采用分步進行消元.首先,對式(11)(取 i=1)進行線性消元,消去 θ1和 θ6,并將其代入式(11)(取 i=2);然后運用 Dixon結(jié)式消去 θ3,θ4,θ5;最后得到關(guān)于 θ2的 ±8次單變量方程.

當(dāng) i=1時,將式(11)看做關(guān)于 c1c6,c1s6,s1c6,s1s6的線性方程組,得

設(shè) F0(t2,t3,t4,t5)=[f1,f2,f3,f4]T,對式(16)構(gòu)造 Dixon結(jié)式如下:

其中

展開式(17)得

其中

矩陣 D6×6為一個只含變量 t2的 6×6方陣,即 Dixon矩陣.

由線性代數(shù)知識可知,式(18)成立的條件為

矩陣 D6×6每行分子分母關(guān)于變量 t2的最高次數(shù)為 1,2,1,1,2,1,其總和為 8.根據(jù)行列式的運算法則可知,展開式(18)后得到只含變量 t2的多項式分子分母最高次數(shù)不會超過 8.由式(19),不需要提取任何公因式,直接可得

其中,si是由已知參數(shù)確定的系數(shù).求解式(20)可得 t2的 16組解.

2.3 其它變量求解

將所求得的 16個 t2解代入下式:

根據(jù) Cramer法測,可線性求解對應(yīng) t2的 t4和 t5.

將所求的 t2,t4和 t5代入式(16)中任一式,即可求得對應(yīng)的 t3.從而求得 θ2,θ3,θ4,θ5,由式(12)～式(15)可求得 θ1和 θ6.

3 數(shù)值算例

采用本文的算法對文獻[8]中的數(shù)值算例求解,其 6R串聯(lián)機構(gòu)的結(jié)構(gòu)參數(shù)和位置參數(shù)如下:

表 1 數(shù)字實例的 4組實數(shù)解 (°)

4 結(jié) 論

1)將倍四元數(shù)的復(fù)指數(shù)形式應(yīng)用于串聯(lián)機構(gòu)位置逆解分析中,得出 8個約束方程,通過線性變換和 Dixon結(jié)式獲得了一元 16次方程及其及全部解析解,求解過程簡潔,易于程序?qū)崿F(xiàn),為 6R串聯(lián)機構(gòu)位置逆解提供了新的算法.

2)倍四元數(shù)把三維空間的轉(zhuǎn)動和移動統(tǒng)一為四維空間的純轉(zhuǎn)動,因此,本文提出的算法可適用于由轉(zhuǎn)動副、移動副及圓柱副三種運動副組成的空間 6自由度串聯(lián)機器人機構(gòu)位置逆解.

References)

[1]Duffy J.Analysis of mechanisms and robot manipulators[M].London:Ed wart Arnold Ltd,1980

[2]廖啟征,梁崇高,張啟先.空間 7R機構(gòu)位移分析的新研究[J].機械工程學(xué)報,1986,22(3):1-5 Liao Qizheng,Liang Chonggao,Zhang Qixian.New research of displacement analysis of spatial 7R mechanism[J].Chinese Journal of Mechanical Engineering,1986,22(3):1-5(in Chinese)

[3]趙杰,王衛(wèi)忠,蔡鶴皋.可重構(gòu)機器人封閉形式的運動學(xué)逆解計算[J].機械工程學(xué)報,2006,42(8):210-214 Zhao Jie,Wang Weizhong,Cai Hegao.Generation of closed-form inverse kinematics for reconfigurable robots[J].Chinese Journal of Mechanical Engineering,2006,42(8):210-214(in Chinese)

[4]Etzel K R,McCarthy JM.A metric for spatial displacement using biquaternions on SO(4)[C]//Proceedings of IEEE International Conference on Robotics and Automation.Piscataway,NJ:IEEE,1996:3185-3190

[5]Geng Q J,Varshney A,Menon JP,et al.Double quaternions for motion interpolation[C]//Proceedingsof the ASME DETC'98.Atlanta:ASME,1998:13-16

[6]McCarthy JM.Mechanismssynthesis theory and the design of robots[C]//Proceedings of IEEE International Conference on Robotics and Automation.San Francisco,CA:IEEE,2000:24-28

[7]喬曙光.6自由度串聯(lián)機械手位置逆解新方法[D].北京:北京郵電大學(xué)自動化學(xué)院,2008 Qiao Shuguang.New algorithm for inverse kinematics of 6 DOF serial manipulator[D].Beijing:School of Automation,Beijing U-niversity of Posts and Telecommunications,2008(in Chinese)

[8]Clifford W K.A prelim inary sketch of biquaternions[G]//Tucker R.Mathematical Papers.London:London Mathematical Society,1882:658

[9]于艷秋,王品,廖啟征.一般 6R機器人位置反解與運動仿真[J].中國機械工程,2003,14(24):2130-2132 Yu Yanqiu,Wang Pin,Liao Qizheng.Forward kinematics analysis of general 6R robot and its simulation[J].China Mechanical Engineering,2003,14(24):2130-2132(in Chinese)

(編 輯:文麗芳)

New algorithm for inverse kinematics of 6R serial robot mechanism

Huang Xiguang

(School of Mechanical and Electrical Engineering,North China University of Technology,Beijing 100041,China)

Liao Qizheng

(School of Automation,Beijing University of Posts and Telecommunications,Beijing 100876,China)

The theory of double quaternions and its application in the inverse kinematics of serial mechanisms was introduced.A new algorithm for the inverse kinematics of 6R mechanisms was presented based on the complex exponent form of double quaternions.Based on double quaternions,a mathematical model of 6R mechanisms was created.Then,a 6×6 resultant matrix was obtained directly by using linear elimination and Dixon resultant method,without factoring outor deriving the greatest common divisor,due to the proposed algorithm used the complex exponent form of double quaternions.A 16th degree univariate equation was achieved from the determinant of the matrix and all 16 closed-form solutions were also obtained.The proposed algorithm is comparably easy and simple to program.It was verified by a numerical example that the obtained roots satisfy the original equations.The research result provides a new method for the inverse kinematics of serial mechanisms.

serial mechanism;double quaternions;inverse kinematics

TH 112

A

1001-5965(2010)03-0295-04

2009-03-05

國家自然科學(xué)基金資助項目(50775012);北京市屬院校人才強教計劃資助項目;北京市特色專業(yè)建設(shè)資助項目;北方工業(yè)大學(xué)?？蒲谢鹳Y助項目

黃昔光(1979-),男,湖南岳陽人,講師,huangxiguang@gmail.com.