Automobile Technology & Material
6061铝合金搅拌摩擦焊接头失效行为仿真方法研究………………………………………………………………朱学武张鹏杨航杨化伟Research on Failure Behavior Simulation Method of 6061 Aluminum Alloy Friction Stir Welding Joint …………………………………………………………………………………………… Zhu Xuewu, Zhang Peng, Yang Hang, Yang Huawei
（1.中国第一汽车股份有限公司研发总院，长春 130013；2.汽车振动噪声与安全控制综合技术国家重点实验室，长春 130013； 3.汽车车身先进设计制造国家重点实验室，长沙 410082）
摘要：针对6061铝合金搅拌摩擦焊接头失效行为仿真方法进行研究。考虑焊缝及热影响区性能降低的规律，提出对接头不同区域进行精确划分的建模方法。对不同区域赋予不同材料力学性能，建立了铝合金搅拌摩擦焊接头拉伸样件的精细模型，并通过试验验证了精细模型的有效性。为保证在整车碰撞模型中的计算效率，提出将精细模型转化为等效简化模型的方法。将精细模型的材料属性、区域划分简化后赋予等效模型，其仿真峰值力与试验结果一致性良好，表明该等效模型可以表征搅拌摩擦焊接头力学性能，满足精度要求。
关键词：搅拌摩擦焊 精细模型 等效模型 仿真精度
中图分类号：U465 文献标识码：B 10.19710/J.cnki.1003-8817.20200320
DOI:
1 前言
随着国内新能源汽车产量的逐年增加[1]，铝合金材料在汽车中应用越来越多。搅拌摩擦焊是一种新型固相连接技术，同传统的熔化焊相比，搅拌摩擦焊有焊缝晶粒小，接头力学性能良好，焊接时无烟尘、无飞溅、焊后残余应力小等优点，在铝合金焊接中得到了广泛应用[2]
。通过工艺仿真，可以得到搅拌摩擦焊焊接过程中的温度场、等效塑性变形场分布以及材料流
作者简介：朱学武（1979—），男，高级工程师，硕士学位，研究方向为车体耐撞性能开发。基金项目：汽车振动噪声与安全控制综合技术国家重点实验室开放基金(1803J)。
参考文献引用格式：
朱学武, 张鹏, 杨航, 等. 6061铝合金搅拌摩擦焊接头失效行为仿真方法研究[J].汽车工艺与材料, 2021（8）：66-70.
ZHU X, ZHANG P, YANG H, et al. Research on Failure Behavior Simulation Method of 6061 Aluminum Alloy Friction Stir Welding Joint [J]. Auto⁃ mobile Technology & Material, 2021(8): 66-70.
动情况，帮助工程师提升对搅拌摩擦焊机理的理
解，同时也能对焊接参数的选择提供依据[3]。然而，对于搅拌摩擦焊在整车碰撞模型中的接头力学性能表征仍有需要研究的内容。同时，整车模型中模拟接头等效单元的大小和范围决定了在碰撞分析中的精度和计算成本。为保证计算效率一般不采用网格尺寸过小的精细模型，因此以精细模型为参照建立等效模型十分必要。
本研究建立搅拌摩擦焊接头拉伸样件的精细模型，该模型对焊接接头搅拌区域进行精确划分，并通过试验验证精细模型的有效性。为保证整车碰撞计算效率，提出一种将精细模型转化为等效简化模型的方法，将精细模型的材料属性、焊缝区和热影响区的区域划分简化后赋予等效模型，保证等效模型能够准确地模拟连接处的强度特性以及接头失效模式。简化模型的仿真峰值力与试验结果一致性较好，表明等效模型在整车碰撞模型中可以准确表征搅拌摩擦焊接头力学性能。
2 精细模型搭建2.1 焊接接头试样试验
3mm 6061-T6
试验采用 厚的 铝合金薄板对接
Nandan R，debroy T
等的研究结果[4]，对
接头。根据
3mm
于 厚铝合金板材，搅拌摩擦焊的合理参数范
800~1 500 r/min
围为搅拌头转速在 ，焊接速度在
300~1 200 mm/min
。本研究采用工艺参数为搅拌
1 200 r/min，焊接速度为800 mm/min。
头转速为获取焊接接头力学性能，采用准静态拉伸
GB/
试验对接头力学性能进行评估，将焊接样板按
T 228.1—2010《金属材料 1
拉伸试验 第 部分：室温试验方法》切割成标准准静态拉伸试样，尺寸如
1 26℃ 5 mm/min
图 所示。在 室温条件下，以 的速
Instron3309
度在 力学试验机上对接头样件进行拉伸试验，获取载荷-位移曲线。25
R
2.2 焊接接头力学性能搅拌摩擦焊是一种热输入很低的焊接方式[5]，最高温度始终低于母材熔点。6
系铝合金属于可热处理材料，焊接过程产生的热量仍会影响到材料性能。不同于熔焊，搅拌摩擦焊过程中还存在材料剧烈变形，故焊接完成后材料性能及组织会发生相应变化。搅拌摩擦焊接头由焊核区、热机影响区、热影响区以及母材组成。可以看出，不同区域的力学
6061-T6
性能均有所下降，这是由于焊接发热导致铝合金中的增强相熔化或粗化所导致。在焊核区由于受到搅拌头搅拌和摩擦作用，温度最高、变形最为剧烈，产生了动态再结晶，晶粒得到细化，从而缓解了热冲击导致的性能下降。热机影响区位于焊核区边缘，相比焊核区，受到的搅拌作用较弱，晶粒细化效果有限，而受到的热冲也比较大，所以性能最低。热影响区仅受到温度的影响，距离焊缝越远，材料性能损失越少并逐渐接近母材性能。由于热机影响区范围较小且对接头力学性能起决定性作用的是材料性能损失最严重的焊核区，所以在仿真模型中将热机影响区与热影响区划分到一起。2.3 精细模型搭建
为了验证搅拌摩擦焊接头模型的力学性能，需要先建立精细模型，在精细模型中对于焊接接头的不同区域进行精确划分和材料分配，有利于获得更为准确的仿真结果。
Hypermesh
根据标准拉伸样件尺寸使用 建立精
0.6 mm
细模型，网格采用 实体单元。由于搅拌摩擦焊成形质量很好，焊缝无余高焊根，所以焊接接头拉伸试样与母材准静态拉伸试样在形状上无差异，只在
2
试样中间划分了焊核区和热影响区域，如图 所示。
铝合金的焊接件力学性能可以由硬度来表征，材料的硬度和抗拉强度正相关，所以焊缝的区域可以根据硬度曲线来划分。由搅拌摩擦焊焊接试样拉伸试验可知，试样均在焊缝处失效，前进侧和后退侧的温度差异对于焊接件的失效位置影响较小，在整车模型中由于网格尺寸要远大于精细模型中的尺寸，前进侧和后退侧的差异可忽略。
3
根据图 硬度曲线可知，焊核区范围约为
10 mm ，材料性能损失区域总宽度约为17 mm。4
从图 焊接接头搅拌区形状可见，搅拌区是一5.32 mm，与硬度曲线对
个碗状区域，最小宽度为应的范围有差异。
5 470 ℃从图 工艺仿真的温度云图可见，温度左右的范围大致与硬度曲线中性能降低区域一致。这是因为搅拌区是搅拌头挤压旋转作用的区域，而搅拌头近似于圆台状，所以接头搅拌区的形状近似于碗状，但是焊核区还包括搅拌区附近的高温区域，搅拌区只是焊核区的一部分。结合硬度曲线和工艺仿真温度云图，对精细模型中热影响区和焊缝区进行划分，根据工艺仿真温度云图将焊缝划分成三部分，分别为焊缝Ⅰ，焊缝Ⅱ，焊缝Ⅲ，如图6
所示。拉伸试验工况较简单，仿真的精确度只与试样尺寸、网格大小、材料参数以及区域划分有关。试样尺寸以及网格大小在建模时可以确定，所以影响结果准确度的就是材料参数以及区域划分。区域划分已由硬度曲线以及工艺仿真温度云图确定，由于焊接时的热输入对于材料的弹性模量几乎没有影响，所以对于不同区域的材料参数只改变塑性段。对于性能降低区域材料参数的确定，通过测量焊接温度场划分各区域的范围，再通过拉伸试验获取不同区域的材料参数。不同区域材料失效参数是在母材失效参数基础上由拉伸试验标定获得。3 精细模型试验验证3.1 精细模型变形失效仿真及试验对标
搅拌摩擦焊是通过搅拌头将两个工件的材料搅动到一起实现连接，同时搅拌头的轴肩会压入材料，所以焊缝的厚度要低于其他区域，同时焊缝区域温度最高，材料性能损失最严重，故搅拌摩擦
焊焊接试样会在焊缝处断裂失效。从图7可看出，试样是沿着焊缝Ⅱ区域断裂，且位于焊缝边缘的热机和热影响区域，仿真与试验的失效模式一致。8
焊接试样等效塑性应变云图如图 所示。当拉t=0.3 ms
伸到 时，试样开始出现塑性变形，最大等效0.004594，位于焊缝Ⅱ处；当t=2.5 ms
塑性应变为 时，试样开始出现明显颈缩，最大等效塑性应变为0.3128；当t=2.9 ms
时，焊缝单元开始失效删除，最大1.086；当t=3.3 ms
等效塑性应变为 时，试样完全断裂。3.2 精细模型拉伸变形力仿真及试验对标
接头在拉伸工况下仿真与试验的载荷-位移曲
9线如图 所示。拉伸工况下仿真模型最大峰值力为
8.22 kn，试验测得最大峰值力为8.25 kn，通过软件
91.4%计算得两条曲线的拟合度为 ，满足精度要求。试验曲线没有明显的弹性段，这是因为接头材料受到热影响后性能会降低，且是连续变化的，但是在仿真中，不同区域的材料性能是不连续的，所
4 等效模型搭建
4.1 等效模型搭建rigid传统焊接接头模拟方法包括 刚体单元、共节点、单根梁单元、多根梁单元、实体焊点单元等。搅拌摩擦焊接头的连接方式是对接，连接区域是一条焊缝而不是多个连接点，所以对于搅拌摩擦焊接头的等效模型，本研究确定的表征方式是采用壳单元并赋予各区域不同材料属性。
在等效模型中，受制于单元类型和网格尺寸，无法实现精细模型中的精确分区，所以在等效模型中将性能相近的区域合并。为避免模型由于简化造成强度高于精细模型，将合并的区域赋予精细模型中相应区域材料性能最弱的参数。根据精细模型仿真结果可知，失效位置位于焊缝处，所以在等效模型中将焊缝Ⅱ、焊缝Ⅲ合并成焊缝区，赋予焊缝Ⅱ的材料属性；将焊缝Ⅰ与热影响区合并成热影响区，赋予焊缝Ⅰ材料。等效模型中的焊缝、热影响区尺寸范围
10 3mm。
如图 所示。等效模型中网格尺寸为4.2 网格尺寸对仿真精度影响
3mm
本研究等效模型的网格尺寸为 ，在此基
2mm 4mm 5mm 6mm作
础上，将网格尺寸设为 、 、 、
11
为对比。不同网格大小的力与位移曲线如图 所
2～ 4mm
示。当网格尺寸为 时，曲线无明显差异，
5mm
当网格尺寸增加到 时，最大峰值力没有明显差异，但是断裂位移有明显增加，故焊缝处网格尺
4mm。
寸不能大于
5 等效模型试验验证5.1 等效模型变形失效仿真及试验对标
等效模型不同于精细模型的详细划分，焊缝区域简化成长方形，所以断裂会在焊缝的中心，如
12
图 所示。段裂位置大致与试验结果一致，均位于焊缝处，只是分别在焊缝的中心和一侧。
13
等效模型的等效塑性应变场如图 所示。当
t=0.3 ms时，试样开始产生塑性应变，位于焊缝，最
大等效塑性应变为0.004 019，时间、应变值均与精
细模型一致；当t=2.3 ms时，试样开始出现明显颈缩，热影响区也开始出现塑性变形，最大等效塑性应变为0.166 4；当t=3 ms时，模型单元开始失效删
除，此时最大等效塑性应变为0.236；当t=3.1 ms，试样完全断裂失效。5.2 等效模型拉伸变形力仿真及试验对标
等效模型仿真结果与精细模型仿真结果、试
14
验结果的力与位移对比曲线如图 所示。峰值力
8.25 kn、8.22 kn、8.27 kn，通过软件算得等
分别为
86.1%，满足精
效模型曲线与试验曲线的拟合度为度要求。等效模型力与位移曲线中部与精细模型和试验都有较大差异，原因是等效模型分区更少，材料性能分布差异更大。
6 结论
根据拉伸试验试样尺寸及工况，建立了搅拌摩擦焊接头失效的精细模型。综合硬度曲线和工艺仿真温度云图，对精细模型中热影响区和焊缝区进行划分，对不同区域赋予合理的材料参数。精细模型的仿真峰值力、失效模式与试验结果的一致良
-位移曲线的拟好，通过软件计算得仿真和试验力
91.4%。为满足整车碰撞工况中计算效率需合度为求，根据精细模型的材料属性和焊缝区、热影响区的划分建立了等效模型，等效模型仿真峰值力与试验结果接近，通过软件计算得仿真和试验力-位移
86.1%
曲线的拟合度为 。研究结果表明，本次建立的等效模型可以表征搅拌摩擦焊接头力学性能，满
CAE
足整车 仿真精度要求。
参考文献：
[1] 马少超,范英.
基于时间序列调整的中国新能源汽车政
策评估[J]. 中国人口·资源与环境,2018(4):117-124. [2] 韩振宇,朱俊,赵北明,等. 6063
铝合金电池盒下壳体搅
拌摩擦焊工艺试验研究[J]. 铝合金加工技术,2019(47): 56-61．
[3] 鄢东洋.铝合金薄壁结构搅拌摩擦焊热-力学过程的研
究及模拟[D]. 北京:清华大学,2010.
[4] Nandan R, Debroy T, Bhadeshia H K D H. Recent advances in friction- stir welding Process, weldment
–
structure and properties[j]. Progress in Materials Science, 2008, 53(6):980-1023.
[5] Mishra R S, Ma Z Y. Friction Stir Welding and Processing II[J]. Materials Science and Engineering Reports, 2005, 50 (1):1-78.