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

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

确定继续浏览么?

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

梯度硝基发射药制备过程安全风险分析  PDF

  • 梁昊 1
  • 王多良 1
  • 孙倩 2
  • 丁亚军 3
  • 李世影 3
  • 李洪伟 1
  • 肖忠良 3
1. 安徽理工大学化工与爆破学院, 安徽 淮南 232001; 2. 中国爆破器材行业协会, 北京 100089; 3. 南京理工大学化学与化工学院, 江苏 南京 210094

中图分类号: TJ55O64

最近更新:2024-08-22

DOI:10.11943/CJEM2024075

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

摘要

梯度硝基发射药(NGDP)作为一种新型发射药,具有高氧平衡、低发射有害现象以及抗迁移的特性。为了分析NGDP制备过程热危险性,利用反应量热仪(RC1mx)、差式扫描量热仪(DSC)和加速绝热量热仪(ARC)测定了双基球扁药脱硝反应过程中的放热以及反应前后物料的热分解行为。结果表明:反应体系的绝热温升(ΔTad)为58.28 ℃;脱硝反应使双基球扁药的初始分解温度从191.33 ℃提高至194.16 ℃;通过脱硝反应产物体系的ARC测试数据,获得了最大反应速度到达时间(TMRad)和温度之间的关系,计算出绝热诱导期为24h时所对应的温度TD24;反应的TD24、技术最高温度MTT与理论最高可达温度MTSR分别为137.2,108 ℃和128.52 ℃。基于精细化工反应安全风险评估规范,该反应在冷却失效的情况下,体系有冲料的风险,但触发二次放热分解反应的风险较低。

0 引 言

发射药能量释放的渐增性对提高身管武器性能起到了至关重要的作

1,为实现燃烧渐增性,一般通过钝感等手段来降低发射药表层燃速,但钝感剂为“惰性”组分,导致发射药负氧平衡加剧,枪口烟、焰与残渣等有害现象严2

梯度硝基发射药(NGDP)是一种新型发射药,其制备过程是采用化学方法将发射药表层含能的硝酸酯基(─O─NO2)梯度水解脱除,由于固液反应的扩散特性,脱硝试剂由表及里渗透进入发射药从而实现表层梯度硝基结构。Li

3-5通过理论计算探究了NGDP表层梯度硝基结构实现能量释放渐增性的原理,得出了表层硝基基团浓度梯度大小与其渐增性强弱直接相关的结论,并通过激光共聚焦显微拉曼光谱测定表层不同深度处硝基拉曼光谱信号强度,验证了经过脱硝处理后发射药表面的梯度硝基结构。肖忠良6指出,与传统发射药相比,NGDP具有贮存稳定性优越、内弹道性能良好、发射有害现象低等特点。崔璠7研究了反应温度、反应液浓度对NGDP脱硝率的影响,借助Avrami模型研究了NGDP制备过程的反应动力学,为NGDP生产工艺条件的建立和优化提供了理论依据。

火炸药生产的过程物料具有高能、高敏感性的特点,对安全生产的要求更

8‑9。研究人员针对不同的火炸药生产过程开展了大量研究,揭示了火炸药制备过程中的安全风险,对火炸药安全制备具有指导意义。周10等利用反应量热仪测定了1,1‑二氨基‑2,2‑二硝基乙烯(FOX‑7)合成过程中的反应热数据,并分析了反应中间体2‑(二硝基亚甲基)‑5,5‑二硝基嘧啶‑4,6‑二酮的热稳定性。陈利11等对二硝基甲苯(DNT)在经硝硫混酸硝化制取三硝基甲苯(TNT)的过程中的热危险性进行了研究,研究表明DNT硝化过程中放热剧烈,且反应过程中的热量有很大一部分来自于副反应放热。梯度硝基发射药制备过程的脱硝反应为典型的放热反应,对其安全性的研究是实现梯度硝基发射药工程化应用所急需且必要的。

反应热危险性研究为合成工艺安全控制提供了重要的理论基础和实验依据,对NGDP制备过程中的反应热危险性开展研究对其安全生产有重要意义。为此,以双基球扁药制备NGDP的过程为研究对象,利用RC1mx、DSC、ARC等多种分析测试手段,对NGDP制备过程的热危险性进行了分析,对比研究双基球扁药与NGDP的初始分解温度、分解峰温等热危险性指标,以期增加对NGDP制备、储存、运输过程安全风险特点的认识,寻求风险控制方法。

1 实验部分

1.1 试剂与仪器

双基球扁药,泸州北方化学工业有限公司提供,硝化纤维素含氮量13.15%,硝化甘油含量10%,弧厚2e1=0.32 mm,火药力f=1000 kJ·kg-1;水合肼(N2H4·H2O),分析纯,成都科隆化学品有限公司提供;去离子水,实验室自制。

仪器:全自动反应量热仪(RC1mx),配有2L常压玻璃反应釜及哈氏合金搅拌桨,梅特勒‑托利多公司,RC1mx可在接近实际工厂生产条件下,进行反应量热测试,记录并分析釜内反应的吸放热情况;控制及安全分析软件iControl;DSC‑HP高压差式扫描量热仪,梅特勒‑托利多公司;加速绝热量热仪(ARC),英国THT公司。

1.2 实验过程

双基球扁药表面脱硝制备NGDP的反应方程如下:

对于硝化纤维素:

C6H7O2OH(3-x)ONO2xn+n N2H4C6H7O2(OH)(4-x)ONO2x-1n+n H2N-NHNO2

对于硝化甘油:

C3H5(ONO2)3+n N2H4C3H5(ONO2)3-n(OH)n+n H2N-NHNO2

1.2.1 发射药脱硝反应的RC1mx测试

研究采用iControl软件编辑反应流程并发送至RC1mx反应量热仪,仪器通过设定的程序执行并记录实验过程中的放热情况。实验过程如下:(1)通过投料口投入200 mL质量分数为16%的水合肼水溶液进入反应釜,设置加热到40 ℃,并以400 r·min-1的速率搅拌10 min,待体系稳定后,执行标准量热校准程序,以确定反应体系的比热与传热系数;(2)校准程序完成后,读取反应釜内溶液体积,并在1 min内将100 g双基球扁药手动投入反应釜;(3)在30 min内将反应体系加热至80 ℃,保温反应2 h,反应结束后执行后校准程序,待体系降温至室温,取反应后固液混合体系以备测试。

1.2.2 发射药脱硝反应的原料和产物的DSC测试

测试温度为50~500 ℃,双基球扁药测试质量为0.52 mg,NGDP测试质量为0.48 mg,样品池均为高压不锈钢坩埚,升温速率为10 ℃·min-1,测试压力为常压0.1 MPa,测试气氛为氮气,流量50 mL·min-1

1.2.3 脱硝反应液及双基球扁药ARC测试

脱硝反应液测试:由于反应后为固液混合体系,为保证取样固液比例与反应釜内大致相符,取固体质量148 mg,液体320 mg,混合后注入ARC球内。测试模式为“加热‑等待‑搜寻”,升温台阶为5 ℃,等待时间为15 min,检测灵敏度0.02 ℃·min-1,检测温度范围为50~400 ℃,所用样品池为不锈钢量热球,球质量23.92 g。

双基球扁药测试:检测温度范围为50~350 ℃,球质量24.07 g,其他反应条件与脱硝反应液相同,具体反应参数如表1所示。

表1  双基球扁药及脱硝反应液的ARC测试条件
Table 1  ARC test conditions of double based gun propellant and denitration reaction fluids
ARC test materialstemperature range / ℃test cell mass / gspecific heat capacity / J·g-1·K-1mass / mg
denitration reaction fluids 50-400 23.92 3.51 508
double based gun propellant 50-350 24.07 1.50 99.7

2 结果与讨论

2.1 发射药脱硝过程的反应量热分析

参照文献[

12-15]中的方法,采用反应量热仪RC1mx测量记录发射药脱硝反应过程,通过夹套中硅油的温度变化结合校准测量的比热容和传热系数,可推导出反应体系的放热情况,实时监测反应过程中的体系放热速率。双基球扁药脱硝反应过程的反应量热曲线如图1所示。

图1  NGDP制备过程脱硝反应量热曲线

Fig. 1  Reaction calorimetry curves of denitrification reaction during NGDP preparation process

图1放热功率曲线Qr可以看出,该反应热释放属于加料控制型。在投入双基球扁药后,反应体系先经历了一段吸热过程,这是由于投入的物料温度低于反应体系的温度,反应放出的热量在此阶段用于将新投入物料加热。在投入双基球扁药后的升温段,体系的放热功率不断升高,在升温至工艺温度Tp时,放热功率达到峰值19.8 W。随着反应的进行,放热功率逐渐降低,反应至120 min时,放热功率趋于零,表明该时刻水合肼溶液向发射药内部扩散的速率较低,发射药表层的硝基几乎不再参与反应,因此可将120 min时视为脱硝反应的终点。

为得到反应过程的放热数据,将体系的实时放热功率Qr关于时间t积分,可得该条件下反应的总放热量ΔH,即

ΔH=Qrdt (1)

式中,ΔH为反应的总放热量,J;Qr为反应体系的实时放热功率,W;t为反应时间,s。

绝热温升ΔTad是反应风险评估的一个重要参数,可以反映出反应体系在冷却失效时,目标反应放热使反应体系升高的温度。ΔTad的值可通过式(2)计算

ΔTad=ΔHMrCp (2)

式中,Mr为反应体系的质量,kg;Cp为反应结束后体系的比热,J·g-1,可通过RC1mx反应量热校准数据得到。

结合反应的工艺温度Tp,利用式(3)计算出脱硝反应理论上可以达到的最高温度MTSR。

MTSR=Tp+ΔTad (3)

将脱硝反应的放热功率曲线关于时间积分,如图1中黄色区域所示,得到反应的总放热量ΔH为61.46 kJ,根据式(2)可计算出反应的绝热温升ΔTad为58.28 ℃,比放热量为614.56 kJ·kg-1(以双基球扁药质量计)。根据式(3),若采用一次性加料方式,反应的MTSR为128.52 ℃,此温度高于反应体系的技术最高温度(MTT),反应在常压下进行,MTT即反应体系的沸点108 ℃。因此在该脱硝反应中,若冷却装置失效,会导致反应体系局部高温,反应液剧烈汽化,进而形成汽液相混合体喷出,造成冲料事故。在工厂化生产中,应在反应釜合理位置设置进料口,并在釜体加装爆破片,爆破片连接管应通向室外安全处,以防物料喷溅伤人。

2.2 双基球扁药与NGDP的差示扫描量热实验分析

原料双基球扁药与NGDP的DSC曲线如图2所示,由于测试中50~100 ℃和300~500 ℃温度区间未扫描到放热峰,故未在图2体现。以DSC曲线上热流斜率最大点的切线与外延基线的交点所对应的温度作为试样的初始分解温度,由图2可以看到双基球扁药的初始分解温度为191.33 ℃,第一阶段的分解峰温为199.54 ℃,第二阶段的分解峰温为205.17 ℃。梯度脱硝发射药的初始分解温度为194.16 ℃,第一阶段的分解峰温为197.93 ℃。第二个分解峰温约为204.16 ℃,与双基球扁药相比,NGDP的第二个分解峰并不明显,这是由于脱硝处理使发射药表面含能的硝酸酯基团被还原为羟基基团,发射药对温度敏感性降低,使第一个放热峰后延,与第二个放热峰重叠。结果表明,脱硝处理使双基球扁药的初始分解温度提高,分解峰温降低,稳定性和安全性有一定程度的提高,与产物梯度硝基发射药相比,原料双基球扁药在贮存阶段的危险程度更高。

图2  双基球扁药和梯度硝基发射药的DSC曲线

Fig.2  DSC curves of double based propellant and nitro gradiently distributed propellant

2.3 双基球扁药绝热量热分析

由于DSC测试中双基球扁药表现出较低的初始分解温度,因此,进一步对双基球扁药进行绝热量热测试,通过逐台阶缓慢升温,能更为精准的确定试样起始分解温

16,结果如图3所示。由图3可以看出,双基球扁药呈现2个明显的放热阶段,这与图2中DSC曲线出现2个放热峰的情况相符,第一个放热阶段起始于140.15 ℃,结束于180.87 ℃,绝热温升为40.72 ℃,该阶段体系温升速率由慢变快,反应体现出明显的自加速分解放热特点。第二个放热阶段,起始于225.10 ℃,结束于240.53 ℃,绝热温升为15.43 ℃,该阶段的放热速率均匀,推测是随着温度的升高,第一个放热阶段的反应产物在密闭的ARC测试球中发生了放热的水煤气变换反2

图3  双基球扁药的绝热量热测试曲线

Fig.3  ARC testing curve of double based gun propellant

结合2个放热阶段得到的绝热温升,计算得双基球扁药的第一阶段的分解热为4177.62 J·g-1,第二阶段的反应放热为1583.12 J·g-1。双基球扁药的起始分解温度为140.15 ℃,高于脱硝反应体系的MTT 108 ℃,在双基球扁药投入反应釜形成脱硝反应体系后,便无法达到其初始分解温度。NGDP的制备过程中,原料在储存运输阶段的温度敏感程度高于反应体系,工业化生产时,应在原料输送管道加装温度实时监测装置,防止原料在密闭管道中产生热积累导致分解。

2.4 脱硝反应液绝热分解分析

对脱硝反应液进行绝热量热测试,结果见表2

表2  脱硝反应液ARC结果
Table 2  Denitration reaction fluids ARC test results
parametersdata
sample mass / g 0.508
specific heat capacity of sample / J·g-1·K-1 3.51
test cell mass / g 23.923
specific heat capacity of test cell /J·g-1·K-1 0.42
thermal correction coefficient Φ 6.72
onset temperature / ℃ 145.64
final temperature / ℃ 157.45
adiabatic temperature rise / ℃ 11.81
adiabatic temperature rise corrected value / ℃ 79.36
max pressure / MPa 87.18
max temp rate / ℃·min-1 0.072
max pressure rate / MPa·min-1 0.078
specific heat release / J·g-1 261.03

对仪器测得的绝热温升进行修正,试样和容器处于绝热状态时,

16

msCsΔTs=msCs+mcCbΔT (4)

式中,msmc分别代表试样、样品池的质量,g;CsCb分别代表试样、样品池的比热容,J·g-1;ΔTs代表物料放热仅加热试样时温度的升高值,℃;ΔT代表物料放热加热试样及容器时的实测温度升高值,℃。

式(4)得到

ΔTs=1+mbCbmsCsΔT (5)

式(5)中的(1+mbCbmsCs)=ΦΦ称为容器的修正系数,对于ARC测试中的实际绝热温升ΔTs,应由绝热温升测量值修正得到,即ΔTs·ΔT,由表2可计算得Φ=6.72,实验测得的绝热温升ΔT为11.81 ℃,经Φ值修正后的绝热温升ΔTs为79.36 ℃。

ARC测试得到的脱硝反应液放热段的温度、温升速率‑时间曲线和压力、压升速率‑温度曲线如图4~5所示。由图4可以看出,脱硝反应液在145.64 ℃开始自放热分解,放热终点位于157.45 ℃。在分解反应进行了268 min后,温升速率达到最大值0.049 ℃·min-1,此时的体系温度为153 ℃,体系压力为1.05 MPa,在脱硝反应液的整个自放热分解过程中,放热速率在148.5 ℃和157.45 ℃ 2个阶段达到高峰,体现出一定的自加速分解特点。由图5可以看出,在整个自分解过程中,体系的压升速率虽有波动,但无明显变化趋势,体系压力近乎线性升高。工业化生产时,即使反应体系发生热失控导致自加速分解,釜内压力骤升的概率也较小,通过设置合理泄压阀和爆破片,可规避反应釜爆炸风险。

图4  脱硝反应液放热段的温度、温升速率与时间的曲线

Fig.4  Temperature, temperature rate and time curve of the denitration reaction solution

图5  脱硝反应液放热段的压力、压升速率与温度的曲线

Fig.5  Pressure, pressure rate and temperature curve of denitration reaction solution

对于绝热条件下的物质反应,一般用阿伦尼乌斯(Arrhenius)公

17描述反应速度:

dcdt=-Acnexp-ERT (6)

式中,c为未反应物料的浓度,mol·L-1A代表指前因子,s-1n为反应级数,E为反应活化能,J·mol-1,R为摩尔气体常数,8.314 J·mol-1·K-1

若分解反应放出的热量无法向外界传递,产生的热量会将反应体系再度加热,有绝热方程:

Hdcdt=-CpdTdt (7)

式中,H为反应热,J·mol-1Cp为体系的平均比热容,J·g-1·K-1

引入下列式(8)、(9)中的边界条

18

t=0时:

c=c0,T=T0 (8)

t=∞时:

c=0,T=Tf (9)

式中,Tf为绝热条件下所达到的最高温度,℃;T0为反应的初始温度,℃;c0为反应物料的初始浓度,mol·L-1

HCpc00dc=-T0TfdT (10)
HCp=Tf-T0c0=ΔTc0 (11)

式中,ΔT为绝热状态下反应体系温度的升高值,℃。

样品的浓度随温度变化,样品分解放热又导致其自身温度升高,所以样品的浓度c与温度T存在以下关系:

c=Tf-TΔTc0 (12)

结合式(6)、(7)、(11)、(12)有

dTdt=ATf-TΔTnΔTc0n-1exp-ERT (13)

k*=dTdtΔTTf-TnΔT-1 (14)

结合式(13)、(14)有

k*=Ac0n-1exp-ERT (15)

式(15)两边同时取对数,则Arrhenius方程的对数形式为:

lnk*=ln(Ac0n-1)-ERT (16)

针对反应过程,选定适合的反应级数n,使得lnk*对1/T的曲线为一条直线,曲线的斜率和截距即为对应反应的EA。通过反应动力学软件计算,获得脱硝反应液的绝热分解动力学参数为:反应级数n=0.59,反应活化能E为2.08×105 J·mol-1,指前因子A=4.2×1022 s-1

最大反应速率到达时间TMR(Time to Maximum Rate),也称TMRad,是指在绝热条件下,试样从任意温度T开始到达最大反应速率对应时刻tm所经历的时间,可表示为:

TMRad=ttmdt=RT2(dT/dt)E-RTm2(dT/dt)mE (17)

式中最大升温速率(dT/dtm远大于(dT/dt),后一项比前一项小得多,忽略后一项后,即为:

TMRad=RT2(dT/dt)E (18)

结合脱硝反应液的绝热分解动力学参数以及ARC测试结果,计算得到的TMRadT曲线如图6所示。化工生产中,常将TMRad为24 h时所对应的温度作为反应物料的初始分解温度,通过外推法获得脱硝反应液的TD24为137.2 ℃,该温度高于反应的MTSR=128.52 ℃,因此在冷却失效的条件下,引发二次反应的风险较低。

图6  脱硝反应液的最大反应速率到达时间与温度关系曲线

Fig.6  TMRadT curve of denitration reaction fluids

2.5 脱硝反应制备NGDP的工艺热危险分析及改进措施

发射药脱硝反应是典型的非均相放热反应,传热系数决定了反应体系和外界交换热量的速率,在规模化生产试验中可采取增加搅拌速率、增加搅拌桨叶层数、加装挡板等方式加快反应体系的流动速率,提高反应体系的热均匀性,避免局部进入绝热状态。

结合反应量热数据可知,在投料时反应体系即开始放热,目前的反应工艺流程为一次性手动投料,投料后即升温,在升温过程中体系的放热功率不断加快,可见脱硝反应的放热功率与温度成正相关,在升温到保温反应温度时,放热功率达到峰值19.8 W。在工业化生产中可采用缓慢升温的方式来降低冷却系统的压力,固体投料也可采用分批次、均匀的自动进料方式,以延长加料时间,防止在加料过程中产生热积累。

此外,由于反应的工艺温度Tp较为接近体系的MTT,且反应体系中水合肼溶液具有一定的挥发性,在进一步的生产实验中应配备冷凝回流装置,以提高反应物料的利用率。

从实验中得到的安全性参数可知,反应的MTSR为128.52 ℃,高于反应体系水合肼溶液的沸点,反应体系在冷却失效的情况下,会先升温至体系沸点,造成冲料,但由于MTSR小于TD24,发生二次放热分解反应的风险较低。在规模化生产试验中,需设置可靠的大功率冷却降温系统,以规避物料冲料的风险。

3 结 论

(1)以双基球扁药为原料,经脱硝反应制备NGDP的反应过程中的比放热量为614.56 kJ·kg-1(以双基球扁药质量计),绝热温升ΔTad为58.28 ℃。

(2)DSC测试结果表明双基球扁药与NGDP的放热均存在着两个阶段,双基球扁药的初始分解温度低于NGDP,DSC测试中,双基球扁药的初始分解温度为191.33 ℃,ARC测试中双基球扁药的初始放热温度为140.15 ℃,这是由两种测试手段的升温方式的差异所导致的,应以ARC测试中的初始放热温度作为物料的初始分解温度。

(3)脱硝反应液的初始分解温度为145.64 ℃,对分解放热数据进行动力学计算,可知分解反应的反应级数n为0.59,反应活化能E为2.08×105 J·mol-1,指前因子A为4.2×1022 s-1,脱硝反应液的TD24为137.2 ℃,结合反应体系的MTT=108 ℃和MTSR=128.52 ℃,推断出该反应在冷却失效的情况下,体系有冲料的风险,但通过TD24=137.2 ℃判断,体系触发二次分解反应的风险较低。

参考文献

1

王泽山. 火药装药设计原理与技术[M].北京北京理工大学出版社2006. [百度学术] 

WANG Ze‑shanet al. Gunpowder Charge Design Principles and Techniques[M].BeijingBeijing Institute of Technology Press2006. [百度学术] 

2

肖忠良. 火炸药导论[M].北京国防工业出版社2019. [百度学术] 

XIAO Zhong‑liang. Introduction to Fire Explosives[M].BeijingNational Defence Industry Press2019. [百度学术] 

3

LI SCHEN HLI Yet al. Surface denitration structure on dynamic combustion performance and muzzle flame of mixed nitrate gun propellant[J]. Propellants, Explosives, Pyrotechnics2023489). [百度学术] 

4

LI STAO ZDING Yet al. Gradient denitration strategy eliminates phthalates associated potential hazards during gun propellant production and application[J]. Propellants, Explosives, Pyrotechnics2020457): 1156-1167. [百度学术] 

5

LI SLI YDING Yet al. One‑step green method to prepare progressive burning gun propellant through gradient denitration strategy[J]. Defence Technology202122135-143. [百度学术] 

6

肖忠良丁亚军李世影. 发射药表层梯度硝基裁剪方法与效应[J]. 中国材料进展2022412): 92-97. [百度学术] 

XIAO Zhong‑liangDING Ya‑junLI Shi‑yinget al. Gradient Tailoring method and effect of nitro gradiently distributed propellant[J]. Materials China2022412): 92-97. [百度学术] 

7

崔璠范红蕾李世影.梯度硝基单基发射药构筑过程的avrami模型[J]. 兵工学报DOI:10.12382/bgxb.2023.1190. [百度学术] 

CUI FanFAN Hong‑leiLI Shi‑yinget al. Avrami model for the construction process of nitro gadiently distributed single‑base propellant[J]. Acta ArmamentariiDOI:10.12382/bgxb.2023.1190. [百度学术] 

8

丁黎. 剃刀边缘的角逐——火炸药工艺安全风险控制难点分析[J]. 火炸药学报20234610): 851-854. [百度学术] 

DING Li. Analysis of difficulties in safety risk control of explosives process[J]. Chinese Journal of Explosives & Propellants20234610): 851-854. [百度学术] 

9

弗朗西斯·施特塞尔陈网桦彭金华. 化工工艺的热安全风险评估与工艺设计[M]. 北京科学出版社2009. [百度学术] 

FRANCIS·StoesselCHEN Wang‑huaPENG Jin‑huaet al. Thermal Safety of Chemical Processes ‑ Risk Assessment and Process Design[M]. BeijingScience Press2009. [百度学术] 

10

周诚朱勇王伯周. FOX‑7合成过程中硝化反应的热危险性[J]. 含能材料20141): 53-56. [百度学术] 

ZHOU ChengZHU YongWANG Bo‑zhouet al.Thermal hazards of nitration reaction in the synthesis of FOX‑7[J].Chinese Journal of Energetic Materials(Hanneng Cailiao)20141): 53-56. [百度学术] 

11

陈利平陈网桦彭金华. 二硝基甲苯硝化反应的热危险性分析[J]. 含能材料2010186): 706-710. [百度学术] 

CHEN Li‑pingCHEN Wang‑huaPENG Jin‑huaet al.Thermal hazard analysis of dinitrotoluene nitration[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao)2010186):706-710. [百度学术] 

12

ZHAO JZHANG WHU Jet al. Research on the risk of thermal runaway in the industrial process of styrene solution polymerization[J]. Organic Process Research & Development2021256): 1366-1374. [百度学术] 

13

常佩周诚王伯周. 二种NTO合成工艺的安全性分析[J]. 高校化学工程学报2018325): 1223-1227. [百度学术] 

CHANG PeiZHOU ChengWANG Bo‑zhouet al.Safety analysis of two NTO synthesis processes[J]. Journal of Chemical Engineering of Chinese Universities2018325): 1223-1227. [百度学术] 

14

朱勇王玉刘建利. ADN合成过程的热效应[J]. 含能材料2012206): 735-738. [百度学术] 

ZHU YongWANG YuLIU Jian‑liet al. Thermal effects of synthesis process of ADN[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao)2012206): 735-738. [百度学术] 

15

吴展华盛敏. 绝热加速量热仪在反应安全风险评估应用中的常见问题[J]. 化工进展2023427): 3374-3382. [百度学术] 

WU Zhan‑huaSHENG Min. Pitfalls of accelerating rate calorimeter for reactivity hazard evaluation and risk assessment[J]. Chemical Industry and Engineering Progress2023427):3374-3382. [百度学术] 

16

WANG YLIU JCHEN Let al. Thermal decomposition characteristics and runaway boundary conditions of hato at adiabatic and high pressure situations[J]. Process Safety and Environmental Protection2022167601-608. [百度学术] 

17

贾昊楠路桂娥江路明. 用加速量热仪研究双基发射药的绝热分解特性[J]. 火炸药学报2013365): 77-81. [百度学术] 

JIA Hao‑nanLU Gui‑eJIANG Lu‑minget al. Study on adiabatic decomposition properties of double‑base gun propellant by accelerating rate calorimeter[J]. Chinese Journal of Explosives & Propellants2013365): 77-81. [百度学术] 

18

冯长根傅智敏钱新明. 绝热方法评价热安定性的改进模型与应用[J].北京理工大学学报200301): 22-25. [百度学术] 

FENG Chang‑genFU Zhi‑minQIAN Xin‑minget al. Improved models and applications of adiabatic method for evaluating the thermal stability of chemicals[J]. Transactions of Beijing Institute of Technology200301): 22-25. [百度学术]