徐大軍，朱云松，湯永康，宮曉睿

（北京航空航天大學(xué)，北京 102206）

關(guān)鍵字：導(dǎo)流格柵；沖擊波；導(dǎo)彈發(fā)射箱；動(dòng)網(wǎng)格

0 引言

在沖擊波開蓋過程中，發(fā)動(dòng)機(jī)破膜瞬間，燃燒室內(nèi)高壓燃?xì)馀c周圍空氣形成初始?jí)毫缑?。隨后，高壓界面在燃?xì)馍淞鞯淖饔孟孪蛲馔七M(jìn)，高壓界面壓縮并加厚周圍空氣層，形成初始沖擊波，同時(shí)伴隨能量的不斷增強(qiáng)[1]。沖擊波作用于后蓋上，在后蓋的反射作用下沿著彈體與發(fā)射箱之間的間隙向前蓋傳播，最終完成導(dǎo)彈發(fā)射箱前蓋的開蓋過程[2-3]。相對(duì)于燃?xì)馐介_蓋與頂破式開蓋，沖擊波開蓋方式能夠有效避免在導(dǎo)彈發(fā)射過程中存在的發(fā)射箱內(nèi)壓力過大、溫度過高的問題，且能夠避免因前蓋碎片與導(dǎo)彈相撞影響導(dǎo)彈飛行軌跡的問題[4]。在沖擊波開蓋方式中，后蓋開蓋壓力的選擇至關(guān)重要。后蓋開蓋過早，發(fā)射箱內(nèi)無法積蓄足夠強(qiáng)度的沖擊波，發(fā)射箱開蓋時(shí)間過晚，形成的沖擊波強(qiáng)度過大，可能會(huì)損傷彈體。在發(fā)射箱開蓋過程中，沖擊波首先作用在后蓋中心位置，可能使得后蓋無法正常碎裂[5]，針對(duì)這一問題，潘登提出在噴管與后蓋之間的區(qū)域放置導(dǎo)流格柵，可有效改善后蓋壓強(qiáng)不均勻的狀況[6]，但是導(dǎo)流格柵對(duì)發(fā)射箱內(nèi)沖擊波強(qiáng)度的影響尚缺乏有效的仿真計(jì)算分析。

本文對(duì)帶有導(dǎo)流格柵的計(jì)算模型進(jìn)行了數(shù)值仿真計(jì)算研究。結(jié)果表明，在相同后蓋開蓋壓力下，導(dǎo)流格柵能夠增強(qiáng)發(fā)射箱內(nèi)的沖擊波強(qiáng)度。

1 仿真模型

1.1 計(jì)算模型

以某型導(dǎo)彈及發(fā)射箱為仿真計(jì)算對(duì)象，導(dǎo)彈在發(fā)射過程中采取傾斜發(fā)射方式，計(jì)算模型主要包括箱體、彈體、發(fā)動(dòng)機(jī)噴管、前后蓋、外流場等?？紤]到計(jì)算資源的有效性以及研究的主次關(guān)系，對(duì)計(jì)算模型進(jìn)行了以下簡化：

1）假定箱內(nèi)氣體為理想氣體，且忽略燃?xì)馍淞髦械幕瘜W(xué)反應(yīng)影響；

2）假定發(fā)射箱箱體與彈體都是剛體，不考慮塑性變形對(duì)流場的影響；

3）忽略彈翼及箱內(nèi)導(dǎo)軌部件的影響。

同時(shí)在發(fā)射箱內(nèi)設(shè)置監(jiān)測點(diǎn)1、2和3用于監(jiān)測發(fā)射箱內(nèi)壓強(qiáng)變化情況。計(jì)算模型如圖1所示。

圖1 計(jì)算模型示意圖Fig.1 Schematic diagram of calculation model

首先，對(duì)計(jì)算模型進(jìn)行網(wǎng)格劃分[7]。箱體內(nèi)部采用結(jié)構(gòu)化網(wǎng)格，外流場區(qū)域采用非結(jié)構(gòu)化網(wǎng)格，網(wǎng)格總數(shù)約為50.2萬。

1.2 控制方程

燃?xì)饬鲌霾捎们蠼夥嵌ǔ＠字Z平均的N-S方程進(jìn)行數(shù)值模擬[8-9]，其控制方程[10]如下。

質(zhì)量守恒方程：

式（1）～（3）中：Sm為用戶定義的源項(xiàng)；t為時(shí)間；ρ為靜壓；gi是重力在i方向上的投影；Fi為外部力在i方向上的投影，F(xiàn)i代表了相關(guān)源項(xiàng)；τij為應(yīng)力張量；T為溫度；k為流體的傳熱系數(shù)；ST為黏性耗散項(xiàng)；cp為比熱容。

湍流模型采用標(biāo)準(zhǔn)k-ε模型[11]，它要求了解湍動(dòng)能以及耗散率方程。通過精確的方程推導(dǎo)能夠得到湍動(dòng)能輸運(yùn)方程，耗散率方程是通過數(shù)學(xué)上模擬相似原形方程以及物理推理得到的。該模型假設(shè)流動(dòng)為完全湍流，忽略分子間黏性的影響。因此，標(biāo)準(zhǔn)k-ε模型只適合完全湍流的流動(dòng)過程模擬。

標(biāo)準(zhǔn)的k-ε模型的湍動(dòng)能k和耗散率ε方程如下：

1.3 邊界條件設(shè)置

仿真計(jì)算過程采用瞬態(tài)求解，對(duì)流場計(jì)算邊界條件做如下設(shè)定。

1）壓力入口。將發(fā)動(dòng)機(jī)噴管入口作為壓力入口，利用UDF函數(shù)定義壓力時(shí)間歷程曲線，并將其加載到瞬態(tài)計(jì)算過程中。圖2給出了發(fā)動(dòng)機(jī)的壓力時(shí)間變化曲線。

圖2 發(fā)動(dòng)機(jī)壓力時(shí)間變化曲線Fig.2 Pressure change curves of the engine

2）壓力出口邊界。將外流場計(jì)算區(qū)域邊界指定為壓力出口，壓力出口邊界壓力值取大氣壓強(qiáng)P=101 325 Pa，溫度取300 K。

3）壁面邊界條件。彈體表面、發(fā)射箱壁面及發(fā)動(dòng)機(jī)外表面采用壁面邊界條件；物面邊界采用無滑移壁面和絕熱壁面邊界條件；近壁面計(jì)算采用標(biāo)準(zhǔn)壁面函數(shù)法處理。

4）后蓋開蓋方式。當(dāng)前關(guān)于燃?xì)馍淞鲉栴}的研究主要通過數(shù)值仿真以及實(shí)驗(yàn)2 種方法進(jìn)行[12-14]。其中，后蓋開蓋方式的選擇對(duì)數(shù)值仿真結(jié)果存在很大影響。研究表明，相比于傳統(tǒng)的燃?xì)馍淞髁鲌雠c后蓋裂片未耦合的仿真方式，采用基于光順和重構(gòu)方法的動(dòng)網(wǎng)格技術(shù)[15-16]，將燃?xì)馍淞髁鲌雠c后蓋裂片運(yùn)動(dòng)相互耦合得到的仿真結(jié)果更加接近實(shí)驗(yàn)數(shù)據(jù)[17]。

圖3 給出了后蓋簡化模型，后蓋在達(dá)到開蓋壓力后逐步打開，在后蓋打開方式上，采用燃?xì)饬髋c后蓋運(yùn)動(dòng)相耦合的開蓋方式，即后蓋達(dá)到開蓋壓力時(shí)，后蓋裂縫（如圖4所示）先打開，延遲1 ms后，4個(gè)后蓋分塊繞著箱體交線旋轉(zhuǎn)至90°，之后將后蓋邊界條件由“壁面”（wall）變?yōu)椤皟?nèi)部流場”（interior），同時(shí)在箱內(nèi)設(shè)置3個(gè)監(jiān)測點(diǎn)用于監(jiān)測沖擊波在發(fā)射箱內(nèi)的傳播規(guī)律。

圖3 發(fā)射箱后蓋簡化模型Fig.3 Simplified model of rear cover of the launch container

圖4 后蓋裂縫Fig.4 Rear cover crack

2 仿真結(jié)果及分析

2.1 無導(dǎo)流格柵的計(jì)算模型仿真結(jié)果

發(fā)動(dòng)機(jī)堵蓋對(duì)發(fā)射箱內(nèi)沖擊波形成及強(qiáng)度具有很大的影響，堵蓋處通常安放點(diǎn)火器，并對(duì)點(diǎn)火后起到增壓的效果[18]。

發(fā)動(dòng)機(jī)點(diǎn)火后，堵蓋達(dá)到破膜壓強(qiáng)1.5 MPa 時(shí)堵蓋破裂。在仿真過程中表示為：堵蓋由原來的“壁面”轉(zhuǎn)換為“內(nèi)部流場”，堵蓋破裂后，沖擊波先于燃?xì)饬髯饔玫桨l(fā)射箱后蓋上，后蓋絕對(duì)壓強(qiáng)達(dá)到2 MPa時(shí)，后蓋打開。圖5～7給出不同時(shí)刻下的壓強(qiáng)及溫度云圖。

圖5 后蓋達(dá)到開蓋壓力時(shí)流場仿真結(jié)果Fig.5 Results of flowfield simulation when the rear cover reaches the opening pressure

圖5中，發(fā)動(dòng)機(jī)點(diǎn)火后堵蓋破裂，沖擊波先于燃?xì)饬髯饔糜诤笊w上，使得后蓋壓力上升。在后蓋達(dá)到開蓋壓力后，后蓋裂縫破裂，沖擊波一方面通過裂縫向箱外傳播，一方面繼續(xù)向前蓋傳播。圖6 給出后蓋打開3 ms后，后蓋在燃?xì)饬髋c沖擊波作用下圍繞各自與箱體軸線做定軸轉(zhuǎn)動(dòng)，在后蓋開蓋過程中，部分沖擊波繼續(xù)沿著箱體與彈體之間的間隙向前蓋傳播。圖7給出后蓋打開6 ms后，后蓋在燃?xì)饬髯饔孟聡@后蓋與箱體軸線逐步打開至90°，此時(shí)后蓋對(duì)沖擊波與燃?xì)饬鞑淮嬖诜瓷渥饔谩?/p>

圖6 后蓋打開3 ms后流場仿真結(jié)果Fig.6 Results of flowfield simulation when the rear cover opens after 3 ms

圖7 后蓋打開6 ms后流場仿真結(jié)果Fig.7 Results of flowfield simulation when the rear cover opens after 6 ms

圖8給出了3個(gè)監(jiān)測點(diǎn)壓強(qiáng)時(shí)間曲線，可以看出：沖擊波在后蓋反射下向前蓋傳播，在傳播過程中，監(jiān)測點(diǎn)1、2 和3 相繼達(dá)到初次峰值；隨著沖擊波在傳播過程中，由于氣體黏性和發(fā)射箱壁面黏性的作用，造成沖擊波能量不斷損耗，監(jiān)測點(diǎn)1、2 和3 的沖擊波初次峰值依次降低。

圖8 3個(gè)監(jiān)測點(diǎn)壓強(qiáng)時(shí)間曲線Fig.8 Pressure-time curve of three monitoring points

從圖9能夠看到，沖擊波作用到前蓋后，在前蓋的阻擋積聚作用下，前蓋壓力不斷上升。本文初定的前蓋開啟壓強(qiáng)為0.17 MPa，對(duì)應(yīng)圖9 的0.016 5 ms 開啟。由于本次仿真過程中，前蓋始終處于未打開狀態(tài)，所以在前蓋壓強(qiáng)達(dá)到開啟壓力后依舊上升。在前蓋作用下，沖擊波向發(fā)射箱后蓋方向發(fā)生反射，此次沖擊波傳播路徑與初次傳播路徑相反。從圖8能夠看出，監(jiān)測點(diǎn)3、2 和1 相繼達(dá)到第2 次壓強(qiáng)峰值，且由于沖擊波在傳播過程中的損耗，監(jiān)測點(diǎn)3、2 和1 的壓強(qiáng)峰值依次降低。

圖9 前蓋壓強(qiáng)時(shí)間曲線Fig.9 Pressure-time curve of front cover

2.2 帶有導(dǎo)流格柵（A位置）的發(fā)射箱模型仿真結(jié)果

在上述發(fā)射箱模型下加入導(dǎo)流格柵，新添加的導(dǎo)流格柵采取“壁面”邊界條件，其他簡化條件與上述仿真過程一致。導(dǎo)流格柵為1 個(gè)多孔矩形平板，指定平板的厚度為5 mm。導(dǎo)流格柵安裝在發(fā)動(dòng)機(jī)噴管與發(fā)射箱后蓋之間的區(qū)域，導(dǎo)流格柵的形狀以及在發(fā)射箱的安裝方式如圖10、11所示。導(dǎo)流格柵分別在A和B位置安裝：A位置距離噴管口較遠(yuǎn)，導(dǎo)流格柵前端面距離噴口的軸向距離為150 mm；B 位置距離噴管口較近，導(dǎo)流格柵前端面距離噴口的軸向距離為100 mm。

圖10 導(dǎo)流格柵Fig.10 Fluid guiding grid

圖11 導(dǎo)流格柵工況示意圖Fig.11 Schematic diagram of fluid guiding grid

通過數(shù)值仿真，得到了有導(dǎo)流格柵工況下導(dǎo)彈發(fā)射箱內(nèi)各監(jiān)測點(diǎn)及前蓋的壓強(qiáng)變化，并與無導(dǎo)流格柵工況下的仿真結(jié)果形成對(duì)比，各監(jiān)測點(diǎn)及前蓋壓強(qiáng)對(duì)比曲線如圖12～15所示。在發(fā)射箱模型中添加導(dǎo)流格柵后，箱內(nèi)各監(jiān)測點(diǎn)以及前蓋的壓強(qiáng)都有所上升，表1給出了詳細(xì)的數(shù)據(jù)對(duì)比。

表1 有無導(dǎo)流格柵下數(shù)據(jù)對(duì)比Tab.1 Comparison of data with and without fluid guiding grid

圖12 監(jiān)測點(diǎn)1壓強(qiáng)對(duì)比曲線Fig.12 Pressure comparison curve of monitor 1

圖13 監(jiān)測點(diǎn)2壓強(qiáng)對(duì)比曲線Fig.13 Pressure comparison curve of monitor 2

圖14 監(jiān)測點(diǎn)3壓強(qiáng)對(duì)比曲線Fig.14 Pressure comparison curve of monitor 3

圖15 前蓋壓強(qiáng)對(duì)比曲線Fig.15 Pressure comparison curve of front cover

從表1可以看出，存在導(dǎo)流格柵工況下，發(fā)射箱內(nèi)各個(gè)監(jiān)測點(diǎn)沖擊波強(qiáng)度都有所增強(qiáng)，而且相同的后蓋開蓋壓力下，前蓋的壓強(qiáng)峰值增幅達(dá)到5%。

如圖16 所示，堵蓋破裂后，初始沖擊波首先沖擊在導(dǎo)流格柵上，一部分沖擊波從導(dǎo)流格柵中的空格穿過，其他部分在導(dǎo)流格柵的阻攔下向前蓋反射，更短的傳播距離使得沖擊波在傳播過程中的損耗也有所降低。

圖16 后蓋達(dá)到開蓋壓力流場仿真結(jié)果Fig.16 Results of flowfield simulation when the rear cover reaches the opening pressure

2.3 帶有導(dǎo)流格柵（B位置）發(fā)射箱模型仿真結(jié)果

調(diào)整導(dǎo)流格柵在發(fā)射箱中的位置，使之更加靠近發(fā)動(dòng)機(jī)噴管，能夠得到發(fā)射箱內(nèi)各個(gè)監(jiān)測點(diǎn)壓強(qiáng)曲線并與A位置下的數(shù)據(jù)做對(duì)比得到下列曲線。

圖17～20分別給出2種工況下各個(gè)監(jiān)測點(diǎn)及前蓋壓強(qiáng)對(duì)比曲線?？梢钥吹剑贐 位置工況下，各個(gè)監(jiān)測點(diǎn)及前蓋的壓強(qiáng)都有所增加，撞擊在B位置的沖擊波反射路徑相對(duì)于A位置更短，沖擊波在傳播過程中損耗也更少。

圖17 監(jiān)測點(diǎn)1壓強(qiáng)對(duì)比曲線Fig.17 Pressure comparison curve of monitor 1

圖18 監(jiān)測點(diǎn)2壓強(qiáng)對(duì)比曲線Fig.18 Pressure comparison curve of monitor 2

圖19 監(jiān)測點(diǎn)3壓強(qiáng)對(duì)比曲線Fig.19 Pressure comparison curve of monitor 3

圖20 監(jiān)測點(diǎn)4壓強(qiáng)對(duì)比曲線Fig.20 Pressure comparison curve of front cover

3 結(jié)論

1）在發(fā)射箱內(nèi)的發(fā)動(dòng)機(jī)噴管與后蓋之間添加導(dǎo)流格柵后，能夠增強(qiáng)導(dǎo)彈發(fā)射箱內(nèi)沖擊波的強(qiáng)度，其中，前蓋的壓強(qiáng)峰值增幅在5%左右。

2）導(dǎo)流格柵在2個(gè)安裝位置下對(duì)比結(jié)果顯示，導(dǎo)流格柵的安裝位置越靠近發(fā)動(dòng)機(jī)噴管，其反射到前蓋的沖擊波壓強(qiáng)越大。本次仿真中，B 位置工況仿真結(jié)果中的前蓋壓強(qiáng)峰值相比于A 位置工況，其增幅在5%左右。

3）導(dǎo)流格柵的安裝及其位置對(duì)發(fā)射箱內(nèi)沖擊波強(qiáng)度存在較大影響，因此在今后的仿真計(jì)算中不能將其忽略，應(yīng)加以考慮。