焊接接头氢扩散数值模拟Ⅰ.pdf

第三届计算机在焊接中的应用技术交流会论文集 数值模拟与仿真 焊接接头氢扩散数值模拟(1) 哈尔滨焊接研究所 张显辉谭长瑛陈佩寅 摘要本文考虑了焊接接头中焊缝、热影响区及母材氢扩敞系数和溶解度的不同，采用ABAQUS有限元分析软件 对氢在辉接接头中扩散过程进行了数值模拟，并在此基础上，计算了预热及捍后后热对焊接接头氢扩散的影响。 关键词：数值模拟、氧、扩散、焊接接头 0前言 焊接氢致裂纹是低合金钢焊接时最容易产生，而且危害最为严重的焊接工艺缺陷，它常常是焊 接结构失效破坏的主要原因。因此，评定焊接氢致裂纹敏感性，预防氢致裂纹的产生一直是低合金 高强度钢焊接性和焊接工艺研究的最重要内容之一。大量研究工作表明，焊接区扩散氢含量、结构 的拘束应力水平及淬硬组织的存在是氢致裂纹产生的三个主要因素。当下述四种条件同时存在时， 焊接氢致裂纹将可能发生…：(1)扩散氢含量达到临界浓度；(2)足够大的应力强度；(3)对氢敏感的组 织结构；(4)温度低于200℃左右。由于氢原子体积小、活性强，即使在较低的温度下氢在金属中也 有较强的扩散能力。因此焊接接头微区中的瞬态氢含量是难以测定的。尽管圜内外学者对焊接接头 微区氢测定技术进行过多种尝试，但是迄今为止尚无成熟涮试技术可以用于焊接接头微区中的瞬态 氢含量测定。由于缺乏焊接接头徽区中瞬态氢含量的数据。目前尚不能提出氢致开裂的准确判据， 也难于深入地认识氢致裂纹产生的机制。如何确定氢致裂纹产生的临界氢浓度，如何准确地预测氢 致裂纹产生的时间和位置仍然是焊接界的一个难题。近年来，随着计算机技术的发展，采用数值分 析的手段对焊接过程中氢的扩散和聚集行为进行分析受到各国焊接学者的普遍重视。先后发表了一 些文章”1，其中大部分对氢的计算分析是基于Fick第二定律。即将焊接接头假定为一个均匀介质， 氢扩散的驱动力是氢的浓度梯度。但是，实际焊接接头中母材、热影响区和焊缝金属的微观组织存 在着显著的差异，同时，焊接接头还存在着焊接残余应力和应变，实际焊接接头不是一个均匀的介 质，Fick第二定律不能解决焊接接头微观氢扩散及聚集问题。本文采用ABAQUS有限元分析软件， 对氢在非均质焊接接头中的扩散过程进行了初步的计算模拟。 1氢扩散基本理论 扩散相在扩散中应满足质量守恒定律，即 … J小 J 式中V表示体积，S表示体积v的表面．13表示垂直表面S的矢量，J为扩散相的浓度流，n．J表示 离开S表面的浓度流。将式(1)中引入活度变量6由，再经过推导变换可得到 物质扩散的主导方程： 148 第三届计算机在焊接中的应用技术交流会论文集 数值模拟与仿真 在不均匀介质中，扩散的驱动力是扩散相的化学势梯度。氢在不均质焊接接头中的扩敞也遵循同样 的规律，即氢的扩散受氢的化学势梯度驱动： J：一一—旦!一．丝一一一一一一一一一一一一一一一一一一一一一一一一一一一一一(3) ax Rf曰一Ozl o+R(e—e 式(3)中，p表示化学势，p：u z)in巾+pVs 其中u。为一常量，R为气体常数，v*为氢摩尔分压。 由于浓度等于活度-q溶解度的乘积．将C--由S及p带入式(3)，经推导便有： 设 b=≯ln≯ 如=≯i若兰西 式(4)可变为： ．，…sD 扩散的“soret效应”系数，o表示温度，e z表示绝对温度零度，kp(C，0，f)表示由应力梯度 驱动扩散的应力系数，p—trace(0)／3为等效应力。f表示其他预定场变量。 将变量c=巾S及式(5)代入式(2)中便可以获得氢扩散本构方程； 式(6)中q一11．J是通过S面进入到物质中的扩散物浓度流。 2．2锄nNiMo钢焊接接头各区域中氢的扩散系数及溶解度测定 焊接接头组织的非均匀性使氢在焊接接头各区域中有不周的扩散性质。为了模拟计算氨在实际 焊接接头的扩散过程，采用气相渗透法对氢在焊接接头各区域中的扩散系数及溶解度进行测定。