搅拌摩擦焊是一种在机械力和摩擦热作用下的固相连接方法，其通过搅拌头的旋转来使两种相同或者不同材料连接在一起。搅拌摩擦焊焊接原 理如图1所示，焊接过程中，搅拌头旋转着缓慢插入被焊工件，搅拌头和被焊接材料之间的摩擦剪切阻力产生了摩擦热，使搅拌头邻近区域的材料接近焊接材料的熔 点；当搅拌头旋转着向前移动时，搅拌头附近的金属材料从搅拌头的前沿向后沿转移，并且在搅拌头轴肩与工件表层摩擦产热和压力共同作用下，形成致密固相连接 接头[1]。

在工件与搅拌头的材料都一定的情况下，温度场主要是通过搅拌头与工件的摩擦、材料的流动和材料的塑性变形形成。材料的流动主要取决于焊接速度、旋 转速度和压力等，这些因素又决定了材料的焊接质量和微观组织。因此，很多学者研究焊接仿真模型去优化搅拌摩擦焊的焊接工艺条件以得到更好的焊接质量。

搅拌摩擦焊的模拟发展概述
搅拌摩擦焊模拟技术的发展经历了十几年的时间，中外很多专家学者对这种技术的发展做出了贡献。由于材料流动对于材料焊接性能的影响十分重要，人们 希望能预测搅拌头附近的流场。最开始人们是想通过某种方法使得搅拌头周围的流场可视化，很多人也对此进行了大量的研究。1999年，K.Colligan 在6061-T6与7075-T6铝合金焊接材料中加入钢珠，使用X射线分析得出搅拌摩擦焊是由搅拌与材料挤压两种现象的组合[2]。同 年，Reynolds等人进行了类似的实验，在焊缝位置把5454铝合金中嵌入到2195铝合金中，观察了不同焊接速度和旋转速度对搅拌头两侧材料流动的 影响[3]。此外也有许多人利用不同材料对流场进行分析。通过这种方法，Murr等人向人们展示了焊接参数变化引起的复杂流场[4]。Guerra等人使 用在6061铝合金的焊缝上放置铜箔的方法，在焊接过程达到稳态后，停止搅拌头旋转，然后使用金相分析，展示了搅拌头附近流场的区域，旋转区与过渡区 [5]。此外还有很多流场分析的方法, 金相分析可以补充X射线图像分析。

搅拌摩擦焊过程的数值模拟的主要问题是使用等效热流密度方程来代替热输入。这个等效热流密度需要通过实验，实际测量来获得。但是如果可以通过热流 模型来预测由于塑性变形和摩擦生成的热流密度的方法，就可以代替试验，于是人们研究开发了很多流体模型。徐和邓等人依据ABAQUS/Explicit开 发了二维和三维有限元模型，他们使用修改的库伦法以及工件和搅拌头之间的滑移界面的摩擦接触，预测的结果与试验结果比较接近[6]。Colegrove使 用CFD软件FLUENT模拟了2024铝合金与7075铝合金焊接过程中各种形状的搅拌头周围的流场[7]。Seidel等人也使用FLUENT2D稳 态方法模拟了圆柱状搅拌头附近的流场。Ulysse使用FIDA软件模拟了的三维热流模型过程考虑了搅拌头螺纹引起的速度[8]。Nikiforakis 使用有限差分法模拟了搅拌摩擦焊过程[9]，当然热输入全部是由塑性变形计算而来。

morfeo原理介绍

基于上述内容，我们看到人们对搅拌摩擦焊的模拟作了很多的研究，但是真正能用于指导生产以及对实际工艺参数起到优化作用的工作并不是很多。 morfeo有限元仿真软件是比利CENAERO公司开发的模拟软件，它可以模拟搅拌摩擦焊过程中的温度场、流动以及应力场。模型采用交错热流耦合分析， 可以使用欧拉法及拉格朗日法；材料本构采用热粘塑性；温度场计算考虑塑性变形生成的热量；流场结果与一系列的实验结果进行了比较；试验采用铜箔与铜线的方 式，嵌入到焊缝周围的各种位置；焊接后在各个截面上进行金相分析。在某实验参数下，可以观察到板厚各个深度下的铜箔分布如图2所示。

1 计算模型

计算模型采用欧拉法，基于质量守恒方程和动量守恒方程。模型在计算过程中，采用了一定程度的简化。首先，因为大多数的变形为粘塑性变形，假定材料 为不可压缩的，在搅拌头附近弹性变形相比粘塑性变形非常小，因此在本模型里面不进行考虑。其次，此计算采用准稳态分析，所以时间的导数在守恒方程中可以删 掉。再次，雷诺数非常小，大约10-4，因此，动量守恒方程的对流项可以被忽略。

这样守恒方程就大大简化了：∇·v=0，∇·σ=0，ρCpv·∇T=∇·k∇+Q。（1）

这些方程中未知量为速度、压力和温度。由质量守恒方程与动量守恒方程组成方程组求解，能量守恒方程单独求解。
质量守恒方程和动量守恒的方程组：∇·v = 0，∇·σ=0 ，（2）

离散格式为∫s：ε^*dΩ-∫p▽•ν^*dΩ-∫f•ν^*dT=0，∫s：ε^*dΩ-∫p▽•νdΩ=0，其中，s是偏应力；ε^∗是测试应变率；v^*为测试速度；p为压力；f为外部力。