金属学报, 2026, 62(1): 148-158 DOI: 10.11900/0412.1961.2025.00240

研究论文

铝合金搅拌摩擦焊热--流耦合仿真及材料迁移轨迹解析

韩鸿华, 史清宇, 杨诚乐, 孔德帅, 陈高强,

清华大学 机械工程系 清洁高效透平动力装备全国重点实验室 北京 100084

Coupled Thermal-Mechanical-Fluid Simulation and Material Migration Trajectory Analysis for Friction Stir Welding of Aluminum Alloy

HAN Honghua, SHI Qingyu, YANG Chengle, KONG Deshuai, CHEN Gaoqiang,

State Key Laboratory of Clean and Efficient Turbomachinery Power Equipment, Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China

通讯作者: 陈高强,cheng1@tsinghua.edu.cn,主要从事搅拌摩擦焊物理过程的建模仿真等研究

收稿日期: 2025-08-20   修回日期: 2025-09-22  

基金资助: 国家自然科学基金项目(52175334)
清洁高效透平动力装备全国重点实验室课题项目(DEC8300CG202428557A-1228215)

Corresponding authors: CHEN Gaoqiang, associate professor, Tel:(010)62773637, E-mail:cheng1@tsinghua.edu.cn

Received: 2025-08-20   Revised: 2025-09-22  

Fund supported: National Natural Science Foundation of China(52175334)
State Key Laboratory of Clean and Efficient Turbomachinery Power Equipment Project(DEC8300CG202428557A-1228215)

作者简介 About authors

韩鸿华,男,1978年生,研究员,博士生

摘要

搅拌摩擦焊(FSW)过程中,材料在数秒时间内迅速完成材料流动-沉积过程。为解析材料流经搅拌头附近的完整轨迹,本工作以2219-T87铝合金FSW过程为研究对象,构建了FSW过程热-力-流全耦合仿真模型,并运用轨迹追踪算法,着重分析了FSW过程中材料流动-沉积的轨迹。结果表明,焊接过程中材料流动-沉积行为受其所处位置影响:靠近轴肩的材料流经搅拌头附近的轨迹呈多圈绕流模式,材料在流动过程中,绕行半径由大变小,并在搅拌针附近发生显著的向下迁移,沉积位置均位于前进侧;远离轴肩的材料呈现出直通流动模式,材料未完成整圈流动后便沉积在搅拌针后方,前进侧材料焊后沉积在前进侧,后退侧材料焊后沉积在后退侧。

关键词: 搅拌摩擦焊; 铝合金; 材料流动; 数值模拟

Abstract

Although friction stir welding (FSW) is widely employed in aerospace and other engineering fields owing to its ability to produce high-quality joints and low residual stresses, the understanding of its physical mechanisms remains limited by the intense material flow during welding, which lags behind engineering practice. This gap motivates researchers to develop thermo-mechanical coupling models to elucidate the underlying physics, aiming to expand its application boundaries. During FSW, the material undergoes a highly rapid flow-deposition process that is completed within a few seconds. In this study, we present a fully coupled thermo-mechanical-fluid model for friction-stir-welded 2219-T87 aluminum alloy to demonstrate the complete material transport path in the vicinity of the tool. A particle-tracking algorithm was employed to quantitatively analyze the flow and subsequent refilling trajectory of the tracer material. The predicted temperature history and flow-zone geometry showed close agreement with experimental data. At the axis Z = 7 mm, the peak downward vertical velocity of the material reached 16.57 mm/s at the advancing side of the pin, whereas the peak upward vertical velocity reached 10.47 mm/s at the retreating side. The analysis revealed that the flow-deposition behavior of the material was highly position dependent. The material adjacent to the shoulder followed a multi-loop spiral trajectory around the tool, during which the loop radius gradually decreased and the material migrated markedly downward in the vicinity of the pin. Ultimately, all deposits accumulated on the advancing side. In contrast, the material far from the shoulder exhibited a through-flow pattern, depositing behind the pin without completing a full revolution. Post-weld deposits originating from the advancing side remained on the advancing side, whereas those from the retreating side settled on the retreating side.

Keywords: friction stir welding; aluminum alloy; material flow; numerical simulation

PDF (2622KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

韩鸿华, 史清宇, 杨诚乐, 孔德帅, 陈高强. 铝合金搅拌摩擦焊热--流耦合仿真及材料迁移轨迹解析[J]. 金属学报, 2026, 62(1): 148-158 DOI:10.11900/0412.1961.2025.00240

HAN Honghua, SHI Qingyu, YANG Chengle, KONG Deshuai, CHEN Gaoqiang. Coupled Thermal-Mechanical-Fluid Simulation and Material Migration Trajectory Analysis for Friction Stir Welding of Aluminum Alloy[J]. Acta Metallurgica Sinica, 2026, 62(1): 148-158 DOI:10.11900/0412.1961.2025.00240

搅拌摩擦焊(friction stir welding,FSW)问世30余年以来,已广泛应用于运载火箭贮箱、高铁列车车体、新能源汽车电池等关键结构的制造中。然而,对FSW物理过程的认识仍大幅落后于其工程应用进程。FSW物理过程具有鲜明特征:材料虽不发生熔化,但在流经搅拌头附近的数秒时间内发生剧烈的塑性流动并沉积焊缝中。目前,针对FSW材料流动的研究大多集中在分析材料流动的速度及其分布[1~8],但由于在实验观测上缺少有效的研究手段,对材料流经搅拌头附近的完整轨迹仍缺少深入认识。

为探究FSW过程中材料流动-沉积行为,研究人员通过实验和数值模拟进行了广泛研究,包括示踪剂追踪技术[9~15]和数值计算方法[16~23],许多基于实验和模拟的研究揭示了材料在搅拌头周围表现出多圈绕流模式和直通流动模式[23~28]。实验可视化研究主要通过引入与母材物性存在差异的材料作为标记物来实现,利用实验表征手段直接呈现材料流动过程[23,29]。研究人员采用铜条[9]、钨球[10]、铜箔[12]、增厚氧化层[13]和异种焊丝[14]等材料作为标记物,追踪材料的流动轨迹和沉积位置。但是采用追踪材料的实验方法也存在一定的局限性。Schimdt等[9]通过焊缝中心纵向、焊缝横向位置设置0.1 mm厚铜条作为追踪材料,对FSW过程的试板进行计算机断层扫描(CT)分析,观察到追踪材料在轴肩部位呈“S”形或“C”形轨迹,并从材料流动速度和沉积位置,判断出材料在沉积前完成多次旋转,由于实验的局限性,无法定量给出某点追踪材料的流动轨迹和沉积位置。Morisada等[10]分别在A1050铝合金和低碳钢搅拌区焊缝中心预制孔中嵌入钨球,搅拌区±30°的位置放置了两组X射线透射系统,通过观察,钨示踪剂材料在铝合金焊缝搅拌头周围存在多圈绕流现象并有沿垂向方向的迁移,而钨示踪剂材料在低碳钢焊缝中直接绕过后退侧沉积在搅拌头后方,呈现出直通绕流模式,实验中,由于钨球嵌入的位置有限,无法全面反映整个焊接区域的材料流动行为。Huang等[12]和Zeng等[13]采用人工增厚氧化层和分层铜箔直接追踪材料水平流动和垂向流动过程的迁移行为,但是未有效获得三维流动模型。因而,目前FSW材料流动-沉积行为的实验表征研究主要集中在搅拌过程中的材料迁移机理,即焊缝宏观结构的形成,关于材料流动过程中三维空间下的流速、应变及应变速率等关键物理参量的实验表征研究仍严重不足[23,29],且缺乏对材料流经搅拌头复杂几何区域周围完整的轨迹认识。

数值模拟方法是一种获取FSW过程中材料流动和迁移行为的有效工具,主要分为计算固体力学(CSM)和计算流体力学(CFD)两种数值计算方法。基于CSM方法,将发生流动的材料视为塑性变形材料,通过任意Lagrange-Euler (ALE)方程或耦合Euler-Lagrange (CEL)方程预测材料的温度场、应力场、速度场和宏观界面形状等[25]。基于CFD的热-力-流耦合模型能够全面揭示焊接过程中发生的复杂物理现象,并与实验结果相结合进行相互验证,从而为材料流动行为提供更深入的理解。在CFD数值计算方法中,主流的边界模型为速度边界模型和摩擦力边界模型。速度边界模型把滑移系数设为常值[17,19]或忽略了实际焊接过程中可能出现的滑移现象[16,20],导致搅拌头肩部外边缘材料流动速率过大或直接忽视与实验数据的对比。摩擦力边界模型减少了人为假设的因素,可实现工具-工件(T/W)界面接触状态的自适应转换[21,22,30,31],能够更加真实地反映材料流动行为,但需要建立准确的剪切应力条件[32]。Liechty和Webb[18]对比了摩擦力边界模型与速度边界模型的模拟结果,认为采用摩擦力边界模型能更好地反映材料流动行为。Chen等[21]和Qiao等[33]建立了基于摩擦力边界的耦合模型,并预测得到水平面上的材料流动模式。Yang等[22]在此基础上准确预测了3 mm薄板材料的垂直流动行为,并定量分析了材料垂向流动的方向和移动距离。

综上,采用实验的方法研究材料流动-沉积行为存在较大局限性,而数值模拟是研究材料流动-沉积行为的有效方法,有必要运用数值模拟方法对FSW过程中材料流经搅拌头附近的完整轨迹进行深入研究。本工作通过FSW焊接实验,测量了焊接过程中的温度,获得了焊接接头的几何形貌,并基于计算流体动力学建立热-力-流全耦合仿真模型,利用实验结果对模拟预测的温度和流动区域的几何形状进行验证,揭示了FSW过程中材料流动的复杂机制,并定量分析了材料的流动模式及其迁移行为。

1 实验方法

FSW实验设备为XA5032型机床,工件材料为2219-T87铝合金,尺寸为400 mm × 100 mm × 8 mm。工件尺寸、焊接方向及搅拌头旋转方向如图1所示,原点O位于搅拌头中心轴与工件背面的交点,其中X轴表示焊接方向,Y轴表示与焊缝垂直的横向方向,Z轴表示工件厚度方向。焊接工艺参数为:焊接速率180 mm/min,搅拌头转速800 r/min,计算模型中的焊接速率和焊接转速与实际焊接参数保持一致。焊接过程中,底部垫板为钢材料,搅拌头材料为H13钢。搅拌头的轴肩直径为20 mm,螺纹搅拌针呈圆台形,搅拌针长度为7.8 mm,大端直径为7.8 mm,小端直径为4.2 mm,搅拌头倾斜角为2.7°。在焊接过程中将K型热电偶(温度测量范围0~1300 ℃)放置于焊缝的前进侧(AS)和后退侧(RS),距离焊缝中心14 mm、深度1 mm的位置处。使用LR8450-01型温度记录仪记录焊接过程中的温度变化曲线,采样频率50 Hz。焊接实验完成后,从焊缝接头处取样。将切取的样品经磨平、抛光处理后,使用Keller试剂(95 mL H2O + 1 mL HF + 1.5 mL HCl + 2.5 mL HNO3)浸泡15 s进行腐蚀处理,利用BX51M型光学显微镜(OM)观察焊接接头的宏观几何形貌。

图1

图1   搅拌摩擦焊(FSW)实验过程示意图

Fig.1   Schematic of the friction stir welding (FSW) experimental process


2 计算模型

2.1 几何模型

本工作采用基于计算流体力学的热-力-流全耦合数值模型,用于模拟2219-T87的焊接过程,几何模型如图2所示。网格划分采用六面体网格,总网格数量为619208个,节点数为651251个。在搅拌头与工件界面及其附近区域,采用更精细的网格划分,其最小尺寸约为0.1 mm,该网格质量确保了计算的收敛性,并能够准确表征材料的剧烈流动行为。在模型中,忽略搅拌头倾斜角度及下压量引起的几何变化,搅拌头在模型中的位置保持固定。热边界条件设置如下:“材料入口”和“材料出口”处的温度均设定为27 ℃,工件与侧壁及垫板的接触传热系数为500 W/(m2·K),工件上表面与空气之间的对流传热系数为50 W/(m2·K)[21,34,35]。工件材料的密度(ρ)为2840 kg/m³;比热容(cp)、热导率(λ)的取值参考文献[36,37]。数值模型采用商业软件ANSYS Fluent进行求解。

图2

图2   几何模型及六面体网格划分

Fig.2   Geometric modeling and hexahedral meshing


2.2 控制方程

本工作中,工件材料被假设为不可压缩非Newton流体,该模型通过求解连续性方程、动量方程和能量方程,描述焊接过程中材料的流动行为以及热量的产生和传递过程。

连续性方程:

ρt+ρv=0

动量守恒方程:

ρvt+ρvv=-p+μv+vT

能量守恒方程:

ρHt+ρvH=λT+Sv

式中, v 为流体速度;t为流动时间;p为压力;μ为材料黏度;T为温度;Sv为工件材料塑性变形产热的源项;H为热焓,可表示为:

H=TrefTcpdT

式中,Tref为环境温度,本工作中为27 ℃。

2.3 摩擦边界条件

将搅拌头施加在工件上产生的摩擦剪应力(τf)表达式为[21, 34]

τf=μfσnvrelvreltanhαvrel

式中,μf为Coulomb摩擦系数,参考文献[21,22],在本工作计算中取0.33;σn为T/W界面上的法向压强;vrel为界面相对速度,可以表达为vrel=vt-vwpvtvwp分别为界面上搅拌头速度和工件速度;α为控制因子,取值50 s/m[21]。在FSW过程中,搅拌针螺纹对材料流动具有重要影响:一方面,螺纹表面与塑性变形区之间的摩擦力驱动螺纹槽内材料作圆周运动;另一方面,搅拌头旋转前进时垂直于螺纹表面的压力推动材料产生垂直方向运动[38,39]。在本模型中,搅拌针螺纹界面速度(vthread)可以表达成vt和螺纹旋转产生的附加垂向速度(vlp)的矢量和:

vthread=vt+vlp

vlp可表达成:

vlp=κω2π

式中,κ为螺纹结构经验系数,本工作中取值为0.002 m/rev;ω为搅拌头角速度。

2.4 产热模型及散热边界条件

FSW过程中,产生的热量主要来源于两种机制:第一种是搅拌头与工件界面之间的相对滑动摩擦所产生的热量,第二种是由于材料流动过程中发生塑性变形而产生的热量。在本工作的数值模拟中,界面摩擦产热热流密度(qf)主要由搅拌头T/W界面产热量进入工件的比例(η,本工作取0.63[21])、τf以及vrel决定:

qf=ητfvrel

塑性流动产热是能量方程中的源项Sv,塑性流动产热热流密度(qbd)由塑性功的热耗散系数(β)、流动应力(σ¯)以及等效应变速率(ε˙)决定:

qbd=βσ¯ε˙

工件上表面、下底面和侧面对流换热边界条件可以用通式表达:

λTnb=hT-T0

式中,nb为边界外法线方向的单位向量,Tnb为边界附近沿法线方向的温度梯度,h为换热系数,T0为初始温度。

2.5 材料本构模型

为了准确描述流体的流动特性,将流体的黏度表示为流动应力与应变(ε)、等效应变速率及温度的函数[16,40]

μ=σ¯ε, ε˙, T3ε˙

按照材料热变形Sellars-Tegart本构模型,σ¯可描述为:

σ¯ε, ε˙, T=σpsinh-1ZT, ε˙A1n
ZT, ε˙=ε˙expQRT
ε˙=23εijεij12
εij=12uixj+ujxi

式(12)~(15)中,σp为材料应力水平参数,取值50 MPa;ZT, ε˙为Zener-Hollomon参数;εij为应变速率张量的分量;R为理想气体常数,取值8.314 J/(mol·K);Q为热变形激活能,取值153300 J/mol;n为应力指数,取值6.91;A为结构因数,取值2.65 × 1010 s-1 [41]uiuj为三维空间中的速度矢量分量,xixj为三维空间中的位置。材料的固相线温度为543 ℃[42,43],在模拟过程中采用经验软化法[44~47],修正固相线软化区域的流变应力。

3 结果与讨论

3.1 焊接过程中的温度分布

图3a给出了FSW过程中工件表面的温度分布预测结果。可以看出,轴肩与搅拌针根部连接处出现局部高温区,最高温度为524.1 ℃,位于搅拌头侧后方。搅拌头前方的温度梯度明显高于后方,符合焊接温度场的典型特征。图3b给出了经过搅拌头旋转中心X = 0 mm截面上的温度分布情况。在搅拌针与工件界面附近的温度分布并不均匀,形成了碗状的局部高温区域,AS的最高温度为513.6 ℃,而RS的最高温度为483.3 ℃,焊接过程中AS的温度高于RS的现象,与其他学者[3,4]的实验和仿真结果一致。为验证温度场预测结果的可靠性,将仿真结果与测温结果进行对比。图3cd给出了图1所示热电偶安装处温度仿真结果与实测结果。距离焊缝中心线Y = -14 mm、Z = 7 mm处AS位置的预测最高温度为301.3 ℃,实测最高温度为304.0 ℃。距离焊缝中心线Y = 14 mm、Z = 7 mm处RS位置的预测最高温度为296.7 ℃,实测最高温度为288.3 ℃。测温点的最高温度预测结果与实验结果之间的误差均在3%以内。这表明,本工作构建的仿真模型在预测FSW的过程温度场方面具有很高的精度。FSW过程中搅拌头附近高温软化的材料发生塑性流动,准确的温度场预测模型为下文分析材料流动奠定了基础。

图3

图3   FSW过程中工件表面的温度分布预测结果,X = 0 mm截面上的温度分布,及热电偶安装处温度仿真结果与实测结果

Fig.3   Simulated temperature distributions of full workpiece surface (a) and X = 0 mm section (b), and simulated and measured temperatures of AS (c) and RS (d) during FSW (Tmax—maximum simulated temperature of the workpiece; TAS_max, TRS_max—maximum temperatures on the AS and RS of the X = 0 section, respectively; Tmax.Sim., Tmax.Exp.—maximum simulated and experimentally measured temperatures at the thermocouple location, respectively)


3.2 焊接过程中的材料流动行为

图4给出了X = 0截面材料流动速度的分布结果。从轴肩外缘到搅拌针根部,材料流动速率逐渐增大,在搅拌针根部Y = -3.9 mm和Y = 3.9 mm处,材料流动速率最大值分别为285和277 mm/s。从搅拌针根部到搅拌针端部,材料流动速率逐渐减小。图5给出了焊缝塑性变形区预测结果和实验结果。图5a为搅拌头附近X = 0 mm截面上应变速率的预测结果。结果表明,塑性变形区呈近似碗状分布,AS塑性变形区边界与RS塑性变形区边界相对于搅拌头中心轴表现出明显的不对称性,RS的塑性变形区宽度明显大于AS的塑性变形区宽度,搅拌针上部的塑性变形区边界尺寸明显大于搅拌针下部的塑性变形区边界尺寸。预测的轴肩塑性变形区的宽度为18.84 mm,小于搅拌头轴肩直径(20 mm),这表明轴肩边缘附近的材料受到未发生塑性流动的固态材料的约束,这种流动范围的限制有助于减少材料在轴肩边缘形成飞边。图5b给出了焊缝宏观形貌。以0.3 s-1应变速率作为热机械影响区(thermo-mechanical affected zone,TMAZ)的边界,预测的塑性变形区形状与实验测得的TMAZ边界具有较高的吻合度。

图4

图4   X = 0 mm截面材料流动速度的分布结果

Fig.4   Velocity magnitude of material flow at the X = 0 mm section


图5

图5   焊缝塑性变形搅拌头附近X = 0 mm截面上应变速率的预测结果及其实验宏观形貌

Fig.5   Predicted strain rate at X = 0 mm section (a) and experimental macromorphology (b) of weld plastic deformation zone (TMZ—thermo-mechanical affected zone)


为进一步探究FSW过程中的材料流动机制,测量并比较了材料在水平面内的流动速度和沿板厚方向的垂向流速,以揭示材料三维流动特征。图6给出了Z = 7 mm和Z = 3 mm两个典型平面位置的材料流动速度分布。结果表明,越靠近轴肩位置,材料流动速率越大。在Z = 7 mm平面处,材料流动主要受轴肩影响,峰值速率为241 mm/s,从搅拌针根部到端部,材料流动速率不断降低,在Z = 3 mm平面处,材料流动峰值速率已降低到74 mm/s。

图6

图6   Z = 7 mm和Z = 3 mm两个平面位置材料的流动速度分布

Fig.6   Velocity magnitudes at Z = 7 mm plane (a) and Z = 3 mm plane (b)


图7给出了Z = 7 mm和Z = 3 mm两个典型平面位置材料在垂直方向上的速度分布。结果表明,材料在搅拌针附近发生不同程度的垂向流动,材料的垂直流动速度呈现出明显的空间分布特征。越靠近轴肩位置,材料垂直方向的流动速率越大。从图7a可以看出,在Z = 7 mm平面位置,材料流动主要受轴肩影响,材料垂直向下的流动速率峰值为16.57 mm/s,位于搅拌针AS;材料垂直向上的流动速率峰值为10.47 mm/s,位于搅拌针RS。这表明在靠近轴肩的区域,材料的垂直向下流动现象更为显著,且主要集中在AS。从图7b可以看出,随着材料向搅拌针端部靠近,材料温度降低使流变应力增加,导致摩擦力不足以驱动材料发生塑性变形和流动,其垂直流动速率也随之减小。

图7

图7   Z = 7 mm和Z = 3 mm两个平面位置材料在垂直方向上的速度分布

Fig.7   Z-direction velocity distributions at the Z = 7 mm plane (a) and Z = 3 mm plane (b)


3.3 材料流动-沉积轨迹

图8给出了Z = 7 mm和Z = 3 mm两个典型平面材料流动轨迹的分布情况。图8a~c给出了Z = 7 mm平面位置的材料流动轨迹。该区域的材料流动主要由轴肩驱动,无论是AS还是RS的材料,均在搅拌头周围经历了多次环绕,形成了多圈绕流区,最后沉积在搅拌头后方。这种绕流区在Schmidt等[9]的实验研究和Morisada等[10]的原位观测中均得到证实。材料的绕流范围始终限制在轴肩区域内,这是由于材料在轴肩范围内的流动宽度小于轴肩尺寸(图5)。图8d~f给出了在Z = 3 mm平面位置的材料流动轨迹。该区域主要由搅拌针驱动,材料无法完全绕搅拌针旋转整圈,材料绕过搅拌针后直接沉积在搅拌头后方。值得注意的是,由于该位置更靠近搅拌针端部,RS材料在绕流不到半圈时便沉积在搅拌针后方。

图8

图8   Z = 7 mm和Z = 3 mm平面材料流动轨迹的分布

Fig.8   Distributions of the material flow trajectory at Z = 7 mm plane (a-c) and Z = 3 mm plane (d-f)

(a, d) overall views (b, e) AS (c, f) RS


3.4 多圈绕流模式流动轨迹的特征

在FSW过程中,材料流动行为对焊缝缺陷的形成、微观组织的演变以及焊接接头的力学性能起着决定性作用[24,25],因此有必要详细讨论材料的流动模式。图910给出了Y = ± 2.5 mm、Z = 7 mm处(即靠近轴肩下方)附近的典型流线的多圈绕流路径及Z向速度分布。从图9b10b可以看出,受轴肩影响,来自AS和RS的材料流入搅拌针高速流动区后,围绕搅拌针完成多个整圈绕流,材料在搅拌头前方绕行半径小,紧紧围绕搅拌针根部,材料在搅拌针后方,绕行半径逐步由大变小,最终均沉积在搅拌针后方AS。从图9c10c可以看出,材料存在明显的垂向迁移行为,在搅拌针附近先向上迁移,再向下迁移,最终分别沉积在初始流入位置下方1.94 mm处(图9c)和1.13 mm处(图10c)。总体而言,靠近轴肩的材料表现为明显的多圈绕流模式,这是由于该处的温度高,剪切阻力小,材料流动速率高,可充分驱动材料发生多圈绕流。

图9

图9   起始于Y = -2.5 mm、Z = 7 mm的多圈绕流路径及其Z向速度分布

Fig.9   Multi-revolution flow paths initiated at Y = -2.5 mm, Z = 7 mm and its Z-direction velocity distribution (ΔYZ) denotes the offset in the Y(Z) direction between the material outflow and initial inflow positions, the same in Fig.10)

(a) 3D view (b) X-Y view (c) X-Z view


图10

图10   起始于Y = 2.5 mm、Z = 7 mm的多圈绕流路径及其Z向速度分布

Fig.10   Multi-revolution flow paths initiated at Y = 2.5 mm, Z = 7 mm and its Z-direction velocity distribution

(a) 3D view (b) X-Y view (c) X-Z view


3.5 直通绕流模式流动轨迹的特征

图1112分别给出了Y = ±1.2 mm、Z = 3 mm处典型流线的直通流动路径及Z向速度分布。从图11b12b可以看出,在远离轴肩位置,来自AS的材料在绕流搅拌针半圈后沉积在搅拌针后方AS,而来自RS的材料在流经搅拌针后快速脱离搅拌针流动区,并沉积在搅拌针后方RS,原因是该位置附近材料温度低(图3b),材料剪切阻力大,材料流动速率显著降低(图6b),无法有效形成整圈绕流。从图11c12c可以看出,在搅拌针端部附近,材料的垂直迁移能力逐步降低,这是材料的垂向流动速率低导致的。

图11

图11   起始于Y = -1.2 mm、Z = 3 mm 的直通流动路径及其Z向速度分布

Fig.11   Straight-through flow paths initiated at Y = -1.2 mm, Z = 3 mm and its Z-direction velocity distribution

(a) 3D view (b) X-Y view (c) X-Z view


图12

图12   起始于Y = 1.2 mm、Z = 3 mm的直通流动路径及其Z向速度分布

Fig.12   Straight-through flow paths initiated at Y = 1.2 mm, Z = 3 mm and its Z-direction velocity distribution

(a) 3D view (b) X-Y view (c) X-Z view


4 结论

(1) 基于CFD建立了8 mm厚2219-T87铝合金板材FSW过程的热-力-流全耦合仿真模型,实现了FSW过程中的温度场和流动场的仿真预测。预测得到的温度曲线和塑性变形区的几何形状均与实验结果吻合良好。

(2) 靠近轴肩的材料流经搅拌头附近的轨迹呈多圈绕流模式。发生多圈绕流的材料在流动过程中,绕行半径由大变小,并在搅拌针附近发生显著的向下迁移后沉积。来自前进侧材料与来自后退侧的材料均可发生多圈绕流,但发生多圈绕流后沉积位置均位于前进侧。

(3) 远离轴肩的材料呈现出直通流动模式,材料未完成整圈后便沉积在搅拌针后方。直通流动模式下,来自前进侧材料焊后沉积在前进侧,而来自后退侧材料沉积在后退侧。

参考文献

Schmidt H, Hattel J, Wert J.

An analytical model for the heat generation in friction stir welding

[J]. Modell. Simul. Mater. Sci. Eng., 2004, 12: 143

[本文引用: 1]

Colligan K J, Mishra R S.

A conceptual model for the process variables related to heat generation in friction stir welding of aluminum

[J]. Scr. Mater., 2008, 58: 327

Chen G Q, Shi Q Y, Li Y J, et al.

Computational fluid dynamics studies on heat generation during friction stir welding of aluminum alloy

[J]. Comput. Mater. Sci., 2013, 79: 540

[本文引用: 1]

Su H, Wu C S, Pittner A, et al.

Thermal energy generation and distribution in friction stir welding of aluminum alloys

[J]. Energy, 2014, 77: 720

[本文引用: 1]

Fratini L, Buffa G, Palmeri D, et al.

Material flow in FSW of AA7075-T6 butt joints: Continuous dynamic recrystallization phenomena

[J]. J. Eng. Mater. Technol., 2006, 128: 428

Nie L, Wu Y X, Gong H.

Prediction of temperature and residual stress distributions in friction stir welding of aluminum alloy

[J]. Int. J. Adv. Manuf. Technol., 2020, 106: 3301

Heidarzadeh A, Mironov S, Kaibyshev R, et al.

Friction stir welding/processing of metals and alloys: A comprehensive review on microstructural evolution

[J]. Prog. Mater. Sci., 2021, 117: 100752

Wu C S, Su H, Shi L.

Numerical simulation of heat generation, heat transfer and material flow in friction stir welding

[J]. Acta Metall. Sin., 2018, 54: 265

[本文引用: 1]

武传松, 宿 浩, 石 磊.

搅拌摩擦焊接产热传热过程与材料流动的数值模拟

[J]. 金属学报, 2018, 54: 265

[本文引用: 1]

Schmidt H N B, Dickerson T L, Hattel J H.

Material flow in butt friction stir welds in AA2024-T3

[J]. Acta Mater., 2006, 54: 1199

[本文引用: 4]

Morisada Y, Imaizumi T, Fujii H, et al.

Three-dimensional visualization of material flow during friction stir welding of steel and aluminum

[J]. J. Mater. Eng. Perform., 2014, 23: 4143

[本文引用: 3]

Morisada Y, Imaizumi T, Fujii H.

Determination of strain rate in Friction Stir Welding by three-dimensional visualization of material flow using X-ray radiography

[J]. Scr. Mater., 2015, 106: 57

Huang Y X, Wang Y B, Wan D, et al.

Material-flow behavior during friction-stir welding of 6082-T6 aluminum alloy

[J]. Int. J. Adv. Manuf. Technol., 2016, 87: 1115

[本文引用: 2]

Zeng X H, Xue P, Wang L, et al.

Material flow and void defect formation in friction stir welding of aluminium alloys

[J]. Sci. Technol. Weld. Joining, 2018, 23: 677

[本文引用: 2]

Chen S J, Han Y, Jiang X Q, et al.

Study on in-situ material flow behaviour during friction stir welding via a novel material tracing technology

[J]. J. Mater. Process. Technol., 2021, 297: 117205

[本文引用: 1]

Kumar R, Lal A K, Bharti R P, et al.

Experimental and computational analyses of material flow characteristics in friction stir welding

[J]. Int. J. Adv. Manuf. Technol., 2021, 115: 3011

[本文引用: 1]

Colegrove P A, Shercliff H R.

3-dimensional CFD modelling of flow round a threaded friction stir welding tool profile

[J]. J. Mater. Process. Technol., 2005, 169: 320

[本文引用: 3]

Bastier A, Maitournam M H, Dang Van K, et al.

Steady state thermomechanical modelling of friction stir welding

[J]. Sci. Technol. Weld. Joining, 2006, 11: 278

[本文引用: 1]

Liechty B C, Webb B W.

Modeling the frictional boundary condition in friction stir welding

[J]. Int. J. Mach. Tools Manuf., 2008, 48: 1474

[本文引用: 1]

Atharifar H, Lin D C, Kovacevic R.

Numerical and experimental investigations on the loads carried by the tool during friction stir welding

[J]. J. Mater. Eng. Perform., 2009, 18: 339

[本文引用: 1]

Yu Z Z, Zhang W, Choo H, et al.

Transient heat and material flow modeling of friction stir processing of magnesium alloy using threaded tool

[J]. Metall. Mater. Trans., 2012, 43A: 724

[本文引用: 1]

Chen G Q, Feng Z L, Zhu Y C, et al.

An alternative frictional boundary condition for computational fluid dynamics simulation of friction stir welding

[J]. J. Mater. Eng. Perform., 2016, 25: 4016

[本文引用: 7]

Yang C L, Chen G Q, Qiao J N, et al.

Material flow during dissimilar friction stir welding of Al/Mg alloys

[J]. Int. J. Mech. Sci., 2024, 272: 109173

[本文引用: 3]

Tamadon A, Pons D J, Chakradhar K, et al.

3D-printed tool shoulder design for the analogue modelling of bobbin friction stir weld joint quality

[J]. Adv. Mater. Sci., 2021, 21: 27

[本文引用: 4]

Schneider J A, Nunes Jr A C.

Characterization of plastic flow and resulting microtextures in a friction stir weld

[J]. Metall. Mater. Trans., 2004, 35B: 777

[本文引用: 1]

Liu F C, Hovanski Y, Miles M P, et al.

A review of friction stir welding of steels: Tool, material flow, microstructure, and properties

[J]. J. Mater. Sci. Technol., 2018, 34: 39

[本文引用: 2]

Liu X C, Sun Y F, Nagira T, et al.

Experimental evaluation of strain and strain rate during rapid cooling friction stir welding of pure copper

[J]. Sci. Technol. Weld. Joining, 2019, 24: 352

Wu C S, Wang T, Su H.

Material flow velocity, strain and strain rate in ultrasonic vibration enhanced friction stir welding of dissimilar Al/Mg alloys

[J]. J. Manuf. Process., 2022, 75: 13

Qiao J N, Shi Q Y, Chen S J, et al.

Modeling and simulation of interface friction and material flow behavior during friction stir welding

[J]. Electr. Weld. Mach., 2023, 53(3): 15

[本文引用: 1]

乔俊楠, 史清宇, 陈树君 .

搅拌摩擦焊界面摩擦及材料流动行为仿真

[J]. 电焊机, 2023, 53(3): 15

[本文引用: 1]

Zhen Y Q, Liu X C, Shen Z K, et al.

State-of-art of experimental characterization of material flow in friction stir welding

[J]. J. Mech. Eng., 2020, 56(6): 184

[本文引用: 2]

甄云乾, 刘小超, 申志康 .

搅拌摩擦焊材料流动的试验表征研究现状

[J]. 机械工程学报, 2020, 56(6): 184

[本文引用: 2]

Wang X, Gao Y F, Liu X, et al.

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

[J]. Acta Mater., 2021, 213: 116969

[本文引用: 1]

Yang C L, Dai Q L, Shi Q Y, et al.

Flow-coupled thermo-mechanical analysis of frictional behaviors at the tool-workpiece interface during friction stir welding

[J]. J. Manuf. Process., 2022, 79: 394

[本文引用: 1]

Xue W, Xiao L Y, Tao T, et al.

A novel method of lateral tool-tilt for optimizing material flow and temperature distribution in friction stir welding

[J]. Mater. Today Commun., 2024, 40: 109962

[本文引用: 1]

Qiao J N, Shi Q Y, Wu C S, et al.

Elucidation of solid-state metal flow behaviors during friction stir welding: Numerical and experimental investigation

[J]. Phys. Fluids, 2023, 35: 123105

[本文引用: 1]

Yu M, Li W Y, Li J L, et al.

Modelling of entire friction stir welding process by explicit finite element method

[J]. Mater. Sci. Technol., 2012, 28: 812

[本文引用: 2]

Sahu M, Paul A, Ganguly S.

Influence of frictional heat spread pattern on the formation of intermetallic layers at the dissimilar FSW joint interface between AA 5083 and HSLA steel

[J]. J. Manuf. Process., 2022, 83: 555

[本文引用: 1]

Xu W F, Liu J H, Zhu H Q.

Numerical simulation of thermal field of friction stir welded 2219 aluminum alloy thick plate

[J]. Trans. China Weld. Inst., 2010, 31(2): 63

[本文引用: 1]

徐韦锋, 刘金合, 朱宏强.

2219铝合金厚板搅拌摩擦焊接温度场数值模拟

[J]. 焊接学报, 2010, 31(2): 63

[本文引用: 1]

Wan S Q, Wu Y X, Gong H, et al.

Thermal mechanical coupling simulation of temperature and residual stress in friction stir welding of 2219 aluminum alloy

[J]. Hot Work. Technol., 2019, 48(13): 159

[本文引用: 1]

万胜强, 吴运新, 龚 海 .

2219铝合金搅拌摩擦焊温度与残余应力热力耦合模拟

[J]. 热加工工艺, 2019, 48(13): 159

[本文引用: 1]

Ke L M, Pan J L, Xing L, et al.

Sucking-extruding theory for the material flow in friction stir welds

[J]. J. Mech. Eng., 2009, 45(4): 89

[本文引用: 1]

柯黎明, 潘际銮, 邢 丽 .

搅拌摩擦焊焊缝金属塑性流动的抽吸-挤压理论

[J]. 机械工程学报, 2009, 45(4): 89

[本文引用: 1]

He C S, Qie M F, Zhang Z Q, et al.

Effect of axial ultrasonic vibration on metal flow behavior during friction stir welding

[J]. Acta Metall. Sin., 2021, 57: 1614

[本文引用: 1]

何长树, 郄默繁, 张志强 .

轴向超声振动对搅拌摩擦焊过程中金属流动行为的影响

[J]. 金属学报, 2021, 57: 1614

[本文引用: 1]

Jain R, Pal S K, Singh S B.

A study on the variation of forces and temperature in a friction stir welding process: A finite element approach

[J]. J. Manuf. Process., 2016, 23: 278

[本文引用: 1]

Ou L, Sun B, Wang Z.

Flow stress of 2219 aluminium alloy during hot compression deformation

[J]. Hot Work. Technol., 2008, 37(2): 42

[本文引用: 1]

欧 玲, 孙 斌, 王 智.

2219铝合金热压缩变形流变应力

[J]. 热加工工艺, 2008, 37(2): 42

[本文引用: 1]

Wang Z T, Tian R Z. Handbook of Aluminum and Aluminum Alloys [M]. 2nd Ed., Changsha: Central South University Press, 2000: 188

[本文引用: 1]

王祝堂, 田荣璋. 铝合金及其加工手册 [M]. 第2版. 长沙: 中南大学出版社, 2000: 188

[本文引用: 1]

Davis J R. Aluminum and Aluminum Alloys [M]. Materials Park: ASM International, 1993: 68

[本文引用: 1]

Colegrove P A, Shercliff H R.

CFD modelling of friction stir welding of thick plate 7449 aluminium alloy

[J]. Sci. Technol. Weld. Joining, 2006, 11: 429

[本文引用: 1]

Su H, Wu C S, Bachmann M, et al.

Numerical modeling for the effect of pin profiles on thermal and material flow characteristics in friction stir welding

[J]. Mater. Des., 2015, 77: 114

Chen G Q, Wang G Q, Shi Q Y, et al.

Three-dimensional thermal-mechanical analysis of retractable pin tool friction stir welding process

[J]. J. Manuf. Process., 2019, 41: 1

Chen G Q, Liu X, Qiao J N, et al.

Improved analytical model for thermal softening in aluminum alloys form room temperature to solidus

[J]. Materials, 2023, 16: 7358

[本文引用: 1]

/