求解高斯热源加载！

ikuangke

我做了一个平板焊接温度场模拟名令流，但运行完之后，所加载的高斯热源却成菱形分布，如图所示：


  
！
正常的话应该是近似于圆形分布，请高手帮我检查下，到底是什么原因，谢谢了！
以下是我做的命令流：
/clear
/filname,seal
/title,temp seal
/units,si
/prep7
lx=0.06
ly=0.03
lz=0.01
lsize=0.002
pi=3.1415926
r=0.002
v=0.005
p=1200
absor=0.4
tm=lx/v
tinc=lsize/v
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=ly/lsize
qm=2*absor*p/(pi*r**2)
et,1,solid70
keyopt,1,2,1
et,2,surf152
r,2
mp,dens,1,7840
mptemp,1,20,100,300,500,700,900,1100
mpdata,kxx,1,1,47,43,38,34,28,25,24
mpdata,c,1,1,472,480,524,615,854,637,595
/view,1,1,1,1
block,0,lx,0,ly,0,lz
lesize,1,,,ly/lsize,0.2
lesize,6,,,ly/lsize,5
lesize,3,,,ly/lsize,5
lesize,8,,,ly/lsize,0.2
lesize,2,,,lx/lsize
lesize,4,,,lx/lsize
lesize,5,,,lx/lsize
lesize,7,,,lx/lsize
lesize,9,,,lz/lsize
lesize,10,,,lz/lsize
lesize,11,,,lz/lsize
lesize,12,,,lz/lsize
type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,1
type,2
asel,s,,,1,6
amesh,all
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*lsize
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qm*exp(-2*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
/solu
antype,trans
trnopt,full
toffst,273
tunif,20
nsel,s,loc,x,0
dsym,symm,x
kbc,1
timint,on
autots,off
esel,s,type,,2
sf,all,conv,65,20
*do,t,1e-6,tm,tinc
asel,s,area,,2
nsla,s,1
sf,all,hflux,%flux2%
allsel,all
time,t
autots,off
deltim,,,,1
kbc,0
outres,all,all
tsres,erase
solve
*enddo
finish

拜托各位大侠了！

ikuangke

高手呢？高手都藏哪了？  出来吧，帮帮忙啊！

huhui

你的高斯热源公式再看看是不是有问题，我印象中的高斯热源e的前面应该是有关半径变量的

ycs2011

你这个没有个注释，我们这种新手很难给你意见啊，看不懂啊！

天晴

我的也是这种分布状态，不知道这个问题你解决了没？？

Vicly2012

新手路过，表示看不懂。求热源加载例子，多谢楼主啦！