梅 勇，馮韶偉

大推力捆綁運載火箭傳力路徑優(yōu)化設計

梅 勇1，馮韶偉2

（1. 中央軍委后勤保障部工程兵科研三所，洛陽，471023；2. 北京宇航系統(tǒng)工程研究所，北京，100076）

以新一代大推力捆綁火箭捆綁裝置設計為研究背景，開展助推器傳力路徑優(yōu)化設計研究。在方案初步優(yōu)選的基礎上，基于iSΙGHT建立大推力捆綁火箭傳力路徑優(yōu)化平臺，針對現(xiàn)有優(yōu)化方法的不足，選用多島遺傳（Multi-Ιsland Genetic Algorithm，MΙGA）+二次規(guī)劃（Sequential Quadratic Programming，SQP）組合優(yōu)化算法進行傳力路徑優(yōu)化設計，從而大幅優(yōu)化了捆綁聯(lián)接結構載荷。該方法可為大推力捆綁運載火箭方案設計提供參考和依據(jù)。

大推力捆綁運載火箭；傳力路徑；遺傳算法；組合優(yōu)化策略

0 引 言

一個國家進入空間的能力很大程度上決定了其空間活動能力以及改建應用開發(fā)水平，而運載火箭的規(guī)模和水平支撐了一個國家進入空間的能力[1]。在載人登月和大規(guī)模深空探測的背景下，大推力的超大型或重型運載火箭逐漸成為研究的熱點和重點[2]。

在運載火箭芯級周圍捆綁助推器能夠有效提高進入空間的能力。助推器和芯級之間采用捆綁聯(lián)接裝置保持連接，捆綁聯(lián)接裝置常成組安裝，構成助推器推力的傳力路徑。傳力路徑設計是運載火箭方案制定、捆綁聯(lián)接裝置選型、結構設計和強度校核的依據(jù)[3]。國內(nèi)外目前常見捆綁火箭芯級與助推器之間構成靜定體系，將助推器傳遞給芯級的軸向載荷與橫向載荷分別由不同的裝置承擔，分別為主承力結構和輔助承力結構，其中前者用于傳遞來自助推器的軸向推力，為主要的承載裝置，限制助推器的平動自由度；后者傳遞來自助推器的剪力、扭轉(zhuǎn)和徑向力，起限制助推器轉(zhuǎn)動自由度的作用[4]。常見主承力結構方案和輔助承力方案分別如圖1、圖2所示[5]。捆綁聯(lián)接裝置中具有分離機構，在接到分離指令后，迅速解鎖，使助推器與芯級完成分離[6]。

圖1 常見的主承力方案

圖2 常見的輔助承力方案

中國新一代大推力運載火箭的助推器推力是現(xiàn)役火箭最大助推器的2倍以上，若采用傳統(tǒng)的傳力路徑和捆綁聯(lián)接裝置，捆綁聯(lián)接裝置和分離機構的結構強度需求將成倍增加，考慮到結構設計難度及目前的火工品設計水平，現(xiàn)有的傳力路徑及捆綁聯(lián)接裝置無法滿足工作要求，因此需要設計新型傳力路徑，從而保證新一代大推力運載火箭助推器的可靠聯(lián)接。

本文首先初步設計兩類新型的超靜定捆綁聯(lián)接形式，通過多種外載荷工況作用下的結構承載分析，初步優(yōu)選新型火箭傳力路徑，基于多學科優(yōu)化設計軟件iSΙGHT，通過集成捆綁火箭有限元分析模型，建立傳力路徑優(yōu)化平臺，采用多島遺傳（Multi-Ιsland Genetic Algorithm，MΙGA）+二次序列規(guī)劃（Sequential Quadratic Programming，SQP）組合優(yōu)化算法進行優(yōu)化設計，最終得到優(yōu)化的傳力路徑。

1 大推力捆綁火箭捆綁載荷分析

研究表明，主承力裝置的受力狀態(tài)不隨捆綁位置的變化而改變，且采用前捆綁主承力形式、助推器支撐芯級的方案能夠避免捆綁聯(lián)接裝置在拉載下破壞[7]。在中國新一代大推力捆綁火箭論證過程中，助推器的規(guī)模和推力是中國現(xiàn)役最大推力火箭助推器的2倍以上。由于現(xiàn)役火箭的捆綁聯(lián)接裝置，尤其是分離機構已接近火工品設計的最高水平，因此在捆綁載荷成倍增加的情況下，現(xiàn)役的助推器及芯級的聯(lián)接部位結構強度難以保證，因此有必要對大推力火箭的傳力路徑進行重新設計。新一代大推力運載火箭構型見圖3。

圖3 新一代大推力運載火箭構型

1.1 三支點超靜定捆綁聯(lián)接結構載荷分析

在原靜定捆綁構型前、后兩捆綁點之間增加一套中捆綁聯(lián)接裝置，從而構成三支點超靜定捆綁方案，如圖3c所示，其中前捆綁點位于芯級和助推前箱的前短殼之間，中捆綁點位于芯一級和助推器的箱間段之間，二者均為輔助承力裝置，采用圖2的方案A；后捆綁點為主承力裝置，位于芯一級和助推器的后過渡段之間，采用圖1的方案1。通過有限元分析獲得多種外載荷工況作用下，火箭捆綁聯(lián)接結構各部件最大受力情況，如表1所示。

表1 三支點超靜定捆綁構型聯(lián)接結構承載

表1的計算結果表明，采用單主承力的超靜定捆綁方案，主承力裝置的受力狀態(tài)與兩點靜定支撐的相似，即其受力仍遠大于現(xiàn)有結構所能承受的最大值。工程中，超靜定捆綁連接方式主要用于改善助推器的局部模態(tài)。

1.2 兩點超靜定捆綁聯(lián)接結構載荷分析

為了改善大推力運載火箭捆綁聯(lián)接結構的受力特點，可采用兩點均為主承力的超靜定捆綁聯(lián)接方式。在捆綁位置固定的假設下，捆綁結構分別采用圖1所示方案1和方案2的組合形式，共4種構型，分別為：

a）構型Ι：前捆綁采用具有捆綁軸安裝角的主承力結構，后捆綁采用主承力+兩連桿的水平承力結構；

b）構型ΙΙ：前捆綁同構型Ι，后捆綁采用具有捆綁軸安裝角的主承力+兩連桿結構；

c）構型ΙΙΙ：前捆綁同構型Ι的后捆綁結構形式，后捆綁同構型Ι的前捆綁結構形式；

d）構型ΙV：前捆綁同構型ΙΙ的后捆綁結構形式，后捆綁同構型ΙΙ的前捆綁結構形式。

4種構型的有限元模型如圖4所示。

圖4 大推力運載火箭新型兩點超靜定聯(lián)接構型

在相同的外載情況下，火箭捆綁聯(lián)接結構各部件最大受力如表2所示。

表2 兩點超靜定捆綁構型聯(lián)接結構承載

以主捆綁軸向力為主要判斷依據(jù)，綜合考慮承力結構各部分最大受力，同時盡量保持結構處于受壓狀態(tài)[8]，得到捆綁構型從優(yōu)到劣的排序為：構型Ι，構型ΙΙΙ，構型ΙΙ，構型ΙV。因此，本文選定構型Ι進行傳力路徑的進一步優(yōu)化。

2 大推力捆綁火箭傳力路徑優(yōu)化設計

2.1 優(yōu)化函數(shù)

以構型Ι為研究對象，研究當整箭在多種外載荷工況作用下的最優(yōu)化傳力路徑問題。以兩主承力軸軸向力、剪力、主承力鉸受力和輔助承力桿軸力的最大值為優(yōu)化目標，研究前捆綁點坐標、后捆綁點坐標以及主捆綁軸安裝角對于優(yōu)化目標的影響。同時傳力路徑的設計需要考慮姿控系統(tǒng)設計的要求，即整箭的橫向彎曲基頻要小于助推器局部的一階振動頻率。設計變量的選擇范圍為：

a）前捆綁點坐標范圍為61.994～62.994 m，捆綁聯(lián)接段長1 m；

b）后捆綁點坐標范圍為92.722～93.369 m，捆綁聯(lián)接段長0.647 m；

c）前捆綁軸安裝角的范圍：11～45°。

捆綁聯(lián)接傳力路徑的多目標優(yōu)化問題數(shù)學描述為

式中 FAxial為主承力軸最大軸向力；FShear為主承力軸最大剪力；FRBE2為主承力點最大合力；FCrod為輔助承力桿最大軸向力；xfront為前捆綁點軸向位置，xbehind為后捆綁點軸向位置；α為前捆綁軸安裝角；freq為火箭的自振頻率。

2.2 優(yōu)化算法

SQP在解決中小規(guī)模的非線性優(yōu)化問題中優(yōu)勢明顯，其收斂速度快、效率高，是國際上公認最有效的優(yōu)化算法之一[9]。SQP在每個迭代點均構造一個二次規(guī)劃子問題，并將該問題的解作為迭代搜索方向進行一維搜索，從而逼近約束優(yōu)化的解[10]。但SQP受初值影響明顯，容易陷入局部最小[11]。

MΙGA是在傳統(tǒng)的遺傳算法基礎上發(fā)展而來的，該算法將待優(yōu)化群體分為若干子群體（稱為“島嶼”），在每個島嶼上獨立地進行遺傳算法操作，通過定期對各島嶼內(nèi)先進個體的遷移，可以在保持群體多樣性的同時，進行高效的全局尋優(yōu)[12]。

為了綜合2種優(yōu)化算法的優(yōu)勢，利用Patran二次開發(fā)語言PCL編寫的參數(shù)化建模文件，通過自動有限元模型的前處理操作、后臺求解及結果后處理功能實現(xiàn)目標函數(shù)的求解。基于iSΙGHT優(yōu)化分析平臺，采用MΙGA+SQP組合優(yōu)化進行優(yōu)化[13]。首先采用多島遺傳算法找到多個局部近似的最優(yōu)點，然后以這些局部最優(yōu)點為初始點，采用SQP二次規(guī)劃法優(yōu)化，從而得到全局最優(yōu)點[14]。組合優(yōu)化流程如圖5所示。

圖5 組合優(yōu)化流程

優(yōu)化仿真系統(tǒng)數(shù)據(jù)流如圖6所示。

圖6 優(yōu)化仿真系統(tǒng)數(shù)據(jù)流

2.3 優(yōu)化結果

基于iSΙGHT平臺的新一代大推力運載火箭傳力路徑受力的優(yōu)化結果見圖7至圖9，輔助承力桿結構受力的優(yōu)化結果如圖10所示，優(yōu)化結果最終趨于收斂。優(yōu)化前后的受力結果如表3所示。

圖7 主承力軸軸向力優(yōu)化歷程

圖8 主承力軸剪力優(yōu)化歷程

圖9 主承力合力優(yōu)化歷程

圖10 輔助承力桿軸向力優(yōu)化歷程

表3 結構受力優(yōu)化結果

結果表明，經(jīng)過傳力路徑優(yōu)化，大推力運載火箭的捆綁聯(lián)接裝置的受力得到了優(yōu)化，其中前捆綁主承力裝置的軸向力降低了33.4%。

優(yōu)化前后全箭的頻率變化情況見表4，結果表明優(yōu)化后，全箭各階模態(tài)頻率有所提升，其中對于一階橫向基頻提升最為明顯，為1.7%。優(yōu)化后，助推器局部模態(tài)小于全箭橫向基頻，有利于提升飛行的可靠性以及降低控制系統(tǒng)的設計難度。荷降低了33.4%，同時優(yōu)化后的全箭動特性第1階橫向基頻有所提高，大于助推器第一階局部模態(tài)，有利于提升飛行的可靠性以及降低控制系統(tǒng)的設計難度。

表4 整箭自振頻率計算結果

3 結 論

本文以大推力運載火箭捆綁方案論證為背景，針對現(xiàn)有傳力路徑和捆綁結構無法滿足使用要求的情況，設計并初步優(yōu)選一類超靜定捆綁的傳力形式。通過在iSΙGHT優(yōu)化平臺上集成MSC.Patran/Nastran，建立了新一代大推力運載火箭有限元參數(shù)化模型。利用多島遺傳算法和序列二次規(guī)劃法組成的優(yōu)化組合策略對火箭傳力路徑進行多目標優(yōu)化研究。

結果表明，經(jīng)過優(yōu)化后，前捆綁主承力裝置的載

[1] 何巍, 劉偉, 龍樂豪. 下一代大推力運載火箭及其應用探討[J]. 導彈與航天運載技術, 2011(1): 1-5.

[2] 龍樂豪, 王小軍, 容易. 我國一次性運載火箭的發(fā)展展望[J]. 中國科學E輯: 技術科學, 2009, 39(3): 460-463.

[3] 李東, 程堂明. 中國新一代運載火箭發(fā)展展望[J]. 中國工程科學, 2006, 8(11): 33-38.

[4] 馮韶偉, 劉竹生, 馬忠輝, 姚世東. 捆綁火箭助推器與芯級間捆綁聯(lián)接技術應用進展[J]. 導彈與航天運載技術, 2012(6): 20-23.

[5] 馮韶偉. 超靜定捆綁火箭力學特性研究[D]. 北京: 中國航天科技集團公司, 2013.

[6] 馮韶偉. 基于拉壓不同剛度的運載火箭捆綁聯(lián)接裝置力學特性研究[J].導彈與航天運載技術, 2013(4): 9-13.

[7] 魯宇, 龍樂豪, 等. 世界航天運載器大全（第二版）[M]. 北京: 中國宇航出版社, 2007.

[8] 汪軼俊. 運載火箭固體捆綁技術研究[D]. 長沙: 國防科學技術大學, 2007.

[9] Boggs P T , Tolle J W. Sequential quadratic programming for large-scale nonlinear optimization[J]. Journal of Computational and Applied Mathematics, 2000, 124(1-2): 123-137.

[10] Chen H , Ooka R , Kato S. Study on optimum design method for pleasant outdoor thermal environment using genetic algorithms and coupled simulation of convection, radiation and conduction[J]. Building and Environment, 2008, 43(1): 18-31.

[11] 韓邦成, 袁倩. 大型磁懸浮CMG轉(zhuǎn)子的組合優(yōu)化策略[J]. 宇航學報, 2012, 33(2): 275-280.

[12] 任利. 基于iSΙGHT的多學科設計優(yōu)化平臺的研究與實現(xiàn)[D]. 青島:山東科技大學, 2006.

[13] 薛定宇, 陳陽泉. 高等應用數(shù)學問題的MATLAB求解[M]. 北京: 清華大學出版社, 2004.

[14] 賴炎連, 賀國平. 最優(yōu)化方法[M]. 北京: 清華大學出版社, 2008.

An Optimization Design for the Force Transmission Path on the High-thrust Strap-on Launch Vehicle

Mei Yong1, Feng Shao-wei2
(1. The Third Engineer Scientific Research Ιnstitute of Logisric Support Department , Luoyang, 471023; 2. Beijing Ιnstitute of Astronautical Systems Engineering, Beijing, 100076)

On the background of the engineering problem that the strap-on linkage equipment of the launch vehicles on active duty can’t satisfy the load design requirements of the new generation high-thrust strap-on launch vehicle and thus the research and development of new linkage equipment technology is needed. This essay studies the force transmission path optimization and designs two hyper-static categories of force transmission path schemes. Through the finite element analysis of the initial configuration of the high-thrust launch vehicle, a optimal scheme of the force transmission path is selected. Design the combinatorial optimization strategy of the force transmission path of high-thrust strap-on launch vehicle based on the iSΙGHT platform: Firstly, use multi-island genetic algorithm to explore the region where the globally optimal solution lies and then approach the exact optimal solution by sequential quadratic programming. The computational results of the structure force and natural frequency suggest the new scheme of force transmission path can substantially optimize the load of strap-on linkage equipment and also satisfy the requirements of the attitude control design.

High-thrust strap-on launch vehicle; Force transmission path; Genetic algorithm; Combination optimization strategy

V421

A

1004-7182(2017)04-0006-05 DOΙ：10.7654/j.issn.1004-7182.20170402

2015-12-01；

2017-06-09

梅 勇（1990-），男，助理研究員，主要研究方向為結構分析與優(yōu)化設計