基于不確定性的飛行器總體設(shè)計應(yīng)用

池元成 劉聞 李晶 王長慶

基于不確定性的飛行器總體設(shè)計應(yīng)用

池元成 劉聞 李晶 王長慶

(中國運載火箭技術(shù)研究院, 北京 100076)

飛行器總體設(shè)計在飛行器研制過程中具有重要地位。由于工程實際中存在的不確定性因素，需要在傳統(tǒng)確定性問題的基礎(chǔ)上，考慮不確定性因素的影響。因此，通過分析確定性優(yōu)化問題和不確定性優(yōu)化問題的流程，結(jié)合飛行器總體設(shè)計，分析了確定性優(yōu)化問題中不確定性因素的影響，并根據(jù)影響結(jié)果，建立了不確定性優(yōu)化問題，在滿足約束條件下，獲得了不確定性條件下最優(yōu)解，顯著提高了滿足約束的概率。

不確定性；飛行器總體設(shè)計；差分進(jìn)化算法；優(yōu)化設(shè)計

0 引言

飛行器總體設(shè)計是一項復(fù)雜系統(tǒng)工程，根據(jù)戰(zhàn)術(shù)技術(shù)要求和技術(shù)發(fā)展實際狀況，對飛行器及其組成飛行器的各系統(tǒng)進(jìn)行綜合、協(xié)調(diào)、研究、設(shè)計與試驗的過程，直接關(guān)系到飛行器的綜合性能、研制周期和研制經(jīng)費等，往往要經(jīng)過多次反復(fù)設(shè)計，才能得到一個綜合性能最佳的飛行器總體技術(shù)方案，因此在飛行器研制過程中具有重要地位[1]。傳統(tǒng)的飛行器總體設(shè)計過程中，往往忽略不確定因素，或者用參數(shù)的較大包絡(luò)去涵蓋不確定性因素的影響，采用確定性問題獲得飛行器初始方案，再通過多次迭代獲得可靠的總體方案。不確定性存在于飛行器設(shè)計的各個階段，按照不確定性的來源可分為隨機不確定性和認(rèn)知不確定性[2]。在飛行器總體設(shè)計過程中，識別不確定性因素并充分考慮不確定性影響，從而降低綜合性能對不確定影響的敏感程度，提高綜合性能的穩(wěn)健性和可靠性，國內(nèi)外在相關(guān)領(lǐng)域已開展了研究應(yīng)用[3-7]。飛行器設(shè)計優(yōu)化是從一組可行方案中選擇最佳可行方案的過程。即對于優(yōu)化問題，通過優(yōu)化算法求解最優(yōu)解。由于工程實際中存在的不確定性，不可能追求滿足要求的最優(yōu)解，而是退而求其次，追求滿足要求的次優(yōu)解。本文對某飛行器總體方案進(jìn)行了優(yōu)化設(shè)計，再利用不確定性方法開展了總體可靠性設(shè)計，并對結(jié)果數(shù)據(jù)進(jìn)行了分析。

1 不確定性優(yōu)化問題

傳統(tǒng)的確定性優(yōu)化問題可描述為

其中，f為目標(biāo)函數(shù)，g為約束條件，x為設(shè)計變量，p為設(shè)計參數(shù)。確定性優(yōu)化流程是由優(yōu)化算法生成設(shè)計變量初始值，由設(shè)計變量計算目標(biāo)函數(shù)和約束條件，在判斷是否滿足停止條件，不滿足則返回，繼續(xù)由優(yōu)化算法生成下一步設(shè)計變量初始值，如此往復(fù)，直到滿足停止條件，選取最優(yōu)解，如圖1所示。

然而，在工程實際中設(shè)計變量和設(shè)計參數(shù)中或多或少的存在不確定性，如飛行載荷、材料屬性、制造工藝等，即數(shù)學(xué)上服從某種分布的隨機變量。此外，由于約束條件的限制，工程實際確定性優(yōu)化問題的優(yōu)化解往往在約束邊界上，即起作用約束的存在，使得確定性優(yōu)化問題的優(yōu)化解容易違反約束，即違反約束的概率較高，導(dǎo)致優(yōu)化解不可靠，這也是實際工程上對優(yōu)化設(shè)計方法不可信的原因之一。

因此，工程中迫切需要研究不確定性優(yōu)化設(shè)計方法，在建模過程中引入不確定性因素，通過不確定性分析，對確定性優(yōu)化解進(jìn)行適當(dāng)調(diào)整，也就是以較高的概率不違反約束條件，從而獲得滿足工程實際需求的較優(yōu)解。當(dāng)考慮設(shè)計變量和設(shè)計參數(shù)的不確定性時，不確定性優(yōu)化問題可描述為

其中，μ為均值，σ為標(biāo)準(zhǔn)差，Pr為概率。不確定性優(yōu)化流程如圖2所示。由圖1和圖2可知，不確定性優(yōu)化流程是雙層嵌套的流程，相比于確定性優(yōu)化流程，不確定性優(yōu)化流程多了不確定性分析的循環(huán)流程，用于計算概率。不確定性分析過程中，對于數(shù)值計算問題，一般采用蒙特卡洛計算。根據(jù)模型的復(fù)雜度，可選擇不同的不確定性量化方法進(jìn)行分析，如對于計算復(fù)雜度高的模型，可使用近似模型代替復(fù)雜模型，提高計算效率。

2 差分進(jìn)化算法

差分進(jìn)化（Differential Evolution，DE）[8,9]算法作為群智能算法的代表算法之一，能夠?qū)崿F(xiàn)全局優(yōu)化。DE不同于遺傳算法（Genetic Algorithm，GA）使用二進(jìn)制編碼，而是直接采用實數(shù)編碼，省去了不必要的麻煩。DE通過父代個體間的簡單數(shù)學(xué)組合來產(chǎn)生不同的候選個體，再通過比較個體間的函數(shù)適應(yīng)值的優(yōu)劣關(guān)系確定候選個體是否取代父代個體?，F(xiàn)已證明DE較其他現(xiàn)有的進(jìn)化算法是有效的、魯棒的。由于DE簡單易實現(xiàn)，可調(diào)參數(shù)少，目前在很多領(lǐng)域的復(fù)雜優(yōu)化問題求解中得到了應(yīng)用。

DE算法基于父代個體間的差異，經(jīng)變異、交叉和選擇產(chǎn)生子代，根據(jù)在搜索空間內(nèi)隨機選取的兩個個體向量的差值與第三個隨機選取的個體向量的加權(quán)求和來實現(xiàn)群體變異，然后通過比較由交叉算子生成的個體和相應(yīng)父代個體的函數(shù)值優(yōu)劣來保存優(yōu)秀個體，如此反復(fù)迭代，不斷進(jìn)化。定義維搜索空間內(nèi)，存在種群數(shù)為的個體i=(x1,x2,…,x) (=1,2, …,)，此時變異個體V=(v1,v2, …,v)通過式(3)產(chǎn)生

式中，1、2和3是異于且互不相同的[1,]之間的整數(shù)；為縮放比例因子；為當(dāng)前迭代步數(shù)。式(3)中，X1稱為基向量，(X2-X3)稱為差分向量。

父代個體X和變異個體V依據(jù)交叉概率生成實驗個體U=(u1,u2, …,u)

式中，()是[0,1]之間的均勻隨機分布數(shù)。

對于最小化問題，交叉后的個體U+1和父代個體X依據(jù)式(5)進(jìn)行選擇操作，生成子代個體X+1

式中，為求解問題的目標(biāo)函數(shù)值。

3 不確定性優(yōu)化設(shè)計方法

結(jié)合上述的內(nèi)容，對于工程實際問題，不必直接采用不確定性分析方法完成優(yōu)化設(shè)計，可采用如下的步驟完成不確定性分析與優(yōu)化。

步驟1：按式(1)選定目標(biāo)函數(shù)和約束條件，確定設(shè)計變量，即設(shè)定確定性優(yōu)化問題；

步驟2：按圖1所示的流程求解確定性優(yōu)化問題，獲得滿足約束的最優(yōu)解；

步驟3：對于最優(yōu)解完成不確定性分析，若滿足則選擇最優(yōu)解完成方案設(shè)計；否則，轉(zhuǎn)步驟4；

步驟4：按式（2）設(shè)定不確定性優(yōu)化問題；

步驟5：按圖2所示的流程求解不確定性優(yōu)化問題，獲得滿足約束的最優(yōu)解，通過分析作為最終方案。

4 優(yōu)化模型與結(jié)果分析

某飛行器總體設(shè)計涉及動力、控制、結(jié)構(gòu)等多個專業(yè)，本文通過結(jié)構(gòu)設(shè)計矩陣梳理了各專業(yè)之間的參數(shù)傳遞關(guān)系與具體傳遞參數(shù)，并搭建了飛行器總體方案設(shè)計流程，其設(shè)計結(jié)構(gòu)矩陣如圖3所示。

圖3 總體設(shè)計結(jié)構(gòu)矩陣

為了達(dá)到設(shè)計要求，以飛行距離（）最大作為優(yōu)化目標(biāo)，構(gòu)建了確定性優(yōu)化問題

（6）

其中，為設(shè)計變量，選取了燃料質(zhì)量等3個關(guān)鍵變量。為簡化研究，文中令三個優(yōu)化設(shè)計變量均服從正態(tài)分布，變異系數(shù)為0.005，如表1所示。

表1 優(yōu)化設(shè)計變量

利用差分進(jìn)化算法完成優(yōu)化求解后，得到了表2所示的最優(yōu)解。

表2 確定性優(yōu)化結(jié)果

由表2可知，最優(yōu)解的約束條件離邊界非常接近，通過分析可知，違反約束條件1的概率為0.46，違反約束條件2的概率為0.50，最優(yōu)解的可靠性相對較低，采用風(fēng)險較大。

圖4 約束條件g1概率分布

圖5 約束條件g2概率分布

為此，將上述問題轉(zhuǎn)化為不確定性優(yōu)化問題

(7)

其中，滿足約束的概率大于90%。

考慮不確定性條件下，利用差分進(jìn)化算法完成優(yōu)化求解后，得到了表3所示的最優(yōu)解。

表3 不確定性優(yōu)化結(jié)果

由表3可知，考慮不確定性的條件下，最優(yōu)解相比確定性優(yōu)化解，稍差0.05%，然而約束條件1和約束條件2的違反概率分別降低為0.0073和0.0063，可知最優(yōu)解可靠性提高顯著，可以采用。

圖6 約束條件g1概率分布

圖7 約束條件g2概率分布

3 結(jié)論

本文以某飛行器總體設(shè)計為對象，利用不確定性優(yōu)化設(shè)計方法對其進(jìn)行了優(yōu)化分析，獲得了滿足約束條件的較優(yōu)解。然而，在工程實際中還需要在以下幾個方面開展深入的研究。

1）飛行器總體設(shè)計過程是一個多輪次迭代的過程，飛行器總體方案由粗到細(xì)，在使用不確定性方法時，還需要進(jìn)一步分析在什么環(huán)節(jié)考慮哪些不確定性因素；

2）飛行器總體設(shè)計各學(xué)科專業(yè)建模時，要充分考慮不確定性，除了隨機不確定性，更要關(guān)注認(rèn)知不確定性，對工程實際問題完整描述；

3）學(xué)科專業(yè)模型的精度或置信度要提高，當(dāng)模型精度較低時，計算結(jié)果會覆蓋不確定性因素的影響；

4）當(dāng)優(yōu)化問題中涉及的不確定性因素較多時，要事先完成不確定性參數(shù)靈敏度分析；當(dāng)模型計算量較大時，要考慮構(gòu)建近似模型。

[1] 侯世明.導(dǎo)彈總體設(shè)計與試驗[M].北京: 宇航出版社, 1996.

[2] 陳小前.飛行器不確定性多學(xué)科設(shè)計優(yōu)化理論與應(yīng)用[M].北京: 科學(xué)出版社, 2013.

[3] 李曉斌, 張為華, 王中偉.固體火箭發(fā)動機裝藥不確定性優(yōu)化設(shè)計[J].固體火箭技術(shù), 2006, 29(4): 269-273.[Li Xiaobin, Zhang Weihua, Wang Zhongwei.Uncertainty optimization design of solid rocket motor grain [J].Journal of Solid Rocket Technology, 2006, 29(4): 269-273.]

[4] 王宇, 余雄慶.考慮不確定性的飛機總體參數(shù)優(yōu)化方法[J].航空學(xué)報, 2009, 30(10): 1883-1888.[Wang Yu, Yu Xiongqing.Optimization method for aircraft conceptual design under uncertainty, 2009, 30(10): 1883-1888.]

[5] Fuchs M, Neumaier A, Girimonte D.Uncertainty modeling in autonomous robust spacecraft system design[J].PAMM, 2007, 7(1):2060041-2060042.

[6] Green L L, Lin HZ, Khalessi M R.Probabilistic methods for uncertainty propagation applied to aircraft design [R].AIAA-2002-3140, 2002.

[7] 沈鳳霞.力學(xué)環(huán)境試驗測量不確定度分析[J].強度與環(huán)境, 2010, 37(1): 58-64.[Shen Fengxia.Measurement uncertainty analysis of dynamic environment test [J].Structure & Environment engineering, 2010, 37(1): 58-64.]

[8] Storn R, Price K.Differential evolution-a simple and efficient adaptive scheme for global optimization over continuous spaces, TR-95-012 [R].Berkeley: International Computer Science Institute, 1995.

[9] Storn R, Price K.Minimizing the real function of the ICEC’96 contest by DE [C].IEEE International Conference on Evolutionary Computation, 1996: 842-844.

Application of Aircraft Overall Design based Uncertainty

CHI Yuan-cheng LIU Wen LI Jing WANG Chang-qing

(China Academy of Launch Vehicle Technology, Beijing, 100076, China)

The overall design of aircraft plays an important role in the process of aircraft development.Due to the uncertain factors in practical engineering, it is necessary to consider the influence of uncertain factors on the basis of traditional deterministic problems.Therefore, by analyzing the process of deterministic optimization problem and uncertain optimization problem, combined with the overall design of aircraft, the influence of uncertain factors in deterministic optimization problem is analyzed, and according to the influence results, the uncertain optimization problem is established.Under the constraint conditions, the optimal solution under the uncertainty conditions is obtained, which significantly improves the probability of meeting the constraints.

Uncertainty; Aircraft overall design; Differential evolution; Optimize design

V211.4

A

1006-3919(2021)06-0045-05

10.19447/j.cnki.11-1773/v.2021.06.007

2021-05-29；

2021-08-03

國家重點研發(fā)計劃項目（2020YFB1709400）

池元成（1981—），男，博士，高級工程師，研究方向：飛行器總體設(shè)計與優(yōu)化；（100076）北京9200信箱38分箱.