渦輪葉片氣膜孔結(jié)構(gòu)化網(wǎng)格生成

張 祥 劉 勇 趙 璐 李亞清

(南京航空航天大學(xué)能源與動力學(xué)院航空發(fā)動機熱環(huán)境與熱結(jié)構(gòu)工業(yè)和信息化部重點實驗室 江蘇 南京 210016)

0 引 言

渦輪是發(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)格。

1 分塊質(zhì)量評價函數(shù)與遺傳算法

1.1 分塊質(zhì)量評價函數(shù)

目前，對于網(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ū)不成立。

1.2 氣膜孔自動分區(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。

2 算例驗證

通過計算渦輪葉片氣膜孔的兩種布置方案來驗證上述分區(qū)方法的可行性。

2.1 平行氣膜孔算例

平行氣膜孔簡化幾何模型如圖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ì)量

2.2 交錯氣膜孔算例

交叉氣膜孔算例簡化幾何模型如圖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)格

3 結(jié) 語

本文以渦輪葉片氣膜孔為結(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ū)的問題。