熱電偶動(dòng)態(tài)特性的溫度預估方法研究

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