ansys 双椭球热源焊接APDL
毕业论文使用ANSYS有限元分析软件对不锈钢与紫铜焊接过程温度场应力场分布模拟
三、由于在焊接过程中不锈钢管和紫铜管应力分布有很大的不同，因此在冷却结束后两管会产生不同的体积变化。这对于焊缝是不利的，同时也会使焊件不符合焊接前的尺寸公差，导致零件报废。建议焊前对焊件采取反应力、反变形措施。
Firstly, because the heat transferring of copper is faster than stainless steel, so temperature field distribution of thin-walled stainless steel and coppertubesduring the welding process is extremely uneven. Thismayleads to molten inconsistency of base metal on both sides and may cause poor weldseam. It is recommended that the welding arc isslightly partial to copper, in case oftheexcessive heating of stainless steel.
Thirdly, due to the large difference of stress distribution between stainless steel tube and copper tube during the welding process, different volume changes are produced in two tubes at the end of cooling. This is detrimental to the welding seam and also leads to unfitness of dimensional tolerances of welding parts, resulting in scrapping of welding parts. It is recommended that anti-stress and anti-distortion measures should be adopted before welding.
ansys软件在焊接领域的应用
ANSYS软件在焊接领域的应用和前景A NSYS是计算机辅助工程(CA E)领域应用最广泛的有限元分析软件,通过对该软件的系统组成、工作流程和工作原理等方面进行分析,对有限元方法FEM在焊接热效应领域的应用进行了研究。
ANSYS 能与其他主流CAD软件双向传递数据,具有多物理场分析能力和便捷的前后数据处理能力,通过基于ANSYS的虚拟试验平台,可以低成本、高效率优化与焊接热效应相关产品设计方案,因而在焊接研究和生产方面有着广阔的应用前景。
焊接是一个涉及许多学科的复杂的物理—化学过程。
由于焊接过程涉及的变量数目繁多，单凭积累工艺试验数据来深入了解和控制焊接过程，既不切实际又成本昂贵和费时费力。
随着计算机技术的发展，通过一组描述焊接基本物理过程的数学方程来模拟焊接过程，采用数值方法求解以获得焊接过程的定量认识，即焊接过程的计算机模拟，成为一种强有力的手段。
计算机模拟方法为焊接科学技术的发展创造了有力的条件。
1993 年，美国能源部组织美国、加拿大、日本、瑞典、英国的 25 位著名专家对 21 世纪焊接科学技术的发展动向做出预测，其中焊接基本现象的模拟与仿真被列为最重要的研究方向之一。
我国国家自然科学基金委员会制定的学科发展战略也将计算机模拟确定为机械热加工领域的发展方向之一。
计算机模拟是使包括焊接在内的热加工工艺研究从“定性”走向“定量”、从“经验”走向“科学”的重要标志。
采用科学的模拟技术和少量的实验验证，以代替过去一切都要通过大量重复实验的方法，不仅可以节省大量的人力和物力，而且还可以通过数值模拟解决一些目前无法在实验室里进行直接研究的复杂问题。
在制造业，计算机模拟与仿真可以增加材料利用率 25 ％，节约生产成本 30 ％，产品设计至实际投产的时间缩短 40 ％。
焊接过程数值模拟中，热源拟合，温度场的模拟是最基本的工作，然后就是应力和变形的模拟。
我们可以看到大量这方面的文章，温度场的模拟起步也较早，也积累了比较丰富的经验，在实际生产中得到了一定的应用。
ANSYS APDL使用指南
第三章 使用参数........................................................................................................................................... 5
3.1 参数 ...................................................................................................................................................... 5 3.2 参数命名规则 ...................................................................................................................................... 5
3.3.1 在运行过程中给参数赋值 ........................................................................................................... 6 3.3.2 在启动时给参数赋值 ................................................................................................................... 7 3.3.3 将 ANSYS 提供的数据赋值给参数.............................................................................................. 7 3.3.4 列表显示参数............................................................................................................................. 10 3.4 删除参数 ............................................................................................................................................ 11 3.5 使用字符参数 .................................................................................................................................... 11 3.6 数值参数值的替换 ............................................................................................................................ 12 3.6.1 禁止替换..................................................................................................................................... 12 3.6.2 字符参数值的替换 ..................................................................................................................... 12 3.7 数字或字符参数的动态替换 ............................................................................................................ 14 3.8 参数表达式 ........................................................................................................................................ 14 3.9 参数函数 ............................................................................................................................................ 15 3.10 保存、恢复和写出参数 .................................................................................................................. 16 3.11 数组参数 .......................................................................................................................................... 17 3.11.1 数组的基础知识 ....................................................................................................................... 17 3.11.2 数组参数举例 ........................................................................................................................... 18 3.11.3 TABLE 类型数组参数................................................................................................................ 19 3.11.4 定义和列表数组参数 ............................................................................................................... 21 3.11.5 给数组元素赋值 ........................................................................................................................ 21 3.11.5.6 插值 ........................................................................................................................................ 28 3.11.6 写出数据文件 ........................................................................................................................... 31 3.11.7 数组参数之间的运算 ............................................................................................................... 33 3.11.8 图形显示数组参数矢量 ........................................................................................................... 40 3.11.9 修改曲线标识字 ....................................................................................................................... 42
Ansys-APDL-基础
APDL 根底
INTRODUCTION TO ANSYS 5.7 - Part 1
APDL 根底 概述
Training Manual
• APDL 是 ANSYS 参数化设计语言的 的缩写,它是一种允许使用参 数并能完成一系列任务的强大的程序语言。
• 使用 APDL, 您可以: • 用参数而不是用数值输入模型尺寸，材料类型等。 • 从 ANSYS 数据库中猎取信息, 比方节点位置或最大应力。 • 在参数中进展数学运算，包括矢量和矩阵运算。 • 把常用的命令或宏定义成缩写形式。 • 建立一个宏使用if-then-else分支和do循环等来执行一系列任务。
用 *SET 看有用参数列表
Training Manual
g=386 massdens=density/g circumf=2*pi*rad area=pi*r**2 dist=sqrt((y2-y1)**2+(x2-x1)**2) slope=(y2-y1)/(x2-x1) theta=atan(slope) jobname=‘proj1’
January 30, 2001 Inventory #001441
12-5
INTRODUCTION TO ANSYS 5.7 - Part 1
INTRODUCTION TO ANSYS 5.7 - Part 1
APDL 根底 ...定义参数
Training Manual
• 以上例子是关于标量参数的, 它只有一个值 — 数字或者字符。
fk,2,fy,-force ! 2号关键点的力
fk,6,fx,force/2
! 6号关键点的力
Training Manual
January 30, 2001 Inventory #001441
连续油管环焊缝温度场模拟
《热加工工艺》2013年4月第42卷第7期连续油管环焊缝温度场模拟雷青娟,李宵,胡莎莎,武岳(西安石油大学材料科学与工程学院,陕西西安710065)摘要:利用有限元软件ANSYS 建立了连续油管对接环焊缝的三维焊接温度场的计算机数值分析模型。
分析中引入了双椭球热源分布模型,研制了相应的APDL 程序,计算出打底焊的温度场曲线。
模拟计算结果和试验结果基本吻合,表明该计算方法能准确地计算出焊接温度场。
关键词:环焊缝;双椭球;温度场;数值分析中图分类号:TG444;TG402文献标识码:A文章编号:1001-3814(2013)07-0175-02Simulation on Ring Welding Temperature Field for Ring WeldingSeam of Coiled TubingLEI Qingjuan,LI Xiao,HU Shasha,WU Yue(School of Material Science and Engineering,Xi'an Shiyou University,Xi'an 710065,China)Abstract :The coiled tubing ring seam 3D welding temperature field numerical analysis model was set up by finite element software ANSYS.The double ellipsoid heat source distribution model was imported,and the corresponding APDL programs were developed.The backing welding temperature field curve was calculated.The simulation results are basically consistent with experiment results,which show that the method can accurately calculate the welding temperature field.Key words :ring welding;double ellipsoid;temperature field;numerical analysis连续油管是相对于常规螺纹连接油管而言的,又称连续管、挠性油管、蛇形管或盘管,具有高强度和高韧性[1]。
ANSYS APDL热分析入门
APDL热分析关键知识及实例一.关键概念（1）λ：热导率，是指当温度垂直向下梯度为1℃/m时，单位时间内通过单位水平截面积所传递的热量。
（2）E: 弹性模量，材料在弹性变形阶段，其应力和应变成正比例关系（即符合胡克定律），其比例系数称为弹性模量。
（3）Β：热胀系数，物体由于温度改变而有胀缩现象。
其变化能力以等压(p一定)下，单位温度变化所导致的长度量值的变化，即热膨胀系数表示。
各物体的热膨胀系数不同，一般金属的热膨胀系数单位为1/度（摄氏）。
（4）μ：泊松比，指材料在单向受拉或受压时，横向正应变与轴向正应变的绝对值的比值，也叫横向变形系数，它是反映材料横向变形的弹性常数。
（5）α：传热系数、膜系数，单位时间通过单位面积传递的热量。
（6）T u: 接触温度，材料与外界接触处温度。
（7）C: 热容，“当一系统由于加给一微小的热量δQ而温度升高dT时，δQ/dT 这个量即是该系统的热容。
”（8）q: 热通量，单位时间内，通过物体单位横截面积上的热量。
（9）ε：发射系数，原子谱线中发射谱线的辐射能量可用一个发射系数来表示，其含义为单位时间单位体积单位立体角内辐射的能量。
●传热三种基本方式：热传导、热对流及热辐射。
●热流率/热流量(Heat flow)表示单位时间内，通过传导，对流，辐射的方式穿过给定表面传输的热量，也称为热流量。
常表示为Φ，国际单位为瓦特(W)。
这是一种热学上荷载，即热量，相当于功率。
如果大于零，表示热量流入，物体获得热量，反之，热量外流。
●热流密度/热通量（Heat Flux）一般用q表示，定义为单位时间内，通过物体单位横截面积上的热量。
二．基本代号（热力学基本符号）●APDL关键缩略写K：关键点L：线E：单元DENS：密度MAT：材料ET：单元类型KXX：热导率EX：弹性模量HF:传热系数PRXY：泊松比ALPX：热胀系数REFT：参考温度GXY：剪切模量MU：摩擦系数REAL：实常数MP：材料类型PRIN：主应力SINT：应力强度SEQV：等效应力IC：初始条件三．主要知识1.常用分析单元：MASS71:Thermal MassLINK31:Radiation LinkLINK33: convection LinkLINK34:conduction LinkPLANE35:Thermal SolidPLANE55:Thermal SolidSHELL131:Themal ShellPLANE223:Coupled_filed Solid2.热分析的类型ANSYS支持两种类型的热分析：(1)稳态热分析确定在稳态的条件下的温度分布及其他热特性，稳态条件指热量随时间的变化可以忽略。
ansys-APDL编程-图文
ansys-APDL编程-图文目录第一章APDL是什么？第二章在工具条上添加命令2.1修改工具条2.2嵌套工具条缩写第三章使用参数3.1参数3.2参数命名规则3.2.1从某STATUS命令中隐藏参数3.3定义参数3.3.1在运行过程中给参数赋值3.3.2在启动时给参数赋值3.3.3赋ANSYS提供的值给参数3.3.3.1某GET命令的用法3.3.3.2内嵌获取函数的用法3.3.4排列显示参数3.4删除参数3.5字符参数的用法3.6数字参数值的置换3.6.1防止置换3.6.2字符参数值的置换3.6.2.1强制置换3.6.2.2字符参数有效的其它地方3.6.2.3字符参数的限制3.7数字或字符参数的动态置换3.8参数公式3.9带参数的函数3.10保存、恢复、写参数3.11数组参数3.11.1数组的基础知识3.11.2数组参数示例3.11.3TABLE类型数组参数3.11.4定义和列表显示数组参数3.11.5给数组元素赋值3.11.5.1给单独的数组元素赋值3.11.5.2填充数组向量12245556666778101111121212121314141515161717181920202021 I3.11.5.3交互式编辑数组3.11.5.4使用某VREAD命令用数据文件填充数组3.11.5.5使用某TREAD命令用数据文件填充TABLE类型数组3.11.6插入值3.11.6.1把获取值存入数组参数或恢复数组参数值3.11.6.2列出数组参数3.11.7写数据文件3.11.7.1数据格式描述符3.11.8对数组参数的运算3.11.8.1对向量的运算3.11.8.2矩阵运算3.11.8.3用于向量和矩阵运算的命令3.11.9用图形表示数组参数向量第四章作为宏语言的APDL4.1什么是APDL宏4.2产生宏4.2.1宏文件命名规则4.2.2宏的搜索路径4.2.3在ANSYS中生成宏4.2.3.1使用某CREATE4.2.3.2使用某CFWRITE4.2.3.3使用UtilityMenu>Macro>CreateMacro4.2.4用文本编辑器生成宏4.2.5使用宏库文件4.3运行宏和宏库文件4.4局部变量4.4.1传递变量到宏4.4.2宏内的局部变量4.4.3宏外部的局部变量4.5在APDL中控制程序流4.5.1宏嵌套：在宏内调用子程序4.5.2无条件分支：GOTO4.5.3条件分支：某IF命令4.5.4重复一个命令4.5.5循环：DO循环4.6控制函数快速参考4.7在宏中使用_STATUS和_RETURN参数4.8在组和组件中使用宏4.9复习宏例子第五章GUI用户界面212223262728292931313335374141414142434343444445454646474747474 84850505052535457II5.1提示用户输入某个参数的值5.2用户提示对话框5.3用宏显示消息5957585.4在宏中生成并维护状态条5.5在宏中进行拾取操作5.6在宏中调用对话框第六章加密宏6.1准备加密宏6.2生成加密宏6.3运行加密宏60616263636364III第一章APDL是什么？APDL即ANSYS参数化设计语言（ANSYSParametricDeignLanguage），它是一种解释性语言，可用来自动完成一些通用性强的任务，也可以用于根据参数来建立模型。
双椭球热源模型加载算法研究
热过程，Goldak 提出了双椭球形热源模型，该模型能比较灵活地处 理电弧挺度对焊接过程的影响。双椭球形热源模型，如图 1 所示。
Flux q z
y b
a２
O
xc
a１
模型的热源分布规律，分别采用三层嵌套循环计算前沿和后沿内 各节点的 HGEN。
3.2 HGEN 自动加载算法
焊接热分析必须实现热源移动，且这一过程通常自动进行。 规则离散后的热源移动通过热源中心的移动参数化实现。确定热 源中心的位置后，将前沿及后沿的 HGEN 分别加载到对应的节 点上，加载求解流程，如图 3 所示。
MB = l2P
（6）
弯变形过程中不同位置上应力的分布情况，将弯曲线划分为四 段，即曲率弹塑性段、弹塑性自由段、纯弹性自由段、无变形自由 段。利用悬臂梁弯曲理论和材料理想弹塑性模型建立了各个分段 的数学方程。利用边界条件采用数值方法可以确定这些方程中的
MA = l1 +1 MB l2 由图 1几何关系得：
线研究和表示压弯，这条曲线为大挠度平面弯曲线。通过分析压 ＰｒｏｃｅｓｓｉｎｇＴｅｃｈｎｏｌｏｇｙ，Ｖｏｌｕｍｅ９９，Ｉｓｓｕｅｓ１～３，１Ｍａｒｃｈ，２０００：１０３～１１２
第 11 期
董克权等：双椭球热源模型加载算法研究 *
61
２ 算法的相关模型
2.1 双椭球热源模型
有关试验研究表明，三维的热源模型计算能真实体现焊接的
（9）
θ = ψ + θ1 + θ2
（10）
４ 结论
参数，从而确定压弯过程中大挠度平面弯曲线。为进一步的分析 和回弹分析奠定基础。
参考文献
１ 韩志仁． 压弯过程的参数辨识：［硕士学位论文］． 沈阳：沈阳航空工业学院， １９８９（２）：１～３
ANSYS APDL命令汇总
ANSYS APDL命令汇总AA,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。
点要依次序输入，输入的顺序会决定面的法线方向。
如果超过四个点，则这些点必须在同一个平面上。
Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPsABBR*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．ABBRESABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．ABBSAVABBSAV，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认ADDadd, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称ADELEAdele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。
ADRAGAdrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成AFILLTAfillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。
双椭球热源参数调整在预测在役焊接熔池尺寸上的应用
郭广飞 ， 王 勇， 韩 涛， 贾鹏 宇 ２ ５ ７ ０ ６ １ ） （ 中 国石 油大 学 （ 华 东 ）机 电工 程 学 院 ， 山东 东 营
摘 要： 双椭 球 热源模 型常用 于在 役 焊接 过程 的数 值模 拟 ， 模 型 中参 数 的合理 选 用 对 于 计算 结 果 的 正确性 具 有决定 性 的作 用 。 以往 进行 热 源参 数 的调 整 ， 常运 用 经验反 复进 行试 算 以确定 热源 参 数 ，
ｗ ｅ ｌ ｄ ｉ ｎ ｇ ． Ａｎ ｅ ｍｐ ｉ ｉ ｒ ｃ ａ ｌ ｆ ｏ ｒ ｍｕ ｌ ａ ｏ ｆ ｐ ｏ ｏ ｌ ｓ ｉ ｚ ｅ ｐ ａ ｒ ａ ｍｅ ｔ ｅ ｒ ｓ ｏ ｆ ｄ ｏ ｕ ｂ ｌ ｅ ｅ ｌ ｌ ｉ ｐ ｓ ｏ ｉ ｄ ｈ ｅ ａ ｔ ｍｏ ｄ ｅ ｌ ｗ ａ ｓ ｅ ｓ ｔ ａ ｂ ｌ ｉ ｓ ｈ ｅ ｄ ｂ ｙ ｔ ｈ ｅ ｍｅ ａ ｎ ｏ ｆ ｃ ｏ ｍｂ ｉ ｎ ａ ｔ ｉ ｏ ｎ ｏ ｆ ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ａ ｎ ｄ ｅ ｘ ｐ ｅ ｉ ｒ ｍｅ ｎ ｔ １． ａ Ａ ｎ ｄ ｔ ｈ ｅ ｒ ｅ ｓ ｕ ｌ ｔ ｓ ｓ ｈ ｏ ｗ ｔ ｈ ａ ｔ ， ｗｉ ｔ ｈ ｔ ｈ ｅ ｍｅ ｔ ｈ ｏ ｄ， ｔ ｈ ｅ ｅ ｉｃ ｆ ｉ ｅ ｎ ｃ ｙ ｔ ｏ ｄ ｅ ｔ ｅ ｒ ｍｉ ｎ ｅ ｔ ｈ ｅ ｈ ｅ ａ ｔ ｐ ａ ｒ ａ ｍｅ ｔ ｅ ｒ ｓ ｉ ｓ ｉ ｍｐ ｒ ｏ ｖ ｅ ｄ ｓ ｉ ｇ ｎ ｉ ｆ ｉ ｃ ａ ｎ ｔ ｌ ｙ ， ａ ｎ ｄ ｔ ｈ ｅ ｓ ｉ ｚ ｅ ｏ ｆ ｉ ｎ —ｓ ｅ ｖｉ ｒ ｃ ｅ ｗｅ ｌ ｄ ｉ ｎ ｇ
（ Ｃ ｏ ｌ ｌ ｅ ｇ ｅ ｏ ｆ Ｍ ｅ ｃ ｈ ａ ｎ ｉ ｃ ａ ｌ ａ ｎ ｄ Ｅ ｌ ｅ ｃ ｔ ｒ ｏ ｎ ｉ ｃ Ｅ ｎ ｇ ｉ ｎ ｅ ｅ ｉ ｒ ｎ ｇ ｉ ｎ Ｃ ｈ ｉ ｎ ａ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ｏ ｆ Ｐ ｅ ｔ ｒ ｏ ｌ ｅ ｕ ｍ， Ｄ ｏ ｎ ｇ ｙ ｉ ｎ ｇ ２ ５ ７ ０ ６ １ ， Ｃ ｈ ｉ ｎ ａ ）
完整word版,ANSYS焊接高斯热源APDL参考.doc
不准/UNITS,SI/CONFIG ,NRES,10000/PREP7ET,1,SOLID70MP,DENS,1,7930MP,C,1,502mptemp,1,20,100,500mpdata,kxx,1,1,12.1,16.3,21.4/VIEW,1,1,1,1!*************定义基板尺寸***********lx=0.1lz=0.1ly=0.006!************* 定义焊接参数***********weld=0.01 ! 每道焊道宽度wheight=0.004 !焊道的高度LSIZE=0.005V=0.008!速度pi=3.1415926U=25I=180Q=U*I!电源有效功率R=0.006! 电源有效热半径effect=0.8Qmax=effect*Q/(pi*R*R) !中心处最大热流tinc=LSIZE/V! 每小段的时间间隔!*****************建立模型，生成网格************** block,0,lx,0,ly,0,lzwpoff,0,ly,lz/2block,0,lx,0,wheight,-weld/2,weld/2wpoff,0,0,-weld/2vsbw,all,,deletewpoff,0,0,weldvsbw,all,,deletevglue,alllsel,s,,,plsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35lsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35lesize,all,0.005,,,0.2,,,,1lsel,s,line,,9,10lsel,a,line,,12lsel,a,line,,38lesize,all,0.005,,,0.2,,,,1lsel,s,line,,2lsel,a,line,,4,5lsel,a,line,,7lsel,a,line,,14lsel,a,line,,16lsel,a,line,,19lsel,a,line,,39lesize,all,0.002,,,1,,,,1lsel,s,line,,23,24lsel,a,line,,41,42lsel,a,line,,45lsel,a,line,,47lsel,s,line,,23,24lsel,a,line,,41,42lsel,a,line,,45lsel,a,line,,47lesize,all,0.002,,,1,,,,1 vmesh,all/PNUM,DEFAEPLOTfini/soluantype,trans,new!nlgeom,ontimint,0,structtimint,1,thermtimint,0,magtimint,0,electtref,25nropt,autoautos,onkbc,0pred,onlnsrch,on!************杀死焊缝去单元***********nsel,s,loc,y,ly,ly+wheightesln,s,1ekill,allallsel,allesel,s,liveeplot!**********施加对流载荷*************esel,s,livensel,s,loc,z,0nsel,a,loc,z,lznsel,a,loc,x,0nsel,a,loc,x,lxsf,all,conv,10,25!************定义数组维数 **********MAX_X=1+lx/LSIZE!*************定义 table 数组 ****************************************** *do,i,1,MAX_X,1local,12,0,(i-1)*lsize,ly,lz/2,,,,,,*del,_FNCNAME*del,_FNCMTID*del,_FNC_C1*del,_FNC_C2*del,_FNC_C3*del,_FNC_C4*del,_FNCCSYS*del,'heatflux'!定义表格各行*set,_FNCNAME,'heatflux'*dim,_FNC_C1,,1*dim,_FNC_C2,,1*dim,_FNC_C3,,1*dim,_FNC_C4,,1!表格各行赋值*set,_FNC_C1(1),effect*set,_FNC_C2(1),U*set,_FNC_C3(1),I*set,_FNC_C4(1),R*set,_FNCCSYS,12*DIM,%_FNCNAME%,TABLE,6,26,1,,,,%_FNCCSYS%!!Begin of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) %_FNCNAME%(0,0,1)= 0.0, -999 %_FNCNAME%(2,0,1)= 0.0%_FNCNAME%(3,0,1)= %_FNC_C1(1)% %_FNCNAME%(4,0,1)= %_FNC_C2(1)% %_FNCNAME%(5,0,1)= %_FNC_C3(1)% %_FNCNAME%(6,0,1)= %_FNC_C4(1)% %_FNCNAME%(0,1,1)= 1.0, -1, 0, 3, 0, 0, 17%_FNCNAME%(0,2,1)= 0.0, -2, 0, 1, -1, 3, 17%_FNCNAME%(0,3,1)=0, -1, 0, 1, -2, 3, 18%_FNCNAME%(0,4,1)= 0.0, -2, 0, 1, -1, 3, 19%_FNCNAME%(0,5,1)= 0.0, -1, 0, 0, 0, 0, 0%_FNCNAME%(0,6,1)= 0.0, -3, 0, 1, 0, 0, -1%_FNCNAME%(0,7,1)= 0.0, -4, 0, 1, -1, 2, -3%_FNCNAME%(0,8,1)= 0.0, -1, 0, 3, 0, 0, -4%_FNCNAME%(0,9,1)= 0.0, -3, 0, 1, -4, 3, -1%_FNCNAME%(0,10,1)= 0.0, -1, 0, 2, 0, 0, 2%_FNCNAME%(0,11,1)= 0.0, -4, 0, 1, 2, 17, -1%_FNCNAME%(0,12,1)= 0.0, -1, 0, 2, 0, 0, 4%_FNCNAME%(0,13,1)= 0.0, -5, 0, 1, 4, 17, -1%_FNCNAME%(0,14,1)= 0.0, -1, 0, 1, -4, 1, -5%_FNCNAME%(0,15,1)= 0.0, -4, 0, 1, -3, 3, -1%_FNCNAME%(0,16,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,17,1)= 0.0, -3, 0, 1, 20, 17, -1%_FNCNAME%(0,18,1)= 0.0, -1, 0, 1, -4, 4, -3%_FNCNAME%(0,19,1)= 0.0, -1, 7, 1, -1, 0, 0%_FNCNAME%(0,20,1)= 0.0, -3, 0, 1, -2, 3, -1%_FNCNAME%(0,21,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,22,1)= 0.0, -2, 0, 1, 20, 17, -1%_FNCNAME%(0,23,1)= 0.0, -1, 0, 3.14159265358979310, 0, 0, -2 %_FNCNAME%(0,24,1)= 0.0, -4, 0, 1, -1, 3, -2%_FNCNAME%(0,25,1)= 0.0, -1, 0, 1, -3, 4, -4%_FNCNAME%(0,26,1)= 0.0, 99, 0, 1, -1, 0, 0!End of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) !-->!**********激活单元 *********esel,s,livensel,s,loc,x,(i-1)*lsize-0.002,(i-1)*lsize+0.002nsel,r,loc,z,-weld/2,weld/2nsel,r,loc,y,0,wheightesln,s,0ealive,allallsel,all!**********施加热流载荷********** esel,s,livensel,s,loc,Y,wheightsf,all,hflux,%heatflux%allsel,alltime,i*tincnsubst,2SOLVE!*************删除热流载荷*********** nsel,s,loc,y,wheightsfdele,all,hfluxallsel,allesel,s,liveeplotOUTRES,ALL,ALL,/PSF,HFLUX,,2/REPLOTOUTRES,ALL,ALL,。
ANSYS二次开发APDL语言
ANSYS二次开发APDL语言参数化程序设计语言（APDL）参数化程序设计语言(APDL:ANSYS Parametric Design Language)实质上由类似于FORTRAN77的程序设计语言部分和1000多条ANSYS命令组成。
其中，程序设计语言部分与其它编程语言一样，具有参数、数组表达式、函数、流程控制（循环与分支）、重复执行命令、缩写、宏以及用户程序等。
标准的AN SYS程序运行是由1000多条命令驱动的，这些命令可以写进程序设计语言编写的程序，命令的参数可以赋确定值，也可以通过表达式的结果或参数的方式进行赋值。
从ANSYS命令的功能上讲，它们分别对应ANSYS分析过程中的定义几何模型、划分单元网格、材料定义、添加载荷和边界条件、控制和执行求解和后处理计算结果等指令。
用户可以利用程序设计语言将ANSYS命令组织起来，编写出参数化的用户程序，从而实现有限元分析的全过程，即建立参数化的CAD 模型、参数化的网格划分与控制、参数化的材料定义、参数化的载荷和边界条件定义、参数化的分析控制和求解以及参数化的后处理。
宏是具有某种特殊功能的命令组合，实质上是参数化的用户小程序，可以当作ANSY S的命令处理，可以有输入参数或没有输入参数。
缩写是某条命令或宏的替代名称，它与被替代命令或宏存在一一对应的关系，在AN SYS中二者是完全等同的，但缩写更符合用户习惯，更易于记忆，减少敲击键盘的次数。
ANSYS工具条就是一个很好的缩写例子。
用户界面设计语言（UIDL）标准ANSYS交互图形界面可以驱动ANSYS命令，提供命令的各类输入参数接口和控制开关，用户在图形驱动的级别上进行有限元分析，整个过程变得直观轻松。
用户图形界面设计语言（UIDL）就是编写或改造ANSYS图形界面的专用设计语言，主要完成以下三种图形界面的设计：主菜单系统及菜单项对话框和拾取对话框帮助系统通过用户界面设计语言（UIDL），用户可以在扩充ANSYS功能的同时建立起对应的图形驱动界面，如在主菜单的某位置增加菜单项，设计对应的对话框、拾取对话框，实现参数的输入和其它程序运行的控制，同时提供相应的联机帮助，使操作者能方便地获取系统帮助。
ANSYS Mechanical在焊接仿真中的应用详细解析
ANSYS Mechanical在焊接仿真中的应用详细解析1 前言焊接作为现代制造业必不可少的工艺，在材料加工领域一直占有重要地位。
焊接是一个涉及到电弧物理、传热、冶金和力学等各学科的复杂过程，其涉及到的传热过程、金属的融化和凝固、冷却时的相变、焊接应力和变形等是企业制造部门和设计人员关心的重点问题。
焊接过程中产生的焊接应力和变形，不仅影响焊接结构的制造过程，而且还影响焊接结构的使用性能。
这些缺陷的产生主要是焊接时不合理的热过程引起的。
由于高能量的集中的瞬时热输入，在焊接过程中和焊后将产生相当大的残余应力和变形，影响结构的加工精度和尺寸的稳定性。
因此对于焊接温度场合应力场的定量分析、预测有重要意义。
传统的焊接温度场和应力测试依赖于设计人员的经验或基于统计基础的半经验公式，但此类方法带有明显的局限性，对于新工艺无法做到前瞻性的预测，从而导致实验成本急剧增加，因此针对焊接采用数值模拟的方式体现出了巨大优势。
ANSYS作为世界知名的通用结构分析软件，提供了完整的分析功能，完备的材料本构关系，为焊接仿真提供了技术保障。
文中以ANSYS为平台，阐述了焊接温度场仿真和热变形、应力仿真的基本理论和仿真流程，为企业设计人员提供了一定的参考。
2 焊接数值模拟理论基础焊接问题中的温度场和应力变形等最终可以归结为求解微分方程组，对于该类方程求解的方式通常为两大类：解析法和数值法。
由于只有在做了大量简化假设，并且问题较为简单的情况下，才可能用解析法得到方程解，因此对于焊接问题的模拟通常采用数值方法。
在焊接分析中，常用的数值方法包括：差分法、有限元法、数值积分法、蒙特卡洛法。
差分法：差分法通过把微分方程转换为差分方程来进行求解。
对于规则的几何特性和均匀的材料特性问题，编程简单，收敛性好。
但该方法往往仅局限于规则的差分网格（正方形、矩形、三角形等），同时差分法只考虑节点的作用，而不考虑节点间单元的贡献，常常用来进行焊接热传导、氢扩散等问题的研究。
(整理)ANSYS焊接高斯热源APDL参考.
不准/UNITS,SI/CONFIG,NRES,10000/PREP7ET,1,SOLID70MP,DENS,1,7930MP,C,1,502mptemp,1,20,100,500mpdata,kxx,1,1,12.1,16.3,21.4/VIEW,1,1,1,1!*************定义基板尺寸***********lx=0.1lz=0.1ly=0.006!*************定义焊接参数***********weld=0.01 !每道焊道宽度wheight=0.004 !焊道的高度LSIZE=0.005V=0.008 !速度pi=3.1415926U=25I=180Q=U*I !电源有效功率R=0.006 !电源有效热半径effect=0.8Qmax=effect*Q/(pi*R*R) !中心处最大热流tinc=LSIZE/V !每小段的时间间隔!*****************建立模型，生成网格************** block,0,lx,0,ly,0,lzwpoff,0,ly,lz/2block,0,lx,0,wheight,-weld/2,weld/2wpoff,0,0,-weld/2vsbw,all,,deletewpoff,0,0,weldvsbw,all,,deletevglue,alllsel,s,,,plsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35lsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35 lesize,all,0.005,,,0.2,,,,1lsel,s,line,,9,10lsel,a,line,,12lsel,a,line,,38 lesize,all,0.005,,,0.2,,,,1lsel,s,line,,2lsel,a,line,,4,5lsel,a,line,,7lsel,a,line,,14lsel,a,line,,16lsel,a,line,,19lsel,a,line,,39 lesize,all,0.002,,,1,,,,1lsel,s,line,,23,24lsel,a,line,,41,42 lsel,a,line,,45lsel,a,line,,47lsel,s,line,,23,24lsel,a,line,,41,42 lsel,a,line,,45lsel,a,line,,47 lesize,all,0.002,,,1,,,,1vmesh,all/PNUM,DEFA EPLOTfini/soluantype,trans,new!nlgeom,ontimint,0,struct timint,1,therm timint,0,magtimint,0,electtref,25nropt,autoautos,onkbc,0pred,onlnsrch,on!************杀死焊缝去单元***********nsel,s,loc,y,ly,ly+wheightesln,s,1ekill,allallsel,allesel,s,liveeplot!**********施加对流载荷*************esel,s,livensel,s,loc,z,0nsel,a,loc,z,lznsel,a,loc,x,0nsel,a,loc,x,lxsf,all,conv,10,25!************定义数组维数**********MAX_X=1+lx/LSIZE!*************定义table数组****************************************** *do,i,1,MAX_X,1local,12,0,(i-1)*lsize,ly,lz/2,,,,,,*del,_FNCNAME*del,_FNCMTID*del,_FNC_C1*del,_FNC_C2*del,_FNC_C3*del,_FNC_C4*del,_FNCCSYS*del,'heatflux'!定义表格各行*set,_FNCNAME,'heatflux'*dim,_FNC_C1,,1*dim,_FNC_C2,,1*dim,_FNC_C3,,1*dim,_FNC_C4,,1!表格各行赋值*set,_FNC_C1(1),effect*set,_FNC_C2(1),U*set,_FNC_C3(1),I*set,_FNC_C4(1),R*set,_FNCCSYS,12*DIM,%_FNCNAME%,TABLE,6,26,1,,,,%_FNCCSYS%!! Begin of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) %_FNCNAME%(0,0,1)= 0.0, -999%_FNCNAME%(2,0,1)= 0.0%_FNCNAME%(3,0,1)= %_FNC_C1(1)%%_FNCNAME%(4,0,1)= %_FNC_C2(1)%%_FNCNAME%(5,0,1)= %_FNC_C3(1)%%_FNCNAME%(6,0,1)= %_FNC_C4(1)%%_FNCNAME%(0,1,1)= 1.0, -1, 0, 3, 0, 0, 17%_FNCNAME%(0,2,1)= 0.0, -2, 0, 1, -1, 3, 17%_FNCNAME%(0,3,1)= 0, -1, 0, 1, -2, 3, 18%_FNCNAME%(0,4,1)= 0.0, -2, 0, 1, -1, 3, 19%_FNCNAME%(0,5,1)= 0.0, -1, 0, 0, 0, 0, 0%_FNCNAME%(0,6,1)= 0.0, -3, 0, 1, 0, 0, -1%_FNCNAME%(0,7,1)= 0.0, -4, 0, 1, -1, 2, -3%_FNCNAME%(0,8,1)= 0.0, -1, 0, 3, 0, 0, -4%_FNCNAME%(0,9,1)= 0.0, -3, 0, 1, -4, 3, -1%_FNCNAME%(0,10,1)= 0.0, -1, 0, 2, 0, 0, 2%_FNCNAME%(0,11,1)= 0.0, -4, 0, 1, 2, 17, -1%_FNCNAME%(0,12,1)= 0.0, -1, 0, 2, 0, 0, 4%_FNCNAME%(0,13,1)= 0.0, -5, 0, 1, 4, 17, -1%_FNCNAME%(0,14,1)= 0.0, -1, 0, 1, -4, 1, -5%_FNCNAME%(0,15,1)= 0.0, -4, 0, 1, -3, 3, -1%_FNCNAME%(0,16,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,17,1)= 0.0, -3, 0, 1, 20, 17, -1%_FNCNAME%(0,18,1)= 0.0, -1, 0, 1, -4, 4, -3%_FNCNAME%(0,19,1)= 0.0, -1, 7, 1, -1, 0, 0%_FNCNAME%(0,20,1)= 0.0, -3, 0, 1, -2, 3, -1%_FNCNAME%(0,21,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,22,1)= 0.0, -2, 0, 1, 20, 17, -1%_FNCNAME%(0,23,1)= 0.0, -1, 0, 3.14159265358979310, 0, 0, -2%_FNCNAME%(0,24,1)= 0.0, -4, 0, 1, -1, 3, -2%_FNCNAME%(0,25,1)= 0.0, -1, 0, 1, -3, 4, -4%_FNCNAME%(0,26,1)= 0.0, 99, 0, 1, -1, 0, 0! End of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) !-->!**********激活单元*********esel,s,livensel,s,loc,x,(i-1)*lsize-0.002,(i-1)*lsize+0.002nsel,r,loc,z,-weld/2,weld/2nsel,r,loc,y,0,wheightesln,s,0ealive,allallsel,all!**********施加热流载荷********** esel,s,livensel,s,loc,Y,wheightsf,all,hflux,%heatflux%allsel,alltime,i*tincnsubst,2SOLVE!*************删除热流载荷*********** nsel,s,loc,y,wheightsfdele,all,hfluxallsel,allesel,s,liveeplotOUTRES,ALL,ALL,/PSF,HFLUX,,2/REPLOTOUTRES,ALL,ALL,*ENDDO。
ANSYS高级分析之APDL基础
ANSYS高级分析之APDL基础ANSYS是一款广泛使用的工程仿真软件，它可以进行各种复杂的物理和工程分析。
其中，ANSYS Parametric Design Language（APDL）是ANSYS的一种基于命令行交互的脚本语言，它可以用于创建和控制各种物理模型，并进行高级分析。
APDL语言主要通过输入一系列的命令来操作ANSYS软件。
在使用APDL进行高级分析之前，我们需要先了解一些基础知识。
APDL中的命令可以分为几个主要的类别，包括几何命令、物理命令、边界条件命令和求解命令等。
几何命令用于创建和修改几何模型，比如绘制线段、圆弧和矩形等。
物理命令用于定义分析的物理性质，比如材料的力学性质、热物性等。
边界条件命令用于设定边界条件，如约束和载荷。
求解命令用于进行数值计算，如求解结构的位移、应力和应变等。
除了常规的命令之外，APDL还提供了一些高级分析的功能。
其中，参数化分析是其中一项重要功能，它可以通过修改输入参数或模型的几何特性，自动执行多个分析，从而得到一系列的结果。
参数化分析可以通过循环和条件语句来实现。
另外，APDL还提供了特殊命令和工具，用于处理大规模模型和复杂的分析问题。
在使用APDL进行高级分析时，需要遵循一些最佳实践。
首先，我们应该仔细设计分析模型，包括选择适当的边界条件和物理参数，并进行合理的离散化。
其次，我们应该对模型进行验证和校准，比较模拟结果与实验数据或已知解进行比较。
最后，我们应该进行后处理，对模拟结果进行分析和解释。
总之，APDL是ANSYS的一种基于命令行交互的脚本语言，它可以用于进行高级分析。
通过使用APDL，我们可以建立复杂的物理模型，并执行各种高级分析。
在使用APDL进行高级分析时，我们应该熟悉APDL的基本命令和语法，合理设计模型和参数，并进行验证和后处理。
只有掌握了APDL的基础知识，我们才能更好地应用ANSYS进行高级分析。
ANSYS中的APDL命令总结
ANSYS中的APDL命令总结在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成一定的ANSYS功能，这些功能一般来说通过菜单操作也能够实现（而那些命令流能够实现，菜单操作实现不了的单个命令比较少见）。
以下命令是结合我自身经验，和前辈们的一些经验而总结出来的，希望对大家有帮助。
（1）．Lsel, type, item, comp, vmin, vmax, vinc, kswp选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元（2）．Nsel, type, item, comp, vmin, vmax, vinc, kabs选择一组节点type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值（3）．Esel, type, item, comp, vmin, vmax, vinc, kabs选择一组单元type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item：Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号（4）. mp, lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c : 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数（5）. 定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg（6）. 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。
关于Ansys workbench热辐射APDL、Enclosure等的一些讨论
关于Ansys workbench热辐射APDL、Enclosure等的一些讨论ANSYS workbench对于热分析初学者来说相对比较容易简单上手，市面上和网络上也有不少的相关教程。
但是初学者在学到workbench中的热辐射一节时，往往会有一些疑惑。
本文将根据作者的一些理解以及大量的验证性例子，为大家进一步剖析一下，希望对初学者有所助益。
初学者学到热辐射这一节时通常会碰到如下几个问题：1、workbench平台中mechanical进行热辐射时，通过点击工具栏的相关快捷按扭仅能完成几何体对空气的热辐射？2、Workbench想要完成两个或者多个几何体之间的热辐射设置，只能通过APDL？APDL里的所谓的“封闭体数量”是指什么？3、Mechanical平台能否通过点击工具栏完成多个体之间的热辐射？4、Enclosure到底怎么设置，指的是什么？下面作者将对以上4个问题逐一进行讲解，并给出相应的计算结果进行验证。
问题一：workbench平台中mechanical进行热辐射时，通过点击工具栏的相关快捷按扭仅能完成几何体对空气的热辐射？讲解：很多初学者在一些参考书里都能看到这样一段话，就是通过点击工具栏仅仅能完成几何体对环境（空气）的热辐射而无法完成2个面或者多个面之间的热辐射。
对此我们先从下边一个例子开始讲起。
当我们模拟一个物体在空气中自然冷却的时候，mechanical中的设置比较简单，只需要设置物体的发射率和环境温度即可。
这里我们假设有一个直径1000mm的半球壳，球壳厚度50mm，初始温度100°，球壳外表面绝热，内表面参与热辐射。
计算它在23°空气中的导热情况，显然最后球体稳定下来的温度应该和环境温度一致，即23°。
相关设置情况和最终结果如下图所示。
仔细的读者可能已经发现在菜单栏radiation热辐射设置中，在correlation的下拉菜单里有一个surface to surface（面对面辐射）的选项。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
finish/clear/BATCH WPSTYLE,,,,,,,,0/FILNAME,Laser Beam Welding,0/TITLE,laser Beam Welding/units,si!!!参数定义L=0.02B1=0.003B2=0.006B3=0.02H=0.0032D=0.001 !焊缝宽度ah=D/2chf=2*ahchb=8*ahbh=0.005rf=2*chf/(chf+chb)rb=2*chb/(chf+chb)N=0.27 !功率有效系数W0=2200 !功率W=W0*N !有效功率/PREP7!定义导热系数ET,1,SOLID70MPTEMP,1,20MPTEMP,2,100MPTEMP,3,200MPTEMP,4,300MPTEMP,5,400MPTEMP,6,500MPTEMP,7,600MPTEMP,8,700MPTEMP,9,800MPTEMP,10,900MPTEMP,11,1000MPTEMP,12,1200MPTEMP,13,1540MPTEMP,14,1700MPTEMP,15,2000MPTEMP,16,3000MPDATA,KXX,1,,6.5 MPDATA,KXX,1,,7 MPDATA,KXX,1,,8 MPDATA,KXX,1,,9 MPDATA,KXX,1,,10 MPDATA,KXX,1,,11 MPDATA,KXX,1,,12 MPDATA,KXX,1,,13 MPDATA,KXX,1,,13.5 MPDATA,KXX,1,,13.8 MPDATA,KXX,1,,13.9 MPDATA,KXX,1,,14 MPDATA,KXX,1,,14.2 MPDATA,KXX,1,,14.4 MPDATA,KXX,1,,14.4 !定义比热容MPTEMP,1,20 MPTEMP,2,100 MPTEMP,3,200 MPTEMP,4,300 MPTEMP,5,400 MPTEMP,6,500 MPTEMP,7,600 MPTEMP,8,700 MPTEMP,9,800 MPTEMP,10,900 MPTEMP,11,1000 MPTEMP,12,1200 MPTEMP,13,1540 MPTEMP,14,1700 MPTEMP,15,2000 MPTEMP,16,3000 MPDATA,C,1,,600 MPDATA,C,1,,620 MPDATA,C,1,,640 MPDATA,C,1,,660 MPDATA,C,1,,680 MPDATA,C,1,,700 MPDATA,C,1,,725 MPDATA,C,1,,750 MPDATA,C,1,,775 MPDATA,C,1,,800 MPDATA,C,1,,805MPDATA,C,1,,830 MPDATA,C,1,,850 MPDATA,C,1,,900 MPDATA,C,1,,900!定义密度MP,DENS,1,4450!!!建立几何模型BLOCK,0,B1,0,L,0,H, BLOCK,B1,B2,0,L,0,H, BLOCK,B2,B3,0,L,0,H, VGLUE,all!!!划分网格esize,h/16,0VSWEEP,1esize,h/2,0VSWEEP,5MSHAPE,1,3DMSHKEY,0VMESH,4FINISH!!!施加载荷STEP=0.002 ! N=L/STEPV=0.004 !速度*DO,I,0,NYY=I*STEPT=(YY+STEP)/V/SOLUANTYPE,4TRNOPT,FULLTUNIF,293,OUTRES,NSOL,LASTTIME,TNROPT,FULL,,ON AUTOTS,1NSUBST,5, , ,1 !!!KBC,1lnsrch,1 !!!TSRES,ERASETINTP,,,,1,,,,,!!!TIMINT,0,structTIMINT,1,THERMTIMINT,0,MAG!!!添加第一个载荷步（焊接阶段，有对流和热源）NSEL,S,EXTNSEL,U,LOC,X,0ASEL,S, , ,5ASEL,U,LOC,X,0LSLA,SNSLL,A,1SF,ALL,CONV,200,293ALLSEL,ALL!NSEL,ALL!*GET,CENTER,NODE,,MNLOC,X!NSEL,S,LOC,X,CENTER!SF,ALL,HFLUX,0!!!双椭球热源的加载!!!前半球节点的选取LOCAL,11,2,0,YY,0, , , ,3,8,NSEL,S,LOC,X,,DNSEL,R,LOC,Y,,90NSEL,R,LOC,Z,,90CM,FRONT,NODE!!!!!!!!!!!!!!!!!!!!!!!!!!!LOCAL,13,0,0,YY,0, , , ,1,1,!!!前半球函数、表格定义*DEL,_FNCNAME*DEL,_FNCMTID*DEL,_FNC_C1*DEL,_FNC_C2*DEL,_FNC_C3*DEL,_FNC_C4*DEL,_FNC_C5*DEL,_FNCCSYS*DEL,'TFRONT' !!!!!!!!*SET,_FNCNAME,'TFRONT'*DIM,_FNC_C1,,1*DIM,_FNC_C2,,1*DIM,_FNC_C3,,1*DIM,_FNC_C4,,1*DIM,_FNC_C5,,1*SET,_FNC_C1(1),rf !!!rf*SET,_FNC_C2(1),w !!!Q*SET,_FNC_C3(1),chf !!!chf*SET,_FNC_C4(1),ah !!!ah*SET,_FNC_C5(1),bh !!!bh*SET,_FNCCSYS,13*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS% *SET,%_FNCNAME%(0,0,1), 0.0, -999*SET,%_FNCNAME%(2,0,1), 0.0*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1 *SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3 *SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1 *SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2 *SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0 *SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0!!!!!!!!!!!!!!!!!!BF,FRONT,HGEN, %TFRONT%CMDELE,FRONTALLSEL,ALL!!!后半球节点的选取LOCAL,12,2,0,YY,0, , , ,6,8,NSEL,S,LOC,X,,DNSEL,R,LOC,Y,-90,0NSEL,R,LOC,Z,,90CM,BEHIND,NODE!!!后半球加载!!!!!!!!!!!!!!!!!!!LOCAL,14,0,0,YY,0, , , ,1,1,!!!后半球函数、表格定义*DEL,_FNCNAME*DEL,_FNCMTID*DEL,_FNC_C1*DEL,_FNC_C2*DEL,_FNC_C3*DEL,_FNC_C4*DEL,_FNC_C5*DEL,_FNCCSYS*DEL,'TBEHIND' !!!!!!!!!*SET,_FNCNAME,'TBEHIND'*DIM,_FNC_C1,,1*DIM,_FNC_C2,,1*DIM,_FNC_C3,,1*DIM,_FNC_C4,,1*DIM,_FNC_C5,,1*SET,_FNC_C1(1),rb !!!rb*SET,_FNC_C2(1),w !!!Q*SET,_FNC_C3(1),chb !!!chb*SET,_FNC_C4(1),ah !!!ah*SET,_FNC_C5(1),bh !!!bh*SET,_FNCCSYS,14 !!! 0*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%*SET,%_FNCNAME%(0,0,1), 0.0, -999*SET,%_FNCNAME%(2,0,1), 0.0*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0 *SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0BF,BEHIND,HGEN, %TBEHIND%CMDELE,BEHINDALLSEL,ALLSOLVEBFDELE,all,HGEN*ENDDOSAVE!!!添加冷却载荷第一步CSYS,0NSEL,S,EXTSF,all,CONV,200,293ALLSEL,ALLTIME,L/V+1NSUBST,5, , ,1!!!DELTIM,0.05,0.01,0.1,1SOLVE!!!添加冷却载荷第二步CSYS,0TIME,L/V+10NSUBST,5, , ,1!!!DELTIM,0.2,0.1,0.4,1SOLVE/post26nsol,6,node(0.000,5*L/10,0.000),temp,,origin!!!在表面沿着远离焊缝中心（+x）方向取点nsol,7,node(0.0005,2*L/4,0),temp,,top1nsol,8,node(0.001,2*L/4,0),temp,,top2nsol,9,node(0.0012,2*L/4,0),temp,,top3nsol,10,node(0.0015,2*L/4,0),temp,,top4plvar,6,7,8,9,10,savefinish。