做焊接仿真这些年连续驱动摩擦焊一直是我觉得特别值得聊的一个方向。它不是靠熔池结晶而是靠摩擦热把界面材料加热到红硬状态后顶锻成型温度、应力、变形全部挤在毫米级范围内相互作用。ABAQUS 里的二维轴对称热力耦合模型是我解决这类问题的第一选择计算量比三维少一到两个数量级却能抓住温度峰值、轴向缩短量、飞边形态和残余应力这些工艺最关心的结果。这套方法适合正在调连续驱动摩擦焊工艺参数的工程师也适合刚接触焊接仿真的学生在前处理、求解、后处理全流程里少踩几个坑。很多人一看“摩擦焊”三个字就想着必须把旋转建出来然后一头扎进三维网格、接触、大变形机器跑一晚上最后还是不知道参数为什么不对。其实连续驱动摩擦焊的二维轴对称模型并不直接模拟旋转运动而是用等效摩擦热源把旋转带来的能量注入界面。这个“等效”不是拍脑袋后面我会把热流公式、边界条件、材料参数和排错方法全部展开。1. 连续驱动摩擦焊仿真的整体思路拆解1.1 为什么二维轴对称模型是前期最优解连续驱动摩擦焊的几何结构本身是绕轴旋转的两个工件同轴对中一个夹在主轴侧旋转一个固定。摩擦界面的热量主要来自周向相对滑动但这段滑动在周向是均匀分布的。只要转速足够高、摩擦过程足够稳定界面温度场在周向几乎不变化这就给二维轴对称建模提供了物理基础。二维轴对称模型把三维圆柱体压缩到 r-z 平面单元只在径向和轴向有自由度。这样做最大的好处是计算规模大幅下降。三维模型要模拟完整的摩擦界面接触至少需要几十万个六面体或四面体单元一个工况跑完动辄数小时而二维轴对称模型往往几千个单元就能覆盖同样尺寸的工件BMS边界矩很小可以在一个晚上反复试多组工艺参数。但这也不是没有代价。二维模型没有切向位移自由度旋转带来的周向剪切和扭转飞边没法直接显式表达。实际操作中我们把它等效成界面热流把旋转摩擦的机械功转化为热载荷。只要你的关注点是温度场、轴向缩短量、飞边形貌和焊后残余应力这层等效完全够用。如果非要研究飞边里的螺旋流线、周向织构那才需要退回去做全三维模型。1.2 完全热力耦合怎么选Standard还是Explicit摩擦焊的界面热量来自摩擦功摩擦功又取决于接触压力、摩擦系数和相对滑动速度同时温度一上来材料软化、屈服强度和导热系数都会变化反过来又改变接触面积和压力分布。这是一个典型的强耦合现象所以不建议用先算温度场再导入结构的顺序耦合做法。完全热力耦合求解器把温度自由度和位移自由度放在同一个增量步里求解能自然反映这种反馈关系。求解器方面我更推荐 ABAQUS/Explicit 而不是 Standard。连续驱动摩擦焊在顶锻阶段会挤出大量飞边材料应变轻松超过 200%接触状态从贴合到滑移再到分离变化非常剧烈。Standard 在这种强非线性下收敛性很差经常出现负特征值或增量步长被砍到零的尴尬局面。Explicit 用中心差分推进不存在整体刚度矩阵反复迭代的问题只要时间增量满足稳定性条件大变形、接触穿透和网格畸变都能被“硬扛”过去。需要注意显式分析的时间步受最小单元尺寸和材料波速限制而摩擦焊过程又包含了秒级的热传导和长时间的保压冷却。所以每次跑显式模型我都会特别关注质量缩放和动能/内能比这一点到第 3.3 节再细说。2. 模型搭建的关键细节2.1 几何与网格轴对称单元怎么选几何建模时直接在 ABAQUS/CAE 的 Part 模块里画矩形截面即可。以半径 10 mm 的实心圆棒对焊为例两个工件在子午面上表现为两个矩形的堆叠界面位于二者交界处。工件长度不用全部建模通常取远离摩擦面的一侧 20~50 mm 就足够因为该区域对界面温度和变形影响已经很小。太长了只会增加网格数量太短了又会让热边界条件过早影响界面一冷一热都会导致结果失真。网格类型建议选择带温度自由度的轴对称缩减积分单元 CAX4RT四节点双线性位移-热耦合单元中心轴线处可以使用三节点三角形 CAX3T 过渡。摩擦界面附近是整个模型的灵魂区域温度梯度和应变梯度都集中在这里网格必须加密。我用过的经验是界面两侧至少各留 10 层厚度约 0.1~0.3 mm 的单元往下离开界面 5 mm 后再逐步放大到 2 mm。宁可远场稀疏一点也不要吝惜界面附近的网格。大变形发生后飞边被挤出靠近界面的单元会被强烈拉长。为了不让单元畸形过度建议在 Step 模块开启 ALE 自适应网格Adaptive Mesh Domain把高应变区选为自适应区域。ALE 会不断重置网格拓扑让材料流动更平滑温度场的连续性也会好很多。二维模型中 ALE 的开销很低属于“开了不亏”的选项。2.2 材料模型与热物性参数摩擦焊涉及高温、大应变、高应变率材料本构不能只给一个常温屈服强度。我通常用 Johnson-Cook 模型描述中等碳钢的高温流动应力σ (A B ε^n)(1 C ln ε̇_r)(1 - T^*^m)其中 A、B、C、n、m 是材料参数T^* 是归一化温度T^* (T - T_r)/(T_m - T_r)。举个常见的中碳钢参数示例A300 MPaB500 MPan0.3C0.02m1.1熔点 T_m1795 K参考温度 T_r293 K。这些参数不是普适的不同批次、不同热处理状态的材料差别很大如果要做精确工程预测必须用高温霍普金森杆试验数据重新标定。热物性参数不能取常数。摩擦界面的温度可能从室温快速升到 1000°C 以上密度、比热和导热系数都会发生明显变化。我用过一个典型碳钢的温度相关热物性表给你参考温度 (°C)密度 (kg/m³)比热 (J/(kg·K))导热系数 (W/(m·K))207850450462007820510444007790580416007750680378007700780321000760090028在塑性选项里把非弹性热分数Inelastic heat fraction设成 0.9也就是塑性功的 90% 转化为热。摩擦热在二维模型中不再通过默认摩擦接触产生而是由等效热流直接注入所以摩擦生热效率通常直接体现在热流幅值上。2.3 接触与等效摩擦热两个工件之间的接触属性要同时定义力学和热学两套行为。力学方面用罚接触法向选择硬接触Hard Contact切向摩擦系数取 0.3 左右并且可以让它随温度下降——高温软化界面的剪切抗力会显著降低。热学方面需要定义“间隙导热”函数当接触间隙为 0 时界面导热系数取一个很大的值比如 1e6 W/(m²·K)当间隙增大到 1 mm 时降到 1e3 W/(m²·K)这样模拟界面的热阻随接触状态变化。现在到最关键的环节连续驱动摩擦焊的二维轴对称等效热源。设工件半径 R角速度 ω摩擦压力 p摩擦因数 μ。在半径 r 处单位面积摩擦热的表达式是q(r) μ p ω r这个式子理解起来很直观摩擦应力 μp 乘以该点的切向线速度 ωr就是单位面积的功率密度。在二维轴对称模型里把这个热流密度作为面热流施加到摩擦界面上就等效替代了旋转摩擦。举个例子R10 mm转速 n1500 rpm角速度 ω2πn/60≈157 rad/s摩擦压力 p2 MPaμ0.3。那么界面最外侧的峰值热流密度约为 0.3×2×10⁶×157×0.01≈9.4×10⁵ W/m²。整个界面的总热功率可以积分得到P ∫₀^R μ p ω r · 2πr dr μ p ω (2πR³/3)代入数值P≈0.3×2×10⁶×157×2.094×10⁻⁶≈197 W虽然总量不大但考虑到只有半径 10 mm 的小截面局部热流密度已经达到 MW/m² 量级完全足够让界面迅速升温。在 ABAQUS 中施加热流密度最常见的方式是写 DFLUX 子程序。下面是一个简化版的示意代码变量名和真实传参之间需要按版本调整SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS, 1 JLTYP,TEMP,PRESS,SNAME) C PARAMETER (R0.010, OMEGA157.0, MU0.3, P2.0E6) RCOORD COORDS(1) FLUX MU * P * OMEGA * RCOORD RETURN END如果你不想碰子程序也有一个“土办法”把摩擦面沿径向等分成 5~10 个环带每个环带施加一个平均半径对应的热流密度。我用过这种做法做快速参数筛选结果和 DFLUX 的误差在 5% 以内但前处理会麻烦不少。3. 边界条件与载荷设置实操3.1 工艺参数与载荷如何施加连续驱动摩擦焊的工艺过程大体分为摩擦加热、顶锻、保压冷却三个主要阶段。以一个半径 10 mm 的中碳钢实心棒摩擦焊为例我常用的一套工艺参数如下工艺阶段轴向压力 (MPa)转速 (rpm)持续时间 (s)摩擦加热215001.5顶锻501.0保压冷却5010.0几何模型里固定端工件远离界面的那侧端面施加全约束也就是轴向位移设为零旋转端工件远离界面的端面施加轴向压力压力方向指向工件内部模拟顶锻力。中心轴线上的节点要约束径向位移为零这是轴对称模型的必然要求。旋转运动本身不建模所以转速并不是一个边界条件而是右上角“等效热流幅值”的一个控制变量。在 DFLUX 子程序中把 ω 按幅值曲线读取让热流从 0 在 0.01~0.05 s 内平滑升高到稳态模拟摩擦焊初始加速阶段。直接一步加到满值名义上省事但会产生剧烈温度冲击真实材料不会这样响应。压力载荷同样建议用幅值曲线。ABAQUS 的 Smooth Step 幅值非常合适它在起点和终点的一阶导数都为零能避免加载瞬间激起高频应力波。摩擦阶段压力从 0 平滑升到 2 MPa顶锻阶段在 0.1 s 内升到 5 MPa同时热流按转速归零的幅值同步关闭。保压冷却阶段压力维持 5 MPa让接头在压力下完成固态连接和收缩。3.2 热边界条件与散热只在界面加热还不够还需要给模型的外部边界设置准确的热损失。自由表面施加对流换热空气自然对流系数取 10~20 W/(m²·K)高温阶段还要考虑热辐射辐射发射率取 0.7环境温度设 293 K20°C。这两个边界看似简单但如果不设置温度场会明显偏高特别是保压冷却阶段温度降不下来残余应力计算就会失真。工件与夹具接触的端面又是另一回事。实际摩擦焊设备中固定端通常夹持在带有水冷结构的夹具里散热能力比空气强一到两个数量级。我们虽然不建夹具实体但可以在固定端面施加等效换热系数取 500~2000 W/(m²·K)同时把环境温度设成冷却水的温度比如 293 K。旋转端与主轴夹头接触的地方也做类似处理。需要说明的是这种“等效换热”只在远离摩擦面的地方使用。摩擦界面附近的散热主要由工件内部导热承担界面热阻已经在接触属性里定义了不要再额外叠加一个强热沉边界否则会把界面温度人为拉低。3.3 求解控制与时间缩放显式求解器有一个绕不开的问题稳定时间增量由最小单元尺寸和材料波速决定。0.1 mm 的钢单元对应的稳定时间步大约只有 2×10⁻⁸ s而完整焊接过程 12.5 s如果老老实实跑需要数十亿个增量步工程上根本不可行。实际仿真中主要通过两种方式解决。第一是质量缩放把稳定时间增量放大到 1×10⁻⁶~1×10⁻⁵ s 量级。代价是惯性力被人为放大所以一定要在后处理里检查动能与内能的比值。我常用的准则是在整个摩擦加热和顶锻阶段动能/内能最好控制在 5% 以下极限不要超过 10%。如果超过说明质量缩放得过狠结果基本不能用来判断飞边形态。第二是适度放松网格把摩擦面附近的最小尺寸控制在 0.2~0.3 mm远场再放大让稳定步长不至于太小。不要随便缩短物理时间。摩擦阶段的热流持续时间和顶锻阶段的时间窗口都是真实工艺参数缩短它们等于改变焊接能量输入。如果实在太慢可以先把保压冷却这种次要阶段用更大的质量缩放因子跑或者导出中间状态改用隐式求解器继续算冷却两种方式我都试过效果都可控。4. 结果分析与后处理技巧4.1 温度场和变形场怎么读计算完成后打开 ODB第一步先看温度场云图。典型的连续驱动摩擦焊温度场会呈现“碗状”分布摩擦界面边缘温度最高中心温度略低因为热流密度 q(r) 本身就随半径增大。界面峰值温度在中碳钢摩擦焊里通常落在 1000~1300°C低于熔点 1450°C 左右。如果云图显示温度明显超过熔点先别急着调整工艺多数情况下是热流计算或网格问题。第二步看变形场和飞边形态。在轴对称剖面里摩擦界面两侧的材料被挤出形成上下对称的一对飞边。输出“轴向缩短量”可以监测焊接过程中的顶锻进给量方法是在旋转端端面上设一个 Reference Point输出该点的轴向位移历史减去初始位置就是总的轴向缩短量。9~12 mm 直径的棒材摩擦和顶锻总缩短量一般在 2~4 mm 范围。缩短量太小说明热输入不足界面还没软透缩短量过大则说明飞边过多母材损耗严重。第三步是残余应力。停止顶锻后继续保压冷却焊后残余应力主要分布在界面两侧 5~10 mm 范围内。如果后续要模拟接头强度可以把 ODB 的最终应力场作为初始状态导入到静力分析里再施加拉伸或弯曲载荷。4.2 怎么判断仿真结果是否合理仿真不是跑出云图就结束必须拿结果和物理事实对表。我常用的合理性判断有这几条界面峰值温度必须低于材料熔点但又要高到足以发生显著塑性流动和氧化物挤出。中碳钢在摩擦焊中峰值温度出现在固相线以下如果模型里出现了大范围液态区域说明热流密度设错或网格太粗。轴向缩短量要和设备上读到的实际进给量比一比差太多就回头查摩擦压力或顶锻力的幅值曲线。飞边应该对称且连续不应出现中间空洞或明显撕裂——这往往是接触设置失误或单元失效参数太激进。最后检查能量模型内部积累的热能应该约等于摩擦输入的总功减去对流和辐射散失的热量。如果热源输入几万焦耳、内能却只有几百焦耳那热量一定在某个奇怪的地方丢掉了。4.3 常见问题与排查技巧速查我整理了一份小表格按我自己的踩坑频率排序现象可能原因处理办法温度瞬间超过熔点热流密度给错单位或网格太粗检查 μpωr 量纲加密界面网格单元畸变导致增量步缩到极小顶锻变形过大ALE 区域没覆盖飞边区扩大自适应区域适当增加界面单元高度动能/内能比值过高质量缩放因子太大降低质量缩放拉长压力加载的斜坡时间界面温度出现锯齿状波动接触间隙热导设置不合理增大 gap conductance缩小间隙阈值飞边不对称或一侧开裂主从面设置不当把较硬的固定端设为主面旋转端设为从面保压冷却阶段温度下降过慢自由表面对流/辐射没设或夹具换热系数太小检查所有外表面边界这六类问题覆盖了我遇到过的绝大多数收敛和结果异常。最容易被忽视的其实是“单位”。ABAQUS 不帮你换单位你如果用 mm 建模却给了 Pa 的压力力密度会差 1000 倍热流密度用 W/m² 还是 W/mm² 也会导致温度场天差地别。我习惯把长度统一成米压力统一成 Pa热流统一成 W/m²这样至少不会在第一步就崩盘。5. 实操心得与后续扩展5.1 参数敏感性速查二维轴对称模型的参数敏感性很有规律。摩擦压力增大接触面积增速加快热流密度里又包含压力项所以温度升温会更快但同时材料被更快挤出飞边也会更大。转速主要影响热流密度分布转速越高边缘和中心的温差越悬殊中心区域的温度反而可能不足。摩擦时间是最直观的旋钮时间太短界面没热透顶锻时压不动时间太长热输入过量飞边过多且热影响区异常粗大。摩擦系数看似简单其实高温下会显著下降如果全程取常数界面温度会被高估。我建议做参数筛选时用“三步走”先固定转速和压力扫描摩擦时间再固定时间和压力扫描转速最后固定时间和转速扫压力。每次只动一个变量出来的信号才干净。二维模型算得快一晚上扫十几组完全不是问题这也是它在工艺优化里的最大价值。5.2 从二维到三维、微观组织扩展二维轴对称模型不包含周向剪切这是它的天然边界。如果要深挖飞边卷曲形态、扭转织构或者摩擦扭矩的时间历程就必须回到三维模型。三维模型可以在 ABAQUS/Explicit 中让旋转端真正转动接触属性里的摩擦生热会自动生成热量但网格量和时间成本会急剧上升。我通常先跑二维筛选出可行的压力、转速、时间窗口再选 2~3 组关键参数用三维验证效率和精度都能兼顾。后续还可以把这套模型的温度历史导给相变模型预测热影响区的晶粒尺寸和硬度分布把摩擦焊仿真从“热场”推进到“组织场”。最后分享一个我保留到现在的习惯每次跑新参数的模型前先用固定摩擦压力做一次 1~2 s 的短时算例盯着截面上最高温度是否平滑上升。如果温度曲线在第一个 0.05 s 内瞬间跳到上千度大概率是热流密度给错单位或者接触根本没有建立。二维轴对称模型的精度很依赖热流密度和材料参数这两个地方核对清楚后面通常不会翻车。