尚永鋒，馮天杰，付為剛

（中國民用航空飛行學(xué)院航空工程學(xué)院，四川 廣漢 618307）

1 引言

輪盤是壓氣機完成功能轉(zhuǎn)換的重要零件，在工作時承受著極大的負(fù)荷。輪盤轉(zhuǎn)速通常都達(dá)到每分鐘數(shù)千轉(zhuǎn)到數(shù)萬轉(zhuǎn)。輪盤除了固定葉片和承受自身的離心力外，還要承受葉片巨大的離心載荷。因為葉片和輪盤產(chǎn)生的離心載荷很大，所以輪盤需要設(shè)計較大的盤心承載，導(dǎo)致輪盤質(zhì)量增加，同時輪盤的體積也會隨之增大[1]。不僅使發(fā)動機氣動外形設(shè)計難度增加，更甚者輪盤發(fā)生破壞，產(chǎn)生的高能碎塊損壞機匣，影響飛行安全[2-5]。因此隨著優(yōu)化技術(shù)的發(fā)展，通過數(shù)值與仿真分析對發(fā)動機盤體尋求一種簡潔、高效的優(yōu)化方式來滿足輪盤輕質(zhì)化設(shè)計的需求也更具有研究意義和價值[6]。

對于輪盤結(jié)構(gòu)優(yōu)化設(shè)計的研究，大部分文獻(xiàn)都是以有限元分析為基礎(chǔ)，或開發(fā)獨立的優(yōu)化計算模塊，或借助現(xiàn)有的軟件優(yōu)化平臺。常見優(yōu)化設(shè)計所涉及的優(yōu)化算法包括了傳統(tǒng)的SUMT法、拉格朗日乘子法等，以及Kriging方法、遺傳算法等現(xiàn)代智能算法。文獻(xiàn)[7-8]使用ANSYS軟件對發(fā)動機高負(fù)荷渦輪盤進行了結(jié)構(gòu)優(yōu)化設(shè)計，將單幅板結(jié)構(gòu)改為雙幅板結(jié)構(gòu)來實現(xiàn)輪盤的減重。文獻(xiàn)[9]選取輪盤的結(jié)構(gòu)尺寸作為優(yōu)化變量，建立了輪盤結(jié)構(gòu)優(yōu)化模型，通過整體結(jié)構(gòu)優(yōu)化提高輪盤設(shè)計的材料利用率，可將輪盤質(zhì)量減輕10%。文獻(xiàn)[10]采用Isight優(yōu)化平臺，將輪盤質(zhì)量設(shè)定為目標(biāo)函數(shù)，以最大周向應(yīng)力和低循環(huán)疲勞壽命為約束，對輪盤結(jié)構(gòu)進行優(yōu)化，并通過對輪盤參數(shù)的研究，分析了各參數(shù)與目標(biāo)變量之間的關(guān)系。

對于輪盤結(jié)構(gòu)優(yōu)化設(shè)計的研究，大部分文獻(xiàn)都是以有限元分析為基礎(chǔ)，或開發(fā)獨立的優(yōu)化計算模塊，或借助現(xiàn)有的軟件優(yōu)化平臺。常見優(yōu)化設(shè)計所涉及的優(yōu)化算法包括了傳統(tǒng)的SUMT法、拉格朗日乘子法等，以及Kriging方法、遺傳算法等現(xiàn)代智能算法。文獻(xiàn)[7-8]使用ANSYS軟件對發(fā)動機高負(fù)荷渦輪盤進行了結(jié)構(gòu)優(yōu)化設(shè)計，將單幅板結(jié)構(gòu)改為雙幅板結(jié)構(gòu)來實現(xiàn)輪盤的減重。文章[9]選取輪盤的結(jié)構(gòu)尺寸作為優(yōu)化變量，建立了輪盤結(jié)構(gòu)優(yōu)化模型，通過整體結(jié)構(gòu)優(yōu)化提高輪盤設(shè)計的材料利用率，可將輪盤質(zhì)量減輕10%。文獻(xiàn)[10]采用Isight優(yōu)化平臺，將輪盤質(zhì)量設(shè)定為目標(biāo)函數(shù)，以最大周向應(yīng)力和低循環(huán)疲勞壽命為約束，對輪盤結(jié)構(gòu)進行優(yōu)化，并通過對輪盤參數(shù)的研究，分析了各參數(shù)與目標(biāo)變量之間的關(guān)系。

利用ANSYS優(yōu)化平臺，以輪盤的體積作為優(yōu)化目標(biāo)，以某型高負(fù)荷壓氣機輪盤為例，考慮輪盤所承受的載荷過大，按照輪盤設(shè)計準(zhǔn)則所規(guī)定的關(guān)鍵應(yīng)力約束作為優(yōu)化限制，確保輪盤與葉片連接處的尺寸在優(yōu)化中保證足夠的應(yīng)力強度，并依照等強度理論對輪盤的幾何尺寸進行約束，提出并建立輪盤結(jié)構(gòu)優(yōu)化設(shè)計數(shù)學(xué)模型及方法，對輪盤優(yōu)化前和優(yōu)化后進行對比分析，以說明壓氣機輪盤結(jié)構(gòu)優(yōu)化設(shè)計的關(guān)鍵和使用ANSYS一階優(yōu)化方法對輪盤優(yōu)化的優(yōu)勢。

2 輪盤簡化設(shè)計方法

2.1 基于等強度理論輪盤優(yōu)化設(shè)計

為了在簡化設(shè)計中確保輪盤整體結(jié)構(gòu)強度滿足規(guī)定的強度要求，可以采用等強度理論對輪盤進行設(shè)計。等強度理論是指讓輪盤各位置的徑向和周向應(yīng)力都等于某一常數(shù)，這樣設(shè)計出的輪盤質(zhì)量最輕[11]。航空發(fā)動機壓氣機盤近似工作在等溫的情況下，滿足等強度盤只有等溫盤的情況下才能實現(xiàn)等強度設(shè)計的要求。

2.2 軸對稱假設(shè)

雖然輪盤結(jié)構(gòu)十分復(fù)雜，但其幾何形狀通常都是對稱的。因此在實際設(shè)計中對輪盤提出了軸對稱假設(shè)，既包括幾何形狀，載荷分布以及溫度場和約束條件情況等都對稱于某軸[12]。

同樣輪盤中所有的應(yīng)力分量、應(yīng)變分量和位移分量也都對稱于這一軸。這種問題稱為空間軸對稱問題。利用軸對稱假設(shè)便可以極大的簡化輪盤設(shè)計。由于盤外緣的應(yīng)力很復(fù)雜，我們只研究盤本身的應(yīng)力，輪盤簡化結(jié)構(gòu)，如圖1所示。這樣輪盤就可以看作是由一個截面如圖2所示繞指定軸旋轉(zhuǎn)而成。那么問題就可以簡化成對一個截面的分析[13]。

圖1 壓氣機輪盤模型Fig.1 Compressor Disk Model

圖2 輪盤區(qū)域分布Fig.2 Disc Area Distribution

3 基于ANSYS輪盤結(jié)構(gòu)優(yōu)化設(shè)計

3.1 輪盤結(jié)構(gòu)優(yōu)化

以EJ200 發(fā)動機壓氣機輪盤工作溫度為213℃為例。其主要結(jié)構(gòu)參數(shù)為：盤外緣半徑為228mm；盤內(nèi)徑為50mm；輪緣寬度為60mm；輪緣厚度初始設(shè)為H=11mm輪盤轉(zhuǎn)速n=15000r/min；盤外緣分布面載σa=12.83MPa；盤材料用TC11鈦合；密度ρ=4480kg/m3；可查出輪盤在T=213℃下的材料參數(shù)：輪盤彈性模量E=105GPa；泊松比ν=0.3；許用應(yīng)力σ=570MPa。為了防止輪盤有害變形和破裂。輪盤需滿足子午面平均周向應(yīng)力和圓柱面最大平均徑向應(yīng)力的約束。

輪盤可分為三個區(qū)域，如圖2所示。即A區(qū)、B區(qū)和C區(qū)，其中A區(qū)是為了消除等強度盤盤心處的應(yīng)力集中所設(shè)計的盤體；B區(qū)即為等強度盤體，是根據(jù)等強度理論設(shè)計的；C區(qū)是用來連接葉片的，由于工作中輪盤和葉片都高速運動，連接處需要承受巨大載荷在優(yōu)化調(diào)整中一般不進行改動。所以優(yōu)化輪盤時只對A區(qū)和B區(qū)的尺寸進行優(yōu)化[14]。輪盤的初始參數(shù)，如表1所示。輪緣處的倒角半徑初步設(shè)定為R=8mm。根據(jù)輪盤初始結(jié)構(gòu)在ANSYS中通過定義變量、建立模型，分化網(wǎng)格，定義邊界條件，計算求解等步驟。對初始輪盤進行強度分析得到其等效應(yīng)力分布，如圖3所示。

表1 輪盤初始設(shè)計參數(shù)Tab.1 Initial Design Parameters of Disc

圖3 初始輪盤等效應(yīng)力云圖Fig.3 Cloud Chart of Equivalent Stress of Initial Disc

3.2 優(yōu)化過程

首先根據(jù)輪盤初始結(jié)構(gòu)在ANSYS中通過定義變量、建立模型，分化網(wǎng)格，定義邊界條件。通過限制所有節(jié)點等效應(yīng)力STR不超過570MPa，將體積V設(shè)定為目標(biāo)函數(shù)，使用一階優(yōu)化方法進行優(yōu)化。一階優(yōu)化方法通過因變量的一階偏導(dǎo)數(shù)確定在設(shè)計空間內(nèi)的搜索方向[15]。

在ANSYS 中讀取所分析的數(shù)據(jù)，輪盤初始體積為V=0.0025344m3。從主菜單中選擇Main menu：General Postproc＞List Results＞Sorted Listing＞Sort Nodes，在打開的對話框中選中圖4-3所示的設(shè)置，讀取等效應(yīng)力，并在命令輸入框中輸入“*GET，STR，SORT，MAX”，讀取最大等效應(yīng)力。最后所讀取的參數(shù)都在參數(shù)欄里，如圖4所示。

圖4 優(yōu)化參數(shù)Fig.4 Optimization Parameters

然后將此時的數(shù)據(jù)保存，打開Design Opts寫入剛才保存的DB文件。讀入后開始對設(shè)計變量進行設(shè)置。并對H1，H2，H3，H4，和R設(shè)定變化區(qū)間。其次再對狀態(tài)變量進行定義，將STR定義在（10～570000000）之間，容差為10。然后對目標(biāo)函數(shù)V 進行設(shè)置，其收斂判據(jù)為0.00001。最后從主菜單Design opt＞Method/tool，選擇First-order，設(shè)置最大迭代次數(shù)為20。然后在主菜單中Main menu：Design opt＞Run 進行優(yōu)化求解。經(jīng)過九次迭代達(dá)到最優(yōu)解。從主菜單中選擇Main menu：Design opt＞Design Sets＞Graphs/Tables.在X 軸列表框選擇Set numbers，在Y 軸列表框中選擇V，單擊Ok，圖像顯示優(yōu)化變量的變化過程，如圖5所示。

圖5 優(yōu)化變量的變化過程Fig.5 Change Process of Optimization Variables

由圖5可以看出在迭代計算的第9次運算，就中斷了優(yōu)化，就表示在第九次的時候已經(jīng)達(dá)到最優(yōu)解了，此時的von mises stress圖，如圖6所示。這樣我們從輪盤優(yōu)化結(jié)果中讀出一組數(shù)據(jù)，如圖7 所示。圖中第九組數(shù)據(jù)可以讀出在滿足最大應(yīng)力562MPa＜570MPa，滿足材料限制，其優(yōu)化后的體積V2為0.0018250m3。

圖6 優(yōu)化后等效應(yīng)力圖Fig.6 Optimized Equivalent Stress Diagram

圖7 輪盤優(yōu)化結(jié)果Fig.7 Optimization Results of Disc

3.3 優(yōu)化結(jié)果分析

優(yōu)化后的周向應(yīng)力云圖，如圖8所示。

圖8 優(yōu)化后的周向應(yīng)力云圖Fig.8 Cloud Chart of Optimized Circumferential Stress

由圖6與圖8可以看出，優(yōu)化后的輪盤有效的改善應(yīng)力集中問題，并用公式（1）將優(yōu)化前與優(yōu)化后輪盤體積進行計算，得到輪盤體積減少了27.99%。

式中：V1，V2—優(yōu)化前后輪盤體積，單位m3。

4 結(jié)論

（1）運用ANSYS對輪盤進行了優(yōu)化設(shè)計的結(jié)果顯示，基于等強度理論優(yōu)化設(shè)計輪盤，應(yīng)力集中現(xiàn)象減少，應(yīng)力分布優(yōu)于傳統(tǒng)設(shè)計。（2）確保輪盤應(yīng)力強度的同時，通過結(jié)構(gòu)優(yōu)化使得輪盤整體體積減少了27.99%，提高輪盤材料利用率的同時，減輕發(fā)動機重量，對提高發(fā)動機推質(zhì)比有極大的參考意義。（3）一階優(yōu)化輪盤應(yīng)力強度由低到高，慢慢貼近約束強度，輪盤體積慢慢減小，顯示了一階優(yōu)化的穩(wěn)定與可靠性。