基于SYSWELD的搭接接头温度场的数值模拟
基于SYSWELD的运行管道在役焊接热循环数值模拟
1.2 热源模型及焊接工艺参数 由于双椭球 模 型 更 为 准 确［8］，焊 接 热 源 采 用 双
椭球模型。采用 焊 条 电 弧 焊 方 法，根 据 不 同 电 流 下 焊接熔池的深度和宽度初步确定双椭球热源模型的 各参 数，然 后 采 用 SYSWELD 的 热 源 拟 合 工 具 进 行 多 次 校 核 ，直 至 模 拟 出 的 熔 池 形 状 和 实 际 接 头 相 符 。 选用的四种热输入及其相关参数如表 1。
工艺参数 编号
A B C D
表 1 焊接工艺参数 table 1 Welding parameters
焊接电流 I/A
电弧电压 焊接速度 热输入
U/V
1 （/ mm·S - 1）E （/ kJ·cm - 1）
课程设计——基于sysweld软件的T型接头焊接仿真模拟
——焊接基于sysweld软件的T型接头焊接仿真模拟姓名：000班级：材料000班学号：00000000指导老师：000日期：2011年09月SYSWELD——法国ESI公司的焊接仿真分析软件，经20多年发展，已成为热处理、焊接和焊接装配过程模拟的领先模拟软件，能够全面考虑材料特性、设计和过程的各种情况。
随着科学技术的发展，机械制造行业也随之不断的革新和进步。
人们对铸件的质量要求也越来越高，而SYSWELD为其提供了一个良好的工具，对提高铸件的质量有未雨绸缪的作用。
SYSWELD热过程模拟软件对铸件的制造起着非常关键的作用，为解决铸件缺陷问题提供了一个平台。
利用SYSWELD软件对焊缝进行计算机仿真模拟来提高焊缝的质量，本文主要对焊接的热过程模拟来分析T形接头焊焊接热过程，主要通过T形建模、热源校核、焊接向导、求解计算及结果后处理的操作步骤对焊接热过程进行数值模拟。
与测试并修正的传统方法相比，SYSWELD使得成本降低、周期缩短。
另外还能够显著减少物理样机，产生高的投资回报率。
界面友好，轻松易学。
SYSWELD 是用于引导工程师发现关于变形、残余应力和塑性应变的影响因素，然后优化过程参数的专业模拟软件。
2011-09-091、T型接头模型的建立1.1创建Points (1)1.2由Points生成Lines (1)1.3由Lines生成Edges (2)1.4由Edges生成Domains (2)1.5离散化操作 (3)1.6划分2D网格 (5)1.7生成Volumes (6)1.8离散Volumes (8)1.9生成体网格 (10)1.10划分换热面 (11)1.11划分1D网格 (12)1.12合并节点 (13)1.13保存模型 (14)1.14组的定义操作 (15)1.15保存 (17)1.16小结 (17)2、焊接热源校核2.1网格的建立 (18)2.2材料的导入及定义 (20)2.3热源过程参数的定义 (20)2.4求解 (21)2.5热源显示 (21)2.6修改参数 (22)2.7热源校核 (22)2.8检查显示结果 (23)2.9保存函数 (24)2.10热源查看 (24)2.11保存热源 (25)2.12小结 (25)3、焊接模拟向导设置3.1材料的导入 (26)3.2热源的导入 (26)3.3材料的定义 (27)3.4焊接过程的定义 (27)3.5热交换的定义 (28)3.6约束条件的定义 (28)3.7焊接过程求解定义 (28)3.8冷却过程求解定义 (29)3.9检查 (29)3.10小结 (31)4、后处理与结果显示分析4.1计算求解 (32)4 .2导入后处理文件 (32)4.3结果显示与分析 (33)4.4小结 (36)1、T型接头模型的建立1.1创建Points根据所设计T型接头模型的规格，选定原点，然后分别计算出各节点的坐标，按照Geom./Mesh.→geometry→point步骤，建立以下13个点：P1（-25,0，-10）、P2(7,0，-10)、P3(10，0，-10)、P4(13,0，-10)、P5(35,0，-10)、P6(35,0,0)、P7(10,0,0)、P8(10,0,30)、P9(0,0,30)、P10(0,0,3)、P11(-1.5,0,1.5)、P12(-3,0,0)、P13(-25,0,0)如下图所示：1.2由Points生成Lines按照Geom./Mesh.→geometry→1Dentities步骤，按照一定的方向性将各点连接成如下图所示的Lines：1.3由Lines生成Edges按照Geom./Mesh.→geometry→EDGE步骤，点击选择各边，依次生成如下图所示各Edges：1.4由Edges生成Domains按照Geom./Mesh.→geometry→Domains步骤，依次生成如下六个Domains：1.5离散化操作离散化操作是针对由Points所生成的Lines而言，由于除了有这些点生成的线以外，软件本身也会自动产生一些辅助的线条，为了方便清晰地对所生成的主要线条进行选取及其他操作，可以通过“隐藏→显示”处理，只显示如下图所示的十八条线：通过以下操作为后面的离散操作做好准备：→通过Meshing→Definition→Discretisation启动离散化操作界面，将L2、L4、L8、L10四条线均匀离散成3段，将其他十四条线非均匀离散，离散单元数为5，系数为3.5。
基于SYSWELD的焊接模拟仿真
拟分析 ， 为实际焊接生产 中控制焊接 变形提供 了 理 论依 据 ．
参考 文献 ：
［ １ ］ 郭艳 ， 凌泽 民， 齐喜岑 ， 等． 基于 Ｓ Ｙ Ｓ ＷＥ Ｌ Ｄ铝 合金管点 焊温度场数值 模拟 ［ Ｊ ］ ． 热加工工艺 ， ２ ０ １ ０ ， １ ９ ： ４ ５ ４９ ．
Ｅ－ ｍａ ｉ ｌ ： ｓ ｕ ｈ ２ ２ １ ＠１ ６ ３ ． １ ２ ０ １ １ １ ．
８ ２
大 连 交 通 大 学 学 报
第３ ４卷 Leabharlann （ ３ ） 应力 分 析
程 中大 体上 不发 生变化 ；
（ ２ ） 由 于焊 缝 的纵 向收 缩 ， 底 板 在 各 道 焊 缝
源模 型． 本 文所 选 构件 采 用 自动 钨 极 氩 弧焊 工 艺
１ 构 件 模 型 分 析
（ １ ） 三维 网格 划分 钛合金 目前应用 广泛 ｒ ｚ ｊ ， 本 文所选构件 材料为 钛合金 Ｔ Ａ １ ５ ， 底 板 为 ２ ． ０ ｍ ｍ， 立筋 ６为 ２ ． ５ ａ ｒ ｍ， 共 ７道焊 缝 ， 采用 自动 钨极 氩 弧 焊工 艺 进行 焊 接 ． 首先对该 构件进 行 了网格 划分 ， 网格 划分 直接影 响
基于 Ｓ Ｙ Ｓ ＷＥ Ｌ Ｄ 的 焊 接 模 拟 仿 真
苏杭 ， 常荣辉 ， 倪 家强
（ 中国航 空工 业集 团 沈 阳飞机工 业（ 集 团） 有 限公 司，辽宁 沈 阳 １ １ ０ ８ ５ ０ ） 米
摘
要： 针对典型焊接构件 ， 采用 Ｓ Ｙ Ｓ ＷＥ Ｌ Ｄ焊接模 拟软件 进行模 拟 ， 建立 了三维模 型并 进行 了 网格 划
第３ ４卷 第 ２期 ２ ０ １ ３年 ４月
基于SYSWELD的T型接头GMAW焊接热过程模拟及其应用
焊接热过程贯穿整个焊接过程，一切焊接物理化学过 程都是在热过程中发生和发展的。焊接温度场决定了焊接 应力场和应变场，还与冶金、结晶、相变过程以及焊缝成 型有着不可分割的联系。因此，焊接热过程是影响焊接质 量和生产效率的主要因素之一，焊接热过程的准确计算和 测量是进行焊接冶金分析、焊接应力应变分析和对焊接过 程控制的前提。
焊接速度 /（mm/s）
4.5 3.6 3.0 2.6 2.25 2.0
焊接线能量 /（kJ/cm）
9.90 12.37 14.84 17.13 19.79 22.26
根据上述方法确定的焊接线能量对热源参数影响规律
如图 4 所示。图 4（a）、（b）分别为焊接线能量对热源深 度 b 和半宽 a 的影响，由焊接速度和半宽 a 就可以根据图 4 （c）、4（d）确定出热源长度方向的尺寸 c1 和 c2。
设计 与 研 究
1
基于 SYSWELD 的 T 型接头 GMAW 焊接热过程模拟及其应用
卢庆亮 1 曹永华 1 杨 云 1 栾守成 1 左增民 2 华 鹏 3 孙俊生 3
（1. 济南重工股份有限公司，济南 250109；2. 菏泽广泰耐磨制品股份有限公司，菏泽 274600； 3. 山东大学 材料学院，济南 250061）
三维双椭球热源模型把熔池设为两个半椭球的组合体， 其尺寸和形状由参数 c1、c2、a 和 b 来限定，如图 2 所示， 而这些参数根据实际焊缝横截面和焊缝表面波纹的实测数 据确定。
而使得熔池的尺寸变小。因此，计算时“设定”熔池的尺
寸应该比实测大。计算和试验测试结果表明，由实际焊缝
横截面和焊缝表面波纹测得的实际熔池尺寸增加 5% ～ 10%
计算试件的尺寸和 T 型接头的坐标系如图 1 所示，材 料为 JB800 贝氏体钢，JB800 钢的化学成分和力学性能如 表 1、表 2 所示。GMAW 焊接电弧以恒定速度 v 从固定坐 标系 O-yxz 的原点 O 开始，沿 y 方向运动；移动坐标系的 原点 O′位于电极的中心线上，并随同电弧同步移动，所 以 ζ=y-vt。
基于Sysweld的T形管焊接件温度及应力应变场数值模拟分析
z z z PsgiolePfrp
z z z PsgiolePfrp
z z z PsgiolePfrp
９４
焊 接 学 报
第 ３７卷
z z z PsgiolePfrp
［上接 第 ８Ｏ页 ］
［６］ 黎超文 ，王 勇 ，李立英 ，等．Ｔ形接头 的焊接温 度场 三维 动态有限元模拟［Ｊ］．焊接学报 ，２０１１，３２（８）：３３—３６．
一 ． 生 产加工 中大多 由焊工 熟 练 的焊 接技 能 来 保证 仿 真 ，模 拟仿 真流程 如 图 １所 示 ．
焊 接质量 ，而焊 接 时的应力 影 响难 以掌 控 ，需 要 大量
离空 要卜 矗驽 的试验¨Ｊ．当前试验测量应力的方法主要有盲孔法
和 ｘ射线衍 射法 ，这 些方 法 需 要 专 门配 套 的 试验 设 备 ，试验周 期 长 ，并 且 很 难 得 到 分 布 复 杂 多 样 的应
图 １ 基 于 Ｓｙｓｗｅｌｄ数 值 仿 真 分 析 流 程
Ｆｉｇ．１ Ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ ｆｌｏｗ ｂａｓｅｄ ｏｎ Ｓｙｓｗｅｌｄ
成本 ．因此采 用计算 机模 拟 方法 模 拟 焊接 应 力具 有
较大 的经济 效益 ． 基于 Ｓｙｓｗｅｌｄ有 限元 分 析 软 件对 焊 接 过 程 中温
z z z PsgiolePfrp
第３ ７卷 第４期
焊 接 学 报
Ｖｏ１．３７ Ｎ。．４
２ ０ １ ６ 年 ４ 月 ＴＲＡＮＳＡＣＴＩＯＮＳ ＯＦ ＴＨＥ ＣＨＩＮＡ ＷＥＬＤＩＮＧ ＩＮＳＴＩＴＵＴＩＯＮ Ａｐｒｉｌ ２ ０ １ ６
基 于 Ｓｙｓｗｅｌｄ的 Ｔ形 管焊 接 件 温 度 及 应 力 应 变 场 数 值 模 拟 分 析
T型接头激光焊接的温度场和应力场的数值模拟
Ｖｏｌ．２４Ｎｏ．４安徽工业大学学报第２４卷第４期Ｏｃｔｏｂｅｒ２００７Ｊ．ｏｆＡｎｈｕｉＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ２００７年１０月文章编号：１６７１－７８７２（２００７）０４－０３８４－０５Ｔ型接头激光焊接的温度场和应力场的数值模拟丁林，周永涛，李明喜（安徽工业大学材料科学与工程学院，安徽马鞍山２４３００２）摘要：基于ＳＹＳＷＥＬＤ的焊接分析功能，采用有限元方法研究激光动态焊接过程中温度场、应力场、应变场的变化情况，应用ＳＹＳＷＥＬＤ软件的校正工具对三维高斯热源进行校核。
考虑各相的热物理性能参数与温度的非线性关系，建立焊接过程的数学模型和物理模型，以不锈钢Ｘ５ＣｒＮｉ１８１０为例，对Ｔ型接头进行三维动态模拟。
结果表明：随焊接速度的减小，热循环在高温时刻停留时间增加，冷却速度减慢；随着远离起始端距离的增加拉应力值逐渐减小转变为压应力，最后趋向零。
关键词：温度场；应力场；应变场；有限元法中图分类号：ＴＧ４０２文献标识码：ＡＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎｏｆＴｅｍｐｅｒａｔｕｒｅＦｉｅｌｄａｎｄＳｔｒｅｓｓＦｉｅｌｄｏｆＴ－ｊｏｉｎｔＤＩＮＧＬｉｎ，ＺＨＯＵＹｏｎｇ－ｔａｏ，ＬＩＭｉｎｇ－ｘｉ（ＳｃｈｏｏｌｏｆＭａｔｅｒｉａｌｓＳｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇ，ＡｎｈｕｉＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｍａ＇ａｎｓｈａｎ２４３００２，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｔｈｅｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｏｆｔｅｍｐｅｒａｔｕｒｅｆｉｅｌｄ，ｓｔｒｅｓｓｆｉｅｌｄａｎｄｓｔｒａｉｎｆｉｅｌｄｄｕｒｉｎｇｌａｓｅｒｗｅｌｄｉｎｇｂａｓｅｄｏｎｔｈｅｗｅｌｄｉｎｇａｎａｌｙｓｉｓｆｕｎｃｔｉｏｎｏｆＳＹＳＷＥＬＤｗｅｒｅｉｎｔｒｏｄｕｃｅｄ．Ｍｏｒｅｏｖｅｒ，ｔｈｅｈｅａｔｓｏｕｒｃｅｉｓｍｏｄｉｆｉｅｄｗｉｔｈｔｈｅｔｏｏｌｓｕｐｐｌｉｅｄｂｙＳＹＳＷＥＬＤｓｏｆｔｗａｒｅ．Ｔｈｅｔｈｅｒｍｏ－ｐｈｙｓｉｃａｌｐｒｏｐｅｒｔｉｅｓｅｓｔａｂｌｉｓｈｅｄａｓｔｈｅｆｕｎｃｔｉｏｎｓｏｆｔｅｍｐｅｒａｔｕｒｅｗｅｒｅｔａｋｅｎｉｎｔｏｃｏｎｓｉｄｅｒａｔｉｏｎ．Ｔ－ｊｏｉｎｔｗｅｌｄｉｎｇｏｆｓｔａｉｎｌｅｓｓｓｔｅｅｌＸ５ＣｒＮｉ１８１０ｗａｓｓｉｍｕｌａｔｅｄ３Ｄｄｙｎａｍｉｃａｌｌｙ．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｃｅａｓｅｉｎｇｔｉｍｅｏｆｈｅａｔ－ｃｙｃｌｅｉｓｉｎｃｒｅａｓｅｄａｎｄｃｏｏｌｉｎｇｖｅｌｏｃｉｔｙｉｓａｌｓｏｄｅｃｒｅａｓｅｄｗｉｔｈｔｈｅｄｅｃｒｅａｓｅｏｆｗｅｌｄｉｎｇｖｅｌｏｃｉｔｙ．Ｗｉｔｈｉｎｃｒｅａｓｅｏｆｄｉｓｔａｎｃｅ，ｔｅｎｓｉｌｅｓｔｒｅｓｓｉｓｔｒａｎｓｌａｔｅｄｉｎｔｏｃｏｍｐｒｅｓｓｉｖｅｓｔｒｅｓｓａｎｄｔｅｎｄｓｔｏｚｅｒｏ．Ｋｅｙｗｏｒｄｓ：ｔｅｍｐｅｒａｔｕｒｅｆｉｅｌｄ；ｓｔｒｅｓｓｆｉｅｌｄ；ｓｔｒａｉｎｆｉｅｌｄ；ｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄ近年来，随着计算机技术和仿真算法的发展、完善，焊接模拟技术变得越来越重要，它不仅能够有效地提高产品的经济效益，还可以节省大量的时间。
一基于SYSWELD的在役焊接径向变形数值模拟
收稿日期:2013-12-11基金项目:国家自然科学基金项目(51074174);山东省自然科学基金项目(ZR2010EM041)作者简介:王勇(1964-),男,教授,博士生导师,主要从事新材料焊接技术㊁金属失效与表面改性及压力容器与管道安全工程的研究工作㊂E-mail:wangyong@㊂ 文章编号:1673⁃5005(2014)03⁃0117⁃06 doi:10.3969/j.issn.1673⁃5005.2014.03.019基于SYSWELD 的在役焊接径向变形数值模拟王 勇,郭广飞,韩 涛,李 浩,孙启平(中国石油大学机电工程学院,山东青岛266580)摘要:采用焊接模拟软件SYSWELD 研究壁厚㊁管径及熔池尺寸等因素对天然气管道在役焊接径向变形的影响㊂结果表明:发生烧穿的临界瞬时最大变形有时间效应;发生临界变形的时间随着壁厚的增大而增加,径向变形量减小;随着管径增加,径向变形量在同壁厚下逐渐增大;壁厚为4.5mm 时,小管径管道易达到临界变形量;壁厚为6mm 时,管径的增大降低了临界变形发生的可能性;当壁厚为7.5mm,管径由254mm 增大到508mm 时,发生临界变形的可能性减小,而管径在508~1016mm 时,管径的增大增加了临界变形发生的可能性;内壁点的径向变形随着热源的靠近而增大,随着壁厚的增大而减小,但当壁厚增大到6mm 后,壁厚的增大对其不再有显著影响;熔池尺寸影响焊接修复同时刻时的径向变形量,达到临界变形量的时间与熔池尺寸成反比,表现出明显的熔池尺寸效应㊂关键词:在役焊接;天然气管线;径向变形;数值模拟中图分类号:TG 402 文献标志码:ANumerical simulation of radial deformation during in⁃servicewelding based on SYSWELDWANG Yong,GUO Guang⁃fei,HAN Tao,LI Hao,SUN Qi⁃ping(College of Mechanical and Electronic Engineering in China University of Petroleum ,Qingdao 266580,China )Abstract :The welding software of SYSWELD was used to simulate the effects of wall thickness,pipe diameter and the pool size on the radial direction deformation during the in⁃service welding of pipeline.The results show that:the critical instanta⁃neous maximum deformation,which will result in burn through pipeline,has time effect;the time reaching critical deforma⁃tion increases with the increase of wall thickness,at some time,the deformation amount reduced.With the increase of pipe diameter,the radial direction deformation will increase at the same wall thickness;when the wall thickness is 4.5mm,the smaller the pipe diameter is,the more easily the critical deformation is reached;when the wall thickness is 6mm,the in⁃crease of pipe diameter will reduce the generation rate of critical deformation;when the wall thickness is 7.5mm,the gener⁃ation rate of critical deformation decreases with the increase of diameter from 254mm to 508mm,and at a pipe diameter range from 508mm to 1016mm,the generation rate increases with the increase of pipe diameter.The radial direction de⁃formation of inner wall increases when the heat source is closer,and it reduces with the increase of wall thickness.And when the wall thickness is greater than 6mm,the wall thickness no longer has significant effect on the radial direction deformation.The size of pool affects the radial direction deformation of weld repair,and the time reaching critical deformation reduces with the increase of pool size,which shows significant pool size effects.Key words :in⁃service welding;gas pipeline;radial deformation;numerical simulation 天然气管道在使用过程中由于腐蚀等因素的影响而局部减薄,需要进行修复,在役焊接具有不停输带压修复的特点,可保证管道输送的连续性,避免管道停输进行修复所带来的损失[1⁃2]㊂焊接电弧热使得焊接处的承载能力下降到一定程度时[3],发生烧穿失稳,带来严重的安全问题[4],所以防止烧穿是2014年 第38卷 中国石油大学学报(自然科学版) Vol.38 No.3 第3期 Journal of China University of Petroleum Jun.2014在役焊接修复所要考虑的首要问题㊂在役焊接时,管道受到内部压力及焊接应力的共同作用,力的作用使得被焊接管道发生一定量的径向变形,焊接热输入则使被焊接管道的强度逐步降低[5],促进了径向变形量的增大,当变形量超过被焊接板厚的0.1倍时[6],极易发生烧穿㊂被焊接修复管道的壁厚及管径对径向变形有显著影响,进而影响烧穿㊂笔者采用焊接模拟软件SYSWELD 进行模拟计算,并结合试验研究壁厚及管径对在役焊接径向变形的影响,揭示这些因素对烧穿的影响㊂1 在役焊接数值模型的建立1.1 几何模型在SYSWELD 自带的Visual⁃Mesh 软件下进行建模,采用SYSWELD 进行焊接模拟计算,焊接修复接头㊁焊道分布及焊接顺序如图1所示㊂图1中1㊁2㊁3㊁4表示在役焊接修复时的4道焊缝,由于烧穿一般发生于第1道焊缝,所以针对焊接接头的第1道焊缝焊接时管道的径向变形进行研究,且由于第1道焊缝对于被修管线的焊接属于表面堆焊,采用图2所示的有限元模型进行数值计算㊂图1 焊接接头示意图Fig.1 Sketch map of welded joint由于管道是轴对称的,采用二分之一模型以及四分之一两种模型,其中管径U r 为254mm 及508mm 的管道采用二分之一模型(图2(a)),管径1016mm 的管道采用四分之一模型(图2(c)),焊缝区的网格划分较密,其余部分划分稀疏,减少了计算量并得到较好的计算精度㊂模型采用4.5㊁6.0和7.5mm 三种厚度,管道半宽均为200mm㊂图2 有限元模型Fig.2 Finite element model1.2 焊接热源模型及参数双椭球热源模型能准确地描述在役焊接过程的热源特性[7],所以采用双椭球模型进行焊接模拟,如图3所示㊂以O 点为原点建立三维坐标系,其中图3 双椭球热源模型Fig.3 Double ellipsoid heat source modela f 和a r 分别为前半椭球和后半椭球的长度参数,b h ㊁c h 分别为椭球宽度和深度参数;采用二氧化碳气体保护焊(GMAW)进行在役焊接修复试验,根据实际焊接过程中产生熔池的深度及宽度确定双椭球热源模型的各参数,再运用SYSWELD 的热源校核工具进行多次校核,直至计算模拟的熔池形状与实际接头符合为止,得到热源模型;采用110A 的电流㊁24V 的电压㊁焊接速度为3mm /s 及热输入为8.8kJ /cm 的焊接参数进行在役焊接试验,管道内部压力为6MPa㊂1.3 换热边界条件及约束管道外表面和空气的换热方式主要考虑热辐射及空气的自然对流换热[8⁃10],总换热系数为α=4.536×10-8(546.3+T 0+T 1)[(273.15+T 0)2+(273.15+T 1)2]+25.(1)式中,T 0㊁T 1分别为环境温度(25℃)和焊接接头与空气接触表面的温度,℃㊂在役焊接接头背面,即焊接管道的内壁和天然气的换热为强迫对流换热形式,其换热系数为α1=0.027λd i Re 0.8Pr 1/3μμæèçöø÷w 0.14.㊃811㊃中国石油大学学报(自然科学版) 2014年6月其中μw =μ0(273.15+T 2)/273.15.(2)式中,λ为气体的导热系数,W /(m㊃K);Re 为气体的雷诺数;Pr 为气体的普朗特数;μ为气体的动力黏度,Pa㊃s;d i 为修复管道的内径,m;μw 为气体在壁温时的动力黏度,Pa㊃s;T 2为管道内壁的温度,℃㊂因为天然气中甲烷占绝大多数,计算中天然气的换热系数采用甲烷的代替,参照文献[9]的数据进行确定㊂通常,天然气的换热系数随着压力的不同也发生着相应的变化,需要针对不同压力的情况,对换热系数α1进行相应的调整,本文中采用6MPa 的恒压,换热系数采用6MPa 下的定值㊂因为要施加内部压力载荷,须对计算模型进行约束㊂由于分别采用了焊道的1/2及1/4进行计算,图2(a)和图2(c)中的模型所受的约束是不同的,令图2(a)中的模型受到沿对称面yz 及面xy 的面约束以及沿y 轴的线约束,而令图2(c)中的模型受到沿对称面xy ㊁xz ㊁zy 的面约束,沿y ㊁z 轴的线约束㊂X70管线钢的导热系数及质量热容比采用文献[11]的公式计算得到,力学性能采用文献[12]提供的数据㊂2 计算结果分析采用254㊁508及1016mm 三种不同管径以及同管径4.5㊁6.0及7.5mm 三种壁厚的模型,用SY⁃SWELD 进行模拟计算,得到各壁厚及管径下150s 各时刻变形场,然后运用软件自带的后处理模块Post processing 得到各时刻最大瞬时径向变形量的计算数据,再进行数据处理,得到相应的计算结果㊂2.1 壁 厚进行模拟计算,得到同管径时壁厚对径向变形量的影响见图4㊂由图4可知,无论何种壁厚以及何种管径,发生临界径向变形量(即壁厚的0.1倍)都需要一定时间,表现出一定的时间效应㊂图4 壁厚对径向变形量的影响Fig.4 Effect of thickness on radial direction deformation 图4(a)中,管径为254mm㊁壁厚4.5mm 的管道在焊接进行到38s 时,变形超过0.45mm;而在焊接修复进行到47s 时壁厚6.0mm 的管道发生0.6mm 的变形量;壁厚7.5mm 的管道发生0.75mm 的径向变形量则是在70s 时㊂图4(b)中壁厚为4.5㊁6.0㊁7.5mm 时,临界变形量分别发生在42㊁58㊁85s 以后㊂同种壁厚时,管径508mm 的管道比管径254mm 的管道发生临界变形量的时间要延后㊂图4(c)中壁厚为4.5㊁6.0㊁7.5mm 时,管道临界烧穿变形发生在42㊁62㊁81s 以后㊂发生时间和管径为508mm 的管道时间相差很小,所以可使用管径508mm 的模型代替管径1016mm 的模型考察壁厚对径向变形的影响㊂由此可见,相同管径下,随着壁厚增大,需要达到相应壁厚的临界径向变形量的时间越来越长,在限定时间内,发生烧穿的可能性越小,可以安全施焊的时间逐渐增加㊂此外,在同管径不同壁厚时,壁厚越小,焊接同时刻的径向变形量越大,并且随着焊接时间的延长,差异性越来越大,呈现出明显的时间效应,但在具体情形中又呈现出不同㊂在图4(a)中,壁厚4.5mm 及6.0mm 的管道在焊接开始一直到40s 之间,径向变形量几乎相同,而壁厚7.5mm 的管道变形量明显较小㊂2.2 管 径使用相同壁厚不同管径的有限元模型分别计算,考察管径对径向变形的影响见图5㊂用U r 254㊁U r 508㊁U r 1016表示管径为254㊁508及1016mm 的管道径向变形㊂㊃911㊃第38卷 第3期 王 勇,等:基于SYSWELD 的在役焊接径向变形数值模拟图5 不同管径的径向变形量Fig.5 Radial direction deformation of different diameters 由图5可知:壁厚为4.5mm时,在焊接修复开始一直到50s期间,3种管径的径向变形量差别很小㊂从50s开始,管径254mm的管道在50s后变形有一个变化率的突变,但是在90~150s期间总体变形变化率是减小的;管径508mm以及管径为1016mm的管道在焊接修复起始一直到120s期间,总体上的变形量几乎相同,而在150s时Ur1016最大,U r508次之, U r254最小㊂壁厚为6mm时,在80s之前,U r254大于U r508及U r1016,而在90s后,U r254的变化将趋缓并有减小趋势,而U r508及U r1016依然呈现增大的趋势㊂在壁厚为7.5mm的管道模型计算结果中,3种管径在前80s的径向变形量差异已经很小,在80s 后,逐渐表现出变化率的不同,其中U r508及U r1016迅速增大㊂分析3种壁厚不同管径的模型计算结果可知,随着管径增大,径向变形量在相同壁厚的情形下不断增加;随着壁厚增大,在相同时刻的变形量逐渐减小,壁厚的增大会显著地降低径向变形,从而降低发生烧穿的可能性㊂管径的大小在前期并不会显著地改变变形量的增加趋势,并且管径较小的管道在焊接修复前期变形量较大㊂随着焊接修复的进行,在80s或者90s 之后会呈现不同的变化,在90s之后,变形量随着管径的增大而增加,呈现出一定的时间效应㊂按照文献[6]的评判标准,在壁厚为4.5mm的模型中,当焊接到40s时,径向变形量达到0.5mm,已经达到了烧穿的临界值;在此壁厚下,管径较小的管道发生烧穿的可能性更大㊂当模型壁厚为6.0mm,在50s时,U r254和U r508均超过0.6mm,而U r1016在65s左右时才达到临界值0.6mm㊂壁厚为7.5mm,在80s时U r1016先达到临界值0.75mm,管径分别为254和508mm的管道在90s 时达到0.75mm的变形量㊂表1 熔池尺寸Table1 Pool size熔池名称a/mm b/mm c/mmA′9.696.753.56B′7.274.691.78C′7.622.811.41综上分析可知,壁厚为4.5mm的管道焊接修复时,管径254mm的管道较508mm及1016mm的管道易发生烧穿,壁厚为6.0mm的管道修复过程中也表现出此种特性;而在7.5mm壁厚的管道修复时,管径为1016mm的管道先于管径为508mm 及254mm的管道达到临界值0.75mm,说明随着壁厚增大,管径增加增大了烧穿的可能性㊂2.3 内壁径向变形以管径为508mm的模型为例进行模拟计算,考察壁厚对管道内壁一点的径向变形影响,内壁点的选取如图2(a)所示,计算结果见图6㊂图6 内壁点径向变形示意图Fig.6 Sketch map of radial direction deformationof node with different thicknesses由图6可见,对于管径508mm的管道,其内壁一点径向变形相对于径向瞬时最大变形(图4(b))呈现如下特点:首先是变形量较之小得多,其在150 s内的最大变形量仅在0.75mm左右;其次,壁厚6㊃021㊃中国石油大学学报(自然科学版) 2014年6月mm 和壁厚7.5mm 管道的内壁点的径向变形不论是在趋势还是量上基本上一致㊂壁厚因素在焊接修复的前期,对内壁点径向变形的影响差别不是很大,但是随着焊接修复时间增加,热源逐渐靠近,壁厚4.5mm 的管道变形量的增大率迅速增大㊂当壁厚增加到6mm 时,壁厚的增加对于内壁点的径向变形量不再有显著影响㊂2.4 径向变形的尺寸效应以往的计算过程中,很少考虑到熔池尺寸的效应对径向变形量的影响,但在计算过程中发现,其对变形影响十分明显㊂采用相同的焊接参数进行焊接,得到3种不同形状的熔池(表1)的计算温度场,对熔池温度场进行后处理,选取同时刻3种熔池的截面,a =a f +a r ,b ㊁c 各含义如图7所示㊂图7 熔池截面示意图Fig.7 Section diagram of pool 3种熔池计算过程中受到的内部压力均为6MPa,焊接修复时间为90s,对3种熔池在焊接过程中的径向瞬时最大变形量进行计算,计算结果见图8㊂图8 不同熔池的径向变形Fig.8 Radial deformation of different pools从图8的计算结果可知,在焊接修复前30s 期间,熔池A′㊁B′产生的变形量基本相同,熔池C′所形成的焊道变形量较小,并且在随后的焊接修复过程中一直表现着较小的变形量,达到临界变形量的时间也明显小于熔池A′㊁B′所产生的焊道;随着焊接修复的不断进行,从30s 开始,熔池A′所产生的径向变形量开始明显大于较小的熔池所产生的径向变形量,表现出较为明显的熔池尺寸效应,并率先达到临界变形量0.45mm,而熔池尺寸较小的修复方案达到临界径向变形量所需要的时间较长,与熔池尺寸成反比㊂熔池尺寸效应的显现也需要一定的时间体现,其与径向变形量的时间效应是吻合的,这是因为焊接修复过程中,热的传导及被焊管道的强度随着热作用的下降程度都需要一个过程时间来显现㊂3 结 论(1)发生烧穿的临界瞬时最大变形有时间效应,发生临界变形的时间随着壁厚增大而增加,变形量减小㊂(2)管径增加,径向变形量在同壁厚下逐渐增大;壁厚为4.5mm 时,小管径管道易达到临界变形量;壁厚为6mm 时,管径的增大降低临界变形发生的可能性;当壁厚为7.5mm 时,管径由254mm 增大到508mm 时,发生临界变形的可能性减小,而在508~1016mm 区间时,管径的增大增加了临界变形发生的可能性㊂(3)内壁点的径向变形随着热源的靠近而增大,随着壁厚的增大而减小,但当壁厚增大到6mm 后,壁厚的增大对其不再有显著影响㊂(4)焊接修复过程中,熔池尺寸影响焊接修复同时刻时的径向变形量,达到临界变形量的时间与熔池尺寸成反比,表现出明显的熔池尺寸效应㊂参考文献:[1] 黎超文,王勇,韩彬,等.高压管线在役焊接烧穿研究进展[J].压力容器,2008,25(8):40⁃45.LI Chao⁃wen,WANG Yong,HAN Bin,et al.Advanceof burn⁃through research on in⁃service welding of high pressure gas pipelines [J].Pressure Vessel Technology,2008,25(8):40⁃45.[2] 陈玉华,王勇,韩彬.X70管线钢在役焊接局部脆化区的组织及精细结构[J].材料热处理学报,2007,28(1):77⁃80.㊃121㊃第38卷 第3期 王 勇,等:基于SYSWELD 的在役焊接径向变形数值模拟CHEN Yu⁃hua,WANG Yong,HAN Bin.Metallurgicalmicrostructure and fine structure in local brittle zone ofin⁃service welding of X70pipeline steel[J].Transactionsof Materails and Heat Treatment,2007,28(1):77⁃80.[3] 边育智,刘洋,王新华.X65㊁X70管线钢高温延塑性的对比分析[J].钢铁研究学报,2007,19(11):51⁃54.BIAN Yu⁃zhi,LIU Yang,WANG Xin⁃parisonand analysis on high temperature ductility of X65and X70pipeline steel[J].Journal of Iron and Steel Research,2007,19(11):51⁃54.[4] OTEGUI J L,CISILINO A,RIVAS A,et al.Influenceof multiple sleeve repairs on the structural integrity of gaspipelines[J].International Journal of Pressure Vesselsand Piping,2002,79(11):759⁃765.[5] 郭广飞,王勇,韩涛.压力管道在役焊接烧穿失稳机制[J].材料科学与工艺,2012,20(3):62⁃66.GUO Guang⁃fei,WANG Yong,HAN Tao.Burn throughfailure mechanism of in⁃service welding of pipeline[J].Materials Science&Technology,2012,20(3):62⁃66.[6] API1104.Welding of pipelines and related facilities,ap⁃pendix B:in⁃service welding[S].USA:American Petro⁃leum Institute,1999.[7] 李培麟,陆皓.双椭球热源参数的敏感性分析及预测[J].焊接学报,2011,22(1):85⁃89.LI Pei⁃lin,LU Hao.Sensitivity analysis and prediction ofdouble ellipsoid heat source parameters[J].Transactionsof the China Welding Institution,2011,22(1):85⁃89.[8] 陈玉华,王勇.基于SYSWELD的运行管道在役焊接热循环数值模拟[J].焊接学报,2007,28(1):85⁃89.CHEN Yu⁃hua,WANG Yong.Numerical simulation ofthermal cycle of in⁃service welding onto active pipelinebased on SYSWELD[J].Transactions of the ChinaWelding Institution,2007,28(1):85⁃89.[9] 任瑛,张弘.传热学[M].东营:石油大学出版社,1988.[10] WATT D F.An algorithm for modeling microstructuredevelopment in weld heat⁃affected zones(part A)[J].Acta Metallurgica,1988,36(11):3029⁃3035. [11] 陈玉华,王勇,何建军.输气管线在役焊接管道内壁变形的数值模拟[J].焊接学报,2010,31(1):110⁃112.CHEN Yu⁃hua,WANG Yong,HE Jian⁃jun.Numericalsimulation on deformation in inner pipe wall of in⁃servicewelding onto gas pipeline[J].Transactions of the ChinaWelding Institution,2010,31(1):110⁃112. [12] BANG IW,SON Y P,OH K H,et al.Numerical simu⁃lation of sleeve repair welding of in⁃service gas pipelines[J].Welding Journal,2002,81(12):273⁃282.(编辑 沈玉英)㊃221㊃中国石油大学学报(自然科学版) 2014年6月。
GTAW焊接热场的三维动态数值模拟
GTAW焊接热场的三维动态数值模拟
李瑞英
【期刊名称】《大庆师范学院学报》
【年(卷),期】2009(29)6
【摘要】基于SYSWELD软件平台,采用双椭球体移动热源模式,利用有限元方法建立了运动电弧作用下GTAW焊接不锈钢薄板焊接热场的三维动态有限元分析模型.利用所建模型,对GTAW焊接热场和熔池形状参数的动态演变进行了预测,达到宏
观准稳态的时间与试验结果基本一致.
【总页数】5页(P90-94)
【作者】李瑞英
【作者单位】大庆师范学院物理与电气信息工程系,黑龙江,大庆,163712
【正文语种】中文
【中图分类】TG444+.74
【相关文献】
1.2219-T87铝合金GTAW焊接接头的局部本构关系测量及单向拉伸数值模拟 [J], 李艳军;李权;吴爱萍;麻宁绪;王国庆;HidekazuMURAKAWA;鄢东洋;吴会强
2.基于ANSYS薄板GTAW焊接温度场数值模拟 [J], 卢振洋;刘建;黄鹏飞;吉玲
3.基于ANSYS薄板GTAW焊接温度场数值模拟 [J], 卢振洋;刘建;黄鹏飞;吉玲
4.2219-T87铝合金GTAW焊接接头的局部本构关系测量及单向拉伸数值模拟 [J], 李艳军;李权;吴爱萍;麻宁绪;王国庆;Hidekazu MURAKAWA;鄢东洋;吴会强;
5.外加变位磁场作用GTAW焊接电弧的数值模拟 [J], 周祥曼;刘练;陈永清;袁有录;田启华;杜义贤;何青松;付君健
因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Ｗ ｕ ｕ ２ ０ ， ｈ ｎ ； ． ｈｎ ｓ ａ ｅ ｆ ｎ ｔ ｕｅｏ ｌｓ ｈ ｓｃ ， ｅ ｅ ２ ０ ０ ， ｈｎ ） ｈ ４ ０ Ｃ ｉ ａ ３ Ｃ ｉｅ ｅＡｃ ｄ ｍｙ ｏ ｓｉ ｔ ｆ ａ ｍａ ｙ ｉｓ Ｈ ｆｉ ３ ０ ０ Ｃ ｉａ １ ０ Ｉ ｔ Ｐ Ｐ
Ａｂ ｔａ ｔ Ｂａｅ ｎｔｅｆ ｉ ｌ ｎ ｏ ｗａｅＳ Ｓ Ｅ Ｄ ｔ ｖ ｒ ｐｊ ｉｔｔｍｐ ｒｔｒ ｅｄ ｆｒ３ ｙ ａａ ｓｒ ｃ： ｓｄｏ ｎｔ ｅｅ ｈ ｉ ｅ ｍｅ ｔｓｆ ｒ Ｙ Ｗ Ｌ ｏ ｏ ｅｌ ｏｎ ｅ ｅａｕｅｆ ｌ ｏ ｄｄ ｎ ￣ｉ ｔ ａ ｉ ｃ
第２ ７卷 第 ６ 期
、０ ． ７ Ｎ Ｏ ，１ ２ ．６
新 乡学院学报 ： 自然科学版
Ｊ ｕ ａ ｆ ｎ ｉｎ ｉ ｅ ｓｔ Ｎａ ｕ ａ ｃ ｅ ｃ ｄ ｔ ｎ ｏ ｒ ｌ ｘ ａ ｇ Ｕｎ ｖ ｒ ｉ ｎ ｏ Ｘｉ ｙ： ｔ ｒ ｌ ｉ ｎ ｅ Ｅ ｉ ｏ Ｓ ｉ
（． ｃｏ ｌ ｆ ａｈｎ ｒ ｎ ｔｍｏ ｉ ｎ ｉｅｒ ｇ Ａｎ ｕ ｏｙｅｈ ｉ Ｕｎｖ ｒｉ ， ｈ ４ ０ ０ Ｃ ｉａ １Ｓ ｈ ｏ Ｍ ｃ ｉｅｙａ ｄＡｕｏ ｂｌ Ｅ ｇｎ ｅｉ ， ｈ ｉ １ｔｃ ｎｃ ｉｅｓ ｙ Ｗｕ ｕ２ １ ０ ， ｈｎ ； ｏ ｅ ｎ Ｐ ｔ
Ｎｕ ｅ ｉ ａ ｉ ｕ ａ ｉ ｎ ｏ ｍ ｐｅ ａ ｕ ｅＦｉ ｌ ｆＯｖ ｒ ａ ｏ ｎ ｓ ｄ ｍ ｒｃ ｌ ｍ ｌ ｔｏ ｆ Ｓ Ｔｅ ｒ ｔ ｒ ｅｄ ｏ ｅ ｌ ｐ Ｊ ｉ ｔ Ｂａ ｅ
ｏ ｈ ＹＳ ＥＬＤ ｎ ｔ ｅＳ Ｗ
Ｚ ＨＡＮＧ ｈ ．ｕ ｎ， ＷＡＮＧ ｈ ｎ － ｅ， ＨＵⅪ ｅｂｎ， ｉ， Ｎ Ｗｅ Ｓ ｕｑ ａ ， ＿ Ｚ ｏ ｇｊ Ｚ － ｕ ．ｉ ＤＡＩ ２ ＷＥ ｉ Ｌ ３
应 力应 变和 减 少焊接 应 力与 变形提供 了参 考依 据 。
关键 词 ：Ｓ Ｓ Ｌ Ｙ ＷＥ Ｄ；搭 接 接 头 温 度 场 ； 焊 接 热 循 环 ；相 演 变
中图 分类 号 ：Ｔ ０ Ｇ４ ２
文 献标 志 码 ：Ａ
文 章编 号 ：１ ７ — ３ ６２ １ ） ６－０ ７ ０ ６４ ３ ２（００ ０ － ６－ ３ ０
摘
要 ：基 于有限元 软件 Ｓ ＷＥ Ｄ 对搭接 接 头温度 场进 行 了 ３维动 态模 拟 ，得 出了瞬 态温度场 分布 图和 ＹＳ Ｌ
特征 点 的热循 环 曲线。 同时通 过对 Ｓ ＷＥ Ｄ 后 处理 分析 ，得 出了焊件 上任 一点 的温度 变化 与相 演 变的 关 ＹＳ Ｌ 系。与 文献 资料 比较表 明 ，所 建立 的数 值模 拟仿 真 模型 可 以较好 地模 拟 焊接 温度 场 ，为研 究 焊接 过程 中的
ｓ ｍｕｌｔｏ ，ｔ ｅ ｔ ａ ｓ ｅ ｔ ｅ ｅ ａ ｕ ｅ ｆｅ ｄ ｄ ｓｒ ｂ ｔｏ ｓ ｗｅ ｌａ ｅ ｔ ｒ ｏ ｎ ｓｔ ｅ ｍ ａ ｙ ｌ ｇ ｃ ｒ ｅ ｉ ｏ ， ｎ ｉ ａ ｉ ｎ ｈ ｒ ｎ ｉ ｎ ｍｐ ｒ ｔ ｒ ｌ ｉ ｔｉ ｕ ｉ ｎ ａ ｌ ｓｆ ａ ｕ ｅ ｐ ｉ ｔ ｈ ｒ ｌｃ ｃ ｉ ｕ ｖ ｓ ｇ ｔ ａ ｄ ｔ ｉ ｎ ｔ ｒ ｕ ｈ ｔ ｅ ｆｎ ｔ ｌｍ ｅ ｔａ ａ ｙ ｉ ， ｈ ｅ ａ ｉ ｎ ｏ ｍ ｐ ｒ ｔ ｒ ｉ ｈ ｉ ｏ ａ ｅ ｎ ａ ｙ ｐ ｉ ｔｏ ｒ ｉ ｃ ｄ ｈ ｏ ｇ ｈ ｉ ｅ ｅ ｎ ｎ ｌ ｓｓ ｔ ｅ ｒ ｌ ｔ ｆｔ ｉ ｅ ｏ ｅ ｅ ａ ｕ ｅ ｗｈ ｃ ｓ ｌ ｃ ｔ ｄ ｉ ｎ ｏ ｎ ｆ ｗｏ ｋ ｐ ｅ ｅ ａ ｎ ｐ ａ ｅ ｔａ ｓ ｏ ａ ｉ ｎ ｉ ｂ ａ ｎ ｄ Ｃｏ ｈ ｓ ｒ ｎ ｆ ｒ ｔｏ ｓｏ ｔ ｉ ｅ ． ｍｐ ｒｎ ｈ ｉ ｌ ｔｏ ｅ ｕ ｔ ｔ ｈ ｉｅ ａ ｕ ｅ ｄ ｔ ， ｔ ｓｄ ｍ ｏ ｓ ｒ ｔ ｄ ｔ ａ ｍ ａ ｉ ｇ ｔ ｅ ｓｍｕ ａ ｉ ｎ ｒ ｓ ｌ ｗｉｈ ｔ ｅ ｌ ｒ ｔ ｒ ａ ａ ｉ ｉ ｅ ｎ ｔａ ｅ ｈ ｔ ｓ ｔ ｔ ｅ ｅ ｔ ｂ ｉ ｈ ｄ ｍｏ ｅ ｆ ｎ ｍｅ ｉ ａ ｉ ｕ ａ ｉ ｎ ｃ ｎ ｐ ｅ ｅ ａ ｌ ｉ ｌｔ ｈ ｅ ｅ ａ ｕ ｅ ｆｅ ｄ ｈ ｓ ａ ｌ ｅ ｄ ｌ ｏ ｕ ｒｃ ｌｓｍ ｌｔｏ ａ ｒ ｆ ｒ ｂ ｙ ｓ ｍｕ ａ ｅ ｔ ｅ ｔ ｍｐ ｒ ｔ ｒ ｌ ，ｗｈ ｃ ｒ ｖ ｄ ｓ ａ ｓ ｉ ｉｈ ｐ ｏ ｉ ｅ ｒ ｆ ｒ ｎ ｅ ｆ ｒａ ａ ｙ ｉ ｇ ｓｒ ｓ ， ｔａｎ ａ ｄ ｄ ｃ ｅ ｓｎ ｈ ｌ ｉ ｇ ｓｒ ｓ ｎ ｅ ｏ ｍａ ｉ ｎ ｉ ｌ ｉ ｇ ｐ ｏ ｅ ｓ ｅ ｅ ｅ ｃ ｏ ｎ ｌ ｚ ｎ ｔｅ ｓ ｓｒ ｉ ｎ ｅ ｒ ａ ｉ ｇ ｔ ｅｗｅ ｄ ｎ ｔ ｅ ｓ ａ ｄ ｄ ｆ ｒ ｔｏ ｗｅ ｄ ｎ ｒ ｃ ｓ ． ｎ
２ Ｄｅ ａ ｔ ｎ ｆＭ ｅ ｈ ｎ ｃ ｌ ｎ ｉ ｅ ｒ ｎ ， ｈ ｉＴ ｃ ｎ ｃ ｌＣｏ ｌ ｇ ｆＭ ｅ ｈ ｎｉ ａ ｎ ｅ ｔｉ ａ ｇ ｎ ｅ ｉ ｇ ． ｐ ｒｍｅ ｔｏ ｃ ａ ｉ ａ ｇ ｎ ｅ ｉ ｇ Ａｎ ｕ ｅ ｈ ｉ ａ ｌ ｅ ｏ ｃ ａ ｃ ｌ ｄ Ｅｌ ｃ ｒｃ ｌ Ｅ ｅ ａ Ｅｎ ｉ ｅ ｒｎ ，Biblioteka ２ １ ０ ０年 １ ２月
Ｄｅ ． ０１ ｃ２ ０
基 于 Ｓ Ｓ Ｌ 的搭接接头温度场 的数值模拟 Ｙ ＷＥ Ｄ
张书权 １ ，王仲 珏 ，朱协 彬 ，代 ， ２ 礼 ２ ，文 伟
（． １安徽 工程 大学 机械 与 汽车 工程 学院 ，安徽 芜湖 ２ １ ０ ；２安徽 机 电职 业技 术 学院 机械 工程 系， ４０ ０ ． 安徽 芜湖 ２ １ ０ ；３中 国科 学院 等 离子 体物 理研 究所 ，安徽 合 肥 ２ ０ ０ ４ ００ ． ３ ０ ０）