王庚祥，劉宏昭, 鄧培生

(西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院，西安 710048)

并聯(lián)機(jī)構(gòu)因承載能力強(qiáng)、剛度大、精度高等優(yōu)點(diǎn)而成為研究熱點(diǎn)。早期并聯(lián)機(jī)構(gòu)為Stewart[1]平臺(tái)。隨并聯(lián)機(jī)構(gòu)與空間機(jī)構(gòu)在航空領(lǐng)域的發(fā)展[2-3]，對(duì)精確預(yù)測(cè)系統(tǒng)動(dòng)力學(xué)行為要求愈加迫切，且相關(guān)理論研究不斷發(fā)展。如并聯(lián)機(jī)構(gòu)運(yùn)動(dòng)學(xué)研究、動(dòng)力學(xué)研究、奇異位形分析、剛度分析及性能評(píng)價(jià)[5-9]。但對(duì)并聯(lián)機(jī)構(gòu)中關(guān)節(jié)間隙研究相對(duì)較少。

運(yùn)動(dòng)副是構(gòu)成并聯(lián)機(jī)構(gòu)多剛體系統(tǒng)的重要組成部分。由于關(guān)節(jié)元素間的磨損，局部變形、瑕疵及制造誤差會(huì)導(dǎo)致運(yùn)動(dòng)副中存在間隙；而間隙為保證關(guān)節(jié)元素相對(duì)運(yùn)動(dòng)與裝配的必要條件。因此，若運(yùn)動(dòng)副中不存在間隙將導(dǎo)致機(jī)構(gòu)關(guān)節(jié)鎖死，無(wú)法裝配。關(guān)節(jié)中間隙的存在不僅改變運(yùn)動(dòng)副構(gòu)件間自由度及受力狀態(tài)，且當(dāng)關(guān)節(jié)處于分離狀態(tài)時(shí)產(chǎn)生的碰撞與沖擊力將造成磨損加劇，噪聲及振動(dòng)會(huì)導(dǎo)致機(jī)構(gòu)喪失精度，降低使用效率及壽命[9]。

關(guān)節(jié)間隙對(duì)機(jī)構(gòu)動(dòng)力學(xué)性能影響研究主要集中于平面機(jī)構(gòu)分析[10-15]。但隨空間機(jī)構(gòu)的應(yīng)用廣泛及迅速發(fā)展，對(duì)關(guān)節(jié)間隙研究轉(zhuǎn)向空間[9,16-21]。本文考慮球面副間隙，基于連續(xù)接觸模型，將球面副間隙視為剛性無(wú)質(zhì)量桿[11,17-18,21]，利用牛頓-歐拉法建立4-SPS/CU并聯(lián)機(jī)構(gòu)[22]含關(guān)節(jié)間隙的動(dòng)力學(xué)模型。為該并聯(lián)機(jī)構(gòu)作為重型作戰(zhàn)車(chē)輛上的導(dǎo)彈發(fā)射跟蹤裝置提供理論依據(jù)，考慮關(guān)節(jié)間隙可改善4-SPS/CU并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)性能，提高其控制精度。

1 4-SPS/CU并聯(lián)機(jī)構(gòu)結(jié)構(gòu)特點(diǎn)

4-SPS/CU機(jī)構(gòu)模型[22]及坐標(biāo)系建立見(jiàn)圖1。設(shè)驅(qū)動(dòng)支鏈A1B1與下平臺(tái)相連的球面副(鉸點(diǎn)B1)含有間隙。

圖1 4-SPS/CU并聯(lián)機(jī)構(gòu)簡(jiǎn)圖

2 含關(guān)節(jié)間隙并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)建模

運(yùn)動(dòng)副間隙描述方法有4種，即動(dòng)量平衡法、接觸變形模型、有限元法及連續(xù)接觸模型[9,23]。由于并聯(lián)機(jī)構(gòu)自身動(dòng)力學(xué)模型具有強(qiáng)非線(xiàn)性，可避開(kāi)運(yùn)動(dòng)副中剛度系數(shù)、阻尼系數(shù)、摩擦系數(shù)及恢復(fù)系數(shù)等物理參數(shù)簡(jiǎn)化含間隙機(jī)構(gòu)動(dòng)力學(xué)建模；由于含間隙關(guān)節(jié)元素分離、碰撞時(shí)間較短，可設(shè)運(yùn)動(dòng)副元素間始終處于連續(xù)接觸狀態(tài)。

2.1 球面副間隙模型

4-SPS/CU并聯(lián)機(jī)構(gòu)的驅(qū)動(dòng)支鏈中含球面副與移動(dòng)副，而約束從動(dòng)支鏈中含萬(wàn)向鉸與圓柱副，而本文只對(duì)B1球面副間隙量進(jìn)行研究，由于① 球面副位于該并聯(lián)機(jī)構(gòu)的驅(qū)動(dòng)支鏈上，相比約束支鏈的運(yùn)動(dòng)副，考慮球面副間隙更能體現(xiàn)關(guān)節(jié)間隙對(duì)該并聯(lián)機(jī)構(gòu)輸入、輸出參數(shù)影響；② 為簡(jiǎn)化原動(dòng)力學(xué)模型而未選擇驅(qū)動(dòng)支鏈與上平臺(tái)相連的球面副。因此本文采用連續(xù)接觸模型視球面副間隙[24]為一剛性無(wú)質(zhì)量桿-間隙桿。

圖2 球關(guān)節(jié)間隙模型

(1)

(2)

2.2 考慮球面副間隙的并聯(lián)機(jī)構(gòu)位置方程

由于4-SPS/CU并聯(lián)機(jī)構(gòu)為無(wú)過(guò)約束機(jī)構(gòu)，且理想運(yùn)動(dòng)副下SPS驅(qū)動(dòng)支鏈[22]為不受約束的6自由度剛體。當(dāng)驅(qū)動(dòng)支鏈B1鉸點(diǎn)處考慮球面副間隙時(shí)，原3個(gè)轉(zhuǎn)動(dòng)自由度的球面副被引入3個(gè)移動(dòng)自由度，即該支鏈在固定坐標(biāo)系OXYZ的螺旋表示為

(3)

式中：p1,q1為下球鉸(SPS支鏈與下平臺(tái)連接的球面副)相對(duì)坐標(biāo)系OXYZ位置；u1,v1,w1為移動(dòng)副軸線(xiàn)方向余弦；l1,m1,n1為上球鉸(SPS支鏈與下平臺(tái)連接的球面副)相對(duì)于坐標(biāo)系OXYZ位置；$1,$2,$3為下球鉸轉(zhuǎn)動(dòng)自由度螺旋表示；$4為移動(dòng)副螺旋表示；$5,$6,$7為上球鉸轉(zhuǎn)動(dòng)自由度螺旋表示；$11,$12,$13為下球鉸考慮間隙新引入3個(gè)移動(dòng)自由度的螺旋表示。

設(shè)式(3)螺旋系反螺旋為

$r=(LrMrNr;PrQrRr)

(4)

基于互易積原理建立方程：

$ij°$r=0, (i=1,2,3,…,7,j=1,2,3)

(5)

對(duì)該線(xiàn)性方程求解后得方程無(wú)解，該螺旋系無(wú)反螺旋，即該含間隙的驅(qū)動(dòng)支鏈與理想支鏈自由度保持一致。4-SPS/CU并聯(lián)機(jī)構(gòu)[22]中無(wú)論考慮任何球面副含有間隙均不改變?cè)摍C(jī)構(gòu)上平臺(tái)自由度。

(6)

式中：i=1,2,3,4。

考慮球面副間隙的4-SPS/CU并聯(lián)機(jī)構(gòu)相比式(6)只有驅(qū)動(dòng)支鏈l1的位置方程發(fā)生改變，其它支鏈位置方程保持不變。驅(qū)動(dòng)支鏈l1位置方程為

(7)

2.3 4-SPS/CU并聯(lián)機(jī)構(gòu)理想運(yùn)動(dòng)學(xué)方程

(8)

(9)

式中：

2.4 理想4-SPS/CU并聯(lián)機(jī)構(gòu)的動(dòng)力學(xué)分析

為簡(jiǎn)化4-SPS/CU并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)模型，將驅(qū)動(dòng)支鏈作整體研究，求驅(qū)動(dòng)支鏈質(zhì)心位置[25]：

(10)

式中：m1為驅(qū)動(dòng)桿質(zhì)量；m2為缸體質(zhì)量；m為驅(qū)動(dòng)支鏈質(zhì)量；lsi為驅(qū)動(dòng)桿質(zhì)心至動(dòng)鉸點(diǎn)Ai距離；lmi為缸體質(zhì)心至固定鉸點(diǎn)Bi距離；li為驅(qū)動(dòng)支鏈總長(zhǎng)(隨時(shí)間變化)；lci為驅(qū)動(dòng)支鏈質(zhì)心距固定鉸點(diǎn)Bi距離；i=1,2,3,4。

2.4.1 驅(qū)動(dòng)支鏈?zhǔn)芰Ψ治?/p>

由于SPS驅(qū)動(dòng)支鏈含局部自由度，為避免其對(duì)該驅(qū)動(dòng)支鏈動(dòng)態(tài)靜力學(xué)平衡方程結(jié)構(gòu)影響[25]。需對(duì)驅(qū)動(dòng)支鏈在固定坐標(biāo)系OXYZ與局部坐標(biāo)系Biuiviwi中方位進(jìn)行分析[25]。坐標(biāo)系OXYZ相對(duì)Biuiviwi的轉(zhuǎn)換矩陣為

R(αiZ,θiX)=[RZαi][RXθi]=

(11)

驅(qū)動(dòng)支鏈在固定坐標(biāo)系OXYZ中具體受力分析見(jiàn)文獻(xiàn)[25]。驅(qū)動(dòng)支鏈在固定坐標(biāo)系OXYZ中繞固定鉸點(diǎn)Bi作純轉(zhuǎn)動(dòng)的角速度為

(12)

ηi=pilciSi

(13)

驅(qū)動(dòng)支鏈處于動(dòng)態(tài)平衡時(shí)，力平衡方程為

(14)

在局部坐標(biāo)系Biuiviwi對(duì)驅(qū)動(dòng)支鏈固定鉸點(diǎn)Bi取矩，力矩平衡方程為

(15)

據(jù)研究結(jié)論[25]，驅(qū)動(dòng)支鏈力矩平衡方程因局部自由度需去掉繞Z軸轉(zhuǎn)動(dòng)力矩平衡方程，則式(15)去掉多余后，在固定坐標(biāo)系OXYZ下力矩平衡方程為

(16)

2.4.2 恰約束從動(dòng)支鏈?zhǔn)芰Ψ治鯷25]

(17)

對(duì)下平臺(tái)質(zhì)心O取矩，力矩平衡方程為

(18)

2.4.3 上平臺(tái)受力分析[25]

上平臺(tái)質(zhì)心OXYZ在固定坐標(biāo)系角速度ω=[ωxωyωz]T，上平臺(tái)在動(dòng)坐標(biāo)系Pxyz的各轉(zhuǎn)動(dòng)慣量為I′，則上平臺(tái)相對(duì)固定坐標(biāo)系OXYZ的轉(zhuǎn)動(dòng)慣量為

I=RI′RT

(19)

式中：

當(dāng)上平臺(tái)處于動(dòng)態(tài)平衡時(shí)，力平衡方程為

(20)

對(duì)上平臺(tái)質(zhì)心O點(diǎn)取矩，力矩平衡方程為

(21)

聯(lián)立式(14)、(16)～(18)、(20)、(21)組成(6×6-4)維方程組Γ中含(6×4+8)個(gè)未知力/力矩。方程組Γ的矩陣表示記為AX=B，A為由4-SPS/CU并聯(lián)機(jī)構(gòu)結(jié)構(gòu)參數(shù)與旋轉(zhuǎn)角組成的系數(shù)矩陣；B為由該并聯(lián)機(jī)構(gòu)重力/力矩與外載荷及系統(tǒng)慣性組成的已知矩陣；X為由該機(jī)構(gòu)所有運(yùn)動(dòng)副反力組成的未知量矩陣X=A-1B。

2.5 含球面副間隙的驅(qū)動(dòng)支鏈動(dòng)力學(xué)分析

(22)

圖3 4-SPS/CU并聯(lián)機(jī)構(gòu)間隙桿動(dòng)態(tài)靜力分析

固定坐標(biāo)系OXYZ下驅(qū)動(dòng)支鏈l1力平衡方程為

(23)

在局部坐標(biāo)系Biuiviwi對(duì)驅(qū)動(dòng)支鏈固定鉸點(diǎn)B1取矩，力矩平衡方程為

(24)

同理，據(jù)式(16)寫(xiě)出式(24)去掉多余后在固定坐標(biāo)系OXYZ下力矩平衡方程為：

(25)

考慮該并聯(lián)機(jī)構(gòu)球面副關(guān)節(jié)間隙后，該機(jī)構(gòu)驅(qū)動(dòng)力為

(26)

2.5.1 含關(guān)節(jié)間隙動(dòng)力學(xué)方程求解

(1) 不考慮球面副間隙時(shí)，方程組Γ記為AX=B。通過(guò)該方程求出4-SPS/CU并聯(lián)機(jī)構(gòu)關(guān)節(jié)約束反力/力矩X=A-1B，記關(guān)節(jié)約束反力為X1。

(27)

(28)

(3) 將(2)中結(jié)論代入式(7)、(6)中并對(duì)該并聯(lián)機(jī)構(gòu)運(yùn)動(dòng)學(xué)進(jìn)行分析，導(dǎo)出所有含關(guān)節(jié)間隙的運(yùn)動(dòng)學(xué)參數(shù)，獲得關(guān)節(jié)間隙對(duì)該并聯(lián)機(jī)構(gòu)運(yùn)動(dòng)學(xué)影響，為動(dòng)力學(xué)分析奠定基礎(chǔ)。

2.6 數(shù)值計(jì)算

4-SPS/CU并聯(lián)機(jī)構(gòu)設(shè)計(jì)參數(shù)、間隙參數(shù)見(jiàn)表1。

表1 結(jié)構(gòu)參數(shù)

設(shè)上平臺(tái)不受外力情況下，質(zhì)心P點(diǎn)運(yùn)動(dòng)規(guī)律為

(29)

圖4 B1處球面副受力分析

4-SPS/CU并聯(lián)機(jī)構(gòu)B1點(diǎn)處球面副考慮關(guān)節(jié)間隙與不考慮關(guān)節(jié)間隙時(shí)對(duì)比分析曲線(xiàn)見(jiàn)圖4。由圖4看出，時(shí)間處于0.1 s，0.95 s附近時(shí)，含間隙球面副關(guān)節(jié)反力小于等于0，說(shuō)明含間隙球面副元素發(fā)生分離；時(shí)間處于0.2 s，0.9 s附近時(shí)，含間隙球面副關(guān)節(jié)反力為理想關(guān)節(jié)副反力的2～3倍，說(shuō)明含間隙球面副元素發(fā)生碰撞；時(shí)間處在0.3～0.8 s時(shí)，考慮間隙與不考慮間隙對(duì)球面副關(guān)節(jié)反力影響不大，說(shuō)明該機(jī)構(gòu)球面副元素處于接觸狀態(tài)。因此，可利用該方法判定B1點(diǎn)處球面副元素接觸分離狀態(tài)。

關(guān)節(jié)間隙對(duì)4-SPS/CU并聯(lián)機(jī)構(gòu)驅(qū)動(dòng)支鏈長(zhǎng)度與驅(qū)動(dòng)速度影響曲線(xiàn)見(jiàn)圖5。由圖5知，關(guān)節(jié)間隙對(duì)驅(qū)動(dòng)支鏈長(zhǎng)度與驅(qū)動(dòng)速度均有影響，且對(duì)驅(qū)動(dòng)速度影響較明顯。由于關(guān)節(jié)間隙位置角為該機(jī)構(gòu)運(yùn)動(dòng)學(xué)參數(shù)函數(shù)，當(dāng)上平臺(tái)速度增加時(shí)，關(guān)節(jié)間隙對(duì)驅(qū)動(dòng)支鏈長(zhǎng)度、驅(qū)動(dòng)速度影響越顯著，圖中數(shù)值分析與理論研究相吻合。

圖5 驅(qū)動(dòng)支鏈運(yùn)動(dòng)學(xué)分析

4-SPS/CU并聯(lián)機(jī)構(gòu)考慮間隙與不考慮間隙時(shí)驅(qū)動(dòng)力對(duì)比分析曲線(xiàn)見(jiàn)圖6。由圖6看出，以無(wú)質(zhì)量桿代替B1處球面副關(guān)節(jié)間隙，對(duì)該機(jī)構(gòu)驅(qū)動(dòng)力無(wú)明顯影響。原因?yàn)楸疚囊暻蛎娓遍g隙為無(wú)質(zhì)量桿避免因間隙導(dǎo)致的接觸碰撞力，忽略間隙接觸體剛度系數(shù)、阻尼系數(shù)、摩擦系數(shù)及碰撞過(guò)程中能量損失。因此，球面副間隙對(duì)該并聯(lián)機(jī)構(gòu)驅(qū)動(dòng)力影響較弱。

3 結(jié) 論

(1) 以無(wú)質(zhì)量桿代替4-SPS/CU并聯(lián)機(jī)構(gòu)中球面副間隙。建立球面副間隙模型及該機(jī)構(gòu)含間隙鉸時(shí)的位置方程。令間隙量為零研究理想機(jī)構(gòu)運(yùn)動(dòng)學(xué)，并基于達(dá)朗貝爾原理建立該機(jī)構(gòu)動(dòng)態(tài)靜力學(xué)方程，求出該機(jī)構(gòu)所有關(guān)節(jié)副反力。設(shè)含間隙運(yùn)動(dòng)副反力位置角與無(wú)間隙運(yùn)動(dòng)副反力位置角近似相等，對(duì)4-SPS/CU并聯(lián)機(jī)構(gòu)中無(wú)質(zhì)量桿進(jìn)行運(yùn)動(dòng)學(xué)分析。

(2) 通過(guò)對(duì)含關(guān)節(jié)間隙驅(qū)動(dòng)支鏈重新受力分析，寫(xiě)出該支鏈含球面副間隙動(dòng)態(tài)平衡方程，結(jié)合理想機(jī)構(gòu)動(dòng)態(tài)平衡方程組裝考慮關(guān)節(jié)間隙時(shí)的動(dòng)力學(xué)方程HX′=K。給出該方程的求解過(guò)程。通過(guò)數(shù)值方法對(duì)比分析該機(jī)構(gòu)在考慮間隙與不考慮間隙時(shí)驅(qū)動(dòng)力與關(guān)節(jié)反力及利用該機(jī)構(gòu)含間隙位置方程研究其運(yùn)動(dòng)學(xué)反解，獲得連續(xù)接觸模型中關(guān)節(jié)間隙對(duì)該機(jī)構(gòu)的驅(qū)動(dòng)力影響不明顯，而對(duì)該機(jī)構(gòu)運(yùn)動(dòng)學(xué)特性影響顯著。

[1]Gough V E, Mech E M I, Whitehall S G,et al.Universal tyre test machine[C]. in: Proceedings of Ninth International Technical Congress,F.I.S.I.T.A.,1962:117-137.

[2]Zhao Yang, Bai Zheng-feng.Dynamics analysis of space robot manipulator with joint clearance[J].ActaAstronautica,2011,68(7-8): 1147-1155.

[3]Li Tuan-jie, Guo Jian, Cao Yu-yan.Dynamic characteristics analysis of deployable space structures considering joint clearance[J].Acta Astronautica,2011,68(7/8):974-983.

[4]Jaime G A ,Agust n R A ,H ctor R G ,et al.Kinematics of an asymmetrical three-legged parallel manipulator by means of the screw theory[J].Mechanism and Machine Theory,2010,45(7):1013-1023.

[5]Staicu S, Zhang D.A novel dynamic modelling approach for parallel mechanisms analysis[J].Robotics and Computer-Integrated Manufacturing,2008,24(1):167-172.

[6]Amine S, Masouleh M T,Caro S,et al.Singularity analysis of 3T2R parallel mechanisms using grassmann-cayley algebra and grassmann geometry[J].Mechanism and Machine Theory,2012,52:326-340.

[7]Pashkevich A, Chablat D, Wenger P.Stiffness analysis of overconstrained parallel manipulators[J].Mechanism and Machine Theory,2009,44(5):966-982.

[8]Wang Jin-song, Wu Chao, Liu Xin-jun.Performance evaluation of parallel manipulators: motion/force transmissibility and its index[J].Mechanism and Machine Theory,2010,45(10):1462-1476.

[9]Shiau T N, Tsai Y J, Tsai M S.Nonlinear dynamic analysis of a parallel mechanism with consideration of joint effects[J].Mechanism and Machine Theory,2008,43(4):491-505.

[10]kErkaya S, Uzmay I.Determining link parameters using genetic algorithm in mechanisms with joint clearance[J]. Mechanism and Machine Theory,2009,44(1):222-234.

[11]Tsai M J,Lai T H.Kinematic sensitivity analysis of linkage with joint clearance based on transmission quality[J].Mechanism and Machine Theory,2004,39(11):1189-1206.

[12]Flores P, Ambrosio J.Revolute joints with clearance in multibod y systems[J].Computers and Structures, 2004, 82(17/19): 1359-1369.

[13]Wang Jin-ge, Zhang Jun-fu, Du Xiao-ping.Hybrid dimension reduction for mechanism reliability analysis with random joint clearances[J].Mechanism and Machine Theory,2011,46(10):1396-1410.

[14]Muvengei O, Kihiu J, Ikua B.Numerical study of parametric effects on the dynamic response of planar multi-bodysystems with differently located frictionless revolute clearance joints[J].Mechanism and Machine Theory,2012,53:30-49.

[15]Bai Zheng-feng, Zhao Yang.Dynamicbehaviour analysis of planar mechanical systems with clearance in revolute joints using a new hybrid contact force model[J].International Journal of Mechanical Sciences,2012,54:190-205.

[16]Huang T, Chetwynd D G, Mei J P,et al.Tolerance design of a 2-DOF overconstrained translational parallel robot[J].IEEE Transactions on Robotics, 2006,22(1):167-172.

[17]Parenti-Castelli V, Venanzi S. Clearance influence analysis on mechanisms[J].Mechanism and Machine Theory,2005, 40(12):1316-1329.

[18]Chebbi A H, Affi Z, Romdhane L.Prediction of the pose errors produced by joints clearance for a 3-UPU parallel robot[J]. Mechanism and Machine Theory,2009,44(9):1768-1783.

[19]Frisoli A, Solazzi M, Pellegrinetti D,et al.A new screw theory method for the estimation of position accuracy in spatial parallel manipulators with revolute joint clearances[J].Mechanism and Machine Theory,2011,46(12):1929-1949.

[20]Chaker A, Mlika A, Laribi M A,et al.Clearance and manufacturing errors' effects on the accuracy of the 3-RCC spherical parallel manipulator[J].European Journal of Mechanics A/Solids,2013,37:86-95.

[21]Bu Wang-hui, Liu Zhen-yu, Tan Jian-rong, et al.Detachment avoidance of joint elements of a robotic manipulator with clearances based on trajectory planning[J].Mechanism and Machine Theory,2010,45(6):925-940.

[22]王庚祥,原大寧,劉宏昭.空間4-SPS/CU并聯(lián)機(jī)構(gòu)運(yùn)動(dòng)學(xué)分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2012,43(3):207-212.

Wang Geng-xiang,Yuan Da-ning,Liu Hong-zhao. Kinematics analysis of spatial 4-SPS/CU parallel mechanism[J].Transactions of the Chinese Society for Agricultural Machinery,2012,43(3):207-212.

[23]閻紹澤,陳鹿民,季林紅,等.含間隙鉸的機(jī)械多體系統(tǒng)動(dòng)力學(xué)模型[J].振動(dòng)工程學(xué)報(bào),2003,26(3):290-294.

YAN Shao-ze,CHEN Lu-min,JI Lin-hong,et al. Dynamie modeling of multibody system with clearanceJoints[J]. Journal of Vibration Engineering, 2003,16(3):290-294.

[24]Tian Qiang, Zhang Yun-qing, Chen Li-ping, et al.Dynamics of spatial flexible multibody systems with clearance and lubricated spherical joints[J].Computers and Structures,2009,87(13-14):913-929.

[25]王庚祥,劉宏昭,原大寧.空間4-SPS/CU并聯(lián)機(jī)構(gòu)的受力分析[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(22):30-38.

WANG Geng-xiang, LIU Hong-zhao, YUAN Da-ning. Force analysis of spatial 4-SPS/CU parallel mechanism[J].Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(22): 30-38.