ansys焊接高斯热源apdl参考
ANSYSAPDL命令汇总(2021年整理精品文档)
(完整版)ANSYSAPDL命令汇总编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整版)ANSYSAPDL命令汇总）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为(完整版)ANSYSAPDL命令汇总的全部内容。
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时低单元点一并删除。
利用APDL语言的焊接温度场参数化分析论文
利用APDL语言的焊接温度场参数化分析论文选择合适的焊接工艺参数，对提高焊接质量和提高生产效率很重要。
长期以来，焊接工艺优化主要依靠阅历积累及试验测试，但试验既繁琐又铺张许多时间，增加了本钱。
随着计算机技术的进展，采纳有限元分析软件进行数值模拟的方法模拟焊接过程为实际焊接供应了理论根据，降低了生产本钱。
ANSYS有限元分析软件就是其中最典型的代表之一。
通过其APDL CANSYS Paramelric Design Language)语言编写的指令流可以建立智能的分析过程，自动完成冗杂的分析计算过程。
本文分析就是采纳APDL语言来编写，以平板对接焊为例。
把不同的焊接工艺参数和平板的儿何参数作为分析变量，分析时只需输入转变的变量值，软件能自动完成整体的平板对接焊温度场数值模拟。
1 实例分析模拟试验中母材为20钢((20钢的高温物理性能从文献中得到)。
模型创建前要确定好参数变量，便于后续指令的编写，避开消失混乱。
分析时会消失所不的参数值输入框，便利修转变量数值。
在ANSYS中进行焊接模拟时，焊接热源模型的选择是重中之重，本文选择生死单元热源加载模型，该模型热源能够模拟焊缝的形成过程。
即在开头计算前，将全部焊缝单元“杀死”，相当于焊前的装配状态。
计算时，按挨次将被“杀死”的单元逐一“激活”，模拟焊缝形成过程，同时给激活的.单元施加生热率，每一步计算完成后，删除该步的生热率，重新进入下一步加载计算a。
山于生死单元不2 结果分析不同的参数值，温度场的改变状况不同。
分别是转变焊接电流((150 A,160 A,170 A,180 A,190 A)、焊接电压(15V,16V,17V,18V、19 V)和坡口角度(300,400,500,600,700)温度场中最高温度的改变状况。
山图中可以看出，随着焊接电流、焊接电压的增加，温度场的最高温度也在增加;随着坡口角度的增加，温度场的最高温度在下降。
3 结语本文通过对焊接温度场的数值模拟，阐述了APDL语言的参数化分析方法。
ANSYS焊接仿真中高斯热源加载算法研究
第２期 ２００７ 年 ２ 月
机械设计与制造 Ｍａｃｈｉｎｅｒｙ Ｄｅｓｉｇｎ ＆ Ｍａｎｕｆａｃｔｕｒｅ
文章编号： １００１－ ３９９７（２００７）０２－ ００８３－ ０２
Ａ Ｎ Ｓ Ｙ Ｓ 焊接仿真中高斯热源加载算法研究 ＊
－ ８３ －
董克权 １， ２ 刘超英 １ 陈英俊 １ （ １ 广东肇庆学院 电子工程系， 肇庆 ５２６０６１） （ ２ 广东工业大学 计算机学院， 广州 ５１００９０）
１ 算法的相关模型
称绝热条件。为减少计算工作量， 只取模型的一半进行分析。
１．１ 高斯热源模型
有关试验研究表明， 高斯热源模型能够表征焊接电弧的热 流分布特征［３］。高斯热源模型如图 １ 所示。
图中， Ｏ 为热源的中心， ｄＨ 为焊接电弧有效加热直径。高斯 热源模型的热流密度公式为：
ｑ＝ｑｍａｘｅｘｐ（－ ｃｒ２）
４．２ 单面双丝焊仿真消耗时间比较
在相同的软硬件环境下， 取大小不同的模型， 焊接方向为 模型的长度方向， 载荷步数与长度方向的单元数相同， 热分析采 用默认 的 求 解 器 、默 认 的 求 解 精 度 ， 用 两 种 算 法 对 ２５ 钢 的 双 丝 焊 接 进 行 了 仿 真 实 验 （材 料 性 能 的 非 线 性 参 数 见 参 考 文 献 ［４］）。 记
Ｋｅｙ ｗｏｒ ｄｓ： Ｗｅｌｄｉｎｇ ｓｉｍｕｌａｔｉｏｎ； Ｇａｕｓｓ ｈｅａｔ ｓｏｕｒ ｃｅ； Ｌｏａｄｉｎｇ ａｒ ｉｔｈｍｅｔｉｃ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
中图分类号： ＴＨ１２ ＴＰ３１１ 文献标识码： Ａ
第２期 ２００７ 年 ２ 月
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
高斯热源3D打印APDL
/clear
jbc=2.5e-3/2 Leabharlann jbk=1.6e-3/2
jbh=0.4e-3
fmc=0.9e-3
fmk=0.6e-3
fmh=0.12e-3
v=0.12
p=195
r=0.05e-3
sh=0.04e-3
sc=0.025e-3
sk=0.025e-3
!************************参数计算莵E
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, 21, 17, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -1, 4, -3
TUNIF,298,
NSUBST,3,0,0
ESEL,S,MAT,,2
EKILL,all
alls
pred,off
alls
solve
QMX=2*0.72*P/3.1415926/R/R*1800
ZDZB=R
/SOLU
*DEL,_FNCNAME
*DEL,_FNCMTID
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,sh, , , , , , ,0
FLST,5,4,4,ORDE,4
FITEM,5,2
FITEM,5,4
FITEM,5,-5
FITEM,5,7
ansys高斯热源施加[整理版]
![ansys高斯热源施加[整理版]](https://img.taocdn.com/s3/m/cb5124d259f5f61fb7360b4c2e3f5727a5e92480.png)
!平板对接开v型坡口多道焊成型(温度场)!采用guass表面热源finish/clear,nostart/filname,1.11thermal_analysis/title,3D_weld_thermal_analysis/units,si !定义国际制单位/prep7!**************************************************!定义焊接尺寸及焊接参数!**************************************************!******************************************构件尺寸len=100e-3 !长0.1mwid=50e-3 !宽0.05mthic=6e-3 !厚0.006m!b=0.002 !根部预留(m)!h=0 !钝边(m)wweld=15e-3 !焊接热影响区宽度lnum=10 !沿焊缝长度方向划分份数!******************************************焊接参数I=200 !焊接电流(A)U=30 !焊接电压(V)v=0.5/60 !焊接速度(m/s)arf=60 !坡口角度（°）pi=3.1415926 !pi为圆周率eta=0.7 !焊接热效率,手工电弧焊0.7,埋弧焊0.8K=1.2e4 !热流集中程度系数(/m^2) ?them0=25 !环境温度及母材初始温度,取室温25℃!**********************************************定义高斯表面移动热源参数!GAUSS热源模型qr=qmax*exp(-k*r*r) !(单位W/m2)!假设热源中心坐标为(a,b,c),a=0,b=每个焊缝的最高点的坐标y,c=v*dt!每个焊缝上表面任一点的坐标为(x,y,z) ,x>0, y>0, z>0!r^2=x^2+(y-b)^2+(z-c)^2，每个焊缝上表面任一点距热源中心距离qmax=0.24*k*eta*I*u/pi !加热斑点中心热流(W/m2) percent=0.99radius=abs(sqrt(1/k*log(1/(1-percent)))) !有效热半径，单位为m!********************************************************************** !定义单元类型及材料属性!********************************************************************** et,1,solid70 !计算单元类型!nothing !定义实常数mptemp,1,0,100,200,300,400,500 !定义材料属性温度范围mptemp,,1200,1250,1350,1400,1450,1500mptemp,,1505,1705,1905,2105,2305,2500mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 !定义材料密度(kg/m3) mpdata,dens,1,,7800,7800,7800,7800,7800,7800mpdata,dens,1,,7800,7800,7800,7800,7800,7800mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 !热传导系数(W/(m*K)) mpdata,kxx,1,,30,30,30,30,30,30mpdata,kxx,1,,30,30,30,30,30,30mpdata,c,1,1,450,469,481,508.5,536,569 !考虑相变潜热的比热容(J/kg*K) mpdata,c,1,,700,2172,5116,6589,8061,9533mpdata,c,1,,9533,7757,5982,4206,2431,700!****************************************************************!建立有限元计算模型!*****************************************************************afun,deg !如果涉及到角度，将使用“度”(系统默认是弧度)K,1,0,0,0K,2,wweld,0,0K,3,wid,0,0K,4,wid,thic,0K,5,thic*tan(arf/2)+wweld,thic,0K,6,thic*tan(arf/2),thic,0K,7,0,0,lenA,1,3,4,6A,1,2,5,6cyl4,0,0,thic/cos(arf/2),arf,0,90cyl4,0,0,thic/cos(arf/2)*sqrt(2)/2,arf,0,90L,1,7 !形成拖拉路径L14 Vdrag,all,,,,,,14 !沿线拖拉面Vptn,all !分割体Vglue,all !粘贴体？Nummrg,all!*************************************************************!划分网格!*************************************************************aadd,24,32 !面相加，以便生成映射网格lccat,40,49lccat,42,52!*************************************沿焊缝z方向划分，均为2mm lsel,s,line,,17,21,2lsel,a,line,,24,26,2lsel,a,line,,35,36lsel,a,line,,14,30,16lesize,all,2e-3!***************************************沿厚度y方向划分1-1.5mm lsel,s,line,,6,25,19lsel,a,line,,2,18,16lesize,all,,,6lsel,s,line,,41,43,2lsel,a,line,,50,53,3lesize,all,1.5e-3!***********************************************沿x方向划分1mm lsel,s,line,,45,47,2 !母材lsel,a,line,,44,46,2lesize,all,,,18,0.8lsel,s,line,,7,27,20 !母材lesize,all,,,12,0.7lsel,a,line,,5,23,18lesize,all,,,12,10/7lsel,s,line,,11,34,23 !焊缝lsel,a,line,,48,51,3lesize,all,,,4!****************************************划分体，生成六面体单元allsmshkey,1 !映射网格mshape,0,3d !六面体形状vsel,allvatt,,,1aslv,saatt,,,1type,1 !单元类型为1vmesh,all!*******************************************存档备份有限元模型allssave,'1.11thermal_analysis_meshed','db'finish!================================================================= !重新载入有限元模型进行加载求解!================================================================= /solu/view,,-3,2,-5 !调整窗口视角!*********************************************************************** ! 杀死焊缝区单元!***********************************************************************v1=5 !焊缝1所在的体v2=7 !焊缝2所在的体a1=18 !焊缝1的上表面a2=31 !焊缝2的上表面esel,s,type,,1 !选取1类单元*get,nemax,elem,,num,max !得到所有1类单元号码的最大值*get,nemin,elem,,num,min!********************************************分段杀死第2道焊缝vsel,s,,,v2eslvekill,allesel,s,liveeplot!***************************************************************!边界条件(暂时没有考虑对流)!***************************************************************nsel,s,loc,x,0 !在x=0处施加对称边界条件dsym,symm,x!*****************************************************************!设置非线性求解选项!*****************************************************************allsantype,trans !瞬态分析trnopt,full !瞬态分析选项，完全分析nropt,full,,on !定义完全牛顿--拉普森方法，激活自适应下降（默认）pred,on !打开预测校正timint,on !时间积分设置tintp,0.005,,,1,0.5,0.2 !定义瞬态综合参数tref,25 !参考温度25摄氏度!***************************************************************!设置载荷步参数!***************************************************************t=0 !求解时间初始值dt0=1e-6 !建立初始条件的小时间段dt1=0.5 !起始加热点对应的时间tinc=len/(lnum*v) !载荷步时间间隔,T =len/v,tinc =T/lnum!tsub=0.2 !求解温度场时子步时间!***************************************************************!稳态分析确定初始温度场!***************************************************************!初始条件time,dt0 !loadstep=1,time=dt0timint,off !时间积分选项，关闭（稳态分析）kbc,1 !阶越加载!nsubst,1 !子步为1步esel,s,type,,1 !选择类型1的单元节点nsle,sic,all,temp,25 !初始状态，温度25摄氏度allsoutres,all,allsolvesave,'1.11thermal_analysis_initial_temp','db'!查看初始温度场，应该全为25摄氏度!***************************************************************!按顺序模拟温度场，以第一道为例!***************************************************************!**************************************************第1道焊缝填充esel,s,type,,1 !在第一类单元中加热流载荷esel,r,liveeplot*do,im,0,lnum,1!****************************选择有效半径内节点为当前节点vsel,s,,,v1 !选择焊缝1上表面单元eslv,rnsleasel,s,,,a1nsla,r,1esln,r,0!热源加载, 将各段后点（0-lnum）为热源中心,先加载后删除热源!即先消除上段所加高斯热源,并将上段的温度值作为下段的初始值!*******************************************包括起始加热点c=V*tinc*im !热源位置b=thic/cos(arf/2)*sqrt(2)/2tm=dt1+im*tinctime,tmantype,4,resttimint,onautots,onkbc,1deltim,0.01,0.01,0.2 !时间步长!***********************************以下为施加热流密度载荷*do,i,nemin,nemax,1*if,esel(i),eq,1,then !如果该单元在上面的集合里，则xsy=centrx(i) !读取该单元的中心坐标ysy=centry(i)zsy=centrz(i)rr=abs(sqrt(xsy*xsy+(ysy-b)*(ysy-b)+(zsy-c)*(zsy-c))) !该单元中心距离热源中心*if,rr,le,radius,then !该单元中心在加热半径范围内qr=qmax*exp(-k*rr*rr) !该单元中心处的热流大小sfe,i,nmface(i),hflux,,qr !在每个单元所指定的面上施加热源载荷*endif*endif*enddoallsel !全选求解outres,all,allsolveesel,s,type,,1esel,r,livevsel,s,,,v1 !选择焊缝1上表面单元eslv,rnsleasel,s,,,a1nsla,r,1esln,r,0输入先暂时停止一下!插入"/post1$plnsol,temp$/psf,hflux,,1,1,on$/replot"的内容进行查看!问题：热流载荷的分布和温度场的分布不一致（该完全一致才对），有点偏移？!若要继续，读入"finish$/solu"*do,i,1,6sfedele,all,i,hflux !删除每个单元六个面上的热载荷*enddo*enddoallsesel,s,liveeplot !查看温度场save,'1.11thermal_analysis_alive1','db'!=================================================================!生成节点温度分布云图动画文件的程序段!=================================================================/post1/seg,dele !允许图形数据存储到局部终端存储器/dscale,1,1.0 !显示位移时，设置放大系数avprin,0,0 !规定主项量和总向量如何被计算avres,1 !规定当powergraphies可用时，结果数据如何被平均/seg,multi,1.11temp-10,0.2 !存动画,分节存储后来的显示,动画名字,每个画面间延迟系数esel,s,live!*************************************************读取初始稳态温度场结果set,,,1,,dt0 !从结果文件读,载荷步,子步,比例因子实部,比例因子虚部,时间plnsol,temp!*****************************************读取第1道焊缝过程中的温度场结果*do,im,0,lnum,1tm=dt1+im*tincset,,,1,,tmplnsol,temp*enddo/seg,off,1.11temp-10,0.2anim,1,1 !*********************************生成动画完毕!****************************************存储第1道焊缝过程中的温度场结果图片esel,s,live*do,im,0,lnum,1*if,im,eq,0,thentm=dt1set,,,1,,tmplnsol,temp/image,save,1.11temp0,jpeg*endiftm=dt1+im*tincset,,,1,,tmplnsol,temp*if,tm,eq,dt1+lnum/2*tinc,then set,,,1,,tmplnsol,temp6.5,jpeg*endif*enddo/image,save,1.11temp12.5,jpeg。
ansys 双椭球热源焊接APDL
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。
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)稳态热分析确定在稳态的条件下的温度分布及其他热特性，稳态条件指热量随时间的变化可以忽略。
APDL-ansys
本文列出了A开头的ansys命令及功能。
学习ansys一定要学学APDL，初学者很多人都疑问，那么多指令怎么去背啊其实很简单，常用的也就那么几个，用多了就知道了，忘记查查就行，更多的时候我们并不是自己去写，而是通过ANSYS产生的日志去稍微修改下就可以了。
但是编比较高级的代码流，还要对fortran有所了解，下面是部分ansys指令可供参考。
1. A，P1，P2，…，P17，P18(以点定义面)2. AADD，NA1，NA2，…NA8，NA9(面相加)3. AATT，MAT，REAL，TYPE，ESYS，SECN(指定面的单元属性)【注】ESYS为坐标系统号、SECN为截面类型号。
4. *ABBR，Abbr，String(定义一个缩略词)5. ABBRES，Lab，Fname，Ext(从文件中读取缩略词)6. ABBSAVE，Lab，Fname，Ext(将当前定义的缩略词写入文件)7. ABS，IR，IA，--，--，Name，--，--，FACTA(取绝对值)【注】*************8. ACCAT，NA1，NA2(连接面)9. ACEL，ACEX，ACEY，ACEZ(定义结构的线性加速度)10. ACLEAR，NA1，NA2，NINC(清除面单元网格)11. ADAMS，NMODES，KSTRESS，KSHELL【注】*************12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC【注】*************13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC(变量加运算)14. ADELE，NA1，NA2，NINC，KSWP(删除面)【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6(将既有线沿一定路径拖拉成面)16. AESIZE，ANUM，SIZE(指定面上划分单元大小)17. AFILLT，NA1，NA1，RAD(两面之间生成倒角面)18. AFSURF，SAREA，TLINE(在既有面单元上生成重叠的表面单元)19. *AFUN, Lab(指定参数表达式中角度单位)20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE(复制面)21. AGLUE，NA1，NA2，…，NA8，NA9(面间相互粘接)22. AINA，NA1，NA2，…，NA8，NA9(被选面的交集)23. AINP，NA1，NA2，…，NA8，NA9(面集两两相交)24. AINV，NA，NV(面体相交)25. AL，L1，L2，…，L9，L10(以线定义面)26. ALIST，NA1，NA2，NINC，Lab(列表显示面的信息)【注】Lab=HPT时，显示面上硬点信息，默认为空。
ANSYS-APDL命令汇总
ANSYS-APDL命令汇总ANSYS APDL命令汇总AA,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。
点要依次序输入，输入的顺序会决定面的法线方向。
如果超过四个点，则这些点必须在同一个平面上。
Menu Paths:MainMenu>Preprocessor>Create>Arbitrary>Throug h 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为半径。
基于ANSYSAPDL的堆焊热过程数值模拟
机械 ２００５年第９期 总第３２卷 计算机应用技术 ・４５・———————————————收稿日期：2005－06－12 作者简介：张光明（1972－），男，讲师，工程师，硕士，研究方向为计算机在材料成型中的应用。
基于ANSYS APDL 的堆焊热过程数值模拟张光明，李志宏，彭显平(四川工程职业技术学院 机电系，四川 德阳 618000)摘要：结合课题研究需要，本文建立了适当的有限元计算模型，利用ANSYS 软件的APDL 语言编写了程序，解决了呈高斯函数分布且移动的焊接热源的施加，实现了焊接温度场的三维动态模拟。
关键词：温度场 ；数值模拟；ANSYS ；APDL中图分类号：TG4 文献标识码：A 文章编号：1006－0316(2005)09－0045－03Numerical simulation of temperature field for build-up welding based on Ansys APDLZHANG Guang-ming ，LI Zhi-hong ，PENG Xian-ping(Sichuan Engineering Technical College ，Deyang 618000，China)Abstract ：For the need of researches, the proper finite element model was built. Based on ANSYS APDL developing platform, the program was designed to solve the GAUSS function heat loads applying and obtain the three dimensional dynamical temperature field for build-up welding.Key words ：temperature field ；numerical simulation ；ANSYS ；APDL焊接温度场的准确计算是焊接冶金分析、焊接应力和变形分析以及焊接质量控制的前提。
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参考
不准/UNITS,SI/CONFIG,NRES,100/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.U=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,,,,1 lsel,s,line,,9,10lsel,a,line,,12lsel,a,line,,38 lesize,all,0.005,,,0.2,,,,1 lsel,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,,,,1 lsel,s,line,,23,24 lsel,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,,,,1 vmesh,all/PNUM,DEFA EPLOTfini/soluantype,trans,new!nlgeom,on timint,0,struct timint,1,therm timint,0,mag timint,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.149310, 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.002 nsel,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,liveeplot OUTRES,ALL,ALL, /PSF,HFLUX,,2 /REPLOT OUTRES,ALL,ALL,。
(整理)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模拟不同参数对20钢焊接温度场的影响
温度 /℃
2500 2250 2000 1750 1500 1250 1000 750 500 2500 4
(c) 10 mm
8 12 16 20 24 28 32 36 40 时间 /s
由图 7 可知，随着板材厚度的增加，焊缝处最上 来的影响。
面两个节点的温度增加， 最下面几个节点的温度降
(2) 电流和电压越大， 焊接过程所能达到的温
122
Hot Working Technology 2011 , Vol.40 , No. 15
上半月出版
Casting·Forging·Welding 金属铸锻焊技术
焊缝对称，采用物理模型的一半进行有限元分析，有 限元模型如图 2 所示。
图 2 有限元模型 Fig.2 Finite element model
部的温度变化相反，发生温度相反的界面处，也是缺 陷最容易出现的地方，这也与事实相符。
3 结论
(1) 通过对三种因素进行分析对比， 电流和电 压对焊缝的影响最大， 板材厚度其次， 板材长度最 小。 对于厚度和长度不同的板材， 只要选取的电流 和电压不同， 就会避免板材厚度和长度对焊缝所带
2.2 板长对温度场的影响 电压、电流和板厚相同，焊接开始、焊接结束和
冷却时间相同。 电压 27 V，电流 160 A，板材尺寸为 120 mm ×(40、60) mm × 6 mm。 在沿着焊缝的方向上， 取融合线上等距离的六个点进行温度场模拟， 温度
《热加工工艺》 2011 年第 40 卷第 15 期
值模拟分析与研究; 电话:15270821292; E-mail:linger_6102@
低，强度适中，塑性好，具有良好的焊接性能。故本文 拟 采 用 ANSYS 软 件 对 20 钢 的 焊 接 温 度 场 进 行 模 拟，分析讨论不同焊接参数对其温度场的影响。从而 为实际焊接工艺的设定和优化提供理论指导， 具有 重要的意义。
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: 待耦合的节点号。
《2024年基于ANSYS的焊接温度场和应力的数值模拟研究》范文
《基于ANSYS的焊接温度场和应力的数值模拟研究》篇一一、引言随着制造业和机械工程领域的不断发展，焊接技术已成为生产制造过程中重要的工艺手段。
在焊接过程中，温度场和应力的变化对焊接质量、产品性能及使用寿命具有重要影响。
因此，对焊接过程中的温度场和应力进行准确预测和控制，对于提高产品质量和优化生产过程具有重要意义。
本文基于ANSYS软件，对焊接过程中的温度场和应力进行数值模拟研究，旨在为实际生产过程中的焊接工艺优化提供理论依据。
二、焊接温度场的数值模拟1. 模型建立首先，根据实际焊接过程，建立三维有限元模型。
模型中应考虑焊缝、母材等关键部分的几何形状和材料属性。
同时，为提高计算效率，可对模型进行合理简化。
2. 材料属性及热源模型在模拟过程中，需要输入材料的热导率、比热容、热扩散率等热物理参数。
此外，选择合适的热源模型也是关键。
本文采用高斯热源模型，该模型能较好地描述焊接过程中的热输入分布。
3. 数值求解利用ANSYS软件的热分析模块，对焊接过程中的温度场进行数值求解。
通过设定合理的初始条件和边界条件，求解出焊接过程中的温度分布。
三、焊接应力的数值模拟1. 模型转换在得到温度场分布后，将热分析结果作为应力分析的初始条件。
将热分析模型转换为应力分析模型，并设定相应的材料属性。
2. 应力分析利用ANSYS的应力分析模块，对焊接过程中的应力进行数值模拟。
考虑焊缝收缩、母材约束等因素对应力的影响。
通过求解，得到焊接过程中的应力分布。
四、结果与讨论1. 温度场分析通过数值模拟，可以得到焊接过程中的温度场分布。
分析温度场的变化规律，可以了解焊接过程中的热输入、热传导及热扩散等情况。
同时，还可以预测焊接过程中的潜在问题，如热裂纹、热变形等。
2. 应力分析根据应力分布结果，可以了解焊接过程中产生的残余应力。
残余应力对产品的性能和使用寿命具有重要影响。
通过分析残余应力的分布和大小，可以为优化焊接工艺提供依据。
此外，还可以考虑采用相应的工艺措施，如焊后热处理、优化焊接顺序等，以降低残余应力。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
nsel,a,loc,x,lx
sf,all,conv,10,25
!************定义数组维数**********
MAX_X=1+lx/LSIZE
!*************定义table数组******************************************
antype,trans,new
!nlgeom,on
timint,0,struct
timint,1,therm
timint,0,mag
timint,0,elect
tref,25
nropt,auto
autos,on
kbc,0
pred,on
lnsrch,on
!************杀死焊缝去单元***********
lsel,a,line,,35
lesize,all,0.005,,,0.2,,,,1
lsel,s,line,,9,10
lsel,a,line,,12
lsel,a,line,,38
lesize,all,0.005,,,0.2,,,,1
lsel,s,line,,2
lsel,a,line,,4,5
%_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
V=0.008 !速度
pi=3.1415926
U=25
I=180
Q=U*I !电源有效功率
R=0.006 !电源有效热半径
effect=0.8
Qmax=effect*Q/(pi*R*R) !中心处最大热流
tinc=LSIZE/V !每小段的时间间隔
lx=0.1
lz=0.1
ly=0.006
!*************定义焊接参数***********
weld=0.01 !每道焊道宽度
wheight=0.004 !焊道的高度
LSIZE=0.005
OUTRES,ALL,ALL,
*ENDDO
%_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)
lsel,a,line,,7
lsel,a,line,,14
lsel,a,line,,16
lsel,a,line,,19
lsel,a,line,,39
lesize,all,0.002,,,1,,,,1
lsel,s,line,,23,24
lsel,a,line,,41,42
lsel,a,line,,45
nsel,s,loc,y,ly,ly+wheight
esln,s,1
ekill,all
allsel,all
esel,s,live
eplot
!**********施加对流载荷*************
esel,s,live
nsel,s,loc,z,0
nsel,a,loc,z,lz
/UNITS,SI
/CONFIG,NRES,10000
/PREP7
ET,1,SOLID70
MP,DENS,1,7930
MP,C,1,502
mptemp,1,20,100,500
mpdata,kxx,1,1,12.1,16.3,21.4
/VIEW,1,1,1,1
!*************定义基板尺寸***********
!-->
!**********激活单元*********
esel,s,live
nsel,s,loc,x,(i-1)*lsize-0.002,(i-1)*lsize+0.002
nsel,r,loc,z,-weld/2,weld/2
nsel,r,loc,y,0,wheight
esln,s,0
ealive,all
allsel,all
!**********施加热流载荷**********
esel,s,live
nsel,s,loc,Y,wheight
sf,all,hflux,%heatflux%
allsel,all
time,i*tinc
nsubst,2
!定义表格各行
*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
!*****************建立模型，生成网格**************
block,0,lx,0,ly,0,lz
wpoff,0,ly,lz/2
block,0,lx,0,wheight,-weld/2,weld/2
wpoff,0,0,-weld/2
vsbw,all,,delete
%_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,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
*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,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
*do,i,1,MAX_X,1
local,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'
%_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,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,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
lsel,a,line,,47
lsel,s,line,,23,24
lsel,a,line,,41,42
lsel,a,line,,45
lsel,a,line,,47
lesize,all,0.002,,,1,,,,1
vmesh,all
/PNUM,DEFA