CN113985001A - 一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法 - Google Patents
一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法 Download PDFInfo
- Publication number
- CN113985001A CN113985001A CN202111158874.6A CN202111158874A CN113985001A CN 113985001 A CN113985001 A CN 113985001A CN 202111158874 A CN202111158874 A CN 202111158874A CN 113985001 A CN113985001 A CN 113985001A
- Authority
- CN
- China
- Prior art keywords
- metal salt
- total nitrogen
- anion metal
- nitrogen anion
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- IJGRMHOSHXDMSA-UHFFFAOYSA-N nitrogen Substances N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 title claims abstract description 221
- 229910052757 nitrogen Inorganic materials 0.000 title claims abstract description 166
- 239000002184 metal Substances 0.000 title claims abstract description 141
- 229910052751 metal Inorganic materials 0.000 title claims abstract description 141
- 150000003839 salts Chemical class 0.000 title claims abstract description 134
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000005474 detonation Methods 0.000 title claims abstract description 41
- 239000002360 explosive Substances 0.000 title abstract description 21
- -1 nitrogen anion Chemical class 0.000 claims abstract description 116
- 238000012360 testing method Methods 0.000 claims abstract description 64
- 238000004088 simulation Methods 0.000 claims abstract description 54
- 230000000977 initiatory effect Effects 0.000 claims abstract description 42
- 230000008569 process Effects 0.000 claims abstract description 30
- 238000011156 evaluation Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims description 26
- 238000000354 decomposition reaction Methods 0.000 claims description 15
- 150000001875 compounds Chemical class 0.000 claims description 9
- 238000009795 derivation Methods 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 abstract description 3
- 239000010949 copper Substances 0.000 description 12
- 238000006243 chemical reaction Methods 0.000 description 9
- 238000012937 correction Methods 0.000 description 9
- 125000004429 atom Chemical group 0.000 description 8
- 238000000329 molecular dynamics simulation Methods 0.000 description 8
- 230000003993 interaction Effects 0.000 description 7
- 238000000197 pyrolysis Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 6
- 125000000129 anionic group Chemical group 0.000 description 4
- 150000001540 azides Chemical class 0.000 description 4
- 239000013078 crystal Substances 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- IVRMZWNICZWHMI-UHFFFAOYSA-N Azide Chemical compound [N-]=[N+]=[N-] IVRMZWNICZWHMI-UHFFFAOYSA-N 0.000 description 2
- PXIPVTKHYLBLMZ-UHFFFAOYSA-N Sodium azide Chemical compound [Na+].[N-]=[N+]=[N-] PXIPVTKHYLBLMZ-UHFFFAOYSA-N 0.000 description 2
- 229910052783 alkali metal Inorganic materials 0.000 description 2
- 150000001340 alkali metals Chemical class 0.000 description 2
- 238000012018 process simulation test Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000002253 acid Substances 0.000 description 1
- 229910052784 alkaline earth metal Inorganic materials 0.000 description 1
- 150000001342 alkaline earth metals Chemical class 0.000 description 1
- 150000001450 anions Chemical class 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000001768 cations Chemical class 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 238000011407 dose test method Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- XLYOFNOQVPJJNP-ZSJDYOACSA-N heavy water Substances [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 description 1
- 239000003999 initiator Substances 0.000 description 1
- 150000002736 metal compounds Chemical class 0.000 description 1
- 150000004767 nitrides Chemical class 0.000 description 1
- 125000004433 nitrogen atom Chemical group N* 0.000 description 1
- 101150013860 oppB gene Proteins 0.000 description 1
- WUHLVXDDBHWHLQ-UHFFFAOYSA-N pentazole Chemical compound N=1N=NNN=1 WUHLVXDDBHWHLQ-UHFFFAOYSA-N 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- DTPQZKZONQKKSU-UHFFFAOYSA-N silver azanide silver Chemical compound [NH2-].[Ag].[Ag].[Ag+] DTPQZKZONQKKSU-UHFFFAOYSA-N 0.000 description 1
- 239000011734 sodium Substances 0.000 description 1
- 229910052723 transition metal Inorganic materials 0.000 description 1
- 150000003624 transition metals Chemical class 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/22—Fuels; Explosives
- G01N33/227—Explosives, e.g. combustive properties thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Computing Systems (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Algebra (AREA)
- Pathology (AREA)
- Food Science & Technology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明提供了一种适用于全氮阴离子金属盐起爆能力的评价方法，通过对全氮阴离子金属盐的起爆过程进行模拟试验，根据模拟试验的结果数据，计算获得T‑Km值，采用T‑Km值对全氮阴离子金属盐的起爆能力进行评价，T‑Km值的绝对值越大，则全氮阴离子金属盐的起爆能力越强。本发明的适用于全氮阴离子金属盐起爆药起爆能力的评价方法，通过对全氮阴离子金属盐的起爆过程进行模拟试验，能够避免危险性实验，获得的评价结果与现有文献中全氮阴离子金属盐起爆能力的实际测试结果相吻合，说明本发明的适用于全氮阴离子金属盐起爆能力的评价方法具有可靠性。因此，本发明的评价方法能够在避免危险性试验的前提下，从理论上可靠预测全氮阴离子金属盐起爆药的起爆能力。
Description
技术领域
本发明属于火炸药技术领域，涉及起爆药起爆能力评价，具体涉及一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法。
背景技术
全氮阴离子金属盐中金属阳离子包含碱金属、碱土金属、过渡金属和部分重金属。常见的全氮阴离子是叠氮根离子，其碱金属盐叠氮化钠不仅在汽车安全气囊中得到应用，而且是制备其他金属叠氮化物的关键原料之一。叠氮化铅是目前得到广泛应用的起爆药之一，其在雷管中的应用目前无可替代。叠氮化铜也因其优异的起爆性能而备受关注。根据现有文献的报道，关于五唑阴离子(N5 -)金属化合物，例如[Na(H2O)(N5)]·2H2O、Co(N5)2(H2O)4·4H2O、[Cu(N5)(N3)]n等，由于N原子组成的五唑阴离子环结构相比叠氮根离子的直线型结构有更大的结构张力，故N5 -阴离子金属盐具有成为起爆药的潜力。
衡量起爆药起爆能力的方法是极限起爆药量测试方法，通过该测试方法，国内外诸多学者已经测试并对比了氮化铅、氮化铜和氮化银的起爆能力。然而，目前国内外文献并没有报道新型的N5 -阴离子金属盐的起爆能力。因此，既无法通过实验验证N5 -阴离子金属盐是否能作为起爆药，又无法对比该新型化合物与传统叠氮金属起爆药的起爆能力。
发明内容
针对现有技术存在的不足，本发明的目的在于，提供一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法，解决现有技术中极限起爆药量测试方法无法评价全氮阴离子金属盐的起爆能力的技术问题。
为了解决上述技术问题，本发明采用如下技术方案予以实现：
一种适用于全氮阴离子金属盐起爆能力的评价方法，通过对全氮阴离子金属盐的起爆过程进行模拟试验，根据模拟试验的结果数据，拟合并计算后获得T-Km值，采用所述的T-Km值对全氮阴离子金属盐的起爆能力进行评价，T-Km值的绝对值越大，则全氮阴离子金属盐的起爆能力越强；
所述的T-Km值为拐点时间T值减去半分解期Km值，所述的半分解期Km值为模拟试验进行至全氮阴离子金属盐分解一半时所需要的时间，所述的拐点时间T值为，类米氏方程的曲率函数的极大值点的横坐标值；
所述的类米氏方程的曲率函数的表达式如下：
R代表氮气分子数量随时间变化的类米氏方程曲率半径；
TR代表模拟试验进行的时间，单位为ps。
本发明还具有如下技术特征：
具体的，该方法包括如下步骤：
步骤一，对全氮阴离子金属盐的起爆过程进行模拟试验，获得模拟试验的结果数据，即氮气分子数量随时间变化的平滑曲线；
步骤二，将步骤一中所述的氮气分子数量随时间变化的平滑曲线代入米氏方程中进行拟合，获得氮气分子数量随时间变化的类米氏方程和半分解期Km值；
步骤三，对步骤二中所述的氮气分子数量随时间变化的类米氏方程的进行求导，获得类米氏方程的曲率函数，再对类米氏方程的曲率函数进行求导，获得类米氏方程的曲率函数的极大值点，从而获得拐点时间T值；
步骤四，用步骤三中所述的拐点时间T值减去步骤二中所述的半分解期Km值，获得T-Km值并对全氮阴离子金属盐的起爆能力进行评价。
具体的，步骤二中，所述的氮气分子数量随时间变化的类米氏方程如下：
TR代表模拟试验进行的时间，单位为ps；
x代表氮气分子数量，单位为mol；
Nmax代表以曲线趋势拟合得到的理论氮气分子数量的极限值，单位为mol；
Km代表为模拟试验进行至全氮阴离子金属盐分解一半时所需要的时间，单位为ps。
具体的，所述的步骤一具体包括如下步骤：
步骤1.1，获得优化全氮阴离子金属盐晶胞：
构建全氮阴离子金属盐晶胞，对全氮阴离子金属盐晶胞的晶胞结构进行优化，获得优化后的全氮阴离子金属盐晶胞；
步骤1.2，构建全氮阴离子金属盐超胞：
选取步骤一中所述的优化后的全氮阴离子金属盐晶胞，构建全氮阴离子金属盐超胞，所述的全氮阴离子金属盐超胞中包含两个以上优化后的全氮阴离子金属盐晶胞；
步骤1.3，进行全氮阴离子金属盐超胞起爆过程的模拟试验：
对步骤二中所述的全氮阴离子金属盐超胞的起爆过程进行模拟试验，在模拟试验过程中，全氮阴离子金属盐超胞会生成氮气分子和金属团簇物；
步骤1.4，获得氮气分子数量随时间变化的平滑曲线：
步骤1.3所述的模拟试验结束后，绘制氮气分子数量随时间的实际变化曲线，将该实际变化曲线进行降噪平滑处理后，获得氮气分子数量随时间变化的平滑曲线。
具体的，步骤1.1中，所述的对全氮阴离子金属盐晶胞的晶胞结构进行优化时，采取的平面波截断能为500eV～1000eV。
具体的，步骤1.3中，所述的对全氮阴离子金属盐超胞的起爆过程进行模拟试验时，采取的平面波截断能的取值为步骤1.1中平面波截断能的取值的80％。
具体的，步骤1.3中，所述的模拟试验的过程为，先将全氮阴离子金属盐超胞在298K的温度下平衡1ps，随后在4000K的起爆温度下，以0.5fs的步长时间模拟5ps时间。
具体的，所述的全氮阴离子金属盐为AgN3、AgN5、CuN3、Cu(N3)2、CuN3N5和Pb(N3)2。
本发明与现有技术相比，具有如下技术效果：
本发明的评价方法能够对全氮阴离子金属盐起爆药的起爆能力进行评价，通过对全氮阴离子金属盐的起爆过程进行模拟试验，能够避免危险性实验；根据模拟试验的结果数据计算获得T-Km值，采用所述的T-Km值对全氮阴离子金属盐的起爆能力进行评价，获得的评价结果与现有文献中全氮阴离子金属盐起爆能力的实际测试结果相吻合，说明本发明的评价方法具有可靠性。本发明的评价方法能够在避免危险性试验的前提下，能够从理论上可靠预测全氮阴离子金属盐起爆药的起爆能力。
附图说明
图1为实施例1的全氮阴离子金属盐的起爆过程模拟试验的函数曲线图，其中，横坐标轴代表模拟试验进行的时间，纵坐标轴代表每摩尔全氮阴离子金属盐的氮气分子数量。
图2为实施例2的全氮阴离子金属盐的起爆过程模拟试验的函数曲线图，其中，横坐标轴代表模拟试验进行的时间，纵坐标轴代表每摩尔全氮阴离子金属盐的氮气分子数量。
图3为实施例3的全氮阴离子金属盐的起爆过程模拟试验的函数曲线图，其中，横坐标轴代表模拟试验进行的时间，纵坐标轴代表每摩尔全氮阴离子金属盐的氮气分子数量。
图4为实施例4的全氮阴离子金属盐的起爆过程模拟试验的函数曲线图，其中，横坐标轴代表模拟试验进行的时间，纵坐标轴代表每摩尔全氮阴离子金属盐的氮气分子数量。
图5为实施例5的全氮阴离子金属盐的起爆过程模拟试验的函数曲线图，其中，横坐标轴代表模拟试验进行的时间，纵坐标轴代表每摩尔全氮阴离子金属盐的氮气分子数量。
图6为实施例6的全氮阴离子金属盐的起爆过程模拟试验的函数曲线图，其中，横坐标轴代表模拟试验进行的时间，纵坐标轴代表每摩尔全氮阴离子金属盐的氮气分子数量。
以下结合实施例对本发明的具体内容作进一步详细解释说明。
具体实施方式
现有技术中，无法评价新型的N5 -阴离子金属盐的起爆能力的主要因素可能有以下几点：第一，极限起爆药量测试有较大风险，大部分科研院所不具备测试条件；第二，测试起爆药极限起爆药量需要一定量的产物，实验室合成的新型化合物通常是毫克量级，药量无法满足测试需求；第三，N5 -阴离子金属本身稳定性差，单质化合物在装填雷管过程中可能存在压爆现象导致无法测试。因此，针对于全氮阴离子金属盐起爆药起爆能力，研究一种既不具有危险性又具有可靠性的评价方法成为一个亟待解决的课题。
本发明中：
第一性原理分子动力学方法采用本领域常用的已知的第一性原理分子动力学方法。
范德华弱相互作用矫正方法DFT-D3法、DFT-D2法或optB88-vdW法均为现有技术中已知的范德华弱相互作用矫正方法；其中，DFT-D2法和DFT-D3法均源于文献《A consistentand accurate ab initio parametrization of density functionaldispersioncorrection(DFT-D)for the 94elements H-Pu》；optB88-vdW法源于文献《Chemicalaccuracy for the van der Waals density functional》。
ps为时间单位皮秒，fs为时间单位飞秒。
遵从上述技术方案，以下给出本发明的具体实施例，需要说明的是本发明并不局限于以下具体实施例，凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。
实施例1：
本实施例提供一种适用于全氮阴离子金属盐起爆能力的评价方法，本实施例中，全氮阴离子金属盐为AgN3。该方法包括如下步骤：
步骤一，采用分子动力学方法，对全氮阴离子金属盐的起爆过程进行模拟试验，获得氮气分子数量随时间变化的平滑曲线；步骤一具体包括如下步骤：
步骤1.1，获得优化全氮阴离子金属盐晶胞：
构建全氮阴离子金属盐晶胞，对全氮阴离子金属盐晶胞的晶胞结构进行优化，获得优化后的全氮阴离子金属盐晶胞；
本实施例中，采取第一性原理分子动力学方法对全氮阴离子金属盐晶胞的晶胞结构进行优化，步骤1.1中的第一性原理分子动力学方法中，平面波截断能为850eV；范德华弱相互作用矫正方法采用DFT-D3法，采用布里渊区划分方式k点个数；本实施例中，AgN3布里渊区划分方式k点采样个数为5×5×5；
本实施例中，第一性原理分子动力学方法是否合适的评判标准为：
第一，平面波截断能是否合适的评判标准为：若升高平面波截断能后分子晶体的能量变化小于0.001eV/atom，则升高前的平面波截断能的取值合适。
第二，范德华弱相互作用矫正方法是否合适的评判标准为：若采取该范德华弱相互作用矫正方法优化后得到的晶胞体积与X射线单晶衍射测得的晶胞体积，二者之间的偏差小于5％，则该范德华弱相互作用矫正方法合适。
第三，k点个数是否合适的评判标准为：若增加k点个数后分子晶体的能量变化小于0.0001eV/atom，则增加前的k点个数的取值合适。
步骤1.2，构建全氮阴离子金属盐超胞：
选取步骤一中的优化后的全氮阴离子金属盐晶胞，构建全氮阴离子金属盐超胞，的全氮阴离子金属盐超胞中包含两个以上优化后的全氮阴离子金属盐晶胞；本实施例中，AgN3超胞的尺寸为2×2×2，共包含128个原子。
步骤1.3，进行全氮阴离子金属盐超胞起爆过程的模拟试验：
对步骤二中的全氮阴离子金属盐超胞的起爆过程进行模拟试验，在模拟试验过程中，全氮阴离子金属盐超胞会生成氮气分子和金属团簇物；的模拟试验的过程为，先将全氮阴离子金属盐超胞在298K(即室温)的温度下平衡1ps，随后在4000K的起爆温度下，以0.5fs的步长时间模拟5ps时间。
本实施例中，采取第一性原理分子动力学方法对步骤二中的全氮阴离子金属盐超胞的起爆过程进行模拟试验，步骤1.3中的第一性原理分子动力学方法具体为：采用Born-Oppenheimer(BOMD)分子动力学方法、正则系综(NVT)和Nosé-Hoover恒温器，用Verlet算法积分牛顿运动方程，平面波截断能的取值为步骤1.1中平面波截断能的取值的80％，k点选用1×1×1的Gamma点，范德华弱相互作用矫正方式与步骤1.1保持一致。
步骤1.4，获得氮气分子数量随时间变化的平滑曲线：
步骤1.3的模拟试验结束后，绘制氮气分子数量随时间的实际变化曲线，如图1所示；将该实际变化曲线进行降噪平滑处理后，获得氮气分子数量随时间变化的平滑曲线，如图1所示。
本实施例中，绘制氮气分子数量随时间的实际变化曲线的过程为：统计全氮阴离子金属盐超胞中氮气分子数量随模拟时长的变化，用氮气分子摩尔数除以超胞中全氮阴离子金属盐的摩尔数，接着绘制每摩尔物质生成氮气数量随时间变化曲线。降噪平滑处理具体为，在Origin 9.0软件中采用快速傅里叶变化滤波器平滑处理得到的降噪曲线。
步骤二，将步骤一中的氮气分子数量随时间变化的平滑曲线代入米氏方程中进行拟合，获得氮气分子数量随时间变化的类米氏方程和半分解期Km值；氮气分子数量随时间变化的类米氏方程如下：
TR代表模拟试验进行的时间，单位为ps；
x代表氮气分子数量，单位为mol；
Nmax代表以曲线趋势拟合得到的理论氮气分子数量的极限值，单位为mol；
Km代表为模拟试验进行至全氮阴离子金属盐分解一半时所需要的时间，单位为ps。本实施例中，类米氏方程的拟合曲线如图1所示。
步骤三，对步骤二中的氮气分子数量随时间变化的类米氏方程的进行求导，获得类米氏方程的曲率函数，类米氏方程的曲率函数曲线如图1所示；再对类米氏方程的曲率函数进行求导，获得类米氏方程的曲率函数的极大值点，该极大值点的横坐标值即为拐点时间T值；类米氏方程的曲率函数的表达式如下：
R代表氮气分子数量随时间变化的类米氏方程曲率半径；
TR代表模拟试验进行的时间，单位为ps。
步骤四，用步骤三中的拐点时间T值减去步骤二中的半分解期Km值，获得T-Km值并对全氮阴离子金属盐的起爆能力进行评价，T-Km值如表1所示。
实施例2：
本实施例提供一种适用于全氮阴离子金属盐起爆能力的评价方法，本实施例中，全氮阴离子金属盐为AgN5。该方法包括如下步骤：
本实施例中，步骤一与实施例1的步骤一基本相同，区别在于，AgN5布里渊区划分方式k点采样个数为5×3×3；AgN5超胞的尺寸为2×1×1，共包含384个原子；氮气分子数量随时间的实际变化曲线和平滑曲线如图2所示。
本实施例中，步骤二与实施例1的步骤二相同。
本实施例中，步骤三与实施例1的步骤三相同。
本实施例中，步骤四与实施例1的步骤四相同，T-Km值如表1所示。
实施例3：
本实施例提供一种适用于全氮阴离子金属盐起爆能力的评价方法，本实施例中，全氮阴离子金属盐为CuN3。该方法包括如下步骤：
本实施例中，步骤一与实施例1的步骤一基本相同，区别在于，CuN3布里渊区划分方式k点采样个数为9×9×5；CuN3超胞的尺寸为2×2×2，共包含256个原子；氮气分子数量随时间的实际变化曲线和平滑曲线如图3所示。
本实施例中，步骤二与实施例1的步骤二相同；类米氏方程的拟合曲线如图3所示。
本实施例中，步骤三与实施例1的步骤三相同；类米氏方程的曲率函数曲线如图3所示。
本实施例中，步骤四与实施例1的步骤四相同，T-Km值如表1所示。
实施例4：
本实施例提供一种适用于全氮阴离子金属盐起爆能力的评价方法，本实施例中，全氮阴离子金属盐为Cu(N3)2。该方法包括如下步骤：
本实施例中，步骤一与实施例1的步骤一基本相同，区别在于，Cu(N3)2布里渊区划分方式k点采样个数为5×21×7；Cu(N3)2超胞的尺寸为1×4×2，共包含224个原子；氮气分子数量随时间的实际变化曲线和平滑曲线如图4所示。
本实施例中，步骤二与实施例1的步骤二相同；类米氏方程的拟合曲线如图4所示。
本实施例中，步骤三与实施例1的步骤三相同；类米氏方程的曲率函数曲线如图4所示。
本实施例中，步骤四与实施例1的步骤四相同，T-Km值如表1所示。
实施例5：
本实施例提供一种适用于全氮阴离子金属盐起爆能力的评价方法，本实施例中，全氮阴离子金属盐为CuN3N5。该方法包括如下步骤：
本实施例中，步骤一与实施例1的步骤一基本相同，区别在于，CuN3N5布里渊区划分方式k点采样个数为5×5×5；CuN3N5超胞的尺寸为2×2×2，共包含288个原子；氮气分子数量随时间的实际变化曲线和平滑曲线如图5所示。
本实施例中，步骤二与实施例1的步骤二相同；类米氏方程的拟合曲线如图5所示。
本实施例中，步骤三与实施例1的步骤三相同；类米氏方程的曲率函数曲线如图5所示。
本实施例中，步骤四与实施例1的步骤四相同，T-Km值如表1所示。
实施例6：
本实施例提供一种适用于全氮阴离子金属盐起爆能力的评价方法，本实施例中，全氮阴离子金属盐为Pb(N3)2。该方法包括如下步骤：
本实施例中，步骤一与实施例1的步骤一基本相同，区别在于，Pb(N3)2布里渊区划分方式k点采样个数为5×5×5；Pb(N3)2超胞的尺寸为2×1×1，共包含168个原子；氮气分子数量随时间的实际变化曲线和平滑曲线如图6所示。
本实施例中，步骤二与实施例1的步骤二相同；类米氏方程的拟合曲线如图6所示。
本实施例中，步骤三与实施例1的步骤三相同；类米氏方程的曲率函数曲线如图6所示。
本实施例中，步骤四与实施例1的步骤四相同，T-Km值如表1所示。
表1.实施例1至实施例6的T-Km值
由实施例1至实施例6可知：
(A)表1中的R2均接近1.0，说明实施例1至实施例6中的六种全氮阴离子金属盐的氮气分子数量随时间变化的平滑曲线拟合米氏方程的相关度较好。
(B)由表1中的Km值可知，实施例1至实施例6中，六种全氮阴离子金属盐快速热解反应前半段的反应速率的关系为：CuN3N5>AgN5>AgN3>Cu(N3)2>CuN3>Pb(N3)2。Km值越小，说明在起爆过程的模拟试验中，快速热解反应前半段的反应速率越快，即更容易分解产生氮气。由上述分析可知，N5 -比N3 -更容易分解产生氮气。
(C)由表1中的拐点时间T值可知，实施例1至实施例6中，六种全氮阴离子金属盐分解进程开始衰减的时间节点的先后顺序为：CuN3<AgN3<Pb(N3)2<AgN5<Cu(N3)2<CuN3N5。拐点时间T越小，说明快速热解反应进入衰减过程越早。
(D)当T-Km值为正值时，说明快速热解反应的进程衰减都出现在其半分解期之后，T-Km值为负值时，说明快速热解反应的进程衰减都出现在其半分解期之前。由表1中的T-Km值可知，实施例1至实施例6中，其中实施例5中的CuN3N5的T-Km值为正值，说明CuN3N5快速热解反应的进程衰减都出现在其半分解期之后；其余实施例中的全氮阴离子金属盐的T-Km值均为负值，说明其余五种全氮阴离子金属盐快速热解反应的进程衰减都出现在其半分解期之前。
T-Km值的绝对值越大，代表全氮阴离子金属盐的起爆能力越强。实施例1至实施例6中，六种全氮阴离子金属盐T-Km值的绝对值的大小按照从大到小的排序为：CuN3N5>Cu(N3)2>AgN5>AgN3>CuN3>Pb(N3)2，因此，实施例1至实施例6中，六种全氮阴离子金属盐的起爆能力按照从强到弱的排序为：CuN3N5>Cu(N3)2>AgN5>AgN3>CuN3>Pb(N3)2。
(E)现有研究中的极限起爆药量测试结果表明，Cu(N3)2的测试极限药量小于5mg，AgN3的测试极限药量大于5mg，Pb(N3)2的测试极限药量大于15mg，测试极限药量越小，则全氮阴离子金属盐的起爆能力越强，因此由上述研究可知，Cu(N3)2、AgN3和Pb(N3)2三者的起爆能力从强到弱的排序为：Cu(N3)2>AgN3>Pb(N3)2。上述测试结果与本发明中模拟试验的结果相符，说明本发明的适用于全氮阴离子金属盐起爆能力的评价方法具有可靠性。
根据本发明的结果能够推测出，实施例1至实施例6中，其中实施例5中的CuN3N5拥有最强的起爆能力，实施例2中的AgN5的起爆能力强于实施例1中的AgN3，但弱于实施例4中的Cu(N3)2，实施例3中CuN3的起爆能力强于实施例6中的Pb(N3)2。
综合上述分析可知，根据本发明的评价方法对全氮阴离子金属盐的起爆能力进行评价，其评价结果与现有文献中全氮阴离子金属盐起爆能力的实际测试结果相吻合，说明本发明的适用于全氮阴离子金属盐起爆能力的评价方法具有可靠性。
Claims (8)
1.一种适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，对全氮阴离子金属盐的起爆过程进行模拟试验后获得T-Km值，采用所述的T-Km值对全氮阴离子金属盐的起爆能力进行评价，T-Km值的绝对值越大，则全氮阴离子金属盐的起爆能力越强；
所述的T-Km值为拐点时间T值减去半分解期Km值，所述的半分解期Km值为模拟试验进行至全氮阴离子金属盐分解一半时所需要的时间，所述的拐点时间T值为，类米氏方程的曲率函数的极大值点的横坐标值；
所述的类米氏方程的曲率函数的表达式如下：
R代表氮气分子数量随时间变化的类米氏方程曲率半径；
TR代表模拟试验进行的时间，单位为ps。
2.如权利要求1所述的适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，该方法包括如下步骤：
步骤一，对全氮阴离子金属盐的起爆过程进行模拟试验，获得氮气分子数量随时间变化的平滑曲线；
步骤二，将步骤一中所述的氮气分子数量随时间变化的平滑曲线代入米氏方程中进行拟合，获得氮气分子数量随时间变化的类米氏方程和半分解期Km值；
步骤三，对步骤二中所述的氮气分子数量随时间变化的类米氏方程的进行求导，获得类米氏方程的曲率函数，再对类米氏方程的曲率函数进行求导，获得类米氏方程的曲率函数的极大值点，从而获得拐点时间T值；
步骤四，用步骤三中所述的拐点时间T值减去步骤二中所述的半分解期Km值，获得T-Km值。
3.如权利要求2所述的适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，步骤二中，所述的氮气分子数量随时间变化的类米氏方程如下：
TR代表模拟试验进行的时间，单位为ps；
x代表氮气分子数量，单位为mol；
Nmax代表以曲线趋势拟合得到的理论氮气分子数量的极限值，单位为mol；
Km代表为模拟试验进行至全氮阴离子金属盐分解一半时所需要的时间，单位为ps。
4.如权利要求2所述的适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，所述的步骤一具体包括如下步骤：
步骤1.1，获得优化全氮阴离子金属盐晶胞：
构建全氮阴离子金属盐晶胞，对全氮阴离子金属盐晶胞的晶胞结构进行优化，获得优化后的全氮阴离子金属盐晶胞；
步骤1.2，构建全氮阴离子金属盐超胞：
选取步骤一中所述的优化后的全氮阴离子金属盐晶胞，构建全氮阴离子金属盐超胞，所述的全氮阴离子金属盐超胞中包含两个以上优化后的全氮阴离子金属盐晶胞；
步骤1.3，进行全氮阴离子金属盐超胞起爆过程的模拟试验：
对步骤二中所述的全氮阴离子金属盐超胞的起爆过程进行模拟试验，在模拟试验过程中，全氮阴离子金属盐超胞会生成氮气分子和金属团簇物；
步骤1.4，获得氮气分子数量随时间变化的平滑曲线：
步骤1.3所述的模拟试验结束后，绘制氮气分子数量随时间的实际变化曲线，将该实际变化曲线进行降噪平滑处理后，获得氮气分子数量随时间变化的平滑曲线。
5.如权利要求4所述的适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，步骤1.1中，所述的对全氮阴离子金属盐晶胞的晶胞结构进行优化时，采取的平面波截断能为500eV～1000eV。
6.如权利要求5所述的适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，步骤1.3中，所述的对全氮阴离子金属盐超胞的起爆过程进行模拟试验时，采取的平面波截断能的取值为步骤1.1中平面波截断能的取值的80％。
7.如权利要求4所述的适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，步骤1.3中，所述的模拟试验的过程为，先将全氮阴离子金属盐超胞在298K的温度下平衡1ps，随后在4000K的起爆温度下，以0.5fs的步长时间模拟5ps时间。
8.如权利要求1所述的适用于全氮阴离子金属盐起爆能力的评价方法，其特征在于，所述的全氮阴离子金属盐为AgN3、AgN5、CuN3、Cu(N3)2、CuN3N5和Pb(N3)2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111158874.6A CN113985001B (zh) | 2021-09-30 | 2021-09-30 | 一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111158874.6A CN113985001B (zh) | 2021-09-30 | 2021-09-30 | 一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113985001A true CN113985001A (zh) | 2022-01-28 |
| CN113985001B CN113985001B (zh) | 2023-09-12 |
Family
ID=79737398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111158874.6A Active CN113985001B (zh) | 2021-09-30 | 2021-09-30 | 一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113985001B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE934694C (de) * | 1951-12-06 | 1955-11-03 | Nitroglycerin Aktiebolaget | Sprengstoffe auf der Grundlage von Nitroform |
| RU2002122484A (ru) * | 2002-08-19 | 2004-03-20 | Федеральное государственное унитарное предпри тие "Муромский приборостроительный завод" | Устройство для определения скорости детонации маломощных детонирующих шнуров типа "волновод" со светопроницаемой оболочкой |
| US20050247109A1 (en) * | 2002-03-28 | 2005-11-10 | Meyer Eric N | System and method for monitoring features of a blast |
| CN103558118A (zh) * | 2013-10-11 | 2014-02-05 | 西安近代化学研究所 | 基于内爆炸试验的温压炸药温压效应定量评价的测试方法 |
| RU141754U1 (ru) * | 2014-01-29 | 2014-06-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Самарский государственный технический университет" | Учебная лабораторная установка для исследования и демонстрации взрывных процессов |
| CN104034859A (zh) * | 2014-04-23 | 2014-09-10 | 西安近代化学研究所 | 炸药爆炸驱动金属能力的测量装置 |
| CN109696538A (zh) * | 2019-03-06 | 2019-04-30 | 西安近代化学研究所 | 一种量化评价多孔炸药起爆能力的方法 |
-
2021
- 2021-09-30 CN CN202111158874.6A patent/CN113985001B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE934694C (de) * | 1951-12-06 | 1955-11-03 | Nitroglycerin Aktiebolaget | Sprengstoffe auf der Grundlage von Nitroform |
| US20050247109A1 (en) * | 2002-03-28 | 2005-11-10 | Meyer Eric N | System and method for monitoring features of a blast |
| RU2002122484A (ru) * | 2002-08-19 | 2004-03-20 | Федеральное государственное унитарное предпри тие "Муромский приборостроительный завод" | Устройство для определения скорости детонации маломощных детонирующих шнуров типа "волновод" со светопроницаемой оболочкой |
| CN103558118A (zh) * | 2013-10-11 | 2014-02-05 | 西安近代化学研究所 | 基于内爆炸试验的温压炸药温压效应定量评价的测试方法 |
| RU141754U1 (ru) * | 2014-01-29 | 2014-06-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Самарский государственный технический университет" | Учебная лабораторная установка для исследования и демонстрации взрывных процессов |
| CN104034859A (zh) * | 2014-04-23 | 2014-09-10 | 西安近代化学研究所 | 炸药爆炸驱动金属能力的测量装置 |
| CN109696538A (zh) * | 2019-03-06 | 2019-04-30 | 西安近代化学研究所 | 一种量化评价多孔炸药起爆能力的方法 |
Non-Patent Citations (9)
| Title |
|---|
| QIAN-YU WANG等: "High-performance primary explosives derived from copper thiolate cluster-assembled materials for micro-initiating device", 《CHEMICAL ENGINEERING JOURNAL》, no. 389, pages 124455 * |
| 刘强, 许志东: "氯乙烯悬浮聚合引发剂的选择及聚合度预测", 齐鲁石油化工, no. 04, pages 45 - 46 * |
| 宋保永;胡振兴;孙;刘朋朋;汤忠斌;索涛;: "膨胀管分离装置分离时间测试研究", 实验力学, no. 02, pages 24 - 29 * |
| 张国英;王虹阳;杨利;: "Cu(N_3)_2及其复合结构的分子动力学模拟", 火工品, no. 01, pages 60 - 62 * |
| 张增志;谷娜;张际飞;: "瓦斯抑爆材料研究进展及瓦斯吸收剂初步研究", 煤矿安全, no. 08, pages 8 - 12 * |
| 张鹏;朱平安;马翰宇;王厚莲;杜保亭;: "炮射云爆战斗部二次起爆时间的试验研究", 火炸药学报, no. 05, pages 63 - 66 * |
| 武晨飞: "引信导爆药降解液释放机构的设计与其过载特性研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, pages 032 - 56 * |
| 肖绍清, 李胜, 林大泽: "采矿业预装药中使用高强度高精度导爆管起爆系统的可靠性和安全性", 安全, no. 04, pages 3 - 7 * |
| 赵宏安;胡荣祖;张晓亮;赵凤起;张海;高红旭;王喜军;冯煜;马海霞;: "非线性等转化率的微、积分法及其在含能材料物理化学研究中的应用Ⅶ. HHTDD的热分解", 西北大学学报(自然科学版), no. 04, pages 77 - 82 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113985001B (zh) | 2023-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wen et al. | Cluster evolution during the early stages of heating explosives and its relationship to sensitivity: a comparative study of TATB, β-HMX and PETN by molecular reactive force field simulations | |
| Tsai et al. | Comparison of the inhibition mechanisms of five types of inhibitors on spontaneous coal combustion | |
| Sobiczewski | Theoretical predictions for the nucleus 118 296 | |
| Muravyev et al. | HP-DSC study of energetic materials. Part I. Overview of pressure influence on thermal behavior | |
| Xiao et al. | Reaction mechanism of embedding oxidizing small molecules in energetic materials to improve the energy by reactive molecular dynamics simulations | |
| CN113985001A (zh) | 一种适用于全氮阴离子金属盐起爆药起爆能力的评价方法 | |
| Junying et al. | Reactive Molecular Dynamics Simulations of the Thermal Decomposition Mechanism of 1, 3, 3‐Trinitroazetidine | |
| Piernas-Muñoz et al. | Effect of temperature on capacity fade in silicon-rich anodes | |
| Ahmad et al. | Properties of Z= 120 nuclei and the α-decay chains of the 292, 304120 isotopes using relativistic and nonrelativistic formalisms | |
| CN106886615B (zh) | 一种rdx基多组分含能复合物的模拟方法 | |
| Xue et al. | Energy Performance and Aging of RDX‐based TiH2, MgH2 Explosive Composites | |
| Huber Jr et al. | The heat of combustion of calcium | |
| CN102539623A (zh) | 一种含硼富燃推进剂燃烧残渣精确定量分析方法 | |
| Wang et al. | Synthesis, characterization, thermal and explosive properties of potassium salts of trinitrophloroglucinol | |
| Yang et al. | Studies on the thermal behavior and safety of a novel thermostable explosive BPTAP | |
| CN106631642A (zh) | 含钛粉高能混合炸药及其制备方法 | |
| Baker et al. | Combined effects aluminized explosives modeling and development | |
| Dmitriev et al. | A study of the quality of preforms from iron-base powders produced by forming combined with sintering by electric current pulses | |
| Wang et al. | Preparation and performance of a novel water gel explosive containing expired propellant grains | |
| Makhov et al. | Method for estimating the acceleration ability of aluminized high explosives | |
| Keshavarz et al. | A simple way to predict electric spark sensitivity of nitramines | |
| Tan et al. | Effect of mixing methods on the thermal stability and detonation characteristics of ammonium nitrate and sodium chloride mixtures | |
| Ding et al. | Molecular dynamics study on the effect of metal ion doping on the performance of HMX | |
| Davydov et al. | Acceleration ability of high explosives and their mixtures with fuel additives. 1. Assessment by the M-40 method | |
| CN112033998A (zh) | 一种基于绝热加速量热法的火炸药物料热稳定性分级方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |