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,,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,,,,1vmesh,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%!!Beginofequation: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!Endofequation: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