3-RRP平面并聯(lián)機構(gòu)的工作空間和靈巧度分析

劉小娟，李 虹，李瑞琴，郝亮亮，汪 輝

(中北大學(xué) 機械與動力工程學(xué)院，山西 太原 030051)

3-RRP平面并聯(lián)機構(gòu)的工作空間和靈巧度分析

劉小娟，李 虹，李瑞琴，郝亮亮，汪 輝

(中北大學(xué) 機械與動力工程學(xué)院，山西 太原 030051)

建立了3-RRP平面并聯(lián)機構(gòu)的正反解方程，利用Matlab/Simulink中的SimMechanics功能模塊建立3-RRP平面并聯(lián)機構(gòu)以及其結(jié)構(gòu)框圖.采用滑桿的行程限制條件編程得出了3-RRP平面并聯(lián)機構(gòu)的工作空間范圍.利用條件數(shù)轉(zhuǎn)化到其倒數(shù)來對3-RRP并聯(lián)機構(gòu)進(jìn)行靈巧度分析，得出了該機構(gòu)的靈巧度的分布等高線，靈巧度在整個工作空間分布均勻，驗證了該機構(gòu)的運動學(xué)處于各向同性，各項性能也較好，研究結(jié)果為3-RRP平面并聯(lián)機構(gòu)優(yōu)化設(shè)計和動力學(xué)研究等奠定理論基礎(chǔ).

Matlab；并聯(lián)機構(gòu)；工作空間；靈巧度

三自由度平面并聯(lián)機構(gòu)是平面并聯(lián)機構(gòu)研究較多的類型之一[1-2]，近年來吸引了眾多學(xué)者對其進(jìn)行研究與開發(fā).文獻(xiàn)[3]運用BP神經(jīng)網(wǎng)絡(luò)建立適合3-RRP平面并聯(lián)機構(gòu)的神經(jīng)網(wǎng)絡(luò)模型，并求解出該機構(gòu)的位置正解.文獻(xiàn)[4]用平面3-PPR并聯(lián)機構(gòu)的位置逆解以及該機構(gòu)的桿長等限制因素，采用Matlab編程得出其工作空間.文獻(xiàn)[5]提出3-PRR并聯(lián)機構(gòu)獲得最大區(qū)域沒有奇異點的圖形算法.文獻(xiàn)[6]利用可操作數(shù)和條件數(shù)對3-RRR并聯(lián)機構(gòu)進(jìn)行了靈巧度分析，得出該機構(gòu)的靈巧性特點.文獻(xiàn)[7]用幾何方法建立平面三自由度并聯(lián)機構(gòu)工作空間與結(jié)構(gòu)參數(shù)的關(guān)系，提出機構(gòu)滿足運動可能性的條件.

目前，國內(nèi)外對3-RRP并聯(lián)機構(gòu)的運動學(xué)分析、動力學(xué)及仿真有一定的研究，而工作空間[8]和靈巧度是并聯(lián)機構(gòu)性能的重要指標(biāo)，機構(gòu)的靈巧度是描述一個機構(gòu)運動失真程度及遠(yuǎn)離奇異點的程度，需對該機構(gòu)靈巧度進(jìn)行分析研究，以改善機構(gòu)的運動性能.本文運用Matlab/Simulink中的SimMec-hanics功能模塊建立3-RRP并聯(lián)機構(gòu)，對3-RRP并聯(lián)機構(gòu)的工作空間進(jìn)行了仿真，以利用條件數(shù)對3-RRP并聯(lián)機構(gòu)進(jìn)行了靈巧度分析，為此構(gòu)型的特性研究及優(yōu)化設(shè)計提供理論基礎(chǔ).

1 3-RRP并聯(lián)機構(gòu)描述及位置分析

3-RRP平面并聯(lián)機構(gòu)的機構(gòu)簡圖如圖1所示. 該機構(gòu)可以實現(xiàn)2個平移和1個轉(zhuǎn)動，由靜平臺、動平臺和3個完全對稱分布的運動支鏈組成，每條支鏈都是由2個轉(zhuǎn)動副和1個移動副組成，移動副為末端支鏈與動平臺相連.

圖1 3-RRP并聯(lián)機構(gòu)簡圖Fig.1 Diagram of 3-RRP parallel mechanism

如圖1所示，建立局部坐標(biāo)系Bi-xiyi， △B1B2B3為正三角形，建立動坐標(biāo)系p-xy，動坐標(biāo)系固接在動平臺上，隨著動平臺的運動而運動，中心點p處坐標(biāo)值為(xp,yp,αp)，其中，xp和yp分別為動坐標(biāo)系在定坐標(biāo)系中沿x軸和y軸方向上的位置坐標(biāo)，αp為動坐標(biāo)系在定坐標(biāo)系的轉(zhuǎn)動角度.設(shè)定轉(zhuǎn)動副的初始轉(zhuǎn)角為α1,α2,α3，3條支鏈中驅(qū)動桿BiCi的長度為H，連接桿CiDi的長度為L，三角形動平臺外接圓半徑為e，定平臺B1B2B3的外接圓半徑為R.

1.1 求解3-RRP并聯(lián)機構(gòu)位置正解

已知各驅(qū)動輸入?yún)?shù)(θ1,θ2,θ3)的大小，求解動平臺輸出位姿參數(shù)(xp,yp,αp).建立如式(1)和(2)所示的位移方程：

H[cos(αi+αp)-cosαi]+L[sinαi-sin(θi+
αi)]+esinαi=xp+cisin(αi+αp)

(1)

H[sinαi-sin(αi+αp)]-L[cos(θi+ai)-
cosαi]+ecosαi=yp-cicos(αi+αp)

(2)

消去式中ci(ci為位移變量)，整理化簡為如式(3)所示的運動方程：

Eixp+Fiyp+Gi=0

(3)

式中：Ei=cos (αi+αp)，F(xiàn)i=sin (ai+αp),Gi=H(cosαp-1)+Lsin (θi-αp)+(L+e)sinαp.

根據(jù)機構(gòu)的幾何關(guān)系：

E1+E2+E3=0

(4)

F1+F2+F3=0

(5)

G1+G2+G3=0

(6)

將式(4)～(6)代入式(3)得到：

(7)

式中：u=3(L+e)-L(cosθ1+cosθ2+cosθ3),v=3H+L(sinθ1+sinθ2+sinθ3),w=3H.

求得的位置角αp代入式(1)和(2)中，求解得：

xp=sinαp[2Lsin(θ2-αp)+(L+e)sinαp+
H(cosαp-1)-Lsin(θ1-αp)]-
Lsin(θ1-αp)cosαp-(L+e)sinαp-
H(cosαp-1)

(8)

yp=cosαp[-2Lsin(θ2-αp)-3(L+e)sinαp-
3H(cosαp-1)-Lsin(θ1-αp)]-
Lsin(θ1-αp)sinαp-(L+e)sinαp-
H(cosαp-1)

(9)

式(7)～(9)聯(lián)立就可得出3-RRP平面并聯(lián)機構(gòu)的位置正解.

1.2 求解3-RRP并聯(lián)機構(gòu)位置反解

求解并聯(lián)機構(gòu)的位置逆解，就是已知該機構(gòu)動平臺運動位置參數(shù)(xp,yp,αp)，求解3個連桿的轉(zhuǎn)動輸入角參數(shù)(θ1,θ2,θ3).根據(jù)式(1)和(2)化簡整理后得到運動學(xué)反解方程為

Risinθi+Sicosθi-Ti=0

(10)

式中：Ri=Lcosαp,Si=sinαp,Ti=H(1-cosαp)-cos (αp+αi)xp-sin (αp+αi)yp-(L+e)sinαp.

化簡可得：

(11)

2 3-RRP并聯(lián)機構(gòu)工作空間分析

已知3-RRP平面并聯(lián)機構(gòu)參數(shù)：L=160 mm，H=80 mm，R=316 mm.本文求解圖1所示的3-RRP并聯(lián)機構(gòu)工作空間：即在滿足桿長尺寸范圍的條件下，求解動平臺中心點p可達(dá)空間點的集合.

圖2 3-RRP并聯(lián)機構(gòu)SimMechanics仿真模型Fig.2 SimMechanics simulation model of 3-RRP parallel mechanism

2.1 構(gòu)建SimMechanics模型

在Matlab/Simulink工具箱建立并聯(lián)機構(gòu)模擬框圖的步驟：(1) 確定剛體的慣性量、自由度、約束；(2) 建立檢測和驅(qū)動模塊組；(3) 以系統(tǒng)的運動形式和驅(qū)動力開始模擬；(4) 在SimMechanics窗口顯示模擬結(jié)果.

根據(jù)3-RRP平面并聯(lián)機構(gòu)的參數(shù)，結(jié)合正反解的方程，建立平面3-RRP并聯(lián)機構(gòu)的模擬框圖和機構(gòu)簡圖，如圖2和3所示.

圖3 3-RRP并聯(lián)機構(gòu)簡圖Fig.3 Diagram of 3-RRP parallel mechanism

由圖2可知，RootPart為靜定平臺，由無自由度Weld模塊連接到RootGround模塊；機構(gòu)的工作環(huán)境模塊Env設(shè)置重力加速度向量和容許誤差等；機構(gòu)的旋轉(zhuǎn)矩陣、剛體位置、速度(或角速度)和加速度(或角加速度)由剛體傳感器Body Sensor模塊測量. 機構(gòu)3個軸向上的位置數(shù)據(jù)由傳感器采集，再經(jīng)過To Work Space模塊傳遞到Matlab工作空間中并為用戶查看及使用.

圖3中的機構(gòu)簡圖與圖1中的是一致的，從而驗證了用Matlab/Simulink建模的正確性.

2.2 求解3-RRP工作空間

設(shè)定滑桿桿長L變化范圍為50～70 mm，利用Matlab可求得3-RRP的工作空間，如圖4所示.

圖4 3-RRP并聯(lián)機構(gòu)工作空間Fig.4 Working space of 3-RRP parallel mechanism

由于驅(qū)動副支鏈的慣性矩存在，使得該支鏈的運動產(chǎn)生偏離，故該3-RRP平面并聯(lián)機構(gòu)的工作空間也繞x軸產(chǎn)生了偏差Δx. 但是從圖4可看出, 該機構(gòu)的工作空間仍呈現(xiàn)對稱性，在整個工作空間內(nèi)結(jié)構(gòu)緊湊，無空洞，說明該機構(gòu)具有良好的工作性能.

3 3-RRP并聯(lián)機構(gòu)的靈巧度分析

定義條件數(shù)為

(12)

式中：‖*‖為任意矩陣的Euclide(歐幾米德)范數(shù).奇異值和條件數(shù)之間的相互關(guān)系[9]為

k(J)=σmax/σmin

(13)

式中：σmax為最大奇異值；σmin為最小奇異值.矩陣的條件數(shù)k取值范圍為1～∞，本文根據(jù)該機構(gòu)的雅克比矩陣可求得線速度和角速度的傳遞矩陣表示為Jv和Jw，求得線速度條件數(shù)k(Jv)和角速度條件數(shù)k(Jw)分別為

(14)

式中：σmaxv和σminv分別表示線速度傳遞矩陣Jv的最大奇異值和最小奇異值；σmaxw和σminw分別表示角速度傳遞矩陣Jw的最大奇異值和最小奇異值；1≤k(Jv)≤∞；1≤k(Jw)≤∞.由于k(Jv)和k(Jw)的變化范圍較大，將條件數(shù)的倒數(shù)作為速度傳遞評價指標(biāo)[10]，如式(15)所示.

(15)

式中：c為機構(gòu)的雅克比矩陣條件數(shù)的倒數(shù)．由此可知，c值越接近1，機構(gòu)的輸出偏差就越小，機構(gòu)的各項性能也較好，并且此時機構(gòu)運動學(xué)處于各向同性.

3-RRP平面并聯(lián)機構(gòu)的c分布如圖5所示.由圖5可知：c截面的變化值從外到內(nèi)逐漸減小，c最大值和最小值分別為0.8和0.3；c在整個工作空間較均勻分布且由中心向四周逐漸減小，c值在工作空間的中心處接近1，則該機構(gòu)的輸出偏差較小，此時機構(gòu)的各項性能也最佳，機構(gòu)運動學(xué)處于各向同性.

圖5 3-RRP并聯(lián)機構(gòu)的c分布圖Fig.5 c distribution of 3-RRP parallel mechanism

4 結(jié) 論

(1) 本文建立了3-RRP平面并聯(lián)機構(gòu)的正反解方程.

(2) 利用Matlab構(gòu)建3-RRP平面并聯(lián)機構(gòu)的SimMechanics模型與機構(gòu)模型一致，說明Matlab/ Simulink建模的正確性.結(jié)合正反解的方程及特定參數(shù)，經(jīng)Matlab編程求得3-RRP平面并聯(lián)機構(gòu)的工作空間，在整個工作空間內(nèi)該機構(gòu)結(jié)構(gòu)緊湊且無空洞，說明該機構(gòu)具有良好的工作性能.用本文方法求解機構(gòu)工作空間，克服了利用極限理論編程運算的復(fù)雜，對高冗余度的機構(gòu)工作空間求解優(yōu)勢更明顯.

(3) 以機構(gòu)的雅克比矩陣條件數(shù)的倒數(shù)作為速度傳遞評價指標(biāo)，繪制了3-RRP并聯(lián)機構(gòu)的c分布圖.在工作空間的中心處c的值接近1，此時機構(gòu)運動學(xué)處于各向同性，各項性能最佳，為該機構(gòu)特性研究及優(yōu)化設(shè)計提供理論基礎(chǔ).

[1] ZEIN M, WENGER P, CHABLAT D. Non-singular assembly-mode changing motions for 3-RPR parallel manipulators [J]. Mechanism and Machine Theory, 2008, 43(4) :480-490.

[2] BINAUD N, CARO S, WENGER P. Sensitivity comparison of planar parallel manipulators [J]. Mechanism and Machine Theory, 2010, 45:1477-1490.

[3] 韓霄,李虹,李瑞琴.基于BP法的3-RRP并聯(lián)機構(gòu)位置正解研究[J]. 組合機床與自動化加工技術(shù), 2015(2): 94-97.

[4] 高曉雪,梅瑛,李瑞琴. 3-PPR平面并聯(lián)機構(gòu)的工作空間分析[J]. 機械傳動, 2014, 38(2): 90-92.

[5] KALOORAZI M H F, MASOULEH M T, CARO S. Deter-mination of the maximal singularity-free workspace of 3-DOF parallel mechanisms with a constructive geometric approach[J]. Mechanism and Machine Theory, 2015, 43(4):25-36.

[6] 杜聿靜, 王雷. 3-RRR并聯(lián)機構(gòu)靈巧度分析[J]. 機械工程與自動化, 2014,183(2):75-77.

[7] 史革盟,魯開講. 并聯(lián)機構(gòu)工作空間與奇異位形分析[J]. 機械傳動, 2014(11): 54-59.

[8] HUSNG M Z, THEBERT J L. A study of workspace and singularity characteristics for design of 3-DOF planar parallel robot[J]. International Journal of Advanced Manufacturing Technology, 2010, 51(5): 789-797.

[9] 吳孟麗,張大衛(wèi),趙興玉. 一種新型非對稱并聯(lián)機構(gòu)的運動學(xué)分析[J]. 中國機械工程,2008(12):1423-1428.

[10] 賈曉輝,劉今越. 3-PRR柔性并聯(lián)機構(gòu)的優(yōu)化設(shè)計[J]. 機械設(shè)計,2014(1):26-29.

Workspace and Dexterity Analysis of 3-RRP Planar Parallel Mechanism

LIUXiao-juan，LIHong，LIRui-qin，HAOLiang-liang，WANGHui

(School of Mechanical and Power Engineering, North University of China, Taiyuan 030051, China )

The positive and inverse solution equations of the 3-RRP planar parallel mechanism are established, and the 3-RRP planar parallel mechanism and its structure block diagram are established by using the SimMechanics function module in Matlab/Simulink. The 3-RRP planar parallel mechanism workspace range is obtained by the sliding rod travel constraint programming.The dexterity of the contour line of the 3-RRP parallel mechanism is obtained by the condition number into the countdown to the dexterity analysis, the distribution of the whole workspace is uniform, and the kinematic of the mechanism is isotropic, and a variety of performance are also good. The results of the study for 3-RRP planar parallel mechanism takes as the theoretical basis for optimization design and dynamics research.

Matlab；parallel mechanism；workspaces；dexterity

1671-0444 (2016)04-0523-04

2015-12-02

國家自然科學(xué)基金資助項目(51275486)

劉小娟(1989—)，女，山西忻州人，碩士研究生，研究方向為機構(gòu)動態(tài)設(shè)計與優(yōu)化. E-mail:994512601@qq.com 李 虹(聯(lián)系人)，女，副教授，E-mail:349824199@qq.com.

TH 124

A