液體火箭發(fā)動機(jī)組件熱真空虛擬試驗技術(shù)

張魏靜，劉占一，劉計武，石曉波，胡錦華，程亞威

(西安航天動力研究所 液體火箭發(fā)動機(jī)技術(shù)重點實驗室，陜西 西安 710100)

0 引言

航天器在軌運行時所處的真空熱環(huán)境復(fù)雜多變，為了保證航天器上各個組件能夠承受這種環(huán)境的影響，并在整個運行過程中正常工作，需要通過熱試驗來考核航天器的熱結(jié)構(gòu)設(shè)計、驗證組件的熱特性。對國外衛(wèi)星地面環(huán)境試驗的故障進(jìn)行統(tǒng)計后得到，在187顆衛(wèi)星故障事例中，熱真空試驗暴露出的故障占47%，可見熱真空試驗對于評判航天器的工作性能、暴露設(shè)計與工藝問題、提升產(chǎn)品可靠度的重要性。此前，熱真空試驗主要針對衛(wèi)星、飛船、空間站等航天器及其組件，但實際上，火箭發(fā)動機(jī)及其零部組件同樣需要進(jìn)行熱真空試驗，尤其是開展月球探測以及深空探測所用的軌控、姿控發(fā)動機(jī)在軌運行時間越來越長，更需要通過發(fā)動機(jī)熱真空試驗驗證其熱設(shè)計的可靠性。熱真空試驗對考核驗證發(fā)動機(jī)組件的性能起到了重要作用，但也存在一些局限性:①熱真空試驗需要的周期較長且花費較高；②試驗溫度測點有限，傳感器布置位置受試件表面形狀和大小限制，通常無法對零部件內(nèi)部的關(guān)注點進(jìn)行溫度測量。

近些年來，隨著數(shù)值模擬技術(shù)的快速發(fā)展，熱真空虛擬試驗已經(jīng)成為可能。目前，國外已經(jīng)廣泛采用熱真空虛擬試驗技術(shù)進(jìn)行產(chǎn)品設(shè)計研究。美國NASA開展研究構(gòu)建全系統(tǒng)虛擬試驗平臺，力圖實現(xiàn)運載火箭系統(tǒng)、發(fā)射系統(tǒng)、回收系統(tǒng)、生保系統(tǒng)等同步設(shè)計。

我國也在逐步開展熱真空虛擬試驗技術(shù)研究，該技術(shù)在衛(wèi)星領(lǐng)域得到較多應(yīng)用。衛(wèi)星虛擬試驗是通過構(gòu)建虛擬熱試驗平臺，并建立虛擬熱試驗環(huán)境，然后利用衛(wèi)星熱物理數(shù)值仿真模型開展試驗，從而獲得衛(wèi)星在外熱流作用下的瞬態(tài)溫度場分布。北京衛(wèi)星環(huán)境工程研究所研發(fā)建立的衛(wèi)星虛擬熱試驗平臺，具備建模工具模塊、加熱系統(tǒng)模塊和熱沉系統(tǒng)模塊等多個相互獨立的模塊，能夠快速直觀的進(jìn)行虛擬試驗。北京強(qiáng)度環(huán)境研究所建立了用于開展飛行器結(jié)構(gòu)熱試驗的虛擬熱試驗平臺，平臺具有結(jié)構(gòu)熱分析、試驗邊界分析等功能，同時具備試驗控制仿真以及電功率系統(tǒng)數(shù)字仿真功能。李濤等利用NX—TMG軟件對某天線熱真空試驗過程進(jìn)行模擬，獲得天線的詳細(xì)溫度場分布，并發(fā)現(xiàn)影響天線溫度均勻性的關(guān)鍵因素。紀(jì)欣言等結(jié)合加熱板在航天器熱試驗中的傳熱特點，在軟件Thermal Desktop中建立數(shù)值仿真模型，分析了加熱板材料、電熱絲數(shù)量和熱沉管路布置方式等因素，對加熱板升溫時間和降溫時間以及溫度分布的影響，建立了模擬星、加熱板和熱沉的仿真模型，分析了加熱板和星表間距對模擬星表面溫度的影響。仇善昌等針對數(shù)字式太陽敏感器，建立了熱真空試驗的仿真模型，開展了太陽敏感器溫度場仿真計算，仿真結(jié)果與真實熱真空試驗數(shù)據(jù)對比表明兩者偏差小于3 K。

國內(nèi)通過開展上述研究工作，積累了較多關(guān)于衛(wèi)星的熱真空虛擬試驗經(jīng)驗。然而，液體火箭發(fā)動機(jī)在飛行過程中，其組件面臨的熱環(huán)境包括內(nèi)熱與外熱兩類：內(nèi)熱主要有發(fā)動機(jī)高溫噴管壁面輻射熱、發(fā)生器壁面輻射熱、渦輪殼體輻射熱，以及推力室出口噴流加熱等；外熱主要有氣動熱、真空、太陽輻射、冷黑背景等復(fù)雜的熱真空環(huán)境。具體來看，發(fā)動機(jī)高溫壁面的輻射熱在工作中一直存在，其他熱流隨飛行高度變化:一級發(fā)動機(jī)工作高度0～80 km，對應(yīng)環(huán)境壓力100 kPa～1 Pa，飛行過程中，發(fā)動機(jī)的某些零部組件可能受到發(fā)動機(jī)噴流加熱;二級及以上的高空發(fā)動機(jī)(不包括軌姿控發(fā)動機(jī))工作高度較高，環(huán)境壓力更接近真空，發(fā)動機(jī)羽流覆蓋面更大，發(fā)動機(jī)零部組件也可能受到噴流加熱，另外由于工作時間可達(dá)數(shù)十分鐘，太陽輻射也不能忽略。軌姿控發(fā)動機(jī)一般工作在高真空環(huán)境，環(huán)境壓力低于1.0×10Pa，而且面臨冷黑背景(4 K)，工作時間達(dá)數(shù)小時至數(shù)天、甚至更長，由于與太陽相對位置的變化，冷熱溫度交變是其面臨的主要熱環(huán)境。

為了保證各類發(fā)動機(jī)的可靠工作，在發(fā)動機(jī)研制過程中，需要在地面營造相應(yīng)的熱真空環(huán)境，并模擬熱流變化過程，對其零部組件甚至整機(jī)進(jìn)行真空熱試驗。

本文基于蒙特卡洛光學(xué)追蹤法進(jìn)行了輻射換熱計算，采用Sinda-fluint熱分析軟件建立了由真空艙、熱沉、石英燈陣和具體發(fā)動機(jī)組件構(gòu)成的熱真空虛擬試驗平臺。

1 發(fā)動機(jī)組件熱真空試驗

1.1 熱真空試驗臺

液體火箭發(fā)動機(jī)組件熱真空試驗一般需要營造三方面的環(huán)境:真空環(huán)境、冷黑背景和熱流模擬(包括太陽輻射和發(fā)動機(jī)自身熱源)。從傳熱角度來看，真空環(huán)境下發(fā)動機(jī)與外界的傳熱主要通過輻射方式進(jìn)行，發(fā)動機(jī)零部組件的表面輻射特性對發(fā)動機(jī)熱狀態(tài)起著重要作用；另外發(fā)動機(jī)中靜態(tài)接觸的部件(如熱控制、熱防護(hù)包覆材料等)，由于材料之間的真空空隙使得接觸熱阻增大，造成地面大氣環(huán)境試驗結(jié)果與真空試驗結(jié)果偏差大。冷黑背景是指宇宙空間相當(dāng)于溫度4 K的黑體，是一個巨大的吸熱裝置，又稱熱沉，這與常溫的地面環(huán)境也有很大差異，在進(jìn)行發(fā)動機(jī)熱設(shè)計時也必須考慮。

某熱真空試驗臺實物如圖1和圖2所示。該實驗臺由試驗艙、熱沉、真空抽氣機(jī)組、大熱流密度加熱系統(tǒng)、液氮系統(tǒng)、控制和數(shù)據(jù)采集處理系統(tǒng)等組成，系統(tǒng)圖如3所示。其中真空艙有效容積12 m，艙內(nèi)壓力低于1.0×10Pa，艙內(nèi)四周的液氮熱沉溫度可低至100 K，艙內(nèi)上部的石英燈加熱器最大加熱功率為100 kW。

圖1 試驗臺外觀照片F(xiàn)ig.1 Appearance photo of test bed

圖2 試驗臺艙內(nèi)照片F(xiàn)ig.2 Interior photo of test bed

圖3 熱真空試驗臺系統(tǒng)圖Fig.3 System diagram of thermalvacuum test bed

1.2 熱真空虛擬試驗

本文以熱真空試驗臺為對象建立虛擬試驗平臺，并開展對照測試對虛擬試驗平臺進(jìn)行驗證，達(dá)到提高發(fā)動機(jī)組件熱真空試驗效率的目的，同時也可以獲取發(fā)動機(jī)組件更詳細(xì)的溫度場分布。虛擬試驗平臺幾何模型如圖4所示，主要包含試驗艙壁、熱沉、石英燈陣及底踏板，下文將對虛擬試驗平臺進(jìn)行詳細(xì)介紹。

圖4 試驗臺模型Fig.4 Model of test bed

1.2.1 數(shù)值分析與模型

本文采用了航天器常用的熱網(wǎng)絡(luò)分析方法進(jìn)行虛擬試驗仿真計算。該方法基于熱傳遞過程與導(dǎo)電過程的相似性，采取離散化的數(shù)學(xué)處理，通過有限差分法求解計算。計算時將發(fā)動機(jī)組件劃分為若干單元，單元之間的換熱過程(包括輻射、傳導(dǎo)和對流)可歸結(jié)為節(jié)點之間的熱流傳遞過程，各節(jié)點通過不同類型的熱阻連接。依照這種方式形成的熱網(wǎng)絡(luò)與電網(wǎng)絡(luò)具有相似性，通過兩者在圖形及數(shù)學(xué)模型上的類比關(guān)系，借用電學(xué)上的基爾霍夫第二定律即可得到各個節(jié)點的熱平衡方程，然后求解得出發(fā)動機(jī)各點的溫度及其變化率。在試驗系統(tǒng)內(nèi)，液體火箭發(fā)動機(jī)組件任意節(jié)點

j

的熱平衡方程為

Q

+

Q

+

Q

+

Q

=

Q

+

Q

(1)

其中

Q

=

ε

E

φ

F

式中:

Q

為節(jié)點

j

表面從熱源吸收的外熱流，W；

Q

為試驗系統(tǒng)內(nèi)各表面發(fā)射的能量被節(jié)點

j

吸收的部分，W；

Q

為節(jié)點

j

內(nèi)熱源，W；

Q

為節(jié)點

j

與其他節(jié)點的導(dǎo)熱傳熱量，W；

Q

為節(jié)點

j

內(nèi)能的變化值，W；

Q

為節(jié)點

j

表面向內(nèi)、外輻射的熱量，W；

ε

為節(jié)點

j

表面發(fā)射率，無量綱；

E

為熱源的平均輻射密度，W/m；

φ

為節(jié)點

j

相對于熱源的幾何角系數(shù)，無量綱；

F

為節(jié)點

j

的外表面積，m；

R

為節(jié)點

i

和節(jié)點

j

之間的熱輻射網(wǎng)絡(luò)傳遞眾數(shù)，m;

σ

為斯忒藩—玻耳茲曼常數(shù)，5.67×10W/(m·K)；

T

為節(jié)點

i

的溫度，K；

m

為與節(jié)點

j

有輻射換熱的節(jié)點數(shù)；

D

為節(jié)點

i

與

j

之間的熱傳導(dǎo)系數(shù)，W/K；

T

為節(jié)點

j

的溫度，K；

n

為與節(jié)點

j

有傳導(dǎo)換熱的節(jié)點數(shù)；

G

為節(jié)點

j

的質(zhì)量，kg；

c

為節(jié)點

j

的比熱容，J/(kg·K)；

τ

為時間，s。

在液體火箭發(fā)動機(jī)熱真空試驗中，存在大量輻射換熱過程，輻射換熱通常會耗費較多時間、并占用大量存儲空間，因而實現(xiàn)準(zhǔn)確高效求解計算對于虛擬試驗平臺來說十分重要。蒙特卡洛法是求解輻射傳熱問題的有效方法，作為一種概率模擬方法，該方法可將輻射熱的傳輸過程分解為發(fā)射、反射、透射、吸收和散射等子過程，把其轉(zhuǎn)化為隨機(jī)問題。令每個單元發(fā)射一定的光束，跟蹤、統(tǒng)計每束光束的歸宿，從而得到該單元輻射能量分配的統(tǒng)計結(jié)果。本文虛擬試驗中的輻射換熱關(guān)系即通過蒙特卡洛法求解。

根據(jù)研究需求與思路，選擇Sinda-fluint的Thermal Desktop模塊作為建模仿真平臺,并基于蒙特卡洛光學(xué)追蹤法進(jìn)行輻射計算，可進(jìn)行各類表面無限次光線反射的模擬，除了可以模擬復(fù)雜工程系統(tǒng)內(nèi)存在的熱輻射傳熱之外，還可利用節(jié)點網(wǎng)絡(luò)法模擬導(dǎo)熱傳熱、流固耦合等。其自帶的前處理模塊，可對導(dǎo)入模型自動進(jìn)行簡化與修復(fù)，幾何模型的修改可自動更新到熱模型中，這些特性方便了仿真模型的建立與修正。

1.2.2 幾何模型

虛擬試驗平臺由真空艙、熱沉、石英燈以及試驗件構(gòu)成，其中真空艙、熱沉、石英燈為固定幾何，每次試驗均包括這些幾何結(jié)構(gòu)。試驗件為可更換結(jié)構(gòu)，根據(jù)試驗需求調(diào)整替換。

1.2.2.1 試驗艙和熱沉模型

部件空間試驗臺的試驗艙分為圓柱段艙體、球形封頭和球形艙門3部分，艙底有不銹鋼踏板。低溫?zé)岢劣刹讳P鋼管緊密排列形成，艙體圓柱段的液氮管束貼合于艙內(nèi)壁布置，圓柱段型面與艙體內(nèi)壁型面相一致，球形艙門和尾端封頭為液氮直管束形成的平面熱沉，可實現(xiàn)大熱流加熱狀態(tài)下的低溫冷黑背景。

由于熱沉低溫管路有兩層，管外焊接有導(dǎo)熱率較高的銅翅片，且相鄰的翅片相互交錯，基本保證試驗艙內(nèi)表面被液氮管路覆蓋。因此在搭建虛擬試驗平臺時，不考慮艙壁的熱影響，只建立艙底踏板和熱沉模型，熱沉由中間的圓柱和兩端的圓盤封頭組成。部件空間試驗臺的熱沉由開式液氮循環(huán)形成，艙內(nèi)熱沉管的壁溫維持在100 K，液氮管路外表面涂覆黑漆，發(fā)射率為0.97以上，因此虛擬試驗平臺將熱沉設(shè)置為定溫邊界100 K，黑度取為1。

1.2.2.2 石英燈陣模型

試驗艙艙頂布置有平板石英燈陣，燈陣尺寸為1 m×1 m，通過調(diào)節(jié)電壓和電流加載不同大小的熱流密度，該試驗臺能達(dá)到的最大熱流密度為100 kW/m。石英燈陣以輻射形式向外散熱，受熱部組件接受的輻射熱流與熱源的方向和位置有關(guān)，即石英燈陣投射到試驗件表面的熱流并不是一個定值。為更真實模擬石英燈陣輻射出的熱流分布，虛擬試驗平臺在艙頂建立一塊1 m×1 m的高溫平板，令高溫平板輻射出的熱流等于試驗加載熱流，平板溫度根據(jù)斯忒藩—玻爾茲曼四次方定律確定。試驗加載熱流

q

=

σεT

(2)

式中

T

為輻射體表面溫度，K。

依據(jù)式(2)反算出平板溫度，然后將此溫度作為確定的熱邊界條件，開展熱真空虛擬試驗時，直接計算高溫平板與受熱體之間的輻射換熱量，這樣可以保證虛擬試驗時石英燈陣對受熱體投射的熱流與真實試驗相一致。

為驗證該方法的有效性，開展測試試驗，并建立相應(yīng)的仿真模型進(jìn)行計算。分析對象為一直徑100 mm，厚度2 mm的不銹鋼圓盤，圓盤表面涂覆黑漆，使發(fā)射率接近1，以排除圓盤表面發(fā)射率對測試過程的干擾。試驗圓盤如圖5所示，圖中右側(cè)為未涂覆黑漆的圓盤。

圖5 圓盤照片F(xiàn)ig.5 Round plates

試驗時圓盤吊裝在距離石英燈200 mm的正下方，給石英燈加載1 000 W/m的熱流，持續(xù)600 s，在圓盤下表面中心位置處及邊緣位置處分別布置熱電偶測量溫度。在熱分析軟件中建立石英燈和圓盤試件模型，石英燈溫度根據(jù)前述方法反算得到，進(jìn)行非穩(wěn)態(tài)傳熱計算。對比計算結(jié)果與實測結(jié)果，曲線如圖6所示，可以看出兩者吻合良好。

圖6 實測結(jié)果與仿真結(jié)果對比Fig.6 Result comparison of test and simulation

1.2.2.3 試驗件模型

在建立虛擬試驗仿真模型時，需要對幾何結(jié)構(gòu)進(jìn)行合理簡化。在具體實施過程中，依據(jù)發(fā)動機(jī)實際物理模型，去掉安裝板和發(fā)動機(jī)面板上的小孔，將復(fù)雜的曲線簡化為直線，復(fù)雜的曲面簡化為多個平面的組合，并保證兩者的面積相等、面法向基本保持一致，在保持幾何形體的表面面積和導(dǎo)熱截面積不變的前提下，將不規(guī)則的幾何體簡化為規(guī)則的幾何形體。將試驗組件劃分為若干個單元，單元體的幾何中心稱為節(jié)點，節(jié)點的溫度和熱物性代表了整個單元體內(nèi)的溫度和平均熱物性參數(shù)值。節(jié)點劃分原則為，在能反映試驗組件主要熱特性和滿足工程設(shè)計要求的前提下，盡可能地減少節(jié)點數(shù)量。通過熱網(wǎng)絡(luò)系數(shù)建立各節(jié)點之間的傳熱關(guān)系，對于不同傳熱方式，熱網(wǎng)絡(luò)系數(shù)的表達(dá)式為:

一般的熱傳導(dǎo)

(3)

式中:

A

為有效導(dǎo)熱面積，m；

L

,為節(jié)點

i

與

j

間的有效導(dǎo)熱距離，m；

λ

為熱導(dǎo)率，W/(m·K)。

接觸熱導(dǎo)

D

,=

A

·

H

(4)

式中:

A

為接觸熱導(dǎo)面積，m；

H

為接觸傳熱系數(shù)，W/(m·K)。

輻射傳遞

R

,=

A

ε

B

,

(5)

式中:

A

為節(jié)點

i

的輻射面積，m；

B

為節(jié)點

i

輻射的能量被節(jié)點

j

所吸收的份額(包括多次反射吸收)，稱之為吸收因子，無量綱。

有交界面的兩個節(jié)點間熱傳導(dǎo)

D

=1

/

[

L

/

(

A

λ

)+1

/

(

A

H

)+

L

/

(

A

λ

)]

(6)

2 實例

液體火箭發(fā)動機(jī)流量調(diào)節(jié)裝置的工作性能對溫度較敏感，在工作過程中受到噴管壁面輻射和噴流的熱影響?？己嗽囼炘诓考臻g試驗臺上進(jìn)行，利用石英燈實現(xiàn)熱源模擬。流量調(diào)節(jié)裝置安裝于石英燈陣正下方200 mm處、固定于安裝支架上，安裝時保持流量調(diào)節(jié)裝置與熱源的相對方位與實際產(chǎn)品一致。封閉試驗艙門后開啟真空泵，當(dāng)艙內(nèi)壓力低于10Pa后給石英燈通電加熱，以5 000 W/m的熱流密度持續(xù)加熱1 200 s，斷電停止加熱后繼續(xù)采集200 s數(shù)據(jù)。

在熱真空虛擬試驗平臺上同步開展考核試驗，試驗時產(chǎn)品安裝方位與真實試驗相同，加載的熱流值和作用時間也與真實試驗一致。由于產(chǎn)品與安裝支架的接觸面積很小，虛擬試驗中不考慮支架的熱影響。試驗件仿真模型如圖7所示。試驗結(jié)束后將結(jié)果與真實試驗進(jìn)行對比，結(jié)果如圖8和圖9所示?？梢钥闯?，兩種試驗得到的組件溫度變化趨勢一致、溫度值接近，測點最高溫度值的偏差小于10%。測點3溫度差異稍大是由于虛擬試驗中采用的三維幾何模型與實際產(chǎn)品略有不同，測點3所在位置為流量調(diào)節(jié)器電機(jī)，實際產(chǎn)品的電機(jī)內(nèi)部存在間隙與空腔，而三維模型為實體結(jié)構(gòu)，故三維模型的熱容較實際產(chǎn)品大，溫升速率相應(yīng)較慢。以上結(jié)果與分析表明，虛擬試驗結(jié)果與真實試驗結(jié)果比較吻合，虛擬試驗平臺可以用作發(fā)動機(jī)組件熱考核。

圖7 試驗件仿真模型Fig.7 Simulation model of test specimen

圖8 t=1 200 s時刻試驗件溫度分布云圖Fig.8 Temperature distribution of test specimen at t=1 200 s

圖9 計算結(jié)果與實測結(jié)果對比Fig.9 Result comparison of test and simulation

3 結(jié)論

針對液體火箭發(fā)動機(jī)組件開展熱真空虛擬試驗技術(shù)研究，采用蒙特卡洛光學(xué)追蹤法進(jìn)行輻射換熱計算，基于Sinda-fluint熱分析軟件建立了由真空艙、熱沉、石英燈陣和具體發(fā)動機(jī)組件構(gòu)成的熱真空虛擬試驗平臺，并針對石英燈陣模型進(jìn)行了單獨試驗驗證。以某型發(fā)動機(jī)的流量調(diào)節(jié)裝置為例，同步開展了實際熱真空實驗和虛擬熱真空試驗，得到結(jié)論如下:

1)將石英燈陣簡化為等溫平板模型能夠有效反映熱源特征。

2)在整個測試過程中，虛擬試驗的試驗件溫度和實測值變化趨勢一致，最高溫度偏差在10%以內(nèi)。

3)該熱真空虛擬試驗平臺能夠用于火箭發(fā)動機(jī)組件的熱環(huán)境考核。