【讨论】激光加热微流体温度场分析

feaxiang

我做了一个激光加热微流体热分析，激光模型是J=0.92*P*exp(-((x-v*t)^2+y^2)/r^2),其中，P，功率，r聚焦激光半径，分析类型是瞬态热分析，
模型尺寸取为2r*1.5r*1.5r,在xy平面网格取均匀网格h=1e-5,在z向变网格
我的有限元程序
  
/title,transient temperature analysis
/prep7
/units,si
et,1,70 !model
!material parameters
mp,c,1,1420
mp,dens,1,1190
mp,kxx,1,0.19
!real number
tmax=1617!降解温度
!time=1.2e-3
r=1.2e-4
h=1e-5
length=2*r
width=1.5*r
  
depth=1.5*r
!model solid
block,0,length,0,width,0,depth
  
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2   
AESIZE,P51X,h,   
LPLOT   
FLST,5,4,4,ORDE,2   
FITEM,5,9   
FITEM,5,-12  
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
LESIZE,_Y1,0.32*h, ,,1.2, , , ,1   
vsel,,,,1
vmesh,all
v=0.2
tau=h/v!time step
numtime=2*r/h+1
qmax=0.92*15/(3.1416*r**2)
ii=0
nsel,s,loc,z,0
esln
*get,nse,elem,0,count
*get,minelem,elem,0,num,min
*get,maxelem,elem,0,num,max
*dim,neorder,array,nse!define an array to receive hflux
!elem order
*do,i,minelem,maxelem
  *if,esel(i),eq,1,then
    ii=ii+1
    neorder(ii)=i!单元排序
  *endif
*enddo
allsel,all
finish
save
/solu
antype,4
TRNOPT,FULL  
nropt,full
ESTIF,0
kbc, 1
autos,1
tunif,20
tref,273
time,tau
deltim,tau
*do,j,1,nse
  sfe,neorder(j),1,hflux,,qmax*exp(-centrx(neorder(j))**2/r**2-centry(neorder(j))**2/r**2)
*enddo
allsel,all
solve
finish
/post1
/view,1,1,1,1
set,last
!plnsol,temp
!prnsol,temp
etable,eras
etable,ll,temp
!pretab,ll
esel,s,etab,ll,tmax!结果大于tmax的为死掉的单元，也就是被加工的单元
*get,minek,elem,0,num,min
*get,maxek,elem,0,num,max
*do,ei,minek,maxek
  *if,esel(ei),eq,1,then
    *do,orgin,1,nse
       *if,neorder(orgin),eq,ei,then
           neorder(orgin)=eladj(ei,6)
       *endif
    *enddo
  *endif
*enddo
  
finish
!the following step load
*do,tm,2,numtime
  
time，tm*tau
/solu
  
ekill,all
esel,s,live
!eplot
tm=2
*do,j,1,nse
  sfedele,neorder(j),1,hflux
  sfe,neorder(j),1,hflux,,qmax*exp(-(centrx(neorder(j))-(tm-1)*h)**2/r**2-centry(neorder(j))**2/r**2)
*enddo  
deltim,tau
!作用时间
!allsel,all
solve
  
/post1
set,last
allsel,all
plnsol,temp
!prnsol,temp
etable,eras
etable,ll,temp
!pretab,ll
esel,r,etab,ll,tmax
*get,minek,elem,0,num,min
*get,maxek,elem,0,num,max
*do,ei,minek,maxek
  *if,esel(ei),eq,1,then
    *do,orgin,1,nse
       *if,neorder(orgin),eq,ei,then
           neorder(orgin)=eladj(ei,6)
       *endif
    *enddo
  *endif
*enddo
  
*enddo
/solu
ekill,all
esel,s,live
eplot
/post1
plnsol,temp
etable,eras
etable,ll,temp
pretab,ll
finish
  
我要求模拟加工过程，也就是当单元温度达到或超过tmax＝1617时，单元就死掉，这样通过循环就得到了流到形状，但是我在time中遇到问题，time是载荷作用时间还是最终时间，我做了发现结果与要求的差别很大，实际最大加工深度为113um，不知问题出在哪里，另外自适应网格怎么做，有没有例子，有没有根据结果来调整网格大小的例子，谢谢，我的email: hfxiang@zju.edu.cn