陳朗，馬欣，黃毅民，伍俊英，常雪梅

(1.北京理工大學(xué) 機(jī)電學(xué)院，北京100081;2.中國(guó)工程物理研究院 化工材料研究所，四川 綿陽(yáng)621900)

0 引言

烤燃實(shí)驗(yàn)是研究炸藥熱安全性主要方法之一。在烤燃實(shí)驗(yàn)中觀測(cè)炸藥溫度變化，能夠獲得炸藥熱反應(yīng)狀態(tài)，從而為炸藥熱安全性分析提供依據(jù)。1982年，Kent 等［1］對(duì)固體推進(jìn)劑進(jìn)行了電加熱下的烤燃實(shí)驗(yàn)，并使用熱電偶測(cè)量了藥柱中心溫度;1984年，Pakulak 等［2］設(shè)計(jì)了小型烤燃彈實(shí)驗(yàn)(簡(jiǎn)稱SCB實(shí)驗(yàn))，采用可以控制加熱速率的電加熱帶，對(duì)烤燃彈殼體加熱和用熱電偶測(cè)量藥柱表面溫度;1991年，Jones 等［3］對(duì)小型烤燃彈實(shí)驗(yàn)裝置進(jìn)行了改進(jìn)，增加了熱電偶測(cè)量炸藥內(nèi)部溫度。在小型烤燃彈實(shí)驗(yàn)中放置熱電偶數(shù)量有限，要全面反映炸藥內(nèi)部溫度變化有技術(shù)困難。2002年，Michael 等［4］設(shè)計(jì)了熱點(diǎn)火實(shí)驗(yàn)裝置，采用在兩塊炸藥之間放置微型熱電偶的方法，對(duì)多個(gè)熱電偶在炸藥中放置位置進(jìn)行精確定位，提高了溫度測(cè)量精度，同時(shí)，也減小了熱電偶對(duì)炸藥溫度的影響。在炸藥烤燃實(shí)驗(yàn)研究的同時(shí)，人們也在對(duì)烤燃炸藥進(jìn)行數(shù)值模擬研究。炸藥烤燃數(shù)值模擬一方面可以得到實(shí)驗(yàn)沒(méi)有測(cè)量的性能參量，另一方面還可以很方便地改變條件，對(duì)不同烤燃情況進(jìn)行預(yù)測(cè)分析。1981年，McGuire 等［5］提出了用多步化學(xué)分解模型和熱傳導(dǎo)模型來(lái)分析炸藥的烤燃過(guò)程，并用化學(xué)傳熱軟件TOPAZ 計(jì)算了炸藥點(diǎn)火時(shí)間。1997年，Chidester 等［6］對(duì)HMX 和TATB為基的幾種炸藥烤燃過(guò)程進(jìn)行了二維數(shù)值模擬計(jì)算。2002年，Michael 等［4］用Coyote 軟件，對(duì)PBXN-109 和PBX-9501 炸藥的熱點(diǎn)火實(shí)驗(yàn)進(jìn)行了三維模擬計(jì)算。2006年，McClelland 等［7］用ALE3D 軟件對(duì)LX-10(HMX94.7%，氟橡膠5.3%)炸藥烤燃實(shí)驗(yàn)進(jìn)行了數(shù)值模擬，采用多步反應(yīng)的化學(xué)動(dòng)力學(xué)模型，并且考慮了鋼殼和炸藥柱之間空氣層的熱傳導(dǎo)和輻射。從目前研究情況看，采用實(shí)驗(yàn)和數(shù)值模擬計(jì)算相結(jié)合的方法，能夠更有效地對(duì)烤燃條件下炸藥熱反應(yīng)規(guī)律進(jìn)行全面和深入的分析。

本文采用多點(diǎn)測(cè)溫的炸藥烤燃實(shí)驗(yàn)裝置，對(duì)PBXC10 炸藥進(jìn)行了不同加熱速率下的烤燃實(shí)驗(yàn)。通過(guò)測(cè)量炸藥內(nèi)部不同位置的溫度變化，分析了炸藥熱反應(yīng)狀態(tài)。建立了炸藥烤燃實(shí)驗(yàn)三維計(jì)算模型，根據(jù)一種加熱速率下PBXC10 炸藥的實(shí)驗(yàn)溫度與時(shí)間曲線，標(biāo)定炸藥指前因子和活化能，對(duì)其它加熱速率下的炸藥熱反應(yīng)過(guò)程進(jìn)行了數(shù)值模擬計(jì)算分析，并對(duì)PBXC10，JB9014 和JOB9003 炸藥在不同溫度下的熱爆炸進(jìn)行了對(duì)比。

1 炸藥烤燃實(shí)驗(yàn)

采用與熱點(diǎn)火實(shí)驗(yàn)裝置［4］相似的多點(diǎn)溫度測(cè)量裝置進(jìn)行炸藥烤燃實(shí)驗(yàn)。圖1為烤燃實(shí)驗(yàn)裝置簡(jiǎn)圖。烤燃裝置的殼體及炸藥均由上、下兩部分組成。藥柱直徑為50 mm，高為50 mm?？救佳b置殼體內(nèi)徑為50 mm，內(nèi)部高為50 mm，外徑為100 mm，外部高為80 mm，分別在炸藥與殼體接觸的頂部和底部留有φ44 mm×2 mm 的空氣槽，當(dāng)炸藥受熱膨脹時(shí)，可緩解炸藥對(duì)殼體的壓力。在殼體上開(kāi)4 個(gè)孔，安裝4 個(gè)鋼質(zhì)6 角頭的螺栓。在下部殼體及炸藥表面開(kāi)7 個(gè)凹槽，凹槽的寬度和深度均為0.5 mm，將7 根直徑為0.5 mm 的K 型熱電偶或7 根直徑為0.08 mm的微型熱電偶嵌入槽中，然后將上下殼體和炸藥扣合，并用螺栓和螺帽固定，并在兩個(gè)殼體間涂抹高溫膠密封粘結(jié)。采用上下兩個(gè)陶瓷加熱套對(duì)烤燃彈加熱。圖2為烤燃彈上下兩部分扣合前的實(shí)物照片。圖3為7 根熱電偶在下半塊炸藥上對(duì)應(yīng)7個(gè)監(jiān)測(cè)點(diǎn)的位置分布示意圖。圖4為烤燃彈上下兩部分扣合后的實(shí)物照片。

圖1 烤燃實(shí)驗(yàn)裝置簡(jiǎn)圖Fig.1 Schematic diagram of cook-off test setup

圖2 烤燃彈上下兩部分扣合前的照片F(xiàn)ig.2 Photograph of upper and lower parts of cook-off bomb before assembling

由控溫儀控制加熱套的升溫速率，熱電偶的電信號(hào)通過(guò)導(dǎo)線傳輸?shù)接涗泝x，記錄炸藥內(nèi)部各點(diǎn)溫度變化。實(shí)驗(yàn)中分別采用9.1 K/min，5.0 K/min，1.1 K/min，0.8 K/min 和0.5 K/min 5 種加熱速率，加熱炸藥直至發(fā)生劇烈反應(yīng)。實(shí)驗(yàn)炸藥為PBXC10炸藥(TATB/HMX)。

2 炸藥烤燃數(shù)值模擬計(jì)算

根據(jù)炸藥烤燃實(shí)驗(yàn)裝置，建立炸藥烤燃三維數(shù)值計(jì)算模型。模型中主要考慮炸藥、殼體和空氣槽，其中殼體外壁為加熱邊界，以替代外部加熱器的作用，炸藥被認(rèn)為是固體，自熱反應(yīng)遵循阿倫尼烏斯定律，不考慮炸藥晶型轉(zhuǎn)變和氣體釋放。

質(zhì)量，動(dòng)量、能量的輸運(yùn)方程都可以用下面的通用形式表達(dá)［8］:

圖3 熱電偶分布Fig.3 Arrangement of thermocouples

圖4 烤燃彈上下兩部分扣合后的照片F(xiàn)ig.4 Photograph of upper and lower parts of cook-off bomb after assembling

式中:φ 是通用的變量，代表質(zhì)量、動(dòng)量、能量等;ρf代表流體密度;Г 是通用的擴(kuò)散系數(shù);S 代表炸藥自熱反應(yīng)放熱源項(xiàng)，用Arrhenius 方程來(lái)表述［9］:

式中:ρe為炸藥元件的密度(kg/m3)，Q 為反應(yīng)熱(J/kg)，Z 為指前因子(s-1)，α 為反應(yīng)分?jǐn)?shù);E 為活化能(J/mol);R 為普適氣體常數(shù)(J/(mol·K)。

炸藥與鋼殼間的空氣區(qū)域采用P1 輻射模型［10］。對(duì)于輻射熱流qr，采用如下方程:

式中:a 為吸收系數(shù);σs為散射系數(shù);G 為入射輻射;C 為線性各相異性相位函數(shù)系數(shù)。

由于烤燃裝置為軸對(duì)稱結(jié)構(gòu)，建立1/4 模型。模型采用尺寸為2 mm 的6 面體網(wǎng)格。采用FLUENT 軟件對(duì)炸藥烤燃過(guò)程進(jìn)行數(shù)值模擬計(jì)算［10］，炸藥自熱反應(yīng)源通過(guò)編寫子程序加載到Fluent 軟件中。

在計(jì)算過(guò)程中，鋼殼側(cè)壁為加熱面;炸藥和鋼殼之間被認(rèn)為是耦合的熱傳導(dǎo)界面，在耦合界面上兩種物質(zhì)的溫度及熱流連續(xù);炸藥與鋼殼間的空氣區(qū)域通過(guò)輻射傳熱，將熱量傳遞給炸藥。計(jì)算中根據(jù)實(shí)驗(yàn)測(cè)量的殼體表面溫度曲線擬合外熱源加熱速率方程，以子函數(shù)的形式加載程序中作為加熱條件。

確定活化能和指前因子是炸藥烤燃數(shù)值模擬計(jì)算的關(guān)鍵所在。計(jì)算中首先以5.0 K/min 加熱速率下，7 號(hào)熱電偶溫度與時(shí)間曲線為基準(zhǔn)，帶入一組指前因子和活化能進(jìn)行計(jì)算，將計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行比較，反復(fù)調(diào)整指前因子和活化能，直至得到與實(shí)驗(yàn)結(jié)果相吻合的計(jì)算結(jié)果。然后把不同加熱速率下計(jì)算的溫度與時(shí)間曲線與實(shí)驗(yàn)曲線進(jìn)行對(duì)比，以驗(yàn)證計(jì)算的正確性。

3 結(jié)果與討論

3.1 實(shí)驗(yàn)結(jié)果

表1 不同加熱速率下炸藥的點(diǎn)火時(shí)間和點(diǎn)火時(shí)的壁面溫度Tab.1 Ignition time of explosive and shell temperature at igniting at different heating rates

圖5分別是3 種加熱速率下炸藥劇烈反應(yīng)后烤燃裝置的照片。從烤燃裝置變形和殘留炸藥反應(yīng)產(chǎn)物看，裝置雖然被沖開(kāi)，但形狀基本完整，殼體內(nèi)殘留少量黑色的炸藥反應(yīng)產(chǎn)物，表明炸藥都是發(fā)生了燃燒反應(yīng)，而沒(méi)有發(fā)生爆轟。在裝置周圍都存在未反應(yīng)的黃色炸藥粉末，顯示炸藥沒(méi)有完全燃燒。

表1是不同加熱速率下炸藥的點(diǎn)火時(shí)間和點(diǎn)火時(shí)的壁面溫度。從表1可以看出，加熱速率為9.1 K/min時(shí)，炸藥點(diǎn)火時(shí)間為1 884 s，壁面溫度為571 K;加熱速率為0.5 K/min 時(shí)，炸藥點(diǎn)火時(shí)間為17 488 s，壁面溫度為494 K。隨著加熱速率的降低，PBXC10 炸藥的點(diǎn)火時(shí)間增長(zhǎng)，炸藥點(diǎn)火時(shí)殼體溫度(點(diǎn)火時(shí)壁面溫度)降低。

表2為在不同加熱速率下，藥柱點(diǎn)火時(shí)測(cè)量點(diǎn)的溫度。加熱速率為9.1 K/min，5.0 K/min 和1.1 K/min時(shí)，炸藥點(diǎn)火時(shí)7 號(hào)熱電偶(距離炸藥中心22 mm)溫度最高;加熱速率為0.8 K/min時(shí)，炸藥點(diǎn)火時(shí)5 號(hào)熱電偶(距離炸藥中心14.42 mm)溫度最高;加熱速率為0.5 K/min 時(shí)，炸藥點(diǎn)火時(shí)1 號(hào)熱電偶(位于炸藥中心)溫度最高。隨著加熱速率的降低，藥柱點(diǎn)火時(shí)的高溫區(qū)域逐漸向藥柱中心移動(dòng)，藥柱點(diǎn)火位置逐漸趨向于藥柱中心。

圖5 3 種加熱速率下炸藥劇烈反應(yīng)后烤燃彈照片F(xiàn)ig.5 Results of cook-off experiments at different heating rates

表2 不同加熱速率下測(cè)量點(diǎn)在點(diǎn)火時(shí)刻的溫度Tab.2 Measuring point temperatures at ignition at different heating rates

圖6是加熱速率0.5 K/min 時(shí)，藥柱內(nèi)部從200 min到250 min 溫度分布。從圖中可以看出，5 號(hào)點(diǎn)、4 號(hào)點(diǎn)、2 號(hào)點(diǎn)和1 號(hào)點(diǎn)依次出現(xiàn)溫度上升速率減慢的溫度平臺(tái)，持續(xù)時(shí)間約為7 min 30 s.5 號(hào)點(diǎn)距離殼體14.42 mm，最先出現(xiàn)溫度平臺(tái);1 號(hào)點(diǎn)位于藥柱中心，溫度上升較慢，最后出現(xiàn)溫度平臺(tái)。7 號(hào)點(diǎn)距離殼體3 mm，受到殼體溫度影響最大，沒(méi)有出現(xiàn)溫度平臺(tái)。圖7是Michael 等對(duì)PBX-9501(95%HMX)炸藥進(jìn)行465 K 恒溫加熱時(shí)測(cè)量炸藥內(nèi)部的溫度與時(shí)間曲線［4］。在約440 K 時(shí)，5 個(gè)測(cè)量點(diǎn)溫度與時(shí)間曲線出現(xiàn)平臺(tái)。他們認(rèn)為這是HMX 發(fā)生晶型轉(zhuǎn)變，從β 相→δ 相，吸熱形成了溫度平臺(tái)。本文測(cè)量的溫度平臺(tái)出現(xiàn)時(shí)，溫度也為440 K 左右。由于PBXC10 炸藥中含有40%HMX，因此，認(rèn)為這是HMX 晶型轉(zhuǎn)變吸熱的結(jié)果，另外，在加熱速率為1.1 K/min 和0.8 K/min 時(shí)也觀測(cè)到了HMX 晶型轉(zhuǎn)變吸熱形成的溫度平臺(tái)。

3.2 計(jì)算結(jié)果

圖8是加熱速率5.0 K/min 時(shí)，炸藥7 號(hào)監(jiān)測(cè)點(diǎn)計(jì)算溫度與時(shí)間曲線和實(shí)驗(yàn)曲線對(duì)比。計(jì)算曲線與實(shí)驗(yàn)曲線能夠很好的吻合。計(jì)算使用的PBXC10炸藥參數(shù)如表3所示。

表3 PBXC10 物性參數(shù)Tab.3 Parameters of PBXC10

圖6 加熱速率0.5 K/min 時(shí)PBXC10 炸藥內(nèi)部溫度分布Fig.6 Internal temperature distribution of PBXC10 explosive at heating rate of 0.5 K/min

圖7 體積膨脹率13.8%和加熱速率0.5 K/min 時(shí)PBX-9501 炸藥內(nèi)部溫度分布Fig.7 Internal temperature distribution of PBX-9501 explosive with 13.8% expansion volume at heating rate of 0.5 K/min

圖8 加熱速率5.0 K/min 時(shí)7 號(hào)監(jiān)測(cè)點(diǎn)計(jì)算與實(shí)驗(yàn)溫度曲線Fig.8 Measured and calculated T-t curves of point 7 at heating rate of 5.0 K/min

圖9和圖10分別是9.1 K/min 和0.5 K/min 加熱速率時(shí)，炸藥4 號(hào)監(jiān)測(cè)點(diǎn)計(jì)算溫度與時(shí)間和實(shí)驗(yàn)曲線對(duì)比。計(jì)算值與實(shí)驗(yàn)值基本符合。在0.5 K/min加熱速率實(shí)驗(yàn)中，為了節(jié)省實(shí)驗(yàn)時(shí)間，先對(duì)炸藥進(jìn)行快速加熱，加熱速率為5.0 K/min，當(dāng)炸藥中心溫度達(dá)到343 K 時(shí)，再以0.5 K/min 的加熱速率進(jìn)行加熱。

圖9 加熱速率9.1 K/min 時(shí)4 號(hào)監(jiān)測(cè)點(diǎn)計(jì)算和實(shí)驗(yàn)溫度曲線Fig.9 Measured and calculated T-t curves of point 4 at heating rate of 9.1 K/min

圖10 加熱速率0.5 K/min 時(shí)4 號(hào)監(jiān)測(cè)點(diǎn)計(jì)算和實(shí)驗(yàn)溫度曲線Fig.10 Measured and calculated T-t curves of point 4 at heating rate of 0.5 K/min

表4 不同加熱條件下實(shí)驗(yàn)和計(jì)算的炸藥點(diǎn)火時(shí)間以及計(jì)算的點(diǎn)火溫度Tab.4 Measured and calculated ignition time，and calculated ignition temperature of explosive at different heating rates

表4是不同加熱速率下實(shí)驗(yàn)和計(jì)算的炸藥點(diǎn)火時(shí)間以及計(jì)算的點(diǎn)火溫度。表中數(shù)據(jù)顯示，計(jì)算點(diǎn)火時(shí)間與實(shí)驗(yàn)值相差在5%以內(nèi)。實(shí)驗(yàn)中沒(méi)有測(cè)量到炸藥點(diǎn)火點(diǎn)的點(diǎn)火溫度，只能通過(guò)計(jì)算得到點(diǎn)火溫度，從計(jì)算結(jié)果可以看出，加熱速率對(duì)炸藥的點(diǎn)火溫度影響不大。以上結(jié)果表明:根據(jù)5.0 K/min 加熱速率下實(shí)驗(yàn)測(cè)量的溫度與時(shí)間曲線標(biāo)定炸藥指前因子和活化能。利用建立炸藥烤燃計(jì)算模型和計(jì)算方法，可以準(zhǔn)確預(yù)測(cè)不同加熱速率下炸藥的烤燃熱反應(yīng)過(guò)程。

除了對(duì)實(shí)驗(yàn)的5 種加熱速率下炸藥烤燃過(guò)程進(jìn)行了計(jì)算外，還對(duì)3.3 K/h 加熱速率下炸藥烤燃進(jìn)行了計(jì)算。圖11是不同加熱速率下炸藥點(diǎn)火時(shí)刻烤燃裝置剖面溫度分布。在9.1 K/min 加熱速率下，炸藥點(diǎn)火區(qū)域在炸藥中部外緣。由于殼體加熱速率較快，直到炸藥即將點(diǎn)火時(shí)，殼體溫度一直高于炸藥內(nèi)部溫度，熱量由炸藥外部向內(nèi)部傳遞。在5.0 K/min 加熱速率下，點(diǎn)火區(qū)域移向炸藥內(nèi)部，殼體溫度仍然高于炸藥內(nèi)部溫度，熱量還是由炸藥外部向內(nèi)部傳遞。在0.8 K/min 加熱速率下，點(diǎn)火區(qū)為以炸藥中心為圓心的環(huán)狀區(qū)域，但點(diǎn)火前炸藥內(nèi)部溫度已高于殼體溫度，熱量會(huì)從炸藥向殼體傳遞，另外，由于炸藥中心溫度還低于高溫區(qū)，熱量還繼續(xù)向炸藥中心傳遞。但隨著加熱速率的進(jìn)一步減小，炸藥環(huán)狀點(diǎn)火區(qū)域的半徑會(huì)逐漸縮小。在0.5 K/min加熱速率下，高溫的點(diǎn)火區(qū)域?yàn)檎ㄋ幹行奶幍那蛐螀^(qū)域，而在3.3 K/h 加熱速率下球形點(diǎn)火區(qū)尺寸進(jìn)一步減小，熱量由炸藥中心向外傳遞。

圖11 不同加熱速率下點(diǎn)火時(shí)刻烤燃裝置剖面溫度分布Fig.11 Temperature distributions on bomb section at different heating rates

在炸藥烤燃研究中，關(guān)于炸藥快速烤燃和慢速烤燃判斷，一直沒(méi)有明確的定義和方法。一般認(rèn)為，加熱速率快就是快速烤燃，加熱速率慢就是慢速烤燃，但如何定義快慢加熱速率的分界點(diǎn)是一難題，并且一直沒(méi)有明確的判斷方法。根據(jù)以上的研究結(jié)果，以炸藥烤燃過(guò)程中熱量傳遞方向?yàn)橐罁?jù)，把烤燃分為快速、中速和慢速烤燃3 種。在炸藥烤燃中，在炸藥點(diǎn)火之前，外界溫度大于炸藥內(nèi)任何一點(diǎn)，熱量總是由外界向炸藥內(nèi)部傳遞，稱為快速烤燃;炸藥任何一定點(diǎn)溫度高于外界溫度，同時(shí)高于炸藥中心點(diǎn)溫度，熱量同時(shí)向外界和炸藥內(nèi)部傳遞，稱為中速烤燃;炸藥中心點(diǎn)溫度高于其他點(diǎn)和外界溫度，熱量從炸藥中心向外部傳遞，稱為慢速烤燃。根據(jù)此定義，對(duì)于本文烤燃實(shí)驗(yàn)，在9.1 K/min 和5.0 K/min 加熱速率為炸藥快速烤燃，1.1 K/min 和0.8 K/min 加熱速率為中速烤燃，而0.5 K/min 和3.3 K/h 加熱速率為慢速烤燃。

利用烤燃實(shí)驗(yàn)計(jì)算模型，分別對(duì)PBXC10 炸藥、JB9014 炸藥和JOB9003 炸藥，在不同環(huán)境溫度下熱反應(yīng)過(guò)程進(jìn)行數(shù)值模擬。JB9014 和JOB9003 炸藥計(jì)算參數(shù)分別來(lái)自文獻(xiàn)［11-12］。圖12是計(jì)算的3 種炸藥點(diǎn)火時(shí)間與環(huán)境溫度倒數(shù)1/Ta的關(guān)系?？梢钥吹?，在相同加熱條件下，JB9014 炸藥點(diǎn)火延遲時(shí)間最長(zhǎng)，JOB9003 最短，PBXC10 居中。PBXC10炸藥熱感度在JB9014 炸藥和JOB9003 炸藥之間。

圖12 計(jì)算3 種炸藥點(diǎn)火時(shí)間t 與環(huán)境溫度倒數(shù)1/Ta 的關(guān)系Fig.12 Relation between ignition time of explosive and 1/Ta

4 結(jié)論

建立了炸藥烤燃實(shí)驗(yàn)多點(diǎn)溫度測(cè)量方法，能夠?qū)φㄋ巸?nèi)部測(cè)量點(diǎn)進(jìn)行準(zhǔn)確定位，記錄從炸藥外沿到炸藥中心不同位置的溫度變化歷程。觀測(cè)到了PBXC10 炸藥中HMX 晶型轉(zhuǎn)變吸熱引起的溫度細(xì)微變化。建立了炸藥烤燃計(jì)算模型和計(jì)算方法，可以準(zhǔn)確預(yù)測(cè)不同加熱條件下炸藥的熱反應(yīng)過(guò)程。提出了快速烤燃，中速烤燃和慢速烤燃，以及按照炸藥內(nèi)部熱量傳遞方向區(qū)分烤燃類型的方法。PBXC10炸藥熱感度在JB9014 炸藥和JOB9003 炸藥之間。

References)

［1］ Kent R，Rat M.Explosion thermique(cook-off)des propergols solides［J］.Propellant，Explosion，Pyrotechnics，1982，7:129-135.

［2］ Pakulak J M.USA small-scale cook-off bomb (SCB)test［C］∥Minutes of 21th Department of Defense Explosives Safety Board Explosives Safety Seminar.Houston:Defense Technical Information Center，1984.

［3］ Jones D A，Parker R P.Heat flow calculations for the small-scale cook-off bomb test，AD-A236829［R］.US:DTIC，1991.

［4］ Kaneshige M J，Renlund A M，Schmitt R G，et al.Cook-off experiments for model validation at sandia national laboratories［C］∥Proceeding of the 12th International Detonation Symposium.Norfolk，VA:Naval Surface Warfare Center and Lawrence Livermore National Laboratory，2002

［5］ McGuire R R，Tarver C M.Chemical decomposition models for the thermal explosion of confined HMX，TATB，RDX and TNT explosives［C］∥The 7th International Detonation Symposium，Anapolis，Maryland，US:Office of Naval Research，1981.

［6］ Chidester S K，Tarver C M，Green L G，et al.On the violence of thermal explosion in solid explosives［J］.Combustion and Flame，1997，110:264-280.

［7］ McCelland M A，Maienschein J L，Howard W M，et al.ALE3D simulation of heating and violence in a fast cookoff experiment with LX-10［C］∥The 13th International Detonation Symposium.Norfolk，VA，US:Naval Surface Warfare Center，2006.

［8］ 王福軍.計(jì)算流體動(dòng)力學(xué)分析［M］.北京:清華大學(xué)出版社，2004.WANG Fu-jun.Analysis of computational fluid dynamics［M］.Beijing:Tsinghua University Press，2004.(in Chinese)

［9］ WANG Pei，CHEN Lang，WANG Yan，et al.Numerical simulation of explosive cook-off at different heating rate［C］∥2007 International Autumn Seminar on Propellants，Explosives，Pyrotechnics.Xi'an:China Ordnance Society and Beijing Institute of Technology，2007.

［10］ Fluent Inc.FLUENT user’s guide［M］.US:Fluent Inc，2006.

［11］ 任艷.炸藥在高溫環(huán)境下的熱反應(yīng)特性研究［D］.北京:北京理工大學(xué)，2008.REN Yan.Investigation on explosive thermal character under high temperature［D］.Beijing:Beijing Institute of Technology，2008，(in Chinese)

［12］ 代曉淦，呂子劍，申春迎，等.兩種尺寸JOB9003 炸藥慢速烤燃試驗(yàn)及數(shù)值模擬研究［C］∥第八屆全國(guó)爆炸力學(xué)學(xué)術(shù)會(huì)議論文集，江西，吉安:中國(guó)力學(xué)學(xué)會(huì)，2007:29-32.DAI Xiao-gan，Lü Zi-jian，SHEN Chun-ying，et al.Slow cookoff test and numerical simulation for two sizes JOB9003 explosive［C］∥The 8th National Explosion Mechanics Symposium.Ji'an，Jiangxi:Chinese Society of Theoretical and Applied Mechanics，2007:29-32.(in Chinese)