高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

炸药覆盖层对爆炸焊接影响的数值模拟

缪广红 胡昱 艾九英 祁俊翔 马宏昊 沈兆武

缪广红, 胡昱, 艾九英, 祁俊翔, 马宏昊, 沈兆武. 炸药覆盖层对爆炸焊接影响的数值模拟[J]. 焊接学报, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002
引用本文: 缪广红, 胡昱, 艾九英, 祁俊翔, 马宏昊, 沈兆武. 炸药覆盖层对爆炸焊接影响的数值模拟[J]. 焊接学报, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002
MIAO Guanghong, HU Yu, AI Jiuying, QI Junxiang, MA Honghao, SHEN Zhaowu. Numerical simulation research on the effect of explosive covering on explosive welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002
Citation: MIAO Guanghong, HU Yu, AI Jiuying, QI Junxiang, MA Honghao, SHEN Zhaowu. Numerical simulation research on the effect of explosive covering on explosive welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002

炸药覆盖层对爆炸焊接影响的数值模拟

doi: 10.12073/j.hjxb.20220121002
基金项目: 国家自然科学基金(11902003, 51874267);安徽省高校自然科学基金重点项目(KJ2017A089);高校优秀青年骨干人才国外访学研修项目(gxgwfx2019017).
详细信息
    作者简介:

    缪广红,1985年出生,博士,副教授,硕士生导师;主要从事含能材料、爆炸复合及爆炸安全等项目研究. Email: miaogh@mail.ustc.edu.cn

  • 中图分类号: TG 456.6

Numerical simulation research on the effect of explosive covering on explosive welding

  • 摘要: 为探究炸药覆盖层厚度对爆炸焊接的影响,采用ANSYS/LS-DYNA软件并结合SPH-FEM耦合算法,对不同覆层厚度下的爆炸焊接试验进行三维数值模拟. 文中采用厚度为 20 mm 的Q235钢和厚度为 2.5 mm 的304不锈钢作为基板和复板. 根据相应的材料参数理论计算了焊接过程中的动态参数,并以此建立爆炸焊接窗口. 仿真结果表明,与无覆盖层爆炸焊接相比,覆盖层厚度为15 mm、 30 mm 和45 mm 时冲击速度分别提高了39.3%, 58.1%和68.8%,碰撞压力分别增大了41.0%, 65.6% 和80.6%. 仿真结果与试验结果基本一致. 利用SPH法进行二维数值模拟,得到了装配炸药覆盖层时复板与基板的复合界面. 仿真结果表明,复合板在覆层厚度为15 mm时具有良好的波形复合界面,且界面波形与试验金相分析结果较为吻合.
  • 图  1  计算模型示意图

    Figure  1.  Schematic diagram of calculation model

    图  2  爆炸焊接窗口

    Figure  2.  Explosive welding window.

    图  3  结合面上的一对特征单元

    Figure  3.  A pair of element on the joint surface.

    图  4  不同覆层厚度下特征单元A与B的速度-时间历程

    Figure  4.  Velocity-time history of element A and B under different covering thicknesses. (a) Colloidal water covering thickness 0 mm; (b) Colloidal water covering thickness 15 mm; (c) Colloidal water covering thickness 30 mm; (d) Colloidal water covering thickness 40 mm

    图  5  碰撞速度随覆层厚度变化的增益曲线

    Figure  5.  Gain curve of impact velocity with variation of covering thickness

    图  6  基板与覆板的碰撞角β

    Figure  6.  Collision angle β of base plate and flyer plate

    图  7  覆层厚度15 mm时复板的压力云图

    Figure  7.  The pressure cloud diagram of the flyer plate when covering thickness is 15 mm. (a) 4.99 μs; (b) 18.99 μs; (c) 25.96 μs; (d) 46.99 μs; (e) 74.99 μs; (f) 90.98 μs

    图  8  15 mm覆层下复板界面结合处的三个特征单元

    Figure  8.  Three element at the junction of the flyer plate under the 15 mm covering

    图  9  三个特征单元的碰撞压力-时间分布图

    Figure  9.  Collision pressure-time distribution diagram of three element

    图  10  不同覆层厚度下特征单元A (元素17911)的碰撞压力分布

    Figure  10.  Impact pressure distribution of element A (element 17911) under different covering thicknesses. (a) Covering thickness 0 mm; (b) Covering thickness 15 mm; (c) Covering thickness 30 mm; (d) Covering thickness 45 mm

    图  11  不同覆层厚度下的位移云图和位移-时间历程图

    Figure  11.  Displacement cloud diagram and displacement-time history diagram under different covering thickness. (a) covering thickness 0 mm; (b) covering thickness 15 mm; (c) covering thickness 30 mm; (d) covering thickness 45 mm

    图  12  覆层厚度15mm复合界面波形对比图

    Figure  12.  Waveform comparison diagram on the composite interface when the covering thickness is 15 mm. (a) metallographic analysis composite interface; (b) 2D simulation composite interface

    表  1  乳化炸药的JWL状态方程参数

    Table  1.   JWL equation-of-state parameters of emulsion explosive

    炸药爆速
    D/(m∙s−1)
    炸药密度
    ρ/(g·cm−3)
    单位体积内能
    E0 /(kJ∙cm−3)
    材料常数1
    AJWL /GPa
    材料常数2
    BJWL/GPa
    材料常数3
    R1
    材料常数4
    R2
    材料常数5
    ω
    3 0271.013.323326.425.808 95.81.560.57
    下载: 导出CSV

    表  2  Q235钢与304不锈钢的Johnson-Cook材料模型参数

    Table  2.   Parameters of Johnson-Cook mode1l of Q235 steel and SUS304 steel

    材料密度
    ρ/(g·cm−3)
    剪切模量
    G/GPa
    初始屈服
    强度A/GPa
    硬度常数
    B/GPa
    硬化指数
    n
    应变率强化
    参数c
    软化指数
    m
    室内温度
    Tr /K
    金属熔点
    Tm /K
    Q2357.8577.00.7920.510.260.0141.03294.01 493
    SUS3047.9324.00.7001.300.750.0210.90294.01 454
    下载: 导出CSV

    表  3  Q235钢与304不锈钢的Gruneisen状态方程参数

    Table  3.   Gruneisen EOS parameters of Q235 steel and SUS304 steel

    材料体积声速
    C/(km·s−1)
    斜率系数
    S
    Gruneisen系数
    γ0
    体积修正系数
    a
    Q2356.01.492.170.46
    SUS3044.51.491.930.50
    下载: 导出CSV

    表  4  胶体水覆层的材料模型与状态方程参数

    Table  4.   Model and EOS parameters of Colloidal water

    密度
    ρ/(g·cm−3)
    截止压力
    PC/Pa
    动态粘度系数
    MU/(10−4 N·s·m−2)
    体积声速
    C/(km·s−1)
    斜率
    系数S
    Gruneisen
    系数γ0
    0.998−10.08.6840.164 71.9210.35
    下载: 导出CSV

    表  5  Q235钢与304不锈钢的材料性能

    Table  5.   The material properties of Q235 steel and SUS304 steel

    材料密度
    ρ/(g·cm−3)
    材料声速
    C0/(km·s−1)
    拉伸强度
    σb/GPa
    维氏硬度
    HV/GPa
    材料熔点
    Tm/℃
    热导率
    κ/(W·m−1·℃−1)
    Q2357.856.000.4051.301 49338
    SUS3047.934.500.5601.701 45422
    下载: 导出CSV

    表  6  不同覆层厚度下的炸药爆速计算值

    Table  6.   Calculation value of explosive detonation velocity under different cladding thickness

    覆层厚度h/mm炸药爆速Vd/(m·s−1)
    02 950
    153 398
    303 530
    453 582
    下载: 导出CSV
  • [1] Fronczek D M, Wojewoda-Budka J, Chulist R, et al. Structural properties of Ti/Al clads manufactured by explosive welding and annealing[J]. Materials & Design, 2016, 91: 80 − 89.
    [2] Mahmood Y, Dai K D, Chen P W, et al. Experimental and Numerical Study on Microstructure and Mechanical Properties of Ti-6Al-4V/Al-1060 Explosive Welding[J]. Metals, 2019, 9(11): 1189. doi:  10.3390/met9111189
    [3] Ren B X, Tao G, Wen P, et al. Study on weldability window and interface morphology of steel tube and tungsten alloy rod welded by explosive welding[J]. International Journal of Refractory Metals and Hard Materials, 2019, 84: 105005. doi:  10.1016/j.ijrmhm.2019.105005
    [4] Liu K Y, Chen P W, Feng J R, et al. Fabrication and characterization of the Mo/CU bimetal with thick Mo layer and high interfacial strength[J]. International Journal of Refractory Metals and Hard Materials, 2021, 94: 105383. doi:  10.1016/j.ijrmhm.2020.105383
    [5] Akbari M S A A, Farhadi S P. Experimental investigation of explosive welding of cp-titanium/AISI 304 stainless steel[J]. Materials & Design, 2009, 30(3): 459 − 468.
    [6] Gladkovsky S V, Kuteneva S V, Sergeev S N. Microstructure and mechanical properties of sandwich copper/steel composites produced by explosive welding[J]. Materials Characterization, 2019, 154: 294 − 304. doi:  10.1016/j.matchar.2019.06.008
    [7] Loureiro A, Mendes R, Ribeiro J B, et al. Effect of explosive mixture on quality of explosive welds of copper to aluminium[J]. Materials & Design, 2016, 95: 256 − 267.
    [8] 邢廷勇, 郑雁军, 崔立山, 等. 时效对爆炸焊接合成TiNi/TiNi合金阻尼行为的影响[J]. 中国有色金属学会学报, 2009, 19(6): 1470 − 1473. doi:  10.1016/S1003-6326(09)60053-4

    Xing Tingyong, Zheng Yanjun, Cui Lishan, et al. Influence of aging on damping behavior of TiNi/TiNi alloys synthesized by explosive welding[J]. Transactions of Nonferrous Metals Society of China, 2009, 19(6): 1470 − 1473. doi:  10.1016/S1003-6326(09)60053-4
    [9] 吴晓明, 史长根, 房中行, 等. 有无夹层的钛铝爆炸焊接能量与界面特性的对比研究[J]. 材料与设计, 2021, 197: 109279.

    Wu Xiaoming, Shi Changgen, Fang Zhonghang, et al. Comparative study on welding energy and Interface characteristics of titanium aluminum explosive composites with and without interlayer[J]. Materials & Design, 2021, 197: 109279.
    [10] Mori D, Kasada R, Konishi S, et al. Underwater explosive welding of tungsten to reduced-activation ferritic steel F82H[J]. Fusion Engineering and Design, 2014, 89(7-8): 1086 − 1090. doi:  10.1016/j.fusengdes.2013.12.038
    [11] Liu W D, Liu K X, Chen Q Y, et al. Metallic glass coating on metals plate by adjusted explosive welding technique[J]. Applied Surface Science, 2009, 255(23): 9343 − 9347. doi:  10.1016/j.apsusc.2009.07.033
    [12] Sui G F, Li J S, Li H W, et al. Investigation on the explosive welding mechanism of corrosion-resisting aluminum and stainless steel tubes through finite element simulation and experiments[J]. International Journal of Minerals Metallurgy and Materials, 2012, 19(2): 151 − 158. doi:  10.1007/s12613-012-0531-4
    [13] 杨明, 马宏昊, 沈兆武, 等. 炸药覆层对不锈钢与钢爆炸焊接的影响研究[J]. 推进剂, 炸药, 烟火, 2019, 44(5): 609 − 616. doi:  10.1002/prep.201800160

    Yang Ming, Ma Hong-hao, Shen Zhao-wu,. et al. Study on the Effects of Explosive Covering on Explosive Welding of Stainless Steel to Steel[J]. Propellants, Explosives, Pyrotechnics, 2019, 44(5): 609 − 616. doi:  10.1002/prep.201800160
    [14] 章冠人, 陈大年. 凝聚炸药起爆动力学[M]. 北京: 国防工业出版社, 1991.

    Zhang Guanren, Chen Danian. Condensed Explosive Detonation Dynamics[M]. Beijing: National Defense Industry Press, 1991.
    [15] 李裕春, 时党勇, 赵远. ANSYS11.0/LS-DYNA基础理论与工程实践[M]. 北京: 中国水利水电出版社, 2008.

    Li Yuchun, Shi Dangyong, Zhao Yuan. ANSYS11.0/LS-DYNA Basic Theory and Engineering Practice[M]. Beijing: China Water&Power Press, 2008.
    [16] 辛春亮, 薛再清, 涂建, 等. 有限元分析常用材料参数手册[M]. 北京: 机械工业出版社, 2019.

    Xin Chunliang, Xue Zaiqing, Tu Jian, et al. Handbook of Common Material Parameters for Finite Element Analysis. Beijing: China Machine Press, 2019.
    [17] Steinberg D J, Cochran S G, Guinan M W. A constitutive model for metals applicable at high-strain rate[J]. Journal of Applied Physics, 1980, 51: 1498 − 1504. doi:  10.1063/1.327799
    [18] Liu G R, Liu M B. Smoothed Particle Hydrodynamics: A Meshfree Particle Method[M]. Changsha: Hunan University Press, 2005.
    [19] Walsh J, Shreffler R, Willig F. Limiting conditions for jet formation in high velocity collisions[J], Journal of Applied Physics, 1953, 24: 349359.
    [20] 郑远谋. 爆炸焊接和爆炸复合材料[M]. 北京: 国防工业出版社, 2017.

    Zheng Yuanmou. Explosive Welding and Explosive Composite Material[M]. Beijing: National Defense Industry Press, 2017.
    [21] Cowan G R, Bergmann O R, Holtzman A H. Mechanism of bond zone wave formation in explosion-clad metals[J]. Metallurgical and Materials Transactions B, 1971, 2: 3145 − 3155. doi:  10.1007/BF02814967
    [22] Shao B H, Zhang K. Explosive welding principle and its application[M]. Dalian: Dalian University of Science and Technology Press, 1987.
    [23] 缪广红, 马雷鸣, 吴建强, 等. 基复板间距对爆炸焊接质量影响的数值模拟[J]. 爆破, 2020, 37(2): 106 − 114. doi:  10.3963/j.issn.1001-487X.2020.02.019

    Miao Guanghong, Ma Leiming, Wu Jianqiang, et al. Numerical Simulation of Influence of Spacing between Base Plate and Flyer Plate on Welding Quality[J]. Blasting, 2020, 37(2): 106 − 114. doi:  10.3963/j.issn.1001-487X.2020.02.019
    [24] Zheng Z M, Yang Z S. Explosion Processing[M]. Beijing: National Defense Industry Press, 1981.
  • [1] 李雪交, 钱敬业, 毕志雄, 张廷赵, 代弦德, 荣凯.  基于自约束结构炸药的T2/Q345爆炸焊接及数值模拟 . 焊接学报, 2023, 44(3): 70-76. doi: 10.12073/j.hjxb.20220419002
    [2] 周进鹏, 马杰, 陆晓峰, 朱晓磊, 王健.  喷雾辅助FSW焊接RAFM钢数值模拟与性能 . 焊接学报, 2022, 43(9): 104-112. doi: 10.12073/j.hjxb.20210922001
    [3] 缪广红, 胡昱, 艾九英, 马秋月, 孙志皓, 马宏昊, 沈兆武.  基于不同算法的金属管棒爆炸焊接模拟 . 焊接学报, 2022, 43(12): 64-71. doi: 10.12073/j.hjxb.20211219002
    [4] 缪广红, 艾九英, 胡昱, 马宏昊, 沈兆武.  基于SPH法的爆炸焊接边界效应二维数值模拟 . 焊接学报, 2021, 42(9): 61-66. doi: 10.12073/j.hjxb.20210203002
    [5] 孙加民, 蔡建鹏, 叶延洪, 邓德安.  电渣焊接头温度场数值模拟 . 焊接学报, 2015, 36(7): 93-96.
    [6] 周广涛, 郭广磊, 方洪渊.  激光诱导焊接温度场的数值模拟 . 焊接学报, 2014, 35(7): 22-26.
    [7] 邓伟, 陆明, 田晓洁.  铝/钛复合管爆炸焊接三维数值模拟 . 焊接学报, 2014, 35(12): 63-66.
    [8] 闫德俊, 刘雪松, 杨建国, 方洪渊, 杨永强.  焊接过程数值模拟材料参量的确定 . 焊接学报, 2012, (9): 24-28.
    [9] 史长根, 汪育, 徐宏.  双立爆炸焊接及防护装置数值模拟和试验 . 焊接学报, 2012, (3): 109-112.
    [10] 邓德安, 童彦刚, 周中玉.  薄壁低碳钢管焊接变形的数值模拟 . 焊接学报, 2011, (2): 81-84.
    [11] 王呼和, 佟铮.  平板金属爆炸焊接过程数值模拟 . 焊接学报, 2010, (9): 101-104,108.
    [12] 张增磊, 史清宇, 刘园, 鄢东洋.  焊接数值模拟材料新模型的建立及应用 . 焊接学报, 2009, (2): 45-48.
    [13] 马贝, 李宏伟, 常辉, 胡锐.  间隙对三层圆管爆炸焊接影响的数值模拟 . 焊接学报, 2009, (9): 33-36,40.
    [14] 张华军, 张广军, 蔡春波, 王俊恒, 吴林.  摆动焊接动态过程温度场数值模拟 . 焊接学报, 2008, (2): 69-72,76.
    [15] 王建民, 朱锡, 刘润泉.  爆炸焊接三维数值模拟 . 焊接学报, 2007, (5): 109-112.
    [16] 赵衍华, 林三宝, 贺紫秋, 吴林.  三维搅拌摩擦焊接过程数值模拟 . 焊接学报, 2006, (6): 75-78.
    [17] 谢飞鸿, 罗冠炜, 廖军生, 汪旭光.  金属爆炸焊接界面应力场数值计算分析 . 焊接学报, 2006, (10): 26-28,93.
    [18] 王希靖, 韩晓辉, 郭瑞杰, 李晶.  搅拌摩擦焊接过程温度场数值模拟 . 焊接学报, 2005, (12): 17-20.
    [19] 曹健, 冯吉才, 李卓然.  Ti/TiAl钎焊接头应力的数值模拟 . 焊接学报, 2005, (1): 5-8.
    [20] 芦凤桂, 唐新华, 李少青, 姚舜, 楼松年.  定点TIG焊接电弧与熔池交互耦合数值模拟 . 焊接学报, 2005, (9): 69-72,76.
  • 加载中
图(12) / 表 (6)
计量
  • 文章访问数:  68
  • HTML全文浏览量:  14
  • PDF下载量:  16
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-01-21
  • 网络出版日期:  2023-02-10
  • 刊出日期:  2023-01-25

炸药覆盖层对爆炸焊接影响的数值模拟

doi: 10.12073/j.hjxb.20220121002
    基金项目:  国家自然科学基金(11902003, 51874267);安徽省高校自然科学基金重点项目(KJ2017A089);高校优秀青年骨干人才国外访学研修项目(gxgwfx2019017).
    作者简介:

    缪广红,1985年出生,博士,副教授,硕士生导师;主要从事含能材料、爆炸复合及爆炸安全等项目研究. Email: miaogh@mail.ustc.edu.cn

  • 中图分类号: TG 456.6

摘要: 为探究炸药覆盖层厚度对爆炸焊接的影响,采用ANSYS/LS-DYNA软件并结合SPH-FEM耦合算法,对不同覆层厚度下的爆炸焊接试验进行三维数值模拟. 文中采用厚度为 20 mm 的Q235钢和厚度为 2.5 mm 的304不锈钢作为基板和复板. 根据相应的材料参数理论计算了焊接过程中的动态参数,并以此建立爆炸焊接窗口. 仿真结果表明,与无覆盖层爆炸焊接相比,覆盖层厚度为15 mm、 30 mm 和45 mm 时冲击速度分别提高了39.3%, 58.1%和68.8%,碰撞压力分别增大了41.0%, 65.6% 和80.6%. 仿真结果与试验结果基本一致. 利用SPH法进行二维数值模拟,得到了装配炸药覆盖层时复板与基板的复合界面. 仿真结果表明,复合板在覆层厚度为15 mm时具有良好的波形复合界面,且界面波形与试验金相分析结果较为吻合.

English Abstract

缪广红, 胡昱, 艾九英, 祁俊翔, 马宏昊, 沈兆武. 炸药覆盖层对爆炸焊接影响的数值模拟[J]. 焊接学报, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002
引用本文: 缪广红, 胡昱, 艾九英, 祁俊翔, 马宏昊, 沈兆武. 炸药覆盖层对爆炸焊接影响的数值模拟[J]. 焊接学报, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002
MIAO Guanghong, HU Yu, AI Jiuying, QI Junxiang, MA Honghao, SHEN Zhaowu. Numerical simulation research on the effect of explosive covering on explosive welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002
Citation: MIAO Guanghong, HU Yu, AI Jiuying, QI Junxiang, MA Honghao, SHEN Zhaowu. Numerical simulation research on the effect of explosive covering on explosive welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(1): 40-48. doi: 10.12073/j.hjxb.20220121002
    • 复合材料的概念最早源于对贝类生物的研究,基于贝壳层次结构的复合金属材料具有单金属材料所不具备的特殊力学性能[1]. 连接两种或多种材料的方法有很多,常见的方法有热压、铸造、电镀、焊接等. 爆炸焊接是一种使用炸药爆炸生成的冲击力加速复层材料,使复层材料在短时间内与基层材料高速碰撞产生高温高压的焊接方法[2]. 爆炸焊接具有其他复合技术无法比拟的优点:当待焊接材料的碰撞速度被控制在600 ~ 3 000 m/s 之间时,碰撞产生的高压能够使绝大部分材料产生塑性变形[3],再结合高达107 K/s 的冷却速率,有效杜绝了结合界面上生成意外相[4]. 这些优点使爆炸焊接成为了一种优秀的金属复合方法,应用前景广阔[5].

      一直以来,关于爆炸焊接的研究主要集中于测试复合材料的物理与金属性能和炸药配比对焊接质量的影响等方面. Gladkovsky等人[6]对铜与低碳钢复合材料的界面组织和力学性能特征进行了研究;Loureiro等人[7]讨论了炸药爆炸比和敏化剂类型对铜铝板爆炸焊缝质量的影响;刑廷勇等人[8]利用爆炸焊接成功制备了2种不同组分比的TiNi/TiNi复合合金,并分析了时效时间对TiNi/TiNi复合合金内摩擦力的影响. 相较于以上研究方向,对炸药能效所开展的研究较少. 事实上,炸药的能量利用率直接关系到材料的复合效果与制备成本,传统爆炸焊接对炸药的能量利用率极低,且在生产过程中会带来粉尘、噪音等次生污染. 针对此类问题,多名学者以不同角度给出了改进方法. 吴晓明等人[9]提出了一种带夹层钛铝复合板的爆炸焊接方法,能够提高近14% 的炸药能量利用率;Mori等人[10]研究了商业纯钨与F82H素体钢板的水下爆炸焊接效果.

      由于爆炸焊接工艺的复杂性、高压性和瞬时性等特点,在试验中直接观察爆炸焊接的微观现象难度较大,因此仿真模拟对于研究爆炸焊接的过程机理有着不可忽视的作用. Liu等人[11]利用自行开发的水力学程序模拟了详细的爆炸焊接过程,并通过数值分析研究了复合板的结合机制;Sui等人[12]利用ANSYS/LS-DYNA平台建立了三维有限元模型,仿真研究了不锈钢管与铝管的爆炸焊接机理. 然而,关于炸药覆盖层对爆炸焊接影响的研究较少. 事实上,提高炸药利用率的有效途径之一就是在炸药顶部合理加装覆盖物. 试验在杨明等人[13]的试验基础上,利用ANSYS/LS-DYNA软件和光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)算法对胶体水包覆炸药的爆炸焊接进行数值模拟,综合考虑了模拟精度与计算效率. 通过将仿真结果与试验结果进行对比,分析了在胶体水包覆炸药的爆炸焊接仿真中该算法的有效性. 此外,二维SPH法模拟的基板和覆板的界面波形也与试验结果吻合较好.

    • 文献[13]提出了一种胶体水包覆炸药的爆炸焊接实验,利用胶体水对炸药施加顶端约束达到提升炸药爆速的效果. 基于此试验利用SPH-FEM(smoothed particle hydrodynamics-finite element method)耦合法并结合ANSYS/LS-DYNA软件建立了带胶体水覆层的爆炸焊接三维模型. 考虑到计算效率及计算精度仅建立1/2模型进行计算,图1为计算模型示意图. 基板、覆板等小变形区域采用有限元网格划分,网格长度取0.8 mm;炸药和胶体水覆层此类大变形区域则使用SPH无网格法生成光滑粒子,粒子大小 Δr 取0.8 mm. 计算模型中基板尺寸为300 mm × 75 mm × 20 mm,覆板尺寸为300 mm × 75 mm × 2.5 mm,炸药厚度为10 mm,基板与覆板间隙为8 mm,胶体水覆层厚度分别为0,15,30,45 mm. 起爆方式为点起爆,忽略空气作用. 模型中所采用的单位制均为 cm,g,μs.

      图  1  计算模型示意图

      Figure 1.  Schematic diagram of calculation model

    • 文中选用乳化炸药作为爆炸材料,炸药的起爆方式设置为点起爆. 在数值计算中,乳化炸药采用高能燃烧模型[14]及Jones-Wilkings-Lee(JWL)[15]状态方程. JWL状态方程表达为

      $$ P = {A_{{\rm{JWL}}}}\left(1 - \frac{\omega }{{{R_1}V}}\right){e^{ - {R_1}V}} + {B_{{\rm{JWL}}}}\left(1 - \frac{\omega }{{{R_2}V}}\right){e^{ - {R_2}V}} + \frac{{\omega {E_0}}}{V} $$ (1)

      式中:P为爆轰产物压力;AJWLBJWLR1R2ω为材料常数;V为爆轰气体产物的相对比容,为无量纲量;E0为初始比内能. 炸药的相关参数具体见表1[16].

      表 1  乳化炸药的JWL状态方程参数

      Table 1.  JWL equation-of-state parameters of emulsion explosive

      炸药爆速
      D/(m∙s−1)
      炸药密度
      ρ/(g·cm−3)
      单位体积内能
      E0 /(kJ∙cm−3)
      材料常数1
      AJWL /GPa
      材料常数2
      BJWL/GPa
      材料常数3
      R1
      材料常数4
      R2
      材料常数5
      ω
      3 0271.013.323326.425.808 95.81.560.57

      基板与覆板均采用Johnson-Cook材料模型和Mie-Gruneisen状态方程[17]. Johnson-Cook材料模型的表达式为

      $$ \sigma = (1 + C{\rm{ln}}\mathop {\varepsilon _p^*}\limits^ \bullet )(1 - {T^*}^m)(A + B\varepsilon _p^n) $$ (2)

      式中:ABCmn为与材料相关的常数;εp为有效塑性应变;$ \mathop {\varepsilon _p^*}\limits^ \bullet = \mathop {\varepsilon _p^{}}\limits^ \bullet /\mathop {\varepsilon _p^0}\limits^ \bullet $为有效塑性应变率;$ \mathop {\varepsilon _p^0}\limits^ \bullet $为参考应变率;T* = (TTr)/(TmTr)为无量纲温度,其中Tm为熔点,Tr为室温. Q235钢与304不锈钢的Johnson-Cook材料模型参数如表2[16]所示.

      表 2  Q235钢与304不锈钢的Johnson-Cook材料模型参数

      Table 2.  Parameters of Johnson-Cook mode1l of Q235 steel and SUS304 steel

      材料密度
      ρ/(g·cm−3)
      剪切模量
      G/GPa
      初始屈服
      强度A/GPa
      硬度常数
      B/GPa
      硬化指数
      n
      应变率强化
      参数c
      软化指数
      m
      室内温度
      Tr /K
      金属熔点
      Tm /K
      Q2357.8577.00.7920.510.260.0141.03294.01 493
      SUS3047.9324.00.7001.300.750.0210.90294.01 454

      Mie-Gruneisen状态方程[18]多用于材料产生大变形的情况,即

      $$ P = \dfrac{{{\rho _0}{C^2}\mu \left[ {1 + \left( {1 - \dfrac{{{\gamma _0}}}{2}} \right)\mu - \dfrac{a}{2}{\mu ^2}} \right]}}{{{{\left[ {1 - \left( {{S_1} - 1} \right)\mu - {S_2}\dfrac{{{\mu ^2}}}{{\mu + 1}} - {S_3}\dfrac{{{\mu ^3}}}{{{{\left( {\mu + 1} \right)}^2}}}} \right]}^2}}} + \left( {{\gamma _0} + a\mu } \right)E $$ (3)

      式中:C为材料体积声速;γ0为Gruneisen系数;aγ0的一阶体积校正系数;S1S2S3为拟合系数;μ = ρ/ρ0 − 1,ρ为当前密度,ρ0为初始密度. Q235与304的Gruneisen状态方程参数见表3[16].

      表 3  Q235钢与304不锈钢的Gruneisen状态方程参数

      Table 3.  Gruneisen EOS parameters of Q235 steel and SUS304 steel

      材料体积声速
      C/(km·s−1)
      斜率系数
      S
      Gruneisen系数
      γ0
      体积修正系数
      a
      Q2356.01.492.170.46
      SUS3044.51.491.930.50

      设置于炸药上部的胶体水覆层采用了Null材料模型,此模型多用于定义流体与气体等无剪切刚度的材料,状态方程则采用了与基覆板相同的Mie-Gruneisen状态方程. 具体参数如表4所示[16].

      表 4  胶体水覆层的材料模型与状态方程参数

      Table 4.  Model and EOS parameters of Colloidal water

      密度
      ρ/(g·cm−3)
      截止压力
      PC/Pa
      动态粘度系数
      MU/(10−4 N·s·m−2)
      体积声速
      C/(km·s−1)
      斜率
      系数S
      Gruneisen
      系数γ0
      0.998−10.08.6840.164 71.9210.35
    • 爆炸焊接中工艺参数的选择对金属材料的复合质量有很大影响,如覆板冲击速度Vp、动态碰撞角β和碰撞点移动速度Vc等. 一般来说,由任意两种或以上参数组合即可构成一个平面区域,此区域能给出不同材料组合的爆炸焊接性能范围,因此称为“爆炸焊接窗口”. 基于前文献[13]中的爆炸焊接试验,选取碰撞点速度Vc和覆板碰撞速度Vp构建爆炸焊接窗口.

    • 碰撞点的射流产生情况是影响爆炸焊接质量的重要因素之一. 通常来说,碰撞点的速度与碰撞点射流的形成关系密切,且碰撞点速度不应超过所选材料体积声速的1.2倍,否则会导致射流无法形成. 因此,一般将材料的体积声速设为碰撞点速度上限,计算公式[19],即

      $$ {V_{{\rm{c}}\max }} = {C_0} $$ (4)

      碰撞点速度上限由公式(4)得Vcmax = 4 500 m/s,又由于平行法爆炸焊接中炸药爆速Vd等于碰撞点速度Vc[20],故选用炸药的爆速不得超过4 500 m/s.

      爆炸焊接中基覆板的碰撞点速度存在一种下限,只有当碰撞点速度大于速度下限时,才能获得高强度的波状结合界面. 碰撞点速度下限计算公式[21]

      $$ {V_{{\rm{c}}\min}} = \sqrt {\frac{{2{R_{\rm{e}}}\left( {H{V_{\rm{b}}} + H{V_{\rm{f}}}} \right)}}{{{\rho _{\rm{b}}} + {\rho _{\rm{f}}}}}} $$ (5)

      式中:ρbρf分别为基板和覆板的密度;Re为雷诺数;HVbHVf为基覆板的维氏硬度. 结合表5中所列出的基覆板性能参数,求得碰撞速度下限为Vcmin = 2 000 m/s. 平行法爆炸焊接中,炸药爆速Vd等于碰撞点速度Vc[20],故炸药的爆速不得低于2 000 m/s.

      表 5  Q235钢与304不锈钢的材料性能

      Table 5.  The material properties of Q235 steel and SUS304 steel

      材料密度
      ρ/(g·cm−3)
      材料声速
      C0/(km·s−1)
      拉伸强度
      σb/GPa
      维氏硬度
      HV/GPa
      材料熔点
      Tm/℃
      热导率
      κ/(W·m−1·℃−1)
      Q2357.856.000.4051.301 49338
      SUS3047.934.500.5601.701 45422
    • 为覆板碰撞速度设定上限能够防止在焊接过程中基覆板间出现连续的界面融化区,式(6)[22]给出了一种覆板碰撞速度上限理论,即

      $$ {V_{p\max }} = \frac{1}{N}\frac{{{{\left( {{T_m}{C_0}} \right)}^{\tfrac{1}{2}}}}}{{{V_c}}}{\left( {\frac{{\kappa Cp{C_0}}}{{\rho h}}} \right)^{\tfrac{1}{4}}} $$ (6)

      式中:κ为热导率;Tm为材料熔点;N为钢的常数取0.062;C0为材料声速;Cp为比热容;ρ为材料密度;h为覆板厚度取2.5 mm. 由式(5)和式(6)结合表5可求得覆板碰撞速度上限Vpmax = 950 m/s.

      由于覆板碰撞速度Vp应保证金属材料的塑性变形,因此以能够使碰撞点处的冲击压力大于金属材料的屈服应力的最小碰撞速度来定义碰撞速度下限Vpmin,可表示为[20]

      $$ {V_{p\min }} = {\left( {\frac{{{\sigma _b}}}{\rho }} \right)^{\tfrac{1}{2}}} $$ (7)

      式中:Vpmin为覆板最小碰撞速度;ρ为覆板密度;σb为拉伸强度. 结合数据计算得出覆板的最小碰撞速度为256 m/s,即覆板的碰撞速度必须大于最小碰撞速度Vpmin = 256 m/s. 由式(4) ~ 式(7)并结合具体参数绘制了爆炸焊接窗口示意图,如图2所示.

      图  2  爆炸焊接窗口

      Figure 2.  Explosive welding window.

    • 基于文献[13]试验中所设置的传感器位置,如图3所示,在基覆板结合面上选取一对特征单元A(元素17911)和B(元素532890).

      图  3  结合面上的一对特征单元

      Figure 3.  A pair of element on the joint surface.

      图4是胶体水覆层厚度分别取0,15,30,45 mm时特征单元A与B的速度-时间分布图,由图4可见,位于基板上的单元B在与覆板碰撞产生负向速度前还具有一个正向速度峰. 出现这种小型的正向速度峰的原因是基覆板在焊接结合时会产生振动位于碰撞点前方的待焊接区域受振动的影响产生了正向翘曲[23]. 结合面上所取单元的最大碰撞速度分别为364.7 m/s,508.1 m/s,576.6 m/s,615.9 m/s,这与文献[13]中所测得的碰撞速度360 m/s,500 m/s,567 m/s,619 m/s 基本一致.

      图  4  不同覆层厚度下特征单元A与B的速度-时间历程

      Figure 4.  Velocity-time history of element A and B under different covering thicknesses. (a) Colloidal water covering thickness 0 mm; (b) Colloidal water covering thickness 15 mm; (c) Colloidal water covering thickness 30 mm; (d) Colloidal water covering thickness 40 mm

      图2所示的爆炸焊接窗口对比可以发现,当胶体水覆层厚度取0 mm 和15 mm 时,覆板碰撞速度位于爆炸窗口之内;当厚度取30 mm 时,碰撞速度已经略超出爆炸窗口的区域;而当厚度到达45 mm 时,碰撞速度已经远超爆炸窗口的区域. 因此,可认为当胶体水覆层厚度取到15 mm 附近时焊接质量较高,当厚度达到乃至超过30 mm 时焊接质量相对较差.

      图5为覆板碰撞速度Vp随覆盖层厚度h增加的变化曲线. 如图5所示,覆板碰撞速度随覆层厚度的增加而增大,这意味着覆盖层能够有效提高炸药对覆板的驱动能力. 与无覆盖层相比,当覆盖厚度为15 mm,30 mm,45 mm 时,碰撞速度分别增大了39.3%,58.1%,68.8%,表明碰撞速度的增大速率呈逐渐降低的趋势.

      图  5  碰撞速度随覆层厚度变化的增益曲线

      Figure 5.  Gain curve of impact velocity with variation of covering thickness

    • 碰撞角是爆炸焊接中覆板和基板斜向碰撞所形成的夹角,其作为动态参数对控制基覆板焊接质量起着重要的作用. 图6为爆炸焊接过程中特征单元A(元素17911)和B(元素532890)结合处的碰撞角示意图,当胶体水覆层厚度为0 mm,15 mm,30 mm,45 mm时,通过ImageJ软件测出碰撞角为7.1°,8.6°,9.4°,9.9°. 这与文献[13]试验所得碰撞角6.9°,9.0°,9.5°,10.3°基本一致,误差幅度在2.9 % ~ 4.8 % 之间.

      图  6  基板与覆板的碰撞角β

      Figure 6.  Collision angle β of base plate and flyer plate

      在平行法爆炸焊接中,碰撞角β、覆板撞击速度Vp和炸药爆速Vd的关系式[24],即

      $$ {V_{\rm{d}}} = \frac{{{V_{\rm{p}}}}}{{sin\beta }} $$ (8)

      将测得的碰撞速度Vp与碰撞角β分别代入式(8),见表6所示,计算出不同厚度覆层对应的炸药爆速Vd. 未安置胶体水覆层时,炸药爆速的计算值与初始爆速基本一致;而当覆层厚度取15,30和45 mm 时,爆速的计算值分别比初始爆速增加了12.3%,16.6%和18.3%. 与覆板碰撞速度随覆层厚度变化的关系相似,炸药爆速同样随覆层厚度的增大而上升,而爆速的增长速率则呈逐渐降低的趋势. 实际上,影响炸药爆速的因素有很多,稀疏波就是其中之一. 这种稀疏波在向爆炸产物内部传播的同时降低了炸药的爆速,而在炸药上层铺设胶体水覆层后,稀疏波对炸药的影响得到了部分控制,相应的炸药爆速也有所提升. 这与文献[13]所提出的观点基本一致.

      表 6  不同覆层厚度下的炸药爆速计算值

      Table 6.  Calculation value of explosive detonation velocity under different cladding thickness

      覆层厚度h/mm炸药爆速Vd/(m·s−1)
      02 950
      153 398
      303 530
      453 582
    • 利用后处理软件LS-Prepost导出胶体水覆层厚度为15 mm 时覆板在整个爆炸焊接过程中的压力分布,如图7所示,能够发现不同时间段中覆板上的最大压力轮廓与炸药爆炸时产生的弧形爆轰波轮廓基本一致. 当25 μs 时,覆板开始与基板相撞并产生了第二道压力轮廓,其与第一道压力轮廓之间的距离保持相对恒定. 在压力轮廓移动的过程中,峰值压力集中于覆板上受爆轰波作用的部分和覆板与基板相撞的部分,而其余区域的压力几乎无明显变化.

      图  7  覆层厚度15 mm时复板的压力云图

      Figure 7.  The pressure cloud diagram of the flyer plate when covering thickness is 15 mm. (a) 4.99 μs; (b) 18.99 μs; (c) 25.96 μs; (d) 46.99 μs; (e) 74.99 μs; (f) 90.98 μs

      图8所示,在覆板上设置3个特征单元(元素8551,元素18201和元素24211),并将各特征单元的碰撞压力分布情况记录于图9中. 可以看出,在焊接的起始阶段,覆板上距离起爆端较近的部分碰撞压力较小. 而随着测点与起始端间距的逐渐增加,覆板的碰撞压力也在不断升高. 这种现象是由炸药在起爆初始阶段能量尚不稳定引起的,被称为边界效应,而在实际生产环节中要合理计算调整静态参数,尽量降低边界效应的作用范围[20].

      图  8  15 mm覆层下复板界面结合处的三个特征单元

      Figure 8.  Three element at the junction of the flyer plate under the 15 mm covering

      图  9  三个特征单元的碰撞压力-时间分布图

      Figure 9.  Collision pressure-time distribution diagram of three element

      爆炸焊接中,碰撞压力与碰撞速度的关系可表示为[20]

      $$ P = \dfrac{{{\rho _1}{v_{s,1}}{v_p}}}{{1 + \dfrac{{{\rho _1}{v_{s,1}}}}{{{\rho _2}{v_{s,2}}}}}} $$ (9)

      式中:P为碰撞压力;ρ1为基板密度,ρ2为覆板密度;vs1为基板材料声速,vs2为覆板材料声速; vp为碰撞速度. 由式(9)可知,在基覆板密度与材料声速不变时碰撞压力与碰撞速度成正比,再结合覆板碰撞速度随覆层厚度的增大而上升,提出碰撞压力与覆层厚度之间也具有相似的变化趋势的猜想.

      为进一步验证胶体水覆层厚度对碰撞压力的影响,在基覆板结合面上选取图3所示的特征单元A(元素17911),图10展示了覆层厚度取0 mm,15 mm,30 mm和45 mm 时特征单元A上的碰撞压力变化情况. 如图10所示,在不同覆层厚度下特征单元A所测得的碰撞压力分别为2.88 GPa,4.06 GPa,4.77 GPa和5.20 GPa. 与无覆层相比,15 mm,30 mm和45 mm 厚度覆层下的碰撞压力分别增大了41.0%,65.6%和80.6%. 可以发现,随着覆层厚度的增大碰撞压力逐渐上升,且这种上升速率呈递减趋势. 这说明胶体水覆层能够在一定程度上增加基覆板的碰撞压力,但这种效果会随覆层厚度的增加而减弱. 同时也验证了所提出的碰撞压力与胶体水覆层厚度关系的猜想.

      图  10  不同覆层厚度下特征单元A (元素17911)的碰撞压力分布

      Figure 10.  Impact pressure distribution of element A (element 17911) under different covering thicknesses. (a) Covering thickness 0 mm; (b) Covering thickness 15 mm; (c) Covering thickness 30 mm; (d) Covering thickness 45 mm

      图11为不同覆层厚度下覆板的位移云图和覆板上3个特征单元(元素3601,元素14221,元素25111)的位移-时间历程图. 可以看出,各组覆板的位移量均接近8 mm,即基覆板的间隙距离. 这表明四组爆炸焊接试验中基覆板均成功复合,并未出现大面积脱焊等严重质量问题,而具体焊接质量则还需要进一步分析.

      图  11  不同覆层厚度下的位移云图和位移-时间历程图

      Figure 11.  Displacement cloud diagram and displacement-time history diagram under different covering thickness. (a) covering thickness 0 mm; (b) covering thickness 15 mm; (c) covering thickness 30 mm; (d) covering thickness 45 mm

      另一方面,各特征单元位移量随覆层厚度的增加略有增大,这是因为碰撞界面的能量大小与覆板碰撞速度成正比关系. 在胶体水覆层的作用下,覆板的碰撞速度得到提升,促使碰撞界面获得了更多能量,因此特征单元的位移量也随之增大.

    • 与传统的拉格朗日法相比,SPH-FEM耦合法在模拟爆炸焊接时可以节省计算时间,提高仿真精度. 然而,基于有限元法分割的金属板材在模拟过程中无法反映复合界面的波形特征,因此选择SPH算法对覆盖层厚度为15 mm 的基覆板的爆炸焊接试验进行二维数值模拟. 板材和炸药均采用SPH粒子化生成,旨在研究胶体水覆盖作用下的基覆板复合界面形貌. 二维SPH法模拟所用状态方程和材料参数与三维建模完全相同,区别在于建模过程中的三维与二维选择上,且三维SPH-FEM耦合法模拟中SPH部分与有限元网格部分之间需要严格设置接触,而二维SPH法模拟无需添加接触选项.

      图12为覆层厚度15 mm复合界面波形对比图. 其中,图12a是文献[13]中试验得到的复合界面金相图,图12b是二维SPH算法模拟得到的复合界面粒子图. 二维模拟中所设置粒子直径均为相同大小Δr = 0.015 cm,利用图像分析软件ImageJ对图12b中的波形粒子数进行分析,计算出二维模拟所得波形的波高约为150 ~ 160 μm,波长约350 ~ 400 μm. 通过对比两图可以发现,试验得到的复合界面波形与模拟得到的复合界面波形相似,说明SPH算法很好地复现了界面波的产生过程,且基覆板复合界面的结合质量较为优良.

      图  12  覆层厚度15mm复合界面波形对比图

      Figure 12.  Waveform comparison diagram on the composite interface when the covering thickness is 15 mm. (a) metallographic analysis composite interface; (b) 2D simulation composite interface

    • (1)利用ANSYS/LS-DYNA软件与光滑粒子流体动力学算法,对胶体水包覆炸药的爆炸焊接试验进行了数值模拟. 覆层厚度设置为0,15,30和45 mm 的情况下,覆板的碰撞速度分别为364.7,508.1 ,576.6和615.9 m/s,碰撞压力分别为2.88 ,4.06,4.77和5.20 GPa. 模拟结果与前期工作中通过试验测得的数据较为吻合. 与无覆层对照组相比,覆层厚度为15,30和45 mm 时,碰撞速度分别提高了39.3%,58.1% 和68.8%,碰撞压力分别提高了41.0%,65.6% 和80.6%. 表明胶体水覆层的设置能够显著提高爆炸焊接中炸药的能量利用率.

      (2)在爆炸焊接过程中,峰值压力集中于覆板上受爆轰波实时作用的区域,且不同时间段里覆板上的最大压力轮廓与爆轰波轮廓基本相同. 在焊接起始阶段,覆板上靠近起爆点部分的碰撞压力较小,而随着与起爆端距离的不断增大,其余区域的碰撞压力也在逐渐增加. 此外,覆层厚度的增大会导致覆板的纵向位移量略微增加. 这是因为覆层的存在一定程度上提升了覆板的碰撞速度,促进了碰撞界面获得更多的能量,覆板的纵向位移量也由此得到了增加.

      (3)在覆层厚度取15 mm的情况下,基于SPH法对爆炸焊接试验进行二维数值模拟,得到了一种连续的波形界面. 其波高约150 ~ 160 μm,波长约350 ~ 400 μm. 此模拟结果与前试验得到的金相图具有一致性. 综合考虑到基覆板的动态参数均位于焊接窗口内,判断此复合界面的结合质量较为优良. 利用SPH算法很好地演示了界面波的产生过程,表明SPH法对于爆炸焊接中结合界面的观测是有效的.

参考文献 (24)

目录

    /

    返回文章
    返回