材料工程  2017, Vol. 45 Issue (4): 77-83   PDF    
http://dx.doi.org/10.11868/j.issn.1001-4381.2016.001212
0

文章信息

刘晓俊, 任会兰, 宁建国
LIU Xiao-jun, REN Hui-lan, NING Jian-guo
不同配比W/Zr活性材料冲击反应实验研究
Experimental Study on Impact Response of W/Zr Reactive Materials with Different Proportions
材料工程, 2017, 45(4): 77-83
Journal of Materials Engineering, 2017, 45(4): 77-83.
http://dx.doi.org/10.11868/j.issn.1001-4381.2016.001212

文章历史

收稿日期: 2016-10-13
修订日期: 2016-12-15
不同配比W/Zr活性材料冲击反应实验研究
刘晓俊, 任会兰 , 宁建国    
北京理工大学 爆炸科学与技术国家重点实验室,北京 100081
摘要: 以W,Zr和ZrH2为原料,采用热压烧结工艺制备不同配比的W/Zr活性材料。通过分离式霍普金森压杆装置 (SHPB) 对W/Zr材料进行动态压缩测试并记录其冲击引发反应过程,同时采用X射线衍射仪 (XRD)、扫描电镜 (SEM)、能谱仪 (EDS) 等对试件反应前后的物相、成分、微观形貌进行测试表征。结果表明:热压烧结制备的W/Zr活性材料结构较为致密,其相对密度≥87.5%;W/Zr活性材料准静态压缩强度≥1022 MPa,破坏应变≤1%,为典型高强脆性材料;3组配比材料动态压缩强度均随应变率增加而提高,当冲击加载达到一定强度时,W/Zr材料破碎并与空气发生剧烈燃烧反应并生成ZrO2
关键词: W/Zr    活性材料    热压    冲击反应   
Experimental Study on Impact Response of W/Zr Reactive Materials with Different Proportions
LIU Xiao-jun, REN Hui-lan , NING Jian-guo    
State Key laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
Abstract: Using W, Zr and ZrH2 as raw materials, W/Zr reactive materials with different proportions were prepared by hot-pressing method. The dynamic behavior of W/Zr materials was investigated by a split Hopkinson pressure bar apparatus (SHPB) and the impact-initiated reaction process was recorded. The characteristics of the samples, such as the phase transformation, component and microstructure before and after reaction, were analyzed by X-ray diffraction (XRD), energy dispersive spectroscopy (EDS) and scanning electronic microscopy (SEM). The results show that the W/Zr material sintered by hot-pressing is more compact and the relative density exceeds 87.5%. W/Zr material is a typical brittle material with high strength, the quasi-static compressive strength is higher than 1022 MPa and the failure strain is less than 1%. The effect of strain rate on the dynamic compressive strength for all three group samples is obvious, strength increases with strain rate. When the impact load is strong enough, W/Zr material crushes, combusts with air violently, and generates ZrO2.
Key words: W/Zr    reactive material    hot-pressing    impact-induced reaction   

活性材料 (reactive materials) 最早由Montgomery[1]在其专利US3961576中提出,是一类实现战斗部高效毁伤的亚稳态含能材料。活性材料通常由两种或更多种非爆炸性固体材料组成,正常处置和存储时都保持惰性,而当受到强冲击加载时会发生化学反应并释放能量,产生燃烧、爆轰等效果,但与常规弹药不同的是活性材料并不能自持爆轰[2, 3]。利用活性材料制备含能毁伤元—活性破片,其结合动能侵彻和内爆效应对目标造成致命打击,能够显著提高战斗部的毁伤效能[4]。美国海军水面作战中心 (Naval Surface Warfare Center,NSWC) 于2002年11月进行的地面静爆威力演示实验表明,活性破片战斗部对典型空中目标摧毁级杀伤半径是惰性破片战斗部的2倍,毁伤威力为惰性破片的5倍,化学潜能更是惰性破片平均动能的12倍[5, 6]。此外,活性破片在侵彻目标时会迅速发生反应而非继续高速飞散,这也能够降低对目标周围环境的附带损害[7]。活性材料因其独特的特性和在军事领域的广阔应用前景而受到越来越多的关注。

活性材料包括铝热剂、金属/聚合物、金属化合物、亚稳态分子化合物和氢化物等[8-10]。目前,研究较多的是金属/聚合物类型活性材料,其具有含能高、加工简单等特性,但Nielson在其专利US0020397A1中发现该类型活性材料的密度和强度普遍偏低,这使其在应用时受到一定的限制[11]。相对而言,金属化合物类型活性材料的密度和强度则普遍要高很多,钨锆合金就是一个典型代表,除含能外还可作为结构材料。王树山、陈伟等对钨锆合金破片和93钨合金的毁伤过程进行了研究,实验结果表明钨锆合金引爆战斗部装药和引燃油箱的能力更强[12, 13]。Zhang等采用弹道枪驱动钨锆合金破片对金属腔体进行侵彻实验研究,研究发现撞击速率对能量释放率有较大影响[14]。Luo等采用同样测试手段对钨锆基复合材料破片进行研究,实验结果表明靶板厚度也会影响反应进程,而且发现钨锆基玻璃合金比钨锆合金含能更高[15]

本工作采用热压烧结工艺制备钨锆材料,测试并分析材料的显微组织形貌、力学性能和反应活性,为其作为一种高密度活性材料的选用和加工提供理论依据。

1 实验材料与方法 1.1 实验材料

W粉,平均粒径3 μm,纯度99.9%;Zr粉,平均粒径40 μm,纯度99%;ZrH2粉,平均粒径3 μm,纯度99%;分析纯乙醇。

1.2 试件制备

本研究设计制备3种不同配比的W/Zr材料,各组分配比如表 1所示。试件的制备步骤主要有湿混、烧结和切割,具体制备流程如图 1所示。实验所用烧结设备为真空气氛压力炉。

表 1 W/Zr试件成分及配比 (质量分数/%) Table 1 Compositions and proportions for W/Zr samples (mass fraction/%)
Sample Zr ZrH2 W
WZ-1 66 0 34
WZ-2 0 66 34
WZ-3 43 0 57
图 1 W/Zr试件制备过程示意图 Fig. 1 Schematic representation of preparation of process for W/Zr samples
1.3 实验方法

准静态和动态力学测试分别是在北京理工大学爆炸科学与技术国家重点实验室的WDW-300型材料试验机和Ф14.5 mm和Ф5 mm分离式霍普金森压杆实验平台上进行,并采用Photron FASTCAM SA5型高速摄影机观测整个动态撞击过程。SHPB装置示意图如图 2所示,本实验在入射杆撞击端加装黄铜垫片进行波形整形,以实现恒应变率加载,并使用超动态信号测试分析系统进行应变数据采集。此外,在试件两端均加装高强钢垫块以防止杆端面受损,并利用塑料外壳罩住整个试件回收反应后的碎片。采用S-4800型场发射扫描电子显微镜进行SEM表征;使用D8 advance型X射线衍射仪进行XRD测试。

图 2 分离式霍普金森压杆装置示意图 Fig. 2 Schematic of a split Hopkinson pressure bar apparatus
2 结果与分析 2.1 试件形貌表征及物相分析

图 3为热压烧结工艺制备的W/Zr试件照片。3种不同配比的烧结块体表面均未呈现金属光泽,如图 3(a)所示,主要是由于表面包覆着一层用于脱模的石墨纸。经线切割加工成尺寸为Ф5 mm×5 mm的试件,如图 3(b)所示,此时明显发现WZ-1~3的色泽依次变暗。

图 3 W/Zr试件图片 (a) 切割前;(b) 切割后 Fig. 3 Photographs of W/Zr samples (a) before cutting; (b) after cutting

在W-Zr二元系相图中显示,1400~1735 ℃下在富Zr区存在化合物W2Zr和β (Zr)[16, 17]图 4是3组W/Zr材料XRD分析结果,它们各自的组成成分为WZ-1:W2Zr,ZrC,Zr;WZ-2:ZrC0.32H1.2,W,Zr;WZ-3:W2Zr,ZrC,W。3组材料均检测到Zr的碳化物,这表明在烧结过程中,用于脱模的石墨纸部分与Zr发生了反应。郭星渠[18]等在用ZrC细化晶粒制取ZrH2块体材料时发现,弥散于Zr基体中的ZrC经900 ℃氢化后会转化为ZrC0.32H1.2。ZrH2于700 ℃温度下脱氢,即可制得Zr粉。WZ-2试件中存在ZrC0.32H1.2表明,受Ar气氛的正压环境影响,ZrH2脱氢后产生的H2没能完全释放到外部而再次参与了反应。WZ-2产物中除了ZrC0.32H1.2之外还含有少量的单质Zr,这是由于石墨纸只有少部分渗入材料内部,C不足未能使Zr完全继续转化为ZrC。

图 4 X射线衍射分析 Fig. 4 Analysis of X-ray diffraction pattern

上述WZ-1~3材料生成机理可由以下方程表示:

(1)
(2)
(3)
(4)

对3组W/Zr材料进行电镜扫描和EDS分析,WZ-1~3试件的微观形貌照片如图 5所示。图 5(a)显示WZ-1试件内部的结构十分致密,经测定其密度为8.34 g/cm3,接近理论密度 (相对密度99.2%)。图 5(a)中呈现3种不同灰度值的区域 (白-浅灰-深灰),图 6给出了A点附近浅灰色区域EDS分析结果,检测出Zr和C元素的质量分数分别为83.7%和16.3%,结合图 4的XRD分析结果可以判定所测区域物质为ZrC。同时测定白色相为化合物W2Zr,少数深灰色区域物质为单质Zr。WZ-2试件的密度可达其理论密度的95.3%。图 5(b)中呈现白色的W颗粒 (直径1~5 μm),均匀地分布在ZrC0.32H1.2基体上。图 5(c)显示WZ-3试件内部结构松散,存在大量的孔洞,其相对密度仅为87.5%。这些与WZ-1试件差异较大,可能是由于WZ-3试件W含量更高,烧结温度不够所致。

图 5 W/Zr试件SEM形貌 (a) WZ-1;(b) WZ-2;(c) WZ-3 Fig. 5 SEM micrographs of W/Zr samples (a) WZ-1; (b) WZ-2; (c) WZ-3
图 6 图 5(a)中A点周围区域灰色相EDS分析结果 Fig. 6 EDS analysis results of gray phase around point A in fig. 5(a)
2.2 试件力学性能分析

对W/Zr材料进行准静态压缩测试,WZ-1~3试件强度分别为1880,1200 MPa和1022 MPa,破坏应变分别为1.0%,0.7%和0.4%,这表明W/Zr材料是一种典型的高强度脆性材料。对WZ-1~3试件进行恒应变率SHPB测试,多次测量取平均后获得的动态压缩强度如图 7所示。

图 7 不同应变率下W/Zr试件压缩强度 Fig. 7 Compressive strength of W/Zr samples under different strain rates

图 7可见,3种材料的压缩强度都表现出一定的正应变率效应。WZ-1试件动态压缩强度明显高于另外两组,应变率为1000 s-1时其值为2690 MPa。WZ-2试件应变率在230~1200 s-1,其压缩强度变化不大 (2186~2340 MPa),仅提高7%。WZ-3试件在应变率为470 s-1时的压缩强度为1060 MPa,与准静态压缩强度1022 MPa差别不大,而当加载应变率为690~ 1120 s-1时,其动态压缩强度则提高至1410~1581 MPa。这主要是由于WZ-3试件内部存在许多微孔洞,这些孔洞与周围脆性材料组成了大量微折板结构,而该折板结构存在一个显著的应力下降现象,在高速撞击时 (高应变率) 强度增加明显,呈现较强的应变率敏感性[19]

2.3 试件冲击反应活性分析

图 8给出WZ-1试件在Ф14.5 mm SHPB冲击加载反应前后照片。通过照片可以看出反应十分剧烈,在t=103.5 ms时火光最为明亮 (图 8(d)),整个火光持续约为300 ms。3组W/Zr材料破坏应变很小,其在受到强冲击载荷作用下均会迅速破碎并发生反应。以火光大小为依据,相比较下WZ-1试件反应最为剧烈,WZ-2试件次之,WZ-3试件最弱。

图 8 HPB测试WZ-1试件反应典型图像 (a) t=0 ms;(b) t=34.5 ms;(c) t=69 ms;(d) t=103.5 ms;(e) t=310 ms;(f) t=550 ms Fig. 8 Typical video frames of WZ-1 sample reaction following a SHPB test (a) t=0 ms; (b) t=34.5 ms; (c) t=69 ms; (d) t=103.5 ms; (e) t=310 ms; (f) t=550 ms

为了进一步研究活性材料撞击反应行为,选用Ф5 mm SHPB实验平台对规格为Ф2 mm×2 mm试件进行加载并采用高速摄影技术进行观测。高速摄影设置中心触发方式,帧率为20000 fps,快门速率为1/50000 s。与以Al/PTFE为代表的金属/聚合物类活性材料撞击产生的爆轰效果不同,W/Zr活性材料冲击反应产生燃烧效果[20]图 9给出了WZ-1试件冲击反应的全过程。图 9(a)为试件加载前状态,在t=0.05 ms时,试件受强撞击已经开始发生破碎 (图 9(b)),而随着入射杆的继续推进,在t=0.1 ms时,试件已经完全压溃,外侧的碎片高速飞溅并伴有大量火星产生 (图 9(c))。对比图 9(b)图 9(c)中方框区域 (局部调高对比度处理后) 发现,在0.05 ms内碎片飞散的最远距离约为8 mm,也就是碎片运动的最高速率可达150 m/s。图 9(d)~(f)给出了碎片继续四处飞散的过程,在此阶段发现,较大的碎片 (如图 9(d)中圆圈所示) 一直保持黑暗,表明其没有反应。对比图 9(d)图 9(e)火星的拖影程度,可以判断碎片的运动速度明显降低。随着入射杆与透射杆分离,在两杆之间逐渐形成由细小碎片和粉末组成的云团,当该云团中活性金属粉末达到一定浓度范围时则十分容易被引燃。t=2.1 ms时,在如图 9(f)中圆圈所示区域出现一个十分微弱的火团并迅速向周围扩散 (图 9(g))。t=4.4 ms时,火团尺寸最大,如图 9(h)所示,在方框内的细小碎片云团浓度下降明显而未燃烧。活性金属燃烧产生的高温可使火团周围稍大的碎片表面融化,如图 9(i)方框中所示,一碎片受高温作用并粘连在入射杆端面。对比前后两组实验发现,在Ф14.5 mm SHPB平台上试件反应的持续时间要远大于在Ф5 mm SHPB平台上试件反应的持续时间。这主要是由于前者试件尺寸更大,撞击产生的碎片数量更多,而绝大部分碎片又被外部瓶体约束在瓶内形成高浓度的细小碎片云团,“燃料”更为充足且反应更加充分所致。

图 9 WZ-1试件冲击反应高速摄影截图 (a) t=0 ms;(b) t=0.05 ms;(c) t=0.10 ms;(d) t=0.35 ms;(e) t=1.25 ms;(f)t=2.10 ms;(g) t=2.80 ms;(h) t=4.40 ms;(i) t=5.70 ms Fig. 9 High-speed photographs of WZ-1 sample impact-induced reaction (a) t=0 ms; (b) t=0.05 ms; (c) t=0.10 ms; (d) t=0.35 ms; (e) t=1.25 ms; (f) t=2.10 ms; (g) t=2.80 ms; (h) t=4.40 ms; (i) t=5.70 ms

图 10给出了3组W/Zr材料在准静态和动态压缩 (Ф14.5 mm SHPB加载平台) 后的碎片照片。通过对比发现3组材料准静态压缩产生的碎片尺寸区别不大,但颜色差异明显,WZ-1~3碎片依次呈现银白色、银灰色和乌黑色。3组材料动态压缩产生的碎片尺寸明显都比准静态压缩产生的碎片小,WZ-1~3碎片尺寸依次减小,其中WZ-3碎片接近原始混合粉体状态。对3组动态碎片进行XRD测试,与原始试件成分对比分析发现W/Zr材料中能够参与反应的物质为Zr,ZrC,ZrC0.32H1.2,这些飞散的细小碎片与空气充分接触并剧烈燃烧生成ZrO2,而W和W2Zr则需要在更极端 (如爆炸) 加载条件下才能反应[21]

图 10 加载测试后W/Zr试件形貌 Fig. 10 Photographs of W/Zr sample after loading tests

图 11为WZ-1~3试件冲击反应后碎片的SEM照片。Zr粉在空气中燃烧产生的温度高达3400 ℃,而ZrO2的熔点约为2700 ℃。WZ-1试件部分直径为1~10 μm的碎片燃烧生成的ZrO2受高温熔化并在表面张力作用下最终呈现球形颗粒,如图 11(a)中圆圈所示[22],而尺寸更小 (直径小于500nm) 的碎片燃烧后形成的ZrO2附着在大碎片表面,呈现绒毛状。WZ-1试件原始粉末中的W全部用于形成结构尺寸更大的W2Zr合金,这使得WZ-1试件内各组分结合更为紧密,致使其动态碎片尺寸较大。WZ-2试件主要成分为ZrC0.32H1.2,而图 11(b)显示WZ-2试件碎片存在大量裂缝,且断口界面棱角分明,这证明由ZrC经氢化形成的ZrC0.32H1.2呈现明显的脆性。在WZ-2试件碎片表面同样可以观察到与WZ-1试件相似的ZrO2粉末和圆形颗粒,但数量明显偏少。图 11(c)显示WZ-3试件碎片的结构松散,与原始试件类似,这些大量缺陷致使WZ-3试件更容易压溃,所以其压缩强度在3组中最低,碎片尺寸最小。在WZ-3试件中,ZrC为主要活性物质,其在700 ℃时能与空气发生燃烧反应并产生2550 ℃高温,但由于未超过ZrO2的熔点,因此在WZ-3试件碎片表面未能观察到同前两组一样的球形ZrO2颗粒[23]

图 11 冲击加载后试件SEM微观形貌 (a) WZ-1;(b) WZ-2;(c) WZ-3 Fig. 11 SEM micrographs of W/Zr residues after dynamic tests (a) WZ-1; (b) WZ-2; (c) WZ-3
3 结论

(1) 采用热压烧结工艺,制备得到了3种配比的W/Zr活性材料,其相对密度为87.5%~99.2%。

(2) W/Zr活性材料是一种高强度脆性材料,其准静态压缩强度为1022~1880 MPa,破坏应变小于1%。动态恒应变率压缩实验时,3种W/Zr材料的压缩强度均呈现出一定的正应变效应。

(3) W/Zr活性材料在受到强冲击载荷作用时会与空气发生剧烈反应并生成ZrO2,与金属/聚合物类活性材料冲击发生的爆轰效果不同,W/Zr活性材料的反应类似燃烧效果。反应剧烈程度和持续时间除受试件原始成分制约外,也会受到破碎程度和碎片浓度影响。WZ-1试件力学性能和活性均优于另外两组。

参考文献(References)
[1] MONTGOMERY H E. Reactive fragment[P].United States Patent: US 3961576, 1976-06-08.
[2] THADHANI N N. Shock-induced and shock-assisted solid-state chemical reactions in powder mixtures[J]. Journal of Applied Physics, 1994, 76 (4): 2129–2138. DOI: 10.1063/1.357624
[3] National Research Council, Committee on Advanced Energetic Materials and Manufacturing Technologies. Advanced Energetic Materials[M]. Washington: National Academies Press, 2004. https://link.springer.com/article/10.1134/S0010508213060014?no-access=true
[4] AMES R G, WAGGENER S S. Reaction efficiencies for impact initiated energetic materials[C]//32nd International Pyrotechnics Seminar. Karlsruhe: International Pyrotechnics Society, 2005: 180/1-180/9.
[5] WILLIAM J F. Reactive fragment warhead for enhanced neutralization of mortar, rocket, and missile threats[R]. PA: De Technologies Inc, 2006.
[6] RITTER J J, BRANT A L, COLBURN J W. Characterization techniques employed to determine the energy release of reactive materials[R]. MD: Aberdeen Proving Ground, 2010. https://www.researchgate.net/publication/235133063_Characterization_Techniques_Employed_to_Determine_the_Energy_Release_of_Reactive_Materials
[7] AYDELOTTE B, BRAITHWAITE C H, MCNESBY K, et al. A study of fragmentation in a Ni+Al structural energetic material[J]. AIP Conference Proceedings, 2012, 1426 : 1097–1100.
[8] AMATO B, GOTZMER C, KIM S. Applications overview of IHDIV NSWC's reactive materials[R]. MD: Naval Surface Warfare Center, 2009. http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADA531456
[9] 张先锋, 赵晓宁. 多功能含能结构材料研究进展[J]. 含能材料, 2009, 17 (6): 731–739. ZHANG X F, ZHAO X N. Review on multifunctional energetic structural materials[J]. Journal of Energetic Materials, 2009, 17 (6): 731–739.
[10] 安亭, 赵凤起, 高红旭, 等. 超级铝热剂的制备及其与双基系推进剂组分的相容性[J]. 材料工程, 2011 (11): 23–28. AN T, ZHAO F Q, GAO H X, et al. Preparation of super thermites and their compatibilities with DB propellants components[J]. Journal of Materials Engineering, 2011 (11): 23–28. DOI: 10.3969/j.issn.1001-4381.2011.11.006
[11] NIELSON D B, TRUITT R M, RASMUSSEN N. Low temperature, extrudable, high-density reactive materials[P]. United States Patent: US 0020397A1, 2004-02-05. http://www.freepatentsonline.com/RE45899.html
[12] 王树山, 李朝君, 马晓飞. 钨合金破片对屏蔽装药撞击起爆的实验研究[J]. 兵工学报, 2001, 22 (2): 189–191. WANG S S, LI C J, MA X F. An experimental study on the initiation of covered charge impacted by tungsten[J]. Journal of China Ordnance, 2001, 22 (2): 189–191.
[13] 陈伟, 赵文天, 王健, 等. 钨锆合金破片毁伤过程研究[J]. 兵器材料科学与工程, 2009, 32 (2): 108–110. CHEN W, ZHAO W T, WANG J, et al. Study of damaging process for W-Zr alloy fragment[J]. Ordnance Material Science and Engineering, 2009, 32 (2): 108–110.
[14] ZHANG X F, SHI A S, QIAO L, et al. Experimental study on impact-initiated characters of multifunctional energetic structural materials[J]. Journal of Applied Physics, 2013, 113 (8): 083508. DOI: 10.1063/1.4793281
[15] LUO P G, WANG Z C, JIANG C L, et al. Experimental study on impact-initiated characters of W/Zr energetic fragments[J]. Materials & Design, 2015, 84 : 72–78.
[16] NAGENDER N S V, RAMA R P. Phase Diagrams of Binary Tungsten Alloys[M]. Calcutta: Indian Institute of Metals, 1991.
[17] OKAMOTO H. Desk Handbook: Phase Diagram for Binary Alloys[M]. American Society for Metals, 2000. https://www.valorebooks.com/textbooks/desk-handbook-phase-diagrams-for-binary-alloys/9780871706829
[18] 郭星渠, 赵彤, 王志英, 等. 用碳化锆细化晶粒制取大块氢化锆[J]. 原子能科学技术, 1980 (5): 641–644. GUO X Q, ZHAO T, WANG Z Y, et al. Preparation of large pieces of zirconium hydride from zirconium carbide by grain refinement[J]. Atomic Energy Science and Technology, 1980 (5): 641–644.
[19] 王鹏飞, 徐松林, 李志斌, 等. 微结构对多孔材料应变率效应影响的机理研究[J]. 爆炸与冲击, 2014, 34 (3): 285–291. WANG P F, XU S L, LI Z B, et al. Effect of micro-structure on the strain rate of cellular materials[J]. Explosion and Shock Waves, 2014, 34 (3): 285–291. DOI: 10.11883/1001-1455(2014)03-0285-07
[20] 任会兰, 李尉, 刘晓俊, 等. 钨颗粒增强铝/聚四氟乙烯材料的冲击反应特性[J]. 兵工学报, 2016 (5): 872–878. REN H L, LI W, LIU X J, et al. Reaction behaviors of Al/PTFE materials enhanced by W particles[J]. Acta Armamentarii, 2016 (5): 872–878.
[21] COVERDILL A, DELANEY C, JENNRICH A, et al. Tungsten combustion in explosively initiated W/Zr mechanical alloys[J]. Journal of Energetic Materials, 2014, 32 (3): 135–145. DOI: 10.1080/07370652.2013.795203
[22] BADIOLA C, DREIZIN E L. Combustion of micron-sized particles of titanium and zirconium[J]. Proceedings of the Combustion Institute, 2013, 34 : 2237–2243. DOI: 10.1016/j.proci.2012.05.089
[23] BAKER F B, STORMS E K, HOLLEY C E. Enthalpy of formation of zirconium carbide[J]. Journal of Chemical and Engineering Data, 1969, 14 (2): 244–246. DOI: 10.1021/je60041a034