航空發動機葉片專用制造裝備設計研究*

劉翠煜，鐘建琳

(北京信息科技大學 機電工程學院，北京 100192)

?

航空發動機葉片專用制造裝備設計研究*

劉翠煜，鐘建琳

(北京信息科技大學 機電工程學院，北京 100192)

航天技術的競爭愈演愈烈，制造業作為國民經濟的基礎力量，是國際間競爭的關鍵性因素。葉片是飛行器發動機的重要組成部分，需具有高的精度和加工質量。文章以五軸聯動葉片加工中心為例，對其結構構件進行了動態特性分析，并對薄弱環節進行了優化設計。最后通過優化前后結果對比分析，驗證了優化的合理性和可行性，為今后的設計提供了一定的參考。

葉片；加工中心；動態特性；優化設計

0　引言

科學技術飛速發展，航空制造領域的先進性成為了影響國家國際地位的又一重要標志。航空發動機的葉片屬于高、精、尖產品，多年來高精度和高表面質量的自由曲面零件的加工一直是制造業的難點。目前，國內、外對葉片制造進行了大量的研究，主要包括以下三個方面：①運用計算機分析手段對葉片進行多學科優化設計[1-4]；②優化葉片生產過程的加工工藝，選取適當加工參數和加工方法；③提高葉片制造設備的性能，保證其加工質量。其中，提高制造裝備研究的主要手段之一就是機床動態特性研究。本文即從改善制造裝備入手，通過應用計算機輔助軟件對某單位自主研發的航空發動機葉片專用生產設備的主要結構進行分析和優化來提高機床的動態特性，改善葉片加工質量。驗證合理的優化設計，不但提高了本臺機床的性能，為今后機床結構的設計尤其是主要受力部位設計也提供了一定的參考。

1　機床設計

研究發現，產品的制造成本超過 70%是在設計階段被決定的，但設計階段的成本僅占整個開發階段的6%，由此可見設計階段對產品的發展十分重要[5]。機床是制造業的基本裝備，其設計可以劃分為產品規劃、方案設計、技術設計和工藝設計四個階段。傳統的設計方法在功能、周期及成本方面已經無法滿足新的需求。現代設計理論采用先進的CAD/CAE技術，結合有限元對產品進行分析、優化，省去了產品樣機反復試制、測試及改進過程，提高產品的動態性能的同時，有效縮短了設計周期，實現了企業低成本快速響應市場、產品批量化生產及大規模生產條件下的個性化需求[6]。

2　葉片加工中心設計

葉片是發動機結構及功能實現的主要組成部件，其形狀、強度以及抗熱疲勞等性能決定了發動機的工作效率和使用壽命。五軸聯動加工中心是目前更高效、高精的數控加工設備，被廣泛應用于復雜曲面的高效、精密、自動化加工，特別是在解決葉輪、葉片等方面具有獨特的優勢[7]。

2.1葉片加工特點及設計要求

應用環境的特殊需求使得發動機葉片一般設計的比較薄，本文加工中心生產的葉片，長約26.5mm，寬約28.8mm，前后緣轉接R0.1～0.3，輪廓度R0.05～0.08，結構如圖1所示。葉片的表面為自由曲面，采用數控銑削方法加工時容易產生較大的加工變形[8]。為滿足發動機高耐熱性、高強度的使用要求，本文加工的葉片即采用性能較好的高溫鎳基合金材料GH4169，其加工需要的切削力大，導致切削溫度較高，易造成切屑與刀具粘連以及機械的磨損，因此，對加工機床的性能提出了更高的要求。

圖1　葉片

2.2葉片加工中心尺寸及參數設計

根據本文葉片的結構特點，國內某單位研制開發了葉片專用五軸聯動加工中心如圖2所示。該加工中心能夠實現X、Y、Z軸及 A、C 軸聯動，X、Y、Z 軸直線移動距離分別為300mm、250mm、300mm，主軸最高轉速可達24000r/min，移動部件盡量采用方形框架結構，整體采用對稱結構，穩定性較高。

圖2　葉片加工中心示意圖

3　葉片加工中心結構分析

在高速切削加工過程中，機床各結構部件的動態性能將直接影響機床的精度、抗振性能、切削效率以及壽命[9]，從而影響到零件表面成形運動軌跡的準確性。因此，對加工中心重要構件的動、靜態特性進行了分析。

3.1床身動態特性分析及優化

3.1.1床身模型簡化與處理

使用三維軟件進行模型簡化，例如刪除小孔、倒角和圓角等特征，對小曲率、小錐度的直線和平面進行處理等。該床身與地面間通過M16的地腳螺栓連接，采用線性彈簧-阻尼單元法進行建模。根據張學良[10]建立的數學模型對螺栓結合面剛度特性進行識別、計算，將結合面特性簡化為一組空間彈簧-阻尼單元。文獻[11]中提出結合面剛度參數識別經驗公式為：

k=αpβΔS

(1)

式中，k—法向剛度kn或切向剛度kτ，單位為N/μm；P—壓力，單位為MPa；α、β—與結合面的材料、加工方式、表面粗糙度和潤滑狀況等有關因素；ΔS—實際接觸面積，單位為m2。根據文獻[12]中數據，法向：α=806256，β=0.8 ；切向： α=132786 ，β=0.37 。計算彈簧剛度為：kn=5.3×107N/mm，kτ=6.0×105N/mm。

3.1.2床身動態特性及靜力分析

床身前3階模態固有頻率與振型描述如表1所述，各階振型如圖3所示。靜力分析結果如圖4所示，中間部位凹陷明顯，最大變形量為9.34 μm。

表1　床身前3階固有頻率及振型

圖3　前三階模態振型圖

圖4　靜力分析結果

3.1.3床身模型優化及改進結果

分析仿真結果可知，在低階時，床身出現整體移動、變形，靜力分析顯示變形部位較為集中、變形量大，且為實際加工中的重要承載位置。因此，需對模型進行優化以提高其穩定性。為盡量避免增加工藝的復雜性，采取的主要措施為加大床身底面的支撐面積，同時針對靜力變形較大的主要承載部位，對其床身底部對應位置的回型筋板進行填充和加固處理。優化前后模型對比如圖5。

圖5　優化前后模型對比

優化后的模型模態分析結果與優化前對比如表2所示，各階固有頻率均有較大提高，且相近固有頻率下床身的變形情況改善明顯，振型圖如圖6所示。靜力分析結果如圖7所示，最大變形量為5.74 μm，與優化前相比，變形量減少約4 μm。

表2　床身模型優化前、后前三階結果對比

圖6　模型優化后前三階模態振型圖

圖7　模型優化后的靜力分析結果

3.2X軸導向滑座動態特性分析及優化

3.2.1X軸導向滑座模型簡化與處理

模型簡化與處理方法參照3.1.1，X軸導向滑座通過M24螺栓與床身連接，計算等效彈簧剛度kn=4.14×1010N/mm，kτ=5.36×108N/mm。

3.2.2X軸導向滑座動態特性及靜力分析

結合X軸導向滑座在機床整體布局中的受力特點，重點分析其在X方向的振型及受力變形情況。前3階模態固有頻率與振型描述如表3所述，各階振型如圖8所示。靜力分析結果如圖9所示，最大變形量為11.6 μm，前側主要受力部位變形明顯。

表3　X軸導向滑座前3階固有頻率及振型

圖8　模型優化后前三階模態振型圖

圖9　靜力分析結果

3.2.3模型優化及改進結果

分析結果可知，低階X軸導向滑座前側變形量較大，會嚴重影響刀具的位置精度。觀察其機構可知，X軸滑座的前側中間位置設計了一段凹槽，如圖9中棕色變形位置中間。結合加工中心整體布局及各部件運動要求，發現此段凹槽并無特殊工藝要求，且考慮到處為承受X軸滑座及主軸重量的主要部位，應具有足夠的強度，盡量減少變形。因此，設計優化時對凹槽位置進行了填補，模型對比如圖10。

圖10　優化前后模型對比

對優化后的模型進行了分析，其模態分析結果與優化前對比如表4所示，各階固有頻率均有較大提高，且相近固有頻率下床身的變形情況改善明顯，振型圖如圖11所示。靜力分析結果如圖12所示，最大變形量為6.9 μm，與優化前相比，減少約4.7μm。

表4　X軸導向滑座模型優化前、后前三階結果對比

圖11　模型優化后前三階模態振型圖

圖12　模型優化后的靜力分析結果

綜上所述，對床身和X軸導向滑座優化后，其固有頻率明顯提高，靜力變形量有所減少。因此，該實例的優化方案合理可行，達到了改善加工中心剛度、穩定性的效果，有效提升了機床整體的加工性能。

4　結束語

航空發動機應用環境的特殊性決定了航空零件制造以及生產設備的特殊性。相比于一般的零件制造，面向航空應用的零件，其制造精度要求更嚴格，工藝過程更為復雜，制造設備的剛度及穩定性要求更高。本文以葉片五軸聯動葉片加工中心為例，對其主要結構進行動態特性及靜力分析，并對薄弱環節進行了優化。優化前、后結果對比，驗證了優化設計的合理性，為之后的設計提供了一定的參考。

[1]MartinBuβmann,JürgenKraus,ErwinBayer.AnIntegratedCost-EffectiveApproachtoBliskManufacturing[C]．InternationalSymposiumonAirBreathingEngines17th.MunichGermany:AmericanInstituteofAeronauticsandAstronautics, 2005:1100-1108.

[2] 于紅英,唐德威,傘紅軍.汽輪機葉片參數化設計關鍵技術研究[J].計算機集成制造系統,2006,12(10):1537-1542.

[3]IdahosaU,GolubevVladimirV,BalabanovVladimirO.AnautomatedoptimaldesignofafanbladeusinganintegratedCFD/MDOcomputere-nvironment[J].Engineeringapplicationofomputationalfluidmechanics,2008,2(2):141-154.

[4]BeirowBernd,KiihhomArnold,SchrapesSven.Adiscretemodeltoconsidertheinfluenceoftheairflowonbladevibrationsofanintergralbliskcompressorrotor[C].Berlin:ProceedingsoftheASMETurboExpo, 2008:38-392.

[5] 馬超民. 產品設計評價方法研究[D]. 長沙:湖南大學,2007.

[6] 張興朝. 基于有限元分析的模塊化數控機床結構動態設計研究[D].天津:天津大學,2001.

[7] 陸啟建, 楮輝生. 高速切削與五軸聯動加工技術[M]. 北京:機械工業出版社, 2003.

[8] 劉維偉, 李杰光, 趙明．航空發動機薄壁葉片加工變形誤差補償技術研究[J].機械設計與制造,2009,(10):175-177.

[9] 徐開元, 徐武彬, 張宏獻, 等.VMC-1000型立式加工中心立柱結構分析與動態設計[J]. 組合機床與自動化加工技術, 2010(4):11-17.

[10] 張學良. 機械結合面動態特性及應用[M]. 北京:中國科學技術出版社, 2002.

[11] 穆希望, 米潔. 基于關鍵螺栓結合面的立主軸模態分析[J]. 制造技術與機床, 2014(8):78-82.

[12] 翁德凱, 程寓, 李奎, 等. 基于結合面參數的機床整機有限元建模與分析[J]. 組合機床與自動化加工技術, 2012(3):29-33.

(編輯李秀敏)

AeroEngineBladeEedicatedManufacturingEquipmentDesignResearch

LIUCui-yu,ZHONGJian-lin

(SchoolofMechanicalandElectricalEngineering,BeijingInformationScienceandTechnologyUniversity,Beijing100192,China)

Thecompetitioninaerospacetechnologyisintensifiedandmanufacturingstrengthasthebasisofthenationaleconomyisakeyfactorintheinternationalcompetition.Bladesasanimportantpartoftheaircraftengineshouldhaveahighprecisionandprocessingquality.Inthispaper,takingthefive-axisbladesmachiningcenterasanexample,thedynamiccharacteristicsofitsStructuralmembersareanalyzedandweaklinksareoptimized.Finally,theresultsbeforeandoptimizedarecomparedtoverifytherationalityandfeasibilityoftheoptimizationandthisconclusionprovidesacertainreferencetothefutureoptimizeddesign.

blades;machiningcenter;dynamiccharacteristics;optimizationdesign

1001-2265(2016)08-0137-03DOI:10.13462/j.cnki.mmtamt.2016.08.037

2015-09-16

國家重大專項資助項目(2013ZX0400-1061)

劉翠煜(1989—)，女，石家莊人，北京信息科技大學碩士研究生，研究方向為機械制造及自動化，(E-mail)liucuiyu0311@163.com。

TH113；TG506

A