牛公杰，錢建平，劉榮忠

（南京理工大學(xué)機械工程學(xué)院，江蘇 南京210094）

引 言

傳統(tǒng)的小型自由裝填式固體火箭發(fā)動機在設(shè)計時，為了控制藥柱的燃燒面積，通常用包覆層將藥柱的部分表面進行包覆［1］。根據(jù)所設(shè)計的裝藥結(jié)構(gòu)幾何尺寸，按照平行層燃燒理論可準(zhǔn)確地計算出燃燒過程中不斷變化的燃面面積，進而計算出火箭發(fā)動機內(nèi)彈道諸元。對于采用星孔裝藥的小型自由裝填式固體火箭發(fā)動機，常用的包覆措施是將藥柱側(cè)面、上、下端面包覆；但也有采取僅將藥柱側(cè)面包覆，而上下端面不包覆的方法，通過藥柱端面和墊片間的間隙使藥柱端面在火箭發(fā)動機工作過程中也參與燃燒，進一步提升火箭發(fā)動機性能［2］。

星孔火箭裝藥的端部初始燃燒面積對火箭發(fā)動機性能有著重要影響。而間隙的大小又決定了端部初始燃面的變化規(guī)律。裝藥端部參與燃燒雖然有助于提升火箭發(fā)動機性能，但也給端部初始燃面的計算帶來困難，無法使用經(jīng)典的燃面計算理論進行描述和計算［3］。本研究通過引入端面燃燒參與度和端面等效燃燒厚度的方法來解決燃面計算問題，并結(jié)合試驗條件推導(dǎo)出了燃面面積經(jīng)驗計算公式，得出溫度對火箭發(fā)動機性能的影響規(guī)律和可用于指導(dǎo)實際工程應(yīng)用的結(jié)論。

1 燃面計算方法

1.1 基本假設(shè)

火箭發(fā)動機裝藥在點火具作用下，藥柱星孔全面瞬時點燃，高溫燃?xì)夂突鹧鏁Z入藥柱端面和墊片之間的間隙并逐漸點燃端面，過程較為復(fù)雜，難以用量化的方式描述，應(yīng)用傳統(tǒng)的燃面計算理論無法描述燃面變化規(guī)律和計算燃面面積，需引入新的方法解決燃面計算問題。

根據(jù)藥柱燃燒特點，可將燃面計算過程分為星孔燃面面積計算和端部燃面面積計算兩部分，并作如下假設(shè)：（1）在點火壓強作用下星孔全面瞬時點燃；（2）忽略藥柱的侵蝕燃燒效應(yīng)；（3）星孔燃燒遵循平行層燃燒定律，燃面沿起始表面各點的法線向內(nèi)部推移。

對于星根圓角化的星孔裝藥，星孔燃燒過程可分為4個階段［4］，即：（1）星根圓角消失前，星邊周邊長為s1，單位：m；（2）星邊消失前，星邊周邊長為s2，單位：m；（3）星邊消失后，星邊周邊長為s3，單位：m；（4）余藥階段燃燒階段，星邊周邊長為s4，單位：m。

圖1為星孔裝藥結(jié)構(gòu)示意圖。其中左側(cè)為星孔藥柱端面燃燒示意圖，右側(cè)為藥柱沿軸向的剖面圖。端部燃燒是個極其復(fù)雜的過程，為簡化計算，作如下4點假設(shè)：（1）藥柱燃燒到t時刻時，一個端面的總面積為Stt，單位：m2，且：Stt=πD2／4-Ap，D為藥柱外徑，單位：m，Ap為藥柱星孔通氣面積，單位：m2；（2）藥柱燃燒到t時刻時，一個端面參與燃燒的面積為Stb，單位：m2，并且假設(shè)燃面Stb平行于藥柱初始端面；（3）引入端面燃燒參與度α，表示一個端面當(dāng)前參與燃燒的面積Stb與當(dāng)前端面總面積Stt之比，即：且0≤α≤1；（4）根據(jù)藥柱質(zhì)量守恒原則，引入端面等效燃燒厚度et，單位：m。

1.2 端面燃燒參與度α

藥柱端面和墊片之間的間隙是影響α變化過程的主要因素。發(fā)動機的零件加工誤差和裝配誤差直接影響藥柱端面與墊片間初始間隙的大小，而溫度變化則會改變初始間隙的大小。一般發(fā)動機殼體的熱膨脹系數(shù)遠(yuǎn)小于藥柱的熱膨脹系數(shù)，而墊片的厚度與發(fā)動機藥室高度和藥柱高度相比很小，故當(dāng)發(fā)動機殼體、藥柱和墊片三者經(jīng)歷相同的溫度變化后，藥柱高度變化量要遠(yuǎn)大于發(fā)動機藥室高度變化量和墊片高度變化量，所以間隙的變化主要取決于藥柱的高度變化量。以發(fā)動機在常溫狀態(tài)下裝配產(chǎn)生的間隙為基準(zhǔn)，當(dāng)發(fā)動機在高溫狀態(tài)下，間隙變小，在低溫狀態(tài)下，間隙變大。

圖1 星孔裝藥結(jié)構(gòu)示意圖Fig.1 Schematic diagram of star shaped rocket grain

從定性的方面分析，藥柱端面與墊片之間的間隙越大、燃燒室壓強p越大、藥柱初始溫度越高、藥柱的長徑比L／D越大，則端面點燃速度越快，端面全面參與燃燒所用的時間越短，反之亦然。因藥柱端面隨著藥柱燃燒時間的增加而逐漸點燃，選取時間t為自變量，考慮到α的變化過程與藥柱燃燒室壓力p、藥柱長度L、初始溫度T、藥柱直徑D等因素密切相關(guān)，α可表示為：

根據(jù)試驗，發(fā)現(xiàn)對于長徑比在1.3左右、特征長度與肉厚比在0.5左右、星角數(shù)為5的星孔裝藥，其端面正常逐漸參與燃燒時，α呈冪函數(shù)形式變化。研究結(jié)合試驗條件和測試結(jié)果，推導(dǎo)出了不同情況下α的經(jīng)驗計算公式。

（1）常溫和低溫情況下，藥柱端面在點火初始時刻即可正常參與燃燒的經(jīng)驗公式為：

式中：t為燃燒時間，s；T為藥柱初始溫度，K；T0為參考基準(zhǔn)溫度，K，T0=273.15K；p0為參考基準(zhǔn)壓力，Pa，p0=106Pa。

（2）高溫情況下，由于藥柱端面與墊片的間隙變小或消失，初期燃燒室內(nèi)壓力相對較小，燃?xì)獠荒茼樌c端面接觸，藥柱端面在初始燃燒階段燃燒參與程度較小，當(dāng)燃燒進行到某一臨界時刻t0（t0＜1），燃燒室內(nèi)壓強升高到某一個臨界值時，燃?xì)饧眲「Z入間隙，迅速全面點燃尚未參加燃燒的端面，α的經(jīng)驗公式為：

式中：t0為臨界時刻，s，其大小根據(jù)具體試驗條件確定。

1.3 端面等效燃燒厚度et

由于端面燃燒參與度系數(shù)α的引入，為保證藥柱燃燒質(zhì)量守恒，需將端面燃燒厚度做等效處理。在端面從起始燃燒到全面燃燒的過程中，即當(dāng)0≤α＜1時，設(shè)在dt時間微元內(nèi)一個端面燃燒掉的質(zhì)量微元為dmt，它應(yīng)該等于當(dāng)前時刻的端面參與燃燒的燃燒面積Stb、燃燒厚度微元de和推進劑密度ρp三者的乘積，即：

式中：mt為端面燃燒掉的質(zhì)量，單位：kg。

將質(zhì)量微元dmt等效由當(dāng)前端面面積為燃面燒去的質(zhì)量，則質(zhì)量微元dmt應(yīng)等于當(dāng)前端面面積Stt、端面等效燃燒厚度微元det和推進劑密度ρp三者的乘積，即：

根據(jù)（4）、（5）兩式相等，可導(dǎo)出端面等效燃燒厚度et的微分方程：

記端面面積全部參與燃燒瞬時時刻為tα1，當(dāng)t≥tα1時，端面面積全部參與燃燒，Stb=Stt，α=1。

式中：a為燃速系數(shù)，m／（s·Pan）；n為燃速指數(shù)。

由（7）式即得到端面等效燃燒厚度計算公式

1.4 總?cè)济嬗嬎?/h3>

總?cè)济婷娣e應(yīng)等于星孔燃面面積加上端部燃面面積。星孔燃面面積根據(jù)平行層燃燒定律計算，藥柱星孔長度應(yīng)等于藥柱初始長度L減去兩個端面的等效燃燒厚度2et，則星孔燃面面積計算公式：

式中：Asb為星孔燃面面積，m2；i=1，2，3，4，對應(yīng)4個燃燒階段。

將Stt=πD2／4-Ap代入上式即得端面燃燒面積計算公式

則總?cè)紵娣e：

式中：Ab為總?cè)济婷娣e，m2。

2 結(jié)果與分析

在相同的室溫條件下完成9 發(fā)火箭發(fā)動機的裝配，放入保溫箱保溫。試驗分常溫、低溫和高溫試驗3組，每組試驗3發(fā)，取3次試驗數(shù)據(jù)的平均值作為每組的試驗結(jié)果。采用零維火箭內(nèi)彈道計算方法［5］，利用計算機采用四階龍哥庫塔法編程求解火箭內(nèi)彈道方程?？紤]到固體火箭發(fā)動機在工作時可能出現(xiàn)兩種端部初始燃面恒定的極端工作狀態(tài)：（1）藥柱端面不參與燃燒，端部初始燃面面積為零；（2）藥柱端面在發(fā)動機工作開始時刻全面瞬時點燃。所以計算時采用3 種燃面計算方法：（1）星孔全面瞬時點燃，上端面、下端面不參與燃燒；（2）星孔、上端面、下端面全面瞬時點燃參與燃燒；（3）星孔全面瞬時點燃、上端面、下端面逐漸點燃參與燃燒。

火箭裝藥為某雙基型推進劑，藥柱密度ρp=1.612g／cm3，火藥力f=929 000J／kg，燃?xì)饨^熱指數(shù)k=1.225，不同溫度條件下的a、n值見表1。

表1 不同溫度條件下的a、n值Table 1 The value of aand nat different temperatures

2.1 常 溫

試驗溫度為+15℃，藥柱初始溫度為288.15K，α采用公式（2）計算。表2為火箭發(fā)動機性能特征量的試驗值和計算值，得到的燃燒室頭部p-t曲線見圖2。

表2 +15℃時發(fā)動機特征量的試驗值與計算值Table 2 Experimental and calculated results at of characteristic values of motor+15℃

圖2 +15℃時燃燒室頭部的p-t曲線Fig.2 The p-t curves of chamber head at+15℃

從表2可以看出，方法1在最大壓強和最大推力方面的相對誤差分別高達(dá)26.97%和27.83%，理論值與試驗值相差甚遠(yuǎn)；而方法2和方法3計算出的各項特征量相對誤差均小于5%。

試驗測試數(shù)據(jù)處理時剔除了試驗測試曲線的初始壓強峰。從圖2可以看出，方法1計算得出的p-t曲線與試驗明顯不符，原因是燃燒過程中只有星孔參與燃燒，星孔的初始燃面面積小，導(dǎo)致初始階段壓強較試驗值小，隨著星孔燃面面積不斷增大，壓強也逐漸增大，壓強呈上升趨勢，與實際情況不符；方法2計算得出的p-t曲線在藥柱初始燃燒階段明顯高于試驗測試曲線，原因是方法2采用的是星孔和兩端面瞬間全部點燃，初始燃面面積較大，致使初始階段燃燒壓力較高；方法3計算得到的p-t曲線僅在0.5～1.0s時間內(nèi)稍低于試驗測試曲線，其余部分與試驗曲線基本吻合，特別是初始燃燒階段吻合較好。

2.2 低 溫

試驗溫度為+15℃，藥柱初始溫度T為233.15K，α采用公式（2）計算。火箭發(fā)動機性能特征量試驗值和計算值見表3，發(fā)動機燃燒室頭部p-t曲線，見圖3。

從表3可以看出，方法1在最大壓強和最大推力方面的相對誤差分別高達(dá)27.53%和28.66%，誤差較大；而方法2和方法3計算出的各項特征量相對誤差均小于5%。

表3 -40℃時發(fā)動機特征量的試驗值與計算值Table 3 Experimental and calculated results of characteristic values of motor at-40℃

圖3 -40℃時燃燒室頭部的p-t曲線Fig.3 The p-t curves of chamber head at-40℃

圖3中，用方法1計算出的p-t曲線與試驗曲線明顯不符；用方法2計算出的p-t曲線在藥柱初始燃燒階段明顯高于試驗測試曲線，而在后期又明顯低于試驗測試曲線；用方法3計算出的p-t曲線在整個燃燒過程中與試驗曲線吻合較好，計算曲線和試驗曲線均比較平滑。

2.3 高 溫

試驗溫度為+50℃，藥柱初始溫度T為328.15K，α采用公式（3）計算，根據(jù)試驗結(jié)果，取t0為0.5s?；鸺l(fā)動機性能特征量試驗測試和理論計算結(jié)果見表4，發(fā)動機燃燒室頭部的p-t曲線見圖4。

表4 +50℃下試驗測試與理論計算結(jié)果Table 4 The experimental test and calculated results of characteristic value of motor at+50℃

圖4 +50℃時燃燒室頭部的p-t曲線Fig.4 The p-t curves of chamber head at+50℃

從表4和圖4可以看出，用方法1和方法2計算出的p-t曲線與試驗曲線相差甚遠(yuǎn)，無法對火箭發(fā)動機性能進行有效預(yù)示，用方法3計算出的曲線在初始階段與試驗吻合較好，但在0.5s以后的曲線上升段和后期吻合不夠理想。在常溫下將火箭藥柱裝入發(fā)動機藥室，令藥柱的一端無約束，在+50℃條件下保溫48h后，藥柱沿軸向伸長0.58mm，火箭藥柱在高溫狀態(tài)下受熱膨脹，致使藥柱端面與墊片間的間隙變小或消失，導(dǎo)致在初始燃燒階段藥柱端面正常參與燃燒，p-t曲線在初始階段上升緩慢；但隨著燃燒進行到大約0.5s，燃燒室壓力升高到11MPa左右時，高溫燃?xì)忾_始急劇竄入藥柱端面和墊片之間的間隙，點燃尚未參加燃燒的端面部分，總?cè)济婷娣e突然增大，致使燃燒室壓強陡然升高。

藥柱在高溫條件下受熱膨脹伸長量若過大，嚴(yán)重時可使藥柱端面、墊片和發(fā)動機壁面三者處于受壓狀態(tài)，間隙消失，使藥柱端面不能正常燃燒，同時又使藥柱內(nèi)部產(chǎn)生應(yīng)力，過大的應(yīng)力又可誘發(fā)藥柱表面或內(nèi)部產(chǎn)生裂紋，甚至發(fā)生破碎，導(dǎo)致火箭發(fā)動機工作異常。燃燒室壓強階躍式升高很可能會破壞火箭發(fā)動機的結(jié)構(gòu)，應(yīng)避免此種情況的發(fā)生［6］。在火箭發(fā)動機裝配時，藥柱和墊片在軸向方向上要保留適當(dāng)?shù)拈g隙或加入間隙調(diào)節(jié)裝置，以保證火箭發(fā)動機在高溫條件下正常工作。

此外，在常溫、低溫和高溫下，用方法3計算的藥柱燃燒質(zhì)量均為1.444kg，與藥柱實際質(zhì)量1.445kg接近，相對誤差僅為-0.07%。

3 結(jié) 論

（1）在常溫和低溫條件下，星孔裝藥端面是以逐漸參與燃燒的方式進行燃燒，試驗中可得到較為平滑的p-t曲線；而裝藥端面不參與燃燒和裝藥端面瞬時全面參與燃燒的方式與實際情況不符。

（2）通過引入端面燃燒參與度和端面等效燃燒厚度的方法，并結(jié)合試驗條件推導(dǎo)出燃面面積經(jīng)驗計算公式，解決了藥柱在燃燒過程中的燃面計算問題。在常溫和低溫條件下，計算值與試驗值的相對誤差均小于5%，計算得到的p-t曲線與試驗曲線基本吻合；高溫條件下稍顯不足。

（3）采用上、下端面均無包覆星孔裝藥的火箭發(fā)動機工作穩(wěn)定性受溫度影響較大。藥柱端面和墊片之間要留有適當(dāng)?shù)拈g隙或加入間隙調(diào)節(jié)裝置，間隙大小可根據(jù)火箭發(fā)動機的具體藥型尺寸和試驗測試結(jié)果調(diào)整。

（4）端面燃燒參與度和端面等效燃燒厚度的引入提供了一種解決端部初始燃面可變的星孔火箭裝藥燃面計算問題的方法，可為同類固體火箭發(fā)動機內(nèi)彈道性能預(yù)測提供參考。

［1］侯林法，楊忠雄.復(fù)合固體火箭推進劑［M］.北京：中國宇航出版社，1994.

［2］夏志全.端面燃燒發(fā)動機研制中應(yīng)注意的幾個問題［J］.固體火箭技術(shù)，2000，23（2）：12-15.

XIA Zhi-quan.Several problems for end burning solid motor development［J］.Journal of Solid Rocket Technology，2000，23（2）：12-15.

［3］許鵬，梁景媛，梁國柱，等.固體火箭發(fā)動機裝藥端面燃燒燃速畸變現(xiàn)象分析［J］.火炸藥學(xué)報，2006，29（5）：32-34.

XU Peng，LIANG Jing-yuan，LIANG Guo-zhu，et al.Preliminary analysis on cone burning of solid propellant in end-burning rocket motors［J］.Chinese Journal of Explosives and Propellants，2006，29（5）：32-34.

［4］周長省，鞠玉濤，朱福亞.火箭彈設(shè)計理論［M］.北京：北京理工大學(xué)出版社，2005.

［5］董師顏，張兆良.固體火箭發(fā)動機原理［M］.北京：北京理工大學(xué)出版社，1996.

［6］邢耀國，董可海，沈偉.固體火箭發(fā)動機使用工程［M］.北京：國防工業(yè)出版社，2010.