移动热源的实现，改进问题

tulongpower

各位大侠，版主
一下是我做的一个移动热源的实现，感觉效果很不明显，而且中间貌似有些错误
望能帮我改进下，发回我邮件tulongpower@qq.com，非常感谢。


/prep7
et,1,70
MP,DENS,1,7800
mptemp,,100,200,300,400,500,600,700,800,1200,1500,2500！定义表格温度值
mpdata,kxx,1,,50.66,48.57,42.28,41.03,37.68,34.33,31.4,27.63,30,30,30！导热系数
mpdata,c,1,,469,481,500,536,569,648,695,700,687,870,700  ！比热
mp,kxx,2,52
mp,c,2,450
mp,dens,2,7800
k,1,0,0,0       !关键点
k,2,24,0,0
k,3,30,0,0
k,4,0,0,24
k,5,0,0,30
k,6,0,0,3
k,7,50,0,0
k,8,50,0,3
lstr,2,3        !直线
lstr,1,7
lstr,7,8
lstr,6,8
lstr,4,5
lstr,1,6
LARC,2,4,1,24,        !弧线
LARC,5,3,1,30
FLST,2,4,4            !面
FITEM,2,7   
FITEM,2,8   
FITEM,2,5   
FITEM,2,1   
AL,P51X
FLST,2,4,4  
FITEM,2,4   
FITEM,2,2   
FITEM,2,6   
FITEM,2,3   
AL,P51X
APLOT   
ASBL,       2,       8  
ADELE,       4, , ,1
wprot,,90,           !园
CYL4,29.84962,3,5   
FLST,3,2,4,ORDE,2   ！切割
FITEM,3,1   
FITEM,3,9   
ASBL,       2,P51X  
ASBL,       5,       8  
ASBL,       6,      11  
FLST,2,3,5,ORDE,3   
FITEM,2,2   
FITEM,2,4   
FITEM,2,7   
ADELE,P51X, , ,1
FLST,2,3,5,ORDE,3           !拉伸
FITEM,2,1   
FITEM,2,3   
FITEM,2,5   
VEXT,P51X, , ,0,100,0,,,,   
vglue,all
EXTO,VSWE,AUTO,OFF        !网格划分?
EXTO,VSWE,AUTO,ON   
FLST,5,3,6,ORDE,2   
FITEM,5,4   
FITEM,5,-6  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
numcmp,all      !压缩编号
/SOLU
AnTYPE,TRANS            ！
定义分析类型
TUNIF,20                    ！
赋予初温20度
kbc,1                         ！定义加载
方式为阶跃载荷
timint,on                       ！打开瞬态
效应
CONV1=110                  ！
定义对流系数
asel,s,area,,8,9,13,16            
！选择散热面
sfa,all,,conv,CONV1,200000          ！
在面上施加对流载荷
allsel                           ！
选择所有
vsel,s,volume,,3              ！
选择编号为3实体
eslv,s,1                        ！
选择实体单元
ekill,all                 
！杀死所有选择单元
allsel                   ！
选择所有
delt=1                     ！
定义步长
*do,i,1,6                  ！
定义循环
time,i*delt                  
！定义变量
vsel,s,volume,,1              
！选择编号为3实体
nslv,s,1                     ！
选择实体单元
nsel,r,loc,y,(i-1)*5-1,(i-1)*0.5+1  ！定义加载温度的区间
esln,s,0
ealive,all                                           ！激活所有选择区间单元
allsel                                              ！选择所有单元
nsel,s,loc,y,(i-1)*5-1,(i-1)*5+1     ！定义加载温度的区间
esln,s,0
esel,r,mat,,1                                         ！选择编号1的材料,esel选择现有的
bfe,all,hgen,1,1.23e+100                                ！热源温度（生热率）,定义单元体载
荷
allsel
*if,i,gt,1,then                                        ！定义循环条件
nsel,s,loc,y,-1,(i-2)*5+1                               ！选择单元区间
esln,s,0                                            ！选择附着与被选中节点的单元
esel,r,mat,,1
BFEDELE,ALL,HGEN                               ！施加冷却约束，删除热生成载荷
allsel
*endif   
deltime,delt                                    ！指定载荷步
autots,on                                             ！打开时间积分
outres,all,all                                          ！定义结果输出方式
solve                                                ！求解
*enddo                                              ！结束循环

wxiaotm

我也正想学习怎样实现热源的移动，谢楼主分享

mochu

我非常想问一个问题啊？我感觉界面操作很简单啊，为什么要写代码来完成呢

cfrfosh

是的模型是很简单的，但只用来调试命令流的，在你做实际模型中，可能是一个很复杂的模型，杀鸡为什么要用牛刀，应为鸡的背后是一个更强大猎物，所以用命令流可以让你迅速调试

wxpand2004

如果楼主已经解决问题的话就麻烦给我也发一份呗！
wxpand2004@163.com
不胜感谢！

wuyp802

能不能 把你想做的模型 写出来
只有程序 不好判断 核心部分  ，只能看见一些硬伤

wuyp802

如果已经改好了  能不能把定稿程序 发一份？谢谢
tianquanrong1986@126.com

wuyp802

如果后面的解释 是你自己加的     也有一些错误

守护蓝天

你好！你解决了吗？可以指点下吗？谢谢！ 5# wxpand2004