劉涌 胡永力 王偉之 張博文

（北京空間機電研究所，北京 100076）

1 引言

航天器在軌運行時，以動量輪為代表的旋轉(zhuǎn)活動部件，由于自身的靜動不平衡，當(dāng)其工作時會輸出擾振力和擾振力矩，使衛(wèi)星發(fā)生輕微的振動。這種振動會引起空間相機內(nèi)部結(jié)構(gòu)零件和光學(xué)系統(tǒng)元件的微位移，進(jìn)而影響成像質(zhì)量，Eyerman 等認(rèn)為飛輪（動量輪等）工作時產(chǎn)生的擾振，是影響有效載荷成像質(zhì)量的主要因素［1］。為了評估動量輪在軌擾振對相機成像質(zhì)量的影響，確定是否須要采取隔振措施，傳統(tǒng)的方法一般是首先分析擾振源對空間相機底座的響應(yīng)，相機研制單位根據(jù)底座響應(yīng)分析出擾振對成像質(zhì)量的影響。這種分析方法將各個學(xué)科、各個子系統(tǒng)分開來考慮，沒有考慮到衛(wèi)星平臺和相機的耦合，從相機的角度來說，相當(dāng)于將衛(wèi)星平臺剛體化了，這樣使得分析出來的結(jié)果相對偏保守。

為了得到更為準(zhǔn)確的分析結(jié)果，國外常常采用集成分析手段，發(fā)展出了諸如集成模型環(huán)境（Integrated Modeling Environment，IME）的專用抖動分析軟件包［2］。集成分析是指綜合地考慮結(jié)構(gòu)、控制、光學(xué)等對成像質(zhì)量有重要影響的子系統(tǒng)，運用各學(xué)科建模工具分別建立分析模型，然后通過分析目標(biāo)在各子系統(tǒng)間的相互關(guān)系進(jìn)行集成，形成一條從擾振源到光學(xué)系統(tǒng)的完整分析鏈的方法［3］。

本文根據(jù)集成分析的思想，建立了干擾源模型、控制模型、結(jié)構(gòu)模型和光學(xué)模型，并給出了各模型間數(shù)據(jù)聯(lián)系的具體方法，編寫了光機集成程序，據(jù)此進(jìn)行了動量輪擾振對成像質(zhì)量的全路徑分析，得到了某型號動量輪工作時對成像質(zhì)量的影響結(jié)果，可為后續(xù)的隔振設(shè)計提供參考。

2 集成分析技術(shù)途徑

一般顫振集成分析的基本技術(shù)途徑如圖1［4］所示。

圖1 成像質(zhì)量集成分析技術(shù)路徑圖Fig．1 Path to integrated analysis of imaging quality

集成分析中的關(guān)鍵，在于合理正確地建立相互聯(lián)系的各子系統(tǒng)的描述模型。其中的干擾源模型和控制模型，可以根據(jù)其自身特點直接集成在結(jié)構(gòu)動力學(xué)模型中，工程上實現(xiàn)起來較為容易。光學(xué)模型由于其學(xué)科跨度比較大，一般須要另外編寫調(diào)用程序來實現(xiàn)與其他模型間的數(shù)據(jù)交換。

3 全路徑分析鏈的建立

3．1 干擾源模型

Liu和Maghami等人對動量輪的振動特性進(jìn)行了深入的試驗研究，其研究表明［5］：動量輪的輸出擾動頻譜主要是以飛輪轉(zhuǎn)速的工頻成分為主，擾動模型可以利用頻率為飛輪轉(zhuǎn)速的正弦載荷來模擬［6］，載荷的幅值應(yīng)以實驗數(shù)據(jù)的均方根值為準(zhǔn)，這樣可以較好地模擬擾動能量。

3．2 結(jié)構(gòu)模型

結(jié)構(gòu)模型采用的是有限元模型，由于主鏡、次鏡、三鏡是重點研究的部件，后續(xù)研究會用到光學(xué)系統(tǒng)鏡面的平移和轉(zhuǎn)動數(shù)據(jù)，因此三個光學(xué)元件均建成了三維實體單元。圖2、圖3給出了相機有限元模型及裝配后的衛(wèi)星-相機有限元模型，圖中坐標(biāo)系為衛(wèi)星本體坐標(biāo)系，X軸為滾動軸，Y軸為俯仰軸，Z軸為偏航軸。表1給出了衛(wèi)星相機在軌模態(tài)的前十階頻率。

圖2 相機單機結(jié)構(gòu)有限模型示意圖Fig．2 Finite element model of the camera

圖3 裝配后的衛(wèi)星-相機有限元模型圖Fig．3 Finite element model of the camera and satellite

表1 相機單機和衛(wèi)星整星在軌有限元模態(tài)表Table 1 Modal analysis results of satellite-camera Hz

3．3 控制系統(tǒng)模型

控制系統(tǒng)在結(jié)構(gòu)系統(tǒng)上增加了閉環(huán)回路，限制了結(jié)構(gòu)系統(tǒng)三個方向轉(zhuǎn)動的剛體自由度，其作用相當(dāng)于增加了一個高通濾波器，使得低頻段的結(jié)構(gòu)響應(yīng)被抑制，有研究指出，控制系統(tǒng)可對1 Hz以下的擾振進(jìn)行補償［7］。假如不考慮控制系統(tǒng)的影響，得出的動力學(xué)響應(yīng)是漂移的，并不能準(zhǔn)確反映擾動對衛(wèi)星的影響。

雖然要考慮控制系統(tǒng)的影響，但是在集成分析中不需要詳細(xì)復(fù)雜的控制系統(tǒng)建模，動力學(xué)分析中可以采用模態(tài)疊加的方法來忽略低階頻率，這種方法簡單有效，已被多篇文獻(xiàn)證明［4，8］。

3．4 光學(xué)模型

光學(xué)模型采用codev軟件建立，codev軟件具有強大的光學(xué)建模及分析能力，支持使用指令語句對光學(xué)模型進(jìn)行修改，有利于編寫程序直接調(diào)用codev程序?qū)懭雽懗?，便于集成。利用codev軟件建立的某型高分辨率相機光學(xué)系統(tǒng)及其坐標(biāo)系示意圖如圖4，X軸為弧矢方向，垂直YZ平面向內(nèi)，Y軸為子午方向，弧矢方向和子午方向的光學(xué)特性會有差別。

圖4 某型高分辨率空間相機光路圖Fig．4 Optical system of space camera

顫振會引起視軸的漂移、抖動，造成光學(xué)系統(tǒng)自身結(jié)構(gòu)性的調(diào)制傳遞函數(shù)（MTF）下降，另外還造成主光學(xué)（chief ray）在焦平面的位移（像移），因此可以以MTF和像移來衡量成像質(zhì)量。由于光學(xué)元件本身的剛度很高（主鏡基頻可以到1700 Hz），振動造成的各鏡面面形變化可不予考慮［9］，即將鏡面視為剛體，主要分析各光學(xué)元件的偏心和傾斜所引起的光學(xué)系統(tǒng)MTF的下降和像移。利用codev軟件的偏心和回歸（DAR）功能來修改光學(xué)模型的參數(shù)，輸出修改后的MTF和像移結(jié)果。

3．5 光機接口程序

有限元動力學(xué)分析只能得出光學(xué)元件某個點的位移，光學(xué)模型需要各光學(xué)元件剛體六自由度位移，另外光學(xué)模型和有限元模型的坐標(biāo)系不一致，為了協(xié)調(diào)光學(xué)模型和有限元模型間的數(shù)據(jù)交換，作者開發(fā)了光機接口程序（JITTER）。JITTER程序可以讀入動力學(xué)分析的文本結(jié)果，計算光學(xué)系統(tǒng)坐標(biāo)系下的光學(xué)件六自由度位移，寫入codev光學(xué)模型，待codev計算結(jié)束后，寫出光學(xué)模型的MTF和像移。

光學(xué)元件六自由度位移的計算方法如圖5 所示，以主鏡為例，取有限元模型主鏡鏡面上＋X、-X、＋Y、-Y軸方向上最外圍的端點，可以利用有限元分析出不同載荷下這四個點的三個方向上的響應(yīng)平動位移，分別用（ΔxA，ΔyA，ΔzA）、（ΔxB，ΔyB，ΔzB）、（ΔxC，ΔyC，ΔzC）、（ΔxD，ΔyD，ΔzD）表示。這樣光學(xué)元件三個方向上的平動位移可以利4四個點在X 軸、Y 軸、Z軸三個方向上的平動位移的平均值來獲取，轉(zhuǎn)動位移為

式中：RX、RY、RZ分別是光學(xué)鏡在X、Y、Z方向上的轉(zhuǎn)動偏移量，單位是弧度，d是鏡子的直徑。

圖5 偏移量計算方法示意圖Fig．5 Method for calculating displacement

4 集成分析算例

按照3．1節(jié)的原則模擬某動量輪的輸出載荷，使用patran／nastran軟件對結(jié)構(gòu)模型進(jìn)行時域響應(yīng)分析。模擬載荷加載頻率為64 Hz，加載時間應(yīng)該長到足以反映動量輪的穩(wěn)態(tài)輸出，這里取20 個周期。計算步長應(yīng)該小到足以反映結(jié)構(gòu)所關(guān)心的最高頻率（本文相機的第十階模態(tài)頻率為150 Hz），這里取0．000 5s，可以精確反映500Hz內(nèi)的結(jié)構(gòu)振動，結(jié)構(gòu)阻尼按照一般經(jīng)驗暫取臨界值0．02，后續(xù)可根據(jù)試驗值來修正。分析結(jié)果經(jīng)JITTER 程序處理后，得到各光學(xué)元件六自由度位移隨時間t的變化如圖6所示。

圖6 動量輪擾振下各光學(xué)元件六自由度位移時域響應(yīng)圖Fig．6 Displacement of optical components under momentum wheel jitter

圖6是動量輪擾振下，控制系統(tǒng)不能控制的光學(xué)元件高頻抖動響應(yīng)。0～0．31s是載荷的加載時間，可以認(rèn)為是動量輪工作的時間。從圖6可以看出，動量輪工作時，空間相機各光學(xué)元件的平動偏移幅值為0．04～0．8μm，轉(zhuǎn)動位移偏移幅值為0．01″～0．08″。轉(zhuǎn)動位移可以反映相機的剛性（光學(xué)元件轉(zhuǎn)動位移相同說明顫振沒有引起相機內(nèi)部運動），從圖6中還可看出，各光學(xué)元件的轉(zhuǎn)動位移并不相同，甚至在X方向，次鏡的轉(zhuǎn)動位移遠(yuǎn)大于主鏡、三鏡的轉(zhuǎn)動位移，證明次鏡的響應(yīng)要遠(yuǎn)大于其他部件，微振動影響到了相機內(nèi)部，以往把相機作為剛體來分析成像質(zhì)量的方法值得商榷。

利用JITTER 程序?qū)?～0．3s內(nèi)的六自由度峰值偏移量寫入codev軟件，調(diào)用其中的分析功能，得到光學(xué)系統(tǒng)奈奎斯特頻率處各視場下的MTF最大變化值，見表2。

取各視場下的MTF 平均值作為指標(biāo)，得出光學(xué)系統(tǒng)自身的 MTF 下降很少，平均值下降0．213%。同時通過調(diào)用codev 中的spot diagram功能，得到相機（像元7μm）的像移量最大值為：俯仰方向1．28μm（0．183 像元），滾動方向0．85μm（0．12像元），偏航方向的偏移計算方法與滾動、俯仰不同，鑒于偏航方向的像移敏感量要遠(yuǎn)小于滾動、俯仰方向［10］，這里不予考慮，考慮到MTF指標(biāo)在相機設(shè)計中的重要性，可以計算像移造成的MTF 下降量。

表2 顫振造成的光學(xué)系統(tǒng)各視場MTF變化值Table 2 Change of optical system MTF by jitter

像移引起的MTF為［11］

式中：f是已知相機的奈奎斯特頻率；s是積分時間內(nèi)的像移。算出最大像移引起的奈奎斯特頻率處的MTF為0．986 2，下降了1．38%。

一般認(rèn)為當(dāng)顫振在焦面引起的TDICCD 像移不超過0．2 像元（對應(yīng)的奈奎斯特頻率處的MTF為0．983 6）時，顫振對成像質(zhì)量的影響是可以接受的［12］。動量輪對該型衛(wèi)星相機造成的最大像移為0．183像元，滿足圖像質(zhì)量要求。

5 結(jié)束語

本文為了分析動量輪顫振對成像質(zhì)量的影響，建立了干擾源模型、衛(wèi)星-相機結(jié)構(gòu)有限元模型、控制模型和光學(xué)模型，并給出了模型間聯(lián)系的具體方法，構(gòu)成了從擾振源到光學(xué)系統(tǒng)的全路徑分析鏈。通過對動量輪微振動的集成分析，發(fā)現(xiàn)在該型號動量輪發(fā)生擾振時，微振動影響到了相機內(nèi)部。另外，本文依據(jù)各光學(xué)元件的偏移得到了動量輪擾振對衛(wèi)星相機成像的具體影響。結(jié)果證明，在動量輪輸出擾振下，衛(wèi)星的成像質(zhì)量基本沒有受到影響，可以不對該動量輪或相機采取隔振措施。此分析方法也可為分析其它擾振源對成像質(zhì)量的影響提供一定的參考。

（References）

［1］Eyeman C E，Shea J F．A systems engineering approach to disturbance minimization for spacecraft utilizing controlled structures technology，MIT SERC Report＃2-90［R］．Boston：MIT，1990

［2］Ston C，Holtery C．The JWST integrated modeling environment［C］／／2004IEEE Aerospace Conference，New York：IEEE，2004：4041-4047

［3］龐世偉，楊雷，曲廣吉．高精度航天器微振動建模與評估技術(shù)最近進(jìn)展［J］．強度與環(huán)境，2007，34（6）：1-9

Pang Shiwei，Yang Lei，Qu Guangji．New development of micro-vibration integrated modeling and assessment technology for high performance spacecraft［J］．Structure ＆Environment Engineering，2007，34（6）：1-9（in Chinese）

［4］龐世偉，潘騰．微振動對圖像質(zhì)量影響評估方法研究［C］／／高分辨率遙感衛(wèi)星結(jié)構(gòu)振動及控制技術(shù)研討會．北京：中國宇航學(xué)會，2011：163-170

Pang Shiwei，Pan Teng．The method to analysis of impact by micro-vibration on imaging［C］／／Structural Vibration and Control of High Resolution Remote Sensing Satellites．Beijing：Chinese Society of Astronautics，2011：163-170（in Chinese）

［5］Liu K，Maghami P，Blaurock C．Recaction wheel disturbance modeling，jitter analysis，and validation tests for solar dynamics observatory，AIAA 2008-7232［R］．Washington D．C．：AIAA，2008

［6］Rudoler S，Hadar O，F(xiàn)isher M，et al．Image resolution limits resulting from mechanical vibrations，Part3：numerical calculation of modulation transfer function［J］．Optical Engineering，1992，31（3）：581-589

［7］Toyoshima Morio，Takashi Jono，Takahashi Nobuhiro．Transfer functions of micro-vibration disturbances on a satellite［R］．Tsukuba ：National Space Development Agency of Japan，2006

［8］史紀(jì)鑫，于偉．基于全柔性模型的高分衛(wèi)星微振動分析方法［C］／／高分辨率遙感衛(wèi)星結(jié)構(gòu)振動及控制技術(shù)研討會論文集．北京：中國宇航學(xué)會，2011：210-218

Shi Jixin，Yu Wei．Method to analysis of micro-vibration based on soft model［C］／／Proc．of Structural Vibration and Control of High Resolution Remote Sensing Satellites．Beijing：Chinese Society of Astronautics，2011：210-218（in Chinese）

［9］張博文，王小勇，胡永力．微振動對高分辨率空間相機成像影響的集成分析［J］．航天返回與遙感，2012，33（2）：60-65

Zhang Bowen，Wang Xiaoyong，Hu Yongli．Integrated analysis on effect of micro-vibration on high resolution space camera imaging［J］．Spacecraft Recovery and Remote Sensing，2012，33（2）：60-65（in Chinese）

［10］樊超，李英才，易紅偉．顫振對TDICCD 相機像質(zhì)的影響分析［J］．光子學(xué)報，2007，36（9）：36-44

Fan Chao，Li Yingcai，Yi Hongwei．Influence analysis of buffeting on image quality of TDICCD camera［J］．Acta Photonica Sinica，2007，36（9）：36-44（in Chinese）

［11］Rudoler S，Hadar O，F(xiàn)isher M．Image resolutions limits resulting from mechanical vibration［J］．Optical Engineering．1991，30（5）：577-589

［12］Holst G C．CCD arrays cameras and displays［J］．SPIE Optical Engineering，1998（2272）：332-331