彭瑞 謝任遠(yuǎn)

【摘 要】為驗(yàn)證某衛(wèi)星轉(zhuǎn)動載荷的大角動量補(bǔ)償功能是否滿足實(shí)際使用需求，設(shè)計(jì)出基于單軸氣浮轉(zhuǎn)臺，利用單軸氣浮臺測角裝置和陀螺組合，對轉(zhuǎn)動載荷在不同質(zhì)量（轉(zhuǎn)動慣量相同）、不同轉(zhuǎn)動角速度的情況下的運(yùn)動特性進(jìn)行測量的地面試驗(yàn)系統(tǒng)，并以某一工況為例，得出載荷轉(zhuǎn)動的時(shí)域與頻域特性。試驗(yàn)結(jié)果為大角動量補(bǔ)償功能的設(shè)計(jì)完善提供了依據(jù)，對衛(wèi)星后續(xù)研制意義重大。

【關(guān)鍵詞】單軸氣浮臺 頻譜分析 運(yùn)動特性 地面試驗(yàn)系統(tǒng)

1引言

隨著衛(wèi)星載荷的不斷發(fā)展和衛(wèi)星平臺技術(shù)的日漸成熟，載荷對衛(wèi)星平臺的要求越來越高，兩者的關(guān)系逐漸從載荷適用平臺發(fā)展為載荷決定平臺，由此給衛(wèi)星平臺的研制工作提出了更高的要求。為了降低研制風(fēng)險(xiǎn)， 需要通過各種類型的大量地面仿真試驗(yàn)對系統(tǒng)方案進(jìn)行驗(yàn)證[1]。

某衛(wèi)星載荷在衛(wèi)星入軌后一直處于勻速轉(zhuǎn)動狀態(tài)，其轉(zhuǎn)動產(chǎn)生的角動量由衛(wèi)星姿軌控分系統(tǒng)進(jìn)行長期補(bǔ)償，因此探測頭部旋轉(zhuǎn)角速度的特性，直接影響衛(wèi)星姿態(tài)穩(wěn)定度指標(biāo)。

本文根據(jù)某衛(wèi)星載荷的運(yùn)動特性，構(gòu)建了一個(gè)基于單軸氣浮臺的對載荷在不同質(zhì)量（轉(zhuǎn)動慣量相同）、不同轉(zhuǎn)動角速度的情況下的運(yùn)動特性進(jìn)行測量的地面試驗(yàn)系統(tǒng)，同時(shí)對測試出的轉(zhuǎn)速進(jìn)行頻譜分析，得到載荷轉(zhuǎn)動的時(shí)域與頻域特性，為優(yōu)化角動量補(bǔ)償控制設(shè)計(jì)提供了依據(jù)。

2 頻譜分析理論及在Matlab中的應(yīng)用

控制系統(tǒng)的設(shè)計(jì)應(yīng)該從分析性能要求開始，而頻譜分析是性能分析的基礎(chǔ)。另外，頻譜分析也是設(shè)計(jì)和分析中常用到的一種數(shù)據(jù)處理手段。因此，設(shè)計(jì)者不僅應(yīng)該知道頻譜的概念，還應(yīng)該掌握頻譜的分析方法[2]。

2.1傅里葉級數(shù)

對于周期為T的周期函數(shù)f（t）：

（1）

如果f（t）滿足狄里赫利條件：在區(qū)間T上有界，且僅有有限個(gè)極大值和極小值，則f（t）可用收斂的傅立葉級數(shù)來表示。

三角函數(shù)形式的傅里葉級數(shù)為

（2）

上式的系數(shù)為

（3）

控制系統(tǒng)中使用復(fù)數(shù)形式的傅里葉級數(shù)

（4）

其中的系數(shù) 為

（5）

為復(fù)數(shù)，一般可表示為如下形式

（6）

從（5）式可知 和 為共軛復(fù)數(shù)，對于其中的每一組可以寫成

（7）

式（7）表明，當(dāng)用復(fù)數(shù)形式表示是，復(fù)系數(shù) 的幅值 表示了第k次諧波的幅值（幅值為2 ），而 的相角 則為該次諧波的相移。這種用復(fù)數(shù)形式來表示的諧波常稱為復(fù)數(shù)正弦。

用傅里葉級數(shù)來表示函數(shù)f（t），無論是實(shí)數(shù)形式還是復(fù)數(shù)形式，都是將f（t）看成是由各次諧波所組成。傅里葉級數(shù)的系數(shù)表示了各次諧波的幅值和相位，這些系數(shù)的集合稱為頻譜。

2.2傅里葉積分和傅里葉變換

控制系統(tǒng)分析中，經(jīng)常遇到非周期函數(shù)，這時(shí)傅里葉級數(shù)就不能應(yīng)用了，但可以用傅里葉積分來處理[2]。

（8）

（9）

式（8）稱為傅里葉積分，式（9）中的 稱為函數(shù)f（t）的傅里葉變換。

設(shè)將 作為頻率的橫坐標(biāo)，這時(shí)式（8）改寫為：

（10）

對于（8）式，按照積分的概念可得：

（11）

或者可以寫為：

（12）

傅里葉積分將一個(gè)非周期函數(shù)f（t）分解為各次諧波，每一個(gè)諧波的幅值由（8）式可知為：

（13）

這幅值為無窮小，所以一般用相對幅值 來表示其頻譜。這就是說，傅立葉變換 表示的是該非周期信號諧波的分布特性，故 稱為信號的頻譜特性，簡稱為頻譜。

傅里葉變換使得我們可以從頻譜的角度來分析信號。

2.3頻譜分析在MATLAB中的應(yīng)用

離散傅里葉變換是進(jìn)行頻譜分析的基本工具，在實(shí)際中得到廣泛應(yīng)用的是快速傅里葉變換（FFT算法）?，F(xiàn)在FFT已可以用MATLAB軟件中的fft函數(shù)來進(jìn)行計(jì)算，這個(gè)函數(shù)做的是下面的變換

（14）

式（14）中 ，N是x的點(diǎn)數(shù)（即數(shù)據(jù)長度）。

實(shí)際使用中需要注意FFT算法與其它函數(shù)一起運(yùn)用。本文中頻譜分析的函數(shù)算法如下：

function[AW，f]=cftbyfft（wt，t，flag）

if nargin==2； flag=1； end

N=length（t）；

T=t（length（t）-t（1））；

dt=T/N；

W0=fft（wt）；

W=dt*W0；

df=1/T；

n=0：1：（N-1）；

if flag==0

n=-N/2：（N/2-1）；

W=fftshift（W）；

end

f=n*df；

AW=abs（W）；

if nargout==0

plot（f，AW）；grid，xlabel（'頻率f（Hz）'）；ylabel（'幅值|W（f）|'）

end

3試驗(yàn)系統(tǒng)設(shè)計(jì)

氣浮臺依靠氣源提供的壓縮空氣在氣浮軸承與軸承座之間形成氣膜， 從而使氣浮臺的轉(zhuǎn)臺浮起， 實(shí)現(xiàn)近似無摩擦的相對運(yùn)動條件，以模擬衛(wèi)星載荷在外層空間所受干擾力矩很小的力學(xué)環(huán)境[2]。某衛(wèi)星載荷在衛(wèi)星入軌后一直處于勻速轉(zhuǎn)動狀態(tài)，為測試其運(yùn)動特性，設(shè)計(jì)的試驗(yàn)系統(tǒng)組成如圖1所示。本試驗(yàn)中氣浮臺為單軸精密氣浮臺；反作用飛輪用于標(biāo)定圖1所示的氣浮臺總的轉(zhuǎn)動慣量；陀螺組合用于測量載荷轉(zhuǎn)動時(shí)氣浮臺的相對轉(zhuǎn)動角速度；載荷探測頭部模擬慣量盤用于模擬載荷探測頭部的轉(zhuǎn)動（質(zhì)量可調(diào)、轉(zhuǎn)動慣量相同），由載荷伺服控制器控制其轉(zhuǎn)動角速度。

圖1試驗(yàn)系統(tǒng)組成

圖1中，載荷探測頭部模擬慣量盤、動量輪通過工裝設(shè)備固連于單軸氣浮臺上，陀螺組合、反作用飛輪地面試驗(yàn)設(shè)備、陀螺組合地面試驗(yàn)設(shè)備直接放置在氣浮臺上；載荷伺服控制器、載荷地檢設(shè)備放置于臺下。載荷伺服控制器通過載荷地檢設(shè)備進(jìn)行供電、控制與相關(guān)通訊；反作用飛輪、陀螺組合分別通過專用測試設(shè)備進(jìn)行供電、控制與相關(guān)通訊；單軸氣浮臺通過配平塊進(jìn)行配平，通過自準(zhǔn)直儀進(jìn)行測角與測角速度；遠(yuǎn)控機(jī)通過網(wǎng)絡(luò)控制飛輪和陀螺的地面軟件，用以發(fā)送指令和存儲數(shù)據(jù)。載荷的臺上臺下連接電纜加工材料采用柔性較好軟線或散線，通過吊裝設(shè)備從氣浮臺中心上方垂下，以減小電纜產(chǎn)生的扭矩對試驗(yàn)的影響。

4試驗(yàn)數(shù)據(jù)分析和處理方法

4.1氣浮臺整體轉(zhuǎn)動慣量J0的計(jì)算

已知反作用飛輪的轉(zhuǎn)動慣量是J1=0.0465kgm?，轉(zhuǎn)速為ω1為飛輪轉(zhuǎn)速；臺體轉(zhuǎn)動慣量設(shè)為J0，轉(zhuǎn)速ω0為氣浮臺自準(zhǔn)儀的測量值；根據(jù)動量守恒公式J0*ω0=J1*ω1，可計(jì)算出氣浮臺體的轉(zhuǎn)動慣量J0=（J1*ω1）/ω0。

試驗(yàn)中共進(jìn)行了4組測試，分別根據(jù)反作用飛輪0rpm->1000rpm、1000rpm->0rpm、0rpm->-1000rpm和-1000rpm->0rpm四種工況下的氣浮臺轉(zhuǎn)速數(shù)據(jù)（△ω01、△ω02，△ω03，△ω04），算出臺體的轉(zhuǎn)動慣量（J01、J02，J03，J04），取平均值作為氣浮臺體的轉(zhuǎn)動慣量。

4.2載荷轉(zhuǎn)速波動量的計(jì)算

載荷探測頭部轉(zhuǎn)動慣量J2取理論值19 kgm?，轉(zhuǎn)速設(shè)為ω2；臺體轉(zhuǎn)動慣量為J0（4.1中計(jì)算值），轉(zhuǎn)速ω0為氣浮臺自準(zhǔn)儀測量值；根據(jù)動量守恒公式J0*ω0=J2*ω2，可計(jì)算出ω2=（J0*ω0）/J2。

試驗(yàn)中分別在模擬慣量盤質(zhì)量60kg和90kg條件下，進(jìn)行了轉(zhuǎn)動角速度108°/s、120°/s、150°/s共六種工況的測試；根據(jù)氣浮臺轉(zhuǎn)速數(shù)據(jù)（陀螺數(shù)據(jù)和自準(zhǔn)儀數(shù)據(jù)），計(jì)算并畫出每一種工況下三個(gè)轉(zhuǎn)動角速度波動量（載荷自身采集，陀螺數(shù)據(jù)換算和自準(zhǔn)儀數(shù)據(jù)換算）的數(shù)值對比曲線，得出每種工況下的載荷探測頭部轉(zhuǎn)動角速度波動量范圍。

為更好的反映運(yùn)動特性，分別對不同工況下的陀螺轉(zhuǎn)速、載荷轉(zhuǎn)速、氣浮臺轉(zhuǎn)速利用FFT進(jìn)行頻譜分析，從頻率角度進(jìn)一步分析其運(yùn)動特性。

5試驗(yàn)結(jié)果分析

試驗(yàn)共進(jìn)行了六種工況的測試，現(xiàn)取其中一種工況（質(zhì)量60kg、轉(zhuǎn)速120°/s）進(jìn)行試驗(yàn)結(jié)果分析。在此種工況下，經(jīng)測試標(biāo)定，氣浮臺整體平均慣量為2951.045 kgm2，載荷探測頭部與氣浮臺的慣量比為142。

此種工況下的轉(zhuǎn)速波動量測試結(jié)果對比曲線如圖2所示。

圖2轉(zhuǎn)速波動量測試結(jié)果對比圖

陀螺轉(zhuǎn)速、氣浮臺轉(zhuǎn)速、載荷轉(zhuǎn)速的時(shí)域與頻域分析曲線對比如圖3所示。

圖3轉(zhuǎn)速頻域分析結(jié)果對比

通過測試曲線對比分析，可以得出試驗(yàn)結(jié)果為：

（1）氣浮臺采集數(shù)據(jù)、陀螺采集數(shù)據(jù)與載荷采集數(shù)據(jù)基本一致，載荷探測頭部穩(wěn)態(tài)運(yùn)行時(shí)的轉(zhuǎn)速波動量基本在±0.2°/s范圍內(nèi)；

（2）氣浮臺、陀螺與載荷數(shù)據(jù)頻譜分析一致。基頻為0.33 Hz，有多倍頻，且基頻分量均遠(yuǎn)大于倍頻分量。

6結(jié)語

本文設(shè)計(jì)出了基于單軸氣浮臺的某衛(wèi)星載荷運(yùn)動特性測試的地面試驗(yàn)系統(tǒng)，并提出了試驗(yàn)數(shù)據(jù)分析和處理方法。

通過對某一工況的測試結(jié)果分析，獲得了載荷時(shí)域和頻域的轉(zhuǎn)動特性，根據(jù)測試和分析結(jié)果進(jìn)一步完善角動量補(bǔ)償控制的設(shè)計(jì)，對載荷的優(yōu)化設(shè)計(jì)提出進(jìn)一步要求。本次試驗(yàn)對后續(xù)研制工作意義重大。

參考文獻(xiàn)：

[1]李季蘇，牟小剛，張錦江.衛(wèi)星控制系統(tǒng)全物理仿真[J].控制工程，2003（1）：20- 26.

[2]王廣雄，何朕著.控制系統(tǒng)設(shè)計(jì).清華大學(xué)出版社，2008.3.

[3]劉飛，董云峰.基于單軸氣浮臺的空間跟瞄地面試驗(yàn)的設(shè)計(jì)與實(shí)現(xiàn).計(jì)算機(jī)測量與控制，2010，18（3）：626-628.