【请教】:传热问题选择三角形单元和矩形单元差别咋那么大?
传热问题,分别采用矩形单元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 编辑 ]
页:
[1]