管道环焊缝接头焊接残余应力的数值模拟论文.pdfVIP

第三届计算机在焊接中韵应用技术交流会论文集 数值模拟与仿真 管道环焊缝接头焊接残余应力的数值模拟 天津大学材料科学与工程学院董俊慧霍立兴张玉凤杨新歧荆洪阳 摘要本文利用ADTHA非线性分析有限元程序，对管道环焊缝接头焊接残余应力进行二维和三维有限元数值模 拟。在热弹塑性分析中考虑了材料熟物理和力学性能依赖于温度变化。通过焊缝单元的。死活”选择，模拟了焊接 熟源的移动过程。结果表明：在管道接头内表面焊缝中心及近缝区轴向和环向残余应力均为拉应力；|；l【离开焊缝距 离的增加，逐渐过渡为压应力。在管道接头外表面捍缝中心处的轴向残余应力为压应力，而环向残余应力为拉应力。 计算预测值与实测值基本一致．表明有限元法能够经济而有效地模拟预测管道环焊缝接头的焊接残余应力。 关键词：管道、环焊缝接头、残余应力、数值模拟 0 前言 管道结构广泛应用于石油、化工和市政建设等重要部门，在其焊接制造和安装过程中，不可避 免地会产生焊接残余应力。它的存在是造成焊接接头或结构脆性断裂、疲劳断裂和应力腐蚀开裂等 失效形式的重要因素。所以，焊接残余应力一直为人们所关注。 焊接臻余应力的大小和分布与焊接热源、接头形式、材料性能等多种因素有关。通常在实践中 采用理论分析和试验相结合的方法研究焊接残余应力的分布规律，以便达到预测、控制和调整它们 的目的。 本文在引入热弹塑性分析数学模型基础上，对管道环焊缝接头进行二维和三维有限元分析，结 果表明：计算值与实测值基本一致，证实了有限元法能够有效地模拟、预测焊接残余应力。 I 数学模型 对于焊接接头进行有限元分析有如下特殊性。(1)无外载荷力的作用，内应力保持平衡：(2) 焊接温度场是载荷条件：(3)材科热物理和力学性能随温度而交化。对于低碳铜管道，常用的弹塑 性本构关系为：材料服从Mises屈服准则，塑性区符合流变法则，并假设各向同性硬化，应力应变 在微小的时问增量内呈线性变化。由此得到弹塑性应力和应变偏增量之间有下列关系式‘“2】 2a。‘(1+去)口．(1+去)3(1十鲁)(¨ 峨瑙‰一蔫+焉一篝，(㈩3-4) 拈*13徘HⅦ删壶)-3拈7】 ：3，-14’ 式中：品—应力偏量；8r一应变偏量；吼广_应力张量；暂一应变张量；op一有效应力；q一屈服 应力；耳—温度；如7二轴向热膨胀增量；G一剪切模量；ⅣL塑性模量；忌一体积模量【K；肋(1-2v)】； 其中B一弹性模量；p一泊松比。材料常数G和K仅随温度而变化。 利用式(1)和(2)可求得应力的变化d嘞为 峨-2G㈣口+南删㈦一‘：3，G2(1Su+d尹,一 c3㈤-4) 州剪帕{。焉8T一．．}Sv-d“K[d(扣3∽ 第三届计算机在焊接中的应用技术交流会论文集 数值模拟与仿真 I?der。|酗￡dV=0 (3-4) (4) 在弹性条件下，括号{)中部分为零。由于无外载荷，虚功原理可写成 式中：耐旬是应变增量幽的变量，应力增量硒由(3)式给出。 公式(3)也可改写成矩阵表示的形式吲： …、 LjqJ 7】 [da】=【D11ds】+ClSl一M【d5【5) 在平面应力条件下，式(5)可写成广义记号法形式。对于未知的结点位移和位移导数方程系 统可由广义记号法形式得到81 【同【胡：【用 ‘翟 式中膳7是单元集合的刚度矩阵．它是由单元刚度矩阵，叼。之和组成。皿即是结点位移矢量。，刀 是由单元的等效结点载荷矢量，R，之和组成，矩阵‘W为应变矩阵。 每个单元的刚度矩阵为