范一清，王 炅

(南京理工大學(xué)，江蘇 南京 210094)

引信傳爆序列殉爆特性數(shù)值模擬方法

范一清，王 炅

(南京理工大學(xué)，江蘇南京210094)

針對(duì)殉爆數(shù)值模擬方法的研究對(duì)象主要為炸藥的問題，提出了引信傳爆序列殉爆特性數(shù)值模擬方法。該方法采用Autodyn軟件對(duì)典型引信傳爆序列殉爆特性進(jìn)行了數(shù)值模擬，得到了主發(fā)傳爆序列傳播到被發(fā)傳爆序列周圍的爆炸沖擊波。數(shù)值模擬結(jié)果表明，由于主發(fā)引信傳爆序列爆炸沖擊波的瞬間沖擊作用，造成被發(fā)引信傳爆序相應(yīng)質(zhì)點(diǎn)壓力劇增達(dá)到臨界爆發(fā)點(diǎn)，最終導(dǎo)致被發(fā)引信傳爆序列發(fā)生殉爆。

引信；數(shù)值模擬；傳爆序列；殉爆

0 引言

引信傳爆序列是爆轟能量傳遞的關(guān)鍵部件。傳爆序列的安全性直接關(guān)系著彈藥的安全性和可靠性[1]。傳爆序列中的含能材料通常具備了較高的機(jī)械感度和沖擊波感度，在外界刺激作用下極易發(fā)生爆炸或爆轟，導(dǎo)致災(zāi)難性事故發(fā)生[2]。

殉爆是含能材料具備的危險(xiǎn)性特征之一。通過引爆殉爆源，觀察在殉爆源爆炸能量作用下，是否引起被測(cè)試物品發(fā)生殉爆現(xiàn)象，從而判定在生產(chǎn)，貯存，運(yùn)輸，銷毀過程中是否存在危險(xiǎn)[3]。

當(dāng)引信采用分裝方式時(shí)，一個(gè)包裝箱中數(shù)發(fā)或十幾發(fā)引信的同時(shí)爆炸也具有極大的危害，還有可能引起附近彈藥的連鎖反應(yīng)。目前，國內(nèi)外對(duì)炸藥殉爆仿真研究較多。研究引信在爆轟環(huán)境中是否發(fā)生殉爆反應(yīng)，可參照彈藥殉爆仿真過程對(duì)引信進(jìn)行仿真分析。陳朗等人采用非線性有限元計(jì)算方法對(duì)炸藥殉爆實(shí)驗(yàn)進(jìn)行了數(shù)值模擬計(jì)算，分析了炸藥殉爆過程中，被發(fā)炸藥爆轟波的成長(zhǎng)歷程[4]。周保順等人對(duì)沖擊波作用下炸藥殉爆過程進(jìn)行了數(shù)值模擬研究[5]。魯忠寶等人采用ANSYS/LS-DYNA軟件建立了典型裝藥水下殉爆的有限元仿真模型[6]。目前尚未見有詳細(xì)描述引信傳爆序列殉爆數(shù)值仿真的文獻(xiàn)。本文針對(duì)此問題，提出了引信傳爆序列殉爆特性數(shù)值模擬方法。

1 計(jì)算模型

1.1 物理模型及算法

引信傳爆序列由雷管、導(dǎo)爆藥、傳爆藥組成。類比彈藥殉爆[7]過程可得，引信殉爆模型應(yīng)包括主發(fā)引信殼體和傳爆序列，被發(fā)引信殼體和傳爆序列，空氣。如圖1所示。

主發(fā)引信傳爆序列含能材料為長(zhǎng)方形物理模型?？紤]到模型的對(duì)稱性，實(shí)際建模為1/2長(zhǎng)方形。圖1中，主發(fā)引信傳爆序列中雷管，導(dǎo)爆藥，傳爆藥長(zhǎng)度依次為10.3 mm,5.15 mm和18 mm。

采用的算法為單點(diǎn)多物質(zhì)流固耦合算法，即計(jì)算模型中用歐拉法描述主發(fā)引信傳爆序列裝藥及周圍空氣介質(zhì)。被發(fā)引信體用拉格朗日法描述。主發(fā)引信傳爆序列中的雷管上端面中心點(diǎn)設(shè)為起爆點(diǎn)。

1.2 材料參數(shù)和物態(tài)方程

雷管常采用分壓多層裝藥，選藥主要考慮簡(jiǎn)化傳爆序列的模型，在接近實(shí)際雷管能量的情況下，選用B炸藥。

導(dǎo)爆藥受其結(jié)構(gòu)的限制，選擇臨界直徑較小的黑索金。

傳爆藥選用具有與特屈兒相同感度的PBX9404。

主發(fā)引信傳爆序列含能材料包括雷管裝藥、導(dǎo)爆藥、傳爆藥。均選用常用的Jones-Wilkins-Lee(JWL)狀態(tài)方程。其基本形式如下：

表1給出了主發(fā)引信傳爆序列中各含能材料的爆轟性能參數(shù)和JWL狀態(tài)方程參數(shù)[8-10]。

表1 JWL模型參數(shù)

被發(fā)引信傳爆序列含能材料選用點(diǎn)火增長(zhǎng)模型。其基本形式如下：

G2(1-λ)eλgpz

表2給出了被發(fā)引信傳爆序列各含能材料的反應(yīng)模型參數(shù)[11-13]。

表2 點(diǎn)火增長(zhǎng)反應(yīng)模型參數(shù)

采用彈塑性模型，用Johnson-cook和Shock方程描述引信殼體的本構(gòu)關(guān)系和沖擊波作用下材料的狀態(tài)變化，用失效應(yīng)變來判斷材料的失效，密度為790 kg/cm3，彈性模量206 GPa，屈服應(yīng)力350 MPa，殼體的失效應(yīng)變?yōu)?.15?？諝獠捎镁€性多項(xiàng)式方程和理想氣體的狀態(tài)方程，密度取標(biāo)準(zhǔn)狀態(tài)下ρ=1.298 kg/m3。

2 模擬結(jié)果分析

上述主、被發(fā)引信傳爆序列數(shù)值模擬結(jié)果如下圖2所示。由圖可知主發(fā)引信傳爆序列發(fā)生可靠作用后，被發(fā)引信發(fā)生了殉爆。根據(jù)被發(fā)引信體中爆轟波傳播方向，選取不同位置單元點(diǎn)觀察爆轟波成長(zhǎng)情況，分別在被發(fā)引信體對(duì)應(yīng)的雷管，導(dǎo)爆藥和傳爆藥位置進(jìn)行標(biāo)記1，2，3。如圖3所示。各節(jié)點(diǎn)的壓力時(shí)程曲線如圖4,圖5,圖6所示。

3 模擬結(jié)果驗(yàn)證

對(duì)于炸藥殉爆來說，沖擊波壓力大于臨界起爆壓力是炸藥殉爆的必要條件。同理，對(duì)于引信殉爆也采取同樣的分析方式。炸藥B所需的起爆壓力至少要達(dá)到4.950 495 MPa，黑索金起爆壓力至少達(dá)到5.630 027 MPa，PBX-9404至少達(dá)到5.543 478 MPa。進(jìn)行試驗(yàn)驗(yàn)證，由表3可知試驗(yàn)中被發(fā)引信體的雷管，導(dǎo)爆藥，傳爆藥即對(duì)應(yīng)的節(jié)點(diǎn)1,2,3處的沖擊波超壓峰值均大于臨界起爆壓力。所以被發(fā)引信發(fā)生殉爆反應(yīng)與數(shù)值模擬結(jié)果一致。

節(jié)點(diǎn)名稱沖擊波超壓峰值大小/GPa13.344321.284431.0885

4 結(jié)論

本文提出了引信傳爆序列殉爆特性數(shù)值模擬方法。該方法利用AUTODYN軟件進(jìn)行引信傳爆序列殉爆仿真。數(shù)值模擬結(jié)果表明，引信傳爆序列在遭遇意外刺激發(fā)生可靠作用時(shí)會(huì)使鄰近引信產(chǎn)生殉爆反應(yīng)。所以為適應(yīng)不敏感引信技術(shù)的要求，急需設(shè)計(jì)新型的引信傳爆序列，含能材料和包裝技術(shù)，減小引信載體(包括人員)的損害程度。

[1]張合, 李豪杰. 引信機(jī)構(gòu)學(xué)[M]. 北京：北京理工大學(xué)出版社, 2014.

[2]張春海. 不敏感彈藥,讓士兵和武器更安全[J]. 現(xiàn)代軍事, 2006(2):54-59.

[3]鄒金龍，雷雅茹. 不敏感彈藥對(duì)引信技術(shù)的要求內(nèi)涵[J]. 探測(cè)與控制學(xué)報(bào)，2016, 38(1):1-6.

[4]陳朗, 王晨, 魯建英,等. 炸藥殉爆實(shí)驗(yàn)和數(shù)值模擬[J]. 北京理工大學(xué)學(xué)報(bào), 2009, 29(6):497-500.

[5]周保順, 王少龍, 徐明利,等. 非均質(zhì)炸藥殉爆試驗(yàn)數(shù)值模擬[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2009, 29(5):145-148.

[6]魯忠寶, 胡宏偉, 劉銳,等. 典型裝藥水下爆炸的殉爆規(guī)律研究[J]. 魚雷技術(shù), 2014, 22(3):230-235.

[7]王翔, 盧校軍, 黃毅民,等. 殉爆試驗(yàn)研究[C]//第九屆全國爆炸與安全技術(shù)學(xué)術(shù)會(huì)議論文集. 北京：中國兵工學(xué)會(huì)，2006.

[8]Balagansky I A, Stepanov A A. Numerical simulation of composition B high explosive charge desensitization in gap test assembly after loading by precursor wave[J]. Shock Waves, 2016, 25(2):109-115.

[9]趙錚, 閆鴻浩, 李曉杰,等. 應(yīng)用LS-DYNA進(jìn)行爆炸壓實(shí)爆轟過程的數(shù)值模擬[C]//第三屆全國計(jì)算爆炸力學(xué)會(huì)議論文集. 北京：中國力學(xué)學(xué)會(huì)爆炸力學(xué)專業(yè)委員會(huì)，2006.

[10]梁增友, 黃風(fēng)雷, 張震宇. PBX-9404的化學(xué)反應(yīng)速率方程及起爆特性[J]. 爆炸與沖擊, 2008, 28(1):38-43.

[11]賈憲振, 楊建, 陳松,等. 帶殼B炸藥在鎢珠撞擊下沖擊起爆的數(shù)值模擬[J]. 火炸藥學(xué)報(bào), 2010, 33(5):43-47.

[12]楊剛, 傅奕軻, 鄭建民,等. 基于FE-SPH自適應(yīng)耦合算法對(duì)長(zhǎng)桿彈侵徹反應(yīng)裝甲的模擬研究[J]. 工程力學(xué), 2016, 33(1):223-231.

[13]楊剛, 傅奕軻, 鄭建民,等. 基于FE-SPH自適應(yīng)耦合算法對(duì)長(zhǎng)桿彈侵徹反應(yīng)裝甲的模擬研究[J]. 工程力學(xué), 2016, 33(1):223-231.

NumericalSimulationofSympatheticDetonationofFuzeExplosiveTrains

FAN Yiqing,WANG Jiong

(Nanjing University of Science & Technology, Nanjing 210094,China)

According to the situation that the study object of numerical simulation of sympathetic detonation was usually explosives,this article present the numerical simulation of sympathetic detonation of explosive trains in fuze. The autodyn was used to research sympathetic detonation property of explosive trains in fuze, and the shock waves produced by the donor fuze explosive trains propagated to the acceptor was cut. Simulation results turned out that under the blast impact of the donor fuze explosive trains, which led to the acceptor’s pressure excursion, reaching the critical point, which finally caused the sympathetic detonation of the fuze explosive train.

fuze; numerical simulation; explosive train; sympathetic detonation

2017-05-06

范一清(1993—)，女，安徽桐城人，博士研究生，研究方向：不敏感引信技術(shù)。E-mail:fyq930222@sina.com。

TJ430;O383

A

1008-1194(2017)05-0007-03