求发一个移动高斯热源的ansys程序

icesblue

刚接触到高斯热源,哪位高手能发一个焊接移动高斯热源的程序让我学习一下,感激不尽!

qiqik0521

/filname,maichong3,0
/title,the thick of the board is 3mm,velocity is 0.5m/min
/prep7  
et,1,solid70
keyopt,1,2,0
keyopt,1,4,0
keyopt,1,7,0
keyopt,1,8,0
!*  
!99.75%fe
mptemp,1,0,20,200,500,800,1100,1500
mpdata,dens,1,1,7880,7880,7710,7588,7320,7278
mpdata,kxx,1,1,53.18,47.74,36.01,32,29.3,25
mpdata,c,1,1,470,470,470,935,1879,1879
mpdata,enth,1,1,7.17e8,14.35e8,28.7e8,43.7e8,52.95e8,82.15e8
block,-0.025,0.025,-0.003,0,0,0.3,
block,-0.005,0.005,-0.003,0,0,0.3,
vsel,s,volu,,1,2
vovlap,all
vglue,all
numcmp,line
numcmp,area
numcmp,volu
/replot
aesize,all,0.001,
lesize,all,0.001, , , ,1, , ,1,
esize,0.001,0,
mshkey,1
mshape,0,3d
vmesh,1
esize,0.01,0,   
aesize,all,0.01,   
lesize,all,0.01, , , ,1, , ,1,
mshkey,1
mshape,0,3d
vmesh,2,3

/view, 1 ,1,1,1
/ang, 1
/rep,fast   
/psf,hflu,,2,0
/pnum,line,1
finish
/solu
antype,trans
tunif,20
timint,on !time integration = on
sfa,1,,conv,50,25,!w/m**2*c
sfa,2,,conv,50,25,
sfa,12,,conv,50,25,
sfa,16,,conv,50,25,
sfa,11,,conv,12.5,25,!bottom
sfa,5,,conv,12.5,25,
sfa,15,,conv,12.5,25,
allsel,all
autots,on
timint,on
kbc,1 !step
time,2
deltim,0.02,0.02,0.1
outres,all,4
nsel,s,loc,x,-0.003,0.003
nsel,r,loc,y,0
nsel,r,loc,z,0,0.3
sf,all,hflu,%feng%
!其中%feng%为
200*24*0.6/{pi}/{pi}/(0.003^2)*exp(-1e6*100*({x}^2+({z}-8.3e-3* {time})^2))
allsel,all
solve

tnanr

上面的程序好用吗？