焊接有限元分析：残余应力预测与焊缝强度评估实战

一、背景与需求

焊接有限元分析是评估焊接残余应力和变形的定量工具。焊接过程中的局部高温加热和不均匀冷却在结构中产生残余应力场，这些应力在服役中与外载荷叠加，可能导致疲劳寿命降低、应力腐蚀开裂或焊接变形超差。焊接有限元分析通过模拟焊接热过程和随后的热应力演化，预测残余应力分布和焊接变形量，为焊接工艺优化和结构设计提供指导。

焊接有限元分析的复杂性来自多物理场耦合。焊接过程涉及传热（热传导+对流+辐射）、冶金相变（奥氏体-铁素体转变伴随体积变化）、力学（热膨胀+相变塑性+材料性能随温度变化）三者的耦合。钢材从室温到熔点（约1500°C），弹性模量从206 GPa降到接近0，屈服强度从345 MPa降到几十MPa，热膨胀系数从1.2×10⁻⁵变到1.6×10⁻⁵/°C。这些参数的剧烈变化使焊接过程的数值模拟成为一个高度非线性问题。

做过一个压力容器接管焊接的残余应力分析项目。容器材质Q345R，壁厚32mm，接管外径219mm，壁厚8mm。焊接工艺：氩弧焊打底+埋弧焊填充盖面，预热温度150°C，层间温度200°C。初始分析用双椭球热源模型，但残余应力与实测偏差达45%——熔池尺寸在模拟中偏大，导致热影响区（HAZ）过宽。问题出在热源参数直接套用文献值而未做标定。重新用焊缝横截面金相实测熔池尺寸（宽12mm、深8mm）标定双椭球热源参数（前后半轴a1=3mm、a2=12mm、宽b=6mm、深c=4mm）后，HAZ宽度模拟值8.2mm与实测8.0mm偏差2.5%，残余应力预测值与X射线衍射实测偏差降至12%以内。

二、核心原理

焊接传热模型

焊接有限元分析的第一步是传热分析。焊接热源在移动过程中将热量输入材料，形成熔池和热影响区。热传导方程：

ρc(∂T/∂t) = ∇·(k∇T) + Q

其中Q是热源项，k是导热系数（随温度变化），ρ是密度，c是比热（随温度变化）。钢的导热系数在20°C时约50 W/(m·K)，在800°C时降到约28 W/(m·K)——高温下传热变慢，热梯度更陡。

边界条件：焊缝表面有对流换热（h≈10-20 W/(m²·K)）和辐射换热（ε≈0.8）。辐射在高温下占主导——800°C以上辐射散热量超过对流。

热源模型

焊接热源的模拟是焊接有限元分析的核心。常用热源模型：

双椭球热源（Goldak模型）：最广泛使用的熔焊热源模型。前半球（焊接方向）热流密度高、尺寸小，后半球热流密度低、尺寸大。需要标定4个几何参数（a1, a2, b, c）和热输入功率Q。

高斯面热源：适用于激光焊、电子束焊等高能束焊接。热流密度按高斯分布衰减，需要标定有效半径和峰值功率。

体积热源：适用于深熔焊（如激光深熔焊的匙孔效应），热源分布在一定体积内而非表面上。

热源标定是焊接仿真中最耗时的环节。标准方法：焊接一道试验焊缝，切取横截面做金相，测量熔池宽度和深度。调整热源参数使模拟的熔池形状与实测吻合。对于多层多道焊，每道焊缝可能需要单独标定——因为坡口形状变化影响热流分布。

顺序耦合热应力分析

焊接有限元分析通常用顺序耦合（而非完全耦合）：先做传热分析得到温度场历程，然后将温度场作为输入做热应力分析。顺序耦合适用于热-力单向耦合（温度场影响应力场，但应力场不反向影响温度场）的情况——对于大多数焊接问题这是合理的近似。

热应力分析中，总应变 = 弹性应变 + 塑性应变 + 热应变 + 相变应变：

ε_total = ε_elastic + ε_plastic + ε_thermal + ε_phase

热应变：ε_thermal = α(T) × ΔT。热膨胀系数α随温度变化，需要输入温度相关的α(T)曲线。

相变应变：奥氏体-铁素体转变伴随体积膨胀（约0.4%），相变塑性在相变温度区间产生额外应变。相变应变的影响在低合金钢中不可忽略——忽略相变应变会使残余应力预测偏差15-25%。

材料性能的温度依赖性

焊接有限元分析需要输入材料性能从室温到熔点的温度曲线。关键参数：

弹性模量E(T)：20°C→206 GPa，600°C→约120 GPa，1000°C→约20 GPa，1500°C→接近0。

屈服强度σ_y(T)：20°C→345 MPa，600°C→约100 MPa，1000°C→约20 MPa。

热膨胀系数α(T)：20°C→1.2×10⁻⁵/°C，800°C→1.6×10⁻⁵/°C。

导热系数k(T)：20°C→50 W/(m·K)，800°C→28 W/(m·K)。

这些参数通常从材料数据库（如JMatPro、ANSYS材料库）或实验获取。如果高温段数据缺失（>1000°C），可以外推但需注意——熔池区域E≈0，如果仍用室温值会导致刚度矩阵病态。

在焊接有限元分析中，材料高温性能数据的准确性是残余应力预测精度的主要决定因素——特别是600-1200°C区间的屈服强度曲线，直接决定冷却后的残余应力水平。

三、关键技术要点

单元生死技术（Element Birth/Death）

焊接过程中焊缝金属是逐步填充的——初始模型中焊缝区域不存在，随焊接进行逐步”激活”。ABAQUS用”Model Change”功能实现：焊缝单元初始处于”死”状态（刚度乘以极小因子），在对应焊接时刻”激活”。

激活时间：按焊接速度计算。如焊接速度5 mm/s，焊缝长100mm，总焊接时间20秒。焊缝沿焊接方向分为20个单元段，每秒激活一段。激活温度设为熔点（1500°C）或稍低（1300°C）。

热源移动模拟

移动热源用DFLUX子程序（ABAQUS）实现。子程序在每个增量步计算热源当前位置，在该位置附近施加双椭球热流。热源移动的精度取决于时间步长——时间步太大热源”跳跃”大，温度场不连续。标准做法：时间步 ≤ 热源半径/焊接速度 × 0.2，如热源半径6mm、速度5mm/s，时间步≤0.24s。

冷却阶段模拟

焊接后冷却至室温是残余应力形成的关键阶段。冷却分两段：(1)焊后缓冷（从层间温度200°C到环境温度20°C），约30-60分钟；(2)自然冷却至室温。

冷却阶段的残余应力演化：当焊缝和HAZ冷却到相变温度（约700°C）时，奥氏体-铁素体相变产生体积膨胀，抵消部分热收缩产生的拉应力。冷却到室温后，焊缝中心通常为拉应力（接近屈服强度），远离焊缝处为压应力——形成残余应力自平衡场。

多层多道焊模拟

实际焊接通常是多层多道焊。每道焊缝的顺序激活模拟：道1→冷却→道2→冷却→…→最后一道→冷却到室温。每道焊缝的热源参数可能不同（坡口形状变化），每道之间的层间温度需要控制。

多层多道焊的计算量随焊道数线性增加——8道焊缝的计算时间约为单道焊的8倍。优化策略：将相邻焊道合并为一个”等效焊道”，减少激活步数。但合并会损失层间温度和道间残余应力的细节。

网格策略

焊接区域的网格密度要求：焊缝中心到HAZ边缘（约5-10mm）需要细网格（0.5-1mm），因为温度梯度和残余应力梯度在此区域最陡。远离焊缝的区域可以用粗网格（5-10mm）。

温度梯度方向上的网格要求：沿板厚方向，焊缝处至少4-6层单元，因为温度在厚度方向衰减快。沿焊接方向，一个热源半径内至少3-5个单元。

四、实操流程

以压力容器接管焊接为例。

第一步：传热分析标定

切取试验焊缝横截面做金相，测量熔池形状：宽12mm、深8mm。调整双椭球热源参数：前半椭球a1=3mm，后半椭球a2=12mm，宽b=6mm，深c=4mm，热输入功率Q=1800W（热效率75%×电弧功率2400W）。

验证：模拟的熔池宽11.8mm、深7.9mm，与实测偏差<3%。HAZ（峰值温度>723°C区域）宽度8.2mm，实测8.0mm。

第二步：热分析

模型：接管与容器相交区域用细网格（焊缝0.5mm，HAZ 1mm，远场5mm），总单元约3.5万。3道焊缝（打底+填充2道），每道焊缝长500mm，焊接速度4 mm/s，每道焊接时间125s。

温度相关材料参数：k(T)从50 W/(m·K)（20°C）到28 W/(m·K)（800°C）到15 W/(m·K)（1500°C）；c(T)从460 J/(kg·K)（20°C）到670 J/(kg·K)（800°C）。

边界条件：空气对流h=15 W/(m²·K)，辐射ε=0.85（Stefan-Boltzmann常数σ=5.67×10⁻⁸ W/(m²·K⁴)）。预热温度150°C施加为初始条件。

时间步：焊接阶段0.2s，冷却阶段逐步放大到10s。总分析时间从焊接开始到冷却至室温，约3小时物理时间。

第三步：热应力分析

将传热分析的温度场历程导入热应力分析。材料力学性能温度曲线：E(T)从206 GPa（20°C）到20 GPa（1000°C）；σ_y(T)从345 MPa（20°C）到15 MPa（1000°C）；α(T)从1.2×10⁻⁵（20°C）到1.6×10⁻⁵（800°C）。

相变模型：用ABAQUS的TRIP（Transformation Induced Plasticity）选项。奥氏体-铁素体相变温度区间800→500°C，体积膨胀0.4%。

边界条件：容器两端面约束轴向位移（防止刚体运动），允许径向膨胀。

第四步：残余应力结果分析

焊缝中心线纵向残余应力（沿焊接方向）：峰值+285 MPa（拉应力），接近Q345R屈服强度的82%。这个拉应力区域沿焊缝中心约8mm宽，在HAZ边缘转变为压应力（-80 MPa），在距焊缝30mm处趋近于零。

焊缝中心线横向残余应力（垂直焊接方向）：峰值+165 MPa（拉应力），低于纵向残余应力——这是焊接残余应力的典型特征，纵向残余应力通常大于横向。

接管与容器交界处的残余应力集中：由于几何不连续（接管开孔削弱了容器壁），残余应力在此处升高至+310 MPa，超过屈服强度的90%。这个区域是疲劳和应力腐蚀的高风险区。

第五步：实验对标

用X射线衍射法测量焊缝表面残余应力。测量点：焊缝中心、HAZ中心（距焊缝4mm）、距焊缝15mm、30mm。

对比结果：

• 焊缝中心：实测+268 MPa，模拟+285 MPa，偏差6.3%
• HAZ中心：实测+82 MPa（压应力），模拟-80 MPa，偏差3%
• 距焊缝15mm：实测-75 MPa，模拟-65 MPa，偏差13%
• 距焊缝30mm：实测-18 MPa，模拟-12 MPa，偏差33%（绝对值小，相对偏差大但影响不大）

四个位置中三个偏差<15%，残余应力场预测可信度较高。

在焊接有限元分析中，X射线衍射或盲孔法实测对标是验证残余应力预测的标准手段——不经验证的残余应力分析结果不能直接用于工程设计决策。

第六步：焊接变形预测

接管焊接后的变形：接管端面翘曲0.35mm（向焊接侧偏转），容器壁局部凹陷0.12mm。变形量在允许范围内（机加工余量0.5mm可以覆盖）。

消除残余应力方案分析：模拟焊后退火（580°C×2h随炉冷却），退火后残余应力峰值从285 MPa降至65 MPa（降幅77%）。但退火变形增大——接管端面翘曲增加到0.48mm，仍可接受。

五、常见问题与排查

残余应力峰值远高于屈服强度

检查高温材料参数输入。如果1000°C以上屈服强度设为室温值（345 MPa而非15 MPa），冷却时塑性变形不足，残余应力被严重高估。正确做法：高温段屈服强度按实际曲线输入，熔池区域（>1400°C）屈服强度设为接近0。

熔池形状与实测不符

热源参数未标定。直接用文献中的双椭球参数通常不适用——不同焊接工艺（电流、电压、速度、保护气体）的热效率差异大。必须用同工艺的试验焊缝金相标定热源几何参数。

相变应变被忽略

许多分析忽略了奥氏体-铁素体相变应变，导致焊缝中心残余应力偏高15-25%。检查是否启用了TRIP模型。对于低碳钢和低合金钢，相变应变不可忽略。对于不锈钢（无相变），可以忽略。

计算时间过长

焊接仿真的计算时间通常以天计。优化策略：(1)用对称面减少模型规模——如果焊缝和结构有对称性，用1/2或1/4模型；(2)简化焊道——将多层多道焊合并为等效焊道（损失精度但减少计算量50-80%）；(3)用显式方法替代隐式方法——对于大变形焊接（薄板焊接变形），显式算法的效率可能高于隐式。

冷却阶段不收敛

高温时材料接近液态（E≈0），刚度矩阵可能病态。解决方法：设定一个最低刚度（如E_min=1 MPa），确保刚度矩阵非奇异。或在高温区域用”annealing”功能——超过指定温度后应力归零，模拟材料重结晶。

六、复盘总结

焊接有限元分析最关键的经验：热源参数标定是分析质量的基石。不做标定直接用文献参数，熔池尺寸可能偏差30-50%，导致热影响区宽度错误，残余应力预测偏差可达40%以上。花半天时间做一道试验焊缝、切金相、标定热源参数，是投入产出比最高的工作。其次是材料高温性能曲线的完整性——600-1500°C区间的E(T)、σ_y(T)、α(T)数据缺一不可，这些数据决定了冷却过程中塑性变形的累积量，进而决定残余应力水平。

方法局限：焊接有限元分析的准确性受限于物理模型的完备性。标准分析通常忽略熔池流体动力学（对流对熔池形状的影响）、金属蒸汽和等离子体效应。对于电子束焊和激光焊，这些效应不可忽略，需要用更复杂的多物理场模型。残余应力测量方法的局限性——X射线衍射只能测量表面残余应力，内部残余应力需要用中子衍射或深孔法，成本高且可获得性有限。对于大型焊接结构（如船体分段），全尺寸焊接仿真的计算量过大，需要用简化模型（如固有应变法）替代热弹塑性分析——固有应变法将焊接残余应力等效为初始应变施加，计算量减少两个数量级，但精度依赖于实验标定的固有应变数据库。