劉明航, 胡峪
(西北工業(yè)大學(xué) 航空學(xué)院, 陜西 西安 710072)
基于自適應(yīng)近似模型的RSCSSO及其應(yīng)用研究
劉明航, 胡峪
(西北工業(yè)大學(xué) 航空學(xué)院, 陜西 西安 710072)
為了更好地解決飛行器設(shè)計中多學(xué)科設(shè)計優(yōu)化問題,在對傳統(tǒng)的基于響應(yīng)面RSCSSO研究的基礎(chǔ)上,提出了改進的基于自適應(yīng)近似模型的RSCSSO。為減少學(xué)科分析的次數(shù),采用均勻?qū)嶒炘O(shè)計代替學(xué)科級優(yōu)化,直接獲得性能優(yōu)良的初始設(shè)計樣本點;在系統(tǒng)級優(yōu)化過程中引入自適應(yīng)近似模型算法,并通過在迭代過程中對兩種近似模型進行對比,選用精度高的近似模型提高系統(tǒng)級優(yōu)化效率。采用改進的RSCSSO解決了一架飛翼布局手拋無人機設(shè)計優(yōu)化問題,并與傳統(tǒng)RSCSSO進行了對比。結(jié)果表明,改進的RSCSSO不但有更小的計算量,而且獲得了更優(yōu)的結(jié)果。
多學(xué)科設(shè)計優(yōu)化; 并行子空間優(yōu)化算法; 神經(jīng)網(wǎng)絡(luò); 遺傳算法
飛行器設(shè)計過程是一個多學(xué)科參與的反復(fù)迭代過程,通常單獨學(xué)科的最優(yōu)解并不是整個系統(tǒng)的最優(yōu)解。因此首要的問題是如何有效地解決各學(xué)科間的信息交流,以及由于學(xué)科分析產(chǎn)生的龐大的計算量問題。
為了更好地解決這類復(fù)雜的系統(tǒng)優(yōu)化問題,本文對傳統(tǒng)的RSCSSO進行了研究和改進,提出了基于自適應(yīng)近似模型的RSCSSO。采用改進的優(yōu)化算法對一架小型飛翼布局手拋無人機進行了設(shè)計優(yōu)化,并與傳統(tǒng)的優(yōu)化算法相比較,驗證了所提算法的有效性和實用性。
目前,為了解決復(fù)雜的耦合問題,提出了許多MDO算法[1]。1988年,Sobieszczanski-Sobieski[2]提出的并行子空間算法(CSSO),在每個子空間分別優(yōu)化一組獨立的設(shè)計變量,分解設(shè)計變量可以起到有效的降維作用,然后利用全局敏度方程(GSE)在系統(tǒng)級對這些設(shè)計變量進行協(xié)調(diào)。CSSO能夠讓各學(xué)科專家根據(jù)本學(xué)科的實際情況選擇適合的分析、計算工具和優(yōu)化方法,設(shè)計自由度高、人工干預(yù)性強;學(xué)科(子空間)劃分與設(shè)計部門劃分一致,便于組織協(xié)調(diào),能夠靈活有效地解決復(fù)雜的系統(tǒng)問題[3]。然而,由于CSSO 算法是基于GSE的線性近似,因而子空間設(shè)計變量的可行域較小。因為在實際設(shè)計問題中,有些設(shè)計變量同時對幾個子空間均有影響,子空間中設(shè)計變量互不重疊的要求不太合理。在求解GSE方程時需要對變量求導(dǎo),無法求解離散問題[4],這些都局限了它的應(yīng)用。1996年,RSCSSO在此基礎(chǔ)上被提出,RSCSSO允許子空間設(shè)計變量共享,同一設(shè)計變量可以在幾個子空間進行優(yōu)化設(shè)計,并且能夠解決離散型問題,適用范圍更廣。
RSCSSO算法流程圖如圖1所示。
圖1 RSCSSO算法流程圖Fig.1 The RSCSSO algorithm structure flowchart
通過對RSCSSO算法的研究分析,可以看出:
(1)收斂準則和優(yōu)化結(jié)果的好壞很大程度上依賴于系統(tǒng)級優(yōu)化,而與學(xué)科級優(yōu)化關(guān)系不大。
(2)系統(tǒng)級優(yōu)化和學(xué)科級優(yōu)化的物理意義是一致的,只是優(yōu)化的設(shè)計變量不完全一致。
(3)根據(jù)數(shù)據(jù)庫建立的近似模型精度直接影響系統(tǒng)級優(yōu)化過程的質(zhì)量。
以上三點表明,學(xué)科級優(yōu)化的目的僅在于提供性能優(yōu)良(更加靠近全局最優(yōu)解)的樣本點,用這些點能夠構(gòu)建更加精確的近似模型,尤其在接近最優(yōu)解的小范圍內(nèi)精確度更高。最后一點表明了近似模型的精度直接決定了收斂過程的快慢和優(yōu)化算法的效率。因此,針對這些特點,本文對該算法進行了改進,提出了基于自適應(yīng)近似模型的RSCSSO。
根據(jù)對傳統(tǒng)RSCSSO的研究,由于近似模型的精度主要由兩個因素決定,一個是初始樣本點的質(zhì)量,另一個是所選近似模型是否合適,因此對RSCSSO進行了如下改進:
(1)在RSCSSO算法中,用均勻?qū)嶒炘O(shè)計代替學(xué)科級優(yōu)化,直接獲得性能優(yōu)良的初始設(shè)計樣本點,學(xué)科級只進行學(xué)科分析,以減小計算量[5]。
(2)二次響應(yīng)面和徑向神經(jīng)網(wǎng)絡(luò)是兩種常用的近似模型[6],由于對不同問題的近似模型的近似精度不同,所以本文提出自適應(yīng)的近似模型方法。
實驗設(shè)計(DOE)是多學(xué)科優(yōu)化設(shè)計的一個重要組成部分,主要目標是用最少的實驗次數(shù),建立設(shè)計變量和狀態(tài)變量之間的關(guān)系。均勻?qū)嶒炈@得的樣本點與正交實驗相比,分布更加均勻、更具有代表性,所以用均勻?qū)嶒炘O(shè)計的方法來直接獲得性能良好的初始設(shè)計樣本點,用以構(gòu)造響應(yīng)面[7]。
(1)
(2)
(3)
改進后的基于自適應(yīng)近似模型的RSCSSO流程如圖2所示。
改進的算法步驟如下:
(1)由均勻?qū)嶒炘O(shè)計(UED)獲得初始設(shè)計樣本點(同時獲得一組測試點)并進行系統(tǒng)分析,得到與之對應(yīng)的一組狀態(tài)變量,構(gòu)成初始數(shù)據(jù)庫;
(3)系統(tǒng)級優(yōu)化根據(jù)步驟(2)中選擇的近似模型進行,對得到的最優(yōu)解進行系統(tǒng)分析并判斷收斂準則。如果收斂則優(yōu)化結(jié)束,否則將本次迭代得到的最優(yōu)解加入數(shù)據(jù)庫然后跳轉(zhuǎn)至步驟(2)。
圖2 基于自適應(yīng)近似模型的RSCSSOFig.2 Improved RSCSSO based on self-adaptive approximation model
利用改進的基于自適應(yīng)近似模型的RSCSSO進行了一架飛翼布局手拋無人機的優(yōu)化設(shè)計。首先,將優(yōu)化設(shè)計問題解耦為一個系統(tǒng)級優(yōu)化和氣動、重量、性能三個學(xué)科級優(yōu)化。在重量學(xué)科中,根據(jù)外形尺寸,將利用幾何積分來計算各部分重量的方法作為該學(xué)科的精確計算模型,完成學(xué)科分析。Tornado (Vortex Lattice Method)用作氣動學(xué)科的精確計算模型[8],根據(jù)文獻[9-10]建立性能學(xué)科精確計算模型,根據(jù)文獻[10]建立電池模型。
將航時T作為優(yōu)化目標尋找其最大值,設(shè)計變量和設(shè)計參數(shù)分別列于表1和表2。
表1 設(shè)計變量Table 1 Design variables
表2 設(shè)計參數(shù)Table 2 Design parameters
其中,VLOF為起飛速度;η為動力系統(tǒng)總效率(包括電機調(diào)速器效率、電機效率和螺旋槳效率);Wload為有效裝載重量;Λ為四分之一弦線后掠角;τ為幾何扭轉(zhuǎn)角;fspan為襟翼展長占機翼展長的比例;fcho為襟翼弦長占當(dāng)?shù)貦C翼弦長的比例;K1為襟翼距翼根的距離占機翼展長的比例;aspan為副翼展長占機翼展長的比例;acho為副翼弦長占當(dāng)?shù)貦C翼弦長的比例;K2為副翼距翼稍的距離占機翼展長的比例。
優(yōu)化問題模型為:
尋找X=[AλSVcruH]T
MaxT=f(X)
fspan+aspan+K1+K2<1
Re≤25 000
解耦后優(yōu)化問題及設(shè)計變量學(xué)科分布情況如圖3所示??梢钥吹?設(shè)計變量被分配到三個學(xué)科中,氣動學(xué)科狀態(tài)變量升阻比L/D同時也是重量學(xué)科和性能學(xué)科分析所必須的設(shè)計參數(shù)。與此同時,起飛重量WTO是重量學(xué)科的狀態(tài)變量,也是性能學(xué)科的輸入量。利用這幾個狀態(tài)變量集中反映學(xué)科間的耦合情況,起到了解耦的作用。
圖3 解耦后優(yōu)化問題及設(shè)計變量學(xué)科分布情況Fig.3 Decoupled optimization problem and design variables distribution
基于梯度的優(yōu)化算法雖然比基于隨機性的優(yōu)化算法更精確一些,但是在求解過程中需要求解最優(yōu)解附近的導(dǎo)數(shù)值容易陷入局部最優(yōu)解?;谔荻葘?dǎo)數(shù)的優(yōu)化算法只有初值取得合適時,才不會陷入局部最優(yōu),找到合理的最優(yōu)解。遺傳算法(GA)是一種基于種群的非梯度隨機性尋優(yōu)算法,并且對于離散型和高階非線性問題均適用,尋優(yōu)效率高,不依賴于初值的取值。本文將GA算法作為每次迭代過程中系統(tǒng)級優(yōu)化的算法。GA算法流程如圖4所示。
圖4 GA算法流程圖Fig.4 Flowchart of GA
針對無人機優(yōu)化設(shè)計問題,本文分別用傳統(tǒng)的RSCSSO和基于自適應(yīng)近似模型的RSCSSO算法進行了計算。自適應(yīng)過程列于表3,計算結(jié)果列于表4。
從表4可以看出,基于自適應(yīng)近似模型的RSCSSO得到了更好的結(jié)果。不但升阻比有略微升高,而且由于翼面積的減少重量明顯減輕,使得航時提升3.67%。
表3 狀態(tài)變量自適應(yīng)過程Table 3 Self-adaptive procedure for the state variables
表4 優(yōu)化結(jié)果比較Table 4 Comparison of the optimization results
兩種算法的迭代收斂過程如圖5所示。可以看出,改進的RSCSSO算法比傳統(tǒng)RSCSSO算法收斂更快、計算量更少。所以,改進的算法更有效、更具實用性。
圖5 迭代收斂過程對比Fig.5 Comparison of the iterative convergence process
本文研究了傳統(tǒng)的RSCSSO算法,在此基礎(chǔ)上提出了基于自適應(yīng)近似模型的RSCSSO,并對一架小型飛翼布局手拋無人機進行了優(yōu)化設(shè)計,得出以下結(jié)論:
(1)由于自適應(yīng)近似模型的引入,采用精度更好的近似模型,使得改進的RSCSSO算法得到了更好的最優(yōu)解。
(2)利用均勻?qū)嶒炘O(shè)計(UED)代替學(xué)科級優(yōu)化過程,減少了學(xué)科分析次數(shù)、降低了計算量,使得所提算法更加有效、收斂更快。
(3)對比優(yōu)化結(jié)果可以發(fā)現(xiàn),改進的優(yōu)化算法得到的優(yōu)化結(jié)果升阻比略微提高、重量減輕,從而使得航時得到提升,學(xué)科間綜合分析完善,更易得到全局最優(yōu)解。
[1] Wit A J D,Keulen F V.Overview of methods for multi-level and/or multi-disciplinary optimization[R].AIAA-2010-2914,2010.
[2] Sobieszczanski-Sobieski J.A linear decomposition method for large optimization problems-blueprint for development [R].NASA-TM-83248,1982.
[3] Sobieszczanski-Sobieski J.Optimization by decomposition:a step from hierarchic to non-hierarchic systems[R].NASA-CP-3031,1988.
[4] 余雄慶.多學(xué)科設(shè)計優(yōu)化算法及其在飛行器設(shè)計中的應(yīng)用[J].航空學(xué)報,2000,21(1):2-5.
[5] 張科施.飛機設(shè)計的多學(xué)科優(yōu)化方法研究[D].西安:西北工業(yè)大學(xué),2006.
[6] 張健,李為吉.飛機多學(xué)科設(shè)計優(yōu)化中的近似方法分析[J].航空計算技術(shù),2005,35(3):1-2.
[7] 方開泰,馬長興.正交與均勻試驗設(shè)計[M].北京:科學(xué)出版社,2001:35-152.
[8] Melin T.A vortex lattice MATLAB implementation for linear aerodynamic wing applications [D].Sweden:Royal Institute of Technology,2000.
[9] Wagner N,Boland S.Powertrain design for hand-launchable long endurance unmanned aerial vehicles[R].AIAA-2011-6047,2011.
[10] Traub L W.Range and endurance estimates for battery-powered aircraft[J].Journal of Aircraft,2011,48(2):703-707.
ResearchontheRSCSSObasedonself-adaptiveapproximationmodelanditsapplication
LIU Ming-hang, HU Yu
(College of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China)
To solve the problem of aircraft design in multidisciplinary design optimization efficiently, the traditional concurrent subspace optimization based on response surface method (RSCSSO) was studied. Then the RSCSSO based on self-adaptive approximation model was developed. In order to reduce the amount of disciplinary analysis, uniform experiment design is introduced to replace the disciplinary optimization to obtain a set of design points directly and a self-adaptive approximation algorithm is introduced in system level optimization. The accuracies of two approximation models are compared in each iteration and the better model will be used. The improved RSCSSO algorithm was validated by a hand-launch UAV test case. In comparison with the traditional RSCSSO optimization, the better results with much less computation cost was found by the improved RSCSSO.
multidisciplinary design optimization; concurrent subspace optimization; neural network; genetic algorithm
V221
A
1002-0853(2013)06-0496-05
2012-12-04;
2013-04-02; < class="emphasis_bold">網(wǎng)絡(luò)出版時間
時間:2013-08-21 16:13
劉明航(1989-),男,陜西咸陽人,碩士研究生,主要研究方向為飛機總體設(shè)計、多學(xué)科設(shè)計優(yōu)化;
胡峪(1974-),男,四川溫江人,副教授,博士,主要研究方向為飛機總體設(shè)計、微小型無人機多學(xué)科設(shè)計優(yōu)化。
(編輯:李怡)