tulongpower 发表于 2010-3-8 22:03:39

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

各位大侠，版主
一下是我做的一个移动热源的实现，感觉效果很不明显，而且中间貌似有些错误
望能帮我改进下，发回我邮件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 发表于 2010-3-8 23:09:26

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

mochu 发表于 2010-3-9 15:07:56

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

cfrfosh 发表于 2010-3-9 21:09:00

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

wxpand2004 发表于 2010-3-10 09:39:26

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

wuyp802 发表于 2010-5-13 12:46:58

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

wuyp802 发表于 2010-5-13 12:48:24

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

wuyp802 发表于 2010-5-13 12:56:28

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

守护蓝天 发表于 2010-5-27 17:11:52

你好！你解决了吗？可以指点下吗？谢谢！ 5# wxpand2004
页: [1]
查看完整版本: 移动热源的实现，改进问题