賈 超, 蔡 杰, 熊朝暉
(中國工程物理研究院計量測試中心,四川綿陽621900)
溫度傳感器作為動態(tài)溫度測量的主要測量手段,其動態(tài)特性是影響動態(tài)溫度測量準(zhǔn)確性的重要參數(shù)[1,2]。為了確定動態(tài)特性參數(shù),國內(nèi)外學(xué)者進(jìn)行了大量研究,形成了投擲法(油浴法、水浴法)[3~5]、激波管法、電加熱法[6]、熱風(fēng)洞法[7,8]和脈沖激光法[9~11]等動態(tài)測試方法,但至今仍未形成統(tǒng)一的動態(tài)特性參數(shù)標(biāo)定方案。究其原因主要是在實際測量過程中,動態(tài)溫度場環(huán)境較復(fù)雜,各種傳熱過程同時發(fā)生,實驗室條件下測量的動態(tài)特性參數(shù)與實際值差別較大,無法完全滿足實際需求[12]。因此,研究不同傳熱方式下溫度傳感器的動態(tài)特性對于溫度傳感器在動態(tài)溫度測量領(lǐng)域的研究和應(yīng)用具有一定的參考價值。
傳感器的動態(tài)響應(yīng)特性是指溫度傳感器的溫度與被測介質(zhì)溫度增量之間的關(guān)系,其特征參數(shù)主要有熱響應(yīng)時間和時間常數(shù)。熱響應(yīng)時間是影響測試結(jié)果的重要因素,而時間常數(shù)是反映其動態(tài)響應(yīng)性能最重要的指標(biāo)[13]。
本文以工業(yè)熱電偶為研究對象,通過時間常數(shù)理論推導(dǎo),分析時間常數(shù)的主要影響因素;根據(jù)時間常數(shù)推導(dǎo)結(jié)果研究不同傳熱方式與動態(tài)特性之間的關(guān)系,并通過理論推導(dǎo)結(jié)合實驗驗證的方式分析了不同傳熱方式下溫度傳感器的動態(tài)響應(yīng)特性,實驗結(jié)果驗證了理論推導(dǎo)的正確性。
假設(shè)將體積V、密度ρ、比熱容c、表面積S的熱電偶放置于溫度為Te的恒定溫場中,動態(tài)平衡時熱電偶的能量方程可表示為:
qh=qc
(1)
式中:qh為對流傳熱量;qc為熱電偶熱結(jié)點儲熱量。
經(jīng)過時間t后,熱電偶指示溫度為T,那么:
qh=Sh(Te-T)
(2)
(3)
式中:h為熱電偶與溫場介質(zhì)之間的表面?zhèn)鳠嵯禂?shù)。
將式(2)、式(3)代入式(1)式得:
(4)
在初始條件t=0,T=T0時,對式(4)積分簡化可得:
(5)
(6)
圖1 熱電偶階躍響應(yīng)曲線上的時間常數(shù)Fig.1 Time constant on the thermocouple step response curve
從表達(dá)式可以看出,時間常數(shù)是一個多影響量參數(shù),除了傳感器自身特性(材料、結(jié)構(gòu))的影響外[14],傳感器與溫場介質(zhì)之間的表面?zhèn)鳠嵯禂?shù)為主要影響量,即傳熱方式。
為了得到正確的時間常數(shù)值,分析不同傳熱方式下熱電偶的動態(tài)響應(yīng)成為關(guān)鍵。目前熱電偶時間常數(shù)的測試方法主要有熱風(fēng)洞法、水浴法、油浴法和脈沖激光法等,按傳熱方式可分為對流傳熱和輻射傳熱[15]。
當(dāng)采用熱風(fēng)洞法、水(油)浴法作為階躍溫度激勵方法時,忽略自身熱傳導(dǎo),熱電偶與外界的主要傳熱方式為對流傳熱。無論對流傳熱過程的具體特性如何,其能量平衡方程可用牛頓冷卻公式表示:
(7)
式中:h1表示熱電偶與硅油之間的表面?zhèn)鳠嵯禂?shù)。
當(dāng)激勵方式為脈沖激光時,激光主要以輻射傳熱作用在熱電偶表面,但隨著溫度的升高,熱電偶會以對流傳熱的方式與周圍空氣發(fā)生熱交換。根據(jù)能量守恒關(guān)系,熱電偶的能量平衡方程可表示為:
Q-(QF+QD)=qc
(8)
式中:Q表示熱電偶吸收激光能量;QF表示熱電偶受激光照射的凈輻射換熱量;QD表示熱電偶與周圍環(huán)境的換熱量。
根據(jù)斯特藩-玻爾茲曼定律可得:
(9)
式中:α為熱電偶吸收率;ε為熱電偶表面發(fā)射率;σ為斯特藩-玻爾茲曼常數(shù);G表示輻照功率;E表示發(fā)射功率。
對熱電偶而言,若忽略自身熱傳導(dǎo),在瞬態(tài)條件下其能量平衡方程可表示為:
(10)
式中:αI為熱電偶對激光的吸收系數(shù);I為激光的功率密度;h2為熱電偶與空氣之間的表面?zhèn)鳠嵯禂?shù)。
式(10)是一個非線性、一階、非齊次常微分方程,無法積分得到嚴(yán)格解,但利用數(shù)值積分方法或有限差分方法可以得到該微分方程解。
為了獲取對流傳熱和輻射傳熱方式下熱電偶的動態(tài)響應(yīng),分別建立了基于恒溫油槽和脈沖激光器的時間常數(shù)測試系統(tǒng),其系統(tǒng)裝置如圖2所示。
圖2 時間常數(shù)測試系統(tǒng)Fig.2 Time constant test system
熱電偶動態(tài)響應(yīng)測試系統(tǒng)激勵源分別選用湖州晶覓試驗設(shè)備有限公司HTS-300A型號恒溫油槽和北京鐳志威光電技術(shù)公司研制的連續(xù)激光器,數(shù)據(jù)采集裝置為Agilent DSO8104A型數(shù)字示波器,采集速率為4 GSa/s。實驗過程中,將示波器設(shè)置為滾動采集模式。采用投擲法將熱電偶測量端突然置于激勵溫場中,確保熱電偶輸入為溫度階躍信號,同時示波器采集得到熱電偶電勢輸出值與響應(yīng)時間曲線。
本文選用重慶川儀一廠WRNK型熱電偶為實驗對象,其有關(guān)參數(shù)為熱結(jié)點直徑D=2 mm,導(dǎo)熱系數(shù)λ=20 W/(m·K),比熱容c=400 J/(kg·K),密度ρ=8 500 kg/m3。
在對流傳熱方式下,采用恒溫槽作為溫度激勵源,介質(zhì)為有機(jī)硅油,其傳熱系數(shù)h=1 100 W/(m2·K),根據(jù)集總熱容應(yīng)用條件[16],熱電偶畢渥數(shù)Bi=hLc/λ=0.03<0.1(Lc稱為定性長度,其值近似等于D/3),采用集總熱容法是極好的近似,因此在對流傳熱方式下研究熱電偶動態(tài)特性時可忽略其自身熱傳導(dǎo)帶來的影響。設(shè)定恒溫油槽溫度100 ℃,溫度波動度小于0.02 ℃/10 min,在穩(wěn)定環(huán)境溫度下,將熱電偶突然置于恒溫槽中,產(chǎn)生一個80 ℃的溫度階躍信號,示波器采集熱電偶動態(tài)響應(yīng)曲線如圖3所示。
圖3 階躍溫度80 ℃時熱電偶輸出響應(yīng)曲線Fig.3 Thermocouple output response curve at step temperature 80 ℃
由對流傳熱方式下熱電偶能量平衡方程可知,80 ℃油浴階躍輸入信號下熱電偶的時間常數(shù)理論值為1.0 s。由圖3可知,相同條件下熱電偶達(dá)到動態(tài)平衡時的時間為6.2 s左右,由時間常數(shù)的定義可得時間常數(shù)實際測量值為1.2 s,與理論值基本一致。
在輻射傳熱方式下,熱電偶與周圍空氣發(fā)生熱交換,其傳熱系數(shù)為h=20 W/(m2·K),同理可得畢渥數(shù)Bi=hLc/λ=2×10-3<0.1,因此可忽略自身熱傳導(dǎo)。該方法采用連續(xù)激光器作為激勵源,輸出波長為1 064 nm,輸出功率0~5 W可調(diào),出光孔光斑<5 mm,輸出功率穩(wěn)定性小于1%。常溫下,設(shè)置激光器輸出功率5 W,并聚焦至熱電偶熱結(jié)點上。數(shù)字示波器采集到的動態(tài)響應(yīng)曲線如圖4所示。由圖4可以看出,達(dá)到動態(tài)平衡時熱電偶熱電勢變化值為11.78 mV,根據(jù)K型熱電偶電勢值-溫度換算公式可得階躍溫度為290 ℃。
圖4 激光輸出功率5 W時熱電偶輸出響應(yīng)曲線Fig.4 Thermocouple output response curve with laser output power of 5 W
根據(jù)式(10),采用數(shù)值積分方法,可得在輻射傳熱方式下熱電偶的理論時間常數(shù)為4.3 s;根據(jù)時間常數(shù)定義,從圖4可以看出達(dá)到溫度階躍量的63.2%時實際測量時間常數(shù)為6.2 s。
針對不同傳熱方式下溫度傳感器的動態(tài)響應(yīng)問題,本文以K型工業(yè)熱電偶為研究對象,從時間常數(shù)的定義出發(fā),分析其影響因素除了傳感器本身外,與溫場介質(zhì)間的傳熱系數(shù)為主要影響量;然后通過理論推導(dǎo)對流傳熱和輻射傳熱方式下熱電偶的能量平衡方程,分析其動態(tài)特性,并開展了兩種傳熱方式下的時間常數(shù)測試實驗,從實驗結(jié)果可以看出,在對流傳熱方式下,熱電偶時間常數(shù)理論值與實際值基本一致;但在輻射傳熱方式下,熱電偶對激光能量的吸收情況決定了時間常數(shù)測試值的大小,影響因素較多,導(dǎo)致時間理論值與實驗值存在較大偏差。因此在研究溫度傳感器動態(tài)特性時,要充分考慮傳熱方式的不同對測量結(jié)果帶來的影響,根據(jù)實際環(huán)境條件確定傳熱過程,以確保測量結(jié)果的準(zhǔn)確性與可靠性。