背壁熱解效應(yīng)對喉襯組件溫度場的影響*

閆寶任,史宏斌,唐 敏,李書良

(西安航天動力技術(shù)研究所,西安 710025)

背壁熱解效應(yīng)對喉襯組件溫度場的影響*

閆寶任,史宏斌,唐 敏,李書良

(西安航天動力技術(shù)研究所,西安 710025)

文中以固體火箭發(fā)動機復(fù)合噴管為研究對象,測試了背壁材料的熱解響應(yīng)特性及隨溫度及壓強變化的接觸熱阻,用有限元法進行了二維軸對稱噴管瞬態(tài)熱分析。然后給出考慮與不考慮背壁熱解效應(yīng)情況下的溫度場計算結(jié)果,并進行了驗證試驗。最終將計算結(jié)果與實驗結(jié)果進行比較,結(jié)果表明,背壁熱解效應(yīng)可以降低喉襯組件的溫度,考慮背壁熱解效應(yīng)的喉襯組件溫度場分析更符合實驗值。

喉襯組件;背壁;溫度場;熱解

0 引言

喉襯組件是維持固體火箭發(fā)動機燃燒室預(yù)定工作壓強,產(chǎn)生推力的關(guān)鍵功能部件之一,其工作環(huán)境惡劣,對發(fā)動機的性能影響明顯,喉襯破壞將直接影響發(fā)動機正常工作;背壁延緩喉襯溫度向外部傳遞的絕熱部件,是確保噴管殼體溫度在許用范圍的重要措施[1-2]。

目前,國內(nèi)外在喉襯溫度場分析方面已經(jīng)取得了較大的成績,國內(nèi),孫冰[3]用有限元法計算了邊界移動的噴管溫度場。亢麗娟[4]利用有限元軟件對復(fù)合噴管全域瞬態(tài)溫度場進行了有限元計算分析。田四朋[5]等人通過計算溫度場和應(yīng)力場,分析了噴管在發(fā)動機工作過程中的結(jié)構(gòu)完整性。國外,Kearney W.J.[6]、JohnW.Edwards[7]、Maw,J.F.[8]和Alhama F[9]等人均對噴管熱分析做過深入研究,特別是JohnW.Edwards對ASRM噴管的傳熱分析進行了深入細(xì)致的工作。R.RameshKumar[10]等人針對噴管的復(fù)合材料結(jié)構(gòu)進行了熱力耦合分析。總體上看,受固體發(fā)動機熱流環(huán)境、邊界條件、材料參數(shù)和熱反應(yīng)模型等復(fù)雜因素的制約,國內(nèi)針對熱載荷下喉襯組件的溫度場計算,誤差達(dá)到30%以上。

綜上可知國內(nèi)外文獻中關(guān)于背壁熱解效應(yīng)對喉襯組件溫度場的影響分析較少,因此文中在接觸熱阻隨溫度及壓強變化的基礎(chǔ)上建立了含熱解喉襯組件的傳熱模型,通過ABAQUS軟件,建立一套工程可用的數(shù)值模擬方法,并進行實驗驗證,對比喉襯組件的溫度場計算結(jié)果與實驗測試結(jié)果,驗證算法精度,研究更為準(zhǔn)確的喉襯組件邊界仿真方法。

1 含熱解喉襯組件的傳熱模型

圖1為喉襯組件燒蝕和傳熱模型,根據(jù)Fourier定律,一維圓柱坐標(biāo)多層復(fù)合結(jié)構(gòu)的非穩(wěn)態(tài)傳熱方程如式(1)所示,多數(shù)文獻中將熱解氣體流過時所帶走的熱量Q1,i忽略,使得溫度場計算誤差較大,文中主要考慮Q1,i對喉襯組件溫度場的影響:

(1)

式中:下標(biāo)i=0,1,2,…,5分別對應(yīng)喉襯燒蝕層、喉襯原始層、背壁碳化層、背壁熱解層(或熱界面)、原始材料層和金屬殼體層;k為導(dǎo)熱系數(shù);ρ為密度;C為比熱;Q1為熱解氣體流過時所帶走(輸入)的熱量;Q2為熱解潛熱。

圖1 喉襯組件燒蝕和傳熱模型

2 理論模型和邊界條件

2.1 理論模型

喉襯組件的軸對稱結(jié)構(gòu)如圖2,喉襯采用石墨材料,背壁采用5-Ⅱ石棉酚醛模壓材料,收斂段采用高硅氧/酚醛材料,擴張段采用C/C材料,金屬殼體采用鋼材料。T1、T2分別為距喉襯內(nèi)邊界20 mm、34 mm計算溫度點。

圖2 喉襯組件的軸對稱結(jié)構(gòu)

2.2 流場特性

由傳熱學(xué)知識可知,發(fā)動機點火后,噴管壁內(nèi)各點的溫度實際上是時間和空間坐標(biāo)的函數(shù)。但由于噴管的幾何形狀是軸對稱的,受載荷邊界條件也是軸對稱的。考慮到所用材料的特點,因材料本身的缺陷所產(chǎn)生的周向熱流是很小的,可以忽略不計。這樣,可將實際上的三維空間問題,簡化為軸對稱二維空間不穩(wěn)定導(dǎo)熱問題[11]。并在分析求解時作如下假設(shè):

1)發(fā)動機穩(wěn)態(tài)工作,燃?xì)鈪?shù)(壓強和溫度等)不隨時間變化;

2)進行溫度計算時,不計噴管內(nèi)壁燒蝕;

3)不考慮輻射傳熱和顆粒接觸傳熱;

4)計算中考慮隨溫度及壓強變化的接觸熱阻。

2.3 燃?xì)鈪?shù)

為便于進行溫度分析,假定燃?xì)饬鲃邮欠€(wěn)態(tài)的;燃燒產(chǎn)物是組分均勻的完全氣體;流動是等熵的。實踐證明,在噴管型面選定后,采用一維等熵流分析噴管流場即可滿足要求。燃?xì)鈱α鲹Q熱可由式(2)表示:

(2)

式中:q表示熱流(W/m2);hg表示對流換熱系數(shù),可由巴茲公式確定(W/(m2·K));T?w表示燃?xì)鉁囟?Twg表示噴管內(nèi)壁溫度。巴茲公式可寫為:

(3)

式中:dt為噴管喉部直徑(m);μ為燃?xì)怵ば韵禂?shù)(kg·s/m2);cpg為燃?xì)舛▔罕葻崛?kcal/(kg·K));Pr為燃?xì)馄绽侍財?shù);c*為燃?xì)馓卣魉俣?m/s);rc為噴管喉部曲率半徑(m);At為噴管喉部面積(m2);A為噴管擴張段某計算截面的面積(m2);σ1為邊界層修正系數(shù)。

文中,取滯止溫度T0為2 751 ℃,可獲得溫度、對流換熱系數(shù)隨噴管軸向的變化情況,如圖3所示。其中,h0為噴管喉部的燃?xì)鈱α鲹Q熱系數(shù)。

圖3 溫度、對流換熱系數(shù)隨噴管軸向變化曲線

2.4 接觸熱阻試驗

本實驗所使用的試件為長方體,邊長為30 mm×30 mm×90 mm。加工M8×0.75的螺紋洞,深度為20 mm。各個孔心保持在一條直線上。本實驗由兩部分組成,一是碳/碳復(fù)合材料和酚醛樹脂的溫度與界面接觸熱阻的關(guān)系試驗,另一組是界面接觸應(yīng)力與界面接觸熱阻的試驗。表1給出了通過靜態(tài)熱流法測得的接觸熱阻隨溫度及壓強變化的實驗數(shù)據(jù)。

表1 接觸熱阻隨溫度及壓強變化實驗數(shù)據(jù)

2.5 5-Ⅱ材料的熱解響應(yīng)特性

圖4為采用NETZSCH DSC404 F3差示掃描量熱儀進行的潛熱測試結(jié)果。試驗在氬氣保護下進行,升溫速率為20 ℃/min,試樣質(zhì)量3.7 mg,采用Al2O3坩堝盛放試樣,最高溫度升至900 ℃。

測試結(jié)果顯示,與熱失重分析的結(jié)果類似,熱解反應(yīng)也分3個階段:

第一階段為放熱反應(yīng),反應(yīng)開始的溫度在355 ℃左右,結(jié)束溫度在390 ℃左右,放出的熱量為19.62 J/g;

第二階段反應(yīng)為吸熱反應(yīng),反應(yīng)開始的溫度在390 ℃左右,結(jié)束溫度在690 ℃左右,吸收的熱量為1 214 J/g;

第三階段反應(yīng)為放熱反應(yīng),反應(yīng)開始的溫度在785 ℃左右,結(jié)束溫度在845 ℃左右,放出的熱量為109.3 J/g。

測試材料受熱時有3個明顯的反應(yīng)區(qū)間,第一次反應(yīng)發(fā)生在400 ℃附近,第2次發(fā)生在600 ℃附近,第3次發(fā)生在800 ℃附近。

圖4 石棉酚醛的DSC曲線

通過ABAQUS STANDER求解器HEAT TRANSFER分析步對喉襯組件進行熱傳導(dǎo)分析,分析共分兩步,分析時間分別為8 s和90 s,用于模擬發(fā)動機的工作過程和其后的自然冷卻過程。噴管的初始溫度為室溫(20 ℃);外界大氣壓為一個標(biāo)準(zhǔn)大氣壓,溫度為室溫;外壁與空氣進行自然對流換熱,對流換熱系數(shù)為5 W/(m2·K),接觸熱阻隨溫度及壓強變化如表一,計算時考慮背壁材料的潛熱,5-Ⅱ材料的熱解主要發(fā)生在390 ℃～690 ℃,潛熱為1 214 J/g。

3 計算結(jié)果及分析

由圖5可以看出,T1點處加潛熱與無潛熱計算溫度隨時間變化趨勢相同,25 s前T1點隨時間增加溫度快速升高,到達(dá)25 s后,溫度變化較小,趨于平緩,55 s時加潛熱與無潛熱溫度達(dá)到最高值,加潛熱最高溫度為712.01 ℃,無潛熱最高溫度為811.08 ℃。無潛熱溫度值高于加潛熱溫度值主要是由背壁中的酚醛成分熱解吸收熱量所致。

由圖6可以看出,T2點處無潛熱計算溫度隨時間增加迅速升高,20 s時溫度達(dá)到最高,最高溫度為499.1 ℃,隨后溫度降低,但變化較為平緩,T2點處加潛熱計算溫度隨時間增加溫度快速升高,到達(dá)20 s時溫度達(dá)到最高,最高溫度為415.51 ℃,隨后溫度突然降低,直至30 s時溫度繼續(xù)降低,但變化較為平緩,30 s后T2點處加潛熱與無潛熱溫度變化趨勢相同。造成20 s時T2點處加潛熱溫度突然降低的原因為背壁材料中大量酚醛樹脂發(fā)生分解反應(yīng),分子鏈裂解產(chǎn)生H2O、CH4、H2、CO和CO2吸熱所致。

圖5 T1處計算結(jié)果

圖6 T2處計算結(jié)果

4 驗證試驗

本實驗在Φ340標(biāo)準(zhǔn)試驗發(fā)動機上進行,發(fā)動機工作時間約為8 s,滯止壓強為10.6 MPa。采用在喉襯上打孔安裝熱電偶的方式來測試噴管溫度,分別在距喉襯內(nèi)邊界20 mm、34 mm處插入兩個熱電偶T1、T2如圖7,測試出T1、T2兩點溫度隨時間變化規(guī)律。試驗中用鉑銠熱電偶進行溫度測試,采用了螺紋卡套連接結(jié)構(gòu)固定熱電偶,使熱電偶與測試部位緊密接觸,減弱了發(fā)動機振動帶來的影響。測試所用熱電偶如圖8所示,測溫孔與熱電偶連接螺紋的密封采用高溫氧化銅膠。

圖7 喉襯、背壁結(jié)構(gòu)及熱電偶布置

圖8 測試用熱電偶及螺紋卡套結(jié)構(gòu)

試驗測試結(jié)果如圖9,時間零點為發(fā)動機點火時刻,發(fā)動機工作時間為8 s,發(fā)動機停止工作后讓其自然對流換熱82 s,最終采集了90 s時長的溫度數(shù)據(jù)。B測點深,靠近喉襯內(nèi)表面,溫升快,在發(fā)動機點火后溫度迅速升高,最高溫度639.68 ℃。由于A測點位置較淺,溫升速率降低,最高溫度出現(xiàn)的時間為20 s,達(dá)到380.77 ℃。

圖9 T1、T2兩點試驗測溫曲線

5 試驗結(jié)果與計算結(jié)果的對比分析

通過計算結(jié)果與實驗結(jié)果對比曲線圖10和圖11可知,T1點加潛熱與無潛熱溫度計算趨勢與試驗實測趨勢相同,加潛熱計算值比無潛熱計算值更符合試驗值。T2點加潛熱溫度計算趨勢與試驗實測趨勢相同,但T2點無潛熱溫度計算值在20 s后趨于平緩,與試驗實測溫度變化曲線差異較明顯,加潛熱計算值比無潛熱計算值更符合試驗值。

圖10 T1處計算結(jié)果與實驗結(jié)果對比圖

圖11 T2處計算結(jié)果與實驗結(jié)果對比圖

綜上可知,T1點處加潛熱最高溫度與實際測量最高溫度相對誤差為11.31%,T2點處加潛熱最高溫度與實際測量最高溫度相對誤差為9.12%。理論計算值高于實驗值的原因:物性參數(shù)的選取會影響計算結(jié)果;實際測量中,熱電偶的慣性會使測量值略低于結(jié)構(gòu)的實際溫度值。

6 結(jié)論

文中基于二維有限元計算模型,分析了喉襯組件的瞬態(tài)溫度場,考慮了加潛熱與不加潛熱喉襯組件溫度變化規(guī)律,得出以下結(jié)論:

1)背壁熱解吸熱減輕了溫度向噴管外壁的傳遞,降低了喉襯組件的溫度,有利于熱防護。

2)加潛熱的溫度場理論計算結(jié)果與實驗結(jié)果的比較說明,加潛熱的計算模型更符合實際測量,結(jié)果更可信。

3)背壁材料的熱解對溫度場的影響是顯著的,計算時需要考慮熱解對喉襯組件溫度場的影響。

[1] 陳汝訓(xùn). 固體火箭發(fā)動機設(shè)計與研究 [M]. 北京: 宇航出版社, 1991: 1-2.

[2] 王錚, 胡永強. 固體火箭發(fā)動機 [M]. 北京: 宇航出版社, 1993: 238-239.

[3] 孫冰, 孫菊芳. 用有限元法計算邊界移動的噴管溫度場 [J]. 推進技術(shù), 1999(5): 54-58.

[4] 亢麗娟. 噴管整體溫度場有限元分析研究 [C]∥中國宇航學(xué)會. 2005年中國固體火箭推進第22屆年會論文集(發(fā)動機分冊), 2005: 141-144.

[5] 田四朋, 唐國金, 李道奎, 等. 固體火箭發(fā)動機噴管結(jié)構(gòu)完整性分析 [J]. 固體火箭技術(shù), 2005, 28(3): 180-183.

[6] Kearney W J. Advanced solid rocket motor nozzle development status, AIAA 1993-2596 [R]. 1993.

[7] John W. Edwards. ASRM nozzle design and develo pment, AIAA 1991-2070 [R]. 1991.

[8] Maw J F, Buttars R L. Thermal analysis of ASRM nozzle, TWR 62159 RerA [R]. 1982.

[9] Alhama F, Campo A. Network simulation of the rapid temperature changes in the composite nozzle wall of an experimental rocket engine during a ground firing test [J]. Applied Thermal Engineering, 2003, 23: 37-47.

[10] Ramesh Kumar G, Vinod S, Renjith G, et al. Harikrishnan. Thermo-structural analysis of composite structures [J]. Materials Science and Engineering, 2005, 412(1): 66-70.

[11] 胡江華. 錐形套式連接C/C噴管擴張段溫度場與應(yīng)力場分析 [J]. 固體火箭技術(shù), 2012, 35(1): 64-68.

The Influence of Pyrolytic Effect of the Back Surface on the Temperature Field of Throat Lining Component

YAN Baoren,SHI Hongbin,TANG Min,LI Shuliang

(Xi’an Institute of Aerospace Propulsion Technology, Xi’an 710025, China)

In this paper, composite nozzle of solid rocket motor was considered as research object, pyrolysis characteristics of back surface material and contact thermal resistance changing with temperature and pressure were tested, finite element method was used for transient thermal analysis of two-dimensional axisymmetric nozzle. The pyrolytic effect of the back surface was considered, then calculation results of the temperature field with expected and unexpected pyrolytic effect of the back surface were given, and validation test was done. Finally, the calculated results were compared with the experimental results, the results show that, pyrolysis effect of the back surface can reduce temperature of the throat lining components, and temperature field analysis with consideration of pyrolysis effect of the back surface of the throat lining components accords more with the experimental value.

throat lining component; back surface; temperature field; pyrolysis

2014-07-29

閆寶任(1989-),男,黑龍江人,碩士研究生,研究方向:固體火箭發(fā)動機工程力學(xué)。

V435.+14

A