董得義，辛宏偉，楊利偉，李志來，關(guān)英俊

（1.中國科學(xué)院 長春光學(xué)精密機械與物理研究所，長春 130033；2.中國科學(xué)院 研究生院，北京 100039；3.長春工業(yè)大學(xué) 機電工程學(xué)院，長春 130012）

為了滿足空間光學(xué)感遙器高分辨、寬覆蓋、大視場［1，2］的要求，反射鏡孔徑越做越大，質(zhì)量越做越重，這給反射鏡支撐結(jié)構(gòu)的設(shè)計帶來的極大的困難。

當(dāng)前，各國空間光學(xué)遙感器的研制過程中，大孔徑反射鏡的支撐技術(shù)都是空間光學(xué)遙感器的關(guān)鍵技術(shù)之一［3，4］。

由于光學(xué)遙感器在運載、發(fā)射以及在軌運行過程中，需要承受各種動力學(xué)環(huán)境以及熱環(huán)境的交互作用［5］，這些環(huán)境約束的存在使得光學(xué)遙感器反射鏡支撐結(jié)構(gòu)的設(shè)計處于一個矛盾的境地：一方面，為了避免在運載、發(fā)射過程中反射鏡組件出現(xiàn)結(jié)構(gòu)變形、微屈服、疲勞破壞等現(xiàn)象，需要反射鏡的支撐結(jié)構(gòu)具有足夠的剛度、強度特性；而另一方面，由于反射鏡組件中各種材料的線脹系數(shù)的不匹配［6］，為了避免熱變形導(dǎo)致鏡面面形產(chǎn)生相對位置以及光學(xué)元件的波前畸變，影響光學(xué)系統(tǒng)的成像質(zhì)量，又需要反射鏡的支撐結(jié)構(gòu)具有一定的柔性，通過支撐結(jié)構(gòu)在熱載荷中產(chǎn)生較大的變形達到卸載和吸收由熱變形引起的應(yīng)變能，從而達到保護鏡面不會產(chǎn)生較大的畸變的目的。由此可見，如何保證反射鏡組件具有足夠的動態(tài)特性、較高的熱尺寸穩(wěn)定性，反射鏡支撐結(jié)構(gòu)的設(shè)計成為解決問題的關(guān)鍵［7］。

如何保證大孔徑反射鏡的支撐結(jié)構(gòu)在具有一定柔性的基礎(chǔ)上，順利通過鑒定級的隨機振動試驗成為了結(jié)構(gòu)設(shè)計人員亟需解決的問題。據(jù)國外統(tǒng)計，衛(wèi)星發(fā)射上天后，第一天所出現(xiàn)的故障，有30% ～60%是由動力學(xué)環(huán)境所引起的［8］，動力學(xué)環(huán)境對星上載荷結(jié)構(gòu)造成的損壞已經(jīng)成為星上載荷結(jié)構(gòu)失效的主要因素之一，因此對星上載荷的支撐結(jié)構(gòu)進行動力學(xué)響應(yīng)分析，根據(jù)分析結(jié)果對結(jié)構(gòu)進行優(yōu)化設(shè)計已經(jīng)成為提高產(chǎn)品可靠性的關(guān)鍵技術(shù)手段［9］。

本文對某光學(xué)遙感器的大孔徑反射鏡組件進行了隨機振動響應(yīng)分析，預(yù)示了反射鏡柔性支撐結(jié)構(gòu)上的應(yīng)力響應(yīng)，并在反射鏡組件的隨機振動試驗中對反射鏡柔性支撐結(jié)構(gòu)上的應(yīng)力響應(yīng)進行了測量，通過分析結(jié)果與試驗結(jié)果的對比，驗證了分析方法的有效性。

1 隨機振動響應(yīng)分析的基本原理

在隨機振動試驗中，任一瞬時各種頻率成分的激勵同時作用在反射鏡組件上，因此無法預(yù)測其瞬時的振動變化規(guī)律，一般用整個加載周期內(nèi)的均值、均方根、累計均方根和功率譜密度等統(tǒng)計量研究反射鏡在隨機振動試驗中的響應(yīng)［10］。

本文采用MSC.NASTRAN進行反射鏡組件的隨機振動響應(yīng)分析，在MSC.NASTRAN中把隨機響應(yīng)分析當(dāng)作頻率響應(yīng)的后處理來進行［11］。其輸入文件不但包括頻率響應(yīng)分析的輸出結(jié)果，而且包括用戶以功率和相關(guān)譜密度形式提供的載荷情況。輸出結(jié)果為響應(yīng)的功率譜密度、自相關(guān)函數(shù)、每單位時間的正斜率方向的零交點個數(shù)、以及響應(yīng)的均方根值（RMS值）等等。

在直接頻率響應(yīng)分析中，通過用復(fù)數(shù)代數(shù)算法求解一系列耦合的矩陣方程，計算離散激勵頻率結(jié)構(gòu)響應(yīng)。對簡諧激勵下有阻尼強迫振動，運動方程為：

對簡諧運動（頻率響應(yīng)分析的基礎(chǔ)），假定一個簡諧形式的解：

其中，{u（ω）}為復(fù)位移向量。

對式（2）求一階及二階導(dǎo)數(shù)得：

將上述表達式代入式（1）得：

除以 eiωt簡化為：

利用復(fù)數(shù)算法，對于每一個輸入激勵頻率的運動方程，可以像靜力學(xué)問題類似地求解。

2 大孔徑反射鏡組件的有限元模型

2.1 大孔徑反射鏡組件的結(jié)構(gòu)

大孔徑反射鏡組件結(jié)構(gòu)的剖面圖如圖1所示，主要由反射鏡、柔性支撐以及三角板等構(gòu)件組成。

柔性支撐結(jié)構(gòu)如圖2所示。

圖1 反射鏡組件結(jié)構(gòu)剖面圖Fig.1 Cross-section view of large mirror component

2.2 大孔徑反射鏡組件的有限元模型

根據(jù)零件結(jié)構(gòu)的特點，采用 8節(jié)點六面體單元對各零件進行有限元網(wǎng)格劃分，個別拐角的地方采用五面體楔形單元過渡，整個組件的有限元模型如圖3所示。

在反射鏡組件的有限元模型中共劃分了124 150個單元，180 020個節(jié)點。在大孔徑反射鏡組件的建模過程中主要考慮了以下幾點：

圖2 反射鏡柔性支撐結(jié)構(gòu)模型Fig.2 The structure model of the flexible support

圖3 反射鏡組件有限元模型Fig.3 The finite element model of the mirror component

（1）有限元模型按照設(shè)計人員提供的三維模型構(gòu)造，以保證有限元模型與設(shè)計結(jié)構(gòu)一致；

（2）在關(guān)鍵的力的傳遞路徑上（比如：柔性支撐結(jié)構(gòu)）：網(wǎng)格劃分較密；非關(guān)鍵部位網(wǎng)格劃分可適當(dāng)稀疏一些，反射鏡柔性支撐的有限元模型如圖4所示。

（3）反射鏡組件中的螺釘連接采用桿單元與RBE2多點約束進行模擬。

圖4 柔性支撐的有限元模型Fig.4 The finite element model of flexible support

有限元模型中所用到的主要材料參數(shù)見表1。

表1 材料主要參數(shù)Tab.1 Main parameters for materials

3 反射鏡組件結(jié)構(gòu)阻尼的測定

采用直接法進行頻響分析過程中，需要輸入結(jié)構(gòu)阻尼。阻尼代表結(jié)構(gòu)中能量耗散的一種數(shù)學(xué)近似。阻尼由多種機制引起，例如：粘性效應(yīng)；外摩擦（結(jié)構(gòu)結(jié)合部之間的滑動）；內(nèi)摩擦（材料類型的特點）；結(jié)構(gòu)的非線性（塑性、間隙）。因為這些影響很難定量，所以阻尼很難準(zhǔn)確模擬，阻尼經(jīng)?；趧恿υ囼灥慕Y(jié)果進行分析計算。

阻尼是許多復(fù)雜機制引起的結(jié)果，在計算響應(yīng)中，阻尼的影響有賴于動力分析的類型及載荷持續(xù)的時間。在短時間載荷作用下，阻尼經(jīng)?？梢院雎?，例如沖擊載荷的作用，因為能量還沒有來得及明顯耗散結(jié)構(gòu)就已經(jīng)達到了它的峰值響應(yīng)。但對長時間載荷，阻尼非常重要，因為這些載荷連續(xù)不斷的向結(jié)構(gòu)加入能量。

在MSC.NASTRAN中，采用直接法進行頻響分析時，其阻尼用阻尼矩陣代表，合并在剛度矩陣中。其阻尼矩陣由兩部分組成，一部分為阻尼單元（粘性阻尼），一部分結(jié)構(gòu)阻尼。因為求解域為復(fù)數(shù)域，允許剛度矩陣為復(fù)數(shù)，所以不必將結(jié)構(gòu)阻尼轉(zhuǎn)換為等效的粘性阻尼形式。對沒有引入阻尼單元的結(jié)構(gòu)，主要的阻尼輸入就是結(jié)構(gòu)阻尼。而對于結(jié)構(gòu)阻尼的輸入，一方面在材料屬性輸入時，會指定用這種材料的所有單元的結(jié)構(gòu)阻尼；另一方面有一個直接輸入的結(jié)構(gòu)阻尼系數(shù)，這個結(jié)構(gòu)阻尼系數(shù)G與結(jié)構(gòu)的動力放大因子Q、結(jié)構(gòu)的阻尼比ζ有如下關(guān)系：

對于結(jié)構(gòu)的阻尼比ζ的計算一般是通過模態(tài)測試試驗獲?。ü逃蓄l率、主振型、阻尼比等），此外還能通過幅頻特性曲線，采用半功率法計算獲取，半功率法計算阻尼的原理［12］如圖5所示。

幅頻特性曲線的共振峰兩側(cè)，可各找到一幅值為0.707Hm的點P1和P2，稱為半功率點。這是因為在這兩點時，系統(tǒng)的振動能量是共振時能量的一半。在小阻尼時，P1和P2對應(yīng)的頻率為：

圖5 半功率法計算阻尼原理圖Fig.5 Acceleration response curve of 0.2 g input

令半功率帶寬為：

故得：

本文中，先對反射鏡組件進行了0.2 g特征掃頻試驗，通過反射鏡組件的幅頻特性曲線，采用半功率法獲取了大孔徑反射鏡組件的結(jié)構(gòu)阻尼。

本文只對水平沿鏡頭長邊方向的分析結(jié)果與試驗結(jié)果進行對比研究（其他兩個方向的情況與之類似），因此，在下面的分析中只對這一個方向進行討論（下同）。

反射鏡組件0.2 g，10 Hz～2 000 Hz頻段的特征掃頻曲線如圖6所示。

圖6 0.2 g掃頻加速度響應(yīng)曲線Fig.6 Acceleration response curve of 0.2 g input

通過半功率法計算得到，反射鏡組件在0.2 g載荷條件下，沿鏡頭長邊方向的結(jié)構(gòu)阻尼為0.015。

4 隨機振動響應(yīng)分析

4.1 基礎(chǔ)隨機激勵載荷譜

空間光學(xué)遙感器隨機振動環(huán)境以基礎(chǔ)隨機激勵的加速度功率譜密度形式給出，對于反射鏡組件隨機振動載荷作用于反射鏡三角板的底部。根據(jù)反射鏡三角板與遙感器平臺對接面上載荷譜的包絡(luò)線和實際的電磁振動臺激勵頻率限制的條件，施加的載荷譜如圖7所示，頻率范圍為10 Hz～2 000 Hz，輸入條件的均方根加速度為6.0 Grms。

圖7 隨機振動試驗條件Fig.7 The test condition of random vibration

4.2 隨機振動加速度響應(yīng)分析結(jié)果

圖8、圖9分別為沿鏡頭長邊方向，反射鏡鏡頭上的加速度功率譜密度響應(yīng)與累積均方根值。加速度功率譜密度均方根響應(yīng)為16.0 Grms，放大倍數(shù)為2.7倍。

4.3 均值應(yīng)力響應(yīng)分析結(jié)果

圖10、圖11分別為沿鏡頭長邊方向激勵下，柔性支撐上的應(yīng)力均方根響應(yīng)與累積均方根值。反射鏡柔性支撐結(jié)構(gòu)上的均值應(yīng)力響應(yīng)分析結(jié)果為34.9 MPa。

5 試驗結(jié)果

5.1 試驗裝置

為了驗證反射鏡支撐結(jié)構(gòu)設(shè)計是否合理，并驗證分析結(jié)果的正確性，對反射鏡組件進行了隨機振動試驗，試驗現(xiàn)場如圖12所示。

圖12 反射鏡組件隨機振動試驗現(xiàn)場Fig.12 Placement of random vibration test

5.2 試驗結(jié)果與分析

反射鏡組件沿鏡頭長邊方向的均方根加速度響應(yīng)曲線如圖13所示。均方根加速度響應(yīng)為16.3 Grms。

反射鏡柔性支撐結(jié)構(gòu)在隨機振動試驗中的應(yīng)變響應(yīng)曲線如圖14所示。

圖13 加速響應(yīng)功率譜密度（PSD）試驗結(jié)果Fig.13 The test result of power spectrum density（PSD）of the acceleration response

圖14 隨機振動柔性支撐結(jié)構(gòu)應(yīng)變響應(yīng)曲線Fig.14 Strain response curves of flexible support under random vibration

從試驗獲得的結(jié)果曲線可以看出，反射鏡組件的均方根加速度響應(yīng)為16.3 Grms，加速度放大倍數(shù)為2.7倍；柔性支撐在隨機振動試驗中的均值應(yīng)力響應(yīng)為31.3 MPa。試驗結(jié)果與分析結(jié)果的對比分析表明：均方根加速度響應(yīng)分析誤差為1.8%，均值應(yīng)力分析誤差為13.2%，滿足工程的精度要求。

由于隨機振動測試結(jié)果存在一定的誤差，會影響分析結(jié)果精度的計算，若要更加準(zhǔn)確的評價分析結(jié)果的精度，須提高測試結(jié)果的精度。關(guān)廣豐等人［13］的研究結(jié)果表明，可通過功率譜再現(xiàn)的自適應(yīng)算法控制振動臺，提高振動臺的控制精度；游偉倩等人［14］提出了采用雙振動臺雙自由度綜合控制的算法，能將隨機振動的的輸出控制譜和參考普的誤差控制在±1 dB以內(nèi)。如何保證測試結(jié)果的精度是驗證分析結(jié)果準(zhǔn)確與否一個關(guān)鍵。

6 結(jié)論

大孔徑反射鏡的支撐結(jié)構(gòu)在具有一定柔性的同時，能夠抵抗發(fā)射段隨機振動的干擾，是空間光學(xué)遙感器入軌后成像質(zhì)量的重要保證。通過有限元分析技術(shù)，對反射鏡組件進行了隨機振動響應(yīng)分析，可以預(yù)示反射鏡柔性支撐結(jié)構(gòu)經(jīng)受動力學(xué)環(huán)境的能力，并驗證大孔徑反射鏡支撐結(jié)構(gòu)設(shè)計是否合理，為反射鏡支撐結(jié)構(gòu)的優(yōu)化設(shè)計提供了理論支持；同時能有效的指導(dǎo)反射鏡組件的隨機振動試驗的進行，預(yù)防在環(huán)境試驗中出現(xiàn)質(zhì)量事故。

本文利用有限元結(jié)構(gòu)分析軟件對某空間光學(xué)遙感器大孔徑的反射鏡組件進行了鑒定級的隨機振動響應(yīng)分析，得到了反射鏡組件的均方根加速度響應(yīng)與柔性支撐結(jié)構(gòu)在隨機振動試驗中的均方根應(yīng)力響應(yīng)，試驗結(jié)果表明：分析結(jié)果與試驗結(jié)果基本吻合，分析結(jié)果對大孔徑反射鏡組件的優(yōu)化設(shè)計具有一定的指導(dǎo)意義。

［1］李志來，薛棟林，張學(xué)軍.長焦距大視場光學(xué)系統(tǒng)的光機結(jié)構(gòu)設(shè)計［J］.光學(xué)精密工程，2008，16（12）：2485 －2490.

［2］韓昌元.高分辨力空間相機的光學(xué)系統(tǒng)研究［J］.光學(xué)精密工程，2008，16（11）：2164 －2172.

［3］鮑 赫，李志來.長條形SiC反射鏡輕量化及支撐結(jié)構(gòu)的設(shè)計［J］.光學(xué)技術(shù)，2008，34（4）：593 －596.

［4］閆 勇，王 棟，金 光.大孔徑SiC反射鏡組件研制技術(shù)［J］.光電工程，2010，37（6）：108 －112.

［5］郭 疆.大孔徑空間遙感相機主反射鏡支撐設(shè)計［J］.光學(xué) 精密工程，2008，16（9）：1642 －1647.

［6］王忠素，翟 巖，梅 貴，等.空間光學(xué)遙感器反射鏡柔性支撐的設(shè)計［J］.光學(xué) 精密工程，2010，18（8）：1833－1841.

［7］孫寶玉.光學(xué)反射鏡柔性支撐結(jié)構(gòu)尺寸穩(wěn)定性分析［J］.光電工程，2009，36（9）：142 －145.

［8］柯受全.衛(wèi)星環(huán)境工程和模擬試驗（下）［M］.北京：宇航出版社，1996，23 －73.

［9］韓 旭，吳清文，董得義，等.室溫硫化膠層建模在透鏡結(jié)構(gòu)分析中的應(yīng)用［J］.光學(xué) 精密工程，2010，18（1）：118－125.

［10］張 軍，諶 勇，張志誼，等.衛(wèi)星隨機試驗的振動響應(yīng)分析［J］.機械強度，2006，28（1）：16－19.

［11］MSC.Software Corporation.MSC.Nastran動力學(xué)培訓(xùn)教程［M］.

［12］劉晶波，杜修力.結(jié)構(gòu)動力學(xué)［M］.北京：機械工業(yè)出版社，2005，60 －61.

［13］游偉倩，崔旭利，陳懷海，等.雙振動臺隨機振動H∞雙自由度綜合控制研究［J］.振動與沖擊，2011，30（5）：179－183，204.

［14］關(guān)廣豐，王海濤，熊 偉.隨機振動功率譜再現(xiàn)自適應(yīng)控制算法研究［J］.振動與沖擊，2011，30（3）：01 －04.