急！激光移动热源怎么总是出现偏离？

418884924

finish
/clear,all
/filname,weldtem
/title,welding
/units,si
/prep7
et,1,solid70
et,2,solid70
keyopt,1,2,1
et,3,surf152
r,3
keyopt,3,4,1
keyopt,3,5,0
keyopt,3,8,3
mptemp,1,0,400,600,800,1000,1500
mpdata,kxx,1,1,18.9,28.35,31.08,32.76,31.5,28.2
mpdata,c,1,1,580,756,840,882,756,797
mpdata,hf,1,1,9.66,59.43,84.32,109.2,180.6,359
mpdata,dens,1,1,7550,7550,7550,7550,7550,7550
pi=3.1415926
i=160
u=23
efficient=0.78
r=0.0045
qm=3*i*u*efficient/(pi*r**2)
kk=0.85
qmax=qm*kk
ly=0.12
lx=0.4
lz=0.012
lsize=0.01
v=0.00185
tm=lx/v
tinc=lsize/v
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=1+ly/lsize
/view,1,1,1,1
hy=0.5*ly
block,0,lx,-hy,hy,0,-lz
block,0,lx,-hy,hy,0,lz
VGLUE,ALL
NUMCMP,ALL
save
esize,lsize
vsel,s,volu,,1
mshkey,1
vatt,1,1,1
vmesh,all
ALLSEL,ALL
vsel,s,volu,,2
mshkey,1
vatt,1,1,2
vmesh,all
ALLSEL,ALL
type,3            !指定单元类型
mat,1             !指定材料类型
real,             !指定实常数类型
esys,0            !指定坐标系
secnum,            !指定梁号
tshape,line        !指定接触目标面为线接触
nsel,all            !选择所有节点
esurf,0    !？？？在体表面生成表面效应单元
allsel,all  !选择所有
/pnum,line,off
/pnum,defa
eplot
finish
*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)=qmax*exp(-3*kk*distance**2/(r**2))
     *else
       flux2(i,j,k)=0
     *endif
     *enddo
   *enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
/solu
antype,trans
tunif,15
kbc,1
timint,on
autots,off   
esel,s,type,,3
nsle,s
sf,all,conv,-1
allsel
*do,t,1e-6,tm,tinc
time,t
autots,0
deltim,,,,1
kbc,0
outres,all,all
tsres,erase
sf,all,hflux,%flux2%
solve
*enddo
finish

——————————————————————————
下面是温度场，我的目的是想让热源沿着y=0，x方向移动，可算出来后发生了偏移

418884924

没人回答么，自己顶一下，别沉了

418884924

没有人帮忙么，急用，斑竹有没有时间帮我看看是怎么回事啊

418884924

都沉到底了，还没有人来看我啊

rizen

根据你的图 好像是沿着长轴X移动呀？ 你的偏移是指向Y方向偏了？

iewg

这个看不懂，现在正在学习移动热源的问题

yulu10

嘎嘎嘎嘎嘎嘎嘎