王昊，張得禮，王宏濤

(南京航空航天大學(xué) 機電學(xué)院，江蘇 南京 210016)

0 引言

工業(yè)機器人由于其工作靈活、可移動和高重復(fù)定位精度等特點被廣泛應(yīng)用于搬運、噴漆、銑削、焊接、制孔等工業(yè)加工領(lǐng)域。工業(yè)機器人與CNC機床一樣具有多軸功能，然而由于工業(yè)機器人獨特的懸臂梁結(jié)構(gòu)，與CNC機床相比剛度相對較差。CNC機床的剛度通常>50 N/μm[1]，而工業(yè)機器人的剛度僅為CNC機床的1/50～1/20，提高工業(yè)機器人的剛度對提高加工質(zhì)量具有重要影響。工業(yè)機器人末端剛度由機器人關(guān)節(jié)剛度和機器人位姿決定[2]，所以對于型號已定的工業(yè)機器人，其位姿對末端剛度有著決定性作用。此外，機器人傳遞力和運動的運動學(xué)性能、規(guī)避奇異點的能力對機器人的加工性能也有較大影響[3]。國內(nèi)外學(xué)者有關(guān)研究表明，可以通過尋優(yōu)機器人加工位姿達到提高機器人末端剛度性能和運動學(xué)性能的目的。CHEN S F等建立了工業(yè)機器人末端剛度和關(guān)節(jié)剛度、機器人姿態(tài)的關(guān)系，并表明末端剛度隨機器人姿態(tài)變化而變化[4]；ZARGARBASHI S H H等人提出一種與任務(wù)相關(guān)的優(yōu)化方法，對機器人姿態(tài)進行尋優(yōu)，從而選取滿足任務(wù)要求的最優(yōu)加工姿態(tài)[5]；KHAN W等人基于機器人運動性能評價指標(biāo)，研究了機器人不同運動性能對機器人加工軌跡的影響[6]。然而上述優(yōu)化方法僅考慮優(yōu)化機器人的剛度性能或機器人的運動性能，并沒有將二者綜合考慮。LIN Y等人在對機器人加工位姿的優(yōu)化時，對機器人的運動性能和剛度性能皆進行了研究，但其選擇優(yōu)化位姿的方式是在性能分布圖中直接選擇機器人關(guān)節(jié)角，未使用優(yōu)化算法[7]。

本文基于機器人性能指標(biāo)并結(jié)合遺傳算法對機器人加工位姿進行優(yōu)化。

1 機器人運動性能對加工位姿的影響

1.1 機器人正運動學(xué)模型

埃斯頓ER30機器人為六自由度關(guān)節(jié)式工業(yè)機器人，且所有關(guān)節(jié)均為旋轉(zhuǎn)副。按照D-H參數(shù)法建立機器人D-H參數(shù)模型如表1所示。表1中，關(guān)節(jié)1-6分別代表機器人的6個關(guān)節(jié)；θi為關(guān)節(jié)i的關(guān)節(jié)角度；di為關(guān)節(jié)i的偏距，即連桿i-1與連桿i的距離；αi為關(guān)節(jié)i-1與關(guān)節(jié)i軸線之間的夾角；ai-1為關(guān)節(jié)i-1與關(guān)節(jié)i軸線之間公法線距離。

表1 ER30機器人D-H參數(shù)

(1)

1.2 機器人運動性能評價指標(biāo)

機器人運動性能是指機器人系統(tǒng)傳遞力和運動的能力，機器人的運動性能可以通過運動性能指標(biāo)來評價。機器人運動性能指標(biāo)除了可以評價機器人的運動性能還可以直觀地衡量機器人當(dāng)前位姿與奇異點距離的遠近。本文選擇基于條件數(shù)[8]的機器人運動性能的評價指標(biāo)，條件數(shù)定義見式(2)，機器人運動性能評價指標(biāo)KM見式(3)。

(2)

KM=1/K(JN)×100%

(3)

式中：n為雅克比矩陣的維數(shù)，本文機器人為6R機器人，所以n=6；tr(JN)為機器人雅克比矩陣的跡；JN為雅克比矩陣的規(guī)范標(biāo)準(zhǔn)型。

1.3 各關(guān)節(jié)角對運動性能的影響

在機器人的6個關(guān)節(jié)角中，θ1、θ6對機器人末端運動性能無影響[7]，所以設(shè)置θ1=0,θ6=0，僅需研究θ2、θ3、θ4、θ5對機器人末端運動性能的影響。式(1)中機器人空間位置px、py、pz主要隨θ2、θ3變化而變化，機器人末端姿態(tài)n、o、a主要隨θ4、θ5變化而變化，所以可通過KM在2、3關(guān)節(jié)空間中的變化，研究機器人空間位置變化對KM的影響。通過KM在4、5關(guān)節(jié)空間中的變化，研究姿態(tài)變化對KM影響。

為研究θ2、θ3對KM的影響，通過蒙特卡洛法繪制KM在2、3關(guān)節(jié)空間的分布如圖1(a)所示。圖1(b)為機器人末端運動性能在笛卡兒空間的分布。在R1、R2、R3三個矩形區(qū)域中機器人的末端運動性能較好。通常情況下，機器人進行加工時，待加工工件放于機器人前方，所以當(dāng)機器人進行加工時，機器人末端通常處在R1區(qū)域，此區(qū)域內(nèi)的機器人加工位姿符合加工實際且機器人運動性能較佳，此時θ2∈[-1,0.5]，θ3∈[-0.8,1.5]。

圖1 機器人運動性能在關(guān)節(jié)空間中的分布

同理，繪制KM在4、5關(guān)節(jié)空間的分布如圖2所示。由圖2可知，KM主要隨θ5變化而變化。在考慮姿態(tài)對KM的影響時，只需優(yōu)化θ5的取值，以避免機器人四軸與六軸平行導(dǎo)致接近奇異點，機器人末端即可達到較好的運動性能。

圖2 機器人運動性能4、5關(guān)節(jié)空間中的分布

綜合上文分析，當(dāng)機器人四軸與六軸不平行遠離奇異時，θ2、θ3變化導(dǎo)致的KM變化范圍為[0,0.7]，θ4、θ5變化導(dǎo)致的KM變化范圍為[0.5,0.7]。機器人末端運動性能主要受機器人位置變化即θ2、θ3變化的影響，機器人姿態(tài)變化即θ4、θ5變化對機器人運動性能影響較小。

2 機器人剛度性能對加工位姿的影響

2.1 機器人靜剛度模型

對于6R串聯(lián)機器人，加工時其末端誤差主要來自于機器人關(guān)節(jié)的受力變形[9]?；诖丝梢越C器人末端笛卡兒剛度矩陣K與關(guān)節(jié)剛度矩陣Kθ的映射關(guān)系[10]：

K=J-TKθJ-1

(4)

式中：J為機器人雅克比矩陣；Kθ=diag([Kθ1,Kθ2,Kθ3,Kθ4,Kθ5,Kθ6])，為機器人關(guān)節(jié)剛度矩陣。

由式(4)可知，機器人的末端剛度與機器人的雅克比矩陣J和關(guān)節(jié)剛度矩陣Kθ有關(guān)。機器人的雅克比矩陣J僅由機器人位姿決定，所以機器人末端剛度具有位姿依賴性。

2.2 機器人剛度評價指標(biāo)

為減少計算難度和計算誤差，引入機器人柔度模型：

(5)

由于機器人末端所受力矩帶來的末端變形遠小于力帶來的末端變形。機器人末端的角位移也遠小于線位移，所以機器人剛度性能評價指標(biāo)Ks可簡化為[11]

(6)

2.3 各關(guān)節(jié)角對剛度性能的影響

由參考文獻[7]可知，θ1、θ6對機器人剛度評價指標(biāo)Ks無影響，設(shè)置θ1=0，θ6=0。機器人剛度性能在2、3關(guān)節(jié)空間的映射如圖3所示。對比圖3(a)和圖3(b)，當(dāng)θ4、θ5變化時，機器人末端剛度性能在2、3關(guān)節(jié)空間的分布基本不變。當(dāng)機器人運動性能較好，即θ2∈[-1,0.5]、θ3∈[-0.8,1.2]時，由圖可知，機器人的剛度性能在θ2→0.5，θ3→-0.8時較好。

圖3 2、3關(guān)節(jié)空間的機器人末端剛度性能

圖4為機器人末端剛度性能在機器人4、5關(guān)節(jié)空間中的映射。由圖4可知，當(dāng)θ2、θ3為隨機值時，機器人剛度性能的分布并無明顯規(guī)律；θ2、θ3為固定值時，機器人剛度性能分布呈現(xiàn)明顯規(guī)律，此時Ks取值區(qū)間為[2,4]。結(jié)合圖3，無論θ4、θ5為何值時，機器人剛度性能分布規(guī)律不變，Ks取值區(qū)間為[1,8]。由此可知θ2、θ3對機器人剛度性能有較大影響，θ4、θ5對機器人剛度性能影響較小，即相較于機器人末端的位置變化，機器人末端的姿態(tài)變化對機器人剛度性能影響較小。

圖4 4、5關(guān)節(jié)空間機器人末端剛度性能

3 基于遺傳算法的機器人加工位姿尋優(yōu)

根據(jù)上文的分析，可以歸納得到機器人加工位姿優(yōu)化問題的數(shù)學(xué)模型，如式(7)所示。這是一個多目標(biāo)、多約束且目標(biāo)數(shù)學(xué)表達式不可導(dǎo)的優(yōu)化問題，解決此優(yōu)化問題需要在變量全局范圍內(nèi)對最優(yōu)解進行搜索。遺傳算法對優(yōu)化問題沒有太多數(shù)學(xué)要求，其較強的全局尋優(yōu)能力的特點符合本文優(yōu)化問題的要求。所以選用遺傳算法實現(xiàn)機器人位姿優(yōu)化。

(7)

1)使用遺傳算法尋優(yōu)機器人加工位姿需先設(shè)計變量和適應(yīng)度函數(shù)，設(shè)置選擇、交叉、變異算子，設(shè)置終止條件，步驟如下。

步驟1：設(shè)計變量和取值范圍。機器人末端運動性能和剛度性能主要受機器人位姿影響，其中θ1、θ6對機器人末端運動性能和剛度性能均無影響，所以待優(yōu)化變量為θ2、θ3、θ4、θ5。結(jié)合上文分析設(shè)置各變量取值范圍為：

θ2∈[-1,0.5],θ3∈[-0.8,1.2],θ4∈[-π,π],

θ5∈[-2.5，-0.5]∪[0.5,2.5]。

步驟2：設(shè)計適應(yīng)度函數(shù)。適應(yīng)度函數(shù)有兩個評價指標(biāo)KM和Ks，以KM和Ks的和作為適應(yīng)度函數(shù)。參考文獻[12]考慮到在加工時機器人的運動性能與剛度性能對加工效果都有明顯影響，設(shè)置權(quán)重ω1=2/3，ω2=1/3。最終的適應(yīng)度函數(shù)如式(8)所示。

Fit(θ2，θ3，θ4，θ5)=0.000 1ω1Ks+ω2KM

(8)

步驟3：設(shè)置選擇、交叉、變異算子。使用輪盤賭選擇法對種群中的個體進行選擇。為快速淘汰低適應(yīng)度值的個體，交叉方法為均勻交叉，交叉算子pc=0.8。變異算子pm通常在[0.001,0.1]區(qū)間取值，本文變異算子pm=0.05。

步驟4：設(shè)置終止條件。遺傳算法終止條件有多種，本文以達到目標(biāo)代數(shù)作為終止條件，當(dāng)?shù)鷶?shù)達到50代時停止優(yōu)化。

2)完成以上設(shè)置后，使用遺傳算法對機器人加工位姿進行優(yōu)化，步驟如下所示。

步驟1：生成初始種群。隨機生成符合取值范圍的25組關(guān)節(jié)角θ2i,θ3i,θ4i,θ5i(i=1～25)，并以二進制編碼的方式將25組關(guān)節(jié)角編譯成25個染色體作為初始種群。

步驟2：計算Fit。計算初始種群中25個染色體Fit(θ2，θ3，θ4，θ5)。

步驟4：判斷是否終止遺傳算法。若代數(shù)g<50，則g+1，進行步驟5；若代數(shù)g>50,則停止優(yōu)化，輸出優(yōu)化結(jié)果。

步驟5：遺傳操作生成下一代子群。根據(jù)當(dāng)前種群中各染色體Fit值，對種群中所有個體進行選擇、復(fù)制、交叉以得到新的種群。進行步驟2。

運行遺傳算法得到的優(yōu)化結(jié)果如表2所示。此優(yōu)化結(jié)果對應(yīng)的機器人位姿如圖5所示。優(yōu)化位姿與隨機位姿的性能指標(biāo)值如表3所示，相較于隨機位姿優(yōu)化后的位姿各項性能得到了較大提升。

表2 遺傳算法優(yōu)化結(jié)果(θ1、θ6=0)

圖5 優(yōu)化后的機器人加工位姿

表3 優(yōu)化位姿與隨機位姿的各評價指標(biāo)值

4 結(jié)語

機器人6個關(guān)節(jié)中，θ2、θ3對機器人運動性能影響較大，θ4、θ5對機器人運動性能影響較小，θ1、θ6對機器人運動性能無影響。當(dāng)θ2∈[-1,0.5],θ3∈[-0.8,1.2],θ4∈[-π,π],θ5∈[-2.5,-0.5]∪[0.5,2.5]時機器人運動性能較好且機器人位姿符合加工實際。

使用遺傳算法對加工位姿進行優(yōu)化，設(shè)計了一種綜合機器人運動性能評價指標(biāo)和機器人剛度性能評價指標(biāo)的適應(yīng)度函數(shù)，并通過系數(shù)設(shè)計調(diào)整兩種性能在優(yōu)化過程中的比重，最終得到了一個符合加工實際且性能較好的機器人加工位姿。