基于ANSYS平台的堆焊热应力动态模拟学习体会
文献概况与研究背景
本文发表于2004年《武汉理工大学学报(交通科学与工程版)》第28卷第6期,由武汉理工大学材料科学与工程学院王强与湖北工学院机械工程学院李冬林合作完成,页码为866至869页。文章的核心工作在于建立适当的有限元模型,对平板堆焊过程中的温度场、应力场瞬态变化以及焊后残余应力和变形进行动态模拟,并利用ANSYS软件的APDL语言编写程序实现移动焊接热源的施加。分类号为TG455,关键词涵盖温度场、应力场、ANSYS及数值模拟。
在2004年前后,国内焊接数值模拟领域正处于从理论探索向工程应用过渡的阶段。堆焊作为一种特殊的焊接工艺,其目的是在基体表面熔敷一层具有特定性能(如耐磨、耐腐蚀、耐高温)的合金层,广泛应用于石油管道内壁防腐、耐磨衬里、阀门密封面修复等工程领域。堆焊过程中,由于多层多道焊的反复热循环,残余应力分布复杂,极易导致堆焊层开裂、剥落或基体变形,严重制约了堆焊件的使用寿命。传统实验方法难以实时获取熔池内部及热影响区的温度与应力信息,而数值模拟技术则为深入理解堆焊热-力耦合过程提供了强有力的工具。
核心观点与技术路线
作者选择ANSYS作为数值模拟平台,其核心思路是利用ANSYS中热-结构耦合分析功能,通过APDL(ANSYS Parametric Design Language)参数化编程实现移动热源的动态施加。这一技术路线的关键在于将焊接热源从静态加载转化为随时间移动的动态载荷,从而真实反映焊接过程中温度场和应力场的瞬态演变规律。
在有限元模型的建立方面,作者采用了合理的网格划分策略,对焊缝区域进行了局部网格加密,以提高温度梯度和应力集中区域的计算精度。热源模型通常采用Goldak双椭圆模型或点热源模型,其中Goldak双椭圆模型能够较好地描述焊接热输入在熔池前方的集中和熔池后方的拖尾效应,参数包括熔池长度、宽度、峰值热流密度等,这些参数的确定直接决定了模拟结果的准确性。
在材料本构关系的选取上,堆焊过程涉及基体材料与堆焊层材料两种不同合金体系,其热物理性能(导热系数、比热容、密度)和力学性能(弹性模量、屈服强度、热膨胀系数)均随温度发生显著变化。作者需要建立温度相关的材料属性数据库,并在APDL程序中通过TABLE命令或BF命令定义这些非线性关系。
温度场与应力场动态模拟技术要点
堆焊过程的热-力耦合分析可分为两个阶段:第一阶段为瞬态热分析,求解焊接过程中温度场的动态变化;第二阶段为热力耦合分析,将热分析结果作为载荷施加到结构分析中,求解应力和变形。ANSYS中实现这一过程的关键技术点包括:
| 技术参数 | 典型取值范围 | 说明 |
|---|---|---|
| 热源峰值热流密度 | 10^7 ~ 10^9 W/m² | 取决于焊接工艺与热输入 |
| 焊接速度 | 5 ~ 50 mm/min | 堆焊通常取较低值以保证熔敷率 |
| 熔池长度 | 10 ~ 40 mm | Goldak模型前椭圆半长 |
| 熔池宽度 | 5 ~ 20 mm | Goldak模型椭圆半宽 |
| 网格尺寸(焊缝区) | 0.5 ~ 2.0 mm | 需满足网格收敛性要求 |
| 时间步长 | 0.1 ~ 1.0 s | 动态模拟需足够小的时间步 |
| 热对流系数 | 5 ~ 25 W/(m²·K) | 自然对流与强制对流 |
| 热辐射系数 | 0.8 ~ 0.95 | 表面发射率 |
在APDL编程实现移动热源方面,作者需要编写宏命令或用户子程序,在每一步计算中根据当前时间确定热源位置,并通过BF命令或BFUNIF命令在对应单元上施加面载荷。这一过程需要精确控制时间步长与移动速度的匹配关系,即每步热源移动距离等于焊接速度乘以时间步长。
残余应力与变形分析
堆焊后的残余应力分布是评价堆焊质量的关键指标。由于堆焊层与基体材料热膨胀系数的差异,冷却过程中堆焊层收缩受到基体的约束,在堆焊层中产生拉应力,在基体中产生压应力。当堆焊层中的拉应力超过材料的屈服强度时,将产生塑性变形,冷却完成后形成残余应力。
多层堆焊的残余应力分布具有明显的层状特征。最外层堆焊层由于最后冷却,其残余拉应力通常最大;而底层堆焊层由于后续层焊接的热回火作用,残余应力有所降低。这种应力分布规律对堆焊层开裂风险具有直接指导意义。
在变形预测方面,堆焊过程产生的变形主要包括纵向收缩变形、横向收缩变形和角变形。对于平板堆焊,角变形是由于焊缝上下表面温度分布不均匀导致的弯曲变形。ANSYS模拟可以预测焊后变形量,为实际生产中预变形补偿提供依据。
工程实践启示与技术反思
从工程实践角度来看,本文的研究具有重要的指导意义。在石油天然气管道内壁防腐堆焊、耐磨衬里堆焊等实际应用中,堆焊层开裂和剥落是长期困扰工程技术人员的质量难题。通过数值模拟预测残余应力分布,可以有针对性地优化堆焊工艺参数,如调整堆焊层数、控制层间温度、优化焊接顺序等,从而有效降低残余应力水平。
然而,2004年的ANSYS版本在热-力耦合分析能力上与现代版本相比仍有较大差距。现代ANSYS已具备更完善的用户子程序接口(如USDFLD、UDEC等),能够更精确地描述焊接过程中的材料相变、塑性流动和损伤演化。此外,自适应网格技术(ALE方法)的引入使得熔池区域的网格能够随温度场动态调整,显著提高了计算精度和效率。
在堆焊工艺优化方面,基于数值模拟的结果可以指导以下工艺改进措施:第一,采用多层多道焊并控制层间温度在200至300摄氏度之间,利用热回火效应降低残余应力;第二,在堆焊层中引入过渡层,减小基体与堆焊层之间的热膨胀系数差异;第三,优化焊接顺序,采用分段退焊或跳焊方式分散热输入。
本文的研究体现了焊接数值模拟从静态分析向动态分析、从单场分析向多场耦合分析的发展趋势。对于从事堆焊工艺研究的工程技术人员而言,掌握ANSYS等有限元分析工具,理解热-力耦合分析的基本原理,是实现堆焊工艺科学优化的必备技能。在当前的智能制造和数字孪生背景下,焊接过程的数字化仿真与实时监测将更加紧密地结合,为堆焊质量的在线控制和预测性维护提供新的技术手段。
卓金管件有限公司