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

极小装药爆炸驱动两级重载荷研究

李鹏飞 杨磊 赵向军 肖川 宋浦

缪广红, 李亮, 江向阳, 刘文震, 李雪交, 汪泉, 余勇, 沈兆武. 双面爆炸焊接的数值模拟[J]. 高压物理学报, 2018, 32(4): 045202. doi: 10.11858/gywlxb.20180513
引用本文: 李鹏飞, 杨磊, 赵向军, 肖川, 宋浦. 极小装药爆炸驱动两级重载荷研究[J]. 高压物理学报, 2018, 32(6): 065102. doi: 10.11858/gywlxb.20180581
MIAO Guanghong, LI Liang, JIANG Xiangyang, LIU Wenzhen, LI Xuejiao, WANG Quan, YU Yong, SHEN Zhaowu. Numerical Simulation of Double-Sided Explosive Welding[J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 045202. doi: 10.11858/gywlxb.20180513
Citation: LI Pengfei, YANG Lei, ZHAO Xiangjun, XIAO Chuan, SONG Pu. Study on Minimal Charge Explosive Driving Two Mass Load[J]. Chinese Journal of High Pressure Physics, 2018, 32(6): 065102. doi: 10.11858/gywlxb.20180581

极小装药爆炸驱动两级重载荷研究

doi: 10.11858/gywlxb.20180581
详细信息
    作者简介:

    李鹏飞(1983-), 男, 硕士, 工程师, 主要从事战斗部设计及毁伤技术研究. E-mail:pengfeili@163.com

    通讯作者:

    宋浦(1973—),男,博士,研究员,主要从事战斗部毁伤技术研究. E-mail:songpu73@163.com

  • 中图分类号: O383

Study on Minimal Charge Explosive Driving Two Mass Load

  • 摘要: 根据瞬时爆轰理论,在极小装填比(C/M)条件下进行了小质量装药驱动大质量载荷运动研究,理论分析了弱约束和强约束装药条件下载荷驱动速度。弱约束装药条件下在格尼方程基础上通过考虑装药壳体径向飞散获得二维炸药驱动载荷运动的计算模型,强约束条件下按照爆轰产物等熵膨胀理论获得大质量载荷驱动速度的计算模型。开展了两种约束条件下爆炸驱动大质量载荷的运动试验研究,弱约束条件下,载荷驱动速度随着装药长径比减小而增大;强约束条件下,载荷驱动速度随着装药量增加而增大,随装药腔体与运动腔体直径比增大而增大。

     

  • 光滑粒子流体动力学(Smoothed Particle Hydrodynamics, SPH)方法中的搜索算法较耗时,即每个时间步都要对领域粒子进行搜索,粒子越多,耗时情况越突出,与有限元法相比,SPH方法的计算效率要低得多。为了解决SPH方法计算效率低的问题,Johnson等[1-2]和Attaway等[3]将有限元与SPH方法相结合,提出了SPH-FEM耦合的算法,即:在小变形区域使用有限元法,大变形区域仍使用SPH方法。该方法不仅提高了计算效率,而且适应性较强。

    目前,采用SPH方法对爆炸焊接进行数值模拟的相关报道较少,而且多采用二维SPH方法。Tanaka[4]采用SPH方法对爆炸焊接的斜碰撞过程进行了数值模拟,成功地模拟出射流、波形和涡旋,波长的模拟结果相对实验结果偏大。李晓杰等[5]采用SPH方法及热塑性流体力学模型对爆炸复合板的斜碰撞过程中出现的界面波进行了数值模拟,模拟结果与张登霞等[6-7]实验结果的一致性较好。刘江等[8]利用AUTODYN软件中的SPH方法模拟了爆炸复合的斜碰撞,结合模拟中有效塑性变形、温度及剪切应力呈现的变化规律发现,爆炸复合的结合机理集塑性变形、熔化和扩散为一体。本研究将采用三维SPH方法对双面爆炸焊接过程进行模拟,将其结果与实验及理论结果进行对比,分析SPH-FEM耦合方法对爆炸焊接模拟的有效性。

    以前期45钢/Q235钢双面爆炸焊接实验[9]为基础,考虑到计算效率,利用LS-DYNA建立如图 1图 2所示的两组双面爆炸焊接SPH-FEM耦合的三维真实计算模型,选用的炸药为乳化炸药(玻璃微球的质量分数为5%),计算模型中基板和复板的材料、尺寸、间隙(δ)及药厚如表 1所示。起爆方式为点起爆。

    图  1  计算模型Ⅰ(10 mm药厚)
    Figure  1.  Calculation model Ⅰ with explosive thickness of 10 mm
    图  2  计算模型Ⅱ(5 mm药厚)
    Figure  2.  Calculation model Ⅱ with explosive thickness of 5 mm
    表  1  计算模型中材料的相关参数
    Table  1.  Related parameters of materials in calculation models
    Calculationmodel Flyer plate Base plate Gap
    δ/mm
    Size of explosive/(mm×mm×mm)
    Material Size/(mm×mm×mm) Material Size/(mm×mm×mm)
    45 steel 300×150×2 Q235 300×150×16 6 300×150×10
    45 steel 300×150×2 Q235 300×150×16 6 300×150×5
    下载: 导出CSV 
    | 显示表格

    基、复板采用3D Solid 164实体单元,单元边长为0.1 cm;炸药划分为光滑粒子,粒子的大小Δr取为0.1 cm。考虑到模型的对称性,为了提高计算效率,采用1/2模型进行计算。单位制为cm-g-μs。

    数值计算中乳化炸药采用高能燃烧模型[10-11]及JWL状态方程[12]。JWL状态方程表达式为

    p=AJWL(1ωR1v)eR1v+BJWL(1ωR2v)eR2v+ωE0v (1)

    式中:AJWLBJWLR1R2ω为材料参数;p为爆轰产物压力,GPa;E0为初始比内能,kJ/cm3v为爆轰气体产物的相对比容,为无量纲量。炸药的相关参数见表 2,其中:ρ为密度,D为炸药爆速。

    表  2  乳化炸药的JWL状态参数[13]
    Table  2.  JWL equation-of-state parameters of emulsion explosives[13]
    ρ/(g·cm-3) D/(m·s-1) AJWL/GPa BJWL/GPa R1 R2 ω E0/(kJ·cm-3)
    1.12 4 510 326.42 5.808 9 5.80 1.56 0.57 3.323
    下载: 导出CSV 
    | 显示表格

    数值计算中,基、复板均采用Mie-Grüneisen状态方程[14]和Johnson-Cook材料模型[15]。Johnson-Cook材料模型的形式如下

    σ=(A+Bεnp)(1+Cln˙εp)(1Tm) (2)

    式中:εp为有效塑性应变;˙εp=˙εp/˙ε0p为有效塑性应变率,其中˙ε0p为参考应变率;ABCmn为与材料相关的常数;无量纲温度T*表示为T*=(T-Tr)/(Tm-Tr),其中Tr为室温, Tm为熔点。45钢选用与Q235钢相同的Johnson-Cook材料模型参数,具体参数如表 3所示。

    表  3  Q235钢的Johnson-Cook模型参数[16]
    Table  3.  Johnson-Cook parameters of Q235 steel[16]
    ρ/(g·cm-3) G/GPa A/GPa B/GPa C n m Tm/K Tr/K
    7.83 77 0.792 0.51 0.014 0.26 1.03 1 793 294
    下载: 导出CSV 
    | 显示表格
    2.1.1   碰撞点位移

    图 3所示是爆炸焊接结束时复板的竖向位移云图。由图 3可看出,复板的位移大致均为6 mm,表明基、复板已完全复合。为了更直观地观察复板单元位移的变化情况,在复板上选择3个特征单元(431 806、437 359、444 788),输出其位移-时间曲线,如图 4所示。由图 4可看出,特征单元的竖向位移均略大于间隙(6 mm),这是由于在爆炸载荷作用下复板有一定程度的减薄率所致。

    图  3  10 mm药厚下爆炸焊接结束时复板的z向位移云图
    Figure  3.  z-direction displacement contour of flyer plate with explosive thickness of 10 mm at the end of explosive welding
    图  4  10 mm药厚下复板上3个特征单元的z向位移-时间历程
    Figure  4.  z-direction displacement histories of 3 characteristic elements with explosive thickness of 10 mm
    2.1.2   复板碰撞速度

    图 5所示是一对分别取自基板与复板结合界面处的特征单元(基板单元:798 751;复板单元:416 251),特征单元的选取与前期实验[9]中金相试样的取样位置一致。

    图  5  10 mm药厚下的一对特征单元
    Figure  5.  A pair of characteristic elements with explosive thickness of 10 mm

    图 6所示是这对特征单元的速度-时间曲线。可以看出,基板在碰撞前有一个正的速度峰;该现象的产生如文献[17]所述,是由于爆轰产物不断堆积以及前碰撞点在待复合区产生的振动能所致。复板上所取单元的最大碰撞速度为897 m/s。

    图  6  一对特征单元(见图 5)的速度-时间曲线
    Figure  6.  Velocity-time curves of the pair of characteristic elements (see Fig. 5)

    图 7所示是在复板结合界面处所选取的3个特征单元(410 476、416 251、420 976)。图 8所示是这3个特征单元的速度-时间曲线。

    图  7  10 mm药厚下复板结合界面处的3个特征单元
    Figure  7.  3 characteristic elements at the bonding interface of flyer plate with explosive thickness of 10 mm
    图  8  10 mm药厚下3个特征单元的速度-时间曲线
    Figure  8.  Velocity-time curves of 3 characteristic elements with explosive thickness of 10 mm

    图 8可以看出,随着距起爆端距离的增加,复板的碰撞速度增大。由文献[17]的结论可知,该现象是由于基板与复板的碰撞在金属板的待复合区产生了强烈振动引起的。

    2.1.3   碰撞点压力分布

    图 9所示是在结合界面处选取的3个特征单元(415 576、418 051、419 776),单元415 576取在复板中心处,与前期实验[7]中取样做金相观察的位置一致。图 10所示是3个特征单元的压力历程。

    图  9  10 mm药厚下复板结合界面处的3个特征单元
    Figure  9.  3 characteristic elements at the bonding interface of flyer plate with explosive thickness of 10 mm
    图  10  10 mm药厚下3个特征单元的压力-时间曲线
    Figure  10.  Pressure-time curves of 3 characteristic elements with explosive thickness of 10 mm

    图 10可以看出,随着距起爆端距离的增加,复板的碰撞压力增大。由文献[17]的结论可知,该现象是爆轰产物不断堆积以及前碰撞点在金属板待复合区振动能不断增加的共同作用结果。

    2.2.1   碰撞点位移

    图 11所示是爆炸焊接结束时复板的竖向位移云图。由图 11可看出,复板的位移大致均为6 mm,表明基、复板已完全复合。为了更加直观地观察复板单元位移的变化情况,在复板上选择3个特征单元(432 182、438 034、443 960),输出其位移-时间曲线,如图 12所示。由图 12可看出,特征单元的竖向位移均略大于6 mm,但较10 mm药厚下的竖向位移小。这是由于5 mm药厚下的爆炸载荷作用比10 mm药厚下小,导致5 mm药厚下的复板减薄率比10 mm药厚下低。

    图  11  5 mm药厚下爆炸复合结束时复板的z向位移云图
    Figure  11.  z-displacement contour of flyer plate with explosive thickness of 5 mm at the end of explosive welding
    图  12  5 mm药厚下复板上特征单元的z向位移-时间历程
    Figure  12.  z-displacement histories of 3 characteristic elements with explosive thickness of 5 mm
    2.2.2   复板碰撞速度

    图 13所示是一对分别取自基板与复板结合界面处的特征单元(基板单元:799 201;复板单元:416 701),特征单元的选取与前期实验[9]中金相试样的取样位置一致。

    图  13  5 mm药厚下的一对特征单元
    Figure  13.  A pair of characteristic elements with explosive thickness of 5 mm

    图 14所示是这对特征单元的速度-时间曲线, 可以看出,基板在碰撞前也有一个正的速度峰。复板上所取单元的最大碰撞速度为565 m/s。

    图  14  一对特征单元(见图 13)的速度-时间曲线
    Figure  14.  Velocity-time curves of the pair of characteristic elements (see Fig. 13)

    图 15所示是在复板结合界面处所选取的3个特征单元(411 976、417 001、423 826)。图 16所示是这3个特征单元的速度-时间曲线。由图 16可以看出,随着距起爆端距离的增加,复板的碰撞速度增大。

    图  15  5 mm药厚下复板结合界面处的3个特征单元
    Figure  15.  3 characteristic elements at the bonding interface of flyer plate with explosive thickness of 5 mm
    图  16  5 mm药厚下3个特征单元的速度-时间曲线
    Figure  16.  Velocity-time curves of 3 characteristic elements with explosive thickness of 5 mm
    2.2.3   碰撞点压力分布

    图 17所示是在结合界面处选取的3个特征单元(416 326、418 801、422 776),单元416 326取在复板中心处,与前期实验[7]中取样做金相观察的位置一致。图 18所示是这3个特征单元的压力历程。由图 18可以看出,随着距起爆端距离的增加,复板的碰撞压力增大。

    图  17  5 mm药厚下复板结合界面处的3个特征单元
    Figure  17.  3 characteristic elements at the bonding interface of flyer plate with explosive thickness of 5 mm
    图  18  5 mm药厚下3个特征单元的压力历程
    Figure  18.  Pressure histories of 3 characteristic elements with explosive thickness of 5 mm

    图 6可以看出,10 mm药厚下复板的最大碰撞速度为897 m/s。由图 14可以看出,5 mm药厚下复板的最大碰撞速度为565 m/s。利用前期工作[18]中提到的3种理论公式(Gurney公式、Aziz公式、Deribas公式)计算了复板的碰撞速度,如表 4表 5所示,并与数值模拟结果进行了比较。由表 4表 5可以看出:Gurney公式和Aziz公式的计算结果均存在较大的偏差;而由Deribas公式计算的两组结果与数值模拟结果较接近,误差均未超过5%,且与前期实验结果较吻合,证明了SPH-FEM耦合算法的可靠性。

    表  4  10 mm药厚下碰撞速度理论计算结果与数值模拟结果的比较
    Table  4.  Comparison of collision velocity between theoretical calculation and numerical simulation with explosive thickness of 10 mm
    Theoreticalformula Massfraction Collision velocity/(m·s-1) Error/%
    Theoretical calculation[18] Simulation
    Gurney 0.75 1 089 897 -21.0
    Aziz 0.75 711 897 20.0
    Deribas 0.75 853 897 4.9
    下载: 导出CSV 
    | 显示表格
    表  5  5 mm药厚下碰撞速度理论计算结果与数值模拟结果的比较
    Table  5.  Comparison of collision velocity between theoretical calculation and numerical simulation with explosive thickness of 5 mm
    Theoreticalformula Massfraction Collision velocity/(m·s-1) Error/%
    Theoretical calculation[18] Simulation
    Gurney 0.45 863 565 -52.7
    Aziz 0.45 480 565 15.0
    Deribas 0.45 576 565 -1.9
    下载: 导出CSV 
    | 显示表格

    图 10可以看出,10 mm药厚下复板单元415 576处的碰撞压力为17.08 GPa。由图 18可以看出,5 mm药厚下复板单元416 326处的碰撞压力为11.25 GPa。

    Ezra等提出的碰撞压力的计算公式为[16]

    p=ρ1vs,1vp1+ρ1vs,1ρ2vs, 2 (3)

    式中:vs, 1vs, 2分别表示复板、基板的声速,m·s-1ρ1ρ2分别表示复板、基板的密度,g·cm-3vp表示复板的碰撞速度,m·s-1

    结合表 4表 5中3种理论公式计算得到的碰撞速度,通过(3)式可得到复板的碰撞压力,表 6表 7为其理论计算值与数值模拟结果的比较。可见:Gurney公式和Aziz公式的计算结果均存在较大的偏差;而由Deribas公式计算的两组结果与数值模拟结果较接近,误差均未超过5%,说明Deribas公式和SPH-FEM耦合方法对双面爆炸焊接具有较好的指导意义。

    表  6  10 mm药厚下碰撞压力理论计算结果与数值模拟结果的比较
    Table  6.  Comparison of collision pressure betweentheoretical calculation and numerical simulationwith explosive thickness of 10 mm
    Theoreticalformula Collision pressure/GPa Error/%
    Calculation Simulation
    Gurney 22.08 17.08 -29.3
    Aziz 14.42 17.08 15.6
    Deribas 17.30 17.08 -1.3
    下载: 导出CSV 
    | 显示表格
    表  7  5 mm药厚下碰撞压力理论计算结果与数值模拟结果的比较
    Table  7.  Comparison of collision pressure betweentheoretical calculation and numerical simulationwith explosive thickness of 5 mm
    Theoreticalformula Collision pressure/GPa Error/%
    Calculation Simulation
    Gurney 17.50 11.25 -55.6
    Aziz 9.73 11.25 13.5
    Deribas 11.68 11.25 -3.8
    下载: 导出CSV 
    | 显示表格

    利用LS-DYNA软件和SPH-FEM耦合方法对前期双面爆炸焊接实验进行了三维数值模拟,并将模拟结果与实验及理论计算结果进行了对比,得到如下结论。

    (1) 10 mm药厚和5 mm药厚下复板位移均略大于间隙值6 mm,这是由于爆轰载荷作用下复板有一定的减薄率所致。

    (2) 10 mm药厚下,复板中部的最大碰撞速度为897 m/s,碰撞压力为17.08 GPa;5 mm药厚下,复板中部的最大碰撞速度为565 m/s,碰撞压力为11.25 GPa。通过与3种理论公式(Gurney公式、Aziz公式、Deribas公式)计算得到的碰撞速度进行比较发现,数值模拟结果与Deribas公式的计算结果较接近,误差较小,且与实验结果较吻合,证明了SPH-FEM耦合方法用于双面爆炸复合模拟的有效性,同时Deribas公式和SPH-FEM耦合方法对双面爆炸复合具有较好的指导意义。

    (3) 10 mm药厚和5 mm药厚下复板的碰撞速度及碰撞压力均随着距起爆端距离的增加而增大,该现象是由于爆轰产物的不断堆积和前碰撞点在金属板待复合区振动能的不断增加共同作用的结果。

  • 图  弱约束条件下极小装药驱动大质量载荷模型

    Figure  1.  Schematic illustration of a small charge-driven mass load model with weak constraints

    图  强约束条件下极小装药驱动大质量载荷模型

    Figure  2.  Schematic illustration of a small charge-driven mass load model with strong constraints

    图  试验布局图

    Figure  3.  Experiment layout

    图  载荷驱动高速摄影图

    Figure  4.  Load driven high speed photogram

    图  弱约束条件下载荷运动速度

    Figure  5.  Load velocity under the weak constraint

    图  装药质量对速度的影响

    Figure  6.  Effect of charge mass on speed

    图  装药和腔体尺寸对运动速度的影响

    Figure  7.  Effect of charge and cavity size on the velocity of movement

  • [1] 孙承纬, 文尚刚, 赵峰.多级炸药爆轰高速驱动技术的Gurney模型优化分析[J].爆炸与冲击, 2004, 24(4):299-303. doi: 10.3321/j.issn:1001-1455.2004.04.002

    SUN C W, WEN S G, ZHAO F.An optimal analysis of multi-stage explosive accelerated high velocity flyers with the improved Gurney model[J]. Explosion and Shock Waves, 2004, 24(4):299-303. doi: 10.3321/j.issn:1001-1455.2004.04.002
    [2] 雷灏, 姚志敏, 尉广军.基于Gurney公式的半预制破片飞散特性研究[J].爆破, 2014, 31(1):38-41. http://d.old.wanfangdata.com.cn/Periodical/bp201401008

    LEI H, YAO Z M, YU G J.Research on dispersion characteristics of pre-formed fragments based on gurney's formula[J]. Blasting, 2014, 31(1):38-41. http://d.old.wanfangdata.com.cn/Periodical/bp201401008
    [3] 李定和.云爆弹云雾引信抛射轨迹的研究[J].兵工学报, 1993, 14(2):74-77. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000001078525

    LI D H.A study on the trajectory of secondary fuze (for initiation of the fillermist) in a fuel-air explosive munition[J]. Acta Armamentarii, 1993, 14(2):74-77. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000001078525
    [4] 何勇, 梁安定, 孙兴昀.一种新型防空反导杀伤增强装置驱动技术研究[J].科学技术与工程, 2016, 16(18):192-195. doi: 10.3969/j.issn.1671-1815.2016.18.036

    HE Y, LIANG A D, SUN X Y.Research on a new drive technology of kill-enhanced device of air defense and anti-missile[J]. Science Technology and Engineering, 2016, 16(18):192-195. doi: 10.3969/j.issn.1671-1815.2016.18.036
    [5] 隋树元, 王树山.终点效应学[M].北京:国防工业出版社, 2000:75-82.

    SUI S Y, WANG S S.Terminal effects[M]. Beijing:National Defense Industry Press, 2000:75-82.
    [6] YADAV H S.Flyer plate motion by thin sheet of explosive[J]. Propellants Explosives Pyrotechnics, 1988, 13(1):17-20. doi: 10.1002/(ISSN)1521-4087
    [7] 孙承纬, 卫玉章, 周之奎.应用爆轰物理[M].北京:国防工业出版社, 2000:273-276.

    SUN C W, WEI Y Z, ZHOU Z K.Application of detonation physics[M]. Beijing:National Defense Industry Press, 2000:273-276.
    [8] 沈飞, 王辉, 袁建飞, 等.炸药格尼系数的一种简易估算法[J].火炸药学报, 2013, 36(6):36-38. doi: 10.3969/j.issn.1007-7812.2013.06.007

    SHEN F, WANG H, YUAN J F, et al.A simple estimation method for gurney coefficient of explosive[J]. Chinese Journal of Explosives & Propellants, 2013, 36(6):36-38. doi: 10.3969/j.issn.1007-7812.2013.06.007
  • 加载中
图(7)
计量
  • 文章访问数:  6424
  • HTML全文浏览量:  3026
  • PDF下载量:  18
出版历程
  • 收稿日期:  2018-06-13
  • 修回日期:  2018-07-20

目录

/

返回文章
返回