金秉寧，劉佩進(jìn)，徐冠宇

（西北工業(yè)大學(xué) 燃燒、熱結(jié)構(gòu)與內(nèi)流場重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710072）

1 引言

固體推進(jìn)劑燃燒響應(yīng)是燃燒不穩(wěn)定的主要影響因素。近年來，從全尺寸發(fā)動(dòng)機(jī)研制過程中發(fā)現(xiàn)，同樣的裝藥結(jié)構(gòu)，推進(jìn)劑配方的微小變化、原材料批次以及工藝等因素均會導(dǎo)致推進(jìn)劑燃燒響應(yīng)特性的變化，產(chǎn)生燃燒不穩(wěn)定問題。由于缺乏對燃燒響應(yīng)機(jī)理的深入認(rèn)識，發(fā)動(dòng)機(jī)研制中往往采用“改配方-測量-再改配方-再測量”的方法，降低推進(jìn)劑燃燒響應(yīng)的增益作用。但這種研究方法效率低、周期長、實(shí)驗(yàn)次數(shù)多、難以快速準(zhǔn)確地提出合理解決方案，從而給研制者們帶來極大的困難。因此，十分有必要開展實(shí)際配方組分對固體推進(jìn)劑燃燒響應(yīng)影響規(guī)律的理論研究，為發(fā)動(dòng)機(jī)研制提供有力的理論支撐和改進(jìn)依據(jù)。

固體推進(jìn)劑燃燒響應(yīng)理論的研究需從穩(wěn)態(tài)燃燒特性和非穩(wěn)態(tài)燃燒特性［1-2］兩個(gè)方面開展。國外開展了大量研究。其中，Culick［3-5］建立了均質(zhì)推進(jìn)劑一維準(zhǔn)穩(wěn)態(tài)氣相燃燒模型（QSHOD模型），并將高氯酸銨（AP）燃燒模型帶到參數(shù)A和B中進(jìn)行修正，獲得AP粒度對響應(yīng)函數(shù)的影響。Rasmussen等［6］以一維多火焰穩(wěn)態(tài)燃燒模型（BDP）［7］為基礎(chǔ)，建立了一維復(fù)合推進(jìn)劑壓強(qiáng)耦合響應(yīng)模型，獲得了AP含量、粒度以及工作壓強(qiáng)等參數(shù)對燃燒響應(yīng)特性的影響，加深了研究者們對非均質(zhì)推進(jìn)劑燃燒響應(yīng)特性的認(rèn)識。Glick等［8］在燃燒過程中考慮了多元、多級配氧化劑的組成，建立了“小系統(tǒng)綜合（PEM）”穩(wěn)態(tài)燃燒模型。Cai等［9-10］建立了二維三明治模型，獲得了AP粒徑、壓強(qiáng)對燃速和火焰高度的影響以及瞬態(tài)壓強(qiáng)振蕩之后推進(jìn)劑的火焰結(jié)構(gòu)和燃 燒響應(yīng)特 性。Jackson 等［11］采用“Random Pack”模型，對復(fù)合推進(jìn)劑進(jìn)行了三維數(shù)值模擬，獲得了多級配AP復(fù)合推進(jìn)劑的壓強(qiáng)耦合響應(yīng)函數(shù)。但該模型需要的計(jì)算資源較大，計(jì)算時(shí)間較長。Renie和Osborn［12］基于小擾動(dòng)理論，分析了氧化劑顆粒大小、燃燒室壓強(qiáng)和平均流速對燃燒響應(yīng)函數(shù)的影響。國內(nèi)主要開展了穩(wěn)態(tài)燃燒模型方面的研究。其中，周志清等［13］采用多步化學(xué)反應(yīng)動(dòng)力學(xué)機(jī)制，建立了三明治二維細(xì)觀燃燒模型。趙瑜基于詳細(xì)化學(xué)動(dòng)力學(xué)機(jī)制，建立了復(fù)合推進(jìn)劑燃燒模型［14］。張煒等建立了AP/黑索今（RDX）兩種氧化劑含量和粒度級配可任意變化的AP/RDX/端羥基聚丁二烯（HTPB）推進(jìn)劑穩(wěn)態(tài)燃燒模型［15］。彭培根等基于以往復(fù)合推進(jìn)劑燃燒模型和大量實(shí)驗(yàn)現(xiàn)象提出了一個(gè)AP/RDX推進(jìn)劑的復(fù)合多火焰燃燒模型［16］。劉佩進(jìn)等從燃燒學(xué)、聲學(xué)和動(dòng)力學(xué)出發(fā)，對復(fù)合推進(jìn)劑的壓強(qiáng)耦合響應(yīng)開展了理論、實(shí)驗(yàn)和數(shù)值計(jì)算研究［17］，對復(fù)合推進(jìn)劑的燃燒響應(yīng)特性有了較全面的認(rèn)識。

綜上所述，國外對于均質(zhì)和非均質(zhì)推進(jìn)劑燃燒響應(yīng)模型開展了大量的研究，但關(guān)于非勻質(zhì)推進(jìn)劑中氧化劑粒度級配及配比對響應(yīng)函數(shù)影響的研究較少；國內(nèi)對于均質(zhì)和非均質(zhì)推進(jìn)劑的穩(wěn)態(tài)燃燒模型開展了大量的研究，而燃燒響應(yīng)模型則研究較少，關(guān)于推進(jìn)劑配方組分對響應(yīng)函數(shù)影響的認(rèn)識大多是基于實(shí)驗(yàn)測量結(jié)果，對燃燒響應(yīng)機(jī)理缺少深入的認(rèn)識，使得工程研制中較難提出降低推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)的抑制方法。

為此，本研究以BDP多火焰燃燒模型為框架，考慮實(shí)際推進(jìn)劑中AP多級配粒度分布，建立AP多粒度級配的AP/HTPB推進(jìn)劑燃燒響應(yīng)模型，對模型進(jìn)行校驗(yàn)。并針對幾種常用的AP粒度，在不同粒度級配和配比條件下開展非穩(wěn)態(tài)計(jì)算，深入分析AP粒度、級配和配比等參數(shù)變化對壓強(qiáng)耦合響應(yīng)函數(shù)的影響規(guī)律。

2 AP多粒度級配燃燒響應(yīng)模型

BDP多火焰模型假設(shè)復(fù)合推進(jìn)劑僅由單一粒徑AP和 HTPB粘合劑組成［4，6］，屬于雙組元推進(jìn)劑。模型包含AP分解焰、擴(kuò)散火焰和初焰的多火焰結(jié)構(gòu)。AP分解焰為AP分解產(chǎn)物NH3與HClO4的反應(yīng)，產(chǎn)生于AP“氣柱”中，向燃面反饋熱量，屬于預(yù)混火焰。初焰為AP分解產(chǎn)物與粘合劑熱解產(chǎn)物之間的化學(xué)反應(yīng)，受擴(kuò)散和化學(xué)反應(yīng)速率影響，屬于擴(kuò)散火焰。終焰為粘合劑熱解產(chǎn)物和AP焰富氧產(chǎn)物間的反應(yīng)，位于AP分解焰外端，屬于擴(kuò)散火焰。與AP分解焰是相連，并向燃面反饋熱量。該模型的主要假設(shè)如下：

（1）沿垂直燃燒表面方向視為一維；

（2）AP和粘合劑表面分解反應(yīng)符合Arrhenius定律；

（3）AP和粘合劑具有不同的表面溫度；

（4）氣相反應(yīng)為簡單均相反應(yīng)，且產(chǎn)物為完全氣體。

2.1 AP多級配模型構(gòu)建

對于AP多級配推進(jìn)劑來說，可將AP多粒度級配AP/HTPB復(fù)合推進(jìn)劑視為由多個(gè)假定AP單粒度AP/HTPB系統(tǒng)組合而成［18］，假設(shè)AP表面的單元火焰是相互獨(dú)立燃燒，并將推進(jìn)劑重新排列成多組各包含一種AP粒度的假想單級配推進(jìn)劑（如圖1所示）。這里考慮的AP四級配推進(jìn)劑，其中AP粒度分別為Ⅰ類、Ⅱ類、Ⅲ類和細(xì)AP。

圖1 AP四級配假想推進(jìn)劑模型Fig.1 The sketch of Quadruplemodal propellant model

重新排列的AP單級配假想推進(jìn)劑火焰構(gòu)型滿足BDP多火焰模型，其中，AP分解火焰和擴(kuò)散終焰受AP粒度尺寸影響較大，需要根據(jù)不同AP粒度尺寸條件分別計(jì)算，圖2所示為相鄰兩個(gè)不同AP粒徑的火焰模型示意圖。

從圖2可以看出，相鄰的兩個(gè)AP粒子，由于粒度不同，相應(yīng)的分解焰xf，AP，xr和終焰 xf各不相同，并且AP表面溫度Ts，AP和熱反饋量QL也各不相同。因此，多級配模型中的火焰?zhèn)€數(shù)取決于AP粒徑種類數(shù)量，如果是AP四級配，則有四個(gè)不同尺寸的四個(gè) AP分解焰xf，AP，n、四個(gè) AP 表面溫度Ts，AP，n；對于粘合劑，為了簡化計(jì)算，其表面溫度按照小尺寸AP顆粒周圍的粘合劑溫度計(jì)算；多級配推進(jìn)劑的燃燒速率由各個(gè)單級配的結(jié)果按照其AP質(zhì)量含量比例所占權(quán)重綜合統(tǒng)計(jì)計(jì)算，且模型中各組分的質(zhì)量通量由AP單級配BDP火焰模型決定。

2.1.1 穩(wěn)態(tài)燃燒模型

根據(jù)實(shí)際推進(jìn)劑配方AP粒度分布，改進(jìn)模型中Ⅰ類 AP 平均粒徑為DAP，1，含量為α1；Ⅱ類 AP 平均粒徑DAP，2，含量α2；Ⅲ類 AP 平均粒徑DAP，3，含量為α3；細(xì)粒度 AP平均粒徑DAP，4，含量為α4；αAP為 AP總含量，如公式（1）所示，而HTPB含量為1-αAP。

圖2 相鄰兩個(gè)不同AP粒徑的火焰模型示意圖Fig.2 The proposed flame structure for two adjacent AP particles

根據(jù)文獻(xiàn)［6］中AP單級配推進(jìn)劑質(zhì)量通量的表達(dá)式，AP多級配推進(jìn)劑中每種AP單級配推進(jìn)劑的質(zhì)量通量Gp，n和多級配推進(jìn)劑總質(zhì)量通量Gp可表示如下：

式中，αAP為 AP（氧化劑）含量，n則代表AP級配數(shù)，GAP，n為各個(gè)粒度的氧化劑質(zhì)量生成速率，Gb為粘合劑的質(zhì)量生成速率，其表達(dá)式如下［6］：

式中，在多級配模型中，不同AP粒子具有不同表面溫度Ts，AP，n，其表達(dá)式如下所示［6］：

而粘合劑溫度采用統(tǒng)一值Ts，b，取最小尺寸AP顆粒周圍的粘合劑溫度［6］：

相應(yīng)的 AP 多級配的 AP 分解焰火焰溫度Tf，AP，n、AP 分解焰火焰高度xf，AP，n，初始擴(kuò)散火焰高度xr，n，擴(kuò)散終焰高度xf，n均與單級配模型表達(dá)式相一致，與AP等效粒徑有關(guān)，其各個(gè)級配AP等效粒徑的表達(dá)式如下［6］：

2.1.2 非穩(wěn)態(tài)燃燒模型

AP多級配非穩(wěn)態(tài)模型計(jì)算，基本上與AP單級配非穩(wěn)態(tài)模型相同。當(dāng)推進(jìn)劑進(jìn)行非穩(wěn)態(tài)燃燒時(shí)，?T/?t≠0，此時(shí)模型中各個(gè)火焰的高度、溫度等參數(shù)均會隨時(shí)間變化，因此需要考慮AP和粘合劑表面能量方程中隨時(shí)間變化的積分項(xiàng)。其非穩(wěn)態(tài)條件下的AP和粘合劑表面溫度［6］：

式中，多級配 AP 顆粒表面溫度Ts，AP，n和粘合劑表面溫度Ts，b須采用非穩(wěn)態(tài)計(jì)算結(jié)果表征非穩(wěn)態(tài)特性。而AP質(zhì)量通量GAP，n、粘合劑質(zhì)量通量Gb和推進(jìn)劑總質(zhì)量通量Gp是Ts，AP，n、Ts，b和（Ts，AP，n，Ts，b）的 函 數(shù) 。 關(guān) 于xf，n、xr，n、xf，AP，n和Tf，AP，n的計(jì)算，則采用穩(wěn)態(tài)計(jì)算結(jié)果進(jìn)行插值計(jì)算出壓強(qiáng)擾動(dòng)條件下不同時(shí)刻壓強(qiáng)所對應(yīng)的值，作為已知參數(shù)來求解非線性方程組。

2.1.3 AP多級配燃燒響應(yīng)函數(shù)構(gòu)造

穩(wěn)態(tài)和非穩(wěn)態(tài)結(jié)果中包含了多個(gè)參變量，如Gp，n、GAP，n、Gb、Tf，n、Ts，AP，n、Ts，b、xf，n、xr，n、xf，AP，n，這 些 參 數(shù) 作 為非穩(wěn)態(tài)計(jì)算工況的初始參數(shù)，用來計(jì)算推進(jìn)劑內(nèi)部溫度分布和不同時(shí)刻下各參數(shù)變化情況，并根據(jù)計(jì)算結(jié)果獲得推進(jìn)劑的壓強(qiáng)耦合響應(yīng)函數(shù)。壓強(qiáng)耦合響應(yīng)函數(shù)中的為壓強(qiáng)波動(dòng)與平均壓強(qiáng)的比值，是模型初始參數(shù)給定的量；因此，構(gòu)建響應(yīng)函數(shù)的關(guān)鍵參數(shù)即通過模型獲得燃速擾動(dòng)量r′和平均量的比值。然而對于非穩(wěn)態(tài)模型計(jì)算結(jié)果可以獲得燃速擾動(dòng)量r′，通過r′進(jìn)行平均化處理方可得到燃速的平均量。燃速平均量的處理方法如下：

式中，p0為平衡壓強(qiáng)；mag為擾動(dòng)幅值（未特殊說明均采用平衡壓強(qiáng)的20%）；f為擾動(dòng)頻率；t為燃燒時(shí)間。通過改變正弦壓強(qiáng)波動(dòng)函數(shù)的頻率進(jìn)行非穩(wěn)態(tài)計(jì)算可以獲得不同擾動(dòng)頻率下的燃速-時(shí)間曲線。

圖3為平衡壓強(qiáng)在10 MPa，振蕩頻率為175 Hz，振蕩幅值為平衡壓強(qiáng)的20%時(shí)的壓強(qiáng)振蕩時(shí)間曲線以及模型中所有參變量隨時(shí)間振蕩變化曲線，模型中AP總質(zhì)量分?jǐn)?shù)為80%，AP采用兩種粒度（50/200 μm），質(zhì)量配比為8∶1，表面溫度為298 K?？梢钥闯?，隨著壓強(qiáng)振蕩輸入量的變化，模型中AP和粘合劑以及推進(jìn)劑總質(zhì)量通量隨時(shí)間呈周期性的振蕩變化，振蕩頻率與壓強(qiáng)振蕩頻率相同。并且推進(jìn)劑燃燒表面的溫度也隨壓強(qiáng)振蕩的變化呈正弦波動(dòng)，波動(dòng)頻率與壓強(qiáng)振蕩頻率相同。

圖3 振蕩條件下模型參數(shù)振蕩曲線Fig.3 Typical model parameters response to oscillation pressure input

AP多級配推進(jìn)劑燃燒模型采用Matlab計(jì)算平臺編程計(jì)算，與AP單級配模型基本上相同，均是分別求解穩(wěn)態(tài)（公式（1）～（8））和非穩(wěn)態(tài)（公式（9）～（11））的非線性方程組。然而，對于不同AP級配來說，其穩(wěn)態(tài)和非穩(wěn)態(tài)的非線性方程組的數(shù)量不同，級配越多，方程數(shù)量也越多，計(jì)算量也越大，計(jì)算時(shí)間更長。AP單級配燃燒響應(yīng)模型計(jì)算的非線性方程數(shù)量是11個(gè)，而本文建立的AP四級配燃燒響應(yīng)模型則由36個(gè)非線性方程組成。該模型仍是AP和HTPB雙組元推進(jìn)劑，其中AP質(zhì)量分?jǐn)?shù)為73%～87%，初始溫度為298 K，計(jì)算壓強(qiáng)2～20 MPa內(nèi)均可以獲得燃燒響應(yīng)函數(shù)。

2.2 模型校驗(yàn)

2.2.1 穩(wěn)態(tài)模型校驗(yàn)

針對建立的AP多組元AP/HTPB推進(jìn)劑燃燒模型，首先根據(jù)國內(nèi)外文獻(xiàn)資料中的實(shí)驗(yàn)結(jié)果對穩(wěn)態(tài)燃燒模型進(jìn)行驗(yàn)證。分別采用文獻(xiàn)［19］中相關(guān)的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行模型校驗(yàn)。其中AP粒度含量約為80%，AP級配分別為雙級配、三級配和四級配，AP粒度分別采用400，200，90，50 μm和20 μm，基本涵蓋了AP各種粒度分布。表1給出了具體的粒度級配和配比參數(shù)。

模型中用到的 AP 和 HTPB 的指前因子As，AP，As，b、活 化 能Es，AP，Es，b以 及 比熱 容cp，g，AP，cp，s，b等 參數(shù) 值 詳 見文獻(xiàn)［20］。模型計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果如圖4所示。

由圖4可知，在高壓強(qiáng)條件下（＞5 MPa），多級配穩(wěn)態(tài)燃燒模型的計(jì)算值與實(shí)驗(yàn)測量值之間的誤差均小于10%，與實(shí)驗(yàn)測量結(jié)果吻合得較好。同時(shí)，從AP粒度配比的涵蓋范圍來看，包含了實(shí)際推進(jìn)劑中采用的多種粒度分布，具有一定的通用性。因此，可以認(rèn)為所建立的多級配多火焰穩(wěn)態(tài)燃燒模型是可行的，能較好地反映推進(jìn)劑配方參數(shù)變化的燃速特性，可以進(jìn)行非穩(wěn)態(tài)壓強(qiáng)耦合響應(yīng)函數(shù)的計(jì)算工作。

表1 推進(jìn)劑配方參數(shù)及計(jì)算工況Table 1 Propellant formulations and detailed conditions for calculation

圖4 5種AP多級配模型計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果比較Fig.4 A comparison of between the modeling results and the experimental data for five propellants

2.2.2 多級配響應(yīng)函數(shù)驗(yàn)證

針對文獻(xiàn)［21］中AP三級配AP/HTPB推進(jìn)劑配方（表2），開展多級配響應(yīng)函數(shù)計(jì)算，并與文獻(xiàn)中實(shí)驗(yàn)結(jié)果進(jìn)行對比，對多級配響應(yīng)函數(shù)模型進(jìn)行驗(yàn)證。文獻(xiàn)采用脈沖觸發(fā)T型燃燒器實(shí)驗(yàn)測量裝置，實(shí)驗(yàn)測量工作壓強(qiáng)約為6.9 MPa，測量頻率為400～1300 Hz。因此，AP多級配推進(jìn)劑響應(yīng)函數(shù)模型的計(jì)算工作壓強(qiáng)約為6.9 MPa，計(jì)算頻率范圍為200～1300 Hz，AP總含量為87.6%，HTPB含量為12.4%。計(jì)算結(jié)果如圖5所示。

由圖5可知，采用粒度分布為400/200/17 μm的AP三級配推進(jìn)劑，其粒度配比分別為23.23∶34.98∶29.15，其響應(yīng)函數(shù)計(jì)算結(jié)果與實(shí)驗(yàn)擬合結(jié)果在200～1300 Hz內(nèi)，基本上吻合得較好，其中振蕩頻率為250～800 Hz內(nèi)，響應(yīng)函數(shù)值的誤差在5%～9%；而振蕩頻率約為1200Hz時(shí)，模型計(jì)算結(jié)果與實(shí)驗(yàn)擬合結(jié)果（黑色虛線）的誤差值相對較大，約為14%；但該模型計(jì)算結(jié)果仍在實(shí)驗(yàn)測量值（黑色方框□）的上下限范圍內(nèi)，因此認(rèn)為模型計(jì)算的誤差是可接受的。

表2 AP三級配推進(jìn)劑配方參數(shù)［21］Table 2 Formulation of AP Triple modalpropellant［21］

圖5 計(jì)算模型與實(shí)驗(yàn)結(jié)果對比驗(yàn)證Fig.5 A comparison of the predicted and measured Rp values

通過與參考文獻(xiàn)中［21］的實(shí)驗(yàn)測量值進(jìn)行比較發(fā)現(xiàn)，AP多級配推進(jìn)劑的壓強(qiáng)耦合響應(yīng)函數(shù)計(jì)算值與實(shí)驗(yàn)測量值之間的誤差均小于14%，與實(shí)驗(yàn)測量結(jié)果吻合較好?？梢哉J(rèn)為，建立的多級配推進(jìn)劑壓強(qiáng)耦合響應(yīng)函數(shù)計(jì)算模型是可行的，能較好地反映出AP多級配推進(jìn)劑響應(yīng)函數(shù)特性，可以開展AP多級配推進(jìn)劑響應(yīng)函數(shù)特性的計(jì)算研究。

3 計(jì)算結(jié)果及影響規(guī)律

根據(jù)實(shí)際推進(jìn)劑中經(jīng)常采用的四種AP顆粒粒度（330，250，110 μm 和 50 μm），通過調(diào)整各個(gè)粒子含量的配比關(guān)系，分別得到AP多粒度級配的AP/HTPB推進(jìn)劑，通過模型開展穩(wěn)態(tài)和非穩(wěn)態(tài)計(jì)算，獲得不同的級配和配比條件下各種粒度顆粒對響應(yīng)函數(shù)特性的影響。模型中推進(jìn)劑仍為雙組元AP/HTPB復(fù)合推進(jìn)劑，其中AP含量為80%，平衡工作壓強(qiáng)為10 MPa，壓強(qiáng)振蕩幅值為平衡壓強(qiáng)的20%，振蕩頻率為25～1000 Hz。

3.1 AP粒度配比對燃燒響應(yīng)的影響規(guī)律

為獲得 330，250，110，50 μm四種 AP粒度中某一粒度含量變化對響應(yīng)函數(shù)影響，粒度配比變化規(guī)則為：AP總含量不變，一次只改變一種AP粒度的含量（a），相應(yīng)其它三種粒度含量同時(shí)改變，計(jì)算結(jié)果如圖6所示。

由圖6可知，不同AP粒度主導(dǎo)變化對響應(yīng)函數(shù)的變化規(guī)律有明顯的不同。其中圖6a中隨著330 μm粒度AP含量的降低，中低頻（＜250 Hz）時(shí)，含量從8∶1∶1∶1降低至1∶1∶1∶1時(shí)，響應(yīng)函數(shù)值快速降低，其中峰值變化最大，從1.7降低至1.35，但峰值頻率變化不大；而含量從1∶1∶1∶1降低至0.125∶1∶1∶1時(shí)，響應(yīng)函數(shù)值基本變化不大。對于中高頻來說（＞250 Hz），隨著含量的降低，響應(yīng)函數(shù)值增大；并且頻率越大，響應(yīng)函數(shù)值變化較大。

圖6b中，250 μm粒度AP含量變化對響應(yīng)函數(shù)的影響基本上與330 μm相同，響應(yīng)函數(shù)變化仍在不同頻率范圍變化不同，即中低頻范圍內(nèi)，隨著a250含量降低，響應(yīng)函數(shù)值下降，峰值變化從1.56降至1.34；在中高頻范圍內(nèi)，隨著a250含量降低，響應(yīng)函數(shù)值上升。單臨界頻率約為400 μm，向高頻移動(dòng)。

圖6c中，在整個(gè)頻率范圍內(nèi)，隨著a110含量降低，響應(yīng)函數(shù)值均下降；但不同頻率范圍內(nèi)下降的趨勢不同，低頻范圍內(nèi)變化趨勢較為緩慢，而高頻范圍內(nèi)變化趨勢較快。并且隨著響應(yīng)函數(shù)值降低，響應(yīng)函數(shù)峰值頻率向低頻移動(dòng)。

圖6d中，以含量1∶1∶1∶1為臨界情況，當(dāng)a110含量＞該臨界條件時(shí)，在整個(gè)頻率范圍內(nèi)，隨著a110含量降低，響應(yīng)函數(shù)呈整體型降低的趨勢變化；而當(dāng)a110含量＜臨界條件時(shí)，隨著a110含量降低，響應(yīng)函數(shù)呈區(qū)域型變化，其中中低頻范圍內(nèi)響應(yīng)函數(shù)上升，而中高頻范圍內(nèi)，響應(yīng)函數(shù)下降。a110含量的降低，對于峰值頻率向低頻移動(dòng)。

圖6 AP四級配單粒度配比變化對響應(yīng)函數(shù)影響（10 MPa）Fig.6 Effect of the change of particle content on response function in Quadruplemodal propellants（10 MPa）

圖7 擴(kuò)散火焰高度與AP、粘合劑表面溫度關(guān)系（10 MPa）Fig.7 Change of diffusion flame height and surface temperature of AP and the binder with particle size of AP（10 MPa）

圖8 AP和粘合劑的質(zhì)量通量以及特征時(shí)間隨AP粒度變化Fig.8 Change of mass flux and characteristic response time of AP and binder with particle size of AP

圖9 四種AP粒度響應(yīng)函數(shù)分布曲線Fig.9 Effect of the particle size of AP on Rp

3.2 AP粒度級配對燃燒響應(yīng)的影響規(guī)律

針對雙級配 110 μm/50 μm、250 μm/50 μm 和330 μm/50 μm三種不同AP粒度的雙級配推進(jìn)劑，采用粗細(xì)粒度配比為4∶1，分別用250 μm或者330 μm替代雙級配中粗粒度AP的50%含量，而小粒度50 μm含量保持不變，形成三種AP三粒度級配推進(jìn)劑，配比為2∶2∶1。AP粒度級配對響應(yīng)函數(shù)結(jié)果如圖10所示。

圖10 AP粒度級配對響應(yīng)函數(shù)的影響Fig.10 Effect of the AP modal on Rp

從圖10可以看出，在110 μm/50 μmAP雙級配推進(jìn)劑中，當(dāng)用 250 μm 或 330 μm AP替代 50% 的 110 μm AP含量時(shí)，響應(yīng)函數(shù)在不同頻率范圍變化不同：小于200 Hz時(shí)，隨著更大粒度的加入，響應(yīng)函數(shù)增大；大于200 Hz時(shí)，則趨勢相反。在 330 μm/50 μm AP雙級配推進(jìn)劑中，當(dāng)用250 μm或者110 μm AP替代50%的330 μm，響應(yīng)函數(shù)仍是在不同頻率范圍變化不同：小于200 Hz時(shí)，隨著更大粒度的加入，響應(yīng)函數(shù)快速降低：大于200 Hz時(shí)，趨勢相反。對于AP粒度組合250 μm/50 μm 推進(jìn)劑來說，當(dāng)用 330 μm或者 110 μm AP替代50%的250 μm AP含量時(shí)，響應(yīng)函數(shù)的變化趨勢完全相反；用大于原有粗粒度AP替代時(shí)，響應(yīng)函數(shù)變化規(guī)律與110 μm/50 μm AP雙級配推進(jìn)劑變化規(guī)律相同；而用小于原有粗粒度AP替代時(shí)，響應(yīng)函數(shù)變化規(guī)律與330 μm/50 μm AP雙級配推進(jìn)劑結(jié)果變化規(guī)律相同。

因此，對于AP雙級配推進(jìn)劑來說，用某一不同粒度AP替代雙級配中粗粒度含量時(shí)，首選需要判斷該替換AP粒度與被替換AP粒度大小關(guān)系，如果替換AP粒度＞被替換AP粒度時(shí)，在中低頻響應(yīng)函數(shù)呈增大的趨勢變化，而在中高頻呈降低的趨勢變化；如果AP粒度＞被替換AP粒度時(shí)，響應(yīng)函數(shù)變化規(guī)律相反，即在中低頻呈下降的趨勢變化，而在中高頻呈增大的趨勢變化。此規(guī)律可以同樣適用于AP三級配和四級配推進(jìn)劑。因此，對于AP粒度級配對響應(yīng)函數(shù)的影響，需要根據(jù)AP粒度替代關(guān)系進(jìn)行具體分析，并且在不同的振蕩頻率條件下，變化規(guī)律完全不同。

4 結(jié)論

建立了AP多級配的AP/HTPB推進(jìn)劑燃燒響應(yīng)模型，針對四種 AP粒度 330，250，110 μm 和 50 μm 組合成的多種AP多級配推進(jìn)劑，分別研究了響應(yīng)函數(shù)特性，獲得了AP粒度配比和粒度級配對響應(yīng)函數(shù)分布的影響規(guī)律。主要結(jié)論如下：

（1）建立的AP多粒度多級配燃燒響應(yīng)模型能較好地反映出AP多粒度多級配推進(jìn)劑燃燒響應(yīng)特性，100～1000 Hz內(nèi)，燃燒響應(yīng)模型計(jì)算結(jié)果與文獻(xiàn)實(shí)驗(yàn)測量結(jié)果吻合的較好，誤差小于9%。

（2）AP粒度對響應(yīng)函數(shù)分布的影響較大，影響規(guī)律基本滿足：小粒度增益中高頻、抑制中低頻，大粒度增益中低頻、抑制中高頻。即大粒度含量增大，中低頻響應(yīng)函數(shù)值增大，中高頻響應(yīng)函數(shù)值降低；而細(xì)粒度含量增大時(shí)，中低頻響應(yīng)函數(shù)值降低，中高頻響應(yīng)函數(shù)值增大。

（3）AP級配對響應(yīng)函數(shù)分布的響應(yīng)影響也較大，用中粒度AP替代部分小粒度AP，可以抑制中高頻振蕩，但同時(shí)會增益中低頻振蕩；用中粒度AP替代大粒度AP時(shí)，則結(jié)果相反。