厚板不锈钢凝固裂纹驱动力的三维数值模拟研究.pdf

厚板不锈钢凝固裂纹驱动力的三维数值模拟 1．哈尔滨工业大学现代焊接生产技术国家重点实验室董志波魏艳红 2．机械科学研究院哈尔滨焊接研究所 刘仁培 董祖珏 摘要：利用三维有限元网格划分技术，对工件进行网格剖分，采用双椭球形热源模型研究了SUS310不 锈钢的三维焊接温度场的分布规律。在此基础上重点研究了SUS310不锈钢在凝固裂纹敏感温度区间内， mm厚板的凝固裂纹驱动力曲线，并利用动态单元再生方法， 焊缝金属应变场的动态演变过程，得到了lO 消除了焊接凝固过程对应变场的影响。模拟计算得到的驱动力曲线与前人实验测量的凝固裂纹阻力曲线进 行了比较，初步预测了实验条件下SUS310不锈钢的凝固裂纹敏感性。 关键词：凝固裂纹，应变场，单元再生方法 前富 长期以来，焊缝金属凝固裂纹一直是焊接结构中重要焊接缺陷之一。凝固裂纹的产生与否取决 于两个方面的因素，材料凝固裂纹阻力和凝固裂纹驱动力IlJ。因此，实现凝固裂纹开裂倾向的预测需 要得到两方面的数据：材料在脆性温度区间内的凝固裂纹阻力曲线和在这一温度区间的驱动力曲线。 材料凝固裂纹阻力可以通过实验方法加以测量。利用横向可调拘束实验，Senda和Arata等【2，3】系统研 究了碳钢和不锈钢等的凝固裂纹敏感性。Goodwin 14J使用Sigmajig实验分析了薄板不锈钢的热裂纹 敏感性。Matsuda等【5】通过MISO方法提供了一系列凝固裂纹现象的基本研究结果。 由于焊接过程的复杂性，很难确定和测量焊接动态应力应变场的分布，因此数值模拟方法对研究 焊缝在凝固过程中所受的应力应变及最终获得凝固裂纹驱动力曲线将具有十分重要的意义。Brooks[6l 纹驱动力曲线的模拟问题进行了初步研究。然而二维的简化对实际的三维构件势必产生一定的误差， 因此三维模型才能真正体现实际的焊接过程。 本文在研究焊接温度场的基础上，着重研究了脆性温度区间内位移场的分布情况，得到了凝固裂 纹驱动力曲线，并与实验测量的凝固裂纹阻力曲线进行了 对比，从而对SUS310不锈钢的凝固裂纹敏感性进行了预测，为进一步深入的研究奠定了基础。 热源模型 三维的热源模型计算才能真实体现焊接的热过程，本文利用Goldak研究并提出的双椭球形热源 模型进行数值模拟[10]o前部分是个1／4椭球，后半部分是另1／4椭球。办、工与双椭圆高斯热源 模型中含义相同，则在前半部分椭球内热源分布为： e-3Z／c} (1) 拈烘加掣犟e-3X2／ae-3Y2／bZ abcl#xl# 同样，对后半部分椭球内热源分布为： 国家自然科学基金项目资助 245 (2) q(x，y，罗)：掣彩．3；V矿箩q，肛e．3：Vd 口bc，霞心霞 在公式(1)、(2)中参数8，b，c．--／有不同的值，它们是相互独立的。 有隈元模型的建立 本文在计算中采用的模型尺寸为loo锄× 100衄×t0眦，材料为SUS310不锈钢。由于对 称性取其一半进行分析，网格划分如图2所示。 网格划分时在焊缝金属及其附近区域进行了细 分，而在远离焊缝的地方网格划得较粗，这样不 仅可以减少模型的计算量，又保证了焊缝附近计 算的精度。 SUS310型单相奥氏体不锈钢材料的力学性 能参数采用参考文献【9】中数据。模型计算过程中 A，电压13．5V， 采用的焊接参数为：电流200 焊接速度2mm／$，焊缝未焊透，采用自由边界条件， 图2有限元网格划分 并在机械模型中充分考虑了凝固过程—凝固收缩，熔池变形的影响。 模拟结果及分析 本文在求解过程中，采用弱耦合解法，在每个时间增量步内，首先进行温度场分析，完成温度 场分析后再根据温度场的结果计算应交场，该方法只进行单向耦合，即只考虑温度场产生的应力应 变场，而忽略应力应变场对温度场的影响，应力应交场对温度场的影响相对于温度场产生的应力应 变场