【求助】有谁做过移动热源的切割或焊接过程温度场分析？

feaxiang

也给我一份好吗，我的信箱hfxiang@zju.edu.cn
下面是我做的
finish
/clear
/title,transient temperature analysis
/prep7
et,1,70 !model
!material parameters
mp,c,1,1.42e+009
!mp,c,1,6563!等效比热法，就是将降解潜热等效到比热中去
mp,dens,1,1.19e-6
mp,kxx,1,190
!real number
tmax=1617!降解温度
!tg=370!降解温度
r=1.2e-1
h=0.5e-2
length=2*r
width=1.2*h
depth=1.2*r
!model solid
block,0,length,0,width,0,depth
!net size
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.3*h, ,,1.28, , , ,1   
vsel,,,,1
vmesh,all
!keyopt,1,8,1
v=200
tau=h/v!time step
numtime=2*r/h
qmax=0.92*15e+006/(3.1416*r**2)
nsel,s,loc,z,0
esln
!select element and sort element define an array neorder
*get,nse,elem,0,count
!*get,minelem,node,0,num,min
!*get,maxelem,node,0,num,max
*dim,neorder,array,nse
!elem order
*do,i,1,nse
  neorder(i)=i
*enddo
allsel,all
finish
save
  
!calculate
/solu
antype,4
nropt,full
time,tau
nsubst,2
kbc, 0
/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST
tunif,20
tref,273
timint,on
autos,1
TSRES,ERASE
!first time step
*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
/post1
plnsol,temp
etable,ll,temp
pretab,ll
fini
/solu
*do,j1,1,nse
    sfedele,neorder(j1),1,hflux
*enddo
!d,all,temp,0
esel,s,etab,ll,tmax
*get,minek,elem,0,num,min
*get,maxek,elem,0,num,max
!sort ekill element,ekilled element is replaced
*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
ekill,all
  esel,s,live
  eplot
!循环加载
*do,tm,2,numtime
  
time,tm*tau
nropt,full
kbc,0
autos,1
timint,on
TSRES,ERASE
*do,j,1,nse
    sfe,neorder(j),1,hflux,,qmax*exp(-(centrx(neorder(j))-(tm-1)*h)**2/r**2-centry(neorder(j))**2/r**2)
*enddo
allsel,all
solve
fini
/post1
plnsol,temp
etable,ll,temp
!pretab,ll
fini
/solu
*do,j1,1,nse
    sfedele,neorder(j1),1,hflux
*enddo
  
!d,all,temp,0
esel,r,etab,ll,tmax
*get,minek,elem,0,num,min
*get,maxek,elem,0,num,max
!sort ekill element,ekilled element is replaced
*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
ekill,all
esel,s,live
eplot
*enddo
  
但是问题出现了，就是每次加载荷是时，作用时间不是从前一次作用时间开始，而是从0开始，这是什么原因

zyp

能给我发一份吗
zyp_820@163.com
谢谢了

Joanna

我也要一份好吗？joanajia◎sohu.com

yltechno

我有全套phoenics上做的程序，不知各位有无用途！！！

sl3913

太好了，能来一份吗！
请发到slcumt@126.com，谢谢！

wpengw123

能给我一份么？
我的信箱为guyue1103@hotmail.com
谢谢你了

wangzhonghui

请给我一份吧。谢谢，chendlut!
   email:wangzhonghui@bipt.edu.cn

zhtbuaa

我也要一份：zhtbuaa@63.com

ncujjp

o偶也要一份。jpjia@yahoo.com.cn        谢谢

ncujjp

我也要.  jp_jia@sohu.com

ncujjp

wo ye yao .         jpjia@126.com         谢谢

Nancy1101

能不能给我一份啊，我的email：wtt_135@emails.bjpu.edu.cn
谢谢！

wangzhonghui

能否发给我一份？谢谢！
email:wangzhonghui@bipt.edu.cn

Nancy1101

我没有收到，麻烦楼主再给我发一份好吗？
e－mail：wtt_135@emails.bjpu.edu.cn
谢谢！

niathjing

给我也发一份吧.我需要这方面的资料.谢谢.
99111wf3@163.com