可重構(gòu)單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)連桿曲線與姿態(tài)的數(shù)值解法

李庠，李瑞琴，李輝，寧峰平

(中北大學(xué) 機(jī)械工程學(xué)院，山西 太原 030051)

0 引言

目前，在多工件多位姿的搬運(yùn)、裝配、包裝等領(lǐng)域主要采用多驅(qū)動(dòng)的并聯(lián)機(jī)構(gòu)[1]、串聯(lián)機(jī)構(gòu)[2]或者連桿機(jī)構(gòu)[3]來(lái)實(shí)現(xiàn)。多驅(qū)動(dòng)的并聯(lián)機(jī)構(gòu)由于其正解復(fù)雜，導(dǎo)致其控制系統(tǒng)復(fù)雜，而單驅(qū)動(dòng)能夠簡(jiǎn)化控制系統(tǒng)，使機(jī)構(gòu)擺脫對(duì)傳感器的依賴，顯著降低連續(xù)驅(qū)動(dòng)時(shí)的能耗。因此，實(shí)現(xiàn)多個(gè)精確位姿的單驅(qū)動(dòng)平面并聯(lián)機(jī)構(gòu)在工程上有重要的應(yīng)用價(jià)值。

3-RRR平面并聯(lián)機(jī)構(gòu)是典型的平面并聯(lián)機(jī)構(gòu)，相關(guān)文獻(xiàn)研究了其運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)性能[4-8]。Li等[9]研究了3-RRR機(jī)構(gòu)的曲柄存在條件，揭示了該機(jī)構(gòu)的可轉(zhuǎn)動(dòng)性。Zhang等[10]采用轉(zhuǎn)換驅(qū)動(dòng)模式的方法，在不改變其拓?fù)浣Y(jié)構(gòu)、不增加冗余的情況下消除平面3-RRR機(jī)構(gòu)的Ⅱ型奇異性。Quintero-Riaza等[11]提出了一種平面并聯(lián)機(jī)構(gòu)的尺寸綜合方法并應(yīng)用于3-RRR機(jī)構(gòu)，能夠優(yōu)化得到預(yù)期的靈巧度好、傳動(dòng)力大、覆蓋特定工作區(qū)間的機(jī)構(gòu)。

有些學(xué)者將平面并聯(lián)機(jī)構(gòu)約束為連桿機(jī)構(gòu)，進(jìn)行機(jī)構(gòu)綜合分析。Soh等[12]提出由一個(gè)平面并聯(lián)機(jī)構(gòu)約束兩個(gè)R-R(轉(zhuǎn)動(dòng)- 轉(zhuǎn)動(dòng))的曲柄綜合方法，設(shè)計(jì)出能夠?qū)崿F(xiàn)5個(gè)指定位姿的8連桿機(jī)構(gòu)。Chen等[13]通過(guò)一個(gè)連桿連接兩個(gè)4連桿，生成一種能夠?qū)崿F(xiàn)11個(gè)指定位姿的8連桿機(jī)構(gòu)。Chung[14]提出將兩個(gè)4桿機(jī)構(gòu)串聯(lián)起來(lái)轉(zhuǎn)化為可以實(shí)現(xiàn)17個(gè)指定位姿的連桿機(jī)構(gòu)，但該機(jī)構(gòu)具有兩個(gè)自由度。Surez-Velsquez等[15]提出了一種1自由度的連桿機(jī)構(gòu)，可以實(shí)現(xiàn)8個(gè)指定的位姿。

Bai等[16]提出的一種通過(guò)平行四邊形機(jī)構(gòu)約束3-RRR平面并聯(lián)機(jī)構(gòu)得到1自由度平面并聯(lián)機(jī)構(gòu)的連桿綜合方法，計(jì)算得到動(dòng)平臺(tái)能夠?qū)崿F(xiàn)10個(gè)位姿。但該文獻(xiàn)沒(méi)有對(duì)輸入- 輸出(IO)方程進(jìn)行求解，必須借助計(jì)算機(jī)輔助設(shè)計(jì)(CAD)模型仿真結(jié)果確定方程的解屬于機(jī)構(gòu)的何種構(gòu)型，因此還需要進(jìn)一步研究機(jī)構(gòu)運(yùn)動(dòng)軌跡與姿態(tài)的求解方法。

在此基礎(chǔ)上，本文研究單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)的IO特性，提出一種IO方程的數(shù)值解法，利用此方法精確生成機(jī)構(gòu)的連桿曲線(即運(yùn)動(dòng)軌跡)和姿態(tài)，能夠判斷在某一輸入下得到的多個(gè)解分別屬于機(jī)構(gòu)的何種構(gòu)型，擺脫構(gòu)型分析對(duì)CAD模型的依賴；提出一種單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)的重構(gòu)方法，計(jì)算機(jī)構(gòu)重構(gòu)后的連桿曲線與姿態(tài)。

1 3-RRR平面并聯(lián)機(jī)構(gòu)的數(shù)學(xué)表示

傳統(tǒng)3自由度3-RRR平面并聯(lián)機(jī)構(gòu)如圖1(a)所示，需要3個(gè)原動(dòng)件同時(shí)驅(qū)動(dòng)，在通過(guò)指定多個(gè)位姿的連續(xù)運(yùn)動(dòng)中，需要傳感器實(shí)時(shí)檢測(cè)機(jī)構(gòu)的運(yùn)動(dòng)位置，驅(qū)動(dòng)電機(jī)還需要頻繁地進(jìn)行正/反轉(zhuǎn)的切換，經(jīng)常進(jìn)行加速、減速運(yùn)動(dòng)，能耗較大，控制系統(tǒng)復(fù)雜。為改善以上問(wèn)題，本文在實(shí)驗(yàn)室已有3自由度3-RRR平面并聯(lián)機(jī)構(gòu)的基礎(chǔ)上，提出一種單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)，在實(shí)現(xiàn)指定多個(gè)位姿的同時(shí)，降低軟件與硬件成本和能耗。

單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)如圖1(b)所示。圖1(b)中，Ai、Bi、Ci(i=1，2，3)表示第i個(gè)支鏈的各轉(zhuǎn)動(dòng)關(guān)節(jié)，G、H為支鏈2約束支鏈1、支鏈3的轉(zhuǎn)動(dòng)關(guān)節(jié)，l1i表示桿BiCi的長(zhǎng)度，l2i表示桿AiCi的長(zhǎng)度，e1、e2表示動(dòng)平臺(tái)A1A2、A2A3的邊長(zhǎng)，θi0表示第i個(gè)支鏈的初始角度，OBXBYB為固定在B2點(diǎn)的靜坐標(biāo)系。3-RRR平面并聯(lián)機(jī)構(gòu)有3個(gè)原動(dòng)件，但桿B1C1與B3C3分別被C1G與C3H約束，形成兩個(gè)平行四邊形連桿B1C1GB2和B3C3HB2. 由于桿B1C1、B3C3的轉(zhuǎn)角θ10、θ30分別與B2C2的轉(zhuǎn)角θ20相差θ1-2=θ10-θ20、θ3-2=θ30-θ20，3個(gè)角度可以相互表示，故機(jī)構(gòu)自由度為1. 使用θ2j作為輸入角θj，其中j表示機(jī)構(gòu)處于第j個(gè)姿態(tài)：當(dāng)機(jī)構(gòu)處于初始姿態(tài)時(shí)j=0；當(dāng)桿B2C2轉(zhuǎn)過(guò)某一角度θ時(shí)，會(huì)帶動(dòng)B1C1和B3C3同時(shí)轉(zhuǎn)過(guò)θ角，此時(shí)θj=θ20+θ.

圖1 3-RRR平面并聯(lián)機(jī)構(gòu)Fig.1 3-RRR planar parallel mechanism

為了便于描述，將靜坐標(biāo)系固定在B2點(diǎn)，動(dòng)坐標(biāo)系固定在動(dòng)平臺(tái)的A2點(diǎn)并隨之移動(dòng)，如圖2所示。圖2中，OAXAYA為固定在A2點(diǎn)的動(dòng)坐標(biāo)系，以支鏈B2C2A2為例，使用向量b2表示B2的位置(由于b2為0向量，故不在圖2中標(biāo)出)，向量a20、c20表示A2、C2的初始位置，向量a120表示初始位姿下A2指向A1的向量。當(dāng)機(jī)構(gòu)處于第j個(gè)位姿時(shí)，用向量rj和角度φj描述動(dòng)平臺(tái)的位置和姿態(tài)，向量a2j、c2j描述A2、C2的位置，a12j表示A2指向A1的向量。支鏈B1C1A1和支鏈B3C3A3同理可表示。

圖2 機(jī)構(gòu)原動(dòng)件的初始姿態(tài)與第j個(gè)姿態(tài)Fig.2 Initial orientation and jth orientation of driving link of the mechanism

假設(shè)連桿A2C2為剛體，可得

(1)

同理可得關(guān)于A1C1與A3C3的方程，故(1)式可寫為

(2)

式中：aij=rj+Qjai20；cij=bi+Rjc′i，c′i=ci0-bi；Qj、Rj為平面旋轉(zhuǎn)矩陣，Qj=Q(φj)、Rj=R(θj).

故(2)式可以寫為

(3)

(4)

式中：dij=rj+eij，eij=Qjai20-bi.

(4)式簡(jiǎn)寫為

Jicosθj+Kisinθj+Li=0，

(5)

式中：Ji、Ki、Li均為系數(shù)，

(6)

(5)式即為機(jī)構(gòu)的綜合方程。

2 單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)的IO方程

建立單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)的IO方程，便于研究機(jī)構(gòu)的連桿曲線與姿態(tài)，判斷機(jī)構(gòu)的構(gòu)型情況。

圖3所示為IO分析的機(jī)構(gòu)運(yùn)動(dòng)簡(jiǎn)圖。圖3中，μ1為桿A1A2與桿A2C2的夾角，μ2為桿A2C2與桿A2A3的夾角，α為桿A1A2與桿A2A3的夾角，ψ1為桿A2C2與線C1C2的夾角，ψ2為桿A2C2與線C2C3的夾角，σ為線C1C2與線C2C3的夾角，d1、d2分別為C2點(diǎn)到C1點(diǎn)、C3點(diǎn)的距離。

圖3 IO分析的機(jī)構(gòu)運(yùn)動(dòng)簡(jiǎn)圖Fig.3 Kinematic diagram of the mechanism for IO analysis

對(duì)圖3所示機(jī)構(gòu)的任一個(gè)位姿，使用Freudenstein方程將平面連桿機(jī)構(gòu)的IO方程寫成無(wú)量綱形式：

k11-k12cosμ1-k13cosψ1+
cosψ1cosμ1-sinψ1sinμ1=0，

(7a)

k21-k22cosμ2-k23cosψ2+
cosψ2cosμ2-sinψ2sinμ2=0，

(7b)

式中：k11、k12、k13為等效4桿機(jī)構(gòu)A1C1C2A2的Freudenstein參數(shù)，

(8)

k21、k22、k23為等效4桿機(jī)構(gòu)A2C2C3A3的Freudenstein參數(shù)，

(9)

為了方便求解，需要利用輸入和輸出角度的幾何關(guān)系減少變量的個(gè)數(shù)，即

μ1=2π-μ2-α,

(10)

ψ1=2π-ψ2-σ,

(11)

式中：當(dāng)給定輸入角θj后，利用幾何關(guān)系計(jì)算可得σ.通過(guò)(10)式、(11)式消除μ1、ψ1兩個(gè)變量后，(7a)式可重寫為

k11-k12cos(μ2+α)-k13cos(ψ2+σ)+
cos(ψ2+σ)cos(μ2+α)-
sin(ψ2+σ)sin(μ2+α)=0,

(12)

(12)式簡(jiǎn)寫為

D1cosψ2+E1sinψ2+F1=0,

(13)

式中：

(14)

同理,(7b)式可寫為

D2cosψ2+E2sinψ2+F2=0,

(15)

式中：

(16)

聯(lián)立(13)式、(15)式，可得

(17)

得到最終的IO方程為

f(μ2)=(E1F2-F1E2)2+(D1F2-D2F1)2-
(D1E2-D2E1)2=0，

(18)

(18)式中只包含一個(gè)變量μ2，令t=tan(μ2/2)，利用萬(wàn)能公式

(19)

將(18)式化為一個(gè)關(guān)于變量t的6次多項(xiàng)式，此時(shí)t是可解的。得到t后，μ2和ψ2也能夠解出。當(dāng)解出μ2和ψ2后，A1、A2和A3都能用μ2和ψ2表示，

(20)

當(dāng)已知?jiǎng)悠脚_(tái)3個(gè)頂點(diǎn)中任意2個(gè)頂點(diǎn)的位姿時(shí)，便可求得rj與φj. 假設(shè)已知A1、A2兩點(diǎn)的坐標(biāo)，由圖2可知，A1、A3可表示為

a1j=rj+Q(φj)(a10-a20)+a20,

(21a)

a3=rj+Q(φj)(a30-a20)+a20,

(21b)

(21a)式和(21b)式轉(zhuǎn)置后相減，得

(a1j-a3j)T=[Q(φj)(a10-a30)]T=
[cosφj(a10-a30)+sinφjE(a10-a30)]T.

(22)

(22)式等號(hào)兩邊同時(shí)右乘E(a10-a30)，得

(23)

(22)式等號(hào)兩邊同時(shí)右乘a10-a30，得

(24)

令u=sinφj，v=cosφj，則φj可表示為

(25)

此時(shí)，rj可表示為

rj=a1j-a20-Q(φj)(a10-a20).

(26)

由于(18)式為6次多項(xiàng)式，最多具有6個(gè)實(shí)數(shù)解，其每個(gè)解都是在機(jī)構(gòu)原動(dòng)件參數(shù)不變情況下的不同構(gòu)型，即機(jī)構(gòu)除了圖1所示的構(gòu)型外，最多還存在5種不同的構(gòu)型。

3 單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)的重構(gòu)方法與IO方程的數(shù)值解法

3.1 機(jī)構(gòu)的重構(gòu)方法

機(jī)構(gòu)重構(gòu)簡(jiǎn)圖如圖4所示。圖4中，θRE為重構(gòu)角度，是B1G1與B1C1的夾角，θRE0為重構(gòu)后的初始重構(gòu)角。在桿B1C1的C1點(diǎn)固定一個(gè)圓心為B1、半徑為B1C1的圓弧形導(dǎo)軌，使G1點(diǎn)能夠沿著弧形導(dǎo)軌滑動(dòng)θRE后固定，以此改變機(jī)構(gòu)的初始重構(gòu)角θRE0，此時(shí)θRE0=θ10+θRE，約束桿B1C1的平行四邊形變?yōu)锽1G1G2B2.θRE每改變一次，則實(shí)現(xiàn)對(duì)機(jī)構(gòu)的一次重構(gòu)，每次重構(gòu)后均可得到一個(gè)新的機(jī)構(gòu)。在不同原動(dòng)件上增加弧形導(dǎo)軌，均能夠?qū)崿F(xiàn)對(duì)機(jī)構(gòu)的重構(gòu)，通過(guò)此類方法可以得到若干個(gè)新的機(jī)構(gòu)。不失一般性，在此以桿B1C1為例，研究機(jī)構(gòu)重構(gòu)后連桿曲線的求解方法。

圖4 重構(gòu)的機(jī)構(gòu)運(yùn)動(dòng)簡(jiǎn)圖Fig.4 Kinematic diagram of reconfigured mechanism

3.2 機(jī)構(gòu)IO方程的數(shù)值解法與重構(gòu)后IO方程的數(shù)值解法

由于建立IO方程時(shí)使用了萬(wàn)能公式，采用解析法求解IO方程各個(gè)參數(shù)時(shí)需要進(jìn)行角度判定，從而增加了解方程的工作量。并且由于(18)式為6次方程，可能存在多個(gè)實(shí)數(shù)解，但無(wú)法確定在不同輸入角θj下每個(gè)實(shí)數(shù)解之間的對(duì)應(yīng)關(guān)系，導(dǎo)致難以求解機(jī)構(gòu)的IO方程，因此本文提出一種IO方程的數(shù)值解法。

機(jī)構(gòu)的初始位置如圖1(b)所示，以點(diǎn)B2為原點(diǎn)，假設(shè)機(jī)構(gòu)各設(shè)計(jì)參數(shù)與運(yùn)動(dòng)參數(shù)已知。計(jì)算(8)式、(9)式，將所得的Freudenstein參數(shù)代入(18)式、(19)式求解t.通過(guò)(19)式計(jì)算得到μ2的所有實(shí)數(shù)解，即得到機(jī)構(gòu)的所有構(gòu)型。利用(17)式解出ψ2，代入(20)式計(jì)算A1、A2、A3在初始相位角θ10時(shí)的值，此時(shí)rj與φj均為0，將其作為MATLAB軟件中fsovle函數(shù)的初值代入(5)式中，計(jì)算當(dāng)θ10轉(zhuǎn)過(guò)一個(gè)微小角度Δθ后θj的rj與φj，并得到此時(shí)A1、A2、A3的值。再將此時(shí)的rj與φj作為fsovle函數(shù)的初值代入(5)式中，計(jì)算θj轉(zhuǎn)過(guò)一個(gè)微小角度Δθ得到新的輸入角θj時(shí)機(jī)構(gòu)的rj與φj，得到此時(shí)A1、A2、A3的值，以此循環(huán)，直到θj=360°時(shí)，即可獲得機(jī)構(gòu)的連桿曲線。

若已知機(jī)構(gòu)某一位姿，則可將該位姿代入不同構(gòu)型的方程中(參數(shù)t、μ2、ψ2以及動(dòng)平臺(tái)初始位姿不同)，檢查計(jì)算是否收斂，以此判斷該位姿屬于何種構(gòu)型，也可判斷不同位姿下t的多個(gè)解的關(guān)系。

當(dāng)機(jī)構(gòu)重構(gòu)時(shí)只改變了機(jī)構(gòu)的初始重構(gòu)角θRE0，故上述方法同樣適用于重構(gòu)后的機(jī)構(gòu)。在求解時(shí)，只需要改變機(jī)構(gòu)的初始重構(gòu)角，重復(fù)上述方法即可得到機(jī)構(gòu)重構(gòu)后的連桿曲線，計(jì)算流程如圖5所示。

圖5 機(jī)構(gòu)連桿曲線計(jì)算流程圖Fig.5 Flowchart for calculating coupler curve of mechanism

4 數(shù)值算例

如圖1(b)所示，以點(diǎn)B2為原點(diǎn)，機(jī)構(gòu)的設(shè)計(jì)參數(shù)如表1所示，運(yùn)動(dòng)參數(shù)如表2所示。

表1 機(jī)構(gòu)的設(shè)計(jì)參數(shù)Tab.1 Design parameters of mechanism

表2 機(jī)構(gòu)的運(yùn)動(dòng)參數(shù)Tab.2 Kinematic parameters of mechanism

通過(guò)(18)式和(19)式計(jì)算得到μ2有2個(gè)實(shí)數(shù)解，如表3所示，表明該機(jī)構(gòu)除了圖1(b)所示的構(gòu)型以外，還存在一種構(gòu)型，如圖6所示，即該機(jī)構(gòu)共有兩種構(gòu)型。

表3 IO方程的兩個(gè)解Tab.3 Two solutions of IO equation

圖6 單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)的構(gòu)型ⅡFig.6 Configuration Ⅱ of single-driven 3-RRR planar parallel mechanism

將圖1(b)所示的機(jī)構(gòu)稱為構(gòu)型Ⅰ，圖6所示的機(jī)構(gòu)稱為構(gòu)型Ⅱ，構(gòu)型Ⅰ與構(gòu)型Ⅱ的設(shè)計(jì)參數(shù)與運(yùn)動(dòng)參數(shù)完全相同，但動(dòng)平臺(tái)A1A2A3的位姿不同。構(gòu)型Ⅰ的ψ2通過(guò)計(jì)算所得的值為155.48°，但圖1(b)中ψ2明顯大于180°，因此ψ2=360°-155.48°=204.52°. 這表明若采用解析法求解IO方程，則還需解決角度判定問(wèn)題，因此本文提出的數(shù)值解法能夠很大程度上化簡(jiǎn)求解IO方程的工作。

構(gòu)型Ⅰ與構(gòu)型Ⅱ的重構(gòu)方法和連桿曲線與姿態(tài)的計(jì)算步驟相同，因此本文僅以構(gòu)型Ⅰ作為算例。按照?qǐng)D5所示計(jì)算流程計(jì)算得到構(gòu)型Ⅰ的連桿曲線與姿態(tài)φj，通過(guò)三維建模SolidWorks軟件仿真得到連桿曲線，將理論計(jì)算曲線與仿真曲線進(jìn)行比對(duì)驗(yàn)證，發(fā)現(xiàn)二者完全吻合(以B2為坐標(biāo)原點(diǎn)、A2為動(dòng)平臺(tái)原點(diǎn))，如圖7所示。構(gòu)型Ⅰ的連桿曲線在x軸上的坐標(biāo)最大值和最小值分別為50.89 mm和-85.01 mm，跨度為135.9 mm；在y軸上的坐標(biāo)最大值和最小值分別為102.1 mm和9.76 mm，跨度為111.86 mm.

圖7 構(gòu)型Ⅰ的連桿曲線Fig.7 Coupler curve of configuration Ⅰ

通過(guò)動(dòng)力學(xué)仿真ADAMS軟件對(duì)其姿態(tài)φj進(jìn)行仿真，其仿真結(jié)果與理論計(jì)算所得姿態(tài)的曲線吻合，最大誤差約為0.3°，如圖8所示。構(gòu)型Ⅰ的φj在一個(gè)轉(zhuǎn)動(dòng)周期內(nèi)的最大值和最小值分別為21.67°和-3.91°，跨度為25.58°.

圖8 構(gòu)型Ⅰ的姿態(tài)Fig.8 Orientation of configuration Ⅰ

通過(guò)仿真結(jié)果和理論計(jì)算結(jié)果的對(duì)比，表明IO方程與求解方法的正確性。本文提出方程的求解思路同樣適用于其他IO方程求解工作復(fù)雜、困難的連桿機(jī)構(gòu)。

機(jī)構(gòu)重構(gòu)的連桿曲線如圖9所示，隨著重構(gòu)角θRE的增加，機(jī)構(gòu)在x軸上的工作范圍越來(lái)越小，在y軸上的工作范圍越來(lái)越大。構(gòu)型Ⅰ在θRE=45°時(shí)，連桿曲線在x軸的坐標(biāo)最大值和最小值分別為25.07 mm和-55.4 mm，跨度為80.47 mm；y軸的坐標(biāo)最大值和最小值分別為119.6 mm和-7.53 mm，跨度為127.13 mm. 相比于原機(jī)構(gòu)，連桿曲線的工作范圍在x軸上減少了40.79%，在y軸上增加了13.65%.由此可得，重構(gòu)主要影響了構(gòu)型Ⅰ的連桿曲線在x軸上的工作范圍，且θRE每增加15°，連桿曲線在x軸上的坐標(biāo)最大值減少9 mm左右，坐標(biāo)最小值增加10 mm左右。

圖9 構(gòu)型Ⅰ重構(gòu)的連桿曲線Fig.9 Coupler curves of reconfiguration of configuration Ⅰ

由圖9(a)可知，重構(gòu)改變了構(gòu)型Ⅰ的連桿曲線，但所有重構(gòu)后的曲線都集中相交于M、N兩點(diǎn)附近，在實(shí)際應(yīng)用中可將這兩點(diǎn)作為機(jī)構(gòu)的特征點(diǎn)，再通過(guò)具體工況選擇合適的連桿曲線。

構(gòu)型Ⅰ重構(gòu)后的姿態(tài)φj如圖10所示。由圖10可見(jiàn)：隨著重構(gòu)角θRE的增加，構(gòu)型Ⅰ姿態(tài)φj的曲線極值點(diǎn)在橫坐標(biāo)上發(fā)生遷移，曲線的大致形狀相似，都存在一個(gè)波峰和一個(gè)波谷，且姿態(tài)φj的最大值和最小值也隨之增大和減小；當(dāng)θRE=0°時(shí)姿態(tài)φj的最大值和最小值分別為21.67°和-14.78°，當(dāng)θRE=45°時(shí)姿態(tài)φj的最大值和最小值分別為32.05°和3.911°，最大值和最小值的跨度分別為10.3°和18.691°.

圖10 構(gòu)型Ⅰ重構(gòu)的姿態(tài)Fig.10 Orientation of reconfiguration of configuration Ⅰ

圖10(a)中，仿真的姿態(tài)φj曲線的個(gè)別地方出現(xiàn)局部的峰值，這是因?yàn)樵诜抡孢^(guò)程中該點(diǎn)接近于機(jī)構(gòu)的奇異位置，從而導(dǎo)致機(jī)構(gòu)運(yùn)動(dòng)的突變。由于奇異位置的存在也限制了仿真過(guò)程中仿真步長(zhǎng)的設(shè)置，最終導(dǎo)致姿態(tài)φj局部峰值的出現(xiàn)，與理論結(jié)果出現(xiàn)差異，因此在之后的研究中還要對(duì)機(jī)構(gòu)的性能進(jìn)行分析與優(yōu)化。

由圖9、圖10可知，構(gòu)型Ⅰ通過(guò)理論計(jì)算得到連桿曲線與姿態(tài)φj都與仿真結(jié)果相符，表明本文所提IO方程數(shù)值求解方法的正確性，并能夠適用于重構(gòu)后機(jī)構(gòu)的計(jì)算。

5 結(jié)論

本文利用Freudenstein方程推導(dǎo)出單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)的IO方程，提出一種數(shù)值求解方法。在此基礎(chǔ)上，提出了一種通過(guò)改變某一驅(qū)動(dòng)轉(zhuǎn)動(dòng)副的初始角度實(shí)現(xiàn)重構(gòu)的方法并通過(guò)數(shù)算例求解驗(yàn)證。得到以下主要結(jié)論：

1)提出一種IO方程的數(shù)值解法，得到機(jī)構(gòu)的所有構(gòu)型及其連桿曲線與姿態(tài)，能夠通過(guò)計(jì)算判斷某一位姿屬于機(jī)構(gòu)的何種構(gòu)型。

2)通過(guò)對(duì)比理論與仿真得到的連桿曲線與姿態(tài)φj的曲線，證明該數(shù)值求解方法的正確性。該方法化簡(jiǎn)之后的構(gòu)型分析工作，可為其他求解IO方程工作復(fù)雜、困難的連桿機(jī)構(gòu)提供參考。

3)使用提出的求解方法計(jì)算了該可重構(gòu)機(jī)構(gòu)的運(yùn)動(dòng)軌跡，理論計(jì)算結(jié)果與仿真軌跡相符，證明了該方法同樣適用于重構(gòu)后的單驅(qū)動(dòng)3-RRR平面并聯(lián)機(jī)構(gòu)。

4)提出的約束方法能為其他多驅(qū)動(dòng)平面并聯(lián)機(jī)構(gòu)轉(zhuǎn)化為單驅(qū)動(dòng)或少驅(qū)動(dòng)機(jī)構(gòu)提供研究基礎(chǔ)；采用的重構(gòu)方法可為單驅(qū)動(dòng)機(jī)構(gòu)改變運(yùn)動(dòng)軌跡的重構(gòu)設(shè)計(jì)提供參考。

參考文獻(xiàn)(References)

[1] 李建鋼，楊桂林，張林，等.基于柔性并聯(lián)機(jī)構(gòu)的變剛度夾持器設(shè)計(jì)分析[J].機(jī)械工程學(xué)報(bào)，2019，55(21):64-72.

LI J G,YANG G L,ZHANG L,et al.Design analysis of a variable stiffness gripper based on flexible parallel mechanism[J].Journal of Mechanical Engineering,2019,55(21):64-72.(in Chinese)

[2] 鄭昱，廣晨漢，李前，等.一種新型物料裝填機(jī)器人的設(shè)計(jì)與軌跡優(yōu)化[J].兵工學(xué)報(bào)，2020，41(4):763-770.

ZHENG Y,GUANG C H,LI Q,et al.Design and trajectory optimization of an object loading robot[J].Acta Armamentarii,2020,41(4):763-770.(in Chinese)

[3] 彭京徽，王德石，張弘弨.艦炮含鉗爪回轉(zhuǎn)機(jī)構(gòu)的伸爪運(yùn)動(dòng)特性研究[J].兵工學(xué)報(bào)，2020，41(5):874-880.

PENG J H,WANG D S,ZHANG H C.Research on the motion characteristics of extendable gripper of rotary conveyor mechanism in naval gun[J].Acta Armamentarii,2020,41(5):874-880.(in Chinese)

[4] CAN E,STACHEL H.A planar parallel 3-RRR robot with synchronously driven cranks[J].Mechanism and Machine Theory,2014,79:29-45.

[6] LIU S,QIU Z C,ZHANG X M.Singularity and path-planning with the working mode conversion of a 3-DOF 3-RRR planar parallel manipulator[J].Mechanism and Machine Theory,2017,107:166-182.

[7] WU X Y,BAI S P.Analytical determination of shape singularities for three types of parallel manipulators[J].Mechanism and Machine Theory,2020,149:103812.

[8] ZHANG X C,ZHANG X M.A comparative study of planar 3-RRR and 4-RRR mechanisms with joint clearances[J].Robotics and Computer-Integrated Manufacturing,2016,40:24-33.

[9] LI R Q,DAI J S.Crank conditions and rotatability of 3-RRR planar parallel mechanisms[J].Science China-Technological Sciences,2009,52(12):3601-3612.

[10] ZHANG Z K,WANG L P,SHAO Z F.Improving the kinematic performance of a planar 3-RRR parallel manipulator through actuation mode conversion[J].Mechanism and Machine Theory,2018,130:86-108.

[12] SOH G S,YING F T.Dimensional synthesis of planar eight-bar linkages based on a parallel robot with a prismatic base joint[C]∥Proceedings of the ASME 2013 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference.Portland,OR,US:ASME,2013.

[13] CHEN C,ANGELES J.A novel family of linkages for advanced motion synthesis[J].Mechanism and Machine Theory,2008,43(7):882-890.

[14] CHUNG W Y.Synthesis of two four-bar in series for body gui-dance[J].MATEC Web of Conferences,2017,95:04006.

[16] BAI S P，LI Z Y，LI R Q.Exact synthesis and input-output analysis of 1-DOF planar linkages for visiting 10 poses[J].Mechanism and Machine Theory,2020,143:103625.