柳曉寧，郭 鵬，申 彬，曹志松,2

（1.北京衛(wèi)星環(huán)境工程研究所；2.可靠性與環(huán)境工程技術(shù)重點(diǎn)實(shí)驗(yàn)室：北京 100094）

0 引言

航天器真空熱試驗(yàn)是航天器研制過程中的重要環(huán)節(jié)，是驗(yàn)證航天器熱控設(shè)計(jì)正確性、暴露有關(guān)設(shè)計(jì)和質(zhì)量缺陷的必要手段。紅外模擬是我國(guó)當(dāng)前航天器真空熱試驗(yàn)中最常用的外熱流模擬手段，其中使用較多的外熱流模擬裝置有紅外加熱籠、紅外燈陣和薄膜加熱片。紅外加熱棒在國(guó)外航天器真空熱試驗(yàn)中有所應(yīng)用，但在我國(guó)航天器真空熱試驗(yàn)中的應(yīng)用尚處于初步研究階段[1]。

本文針對(duì)新研紅外加熱棒式外熱流模擬器進(jìn)行結(jié)構(gòu)建模，對(duì)其外熱流模擬效果進(jìn)行數(shù)值計(jì)算分 析，并通過真空低溫試驗(yàn)對(duì)模擬器加熱能力進(jìn)行驗(yàn)證。研究結(jié)果將對(duì)紅外加熱棒式外熱流模擬器在航天器真空熱試驗(yàn)中的實(shí)際應(yīng)用具有一定指導(dǎo)意義。

1 紅外加熱棒特性分析

1.1 結(jié)構(gòu)及材料屬性

本文選擇技術(shù)成熟度高的金屬紅外加熱棒為研究對(duì)象，其結(jié)構(gòu)為在無縫金屬管內(nèi)放入電熱絲，并在空隙部分緊密填充具有良好導(dǎo)熱性和電絕緣性的晶體氧化鎂。金屬紅外加熱棒具有結(jié)構(gòu)簡(jiǎn)單、設(shè)計(jì)靈活、加熱效率高、機(jī)械強(qiáng)度大、壽命長(zhǎng)、 以及使用安全等優(yōu)點(diǎn)。

常見紅外加熱棒的電極引線一般有雙端出線和單端出線2 種方式，本文所選加熱棒為單端出線。加熱棒最外層套管為不銹鋼材料，內(nèi)部填充絕緣材料為晶體氧化鎂，電熱絲為鎳鉻絲。加熱棒的截面如圖1所示，組成材料物性見表1。

圖1 加熱棒截面示意 Fig.1 Cross section of the infrared heating rod

表1 紅外加熱棒組成材料物性參數(shù) Table1 Physical properties of the infrared heating rod material

1.2 穩(wěn)態(tài)溫度分析

紅外加熱棒工作時(shí)，電熱絲產(chǎn)生的熱量通過絕緣層和不銹鋼套管，以輻射、對(duì)流或?qū)岬姆绞绞┘拥奖患訜嵛矬w，在達(dá)到溫度穩(wěn)定狀態(tài)后，內(nèi)部溫度將不再變化。對(duì)于電熱絲、絕緣層和不銹鋼套管，其內(nèi)部導(dǎo)熱微分方程[2]均可表示為

式中：t為加熱棒內(nèi)部某點(diǎn)(r,φ,z)的溫度，r,φ,z分別為內(nèi)部某一點(diǎn)的徑向坐標(biāo)、方位角和高度；為內(nèi)熱源強(qiáng)度；λ為材料導(dǎo)熱系數(shù)。紅外加熱棒的軸向尺寸遠(yuǎn)大于其徑向尺寸，且氧化鎂絕緣層導(dǎo)熱性能為各向同性，因此為了簡(jiǎn)化分析，可以忽略加熱棒內(nèi)部溫度隨方位角和高度的變化，將式(1)簡(jiǎn)化為一維導(dǎo)熱方程處理。

下面針對(duì)電熱絲、絕緣層和不銹鋼套管，分別在加熱棒穩(wěn)態(tài)狀態(tài)下對(duì)式(1)進(jìn)行簡(jiǎn)化并求解，得出相應(yīng)材料的穩(wěn)態(tài)溫度分布。假定電熱絲、絕緣層和不銹鋼套管的導(dǎo)熱系數(shù)均為常數(shù)，分別用λ1,λ2,λ3表示，各層外表面溫度分別用t1,t2,t3表示。

1）對(duì)于電熱絲，其穩(wěn)態(tài)導(dǎo)熱微分方程可簡(jiǎn)化為

對(duì)式(2)進(jìn)行積分求解，可得電熱絲的溫度分布表達(dá)式為

2）對(duì)于絕緣層，其穩(wěn)態(tài)導(dǎo)熱微分方程可簡(jiǎn)化為

對(duì)式(4)進(jìn)行積分求解，可得絕緣層的溫度分布表達(dá)式為

3）對(duì)于不銹鋼套管，其穩(wěn)態(tài)導(dǎo)熱微分方程可簡(jiǎn)化為

對(duì)式(6)進(jìn)行積分求解，可得不銹鋼套管的溫度分布表達(dá)式為

參考多層圓柱模型，紅外加熱棒的導(dǎo)熱熱流可表示為（設(shè)加熱棒長(zhǎng)度為l）

在穩(wěn)態(tài)情況下，通過式(3)、(5)、(7)可對(duì)紅外加熱棒內(nèi)部溫度分布進(jìn)行定性分析。加熱棒內(nèi)部電熱絲有最高工作溫度限制，因此實(shí)際使用時(shí)對(duì)加熱棒表面的最高功率密度會(huì)有限制要求[3]。在航天器真空熱試驗(yàn)中使用時(shí)應(yīng)對(duì)加熱棒外表面溫度進(jìn)行監(jiān)測(cè)，避免由于溫度過高而損壞加熱棒。

1.3 應(yīng)用優(yōu)勢(shì)

單位長(zhǎng)度紅外加熱棒的阻值可以按需要設(shè)計(jì)，與傳統(tǒng)鎳鉻加熱帶條相比，在相同覆蓋率條件下，在同一加熱區(qū)域內(nèi)使用紅外加熱棒可以更加方便地變更設(shè)計(jì)阻值，從而提高外熱流模擬器的加熱靈活性。對(duì)于高溫端溫度偏高、低溫端溫度偏低，且比熱容大、降溫困難的航天器部組件熱試驗(yàn)，要求外熱流模擬器不僅覆蓋率小，同時(shí)還須具有足夠的加熱能力。而紅外加熱棒的加熱靈活性使其在航天器真空熱試驗(yàn)中具有明顯的應(yīng)用優(yōu)勢(shì)。

2 外熱流模擬效果分析

為對(duì)紅外加熱棒式和加熱帶式外熱流模擬器 在真空熱試驗(yàn)條件下的外熱流模擬效果進(jìn)行比較分析，針對(duì)這2 種外熱流模擬器分別建立了虛擬試驗(yàn)計(jì)算模型，利用虛擬熱試驗(yàn)計(jì)算平臺(tái)[4]，比較分析了兩者的外熱流模擬效果。

2.1 建模及計(jì)算方法

進(jìn)行表面輻射換熱數(shù)值計(jì)算常用的方法有凈熱量法、網(wǎng)絡(luò)法、區(qū)域法、蒙特卡羅法等[3]。在航天器真空熱試驗(yàn)中，外熱流模擬器與試驗(yàn)件之間參與輻射換熱的表面較多，形狀也常常比較復(fù)雜，采用以上單一理論方法進(jìn)行數(shù)值分析的計(jì)算量很大，求解較為困難[5]。下面采用以蒙特卡羅法為基礎(chǔ)建立熱試驗(yàn)計(jì)算模型來求解，并結(jié)合區(qū)域法難點(diǎn)分離的思想對(duì)計(jì)算過程進(jìn)行簡(jiǎn)化[6]。

蒙特卡羅法進(jìn)行輻射換熱計(jì)算的基本思想是：將一個(gè)表面發(fā)射的輻射能看作是由大量獨(dú)立的能束所組成，并假設(shè)每個(gè)能束具有相同的能量，按照蘭貝特定律的概率分布來描述粒子的輻射方向[7]。實(shí)際計(jì)算的基本步驟為：1）確定發(fā)射點(diǎn)的位置；2）確定發(fā)射能束方向；3）判斷發(fā)射能束的落點(diǎn)；4）統(tǒng)計(jì)到達(dá)吸收表面的能束數(shù)量。

假設(shè)外熱流模擬器上某一加熱棒或加熱帶條表面的發(fā)射點(diǎn)A，在其自身坐標(biāo)系下的坐標(biāo)位置為(x',y',z')，在計(jì)算模型系統(tǒng)坐標(biāo)系下的坐標(biāo)為(x,y,z)，模擬器自身坐標(biāo)的原點(diǎn)在系統(tǒng)坐標(biāo)系中的位置 為(x0′,y0′,z0′ )，則有

式中：α1,α2,α3分別為模擬器自身坐標(biāo)系中的x'軸與系統(tǒng)坐標(biāo)系中x,y,z三個(gè)坐標(biāo)軸的方向夾角。同理，β1,β2,β3和γ1,γ2,γ3分別為y'軸、z'軸與系統(tǒng)坐標(biāo)系中各坐標(biāo)軸的方向夾角。

假設(shè)隨機(jī)能束在模擬器自身坐標(biāo)系中的圓周角為θ，俯仰角為φ，則發(fā)射能束在模擬器自身坐標(biāo)系下的方向向量為(sinθcosφ,sinθsinφ,cosθ)，

利用幾何關(guān)系可計(jì)算出發(fā)射能束在系統(tǒng)坐標(biāo)系下的方向向量為(ax,ay,az)，從而可得到能束發(fā)射方向在系統(tǒng)坐標(biāo)系下的直線方程為[8]

假設(shè)被加熱試驗(yàn)件表面在系統(tǒng)坐標(biāo)系下的曲面方程為f(x,y,z)，通過求解f(x,y,z)和方程(10)可知發(fā)射方向與接收曲面是否存在交點(diǎn)，從而確定發(fā)射能束是否為被加熱表面吸收[9]。當(dāng)發(fā)射的能束足夠多時(shí)，就可以用落在試驗(yàn)件表面的能束數(shù)量來表示試驗(yàn)件表面的熱流密度分布情況。

2.2 虛擬試驗(yàn)計(jì)算模型

對(duì)于紅外加熱棒式和加熱帶式外熱流模擬器的虛擬試驗(yàn)計(jì)算模型如圖2所示，其中外熱流模擬器依據(jù)實(shí)物建模，其外形尺寸為長(zhǎng)0.6 m、寬0.4 m、高0.4 m；試驗(yàn)件為長(zhǎng)400 mm、寬150 mm、厚2 mm的鋁板。2 種外熱流模擬器的加熱棒和加熱帶的覆蓋率均取20%。

圖2 虛擬試驗(yàn)計(jì)算模型 Fig.2 Computational model for virtual thermal test

2.3 計(jì)算結(jié)果及分析

分別在高溫和低溫工況下，對(duì)2 種外熱流模擬器的加熱效果進(jìn)行比對(duì)分析。外熱流模擬器的加熱器表面及試驗(yàn)鋁板表面均按發(fā)黑狀態(tài)處理，表面發(fā)射率取0.9，輻射背景環(huán)境溫度設(shè)置為100 K。

1）高溫工況計(jì)算

在高溫工況下，為保證高溫端鋁板表面到達(dá)熱流密度大于1200 W/m2，2 種外熱流模擬加熱器的總加熱功率均取為800 W，各面輻射熱流密度按總功率平均處理，外熱流模擬器計(jì)算結(jié)果如圖3所示。

圖3 高溫工況計(jì)算結(jié)果 Fig.3 Simulation results in the high-temperature case

可以看出，高溫工況加熱條件下，2 種外熱流 模擬器在試驗(yàn)鋁板表面施加的有效到達(dá)熱流密度分布情況相似，加熱帶式外熱流模擬器的熱流密度略高于加熱棒式的，平均偏高約3.5%。

2）低溫工況計(jì)算

在低溫工況下，鋁板表面到達(dá)熱流密度計(jì)算目標(biāo)值取接近于30～50 W/m2的空間環(huán)境模擬容器背景熱流大小，為此2 種外熱流模擬加熱器的總加熱功率均取為20 W，各面輻射熱流密度按總功率平均處理，外熱流模擬器計(jì)算結(jié)果如圖4所示。

圖4 低溫工況計(jì)算結(jié)果 Fig.4 Simulation results in the low-temperature case

可以看出，低溫工況加熱條件下，2 種外熱流模擬加熱器在試驗(yàn)鋁板表面施加的有效到達(dá)熱流密度分布情況差別不大，同樣是加熱帶式外熱流模擬器略高于加熱棒式的，平均偏高約7.5%。

2.4 小結(jié)

通過虛擬試驗(yàn)計(jì)算分析可以看出，在高低溫工況對(duì)應(yīng)的外熱流條件下，加熱棒式外熱流模擬加熱器和加熱帶式外熱流模擬加熱器對(duì)試件表面形成的有效到達(dá)熱流密度分布方式基本相同，前者可以達(dá)到與后者相似的外熱流模擬效果，通過合理的結(jié)構(gòu)設(shè)計(jì)和熱設(shè)計(jì)，加熱棒式外熱流模擬器的模擬效果可以滿足航天器真空熱試驗(yàn)使用需求。

3 試驗(yàn)驗(yàn)證

3.1 驗(yàn)證試驗(yàn)設(shè)計(jì)

為方便進(jìn)行紅外加熱棒在真空低溫環(huán)境條件下的性能驗(yàn)證，設(shè)計(jì)了各面加熱棒覆蓋率均小于0.1 的長(zhǎng)方體結(jié)構(gòu)外熱流模擬器，并在直徑2 m、長(zhǎng)2.5 m 的空間環(huán)境模擬器內(nèi)進(jìn)行了驗(yàn)證試驗(yàn)。驗(yàn)證試驗(yàn)布局如圖5所示，在外熱流模擬器內(nèi)安裝了一塊鋁板作為控溫對(duì)象。參考航天器大部件真空熱試驗(yàn)的實(shí)際需要，在驗(yàn)證試驗(yàn)中設(shè)計(jì)了單一電流施加和目標(biāo)溫度控制2 種測(cè)試工況。

圖5 驗(yàn)證試驗(yàn)布局示意 Fig.5 Layout of the validation test

3.2 單一電流施加

為充分考核外熱流模擬器在真空低溫環(huán)境中施加不同電流時(shí)的加熱能力，并測(cè)試不同電流條件下加熱棒表面溫度，按照預(yù)設(shè)電流進(jìn)行了單一電流加熱測(cè)試。在試驗(yàn)調(diào)試階段，發(fā)現(xiàn)施加3 A 電流時(shí)，加熱棒表面溫度已經(jīng)接近空間模擬容器溫度測(cè)量系統(tǒng)所用T 型熱電偶溫度高限，因此在正式測(cè)試中將最大測(cè)試電流限定為3 A。

依次從低溫端施加1、2、3 A 的電流，進(jìn)行3個(gè)循環(huán)周期測(cè)試。每個(gè)循環(huán)開始前，首先使鋁板降溫4 h，之后開始施加單一電流，在4 h 內(nèi)完成1個(gè)測(cè)試周期。圖6和圖7分別給出了3 種電流下，紅外加熱棒外表面溫度和溫度變化速率曲線。由圖可以看出，加熱電流越大，其表面溫度達(dá)到平衡的速度越快。表2給出了3 種電流下的加熱棒表面溫度變化的比較數(shù)據(jù)。

圖6 不同電流下紅外加熱棒表面溫度變化曲線 Fig.6 Temperature variations on the infrared heating rod surface

圖7 不同電流下紅外加熱棒表面溫度變化速率 Fig.7 Temperature variations on the infrared heating rod surface

表2 不同電流下紅外加熱棒表面溫度變化 Table2 Temperature variations on the infrared heating rod surface

3.3 目標(biāo)溫度控制

試驗(yàn)驗(yàn)證外熱流模擬器控制既定目標(biāo)溫度的準(zhǔn)確性和穩(wěn)定性，鋁板高溫及低溫端控制目標(biāo)溫度分別設(shè)定為-10 ℃和-80 ℃。在航天器部組件真空熱試驗(yàn)中，低溫端溫度控制難度較大[8]。為更好驗(yàn)證外熱流模擬器的低溫控制效果，延長(zhǎng)了低溫端的驗(yàn)證時(shí)間。圖8為利用PID 控制算法根據(jù)試件溫度變化自動(dòng)調(diào)整加熱器電流實(shí)現(xiàn)的溫度控制曲線，可以看出紅外加熱棒式外熱流模擬器在高低溫端的溫度控制均比較穩(wěn)定，平衡溫度波動(dòng)均小于1 ℃。

圖8 鋁板表面平均溫度變化曲線 Fig.8 Curve of the average temperature on the aluminium surface

4 結(jié)束語

通過對(duì)紅外加熱棒式外熱流模擬器的使用方式、外熱流模擬效果計(jì)算及驗(yàn)證試驗(yàn)分析，認(rèn)為：

1）表面覆蓋率、輻射狀態(tài)、加熱功率相同時(shí)，紅外加熱棒式外熱流模擬加熱器可在試件表面形成與紅外加熱帶式外熱流模擬加熱器接近的到達(dá)熱流分布，兩者高溫端到達(dá)熱流密度相對(duì)平均偏差約為3.5%，低溫端約為7.5%；通過合理的結(jié)構(gòu)設(shè)計(jì)，可以保證外熱流模擬不均勻度優(yōu)于±5%[10]。

2）紅外加熱棒式外熱流模擬器在真空低溫環(huán)境下對(duì)試驗(yàn)件高低溫端溫度的控制準(zhǔn)確性及穩(wěn)定性良好，平衡溫度波動(dòng)均小于1 ℃。

（References）

[1] 楊曉寧,任德鵬,賈陽.紅外加熱棒式空間外熱流模擬器的數(shù)值研究[J].中國(guó)空間科學(xué)技術(shù),2009,29(4): 54-60 Yang Xiaoning,Ren Depeng,Jia Yang.Numerical study on space outer-heat simulator manner with infrared thermal club[J].Chinese Space Science and Technology,2009,29(4): 54-60

[2] 楊世銘,陶文銓.傳熱學(xué)[M].北京: 高等教育出版社,2006: 41-44

[3] 孔寅,姜少飛,李吉泉.電熱變模溫高光模具加熱棒傳熱過程數(shù)值模擬研究[J].輕工機(jī)械.2011,29(6): 34-39 Kong Yin,Jiang Shaofei,Li Jiquan.Numerical value simulation analysis of heat transfer of heat rod in rapid electric heating cycle injection molding[J].Light Industry Machinery,2011,29(6): 34-39

[4] 曹志松,劉紹然,裴一飛.衛(wèi)星虛擬熱試驗(yàn)平臺(tái)建模工具模塊研究[J].航天器環(huán)境工程,2012,29(1): 42-45 Cao Zhisong,Liu Shaoran,Pei Yifei.Tool module for modeling satellite virtual thermal test platform[J].Spacecraft Environment Engineering,2012,29(1): 42-45

[5] 孫玉瑋,楊曉寧,李春楊.圓臺(tái)形紅外籠仿真優(yōu)化研究[J].航天器環(huán)境工程,2011,28(3): 222-227 Sun Yuwei,Yang Xiaoning,Li Chunyang.Simulation study for the optimal design of conical infrared heating cage[J].Spacecraft Environment Engineering,2011,28(3): 222-227

[6] 楊曉寧,孫玉瑋,李春楊.圓柱形紅外加熱籠優(yōu)化設(shè)計(jì)研究[J].航天器環(huán)境工程,2009,26(6): 514-518 Yang Xiaoning,Sun Yuwei,Li Chunyang.The optimization design of cylindrical infrared heating cage[J].Spacecraft Environment Engineering,2009,26(6): 514-518

[7] 劉守文,尹曉芳,裴一飛,等.基于蒙特卡羅方法的紅外燈熱流分布研究[J].宇航學(xué)報(bào),2010,31(2): 608-614 Liu Shouwen,Yin Xiaofang,Pei Yifei,et al.The study of heat flux distribution for infrared lamp based on Monte Carlo method[J].Journal of Astronautics,2010,31(2): 608-614

[8] 尹曉芳,劉守文,郄殿福.紅外燈熱流分布試驗(yàn)研究[J].航天器環(huán)境工程,2010,27(1): 63-66 Yin Xiaofang,Liu Shouwen,Qie Dianfu.Experimental research of the radiant heat flux produced by an infrared lamp[J].Spacecraft Environment Engineering,2010,27(1): 63-66

[9] 楊曉寧,孫玉瑋.利用紅外加熱籠進(jìn)行低熱流模擬的設(shè)計(jì)方法研究[J].航天器環(huán)境工程,2007,24(6): 374-376 Yang Xiaoning,Sun Yuwei.Simulating low flux on spacecraft surface by using infrared heating cage[J].Spacecraft Environment Engineering,2007,24(6): 374-376

[10] 黃本誠(chéng),馬有禮.航天器空間環(huán)境試驗(yàn)技術(shù)[M].北京: 國(guó)防工業(yè)出版社,2002: 62-63