【求助】请教关于3维热分析移动面载的寿命

zlj

最近作料一个程序，金属表面加载hflux匀速直线移动载荷，结果发现在第n load step  时，1，2，3。。。。，n-1  load  step  的 hflux 仍然存在，即所有过程中每一步的载荷，在该载荷步开始的时刻被加上，在该载荷步结束时刻以后并不消失，而是一直持续到热源到达终点的时刻。
在论坛里，wuhanwangtao 曾提及surface loads的寿命有所不同，可惜没有再细讲。 请大家多多指教。

zlj

/prep7  
/UNITS,si
   
et,1,solid70
mp,c,1,754   
mp,kxx,1,30  
mp,dens,1,8500
  
block,0,20,0,40,0,4  
  
esize,1  
vmesh,all  
FINISH  
/solu   
antype,trans
tunif,20
kbc,1
   
*SET,delt,1  
tim=0
  
*do,i,1,40,1 ！！！！！！
tim=tim+0.1   
time,tim
timint,on  
deltim,0.1,0.01,1
autots,on
  
nsel,s,loc,z,2.5,4.5
nsel,r,loc,x,18.5,20.5   
nsel,r,loc,y,i-1.5,i+0.5
esln,s,1
  
sfe,all,6,hflux,,10e8
  
  
allsel
  
outres,all,all  
solve  
save
   
*enddo   
finish  
save

zlj

上述命令运行后最后一步得到下图

zlj

把*do,i,1,40,1
改为*do,i,1,80,1
运行后最后一步得到下图

zlj

进一步修改后，（如下）
  
/prep7  
/UNITS,si
   
et,1,solid70
mp,c,1,754   
mp,kxx,1,30  
mp,dens,1,8500
  
block,0,20,0,40,0,4  
  
esize,1  
vmesh,all  
FINISH  
/solu   
antype,trans
tunif,20
kbc,1
   
*SET,delt,1  
tim=0
  
*do,i,1,40,1
tim=tim+0.1   
time,tim
timint,on  
deltim,0.1,0.01,1
autots,on
  
nsel,s,loc,z,2.5,4.5
nsel,r,loc,x,18.5,20.5   
nsel,r,loc,y,i-1.5,i+0.5
esln,s,1
  
sfe,all,6,hflux,,10e8
  
  
allsel
  
outres,all,all  
solve  
save
！卸掉前一步的hflux
sfcum,hflux
nsel,s,loc,z,2.5,4.5
nsel,r,loc,x,18.5,20.5   
nsel,r,loc,y,i-1.5,i+0.5
esln,s,1
  
sfe,all,6,hflux,,1e-10
  
allsel
outres,all,all  
   
*enddo
  
allsel
outres,all,all   
  
finish
运行后最后一步得到下图

zlj

这个结果就相对合理多了。
我想知道：
1.把*do,i,1,40,1  
改为*do,i,1,80,1  
运行时程序为什麽不提醒出错。
2.对移动的面载有没有其他的处理方法。
请多多指教，十分感谢！

zlj

物理模型示意简图

zlj

指点一下，谢谢！

jinchee

下面這方法你參考一下吧...   細節不用管,
你把前一個 load step 的 heat flux 殺了,
不就好了.
  
==========
  
j=1  
  
*DO,i,200,80,-5
  
  TIME,j   
  DELTIM,0.11,0.11,0.11,  
  
  SFEDELE,i+5,2,HFLUX                          !delete heat flux of previous step
  SFE,i,2,HFLUX, ,5/0.002/0.002, , ,         !apply heat flux, face 2 , 5 W
  solve
  
  j=j+1
  
*ENDDO

zlj

谢谢jinchee,我已经试过了,效果很好。

laokaner

学习学习