楊風(fēng)波，馬大為，任 杰，樂(lè)貴高，蔡德詠

(1. 南京理工大學(xué) 機(jī)械工程學(xué)院，南京 210094；2. 72465部隊(duì)，濟(jì)南 250022)

?

新型同心筒自力發(fā)射熱環(huán)境優(yōu)化設(shè)計(jì)①

楊風(fēng)波1，馬大為1，任 杰1，樂(lè)貴高1，蔡德詠2

(1. 南京理工大學(xué) 機(jī)械工程學(xué)院，南京 210094；2. 72465部隊(duì)，濟(jì)南 250022)

針對(duì)新型路基同心筒自力發(fā)射熱環(huán)境評(píng)估與優(yōu)化設(shè)計(jì)問(wèn)題，依托彈性變形和域動(dòng)分層結(jié)合的動(dòng)網(wǎng)格技術(shù)，求解了二維軸對(duì)稱(chēng)N-S方程，分析了“中段導(dǎo)流”同心筒動(dòng)態(tài)熱環(huán)境特性，確定了熱環(huán)境評(píng)價(jià)指標(biāo)；通過(guò)建立以?xún)?yōu)化拉丁超立方實(shí)驗(yàn)設(shè)計(jì)和四階響應(yīng)面為理論基礎(chǔ)的近似數(shù)學(xué)模型，解決了CFD自動(dòng)建模困難、直接尋優(yōu)計(jì)算量大的難點(diǎn)；利用多島遺傳和梯度優(yōu)化算法搭建組合優(yōu)化策略平臺(tái)，克服了流場(chǎng)在不同熱結(jié)構(gòu)條件下的強(qiáng)非線性問(wèn)題，并構(gòu)建了支持近似數(shù)學(xué)模型的熱環(huán)境優(yōu)化構(gòu)架。對(duì)比數(shù)值結(jié)果表明，倒吸進(jìn)入內(nèi)筒的低溫氣體有力改善了同心筒熱環(huán)境；建立的近似數(shù)學(xué)模型精度較高，滿(mǎn)足工程需求；優(yōu)化后，熱環(huán)境特性發(fā)生良性變化，導(dǎo)彈總體熱環(huán)境得到顯著改善。

同心筒發(fā)射裝置；燃?xì)馍淞鳎粺岘h(huán)境特性；近似數(shù)學(xué)模型；優(yōu)化設(shè)計(jì)

0 引言

同心筒以其結(jié)構(gòu)簡(jiǎn)單、快速反應(yīng)、全方位覆蓋、通用性好等優(yōu)點(diǎn)，滿(mǎn)足了現(xiàn)代戰(zhàn)爭(zhēng)對(duì)武器裝備綜合作戰(zhàn)能力的要求，備受各國(guó)海軍青睞[1-2]。近年來(lái)，為完善三位一體戰(zhàn)略威懾力量品種，路基車(chē)載同心筒自力發(fā)射成為一種新的發(fā)展方向。針對(duì)艦載和潛載同心筒，國(guó)內(nèi)學(xué)者跟進(jìn)國(guó)外研究，做了很多工作[3-6]，致力于從拓?fù)浣Y(jié)構(gòu)和氣液多相流兩方面來(lái)改善熱環(huán)境；同時(shí)，國(guó)內(nèi)開(kāi)始通過(guò)窮舉法[7]分析不同熱結(jié)構(gòu)條件下導(dǎo)彈的熱環(huán)境，而針對(duì)同心筒熱環(huán)境的優(yōu)化設(shè)計(jì)鮮見(jiàn)報(bào)道。同心筒的熱環(huán)境直接關(guān)系到導(dǎo)彈及發(fā)射裝置的熱安全，對(duì)路基新型同心筒進(jìn)行熱環(huán)境優(yōu)化設(shè)計(jì)具有重要的工程應(yīng)用價(jià)值。

目前，多學(xué)科優(yōu)化設(shè)計(jì)(Multidisciplinary design optimization, MDO)已經(jīng)在諸多領(lǐng)域得到應(yīng)用，其中應(yīng)用到工程實(shí)踐中的優(yōu)化算法也有多種。成沉等[8]結(jié)合響應(yīng)面近似模型，采用序列二次規(guī)劃(NLPQL)算法，優(yōu)化了喉栓式噴管的氣動(dòng)性能；陳偉等[9]利用內(nèi)點(diǎn)法和NLPQL算法，對(duì)噴管型面進(jìn)行了優(yōu)化設(shè)計(jì)；粟華等[10]采用響應(yīng)面與梯度算法組合的搜索方式，實(shí)現(xiàn)了再入飛行器的多學(xué)科優(yōu)化設(shè)計(jì)；邱東海等[11]利用多島遺傳算法和二次規(guī)劃法組合優(yōu)化策略，對(duì)無(wú)人機(jī)地面隔振系統(tǒng)的隔振性能進(jìn)行了優(yōu)化設(shè)計(jì)。根據(jù)公開(kāi)的報(bào)道來(lái)看，搭配全局優(yōu)化算法和局部?jī)?yōu)化算法的組合優(yōu)化策略能充分發(fā)揮2種算法的特點(diǎn)，形成互補(bǔ)優(yōu)勢(shì)，成為研究熱點(diǎn)。

基于此，本文采用組合優(yōu)化策略，對(duì)同心筒動(dòng)態(tài)熱環(huán)境進(jìn)行優(yōu)化設(shè)計(jì)。首先，結(jié)合域動(dòng)分層動(dòng)網(wǎng)格技術(shù)，通過(guò)求解N-S方程研究新型同心筒的熱環(huán)境特性，確定評(píng)價(jià)導(dǎo)彈熱環(huán)境品質(zhì)的技術(shù)指標(biāo)；其次，結(jié)合試驗(yàn)設(shè)計(jì)采樣和和四階響應(yīng)面方法，建立精度可靠的近似數(shù)學(xué)模型；最后，搭建基于多島遺傳和二次規(guī)劃算法組合策略的優(yōu)化平臺(tái)，對(duì)導(dǎo)彈熱環(huán)境進(jìn)行優(yōu)化設(shè)計(jì)。

1 流體模型與數(shù)值計(jì)算方法

1.1 流體基本控制方程與湍流模型

本文進(jìn)行動(dòng)態(tài)數(shù)值計(jì)算，屬于流固耦合范疇，求解過(guò)程伴隨著控制體的運(yùn)動(dòng)、消失與生成，對(duì)流體運(yùn)動(dòng)描述宜采用任意拉格朗日歐拉方法(ALE)。ALE形式軸對(duì)稱(chēng)Navier-Stokes方程組的守恒形式可表述為

(1)

式中 各項(xiàng)的物理意義和文獻(xiàn)[12]一致。

1.2 計(jì)算方法與邊界條件

邊界條件、湍流模型和文獻(xiàn)[12]一致。圖1給出了噴管的入口條件，包括總壓和總溫度變化曲線。

2 模型描述與邊界條件

本文以某新型同心筒為基準(zhǔn)模型，針對(duì)主要熱結(jié)構(gòu)參數(shù)，對(duì)導(dǎo)彈熱環(huán)境進(jìn)行優(yōu)化設(shè)計(jì)。

圖2分別給出了傳統(tǒng)同心筒優(yōu)化結(jié)構(gòu)和新型同心筒結(jié)構(gòu)方案。新型方案中，內(nèi)、筒上部結(jié)構(gòu)做成一體，中部有導(dǎo)流板，該方案具有完全軸對(duì)稱(chēng)性質(zhì)，取二維軸對(duì)稱(chēng)流場(chǎng)為研究對(duì)象。

(a) 總壓 (b) 總溫

(a) 傳統(tǒng)同心筒

(b) 新型同心筒

3 熱環(huán)境特性與熱環(huán)境評(píng)價(jià)指標(biāo)

3.1 熱環(huán)境特性

文獻(xiàn)[12]已經(jīng)對(duì)同心筒的流場(chǎng)機(jī)理與熱環(huán)境特性進(jìn)行了詳細(xì)分析，得到了有益結(jié)論。本文在文獻(xiàn)[12]研究工作的基礎(chǔ)上，考慮了新型導(dǎo)流結(jié)構(gòu)和發(fā)射車(chē)的高度匹配，其基準(zhǔn)模型如圖2所示。

為揭示同心筒這一特殊發(fā)射裝置的流場(chǎng)特征與熱環(huán)境機(jī)理，且考慮到傳統(tǒng)優(yōu)化方案流動(dòng)現(xiàn)象更明顯，圖3給出了傳統(tǒng)優(yōu)化方案發(fā)射前期和后期典型的速度矢量圖。從圖3(a)可看出，在前期的0.058 s時(shí)刻，內(nèi)、外筒燃?xì)饩霈F(xiàn)了向上流動(dòng)的特征?？蓮囊韵路矫娼沂荆喊l(fā)動(dòng)機(jī)點(diǎn)火，燃?xì)饬饕环矫嬖趯?dǎo)流器排導(dǎo)作用下，順著外筒流出；另一方面，大量燃?xì)庠谕驳撞荒苎杆倥懦?，發(fā)生了“壅塞現(xiàn)象”，導(dǎo)致筒底壓力大于環(huán)境壓力，在壓差作用下，部分燃?xì)饨?jīng)過(guò)內(nèi)筒向上排出，導(dǎo)致觀測(cè)1面的溫度急劇上升，如圖4(a)所示；在導(dǎo)彈啟動(dòng)的后期，燃?xì)饬魑惭孚呌诜€(wěn)定，筒底導(dǎo)流板遮擋反濺流，導(dǎo)流器實(shí)現(xiàn)了高溫氣流的有效排導(dǎo)，且氣流速度很大，最終使得筒底壓力小于外界壓力，在負(fù)壓的作用下，內(nèi)筒出現(xiàn)“倒吸現(xiàn)象”。內(nèi)筒初始低溫氣體使得觀測(cè)1面溫度急劇下降；隨后，部分氣流從筒口經(jīng)過(guò)內(nèi)筒流向彈底(如圖3(b)、 (c)速度矢量圖所示)，但不同結(jié)構(gòu)同心筒的筒口氣流熱環(huán)境不一樣，使得觀測(cè)1面的熱環(huán)境不同，由于本文方案中部導(dǎo)流板結(jié)構(gòu)和發(fā)射車(chē)的高度相匹配，后期筒口熱環(huán)境最為優(yōu)良，倒吸入內(nèi)筒的均為低溫氣體，如圖4(b)所示。

(a) 傳統(tǒng)方案0.058 s時(shí)刻 (b) 傳統(tǒng)方案0.126 s時(shí)刻

(c) 傳統(tǒng)方案0.126 s時(shí)刻(局部放大)

(a) 0.034 s (b) 0.5 s

從圖5可看出，觀測(cè)1面的溫度時(shí)程曲線和流場(chǎng)特征及熱環(huán)境機(jī)理吻合良好，但不同熱結(jié)構(gòu)其變化細(xì)節(jié)確有不同。傳統(tǒng)方案和導(dǎo)彈底部均出現(xiàn)了溫度先上升、后下降、然后再上升的燃?xì)鉄g過(guò)程，而本文方案中導(dǎo)彈底部在后期熱環(huán)境友好，溫度沒(méi)有反彈。

3.2 熱環(huán)境評(píng)價(jià)指標(biāo)

本文根據(jù)路基新型同心筒的技術(shù)要求，對(duì)圖6中的3個(gè)參數(shù)進(jìn)行優(yōu)化，使得導(dǎo)彈的熱環(huán)境達(dá)到最優(yōu)。而確定評(píng)價(jià)導(dǎo)彈熱環(huán)境的技術(shù)指標(biāo)是進(jìn)行熱環(huán)境優(yōu)化設(shè)計(jì)的基礎(chǔ)。

圖5 2種方案觀測(cè)1溫度時(shí)程曲線Fig.5 Temperature curves of the observation plane 1 in the two scheme

圖6 新型同心筒優(yōu)化參數(shù)Fig.6 Optimization parameters of a new CCL scheme

從3.1節(jié)的分析可看出，在0.2 s之前，倒吸進(jìn)入內(nèi)筒的冷氣體使得觀測(cè)1面的溫度已經(jīng)基本穩(wěn)定在低溫狀態(tài)。為摸清導(dǎo)彈熱環(huán)境特性，本文隨機(jī)取3組L1、L2、Ф組合建模，并進(jìn)行了動(dòng)態(tài)計(jì)算。圖7給出了3種方案觀測(cè)1面的溫度時(shí)程曲線，溫度變化規(guī)律和3.1的分析基本一致。

圖7 3種隨機(jī)方案觀測(cè)1面溫度時(shí)程曲線Fig.7 Temperature curves of the observation plane 1 for three random schemes

4 導(dǎo)彈熱環(huán)境組合優(yōu)化設(shè)計(jì)

新型同心筒自力發(fā)射動(dòng)態(tài)熱環(huán)境優(yōu)化涉及到幾何造型、網(wǎng)格劃分、CFD動(dòng)態(tài)計(jì)算、優(yōu)化目標(biāo)的確定、試驗(yàn)設(shè)計(jì)方案的選擇、近似模型的建立，優(yōu)化策略的確定等，是一個(gè)典型的多學(xué)科協(xié)同問(wèn)題，工作平臺(tái)的搭建是基石，原理流程如圖8所示。

表1 3種隨機(jī)方案觀測(cè)1面相關(guān)參量Table1 Related parameters for the observation plane 1 of three random schemes

圖8 優(yōu)化平臺(tái)原理流程圖Fig.8 Principle flowchart of optimization platform

4.1 試驗(yàn)設(shè)計(jì)

試驗(yàn)設(shè)計(jì)(design of experiment, DOE)提供的樣本空間是構(gòu)造近似數(shù)學(xué)模型的基礎(chǔ)。樣本空間包含了所要構(gòu)造的數(shù)學(xué)模型在不同空間位置的信息，如果試驗(yàn)設(shè)計(jì)不合適，數(shù)學(xué)模型中輸入與輸出的響應(yīng)精度將顯著降低，熱環(huán)境優(yōu)化設(shè)計(jì)也就失去了意義。目前，主要的試驗(yàn)設(shè)計(jì)方法包括部分因子設(shè)計(jì)、全因子設(shè)計(jì)、中心組合設(shè)計(jì)、拉丁超立方設(shè)計(jì)、優(yōu)化拉丁超立方設(shè)計(jì)(Optimal Latin hypercube design, Opt LHD)等。其中，優(yōu)化拉丁超立方設(shè)計(jì)其正交性良好、抽樣均勻、充滿(mǎn)度好等優(yōu)點(diǎn)，使其具有較強(qiáng)優(yōu)勢(shì)，其兩因素25抽樣點(diǎn)效果如圖9所示。本文最后滿(mǎn)足近似數(shù)學(xué)模型精度的Opt LHD樣本空間如表2所示，圖10給出了對(duì)應(yīng)的樣本空間分布圖。圖中顯示，Opt LHD抽樣具有很好滿(mǎn)意度。

圖9 兩因素25抽樣點(diǎn)OLHDFig.9 OLHD of 25 sampling points for two factors

圖10 三因素56抽樣點(diǎn)OLHDFig.10 OLHD of 56 sampling points for three factors

表2 Opt LHD樣本空間Table 2 Sample space of Opt LHD

4.2 近似數(shù)學(xué)模型

新型同心筒自力發(fā)射流場(chǎng)計(jì)算涉及到動(dòng)態(tài)網(wǎng)格，伴隨著控制體的運(yùn)動(dòng)、生成與消失，計(jì)算量很大，對(duì)其熱環(huán)境進(jìn)行優(yōu)化設(shè)計(jì)時(shí)，尋優(yōu)過(guò)程要耗費(fèi)巨大計(jì)算成本。基于合理的試驗(yàn)設(shè)計(jì)方案，建立精度可靠的響應(yīng)面近似數(shù)學(xué)模型，可消除優(yōu)化設(shè)計(jì)過(guò)程中的數(shù)值噪聲[14]，顯著降低計(jì)算成本。針對(duì)導(dǎo)彈熱環(huán)境評(píng)價(jià)指標(biāo)存在強(qiáng)非線性的問(wèn)題，本文采用四階響應(yīng)面構(gòu)造近似數(shù)學(xué)模型，不考慮三次及以上交叉項(xiàng)，其函數(shù)構(gòu)造形式為

(2)

式中n為設(shè)計(jì)變量個(gè)數(shù)；a0為常數(shù)項(xiàng)系數(shù)；ai、aii、aiii、aiiii分別為一次項(xiàng)到四次項(xiàng)系數(shù)；aij為二次交叉項(xiàng)系數(shù)。

式(2)中，各項(xiàng)系數(shù)通過(guò)最小二乘回歸分析求解，使得樣本值與預(yù)測(cè)值誤差平方和最小，則四階響應(yīng)面系數(shù)向量可表述為

a=(xTx)-1xTF

(3)

式中x為優(yōu)化拉丁超立方試驗(yàn)設(shè)計(jì)樣本空間設(shè)計(jì)變量向量，即x=[L1,L2,Φ]；F為樣本空間目標(biāo)向量，即F=[F1,F2]。

為確保后續(xù)優(yōu)化工作可靠，響應(yīng)面擬合精度必須滿(mǎn)足要求。評(píng)價(jià)響應(yīng)面擬合精度的指標(biāo)為R2：

(4)

經(jīng)過(guò)優(yōu)化平臺(tái)多次采樣、計(jì)算與擬合，近似數(shù)學(xué)模型得到了比較滿(mǎn)意的結(jié)果。圖11為觀測(cè)1面和觀測(cè)2面的近似數(shù)學(xué)模型擬合圖，其R2分別為0.975 6、0.967 8，且表3給出了觀測(cè)1面和2面近似數(shù)學(xué)模型的各項(xiàng)系數(shù)。

為驗(yàn)證近似數(shù)學(xué)模型精度，對(duì)觀測(cè)1面和觀測(cè)2面各隨機(jī)抽取5組模型進(jìn)行誤差分析，如表4所示。5組隨機(jī)誤差分析顯示，觀測(cè)1面、2面近似模型誤差較低，滿(mǎn)足進(jìn)一步優(yōu)化分析要求。

4.3 優(yōu)化策略

對(duì)于路基新型同心筒自力發(fā)射導(dǎo)彈熱環(huán)境優(yōu)化這種高度非線性的物理問(wèn)題，本文采用基于近似模型的多島遺傳算法(multi-island genetic algorithm MIGA)和序列二次規(guī)劃梯度算法相結(jié)合的組合優(yōu)化策略，如圖8所示。多島遺傳算法作為一種健壯穩(wěn)定的全局優(yōu)化算法能有效克服高度非線性難題，并對(duì)設(shè)計(jì)空間遍歷進(jìn)行全局搜索以定位敏感區(qū)域；作為一種成熟的梯度優(yōu)化算法，序列二次規(guī)劃[15]能充分發(fā)揮對(duì)敏感區(qū)域局部快速尋優(yōu)的優(yōu)勢(shì)，反復(fù)迭代，直到逼近全局最優(yōu)解。

(a) 觀測(cè)1面

(b) 觀測(cè)2面

表3 近似模型各項(xiàng)系數(shù)Table 3 Coefficients of the approximate model

多島遺傳算法通過(guò)將大種群分成若干個(gè)子種群(島)，在每個(gè)島上進(jìn)行基于傳統(tǒng)遺傳算法的選擇、交叉、變異等優(yōu)化操作，對(duì)子種群進(jìn)行優(yōu)化，提高優(yōu)化效率，便于找到敏感區(qū)域；通過(guò)遷移條件周期性執(zhí)行遷移操作，利于種群的多樣性，確保了全局最優(yōu)解的可靠性，其優(yōu)化歷程如圖12所示。

表4 觀測(cè)1、2面近似模型隨機(jī)誤差分析Table4 Random error analysis for the approximate model of observation plane 1 and plane 2

圖12 多島遺傳算法流程圖Fig.12 Principle flowchart of MIGA

5 優(yōu)化結(jié)果分析

圖 13給出了熱環(huán)境評(píng)價(jià)指標(biāo)的優(yōu)化歷程，從圖13(a)中多島遺傳算法的優(yōu)化歷程可看出，多島遺傳算法對(duì)熱環(huán)境指標(biāo)的敏感區(qū)域進(jìn)行快速局部定位，并呈現(xiàn)出局部震蕩到最優(yōu)解的趨勢(shì)，但無(wú)法對(duì)最優(yōu)解進(jìn)行最優(yōu)定位，主要發(fā)揮最優(yōu)區(qū)域定位的功能； NLPQL算法對(duì)最優(yōu)解非常敏感，具有很好的定位捕捉能力，但容易陷入局部最優(yōu)解；隨著多島遺傳算法對(duì)最優(yōu)區(qū)域的逐步縮小，NLPQL算法充分發(fā)揮了最優(yōu)解快速精確定位的特點(diǎn)，最終收斂到最優(yōu)解，充分體現(xiàn)了“組合優(yōu)勢(shì)”。

限于篇幅 ，各變量的優(yōu)化歷程沒(méi)有給出。表5給出了優(yōu)化前、后新型同心筒的熱結(jié)構(gòu)參數(shù)，以及熱環(huán)境優(yōu)化目標(biāo)值。表5顯示，從熱環(huán)境評(píng)價(jià)目標(biāo)的數(shù)值來(lái)看，優(yōu)化后的熱環(huán)境品質(zhì)提高了96.4%。

(a) MIGA優(yōu)化歷程

(b) NLPQL優(yōu)化歷程

表5 優(yōu)化結(jié)果對(duì)比Table 5 Comparison optical calculation results

圖14給出了優(yōu)化后新型同心筒優(yōu)化方案和其他2種方案觀測(cè)1面(距彈底1 m)的溫度時(shí)程曲線，各方案在0.085 s左右達(dá)到最高值，“引射效應(yīng)”結(jié)束，伴隨著內(nèi)筒冷氣體“倒吸效應(yīng)”的到來(lái)。為便于分析，圖15給出了新型優(yōu)化方案、新型基準(zhǔn)方案和新型隨機(jī)方案在此時(shí)刻彈壁面的溫度分布。圖16給出了新型優(yōu)化方案筒底部此時(shí)刻速度矢量流線圖。圖15顯示，從熱沖擊機(jī)理來(lái)看，優(yōu)化后的觀測(cè)1面(距彈底部1)溫度峰值320 K左右，幾乎沒(méi)被高溫燃?xì)饬髦苯記_擊，只有部分?jǐn)_動(dòng)氣流，導(dǎo)彈熱環(huán)境得到顯著改善；同時(shí)，從圖16中速度矢量流線圖可看出，彈底部的渦旋已經(jīng)無(wú)法影響導(dǎo)彈壁面，冷氣流的“倒吸效應(yīng)”已經(jīng)開(kāi)始。所以，燃?xì)饬骱诵牧鳠o(wú)法影響評(píng)價(jià)導(dǎo)彈熱環(huán)境指標(biāo)的觀測(cè)面1，相對(duì)于隨機(jī)方案和基準(zhǔn)方案，優(yōu)化后導(dǎo)彈的熱環(huán)境得到顯著改善。從圖16可看出，優(yōu)化后的熱結(jié)構(gòu)對(duì)筒底反濺流與渦流均具有很好的規(guī)避效果。

圖14 3種方案觀測(cè)1溫度時(shí)程曲線Fig.14 Temperature curves of the observation plane 1 in the three scheme

圖15 3種方案彈壁面0.085 s時(shí)刻溫度分布Fig.15 Temperature distribution of missile wall in three scheme at 0.085 s

圖16 優(yōu)化方案速度矢量Fig.16 Velocity vector of optimization schemes

6 結(jié)論

(1) 研究了新型路基同心筒熱環(huán)境特性，發(fā)射前期，內(nèi)筒出現(xiàn)“壅塞現(xiàn)象”導(dǎo)致的高溫燃?xì)庋刂鴱棻谙蛏线\(yùn)動(dòng)的現(xiàn)象，發(fā)射后期出現(xiàn)“倒吸效應(yīng)”；結(jié)合導(dǎo)彈壁面溫度“先上升、后下降、最后穩(wěn)定”的熱環(huán)境特性，確定以導(dǎo)彈觀測(cè)面溫差(相對(duì)于環(huán)境溫度)對(duì)時(shí)間的積分為目標(biāo)，來(lái)評(píng)價(jià)導(dǎo)彈熱環(huán)境品質(zhì)。

(2) 以?xún)?yōu)化拉丁超立方實(shí)驗(yàn)設(shè)計(jì)為理論基礎(chǔ)，建立了四階響應(yīng)面近似數(shù)學(xué)模型，避免了直接通過(guò)CFD計(jì)算進(jìn)行優(yōu)化搜索帶來(lái)的巨大計(jì)算量問(wèn)題，R2分析和隨機(jī)誤差分析顯示，代理數(shù)學(xué)模型精度較高，為導(dǎo)彈熱環(huán)境優(yōu)化設(shè)計(jì)提供了堅(jiān)實(shí)的技術(shù)支撐。

(3) 針對(duì)導(dǎo)彈熱環(huán)境在不同熱結(jié)構(gòu)條件下的強(qiáng)非線性問(wèn)題，基于多島遺傳全局優(yōu)化算法和序列二次規(guī)劃梯度算法，搭建了支持近似數(shù)學(xué)模型的導(dǎo)彈熱環(huán)境組合優(yōu)化構(gòu)架，為新型同心筒自力發(fā)射導(dǎo)彈熱環(huán)境友好性研究提供了一種新思路。

(4) 優(yōu)化后，熱環(huán)境評(píng)價(jià)目標(biāo)顯示，熱環(huán)境品質(zhì)提高了96.4%；發(fā)射前期，“壅塞現(xiàn)象”導(dǎo)致的高溫燃?xì)饬餮刂鴱棻谙蛏线\(yùn)動(dòng)的現(xiàn)象得到明顯抑制，發(fā)射后期，在冷氣流流入內(nèi)筒的“倒吸效應(yīng)”發(fā)生之前，評(píng)價(jià)導(dǎo)彈熱環(huán)境的觀測(cè)面幾乎不受燃?xì)饬骱诵牧鞯臒釠_擊，導(dǎo)彈熱沖擊機(jī)理發(fā)生良性改變，導(dǎo)彈總體熱環(huán)境得到顯著改善。

[1] 邵立武,姜毅,馬艷麗,等.新型艦載同心筒發(fā)射過(guò)程流場(chǎng)研究[J].導(dǎo)彈與航天運(yùn)載技術(shù),2011(4):54-58.

[2] 趙汝巖,黃志勇,周紅梅.潛載導(dǎo)彈近筒口點(diǎn)火數(shù)值仿真[J].固體火箭技術(shù),2012,35(2):161-165.

[3] 姜毅,郝繼光,劉群.同心筒垂直發(fā)射裝置排導(dǎo)燃?xì)饬鞯母倪M(jìn)[J].北京理工大學(xué)學(xué)報(bào),2007,27(2):95-98.

[4] 姜毅,郝繼光.新型“引射同心筒”垂直發(fā)射裝置理論及試驗(yàn)研究[J].宇航學(xué)報(bào),2008,29(1):236-241.

[5] 馬艷麗,姜毅,王偉臣,等.濕式同心筒自力垂直熱發(fā)射技術(shù)降溫效果研究[J].彈道學(xué)報(bào),2010,22(4):89-93.

[6] 于勇,母云濤.新型變截面同心筒發(fā)射裝置及其熱環(huán)境氣動(dòng)原理研究[J].宇航學(xué)報(bào),2013,34(9):1281-1287.

[7] 姜毅,郝繼光,劉群.同心筒垂直發(fā)射裝置排導(dǎo)燃?xì)饬鞯母倪M(jìn)[J].北京理工大學(xué)學(xué)報(bào),2007,27(2):95-98.

[8] 成沉,鮑福廷,劉旸.基于響應(yīng)面法的侯栓式噴管型面優(yōu)化設(shè)計(jì)[J].固體火箭技術(shù),2014,37(1):47-51.

[9] 陳偉,梁國(guó)柱.基于CFD的多級(jí)推力固體火箭發(fā)動(dòng)機(jī)軸對(duì)稱(chēng)噴管型面優(yōu)化與高精度性能預(yù)估[J].固體火箭技術(shù),2014,37(1):30-56.

[10] 粟華,谷良賢,龔春林.基于高擬真度模型的再入飛行器多學(xué)科優(yōu)化[J].西北工業(yè)大學(xué)學(xué)報(bào),2013,31(3):339-344.

[11] 邱東海,賈宏光,馬伍元,等.無(wú)人機(jī)地面隔振系統(tǒng)動(dòng)力學(xué)建模與優(yōu)化設(shè)計(jì)[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,42(2):121-126.

[12] 楊風(fēng)波,馬大為,任杰,等.新型車(chē)載同心筒流場(chǎng)機(jī)理與熱環(huán)境研究 [J].固體火箭技術(shù),2014,37(3):301-306.

[13] 李宗瑞.物料加熱與干燥過(guò)程的反向研究-傳熱與傳質(zhì)過(guò)程的最優(yōu)化 [D].沈陽(yáng):東北大學(xué),1991.

[14] 虞跨海,楊茜,倪俊,等.基于響應(yīng)面的渦輪葉片冷卻通道設(shè)計(jì)優(yōu)化 [J].航空學(xué)報(bào),2009,30(9):1630-1634.

[15] Gill P E,Murray W,Saunders M A.SNOPT: an SQP algorithm for large-scale constrained optimization [J].Society for Industrial and Applied Mathematics,2005,47(1):99-131.

(編輯：呂耀輝)

Optimization design for thermal environment of a new roadbed concentric canister launcher

YANG Feng-bo1，MA Da-wei1，REN Jie1，LE Gui-gao1，CAI De-yong2

(1. School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China; 2. 72465 Army,Jinan 250022, China)

Based on dynamic mesh technology with spring based smoothing method,and laying based zone moving method, the axisymmetric N-S equations were solved numerically, the dynamic thermal environment characteristics were obtained to deal with the thermal environment evaluating and optimization design problem of the new “middle diversion” concentric canister launcher (CCL),and evaluation index of thermal environment was also determined.The approximate mathematic model was established by optimal latin hypercube design and fourth-order response surface method to solve the automatic modeling problem of CFD and compensate for the shortcoming of large amount of calculation for direct optimization. A combinatorial optimization strategy platform based on Multi-Island Genetic Algorithm and Sequential Quadratic Programming was established to overcome the problem of strong nonlinear characteristic of the flow field parameters under different thermal structure conditions,and the optimization design system of thermal environment for missile which supports approximate mathematic model was also built. The comparison of the numerical results show that thermal environment of internal canister and external cylinder are improved by the cryogenic gas came from the cylinder port;the accuracy of approximate mathematic model is high enough to meet engineering requirements;after optimization, the thermal environment characteristics change positively,and the thermal environment of missile has been improved significantly.

concentric canister launcher(CCL);missile jet; thermal environment characteristics;approximate mathematic model;optimization design

2014-07-01；

：2014-09-20。

國(guó)防基礎(chǔ)科研資助項(xiàng)目(B2620110005)；國(guó)防預(yù)研項(xiàng)目(403050102)。

楊風(fēng)波(1987—)，男，博士生，研究方向?yàn)楸靼l(fā)射理論與技術(shù)。E-mail:yangfengbo.cool@163.com

任杰(1982—)，男，碩士生導(dǎo)師，研究方向?yàn)楸靼l(fā)射動(dòng)力學(xué)。E-mail:renjie@njust.edu.cn

V553.1

A

1006-2793(2015)02-0172-07

10.7673/j.issn.1006-2793.2015.02.005