下天雄 发表于 2007-9-26 22:37:26

瞬态问题求助

模型很简单
一个简单的细直长方体
先用稳态求解一个温度分布
然后,用这个温度分布作为初始条件,再进行瞬态求解
求解时,下端固定为4.5K,上端有一股热流引进
想求解什么时候,上端温度达到523K

在求解的过程中,老是出现问题
在长方体的中间会出现0.2K的温度
不知道是什么问题造成的?


以下是我的命令流
请大家帮忙看看!
先谢谢了
finish
/clear
heat=15
area=4.3e-3*1.84e-3*2/3
/prep7
et,1,solid70
block,0,4.3e-3,0,1.84e-3,0,450e-3
!Ag/Au HTS current
MPTEMP,,,,,,,,
MPTEMP,1,4,9,14,19,24,29,
MPTEMP,7,34,39,44,49,54,
MPDATA,KXX,1,1,10.095,16.87,23.89,30.992,37.993,44.724
MPDATA,KXX,1,7,51.011,56.68,61.56,65.48,68.26
MPTEMP,,,,,,,,
MPTEMP,1,4.2,10,20,40,60,80
MPTEMP,7,100,120,140,160,180,200
MPTEMP,13,220,240,260,273,293,300
MPTEMP,19,320,340,360,373,400,600,
MPDATA,c,1,1,0.142,2,15.75,62.75,108.5,132
MPDATA,c,1,7,148.5,158.5,164.5,168.5,171.5,174.0
MPDATA,c,1,13,176,178,179.5,180,180.5,180.5
MPDATA,c,1,19,181,181,181,181,181,181
mp,dens,1,15000
lsel,s,length,,1.84e-3
lesize,all,,,3
lsel,s,length,,4.3e-3
lesize,all,,,5
lsel,s,length,,450e-3
lesize,all,,,100
vatt,1
vsweep,all
finish
/solu
antype,static
asel,s,loc,z,0
nsla,s,1
d,all,temp,45
asel,s,loc,z,450e-3
nsla,s,1
d,all,temp,5
allsel
outpr,basic,all
solve
finish
/POST1
SET,LAST
PLNSOL, TEMP,, 0
/PREP7
allsel,all
LSCLEAR,fe      
LSCLEAR,solid   

/solu
ANTYPE,4
ldread,temp,,,,,,rth
outres,,all
solcon,on   
NEQIT,25
LNSRCH,1
timestep=0.5
KBC,1   

!DELTIM,timestep,timestep,timestep
calcno=500
*do,i,1,calcno
/solu
TIME,timestep*i

asel,s,loc,z,450e-3
nsla,s,1
d,all,temp,5

asel,s,loc,z,0
sfa,all,,hflux,heat/area
allsel
solve
/POST1
SET,LAST
PLNSOL, TEMP,, 0
*get,T,node,1004,temp
*if,T,gt,523,then
*exit
*endif
*enddo
/POST1
SET,LAST
PLNSOL, TEMP,, 0

下天雄 发表于 2007-9-27 19:14:10

大家帮帮忙吧
页: [1]
查看完整版本: 瞬态问题求助