Loading [MathJax]/jax/output/SVG/jax.js

冲击波作用下引信传爆序列殉爆的数值模拟

肖向东 肖有才 蒋海燕 范晨阳 王志军

肖向东, 肖有才, 蒋海燕, 范晨阳, 王志军. 冲击波作用下引信传爆序列殉爆的数值模拟[J]. 高压物理学报, 2021, 35(5): 054202. doi: 10.11858/gywlxb.20210706
引用本文: 肖向东, 肖有才, 蒋海燕, 范晨阳, 王志军. 冲击波作用下引信传爆序列殉爆的数值模拟[J]. 高压物理学报, 2021, 35(5): 054202. doi: 10.11858/gywlxb.20210706
XIAO Xiangdong, XIAO Youcai, JIANG Haiyan, FAN Chenyang, WANG Zhijun. Numerical Simulation and Analysis of Fuze Explosive Trains under Shock Waves[J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 054202. doi: 10.11858/gywlxb.20210706
Citation: XIAO Xiangdong, XIAO Youcai, JIANG Haiyan, FAN Chenyang, WANG Zhijun. Numerical Simulation and Analysis of Fuze Explosive Trains under Shock Waves[J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 054202. doi: 10.11858/gywlxb.20210706

冲击波作用下引信传爆序列殉爆的数值模拟

doi: 10.11858/gywlxb.20210706
基金项目: 国家自然科学基金(11802273);山西省青年科学基金(201901D211279);国防科工局技术基础科研项目(JCKY2XXXXX7B0XX)
详细信息
    作者简介:

    肖向东(1994-),男,硕士研究生,主要从事爆炸与冲击相关问题研究. E-mail:theexpectedfuture@126.com

    通讯作者:

    肖有才(1988-),男,博士,副教授,主要从事材料动态力学、损伤力学、爆炸与冲击等相关研究. E-mail:xiaoyoucai@nuc.edu.cn

  • 中图分类号: O381

Numerical Simulation and Analysis of Fuze Explosive Trains under Shock Waves

  • 摘要: 为研究不敏感弹药或不敏感引信在战备和后勤贮存过程中的殉爆现象,利用装填JH-14C传爆药的某引信传爆序列,开展了冲击波作用下的殉爆数值模拟研究,获得了引信传爆序列爆轰波成长历程、传播规律以及临界殉爆距离,建立了冲击波能量判据,并给出了殉爆条件。结果表明:冲击波由传爆药左上向右下传播,在最右下端面处起爆,引信传爆序列的临界殉爆距离为9.7 mm;当作用冲击波能量大于临界起爆能量时,引信传爆序列发生殉爆。

     

  • 引信传爆序列是爆轰能量传播的重要构件,序列中的含能材料通常具有较高的机械感度和冲击波感度,在外界作用下易发生爆炸或爆轰,导致殉爆现象发生[13]

    Kim等[45]通过建立流体力学模型对PBXN-9炸药开展了殉爆数值模拟研究,并将模拟结果与实验结果进行对比,得到了在一对多殉爆情况下,PBXN-9炸药殉爆影响因素。Mostafa等[6]利用一种轻型泡沫塑料对裸装炸药殉爆时的爆轰波历程进行了研究。陈朗等[78]计算了不同距离、不同工况下GHL炸药的殉爆距离,得出了殉爆过程中爆轰波的成长规律,但所研究的爆轰波传播规律仅针对单个被发炸药,并没有涉及一对多以及被发炸药相互作用下爆轰波传播规律分析。

    目前关于殉爆的研究很多,但主要针对裸装炸药或弹药殉爆,对不敏感引信传爆序列殉爆研究较少。本研究基于装填JH-14C传爆药的某引信传爆序列结构,拟建立一种引信传爆序列殉爆特性数值模拟方法,结合见证板实验,验证模型参数,并利用冲击波超压得到冲击波能量方程。

    起爆药的未反应和反应产物均采用Jones-Wilkins-Lee(JWL)状态方程描述

    p=A0eR1V+B0eR2V+ωCVTV (1)

    式中:p为压力,V为相对体积,T为温度,A0B0R1R2ωCV为常数。

    反应速率方程表示为

    dλdt=I(1λ)b(ρρ01a)x+G1(1λ)cλdpy0+G2(1λ)eλgpz0 (2)

    式中:λ为炸药反应度,t为时间,ρ为密度,ρ0为初始密度,p0为反应压强,IG1G2a、b、x、c、d、y、e、gz为拟合系数。

    传爆药的冲击起爆过程模拟分为3个阶段:(1)热点形成及其加热区点火阶段,(2)孤立热点向内或向外低速增长阶段,(3)反应快速完成阶段。

    Grüneisen状态方程是以Hugoniot曲线作为参考曲线,由S1S2S3这3个参数拟合而成的三次多项式,其压缩、拉伸材料方程分别表示为

    p=ρ0C20μ[1+(1γ02)μa2μ2]1(S11)μS2μ2μ+1S3μ3(μ+1)2+(γ0αμ)E (3)
    p=ρ0C20μ+(γ0+αμ)E (4)

    式中:EC0分别为初始内能和曲线的截距,μ=ρ/ρ01S1S2S3为冲击波速度-粒子速度(us-up)曲线斜率的系数,γ0为Grüneisen系数,αγ0的一阶体积修正。

    采用非线性有限元软件LS-DYNA3D,对冲击波作用下以JH-14C为传爆药的引信传爆序列进行一对多殉爆数值模拟。起爆点设置于起爆药上端中心点处,考虑爆炸冲击波对引信传爆序列的冲击作用,引信传爆序列结构如图1所示,其中T1T8为选取的高斯点,传爆序列中JH-14C的尺寸为38.0 mm × 15.5 mm。

    图  1  结构简图
    Figure  1.  Structure diagram

    图2(a)为使用True Grid软件建立的有限元模型,采用cm-g-μs单位制建模,模型中主发引信传爆序列和周围空气介质采用欧拉法描述,网格尺寸2 mm;为便于观察爆炸和变形情况,被发引信传爆序列A1、A2、A3和见证板采用拉格朗日法描述,网格尺寸1 mm,位置分布如图2(b)所示。

    图  2  位置分布与有限元计算模型
    Figure  2.  Location distribution and finite element model

    传爆序列殉爆数值模拟中,起爆药和传爆药分别为 Comp B、JH-14C,具体参数如表1[9-10]表2[11]所示,其中: G为剪切模量,σy为屈服强度,p0为炸药的单位体积爆轰压力,D为爆轰波速,FigmaxFG1maxFG2min分别为炸药反应点火时、炸药反应增长时和炸药反应完成时的最大反应分数,I为点火常数,G1G2分别为炸药反应增长和完成时的压力释放率。

    表  1  Comp B炸药的High Explosive Burn本构模型参数和JWL状态方程参数[910]
    Table  1.  High Explosive Burn constitutive model parameters and JWL equation parameters for Comp B[910]
    ρ/(g·cm−3)S1γ0A0/GPaB0/GPaR1R2ωp0/GPa
    1.7170.7980.33 524.237674.21.10.348.5
    下载: 导出CSV 
    | 显示表格
    表  2  JH-14C的Lee-Tarve参数[11]
    Table  2.  Lee-Tarve parameters of JH-14C[11]
    JWL equation parameters for unreacted explosive and reaction product
    StateA0/GPaB0/GPaR1R2ωG/GPaσy/MPaρ/(g·cm−3)p0/GPaD/(km·s−1)p/GPa
    Unreacted
    explosive
    7 781.0−5.3111.31.130.89383.5 2001.65
    Reaction
    product
    592.710.51 4.41.200.330 0 11.568.1927.67
    Rate of reaction
    a b c d e g x y z
    0.02 0.667 0.667 0.350 0.667 0.667 7.0 2.0 3.0
    Figmax FG1max FG2min I/μs−1 G1/(MPa·μs−1) G2/(MPa·μs−1)
    0.022 1 0 4 × 106 1.4 × 107 7 × 108
    下载: 导出CSV 
    | 显示表格

    雷管座和隔爆板材料为2024铝,保险机构为紫铜,见证板为4340钢,3种材料均利用Johnson-Cook本构方程和Grüneisen状态方程,参数如表3[12-13]表4[12-13]所示,其中:ABncm为材料常数;Tm为材料熔化温度;Tr为参考温度,通常取为室温;Cp为比定压热容,具体参数均源于Autodyn材料库。

    表  3  紫铜、2024铝与4340钢的Johnson-Cook本构模型参数[1213]
    Table  3.  Johnson-Cook constitutive model parameters for copper, 2024 aluminum and 4340 steel[1213]
    Materialρ/(g·cm−3)G/GPaA/MPaB/MPancmTm/KTr/KCp/(J·kg−1·K−1)
    Copper8.96046.0 902920.310.0251.091356300383
    2024 aluminum2.78528.62654260.340.0151.00 445300875
    4340 steel7.83077.07925100.260.0141.031793300477
    下载: 导出CSV 
    | 显示表格
    表  4  紫铜、2024铝、4340钢的Grüneisen状态方程参数[1213]
    Table  4.  Grüneisen equation of state parameters of copper, 2024 aluminum and 4340 steel[1213]
    MaterialC/(km·s−1)S1γ0
    Copper3.9401.4891.99
    2024 aluminum5.3281.3382.00
    4340 steel4.5691.4902.17
    下载: 导出CSV 
    | 显示表格

    为验证上述有限元模型以及参数的正确性,对单个引信传爆序列进行数值模拟,观察见证板的变形情况。有限元模型与实验装置如图3所示。

    图  3  有限元模型与实验装置
    Figure  3.  Finite element model and experimental apparatus

    4次实验中见证板的变形情况如图4所示。实验中采用雷管从引信传爆序列顶部引爆,见证板上凹坑深度分别为5.293、5.134、5.595、5.499 mm。图5为数值模拟结果与实验数据的对比,其中数值模拟得到的凹坑深度为5.355 mm。由图5可知,模拟结果与实验结果基本符合,凹坑深度的实验平均值与数值模拟结果的相对偏差约为0.47%,说明数值模型较好地描述了引信传爆序列爆轰行为。

    图  4  见证板变形
    Figure  4.  Witness plate deformation
    图  5  凹坑深度
    Figure  5.  Cave depths curve

    表5为不同距离下殉爆数值模拟结果,表格中距离如图2(b)所示。结果显示:当距离小于9.7 mm时,所有被发引信传爆序列完全爆轰;距离在9.7~10.3 mm之间时,A1和A2传爆药部分反应,A3未反应;距离大于10.3 mm时,被发引信传爆序列均未发生爆轰。

    表  5  不同距离下殉爆数值模拟结果
    Table  5.  Numerical simulation results at different distances
    Distance/mmNumerical simulation results
    A1A2A3
    8.5DetonationDetonationDetonation
    9.0DetonationDetonationDetonation
    9.7DetonationDetonationDetonation
    10.0Partial reactionPartial reactionNo detonation
    10.3No detonationNo detonationNo detonation
    10.5No detonationNo detonationNo detonation
    下载: 导出CSV 
    | 显示表格

    图6为距离9.0 mm时被发引信传爆序列A1殉爆压力云图。t = 10.0 μs时,主发引信传爆序列爆轰波作用于A1中的JH-14C,并从左上端开始向右下传播[14]。在此过程中,由于爆轰波的压缩作用,JH-14C内部的空气隙或气泡形成局部热点。t = 21.0 μs时,由于JH-14C内部热点汇聚,压力上升,在右下端面发生爆炸,产生的爆轰波从右下端向左上端传播,在极短时间内引爆整个引信传爆序列。

    图  6  距离9.0 mm 时A1殉爆压力云图
    Figure  6.  Pressure nephogram in the A1 at 9.0 mm

    在距离9.7 mm时,不同时刻被发引信传爆序列A3中传爆药变形和等压线如图7所示。t = 12.0 μs时,主发冲击波作用于A3;t = 23.0 μs时,由于A1和A2发生殉爆产生的冲击波作用于A3,因此从图6(c)中看出,A3在侧面发生殉爆,形成爆轰波,并在极短时间内赶上主发冲击波;t = 26.0 μs时,主发产生的冲击波与A3产生的爆轰波在传爆药右部发生叠加继续往右传播。

    图  7  距离9.7 mm时不同时刻A3中传爆药的变形和等压线
    Figure  7.  Deformation and isobaric lines of the explosive in the A3 at different times from 9.7 mm

    图8为距离为9.7 mm时,A1和A3的压力历程曲线。由图8(a)可知:t = 21.4 μs时,A1发生反应,高斯点T8处峰值压力约为6.4 GPa;随着反应进行,从t = 24.8 μs开始,压力值逐渐升高,并发展成稳定爆轰波。由图8(b)可知:t = 12.0 μs 时,主发引信传爆序列爆炸产生的冲击波开始作用于A3,此时压力为0.8 GPa,不足以使引信传爆序列发生反应,但A1和A2在主发引信传爆序列爆轰波作用下发生殉爆,产生的爆轰波作用于A3;t = 24.0 μs时,高斯点T1T2T3处压力发生明显的突跃,最高达到6.8 GPa,导致传爆药中热点增多,并产生局部高温区,A3发生爆炸反应,并且逐渐发展为稳定爆轰。

    图  8  距离为9.7 mm时A1和A3中的压力曲线
    Figure  8.  Column of pressure curves in the A1 and A3 at 9.7 mm

    图9(a)为距离10.0 mm时A1的压力历程曲线。由图9(a)可知,T8点峰值压力约为2.8 GPa,T7点峰值压力略有降低,T4点峰值压力接近8.5 GPa,而T3T2T1点压力呈下降趋势,说明部分炸药已发生反应,但没有形成稳定爆轰波。图9(b)为距离10.5 mm时A1的压力历程曲线。由图9(b)可知,t = 22.0 μs时,A1中传爆药发生反应,T8点压力值为2.0 GPa,随着反应进行,压力值并无上升反而呈现衰减趋势。

    图  9  A1中的压力曲线
    Figure  9.  Column of pressure curves in the A1

    本研究引入 Foan等[15]提出的衰减压力脉冲积分能量判据

    E=tw0p2(tw)dtw (5)

    由守恒方程和雨贡纽曲线关系

    Ec=p2twρus (6)

    式中:Ec为起爆能量,tw为脉冲宽度,us为压力p时的冲击波速度。

    入射冲击波压力达到一定强度时爆轰才可能发生,低于这一压力,爆轰不会发生,该压力称为临界起爆压力。采用作用于装药的临界起爆超压值Δp作为发生殉爆的判据

    Δp=p (7)

    则当外部爆轰波作用在引信传爆序列上的压力pout > Δp时,被发药发生殉爆。

    因而由式(6)、式(7)可得

    Ec=(Δp)2twρus (8)

    超压值Δp[16]

    Δp=1.02(3WR)+3.99(3WR)2+12.6(3WR)3 (9)

    其适用范围

    ¯H=H3W<0.35,1R3W1015 (10)

    式中:W为炸药当量;R为距装药中心距离,m;H为离地面高度,m。

    由式(8)、式(9)可得

    Ec=(0.37R2+0.53R+0.6R3)2twρus (11)

    由式(11)计算得到JH-14C的临界起爆能量Ec = 4.93 × 105 J/m2

    图10为引信传爆序列中高斯点的能量曲线。由图10(a)可知,在距离9.7 mm时,A1中高斯点T4T5处的能量高于临界起爆能量,发生殉爆。而在24.0 μs之前,A3中高斯点T4T5处的能量低于临界起爆能量,但A1与A2产生的冲击波作用使压力发生突跃,内部能量也相应上升,并迅速超过临界起爆能量,致使A3发生殉爆。由图10(b)可知,在距离为10.5 mm时, A1与A3中的高斯点T4T5处的能量低于临界起爆能量。

    图  10  引信传爆序列中高斯点能量
    Figure  10.  Gauss points energy of the fuze explosive train

    (1)Comp B为起爆药、JH-14C为传爆药的引信传爆序列起爆实验中,见证板的凹坑平均深度约为5.38 mm,与模拟计算结果5.35 mm基本吻合,由此验证了引信传爆序列模型参数的准确性。在此实际条件下,获得了引信传爆序列的临界殉爆距离为9.7 mm,殉爆安全距离为10.3 mm。

    (2)以JH-14C为传爆药的引信传爆序列殉爆过程中,冲击波由传爆药左端向右下传播,并在右下端面处起爆,之后爆轰波开始向左上传播,形成稳定爆轰波。

    (3)建立了以JH-14C(38.0 mm × 15.5 mm)为传爆药的引信传爆序列冲击波能量判据,获得了JH-14C的临界起爆能量Ec为4.93 × 105 J/m2。当作用冲击波能量大于Ec时,引信传爆序列发生殉爆反应。

  • 图  结构简图

    Figure  1.  Structure diagram

    图  位置分布与有限元计算模型

    Figure  2.  Location distribution and finite element model

    图  有限元模型与实验装置

    Figure  3.  Finite element model and experimental apparatus

    图  见证板变形

    Figure  4.  Witness plate deformation

    图  凹坑深度

    Figure  5.  Cave depths curve

    图  距离9.0 mm 时A1殉爆压力云图

    Figure  6.  Pressure nephogram in the A1 at 9.0 mm

    图  距离9.7 mm时不同时刻A3中传爆药的变形和等压线

    Figure  7.  Deformation and isobaric lines of the explosive in the A3 at different times from 9.7 mm

    图  距离为9.7 mm时A1和A3中的压力曲线

    Figure  8.  Column of pressure curves in the A1 and A3 at 9.7 mm

    图  A1中的压力曲线

    Figure  9.  Column of pressure curves in the A1

    图  10  引信传爆序列中高斯点能量

    Figure  10.  Gauss points energy of the fuze explosive train

    表  1  Comp B炸药的High Explosive Burn本构模型参数和JWL状态方程参数[910]

    Table  1.   High Explosive Burn constitutive model parameters and JWL equation parameters for Comp B[910]

    ρ/(g·cm−3)S1γ0A0/GPaB0/GPaR1R2ωp0/GPa
    1.7170.7980.33 524.237674.21.10.348.5
    下载: 导出CSV

    表  2  JH-14C的Lee-Tarve参数[11]

    Table  2.   Lee-Tarve parameters of JH-14C[11]

    JWL equation parameters for unreacted explosive and reaction product
    StateA0/GPaB0/GPaR1R2ωG/GPaσy/MPaρ/(g·cm−3)p0/GPaD/(km·s−1)p/GPa
    Unreacted
    explosive
    7 781.0−5.3111.31.130.89383.5 2001.65
    Reaction
    product
    592.710.51 4.41.200.330 0 11.568.1927.67
    Rate of reaction
    a b c d e g x y z
    0.02 0.667 0.667 0.350 0.667 0.667 7.0 2.0 3.0
    Figmax FG1max FG2min I/μs−1 G1/(MPa·μs−1) G2/(MPa·μs−1)
    0.022 1 0 4 × 106 1.4 × 107 7 × 108
    下载: 导出CSV

    表  3  紫铜、2024铝与4340钢的Johnson-Cook本构模型参数[1213]

    Table  3.   Johnson-Cook constitutive model parameters for copper, 2024 aluminum and 4340 steel[1213]

    Materialρ/(g·cm−3)G/GPaA/MPaB/MPancmTm/KTr/KCp/(J·kg−1·K−1)
    Copper8.96046.0 902920.310.0251.091356300383
    2024 aluminum2.78528.62654260.340.0151.00 445300875
    4340 steel7.83077.07925100.260.0141.031793300477
    下载: 导出CSV

    表  4  紫铜、2024铝、4340钢的Grüneisen状态方程参数[1213]

    Table  4.   Grüneisen equation of state parameters of copper, 2024 aluminum and 4340 steel[1213]

    MaterialC/(km·s−1)S1γ0
    Copper3.9401.4891.99
    2024 aluminum5.3281.3382.00
    4340 steel4.5691.4902.17
    下载: 导出CSV

    表  5  不同距离下殉爆数值模拟结果

    Table  5.   Numerical simulation results at different distances

    Distance/mmNumerical simulation results
    A1A2A3
    8.5DetonationDetonationDetonation
    9.0DetonationDetonationDetonation
    9.7DetonationDetonationDetonation
    10.0Partial reactionPartial reactionNo detonation
    10.3No detonationNo detonationNo detonation
    10.5No detonationNo detonationNo detonation
    下载: 导出CSV
  • [1] WANG F J, CHEN H M, MA C, et al. Construction of backscattering echo caused by cloud in laser fuze [J]. Optik, 2018, 171: 153–160. doi: 10.1016/j.ijleo.2018.06.028
    [2] SHARP A, ANDRADE J, RUFFINI N. Design for reliability for the high reliability fuze [J]. Reliability Engineering & System Safety, 2018, 181: 54–61.
    [3] 韩炎晖, 娄文忠, 冯跃, 等. 慢速烤燃环境下引信热响应特性测试与仿真 [J]. 兵工学报, 2019, 40(5): 946–953.

    HAN Y H, LUO W Z, YUE F, et al. Measurement and simulation of thermal response characteristics of fuze in slow cook-off test [J]. Acta Armamentarii, 2019, 40(5): 946–953.
    [4] KIM B, PARK J, LEE K C, et al. A reactive flow model for heavily aluminized cyclotrimethylene-trinitramine [J]. Journal of Applied Physics, 2014, 116(2): 023512. doi: 10.1063/1.4887811
    [5] KIM B, KIM M, SUN T, et al. Simulating sympathetic detonation using the hydrodynamic models and constitutive equations [J]. Journal of Mechanical Science and Technology, 2016, 30(12): 5491–5502. doi: 10.1007/s12206-016-1117-2
    [6] MOSTAFA H E, MEKKY W F, EL-DAKHAKHNI W W. Sympathetic detonation wave attenuation using polyurethane foam [J]. Journal of Materials in Civil Engineering, 2014, 26(8): 04014046. doi: 10.1061/(ASCE)MT.1943-5533.0000934
    [7] 陈朗, 王晨, 鲁建英, 等. 炸药殉爆实验和数值模拟 [J]. 北京理工大学学报, 2009, 29(6): 497–500, 524.

    CHEN L, WANG C, LU J Y, et al. Experiment & simulation of sympathetic detonation tests [J]. Transactions of Beijing Institute of Technology, 2009, 29(6): 497–500, 524.
    [8] CHEN L, WANG C, FENG C G, et al. Study on random initiation phenomenon for sympathetic detonation of explosive [J]. Defence Technology, 2013, 9(4): 224–228. doi: 10.1016/j.dt.2013.12.002
    [9] 刘鹏飞. 破片特性对冲击起爆B炸药比动能阈值的影响[D]. 太原: 中北大学, 2017.

    LIU P F. Influence of fragment characteristics on the threshold specific kinetic energy of impacting on covered Comp B [D]. Taiyuan: North University of China, 2017.
    [10] SCHWER L E. Impact and detonation of COMP-B: an example using the LS-DYNA EOS: ignition and growth of reaction in high explosives [C]//12th International LS-DYNA User Conference. Detroit, USA: 2012.
    [11] 李硕, 袁俊明, 刘玉存, 等. 聚黑-14C 的传爆装置冲击起爆实验及数值模拟 [J]. 火炸药学报, 2016, 39(6): 63–68, 79.

    LI S, YUAN J M, LIU Y C, et al. Experiment and numerical simulation of shock initiation of JH-14C detonation device [J]. Chinese Journal of Explosive & Propellants, 2016, 39(6): 63–68, 79.
    [12] BUYUK M, KURTARAN H, MARZOUGUI D, et al. Automated design of threats and shields under hypervelocity impacts by using successive optimization methodology [J]. International Journal of Impact Engineering, 2008, 35(12): 1449–1458. doi: 10.1016/j.ijimpeng.2008.07.057
    [13] 时党勇, 李裕春, 张胜民. 基于ANSYS/LS-DYNA 8.1进行显式动力分析 [M]. 北京: 清华大学出版社, 2005.

    SHI D Y, LI Y C, ZHANG S M. Explicit dynamics based on ANSYS/LS-DYNA8.1 force analysis [M]. Beijing: Tsinghua University Press, 2009.
    [14] 袁俊明, 李硕, 刘玉存, 等. 聚奥-9C装药的传爆管殉爆 [J]. 爆炸与冲击, 2018, 38(3): 632–638.

    YUAN J M, LI S, LIU Y C, et al. Sympathetic detonation of booster pipe with JO-9C charge [J]. Explosion and Shock Waves, 2018, 38(3): 632–638.
    [15] FOAN G C M, GOLEY G D. Shock initiation in gap test configurations [C]//7th Symposium on Detonation. Annapolis, USA: Naval Surface Weapons Center, 1981.
    [16] 王泽溥, 郑志良. 爆炸及其防护[M]. 北京: 兵器工业出版社, 2008.

    WANG Z P, ZHENG Z L. Explosion and protection [M]. Beijing: Ordnance Industry Press, 2008.
  • 加载中
图(10) / 表(5)
计量
  • 文章访问数:  1572
  • HTML全文浏览量:  683
  • PDF下载量:  48
出版历程
  • 收稿日期:  2021-01-10
  • 修回日期:  2021-02-02

目录

/

返回文章
返回