張 祥 劉 勇 趙 璐 李亞清
(南京航空航天大學(xué)能源與動力學(xué)院航空發(fā)動機熱環(huán)境與熱結(jié)構(gòu)工業(yè)和信息化部重點實驗室 江蘇 南京 210016)
渦輪是發(fā)動機中動力負荷與熱負荷最大的部件,采用有效的冷卻措施是發(fā)動機安全可靠工作的有效保證,也是降低高溫材料成本的措施[1]。隨著航空發(fā)動機性能發(fā)展的需求,渦輪葉片的冷卻結(jié)構(gòu)越來越復(fù)雜[2]。
工程上通常采用數(shù)值模擬技術(shù)預(yù)測渦輪葉片表面溫度分布以計算分析熱應(yīng)力,從而預(yù)估其壽命[3]。在計算域的離散方面,結(jié)構(gòu)化網(wǎng)格以其計算資源占用小、收斂速度快、求解精度高、適用于大規(guī)模并行計算等優(yōu)點,在計算流體力學(xué)CFD(Computational Fluid Dynamics)領(lǐng)域具有不可替代的地位[4]。然而,由于渦輪葉型復(fù)雜的幾何結(jié)構(gòu),其計算域結(jié)構(gòu)化網(wǎng)格生成一直是數(shù)值模擬中的難點。在傳統(tǒng)的網(wǎng)格生成過程中,通常采用現(xiàn)有商業(yè)軟件生成非結(jié)構(gòu)化網(wǎng)格或非結(jié)構(gòu)/結(jié)構(gòu)化混合網(wǎng)格,如Gambit、ICEM、HyperMesh等。即使采用全結(jié)構(gòu)化網(wǎng)格,網(wǎng)格生成周期也十分漫長,而且能否生成結(jié)構(gòu)化網(wǎng)格十分依賴工程技術(shù)人員經(jīng)驗。
為提高渦輪葉片網(wǎng)格結(jié)構(gòu)化網(wǎng)格生成效率與自動化程度,縮短設(shè)計周期,本文在采用分塊拼接網(wǎng)格策略[5-6]與遺傳算法[7]的基礎(chǔ)上,提出了一種針對渦輪葉片表面氣膜孔的結(jié)構(gòu)化網(wǎng)格生成方案。首先建立了分塊質(zhì)量評價函數(shù),進而采用遺傳算法尋求分區(qū)的最優(yōu)方案,最終生成結(jié)構(gòu)化網(wǎng)格。
目前,對于網(wǎng)格分區(qū)的質(zhì)量評價大多依賴經(jīng)驗,常見的經(jīng)驗分區(qū)有針對圓孔特征的“O型分割”、針對三角特征的“Y型分割”(如圖1所示)等。
(a) O型分割 (b) Y型分割 圖1 常見結(jié)構(gòu)化網(wǎng)格分區(qū)方案
分區(qū)結(jié)果的質(zhì)量評價函數(shù)是分區(qū)過程自動化的基礎(chǔ)。分塊質(zhì)量的好壞體現(xiàn)在對其生成的網(wǎng)格質(zhì)量的影響上,而網(wǎng)格的質(zhì)量體現(xiàn)在對CFD計算結(jié)果的影響上。通過理論分析及經(jīng)驗表明,CFD對網(wǎng)格的一般性要求為雅可比矩陣行列式的值大于0,其他指標(biāo)包括網(wǎng)格扭曲、偏斜、翹曲等。本文在總結(jié)和借鑒眾多現(xiàn)有商業(yè)網(wǎng)格生成軟件網(wǎng)格質(zhì)量評價模塊的基礎(chǔ)上,提出了拼接多塊的質(zhì)量評價指標(biāo)。顯然,當(dāng)子分區(qū)為立方體時,所生成的網(wǎng)格質(zhì)量最優(yōu),則該分區(qū)質(zhì)量最好。所以,分塊質(zhì)量評價的根本目的就是確定分區(qū)與立方體的相似程度。依據(jù)評價指標(biāo)計算是在表面面上進行還是在立方體上進行,將評價指標(biāo)分為“2D評價指標(biāo)”與“3D評價指標(biāo)”。
1.1.1 2D評價指標(biāo)
2D評價指標(biāo)將在分區(qū)的各個表面上進行,主要計算指標(biāo)有:翹曲角、偏斜度、扭曲度。
偏斜度(s)表征了該分區(qū)表面的尖銳程度,用四邊形表面對邊中點連線的夾角與直角的偏差來表示,如圖2所示。歸一化后為:
(1)
圖2 偏斜度示意圖
翹曲角(w)定義為:四邊形表面對角線分割成的兩個三角形的夾角的最大值,如圖3所示。歸一化后為:
(2)
圖3 翹曲角示意圖
扭曲度(t)表征了分區(qū)表面的扭曲程度,用四邊形最大內(nèi)角與最小內(nèi)角中與直角偏差較大的角表示,如圖4所示。歸一化后為:
(3)
圖4 扭曲角示意圖
1.1.2 3D評價指標(biāo)
面積比(a)表征了立方體與四面體的相似程度用相對表面面積比的最小值表示。a越小,該分區(qū)質(zhì)量越差。
分區(qū)的雅可比矩陣行列式的值(J)反映了該分區(qū)與立方體的相似程度。如圖5所示,根據(jù)文獻[8], 假設(shè)Xm∈R3為六面體分區(qū)中的第m個頂點,Xm=[xm,ym,zm]T,m= 1,2,…,8。Xm,i∈R3(i=1,2,3)為與Xm共邊的頂點,em,i=Xm,i-Xm(i=1,2,3)為第m個頂點的邊向量的單位向量,則該分區(qū)某一頂點的雅可比矩陣Jm=[em,1,em,2,em,3]。
J=min[det(Jm)]m=1,…,8
(4)
圖5 雅可比矩陣列向量示意圖
式(4)中:J∈[-1,1],J越大,則表明該分區(qū)與立方體相似程度越大。當(dāng)J≤0,則認為該分區(qū)不成立。
針對氣膜孔網(wǎng)格生成的分區(qū)問題,采用遺傳算法來對分區(qū)結(jié)果尋優(yōu)。得到分區(qū)結(jié)果后,采用超限插值算法生成結(jié)構(gòu)化網(wǎng)格。根據(jù)文獻[7], 遺傳算法可以表示為:
SGA=(C,E,P0,M,Φ,Γ,Ψ,T)
(5)
式中:C為個體編碼方法;E為個體適應(yīng)度評價函數(shù);P0為初代種群;M為種群大?。沪禐檫x擇算子;Γ為交叉算子;Ψ為變異算子;T為遺傳運算終止條件。
遺傳算法在眾多領(lǐng)域應(yīng)用廣泛,關(guān)于遺傳算法的原理及流程本文不再贅述。需要指出的是,在染色體的編碼與解碼上,種群中每一個個體代表一種分區(qū)策略,用一個0到1之間的雙浮點精度數(shù)的一維數(shù)組R表示分割點所在的位置,R=[r1,r2,…,rn],rj(j=1,2,…,n)。R代表了該分割點在其所在邊上等分線段上的長度比例。以一個長度為32位的十進制整數(shù)的一維數(shù)組編碼串對R進行編碼。這樣,就可以用一個十進制數(shù)組的集合完成對分區(qū)策略的編碼。
例如:在長度為100的邊上,染色體的表現(xiàn)型為R=[0.5,0.5,0.5,0.5],則該邊分割形式如圖6所示。
圖6 染色體表現(xiàn)型示意圖
針對自動分區(qū)問題,個體適應(yīng)度函數(shù)取分區(qū)子塊質(zhì)量的最小值。這樣,能夠?qū)ι鲜鲈u價指標(biāo)逐一優(yōu)化。
E=min[f(Bi)]
(6)
式中:
f(Bi)=min[w(Bi),s(Bi),t(Bi),a(Bi),J(Bi)]
i=1,2,…,nf(Bi)∈[-1,1]
(7)
終止條件設(shè)定為:當(dāng)種群最大適應(yīng)度與平均適應(yīng)度之差小于10-5且各個分區(qū)的雅可比大于0時,終止計算。
終止條件可以表示為:
T=(Emax-Eavg)≤10-5∩J(Bi)>0
(8)
交叉概率與變異概率分別設(shè)置為0.75和0.1。
通過計算渦輪葉片氣膜孔的兩種布置方案來驗證上述分區(qū)方法的可行性。
平行氣膜孔簡化幾何模型如圖7所示,兩直徑為10的平行氣膜孔布置在200×100×100的平板上,各邊上的數(shù)字代表了邊上的上下兩表面對應(yīng)邊上的分割點的數(shù)目,兩個圓上的分割點數(shù)均為4。
圖7 兩平行氣膜孔幾何模型
采用上述分區(qū)方法,種群個體數(shù)目為300,種群最大適應(yīng)度收斂于0.439,遺傳算法收斂曲線如圖8所示。得到分區(qū)結(jié)果如圖9所示。根據(jù)該結(jié)果,采用超限插值算法[9]生成結(jié)構(gòu)化網(wǎng)格,網(wǎng)格尺寸設(shè)置為3.0,最終生成的網(wǎng)格結(jié)果如圖10所示,網(wǎng)格數(shù)量為39 520。圖11為z=50 mm平面上網(wǎng)格結(jié)果。
圖8 平行氣膜孔算例適應(yīng)度收斂曲線
圖9 平行氣膜孔分區(qū)結(jié)果
圖10 平行氣膜孔結(jié)構(gòu)化網(wǎng)格
圖11 z=50中心截面網(wǎng)格
將該網(wǎng)格結(jié)果導(dǎo)入ICEM進行網(wǎng)格質(zhì)量評價。網(wǎng)格質(zhì)量如表1所示。該分區(qū)結(jié)果表明:根據(jù)遺傳算法尋優(yōu)得到的分區(qū)結(jié)果生成的網(wǎng)格滿足CFD計算需求。
表1 平行孔算例網(wǎng)格質(zhì)量
交叉氣膜孔算例簡化幾何模型如圖12所示。兩直徑為10的氣膜孔交叉布置在200×200×100的平板上,氣膜孔法向分別為(1,0,1)、(0,1,1)。
圖12 交錯氣膜孔簡化模型
采用上述遺傳算法對該模型進行分區(qū),各邊分割點分布與平行氣膜孔算例相同。種群數(shù)量為500,遺傳算法進化374代,適應(yīng)度收斂于0.209。遺傳算法收斂曲線如圖13所示。得到的分區(qū)結(jié)果如圖14所示。同樣采用超限插值法生成結(jié)構(gòu)化網(wǎng)格,網(wǎng)格尺寸為3.0,網(wǎng)格數(shù)量為146 664,網(wǎng)格結(jié)果如圖15所示。
圖13 交錯氣膜孔算例收斂曲線
圖15 交錯氣膜孔網(wǎng)格
將該網(wǎng)格結(jié)果導(dǎo)入ICEM進行質(zhì)量評價。網(wǎng)格質(zhì)量見表2。該網(wǎng)格質(zhì)量滿足CFD計算需求。
表2 交錯孔算例網(wǎng)格質(zhì)量
最后,對渦輪葉片整體生成結(jié)構(gòu)化網(wǎng)格以驗證該方法的可行性。整體葉片網(wǎng)格與氣膜孔局部放大如圖16所示。
圖16 渦輪葉片整體網(wǎng)格
本文以渦輪葉片氣膜孔為結(jié)構(gòu)化網(wǎng)格生成的研究對象,采用多塊拼接網(wǎng)格,結(jié)合遺傳算法,并采用超限插值算法生成結(jié)構(gòu)化網(wǎng)格。通過對渦輪葉片整體生成結(jié)構(gòu)化網(wǎng)格以驗證該方法的可行性。結(jié)果表明:本文提出的基于遺傳算法的渦輪葉片氣膜孔自動分區(qū)方法可以得到滿足CFD計算需求的結(jié)構(gòu)化網(wǎng)格。本文工作總結(jié)如下:
(1) 經(jīng)檢驗,本文建立的分區(qū)質(zhì)量評價函數(shù)能夠反映分區(qū)的質(zhì)量。網(wǎng)格質(zhì)量與分區(qū)質(zhì)量相關(guān)度較高。
(2) 遺傳算法能夠有效處理結(jié)構(gòu)化網(wǎng)格分區(qū)問題。與傳統(tǒng)手工劃分網(wǎng)格方法相比,極大提高了網(wǎng)格生成的成功率與自動化程度。
(3) 在采用遺傳算法分區(qū)前,需要預(yù)先設(shè)定對應(yīng)面上的分割點數(shù)目,尚未解決全自動處理復(fù)雜幾何結(jié)構(gòu)分區(qū)的問題。