單晶高溫合金行業現狀及其數值模擬應用綜述

梁琳++李媛智

摘要： 了解單晶高溫合金的行業現狀及其數值模擬應用，對于了解產品的市場動向以及提高產品的質量和研發效率具有現實意義.概述國內外單晶高溫合金的行業現狀及其數值模擬應用，重點總結我國單晶高溫合金的研發狀況、數值模擬技術在航空發動機用單晶高溫合金制備工藝、性能評價等方面的應用，并提出改善我國單晶高溫合金生長工藝與質量控制的努力方向.

關鍵詞： 航空發動機葉片； 單晶高溫合金； 制備工藝； 數值模擬

中圖分類號： V263文獻標志碼： A

0引言

隨著科學技術與社會經濟的發展，航空航天行業呈現產品軍民多樣化、市場競爭日益激烈的態勢.渦輪發動機葉片的質量是衡量航空航天產品性能的重要指標.單晶高溫合金以其優異的高溫抗疲勞、抗蠕變性能，已成為制造渦輪發動機葉片的重要材料.單晶高溫合金的生產制備工藝、缺陷預測與控制是影響其工程實際應用的關鍵因素.傳統一般采用實驗方法研究上述問題，但實驗往往成本高、周期長.隨著計算機軟、硬件技術的不斷進步以及計算數學、計算力學等學科的交叉滲透，數值模擬技術作為實驗法的輔助手段，已日益在渦輪發動機葉片等高端裝備制造領域中得到廣泛的應用.因此，了解單晶高溫合金的行業現狀及其數值模擬應用，對于了解產品的市場動向以及提高產品的質量與研發效率具有現實意義.

本文從國內外單晶高溫合金的行業現狀出發，結合單晶高溫合金的制備工藝和性能評價等問題，介紹數值模擬技術在航空渦輪發動機葉片單晶高溫合金中的應用，并提出改善發動機葉片單晶高溫合金生長工藝與質量的數值模擬研究方向.

1國外單晶高溫合金行業概況

目前，全球高溫合金材料被廣泛應用于各個工程領域，但總體而言，核心制造技術主要被為數不多的企業所占據.

全球范圍內能夠生產航空航天用高溫合金的企業不超過50家，主要集中在美、英、日等國.[1]發達國家將航空航天領域的高溫合金產品視為其長期發展的戰略物資，嚴格控制對外出口.美國航空航天單晶高溫合金的年產量約為5萬t，其中約60%用于民用工業，其研發工藝和應用技術一直處于國際領先水平.主要高溫合金企業的營業收入和營業利潤均在波動中上升.在美國，生產高溫合金水平較高的公司主要有GE，PrattWhitney，CannonMuskegon，Inco和Carpenter，其中GE公司能夠自主研發與生產航空發動機用的高溫合金[1].

英國是歐盟國家中高溫合金研發和生產的主要代表.英國研發高溫合金起步較早，其鑄造合金技術非常先進，具有代表性的是Mond Nickel公司生產的Nimocast合金.后來該國RollsRoyceplc公司研制SRR2000和SRR2060等航空發動機用定向凝固單晶合金.該公司的主營業務收入和利潤一直處于提升的狀態.

日本近年來的單晶高溫合金工藝發展迅速，在鎳基單晶高溫合金、鎳基超塑性高溫合金等方面取得較大成果，已成功開發出在1 200 ℃高溫下能保持足夠強度的新型耐高溫合金.日本高溫合金生產企業主要有IHI，JFE，新日鐵和神戶制鋼等公司.

總體而言，發達國家美、英、日等國在研發與制備航空發動機用高溫合金方面具有雄厚的技術實力，其核心技術領跑全球，值得國內單晶高溫合金研發與生產人員參考借鑒.

2國內單晶高溫合金行業現狀

20世紀80年代初，中航工業北京航空材料研究院在國內率先開展單晶合金及葉片技術的研究，成功研制出我國第一代單晶高溫合金DD3.20世紀90年代又成功研制出綜合性能優異的第二代單晶高溫合金DD6.DD6拉伸、持久、蠕變、疲勞、抗氧化及耐熱耐腐蝕性能等均達到國外廣泛應用的第二代單晶合金的性能水平.自成功研制出我國首件單晶空心渦輪葉片以來，中航工業北京航空材料研究院為多種型號先進航空發動機提供數萬件單晶葉片，某些裝配單晶渦輪葉片的發動機已隨機翱翔于藍天.

國內具有代表性的高溫合金研發和生產企業見表1.

表 1國內高溫合金研發和生產的主要企業

Tab.1Superalloy research and production enterprises in China單位高溫合金業務簡介鋼研

高納擁有技術、資金、市場和規模等優勢；產品涵蓋鑄造、變形、新型高溫合金等；鑄造高溫合金的產能近2 000 t，變形高溫合金產能超過150 t，粉末高溫合金產能近100 t，ODS合金、核電等新型高溫合金產能近100 t撫順

特鋼國內老牌高溫合金生產企業，產品主要以變形高溫合金為主，規模國內最大；高溫合金產能近3 000 t，民用產品約占40%，高溫合金營業收入近6億元寶鋼

特鋼國內老牌的高溫合金生產企業，民用產品占比高，但也生產GH4169和GH738等航空航天用變形高溫合金；其高溫合金產量超過1 000 t長城

特鋼以一般的變形高溫合金產品為主齊齊哈

爾特鋼以生產一般的民用變形高溫合金產品為主北京航

材院技術實力出眾，主要從事飛機、發動機和直升機用材料、工藝、檢測技術研究，具有高性能材料小批量生產和高難度重要部件的研發能力，在高溫合金單晶葉片領域具備較強科研與生產實力沈陽金

屬所以國家重大攻關項目研究為主，涉及高溫合金等高端裝備材料研發等領域

國內高校及相關企業在耐高溫單晶合金方面也開展相關工作.

燕平等[2]對DD402（CMSX2）單晶合金標準熱處理狀態組織及850 ℃下500，1 500和3 000 h時效的組織進行顯微觀察和比較，測定長期時效后的高溫持久性能.DD402單晶合金屬于第一代單晶合金，已成功應用于某發動機的Ⅰ級渦輪葉片.為提高我國航空發動機的質量，鋼鐵研究總院和南方航空動力機械公司共同開展DD402單晶合金及其Ⅰ級單晶渦輪葉片的研制工作.endprint

北京航空材料研究所在單晶高溫合金的研制方面也開展很多工作.陳榮章[3]總結分析單晶高溫合金的發展現狀，認為從20世紀80年代初第一代單晶高溫合金研制成功以來，單晶合金的發展甚為迅速，第二代、第三代單晶合金相繼出現和應用，為航空發動機的性能提高作出重要貢獻.單晶高溫合金在先進渦輪發動機葉片中的應用日益擴大，但客觀而言，我國在單晶合金工藝研究方面仍落后于當前國際先進水平.張勇[4]介紹一種加礦化劑的單晶殼型.該殼型具有高強度、壁薄且均勻的特點，已成功應用于澆注單晶葉片.張兵等[5]基于近年來國內外對于定向凝固和單晶高溫合金再結晶行為的研究，總結定向凝固和單晶高溫合金再結晶的主要影響因素，包括熱處理溫度與時間、第二相粒子、變形溫度等.李影等[6]總結鎳基單晶高溫合金的反常屈服行為，并討論導致產生這種反常屈服行為的變形機制以及這些變形機制的發展過程.

中國科學院金屬研究所的劉春廷等[7]采用低壓氣相沉積法，在鎳基高溫合金DD32上制備鋁化物涂層.在900和1 000 ℃氧化500 h后，表面氧化膜為致密的αAl2O3和針狀的θAl2O3.氧化后的鋁化物涂層外層為βNi3Al，內層（擴散層）母體為Ni3Al.

朱鷗等[8]總結國內外航空發動機用單晶高溫合金的熱處理工藝，介紹國外第一代至第三代單晶合金成分與熱處理的關系以及最優的熱處理條件.研究結果表明：升溫和冷卻速率對于單晶高溫合金的使用性能影響顯著.

丁智平[9]以鎳基單晶合金渦輪葉片的壽命分析為工程背景，研究鎳基單晶合金的熱彈塑性行為和高溫低周疲勞性能，并利用有限元分析軟件預測某發動機渦輪盤鎳基單晶合金葉片的三維非線性循環應力應變以及葉片的低周疲勞壽命.趙萍[10]將晶體塑性理論和非線性隨動硬化規律應用于低周疲勞研究，基于國內第一代鎳基單晶高溫合金DD3，對其單軸、多軸疲勞開展研究，并建立預測航空發動機單晶葉片低周疲勞的本構模型.

孟杰等[11]研究熱處理溫度與時間、變形程度及合金成分等多種因素對鎳基單晶高溫合金再結晶的影響規律，分析再結晶對其蠕變和疲勞性能的影響.

卿華等[12]采用實驗研究與有限元數值模擬方法相結合的方法對比分析帶孔和不帶孔的某第二代鎳基單晶合金平板試樣的蠕變性能.研究結果表明：在高溫蠕變下，平板試樣的晶體取向與是否開孔對蠕變壽命有明顯的影響.

綜上所述，國內對單晶高溫合金的研究工作主要側重于單晶高溫合金蠕變、疲勞等力學性能及相關的熱處理工藝的研究，對于單晶高溫合金生長過程中的缺陷預測、控制與優化的研究鮮有報道.另外，對單晶高溫合金生長機理的研究幾乎空白.

3數值模擬技術在航空發動機單晶高溫合金葉片領域中的應用現狀用于制造航空發動機渦輪葉片單晶高溫合金的定向凝固鑄件生產工藝復雜、控制要求高，因而通過實驗研究單晶高溫合金葉片的成本較高且研發周期長.隨著計算機軟硬件技術的進步，數值模擬技術發展迅速，在工業領域已得到廣泛應用.利用數值模擬技術，可以有效模擬單晶高溫合金的定向凝固生長過程，預測、控制與優化其微觀組織與缺陷情況，提高單晶高溫合金質量、降低研發成本并縮短研發時間.

薛明等[13]模擬分析空心渦輪葉片在定向凝固過程中陶瓷芯內部的溫度分布，研究陶瓷芯定位及型殼熱物性參數匹配的影響.李嘉榮等[14]采用有限元模擬軟件計算單晶合金定向凝固過程中的溫度場.楊亮等[15]針對單晶高溫合金精鑄薄壁件制備困難的問題，采用數值模擬方法分析DD6單晶高溫合金的定向凝固生長過程.通過研究試件幾何形狀以及工藝參數對定向凝固生長過程中溫度梯度和糊狀區的影響，結果表明：幾何形狀對單晶高溫合金試樣的定向凝固生長有重要影響，提高澆注溫度或降低軸拉速度有助于提高試樣固液界面前沿液相溫度梯度或減小糊狀區的寬度.

于靖等[16]基于CA模型、枝晶生長機理和熱量、溶質傳輸方程，建立定向凝固單晶高溫合金葉片的三維數值分析模型，考慮多葉片之間以及加熱爐之間復雜的輻射傳熱，模擬不同抽拉速度下葉片內部的溫度分布，并采用分層算法模擬定向凝固單晶高溫合金葉片的微觀組織演變.PAN等[17]基于有限差分法數值模擬航空發動機用單晶的定向凝固生長.研究結果表明：通過適當改變回退率，可以既提高生產能力，又避免產生晶粒缺陷.張航等[18]采用元胞自動機以及有限差分法對DD6高溫合金三維定向凝固枝晶生長開展數值模擬研究.

在國外，單晶高溫合金研發與生產單位非常重視晶體生長過程的全程整體預測以及晶體生長缺陷的預測與優化.為適應科研與生產的實際需要，結合傳熱與流體計算力學，開展單晶高溫合金生長數值模擬分析軟件的研發工作，并在數值分析軟件的算法、預測精度與求解效率方面進行很多有益的探索研究，其中，DUPRET等[1920]開展諸多研究工作，并提出卓有成效的晶體生長全局數值模擬算法.基于這些穩健、高效的算法，DUPRET教授主持開發晶體生長數值模擬軟件FEMAG.該軟件可模擬提拉法、泡生法、區熔法、垂直布里奇曼法、垂直梯度凝固法、定向凝固法、熱交換法和物理氣相傳輸法等晶體生長工藝以及多物理場耦合仿真問題.其中，FEMAG軟件的定向凝固法生長工藝模擬功能頗具特色，通過集成定向凝固法生長工藝條件的設置，不僅可以自動捕捉晶體定向凝固過程中的固液界面形狀與位置，還可以計算熔體對流和輻射傳熱、預測加熱器的功率以及計算氧雜質濃度、控制晶體生長的缺陷.

利用FEMAG軟件可以有效模擬包括單晶高溫合金在內的晶體定向凝固過程以及預測與控制晶體生長的質量，其算法穩健，計算精度高，求解效率高.而相比而言，國內尚未研制出同等高效、高精度模擬預測晶體定向凝固法生長的數值分析軟件，這方面的工作尚需國內的研發人員加以重視，并努力提高晶體生長模擬的研發水平，加快自主研發步伐，以縮小與國外的差距.

4結束語

作為高端裝備制造領域航空發動機渦輪葉片的重要材料，單晶高溫合金的研發與生產應用已日益深入.從國內外單晶高溫合金的行業狀況看，單晶高溫合金的核心制造技術仍然為國外某些發達國家所掌握，研發與生產技術水平仍明顯優于國內同行的工作.從研究方法和狀況看，國內外對單晶高溫合金的力學性能（蠕變、疲勞等）研究較多，采用的手段多為有限元數值模擬方法，在工藝方面也主要側重于對定向凝固法生長單晶高溫合金的研究，對于有限元法以外的其他數值模擬方法、定向凝固法以外的其他生長工藝方法研究頗少.對于單晶高溫合金生長機理的理論研究則更少，且難有重大的理論進展突破.從單晶高溫合金生長模擬分析軟件的研制工作看，國內較國外而言，差距仍很顯著.endprint

為此，以下幾方面值得國內從事單晶高溫合金研發與生產的相關人員參考.

1）針對國內在單晶高溫合金理論、工藝以及數值模擬分析方面研發相對落后的局面，國內的相關工作者需要予以重視，并努力提高單晶高溫合金工藝及其實驗、數值預測的研發水平，加快自主產權軟件的研發，努力縮小與國外同行的差距.

2）國內關于單晶高溫合金的研發與生產多為獨立進行，多個單位與行業的合作較少，資源的配置利用較為分散.開展多方合作研發，集中多方的研發優勢，對于核心技術的攻關、資源的高效配置以及研發效率的提高具有重要的促進作用.

3）引進、吸收和消化國外先進的研發經驗與技術，開展擁有自主知識產權的晶體生長數值分析軟件的研發工作，對于提升我國在單晶高溫合金數值模擬分析方面的技術水平與核心競爭力具有現實意義.參考文獻：

[1]徐志國. 高端裝備制造關鍵材料遙看高溫合金兩千億市場空間[N]. 中國航空報，20141206（D04）.

[2]燕平， 單熙， 趙京晨， 等. DD402單晶合金及其Ⅰ級單晶渦輪葉片的組織穩定性[J]. 航空材料學報， 1999， 19（1）： 212.

YAN Ping， SHAN Xi， ZHAO Jingchen， et al. Microstructural stability of single crystal superalloy DD402 and its 1st stage turbine blades[J]. J Aeronautical Mat， 1999， 19（1）： 212.

[3]陳榮章. 單晶高溫合金發展現狀[J]. 材料工程， 1995， 23（8）： 212.

CHEN Rongzhang. Development status of single crystal superalloys[J]. J Mat Eng， 1995，23（8）： 212.

[4]張勇. 單晶高溫合金葉片用熔模殼型的研制[J]. 航空制造工程， 1998（4）： 2022.

ZHANG Yong. Study of mould shells for single crystal blades[J]. Aeronautical Manufacturing Eng， 1998（4）： 2022.

[5]張兵， 李田華， 周靜怡. 定向凝固和單晶高溫合金再結晶影響因素[J]. 航空科學技術， 2014， 25（5）： 109113.

ZHANG Bing， Li Tianhua， ZHOU Jingyi. Inlfuencing factors of recrystallization of directionally solidiifed and single crystal superalloys[J]. Aeronautical Sci & Technol， 2014， 25（5）： 109113.

[6]李影， 蘇彬. 鎳基單晶高溫合金的反常屈服行為與變形機制[J]. 材料工程， 2004， 32（3）： 4548.

LI Ying， SU Bin. Abnormal yield behavior and deformation mechanism of nickel base single crystal superalloy[J]. J Mat Eng， 2004， 32（3）： 4548.

[7]劉春廷， 孫曉峰， 管恒榮. 含錸單晶高溫合金鋁化物涂層的高溫氧化行為[J]. 航空材料學報， 2003， 23（S1）： 118122.

LIU Chunting， SUN Xiaofeng， GUAN Hengrong. High temperature oxidation behavior for aluminide coating on the single crystal superalloy containing rhenium[J]. J Aeronautical Mat， 2003， 23（S1）： 118122.

[8]朱鷗， 李玉龍， 張燕， 等. 航空發動機用單晶鑄造高溫合金熱處理工藝[J]. 鑄造技術， 2013， 34（9）： 11371140.

ZHU Ou， LI Yulong， ZHANG Yan， et al. Heat treatment process for singlecrystal superalloys used in aeroengines[J]. Foundry Technol， 2013， 34（9）： 11371140.

[9]丁智平. 復雜應力狀態鎳基單晶高溫合金低周疲勞損傷研究[D]. 長沙： 中南大學， 2005.

[10]趙萍. 航空發動機單晶葉片的多軸低周疲勞研究[D]. 長沙： 中南大學， 2011.

[11]孟杰， 金濤. 鎳基單晶高溫合金的再結晶[J]. 材料工程， 2011， 39（6）： 9298.

MENG Jie， JIN Tao. Recrystallization of single crystal nickelbased superalloys[J]. J Mat Eng， 2011， 39（6）： 9298.

[12]卿華， 江和甫， 溫衛東， 等. 鎳基單晶合金氣冷葉片模擬試樣的蠕變性能研究[J]. 航空動力學報， 2007， 22（5）： 773778.

QING Hua， JIANG Hefu， WEN Weidong， et al. Study on the creep behavior of model specimens of nickelbased single crystal aircooled blades[J]. J Aerospace Power， 2007， 22（5）： 773778.endprint

[13]薛明， 曹臘梅， 劉世忠， 等. 定向凝固過程中型芯、型殼溫度場數值模擬[J]. 熱加工工藝， 2006， 35（21）： 8688.

XUE Ming， CAO Lamei， LIU Shizhong， et al. Simulation on temperature fields of ceramic core and ceramic shell during directional solidification process[J]. Hot Working Technol， 2006， 35（21）： 8688.

[14]李嘉榮， 劉世忠， 史振學， 等. 第三代單晶高溫合金DD9[J]. 鋼鐵研究學報， 2011， 23（S2）： 337340.

LI Jiarong， LIU Shizhong， SHI Zhenxu， et al. Third generation single crystal superalloy DD9[J]. J Iron & Steel Res， 2011， 23（S2）： 337340.

[15]楊亮， 李嘉榮， 金海鵬， 等. DD6單晶精鑄薄壁試樣定向凝固過程數值模擬[J]. 材料工程， 2014， 42（11）： 1522.

YANG Liang， LI Jiarong， JIN Haipeng， et al. Numerical simulation on directional solidification process of DD6 single crystal superalloy thinwalled specimen[J]. J Mat Eng， 2014， 42（11）： 1522.

[16]于靖， 許慶彥， 柳百成. 鑄鋼件熔模精密鑄造凝固過程數值模擬[J]. 鑄造， 2006， 55（5）： 473476.

YU Jing， XU Qingyan， LIU Baicheng. Numerical simulation of solidification process of cast steel investment castings[J]. Foundry， 2006， 55（5）： 473476.

[17]PAN D， XU Q， YU J， et al. Numerical simulation of directional solidification of single crystal turbine blade casting[J]. Int J Cast Metals Res， 2008， 21（14）： 308312.

[18]張航， 許慶彥， 史振學， 等. DD6高溫合金定向凝固枝晶生長的數值模擬研究[J]. 金屬學報， 2014， 50（3）： 345354.

ZHANG Hang，XU Qingyan，SHI Zhenxue， et al. Numerical simulation of dendrite grain growth of DD6 superalloy during directional solidification process[J]. Acta Metallurgica Sinica， 2014， 50（3）： 345354.

[19]DUPRET F， van den BOGAERT N. Modelling bridgman and czochralski growth[EB/OL]. （20150801）[19981231]. http：//citeseerx.ist.psu.edu/showciting？cid=7866209.

[20]DUPRET F， NICODEME P， RYCKMANS Y， et al. Global modeling of heat transfer in crystal growth furnaces[J]. Int J Heat Mass Transfer， 1990， 33（9）： 18491871.（編輯武曉英）

（上接第17頁）

[7]臺冰豐， 吳杰靈. 清掃車吸塵口結構優化設計[J]. 專用汽車， 2014（4）： 1624.

TAI Bingfeng， WU Jieling. Optimized design of suction mouth of sweeper[J]. Spec Purpose Vehicle， 2014（4）： 1624.

[8]WU B， MEN J， CHEN J. Improving the design of a pickup head for particle removal using computational fluid dynamics[J]. J Mech Eng Sci， 2011， 255（4）： 938948.

[9]SHIH T H， LIOU W W， SHABBIR A， et al. A new kε eddy viscosity model for high Reynolds number turbulent flowsmodel development and validation[J]. Comput & Fluids， 1995， 24（3）： 227238.

[10]王福軍. 計算流體動力學分析[M]. 北京：清華大學出版社， 2004： 120122.

[11]張波， 文建軍. 旋風除塵器流場仿真分析[J]. 計算機輔助工程， 2011， 20（2）： 100103.

ZHANG Bo， WEN JIanjun. Simulation and analysis on flow field of cyclone separator[J]. Comput Aided Eng， 2011， 20（2）： 100103.

[12]MODLLING A. Tracer particles and seeding for particle image velocimetry[J]. Measurement Sci Technol， 1997， 8（12）： 14061416.

[13]HAIDER A， LEVENSPIEL O. Drag coefficient and terminal velocity of spherical and nonspherical particles[J]. Powder Technol， 1989， 58（1）： 6370.

[14]SAFFMAN P G. The lift on a small sphere in a slow shear flow[J]. J Fluid Mech， 1965， 22： 385400.（編輯武曉英）第24卷 第6期2015年12月計 算 機 輔 助 工 程Computer Aided EngineeringVol.24 No.6Dec. 2015endprint