CHINESE JOURNAL OF ENERGETIC MATERIALS
+高级检索
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于X射线显微CT的PBX热冲击损伤特征  PDF

  • 杨存丰
  • 田勇
  • 张伟斌
  • 张催
  • 肖盼
中国工程物理研究院化工材料研究所, 四川 绵阳 621900

中图分类号: TJ55

最近更新:2022-08-24

DOI:10.11943/CJEM2021247

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

为研究小尺寸高聚物黏结炸药(polymer bonded explosive,PBX)半球的热冲击损伤特征,针对Φ10 mm TATB基和HNS基PBX半球样品开展了0~100 ℃水浴温度冲击实验,利用X射线显微层析成像(Micro‑computed tomography,μCT)技术研究了样品损伤的三维形态和分布特征,并采用热弹塑性二维轴对称模型对样品温度冲击过程的热传导和热应力进行了模拟分析。CT结果表明温度冲击后2种样品均由边角起裂,其中TATB基PBX半球沿边部环向扩展,形态曲折,具有撕裂和脆断特征;HNS基PBX半球沿轴向基本贯穿,形态平直,具有脆断特征。数字模拟结果表明温度冲击过程中TATB基PBX半球内部产生较强的拉应力,由边角区域至中心区域拉应力先后超过其拉伸强度,导致主裂纹由半球边角区域萌生并向内扩展。样品损伤形态与温度冲击作用下的应力分布特征、黏结剂温度特性吻合。本研究为TATB基和HNS基PBX温度冲击损伤机理分析奠定了基础。

图文摘要

In order to study the damage characteristics of small scale polymer bonded explosive(PBX) hemisphere under thermal shock, water‑bathed thermal shock test from 0 ℃ to 100 ℃ was carried out for Φ10 mm TATB based and HNS based PBX hemispheres respectively. The thermal elastic‑plastic response of PBX hemisphere was calculated by 2D axisymmetric finite element model which contained temperature related material properties, and 3D morphology and distribution characteristics of inner cracks for PBX hemisphere were studied by X‑ray micro‑computed tomography. The correlation between the damage characteristics and elastic‑plastic response for PBX hemisphere was discussed.

0 引 言

高聚物黏接炸药(polymer bonded explosive,PBX)是由单质炸药晶体和高聚物黏结剂等组成的混合炸药,具有颗粒高度填充、非均相的结构特征。PBX广泛应用于各类武器,具有结构承载、爆轰破坏等重要功能。PBX制造过程中会产生微裂纹、微空隙等初始损

1-3,同时在复杂的服役环境条件下,PBX承受力、热等载荷作用,初始损伤进一步扩展并产生新的损伤。损伤使PBX的力学性能劣化、安全性能和爆轰性能降4-5,影响武器可靠性。

为掌握PBX在温度载荷作用下的损伤机制,研究人员综合运用数值模

6-7以及声发6-7、超8-10、光学显微镜和扫描电11-13、工业CT11-14等技术手段,研究了PBX的热损伤行为,对热损伤机理研究有了一定的认识。其中声发射和超声方法通过PBX损伤的声发射和超声特征信号表征损伤,难以直观地表征损伤形态特征。光学显微镜和扫描电镜方法主要表征PBX的表面损伤,观测视野小且取样具有破坏性。X射线显微CT具有微米尺度分辨、三维层析成像的技术优势,可以无损地透视样品内部结构并精确地定量表征其结构与损伤的三维形态特征,为PBX损伤研究提供了重要的技术支15

热应力是导致PBX损伤的最常见的和最主要的形式,其中温度冲击最容易使PBX产生损

12,由此抗热冲击性能在PBX的设计和应用中受到广泛关注。但目前相关研究主要采用PBX柱形结构作为研究对象,而对PBX炸药半球形结构研究较少,PBX柱形结构和球形结构的温度冲击响应存在明显差异,同时PBX构型尺寸影响其温度冲击过程的热应力分布。为掌握小尺寸PBX半球形结构的温度冲击损伤特征,本研究针对Φ10 mm的TATB(1,3,5‑三氨基‑2,4,6‑三硝基苯)基和HNS(六硝基芪)基PBX半球样品,采用水浴对流热交换法对样品进行温度冲击实验以激发样品产生损伤,利用X射线显微CT和三维图像处理技术表征了2种PBX半球样品的温度冲击损伤的形态特征,采用热弹塑性二维轴对称模型对TATB基PBX样品温度冲击过程的热传导和热应力进行了模拟分析,探讨了温度冲击损伤特征与温度分布、热应力分布以及黏结剂温度特性的关联性,为PBX半球的温度冲击损伤规律研究提供了指导。

1 实验部分

1.1 实验样品

TATB、HNS,化工材料研究所;黏结剂F2311、F2314为工业级,晨光化工研究院。

分别采用TATB和F2311、HNS和F2314材料组合,经水悬浮造粒后压制成TATB基和HNS基PBX半球形样品,数量各4发,编号分别为A1~A4、B1~B4。样品直径为Φ10 mm,样品密度约1.60 g·cm-3。F2311、F2314的玻璃化转变温度分别为0 ℃、49.6 ℃,软化点温度分别约80 ℃、120 ℃。

1.2 仪器设备

X射线显微CT,主要由微焦点X射线源、多轴运动控制精密机械系统、光耦探测器组成,射线源最大管电压150 kV,探测器总像素阵列2048×2048,像素尺寸67.5 μm。

1.3 实验

1.3.1 温度冲击实验

采用水浴对流热交换方法开展温度冲击实验,利用高热流密度冲击激发样品产生损伤,实验方法为:在室温(20 ℃)条件下,样品经CT检测后,依次浸入0 ℃冰水混合物、100 ℃热水中各保持30 min,2种温度间的转换时间控制在2 s,样品由100 ℃热水取出30 min后再次进行CT检测。

1.3.2 显微CT成像与损伤表征

采用X射线显微CT对样品进行三维CT成像,X射线显微CT检测示意图见图1,检测参数为:射线源的管电压50 kV、功率10 W,单幅投影图采样曝光时间2.5 s,样品旋转360°,投影图采集幅数4501,采集的图像像素阵列1024×1024,重建后的样品三维CT图像像素尺寸为11.9 μm。

图1  X射线显微CT检测示意图

Fig.1  Schematic diagram of X‑ray micro‑CT

所得样品CT图像为灰度值fxy)的空间分布图像,fxy)代表样品位置(xy)上材料的X射线衰减系数,PBX基体和裂纹因X射线衰减系数不同导致影像灰度存在差异。基于图像的不连续性和相似性基本原则,利用sobel边缘检测改进的灰度直方图阈值分割方

16,提取PBX内部裂纹的三维形貌,观察裂纹形态特征,评估损伤的产生与扩展方向,测量裂纹的体积和表面积评估PBX损伤程度,测量垂直于PBX半球中轴线的不同位置截面的损伤面积比分析损伤与位置的关系。

1.4 数字模拟

PBX半球温度冲击过程中,结构及载荷具有轴对称性,因此将三维空间问题简化为二维轴对称平面问

6-7,采用Abaqus有限元分析软17进行热传导和热应力分析。几何模型选取SR 5 mm半球的1/2中轴面建立轴对称模型,按四结点轴对称热耦合四边形单元CAX4T进行网格划分。TATB基PBX样品的力学性能和热性能与温度相关,其数字模拟采用的力学性能参数引用文献[18],热物理性能参数引用实测数据,其中线膨胀系数在温度243~303 K时为4.18×10-5 K-1、303~343 K时为5.81×10-5 K-1,比热为1.08×103 J∙kg-1K-1,导热系数为0.247 W∙m-1K-1

2 结果与讨论

2.1 损伤的三维形态特征

温度冲击实验前后对所有样品A1~A4、B1~B4均进行了CT检测,由于同类样品的损伤形态特征相同,研究选用样品A1和B1作损伤特征分析。温度冲击实验前,样品A1的三维CT图像、样品B1的内部裂纹三维形貌分别见图2a和2b,由图2a和2b可见,样品A1内部未见裂纹,样品B1距顶部约1.3 mm处存在环状初始裂纹。由图2c和2d温度冲击实验后样品A1和样品B1的三维CT图像可见,2种样品均产生明显开裂,开裂宽度由样品边角向内逐渐减小,呈现由边角部位起裂并向内扩展的特征;2种样品的主裂纹扩展形态明显不同,样品A1主裂纹沿其边部环向扩展,样品B1主裂纹沿轴向平直扩展。

图2  样品A1与B1温度冲击前后形貌(a,b. 温度冲击前样品A1,B1的三维CT图像; c,d. 温度冲击后样品A1,B1的三维CT图像; e,f. 温度冲击后样品A1,B1内部裂纹三维形貌)

Fig.2  morphology of sample A1 and B1 before and after temperature shock(a,b. 3D CT image for sample A1 and B1 before temperature shock; c,d. 3D CT image for sample A1 and B1 after temperature shock; e,f. 3D morphology of inner cracks for sample A1 and B1 after temperature shock

2e和2f分别为温度冲击后样品A1和B1的内部裂纹形貌,裂纹形貌图反映了裂纹宽度和边缘梯度的变化信息,垂直于梯度变化的方向为裂纹的扩展方向(如图2e和2f箭头所示)。由图2e可见,样品A1边角部位产生多条裂纹并沿边部扩展,导致裂纹三维形貌整体呈环状特征,且多数裂纹的梯度纹理形态不规则、方向不一致,具有撕裂特征。由图2f可见,样品B1内部新生了严重裂纹,其中主裂纹(如图2f中位置1所示)沿边部起裂,轴向扩展至微裂纹(图2f中位置4所示)断裂面,形态平直,具有脆断特征;次裂纹(图2f中位置2、3、5所示)轴向扩展至主裂纹面。样品B1内部的初始裂纹沿球径、环向及轴向3个方向轻微扩展。

基于图像分割处理得到的2种样品内部裂纹三维形貌特征,由裂纹所有体素的体积和计算裂纹体积,由裂纹与背景相邻的像素面积和计算裂纹表面积,分析得到温度冲击前后样品内部裂纹的体积及表面积,结果见表1。由表1可知,样品的损伤程度不相同,其中样品A1~A4损伤的体积比为0.210%~0.45%,表面积为226~601 mm2;样品B1~B4损伤的体积比为0.150%~0.34%,表面积为264~403 mm2。2种样品相比较,TATB基PBX样品损伤的体积比分布和表面积分布大,表明其损伤扩展程度大,其原因是裂纹环向扩展导致长度大、裂纹纹理不规则导致损伤体积和表面积增大。

表1  裂纹的体积和表面积
Table 1  Volume and surface area of cracks

sample

number

before temperature shockafter temperature shock
Volume/mm3crack volume ratio/%surface/mm2volume/mm3crack volume ratio%surface/mm2
A1 - - - 7.85 0.450 345.47
A2 - - - 6.82 0.390 601.73
A3 - - - 4.21 0.240 248.24
A4 - - - 3.70 0.210 226.71
B1 0.08 0.005 17.69 5.92 0.340 403.61
B2 0.11 0.010 14.16 5.69 0.320 272.49
B3 0.04 0.002 5.06 5.31 0.300 236.26
B4 0.14 0.010 12.39 2.71 0.150 264.58

基于CT检测获得的裂纹特征,可以看出在0~100 ℃温度冲击作用下,样品A1和B1的主裂纹由边角部位起裂并向内部扩展,表明温度冲击过程中样品边角部位的热应力先超过了其拉伸强度。样品A1和B1分别采用了F2311、F2314黏结剂,其温度特性存在明显差别,其中F2311、F2314的玻璃化转变温度分别为0 ℃、49.6 ℃,软化点温度分别约为80 ℃、120 ℃。在0~100 ℃温度冲击过程中,F2311由高弹态向黏流态转变,样品A1开裂特征由脆断向撕裂转变;同时随着温度升高,样品内部的黏结剂F2311逐步软化,炸药与黏结剂之间的界面作用相应减弱,造成结构强度降低,在半球内部周向热应力作用下产生环状开裂;F2314由玻璃态向高弹态转变,样品B1仍保持一定的脆性和结构强度,在周向热应力作用下产生轴向脆性开裂。

2.2 损伤的分布特征

样品A1和B1端面及距端面1785 μm、3570 μm位置的径向CT切片见图3。由图3a~3c可见,样品A1的3个切片位置上裂纹形态复杂,呈现撕裂和脆断2种特征。样品端面截面(图3a)上裂纹1和2位于半球边角区域,边角处开裂宽度大,分别约144 μm和101 μm,显示裂纹由边部向内部扩展;裂纹3和4由数量较多的微裂纹构成,形态不规则,具有撕裂特征。距端面1785 μm截面(图3b)上,裂纹1、3、4基本融合,裂纹2、裂纹1及其相连的裂纹3处开裂宽度大,同样显示出裂纹由边部向内部扩展的特征。距端面3570 μm截面(图3c)上,裂纹3和4融合,形态较平直,具有一定的脆断特征。由图3b、3c可见,裂纹5产生于样品中部并偏离中轴方向扩展。综合以上裂纹特征,可以看出:主裂纹由样品边角部位产生并向内扩展;样品的开裂特征主要与黏结剂的温度特性相关,端面处温度高导致黏结剂迅速软化,由此出现明显的撕裂特征,随着热由样品表面向内部传递,样品内部温度逐渐升高,但在样品开裂时样品中部区域黏结剂软化尚不明显,由此开裂呈现一定的脆断特征。

图3  温度冲击后样品A1和B1径向CT切片

Fig.3  Radial CT slices of sample A1 and B1 after temperature shock

图3d~f可见,样品B1的3个切片位置上主裂纹(裂纹1)形态平直、轴向基本贯穿,开裂宽度边部约为121 μm、向内逐渐减小。端面位置(图3d)主裂纹扩展于裂纹4处停止,次裂纹(裂纹2、5)形态较平直,但存在部分不连续。距端面1785 μm位置(图3e)次裂纹(裂纹2、5)形态平直,其扩展于主裂纹处截止。距端面3570 μm位置(图3f),次裂纹(裂纹3、5)形态平直,初始裂纹(裂纹6)出现环向和轴向扩展特征。总体上样品B1内部裂纹形态整体上较为规则,呈现明显的脆断特征;主裂纹由边部向内部扩展。

对样品所有径向切片,由样品端面至顶部,依次测量其损伤面积比(裂纹面积与总面积之比),得到样品A1~A4和样品B1~B4损伤面积比随切片位置的分布曲线见图4,由图4可见,样品A2~A4在0~300切片位置的损伤面积比缓慢下降、样品B1~B4的损伤面积比在0~320切片位置的损伤面积缓慢增大,产生该分布规律的主要原因是样品A2~A4开裂沿边部扩展导致顶部区域损伤低、样品B1~B4轴向扩展导致各截面损伤特征基本相同。

a.  TATB based PBX samples

b.  HNS based PBX samples

图4 不同样品的径向截面损伤面积比

Fig.4 Curves of damage area ratio at different radial slices for different samples

2.3 温度冲击过程的热传导和热应力有限元分析

由于TATB基和HNS基PBX样品均呈现由边角部位起裂向内扩展的形态特征,本研究选用TATB基PBX样品作热传导和热应力分析。TATB基PBX样品受0~100 ℃温度冲击1.3 s时的温度分布和应力分布见图5,由图5可见,半球内部为拉应力,其最大值为7.49 MPa,拉应力超过其拉伸强度,使PBX内部产生开裂;半球外表面为压应力,端面边部压应力最大,最大值为3.2 MPa,远低于其压缩强度,压应力不会使PBX产生开裂。

a.  Temperature distribution

b.  Stress distribution

图5 温度冲击1.3 s时的温度分布和应力分布

Fig.5 Temperature distribution and stress distribution at high temperature shock 1.3 s

为研究TATB基PBX样品内部温度和应力分布规律,在应力最大的方向,由半球中轴至边角部位依次选取4点N1~N4(图5b),采用有限元方法分析得到N1、N2、N3、N4点的温度和应力历史,结果如图6所示。由图6a可见,温度冲击过程中,由边角至中部(N4、N3、N2、N1点)温升梯度依次减小;由图6b可见,温度冲击数秒内,N4、N3、N2、N1各点的拉应力迅速增大,拉应力升高梯度基本相同、应力峰值依次增大,应力峰值分别为7.46、11.43、12.90、13.63 MPa。结合不同温度下的拉伸强度数

18,样品由边角至中部的拉应力先后超过其拉伸强度,这与样品沿边角处起裂向向内扩展的形态吻合。

a.  Temperature history

b.  Stress history

图6 N1~N4点的温度和应力历史

Fig.6 Temperature and stress history of point N1-N4

另外,TATB基PBX样品由0 ℃转至100 ℃环境后,由表及里温度逐渐升高,TATB基PBX样品中的黏结剂F2311软化,导致其脆性降低、韧性增强,在拉应力作用下样品开裂由脆断向撕裂转变,尤其半球近表面撕裂特征更明显。根据TATB基PBX半球样品温度冲击后的损伤形态CT检测结果,主裂纹由边角区域起裂后,在拉应力的作用下,裂纹偏向于沿拉应力超过其拉伸强度、相对的软化区和非软化区中间的强度薄弱部位扩展。HNS基PBX半球中黏结剂F2314在半球内部起裂时处于玻璃态,随着最大拉应力由半球边角向中心演变,裂纹迅速向内扩展,形成了垂直于端面的、形态平直的主裂纹,同时产生部分形态平直的次裂纹等,裂纹具有脆断特征。

3 结 论

(1)0~100 ℃高温冲击下,PBX半球内部拉应力由边角到中心区域先后超过其拉伸强度,导致半球内部产生由边角区域起裂并向内扩展的主裂纹。

(2)由0 ℃转至100 ℃高温环境,PBX半球内部温度由表及里温度依次升高,TATB基PBX半球中的F2311黏结剂逐步软化,导致其开裂呈脆断和撕裂特征,裂纹偏向于沿半球边部环状扩展。HNS‑PBX半球中的F2314黏结剂处于玻璃态,导致其开裂呈脆断特征,裂纹偏向于沿半球轴向扩展。

(3)HNS‑PBX半球顶部的初始裂纹在高温冲击作用下沿径向及环向、轴向扩展,扩展程度远低于半球内部新生的主裂纹和次裂纹。

参考文献

1

梁华琼雍炼唐常良.RDX为基的PBX炸药压制过程损伤形成研究[J]. 含能材料2009176):713-716. [百度学术] 

LIANG Hua‑qiongYONG LianTANG Chang‑liang. Pressing damage of RDX‑based polymer bonded explosive[J]. Chinese Journal of Energetic Materials (Hanneng Cailiao)2009176): 713-716. [百度学术] 

2

梁华琼雍炼唐常良.压制过程中PBX炸药颗粒的破碎及损伤[J]. 火炸药学报2010331): 27-30. [百度学术] 

LIANG Hua‑qiongYONG LianTANG Chang‑lianget al. Crack and damage of PBX during pressing[J]. China Journal of Explosive and Propellants2010331): 27-30. [百度学术] 

3

张伟斌李敬明杨雪海.TATB颗粒温压成形PBX的初始细观损伤[J].含能材料2015232):202-204 [百度学术] 

ZHANG Wei‑binLI Jing‑mingYANG Xue‑haiet al. Initial mesoscopic damage of TATB based PBX pressed by warm compaction[J]. Chinese Journal of Energetic Materials (Hanneng Cailiao)2015232):202-204. [百度学术] 

4

陈鹏万黄风雷.含能材料损伤理论及应用[M]. 北京北京理工大学出版社2006. [百度学术] 

CHEN Peng‑wanHUANG Feng‑Iei. Damage theory and applications of energetic materials[M]. BeijingBeijing Institute of Technology Press2006. [百度学术] 

5

唐明峰温茂萍涂晓珍.高温及机械应力对PBX力学行为的影响规律及机理分析[J].含能材料2018262):150-155 [百度学术] 

TANG Ming‑fengWEN Mao‑pingTU Xiao‑zhenet al. Influence and mechanism of high temperature and mechanical stress on the mechanical behaviors of PBXs[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao)2018262): 150-155. [百度学术] 

6

韦兴文吴束力唐兴. HMX基PBX炸药热损伤的数值计算与实验研究[J].火炸药学报2014374):9-13. [百度学术] 

WEI Xing‑wenWU Shu‑liTANG Xing. Numerical calculation and experimental study on thermal damage of HMX based polymer bonded explosive[J]. Chinese Journal of Explosives & Propellants2014374):9-13. [百度学术] 

7

唐明峰甘海啸温茂萍. 含缺口PBX药柱热冲击响应的数值模拟及试验[J]. 含能材料2021291):41-47. [百度学术] 

TANG Ming‑fengGAN Hai‑xiaoWEN Mao‑pinget al. Simulation and experimental study on the thermal shock behavior of notched PBX cylinders[J].Chinese Journal of Energetic Materials(Hanneng Cailiao)2021291):41-47. [百度学术] 

8

田勇张伟斌郝莹. 炸药热冲击损伤破坏及超声波特性参量检测[J].火炸药学报2000234): 13-15. [百度学术] 

TIAN YongZHANG Wei‑binHAO Yinget al.Thermal shock damage of explosive and its ultrasonic characterization[J]. Chinese Journal of Explosives and Propellants2000234):13-15. [百度学术] 

9

田勇罗顺火张伟斌. JOB‑9003“激热”冲击损伤破坏及超声特征[J].火炸药学报2002253): 17-19. [百度学术] 

TIAN YongLUO Shun‑huoZHANG Wei‑binet al. Water-bathed thermal shock damage of PBX JOB⁃9003 and its ultrasonic characteristics[J]. Chinese Journal of Explosives & Propellants2002253):17-19. [百度学术] 

10

田勇张伟斌温茂萍. JOB‑9003高聚物粘结炸药热冲击损伤破坏相关性研究[J]. 含能材料2004123):174-177. [百度学术] 

TIAN YongZHANG Wei‑binWEN Mao‑pinget al. Research on correlation of thermal shock damage of PBX JOB‑9003[J]. Chinese Journal of Energetic Materials (Hanneng Cailiao)2004123): 174-177. [百度学术] 

11

邵珠格刘如沁吴艳青. 高温作用后HMX基PBX热损伤表征试验研究[J].火炸药学报2020434): 406-412 [百度学术] 

SHAO Zhu‑geLIU Ru‑qinWU Yan‑qinget al. Experimental study on thermal damage characterization of HMX based PBX after high‑temperature treatment[J]. Chinese Journal of Explosives & Propellants2020434):406-412 [百度学术] 

12

尹俊婷袁宝慧牛鹏俊.炸药损伤及损伤炸药环境适应性的实验研究[J].火炸药学报2008312):78-80. [百度学术] 

YIN Jun‑tingYUAN Bao‑huiNIU Peng‑junet al. Experimental study on the explosive damage and environment adaptability of damaged explosive[J]. Chinese Journal of Explosives & Propellants2008312): 78-80. [百度学术] 

13

张冬梅韩芳贾林.温度循环载荷作用下压装A‑IX‑II装药的裂纹机理研究[J].含能材料2016391):89-93 [百度学术] 

ZHANG Dong‑meiHAN FangJIA Linet al. Study on crack formation mechanism of pressed charge A‑IX‑II under temperature cyclic load[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao)2016391):89-93 [百度学术] 

14

刘瑞鹏王红星王浩.高低温环境对不同炸药结构影响研究[J].含能材料20051320):20-22 [百度学术] 

LIU Rui‑pengWANG Hong‑xingWANG Haoet al. Effect of environment temperature on different explosive structures[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao)20051320):2 0-22 [百度学术] 

15

刘本德陈华张伟斌. 基于CT图像序列的围压TATB基PBX冲击损伤特性[J]. 含能材料2019273):196-201. [百度学术] 

LIU Ben‑deCHEN HuaZHANG Wei‑binet al. Impact damage characteristics of TATB based polymer bonded explosive under confining pressure based on the CT image sequences[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao)2019273):196-201. [百度学术] 

16

RAFAEL C GonzalezRICHARD E Woods. 数字图像处理[M].北京电子工业出版社2017443-492. [百度学术] 

Rafael C GonzalezRichard E Woods. Digital Image Processing[M]. BeijingPublishing House of Electronics Industry2017443-492. [百度学术] 

17

岳小亮.温度冲击载荷下药柱的力学响应研究[D].南京南京理工大学201318-24. [百度学术] 

YUE Xiao‑liang. Study on the mechanical response of grain under temperature impact load[D].NanjingNanjing University of Science and Technology201318-24. [百度学术] 

18

涂小珍张波韦兴文.TATB基高聚物粘结炸药高温力学性能[J]. 含能材料2013213):306-309. [百度学术] 

TU Xiao‑zhenZHANG BoWEI Xing‑wenet al. Mechnical properties of TATB based PBX at high temperature[J].Chinese Journal of Energetic Materials(Hanneng Cailiao)2013213):306-309. [百度学术]