命令流:
/PREP7
ET,1,SOLID70
!定义混凝土导热性,比热,密度等
MPTEMP,,20,300,400,500,600,700
MPDATA,KXX,1,,1.34,1.34,1.22,1.1,0.97,0.85
MPTEMP,,800
MPDATA,KXX,1,,0.72
MPTEMP !清除温度场
MPTEMP,,20,500,600,700,800
MPDATA,C,1,,1115.7,1115.7,1812.7,6757.7,1115.7
MP,DENS,1,2300
!建立模型
k,1,0,0,0
k,2,0,150,0
k,3,0,150,300
k,4,0,0,300
k,5,1600,0,0
k,6,1600,150,0
k,7,1400,150,300
k,8,1400,0,300
k,9,1200,0,0
k,10,1200,150,0
k,11,1200,150,300
k,12,1200,0,300
V,1, 2, 3, 4,9,10,11,12
V,9,10,11,12,5,6,7,8
vplot
VGLUE,all
/USER, 1
/VIEW, 1, -0.266795396550 , -0.672493524827 , 0.690342433467
/ANG, 1, 5.02290010730
/REPLO
!划分单元
lesize,1,,,5
lesize,4,,,10
lesize,5,,,40
MSHAPE,0,3D
MSHKEY,1
vmesh,1
lesize,6,,,5
lesize,12,,,10
lesize,17,,,5
lesize,13,,,10
MSHAPE,1,3D
MSHKEY,0
vmesh,2
!定义温度变化和对流系数变化
*DIM,wendu,TABLE,18,1,1,TIME,TEMP,
!*
*SET,WENDU(1,0,1) , 0
*SET,WENDU(1,1,1) , 20
*SET,WENDU(2,0,1) , 240
*SET,WENDU(2,1,1) , 60
*SET,WENDU(3,0,1) , 480
*SET,WENDU(3,1,1) , 100
*SET,WENDU(4,0,1) , 720
*SET,WENDU(4,1,1) , 140
*SET,WENDU(5,0,1) , 960
*SET,WENDU(5,1,1) , 180
*SET,WENDU(6,0,1) , 1200
*SET,WENDU(6,1,1) , 220
*SET,WENDU(7,0,1) , 1440
*SET,WENDU(7,1,1) , 260
*SET,WENDU(8,0,1) , 1680
*SET,WENDU(8,1,1) , 300
*SET,WENDU(9,0,1) , 1920
*SET,WENDU(9,1,1) , 340
*SET,WENDU(10,0,1) , 2160
*SET,WENDU(10,1,1) , 380
*SET,WENDU(11,0,1) , 2400
*SET,WENDU(11,1,1) , 420
*SET,WENDU(12,0,1) , 2640
*SET,WENDU(12,1,1) , 460
*SET,WENDU(13,0,1) , 2880
*SET,WENDU(13,1,1) , 500
*SET,WENDU(14,0,1) , 3120
*SET,WENDU(14,1,1) , 540
*SET,WENDU(15,0,1) , 3360
*SET,WENDU(15,1,1) , 580
*SET,WENDU(16,0,1) , 3600
*SET,WENDU(16,1,1) , 620
*SET,WENDU(17,0,1) , 3840
*SET,WENDU(17,1,1) , 660
*SET,WENDU(18,0,1) , 4080
*SET,WENDU(18,1,1) , 700
*DIM,DUILIU,TABLE,6,1,1,TIME,,
!*
*SET,DUILIU(1,0,1) , 0
*SET,DUILIU(1,1,1) , 11.63
*SET,DUILIU(2,0,1) , 1080
*SET,DUILIU(2,1,1) , 11.64
*SET,DUILIU(3,0,1) , 2280
*SET,DUILIU(3,1,1) , 17.45
*SET,DUILIU(4,0,1) , 2880
*SET,DUILIU(4,1,1) , 23.26
*SET,DUILIU(5,0,1) , 3480
*SET,DUILIU(5,1,1) , 34.89
*SET,DUILIU(6,0,1) , 4080
*SET,DUILIU(6,1,1) , 46.52
FINISH
/SOLU
ANTYPE,TRANS
TRNOPT,FULL
LUMPM,0
allsel,all
IC,all,TEMP,20 !定义初始温度
SFA,2,1,CONV, %DUILIU% , %WENDU%
SFA,3,1,CONV, %DUILIU% , %WENDU%
SFA,5,1,CONV, %DUILIU% , %WENDU%
SFA,7,1,CONV,15,20
SFA,8,1,CONV,15,20
SFA,9,1,CONV,15,20
SFA,10,1,CONV,15,20
SFA,11,1,CONV,15,20
SFA,4,1,CONV,15,20
OUTRES,ALL,ALL !输出所有项
TIME,4080
AUTOTS,1
DELTIM,60, , ,1
KBC,0
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
neqit,50 !每一子步中方程的迭代次数限值
pred,on !打开预测器
save
solve
finish
按照这个计算的时候,最后的最大温度才20多度,
可是我在施加对流边界条件的时候,环境温度最高设为700度,
不知道哪里出了问题?
请大家帮忙 看看 谢谢了 |