徐春雷，薛秀生，樸成杰，楊寶興，劉自海

（中航工業(yè)沈陽發(fā)動機設計研究所，沈陽 110015）

0 引言

航空發(fā)動機的發(fā)展以高溫為主要特征，提高渦輪工作溫度是發(fā)展先進航空發(fā)動機的重要途徑。如果能夠掌握在發(fā)動機真實工作狀態(tài)下高壓渦輪葉片的表面溫度場分布，將對發(fā)動機性能研究和渦輪氣冷葉片的設計與優(yōu)化提供很大幫助[1]。因此，渦輪葉片表面溫度測量成為1個研究熱點。

紅外測溫技術(shù)作為1種先進的測溫方法，在國外航空發(fā)動機研制過程中得到了大量應用[2]。紅外測溫法屬于非接觸式測溫，相比傳統(tǒng)測試方法具有以下優(yōu)點[3]：

（a）不會破壞被測物體的表面；

（b）響應時間短，能實時測量不斷變化的溫度值；

（c）測溫范圍廣，靈敏度高。

因此，紅外測溫法在發(fā)動機熱端部件測溫領(lǐng)域具有十分廣闊的應用前景。但紅外測溫實際測量的是物體的輻射能量，需要已知表面發(fā)射率才能得到溫度值。RotamapⅡ紅外測溫系統(tǒng)在測試前需要輸入葉片平均發(fā)射率或者每片葉片的發(fā)射率數(shù)值[4]，所以提前獲知被測物體的表面發(fā)射率及其變化趨勢能夠有效地減少由于發(fā)射率所帶來的測溫誤差[5]。

1 發(fā)射率標定原理及公式

任何物體只要其溫度高于絕對零度，都會因為分子的熱運動而產(chǎn)生紅外輻射，而且輻射的能量隨溫度的變化而變化[6]。在光學研究中一般將物體分為黑體與灰體。絕對黑體的單色輻射能量E 隨波長變化的規(guī)律如普朗克定律所示[7]

式中：E 為物體的輻射能量，W/cm2·μm；C1、C2為普朗克第1、2輻射常數(shù)；λ 為波長，μm；T 為物體的絕對溫度，K。

然而，工程上遇到的物體一般都是灰體。對于灰體，其單色輻射能量為

式中：ε 為灰體的發(fā)射率，或者叫灰度系數(shù)。

通常紅外測溫儀得到的是灰體的單色輻射能量Eε，該儀器指示的溫度即為絕對黑體在單色輻射能量正好等于Eε時的溫度，稱為該灰體的亮度溫度，用符號Tl表示[8]。

此次測量所采用的紅外測溫儀波長λ=1μm，測溫范圍為600～1000℃，當λT=3120μm·K時，式（1）可簡化為[9]

根據(jù)亮度溫度的定義可知

由此可得

于是，灰體的表面溫度T 便可根據(jù)Tl用上式求算。同時可以得到灰體表面發(fā)射率

若已知物體的表面溫度T 和該輻射能量下黑體的亮度溫度Tl，便可得到物體的表面發(fā)射率。

2 試驗設備與測試方法

2.1 試驗設備

目前，渦輪葉片的表面發(fā)射率是紅外測溫領(lǐng)域的難點之一。國內(nèi)對發(fā)射率的測量實驗多在實驗室高溫爐內(nèi)進行。但由于高溫爐爐腔內(nèi)壁本身會產(chǎn)生較大輻射能量，會給被測件發(fā)射率的測量帶來極大誤差，若采用短腔高溫爐又無法實現(xiàn)大溫度范圍內(nèi)的測量[10-13]。作為創(chuàng)新點之一，本次試驗在中航工業(yè)動力所的熱機復合疲勞試驗器上進行，試驗設備如圖1所示。該試驗器主要由感應加溫系統(tǒng)、冷卻水系統(tǒng)、電氣控制系統(tǒng)以及測試系統(tǒng)構(gòu)成。通過采用電渦流加熱的方式可以實現(xiàn)大溫度范圍內(nèi)的發(fā)射率測量。而且感應銅管內(nèi)部有冷卻水通過，因此，銅管本身的溫度并不高，加熱裝置不會給測試結(jié)果帶來較大誤差。

圖1 熱機疲勞試驗器

2.2 試驗方法

RotamapⅡ紅外測溫系統(tǒng)的測溫位置為每片高壓渦輪葉片葉背區(qū)域，該區(qū)域比較平整所以發(fā)射率數(shù)值基本一致，較有代表性。因此，本次發(fā)射率測量位置為高壓渦輪葉片葉背區(qū)域。葉片測試改裝及電渦流線圈位置如圖2所示。圖2中的1為高壓渦輪葉片，2為電渦流線圈。通過調(diào)整葉片夾具和電渦流線圈的高度使加熱區(qū)域與渦輪葉片考核面對齊?？己嗣嫖挥谌~片葉背較平整區(qū)域，沿水平方向分布有2支熱點偶（GH3030，尺寸為φ1）測溫點。電渦流加熱渦輪葉片時，沿葉片垂直方向溫度差異較大，但在水平方向則溫度差異較小。2支熱電偶測溫點之間的水平距離為5mm，取2支熱電偶的平均溫度作為測試點（圖2中位置3）的溫度，即為渦輪葉片的表面溫度T。這樣處理既能減小葉片水平方向的溫度差異，又沒有破壞葉片的真實表面。RotamapⅡ系統(tǒng)使用條件較復雜，需提供氣源和發(fā)動機轉(zhuǎn)速信號。因此使用與RotamapⅡ紅外系統(tǒng)波長相同的紅外測溫儀Cyclops100測量位置3的溫度。初始設定紅外測溫儀發(fā)射率為1，可以得到該輻射能量下的亮度溫度Tl。在測試過程中，需注意紅外測溫儀垂直葉背平面測量，光線入射角需在±40°之間，紅外測溫儀用三角架固定，確保每次測試的位置不變。通過調(diào)整Cyclops100紅外測溫儀的焦距使被測點位置3在其視場中清晰可見為止。對新舊2片高壓渦輪葉片采用同樣的方法進行試驗。

圖2 改裝的高壓渦輪葉片及電渦流線圈位置

葉片表面的溫度加熱曲線如圖3所示?？拷~前緣的熱電偶1為溫度監(jiān)測點1，靠近葉尾緣的熱電偶2為溫度監(jiān)測點2。以靠近前緣的熱電偶作為溫度控制點。溫度穩(wěn)定時間為120s，在此期間使用紅外測溫儀測量2支電偶中間測溫點的溫度，并記錄相關(guān)數(shù)據(jù)。

圖3 電渦流加溫曲線

3 試驗結(jié)果分析

新、舊葉片的發(fā)射率標定結(jié)果見表1、2。其中表1為經(jīng)過長時間試車的舊葉片的試驗結(jié)果，表2為新葉片的試驗結(jié)果。T1、T2分別為熱電偶1、2測得的溫度值，T 為二者的平均溫度。亮度溫度Tl為采用紅外測溫儀Cyclops100在發(fā)射率設定為1時（被測表面真實發(fā)射率小于1）測得的溫度。

表1 舊葉片的發(fā)射率標定結(jié)果 ℃

表2 新葉片的發(fā)射率標定結(jié)果 ℃

通過表1、2中的數(shù)據(jù)和式（6）可以計算出各考核溫度下舊、新渦輪葉片的表面發(fā)射率數(shù)值，分別見表3、4。

通過表3、4中的數(shù)據(jù)可得舊、新葉片的表面發(fā)射率 平 均值0.864。通過2組對比試驗可以得到舊、新葉片在各溫度下發(fā)射率數(shù)值曲線，如圖4所示。其變化趨勢如圖中趨勢線所示。

表4 各溫度下新葉片的表面發(fā)射率

通過分析試驗數(shù)據(jù)可知，新、舊渦輪葉片的表面發(fā)射率會隨溫度的變化而改變。對于舊葉片，若取 εˉ1=0.893作為其發(fā)射率系數(shù)，則在600℃時引起的誤差較大，在其他溫度下的發(fā)射率誤差為±0.02，趨勢線斜率較小，數(shù)值變化不明顯。對于新葉片在各溫度點下發(fā)射率的偏差較大，且發(fā)射率明顯隨溫度升高而增大。當新渦輪葉片的表面溫度達到1000℃時，其表面發(fā)射率與舊葉片的表面發(fā)射率一致。因此，采用RotamapⅡ紅外系統(tǒng)時，若不知渦輪葉片溫度可使用0.893作為高壓渦輪葉片的表面發(fā)射率系數(shù)；若能估計出大致溫度，則可根據(jù)溫度值選擇發(fā)射率系數(shù)從而減小誤差。

英國RotamapⅡ紅外測試系統(tǒng)是1套專門用于測試發(fā)動機高壓渦輪葉片表面溫度場的設備，其技術(shù)指標見表5。

表5中的測試精度為不考慮葉片表面發(fā)射率誤差的情況。在使用該設備時可以選擇輸入平均發(fā)射率或者是每片葉片的發(fā)射率系數(shù)。通過試驗數(shù)據(jù)和普朗克方程可以得出，若使用平均發(fā)射率εˉ1=0.893，所造成的溫度誤差曲線，如圖5所示。在600℃下使用平均發(fā)射率造成的測溫誤差小于±5℃，在1200℃下使用平均發(fā)射率造成的誤差不超過±10℃，小于被測物體溫度的±1%，滿足測溫精度要求。若根據(jù)經(jīng)驗能估計出葉片溫度范圍，則可以結(jié)合如圖5所示發(fā)射率導致的測溫誤差數(shù)據(jù)對紅外測溫結(jié)果進行修正，在整級渦輪溫場測試中提高絕對溫度測試精度。

圖4 渦輪葉片表面發(fā)射率標定曲線

表5 RotamapⅡ紅外測試系統(tǒng)技術(shù)指標

圖5 發(fā)射率誤差導致的溫度誤差曲線

4 結(jié)束語

由于渦輪葉片表面發(fā)射率跟葉片材料、表面狀況、溫度、形狀有關(guān)，因此，必須使用真實的葉片進行發(fā)射率測量。如果條件允許，應當在發(fā)動機試車前后各進行1次發(fā)射率測試工作。此次發(fā)射率標定試驗主要測量了高壓渦輪葉片的表面發(fā)射率，研究了新、舊葉片間的發(fā)射率數(shù)值及變化趨勢。對提高紅外測溫精度有一定意義。同時為測試發(fā)動機零部件表面發(fā)射率提供了1種有效的方法。

本次試驗未對表面涂有隔溫涂層的葉片進行研究，需要在后續(xù)的試驗過程中逐步探索。對于渦輪葉片在燃燒過程中發(fā)生積炭等情況，需要在熱風洞試驗器上進行后續(xù)試驗[14-15]。

[1]劉大響，程榮輝.世界航空動力技術(shù)的現(xiàn)狀及發(fā)展動向[J].北京航空航天大學學報，2002，28（5）：490-493.LIU Daxiang，CHENG Ronghui.Current status and development direction of aircraft power technology in the world[J].Journal of Beijing University of Aeronautics and Astronautics，2002，28（5）：490-493.（in Chinese）

[2]王如根，高坤華.航空發(fā)動機新技術(shù)[M].北京：航空工業(yè)出版社，2003：10-85.WANG Rugen，GAO Kunhua.The New technology of aeroengine[M].Beiing：Aviation Industry Press，2003：10-85.（in Chinese）

[3]戴景民.輻射測溫的發(fā)展現(xiàn)狀與展望[J].自動化技術(shù)與應用，2004，57（3）：1-5.DAI Jingmin.Survey of radiation thermometry[J].Techniques of Automation and Applications，2004，57（3）：1-5.（in Chinese）

[4]熊兵，石小江，陳洪敏，等.航空發(fā)動機渦輪葉片發(fā)射率測量[J].燃氣渦輪試驗與研究，2011，24（2）：45-48.XIONG Bing，SHI Xiaojiang，CHEN Minhong，et al.Emissivity measurement of aeroengine turbine blade[J].Gas Turbine Experiment and Research，2011，24（2）：45-48.（in Chinese）

[5]Incropera F P，Dewit D P，Bergman T L，等.傳熱和傳質(zhì)基本原理[M].北京：化學工業(yè)出版社，2007：97-190.Incropera F P，Dewit D P，Bergman T L，et al.The basic principles of heat and mass transfer[M].Beijing：Chemical Industry Press，2007：97-190.（in Chinese）

[6]吳燕燕，羅鐵茍，黃杰.基于紅外熱像儀測溫原理的物體表面發(fā)射率計算[J].直升機技術(shù)，2011，169（4）：25-29.WU Yanyan，LUO Tiegou，HUANG Jie.Calculation of object surface emissivity based on measuring temperature principle of thermal imager[J].Helicopter Technique，2011，169（4）：25-29.（in Chinese）

[7]林其勛，游紹堃.航空發(fā)動機氣動參數(shù)測量[M].北京：國防工業(yè)出版社，1988：247-256.LIN Qixun，YOU Shaokun，et al.Aeroengine aerodynamic parameters measuring[M].Beijing：Defense Industry Press，1988：247-256.（in Chinese）

[8]Markham J R，Latvakoski H M，F(xiàn)rank S L F，et al.Simultaneous short and long wavelength infrared pyrometer measurements in a heavy duty gas turbine[R].ASME 2001-GT-0026.

[9]Frank S L F，Holt T O，Eisenlohr H，et al.Application of a high resolution turbine pyrometer to heavy duty gas turbines[R].ASME 2001-GT-0577.

[10]Inagaki T，Okamoto Y.Surface temperature measurement near ambient conditions using infrared radiometers with different detection wavelength bands by applying a grey-body approximation:estimation of radiative properties for non-matal surfaces[J].NDT & E International，1996（6）：363-369.

[11]鄧洪偉，邵萬仁，周勝田，等.某型航空發(fā)動機噴管紅外輻射特征數(shù)值模擬和試驗研究[J].航空發(fā)動機，2010，36（1）：45-48.DENG Hongwei，SHAO Wanren，ZHOU Shengtian，et al.Numerical simulation and experimental investigation of infrared radiation characteristics for an aeroengine exhaust nozzle[J].Aeroengine，2010，36（1）：45-48.（in Chinese）

[12]Char J M，YEH J H.The measurement of open propane flame temperature using infrared technique[J].Journal of Quantitative Spectroscopy and Radiative Transfer，1996，（1）:133-144.

[13]John D.Wrbanek W.Development of thin film ceramic thermocouples for high temperature environments[C].the MAA，ASME，SAE，and ASEE Fort Lauderdale，F(xiàn)lorida，2004.11-14.

[14]馬春武，劉忠奎，薛秀生，等.航空發(fā)動機表面發(fā)射率測量技術(shù)研究[J].航空發(fā)動機，2012，38（6）：58-62.MA Chunwu，LIU Zhongkui，XUE Xiusheng，et al.Investigation on surface emissivity measurement technologies for aeroengine[J].Aeroengine，2012，38（6）：58-62.（in Chinese）

[15]閆久坤.燃氣輪機渦輪出口溫度受感部故障分析與改進設計[J].航空發(fā)動機，2011，37（5）：33-37.YAN Jiukun.Failure analysis and improved design of turbine outlet temperature probe for gas turbine[J].Aeroengine，2011，37（5）：33-37.（in Chinese）