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

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

（1.安徽理工大学深部煤矿采动响应与灾害防控国家重点实验室，安徽 淮南 232001；2.安徽理工大学力学与光电物理学院，安徽 淮南 232001；3.安徽理工大学土木建筑学院，安徽 淮南 232001；4.安徽理工大学化学工程学院，安徽 淮南 232001；5.中国科学院材料力学行为和设计重点实验室，中国科学技术大学，安徽 合肥 230027）

由于爆炸焊接过程具有瞬时性、复杂性，人们在通过实验及相关设备来观察、研究爆炸焊接过程时，效果十分有限。随着计算机技术的发展，数值模拟技术的突破为爆炸焊接研究提供了帮助。迄今为止，国内外诸多学者利用数值模拟软件对爆炸焊接过程及参数开展了大量研究，取得了较理想的研究成果。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法进行二维建模，研究炸药与爆炸焊接界面波的关系及波的形成机理。

1 计算模型

1.1 爆炸几何模型

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 计算模型Fig.1 Computational model

1.2 材料模型和状态方程

硝铵炸药选用 Jones-Wilkins-Lee（JWL）[12]状态方程（EOS），计算公式为

式中：AJ、BJ、R1、R2和ω为材料系数；p为爆轰产物压力，GPa；E0为初始比内能，kJ·cm-3；V为爆轰气体产物的相对比容。炸药的具体JWL状态参数见表1，其中ρ为密度，D为爆速。

表1 硝铵炸药的JWL状态方程参数Table 1 JWL EOS parameters of ammonium nitrate explosive

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

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

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

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

式中：εp为有效塑性应变；为 有效塑性应变率，为参考应变率；A、B、C、m、n是与材料相关的常数；Tr为室温；Tm为熔点。Cu和Q235钢的Johnson-Cook材料模型参数列于表2，Grüneisen EOS参数列于表3，其中：c为声速，km/s；S1为拟合系数；a为Grüneisen系数Γ0的一阶体积修正，G为剪切模量。

表2 Cu和Q235钢的Johnson-Cook模型参数Table 2 Parameters of Johnson-Cook model of Cu and Q235 steel

表3 Cu和Q235钢的Grüneisen方程参数Table 3 Grüneisen EOS parameters of Cu and Q235 steel

1.3 爆炸复合窗口

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

图2 爆炸焊接窗口Fig.2 Explosive welding window

2 模拟结果分析

2.1 均匀布药

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

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

图3 均匀布药模拟效果图Fig.3 Simulation effect diagram of uniform charge

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

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

图4 关键点取样示意图Fig.4 Schematic diagram of key point sampling

表4 均匀布药方案关键点碰撞压力Table 4 Collision pressure of key points in uniform charge scheme

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

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

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

图5 关键点碰撞压力值折线Fig.5 Line diagram of collision pressure of key points

从图3可以看出，复板上部均出现了打伤现象，这是由于前期炸药和板材划分的网格尺寸不同所致。R1和R2两组模拟均得到了波状焊接界面，这被认为是焊接良好的标志[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是均匀布药下各个关键点的碰撞压力折线。考虑到边界效应的影响，去掉首尾两个关键点后，折线显示沿着爆轰方向焊接界面的碰撞压力逐渐增大，由于碰撞压力过大会产生熔化层，产生的稀疏波可能将已结合的复合板拉开，降低焊接质量，因此均匀布药不是最终选择。

2.2 梯形布药

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

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

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

图6 梯形装药结构示意图Fig.6 Schematic diagram of ladder charge structure

表5 梯形布药方案Table 5 Ladder charging scheme

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

图7 关键点碰撞压力值折线图Fig.7 Line diagram of collision pressure of key points

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

表6 梯形布药方案关键点碰撞压力Table 6 Collision pressure of key points of ladder charge scheme

3 界面波形成机理

爆炸焊接过程具有瞬时性，因此很难用实验方法来研究界面波的形成与焊接参数之间的关系。虽然有关界面波形成机理的研究已经进行了数十年，但迄今为止仍然没有达成一个统一的认识[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 波形成示意图Fig.8 Illustration of wave formation

4 结 论

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

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

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

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