毛冰滟 謝志江 吳小勇 袁岳軍

(重慶大學機械傳動國家重點實驗室， 重慶 400044)

基于引導人工蜂群算法的3-RPS并聯(lián)機構(gòu)正解優(yōu)化

毛冰滟 謝志江 吳小勇 袁岳軍

(重慶大學機械傳動國家重點實驗室， 重慶 400044)

提出了一種求解3-RPS并聯(lián)機構(gòu)正解優(yōu)化的快速數(shù)值解算法，采用數(shù)值優(yōu)化方式將正解轉(zhuǎn)換為最小化問題。針對3-RPS并聯(lián)機構(gòu)位姿正解優(yōu)化，采用數(shù)值法建立該并聯(lián)運動平臺正解方程，其實質(zhì)為帶約束的多目標非線性方程組。建立判定方程實現(xiàn)方程組單目標優(yōu)化，采用引導人工蜂群算法進行最小值優(yōu)化求解。該算法既利用了基本人工蜂群算法中鄰近蜜蜂交換蜜源信息的方式，又采用全局最優(yōu)蜜蜂引導所有蜜蜂往更優(yōu)蜜源處移動，更快速地搜索到最優(yōu)蜜源位置。通過求解3-RPS正解數(shù)值解優(yōu)化仿真算例分析，對比改進蟻群算法和基本人工蜂群算法，結(jié)果表明引導人工蜂群算法是高精度高速求解并聯(lián)運動平臺正解優(yōu)化的一種有效方法。

并聯(lián)機構(gòu)； 人工蜂群算法； 運動學正解

引言

與串聯(lián)機器人相比，并聯(lián)機構(gòu)具有高承載能力、高剛度、高精度等特點[1-2]。為達到運動精度高、可達精度小，機構(gòu)的綜合問題分析在于運動學求解優(yōu)化，并聯(lián)機構(gòu)的運動學優(yōu)化關(guān)鍵在于機構(gòu)正解，即已知驅(qū)動得到位置和姿態(tài)[3-5]。

并聯(lián)機構(gòu)的運動學正解包含一組復雜的多目標耦合非線性方程組，傳統(tǒng)采用的方式有數(shù)值法和解析法。任文博等[6-7]得到并聯(lián)機構(gòu)的封閉解析解，需要根據(jù)機構(gòu)的結(jié)構(gòu)特點進行判斷，效率低，計算量大。數(shù)值法通常采用Newton-Raphson法進行迭代求解[8-9]，需計算雅可比矩陣及其逆矩陣，計算工作量大，結(jié)構(gòu)受初值影響，計算精度達不到要求[8-11]。

近年來，智能優(yōu)化算法發(fā)展迅速，如粒子群算法(Particle swarm optimization, PSO)[12-15]、遺傳算法(Genetic algorithm, GA)[16]、連續(xù)蟻群算法(Ant colony algorithm, ACO)[17]等在機器人運動學求解上得到了應(yīng)用。粒子群算法運算速度快，算法易于應(yīng)用，但全局尋優(yōu)性差，容易陷入局部最優(yōu)值；遺傳算法全局并行搜索，但常規(guī)的算法有收斂速度慢及早熟的缺點；連續(xù)蟻群算法為群智能算法模式，但運行速度慢。人工蜂群算法(Artificial bee colony, ABC)[18-25]仿生于蜜蜂群體采蜜時各個蜂種之間的配合，以找到最優(yōu)采蜜位置。算法控制參數(shù)少，全局尋優(yōu)性好，收斂速度快，具有非常好的優(yōu)化搜索特性，可以有效地解決實際工程問題。利用改進的引導人工蜂群算法(Global-best artificial bee colony, GABC)，通過利用全局最優(yōu)蜜蜂進行引導學習，帶領(lǐng)全部蜜蜂向更優(yōu)蜜源移動，加快算法收斂速度。

本文對3-RPS并聯(lián)運動平臺進行運動學分析，將運動學正解多目標非線性方程組規(guī)劃問題轉(zhuǎn)換為單目標優(yōu)化，建立判定方程。采用人工蜂群算法以及改進的引導人工蜂群算法進行評定方程優(yōu)化求解。

1 并聯(lián)機構(gòu)運動學分析

1.1 運動學逆解

圖1為3-RPS并聯(lián)機構(gòu)結(jié)構(gòu)簡圖。該機構(gòu)是由上、下平臺以及連接兩平臺的3個支鏈組成。上平臺為運動平臺，能實現(xiàn)3-RPS的3個自由度，下平臺為固定平臺，位置保持不變。支鏈從下往上由轉(zhuǎn)動副(R副)、移動副(P副)和球副(S副)組成，轉(zhuǎn)動副與固定平臺連接且采用等邊三角形的布置方式，球副與運動平臺連接，布置方式與轉(zhuǎn)動副一一對應(yīng)。3條支鏈的長度可通過移動副的伸長縮短進行改變，從而使得運動平臺位置和姿態(tài)改變。

圖1 3-RPS并聯(lián)機構(gòu)結(jié)構(gòu)簡圖Fig.1 Structure sketch of 3-RPS parallel manipulator

以右手定則在固定平臺中心建立固定坐標系OXYZ，X軸指向轉(zhuǎn)動副中心A1，Z軸垂直固定平臺平面向上。在運動平臺上相應(yīng)的建立動坐標系O1X1Y1Z1，X1軸指向球副中心B1，Z1軸垂直運動平臺向上。固定坐標系不變，該平臺可實現(xiàn)3個自由度為繞X軸和Y軸方向的轉(zhuǎn)動，以及沿Z軸方向的移動。

由機構(gòu)的封閉矢量可得

LAiBi=RLO1Bi+P-LOAi(i=1,2,3)

(1)

其中

(2)

式中R——動坐標系到定坐標系的Z-Y-X歐拉旋轉(zhuǎn)變換矩陣

P——動坐標系到定坐標系的位置變換矩陣

α、β、γ——繞X軸、Y軸、Z軸的轉(zhuǎn)動角度

xe、ye、ze——繞X軸、Y軸、Z軸的移動距離

其中sα=sinα，cα=cosα，以此類推。

根據(jù)機構(gòu)的幾何關(guān)系可以確定A點和B點分別在定坐標系和動坐標系系中的位置，即

(3)

其中

式中ra——轉(zhuǎn)動副所在平臺半徑rb——球副所在平臺半徑

根據(jù)并聯(lián)機構(gòu)實際的裝配關(guān)系可知，轉(zhuǎn)動副R的軸線ei與移動副P的移動方向保持垂直，則可以得到運動學約束方程

(RLO1Bi+P-LOAi)Tei=0

(4)

則可以得到3個轉(zhuǎn)動副R的軸線向量為

(5)

通過式(3)～(5)可得到xe、ye、γ關(guān)于ze、α、β的表達式

(6)

式(6)為該并聯(lián)機構(gòu)的運動學約束方程。在帶約束情況下求該三自由度并聯(lián)機構(gòu)的位置正解，能夠得到該機構(gòu)的唯一驅(qū)動。將式(6)代入式(1)中，使得在逆解的非線性方程組中只含有代表機構(gòu)位置和姿態(tài)的3個自由度ze、α、β。從而可以得到ze、α、β表示向量LAiBi。

通過式(7)～(9)可得，機構(gòu)實現(xiàn)3個自由度的驅(qū)動為移動副的驅(qū)動位移，則機構(gòu)逆解為

(7)

(8)

(9)

(10)

1.2 運動學正解多目標非線性模型優(yōu)化

機構(gòu)正解為已知驅(qū)動位移di的情況下得到機構(gòu)的位姿α、β、ze。在得到機構(gòu)逆解的情況下，正解的本質(zhì)在于求解一組隱式非線性方程組，即

di-fi(α,β,ze)=0 (i=1,2,3)

(11)

采用權(quán)和法將該非線性方程組轉(zhuǎn)換為單目標線性規(guī)劃問題，其實質(zhì)為帶約束的多目標非線性方程優(yōu)化，轉(zhuǎn)換為最小化問題，則

(12)

因此得到該并聯(lián)運動平臺數(shù)值正解的判定方程為

(13)

2 人工蜂群算法

2.1 基本人工蜂群算法

蜜蜂為群居動物，在群體中不同的蜜蜂有不同的角色和任務(wù)。單個蜜蜂行為簡單，而群體運動時能夠適應(yīng)環(huán)境的改變，通過交流蜜源信息來尋求最優(yōu)蜜源。人工蜂群算法的搜索模型包括蜜源、雇傭蜂、觀察蜂和偵察蜂，不同的蜜蜂以不同的方式尋找蜜源，在采蜜過程中進行交流交換信息，最后得到最優(yōu)蜜源。

算法的具體步驟如下：

(1)初始化蜂群

將蜜源抽象為n維搜索空間內(nèi)的點，且參與尋找蜜源的蜜蜂個數(shù)為SN。則在搜索范圍內(nèi)隨機初始化的第i個蜜源為

Xi={xi,1,xi,2,…,xi,n}

(14)

xi,j=xmin,j+rand(0,1)(xmax,j-xmin,j)

(15)

其中i=1,2,…,SN；j=1,2,…,n；Xmin,j和Xmax,j分別為搜索范圍的最小值和最大值。

通過式(14)和式(15)在搜索空間內(nèi)隨機生成SN個蜜源，分別對應(yīng)SN個雇傭蜂的初始搜尋位置，且其判定值fi可以根據(jù)判定方程得出。

(2)雇傭蜂運動

每個雇傭蜂在執(zhí)行一次采蜜任務(wù)后尋找新的蜜源，單個蜜蜂會與隨機選擇的附近的蜜蜂交換蜜源信息，通過信息的交換得到新的蜜源位置。蜜源更新方程為

vi,j=xi,j+φi,j(xi,j-xk,j)

(16)

其中k∈{1,2,…,SN}，j∈{1,2,…,n}，且k與j隨機從范圍內(nèi)選取。φi,j從[-1,1] 中隨機選取。

得到新的蜜源Vi后，每個蜜蜂會與自己之前的蜜源Xi進行比較。根據(jù)判定方程的判定值fi，通過貪婪法則進行比較。如果新的蜜源Vi優(yōu)于Xi，則Vi取代Xi成為新的蜜源，如果Vi并沒有優(yōu)于Xi，則該蜜蜂的蜜源保持為Xi。

(3)選擇概率

所有的雇傭蜂在完成搜索蜜源后，在舞蹈區(qū)域與觀察蜂交換蜜源信息。每只觀察蜂的蜜源信息通過所有的雇傭蜂獲得，且蜜源的位置會根據(jù)花蜜的數(shù)量得到選擇概率。這個概率的選擇是由每個蜜源的適應(yīng)度決定，概率為

(17)

(18)

每個蜜源的適應(yīng)度fiti由判定方程的判定值得出，適應(yīng)度越高的蜜源位置越好。

(4)觀察蜂運動

根據(jù)選擇概率，每只觀察蜂從所有雇傭蜂處獲得蜜源位置信息。在選擇蜜源位置Xi后，依據(jù)式(16)更新自己的蜜源位置。運動到新的蜜源位置后，通過貪婪法則進行選擇，如果新位置優(yōu)于雇傭蜂的蜜源位置，則替代雇傭蜂進入蜂群。

(5)偵察蜂運動

偵察蜂的任務(wù)在于偵察蜜源位置的更新情況。某個雇傭蜂的蜜源位置被更換的次數(shù)超過限定值，則該雇傭蜂需要被淘汰。偵察蜂會根據(jù)式(15)在搜索范圍內(nèi)隨機產(chǎn)生一個位置，替代該雇傭蜂進入蜂群尋找蜜源。

2.2 引導人工蜂群算法

基本的算法中雇傭蜂和觀察蜂都采用式(16)來更新蜜源，即將自己蜜源信息中的一部分與鄰近蜜蜂的蜜源信息進行隨機交換。這樣交換信息更新的方式隨機性大，不易陷入局部最優(yōu)位置。然而，鄰近的蜜蜂是隨機選取，其蜜源優(yōu)劣程度未知，則在一定程度上會引導該蜜蜂轉(zhuǎn)向相對差的蜜源處。

引導人工蜂群算法(GABC)在更新蜜源位置時與基本人工蜂群算法不同，該算法選出全局最優(yōu)的蜜蜂及其蜜源位置，且隨著蜂群的運動不斷更新全局最優(yōu)。全局最優(yōu)的蜜蜂找到了此時蜂群的最佳蜜源，則有責任引導和帶領(lǐng)其余蜜蜂向好的蜜源移動。在更新蜜源時采用

vi,j=xi,j+φi,j(xi,j-xk,j)+φi,j(xGbest,j-xk,j)

(19)

式中xGbest,j——全局最優(yōu)蜜源位置

在引導人工蜂群算法中，雇傭蜂與觀察蜂均采用式(19)的方式更新蜜源位置。這樣的更新方式，一部分的信息與鄰近的蜜蜂進行交換更新，一部分的信息與全局最優(yōu)蜜蜂交換更新，即可以保證搜索的全局性，不會被局部最優(yōu)的蜜源干擾，又可以在全局最優(yōu)蜜蜂的引導下，向更好的蜜源處移動。

3 算例

根據(jù)如圖1所示的3-RPS機構(gòu)，以坐標系中心為原點，轉(zhuǎn)動副所在平臺半徑為ra=274 mm，球副所在平臺半徑為rb=158 mm。該運動平臺的工作空間如表1所示，繞X軸的轉(zhuǎn)動為α，繞Y軸的轉(zhuǎn)動為β。

表1 運動平臺工作空間
Tab.1 Workspace of moving platform

參數(shù)XYZ位移/mm00[400,700]轉(zhuǎn)角/(°)[-20,20][-20,20]0

以求解α=9°,β=13°,ze=622 mm為例分別采用連續(xù)蟻群算法、人工蜂群算法和引導人工蜂群算法進行正解運算。3-RPS并聯(lián)運動平臺的正解數(shù)值解運算為已知驅(qū)動連桿長度d1=598.849 5 mm，d2=671.584 7 mm，d3=629.495 5 mm，求解得到運動平臺的位置和姿態(tài)。本文中采用Matlab R2012b進行仿真，運行電腦配置為Intel? Core i7-4510U CPU @ 2.00 GHz 2.60 GHz RAM 7.71 GB。

采用3種算法均為群智能最優(yōu)化算法，采用仿生方式在種群的基礎(chǔ)上進行尋優(yōu)。3種算法種群數(shù)量設(shè)置為10，迭代次數(shù)為600次，獨立運行30次后得到結(jié)果進行比較。

從表2中可以看出，人工蜂群算法和引導人工蜂群算法得到的計算值相差不大，誤差都小于10-10，比較而言連續(xù)蟻群算法計算值的誤差較大。圖2為3種算法判定值誤差曲線。從圖2可以看出，連續(xù)蟻群算法曲線平緩，收斂速度慢。人工蜂群算法和引導人工蜂群算法分別在迭代170次左右和迭代90次左右收斂，達到最小值，并且引導人工蜂群算法的收斂速度最快。

表3為獨立運行30次后，判定方程誤差的平均值，標準差以及運行平均時間的比較。從表3可以看出，人工蜂群算法及引導人工蜂群算法的單次運行時間比蟻群算法長，相比較而言，引導人工蜂群算法的單次運行時間比人工蜂群算法短。但是，蟻群算法得到的誤差最大，且與人工蜂群算法和引導人工蜂群算法的誤差相差非常大，后2種算法得到的誤差非常小，優(yōu)勢明顯。

表2 各個算法實際值與計算值結(jié)果
Tab.2 Results of actual value and calculated value for each algorithm

算法α/(°)β/(°)ze/mm目標值計算值絕對誤差目標值計算值絕對誤差目標值計算值絕對誤差ACO98.859831.40165×10-11313.00181.84503×10-3622621.9995.38574×10-4ABC99.000004.40234×10-111313.00007.75024×10-12622622.0005.12727×10-10GABC99.000001.07034×10-101313.00006.02700×10-10622622.0007.27595×10-12

圖2 各個算法的誤差曲線Fig.2 Fitness error curves of each algorithm

算法平均值標準差時間/msACO5.69082×10-34.36393×10-3250.567ABC5.41289×10-172.68715×10-17362.400GABC4.10317×10-172.11369×10-17361.333

根據(jù)工程需要，在實際的運動控制中，3-RPS并聯(lián)機構(gòu)運動平臺絕對誤差達到10-4即可滿足需求。人工蜂群算法和引導人工蜂群算法都達到要求，且精度很高。實際控制要求中，精度滿足要求是最基本的條件，在運動過程中需要實時計算出機構(gòu)正解，以達到運動平臺的連續(xù)運動，因此運行時間直接影響到控制的實時性，此工程的實時控制時間要求為100 ms以下。由此，在3種算法計算精度滿足工程實際要求情況下，對運算時間進行比較，在最高迭代次數(shù)為600次的基礎(chǔ)下，設(shè)置判定方程的運算精度為10-5，得到最后的計算值的絕對誤差精度在10-4以上，具體的計算值、絕對誤差、判定誤差以及單次運算時間如表4和表5所示。從這2個表中可以看出，連續(xù)蟻群算法在迭代次數(shù)內(nèi)無法達到工程應(yīng)用要求，且運算時間也大大超出了工程應(yīng)用實時控制運算時間要求。人工蜂群算法和引導人工蜂群算法的計算值和絕對誤差都滿足工程應(yīng)用要求，且運算時間也在實時控制要求以內(nèi)。引導人工蜂群的收斂速度快，運算的時間最短，非常符合工程應(yīng)用要求。

表4 各個算法滿足精度的實際值與計算值結(jié)果
Tab.4 Results of actual value and calculated value met precision requirement for each algorithm

算法α/(°)β/(°)ze/mm目標值計算值絕對誤差目標值計算值絕對誤差目標值計算值絕對誤差ACO98.859831.40165×10-11313.00181.84503×10-3622621.9995.38575×10-4ABC99.000333.39943×10-41312.99945.18093×10-4622622.0001.60948×10-4GABC99.000121.26283×10-41313.00003.14427×10-5622622.0002.23827×10-4

表5 各個算法滿足精度的判定誤差比較
Tab.5 Comparison of fitness error met precision requirement for each algorithm

算法平均值標準差時間/msACO5.69082×10-34.36393×10-3250.567ABC5.72884×10-53.02445×10-569.800GABC5.73793×10-53.35311×10-538.300

具體判定方程的誤差值收斂情況如圖3所示，此誤差達到了10-4以下，且收斂速度非?？欤龑斯し淙簝?yōu)勢明顯。

圖3 各個算法滿足精度要求誤差曲線Fig.3 Fitness error curves met precision requirement of each algorithm

該3-RPS并聯(lián)機構(gòu)運動平臺在試驗過程中需到達工作空間內(nèi)任意位置，要求控制連續(xù)運動，實時得到機構(gòu)的運動學正解。如表6所示，列出了5種不同情況的位姿，采用引導人工蜂群算法進行求解，求解精度和求解時間均滿足工程應(yīng)用運動控制要求。

表6 不同情況下引導人工蜂群計算結(jié)果
Tab.6 Results of calculated value of GABC in different situations

序號參數(shù)真實值計算值絕對誤差適應(yīng)誤差時間/msα/(°)99.000121.26283×10-41β/(°)1313.00003.14427×10-55.73793×10-538.3000ze/mm622622.0002.23826×10-4α/(°)-17-16.99972.55598×10-42β/(°)88.000171.68791×10-45.83262×10-537.9000ze/mm436436.0002.77078×10-4α/(°)44.000141.40896×10-43β/(°)7-7.000131.31228×10-44.95234×10-533.8333ze/mm509509.0001.64332×10-4α/(°)-13-13.00001.37617×10-44β/(°)-19-18.99992.95497×10-44.76015×10-536.7667ze/mm574573.9999.32082×10-5α/(°)1818.00033.41270×10-45β/(°)1110.99963.08625×10-45.97693×10-534.6000ze/mm695694.9992.65672×10-5

4 結(jié)束語

對3-RPS并聯(lián)機構(gòu)運動平臺進行運動學分析，建立了帶約束條件下運動學正解非線性方程組，將多目標非線性方程轉(zhuǎn)換為單目標判定方程優(yōu)化。采用人工蜂群算法，通過全局最優(yōu)進行引導的方式，得到高精度低運算時間的并聯(lián)運動平臺位姿參數(shù)，達到工作空間內(nèi)控制策略的連續(xù)求解要求。實際仿真算例驗證了該算法在并聯(lián)機構(gòu)運動學正解上的有效性和優(yōu)勢。

1 孔憲文, 黃真. 3-RPS 控制位置用并聯(lián)機器人機構(gòu)的反解[J]. 機械科學與技術(shù), 1999, 18(3): 424-426. KONG Xianwen, HUANG Zhen. Inverse displacement analysis of 3-RPS regional parallel manipulators[J]. Mechanical Science and Technology, 1999, 18(3): 424-426. (in Chinese)

2 魯開講, 牛祿峰, 劉亞茹, 等. 3-RPS 并聯(lián)機構(gòu)奇異位形及工作空間研究[J]. 農(nóng)業(yè)機械學報, 2007, 38(5): 143-146. LU Kaijiang, NIU Lufeng, LIU Yaru, et al. Research on singular configuration and workspace of 3-RPS parallel mechanism[J]. Transactions of the Chinese Society for Agricultural Machinery, 2007, 38(5): 143-146. (in Chinese)

3 陳修龍, 孫先洋, 鄧昱. 4-UPS-RPS 空間5自由度并聯(lián)機構(gòu)運動學分析 [J/OL]. 農(nóng)業(yè)機械學報, 2013, 44(8): 257-261. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20130844&flag=1. DOI：10.6041/j.issn.1000-1298.2013.08.044. CHEN Xiulong, SUN Xianyang, DENG Yu. Kinematics analysis of 4-UPS-RPS spatial 5-DOF parallel mechanism[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(8): 257-261. (in Chinese)

4 GALLARDO J, OROZEO H, RICO J M. Kinematics of 3-RPS parallel manipulators by means of screw theory[J]. The International Journal of Advanced Manufacturing Technology, 2008, 36(5-6): 598-605.

5 RAO P S, RAO N M. Position analysis of spatial 3-RPS parallel manipulator[J]. International Journal of Mechanical Engineering and Robotics Research, 2013, 2(2): 80-90.

6 任文博, 顏兵兵, 殷寶麟, 等. 3-RPS 型并聯(lián)機構(gòu)運動正解的研究[J]. 佳木斯大學學報: 自然科學版, 2014, 32(3): 412-414. REN Wenbo, YAN Bingbing, YIN Baolin, et al. Research on direct kinematics for 3-RPS parallel mechanism[J]. Journal of Jiamusi University Natural Science Edition, 2014, 32(3): 412-414. (in Chinese)

7 牛祿峰. 利用特征列方法求解并聯(lián)機構(gòu)位置正解[J]. 機械設(shè)計, 2007, 24(8): 28-30. NIU Lufeng. Using the characteristic row method to solve the positional normal solution of parallel mechanism [J]. Journal of Machine Design, 2007, 24(8): 28-30. (in Chinese)

8 YANG C F, ZHENG S T, JIN J, et al.Forward kinematics analysis of a 3-PRS parallel manipulator[J]. World Academy of Science, Engineering and Technology, 2010, 17: 1264-1270.

9 謝志江, 梁歡, 宋代平. 基于連續(xù)蟻群算法的 3-RPS 并聯(lián)機構(gòu)正解[J]. 中國機械工程, 2015, 26(6): 799-803. XIE Zhijiang, LIANG Huan, SONG Daiping. Forward kinematics of 3-RPS parallel mechanism based on a continuous ant colony algorithm[J]. China Mechanical Engineering, 2015, 26(6): 799-803. (in Chinese)

10 李樹軍, 王陰, 王曉光. 3-RPS 并聯(lián)機器人機構(gòu)位置正解的桿長逼近法[J]. 東北大學學報: 自然科學版, 2001, 22(3): 205-207. LI Shujun, WANG Yin, WANG Xiaoguang. Forward position analysis of 3-RPS in-parallel manipulator using self-modified successive approximation method[J]. Journal of Northeastern University: Natural Science, 2001, 22(3): 205-207. (in Chinese)

11 韓方元, 趙丁選, 李天宇. 3-RPS 并聯(lián)機構(gòu)正解快速數(shù)值算法[J]. 農(nóng)業(yè)機械學報, 2011, 42(4): 229-233. HAN Fangyuan, ZHAO Dingxuan, LI Tianyu. A fast forward algorithm for 3-RPS parallel mechanism [J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(4): 229-233. (in Chinese)

12 陳長憶, 車林仙. 應(yīng)用粒子群算法的 3-RPS 并聯(lián)機器人機構(gòu)位置正解[J]. 現(xiàn)代制造工程, 2006(5): 77-79. CHEN Changyi, CHE Linxian. Forward position analysis of 3-RPS in-parallel manipulators based on particle swarm optimization[J]. Modern Manufacturing Engineering, 2006(5): 77-79. (in Chinese)

13 ZHANG Hongli, REN Tiantian, PAZILAI Mahemuti. Forward position solution of 3-RPS in-parallel manipulator based on particle swarm optimization[C]∥The 26th Chinese Control and Decision Conference (2014 CCDC), IEEE, 2014: 4171-4177.

14 杜義浩, 謝平, 田培濤, 等. 基于改進粒子群算法的并聯(lián)機器人運動學精度提高新方法[J]. 中國機械工程, 2012, 23(16): 1938-1942. DU Yihao, XIE Ping, TIAN Peitao, et al. A new method for improving kinematics accuracy of parallel robot based on improved PSO[J]. China Mechanical Engineering, 2012, 23(16): 1938-1942. (in Chinese)

15 房立金, 黨鵬飛. 基于量子粒子群優(yōu)化算法的機器人運動學標定方法[J]. 機械工程學報, 2016, 52(7): 23-30. FANG Lijin, DANG Pengfei. Kinematic calibration method of robots based on quantum-behaved particle swarm optimization[J]. Journal of Mechanical Engineering, 2016, 52(7): 23-30. (in Chinese)

16 RAO N M, RAO K M. Dimensional synthesis of a spatial 3-RPS parallel manipulator for a prescribed range of motion of spherical joints[J]. Mechanism and Machine Theory, 2009, 44(2): 477-486.

17 吳小勇, 謝志江, 宋代平, 等. 基于改進蟻群算法的 3-PPR 并聯(lián)機構(gòu)位置正解研究[J/OL]. 農(nóng)業(yè)機械學報, 2015, 46(7): 339-344. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20150748&flag=1. DOI：10.6041/j.issn.1000-1298.2015.07.048. WU Xiaoyong, XIE Zhijiang, SONG Daiping, et al. Forward kinematics of 3-PPR parallel mechanism based on the improved ant colony algorithm[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(7): 339-344. (in Chinese)

18 GAO W, LIU S. A modified artificial bee colony algorithm[J]. Computers & Operations Research, 2012, 39(3): 687-697.

19 JADHAV H T, ROY R. Gbest guided artificial bee colony algorithm for environmental/economic dispatch considering wind power[J]. Expert Systems with Applications, 2013, 40(16): 6385-6399.

20 ZHU G, Kwong S. Gbest-guided artificial bee colony algorithm for numerical function optimization[J]. Applied Mathematics and Computation, 2010, 217(7): 3166-3173.

21 王慧穎, 劉建軍, 王全洲. 改進的人工蜂群算法在函數(shù)優(yōu)化問題中的應(yīng)用[J]. 計算機工程與應(yīng)用, 2012, 48(19): 36-39. WANG Huiying, LIU Jianjun, WANG Quanzhou. Modified artificial bee colony algorithm for numerical function optimization[J]. Computer Engineering and Applications, 2012, 48(19)：36-39. (in Chinese)

22 畢曉君, 王艷嬌. 改進人工蜂群算法[J]. 哈爾濱工程大學學報, 2012, 33(1): 117-123. BI Xiaojun, WANG Yanjiao. A modified artificial bee colony algorithm and its application[J]. Journal of Harbin Engineering University, 2012, 33(1): 117-123. (in Chinese)

23 任子武, 王振華, 孫立寧. 基于改進人工蜂群算法的并聯(lián)機器人正運動學解[J]. 機械工程學報, 2013, 49(13): 48-55. REN Ziwu, WANG Zhenhua, SUN Lining. Forward kinematics solution for parallel manipulators based on improved artificial bee colony algorithm[J]. Journal of Mechanical Engineering, 2013, 49(13): 48-55. (in Chinese)

24 畢曉君, 王艷嬌. 加速收斂的人工蜂群算法[J]. 系統(tǒng)工程與電子技術(shù), 2011, 33(12): 2755-2761. BI Xiaojun, WANG Yanjiao. Artificial bee colony algorithm with fast convergence[J]. Systems Engineering and Electronics, 2011, 33(12): 2755-2761. (in Chinese)

25 秦全德, 程適, 李麗, 等. 人工蜂群算法研究綜述[J]. 智能系統(tǒng)學報, 2014, 9(2): 377-385. QIN Quande, CHENG Shi, LI Li, et al. Artificial bee colony algorithm: a survey[J]. CAAI Transactions on Intelligent Systems, 2014, 9(2): 377-385. (in Chinese)

Forward Kinematics Optimization of 3-RPS Parallel Manipulator Based on Global-best Artificial Bee Colony Algorithm

MAO Bingyan XIE Zhijiang WU Xiaoyong YUAN Yuejun

(StateKeyLaboratoryofMechanicalTransmission,ChongqingUniversity,Chongqing400044,China)

A fast numerical algorithm for the forward kinematics optimization of 3-RPS parallel manipulator was presented. Through numerical optimization methods, the problem of forward kinematics of parallel robot was inverted to minimization problem. According to the optimization of forward kinematics for 3-RPS parallel manipulator, the forward kinematic functions, which were the multi-objective nonlinear equations with constraint, were built based on numerical method. The fitness function was established to achieve the optimization of single objective, and the minimum optimization was solved based on the global-best artificial bee colony algorithm. The information carried by employed bees was exchanged among the neighboring bees. The position information was exchanged randomly and the neighborhood bees were chosen randomly as the standard artificial bee colony algorithm. In addition, for the global-best artificial bee colony algorithm, the global best bee which carried the best position information would guide all bees to move to the better position, and convergence fast to the best position. As the simulation result of 3-RPS forward numerical kinematics showed, compared with continuous ant colony algorithm and standard artificial bee colony algorithm, the global-best artificial bee colony algorithm was an effective way, which with high precision and high speed, to solve forward kinematics of parallel manipulators.

parallel manipulator; artificial bee colony algorithm; forward kinematics

10.6041/j.issn.1000-1298.2017.01.045

2016-05-27

2016-07-25

國家自然科學基金項目(U1530138)

毛冰滟(1991—)，女，博士生，主要從事機器人技術(shù)和優(yōu)化算法研究，E-mail: sucaogen@163.com

謝志江(1962—)，男，教授，博士生導師，主要從事機械創(chuàng)新設(shè)計、機電一體化及設(shè)備故障診斷研究，E-mail: xie@cqu.edu.cn

TH112

A

1000-1298(2017)01-0339-07