余德文,趙強(qiáng)強(qiáng),陳飛飛,郭俊康,洪軍

(1.西安交通大學(xué)現(xiàn)代設(shè)計(jì)及轉(zhuǎn)子軸承系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室,710049,西安;2.上海宇航系統(tǒng)工程研究所,201109,上海)

由于運(yùn)載火箭容積限制以及航天領(lǐng)域?qū)臻g天線性能要求的提高,大口徑、高精度的空間可展機(jī)構(gòu)一直是空間技術(shù)的研究熱點(diǎn)。對于衛(wèi)星天線來說,裝配精度很大程度上決定著工作性能。在實(shí)際工程中,由于受到桿件加工誤差、鉸鏈裝配誤差的影響,空間可展機(jī)構(gòu)裝配精度往往達(dá)不到設(shè)計(jì)要求,為此工程人員經(jīng)常在地面裝配過程中對其支撐桿件進(jìn)行尺寸調(diào)整。然而,由于不清楚誤差源與裝配精度之間的關(guān)系,航天企業(yè)目前采用測量、調(diào)整、測量的純經(jīng)驗(yàn)盲調(diào)模式。因此,從提高裝配精度和實(shí)現(xiàn)定量裝調(diào)的角度出發(fā),一套行之有效的裝配精度分析方法,對于保障空間可展機(jī)構(gòu)裝配質(zhì)量和提高衛(wèi)星天線服役性能顯得尤為重要。

作為空間可展機(jī)構(gòu)的典型代表,星載SAR天線空間可展機(jī)構(gòu)備受學(xué)者關(guān)注。早期關(guān)于該類空間可展機(jī)構(gòu)的研究主要集中在構(gòu)型優(yōu)化設(shè)計(jì)[1-5]和運(yùn)動(dòng)機(jī)理分析[6-11]上。繼Campbell等發(fā)明可展支撐桁架結(jié)構(gòu)[12]后,又經(jīng)設(shè)計(jì)人員的不斷改進(jìn)[13-15],星載SAR天線空間展開機(jī)構(gòu)已成功應(yīng)用于美國SEASAT衛(wèi)星[16]、歐空局ERS系列衛(wèi)星[17]、加拿大RADARSAT系列衛(wèi)星[18]、日本ALOS衛(wèi)星[19]以及我國高分三號衛(wèi)星[20]。

隨著運(yùn)動(dòng)控制趨于成熟,該類可展機(jī)構(gòu)的拓?fù)浣Y(jié)構(gòu)也逐漸定型,具有桿件與環(huán)約束多、環(huán)約束互相耦合等特點(diǎn),給幾何精度分析帶來了極大困難,目前只有少數(shù)學(xué)者進(jìn)行了相關(guān)探索。Hedgepeth提出以天線平面節(jié)點(diǎn)偏差均方根為評價(jià)指標(biāo)進(jìn)行型面精度分析[21];Greene建立了一種基于有限元分析的評估桿件單元長度不確定性對天線性能影響的方法[22];Mobrem進(jìn)一步總結(jié)出分析空間可展機(jī)構(gòu)桿長對天線性能影響的蒙特卡羅法、逆頻率平方法、直接法與正則模法[23];受此啟發(fā),吳建云等利用蒙特卡洛法,研究了可展機(jī)構(gòu)中鉸鏈鎖定位置偏差對展開精度的影響[24]。

上述精度分析方法是建立在關(guān)鍵設(shè)計(jì)參數(shù)隨機(jī)分布的假設(shè)之上,結(jié)果無法指導(dǎo)具體產(chǎn)品的裝調(diào)[21-24]。另外,由于星載SAR天線空間可展機(jī)構(gòu)在理想情況下,構(gòu)型關(guān)于展開方向與天線板垂線方向所構(gòu)成的平面對稱,因此一些研究者在分析其幾何精度時(shí)按平面連桿機(jī)構(gòu)簡化處理。例如,文獻(xiàn)[25-26]構(gòu)建了單桿件固定、雙桿件連接和多余桿件連接3種裝配單元的誤差模型,進(jìn)而提出了多環(huán)閉鏈可展機(jī)構(gòu)的精度預(yù)測算法,此模型對三維構(gòu)型裝配誤差分析具有借鑒意義。文獻(xiàn)[27]從幾何誤差角度求解了極惡劣工況下的指向誤差。此外,張武翔等以結(jié)構(gòu)勢能最小為目標(biāo)函數(shù),建立了空間可展機(jī)構(gòu)尺寸調(diào)整的數(shù)學(xué)模型[28],然而定義的綜合調(diào)整量沒有真實(shí)的物理含義,并且去掉了可展機(jī)構(gòu)中的部分支撐桿件,減少了實(shí)際存在的環(huán)約束,因此工程應(yīng)用具有較大局限性。

本文提出了一種空間可展機(jī)構(gòu)裝配誤差建模和靈敏度分析方法。首先,通過空間可展機(jī)構(gòu)的等效轉(zhuǎn)換將其分為兩部分,基于閉環(huán)矢量與虛位移法,分別建立這兩部分結(jié)構(gòu)的裝配誤差模型,經(jīng)線性疊加后,全面表征出天線板位姿誤差與誤差源的函數(shù)關(guān)系,實(shí)現(xiàn)整個(gè)空間可展機(jī)構(gòu)的幾何精度預(yù)測。在此基礎(chǔ)上,運(yùn)用偏微分法分析裝配精度受誤差源的影響情況,從而識別出關(guān)鍵誤差源,并通過具體數(shù)值算例驗(yàn)證所提方法的有效性。

1 問題描述

圖1給出了星載SAR天線空間可展機(jī)構(gòu)的展開構(gòu)型,主要由天線內(nèi)板、天線外板、支撐桿件、90°鎖定鉸鏈和180°鎖定鉸鏈組成。天線內(nèi)板通過兩個(gè)90°鎖定鉸鏈與星體連接,天線外板通過兩個(gè)180°鎖定鉸鏈與天線內(nèi)板連接,天線內(nèi)外板依靠支撐桿件增強(qiáng)與星體的連接剛度。可以看出,空間可展機(jī)構(gòu)實(shí)質(zhì)上是一個(gè)超靜定過約束機(jī)構(gòu)。

圖1 星載SAR天線空間可展機(jī)構(gòu)

在實(shí)際裝配過程中,空間可展機(jī)構(gòu)調(diào)不準(zhǔn)的頑疾至今仍未根除。因此,本文所需解決的問題是:如何從幾何角度出發(fā),構(gòu)建空間可展機(jī)構(gòu)裝配精度與誤差源的函數(shù)關(guān)系,并找出影響裝配精度的關(guān)鍵桿件,以實(shí)現(xiàn)空間可展機(jī)構(gòu)的裝配誤差預(yù)測和精準(zhǔn)定量裝調(diào)。

2 裝配誤差建模

2.1 結(jié)構(gòu)等效轉(zhuǎn)換

如前文所述,空間可展機(jī)構(gòu)的展開構(gòu)型屬于典型超靜定結(jié)構(gòu)。在不考慮變形且不影響裝配精度的前提下,必須進(jìn)行結(jié)構(gòu)等效轉(zhuǎn)換以消除過約束。為此,將空間可展機(jī)構(gòu)分為圖2所示的兩部分。第一部分由外撐桿、中撐桿、內(nèi)撐桿及天線面板組成的不含星體連接桿的支撐結(jié)構(gòu);第二部分由星體連接桿與虛擬桿件(連接撐桿交匯中心O與90°鎖定鉸鏈幾何對稱中心B的假想桿件,其簡要表征不含星體連接桿的支撐結(jié)構(gòu))組成的雙桿件連接結(jié)構(gòu)。

上述轉(zhuǎn)換的理由如下:由于SAR天線內(nèi)外板間180°鎖定鉸鏈的間隙十分微小,并且鉸鏈鎖定后的結(jié)構(gòu)剛度較大,因此兩板原有的平面度在整個(gè)裝調(diào)過程中保持不變,從而可把內(nèi)外板視為一塊整板。不同于180°鎖定鉸鏈的是,連接星體與天線內(nèi)板的90°鎖定鉸鏈間隙相對較大,容許天線面板一定范圍內(nèi)的裝配誤差。綜上所述,將空間可展機(jī)構(gòu)分成兩部分結(jié)構(gòu)具有合理性。

圖2 空間可展機(jī)構(gòu)等效轉(zhuǎn)換模型

2.2 不含星體連接桿的支撐結(jié)構(gòu)誤差建模

將支撐結(jié)構(gòu)第一部分視作6-SPS(S-球鉸副,P-移動(dòng)副)結(jié)構(gòu)。一方面,桿長調(diào)整(加減墊片)的過程類似于移動(dòng)副的直線運(yùn)動(dòng),實(shí)質(zhì)均是沿特定方向的特征尺寸發(fā)生改變;另一方面,由于天線面板位姿變化空間十分微小,且鉸鏈轉(zhuǎn)動(dòng)副間存在徑向間隙,因此球鉸副代替轉(zhuǎn)動(dòng)副為工程所接受。

如圖3所示,在6桿交匯中心的位置O建立局部坐標(biāo)系O1-x1y1z1,x1軸正向沿天線面板長度方向并遠(yuǎn)離星體,y1軸正向沿天線面板寬度方向且豎直向下,z1軸正向由右手法則確定。此外,在天線面板的幾何中心A建立局部坐標(biāo)系O2-x2y2z2,該坐標(biāo)系與坐標(biāo)系O1-x1y1z1平行。

圖3 不含星體連接桿的支撐結(jié)構(gòu)

在第i個(gè)封閉環(huán)O-A-Ai-O(i=1,2,3,4,5,6)中,滿足閉環(huán)矢量方程[29]

liui=r+Rai

(1)

式中:ui為由O指向Ai的單位矢量;li為桿件OAi名義長度;r為坐標(biāo)系O1-x1y1z1下由O指向A的矢量;R為坐標(biāo)系O2-x2y2z2到坐標(biāo)系O1-x1y1z1的齊次變換矩陣;ai是Ai在O2-x2y2z2坐標(biāo)系下的坐標(biāo)。對式(1)兩邊全微分,有

dliui+lidui=dr+dRai+Rdai

(2)

(3)

(4)

(5)

在式(3)中,與名義桿長相比,桿長誤差均可視作無窮小量,即dli=Δli,dr=Δr,dai=Δai。接著,把式(4)帶入式(3)中,整理后可得

(6)

式(6)的矩陣表達(dá)形式如下

(7)

6個(gè)封閉環(huán)O-A-Ai-O均滿足等式(7),將其完全組合有

ΔL=JvΔV+JbΔB

(8)

式中:ΔL=[Δl1,Δl2,Δl3,Δl4,Δl5,Δl6]T,分別對應(yīng)6桿的長度誤差;ΔB=[Δa1,Δa2,Δa3,Δa4,Δa5,Δa6]T,分別為6個(gè)鉸鏈的安裝位置誤差,有Δai=[Δxi,Δyi,Δzi];ΔV=[Δr,Δθ]T,Δr、Δθ分別為天線板相對于坐標(biāo)系O1-x1y1z1的位置誤差、姿態(tài)誤差;Jv、Jb為ΔV、ΔB相應(yīng)的系數(shù)矩陣,其具體表達(dá)式如下

(9)

由式(8)可得到不含星體連接桿的支撐結(jié)構(gòu)裝配誤差的正解方程

(10)

至此,式(10)表達(dá)了天線面板位姿誤差與6個(gè)桿長誤差、6個(gè)鉸鏈安裝位置誤差的函數(shù)關(guān)系,而該位姿誤差是在坐標(biāo)系O1-x1y1z1下衡量的。只要進(jìn)一步得到坐標(biāo)系O1-x1y1z1與全局坐標(biāo)系的相對位姿關(guān)系,便能在全局坐標(biāo)系下求解出天線面板的裝配誤差。

2.3 含星體連接桿的支撐結(jié)構(gòu)誤差建模

如圖4所示,全局坐標(biāo)系Og-xyz的原點(diǎn)為兩個(gè)90°鎖定鉸鏈理論安裝位置的幾何對稱中心B,x軸正向?yàn)檠靥炀€面板長度方向并遠(yuǎn)離星體,y軸正向?yàn)檠靥炀€面板寬度方向且豎直向下,z軸正向根據(jù)右手法則確定。在全局坐標(biāo)系Og-xyz中,局部坐標(biāo)系O1-x1y1z1坐標(biāo)原點(diǎn)O(O1)受星體連接桿影響,名義位置與實(shí)際位置之間存在裝配誤差。因此,下面將構(gòu)建雙桿件連接裝配誤差模型,并確定O點(diǎn)位置誤差對空間可展機(jī)構(gòu)幾何裝配精度的影響。

圖4 雙桿件連接結(jié)構(gòu)及其平面表征

由于全局坐標(biāo)系Og-xyz和局部坐標(biāo)系O1-x1y1z1完全平行,星體連接桿CO和虛擬桿BO都在xz平面,因而投影至平面表征出如圖4所示的雙桿件連接結(jié)構(gòu)。在圖4情況下,O點(diǎn)的裝配位置誤差由桿件BO、CO的桿長誤差與鉸鏈連接處B、C的位置誤差造成。因此,O點(diǎn)位置誤差可以表示為

ΔO=μ1Δl+μ2ΔB+μ3ΔC

(11)

式中:ΔO表示O點(diǎn)位置誤差;Δl=[Δl02,Δl0]T,Δl02、Δl0分別為BO、CO的桿長誤差;ΔB=[ΔzB,ΔxB]T,ΔzB、ΔxB分別表示B在全局坐標(biāo)系下沿z軸、x軸的位置誤差;與ΔB類似,ΔC=[ΔzC,ΔxC]T表示鉸鏈C處的位置誤差;μ1、μ2、μ3分別為相應(yīng)誤差系數(shù)矩陣,且定義如下

(12)

圖5 系數(shù)矩陣部分表達(dá)式計(jì)算原理圖

(13)

通過構(gòu)造與圖5相似的結(jié)構(gòu)可以求解系數(shù)矩陣的其它比值,由此得到的3個(gè)系數(shù)矩陣具體為

(14a)

μ2=

(14b)

μ3=

(14c)

式中:α、β、γ為名義尺寸下各鉸接點(diǎn)所組成三角形的內(nèi)角;ε為BC連線與z軸負(fù)方向的夾角。

需要特別指出的是,利用式(11)求解O點(diǎn)的裝配位置誤差時(shí),兩個(gè)90°鎖定鉸鏈理論安裝位置的幾何對稱中心B點(diǎn)作為基準(zhǔn)點(diǎn),其裝配誤差為0。另外,虛擬桿件BO在理想條件下的矢量表達(dá)式為

l02=r+l01

(15)

(16)

Δl02=((rx+Δx-l01)2+(ry+Δy-δzl01)2+

(17)

將式(17)(14)帶入式(11),可計(jì)算出O點(diǎn)的裝配誤差。注意到ΔO只在xz平面,不影響繞x方向和z方向的轉(zhuǎn)動(dòng)偏差,而所引起空間可展機(jī)構(gòu)繞y軸的轉(zhuǎn)動(dòng)偏差為

ξy=‖ΔO‖/l02

(18)

利用線性疊加,綜合式(10)(18),得到考慮所有桿長誤差與鉸鏈安裝位置誤差的空間可展機(jī)構(gòu)幾何裝配精度模型

(19)

式中:ΔP=[Tx,Ty,Tz,Rx,Ry,Rz]T,為全局坐標(biāo)系下空間可展機(jī)構(gòu)的位姿誤差。

3 靈敏度分析

由于在實(shí)際裝調(diào)過程中,可展機(jī)構(gòu)繞y軸的偏擺誤差通常更難以滿足設(shè)計(jì)要求。因此,下面著重分析各誤差因素對可展機(jī)構(gòu)y軸姿態(tài)誤差的影響,而其余的靈敏度分析可據(jù)此類推。

3.1 桿長誤差影響系數(shù)

基于空間可展機(jī)構(gòu)幾何裝配精度模型,由式(19)的位姿誤差與誤差源的函數(shù)關(guān)系,可分離出到可展機(jī)構(gòu)y軸偏擺誤差與各項(xiàng)誤差的關(guān)系,即

Ry=δy+ξy

(20)

(21)

(22)

式中:J1i,…,J6i為誤差傳遞矩陣J各行第i列的值。

另外,從式(10)可知,ΔV=[Δx,Δy,Δz,δx,δy,δz]T均是桿長誤差Δli(i=1,2,3,4,5,6)的函數(shù)。從式(17)可以看出Δl02是ΔV的函數(shù)。以此推出,Δl02也是Δli(i=1,2,3,4,5,6)的函數(shù)。于是,式(17)兩邊對桿長誤差Δli求偏導(dǎo)有

(23)

(24a)

i=1,2,3,4,5,6

(24b)

令

(25)

綜合式(23)(24)(25),各項(xiàng)桿長誤差對天線面板y軸偏擺誤差影響的靈敏度表達(dá)式如下

(26)

3.2 鉸鏈安裝位置誤差影響系數(shù)

類似地,按照推導(dǎo)桿長誤差影響系數(shù)的思路,利用式(11)和式(20),可以得出C處鉸鏈安裝位置誤差對天線面板繞y軸偏擺誤差影響的靈敏度表達(dá)式如下

(27a)

(27b)

除C處鉸鏈外,天線面板上Aj(j=1,2,3,4,5,6)處鉸鏈同樣存在安裝位置誤差。為簡化表達(dá),以Δajk(j=1,2,3,4,5,6;k=x,y,z)表示Aj處鉸鏈在k方向的安裝位置誤差。由式(10)(17)可知,Δl02也是Δajk的函數(shù),進(jìn)而對表達(dá)式(17)進(jìn)行相應(yīng)鉸鏈安裝位置誤差取偏微分,可以得到

(28)

在此基礎(chǔ)上,構(gòu)建映射關(guān)系:當(dāng)k=x時(shí),n=1;當(dāng)k=y時(shí),n=2;當(dāng)k=z時(shí),n=3。令m=3(j+1)+n(j=1,2,3,4,5,6;n=1,2,3),利用式(20),推導(dǎo)出天線面板處鉸鏈安裝位置誤差對天線面板繞y軸偏擺誤差影響的靈敏度表達(dá)式如下

(q3J5m-q4J6m+q5J1m+q6J2m+q7J3m)

(29)

式中:J1m,…,J6m為誤差傳遞矩陣J各行第m列的值。

至此,式(27)和式(29)完全表達(dá)了各鉸鏈安裝位置誤差對天線面板y軸偏擺誤差影響的靈敏度。

4 案例分析

某星載SAR天線空間可展機(jī)構(gòu)的具體尺寸如表1所示,天線陣面長為2.22 m,寬為1.65 m。為保障衛(wèi)星正常服役所需的電性能,在地面裝配中,特別要求天線面板俯仰角與方位角均優(yōu)于0.020°,該精度指標(biāo)對應(yīng)本模型中的y軸偏角(Ry)與x軸偏角(Rx)。

該機(jī)構(gòu)屬于精密桿系機(jī)構(gòu),桿件加工精度等級取11級,中撐桿的半公差帶寬為0.66 mm,其余桿件的半公差帶寬均為1.10 mm。實(shí)際裝配過程中,鉸鏈安裝位置誤差在3個(gè)坐標(biāo)方向上均不超過0.20 mm。在上述相應(yīng)誤差范圍內(nèi),給定5組如表2所示的桿長誤差和鉸鏈安裝位置誤差。運(yùn)用MATLAB編寫幾何裝配精度分析程序,根據(jù)式(19)計(jì)算出裝配誤差,結(jié)果詳見表3。

表2 桿長與鉸鏈安裝位置誤差值

表3 空間可展機(jī)構(gòu)位姿誤差

從表2與表3可以看出,在誤差組1、2、3作用下,天線面板繞y軸方向的偏擺誤差(Ry)均大于0.02°。這說明,即使桿長誤差滿足設(shè)計(jì)要求,鉸鏈安裝符合工藝要求,依舊無法完全保證天線面板的指向精度。為進(jìn)一步反映桿件動(dòng)態(tài)調(diào)整過程,在誤差組1的基礎(chǔ)上,只調(diào)整支撐桿件長度而保持鉸鏈安裝位置誤差不變,依據(jù)誤差組1、3、5的計(jì)算結(jié)果,注意到Rx始終在允許范圍內(nèi),而Ry從0.072 3°到0.021 4°再到0.007 5°,因此逐漸滿足裝配精度要求。較誤差組3而言,誤差組4全面調(diào)整了安裝位置誤差,誤差組5僅微調(diào)內(nèi)撐桿Δl5與星體桿Δl0的尺寸,兩種方法均獲得了理想的裝配精度。由此可見,在滿足桿件加工精度和鉸鏈安裝精度條件下,盡管調(diào)整鉸鏈安裝位置也能限制天線面板的裝配誤差,但考慮到裝調(diào)時(shí)間和裝調(diào)難度,調(diào)整桿長不失為更合理的選擇。另外,從表3還可以看出,上述誤差作用下的天線面板沿y軸方向的位移偏差(Ty)幾乎為0,且繞z軸方向偏擺誤差(Rz)的絕對值不超過0.001°,這符合平面天線空間展開機(jī)構(gòu)完全展開狀態(tài)下的結(jié)構(gòu)特征。

正如上文分析結(jié)果,在實(shí)際裝調(diào)過程中,工廠也只選擇對桿長加減墊片進(jìn)行空間可展機(jī)構(gòu)裝配精度的調(diào)整。為進(jìn)一步識別出關(guān)鍵桿件以降低裝配難度,依據(jù)式(26),計(jì)算出各桿長誤差對天線面板繞y軸偏擺誤差影響的靈敏度值,詳見表4。

表4 各桿長誤差對y軸偏擺誤差影響的靈敏度

為了更好地表達(dá)誤差參數(shù)靈敏度,進(jìn)行歸一化處理,即

(30)

式中:Si為各誤差源誤差對天線面板y軸偏擺誤差影響的靈敏度,在該數(shù)值案例中,即特指桿長誤差對天線面板繞y軸偏擺角度誤差影響的靈敏度。根據(jù)表4計(jì)算值,對誤差參數(shù)靈敏度進(jìn)行歸一化處理,結(jié)果詳見圖6。

圖6 各參數(shù)靈敏度歸一化處理結(jié)果

從圖5可知,結(jié)構(gòu)對稱的兩撐桿對天線面板繞y軸偏擺誤差影響程度相同。此外,星體連接桿的長度誤差對天線面板繞y軸偏擺誤差影響最大,兩內(nèi)支撐桿長度誤差對天線面板繞y軸偏擺誤差影響最小,前者的歸一化敏感度是后者的35.47倍,這意味著調(diào)整星體連接桿長度所致使天線裝配精度的變動(dòng)效果更為顯著。

5 結(jié) 論

(1)建立了空間可展機(jī)構(gòu)的幾何裝配精度模型。從空間構(gòu)型出發(fā),提出了將可展機(jī)構(gòu)等效為6-SPS與雙桿件連接結(jié)構(gòu)的轉(zhuǎn)化方法,并推導(dǎo)了兩部分結(jié)構(gòu)的裝配誤差表達(dá)式。在此基礎(chǔ)上,結(jié)合幾何約束,構(gòu)建了所有桿長誤差與鉸鏈安裝位置誤差等誤差源與天線面板位姿精度的顯式函數(shù)關(guān)系。

(2)基于可展機(jī)構(gòu)幾何裝配精度模型,利用偏微分法分析了桿長誤差和鉸鏈安裝位置誤差對可展機(jī)構(gòu)y軸偏擺誤差影響的靈敏度,從而能夠有效識別出關(guān)鍵誤差,有助于降低裝配難度。

(3)通過星載SAR天線空間可展機(jī)構(gòu)的分析實(shí)例,所預(yù)測的裝配精度能滿足工程要求,確定星體連接桿為關(guān)鍵桿,驗(yàn)證了空間可展機(jī)構(gòu)裝配誤差建模和靈敏度分析方法的有效性。本文的研究成果對空間可展機(jī)構(gòu)的裝調(diào)具有指導(dǎo)意義和參考價(jià)值。