卓金管件有限公司厚板镁合金MIG焊焊接残余应力场数值模拟的学习体会
研究背景与镁合金焊接的特殊性
吴志生等人于2012年在《焊接技术》发表的这篇研究，针对10 mm厚AZ40M镁合金板进行了MIG焊焊接残余应力场的有限元数值模拟。镁合金以其优异的比强度、比刚度和电磁屏蔽性能，在航空航天、3C电子、汽车轻量化等领域具有广泛应用前景。然而，镁合金焊接面临诸多技术挑战：
- 低熔点与高导热性：镁合金熔点约650℃，导热系数高，焊接热输入窗口窄，易产生气孔和热裂纹。
- 化学活性强：高温下镁与氧、氮反应剧烈，保护气氛要求高，否则焊缝氧化严重。
- 热膨胀系数大：约为钢的2倍，焊后残余应力和变形显著。
- 应力腐蚀敏感性：残余拉应力与腐蚀介质协同作用，加速应力腐蚀开裂（SCC）。
因此，焊接残余应力的精确预测和控制对镁合金焊接结构的服役安全性至关重要。
数值模拟方法与模型构建
文献采用有限元方法（FEM）进行焊接温度场和残余应力场的耦合计算，关键技术要点如下：
1. 热源模型：采用双椭球移动热源（Double Ellipsoid Moving Heat Source）模拟MIG焊电弧热输入。双椭球热源模型将电弧热源分为前椭球（熔化区）和后椭球（热影响区），通过调整椭球参数（a、b、c、f、r）来模拟电弧热输入的不对称分布特征。
2. 材料模型：考虑材料物理性能参数（密度、比热、导热系数、弹性模量、热膨胀系数）随温度的非线性变化。AZ40M镁合金在室温至熔化温度范围内，各项物理性能均呈现显著的温度依赖性，忽略这种非线性将导致计算结果严重偏离实际。
3. 单元技术：采用生死单元技术（Element Birth and Death）模拟焊接过程中材料的熔化和凝固。当单元温度超过熔点时"杀死"该单元（移除刚度矩阵），温度降至凝固温度以下时"激活"该单元，以此模拟焊缝金属的填充过程。
4. 耦合方法：采用间接法（Sequential Coupling）进行温度场-应力场的耦合计算。先独立求解温度场，再将温度场结果作为热载荷施加到应力场计算中，通过热膨胀和热弹塑性本构关系计算残余应力。
残余应力分布特征分析
文献给出的残余应力分布结果具有重要的工程指导意义：
| 应力方向 | 最大拉应力位置 | 最大压应力位置 | 应力特征 |
|---|---|---|---|
| 横向（垂直于焊缝） | 距焊缝中心线30 mm处 | 焊缝两端 | 焊缝区为压应力 |
| 纵向（沿焊缝方向） | 距焊缝中心线15~35 mm处 | 距焊缝中心线45~70 mm处 | 近焊缝区拉应力，远场压应力 |
从残余应力分布规律来看，横向最大拉应力出现在距焊缝中心线30 mm处，而非焊缝表面。这是因为焊缝区在冷却收缩过程中受到周围母材的约束，焊缝本身处于压应力状态，而约束区域（距焊缝一定距离处）承受拉应力。这一分布特征与经典的焊接残余应力理论一致，但具体数值和位置受板厚、焊接参数和拘束条件的影响。
纵向最大拉应力出现在距焊缝中心线15~35 mm的带状区域，反映了焊接过程中热影响区的塑性变形和约束效应。纵向最大压应力出现在距焊缝中心线45~70 mm处，这是焊接热循环对远场材料产生的热膨胀和松弛效应所致。
验证与可靠性分析
文献对最大残余应力区域的关键点进行了X射线衍射法（XRD）或盲孔法（Hole Drilling）的实测验证，测量结果与计算结果相符，残余应力循环曲线与实际理论相符。这一验证步骤对于数值模拟的可信度至关重要。
从数值模拟的可靠性角度，以下因素需要重点关注：
- 网格密度：焊缝区及热影响区需要足够细密的网格以捕捉应力梯度，但网格过密会增加计算成本。通常焊缝区网格尺寸应小于焊丝直径的1/2。
- 材料参数精度：温度相关物理性能参数的准确性直接影响计算结果。镁合金的高温物理性能数据相对匮乏，需要实验测定或可靠文献引用。
- 本构模型选择：弹性-理想塑性模型适用于残余应力的初步预测，但对于大变形和复杂载荷状态，弹塑性本构模型更为准确。
- 边界条件：夹具约束、对称边界和远场自由的边界条件设定需与实际焊接条件一致。
工程应用与残余应力控制策略
基于数值模拟结果，镁合金厚板焊接的残余应力控制可采用以下策略：
| 控制策略 | 原理 | 适用场景 | 效果评估 |
|---|---|---|---|
| 焊前预热 | 降低温度梯度，减少热应力 | 厚板、拘束度大的结构 | 降低残余应力10%~30% |
| 多层多道焊 | 后续焊道对前序焊道产生松弛效应 | 厚板对接焊 | 降低残余应力20%~40% |
| 锤击焊缝 | 引入塑性变形松弛应力 | 表面残余拉应力区 | 降低残余应力30%~50% |
| 焊后热处理 | 再结晶消除残余应力 | 可热处理的材料 | 降低残余应力50%~80% |
| 振动时效 | 低频振动使材料产生微塑性变形 | 大型结构件 | 降低残余应力30%~60% |
对于AZ40M镁合金，焊后热处理需特别谨慎，因为镁合金在高温下强度显著下降，热处理温度过高可能导致力学性能退化。振动时效和锤击法是不涉及高温处理的残余应力控制手段，更适合镁合金焊接结构。
个人反思与技术展望
这篇文献的数值模拟方法在焊接残余应力预测领域具有代表性。有限元数值模拟已成为焊接工艺开发和质量控制的重要工具，能够预测实际焊接难以直接测量的残余应力分布，为焊接工艺优化和结构设计提供指导。
然而，数值模拟结果的可信度取决于模型假设的合理性和输入参数的准确性。在工程实践中，数值模拟结果应与实验测量数据进行交叉验证，避免过度依赖模拟结果。此外，焊接残余应力的测量方法（XRD、盲孔法、中子衍射、超声法）各有优缺点，选择合适的测量方法对验证模拟结果至关重要。
从技术发展趋势看，焊接残余应力的预测和控制正朝着多物理场耦合、实时在线监测和智能优化的方向发展。结合数据分析等数据驱动方法，基于历史焊接数据和实测残余应力建立预测模型，可能实现残余应力的快速预测和工艺参数的智能优化。
总体而言，该文献以严谨的数值模拟方法揭示了AZ40M镁合金厚板MIG焊的残余应力分布规律，并通过实验验证了模拟结果的可靠性。对于从事镁合金焊接结构设计和制造的技术人员，该文献提供了重要的残余应力预测和控制依据，具有直接的工程参考价值。