熱電偶動態(tài)特性的溫度預估方法實驗

摘要:針對高溫高頻和非均勻熱流密度條件下溫度測量問題,采用激光連續(xù)加熱熱電偶的方法,研究了不同熱流密度條件下、不同直徑熱電偶的動態(tài)響應(yīng)特性,分析了熱流密度對熱電偶動態(tài)響應(yīng)時間和電壓變化率的影響,得到熱電偶動態(tài)響應(yīng)時間、電壓變化率與熱流密度的函數(shù)關(guān)系,并提出了一種基于熱電偶動態(tài)響應(yīng)特性溫度預估方法。結(jié)果表明:熱電偶動態(tài)特性受熱流密度影響較大,在相同熱流密度條件下,動態(tài)響應(yīng)時間與電壓變化率的指數(shù)函數(shù)呈線性關(guān)系;文中提出的溫度預估方法在超出熱電偶測溫量程27.3%的范圍內(nèi)具有較小誤差,能滿足工程溫度測量的需要。
0引言
　　在火箭發(fā)動機的研制和試驗過程中,發(fā)動機燃燒室燃氣溫度的測量對于評估推進劑能量水平、發(fā)動機內(nèi)流場特性以及熱防護系統(tǒng)的設(shè)計都具有重要的意義。目前,在工程中主要應(yīng)用金屬熱電偶絲為測溫元件的接觸式測溫方法,由于受到熱電偶絲熔點的限制,接觸式測溫方法的測溫范圍受到一定的限制,在測量過程中經(jīng)常出現(xiàn)熱電偶結(jié)點熔斷或在發(fā)動機稠密兩相燃氣流中沖蝕破壞的現(xiàn)象，導致無法獲得發(fā)動機.穩(wěn)定工作時燃燒室的溫度數(shù)據(jù)。
　　本文針對火箭發(fā)動機接觸式測溫中經(jīng)常出現(xiàn)的熱電偶絲熔斷和沖蝕破壞現(xiàn)象，利用熱電偶在破壞前所采集到的溫度數(shù)據(jù),并結(jié)合熱電偶的動態(tài)響應(yīng)特性,提出了一種基于熱電偶動態(tài)響應(yīng)特性的溫度預估方法。熱電偶的動態(tài)響應(yīng)時間是動態(tài)溫度測量的一個重要參數(shù),也是衡量溫度傳感器動態(tài)響應(yīng)性能的重要指標。熱電偶動態(tài)響應(yīng)時間是指當熱電偶從某-初始溫度快速放人溫度恒定的介質(zhì)中后,熱電偶的溫度變化量,由零達到介質(zhì)溫度與溫度傳感器初始溫度之差的63.2%所用的時間。國內(nèi)外研究表明(7-8]，激光作為加熱源具有能量集中,功率穩(wěn)定、熱流密度可調(diào)、延遲短、安全等優(yōu)點,能得到較精確的研究結(jié)果。
　　本文利用激光器作加熱源,研究了熱流密度對熱電偶動態(tài)響應(yīng)特性的影響，建立了熱電偶升溫階段電壓變化率與響應(yīng)時間的關(guān)系，提出基于熱電偶動態(tài)特性的溫度預估方法,通過測量常溫常壓下3種典型推進劑的火焰溫度來驗證該溫度預估方法的可靠性，為火箭發(fā)動機高溫高頻瞬態(tài)溫度的測量及預估提供技術(shù)支持。
1實驗裝置
　　為了研究不同熱流密度下、不同結(jié)點直徑熱電偶的動態(tài)響應(yīng)特性，文中自制了結(jié)點直徑分別為0.24、0.34.0.42mm的鎳鉻-鎳硅K型熱電偶傳感器,熱電偶絲直徑為0.1mm。實驗前，用西派埃公司的WJL-1A型熱電偶校驗爐對熱電偶進行溫度標定。熱電偶動態(tài)響應(yīng)測試系統(tǒng)如圖1所示，包括CO2激光器、光電傳感器、平面鍍金全反射鏡、聚焦鏡、密閉燃燒器、K型鎳鉻鎳硅熱電偶、信號采集系統(tǒng)及工控機。激光器作為加熱源提供不同熱流密度的激光輻射能,激光器發(fā)出的光束經(jīng)平面鍍金全反射鏡人射到聚焦鏡，兩鏡面距離為80mm,平面反射鏡與輸人激光光束成45°夾角。經(jīng)聚焦鏡后,激光光束直接照射到熱電偶結(jié)點上，激光光斑直徑為3mm。熱電偶通過補償導線與采集系統(tǒng)連接,采集頻率為1kHz。用光電傳感器測得的激光出光信號作為系統(tǒng)的觸發(fā)信號,光電傳感器測得激光信號和熱電偶測得的溫度響應(yīng),在采集系統(tǒng)中以電壓形式輸出,采集信號最后經(jīng)濾波處理。激光器的輸出功率調(diào)整及加熱時間的由工控機控制。
 
2熱電偶動態(tài)響應(yīng)特性
2.1熱電偶校驗
　　為了將熱電偶電壓信號準確地轉(zhuǎn)換為溫度值,需對自制的熱電偶進行校驗。熱電偶的原位校驗熱電偶校驗爐中進行，溫度范圍從25~1200℃,溫度增量25℃,熱電偶在標準溫度下的恒定時間為5min。
　　熱電偶校驗結(jié)果如圖2所示，并對校驗數(shù)據(jù)進行三階多項式擬合,對于結(jié)點直徑分別為0.24、0.34、0.42mm的熱電偶，溫度與電壓之間符合如下關(guān)系:
 
　　另外，由圖2中校驗曲線可看出,不同結(jié)點直徑熱電偶的校驗結(jié)果之間存在一定誤差,這主要由熱電偶焊接工藝和焊接質(zhì)量造成的。
2.2不同熱流密度下熱電偶典型響應(yīng)曲線
　　利用激光連續(xù)加熱方法,研究了不同熱流密度條件下不同結(jié)點直徑熱電偶的動態(tài)響應(yīng)特性。實驗中，熱電偶結(jié)點直徑分別為0.24、0.340.42mm,激光加熱時間為4s,激光熱流密度分別為425.637.778、849.948W/cm²。
 
　　圖3所示為不同熱流密度下不同結(jié)點直徑熱電偶的典型響應(yīng)曲線。如圖3(a)所示,根據(jù)電壓隨時間的變化特點,可將熱電偶響應(yīng)曲線分為3個階段:升溫段I、熱平衡段II和冷卻段M。升溫段I:在激光光束的輻射加熱作用下,熱電偶結(jié)點吸收的熱量大于向外輻射和自然對流損失的熱量,結(jié)點溫度逐漸升高;熱平衡段II:當熱電偶結(jié)點吸收的激光光束熱量等于損失的能量時,熱電偶結(jié)點溫度保持在某一恒值;冷卻段皿:當激光光束停止加熱后，熱電偶結(jié)點溫度逐漸降低。
　　圖3(a)~(c)所示為激光熱流密度分別為425、849、948W/cm²時不同直徑熱電偶的響應(yīng)曲線。由圖3可知,激光熱流密度越小,升溫段曲線的斜率越小，平衡電壓越小。當激光熱流密度增大時,直徑為0.24.mm的熱電偶的電壓變化率最為敏感,這說明熱電偶結(jié)點直徑越小，響應(yīng)越快,響應(yīng)時間越短。
　　實驗中發(fā)現(xiàn),在較高的激光熱流密度情況下,熱電偶結(jié)點出現(xiàn)燒蝕現(xiàn)象[9],燒蝕程度-一般為納米或微米級。當激光熱流密度超過某一特定值時,熱電偶結(jié)點會被燒壞而發(fā)生熔斷，這主要與熱電偶絲熔點有關(guān)。
2.3熱流密度對熱電偶動態(tài)特性的影響
　　研究激光熱流密度對不同直徑熱電偶的動態(tài)響應(yīng)特性的影響。不同直徑熱電偶輸出電壓分別通過對應(yīng)的關(guān)系式(8)~式(10)轉(zhuǎn)換為溫度值。圖4為不同直徑熱電偶響應(yīng)時間隨熱流密度q的變化曲線。由圖4可知，熱流密度對熱電偶響應(yīng)時間的影響較大。對于結(jié)點直徑相同的熱電偶，熱流密度越大,響應(yīng)時間越短;熱電偶結(jié)點直徑越大，響應(yīng)時間越長。
 
　　圖5所示為不同直徑熱電偶在升溫段的電壓變化率隨熱流密度q的變化曲線。由圖5可知，熱流密度對升溫段電壓的變化率有較大影響,對于不同直徑的熱電偶,熱流密度對電壓變化率的影響趨勢基本是一致的。相同直徑的熱電偶，升溫段的電壓變化率隨著熱流密度的增大而增大。熱電偶結(jié)點直徑越小,熱流密度對電壓變化率的影響越大。
 
2.4響應(yīng)時間1.電壓變化率dU/dt與熱流密度q的相關(guān)關(guān)系
　　為了研究響應(yīng)時間、電壓變化率與熱流密度之間的關(guān)系,利用最小二乘法對圖4和圖5中的變化曲線進行擬合,擬合結(jié)果分別見圖4和圖5中的擬合曲線，相關(guān)系數(shù)都在0.99以上。不同直徑熱電偶響應(yīng)時間和電壓變化率與熱流密度的相關(guān)關(guān)系分別為
 
　　由式(4)~式(9)可知,電壓變化率與熱流密度呈二階多項式函數(shù)關(guān)系,隨著結(jié)點直徑的增大，電壓變化率與熱流密度擬合函數(shù)的常數(shù)項逐漸增大,-次項和二次項系數(shù)逐漸減小;響應(yīng)時間與熱流密度的指數(shù)呈線性函數(shù)關(guān)系,隨著結(jié)點直徑的增大，擬合函數(shù)的常數(shù)項逐漸減小,而熱流密度的系數(shù)逐漸增大。
3溫度預估方法的建立及驗證
3.1溫度預估方法建立
　　使用熱電偶進行溫度測量時,由于受到熱電偶絲熔點的限制,測溫范圍有限。因此,在測量溫度超過熱電偶測溫極限的溫度場時，會出現(xiàn)熱電偶絲熔斷現(xiàn)象，這種現(xiàn)象在固體火箭發(fā)動機實驗中經(jīng)常出現(xiàn)。為了能依據(jù)熱電偶破壞前的數(shù)據(jù)來預測被測溫度場的溫度,本文基于熱電偶的動態(tài)響應(yīng)特性提出一種溫度預估方法。
　　假設(shè)在超過熱電偶測溫極限-定范圍時,測量溫度與熱電偶響應(yīng)電壓之間仍符合熱電偶標定結(jié)果。根據(jù)實驗研究結(jié)果得到的熱電偶電壓變化率、動態(tài)響應(yīng)時間與熱流密度的關(guān)系,消去中間變量q,得到熱電偶動態(tài)響應(yīng)時間與電壓變化率的函數(shù)關(guān)系,即
 
　　表1中給出了對應(yīng)于結(jié)點直徑分別為0.24、0.34、0.42mm熱電偶的具體函數(shù)關(guān)系式。熱電偶測量溫度時,通常得到的是熱電偶響應(yīng)電壓隨時間的變化關(guān)系。然后,通過熱電偶標定結(jié)果將電壓轉(zhuǎn)換為溫度。在文中進行溫度預估時，針對某--特定結(jié)點直徑的熱電偶，利用熱電偶絲熔斷前采集到的響應(yīng)電壓隨時間的變化數(shù)據(jù),得到本次溫度測量過程中熱電偶的電壓變化率，由式(10)可計算得到對應(yīng)結(jié)點直徑的熱電偶動態(tài)響應(yīng)時間tc同時由測量數(shù)據(jù)可得到對應(yīng)的電壓幅值U。根據(jù)熱電偶動態(tài)響應(yīng)時間的定義,有
 
　　運用反演法思想，利用得到的Uc代入式(11)反推出熱電偶的平衡電壓。然后,根據(jù)熱電偶標定結(jié)果計算出平衡溫度。
 
3.2試驗驗證
　　為了驗證提出的溫度預估方法的可靠性,利用自制鎳鉻K型熱電偶和鎢錸G型熱電偶在圖1所示裝置中同時對典型雙基推進劑.改性雙基推進劑和復合推進劑室溫常壓下的燃燒火焰溫度進行測量,所用鎳銘熱電偶結(jié)點直徑為0.42mm,鎢錸熱電偶結(jié)點與鎳鉻熱電偶結(jié)點間距離為1mm,測量點距離推進劑表面的高度通過預先實驗確定,分別為8mm(雙基推進劑和改性雙基推進劑)、10mm(復合推進劑)。溫度測量結(jié)果以鎢錸熱電偶所測溫度為參考值,每種推進劑進行5次重復實驗,溫度值取5次試驗的平均值。
　　推進劑燃燒過程中,不同火焰區(qū)的燃燒溫度各不相同,為了使熱電偶所測溫度為某一恒定值，在測量過程中,熱電偶需一直處于同一火焰區(qū)。但在藥柱燃燒過程中,燃面發(fā)生退移,燃面與熱電偶結(jié)點之間的測量高度會發(fā)生變化。為了保證熱電偶結(jié)點與推進劑燃面之間保持在預定的測量距離,本文測量時采用了圖6所示的藥柱燃面調(diào)節(jié)裝置,主要由外殼、推進劑藥柱、托板.彈簧和螺蓋組成。實驗時,將熱電偶調(diào)整至預定高度并固定。用激光器點燃推進劑，然后迅速將圖6所示的裝置平行推至預定位置。推進劑燃燒時,燃面發(fā)生退移。與此同時,推進劑藥柱底部在彈簧力的作用下向上移動,使得燃面與熱電偶結(jié)點之間的距離基本保持不變,從而保證熱電偶所測的溫度為同--火焰區(qū)的溫度值,熱電偶達到平衡時能夠穩(wěn)定在某一恒定溫度值。
 
　　圖7為常壓.室溫時所測得的3種不同類型推進劑的燃燒火焰溫度變化曲線。其中,圖7(a)為雙基推進劑火焰溫度變化曲線，平衡段溫度平均值為1075.2℃;圖7(b)為改性雙基推進劑火焰溫度變化曲線，平衡段溫度平均值為1391.3℃;圖7(c)為復合推進劑火焰溫度變化曲線,平衡段溫度平均值1527.8℃。
 
　　由圖7可知,雙基推進劑室溫常壓下的燃燒溫度最低,溫度震蕩最小;改型雙基推進劑燃燒溫度明顯提高，且溫度曲線伴隨輕微的震蕩;復合推進劑燃燒溫度最高,溫度曲線震蕩最大。這是因為燃燒時，推進劑燃面處有大量的金屬顆粒向外噴射,火焰不均勻,引起溫度震蕩較大。另外,圖7(a)中鎳鉻熱電偶和鎢錸熱電偶所測火焰溫度基本一致,說明熱電偶結(jié)點間距對溫度場測量干擾較小,間距適當;圖7(b)和(c)中鎳鉻熱電偶在火焰溫度大于1250℃時被燒斷,這主要是由鎳鉻熱電偶絲的熔點決定的。
　　依據(jù)文中建立的溫度預估方法,利用鎳鉻熱電偶測得的不完整溫度數(shù)據(jù),對3種推進劑的燃燒火焰溫度進行了預估,結(jié)果見表2。由表2不難看出，預測結(jié).果均小于參考溫度。這主要是因為根據(jù)鎳鉻熱電偶熔斷前數(shù)據(jù)中得到的電壓變化率小于實際電壓變化率,從而引起計算得到的動態(tài)響應(yīng)時間偏小,最終導致預估的平衡電壓值偏小，預估溫度偏小。另外,表2中數(shù)據(jù)表明,在超出熱電偶測溫極限15.9%和27.3%時，預估值誤差分別為5.7%和7.2%，均小于10%，說明文中假設(shè)的合理性。
 
4結(jié)論
(1)熱電偶動態(tài)特性受熱流密度影響較大;在相同熱流密度條件下，熱電偶動態(tài)響應(yīng)時間與升溫段電壓變化率的指數(shù)函數(shù)呈線性函數(shù)關(guān)系。
(2)文中基于熱電偶動態(tài)特性提出的溫度預估方法在超出熱電偶測溫極限27.3%的范圍時具有較小.誤差,能滿足工程溫度測量的需要。