建筑工程钢结构焊接过程模拟与焊接变形、焊接ansys应力有限元分析(详细图解分析)

基于ANSYS的焊接温度场和应力的数值模拟研究一、本文概述随着现代工业技术的飞速发展，焊接作为一种重要的连接工艺，在航空、汽车、船舶、石油化工等领域的应用日益广泛。

然而，焊接过程中产生的温度场和应力场对焊接结构的性能有着至关重要的影响。

为了深入理解焊接过程中的热-力行为，预测焊接结构的变形和残余应力，进而优化焊接工艺参数和提高产品质量，本文旨在利用ANSYS有限元分析软件，对焊接过程中的温度场和应力场进行数值模拟研究。

本文首先简要介绍了焊接数值模拟的意义和现状，包括焊接数值模拟的重要性、国内外研究现状和存在的问题等。

随后，详细阐述了ANSYS 软件在焊接数值模拟中的应用，包括其基本原理、分析流程、模型建立、参数设置等方面。

在此基础上，本文以某典型焊接结构为例，详细阐述了焊接温度场和应力场的数值模拟过程，包括模型的建立、边界条件的设定、求解参数的选择、结果的后处理等。

对模拟结果进行了详细的分析和讨论，验证了数值模拟方法的准确性和可靠性，为实际工程应用提供了有益的参考。

本文的研究不仅有助于深入理解焊接过程中的热-力行为，为优化焊接工艺参数和提高产品质量提供理论支持，同时也为ANSYS软件在焊接数值模拟领域的应用推广和进一步发展奠定了基础。

二、焊接理论基础焊接是一种通过加热、加压或两者并用，使两块或多块金属在原子层面结合形成永久性连接的工艺过程。

焊接过程涉及复杂的物理和化学变化，包括金属的熔化、凝固、相变以及应力和变形的产生等。

因此，深入了解焊接过程的理论基础对于准确模拟焊接过程中的温度场和应力分布至关重要。

焊接过程中，热源将能量传递给工件，导致工件局部快速升温并熔化。

熔池形成后，随着热源的移动，熔池中的液态金属逐渐凝固形成焊缝。

焊接热源的类型和移动速度、工件的材质和厚度等因素都会影响焊接过程的温度场分布。

为了准确模拟这一过程，需要了解各种热源模型（如移动热源模型、体积热源模型等）及其适用范围，并选择合适的模型进行数值模拟。

有限元分析在钢结构工程施工中的应用摘要：现阶段，我国的综合国力不断地提高，人们的生活水平也越来越高。

为满足人们文化及精神生活的需求，各种大型建筑应运而生。

尤其，近年来各种空间钢结构不断涌现，如网架结构、桁架结构、网壳结构等广泛应用于实际工程中。

对大跨空间钢结构而言，由于其结构施工过程复杂，施工方法和施工工艺繁琐，在施工阶段出现风险的概率要比其他结构高。

运用有限元分析，可以在钢结构施工过程中进行计算机模拟跟踪计算，为施工过程提供安全精确的数值分析结果和动态模拟。

关键词：有限元分析；钢结构工程施工；应用引言随着我国经济的发展，大跨度空间钢结构的形式也日趋复杂，施工过程对结构的影响不能忽略。

用施工力学的方法对施工过程进行预分析，不仅可以优选结构施工方案，而且保证施工过程中结构的安全性以及竣工状态结构的内力和位形满足设计要求。

本文基于ANSYS、MARC等大型有限元平台上，并充分考虑施工步骤，等的影响，对结构施工进行跟踪模拟分析。

1有限元方法及软件介绍有限元法可以称为有限单元法或有限元素法，基本思想是将物体（即连续求解域）离散成有限个且按一定方式相互连接在一起的单元组合，来模拟和逼近原来的物体，从而将一个连续的无限自由度问题简化为离散的有限自由度问题求解的数值分析法。

结构在施工过程中是逐层承受荷载的，并引起结构相应的内力和变形，每次对结构施加荷载时，结构便形成刚度，便产生内力与变形。

当增加下一结构时，所施加的荷载与原来形成的荷载一起影响结构的变形与内力，这样不停地变化，内力与变形也在不停地发生变化，每次形成矩阵不断地迭代求解，有限元则是采用单元生死技术来控制结构的先后顺序，模拟变形，得到所需要的结果。

结构施工建模步骤如下：（1）建立构件三维空间有限元模型，形成结构整体刚度矩阵；根据施工步骤划分施工阶段，分阶段建模。

（2）利用有限元软件ANSYS的单元生死技术钝化所有施工步（包括构件及其相应的边界条件、荷载和约束），先将整体结构建模，按照施工的顺序，将未建造结构单元的刚度矩阵乘以一个很小的缩减因子，即单元生死系数，这样单元就处于失效的状态下；（3）将单元载荷、质量、应变和刚度设为0值，未建结构单元的质量、刚度对已建结构不产生任何影响。

《有限元法及其应用》课程作业ANSYS应用分析学号：姓名：专业：建筑与土木工程角托架的有限元建模与分析一 、模型介绍本模型是关于一个角托架的简单加载，线性静态结构分析问题，托架的具体形状和尺寸如图所示。

托架左上方的销孔被焊接完全固定，其右下角的销孔受到锥形压力载荷，角托架材料为Q235A 优质钢。

角托架材料参数为：弹性模量366E e psi =；泊松比0.27ν=托架图（厚度：0.5）二、问题分析因为角托架在Z 方向尺寸相对于其在X,Y 方向的尺寸来说很小，并且压力荷载仅作用在X,Y 平面上，因此可以认为这个分析为平面应力状态。

三、模型建立3.1 指定工作文件名和分析标题（1）选择菜单栏Utility Menu →File →Jobname 命令.系统将弹出Jobname(修改文件名)对话框,输入bracket（2）定义分析标题GUI ：Utility Menu>Preprocess>Element Type>Add/Edit/Delete 执行命令后，弹出对话框，输入stress in a bracket 作为ANSYS 图形显示时的标题。

3.2设置计算类型Main Menu: Preferences … →select Structural → OK3.3定义单元类型PLANE82 GUI ：Main Menu →Preprocessor →Element Type →Add/Edit/Delete 命令，系统将弹出Element Types 对话框。

单击Add 按钮，在对话框左边的下拉列表中单击Structural Solid →Quad 8node 82，选择8节点平面单元PLANE82。

单击ok ，Element Types 对话框，单击Option ，在Element behavior 后面窗口中选取Plane strs w/thk 后单击ok 完成定义单元类型。

∙基于ANSYS焊接残余应力有限元分析技术研究∙以岭澳核电站控制棒驱动机构耐压壳Ω环焊接修复为例，应用ANSYS有限元生死单元技术模拟焊接流程，计算出焊接后残余应力的分布，绘制出残余应力分布曲线，并与美国WSI公司的计算结果进行对比分析。

结果表明，本课题的计算结果与美国焊接公司<WSI公司）一致。

因此，焊接残余应力有限元分析技术可以用于反应堆耐压壳焊接修复评价。

1 前言焊接在工业中的应用是不言而喻的，但同时焊接过程中产生的残余应力往往又会导致焊接失效。

因此，在工业中一般都要对残余应力进行消除，但这种消应力处理往往在实际结构或环境中难以实现，就必须进行破坏性分析。

随着我国核反应堆的建设及运行，核级设备及管道会出现较多的缺陷，有的缺陷必须进行打磨后焊接修复，同时要进行力学分析评价，此时，力学分析就必须考虑由焊接而产生的残余应力。

对于焊接后结构中的残余应力大小及分布，会因结构形式、焊接方式及材料特性的不同而不同。

某核电站控制棒驱动机构(CRDM >耐压壳上部Ω环连续两年都出现了泄漏，并在检修期间进行焊接修复。

焊接公司委托美国公司对修复后的结构进行了力学分析和评定。

焊接残余应力的有限元计算是关键技术之一，也是难点。

通过本课题的研究，掌握有限元模拟焊接过程及残余应力计算，能够提高我国焊接修复工程缺陷的分析能力，优化不符合项的处理程序，达到既节约时间和资金又满足工作性能和安全性能的目的。

因此，进行焊接残余应力有限元分析技术的研究是非常有必要的。

2 焊接实例本文以某核电站CRDM 耐压壳Ω焊接为研究对象，分析研究焊接后的残余应力分布。

CRDM 耐压壳包括上段是驱动杆行程套管和下段的密封壳。

驱动杆行程套管与密封壳采用螺纹连接，Ω焊接密封的结构进行连接和密封。

驱动杆行程套管的上端采用端塞，通过螺纹连接，Ω焊接密封的结构进行密封。

CRDM 耐压壳采用的这种密封结构形式是一种便于拆装的焊接密封结构，由于其内力的整体平衡主要由连接螺纹承担，Ω焊缝功能上主要起密封作用。

焊接过程模拟与焊接变形、焊接Ansys应力有限元分析1.1 焊接变形与焊接应力焊接时，加热和冷却循环总会导致一定程度的变形，焊接变形对尺寸稳定性以及结构力学性能都有很大的影响，控制焊接变形在焊接加工中是一个关键的任务。

在钢结构焊接中，焊接工艺会使构件温度场产生不均匀变化，从而在构件中产生复杂的残余应力分布。

残余应力是一种自相平衡的力系，当构件承受荷载时，如受拉、受压等，荷载引起的应力将与截面残余应力相叠加，从而使构件某些部位提前达到屈服强度，并发生塑性变形，故会严重降低构件的刚度和稳定性以及结构疲劳强度。

对构件进行焊接，在焊件上产生局部高温的不均匀温度场，焊接中心处温度可达1600℃，高温区的钢材会发生较大程度的膨胀伸长，但受到相邻钢材的约束，从而在焊件内引起较高的温度应力，并在焊接过程中，随时间和温度而不断变化，称其为焊接应力。

焊接应力较高的部位，甚至将达到钢材的屈服强度而发生塑性变形，因而钢材冷却后将有残存于焊件内的应力，称为焊接残余应力。

并且在冷却过程中，钢材由于不能自由收缩，而受到拉伸，于是焊件中出现了一个与焊件加热方向大致相反的内应力场。

1.2 Ansys有限元焊接分析为通过对焊接过程的三维有限元模拟分析以及焊接后构件变形及残余应力分布分析，为评估焊接对焊件的影响提供更加合理、有效、可靠的分析数据，并为焊接工艺提供一定的指导，为采用的焊接过程提供一定的分析依据，采用大型有限元计算软件Ansys作为分析工具对焊接过程与焊件的变形与残余应力进行了分析。

ANSYS有2种方式来考虑热分析与力学分析之间的耦合,即直接耦合和间接耦合。

间接耦合法的处理思路为先进行温度场的模拟,然后将求出的结点温度作为体载荷施加在结构中,计算焊接残余应力与变形。

即：(1)使用热分析的手段进行热分析，根据需要可采用瞬态分析与稳态分析模型，此处为瞬态分析。

(2)重新进入前处理中,将热分析单元转换为相应的结构分析单元，设置结构分析中材料属性,如弹性模量、泊松比、热膨胀系数等。

焊接温度场仿真和热变形、应力仿真的基本理论和仿真流程1 前言焊接作为现代制造业必不可少的工艺，在材料加工领域一直占有重要地位。

焊接是一个涉及到电弧物理、传热、冶金和力学等各学科的复杂过程，其涉及到的传热过程、金属的融化和凝固、冷却时的相变、焊接应力和变形等是企业制造部门和设计人员关心的重点问题。

焊接过程中产生的焊接应力和变形，不仅影响焊接结构的制造过程，而且还影响焊接结构的使用性能。

这些缺陷的产生主要是焊接时不合理的热过程引起的。

由于高能量的集中的瞬时热输入，在焊接过程中和焊后将产生相当大的残余应力和变形，影响结构的加工精度和尺寸的稳定性。

因此对于焊接温度场合应力场的定量分析、预测有重要意义。

传统的焊接温度场和应力测试依赖于设计人员的经验或基于统计基础的半经验公式，但此类方法带有明显的局限性，对于新工艺无法做到前瞻性的预测，从而导致实验成本急剧增加，因此针对焊接采用数值模拟的方式体现出了巨大优势。

ANSYS作为世界知名的通用结构分析软件，提供了完整的分析功能，完备的材料本构关系，为焊接仿真提供了技术保障。

文中以ANSYS为平台，阐述了焊接温度场仿真和热变形、应力仿真的基本理论和仿真流程，为企业设计人员提供了一定的参考。

2 焊接数值模拟理论基础焊接问题中的温度场和应力变形等最终可以归结为求解微分方程组，对于该类方程求解的方式通常为两大类：解析法和数值法。

由于只有在做了大量简化假设，并且问题较为简单的情况下，才可能用解析法得到方程解，因此对于焊接问题的模拟通常采用数值方法。

在焊接分析中，常用的数值方法包括：差分法、有限元法、数值积分法、蒙特卡洛法。

差分法：差分法通过把微分方程转换为差分方程来进行求解。

对于规则的几何特性和均匀的材料特性问题，编程简单，收敛性好。

但该方法往往仅局限于规则的差分网格（正方形、矩形、三角形等），同时差分法只考虑节点的作用，而不考虑节点间单元的贡献，常常用来进行焊接热传导、氢扩散等问题的研究。

《基于ANSYS的焊接温度场和应力的数值模拟研究》篇一一、引言随着科技的发展，焊接技术作为制造行业中的关键工艺之一，其质量和效率直接关系到产品的性能和寿命。

因此，对焊接过程中的温度场和应力分布进行精确的数值模拟显得尤为重要。

ANSYS作为一种功能强大的工程仿真软件，被广泛应用于焊接过程的数值模拟。

本文将基于ANSYS，对焊接温度场和应力进行数值模拟研究，以期为实际生产提供理论依据。

二、焊接温度场的数值模拟1. 模型建立在ANSYS中建立焊接过程的有限元模型，包括焊件、焊缝、热源等部分。

其中，焊件采用实体单元进行建模，焊缝则通过线单元进行描述。

热源模型的选择对于模拟结果的准确性至关重要，应根据具体的焊接工艺选择合适的热源模型。

2. 材料属性及边界条件根据实际材料，设定焊件和焊缝的热导率、比热容、热扩散率等物理参数。

同时，设定初始温度、环境温度等边界条件。

3. 数值模拟过程根据焊接过程的实际情况，设定加载步和时间步长，模拟焊接过程中的温度变化。

通过ANSYS的热分析模块，得到焊接过程中的温度场分布。

三、焊接应力的数值模拟1. 耦合分析焊接过程中，温度场的变化会导致应力的产生。

因此，在ANSYS中，需要将在热分析中得到的温度场结果作为应力分析的输入条件，进行热-结构耦合分析。

2. 本构关系与材料模型根据材料的本构关系和力学性能，设定材料的弹性模量、泊松比、热膨胀系数等参数。

同时，选择合适的材料模型，如各向同性模型或各向异性模型。

3. 应力分析通过ANSYS的结构分析模块，结合耦合后的温度场结果，进行应力分析。

得到焊接过程中的应力分布和变化情况。

四、结果与讨论1. 温度场结果分析通过ANSYS的后处理功能，可以得到焊接过程中的温度场分布图。

分析温度场的分布情况，可以了解焊接过程中的热传导和热扩散情况，为优化焊接工艺提供依据。

2. 应力结果分析同样，通过后处理功能可以得到焊接过程中的应力分布图。

分析应力的分布和变化情况，可以了解焊接过程中产生的残余应力和变形情况。

焊接模拟ansys实例!下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进!行焊接过程仿真, 计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。

finish/clear/filnam,1-2D element birth and death/title,Weld Analysis by "Element Birth and Death"/prep7/unit,si !采用国际单位制!******************************************************et,1,13,4 !13号二维耦合单元, 同时具有温度和位移自由度et,2,13,4!1号材料是钢!2号材料是铝!3号材料是铜!铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来!下面是在几个温度点下, 各材料的弹性模量mptemp,1,20,500,1000,1500,2000mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11!假设各材料都是双线性随动硬化弹塑性本构关系!下面是各材料在各温度点下的屈服应力和屈服后的弹性模量tb,bkin,1,5tbtemp,20,1tbdata,1,1200e6,0.193e11tbtemp,500,2tbdata,1, 933e6,0.150e11tbtemp,1000,3tbdata,1, 435e6,0.070e11tbtemp,1500,4tbdata,1, 70e6,0.010e11tbtemp,2000,5tbdata,1, 7e6,0.001e11!tb,bkin,2,5tbtemp,20,1tbdata,1,800e6,0.102e11tbtemp,500,2tbdata,1,400e6,0.050e11tbtemp,1000,3tbdata,1, 70e6,0.008e11tbdata,1, 1e6,0.0001e11tbtemp,2000,5tbdata,1,0.1e6,0.00001e11!tb,bkin,3,5tbtemp,20,1tbdata,1,900e6,0.117e11tbtemp,500,2tbdata,1,700e6,0.090e11tbtemp,1000,3tbdata,1,230e6,0.030e11tbtemp,1500,4tbdata,1, 40e6,0.005e11tbtemp,2000,5tbdata,1, 4e6,0.0005e11!!材料密度(假设为常值)mp,dens,1,8030mp,dens,2,4850mp,dens,3,8900! 热膨胀系数(假设为常值)mp,alpx,1,1.78e-5mp,alpx,2,9.36e-6mp,alpx,3,1.66e-5!泊松比(假设为常值)mp,nuxy,1,0.29mp,nuxy,2,0.30mp,nuxy,3,0.30!热传导系数(假设为常值)mp,kxx,1,16.3mp,kxx,2,7.44mp,kxx,3,393!比热(假设为常值)mp,c,1,502mp,c,2,544mp,c,3,385!热膨胀系数(假设为常值)!由于该13号单元还有磁自由度, 此处假设一磁特性, 但并不影响我们所关心的结果mp,murx,1,1mp,murx,2,1mp,murx,3,1!假设焊料(铝)焊上去后的初始温度是1500℃mp,reft,2,1500mp,reft,3,20!******************************************************!下面建立几何模型csys,0k,1,0,0,0k,2,0.5,0,0k,3,1,0,0 !长1米k,4,0,0.3,0 !厚度0.3米（二维中叫做宽度）k,5,0.35,0.3,0k,6,0.65,0.3,0k,7,1,0.3,0a,1,2,5,4a,2,6,5a,2,3,7,6!!!!!!!!!!!!!!!!!!!!!!!!!!划分网格!!!!!!!!!!!!!!!!!esize,0.025type,2mat,2amesh,2!esize,0.05 !网格划分出现问题type,1mat,1amesh,1!mat,3amesh,3eplot!/soluantype,4 ! 瞬态分析trnopt,full!!!!!!!!!!!!!!!!!!!!!!!!在模型的左边界加位移约束!!!!!!!!!!!!!!!!!!!!!!!!!!!nsel,all*get,minx,node,,mnloc,xnsel,s,loc,x,minxd,all,ux,0*get,miny,node,,mnloc,ynsel,r,loc,y,minyd,all,uy,0!*****假设模型的左右边界处温度始终保持在20摄氏度左右*****!其他边界条件如对流和辐射等均可施加，此处因为只是示意而已，故只施加恒温边界条件nsel,all*get,minx,node,,mnloc,xnsel,s,loc,x,minxd,all,temp,20nsel,all*get,maxx,node,,mxloc,xnsel,s,loc,x,maxxd,all,temp,20!!!!!!!!!!!!!!!由于第2个面是焊接所在区域，因此首先将该区域的单元“杀死”!!!!!!!!!!!!!!!!nna=2esel,all*get,emax,elem,,num,maxasel,s,area,,nnaesla*get,nse,elem,,count*dim,ne,,nse*dim,nex,,nse*dim,ney,,nse*dim,neorder,,nsemine=0!**********************************************!下面的do循环用于将焊料区的单元按其形心y坐标排序!以便后面模拟焊料由下向上逐步“生长”过程*do,i1,1,nseesel,u,elem,,mine*get,nse1,elem,,countii=0*do,i,1,emax*if,esel(i),eq,1,thenii=ii+1ne(ii)=i*endif*enddo*do,i,1,nse1*get,ney(i),elem,ne(i),cent,y*get,nex(i),elem,ne(i),cent,x*enddominy=1e20minx=1e20*do,i,1,nse1*if,ney(i),lt,miny,thenminy=ney(i)minx=nex(i)mine=ne(i)*else*if,ney(i),eq,miny,then*if,nex(i),lt,minx,thenminy=ney(i)minx=nex(i)mine=ne(i)*endif*endif*endif*enddoneorder(i1)=mine*enddo!************************************************************** max_tem=1500 !按照前面假设，焊料的初始温度为1500℃dt1=1e-3 !用于建立初始条件的一个很小的时间段dt=5 !焊接一个单元所需的时间t=0 !起始时间esel,alleplot/auto,1/replot*do,i,1,nseekill,neorder(i)esel,s,liveeplot*enddoallsel,alloutres,all,allic,all,temp,20kbc,1timint,0,structtimint,1,thermtimint,0,magtintp,0.005,,,1,0.5,0.2!nsub1=2nsub2=40!**************************************************do,i,1,nseealive,neorder(i)esel,s,liveeplotesel,all!******下面的求解用于建立温度的初始条件******t=t+dt1time,tnsubst,1*do,j,1,4d,nelem(neorder(i),j),temp,max_tem*enddosolve!****下面的求解用于保证初始的升温速度为零****t=t+dt1time,tsolve!*********下面的步骤用于求解温度分布***********do,j,1,4ddele,nelem(neorder(i),j),temp*enddot=t+dt-2*dt1time,tnsubst,nsub1solve*enddot=t+50000 !*********下面的步骤用于冷却过程求解***** time,tnsubst,nsub2solvesavefinish!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!后处理过程!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/post1!**************下面的一系列命令用于生成应力的动画文件******* /seg,dele/cont,1,15,0,1200e6/16,1200e6/dscale,1,1.0avprin,0,0avres,1/seg,multi,stress1,0.1esel,all*do,i,1,nseesel,u,elem,,neorder(i)*enddo*do,i,1,nseesel,a,elem,,neorder(i)set,(i-1)*3+1,1plnsol,s,eqv*do,j,1,nsub1set,(i-1)*3+3,jplnsol,s,eqv*enddo*enddo*do,i,1,nsub2set,(nse-1)*3+4,iplnsol,s,eqv*enddo/seg,off,stress1,0.1/anfile,save,stress1,avi!**********下面的一系列命令用于生成温度的动画文件************ /seg,dele/cont,1,15,0,1500/16,1500/dscale,1,1.0avprin,0,0avres,1/seg,multi,temp1,0.1esel,all*do,i,1,nseesel,u,elem,,neorder(i)*enddo*do,i,1,nseesel,a,elem,,neorder(i)set,(i-1)*3+1,1plnsol,temp*do,j,1,nsub1set,(i-1)*3+3,jplnsol,temp*enddo*enddo*do,i,1,nsub2set,(nse-1)*3+4,iplnsol,temp*enddo/seg,off,temp1,0.1/anfile,save,temp1,avifinish。

基于ANSYS的某焊接件两焊缝在顺序焊接过程中的分析（生死单元应用案例）焊接几何模型如下图所示，左右两侧90度扇区为焊接材料，其余为钢板材料。

其他更多已知条件请参考命令流，这里不再赘述。

网格单元本实例中顺序焊接分为如下步骤：第一步0-1秒：右侧焊接稳态分析（杀死左焊缝，施加右焊缝温度和焊接件参考温度）第二步1-100秒：相变分析（删除温度载荷，施加对流热传导）第三步100-1000秒：右侧焊缝凝固分析第四步1000-1001秒：激活左侧焊缝单元进行稳态分析（施加左焊缝温度）第五步1001-1100秒：左焊缝相变分析第六步1100-2000秒：左侧焊缝凝固分析第七步：结果后处理ANSYS命令流：FINISH/FILNAME,Exercise ! 定义隐式热分析文件名/PREP7 ! 进入前处理器ET,1,SOLID70 ! 选择8节点实体热分析单元MP,KXX,1,.5e-3MP,C,1,.2MP,DENS,1,.2833MPTEMP,1,0,2643,2750,2875,3000MPDATA,ENTH,1,1,0,128.1,163.8,174.2,184.6 ! 定义右焊缝材料热物理性能MP,KXX,2,.5e-3MP,C,2,.2MP,DENS,2,.2833MP,KXX,3,0.5e-3 ! 定义两块钢板的热物理性能MP,DENS,3,.2833MPTEMP,1,0,2643,2750,2875,3000MPDATA,ENTH,3,1,0,128.1,163.8,174.2,184.6 ! 定义左焊缝材料热物理性能BLOCK,-0.17,0.17,0,0.34,0,1.2BLOCK,0.17,0.34,0,0.34,0,1.2BLOCK,0.34,1,0,0.34,0,1.2BLOCK,-0.17,0.17,0.34,0.51,0,1.2BLOCK,-0.17,0.17,0.51,1.34,0,1.2WPAVE,0.17,0.34,0CYLIND,0.17,0,0,1.2,0,90WPAVE,0,0,0CSYS,0FLST,3,3,6,ORDE,3FITEM,3,2FITEM,3,-3FITEM,3,6VSYMM,X,P51X, , , ,0,0 ! 建立焊接件的几何模型VGLUE,ALL ! 粘接各体VSEL,S,,,10VATT,1,1,1 ! 附于右焊缝的材料属性VSEL,S,,,1VSEL,A,,,12,17,1VATT,2,1,1 ! 附于两块钢板的材料属性VSEL,S,,,11VATT,3,1,1 ! 附于左焊缝的材料属性ALLSEL,ALLESIZE,0.05 ! 定义单元划分尺寸VSWEEP,ALL ! 划分单元ESEL,S,MAT,,3TOFFST,460 ! 定义温度偏移量!第一步：稳态分析EKILL,ALL ! 杀死左焊缝单元ALLSEL,ALL/SOLUANTYPE,TRANS ! 定义瞬态分析类型TIMINT,OFF ! 关闭时间积分ESEL,S,MAT,,1NSLE ! 选择右焊缝节点D,ALL,TEMP,3000 ! 施加右焊缝初始温度载荷NSEL,INVE ! 选择其它节点D,ALL,TEMP,70 ! 施加初始温度载荷TIME,1 ! 定义求解时间KBC,0 ! 设置为斜坡载荷SOLVE ! 求解!第二步：右侧焊缝相变分析（1到100秒)DDELE,ALL,TEMP ! 删除温度载荷TIMINT,ON ! 打开时间积分TINTP,,,,1 ! 定义瞬态积分参数TIME,100 ! 定义求解时间DELTIME,1,.5,10 ! 定义时间子步AUTOTS,ON ! 打开自动时间开关KBC,1 ! 设置为阶越载荷OUTRES,ERASEOUTRES,ALL,ALL ! 设置结果输出ASEL,S,EXTASEL,U,LOC,Y,0SFA,ALL,,CONV,5E-5,70 ! 施加对流换热载荷ALLSEL,ALLSOLVE ! 求解!第三步：右侧焊缝凝固分析（100到1000秒)TIME,1000 ! 定义求解时间DELTIME,50,10,100 ! 定义时间子步AUTOTS,ON ! 打开自动时间开关SOLVE ! 求解!第四步：激活左侧焊缝单元进行分析（1000到1001秒) EALIVE,ALL ! 激活左侧焊缝单元ALLSEL,allESEL,S,MAT,,3NSLE ! 选择左焊缝节点D,ALL,TEMP,3000 ! 施加左焊缝初始温度载荷TIME,1001 ! 定义求解时间DELTIME,1,1,1 ! 定义时间子步ALLSEL,ALLSOLVE ! 求解!第五步：左侧焊缝相变分析（1001到1100秒) DDELE,ALL,TEMP ! 删除温度载荷TIME,1100 ! 定义求解时间DELTIME,1,.5,10 ! 定义时间子步SOLVE ! 求解!第六步：左侧焊缝凝固分析（1100到2000秒) TIME,2000 ! 定义求解时间DELTIME,100,10,200 ! 定义时间子步SOLVE ! 求解!第七步：后处理/POST1 ! 进入通用后处理器SET,,,,,1, , ! 读取1秒中分析结果PLNSOL, TEMP,, 0 ! 显示1秒钟后焊接件的温度分布SET,,,,,100, , ! 读取100秒中分析结果PLNSOL, TEMP,, 0 ! 显示100秒钟后焊接件的温度分布SET,,,,,1000, , ! 读取1000秒中分析结果PLNSOL, TEMP,, 0 ! 显示1000秒钟后焊接件的温度分布SET,,,,,1001, , ! 读取1001秒中分析结果PLNSOL, TEMP,, 0 ! 显示1001秒钟后焊接件的温度分布SET,,,,,1100, , ! 读取1100秒中分析结果PLNSOL, TEMP,, 0 ! 显示1100秒钟后焊接件的温度分布SET,,,,,2000, , ! 读取2000秒中分析结果PLNSOL, TEMP,, 0 ! 显示2000秒钟后焊接件的温度分布/POST26 ! 进入时间历程后处理器NSOL,2,4727,TEMP,, TEMP_2STORE,MERGENSOL,3,4752,TEMP,,TEMP_3STORE,MERGENSOL,4,4808,TEMP,,TEMP_4STORE,MERGENSOL,5,4833,TEMP,,TEMP_5STORE,MERGENSOL,6,4883,TEMP,,TEMP_6STORE,MERGENSOL,7,4908,TEMP,,TEMP_7STORE,MERGENSOL,8,5088,TEMP,,TEMP_8STORE,MERGENSOL,9,5308,TEMP,, TEMP_9STORE,MERGE ! 定义焊接件某些位置8个节点的时间温度变量/AXLAB,X,TIME/AXLAB,Y,TEMPERATURE ! 更改坐标轴标识/XRANGE,0,2000 ! 设定横坐标轴范围PLVAR,2,3,4,5,6,7,8,9, ! 绘制8节点温度随时间的变化曲线温度结果显示1秒时：100秒时：1000秒时：1001秒时：1100秒时：2000秒时：关键点温度变化曲线关键点位置分布：关键点与曲线对应情况：NSOL,2,4727,TEMP,,TEMP_2NSOL,3,4752,TEMP,,TEMP_3NSOL,4,4808,TEMP,,TEMP_4NSOL,5,4833,TEMP,,TEMP_5NSOL,6,4883,TEMP,,TEMP_6NSOL,7,4908,TEMP,,TEMP_7NSOL,8,5088,TEMP,,TEMP_8NSOL,9,5308,TEMP,,TEMP_9关键点温度变化曲线：从该图中可以明显看到：离焊接热影响区域距离较远的关键点温升较小，距离较近的关键点温升较大。