羅 丹,鄭 宇,劉曉軍,王 勇,劉 俊

(1 南京理工大學(xué)機械工程學(xué)院,南京 210094;2 遼沈工業(yè)集團有限公司,沈陽 110045;3 北京礦冶研究總院,北京 100015)

?

新型炸藥ANPyO熱刺激響應(yīng)仿真研究

羅 丹1,鄭 宇1,劉曉軍2,王 勇3,劉 俊1

(1 南京理工大學(xué)機械工程學(xué)院,南京 210094;2 遼沈工業(yè)集團有限公司,沈陽 110045;3 北京礦冶研究總院,北京 100015)

為分析新型鈍感炸藥2,6-二氨基-3,5-二硝基吡啶-1-氧化物(ANPyO)的熱響應(yīng)規(guī)律,利用Fluent軟件建立了帶殼炸藥烤燃模型,對ANPyO在不同升溫速率下(1 ℃/min、2 ℃/min、10 ℃/min、20 ℃/min、80 ℃/min、140 ℃/min)的烤燃情況進行了數(shù)值仿真模擬,并與鈍化黑索今(RDX)進行了比較。結(jié)果表明:隨著升溫速率增大,炸藥點火時間縮短;在同一個升溫速率下,ANPyO的點火時間是鈍化RDX的1.6倍左右,點火溫度比鈍化RDX高出145 ℃左右,可見ANPyO熱刺激鈍感于鈍化RDX。

鈍感炸藥;ANPyO炸藥;烤燃;點火時間;數(shù)值仿真

0 引言

在意外火災(zāi)等異常熱環(huán)境下,炸藥可能發(fā)生燃燒、爆炸等重大事故,從而造成不可挽回的重大損失。為模擬炸藥在熱環(huán)境下的安全性能,國內(nèi)外開展了大量研究工作,如智小琦等[1]設(shè)計了密封狀態(tài)不同的多種烤燃彈結(jié)構(gòu),用快速烤燃實驗方法,研究了殼體密封條件對RDX基傳爆藥的烤燃響應(yīng)特性。在進行試驗研究的同時,使用數(shù)值模擬方法評估含能材料的熱安全性已成為新的研究途徑,如Victor等[2]對烤燃試驗進行了數(shù)值模擬,認為含能材料熱分解和熱傳遞符合Frank-Kamenetskii方程。

近年,新型鈍感炸藥ANPyO因其具有起爆感度低、耐高溫的優(yōu)點受到國內(nèi)外廣泛關(guān)注[3-5],成健等[6-7]研究了ANPyO炸藥的合成新方法與性能。何志偉等[8]提出了ANPyO炸藥的精制工藝方法,研究了包覆對ANPyO炸藥的降感作用。劉華寧等[9]分析了ANPyO的射流沖擊感度,程波等[10]研究了約束條件對ANPyO快烤響應(yīng)的影響,對5種不同約束強度下的ANPyO烤燃彈進行了外部火燒試驗研究。

文中通過烤燃數(shù)值仿真評估其在不同升溫速率下的響應(yīng)規(guī)律,對于分析帶殼ANPyO炸藥在熱刺激不敏感度上的優(yōu)勢與應(yīng)用前景均具有重要意義。

1 理論分析

環(huán)境中熱量傳入反應(yīng)系統(tǒng),炸藥發(fā)生自熱反應(yīng)分解放熱,當(dāng)反應(yīng)系統(tǒng)產(chǎn)生的熱量不能及時散失,系統(tǒng)溫度上升,使炸藥分解加速,系統(tǒng)形成自熱,溫度更快的上升,最終導(dǎo)致燃燒或爆炸。用于研究炸藥熱反應(yīng)的經(jīng)典模型主要是Semenov模型和Frank-Kamenetskii模型,前者主要用于描述均溫系統(tǒng),后者則描述非均溫系統(tǒng),他們構(gòu)成了炸藥熱反應(yīng)理論的基礎(chǔ)[11]。

1.2 數(shù)值仿真理論方程

為降低戰(zhàn)斗部熱易損性提供可靠的技術(shù)支持,建立烤燃過程的物理模型來研究戰(zhàn)斗部在熱刺激作用下的響應(yīng)機制成為必不可少的環(huán)節(jié)。考慮到項目研究內(nèi)容多,試驗成本較高,因此通過數(shù)值仿真研究炸藥的烤燃響應(yīng)規(guī)律也是一種重要途徑。利用計算流體力學(xué)軟件Fluent對ANPyO和鈍化RDX炸藥進行了烤燃數(shù)值模擬研究,計算時,對炸藥的烤燃過程作如下假設(shè)[12]:

1)炸藥為均質(zhì)固相,化學(xué)反應(yīng)是零級放熱反應(yīng),炸藥不發(fā)生相變,反應(yīng)物沒有消耗;

2)炸藥和彈體間無空隙,反應(yīng)區(qū)內(nèi)僅有熱傳導(dǎo),反應(yīng)物質(zhì)無運動,因此無對流傳熱;

3)炸藥的自熱反應(yīng)遵循Arrhenius方程;

4)炸藥和鋼殼的物理化學(xué)參數(shù)為常數(shù)的,在反應(yīng)過程中保持不變。

采用SPSS21.0軟件處理數(shù)據(jù)，計數(shù)資料用率表示，比較采用χ2檢驗，計量資料采用（± s）形式表示，對術(shù)后患者尿潴留與獨立變量之間的相關(guān)性采用多因素回歸分析模型進行研究，P＜0.05，差異具有統(tǒng)計學(xué)意義。

對于非均溫系統(tǒng),計算中藥柱的熱反應(yīng)及熱傳導(dǎo)遵循Frank-Kamenetskii方程[13]:

(1)

上式的物理含義為:左邊是單位體積炸藥在單位時間內(nèi)升高溫度所要求的熱量,右邊第一項為由熱傳導(dǎo)流入或流出的熱量,第二項為單位體積內(nèi)、單位時間內(nèi)炸藥熱分解釋放的能量,即炸藥自熱反應(yīng)源項S,用Arrhenius方程表述為:

(2)

上兩式中:ρ為反應(yīng)物的密度(kg/m3);cv為比熱容(J/(kg·K));λ為熱導(dǎo)率(J/(m·K·s));a為反應(yīng)分數(shù);Q為反應(yīng)熱(J/kg);A為指前因子(1/s);E為活化能(J/mol);R為氣體常數(shù)(J/(mol·K))。

2 數(shù)值仿真及結(jié)果分析

2.1 模型建立

烤燃物理模型如圖1所示,45#鋼殼體圓柱部、上下底端厚度均為3 mm,藥柱直徑30 mm,長120 mm,因為藥柱為圓柱型,為減少計算量,采用二分之一模型計算。由于仿真參數(shù)時間步長設(shè)定為1 s,模型在1 s內(nèi)溫度有一定變化,所得點火時對稱面溫度等高線圖與實際點火時模型溫度有一定誤差,但規(guī)律一致。

仿真計算中,采用Tet/Hybrid網(wǎng)格劃分,炸藥劃分Spacing設(shè)置為1.5,殼體劃分Spacing設(shè)置為2,設(shè)殼體外壁為加熱邊界,殼體與炸藥接觸熱條件為Coupled,殼體升溫速率設(shè)置和炸藥自熱反應(yīng)源項S通過C語言編寫為UDF子程序加載到Fluent軟件中進行計算。計算的初始溫度為293 K(20 ℃)。采用壓力基隱式求解器及非定常流動計算,選擇標(biāo)準(zhǔn)k-epsilon湍流模型,并啟動能量方程。

圖1 烤燃物理模型圖

2.2 材料參數(shù)

ANPyO高聚物粘結(jié)炸藥,其粘結(jié)劑氟橡膠F2311的含量為5%,數(shù)值仿真計算參數(shù)如表1所示。

表1 材料計算參數(shù)

2.3 不同升溫速率下ANPyO的響應(yīng)規(guī)律

模擬計算得到ANPyO在6種不同升溫速率下炸藥的熱響應(yīng)規(guī)律,點火時模型的對稱面溫度等高線見圖2所示。

由圖2可見,隨著升溫速率增大,炸藥的點火時間縮短,炸藥點火位置由藥柱中心移向藥柱邊緣。

ANPyO在不同升溫速率下炸藥的點火溫度Ti和點火時藥柱的中心溫度Tc,以及由點火時間和升溫速率計算出來得到的殼體溫度Ts,也就是點火時的環(huán)境溫度,這3種溫度的變化趨勢見圖3所示。

由圖3可見,升溫速率對Ti影響不大;隨著升溫速率的增加,Ts升高,也就是點火環(huán)境溫度升高,說明在慢烤條件下炸藥點火時的環(huán)境溫度低,危險性更大;隨著升溫速率的增加,Tc降低,且降低的趨勢變平緩,這是因為在升溫速率較大的情況下,藥柱溫度來不及傳到中心已經(jīng)達到點火溫度,此時點火位置在藥柱邊緣,而藥柱中心溫度升高的熱量基本上都來自于熱傳導(dǎo)項,升溫不明顯。

圖2 ANPyO點火時模型截面溫度分布圖

圖3 溫度Ti、Ts和Tc隨升溫速率變化曲線

2.4 與鈍化RDX熱響應(yīng)規(guī)律的比較

通過數(shù)值模擬計算得到鈍化RDX在6種不同升溫速率下的熱響應(yīng)規(guī)律,在同一升溫速率下,ANPyO的點火時間是鈍化RDX的1.6倍左右,如表2所示,且升溫速率越低,點火時間比值越大,說明ANPyO在慢速烤燃情況下熱不敏感性更好。

表2 不同升溫速率下兩種炸藥點火時間比較

隨著升溫速率的增加,鈍化RDX點火時間的變化規(guī)律與ANPyO一致,兩者變化曲線見圖4所示,可見在慢速烤燃情況下,炸藥點火時間曲線變化劇烈,隨升溫速率增加趨于平緩,這是因為殼體在緩慢升溫條件下以傳導(dǎo)傳熱和輻射傳熱的方式將熱量傳給藥柱,藥柱以傳導(dǎo)傳熱的方式將表層的熱量傳遞到藥柱內(nèi)部,同時藥柱發(fā)生自熱反應(yīng)生熱,自熱源項S隨溫度T呈指數(shù)上升,在緩慢升溫情況下,反應(yīng)生熱少,殼體傳熱慢,于是藥柱中心熱量積累過程長,點火時間長。升溫速率稍微增大時,殼體熱傳導(dǎo)量增大,炸藥自熱反應(yīng)源項S明顯增大,整個熱積累過程會明顯縮短,從而點火時間明顯縮短,而在快速升溫速率條件下,藥柱邊緣很快達到點火溫度,相應(yīng)點火時間很短。

圖4 不同升溫速率下兩種炸藥的點火時間變化曲線

隨著升溫速率增大,兩種炸藥點火位置變化規(guī)律一致,炸藥點火位置由藥柱中心移向藥柱邊緣,且點火時的殼體溫度Ts逐漸接近于點火溫度Ti,但兩種炸藥的溫度Ti和Ts相差較大,見圖5所示。

由圖5可見,在同一升溫速率下,ANPyO點火溫度Ti比鈍化RDX高出145 ℃左右,且ANPyO點火時殼體溫度Ts比鈍化RDX高出135 ℃左右,即ANPyO點火環(huán)境溫度較高,說明ANPyO更不容易被點燃,熱不敏感性較鈍化RDX要好。

圖5 兩種炸藥的溫度Ti和Ts隨升溫速率變化曲線

3 結(jié)論

1)升溫速率對點火溫度影響不大;隨著升溫速率的增加,炸藥點火時間縮短,炸藥點火位置由藥柱中心移向藥柱邊緣,炸藥點火時的中心溫度降低;隨著升溫速率的增加,殼體溫度升高,也就是點火時的環(huán)境溫度升高,說明在慢速烤燃條件下炸藥點火時的環(huán)境溫度低,危險性更大。

2)在同一升溫速率下,ANPyO的點火時間是鈍化RDX的1.6倍左右,且升溫速率越低,比值越大,說明ANPyO在慢速烤燃情況下熱鈍感性更好。

3)在同一升溫速率下,兩種炸藥的點火溫度和點火時的殼體溫度(環(huán)境溫度)相差較大,ANPyO點火溫度比鈍化RDX高出145 ℃左右,且ANPyO點火時的殼體溫度比鈍化RDX高出135 ℃左右,即ANPyO點火環(huán)境溫度較高,說明ANPyO更不容易被點燃,熱不敏感性較鈍化RDX要好。

[1] 智小琦, 胡雙啟, 肖志華, 等. 密閉條件對鈍化RDX快速烤燃響應(yīng)特性的影響 [J]. 火炸藥學(xué)報, 2010, 33

(1): 31-34.

[2] Victor A C. Simple calculation methods for munitions cook-off times and temperature [J]. Propellant, Explosion, Pyrotechnical, 1995, 20(5): 252-259.

[3] 李金山, 黃奕剛, 董海山, 等. 多硝基吡啶的密度泛函理論研究 [J]. 含能材料, 2003, 11(4): 178-181.

[4] Hollins R A, Merwin L M, Nissan R A. Aminoni-tropyridines and their N-O xides [J]. Heterocycl Chem, 1996, 33: 895-904.

[5] Licht H H. Performance and sensitivity of explosives [J]. Propellants, Explosives Pyrotechnics, 2000, 25(3): 126-132.

[6] 成健, 姚其正, 劉祖亮, 等. 2, 6-二氨基-3, 5-二硝基吡啶-1-氧化物的合成與性能 [J]. 含能材料, 2008, 16(6): 672-675.

[7] 成健, 姚其正, 劉祖亮. 2, 6-二氨基-3, 5-二硝基吡啶-1-氧化物的合成新方法 [J]. 含能材料, 2009, 17(2): 166-168.

[8] 何志偉. 多氨基多硝基吡啶氮氧化物及其配方的性能研究 [D]. 南京: 南京理工大學(xué), 2010.

[9] 劉華寧, 鄭宇, 邱從禮, 等. 新型炸藥2, 6-二氨基-3, 5-二硝基吡啶-1-氧化物的射流沖擊感度實驗研究 [J]. 含能材料, 2014, 22(3): 337-342.

[10] 程波, 李文彬, 鄭宇, 等. 不同約束條件下ANPyO炸藥快烤試驗研究 [J]. 爆破器材, 2013, 42(5): 53-56.

[11] 荊松吉. 凝聚炸藥烤燃機理研究及二維數(shù)值模擬 [D]. 長沙: 國防科技大學(xué), 2004.

[12] 張亞坤, 智小琦, 李強, 等. 烤燃溫度對凝聚炸藥熱起爆臨界溫度影響的研究 [J]. 彈箭與制導(dǎo)學(xué)報, 2014, 34(1): 69-72.

[13] 馮長根. 熱爆炸理論 [M]. 北京: 科學(xué)出版社, 1988.

[14] 張亞坤, 智小琦, 李強, 等. RDX基炸藥熱起爆臨界溫度的測試及數(shù)值計算 [J]. 火炸藥學(xué)報, 2014, 37(1): 39-43.

The Simulation Research of Thermal Stimulus Response of New Explosive ANPyO

LUO Dan1,ZHENG Yu1,LIU Xiaojun2,WANG Yong3,LIU Jun1

(1 School of Mechanical Engineering, NUST, Nanjing 210094, China; 2 Liaosheg Industries Group Co. Ltd, Shenyang 110045, China; 3 Beijing General Research Institute of Mining and Metallurgy, Beijing 100015, China)

In order to analyze new passive explosive ANPyO’s thermal response regularity, the cook-off model of shelled explosive was established based on software Fluent and it was used for simulating the explosive ANPyO’s cook-off situation at different heating rate (1 ℃/min, 2 ℃/min, 10 ℃/min, 20 ℃/min, 80 ℃/min, 140 ℃/min), and compared with passive RDX. The result shows that:with the heating rate increases, explosive’s ignition time decreases. Under the same heating rate, ANPyO’s ignition time is 1.6 times of passive RDX, and the ignition temperature of ANPyO is higher than passive RDX about 145 ℃. Therefore, ANPyO’s thermal stimulus sensitivity is insensitive to passive RDX.

insensitive explosive; ANPyO explosive; cook-off; ignition time; numerical simulation

2014-12-16

國家自然科學(xué)基金(11302108)資助

羅丹(1992-),女,湖南常德人,碩士研究生,研究方向:不敏感彈藥技術(shù)。

TJ55;O64;O389

A