【请教】:传热问题选择三角形单元和矩形单元差别咋那么大?

ParkerZeng

传热问题,分别采用矩形单元plane55和三角形单元plane35,其它一切不变,代码如下:

/prep7
et,1,35,,,1                                                                               !采用矩形单元为；et,1,55,,,1
et,2,35,,,1                                                                               !采用矩形单元为；et,2,35,,,1
mp,  dens,  1,   7750                                                                !密度
mp,  kxx,    1,   14                                                                   !导热系数
mp,  C,       1,   510                                                                 !比热容
mp,  dens,  2,  7860                                                                !密度  
mp,  kxx,    2,  49.8                                                                 !导热系数      
mp,  c,       2,  465                                                                   !比热容
k, 1,   0,           0
k, 2,   0.15,      0
k, 3,   0.15,      0.15
k, 4,   0,           0.15
k, 5,   0,           0.16
k, 6,   0.03,       0.16
k, 7,   0.03,       0.15
a,1,2,3,4
a,4,5,6,7
aglue,all
numcmp,area
type,1
mat,1
SMRTSIZE,1
amesh,1
type,2
mat,2
SMRTSIZE,1
amesh,2

/solu
antype,trans                                                          !瞬态分析
!********************先做一步稳态分析初始化温度场***********************      
TIMINT,off                                                           !打开瞬态时间积分选项
nsubst,1                                                                 !子步数为1
time,1e-5                                                               !时间为0.0011秒
deltim,1e-5                                                           ! run one substep to initialize temperature               
outres,all,all                                                           !设定将所有求解信息写入输出数据库
esel,s,mat,,2
nsle,s,1
CM,NS1,Node
d,NS1,temp,298                                                      !设定节点处的温度为298K
esel,s,mat,,1
nsle,s,1
CM,NS2,Node
d,NS2,temp,1000                                               !设定节点处的温度为1000K  
allsel,all
solve
!************进入瞬态分析*******************
ddele,all,temp                                             ! delete initial temperature specs
timint,on                                                   ! turn on time integration
tintp,,,,1                                                  ! specify use EULER bakward diff. integration
lnsrch,on                                                ! turn on linesearch
time,60                                             ! set final time to 60 seconds
nsubst,200
solve      

下面两张图分别是采用矩形单元和三角形单元在第二步、第一子步时间为0.30001秒的温度云图。采用矩形单元的模型应该没什么问题，但是采用三角形单元的模型其最低温度是290度。模型中没有施加任何边界，最低温度怎么会比298K还低呢？请前辈们多多指教！！


[ 本帖最后由 ParkerZeng 于 2007-3-21 09:36 编辑 ]