溫貽芳 孫立寧 徐 朋
(①蘇州工業職業技術學院,江蘇 蘇州 215104;②天津大學機械工程學院,天津 300072;③蘇州大學機電學院,江蘇 蘇州 215100; ④南京航空航天大學機電學院,江蘇 南京 210016)
隨著科技進步和智能制造業的飛速發展,對制造件的表面性能要求越來越高[1]。等離子噴涂是一種比較成熟的表面改性方法,等離子體表面改性技術以等離子體噴槍為載體,無需低壓真空設備和隨之而來的高昂維護費,還可使等離子體產生區和應用區分離,是目前國際上等離子體科學和技術領域的研究熱點之一。具有低孔隙率、高結合強度以及低雜質的優點,被廣泛應用于航空航天[2]、生物醫學[3],陶瓷制造[4]以及再制造[5]等方面。傳統的等離子噴涂是通過人工操作實現,定位精度低,靈活性差,產生的化學物質、熱輻射、噪音等對人體也會產生危害。所以將工業機器人應用于等離子噴涂中,不但容易適應操作環境,還能提高噴涂成型的工藝質量[6]。通常等離子噴涂所面對的產品多為形狀復雜,尺寸變化大,需要噴槍具有更高的運動靈活性。為此通過增加機器人自由度實現噴涂運動靈活性,這里簡稱為冗余等離子噴涂機器人[7-9]。但是冗余自由度的增加也同時使得機器人操作性能指標問題變得更加復雜[10]。所以為了更好地評價冗余等離子噴涂機器人的操作性能,研究新的復合性能指標具有重要的意義。
等離子噴涂機器人操作性能指標主要有基于雅可比矩陣的條件數和可操作度等度量指標。條件數本質上是雅可比矩陣變換均一性在機器人各個運動方向上的表現形式[11-12]。條件數越小,機器人運動靈活性越好,速度也會越均勻,同時在各個方向上的運動能力也就越一致。可操作度主要表示機器人整體運動的靈活性,本質上是機器人運動速度在各個方向上的綜合度量。
在機器人條件數和可操作度性能指標研究方面,SALISBURY等[13]基于雅克比矩陣提出了條件數概念,并利用條件數方法評價了機器人的末端運動速度工作空間。YOSHIKAWA[14]利用雅克比矩陣行列式方法提出機器人的可操作度概念,并利用可操作度橢球對其幾何意義進行了相關描述。姚建初等[15]提出基于運動速度方向的可操作度指標,增加了特定方向的運動靈活性,缺點是該指標忽略了條件數對機器人性能的影響。張建富等[16]在可操作度基礎上利用蒙特卡洛法對并聯機器人的隨機點分布不均勻問題進行了研究,得到了空間點密度與雅克比行列式的對應關系式。張博等[17]在傳統可操作度基礎上,系統分析了機器人基座與輔助臂的耦合效應對末端任務空間的影響。謝碧云等[18]綜合考慮了條件數和方向可操作度對機器人運動操作速度的影響,提出基于條件數的方向可操作度指標,但是沒有明確分析各指標對機器人操作性能的影響程度,工作過程中不能對機器人末端運動操作性能進行有效控制。因此,為了評價機器人操作性能在某一特定位形下的運動能力,本文在條件數和速度方向可操作度基礎上提出利用Adaboost加權尋優算法的復合評價指標方法。
Adaboost算法的基本思想是通過不斷迭代將多個“弱”性能指標進行融合產生更加有效的綜合性能指標[19-20]。目前廣泛應用于人臉識別、分類與預測領域。本文將條件數和速度方向可操作度作為弱性能指標,通過迭代計算其占有的不同比重,然后將上述弱性能指標與權重做積加權,得到強性能指標[21-22]。最后利用強性能評價指標對冗余等離子噴涂機器人操作性能進行仿真分析。仿真結果表明:復合性能指標相比條件數和方向可操作度指標具有更強的可操作性,可以適應更復雜的等離子表面改性工作。
圖1為冗余等離子噴涂機器人工作模型。由六自由度ABB機械手、單自由度滑軌及旋轉殼體組成。圖2所示為冗余等離子噴涂機器人在初始工作狀態時的運動旋量模型。{S}表示慣性坐標系;{T}表示工具坐標系;a1-a6分別為機械手的尺寸參數;d1為冗余關節變量值;r1~r5為各自軸線上的點;ξ1~ξ7為各自軸線上的單位運動旋量。


冗余機器人末端運動速度評價指標主要有條件數和速度方向可操作度兩種,分別如式(1)、 (2)所示,J表示雅克比矩陣,v表示等離子噴涂機器人末端運動速度。
k=‖J‖‖J-1‖
(1)
u=(v(JJT)-1v)-1
(2)
當條件數過大時,冗余等離子噴涂機器人的運動方程就容易處于病態,而且會嚴重影響機器人的運動學逆解求解精度,導致表面改性質量變差,產生過大孔隙率、過熔等不良現象。速度方向可操作度只考慮了等離子噴涂機器人在末端運動方向的速度轉換能力,沒有考慮機器人其他運動方向的速度變換。因此用該方法對等離子噴涂機器人進行運動操作速度優化的時候,機器人容易出現奇異位形或條件數變大的位形,這樣機器人會失去其運動靈活性,從而影響等離子表面改性質量。
基于以上問題,本文在條件數和速度方向可操作度基礎上提出一種基于Adaboost加權尋優的復合性能評價方法。該方法不僅包含條件數和速度方向可操作度,而且通過對加權系數的設置,使得條件數永遠小于任務要求的最大條件數,在此基礎上使得速度方向可操作度盡量取最大值,最終使得等離子噴涂機器人具有更好的操作靈活性。
AdaBoost算法是基于Boosting算法的一種迭代算法,其基本思想是通過不斷迭代將多個“弱”性能指標進行合并而產生更加有效的綜合性能指標。主要步驟為:首先給定樣本空間(x,y),將訓練樣本分為m組,每組樣本占有1/m;然后通過boosting算法進行迭代,依據組合結果增加組合性能差的訓練個體權重從而對訓練樣本權重分布進行更新。最終得到一組綜合性能指標為P1,P2,P3,…,Pn,其中各指標占有不同比重,性能指標優的所占比重大。然后將上述指標與權重做積加權,得到強性能指標。具體實現步驟如下:
根據泛函理論,最小化處理經驗風險泛函如下:
(3)
具體步驟如下:
(1)對風險泛函進行迭代,第k次時構造函數:
(4)
其中:k為迭代次數;φr(x),r=1,…,l為指示函數集;βk=(d1,d2,…,dk)是一個k維向量。
(2) 假設第k次迭代時經驗風險值為:
(5)
則在k+1次迭代時,對下面式(6)進行最小化:
f(x,βk+1)=f(x,βk)+dk+1φk+1(x)
(6)
得到以下經驗風險值:
(7)
為了將式(7)最小化,定義如下參數值:
(8)


(9)
這是由下面的事實所決定的:yiφk+1(xi)∈{1,-1}, 式(7)在最優點dk+1上取得最小值,即:

(10)
(3)以此類推,第k+i次時有下面的等式成立:
(11)

(12)
根據泛函理論,最小化處理經驗風險泛函如下:
(13)
其中:
(14)
式(13)就是Adaboost決策規則。其具體算法如下:
(1)樣本選擇和網絡初始化。首先在樣本空間中隨機選擇m組樣本數據,對樣本數據中的每一個樣本的分布權值初始化為Di(i)=1/m。
(2)弱分類器預測。訓練第t個弱分類器時,通過綜合性能評價的結果和上一個性能指標評價結果進行比較做差從而得到誤差ei:
(15)
(3)計算預測序列權重。根據綜合性能評價誤差ei得到權重計算公式at為:
(16)
(4)測試樣本權重調整:
(17)
其中Bt定義為歸一化因子,一般在權重比例不變的情況下所有局部性能指標對應的分布權值為1。
(5)得到強分類器函數。多次迭代結束后,通過加權組合得到強綜合性能公式為:
(18)
基于條件數和速度方向可操作度的加權復合性能指標如式(19)所示。
γNew=ω1δ‖J‖J-1‖+ω2δ(v(JJT)-1v)-1=ω1γ1+ω2γ2
(19)
(20)
式(19)中:ω1、ω2分別代表加權系數,δ根據不同情況來定,一般1.2≤δ≤5。由于條件數和方向可操作度有著相反的變化趨勢,所以確定最優的權系數ω1、ω2成為衡量復合性能指標γNew優劣的重要標準。
根據冗余等離子噴涂機器人的工作特性,本文所討論性能指標同時受到機器人條件數及速度方向可操作度的約束,所以在優先保證性能指標最大條件數小于任務規定數值的同時,機器人的方向可操作度指標取最大值。
式(20)中γ1和γ2被定義為虛擬靈活性指標1和2。那么上述問題就簡化為在滿足較小條件數時,方向可操作度盡可能取最大值時,確定一組最優的w1和w2參數。
具體步驟為:首先將靈活性指標γ1到靈活性指標γn看作弱分類器,然后通過上述Adaboost算法,將靈活性指標進行組合,得到強靈活性指標。原理如圖3所示。
對于本文,只有兩個靈活性指標,對應于文中的虛擬靈活性指標1和2。通過Adaboost算法對上述虛擬靈活性指標進行性能提升。以條件數和速度方向操作度作為迭代條件,得到一組最優w1和w2參數,使得靈活性指標隨時間變化時能保證較小的條件數和較大的方向操作度。最終w1和w2參數分別為0.235、0.635,其中過程迭代如圖4所示。


確定了加權系數的數值,就可以利用新的復合評價指標γNew對冗余等離子機器人運動學性能進行仿真分析,并且與條件數和速度方向可操作度進行對比分析。根據圖2所示冗余等離子噴涂機器人的結構參數,末端運動速度取v=(0.04 0.02 0)Tm/s;初始運動位置定義為:

(21)
仿真時間定義為t=5 s,采樣周期取dt=50 ms,任務規定的最大條件數取18,取δ=3。
(22)
仿真結果如圖5、6所示,分別表示優化指標γNew、γ1、γ2基于條件數和方向可操作度的變化規律。由圖5可知,在1.5 s 本文為了提高冗余等離子噴涂機器人的操作靈活性和穩定性,有效提高等離子表面改性質量,在條件數和速度方向可操作度基礎上,提出一種分析冗余等離子噴涂機器人操作靈活性的復合評價指標,并利用Adaboost加權尋優算法對復合評價指標的加權系數進行確定。最后通過仿真分析證明本文所提方法在冗余等離子噴涂機器人操作靈活性和穩定性方面,比單一條件數法和速度方向可操作度法更加合理有效。

4 結語