基于CuSi3焊丝的镀锌钢板激光钎焊温度场及接头组织分析

基于CuSi3焊丝的镀锌钢板激光钎焊温度场及接头组织分析 Temperature Field Analysis and Joint Microstructure of Laser Brazing of Galvanized Steel Sheet Based on CuSi3 Welding Wire first-author 刘海江（1967—），男，博士，教授，主要从事复杂装备数字化设计与智能制造技术的研究。 刘海江（1967—），男，博士，教授，主要从事复杂装备数字化设计与智能制造技术的研究。 同济大学 机械与能源工程学院，上海 201804 同济大学 机械与能源工程学院，上海 201804 College of Mechanical Engineering, Tongji University, Shanghai 201804, China College of Mechanical Engineering, Tongji University, Shanghai 201804, China 采用焊接试验与焊接数值仿真相结合的方法对镀锌钢板CuSi3焊丝激光钎焊焊接温度场进行研究，使用高斯面热源和旋转高斯体热模型构成的组合热源模型对温度场数值模型进行计算，根据试验和数值计算结果分析激光钎焊过程中焊接线能量对接头显微组织及温度场分布的影响规律。焊接接头母材区主要为尺寸均匀的细密等轴晶，而焊接热影响区主要为垂直于温度场等温线并朝向热源方向生长的柱状晶，柱状晶的宽度随着靠近热源在不断生长；随着焊接线能量的提升，焊接接头截面的钎缝深宽比也展现出增大的规律；表面焊缝中心温度场分布梯度较小，焊接线能量对焊缝中心峰值温度影响不大；焊缝边缘贴近母材区域温度在四种试验参数下均高于锌镀层沸点，从而导致焊接试件表面产生气孔。 The welding temperature field of CuSi3 welding wire laser brazing of galvanized steel sheet was studied by combining welding experiment and welding numerical simulation. The combined heat source model composed of Gaussian surface heat source and rotating Gaussian body heat model was used to calculate the numerical model of the temperature field, and the effect of the welding line energy on the microstructure of the welded joint and the effect of temperature field distribution during the laser brazing process was analyzed according to the experimental and numerical calculation results. The experimental results show that the base metal BM (basic mental, BM) of the welded joint is mainly distributed as fine equiaxed grains with uniform size, while the HAZ (heat affect zone, HAZ) area is mainly distributed columnar grains that the growth direction is perpendicular to the temperature field isotherm and towards the heat source; with the increase of welding line energy, the aspect ratio of the welded joint also shows an increasing law; the temperature field distribution gradient in the center of the surface weld is small, and the welding line energy has little effect on the peak temperature of the weld center; the temperature of the weld edge close to the base metal is higher than the boiling point of the zinc coating under all parameters, which leads to the formation of holes on the surface of specimens. 镀锌钢板 CuSi3焊丝 激光钎焊 复合热源模型 温度场 显微组织 galvanized steel sheet CuSi3 welding wire laser brazing composite heat source model temperature field Microstructure 前言 1 激光钎焊以功率密度高和方向性优良的激光束作为热源来熔化钎料，使熔融金属液与母材完成润湿铺展，自然冷却后实现工件的连接。与其他传统焊接方式相比，激光钎焊具有焊接速度快、能量密度大、焊件变形小、焊缝深宽比大、易实现自动控制等优势，被广泛应用于高强钢和有镀层金属薄板的连接。由于激光钎焊不直接熔化母材，可解决激光深熔焊过程中锌蒸汽排除困难而导致焊缝气孔过多、成形不好的问题 ［ 1 1 - 2 2 1-2 ］ ，因此激光钎焊技术广泛应用于汽车白车身车顶的焊接。近些年来，铜基材料作为填丝金属大量使用于镀锌钢板的激光钎焊加工，而CuSi3焊丝凭借其优秀的机械性能（包括良好的流动性和刚度） ［ 3 3 ］ 被广泛采用，但激光钎焊过程中的温度梯度很大，温度场的分布将直接对CuSi3焊丝在镀锌钢板表面的熔化、湿润和铺展过程造成影响。因此分析研究基于CuSi3焊丝的镀锌钢板激光钎焊的温度场分布，对研究和调控焊接接头显微组织转变、焊后接头形变等关键参数有着重要的工程指导意义。 由于激光加热速率快，板材温度变化剧烈，采用传统的热电偶法难以测量靠近熔池区域的温度。有限元数值模拟是研究焊接温度场的一种高效方法，但激光钎焊涉及激光束的移动、焊丝的熔融铺展，物理过程比较复杂，因此相关的研究并不多见。Jeon ［ 4 4 - 5 5 4-5 ］ 等人基于ABAQUS软件，将激光作为表面热源，把焊接过程分为六个典型阶段，模拟了304不锈钢和5052铝合金板的激光钎焊温度场。赵振龙 ［ 6 6 ］ 采用Ansys Workbench软件，选用高斯移动热源模型，对汽车顶盖和侧围进行了激光钎焊温度场仿真，并与传统点焊对比得出了激光钎焊工艺的优越性。宋永刚 ［ 7 7 ］ 采用ANSYSY有限元分析软件，研究了焊接速度对激光钎焊CBN的影响。以上研究建立的激光钎焊轴对称模型相对简单，没有考虑钎料的动态填充过程。封小松 ［ 8 8 ］ 采用生死单元网格模拟钎料填充过程，分析了镀锌钢板单、双光束温度场的不同特点，但仅对激光相关参数进行了研究。 本文采用实验与数值仿真相结合，对0.7 mm镀锌钢板特殊结构件激光钎焊进行焊接和数值建模，从而模拟镀锌板在车顶焊接过程中的实际工程应用场景。在深入分析激光钎焊加热过程的基础上，应用ABAQUS有限元分析软件建立镀锌钢板搭接接头激光钎焊温度场模型，对焊接过程中钎缝、母材的温度变化及激光功率、焊接速度、光斑直径等焊接参数对温度场的影响规律进行深入研究，并以工艺实验数据为检验标准验证仿真模型的可靠性。 试验材料及方法 1 试验焊接母材是牌号为GMW2M-ST-SCR2HD60G60GE的汽车高强度钢，双面均有镀锌层，镀锌厚度为上表面65 g/m 2 、下表面66 g/m 2 ，母材化学成分如 表1 表1 所示；激光钎焊使用的填丝材料为CuSi3焊丝，其化学成分如 表2 表2 所示。板材尺寸为200 mm×75 mm×0.7 mm，采用两板搭接的形式进行激光钎焊，下部焊接试件为标准平板，上部为标准平板弯成90°直角，弯折边15 mm，翻边圆角半径为2 mm，具体搭接焊接情况如 图1 图1 所示。 C Si Mn P S Al Ti 0.02 0.01 0.015 0.01 0.007 0.04 0.05 Cu Si Zn Mn Fe Sn Al P Base 3.5 0.18 1.2 0.02 0.11 0.009 0.01 焊接光源采用德国通快公司生产的TruDisk-5006型号的碟片式激光器，使用威尔得WPC-600送丝机构进行钎焊送丝，采用KUKA O2ORB-100六轴机器人进行焊接操作，实际焊接工艺参数如 表1 表1 所示。焊后使用奥林巴斯BX53M正立金相光学显微镜观察焊缝形貌及微观组织。 试件 编号 激光功率/W 焊接速度 /（mm·s -1 ） 线能量 /（J·mm -1 ） 送丝速度 /（mm·s -1 ） 光斑直径/mm 离焦量 /mm S1 4 500 70 64.3 70 4.8 0 S2 4 100 60 68.3 65 4.8 0 S3 3 700 52 71.2 60 4.8 0 S4 3 200 43 74.4 43 4.8 0 温度场数值模型 1 几何模型定义和网格划分 2 在有限元仿真和分析环节，根据实际激光钎焊工艺实验情况简化焊接模型，离热源较远的母材几乎不会受到温度变化影响，为加快仿真计算速度在有限元分析软件ABAQUS中建立三维模型如下：平板尺寸100 mm×30 mm×0.7 mm；下部焊接试件为标准平板，上部焊接试件为标准平板弯成90°直角，弯折边长为15 mm，翻边圆角半径为2 mm。考虑到母材与焊丝之间的热效应问题，将模型分为母材基体和焊丝两部分。激光钎焊数值模型中母材镀锌钢板的密度7 840 kg/m 3 、焊丝CuSi3密度8 500 kg/m 3 。镀锌钢板和CuSi3焊丝的基本物理参数如 表4 表4 所示，热力学参数如 表5 表5 所示，室温为20 °C。 牌号 密度 /（kg·m -3 ） 固相线温度/°C 液相线温度/°C 潜热 /（J·kg -1 ） 热对流系数 /［W·（m 2 ·K） -1 ］ 镀锌钢板 7 840 1 025 1 300 225 000 33.5 CuSi3焊丝 8 500 970 1 025 207 000 47.0 温度/°C 导热系数/［W·（m·K） -1 ］ 比热容/［J·（g·K） -1 ］ 28 46.64 0.511 100 43.62 0.559 200 43.10 0.597 400 42.16 0.717 500 39.98 0.773 700 28.03 0.820 780 29.21 0.807 800 29.55 0.800 860 30.38 0.773 900 31.42 0.754 1 000 30.74 0.725 1 060 27.22 0.727 激光钎焊过程中，钎料随着激光热源的移动而熔化流入焊接区域，随着钎料的冷却形成焊缝，该过程对焊接的温度场有较大影响，故在仿真计算时采用生死单元来模拟钎焊过程，即在钎料依次熔化冷却形成焊缝的过程中逐步激活焊缝的单元。为增加钎焊仿真的精准性，将钎缝区域及热源经过路径附近的网格尺寸细化至0.14 mm；同时为节约计算时间并兼顾计算精度与效率，将网格尺寸随着距焊缝处距离的增加而逐渐增大至1 mm，细化后模型的网格分布如 图2 图2 所示。网格类型采用DC3D20型单元进行划分，共生成74 900个单元。 热源模型 2 在激光焊接过程中，镀锌板试件表面温度会因为高功率密度的激光而迅速升高，试件表面材料因此快速熔化、汽化产生蒸汽，从而导致汽化膨胀压力的产生进而出现小孔效应。为了反映出激光焊接过程中小孔效应的作用，体现热源在深度方向上的能量吸收，激光焊接的数值仿真热源模型主要使用三维热源模型或组合热源模型 ［ 10 10 ］ 。由于激光光斑的热输入分布与高斯分布函数相吻合，且激光钎焊过程中没有形成深的熔池，只是温度的骤升和骤降 ［ 11 11 ］ 且激光光源不受气体物质的影响 ［ 12 12 - 13 13 12-13 ］ 。又通过分析认为该激光钎焊模型是介于深熔焊和热导焊之间的焊接模型，所以本文采用旋转高斯体热源模型和平面高斯热源模型构成的组合热源模型。 平面高斯热源模型的数学表达式为 <math display="block" id="M1"><mi>q</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><mi>η</mi><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">π</mi><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac><msup><mrow><mi>e</mi></mrow><mrow><mo stretchy="false">(</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003483&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003479&type= 27.43199921 9.31333351 式中　 <math id="M2"><mi>q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003561&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003557&type= 1.86266661 3.38666677 为热源密度； <math id="M3"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003500&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003496&type= 3.55599999 4.48733330 为最大输入能量； <math id="M4"><mi>η</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003507&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003503&type= 1.86266661 3.38666677 为激光功率有效利用系数； <math id="M5"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003586&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003582&type= 2.62466669 3.72533321 为热源的有效半径； r 为焊接试件表面任意一点到光源中心的距离，其计算公式为 <math id="M6"><mi>r</mi><mo>=</mo><mroot><mrow><msup><mrow><mo stretchy="false">(</mo><mi>x</mi><mo>-</mo><mi>v</mi><mi>t</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003522&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003519&type= 26.41600037 5.33400011 ，其中 <math id="M7"><mi>v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003529&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003525&type= 1.69333339 2.62466669 为实际焊接速度。平面高斯热源模型示意如 图3 图3 所示。 旋转高斯体热源模型的数学表达式为 <math display="block" id="M8"><mi>q</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">9</mn><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">π</mi><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mi>H</mi><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">3</mn></mrow></msup><mo stretchy="false">)</mo></mrow></mfrac><msup><mrow><mi>e</mi></mrow><mrow><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mo>-</mo><mn mathvariant="normal">9</mn><mo stretchy="false">(</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo></mrow><mrow><msubsup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>H</mi><mo>/</mo><mi>z</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></mfenced></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003555&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003549&type= 49.78400040 10.07533360 式中　 <math id="M9"><mi>q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003561&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003557&type= 1.86266661 3.38666677 为热源密度； <math id="M10"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003571&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003567&type= 3.55599999 4.48733330 为最大输入能量； <math id="M11"><mi>H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003578&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003575&type= 3.13266683 2.62466669 为热流的有效作用深度； <math id="M12"><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003586&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003582&type= 2.62466669 3.72533321 为热源上表面半径。旋转高斯热源模型示意如 图4 图4 所示。 边界条件 2 激光钎焊的焊接过程为典型的非线性瞬态热传导问题，其过程方程为： <math display="block" id="M13"><mi>ρ</mi><mi>c</mi><mfrac><mrow><mo>∂</mo><mi>T</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><msub><mrow><mi>K</mi></mrow><mrow><mi>x</mi></mrow></msub><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>T</mi></mrow><mrow><mo>∂</mo><msup><mrow><mi>X</mi></mrow><mrow><msup><mrow/><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></msup></mrow></mfrac><mo>+</mo><msub><mrow><mi>K</mi></mrow><mrow><mi>y</mi></mrow></msub><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>T</mi></mrow><mrow><mo>∂</mo><msup><mrow><mi>Y</mi></mrow><mrow><msup><mrow/><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></msup></mrow></mfrac><mo>+</mo><msub><mrow><mi>K</mi></mrow><mrow><mi>z</mi></mrow></msub><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>T</mi></mrow><mrow><mo>∂</mo><msup><mrow><mi>Z</mi></mrow><mrow><msup><mrow/><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003602&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003599&type= 50.63066483 7.70466709 式中　 <math id="M14"><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003610&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003605&type= 2.03200006 3.38666677 为密度； <math id="M15"><mi>c</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003616&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003613&type= 1.69333339 2.62466669 为比热容； <math id="M16"><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003623&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003620&type= 2.45533323 2.62466669 为温度； <math id="M17"><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003629&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003626&type= 1.10066664 2.62466669 为时间； <math id="M18"><msub><mrow><mi>K</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003635&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003632&type= 3.64066648 3.72533321 、 <math id="M19"><msub><mrow><mi>K</mi></mrow><mrow><mi>y</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003645&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003640&type= 3.64066648 3.72533321 、 <math id="M20"><msub><mrow><mi>K</mi></mrow><mrow><mi>z</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003654&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003650&type= 3.55599999 3.72533321 为各方向的导热系数。 对整个镀锌板激光钎焊工艺参数实验数值模型而言，试件模型的热输入 <math id="M21"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003663&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003660&type= 3.55599999 4.48733330 由旋转高斯体热源和平面高斯热源组合热源模型所提供；热量耗散来自于试件整体与环境的热辐射耗散 <math id="M22"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003667&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003665&type= 4.57200003 4.48733330 、试件整体与环境的对流耗散 <math id="M23"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">v</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003671&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003669&type= 5.92666674 4.48733330 、焊接过程中的液-气相变潜热耗散 <math id="M24"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">v</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003675&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003673&type= 4.91066647 4.48733330 三部分。 数值模型的热辐射边界条件可表示为 <math display="block" id="M25"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi></mrow></msub><mo>=</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub><mi>ε</mi><mo stretchy="false">(</mo><msup><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>-</mo><msubsup><mrow><mi>T</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msubsup><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003680&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003677&type= 25.65399933 4.82600021 式中　 <math id="M26"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003684&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003682&type= 2.79399991 3.72533321 为玻尔兹曼常数； <math id="M27"><mi>ε</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003687&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003686&type= 1.52400005 2.62466669 为表面辐射系数； <math id="M28"><msub><mrow><mi>T</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003713&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003710&type= 3.64066648 3.72533321 为室温。 试件整体表面的对流耗散边界条件可表示为 <math display="block" id="M29"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">v</mi></mrow></msub><mo>=</mo><msub><mrow><mi>h</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo stretchy="false">(</mo><mi>T</mi><mo>-</mo><msub><mrow><mi>T</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003699&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003695&type= 24.29933357 4.48733330 式中　 <math id="M30"><msub><mrow><mi>h</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003706&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003702&type= 2.96333337 3.72533321 试件整体与室温环境的换热系数； <math id="M31"><msub><mrow><mi>T</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003713&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003710&type= 3.64066648 3.72533321 为室温。 母材和钎料因高功率激光表照射而汽化蒸发所造成的液-气相变潜热边界条件可表示为 <math display="block" id="M32"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">v</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">p</mi></mrow></msub><mo>=</mo><mi>W</mi><mi>L</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003718&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003716&type= 13.29266548 4.48733330 式中　 <math id="M33"><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003723&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003721&type= 3.55599999 2.62466669 为蒸发率； <math id="M34"><mi>L</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003727&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003725&type= 2.20133328 2.62466669 为材料的液-气相变潜热。 即该焊接数值模型的整体净热输入可表示为 <math display="block" id="M35"><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">n</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">t</mi></mrow></msub><mo>=</mo><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><mo>-</mo><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi></mrow></msub><mo>-</mo><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">v</mi></mrow></msub><mo>-</mo><msub><mrow><mi>q</mi></mrow><mrow><mi mathvariant="normal">v</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003731&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39003729&type= 36.40666580 4.48733330 结果与分析 1 显微组织 2 分别对四种不同焊接参数的焊接接头显微组织和钎缝形貌进行观察和记录。选取出接头截面中钎缝两侧区域进行分析，分析区域包含钎料、热影响区以及母材区，具体选取区域如 图5 图5 所示。 左侧观测区域显微组织如 图6 图6 所示。由图可知，近似线能量的焊接接头显微组织类型基本相同，但分布和比例随着焊接线能量的增加变化较为显著。显微组织基本由分布于母材区（BM）的等轴晶和靠近热源处的热影响区（HAZ）的柱状晶组成，且随着线能量的增加柱状晶的生长趋势愈发明显。由图 6a 6a 、 6b 6b 可知，在HAZ和BM的交接区域，显微组织分界明显：BM区域以细密的等轴晶为主，其平均尺寸不足10 μm；HAZ则大部分都已经经过高温转变为细长的柱状晶，并夹杂少量尺寸较为粗大的等轴晶，柱状晶平均长度约为35 μm。HAZ区域的柱状晶生长方向明显，垂直于温度场等温线并朝向热源，且柱状晶的宽度也随着靠近钎缝而增大。 1,2;3,4; 右侧观测区域的显微组织如 图7 图7 所示。右侧显微组织分布趋势与左侧大致相同，BM区域主要分布为细密的等轴晶，HAZ主要分布为柱状晶及部分粗大的等轴晶。但与左侧显微组织不同的是右侧HAZ区域夹杂的柱状晶含量减少、尺寸较大的等轴晶数量比例较左侧有所提高，该现象在焊接线能量较高时更加显著；且BM区与HAZ的熔合线明显程度随线能量的降低而降低。 1,2;3,4; 焊缝形状是影响钎焊焊缝质量的关键因素之一，本文引用熔焊中的深宽比概念对钎焊焊缝形状加以表征，来反映镀锌板钎焊过程中不同焊接线能量对钎焊焊缝深宽比的影响结果，具体钎缝宽度、深度取值如 图8 图8 所示。 在四组钎焊试验中，焊接线能量对钎缝深宽比的影响如 图9 图9 所示。四组焊接试件的线能量大小从64.3 J/mm增加至74.4 J/mm，其钎缝深宽比从0.273增大至0.394，反映出焊接热输入的增加对焊缝深宽比影响明显。 温度场计算结果 2 采用焊接接头热影响区分布情况验证仿真结果， 图10 图10 为S1试件的实际焊接热影响区分布与相同参数下热影响区形貌分布对照结果，对照显示分布结果吻合良好、分布基本一致，验证了该温度场有限元模型的准确性。 以焊接试件S1的温度场为例可得到整体数值模型焊接过程：钎料熔化形成的焊缝中心在受热约0.18 s后温度开始急剧上升，在0.22 s即可到达最高温度1 384 ℃，固态钎料随之熔融为液态金属流入母材间隙完成润湿铺展，此时钎料快速传热给母材，温度降至固相线以下后形成钎缝。随着激光源的离去、输入热源消失，钎料快速冷却。随后冷却速度减缓并一直持续到焊接结束，焊缝中心平均温度保持在200 ℃左右至少15 s后并最终降至室温。 焊接温度场分布计算结果如 图11 图11 所示，可以推断，镀锌板激光钎焊模型总体温度梯度不大，但随着焊接线能量的提升温度场的梯度随之提升，且形成的温度场等温线越宽。随着冷却过程的发生，焊接接头钎料部分的温度分布逐渐传导至镀锌板母材，由于钎料CuSi3的热传导系数小于母材，即其在降温过程中的热传导速度和对流换热速度都受制于钎料的影响，造成焊缝中心区域温度很高但温度梯度较小的情况。 1,2;3,4; 焊缝区域的焊接热循环曲线如 图12 图12 所示。由图可知，焊接实验过程中，焊接线能量对焊缝中心温度影响并不明显。不同的焊接线能量其焊件中心点P center 的最高温度略有不同，P center 温度峰值随着焊接功率的提升有一定增大，但中心最高温度基本分布在1 400 ℃附近（见 图12a） 图12a） 。 图12b 图12b 为焊缝边缘处贴近母材表面一点P edge 的焊接热循环曲线，其温度变化趋势与P center 的焊接热循环曲线趋势基本一致，可以看到边缘P edge 处试件S1~S4的峰值温度均在1 100 ℃左右，都已超过锌镀层的沸点908 ℃，因而钎料边缘处的镀锌层在激光钎焊过程中将被汽化烧损。 1,2; 结论与展望 1 本文对镀锌钢板激光钎焊的受热过程进行分析，并在此基础上建立有限元数值模型，结合激光钎焊实验对镀锌钢板CuSi3焊丝激光钎焊过程进行研究，得到以下结论： （1）采用旋转高斯体热源模型和平面高斯热源模型构成的组合热源模型对基于CuSi3焊丝的镀锌钢板激光钎焊过程中的温度场进行数值仿真，通过实验焊接接头热影响区的分布验证了数值仿真模型的准确性。 （2）焊接接头显微组织随着靠近钎料而逐渐由细密的晶粒转变为粗大的晶粒，BM区主要分布为尺寸不足10 μm的等轴晶，HAZ主要分布为尺寸约30 μm的粗大等轴晶，熔合线附近主要分布为垂直于温度场等温线方向生长的柱状晶，柱状晶宽度随着靠近热源也在不断生长。随着焊接线能量的提升，焊接接头的深宽比也呈显著增大的趋势。 （3）由于CuSi3焊丝导热系数远小于镀锌板母材，使得表面焊缝中心温度场分布梯度较小，焊接线能量对焊缝中心峰值温度影响不大。焊缝边缘贴近母材区域温度在四组参数下均高于锌镀层沸点，从而导致焊接试件表面产生气孔。 （4）后续研究将致力于建立焊缝区域钎料熔化过程的热-流耦合数值模型，并通过仿真和实验对钎料湿润、铺展的全过程进行研究，以帮助阐述焊接工艺参数对钎缝深宽比的影响机理。 参考文献 张屹 ， 李时春 ， 金湘中 ， 等 . 镀锌钢板激光焊接关键技术研究 ［J］. 激光与光电子学进展 ， 2010 （ 07 ）： 33 - 41 . 镀锌钢板激光焊接关键技术研究 ZHANG Y ， LI S C ， JIN X Z ， et al . Key Technology of Laser Welding of Galvanized Steel ［J］. Laser & Optoelectronics Progress ， 2010 （ 07 ）： 33 - 41 . Key Technology of Laser Welding of Galvanized Steel 李斌 ， 国莉 ， 郭涟 . 激光钎焊在白车身中的应用 ［C］. 中国汽车工程学会. 2013中国汽车工程学会年会论文集 . 北京 ： 北京理工大学出版社 ， 2013 ： 3 . 激光钎焊在白车身中的应用 马凯 . 激光钎焊条件下CuSi3钎料与镀锌钢板界面传质行为的研究 ［D］. 上海 ： 上海工程技术大学 ， 2015 . 激光钎焊条件下CuSi3钎料与镀锌钢板界面传质行为的研究 MA K . Research of interface mass transfer behavior under the condition of laser brazing of Zinc-coated steel with CuSi3 filler metal ［D］. Shanghai ： Shanghai University of Engineering Science ， 2015 . Research of interface mass transfer behavior under the condition of laser brazing of Zinc-coated steel with CuSi3 filler metal Min-Kyu Jeon ， Weon-Bae Kim ， Guk-Chan Han ， et al . A Study on Heat Flow and Temperature Monitoring in the Laser Brazing of a Pin-to-plate Joint ［J］. Journal of Material Processing Technology ， 1998 ， 82 （ 1-3 ）： 53 - 60 . A Study on Heat Flow and Temperature Monitoring in the Laser Brazing of a Pin-to-plate Joint J. S. Park ， S. J. Na . Heat Transfer in a Stud-to-plate Laser Braze Considering Filler Metal Movement ［J］. Wel-ding Journal . 1998 ， 77 （ 4 ）： 155 - 163 . Heat Transfer in a Stud-to-plate Laser Braze Considering Filler Metal Movement 赵振龙 . 轿车顶盖与侧围激光填丝焊接技术研究 ［D］. 吉林 ： 吉林大学 ， 2012 . 轿车顶盖与侧围激光填丝焊接技术研究 ZHAO Z L . Research on the laser filler wire welding technology of the car roof and side panel ［D］. Jilin ： Jilin University ， 2012 . Research on the laser filler wire welding technology of the car roof and side panel 宋永刚 ， 梁国星 ， 吕明 . 激光钎焊的温度场有限元分析 ［J］. 机械工程与自动化 ， 2009 （ 03 ）： 76 - 78 . 激光钎焊的温度场有限元分析 SONG Y G ， LIANG G X ， LV M . Finite Element Analysis of Temperature Field in Laser Brazing ［J］. Mechanical Engineering & Automation ， 2009 （ 03 ）： 76 - 78 . Finite Element Analysis of Temperature Field in Laser Brazing 封小松 ， 陈彦宾 ， 李俐群 . 镀锌板激光钎焊温度场的数值模拟 ［J］. 金属学报 ， 2006 （ 08 ）： 882 - 886 . 镀锌板激光钎焊温度场的数值模拟 FENG X S ， CHEN Y B ， LI L Q . Temperature Filed Simulation of Laser Brazing for Galvanized Steel Sh-eets ［J］. Acta Metallurgica Sinica ， 2006 （ 08 ）： 882 - 886 . Temperature Filed Simulation of Laser Brazing for Galvanized Steel Sh-eets 冯杰 . 镀锌板激光钎焊温度场及应力应变场分析 ［D］. 黑龙江 ： 哈尔滨工业大学 ， 2005 . 镀锌板激光钎焊温度场及应力应变场分析 FENG J . The Analysis of Temperature field and Stress and Strain field in Laser Brazing of Galvanized Steel Sheet ［D］. Heilongjiang ： Harbin Institute of Technology ， 2005 . The Analysis of Temperature field and Stress and Strain field in Laser Brazing of Galvanized Steel Sheet 金湘中 ， 李力钧 . 激光深熔焊接过程中小孔效应的试验研究 ［J］. 应用激光 ， 1999 （ 05 ）： 293 - 295，322 . 激光深熔焊接过程中小孔效应的试验研究 JIN X Z ， LI L J . An Experimental Study on Keyhole Effect in Deep Penetration Laser Welding ［J］. Applied Laser ， 1999 （ 05 ）： 293 - 295，322 . An Experimental Study on Keyhole Effect in Deep Penetration Laser Welding 席明哲 ， 虞钢 . 连续移动三维瞬态激光熔池温度场数值模拟 ［J］. 中国激光 ， 2004 ， 31 （ 12 ）： 1527 - 1532 . 连续移动三维瞬态激光熔池温度场数值模拟 XI M Z ， YU G . Numerical Simulation for the Transient Temperature Field of 3D Moving Laser Molten Pool ［J］. Chinese Journal of Lasers ， 2004 ， 31 （ 12 ）： 1527 - 1532 . Numerical Simulation for the Transient Temperature Field of 3D Moving Laser Molten Pool Reisgen U ， Schleser M ， Mokrov M ， et al . Optimization of laser welding of DP/TRIP steel sheets using statistical approach ［J］. Optics and Laser Technology ， 2012 ， 44 （ 1 ）： 255 - 262 . Optimization of laser welding of DP/TRIP steel sheets using statistical approach LI L Q ， FENG X S ， CHEN Y B . Influence of laser energy input mode on joint interface characteristics in laser brazing with Cu-base filler metal ［J］. Transactions of Nonferrous Metals Society of China ， 2008 （ 05 ）： 1065 - 1070 . Influence of laser energy input mode on joint interface characteristics in laser brazing with Cu-base filler metal 10.7512/j.issn.1001-2303.2022.11.03 TG402 A 1001-2303（2022）11-0014-08 刘海江，汪乾，张恒，等.基于CuSi3焊丝的镀锌钢板激光钎焊温度场及接头组织分析［J］.电焊机，2022，52（11）：14-21. 刘海江，汪乾，张恒，等.基于CuSi3焊丝的镀锌钢板激光钎焊温度场及接头组织分析［J］.电焊机，2022，52（11）：14-21. LIU Haijiang, WANG Qian, ZHANG Heng, et al.Temperature Field Analysis and Joint Microstructure of Laser Brazing of Galvanized Steel Sheet Based on CuSi3 Welding Wire[J].Electric Welding Machine, 2022, 52(11): 14-21. LIU Haijiang, WANG Qian, ZHANG Heng, et al.Temperature Field Analysis and Joint Microstructure of Laser Brazing of Galvanized Steel Sheet Based on CuSi3 Welding Wire[J].Electric Welding Machine, 2022, 52(11): 14-21. 2022-03-24 2022-09-24 2022-11-25 2023-01-17 电焊机 11 52