王 濤

（南昌理工學(xué)院,南昌 330013）

1 機(jī)翼氣動(dòng)優(yōu)化平臺(tái)的建立

根據(jù)機(jī)翼氣動(dòng)優(yōu)化分析的需要，要求完成三維機(jī)翼模型的自動(dòng)建立、機(jī)翼網(wǎng)格的自動(dòng)劃分、流體計(jì)算以及系統(tǒng)優(yōu)化四個(gè)方面的任務(wù)。優(yōu)化平臺(tái)的設(shè)計(jì)簡(jiǎn)化框圖如圖1所示：

機(jī)翼模型的構(gòu)建機(jī)翼的網(wǎng)格劃分機(jī)翼流體力學(xué)計(jì)算Isight控制數(shù)據(jù)流的傳遞。

圖1 機(jī)翼氣動(dòng)優(yōu)化平臺(tái)簡(jiǎn)化框圖

機(jī)翼模型的構(gòu)建模塊：主要是完成對(duì)三維機(jī)翼的自動(dòng)修形，即根據(jù)Isight優(yōu)化算法的要求，自動(dòng)的調(diào)整諸如弦長(zhǎng)、迎角之類的外形參數(shù)，構(gòu)建新的機(jī)翼外形。

機(jī)翼網(wǎng)格劃分模塊：主要負(fù)責(zé)對(duì)機(jī)翼外形的網(wǎng)格劃分。

機(jī)翼流體力學(xué)計(jì)算模塊：負(fù)責(zé)采用合適的數(shù)值模擬方法計(jì)算出各項(xiàng)氣動(dòng)特性參數(shù)。

Isight模塊：負(fù)責(zé)各項(xiàng)任務(wù)銜接和引導(dǎo)，并根據(jù)流體計(jì)算結(jié)果和一定的優(yōu)化搜索策略，調(diào)整機(jī)翼的外形參數(shù)，到達(dá)自動(dòng)尋優(yōu)的目的。

鑒于以上的分析，在Isight軟件里建立了如下的優(yōu)化平臺(tái)。

圖2 Isight中的機(jī)翼氣動(dòng)優(yōu)化平臺(tái)

從圖2中可以看出，goug是完成參數(shù)化建模的，gogambit用來(lái)實(shí)現(xiàn)網(wǎng)格劃分，流體計(jì)算則是由gofluent來(lái)完成、Calculatvor用來(lái)完成對(duì)研究目標(biāo)的二次計(jì)算，如計(jì)算升阻比CL/CD等。DOE用來(lái)完成實(shí)驗(yàn)設(shè)計(jì)，根據(jù)實(shí)驗(yàn)設(shè)計(jì)的結(jié)果，選擇合適的設(shè)計(jì)點(diǎn)，再進(jìn)行優(yōu)化設(shè)計(jì)。

機(jī)翼氣動(dòng)優(yōu)化平臺(tái)建立的具體思路步驟見(jiàn)圖2。在該圖中，包含了參數(shù)化建模、網(wǎng)格劃分到氣動(dòng)分析所需的各類信息，同時(shí)給出了各模塊在iSIGHT上的具體的輸入輸出參數(shù)文件。其中iSIGHT主要是實(shí)現(xiàn)對(duì)數(shù)據(jù)的自動(dòng)調(diào)用，按照某種預(yù)先設(shè)計(jì)的算法，完成對(duì)數(shù)據(jù)的流的反復(fù)導(dǎo)入，當(dāng)達(dá)到最終設(shè)定的優(yōu)化目標(biāo)時(shí)，則停止運(yùn)算，并輸出最優(yōu)結(jié)果。

Feiji.prtFeiji.expFeiji.x_tGambit.jouFluent.jouFeiji.mshFeiji.datUG構(gòu)建機(jī)翼模型Gambit生成翼型網(wǎng)格Fluent計(jì)算得到Cl，CdiSIGHT控制數(shù)據(jù)流程Ug二次開(kāi)發(fā)程序

圖3 優(yōu)化平臺(tái)建立的具體思路步驟

2 目標(biāo)優(yōu)化

為了驗(yàn)證機(jī)翼氣動(dòng)優(yōu)化的效果，我們選取了典型的joukowski翼型作為研究對(duì)象，首先通過(guò)對(duì)某類翼型的參數(shù)化，然后由ISIGHT負(fù)責(zé)自動(dòng)計(jì)算，實(shí)現(xiàn)參數(shù)化建模的自動(dòng)化。

設(shè)計(jì)變量：提取根弦長(zhǎng)250mm、梢弦長(zhǎng)50mm、半展長(zhǎng)1000mm和前緣后掠角11.5°作為設(shè)計(jì)變量。

其中來(lái)流速度V0=50m/s，攻角為θ=5°；假設(shè)流體為定常的；來(lái)流空氣物理特性取標(biāo)準(zhǔn)海平面的值：壓強(qiáng)為101,325 Pa ，密度= 1.2250 kg/m3，溫度為288.16 K，運(yùn)動(dòng)學(xué)粘性系數(shù)v = 1.4607e-5 m2/s。

目標(biāo)函數(shù)：升阻比CL/CD最大。

優(yōu)化求解：在本算例中，試驗(yàn)設(shè)計(jì)采用拉丁超立方法獲得100個(gè)樣本點(diǎn)。該方法允許實(shí)驗(yàn)樣本點(diǎn)在設(shè)計(jì)有效域內(nèi)進(jìn)行均勻化的排布。該算法和ISIGHT的結(jié)合將可以較好的進(jìn)行反復(fù)的目標(biāo)函數(shù)測(cè)算，并且可以有效減少因?yàn)閷?shí)驗(yàn)樣本點(diǎn)選擇不合理帶來(lái)的優(yōu)化誤差。由于選擇了100個(gè)樣本點(diǎn)，因此基于ISIGHT的集成系統(tǒng)會(huì)對(duì)100多個(gè)不同參數(shù)的模型進(jìn)行測(cè)算。并對(duì)每一個(gè)模型進(jìn)行網(wǎng)格化（通過(guò)調(diào)用Gambit軟件），最后進(jìn)行流體計(jì)算（通過(guò)調(diào)入Fluent軟件）。最后得出100組升阻比。

優(yōu)化結(jié)果如表所示。

在上表中，只給出最終的優(yōu)化結(jié)果，從表中的數(shù)據(jù)可以看出，通過(guò)對(duì)三維機(jī)翼的優(yōu)化，機(jī)翼的升阻比有一定的提高。在最優(yōu)升阻比數(shù)據(jù)對(duì)應(yīng)的三維模型，也即我們想要得到的模型。

優(yōu)化結(jié)果

3 總結(jié)

本章應(yīng)用CAD二次開(kāi)發(fā)技術(shù)和參數(shù)化建模方法，充分利用了CAD軟件豐富的建模功能，通過(guò)和ISIGHT集成軟件的并行工作，最終提升了設(shè)計(jì)效率。

[1]王濤，彭如.多學(xué)科協(xié)同優(yōu)化設(shè)計(jì)技術(shù)的研究初探[J].電子制作，2014(01).