装药方式对铜/钢爆炸焊接界面波的影响及波形成机理

缪广红 马雷鸣 李雪交 艾九英 赵文慧 马宏昊 沈兆武

引用本文:
Citation:

装药方式对铜/钢爆炸焊接界面波的影响及波形成机理

    作者简介: 缪广红(1985-),男,博士,副教授,主要从事含能材料、爆炸复合及爆炸力学相关研究. E-mail:miaogh@mail.ustc.edu.cn;
  • 中图分类号: O389

Effect of Charge Mode on Interface Wave of Copper/Steel Explosive Welding and Wave Formation Mechanism

  • CLC number: O389

  • 摘要: 为了改善爆炸焊接质量,解决高噪低效的问题,选取Cu为复板、Q235钢为基板,采用LS-DYNA软件和光滑粒子流体动力学方法分别设计了均匀布药和梯形布药方案,研究了硝铵炸药对爆炸焊接界面波的影响。均匀布药结果显示:沿着爆轰方向碰撞压力逐渐增大;炸药量越多,碰撞压力越大,界面波波形越大。梯形布药方案中,通过改变炸药起爆端和末端的高度,设计了4种方案,结果显示:梯形布药可以消除爆炸焊接界面波不均匀现象,使界面波形尺寸基本保持一致,而且节省了炸药用量;当起爆端和末端的高度分别为67.2 mm和42.0 mm时,波形效果最好。通过研究界面波的形成过程可知,SPH法模拟的界面波形成过程与复板流侵彻机理的一致性较好,证明了复板流侵彻机理解释界面波形成过程的有效性。
  • 图 1  计算模型

    Figure 1.  Computational model

    图 2  爆炸焊接窗口

    Figure 2.  Explosive welding window

    图 3  均匀布药模拟效果图

    Figure 3.  Simulation effect diagram of uniform charge

    图 4  关键点取样示意图

    Figure 4.  Schematic diagram of key point sampling

    图 5  关键点碰撞压力值折线

    Figure 5.  Line diagram of collision pressure of key points

    图 6  梯形装药结构示意图

    Figure 6.  Schematic diagram of ladder charge structure

    图 7  关键点碰撞压力值折线图

    Figure 7.  Line diagram of collision pressure of key points

    图 8  波形成示意图

    Figure 8.  Illustration of wave formation

    表 1  硝铵炸药的JWL状态方程参数

    Table 1.  JWL EOS parameters of ammonium nitrate explosive

    ρ/(kg·m−3)D/(m·s−1)AJ/GPaBJ/GPaR1R2ω
    8002 800132.750.4235.31.20.21
    下载: 导出CSV

    表 2  Cu和Q235钢的Johnson-Cook模型参数

    Table 2.  Parameters of Johnson-Cook model of Cu and Q235 steel

    Materialρ/(g·cm−3)G/GPaA/GPaB/GPanCmTm/KTr/K
    Cu8.96460.0900.2920.310.0251.091 356294
    Q2357.83770.7920.5100.260.0141.031 793294
    下载: 导出CSV

    表 3  Cu和Q235钢的Grüneisen方程参数

    Table 3.  Grüneisen EOS parameters of Cu and Q235 steel

    Materialc/(km·s−1)S1Γ0a
    Cu3.9401.4892.020.47
    Q2354.5691.4902.170.46
    下载: 导出CSV

    表 4  均匀布药方案关键点碰撞压力

    Table 4.  Collision pressure of key points in uniform charge scheme

    Key pointPressure/GPaKey pointPressure/GPa
    R1 = 1.0R2 = 1.5R1 = 1.0R2 = 1.5
    A10.5811.602A54.9997.086
    A21.6722.002A65.7548.576
    A34.5075.289A70.5261.031
    A44.6546.191
    下载: 导出CSV

    表 5  梯形布药方案

    Table 5.  Ladder charging scheme

    Schemea/mmb/mm
    67.258.8
    67.250.4
    67.242.0
    67.233.6
    下载: 导出CSV

    表 6  梯形布药方案关键点碰撞压力

    Table 6.  Collision pressure of key points of ladder charge scheme

    Key pointPressure/GPa
    Scheme ⅠScheme ⅡScheme ⅢScheme Ⅳ
    A11.7381.1581.6940.351
    A22.5492.5482.1411.936
    A33.8393.8303.6212.776
    A45.2307.5475.4395.665
    A57.8908.3303.4134.537
    A66.4233.5053.7312.064
    A70.7800.4310.3970.435
    下载: 导出CSV
  • [1] NASSIRI A, KINSEY B. Numerical studies on high-velocity impact welding: smoothed particle hydrodynamics (SPH) and arbitrary Lagrangian-Eulerian (ALE) [J]. Journal of Manufacturing Processes, 2016, 24: 376–381. doi: 10.1016/j.jmapro.2016.06.017
    [2] ABE A. Numerical study of the mechanism of wavy interface generation in explosive welding [J]. JSME International Journal Series B—Fluids and Thermal Engineering, 1997, 40: 395–401. doi: 10.1299/jsmeb.40.395
    [3] YUAN X, WANG W, CAO X, et al. Numerical study on the interfacial behavior of Mg/Al plate in explosive/impact welding [J]. Science & Engineering of Composite Materials, 2017, 24(6): 833–843.
    [4] TABBATAEE M, MAHMOUDI J. Finite element simulation of explosive welding [J]. Journal of Applied Physics, 2014, 24(3): 349–359.
    [5] MOUSAVI A A A, BURLEY S J, AL-HASSANI S T S. Simulation of explosive welding using the Williamsburg equation of state to model low detonation velocity explosives [J]. International Journal of Impact Engineering, 2005, 31(6): 719–734. doi: 10.1016/j.ijimpeng.2004.03.003
    [6] MOUSAVI A A A, AL-HASSANI S T S. Simulation of wave and jet formations in explosive/impact welding [C]//ASME 7th Biennial Conference on Engineering Systems Design and Analysis. Manchester, England, 2004: 265–274.
    [7] 王宇新, 李晓杰, 孙国, 等. 无网格MPM法三维爆炸焊接数值模拟 [J]. 计算力学学报, 2013, 30(1): 34–38. doi: 10.7511/jslx201301006
    WANG Y X, LI X J, SUN G, et al. Three dimensional simulation of the explosive welding by using of the MPM [J]. Chinese Journal of Computational Mechanics, 2013, 30(1): 34–38. doi: 10.7511/jslx201301006
    [8] 刘江, 郑远远, 沈宗宝, 等. 基于SPH方法的爆炸焊接过程模拟 [J]. 焊接技术, 2013, 42(12): 17–20.
    LIU J, ZHENG Y Y, SHEN Z B, et al. Simulation of explosive welding process based on SPH method [J]. Welding Technology, 2013, 42(12): 17–20.
    [9] 周春华, 史长根, 蔡立艮, 等. 爆炸焊接布药工艺的研究 [J]. 焊接技术, 2002, 31(6): 17–18. doi: 10.3969/j.issn.1002-025X.2002.06.008
    ZHOU C H, SHI C G, CAI L G, et al. Research on dynamite-distributing technology of explosive welding [J]. Welding Technology, 2002, 31(6): 17–18. doi: 10.3969/j.issn.1002-025X.2002.06.008
    [10] 董刚, 周春华, 史长根, 等. 爆炸焊接不等厚度布药工艺 [J]. 焊接, 2004(6): 35–38. doi: 10.3969/j.issn.1001-1382.2004.06.010
    DONG G, ZHOU C H, SHI C G, et al. Unequal thickness arranging explosive technology of explosive welding [J]. Welding, 2004(6): 35–38. doi: 10.3969/j.issn.1001-1382.2004.06.010
    [11] 缪广红, 李亮, 江向阳, 等. 双面爆炸焊接的数值模拟 [J]. 高压物理学报, 2018, 32(4): 1–8. doi: 10.11858/gywlxb.20180513
    MIAO G H, LI L, JIANG X Y, et al. Numerical simulation of double sided explosive welding [J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 1–8. doi: 10.11858/gywlxb.20180513
    [12] LEE E, FINGER M, COLLINS W. JWL equation of state coefficients for high explosives [R]. Livermore, CA, USA: Lawrance Livermore National Laboratory, 1973.
    [13] LIU G R, LIU M B. 光滑粒子流体动力学——一种无网格粒子法 [M]. 韩旭, 译. 长沙: 湖南大学出版社, 2005.
    [14] 程国强, 李守新. 金属材料在高应变率下的热粘塑性本构模型 [J]. 弹道学报, 2004, 11(6): 18–22.
    CHENG G Q, LI S X. Thermal viscoplastic constitutive model of metallic materials at high strain rate [J]. Journal of Ballistics, 2004, 11(6): 18–22.
    [15] 张振逵, 吴绍尧. 用半圆柱法测定铜-钢爆炸焊接窗口及合理药量 [J]. 焊接学报, 1980(3): 17–30, 67.
    ZHANG Z K, WU S Y. Determination of explosive welding window and reasonable charge content of copper-steel by semi-cylindrical method [J]. Transactions of the China Welding Institution, 1980(3): 17–30, 67.
    [16] SUI G F, LI J S, SUN F, et al. 3D finite element simulation of explosive welding of three-layer plates [J]. Science China-Physics Mechanics & Astronomy, 2011, 54(5): 890–896.
    [17] 孙锦山, 朱建士. 理论爆轰物理 [M]. 北京: 国防工业出版社, 1995: 356–418.
    [18] MOUSAVI A A A, AL-HASSANI S T S. Finite element simulation of explosively-driven plate impact with application to explosive welding [J]. Materials & Design, 2008, 29(1): 1–19. doi: 10.1016/j.matdes.2006.12.012
    [19] 蔡立艮, 卢红标, 周春华, 等. 爆炸焊接布药工艺与微观结合界面形貌分析 [J]. 爆破, 2010, 27(1): 78–81. doi: 10.3963/j.issn.1001-487X.2010.01.021
    CAI L G, LU H B, ZHOU C H, et al. Arranging explosive technology of explosive welding and microanalysis of bonging interfaces [J]. Blasting, 2010, 27(1): 78–81. doi: 10.3963/j.issn.1001-487X.2010.01.021
    [20] 王克鸿, 张德库, 张文军. 爆炸焊接技术研究进展 [J]. 机械制造与自动化, 2011, 40(2): 1–5. doi: 10.3969/j.issn.1671-5276.2011.02.001
    WANG K H, ZHANG D K, ZHANG W J. Research progress of explosive welding technology [J]. Mechanical Manufacturing and Automation, 2011, 40(2): 1–5. doi: 10.3969/j.issn.1671-5276.2011.02.001
    [21] FINDIK F. Recent developments in explosive welding [J]. Materials & Design, 2011, 32(3): 1081–1082.
    [22] 袁晓丹. 铝-镁合金爆炸焊接层状复合界面形成机制及数值模拟 [D]. 太原: 太原理工大学, 2016.
    YUAN X D. Formation mechanism and numerical simulation of layered composite interface in explosive welding of Al-Mg alloy [D]. Taiyuan: Taiyuan University of Technology, 2016.
    [23] LI Y, WU Z. Microstructural characteristics and mechanical properties of 2205/AZ31B laminates fabricated by explosive welding [J]. Metals, 2017, 7(4): 125. doi: 10.3390/met7040125
    [24] 郑远谋. 爆炸焊接和爆炸复合材料 [M]. 北京: 国防工业出版社, 2017: 13–14.
    [25] 缪广红. 蜂窝结构炸药与双面爆炸复合的研究 [D]. 合肥: 中国科学技术大学, 2015.
    MIAO G H. Research on honeycomb structure explosives and double sided explosive cladding [D]. Hefei: University of Science and Technology of China, 2015.
    [26] 王耀华. 金属板材爆炸焊接研究与实践 [M]. 北京: 国防工业出版社, 2007.
    [27] 缪广红, 马宏昊, 沈兆武, 等. 不锈钢-普碳钢的双面爆炸复合 [J]. 爆炸与冲击, 2015, 35(4): 536–540. doi: 10.11883/1001-1455(2015)04-0536-05
    MIAO G H, MA H H, SHEN Z W, et al. Double-sided explosive recombination of stainless steel and plain carbon steel [J]. Explosion and Shock Waves, 2015, 35(4): 536–540. doi: 10.11883/1001-1455(2015)04-0536-05
    [28] WANG X, ZHENG Y, LIU H, et al. Numerical study of the mechanism of explosive/impact welding using smoothed particle hydrodynamics method [J]. Materials & Design, 2012, 35: 210–219. doi: 10.1016/j.matdes.2011.09.047
    [29] BAHRANI A S, BLACK T J, CROSSLAND B. The mechanics of wave formation in explosive welding [J]. Proceedings of the Royal Society of London, Series A, Mathematical and Physical Sciences, 1967, 296(1445): 123–136.
    [30] 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(11): 3145–3155. doi: 10.1007/BF02814967
    [31] COWAN G R, HOLTZMAN A H. Flow configurations in colliding plates: explosive bonding [J]. Journal of Applied Physics, 1963, 34(4): 928–939. doi: 10.1063/1.1729565
    [32] KOWALICK J F, HAY D R. A mechanism of explosive bonding [J]. Metallurgical and Materials Transactions B, 1971, 2(7): 1953–1958.
    [33] REID S R, SHERIFF N H S. Prediction of the wave length of interface waves in symmetric explosive welding [J]. Journal of Mechanical Engineering Science, 1980, 18(2): 87–94.
    [34] GODUNOV S K, DERIBAS A A, ZABRADINA V. Hydrodynamic effect in colliding solids [J]. Computational Physics, 1970, 5: 517–539. doi: 10.1016/0021-9991(70)90078-1
  • [1] 李晓杰莫非闫鸿浩张程娇 . 爆炸焊接斜碰撞过程的数值模拟. 高压物理学报, 2011, 25(2): 173-176 . doi: 10.11858/gywlxb.2011.02.014
    [2] 任会兰宁建国许香照 . 不同炸药量在工事中爆炸的三维数值模拟. 高压物理学报, 2013, 27(2): 216-222. doi: 10.11858/gywlxb.2013.02.008
    [3] 李金河赵继波谭多望王彦平张远平 . 不同起爆方式对含铝炸药水中爆炸近场冲击波压力的影响. 高压物理学报, 2012, 26(3): 289-293. doi: 10.11858/gywlxb.2012.03.007
    [4] 哈跃管公顺迟润强庞宝君 . 铝合金弹丸超高速撞击玄武岩纤维布剩余速度分析. 高压物理学报, 2012, 26(3): 273-280. doi: 10.11858/gywlxb.2012.03.005
    [5] 哈跃刘志勇管公顺庞宝君 . 玄武岩纤维布超高速撞击损伤分析. 高压物理学报, 2012, 26(5): 557-563. doi: 10.11858/gywlxb.2012.05.012
    [6] 张宝玺哈跃邓云飞庞宝君 . 超高速撞击Kevlar纤维布填充防护结构研究. 高压物理学报, 2013, 27(1): 105-112. doi: 10.11858/gywlxb.2013.01.015
    [7] 郑康陈力方秦高飞 . CFRP布/树脂胶复合条带的抗拉力学性能. 高压物理学报, 2017, 31(6): 794-802. doi: 10.11858/gywlxb.2017.06.015
    [8] 苘廉洁苑洪胜秦礼萍张莉 . 下地幔温压条件下碳对(Mg,Fe)SiO3布里奇曼石的影响. 高压物理学报, 2019, 33(6): 060102-1-060102-10. doi: 10.11858/gywlxb.20190788
    [9] 谭华 . 金属的冲击波温度测量(Ⅱ)──界面卸载近似. 高压物理学报, 1996, 10(3): 161-169 . doi: 10.11858/gywlxb.1996.03.001
    [10] 倪小军马宏昊沈兆武蒋耀港李磊 . 流-固界面冲击波入射压力的测试及数值模拟. 高压物理学报, 2013, 27(4): 528-534. doi: 10.11858/gywlxb.2013.04.010
    [11] 刘永贵沈玲燕 . 固定温度界面对相变波传播规律的影响. 高压物理学报, 2018, 32(4): 042301-1-042301-7. doi: 10.11858/gywlxb.20170559
    [12] 胡金彪刘吉平谭华杨嘉陵唐志平 . 冲击波作用下CHBr3/NaCl界面热弛豫特性研究. 高压物理学报, 2001, 15(2): 134-140 . doi: 10.11858/gywlxb.2001.02.010
    [13] 徐森刘大斌彭金华王建灵郭伟金朋刚贾宪政 . 药柱冲击波在有机玻璃中的衰减特性研究. 高压物理学报, 2010, 24(6): 431-437 . doi: 10.11858/gywlxb.2010.06.005
    [14] 谭华 . 金属的冲击波温度测量(Ⅰ)──高温计的标定和界面温度的确定. 高压物理学报, 1994, 8(4): 254-263 . doi: 10.11858/gywlxb.1994.04.003
    [15] 谭华 . 金属的冲击波温度测量(Ⅲ)基板/样品界面间隙对辐射法测量冲击波温度的影响. 高压物理学报, 1999, 13(3): 161-168 . doi: 10.11858/gywlxb.1999.03.001
    [16] 赵同虎韩立石何智 . 猛炸药RHT-901和钝感炸药IHE-2的爆轰波直角绕射. 高压物理学报, 1998, 12(4): 264-270 . doi: 10.11858/gywlxb.1998.04.005
    [17] 谭华韩钧万王晓江苏林祥刘利刘江崔玲 . 炸药爆炸冲击波合成纤锌矿型氮化硼. 高压物理学报, 1991, 5(4): 241-253 . doi: 10.11858/gywlxb.1991.04.001
    [18] 赵继波谭多望李金河龚晏青孙永强 . 含铝炸药水中爆炸冲击波相似律适应性探索. 高压物理学报, 2010, 24(5): 388-394 . doi: 10.11858/gywlxb.2010.05.012
    [19] 姜洋李蕾王涛钟敏李平柏劲松 . 180圆弧形钝感炸药的爆轰波传播. 高压物理学报, 2013, 27(1): 132-136. doi: 10.11858/gywlxb.2013.01.019
    [20] 黄彬彬张旭谭多望 . 谱元法数值模拟约束结构内炸药中应力波的传播. 高压物理学报, 2014, 28(6): 716-722. doi: 10.11858/gywlxb.2014.06.012
  • 加载中
图(8)表(6)
计量
  • 文章访问数:  189
  • 阅读全文浏览量:  163
  • PDF下载量:  1
出版历程
  • 收稿日期:  2019-10-08
  • 录用日期:  2019-11-11
  • 网络出版日期:  2020-03-23
  • 刊出日期:  2020-04-01

装药方式对铜/钢爆炸焊接界面波的影响及波形成机理

    作者简介:缪广红(1985-),男,博士,副教授,主要从事含能材料、爆炸复合及爆炸力学相关研究. E-mail:miaogh@mail.ustc.edu.cn
  • 1. 安徽理工大学深部煤矿采动响应与灾害防控国家重点实验室,安徽 淮南 232001
  • 2. 安徽理工大学力学与光电物理学院,安徽 淮南 232001
  • 3. 安徽理工大学土木建筑学院,安徽 淮南 232001
  • 4. 安徽理工大学化学工程学院,安徽 淮南 232001
  • 5. 中国科学院材料力学行为和设计重点实验室,中国科学技术大学,安徽 合肥 230027

摘要: 为了改善爆炸焊接质量,解决高噪低效的问题,选取Cu为复板、Q235钢为基板,采用LS-DYNA软件和光滑粒子流体动力学方法分别设计了均匀布药和梯形布药方案,研究了硝铵炸药对爆炸焊接界面波的影响。均匀布药结果显示:沿着爆轰方向碰撞压力逐渐增大;炸药量越多,碰撞压力越大,界面波波形越大。梯形布药方案中,通过改变炸药起爆端和末端的高度,设计了4种方案,结果显示:梯形布药可以消除爆炸焊接界面波不均匀现象,使界面波形尺寸基本保持一致,而且节省了炸药用量;当起爆端和末端的高度分别为67.2 mm和42.0 mm时,波形效果最好。通过研究界面波的形成过程可知,SPH法模拟的界面波形成过程与复板流侵彻机理的一致性较好,证明了复板流侵彻机理解释界面波形成过程的有效性。

English Abstract

  • 由于爆炸焊接过程具有瞬时性、复杂性,人们在通过实验及相关设备来观察、研究爆炸焊接过程时,效果十分有限。随着计算机技术的发展,数值模拟技术的突破为爆炸焊接研究提供了帮助。迄今为止,国内外诸多学者利用数值模拟软件对爆炸焊接过程及参数开展了大量研究,取得了较理想的研究成果。Nassiri等[1]分别利用光滑粒子流体动力学(Smoothed particle hydrodynamics,SPH)法、任意拉格朗日-欧拉(Arbitrary Lagrange-Euler,ALE)法模拟了爆炸焊接的高速撞击过程,并研究了结合界面的剪切应力和速度分布,结果显示只有SPH法再现了碰撞时的射流现象。Abe[2]使用二维有限差分法研究了爆炸焊接波的形成机理,计算结果与模拟结果保持一致。Yuan等[3]采用SPH法研究了6061Al/AZ31B的爆炸焊接,模拟得到波形界面和射流现象,但是并没有详细指出波的形成机理。Tabbataee等[4]使用有限单元法模拟了射流。Mousavi等[5-6]采用欧拉法模拟了射流和类似流体状的波状界面,并成功预测了波形尺寸和射流速度。王宇新等[7]应用无网格物质点法(Material point method,MPM),对爆炸荷载作用下基板和复板的动态变形过程进行了三维数值模拟,并分析了碰撞点压力和速度的变化情况。刘江等[8]用SPH法进行二维建模,重现了实验中观察到的射流和界面波现象,并揭示了碰撞压力、有效塑性应变、温度等参数随时间的变化关系。周春华等[9-10]研究了均匀布药下结合界面波的变化规律后,提出了新的不等药厚布药工艺。

    从众多学者的研究成果来看,数值模拟软件能够帮助人们全面、系统地了解爆炸焊接过程,降低实验成本,优化实验参数。在平板爆炸焊接中,通过炸药爆炸产生的爆轰能量推动复板与基板发生斜碰撞,碰撞界面在微秒级的时间内产生数千兆帕高压,发生塑性变形,从而实现基板和复板的固相焊接。由此可以看出,炸药是整个爆炸焊接过程的能量来源,因此对炸药药量和布药方式进行控制具有重要的研究意义。本研究应用LS-DYNA软件,采用SPH法进行二维建模,研究炸药与爆炸焊接界面波的关系及波的形成机理。

    • SPH法是一种无网格法,适合处理爆炸、冲击等大变形问题,但是其搜索算法较为耗时,粒子越多,该问题越突出,计算效率越低[11]。为了提高计算结果的精确性,同时兼顾计算效率,炸药采用粗网格划分,基板和复板采用细网格划分。先采用LS-DYNA软件建立二维平面计算模型,再利用LS-prepost软件将网格全部生成SPH粒子(见图1)。复板和基板材料选取Cu和Q235钢,爆炸焊接选用硝铵炸药。模型尺寸:复板150 mm × 4 mm,基板150 mm × 15 mm,间隙6 mm,基板和复板网格尺寸为0.2 mm × 0.2 mm,炸药网格尺寸为0.4 mm × 0.4 mm。采用cm-g-μs单位制。

      图  1  计算模型

      Figure 1.  Computational model

    • 硝铵炸药选用Jones-Wilkins-Lee(JWL)[12]状态方程(EOS),计算公式为

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

      式中:AJBJR1R2ω为材料系数;p为爆轰产物压力,GPa;E0为初始比内能,kJ·cm–3V为爆轰气体产物的相对比容。炸药的具体JWL状态参数见表1,其中ρ为密度,D为爆速。

      ρ/(kg·m−3)D/(m·s−1)AJ/GPaBJ/GPaR1R2ω
      8002 800132.750.4235.31.20.21

      表 1  硝铵炸药的JWL状态方程参数

      Table 1.  JWL EOS parameters of ammonium nitrate explosive

      基板、复板均采用 Mie-Grüneisen状态方程[13]和Johnson-Cook材料模型[14]。Grüneisen状态方程被广泛应用于解决材料的大变形问题,它描述了粒子速度up和冲击速度D之间的基本关系,计算公式如下

      $ \rho = {\rho _{\rm{H}}} + {\varGamma _\rho }\left( {e - {e_{\rm{H}}}} \right) $

      其中

      $ {\varGamma _\rho } = {\varGamma _0}{\rho _0} = {\rm{Constant}} $

      $ {\rho _{\rm{H}}} = \frac{{{\rho _0}{c_0}\mu \left( {1 + \mu } \right)}}{{{{\left[ {1 - \left( {s - 1} \right)\mu } \right]}^2}}} $

      $ {e_{\rm{H}}} = \frac{1}{2}\frac{{{\rho _{\rm{H}}}}}{{{\rho _0}}} {\frac{\mu }{{1 + \mu }}} $

      $ \mu = {\frac{{{\rho _{}}}}{{{\rho _0}}}} - 1 $

      式中:Γ0为Grüneisen系数,ρ为当前密度,ρ0为初始密度,c0为体积声速,s为压缩比。

      Johnson-Cook材料模型经常用来反映金属的力学行为,计算公式如下

      $ \sigma = \left( {A + B\varepsilon _{\rm p}^n} \right)\left( {1 + C\; {\rm ln}\;{\dot\varepsilon} _{\rm p}^*} \right)\left( {1 - {T^{*m}}} \right) $

      其中无量纲温度T*表示为

      $ {T^*} = \left( {T - {T_{\rm{r}}}} \right)/\left( {{T_{\rm{m}}} - {T_{\rm{r}}}} \right) $

      式中:εp为有效塑性应变;$\dot \varepsilon _p^* = {\dot \varepsilon _{\rm p}}/\dot \varepsilon _{\rm p}^0$为有效塑性应变率,$\dot \varepsilon _{\rm p}^0$为参考应变率;ABCmn是与材料相关的常数;Tr为室温;Tm为熔点。Cu和Q235钢的Johnson-Cook材料模型参数列于表2,Grüneisen EOS参数列于表3,其中:c为声速,km/s; S1为拟合系数;a为Grüneisen系数Γ0的一阶体积修正,G为剪切模量。

      Materialρ/(g·cm−3)G/GPaA/GPaB/GPanCmTm/KTr/K
      Cu8.96460.0900.2920.310.0251.091 356294
      Q2357.83770.7920.5100.260.0141.031 793294

      表 2  Cu和Q235钢的Johnson-Cook模型参数

      Table 2.  Parameters of Johnson-Cook model of Cu and Q235 steel

      Materialc/(km·s−1)S1Γ0a
      Cu3.9401.4892.020.47
      Q2354.5691.4902.170.46

      表 3  Cu和Q235钢的Grüneisen方程参数

      Table 3.  Grüneisen EOS parameters of Cu and Q235 steel

    • 爆炸焊接质量的好坏极其依赖选取的爆炸焊接参数,如药量比、间距、碰撞速度、碰撞角等。参考文献[15]中的爆炸复合窗口,选取复板厚度δ1 = 4 mm时药量比(单位面积炸药与复板质量之比R)分别为R1 = 1.0、R2 = 1.5的两组参数开展研究。经计算得到对应的炸药厚度分别为44.8 mm和67.2 mm。由图2可知,R1靠近可焊窗口下限,R2位于可焊窗口内部,两组参数均位于爆炸复合窗口内部。

      图  2  爆炸焊接窗口

      Figure 2.  Explosive welding window

    • 为了研究炸药量与焊接界面波形的关系,首先进行均匀布药模拟。分别建立炸药厚度为44.8 mm和67.2 mm的计算模型,输出模拟图像,如图3所示。

      图  3  均匀布药模拟效果图

      Figure 3.  Simulation effect diagram of uniform charge

      为了更方便地了解爆炸焊接过程参数的变化规律,在结合面上等间距区域选取7个关键点,选取位置如图4所示。输出关键点的碰撞压力值,如表4所示。硝铵炸药的爆轰压力计算公式为

      图  4  关键点取样示意图

      Figure 4.  Schematic diagram of key point sampling

      Key pointPressure/GPaKey pointPressure/GPa
      R1 = 1.0R2 = 1.5R1 = 1.0R2 = 1.5
      A10.5811.602A54.9997.086
      A21.6722.002A65.7548.576
      A34.5075.289A70.5261.031
      A44.6546.191

      表 4  均匀布药方案关键点碰撞压力

      Table 4.  Collision pressure of key points in uniform charge scheme

      $ p = \frac{{{\rho _0}{D^2}}}{{1 + \gamma }} $

      式中:γ为硝铵炸药的多方指数,ρ0为炸药密度,D为炸药爆速。

      由式(9)可以看出,模拟中使用的硝铵炸药的密度为常数,所以爆轰压力p只与炸药爆速D有关。由于炸药的爆速与单位面积的药量成正比[16],且R1 < R2,因此R2对应的爆轰压力要高于R1的爆轰压力。炸药的爆轰压力作用于复板,复板又传递到焊接界面,最终导致R2对应的基板、复板碰撞压力高于R1。根据图5的模拟结果可知,理论分析结果与模拟结果一致。

      图  5  关键点碰撞压力值折线

      Figure 5.  Line diagram of collision pressure of key points

      焊接界面波长的计算公式[17]

      $ \lambda = \frac{{{\text{π}}{v_{\rm{p}}}d}}{{2c}}{\left( {\frac{c}{{{v_{\rm{c}}}}} - 1} \right)^2} $

      式中:vp为碰撞速度;d为复板厚度;c为材料声速;vc为碰撞点移动速度,等于炸药的爆轰速度。

      从式(10)可以看出,给定金属板材的材料参数后,复板厚度d和材料的声速c为常数,碰撞速度vp和碰撞点移动速度vc为变量,界面波的波长λvp成正比,与vc成反比。

      图3可以看出,复板上部均出现了打伤现象,这是由于前期炸药和板材划分的网格尺寸不同所致。R1R2两组模拟均得到了波状焊接界面,这被认为是焊接良好的标志[18],模拟结果与文献[19]的金相实验结果吻合较好。从图3(a)图3(b)看出,焊接界面在起爆端区域没有出现波状界面,随着距起爆端距离的不断增大,界面先后出现了无波到微波、微波到大波的转变。图3(a)中射流现象最先出现在4.299 6 μs时刻,而图3(b)中的射流最先出现在4.198 2 μs,早于前者,这是由于药量比越大,爆轰能量越大,导致射流出现较早。焊接界面的放大图像显示,波状界面先于射流出现,可以认为波的形成不一定需要射流,印证了Godunov提出的应力波机理[20]。从图3(c)图3(d)可以清晰地看出,碰撞点前方在基板和复板碰撞作用下,板材表面喷射出大量射流,可以清洗掉金属表面的杂质和氧化层,提高焊接质量[21]。喷射出的射流中,基板的成分明显高于复板,基板材料是射流成分的主要来源,这是由于Cu的密度高于Q235钢,密度较小的材料更容易产生射流[22]图3(e)图3(f)为模拟图像,表明基板与复板已经完全复合。从复合板中央部分的取样结果来看,图3(f)的波状尺寸大于图3(e),由于小波增加了焊接面积,有助于提高焊接质量,波形越大,越容易出现裂缝和孔洞等缺陷,因此可以认为药量比R1优于R2的焊接质量,药量比尽量取在可焊窗口的下限。

      图5是均匀布药下各个关键点的碰撞压力折线。考虑到边界效应的影响,去掉首尾两个关键点后,折线显示沿着爆轰方向焊接界面的碰撞压力逐渐增大,由于碰撞压力过大会产生熔化层,产生的稀疏波可能将已结合的复合板拉开,降低焊接质量,因此均匀布药不是最终选择。

    • 根据均匀布药方式下界面波的分布特点,起爆端区域的结合呈平直界面,这是由于炸药起爆时能量不稳定[23],碰撞压力小。在爆炸焊接中,碰撞压力过小会导致界面未复合[16],即起爆端附近极易出现未复合现象,降低金属板材的利用率,这种现象被称为边界效应,在实验和工程实践中应尽量避免。碰撞压力要超过某个临界值,复合板才能获得良好的焊接质量,这也证明了爆炸焊接是一种压力焊[24]。沿着爆轰方向,结合界面逐渐由无波发展成微波、大波界面,这是爆炸产物的堆集以及基板和复板碰撞引起的振动能在待复合区叠加并共同作用的结果[25-26]。结合式(9)、式(10)可知,波状界面的尺寸与碰撞点的移动速度有关,碰撞压力与炸药爆速成正比,在平行法爆炸焊接装置中爆速等于碰撞点的移动速度[27],由此得出,炸药的爆轰能量赋予了复板巨大的压力,碰撞压力决定界面波形。因此可以通过调整炸药厚度来控制爆轰能量,使界面波形维持在微波状,进而提高复合板的焊接质量。

      设计如图6所示的呈梯形装药结构,通过调整炸药起爆端和末端厚度ab的值,来确定一个合适的倾斜装药角度,沿着爆轰方向控制碰撞点的压力基本保持不变,避免因大波状界面而产生涡旋和裂缝,提高复合板的结合强度,使整个焊接界面都保持在微波状结合状态。布药方案见表5

      Schemea/mmb/mm
      67.258.8
      67.250.4
      67.242.0
      67.233.6

      表 5  梯形布药方案

      Table 5.  Ladder charging scheme

      图  6  梯形装药结构示意图

      Figure 6.  Schematic diagram of ladder charge structure

      由于网格划分较大,导致波形尺寸不易测量。当网格划分较小时,可以较清晰地观察界面的形貌特征。由前述分析得知,碰撞压力决定界面波形,因此将界面上的碰撞压力作为研究对象,使界面上各点处碰撞压力基本保持一致。

      不同梯形布药方案下界面关键点的碰撞压力如图7所示,具体压力值列于表6。由图7中折线的波动情况可知,当采用梯形布药方式后,基板、复板之间的碰撞压力不再呈均匀布药下逐渐增长的趋势,而是呈较小幅度的上下波动,因此梯形布药方式可以消除不利因素,使焊接界面的碰撞压力基本保持不变,以获得波长较小、波幅较小的微波界面。观察图7发现,起爆端和末端的碰撞压力最小,这是由稀疏波引发的边界效应所致。在不考虑首尾两端的压力时,方案Ⅲ的碰撞压力波动最小,因此最终焊接界面的波形状态好于其他方案。而方案Ⅱ中界面上的碰撞压力起伏最大,最终可能导致基板与复板未复合。因此确定合适的首尾装药高度是梯形布药方案的关键。

      Key pointPressure/GPa
      Scheme ⅠScheme ⅡScheme ⅢScheme Ⅳ
      A11.7381.1581.6940.351
      A22.5492.5482.1411.936
      A33.8393.8303.6212.776
      A45.2307.5475.4395.665
      A57.8908.3303.4134.537
      A66.4233.5053.7312.064
      A70.7800.4310.3970.435

      表 6  梯形布药方案关键点碰撞压力

      Table 6.  Collision pressure of key points of ladder charge scheme

      图  7  关键点碰撞压力值折线图

      Figure 7.  Line diagram of collision pressure of key points

      综上所述,采用梯形布药方式可以消除焊接界面波形不均匀现象,而且减少了炸药用量,同时也可以减小爆炸焊接实验现场引发的噪声。通过分析各个关键点的碰撞压力可知,方案Ⅲ的复合效果最好,即炸药起爆端和末端的厚度分别为67.2 mm和42.0 mm时装药倾斜角度最佳。在实际的工程实践中,梯形布药方式的应用并不广泛,原因是其影响因素较多,参数确定比较复杂。尽管减小炸药厚度会使碰撞压力减小,但碰撞压力的变化与炸药厚度并非成严格的正比关系,所以才会得到图7的模拟结果。如何优化梯形布药参数使界面波波形基本保持不变,还需要大量深入的研究。

    • 爆炸焊接过程具有瞬时性,因此很难用实验方法来研究界面波的形成与焊接参数之间的关系。虽然有关界面波形成机理的研究已经进行了数十年,但迄今为止仍然没有达成一个统一的认识[28]。期间众多学者提出了一些波形形成机制,主要有复板流侵彻机理[29]、Helm-holtz不稳定流机理[30-31]、涡脱落机理[32-33]和应力波机理[34]。随着数值模拟软件的发展,SPH法能够清晰地模拟出波形成过程,为解释界面波的形成机理提供了帮助。

      图8展示了采用SPH法得到的爆炸焊接界面波形成过程。基板与复板在P点高速碰撞,此时P点可看作低黏塑性流体状态,巨大的碰撞压力使P点产生塑性变形。如图8(a)所示,复板与基板碰撞后形成初始射流和再入射流,再入射流在前方复板的阻碍下向斜上方喷出。沿着焊接方向向前发展,碰撞点P转移到波峰上,见图8(b),先前的再入射流被基板捕获,形成在波峰的前后方,此处极易形成漩涡、空洞等缺陷,从而降低复合板的焊接质量。碰撞点P转移到波峰上后,再入射流与波峰碰撞,沿斜下方喷出。图8(c)中碰撞点P越过波峰后形成新的碰撞点,在巨大的爆轰压力下基板再次发生塑性变形,形成新的波形,由此产生周期性的波状界面。由图8可见,SPH法模拟的波形成过程与复板流侵彻机理的一致性较好,证明了SPH法在解释爆炸焊接界面波形成过程中复板流侵彻机理的有效性。

      图  8  波形成示意图

      Figure 8.  Illustration of wave formation

    • (1)采用LS-DYNA软件和SPH法,模拟了整个爆炸焊接过程,得到典型的波状界面和射流现象,理论分析与模拟结果具有一致性。

      (2)在两种药量比、均匀布药方案中,忽略边界效应后的结果显示,沿着爆轰方向基板与复板之间的碰撞压力均逐渐增大,而且药量比越大,碰撞压力越大,界面波形越大。

      (3)在梯形布药方案中,通过改变炸药起爆端和末端的高度设计了4种方案,结果显示梯形布药可以有效消除爆炸焊接界面波的不均匀现象,使界面波形尺寸基本保持一致,而且节省了炸药用量。当起爆端和末端的高度分别为67.2 mm和42.0 mm时效果最佳。

      (4)介绍了几种波形成机理,展示并研究了界面波形成过程,模拟的波形成过程与复板流侵彻机理的一致性较好,证明了复板流侵彻机理解释界面波形成过程的有效性。

参考文献 (34)

目录

    /

    返回文章
    返回