-
铜和铝具有良好的导电性和导热性,在电气和电子工业中应用广泛. 将铜/铝及其合金连接形成可靠的复合结构,既可以实现结构轻量化,又能降低成本. 铝和铜的熔点、线膨胀系数等性能存在差异,使其焊接较为困难. 有学者研究铜铝异种金属的熔化焊、固相焊及超声焊接,发现存在金属间化合物和氧化颗粒残留等缺陷[1-4]. 磁脉冲半固态辅助钎焊工艺[5]结合电磁成形、半固态成形和钎焊的优点,利用洛伦兹力驱动母材挤压半固态钎料使其产生强烈的剪切流动,破坏并去除氧化膜,细化钎缝组织[6],提高焊接接头的综合性能.王振东等人[7]研究了Zn-15Al钎料的固相率对钎焊接头的影响,结果发现,固相率为60%的界面结合情况较好. 黄海川等人[8]研究了放电电压对钎焊的影响,研究表明,第二次放电电压为7 kV时能实现良好的冶金结合,且铝侧形成α-Al和CuZn5金属间化合物,钎料层出现α-Al、富锌相及CuZn5,铜侧形成锯齿状扩散层. 高远等人[6]研究温度条件对铜/铝磁脉冲半固态钎焊的影响,结果发现,温度升高会导致α相形状变得不规则,并且会有粗大树枝状组织产生的倾向,同时钎缝中铜含量升高. 邓凌波等人[9]研究了铜/铝管磁脉冲半固态钎焊钎料的流变行为,发现带锥角集磁器能有效细化铜/铝接头搭接区域中部钎缝组织, 并最终提高钎焊接头的抗剪强度.
目前关于铜/铝磁脉冲钎焊的研究侧重于工艺优化、钎料优化、钎缝组织及连接强度的研究,而分子动力学(molecular dynamics,MD)在研究结构演化和微观扩散机制方面有一定优势[10]. Zhang等人[11]借助MD分别研究了6061铝合金与AZ31B镁合金爆炸焊中冲击速度和冲击角度对扩散系数的影响,结果表明,碰撞角保持不变时,扩散系数随着碰撞速度的增大呈线性增大;碰撞速度保持440 m/s不变时,扩散系数随着碰撞角的减小呈线性减小. Chen等人[12]利用MD研究铜/铝爆炸焊,结果表明,当未加横向速度时,扩散系数与纵向速度成正比. 当纵向速度固定时,扩散系数与横向速度的平方成正比. Zhang等人[13]利用MD研究不锈钢和纯镍的扩散焊过程中原子的扩散行为,结果表明,表面粗糙度对扩散行为有显著影响,并且扩散焊的压力效应小于温度效应. Mao等人[14]利用MD研究铜/铝固液界面原子的扩散行为,研究表明,Cu原子的扩散深度与体系温度呈线性关系,扩散层厚度与扩散焊时间呈抛物线关系. Yang,Ma和Samanta等人[10, 15-16]利用MD研究铜/铝超声波焊接中焊缝的动态演变,以及温度、压力和超声波速度对铜/铝扩散的影响,结果表明,在超声波焊接中铜/铝出现了不对称扩散,扩散层厚度随焊接时间的延长而增加.袁晓静等人[17]采用MD对微弧等离子增材制造NiCr合金构件生长过程中温度场变化及等轴晶生长过程进行模拟,研究表明,冷却速率为3.38 K/ps和0.675 K/ps时,Ni-Cr体系呈现非晶凝固,冷却速率0.0775 K/ps下,Ni-Cr体系自发形核长大.
通过建立的铜−钎料−铝体系的MD模型,模拟铜/铝磁脉冲半固态钎焊过程原子的扩散行为,研究了不同冲击速度对冶金结合扩散层厚度的影响,找到了模拟扩散层厚度与试验扩散层厚度之间的关系.
-
采用的试验材料为T2紫铜、Zn15Al和1060铝合金,其中1060铝合金为退火态,试验材料化学成分如表1所示. Zn-15Al钎料在固相率为60%时可以实现较好的焊接效果[7],故焊接试验在410 ℃(对应固相率为60%)下进行.
表 1 试验材料的化学成分(质量分数,%)
Table 1. Chemical compositions of experimental materials
材料 Al Zn Cu Fe Ag Si Mg Mn Ni P Pb T2 — 0.003 7 余量 0.004 2 0.051 — — — 0.003 8 0.003 9 0.0032 1060 余量 0.05 0.05 0.35 — 0.2 — 0 — — — Zn15Al 14.95 余量 0.521 0.012 — 0.031 0.01 — — — 0.002 放电设备使用WG-IV60kJ型磁脉冲成形机,最大充电电压为11 kV,试验电容为550 μF. 试验工装主要由磁脉冲成形线圈、集磁器、感应加热线圈、热电偶、外管、钎料和内管组成,如图1所示. 工艺过程有两次放电过程(图2),一是初装定位管件后,在常温条件下通过脉冲放电产生的电磁力对内、外管和钎料进行装夹;二是设备充电后再通过感应线圈加热,当钎料温度达到410 ℃时,停止加热并放电,在室温下冷却,获得铜/铝管接头试样.
-
MD模型使用LAMMPS软件进行模拟. 模拟模型由铝基体、Zn15Al钎料和铜基体组成,如图3所示. 构建3个不同的区域,左中右3层(z方向)分别放置铝原子、钎料原子和铜原子. 铝和铜都是体心立方晶格,铝和铜的晶格常数分别为4.05 Å和3.67 Å;钎料层以锌的密排六方构型为基础随机将33%的Zn原子替换成Al原子,晶格常数为2.67 Å,模拟Zn15Al钎料. 模拟体系x,y和z方向分别对应[100],[010]和[001]晶向,尺寸为2.2 nm × 2.2 nm × 34 nm,其中铝、钎料和铜区域的长度均为10 nm,铝、钎料和铜区域创建的原子数分别为3 025,4 725个和4 725个.3个区域放在同一个模拟盒子,其中铝和钎料、钎料和铜的冲击接触面之间距离均为2 nm.x和y方向设置为周期边界条件;在冲击方向(z方向)上,铝区域的最左端和铜区域的最右端各取10 Å 距离的原子作为固定边界,固定边界原子在整个模拟过程中不受力的作用. 模拟选用的势函数是采用Dickel等人[18]拟合的MgAlZn的meam势函数、Liu等人[19]开发的CuAl的EAM势函数和苏文辉等人[20]开发拟合的CuZn的Morse势函数参数. 系统温度为300 K,时间为30 ps的NVT系综进行弛豫,使模拟体系初始化. 整个模拟体系是金属体系,所以选择0.001 ps作为体系的时间步长. 升温后在模拟冲击钎焊过程中,铜侧边界固定以防止模拟中体系发生移动,在铝侧边界向右施加不同的速度,模拟不同电压下的冲击钎焊.
初始化后使用NPT系综,将x和y方向上的压强设置为1.01 × 105 Pa,压力调节阻尼系数(pdamp系数)设置为0.1,可以更快的稳定体系压强. 将体系的温度从300 K升高至770 K(此时钎料模拟体系的固相率为60%),对铝侧边界原子施加不同的冲击速度,使用NVE系综来计算体系的运动,冲击时间为100 ps,去除铝侧边界原子的运动速度,并继续使用NVE系综模拟冲击后的钎焊扩散混合过程. 当NVE系综的温度达到稳定则表示混合过程结束,混合模拟时间为1 000 ps. 钎焊过程中各原子的运动情况均可在OVITO软件中观察到,从而揭示了钎焊界面各原子混合的机理,为在试验中得到更好的钎焊接头提供理论指导.
-
图4为无冲击速度情况下铝基体/钎料界面扩散过程中运动最远的原子位置变化曲线. 在没有冲击的情况下,扩散界面中Al,Zn和Cu原子的移动驱动力理论上差距不大,但实际上各个原子的移动表现并不同. 从图4可以看出,初始阶段原子移动很快,随着体系逐渐平衡,移动速度降低,最后在一个位置震动. 图中铝基体原子运动的最远距离为0.6 nm,而钎料中Al原子和Zn原子运动最远距离分别为1.6和1.4 nm. 钎料中Zn原子和Al原子的运动距离远大于铝基体中的Al原子,铝基体和钎料中的Zn原子和Al原子随着边界上原子排列无序区域的增加而扩散混合,说明钎料的无序化程度比铝基体无序化程度更大.
-
图5为不同冲击速度下在NVE系综仿真进行到100 ps时原子位形图. 从图5a可以看出,铝基体与钎料的结合界面有明显的相互扩散混合,而铜基体与钎料的结合界面未发现明显的原子相互扩散. 随着冲击速度的增加,钎料和铝基体原子进一步相互扩散,铝基体和钎料之间的扩散层厚度也随之增加;铜基体与钎料之间也发生轻微的相互扩散,在冲击速度100 m/s下,模拟扩散层最厚为0.6 nm左右. 同时发现,铝基体原子与钎料原子发生混合时,钎料中的Zn原子和Al原子向无序排列的铝基体中的运动趋势更加明显. 而铝基体原子进入扩散层后向钎料中的扩散较为困难,铝基体原子的运动的最远位置离平衡位置不到1 nm,运动距离非常短.
图 5 不同冲击速度下原子位形图
Figure 5. Shape and position of atoms under different impact velocity. (a) v = 0 m/s; (b) v = 50 m/s; (c) v = 100 m/s
为了更清楚地了解界面扩散行为,采用NVE系综模拟了冲击速度100 m/s下模拟体系中各原子随时间的位形变化. 图6为模拟体系各原子在0,500 ps和1 000 ps时刻的原子位形图. 开始阶段,钎料中仅有少量的Zn和Al原子向铝基体和铜基体扩散;随着扩散过程的进行,更多的Zn原子和Al原子从钎料中扩散至铜基体和铝基体. 在NVE系综模拟的最终时刻,铝基体和钎料界面形成了Zn和Al元素相互混合且较厚的扩散层,铜基体和钎料的界面处形成较薄的扩散层.
-
利用NVE系综分别模拟冲击速度分别为0,50,60,70,80 m/s和100 m/s,冲击时间为100 ps时,铝基体侧界面附近的原子含量情况,模拟计算到1 000 ps时的分布曲线如图7所示. 随着冲击速度的增大,界面处原子含量变化趋势逐渐缓和过渡,界面扩散层厚度也不断增加.
图 7 从铝基体至钎料界面扩散层原子含量分布
Figure 7. Atom content distribution in the diffusion layer from aluminum matrix to solder interface. (a) v = 0 m/s; (b) v = 50 m/s; (c) v = 60 m/s; (d) v = 70 m/s; (e)v = 80 m/s; (f) v = 100 m/s
图7中Zn原子和Al原子的含量分布曲线的交点更靠近钎料侧,说明铝基体原子和钎料含有的锌原子铝原子的相互混合主要发生在铝基体侧. 为了研究其原因,利用NVE系综模拟冲击速度80 m/s下不同时刻下原子含量分布,如图8所示. 开始时,交点靠近钎料一侧,进一步验证了钎料原子向铝基体扩散速度大于铝原子向钎料扩散速度. 但随着时间增加,交点位置向铝基体侧移动,直至达到中间位置附近,这是铝基体原子向钎料扩散比例不断提高的结果. 从图8a和图8b可以看出,Zn和Al含量变化的区间宽度有明显增大,而从图8b~图8e可以发现,Zn和Al原子扩散层厚度基本没有变化,这说明原子在模拟初期时扩散较快,随后达到稳定,扩散速度降低. 而铜侧则是陡变的曲线,可以判断出在冲击速度为80 m/s、冲击时间为100 ps和模拟时间为1000 ps时,铜侧未形成稳定的界面扩散层.
-
图9和图10分别为冲击速度分别为57,61,67,71 m/s和75 m/s焊接接头的微观形貌和能谱分析图. 能谱扫描点从钎料开始逐渐进入界面层和铝基体. 进入界面过渡层后,Al元素含量单调增加直到100%(即完全进入铝基体),而Zn元素含量单调递减直到0. 不同冲击速度57,61,67,71 m/s和75 m/s下原子含量梯度变化的扩散层厚度分别为3.5,3.8,4.5,5.1 μm和5.6 μm,说明随着冲击速度增加,扩散层厚度随之增加.
图 9 不同冲击速度下铝侧界面微观组织
Figure 9. Microstructure of aluminum side interface at different impact speeds. (a) v = 57 m/s; (b) v = 61 m/s; (c) v = 67 m/s; (d) v = 71 m/s; (e) v = 75 m/s
图 10 不同冲击速度下铝侧界面EDS曲线
Figure 10. EDS curves of aluminum side interface at different impact speeds. (a) v = 57 m/s; (b) v = 61 m/s; (c) v = 67 m/s; (d) v = 71 m/s; (e) v = 75 m/s
图11和图12分别为不同放电电压下钎焊时铜基体与钎料结合界面处的微观形貌和EDS曲线,能谱扫描点从铜基体开始逐渐进入界面层和钎料. Cu元素含量从100%降到60%左右,Zn元素从0增加到20%左右,Al元素含量从0增加到30%左右. 不同放电电压0,4.5,5.0,5.5 kV和6.0 kV下扩散层厚度分别为1.9,1.8,2.0 μm和2.5 μm. 在一定的速度范围内,钎料和铜基体之间的相互扩散受冲击速度影响不大. 从EDS曲线可以看到,在界面扩散层中,铜含量的分布规律与铝、锌有区别,这与铜、锌在扩散混合过程中形成了铜锌化合物有关.
图 11 不同放电电压下铜侧界面微观形貌
Figure 11. Microstructure of the aluminum side interface at different discharge voltages. (a) U = 4.0 kV; (b) U = 4.5 kV; (c) U = 5.0 kV; (d) U = 5.5 kV; (e) U = 6.0 kV
图 12 不同放电电压下铜侧界面EDS曲线
Figure 12. EDS curves of copper side interface at different discharge voltages. (a) U = 4.0 kV; (b) U = 4.5 kV; (c) U = 5.0 kV; (d) U = 5.5 kV; (e) U = 6.0 kV
为了获得模拟结果与试验结果之间的关系,对模拟得到的界面扩散层厚度和试验得到的界面扩散层厚度一一对应. 由于模拟时铜侧未形成明显的界面扩散层,这里主要分析钎料和铝基体侧的界面扩散层厚度变化. 图13为铝侧界面扩散层厚度和冲击速度的关系.图13中显示扩散层厚度随冲击速度的变化接近一条直线,说明界面扩散层厚度与冲击速度成线性关系,在冲击速度为57,61,67,71 m/s和75 m/s下模拟所得的界面扩散层厚度分别为1.475,1.575,1.725,1.825 μm和1.925 μm. 用上述模拟计算出的铝基体侧界面扩散层厚度与在相同冲击速度下试验得到的铝基体侧界面扩散层厚度作出对应的关系. 从图14可以看出,模拟的铝侧界面扩散层厚度与试验得到的铝侧界面扩散层厚度可拟合成线性关系. 若将不同冲击速度下模拟得到的界面扩散层厚度代入此曲线即可预测试验扩散层厚度,如表2所示,试验与模拟预测结果最大误差为2.8%,说明用模拟方法能较好地预测试验扩散层厚度.
图 13 不同冲击速度下铝侧界面模拟扩散层厚度
Figure 13. Simulated diffusion layer thickness of interface at aluminum side under different impact speeds
图 14 模拟扩散层厚度与试验扩散层厚度的关系
Figure 14. Relationship between the simulated thickness of the diffusion layer and the measured thickness in the experiment
表 2 不同冲击速度下界面模拟和试验扩散层厚度
Table 2. Simulation and experimental diffusion layer thickness of interface under different impact speeds
冲击速度
v/(m·s−1)扩散层厚度δ/μm 误差e(%) 试验值 模拟值 57 3.5 3.43 2.0 61 3.8 3.91 2.8 67 4.5 4.62 2.6 71 5.1 5.12 0.1 75 5.6 5.58 0.3 -
(1) 在温度为770 K、无冲击的情况下,钎料中的原子进入铝基体侧的界面扩散层后扩散容易;而铝基体原子进入界面扩散层后向钎料中扩散较为困难.
(2) 随着冲击速度从57 ~ 75 m/s不断增大,铝侧模拟和试验扩散层厚度不断增加;而铜侧模拟扩散层厚度增加不明显,试验中铜侧界面扩散层厚度稳定在 2 μm左右,说明冲击速度对其影响相对较小,界面上无序排列的原子较少,故界面扩散层厚度较小.
(3)根据模拟的情况,建立了在相同放电电压下模拟的界面扩散层厚度和试验界面扩散层厚度之间的关系,由模拟计算出的界面扩散层厚度与试验界面扩散层厚度相比较,最大误差为2.8%.
Interfacial diffusion process of Cu/Al magnetic pulse semi-solid assisted brazing
-
摘要: 结合分子动力学模拟和试验研究,对铜/铝管磁脉冲半固态辅助钎焊界面原子扩散过程进行了研究. 结果表明,在铝侧扩散界面原子主要在界面的无序原子层中相互扩散,且各元素的扩散行为不均匀,铝基体原子向钎料的扩散速度远小于钎料原子向铝基体扩散的速度;铜侧界面在模拟设置的冲击速度下,扩散层很薄且厚度变化并不明显. 测量不同冲击速度下铝侧界面的扩散层厚度发现,随着冲击速度增加,模拟的扩散层厚度呈线性增加,与试验结果相符. 根据模拟与试验结果建立了在相同或相近的冲击速度下,模拟界面扩散层厚度与试验界面扩散层厚度之间的关系,模拟结果能够较好地预测试验界面扩散层厚度, 最大误差为2.8%.Abstract: The atomic diffusion process of the magnetic pulse semi-solid assisted brazing interface of copper/aluminum tube was studied by molecular dynamics simulation and experiment. The results show that the atoms at the aluminum side diffusion interface mainly diffused with each other in the disordered atomic layer at the interface, and the diffusion behavior of the elements was not uniform. The diffusion speed of the aluminum matrix atoms to the filler metal was much lower than that of the filler metal atoms to the aluminum matrix. At the copper side interface, the diffusion layer is very thin (about two atomic layers thickness) and the thickness changes are not obvious at the simulated impact speeds. It is found that the simulated diffusion layer thickness increases linearly with the increase of impact speed, which is consistent with the experimental results. According to the simulation and experimental results, the relationship between the simulated interface diffusion layer thickness and the test interface diffusion layer thickness is established under the same or similar impact velocity. The simulation results can well predict the test interface diffusion layer thickness, with a maximum error rate of 2.8%.
-
Key words:
- Cu/Al tube /
- brazing /
- magnetic pulse /
- semi-solid /
- molecular dynamics
-
表 1 试验材料的化学成分(质量分数,%)
Table 1. Chemical compositions of experimental materials
材料 Al Zn Cu Fe Ag Si Mg Mn Ni P Pb T2 — 0.003 7 余量 0.004 2 0.051 — — — 0.003 8 0.003 9 0.0032 1060 余量 0.05 0.05 0.35 — 0.2 — 0 — — — Zn15Al 14.95 余量 0.521 0.012 — 0.031 0.01 — — — 0.002 表 2 不同冲击速度下界面模拟和试验扩散层厚度
Table 2. Simulation and experimental diffusion layer thickness of interface under different impact speeds
冲击速度
v/(m·s−1)扩散层厚度δ/μm 误差e(%) 试验值 模拟值 57 3.5 3.43 2.0 61 3.8 3.91 2.8 67 4.5 4.62 2.6 71 5.1 5.12 0.1 75 5.6 5.58 0.3 -
[1] Mubiayi M P, Akinlabi E. Friction stir welding of dissimilar materials between aluminium alloys and copper−An overview[C]//Lecture Notes in Engineering & Computer Science, London,UK, 2013, 1990 − 1996. [2] Bergmann J P, Petzoldt F, Schürer R, et al. Solid-state welding of aluminum to copper − case studies[J]. Welding in the World, 2013, 57(4): 541 − 550. doi: 10.1007/s40194-013-0049-z [3] Xu H, Liu C, Silberschmidt V V, et al. Behavior of aluminum oxide, intermetallics and voids in Cu-Al wire bonds[J]. Acta Materialia, 2011, 59(14): 5661 − 5673. doi: 10.1016/j.actamat.2011.05.041 [4] Wang C, Liu S, Zhu H, et al. Effect of process parameters on interfacial microstructure and mechanical properties of Al/Cu friction stir lap welding joints[J]. China Welding, 2022, 31(4): 48 − 58. [5] 黄尚宇, 黄海川, 雷雨, 等. 一种半固态无钎剂辅助的异种金属钎焊装置及方法: 中国, CN108941828B[P]. 2019-11-26. Huang Shangyu, Huang Haichuan, Lei Yu, et al. Semi-solid brazing device and method for dissimilar metals without the aid of brazing flux: China, CN108941828B[P]. 2019-11-26. [6] 高远, 黄尚宇, 邓凌波, 等. 温度条件对铜管/铝管磁脉冲−半固态复合辅助钎焊接头微观组织的影响[J]. 锻压技术, 2019, 44(10): 169 − 175. doi: 10.13330/j.issn.1000-3940.2019.10.030 Gao Yuan, Huang Shangyu, Deng Lingbo, et al. Influence of temperature condition on microstructure of Cu tube/Al tube brazed joint produced by magnetic pulse and semi-solid composite auxiliary technology[J]. Forging & Stamping Technology, 2019, 44(10): 169 − 175. doi: 10.13330/j.issn.1000-3940.2019.10.030 [7] 王振东, 黄尚宇, 李佳琪, 等. Zn-Al钎料固相率及组分对Cu/Al管磁脉冲−半固态复合辅助钎焊接头质量的影响初探[J]. 材料科学与工艺, 2020, 28(4): 1 − 7. doi: 10.11951/j.issn.1005-0299.20190061 Wang Zhendong, Huang Shangyu, Li Jiaqi, et al. Preliminary study on the influence of solder solid phase rate and composition of Zn-Al filler metal on the quality of magnetic pulse-semisolid hybrid assisted soldering of Cu/Al tubes[J]. Materials Science and Technology, 2020, 28(4): 1 − 7. doi: 10.11951/j.issn.1005-0299.20190061 [8] 黄海川, 黄尚宇, 李清宁, 等. 放电电压对Cu/Al管磁脉冲-半固态复合辅助钎焊质量的影响[J]. 塑性工程学报, 2020, 27(2): 60 − 67. doi: 10.3969/j.issn.1007-2012.2020.02.008 Huang Haichuan, Huang Shangyu, Li Qingning, et al. Effect of discharge voltage on quality of Cu/Al tube magnetic pulse-semisolid hybrid assisted brazing[J]. Journal of Plasticity Engineering, 2020, 27(2): 60 − 67. doi: 10.3969/j.issn.1007-2012.2020.02.008 [9] 邓凌波, 黄尚宇, 杨梅, 等. 集磁器结构对磁脉冲辅助半固态钎焊接头组织和性能的影响[J]. 焊接学报, 2022, 43(1): 48 − 54. Deng Lingbo, Huang Shangyu, Yang Mei, et al. Effect of magnetic collector structure on microstructure and properties of semi-solid brazing joint assisted by magnetic pulse[J]. Transactions of the China Welding Institution, 2022, 43(1): 48 − 54. [10] Yang J, Zhang J, Qiao J. Molecular dynamics simulations of atomic diffusion during the Al-Cu ultrasonic welding process[J]. Materials, 2019, 12(14): 2306. doi: 10.3390/ma12142306 [11] Zhang Tingting, Wang Wenxian, Zhou Jun, et al. Molecular dynamics simulations and experimental investigations of atomic diffusion behavior at bonding interface in an explosively welded Al/Mg alloy composite plate[J]. Acta Metallurgica Sinica (English Letters), 2017, 30(10): 983 − 991. doi: 10.1007/s40195-017-0628-x [12] Chen S Y, Wu Z W, Liu K X, et al. Atomic diffusion behavior in Cu-Al explosive welding process[J]. Journal of Applied Physics, 2013, 113(4): 044901. doi: 10.1063/1.4775788 [13] Zhang Yanqiu, Jiang Shuyong. Atomistic investigation on diffusion welding between stainless steel and pure Ni based on molecular dynamics simulation[J]. Materials, 2018, 11(10): 1957. doi: 10.3390/ma11101957 [14] Mao Aixia, Zhang Jinping, Yao Shichang, et al. The diffusion behaviors at the Cu-Al solid-liquid interface: A molecular dynamics study[J]. Results in Physics, 2020, 16: 102998. doi: 10.1016/j.rinp.2020.102998 [15] Ma Qiuchen, Song Cheng, Zhou Jianli, et al. Dynamic weld evolution during ultrasonic welding of Cu-Al joints[J]. Materials Science and Engineering: A, 2021, 823: 141724. doi: 10.1016/j.msea.2021.141724 [16] Samanta A, Xiao S, Shen N, et al. Atomistic simulation of diffusion bonding of dissimilar materials undergoing ultrasonic welding[J]. The International Journal of Advanced Manufacturing Technology, 2019, 103: 879 − 890. doi: 10.1007/s00170-019-03582-9 [17] 袁晓静, 郭晓辉, 关宁, 等. 微弧等离子增材制造NiCr合金的分子动力学数值模拟[J]. 焊接学报, 2021, 42(8): 25 − 32. doi: 10.12073/j.hjxb.20210131001 Yuan Xiaojing, Guo Xiaohui, Guan Ning, et al. Numerical simulation of molecular dynamics for fabrication of NiCr alloy by microarc plasma additive[J]. Transactions of the China Welding Institution, 2021, 42(8): 25 − 32. doi: 10.12073/j.hjxb.20210131001 [18] Dickel D E, Baskes M I, Aslam I, et al. New interatomic potential for Mg-Al-Zn alloys with specific application to dilute Mg-based alloys[J]. Modelling and Simulation in Materials Science and Engineering, 2018, 26(4): 045010. doi: 10.1088/1361-651X/aabaad [19] Liu X Y, Liu C L, Borucki L J. A new investigation of copper's role in enhancing Al-Cu interconnect electromigration resistance from an atomistic view[J]. Acta Materialia, 1999, 47(11): 3227 − 3231. doi: 10.1016/S1359-6454(99)00186-X [20] 苏文辉, 刘春林, 刘承蕙. 锌及 β 黄铜的原子间相互作用势函数与高压状态方程[J]. 吉林大学自然科学学报, 1980(1): 69 − 79. Su Wenhui, Liu Chunlin, Liu Chenghui. The potential functions of interatomic interaction and the high-pressure equations of state for zinc and β-brass[J]. Journal of Natural Science of Jilin University, 1980(1): 69 − 79. -