梁偉，秦開宇，蔣仁奎，劉盾

(1.電子科技大學航空航天學院 成都 611731；2.中國科學院光電技術研究所 成都 610209)

巡天望遠鏡主鏡口徑達到2m，是中國第一個大口徑、大視場空間天文望遠鏡。巡天相機是巡天望遠鏡最重要的科學終端載荷之一。巡天相機焦面尺寸達到610mm×590mm，采用機械快門對焦平面進行曝光控制。在軌工作10 年，曝光次數(shù)約為100 萬次，壽命長、可靠性高?？扉T采用結(jié)構(gòu)簡單、電機直驅(qū)的旋轉(zhuǎn)對開結(jié)構(gòu)形式。機械快門運行的旋轉(zhuǎn)運動將形成微振動，產(chǎn)生激振力。同時方向相反的兩旋轉(zhuǎn)運動帶來的激振力在轉(zhuǎn)軸相對方向的分量抵消，在兩轉(zhuǎn)軸平面法向的分量加倍，成為巡天相機的主要擾動源之一，降低空間光學望遠鏡的成像質(zhì)量，需對其進行抑制[1-2]。

類似旋轉(zhuǎn)對開式機構(gòu)在空間望遠鏡上應用較多，如哈勃望遠鏡(Hubble space telescope)[3-4]、錢德拉X 射線太空望遠鏡(Chandra X-ray observatory)[5-7]、倫琴X 射線天文臺(Rosat)[8-9]和XMM 牛頓衛(wèi)星(X-ray multi-mirror mission)[10-11]的遮光蓋系統(tǒng)。這些遮光系統(tǒng)的功能主要是避免望遠鏡內(nèi)部器件受到強光照射，工作過程中望遠鏡焦面組件未工作，運行中的激振力不影響望遠鏡成像質(zhì)量。歐幾里得空間望遠鏡(Euclid space mission)[12-14]中可見光成像儀(Euclid-VIS)采用的快門形式與巡天相機機械快門類似，通過在葉片后端進行嚴格的靜動平衡設計，抑制激振力，取得了很好的抑制效果。巡天相機機械快門由于空間包絡約束限制難以采取靜、動平衡措施減小激振力?？扉T的運行頻率f小于1Hz，采用傳統(tǒng)的隔振措施，隔振頻率小于0.707f時才會發(fā)揮隔振效果[15-16]，這將削弱快門的支撐剛度，穩(wěn)定性差。

針對旋轉(zhuǎn)對開式機械快門的激振力問題，建立轉(zhuǎn)子及機械快門運動數(shù)學模型，得到快門運動的激振力數(shù)值解。在此基礎上，提出一種動態(tài)補償方法用于激振力抑制：結(jié)合凸輪運動動態(tài)性和杠桿反向補償、放大原理，且杠桿慣性力與快門激振力相等，實現(xiàn)激振力動態(tài)補償抑制。采用單純形法進行凸輪曲線優(yōu)化設計。最后進行運動學仿真對激振力的動態(tài)補償抑制方法進行驗證。

1 巡天相機機械快門及其布局

如圖1 所示，機械快門位于焦面組件前端，與焦面組件一同安裝在主體結(jié)構(gòu)上，快門兩葉片同時開合實現(xiàn)對焦平面的曝光控制。

圖1 巡天相機及機械快門布局

機械快門由正Y側(cè)、負Y側(cè)以及支座組成，正Y側(cè)和負Y側(cè)快門兩葉片同步反向轉(zhuǎn)動，完成快門的開啟或關閉功能。

2 數(shù)學建模

2.1 旋轉(zhuǎn)轉(zhuǎn)子數(shù)學建模

轉(zhuǎn)子繞轉(zhuǎn)軸旋轉(zhuǎn)，假設其滿足如下條件：

1)轉(zhuǎn)子與轉(zhuǎn)軸均為剛性體，為薄圓盤轉(zhuǎn)子；

2)在均質(zhì)圓盤上附加偏心質(zhì)量為m，偏心距為r。

為建立轉(zhuǎn)子數(shù)學模型，采用如圖2 所示的直角坐標系oxyz，轉(zhuǎn)軸與坐標系x軸重合。

圖2 圓盤轉(zhuǎn)子模型

轉(zhuǎn)子圓盤形心為O點，質(zhì)心為C點，其坐標關系為：

式中，質(zhì)心坐標為(xc，yc)，對上式求導數(shù)，得到偏心質(zhì)量的速度與加速度如下：

由于存在偏心質(zhì)量m，轉(zhuǎn)子的激振力如下：

2.2 機械快門運動數(shù)學模型

機械快門簡化模型如圖3 所示，規(guī)定葉片與水平面夾角θ 為葉片轉(zhuǎn)角。當θ=0°時，快門完全關閉，當θ=90°時，快門完全打開。

圖3 巡天相機機械快門簡化模型

將葉片偏心簡化為質(zhì)點，附加質(zhì)量為m，質(zhì)心與轉(zhuǎn)軸的距離為r。

據(jù)圖3，負Y側(cè)快門的數(shù)學模型如下：

同樣，正Y側(cè)快門的數(shù)學模型：

則機械快門的整體振動受力為：

即旋轉(zhuǎn)對開式機械快門在Y向激振力相互抵消，為0；在Z向的激振力為單側(cè)旋轉(zhuǎn)運動激振力的兩倍。

2.3 機械快門運動分析結(jié)果

據(jù)巡天任務規(guī)劃，快門典型運行過程：1)打開和關閉過程時間為1.3s；2)保持打開持續(xù)時間為150～200s。為保證快門開啟和關閉時刻的平穩(wěn)性，葉片采用正弦波加減速驅(qū)動模式，即：

葉片運動曲線如圖4 所示。葉片相對轉(zhuǎn)軸附加質(zhì)量m為1.4kg，葉片質(zhì)心與轉(zhuǎn)軸的距離r為144.75mm，將式(10)～式(12)代入式(9)，得到一個曝光周期下快門的Z向激振力，激振力幅值達到2.38N。

圖4 快門葉片驅(qū)動模式

3 法向激振力補償抑制

3.1 補償抑制方法

杠桿工作原理如圖5 所示。

圖5 杠桿工作原理圖

當質(zhì)點m1沿正Y向以加速度a1運動時，受到的慣性力F1為：

根據(jù)杠桿原理，質(zhì)點m2沿負Y向以加速度a2運動，受到的慣性力F2，兩者關系為：

由式(13)、式(15)知，杠桿加減速過程中，兩側(cè)角加速度大小相等，方向相反。杠桿兩側(cè)端點處線加速度方向仍相反，大小與支點距離成正比。

由式(14)、式(15)知，當快門轉(zhuǎn)軸帶動質(zhì)點m1向上運動時，杠桿另一側(cè)質(zhì)點m2向下運動，慣性力F2與F1方向相反，可以補償快門在F1方向上產(chǎn)生的激振力。杠桿在放大質(zhì)點m2線加速度的同時，取m1<>abs(F1)，可以用小的質(zhì)量代價獲取大的慣性力F2補償快門的激振力，適用于航天任務。

3.2X 向附加分量

由圖5 所示，杠桿在補償機械快門Y向激振力的同時，將產(chǎn)生X向的水平分力。

其合力為：

式中，m1<

考慮使用兩套機構(gòu)，對稱布置在旋轉(zhuǎn)軸下方，如圖6 所示，不占用巡天相機的其他空間。

圖6 補償抑制機構(gòu)布局

當兩機構(gòu)完全一致時，負X側(cè)的機構(gòu)產(chǎn)生正X向的附加分量，正X側(cè)的機構(gòu)產(chǎn)生負X向的附加分量，兩者大小一致，相互抵消，克服補償抑制X向附加分量不為零的問題。

4 補償抑制機構(gòu)

4.1 機構(gòu)設計

兩組凸輪杠桿機構(gòu)位于單側(cè)快門基座下方，如圖7 所示，每套機構(gòu)均由平衡塊、杠桿、轉(zhuǎn)軸、凸輪和導向桿組成。快門軸系與凸輪相連，無需額外驅(qū)動，軸系帶動凸輪轉(zhuǎn)動，凸輪驅(qū)動導向桿移動，導向桿推動杠桿一端使其繞轉(zhuǎn)軸轉(zhuǎn)動，杠桿另一端與平衡塊相連，杠桿轉(zhuǎn)動同時帶動平衡塊實現(xiàn)相應運動。

4.2 凸輪曲線優(yōu)化

據(jù)圖7 中激振力抑制機構(gòu)的零件幾何位置關系，假設轉(zhuǎn)軸中心到導桿中心線的距離為A，到平衡塊質(zhì)心的距離為B；平衡塊質(zhì)量為m0；凸輪對導向桿的驅(qū)動位移量為k(t)。

圖7 激振力抑制機構(gòu)剖視圖

凸輪驅(qū)動導桿移動時，k(t)隨時間變化，杠桿另一端的平衡塊帶來的慣性力為：

式中，m0=2kg；A=35mm；B=180mm。

為了補償快門Z向的激振力，需滿足：

將式(9)～式(12)、式(18)代入式(19)，即可形成凸輪驅(qū)動位移與時間的關系式，但該關系式復雜，無法求出k(t)的顯性解析式，采用數(shù)值優(yōu)化方式對其進行求解。

將k(t)表示成多項式：

k(t)的約束條件為：

3 個約束條件分別表示：1)快門處于關閉位置時，杠桿平行于快門軸系；2)快門啟動時刻凸輪驅(qū)動導向桿的速度為零；3)快門關閉時凸輪驅(qū)動導向桿的速度為零。

將快門激振力和平衡塊慣性力最大差異的絕對值最小作為優(yōu)化目標，即：

由于優(yōu)化目標值無法求導，常規(guī)的梯度優(yōu)化方法如坐標輪換法、最速下降法、牛頓法、共軛梯度法和變尺度法均不適用[17]。單純形法作為非梯度優(yōu)化方法[18]，收斂速度快，穩(wěn)定性好，變量個數(shù)不受限，可通過變換初始頂點避免局部最優(yōu)，適用于凸輪曲線的優(yōu)化。為避免陷入局部最優(yōu)解，在可行域內(nèi)隨機初始頂點，并多次計算，選擇所有計算結(jié)果中的最優(yōu)值。

將式(9)、式(18)和式(20)代入式(22)中，采用單純形法對多項式中各個參數(shù)進行優(yōu)化，使M最小。隨著優(yōu)化變量個數(shù)增多，k(t)自由度增加，目標函數(shù)最優(yōu)解越來越小，迭代次數(shù)和計算時間增加。當目標值變化率趨近于零，停止優(yōu)化迭代，迭代過程如圖8 所示，其中優(yōu)化后目標值未超過0.26N，優(yōu)化變量個數(shù)為8；將式(21)代入式(20)得到3 個方程可以解得3 項系數(shù)值，這樣式(20)共11 個系數(shù)項，次數(shù)n為10。

圖8 目標優(yōu)化過程

優(yōu)化所得多項式的各系數(shù)值如表1 所示，即凸輪驅(qū)動曲線是時間變量的10 次多項式。

表1 多項式系數(shù)

4.3 補償抑制誤差

將各變量代入式(20)，聯(lián)合式(9)和式(12)，得到凸輪轉(zhuǎn)動角度與驅(qū)動位移量的關系，如圖9所示。

圖9 凸輪曲線

得到凸輪驅(qū)動位移曲線k(t)和線加速度曲線(t)后，可以計算慣性力與激振力之間補償誤差如下：

補償抑制誤差曲線如圖10 所示，最大值為0.26N，為最大激振力的10.9%。

圖10 補償誤差曲線

4.4 工程可行性分析

補償抑制機構(gòu)是一被動運動機構(gòu)，不需要增加額外的電學控制系統(tǒng)，但會增加結(jié)構(gòu)的復雜性：包括凸輪和轉(zhuǎn)軸等，均為工程上的常用機構(gòu)，且這些機構(gòu)在航空航天任務中均得到應用[19]，是可行的。對補償抑制機構(gòu)進行質(zhì)量統(tǒng)計：單側(cè)快門兩平衡塊質(zhì)量小計2.0kg、機構(gòu)其他部分質(zhì)量為0.58kg，所以機械快門兩側(cè)機構(gòu)共計5.16kg。而機械快門本體質(zhì)量為17.8kg。增加補償抑制機構(gòu)后的機械快門質(zhì)量共計22.96kg，增加了29.0%。其中任務總體分配給機械快門的質(zhì)量要求不超過25kg，較機械快門質(zhì)量仍有余量2.04kg，滿足任務指標要求。因此，采用補償抑制機構(gòu)來解決機械快門微振動抑制問題，除了增加質(zhì)量代價外，沒有對望遠鏡正常工作產(chǎn)生其他影響，在工程上是可行的。

5 運動學分析驗證

建立快門運動學模型對其進行瞬態(tài)分析，模擬快門運行過程。使用旋轉(zhuǎn)副模擬軸系軸承，采用接觸力模擬凸輪與導桿的驅(qū)動關系，移動副約束導桿的自由度。在快門支座下建立基準，將基準與支座下表面耦合，固定基準即可對整個模型進行約束，通過讀取該基準的支反力獲得快門運行時的激振力。

在激振力補償抑制機構(gòu)使用前后的Z向激振力響應如圖11 所示。抑制后激振力幅值最大值由2.38N 降至0.29N(與優(yōu)化設計得到的最大值0.26 N 基本一致)，抑制率87.8%，驗證了本文提出的補償抑制方法。

圖11 抑制前和抑制后的Z 向激振力

6 結(jié)束語

針對巡天相機旋轉(zhuǎn)對開式機械快門提出一種激振力動態(tài)補償抑制方法。建立快門的運動數(shù)學模型，設計激振力補償抑制裝置，凸輪和杠桿機構(gòu)相結(jié)合實現(xiàn)杠桿質(zhì)量點慣性力對激振力的補償抑制。由單純形法優(yōu)化獲取凸輪曲線，最后進行機械快門運動學仿真，結(jié)果表明：在激振力抑制裝置采用前后，快門激振力分別為2.38N 和0.29N，抑制率87.8%，效果明顯，解決了在尺寸約束下的快門激振力抑制問題。