于榮榮 李永亮 董禮港 呂 寵 王國欣

(1.北京衛(wèi)星制造廠有限公司， 北京 100094； 2.北京空間機電研究所， 北京 100094)

0 引言

在大型航天相機的輻射定標(biāo)測試環(huán)節(jié)，需要對相機進行空間六自由度的姿態(tài)精確調(diào)整。目前，轉(zhuǎn)運多采用吊車加專用轉(zhuǎn)運工裝方式，利用激光跟蹤儀測量當(dāng)前姿態(tài)和目標(biāo)姿態(tài)，由人工對與調(diào)姿面連接的多個絲杠的高度進行反復(fù)手動調(diào)姿，直至調(diào)整至誤差允許范圍之內(nèi)。這種調(diào)姿方法工作量大、工作效率低、操作繁瑣，并且調(diào)姿精度低、產(chǎn)品風(fēng)險高。隨著航天相機功能的多樣化，其結(jié)構(gòu)越來越復(fù)雜，質(zhì)量和尺寸越來越大，傳統(tǒng)的調(diào)姿方法很難滿足超大、超重載荷精確、快速和在有限空間內(nèi)的柔性調(diào)姿要求。因此，需研發(fā)一種大型航天相機全自動六維姿態(tài)調(diào)整、轉(zhuǎn)運一體化設(shè)備，這對提高大型相機轉(zhuǎn)運及定標(biāo)測試效率、提高產(chǎn)品安全性具有重要意義。

近年來，學(xué)者們對多自由度調(diào)姿方法進行了諸多研究：SHANG等[1]將并聯(lián)機構(gòu)(Parallel kinematics machine，PKM)引入飛機機翼柔性裝配中。宋晨等[2]研發(fā)了一種具有串并聯(lián)特征的六自由度調(diào)姿機構(gòu)，用于小型噴嘴包調(diào)姿。尤晶晶等[3]進行了Stewart冗余并聯(lián)機構(gòu)正向運動學(xué)研究。敬興久等[4]研發(fā)了一種基于氣缸驅(qū)動的Stewart平臺并聯(lián)調(diào)姿機構(gòu)，用于大型構(gòu)件裝配。顧營迎等[5]設(shè)計了六足平臺，用于衛(wèi)星平臺微振動環(huán)境的地面測試。耿明超等[6]設(shè)計了6 UPRRUS 折疊式并聯(lián)機構(gòu)。文獻[7-13]對3-RPS并聯(lián)機構(gòu)的結(jié)構(gòu)、運動學(xué)、力學(xué)特性、標(biāo)定方法等進行了研究，實現(xiàn)了柔順控制，可用于踝關(guān)節(jié)康復(fù)訓(xùn)練等。文獻[14-15]設(shè)計了航空發(fā)動機多自由度調(diào)姿裝置。徐德勇等[16]設(shè)計了一種六自由度調(diào)姿平臺。文獻[17-19]對3T1R、3-RRPaR等多種并聯(lián)機構(gòu)進行了運動學(xué)分析。多數(shù)學(xué)者針對并聯(lián)機構(gòu)自身的結(jié)構(gòu)特點和力學(xué)特性進行運動學(xué)分析、控制方法等方面的研究，而對復(fù)雜工況下任意姿態(tài)到目標(biāo)姿態(tài)的量化及精確調(diào)姿方法研究卻鮮見報道。本文針對使用需求，提出基于全向移動平臺和3-RPS并聯(lián)機構(gòu)的移動并聯(lián)式轉(zhuǎn)運調(diào)姿方法，其中全向移動平臺可實現(xiàn)相機長距離轉(zhuǎn)運，并可完成相機在XOY平面內(nèi)任意方向的運動和繞Z軸的旋轉(zhuǎn)，并聯(lián)機構(gòu)可實現(xiàn)繞X、Y軸的轉(zhuǎn)動和沿Z軸的平動。

1 系統(tǒng)概述

如圖1所示，建立的坐標(biāo)系1～6分別為：檢測點坐標(biāo)系O1X1Y1Z1、目標(biāo)位置坐標(biāo)系O2X2Y2Z2、相機端面坐標(biāo)系O3X3Y3Z3、上平臺形心坐標(biāo)系O4X4Y4Z4、下平臺形心坐標(biāo)系O5X5Y5Z5、全向移動平臺端面坐標(biāo)系O6X6Y6Z6。

為了實現(xiàn)總裝測試環(huán)節(jié)中，航天相機相對目標(biāo)設(shè)備的精確瞄準(zhǔn)，需要對相機進行六維位姿調(diào)整，本文將實際應(yīng)用中基于精測的空間位姿調(diào)整問題，簡化為目標(biāo)坐標(biāo)系和相機端面坐標(biāo)系的位姿匹配問題，即通過調(diào)整全向移動平臺和并聯(lián)調(diào)姿機構(gòu)的六維位姿，使坐標(biāo)系3運動至與坐標(biāo)系2重合，實現(xiàn)相機當(dāng)前位姿向目標(biāo)位姿的自動調(diào)整。系統(tǒng)工作流程為：①通過激光測距技術(shù)，由激光跟蹤儀測得當(dāng)前位姿下，相機端面相對于檢測點的位姿、全向移動平臺相對于檢測點的位姿，實現(xiàn)相機端面、全向移動平臺端面在全局坐標(biāo)系下位姿參數(shù)的獲取。②根據(jù)當(dāng)前位姿、全向移動平臺和并聯(lián)機構(gòu)的結(jié)構(gòu)特點，建立位姿模型，利用矩陣變換理論解算各個坐標(biāo)系之間的位姿變換矩陣，并獲得相機當(dāng)前位姿到目標(biāo)位姿的位姿變換矩陣。③通過逆向運動學(xué)模型解算，將相機的位姿矩陣解耦為全向移動平臺及并聯(lián)調(diào)姿機構(gòu)的位姿變換矩陣。④利用該矩陣反算得到全向移動平臺和并聯(lián)調(diào)姿機構(gòu)各軸的運動參數(shù)，以該參數(shù)作為運動輸入，解算各軸驅(qū)動量，通過多軸組合運動實現(xiàn)相機當(dāng)前位置向目標(biāo)位置的自動調(diào)整。

2 位姿矩陣解算

2.1 當(dāng)前位姿解算

On1表示坐標(biāo)系n相對于坐標(biāo)系1(基準(zhǔn)坐標(biāo)系)的位姿矩陣，Gab表示坐標(biāo)系a相對于坐標(biāo)系b的位姿矩陣。激光跟蹤儀分別測得坐標(biāo)系3、6相對于基準(zhǔn)坐標(biāo)系的位姿矩陣O31和O61，再由坐標(biāo)系1、2之間，3、4之間， 5、6之間的相對平移關(guān)系，可知其兩兩之間的位姿矩陣為

(1)

(2)

(3)

借助機器人運動學(xué)理論[20]，按照右乘法則，解算當(dāng)前姿態(tài)下各坐標(biāo)系在基準(zhǔn)坐標(biāo)系下的位姿矩陣分別為

(4)

則通過矩陣變換，獲取的坐標(biāo)系2～5相對于基準(zhǔn)坐標(biāo)系的位姿矩陣，即為各坐標(biāo)系的當(dāng)前位姿。

2.2 當(dāng)前姿態(tài)相對于目標(biāo)的姿態(tài)矩陣解算

如圖2所示，姿態(tài)調(diào)整的目標(biāo)姿態(tài)為：坐標(biāo)系3借助調(diào)姿機構(gòu)及全向移動平臺的運動，與坐標(biāo)系2重合。

假設(shè)此過程中，O3運動至O′3，O4運動至O′4，O5運動至O′5，O6運動至O′6。其中O3、O2重合，GO′3O2=I，GO′3O′4=GO3O4，GO′5O′6=GO5O6。則有

GO′3O2=I=GO′3O′4GO′4O′5GO′5O′6OO′6O2=GO3O4GO′4O′5GO5O6(GO′6O6OO6O2)

(5)

解算得

(6)

其中

(7)

(8)

式中GO′4O′5——并聯(lián)機構(gòu)的當(dāng)前姿態(tài)與目標(biāo)姿態(tài)間的位姿變換矩陣

GO′6O6——全向移動平臺的當(dāng)前姿態(tài)與目標(biāo)姿態(tài)間的位姿變換矩陣

并聯(lián)機構(gòu)運動到目標(biāo)位置的過程為：先繞Y軸旋轉(zhuǎn)β(俯仰)，再繞X軸旋轉(zhuǎn)α(滾轉(zhuǎn))，然后升高到z，且俯仰和滾轉(zhuǎn)運動中沿X軸方向和Y軸方向的伴隨運動分別為rcos(2α)(1-cosβ)/2、-rsin(2α)(1-cosβ)/2[21]，其中r為動平臺3個支點外接圓半徑。全向移動平臺運動到目標(biāo)位置的過程為沿X和Y軸方向平移xveh和yveh、繞Z軸旋轉(zhuǎn)γ。

將式(7)、(8)代入式(6)，根據(jù)實測的姿態(tài)矩陣和各坐標(biāo)系之間的固有關(guān)系，可分別解算出全向移動平臺及并聯(lián)調(diào)姿機構(gòu)的運動參數(shù)(xveh，yveh，γ)和(α，β，z)。

3 調(diào)姿運動實現(xiàn)

3.1 全向移動平臺運動

為實現(xiàn)上述全向移動平臺由當(dāng)前姿態(tài)到目標(biāo)姿態(tài)的各軸運動參數(shù)xveh、yveh、γ的自動調(diào)整，需建立全向移動平臺的4個麥克納姆輪全向輪組的運動學(xué)模型，解算各輪組運動和平臺運動的關(guān)系，通過對輪組的協(xié)同控制，實現(xiàn)平臺在平面內(nèi)連續(xù)、高精度位置及角度調(diào)整,三維模型如圖3a所示。

對建立的運動學(xué)模型構(gòu)建如圖3b所示的坐標(biāo)系：參考坐標(biāo)系O6-0X6-0Y6-0Z6-0以全向移動平臺的位置中心為原點，坐標(biāo)系O6-iX6-iY6-iZ6-i(i=1，2，3，4)以各個麥克納姆輪位置中心為原點，各坐標(biāo)系與全向移動平臺端面坐標(biāo)系O6X6Y6Z6均為平行系。假設(shè)坐標(biāo)原點O6-0與O6-i在X軸方向距離為L，Y軸方向距離為l，車輪轉(zhuǎn)動時線速度為viw，輥子和地面接觸時的線速度為vir，車輪半徑Rveh，輥子軸線與輪轂軸線夾角為α0，全向移動平臺平移速度v=(vx，vy)，繞Z軸旋轉(zhuǎn)角速度為ωz。根據(jù)運動學(xué)模型和逆向運動學(xué)分解，各個車輪轉(zhuǎn)動角速度ωi(i=1，2，3，4)分別為

(9)

根據(jù)全向移動平臺的各軸運動量，設(shè)定運動時間后，已知平臺整體平移速度和繞中心的轉(zhuǎn)動角速度，由式(9)可解算平臺各個車輪轉(zhuǎn)速，進而基于多軸協(xié)同控制實現(xiàn)當(dāng)前姿態(tài)到目標(biāo)姿態(tài)的3自由度自動調(diào)整。

3.2 并聯(lián)調(diào)姿機構(gòu)運動

為了實現(xiàn)并聯(lián)調(diào)姿機構(gòu)由當(dāng)前姿態(tài)到目標(biāo)姿態(tài)(α，β，z)的自動調(diào)整，對機構(gòu)分析如下：并聯(lián)機構(gòu)采用3-RPS機構(gòu)，如圖4所示，其上下平臺都是正三角形，兩平臺以3個分支相聯(lián)，每個分支由3個運動關(guān)節(jié)連接，上平臺各個支點分別連接一個球關(guān)節(jié)，下平臺的每個支點各自連接一個轉(zhuǎn)動關(guān)節(jié)，球關(guān)節(jié)和轉(zhuǎn)動關(guān)節(jié)通過移動關(guān)節(jié)相連[22]。其中移動關(guān)節(jié)采用三級液壓缸，電機驅(qū)動液壓缸伸縮，通過桿長變化實現(xiàn)上平臺俯仰、滾轉(zhuǎn)和升降[23-25]。并聯(lián)調(diào)姿機構(gòu)的下平臺外接圓半徑R為525 mm，上平臺外接圓半徑r為425 mm，液壓缸收縮態(tài)長度為440 mm，伸出態(tài)長度達1 050 mm，Z向可升高616.7 mm，滿足大展收比的需求。

初始姿態(tài)下，以并聯(lián)機構(gòu)上下平臺的形心為坐標(biāo)原點建立的坐標(biāo)系分別為O4X4Y4Z4和O5X5Y5Z5，目標(biāo)姿態(tài)下，坐標(biāo)系分別記為O′4X′4Y′4Z′4和O′5X′5Y′5Z′5。并聯(lián)機構(gòu)的結(jié)構(gòu)特點確定了下平臺支點Bi(i=1,2,3)在坐標(biāo)系O′5X′5Y′5Z′5下的坐標(biāo)(bXi，bYi，bZi)，再根據(jù)上平臺支點Ai在上平臺形心坐標(biāo)系O′4X′4Y′4Z′4中的坐標(biāo)(axi，ayi，azi)，將并聯(lián)調(diào)姿機構(gòu)的目標(biāo)姿態(tài)代入式(5)可得目標(biāo)姿態(tài)上、下平臺形心坐標(biāo)系的位姿關(guān)系GO′4O′5，則下平臺各支點Bi在下平臺形心坐標(biāo)系下的坐標(biāo)(aXi，aYi，aZi)為

[aXiaYiaZi1]T=GO′4O′5[axiayiazi1]T

(10)

可得目標(biāo)姿態(tài)下，各缸長度為

(11)

通過控制3桿的長度進行所需目標(biāo)姿態(tài)的精確調(diào)整。

4 調(diào)姿算法試驗

設(shè)備實物如圖5所示。在已知姿態(tài)下對調(diào)姿算法中所需的姿態(tài)矩陣進行測試，將實測姿態(tài)數(shù)據(jù)及各矩陣間姿態(tài)關(guān)系作為輸入條件，調(diào)姿控制系統(tǒng)自動解算各軸運動量，驗證移動并聯(lián)式六自由度調(diào)姿技術(shù)的可行性。

為了便于測量，將目標(biāo)位置設(shè)置為檢測點，則GO2O1=I。且根據(jù)設(shè)備本身結(jié)構(gòu)特點，姿態(tài)關(guān)系為

試驗數(shù)據(jù)如表1所示。其中，初始姿態(tài)下相機端面坐標(biāo)系O3X3Y3Z3、全向移動平臺端面坐標(biāo)系O3X3Y3Z3均與監(jiān)測點坐標(biāo)系平行，其他當(dāng)前姿態(tài)為以初始姿態(tài)為基礎(chǔ)進行的姿態(tài)轉(zhuǎn)變，各軸運動量解算結(jié)果中，xveh與初始姿態(tài)下O31(1,4)對應(yīng)，yveh與初始姿態(tài)下O31(2,4)對應(yīng)，γ與對應(yīng)的當(dāng)前姿態(tài)下全向移動平臺旋轉(zhuǎn)角對應(yīng)，l1、l2、l3與目標(biāo)姿態(tài)下坐標(biāo)系4和坐標(biāo)系5的相對姿態(tài)有關(guān)。為了便于驗證，表中目標(biāo)姿態(tài)下坐標(biāo)系4和5平行，高度差d為500 mm，因此并聯(lián)調(diào)姿機構(gòu)3桿長度應(yīng)相同，均為

表1 調(diào)姿算法試驗數(shù)據(jù)

經(jīng)試驗驗證，基于移動并聯(lián)式六自由度調(diào)姿技術(shù)將多軸耦合的空間六自由度完全解耦，實現(xiàn)全向移動平臺和并聯(lián)調(diào)姿機構(gòu)分別調(diào)姿，姿態(tài)互補，調(diào)姿設(shè)備能夠快速精準(zhǔn)地實現(xiàn)大型相機任意姿態(tài)向目標(biāo)姿態(tài)的自動化柔性姿態(tài)調(diào)整，全向移動平臺移動精度優(yōu)于0.3 mm，旋轉(zhuǎn)精度優(yōu)于0.05°，桿長精度優(yōu)于0.5 mm，驗證了調(diào)姿算法的可行性和準(zhǔn)確性。

5 結(jié)束語

基于機器人運動學(xué)理論，對航天相機輻射定標(biāo)測試階段的自動調(diào)姿過程及快速調(diào)整算法進行了研究，提出一種基于激光測量和多姿態(tài)矩陣自動解算的移動并聯(lián)式六自由度自動解耦算法。利用測得的實際姿態(tài)和各坐標(biāo)系位姿關(guān)系，自動解算全向移動平臺和并聯(lián)調(diào)姿機構(gòu)當(dāng)前姿態(tài)到目標(biāo)姿態(tài)的位姿變換矩陣，然后對各自由度解耦反算各軸的運動參數(shù)，通過控制全向移動平臺麥克納姆輪4個輪組的轉(zhuǎn)動速度、方向和時間，控制并聯(lián)調(diào)姿機構(gòu)3桿長度，從而實現(xiàn)到目標(biāo)姿態(tài)的自動調(diào)整。經(jīng)試驗驗證，全向移動平臺移動精度優(yōu)于0.3 mm，旋轉(zhuǎn)精度優(yōu)于0.05°，桿長精度優(yōu)于0.5 mm。