0.   序言

激光焊接由于其焊接速度快、热输入量低、焊接变形小等优势，已经广泛应用于汽车制造、航空航天、能源装备及船舶制造等工业生产焊接制造过程中. 激光焊接过程涉及到材料对激光能量的吸收过程、材料的熔化与蒸发、等离子体(羽烟)的形成、熔池与匙孔的形成、凝固与相变等复杂的多物理场问题，建立合理的反映激光焊接热过程特点的数值模型，定量描述和澄清激光焊接过程中的熔化、蒸发、凝固等物理现象，对于更好地理解和掌握激光热过程对焊缝成形、组织及应力演化的影响规律及内在机理十分重要，故数值模拟方法是焊接关键问题科学研究的强有力分析手段^[1].

激光焊接自身快速熔凝的特点导致其焊接过程中更易产生驼峰、气孔、飞溅、侧壁未熔合等焊接缺陷^[2-5]，直接影响激光焊缝成形质量，另外激光焊缝组织梯度大及残余应力等问题也成为影响激光焊接构件服役性能的重要因素.

文中从影响激光焊缝成形质量及服役性能的焊接缺陷、组织演变、残余应力几个方面入手，阐述数值模拟技术在激光焊接基础问题方面应用及发展.

1.   激光焊缝成形问题的数值模拟

激光焊缝成形过程及相关缺陷的产生与熔池的流动行为、匙孔的动态波动及熔池的凝固过程密切相关. 此部分主要是针对熔池流动及匙孔行为进行阐述.

1.1   焊缝表面成形

在激光高速焊接时，焊缝表面成形主要涉及驼峰、咬边、未熔合、焊穿等问题，这与熔池流动及匙孔行为密切相关. 在表面驼峰成形方面，Otto等人^[6]模拟研究了激光焊接表面驼峰的形成过程，发现金属蒸气的剪切力导致熔化金属向后加速，形成的熔体支流沿侧壁向后铺展，两侧支流汇聚形成驼峰.

在激光熔透焊接时，根部焊缝成形控制非常重要. Zhang等人^[7]利用自行研发的三维激光焊接模型，利用锐利界面方法处理了影响匙孔动态行为的关键物理因素，发现小孔根部的流体受到强烈的蒸气反作用力导致根部驼峰的形成，如图1所示. Bachmann等人^[8-9]研究表明激光熔透焊接过程中电磁力可以改变熔池的流动方向，进而阻止熔池下坠，利于熔透焊缝根部成形.

图  1  背部驼峰形成时的小孔行为

Figure  1.  Behavior of keyhole when the root humping is formed

下载: 全尺寸图片 幻灯片

在激光焊接咬边、焊穿等成形方面，Cho等人^[10]发现增加激光束震动频率改变了能量的传输及熔池的流动行为，从而明显改善焊缝成形质量，咬边和焊穿情况也得到明显改善，如图2所示.

图  2  不同震动频率下的激光焊缝成形情况

Figure  2.  Laser weld formation under different vibration frequency

下载: 全尺寸图片 幻灯片

1.2   焊接气孔

在焊缝气孔缺陷形成方面，许多学者分别在气泡形成、影响因素及气孔的演化等方面做了大量的数值模拟工作来研究气孔形成过程及其抑制机理. Zhao等人^[11]认为当匙孔深度突然变大时容易造成匙孔坍塌，气体进入熔池形成气泡. Pang等人^[12-13]自行开发了考虑匙孔、金属蒸气羽烟及熔池动态行为的三维多相自洽模型，发现蒸气羽烟摆动频率与匙孔波动频率近乎一致，阐明了蒸气羽烟形成的负压导致了匙孔后壁形成涡流，从而诱发环境气体卷入匙孔，促进了气泡的形成，如图3所示.

图  3  匙孔内蒸气羽烟行为及压力分布

Figure  3.  Evolution of vapor plume velocity inside the transient keyhole and relative pressure

下载: 全尺寸图片 幻灯片

气泡能否演化为气孔涉及到熔池流动行为及凝固过程，Lu等人^[14-15]模拟发现气泡迁移过程中存在消失、合并等现象，且被凝固界面捕获后才会演变成气孔，如图4所示. 另外低沸点元素蒸发对匙孔动态行为有着重要的影响，Huang等人^[16-17]模拟考虑了Mg元素含量对激光焊接铝合金气孔生成的影响，认为Mg元素含量增加使得蒸气反作用力增加，从而使得匙孔更加稳定，降低了匙孔型气孔的生成. 激光震荡在消除气孔方面的作用机理也得到了模拟澄清^[18-19]，其根本原因除了激光震荡过程可以提高气泡逃逸速度外，模拟还发现气泡在震荡过程中被匙孔合并导致气泡消失，进而消除气孔的机理.

图  4  气泡向气孔的演变过程

Figure  4.  Evolution process from bubble to pore. (a) 0.20 s; (b) 0.25 s; (c) 0.30 s; (d) 0.35 s; (e) 0.40 s; (f) 0.50 s

下载: 全尺寸图片 幻灯片

1.3   焊接飞溅

飞溅的产生会带来金属的元素损失，并影响焊缝成形质量. Hugger等人^[20]模拟了钢、铝激光焊接过程，发现熔体首先在匙孔的边缘产生驼峰，在表面张力作用下颈缩，在冲击力的作用下离开表面形成飞溅，如图5所示. Chang等人^[21]研究认为铝合金激光焊接时飞溅的产生与熔池的流动速度和涡流密切相关. Wu等人^[22]模拟发现匙孔后壁的熔池沿熔合线向上流动，导致上表面飞溅产生.由于飞溅产生的影响因素较为复杂，考虑金属蒸气行为对于理解飞溅产生的过程十分必要.

图  5  铝合金激光焊接飞溅产生过程

Figure  5.  Formation process of spatter

下载: 全尺寸图片 幻灯片

激光焊接过程中低沸点元素蒸发诱发飞溅产生的情况，也引起了学者的兴趣. Hao等人^[23-25]模拟研究了锌蒸气对匙孔与熔池动态行为的影响，分析了飞溅产生的过程，如图6所示，模拟还发现激光折线扫描路径下，由于锌蒸气分阶段进入匙孔，冲击力减小，使得匙孔和熔池更加稳定，明显降低了飞溅情况. Qi等人^[26]模拟环形激光热源作用下镀锌钢板的飞溅情况，认为产生飞溅的主要原因是熔池中的回流和涡流导致，同轴环形和高斯热源的匹配增加了回流的阻力，进而控制了飞溅的产生.

图  6  锌蒸气在匙孔中的流动行为

Figure  6.  Fluid flow of zinc vapor

下载: 全尺寸图片 幻灯片

2.   激光焊缝组织演化数值模拟

宏、微观组织的模拟也是激光焊接数值模拟的重要组成部分，近年来，各国学者已对激光焊缝组织的仿真进行了相关尝试，本节主要针对这方面的研究进展进行介绍.

2.1   宏观偏析

目前，针对焊缝宏观组织的数值模拟主要集中于异种材料因元素混合不均匀引起的宏观偏析方面. Gu等人^[27]利用元素分离法实现了不同材料在激光选区熔化过程中的界面相迁移行为的模拟. Yao等人^[28]利用该方法，并结合微尺度流体动力学和纳秒级别的热扩散分析，实现了316L不锈钢/Inconel 718异种材料激光熔池中鱼鳞状宏观偏析组织的模拟，如图7所示，这种宏观偏析组织的仿真重点在于界面参数的优化.

图  7  单道激光扫描焊缝截面元素分布模拟结果

Figure  7.  Alloy distribution for single-track laser scan. (a) experimental result of element distribution; (b) simulation result of element distribution

下载: 全尺寸图片 幻灯片

2.2   微观组织

焊缝微观组织的数值模拟的研究主要有以下3种方法：相场法、元胞自动机法以及蒙特卡罗法，且每种方法的侧重点不同.

Fallah等人^[29]采用相场法，针对不同条件下Ti-Nb合金激光熔池中的一次枝晶臂间距进行了模拟计算. Mi等人^[30]同样利用这一方法对Al-Cu合金激光焊熔池中的竞争生长行为进行了研究，发现在非最优取向的晶粒会影响到相邻晶粒的生长，晶界处的枝晶更易于在横向生长，通过阻碍其它方向枝晶的生长路径，进而改变不同取向晶粒的尺寸. Geng等人^[31-32]发现Al-Mg较Al-Cu合金会在凝固过程中更早地发生广泛的枝晶连接现象，从而使其裂纹敏感性更低，这可以归因于其一次枝晶较低的溶质偏析，如图8a，8b所示.

图  8  溶质分布及晶粒形貌模拟结果

Figure  8.  Simulation results of solute distribution and microscopic morphologies. (a) solute distribution in columnar growth of Al-4.0%Cu alloy; (b) solute distribution in columnar growth of Al-4.0%Mg alloy; (c) microscopic morphologies with different pre-heating temperatures in a single track

下载: 全尺寸图片 幻灯片

相较于相场法，元胞自动机法以其计算效率高、所需的模型参数少、试验匹配度较好的特点，得到了越来越多的关注. Ao等人^[33]利用该方法对激光熔池的微观组织进行了模拟，结果如图8c所示，可见等轴晶比例会随着预热温度的增加和激光扫描速度的降低而显著增加，而等轴晶形貌会随着相邻两道激光距离的增加变得更加长而窄. Shi等人^[34]利用该模型结合有限元的方法，实现了含气孔的不同柱状晶和等轴晶比例的熔池的仿真，并发现熔池晶粒生长形态取决于熔池的宽度、深度以及过冷度. 近期，Liu等人^[35]提出了结合三维元胞自动机法和一维相场法的新模型，通过结合元胞自动机法的高效性和相场法的准确性，成功实现了大尺度铝合金凝固前沿枝晶形貌的预测.

对于蒙特卡罗法，该模型并不考虑凝固过程中的枝晶形貌、元素偏聚和过冷等因素，目前被广泛地用于晶粒长大、再结晶等晶粒演化行为的数值模拟. Zhang等人^[36]利用该模型成功实现了12%Cr铁素体不锈钢在激光-电弧复合焊条件下热影响区晶粒的模拟，并发现顶部区域热影响区的晶粒尺寸要大于根部. Gleason等人^[37-38]利用该方法针对激光冲击焊接过程中1100铝合金和304不锈钢的不均匀的界面微观组织进行了仿真，如图9所示，实现了对不同屈服面、层错能以及晶界滑动条件下的晶粒伸长行为的预测.

图  9  有、无晶粒模型的铝/钢接触界面的温度与等效应变分布情况

Figure  9.  Comparison of the localized temperature and equivalent plastic strain sampled along the contact interface in both inhomogeneous (grain) and homogeneous (no grain) models

下载: 全尺寸图片 幻灯片

3.   焊接残余应力数值模拟

与传统的弧焊相比，激光焊对焊接残余应力和变形的影响主要体现在制造工艺方面^[39]，以下将对近几年来在激光焊接残余应力方面取得的进展进行介绍. Sun等人^[40]以Q235钢平板对接接头为研究对象，采用试验手段和热-弹-塑性有限元方法研究了激光焊和电弧焊的温度场、残余应力与焊接变形. 图10是激光焊(Case A)与电弧焊(Case C)得到的焊缝中央位置的温度循环曲线，可见激光焊的熔化面积小，而且熔深贯穿整个板厚，而电弧焊尽管熔化面积更大，激光焊的加热速度更快，峰值温度更高，高温停留时间更短而且冷却速度也更快.

图  10  激光焊与电弧焊焊缝中央位置温度循环曲线

Figure  10.  Thermal cycle at the weld center of laser and arc welding process

下载: 全尺寸图片 幻灯片

图11是薄板接头中央断面上的纵向残余应力与横向残余应力的数值模拟结果，可见在激光焊条件下接头的上、下表面的纵向残余应力的大小与分布几乎完全一致，这是因为激光焊产生的熔化区域在板厚方向分布较均匀所致，拉伸残余应力峰值与材料常温屈服强度基本一致. 在电弧焊条件下的横向残余应力的峰值要远高于激光焊，单就焊接残余应力的峰值大小而言，激光焊并不能有效减缓残余应力，该研究结果也得到了其它学者的支持^[41-42].

图  11  中央截面上的残余应力分布

Figure  11.  Residual stress distribution on the central section. (a) longitudinal residual stress; (b) transverse residual stress

下载: 全尺寸图片 幻灯片

在激光厚板焊接残余应力方面，Xu等人^[43]模拟比较了板厚为12 mm的Q460钢在激光焊和多层多道气体保护焊条件下接头的残余应力，发现激光焊接头的高拉伸纵向残余应力范围要明显小于气体保护焊接头，但峰值应力没有明显的差异. Yan等人^[44]研究了板厚为10 mm 的316L不锈钢在激光条件下的残余应力，如图12所示，由于板厚较厚，横向残余应力的大小与薄板接头相比有显著增加. Deng等人^[45]采用考虑加工硬化及退火软化的热-弹-塑性有限元方法研究了SUS304管-管激光焊对接接头的残余应力，并讨论了由固溶处理引起的初期残余应力对焊接残余应力的影响.

图  12  316L对接接头的焊接残余应力分布

Figure  12.  Residual stress distribution of welded joint. (a) mises stress; (b) transverse stress; (c) longitudinal stress

下载: 全尺寸图片 幻灯片

Elmesalamy等人^[46]以板厚为10和20 mm 316L钢对接接头为研究对象，采用轮廓法测量了多道窄间隙激光焊对接接头和多层多道TIG焊对接接头的残余应力. 图13是距焊接接头上表面1.5 mm位置的纵向残余应力分布，可见窄间隙激光焊和多层多道TIG焊产生的峰值应力前者为310 MPa，而后者为520 MPa，且激光焊纵向高拉伸残余应力区域明显较窄，故窄间隙激光焊方法可能是一个控制残余应力的有效途径，且多道窄间隙激光焊接方法将是用于控制加工硬化显著材料厚板接头焊接残余应力一种有效的方法.

图  13  距上表面1.5 mm处的纵向残余应力分布

Figure  13.  Longitudinal residual stress distribution at the location of 1.5 mm away from the upper surface

下载: 全尺寸图片 幻灯片

近年来，激光焊残余应力的报道在铝合金、镁合金、钛合金^[47]以及异种金属激光熔钎焊^[48]方面也越来越多.

4.   结束语与展望

(1) 在激光焊缝成形模拟方面，匙孔的动态行为决定了焊接过程的稳定性，现有数值模型已经实现了匙孔壁与激光能量的实时耦合，但是激光束与等离子体(羽烟)的耦合行为很少涉及，考虑金属蒸气流动行为与激光热源的相互作用，是激光焊接过程建模发展的方向.

(2) 在激光焊缝组织模拟方面，熔池的流动对元素宏观偏析行为方面的模拟结果取决于界面参数的选取与优化；而微观组织方面，进一步完善3种模型需要准确获得新晶粒形核的相关参数十分重要，实现宏观与微观组织多尺度模拟还具有很高的挑战性.

(3) 在激光焊接残余应力方面，对于低合金高强钢及超高强钢而言，如何建立完备的“热-组织-力学”多场耦合模型，高精度的模拟激光焊条件下残余应力将是一个挑战性的课题. 此外，对加工硬化较显著材料优化工艺降低残余应力、降低累积塑性应变和减缓敏化程度是值得深入探索的方向.