HMX基PBX炸藥熱損傷的數值計算與實驗研究

韋興文，吳束力，唐 興

（中國工程物理研究院化工材料研究所，四川綿陽621900）

HMX基PBX炸藥熱損傷的數值計算與實驗研究

韋興文，吳束力，唐 興

（中國工程物理研究院化工材料研究所，四川綿陽621900）

基于二維平面軸對稱模型，采用有限元軟件ANSYS，對HMX基PBX在溫度沖擊載荷作用下的溫度場和應力場進行數值計算，為了驗證有限元模擬分析的可靠性，對HMX基PBX藥柱進行了溫度沖擊試驗，用熱電偶和聲發射技術測量藥柱表面的溫度及熱損傷。數值計算結果表明，在降溫階段，藥柱側面中部受到較大軸向拉應力作用而產生熱損傷。實驗結果與計算結果吻合較好，說明HMX基PBX熱損傷破壞方式為拉應力破壞，抗拉強度可以較好地反映材料的耐熱損傷能力。

材料科學；HMX；高聚物黏結炸藥；PBX；熱損傷；數值模擬；有限元分析

引 言

PBX是炸藥晶體與黏結劑組成的二元復合材料，炸藥晶體與黏結劑之間存在黏結界面，在貯存、運輸和使用時會經歷各種溫度環境。研究表明［1-3］，在熱載荷作用下，尤其是降溫過程中，PBX容易發生界面損傷，首先是材料的微觀界面結構發生變化，產生微小的缺陷，隨著外界載荷的變化，缺陷會進一步擴大，最終可能導致材料的宏觀破壞。尤其是熱沖擊載荷產生的熱應力容易使PBX內部產生微裂紋，甚至開裂失效［4-11］。熱沖擊載荷條件是最危險，也是最難研究的典型工況之一，因此研究PBX的熱沖擊損傷具有重要意義。

本研究從炸藥裝藥的實際使用環境條件出發，利用電阻應變片測量、超聲波檢測和聲發射技術測試了HMX基PBX藥柱在溫度沖擊下的熱應變及熱損傷，采用有限元方法研究溫度沖擊載荷作用下HMX基PBX熱損傷規律，為HMX基PBX的工程應用提供參考。

1 計算模型

1.1 二維瞬態熱傳導微分方程

由于藥柱的形狀是軸對稱的圓柱形，可以將實際的三維空間問題簡化為二維平面軸對稱問題，溫度沖擊作用可以認為是一個均勻、各向同性且無內熱源的二維瞬態熱傳導過程，其偏微分方程無量綱形式為：

要確定上述微分方程的解，還必須給定初始條件和邊界條件。在溫度沖擊試驗中，藥柱初始溫度為T0，且分布均勻。

藥柱只與周圍空氣之間進行對流換熱，其微分方程為：

式中：h為空氣對流換熱系數；Tf為周圍空氣溫度，且為恒定值。

1.2 熱應力的計算模型

采用間接法計算熱應力：先進行傳熱分析，然后將求得的節點溫度作為載荷進行結構靜應力分析計算熱應力。在結構靜應力分析當中，有限元本構方程可簡化成如下形式：

式中：｛K｝為結構剛度矩陣；｛U｝為位移向量；｛F｝為力向量，包括集中力、溫度載荷及慣性載荷等。

采用ANSYS中Plane77平面單元作為熱分析單元，Plane82平面單元進行熱應力分析。在每一個時間步長或載荷步長內，當溫度場被計算出來后，將溫度場作為熱載荷來計算模型的熱應力。

2 結果與討論

2.1 溫度分布

用有限元方法計算得到高溫沖擊過程中藥柱中心和表面溫度隨時間變化的曲線，見圖1。

圖1 高溫沖擊過程藥柱溫度隨時間的變化曲線Fig.1 The temperature variation with the time during high temperature impact test

由圖1可見，在高溫沖擊初期，藥柱表面溫度上升較快，大約在300s時表面溫度（Ts）即達到33℃，內部溫度梯度很大，表面溫度與中心溫度（Tc）的差值達到最大值，約為12℃；隨著時間的增加，由于內部熱量的傳遞，藥柱表面的升溫速率明顯降低，溫度梯度逐漸減小，約7 200s后藥柱中心溫度達到平衡溫度60℃。高溫沖擊300s時藥柱溫度分布見圖2。

圖2 高溫沖擊300s時藥柱溫度分布圖Fig.2 Contour of temperature distribution of PBX cylinder in high temperature impact test（t＝300s）

由圖2可見，在300s時，藥柱端面邊緣點溫度達到33℃，而藥柱中心的溫度為21℃，溫差達12℃，最大溫度梯度出現在藥柱側面中部。

低溫沖擊過程中藥柱中心和表面溫度隨時間的變化曲線見圖3。

圖3 低溫沖擊過程中藥柱溫度隨時間的變化曲線Fig.3 The temperature variation with the time during low temperature impact test

由圖3可見，在低溫沖擊初期，藥柱表面溫度下降較快，約在600s時表面溫度即降到41℃，藥柱內部溫度梯度很大，表面溫度與中心的溫差達最大值，約為14℃；隨著時間的增加，由于藥柱內部熱量的傳遞，藥柱表面的降溫速率明顯降低，溫度梯度逐漸減少，約7 200s后藥柱中心溫度接近平衡溫度20℃。低溫沖擊600s時藥柱溫度分布見圖4。

由圖4可見，在600s時，藥柱端面邊緣點溫度達到41℃，而藥柱中心的溫度為55℃，溫差達14℃，最大溫度梯度出現在藥柱側面中部，側面中部溫度為50℃。

圖4 低溫沖擊600s時藥柱溫度分布圖Fig.4 Contour of temperature distribution of PBX cylinder in low temperature impact test（t＝600s）

2.2 熱應力分布

采用二維有限元模型計算得到藥柱在高溫沖擊300s時藥柱的熱應力分布，見圖5。

圖5 高溫沖擊300s時藥柱的熱應力分布Fig.5 Contour of major principal stress distribution of PBX cylinder in high temperature impact test（t＝300s）

由圖5可見，藥柱在高溫沖擊過程中表面為壓應力，靠近中心區域為拉應力。最大壓應力出現在藥柱外表面，約為0.79 MPa；最大拉應力出現在藥柱中心部位，約為3.3 MPa，小于相應溫度（21℃）下HMX基PBX的拉伸強度（6.08 MPa），應力水平約為50%，不容易產生裂紋和破壞。

藥柱在低溫沖擊600s時的非穩態熱應力場分布見圖6。

由圖6可見，在低溫沖擊過程中，藥柱表面為拉應力，靠近中心區域為壓應力。最大壓應力出現在藥柱中心部位，約為1.5 MPa；而最大拉應力則出現在藥柱側面中部，約為3.3 MPa。相應溫度（50℃）下HMX基PBX的拉伸強度約為4.18 MPa，應力水平接近80%。研究表明［12］，雖然拉伸應力沒有達到HMX基PBX的破壞強度，但較高水平拉伸應力長期作用仍然會導致HMX基PBX損傷，最終斷裂。因此，在低溫沖擊過程中，可能在藥柱中部產生裂紋，裂紋擴展導致藥柱開裂。

圖6 低溫沖擊600s時藥柱的熱應力分布Fig.6 Contour of major principal stress distribution of PBX cylinder in low temperature impact test（t＝600s）

3 實 驗

3.1 實驗條件及性能測試

為了驗證有限元模擬分析的可靠性，對HMX基PBX藥柱進行溫度沖擊試驗。將HMX基PBX粉末壓制成Φ50 mm×76 mm的藥柱，密度為1.850 g／cm3，HMX基PBX的主要成分為HMX，黏結劑為三元乙丙共聚物。

在20～60℃進行熱沖擊損傷試驗，高溫和低溫轉換時間約為5 min，并在20℃和60℃恒溫2h使樣品溫度均勻。

采用熱電偶測量藥柱表面溫度，用Net DAQ 2640溫度巡檢儀采集數據，分辨率為±0.02℃；準確度為±0.35℃。

采用WS-3811-J24數字式靜態應變采集儀進行應變測量，應變儀的精度為±0.1%，分辨率為10－6，在應變測量中采用溫度補償片補償由電阻片環境溫度引起的附加應變。

采用的聲發射監測系統是一個綜合的、全波形數據采集處理系統，能夠對聲發射信號連續進行實時監測，對整個實驗過程中所產生的AE事件總數進行統計，并對其源參數進行計算，給出實時的聲信號計數和強度值。聲發射儀器有24個通道，采樣頻率為10 MHz，其門檻值設定為100 m V。各測量傳感器布局如圖7所示。

圖7 各測量傳感器布局示意圖Fig.7 The position of sensing element

3.2 實驗結果與模擬結果的比較

溫度沖擊實驗中藥柱表面溫度和應變的測量結果分別見圖8（a）和（b）。

圖8 溫度沖擊試驗中藥柱表面各點溫度及應變隨時間的變化曲線Fig.8 The temperature and strain variation with the time in temperature impact test

由圖8（a）可見，藥柱表面各點的溫度變化趨勢與計算結果相似。高溫沖擊階段，藥柱表面溫度初期變化較快，約20 min后趨于緩慢，約80 min表面溫度趨向于平衡；低溫沖擊階段，約20 min后趨于緩慢，約80 min表面溫度趨向于平衡。由圖8（b）可見，藥柱軸向應變和徑向應變的變化趨勢相似，軸向應變相對大一些，這與藥柱的規格有關。藥柱應變變化趨勢與溫度密切相關，主要是溫度變化過程熱脹冷縮引起。溫度沖擊過程中聲發射強度和計數監測結果見圖9。

圖9 溫度沖擊聲發射信號強度和計數監測結果Fig.9 Acoustic emission signal intension and event count in temperature impact test

由圖9可見，在高溫沖擊3～10 min，出現較強的信號，最高強度63d B，但信號較少，持續時間較短。在低溫沖擊初始階段，又監測到了較強的聲信號，最高達到87d B，在20℃保溫階段，出現大量的聲信號，信號強度略有減小。聲發射信號與炸藥內部出現損傷裂紋有關，說明藥柱的損傷主要發生在低溫沖擊階段，這與前面數值模擬結果一致。

藥柱經溫度沖擊試驗后裂紋及斷口形貌照片見圖10。

圖10 溫度沖擊試驗后藥柱裂紋和斷口照片Fig.10 The crack and fracture of PBX cylinder after temperature impact test

由圖10（a）可見，1號藥柱在側面中部有一條貫穿性裂紋，2號藥柱在相同的位置也出現了局部裂紋，進一步說明溫度沖擊過程中HMX基PBX藥柱破壞方式為拉應力破壞。與普通復合材料不同，HMX基PBX為顆粒黏結黏彈性材料，其拉伸強度遠小于壓縮強度，因此，HMX基PBX的破壞一般都是拉應力破壞所致。對于高溫沖擊過程，由于拉應力出現在藥柱中心位置，應力水平較低，且溫度持續升高使得黏結劑軟化，拉應力也迅速減小，這種塑性變形本身就有閉合裂紋的功能，故裂紋萌生和擴展比較困難［13］。但在低溫沖擊過程中，最大拉應力出現在藥柱外表面，應力水平較高，且隨著溫度的進一步降低，黏結劑塑性下降，容易產生裂紋，當藥柱中心溫度繼續下降時，裂紋向中心擴展導致藥柱開裂。由圖10（b）可見，藥柱斷口外周比較平整，為裂紋萌生部位。

4 結 論

（1）采用二維有限元分析模型，對HMX基PBX在溫度沖擊載荷作用下的溫度場和應力場進行了數值計算，計算時考慮了材料熱力參數隨溫度變化的非線性因素。計算結果表明，在降溫開始階段，藥柱表面溫度下降較快，內部溫度梯度很大，藥柱外表面受拉，中心受壓，側面中部受到較大軸向拉應力作用容易導致裂紋的萌生和擴展。

（2）對HMX基PBX藥柱進行了溫度沖擊試驗，并利用熱電偶、電阻應變片和聲發射技術分別對藥柱表面溫度、應變力及熱損傷進行測量，溫度測量值與數值計算結果趨于一致，聲發射技術信號強度和計數表明裂紋萌生和擴展主要出現在低溫沖擊階段，裂紋位置和斷口特征說明了HMX基PBX熱損傷破壞方式為拉應力破壞。

（3）數值計算與實驗結果表明，HMX基PBX熱損傷破壞方式為拉應力破壞，抗拉強度可以較好地反映材料的耐熱損傷能力。

［1］ 韋興文，周筱雨，涂小珍，等.HMX基PBX的溫度環境適應性［J］.火炸藥學報，2012，35（1）：15-18.

WEI Xing-wen，ZHOU Xiao-yu，TU Xiao-zhen，et al. Thermal environment adaptability of H MX based PBX［J］.Chinese Journal of Explosives and Propellants，2012，35（1）：15-18.

［2］ 馮國富，楊朝卿，于蔭林.高溫暴露試驗對PBX藥柱的影響［J］.火炸藥學報，1997，20（3）：17-19.

FENG Guo-fu，YANG Chao-qing，YU Yin-lin.Influence of high temperature cycling test on PBX columns［J］.Chinese Journal of Explosives and Propellants，1997，20（3）：17-19.

［3］ Willey T M，Lauderbach L，Gagliardi F，et al.Comprehensive characterization of voids and microstructure in TATB-based explosives from 10 nm to 1 cm：effects of temperature cycling and compressive creep［C］∥The 14th International Detonation Symposium.Coeur d＇Alene：ID，2010.

［4］ 顏熹琳，周陽，周筱雨.JOB-9003炸藥的載荷環境試驗［J］.火炸藥學報，2009，32（1）：36-39.

YAN Xi-lin，ZHOU Yang，ZHOU Xiao-yu.Load environmental test of explosive JOB-9003［J］.Chinese Journal of Explosives and Propellants，2009，32（1）：36-39.

［5］ Sandusky H W，Bernerker R R.Influence of fresh damage on the shock reactivity and sensitivity of several energetic materials［C］∥Proceedings of the 10th Symposium on Detonation.Massachusetts：［s.n.］，1993：490-498.

［6］ Meuken D，Scholtes G.Quantification of thermal and mechanical damage in PBX［C］∥Insensitive Munitions and Energetic Materials Technology Symposium. Brislol：［s.n.］，2006.

［7］ 姚惠生，黃風雷，張寶平.炸藥沖擊損傷及損傷炸藥沖擊起爆實驗研究［J］.北京理工大學學報，2007，27（6）：487-490.

YAO Hui-sheng，HUANG Feng-lei，ZHANG Baoping.Experimental study on the impact damage and shock initiation of damaged explosives［J］.Transactions of Beijing Institute of Technology，2007，27（6）：487-490.

［8］ 尹俊婷，袁寶慧，牛鵬俊.炸藥損傷及損傷炸藥環境適應性的實驗研究［J］.火炸藥學報，2008，31（2）：78-80.

YIN Jun-ting，YUAN Bao-hui，NIU Peng-jun.Experimental study on the explosive damage and environment adaptability of damaged explosive［J］.Chinese Journal of Explosives and Propellants，2008，31（2）：78-80.

［9］ 田勇，張偉斌，郝瑩.炸藥熱沖擊損傷破壞及超聲波特性參量檢測［J］.火炸藥學報，2000，23（4）：13-15.

TIAN Yong，ZHANG Wei-bin，HAO Ying.Thermal shock damage of explosives and its ultrasonic characteration［J］.Chinese Journal of Explosives and Propellants，2000，23（4）：13-15.

［10］田勇，羅順火，張偉斌.HMX基PBX“激熱”沖擊損傷破壞及超聲特征［J］.火炸藥學報，2002，25（3）：17-19.

TIAN Yong，LUO Shun-huo，ZHANG Wei-bin.Water-bathed thermal shock damage of PBX JOB-9003 and its ultrasonic characteristics［J］.Chinese Journal of Explosives and Propellants，2002，25（3）：17-19.

［11］韋興文，周筱雨，王培，等.溫度對HMX基PBX炸藥熱膨脹系數和熱導率的影響［J］.火炸藥學報，2012，35（3）：33-37.

WEI Xing-wen，ZHOU Xiao-yu，WANG Pei，et al. Influence of temperature on thermal expansion coefficient and thermal conductivity of HMX based polymer bonded explosive［J］.Chinese Journal of Explosives and Propellants，2012，35（3）：33-37.

［12］周紅萍，何強，李明，等.低拉伸應力下PBX的老化實驗研究［J］.火炸藥學報，2009，32（5）：8-10.

ZHOU Hong-ping，HE Qiang，LI Ming，et al.Experimental study on aging of PBX under low tensile stress［J］.Chinese Journal of Explosives and Propellants，2009，32（5）：8-10.

［13］蘭瓊，戴斌，楊白鳳，等.溫壓時效處理PBX內部裂紋愈合現象研究［J］.含能材料，2013，21（2）：205-208.

LAN Qiong，DAI Bin，YANG Bai-Feng，et al.Healing of cracks in PBX by thermal pressure aging treatment［J］.Energetic Material，2013，21（2）：205-208.

Numerical Calculation and Experimental Study on Thermal Damage of HMX Based Polymer Bonded Explosive

WEI Xing-wen，WU Shu-li，TANG Xing
（Institute of Chemical Materials，CAEP，Mianyang Sichuan 621900，China）

The temperature field and the thermal stress field of HMX based polymer bonded explosive（PBX）during thermal impact process was simulated by the finite element software ANSYS based on a two-dimentional axis symmetric model. The experiment of thermal shock was carried out by using thermocouples and acoustic emission technique to measure temperature and thermal damage respectively.The calculated results show that thermal stress leads to the thermal damage in cooling stage，and the largest tensile stress takes place in the lateral centre of the cylindrical sample.The experiment results agree well with calculated results，indicating that thermal damage of HMX based PBX is caused by tensile stress，and the tensile strength can reflect the ability of heat shock resistance for HMX based PBX preferably.

material science；HMX；polymer bonded explosive；PBX；thermal damage；numerical simulation；finite element analysis

TJ55；X93

A

1007-7812（2014）04-0009-05

2014-01-22；

2014-06-15

韋興文（1977－），男，副研究員，碩士，從事含能材料老化研究工作。