各位大侠,遇到一个奇怪的现象,明明加载的是正的生热率,出来却是负的温度,已经被折磨了好几天了,请大家帮帮忙!!先谢过了,程序如下:
/prep7
et,1,SOLID90 ! thermal brick element
mp,kxx,1,2.1 ! thermal conductivity, W/m-C
mp,c,1,420 ! specific heat, J/kg-C
block,0,0.030,0,0.02,0,0.0004 ! create and mesh model
FLST,5,4,4,ORDE,4
FITEM,5,2
FITEM,5,4
FITEM,5,-5
FITEM,5,7
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,30, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,1
FITEM,5,3
FITEM,5,6
FITEM,5,8
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,20, , , , ,1
FLST,5,3,4,ORDE,3
FITEM,5,9
FITEM,5,11
FITEM,5,-12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,1, , , , ,1
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
*dim,pa,,4453,1 ! 定义数组读数组
*vread,pa(1,1),vector,txt,,jik,1,4453
(F20.3)
/SOL
antype,static ! steady-state heat transfer analysis
time,1.0
*get,nmax,node,,num,max
*do,i,1,nmax
BF,i,HGEN,pa(i)
*ENDDO
nsel,s,ext ! exterior beam nodes
sf,all,conv,6.0,0 ! hf=6.0 W/m^2-C and Tbulk=0 degrees C
nsel,all
save
solve
fini
/post1
set,last
plnsol,temp |