基于二級多點逼近算法的航天器結(jié)構(gòu)優(yōu)化設(shè)計

陳珅艷 黃 海

(北京航空航天大學(xué) 宇航學(xué)院,北京 100191)

基于二級多點逼近算法的航天器結(jié)構(gòu)優(yōu)化設(shè)計

陳珅艷 黃 海

(北京航空航天大學(xué) 宇航學(xué)院,北京 100191)

在某新型衛(wèi)星結(jié)構(gòu)研制的初樣階段,建立了衛(wèi)星初始設(shè)計的有限元模型.以主承力桁架結(jié)構(gòu)中各梁的截面尺寸為設(shè)計變量,考慮整星模態(tài)頻率和強度約束,建立了以結(jié)構(gòu)重量最輕為目標(biāo)的優(yōu)化模型.應(yīng)用二級多點逼近優(yōu)化算法進行了結(jié)構(gòu)優(yōu)化計算,在每一個計算周期中,原結(jié)構(gòu)優(yōu)化問題先轉(zhuǎn)化為具有較高精度的第 1級多點近似問題,該問題則通過可由對偶法快速求解的第 2級近似問題逼近.經(jīng)過初步試算、設(shè)計改進和再優(yōu)化 3個階段,設(shè)置并進行了一系列的優(yōu)化計算,逐步明確了設(shè)計的方向和各結(jié)構(gòu)參數(shù)的取值范圍,得到了合理的可行方案,為該衛(wèi)星平臺結(jié)構(gòu)初樣的詳細(xì)設(shè)計提供了參考,同時表明所采用的優(yōu)化方法適用于工程結(jié)構(gòu)優(yōu)化問題.

結(jié)構(gòu)優(yōu)化;衛(wèi)星;工程;應(yīng)用

航天器的尺寸、重量和復(fù)雜性不斷提高,由于火箭的運載能力有限,航天器的結(jié)構(gòu)重量成為一個關(guān)鍵的設(shè)計因素.在衛(wèi)星結(jié)構(gòu)研制初樣階段中,欲實現(xiàn)結(jié)構(gòu)重量最小且同時滿足多個設(shè)計約束,需首先進行反復(fù)的初樣設(shè)計與初樣分析,直到滿足總體的設(shè)計要求.將結(jié)構(gòu)優(yōu)化技術(shù)應(yīng)用于初樣的設(shè)計與分析,可以有效提高設(shè)計效率和水平.在國外的衛(wèi)星結(jié)構(gòu)設(shè)計中已有不少成功案例,如美國 TRW公司在 GRO和 OMV的初樣設(shè)計中應(yīng)用結(jié)構(gòu)優(yōu)化技術(shù),其結(jié)構(gòu)質(zhì)量分別減少 30%和20%[1];在波音等公司,結(jié)構(gòu)優(yōu)化已成為結(jié)構(gòu)設(shè)計中的規(guī)定環(huán)節(jié).我國的衛(wèi)星結(jié)構(gòu)設(shè)計中,優(yōu)化技術(shù)也逐步得到了重視和應(yīng)用,已在某些衛(wèi)星的部件及整星設(shè)計中實施[2-3].由于衛(wèi)星結(jié)構(gòu)中涉及大量內(nèi)部組件和有效載荷的裝配,故應(yīng)用形狀優(yōu)化和拓?fù)鋬?yōu)化的可能性不大.目前,我國的衛(wèi)星結(jié)構(gòu)優(yōu)化設(shè)計中以截面尺寸優(yōu)化為主.然而,截面尺寸優(yōu)化的應(yīng)用過程中,優(yōu)化問題建模的合理性、優(yōu)化系統(tǒng)的可行性、優(yōu)化算法的有效性和計算效率等仍是經(jīng)常遇到的問題.

本文在某新型衛(wèi)星結(jié)構(gòu)研制的初樣階段,建立了衛(wèi)星初始方案的有限元模型.以主承力桁架中各梁的截面尺寸為設(shè)計變量,考慮整星橫向模態(tài)頻率和強度約束,建立了以結(jié)構(gòu)重量最輕為目標(biāo)的結(jié)構(gòu)優(yōu)化模型.應(yīng)用基于二級多點逼近優(yōu)化算法[4]的結(jié)構(gòu)優(yōu)化系統(tǒng) ESSOSII[5]進行優(yōu)化計算,旨在為該衛(wèi)星平臺結(jié)構(gòu)初樣的詳細(xì)設(shè)計提供參考.

1 結(jié)構(gòu)優(yōu)化問題描述

某衛(wèi)星由主結(jié)構(gòu)平臺和有效載荷艙構(gòu)成,設(shè)計總重為 7000 kg.為了提高結(jié)構(gòu)的承載效率,其主結(jié)構(gòu)擬采用梁板復(fù)合結(jié)構(gòu)[6],即由梁和板共用空間、共同承受載荷的結(jié)構(gòu)形式.在初樣階段的初期,應(yīng)根據(jù)工程要求和工藝條件,對初始方案進行分析和優(yōu)化設(shè)計,從而為主結(jié)構(gòu)平臺中承力梁構(gòu)形的拓?fù)湔{(diào)整、截面尺寸和關(guān)鍵結(jié)構(gòu)參數(shù)(如整星質(zhì)心高度等)確定提供設(shè)計依據(jù),同時減輕重量.整星剛度滿足一階橫向和縱向頻率分別不小于 11Hz和 30Hz.同時滿足發(fā)射狀態(tài)下,整星在準(zhǔn)靜態(tài)慣性載荷(見表 1)作用下的強度要求和穩(wěn)定性要求.

表 1 準(zhǔn)靜態(tài)慣性載荷條件

2 衛(wèi)星結(jié)構(gòu)有限元分析模型

根據(jù)方案階段確定的衛(wèi)星主結(jié)構(gòu)的傳力特點和初始設(shè)計構(gòu)形,經(jīng)簡化后建立了結(jié)構(gòu)的有限元分析模型,主要由殼、梁和桿單元構(gòu)成.梁板復(fù)合結(jié)構(gòu)中梁單元與殼單元邊界節(jié)點重合.在結(jié)構(gòu)動力學(xué)特性分析中,考慮結(jié)構(gòu)平臺附屬結(jié)構(gòu)和載荷艙中有效載荷的質(zhì)量分布狀況,將這些質(zhì)量以非結(jié)構(gòu)質(zhì)量的方式施加于各相關(guān)的殼和梁桿單元上,或以集中質(zhì)量點的形式用剛性連接單元RBE2固連于擬安裝的部位.整星有限元模型如圖 1所示,包括 24 669個節(jié)點,26844個單元,圖中去掉了一塊側(cè)板.

圖1 整星有限元模型

根據(jù)星箭的連接狀況,模型的邊界條件定為連接環(huán)底部固支.對其進行模態(tài)分析,并根據(jù)工程需求對其進行修正,使整星模態(tài)頻率分布和振型合理、傳力路線正確,從而可將其作為優(yōu)化設(shè)計的結(jié)構(gòu)分析模型.為了確定合理的整星質(zhì)心高度,通過調(diào)整非結(jié)構(gòu)質(zhì)量的分布,得到質(zhì)心高度分別為1.8m和 1.6m的兩個初始模型,對應(yīng)的整星第 1階模態(tài)頻率分別為 10.41Hz和 11.62Hz.

3 結(jié)構(gòu)優(yōu)化設(shè)計模型

該結(jié)構(gòu)優(yōu)化問題可表述為以下非線性數(shù)學(xué)規(guī)劃問題:

4 二級多點逼近優(yōu)化算法

采用基于二級多點逼近優(yōu)化算法的結(jié)構(gòu)優(yōu)化系統(tǒng) ESSOSII對該問題進行優(yōu)化前處理和尋優(yōu)計算.在每一個優(yōu)化迭代步中,首先進行當(dāng)前結(jié)構(gòu)的模態(tài)分析和 3種過載工況下的靜力分析,然后對約束條件進行評估和臨時刪除,并計算臨界約束和目標(biāo)函數(shù)(結(jié)構(gòu)重量)對設(shè)計變量的一階導(dǎo)數(shù).其中,結(jié)構(gòu)響應(yīng)和相應(yīng)的導(dǎo)數(shù)信息由航天部門廣泛使用的通用有限元程序 Nastran計算得到.

4.1 目標(biāo)函數(shù)和約束條件近似

算法對原優(yōu)化問題實行兩級近似的求解策略.

1)第 1級多點近似函數(shù).算法基于約束函數(shù)在多個點上的函數(shù)值及其導(dǎo)數(shù),構(gòu)造一種多元逼近函數(shù),第 p個結(jié)構(gòu)優(yōu)化周期內(nèi),形成逼近原結(jié)構(gòu)優(yōu)化問題的第 1級序列顯式近似問題為

式中,Δ gj(Xt)為 gj(X)在 Xt的導(dǎo)數(shù);Tt(Xt)為一個對角矩陣;ht(X)為第 t個點的影響因子;指數(shù) rt(t=1,2,…,H)是控制函數(shù)非線性程度的自適應(yīng)參數(shù);H為構(gòu)造近似函數(shù)所用的設(shè)計點數(shù)目上限(取 H=5),當(dāng)已知點的數(shù)目超過 H時,只保留最后得到的 H個點.在已知設(shè)計點上,(X)的數(shù)值及 1階導(dǎo)數(shù)與原函數(shù)相等.在優(yōu)化初始階段,只有一個已知設(shè)計點時,式(3)退化為1階 Taylor展開.隨著優(yōu)化過程的進行,近似問題的有效范圍自動增大,從而具有較高的尋優(yōu)效率和適應(yīng)性.

2)第 2級近似問題.第 1級近似問題(2)中約束函數(shù)的數(shù)目通常遠(yuǎn)小于設(shè)計變量數(shù)目,因此用對偶方法求解具有很高的計算效率.但對偶方法所能求解的規(guī)劃問題要求對偶變量與原變量之間具有顯式的函數(shù)關(guān)系,而問題(2)仍為復(fù)雜非線性,不能利用對偶方法直接求解.因此,問題(2)的目標(biāo)函數(shù)和約束函數(shù)分別在設(shè)計變量空間X和它的倒變量空間進行線性泰勒展開,形成第2級序列近似問題,以逼近問題(2).

4.2 對偶法求解近似問題

第 2級近似問題的解通過求解其對偶問題獲得,詳細(xì)算法見文獻(xiàn)[4].當(dāng)滿足第 2級近似問題的收斂精度后,獲得問題(2)的解.之后,修改設(shè)計變量對應(yīng)的結(jié)構(gòu)參數(shù),進行下一次全結(jié)構(gòu)分析和敏度分析,轉(zhuǎn)入下一個設(shè)計周期,直至達(dá)到設(shè)定的第 1級近似問題的收斂精度后停止迭代.

5 計算結(jié)果

計算過程中,第 1級近似問題和第 2級近似問題的收斂精度分別取 0.005和 0.0005.

5.1 初步試算

首先,對質(zhì)心高度 H為 1.8m的整星模型進行優(yōu)化.由于設(shè)計變量的初值和上下限尚不明確,需在試算的過程中給以初步確定,故梁板復(fù)合結(jié)構(gòu)的內(nèi)埋梁暫取薄壁圓管截面,設(shè)內(nèi)徑 R=0.04m不變,以壁厚 t作為設(shè)計變量.主隔板內(nèi)埋梁左右對稱,分為 5組(如圖 2a所示),分別對應(yīng)第 1～5號設(shè)計變量;次隔板內(nèi)埋梁壁厚對應(yīng)第 6號設(shè)計變量.第 7～11號設(shè)計變量分別為主支架圓管梁壁厚、對接環(huán)法蘭寬度、儲箱拉桿寬度和上下框架矩形截面梁壁厚.設(shè)計變量的初值取經(jīng)驗值,對上下限暫設(shè)置了較大的變化范圍,以考察其變化趨勢.整星基頻約束 ω1分別為 11.0Hz和10.5Hz時的優(yōu)化結(jié)果見表 2.可見,整星模態(tài)基頻是臨界約束;質(zhì)心高度為 1.8m時,欲使基頻達(dá)到11.0Hz將導(dǎo)致重量大幅增加.

圖2 主隔板內(nèi)埋梁

表 2 初步試算優(yōu)化結(jié)果 mm

接著,對質(zhì)心高度 1.6m的整星模型進行優(yōu)化.設(shè)計變量的初值和上下限不變,整星基頻約束分別取 11.8Hz,11.6Hz和 11.4Hz,計算結(jié)果見表 1.可見,質(zhì)心高度 1.6m時,隨著模態(tài)基頻下限的變化,各結(jié)構(gòu)參數(shù)的變化趨勢較為明顯.基頻下限取 11.4Hz時,6組內(nèi)埋梁的壁厚中,次隔板內(nèi)埋梁壁厚最大,為 2.91mm,主隔板內(nèi)埋梁的壁厚均有較大幅度的下降;上下框架的內(nèi)梁壁厚都接近下限值;對接環(huán)上法蘭的寬度增至 135mm.當(dāng)基頻下限取 11.8Hz時,結(jié)構(gòu)質(zhì)量上升.

5.2 改進設(shè)計

經(jīng)過以上試算,該主承力桁架結(jié)構(gòu)設(shè)計方向逐漸明確:①整星質(zhì)心高度為 1.8m時,模態(tài)基頻很難達(dá)到 11.0Hz,整星質(zhì)心高度調(diào)整至 1.6m較為合理;②主支架圓管梁壁厚對結(jié)構(gòu)剛度的影響大;③內(nèi)埋梁中,次隔板的內(nèi)埋梁對結(jié)構(gòu)剛度影響最大,主隔板內(nèi)埋梁的尺寸應(yīng)有所不同;④應(yīng)增大對接環(huán)上法蘭的寬度(設(shè)計變量 8);⑤上下框架內(nèi)梁的壁厚(設(shè)計變量 10,11)可減小.

由此,提出以下改進設(shè)計:①平臺底板由初始的中心開口加下框架的設(shè)計改為封閉結(jié)構(gòu),從而去掉底板處的主隔板內(nèi)埋梁 4,5和下框架;②設(shè)計變量 8由 100mm增至 130mm,不再參與優(yōu)化.主隔板和次隔板的內(nèi)埋梁由圓管梁改為接近實際結(jié)構(gòu)的薄壁矩形梁,取其寬和高分別為 70mm和28mm,對應(yīng)的壁厚分別為 t1=2 mm和 t2=0.7mm,其中 t1作為設(shè)計變量,分組見圖 2b.

5.3 再優(yōu)化

以 11.0Hz為基頻約束,再次進行優(yōu)化計算.此時,為了獲得合理的結(jié)構(gòu)參數(shù),設(shè)計變量的尋優(yōu)區(qū)間縮小,其上下限和優(yōu)化結(jié)果見表 3.目標(biāo)函數(shù)(結(jié)構(gòu)質(zhì)量)的優(yōu)化迭代過程見圖 3,其中 W為可優(yōu)化部分的結(jié)構(gòu)質(zhì)量,初值為 372.8kg.

表 3 改進方案優(yōu)化結(jié)果 mm

可見,改進方案優(yōu)化后基頻由 10.79Hz提高到 11Hz,重量則減輕了 3.2kg,顯示了優(yōu)化的效果.優(yōu)化中,結(jié)構(gòu)模態(tài)基頻仍是臨界約束,第 2種過載工況下的應(yīng)力峰值最大,為 149MPa,遠(yuǎn)小于材料的強度極限.對優(yōu)化后的結(jié)構(gòu)進行校核,縱向模態(tài)頻率為 34.27Hz(＞30Hz),第 2種過載工況下的最小臨界穩(wěn)定性系數(shù)由初始的 2.02增至2.21(＞1.5).可見,優(yōu)化后結(jié)構(gòu)的剛度、強度和穩(wěn)定性均滿足設(shè)計要求,可作為衛(wèi)星初樣詳細(xì)設(shè)計的參考.

圖3 目標(biāo)函數(shù)(結(jié)構(gòu)質(zhì)量)迭代曲線

6 結(jié)束語

本文基于二級多點逼近優(yōu)化算法,對某新型衛(wèi)星主承力桁架結(jié)構(gòu)初樣方案進行了優(yōu)化設(shè)計.經(jīng)過初步試算、設(shè)計改進和再優(yōu)化 3個階段,逐步明確了設(shè)計的方向和各結(jié)構(gòu)參數(shù)的取值范圍,得到了合理的可行方案,為該衛(wèi)星平臺結(jié)構(gòu)初樣的詳細(xì)設(shè)計提供了參考;同時也表明了該優(yōu)化算法的可行性和有效性.

References)

[1]Woo T H.Structural optimization of large spacecraft[R].AIAA 1992-1227,1992

[2]陳珅艷,黃海,王利民.人機交互的航天器結(jié)構(gòu)優(yōu)化設(shè)計[J].計算力學(xué)學(xué)報,2004,21(4):510-512 Chen Shenyan,Huang Hai,Wang Liming.Structural optimization of spacecraft with man-machine interaction[J].Chinese Journal of Computational Mechanics,2004,21(4):510-512(in Chinese)

[3]顧元憲,亢戰(zhàn),趙國忠,等,衛(wèi)星承力筒復(fù)合材料結(jié)構(gòu)的優(yōu)化設(shè)計[J].宇航學(xué)報,2003,24(1):88-91 Gu Yuanxian,Kang Zhan,Zhao Guozhong,et al.Optimal design of composite structure of bearing cylinder[J].Journal of Astronautics,2003,24(1):88-91(in Chinese)

[4]Huang H,Xia RW.Two-level multipoint constraint approximation concept for structural optim ization[J].Structural Optim ization,1995,9(1):38-45

[5]陳珅艷,袁家軍,黃海.基于 MSC.PATRAN/NASTRAN的結(jié)構(gòu)優(yōu)化系統(tǒng)[J].宇航學(xué)報,2005,26(4):514-518 Chen Shenyan,Yuan Jiajun,Huang Hai.Structural optimization system based on MSC.PATRAN/NASTRAN[J].Journal of Astronautics,2005,26(4):514-518(in Chinese)

[6]盛聰,柴洪友,高峰,等.一種梁板復(fù)合結(jié)構(gòu)的設(shè)計[J].航天器工程,2009,18(2):62-67 Sheng Cong,Chai Hongyou,Gao Feng,et al.Design of A kind of combined structure of beam and plate[J].Spacecraft Engineering,2009,18(2):62-67(in Chinese)

(編 輯 :張 嶸)

Structural optimization o f spacecraft based on two-level multi-point approximate method

Chen Shenyan Huang Hai

(School of Astronautics,Beijing University of Aeronautics and Astronautics,Beijing 100191,China)

At the primal design stage of a new satellite,finite element model of the initial structure was established.Taking thicknesses of pipe beams and rectangle thin-wall beams as design variables,and considering constraints of natural frequency and strength,a structural optimization model was developed to approach minimum mass.Two-level multi-point approximation method was applied in structural optimization.In each calculation cycle,the original optimization problem was transformed to a first level multi-point approximate problem,and it was also approached by a second approximate problem which can be quickly solved by dual method.Three stages were implemented in the whole optimization process:initial optimum searching tests,design revision and re-optimization.Among the first stages,five runs were set with different constraints,design variables bounds and structural centroid heights,through which the design direction and interval of various structural parameters were obtained.Then the original structure was revised according to the results of the first stage.A new optimization model with fewer design variables and more reasonable variables bounds is established in the third stage and a feasible design is achieved,which provide references to detail design of the satellite structure.It indicates that the app lied design procedure and optimization method are effective in engineering structural optimization problem.

structural optimization;satellites;engineering;applications

V 414.19

A

1001-5965(2011)02-0237-04

2009-12-16

陳珅艷(1974-),女,浙江東陽人,講師,chenshenyan@buaa.edu.cn.