caicai1959 发表于 2009-4-11 20:57:27

求助！请高手帮我看看温度出了什么问题！

请高手看看这是出了什么问题，我明明把初始温度设为了第一层140、下层45度啊，为什么我输入命令后出来的结果，初始0.01s的稳态温度是-40度！想不通啊，麻烦指教一下！
finish
/clear
!热分析
/filename,themal
/title,041102

Ha=0.06               !第一层厚度
k1=1.28+Ha
Tair=20                        !空气温度
Ta=140
Tdown=44.5                !下层温度
kxx1=2
kxx2=1.9
kxx3=1.75
kxx4=1.6
kxx5=1.45
con=16.61

/PREP7               
!定义单元类型
ET,1,SOLID70
et,2,surf152
KEYOPT,2,3,0
KEYOPT,2,4,1   
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,1   
KEYOPT,2,9,0


!设置材料属性
MPTEMP,,60,80,100,120,140 !建立温度表
MPDATA,KXX,1,,kxx1,kxx2,kxx3,kxx4,kxx5
MPDATA,C,1,,920
MPDATA,DENS,1,,2565
MPDATA,KXX,2,,1.5
MPDATA,C,2,,920
MPDATA,DENS,2,,2550
MPDATA,KXX,3,,1.5
MPDATA,C,3,,920
MPDATA,DENS,3,,2540
MPDATA,KXX,4,,2.44
MPDATA,C,4,,860
MPDATA,DENS,4,,2750
MPDATA,KXX,5,,2.2   
MPDATA,C,5,,860
MPDATA,DENS,5,,2650
MPDATA,KXX,6,,1.7   
MPDATA,C,6,,860
MPDATA,DENS,6,,1870
MPDATA,HF,7,,con   !空气对流系数
Save

!创建几何模型
k,1,-2,0,0,,,,   
k,2,2,0,0,,,,
k,3,2,0.5,0,,,,
k,4,-2,0.5,0,,,,   
k,5,-2,0.8,0,,,,
k,6,2,0.8,0,,,,   
k,7,2,1.1,0,,,,
k,8,-2,1.1,0,,,,
k,9,-2,1.19,0,,,,
k,10,2,1.19,0,,,,
k,11,2,1.28,0,,,,
k,12,-2,1.28,0,,,,   
k,13,-2,K1,0,,,,   
k,14,2,K1,0,,,,
FLST,2,4,3
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
A,P51X
FLST,2,4,3
FITEM,2,4   
FITEM,2,3   
FITEM,2,6   
FITEM,2,5   
A,P51X
FLST,2,4,3
FITEM,2,5   
FITEM,2,6   
FITEM,2,7   
FITEM,2,8   
A,P51X
FLST,2,4,3
FITEM,2,8   
FITEM,2,7   
FITEM,2,10
FITEM,2,9   
A,P51X
FLST,2,4,3
FITEM,2,9   
FITEM,2,10   
FITEM,2,11   
FITEM,2,12   
A,P51X
FLST,2,4,3
FITEM,2,12   
FITEM,2,11   
FITEM,2,14   
FITEM,2,13   
A,P51X
FLST,2,6,5,ORDE,2   
FITEM,2,1   
FITEM,2,-6
VEXT,P51X, , ,0,0,29.88,,,,
FLST,2,6,6,ORDE,2   
FITEM,2,1   
FITEM,2,-6
VGLUE,P51X   
Save

!定义各层属性
!*
CM,_Y,VOLU
VSEL, , , ,       6
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , ,       5
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       2, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , ,       4
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       3, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , ,       3
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       4, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , ,       2
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       5, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       6, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*

!sweep划分：长100份、宽20份、厚度3、3、3、6、6、10份
FLST,5,1,4,ORDE,1   
FITEM,5,25
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,100, , , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,22
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,20, , , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,48
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,3, , , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,43
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,3, , , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,38
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,3, , , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,33
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,6, , , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,28
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*
LESIZE,_Y1, , ,6, , , , ,1
!*
FLST,5,1,4,ORDE,1   
FITEM,5,21
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*LESIZE,_Y1, , ,10, , , , ,1
!*
FLST,5,6,6,ORDE,2   
FITEM,5,1   
FITEM,5,-6
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*

asel,s,loc,y,k1   !对流边界条件施加在顶面
type,2         !定义属性
amesh,all
esel,s,type,,2      !选择表面效应单元
sf,all,conv,-7,Tair !确定对流边界条件
asel,s,loc,x,-2
da,all,,temp,Tdown
asel,s,loc,x,2
da,all,,temp,Tdown
asel,s,loc,y,o
da,all,,temp,Tdown
asel,s,loc,z,o
da,all,,temp,Tdown
asel,s,loc,z,29.88
da,all,,temp,Tdown   !两边和底面绝热

！进入加载求解
/solu
Antype,trans!设定分析类型为瞬态分析
TOFFST,273! 设定绝对零度
timint,off    !先作稳态分析,确定初始条件
time,0.01   !设定只有一个子步的时间很小的载荷步
deltim,0.01
nsel,s,loc,y,1.28,k1
IC,all,TEMP,Ta, ！施加初始第一层温度
nsel,s,loc,y,0,1.28
IC,all,TEMP,Tdown, ！施加初始下层温度
ALLSEL
timint,off    !先作稳态分析,确定初始条件
time,0.01   !设定只有一个子步的时间很小的载荷步
deltim,0.01
solve!得到初始温度分布


!进行瞬态分析
time,60
timint,on !打开瞬态功能
autots,on!打开自动时间步长
deltim,1!设置时间步长，最大及最小时间步长
tim=0
*do,i,0,29.88,1.66
tim=tim+2
time,tim
nsel,s,loc,y,k1
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-0.99,0.99
sf,all,hflux,-4600
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo

solve
finish
save

/post1
Set,last
Plnsol,temp
Allsel
Finish

caicai1959 发表于 2009-4-11 21:13:02

有人吗?

caicai1959 发表于 2009-4-12 18:00:10

怎么没人帮我看看啊。。急啊。。

reltek 发表于 2010-3-2 10:27:29

KEYOPT,2,3,0
KEYOPT,2,4,1   
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,1   
KEYOPT,2,9,0


请问以上命令是什么意思？（别告我去看help，英文很烂）
页: [1]
查看完整版本: 求助！请高手帮我看看温度出了什么问题！