热分析

写的命令流如下
在求解出温度分布之后，求解温度应力时，提示：rth文件没有，请问如何生成呢
命令流如下：
finish
/clear
/filname,s
/title, temperature analysis
!c35砼物理参数
u=0.15       !泊松比
alpx=1e-5    !热膨胀系数,1/oC
kxx=2.49     !导热系数,W/m.oC
dens=2420    !密度,kg/m3
c=0.980      !比热,J/Kg.oC
cc=959.08    !对流系数kJ/m3. 侧面
cf=369.4     !对流系数kJ/m3. 顶面
tt=21.2      !浇筑/入模温度oC
te=18        !环境温度oC  
!c20砼垫层物理参数
u2=0.2        !泊松比
alpx2=1e-5    !热膨胀系数,1/oC
kxx2=65.28     !导热系数,W/m.oC
dens2=2400    !密度,kg/m3
c2=0.920      !比热,J/Kg.oC
tt2=15         !地基温度oC
/prep7
ET,1,SOLID70
!c20砼 基础
MP,ALPX,1,alpx2   
MP,KXX,1,kxx2      
MP,DENS,1,dens2   
MP,C,1,c2
!c35砼材料特性2~4号
MP,ALPX,2,alpx
MP,KXX,2,kxx      !定义上面层2.7材料属性编号2
MP,DENS,2,dens   
MP,C,2,c
MP,ALPX,3,alpx   
MP,KXX,3,kxx    !定义上面层2.7材料属性编号3   
MP,DENS,3,dens   
MP,C,3,c
MP,ALPX,4,alpx     
MP,KXX,4,kxx      !定义上面层2.6材料属性编号4
MP,DENS,4,dens   
MP,C,4,c
block,-41.85,-11.85,-11.75,11.75,-1,8
block,11.85,41.85,-11.75,11.75,-1,8
k,17,-11.85,-7.5,-1
k,18,-11.85,-7.5,8
k,19,-11.85,7.5,-1
k,20,-11.85,7.5,8
k,21,-9.35,-5,-1
k,22,-9.35,-5,8
k,23,-9.35,5,-1
k,24,-9.35,5,8
a,17,18,20,19
a,19,20,24,23
a,21,22,24,23
a,17,18,22,21
a,17,19,23,21
a,18,20,24,22
va,13,14,15,16,17,18
k,25,11.85,-7.5,-1
k,26,11.85,-7.5,8
k,27,11.85,7.5,-1
k,28,11.85,7.5,8
k,29,9.35,-5,-1
k,30,9.35,-5,8
k,31,9.35,5,-1
k,32,9.35,5,8
a,25,26,28,27
a,27,28,32,31
a,29,30,32,31
a,25,26,30,29
a,25,27,31,29
a,26,28,32,30
va,19,20,21,22,23,24
a,21,22,30,29
a,23,24,32,31
a,22,24,32,30
a,21,23,31,29
va,25,26,27,28,21,15
vglue,all
vsbw,all
wpoffs,,,2.7
vsbw,all
wpoffs,,,2.7
vsbw,all
wprota,,90
wpoffs,,,7.5
vsbw,all
wpoffs,,,-15
vsbw,all
wpcsys
allsel
lsel,s,loc,x,-41.85,
lsel,a,loc,x,41.85
lsel,a,loc,x,11.85
lsel,a,loc,x,-11.85
lsel,r,loc,y,-11.74,-7.51
lesize,all,,,5
allsel
lsel,s,loc,x,-41.85,
lsel,a,loc,x,41.85
lsel,a,loc,x,11.85
lsel,a,loc,x,-11.85
lsel,r,loc,y,7.51,11.74
lesize,all,,,5
allsel
lsel,s,loc,x,-41.85,
lsel,a,loc,x,41.85
lsel,a,loc,x,11.85
lsel,a,loc,x,-11.85
lsel,r,loc,y,-7.49,7.49
lesize,all,,,20
allsel,
lsel,s,loc,y,11.75,
lsel,a,loc,y,-11.75
lsel,r,loc,x,-41.84,-11.86
lesize,all,,,30
allsel,
lsel,s,loc,y,11.75,
lsel,a,loc,y,-11.75
lsel,r,loc,x,11.86,41.84
lesize,all,,,30
allsel,
lsel,s,loc,y,7.5,
lsel,a,loc,y,-7.5
lsel,r,loc,x,11.86,41.84
lesize,all,,,30
allsel,
lsel,s,loc,y,7.5,
lsel,a,loc,y,-7.5
lsel,r,loc,x,-41.84,-11.86
lesize,all,,,30
allsel
lsel,s,loc,x,-11.84,-9.36
lsel,a,loc,x,9.36,11.84
lesize,all,,,5
allsel
lsel,s,loc,x,-9.34,9.34
lesize,all,,,20
allsel
lsel,s,loc,x,-9.35
lsel,a,loc,x,9.35
lsel,r,loc,y,-4.99,4.99
lesize,all,,,20
allsel
lsel,s,loc,x,-41.85
lsel,r,loc,y,-11.75
lesize,all,,,4
allsel
lsel,s,loc,x,41.85
lsel,r,loc,y,11.75
lesize,all,,,4
allsel
lsel,s,loc,x,-41.85
lsel,r,loc,y,11.75
lesize,all,,,4
allsel
lsel,s,loc,x,41.85
lsel,r,loc,y,-11.75
lesize,all,,,4
allsel
lsel,s,loc,x,-11.85
lsel,r,loc,y,-7.5
lesize,all,,,4
allsel
lsel,s,loc,x,11.85
lsel,r,loc,y,-7.5
lesize,all,,,4
allsel
lsel,s,loc,x,-11.85
lsel,r,loc,y,7.5
lesize,all,,,4
allsel
lsel,s,loc,x,11.85
lsel,r,loc,y,7.5
lesize,all,,,4
allsel
lsel,s,loc,x,-9.35
lsel,r,loc,y,-5
lesize,all,,,4
allsel
lsel,s,loc,x,9.35
lsel,r,loc,y,5
lesize,all,,,4
allsel
lsel,s,loc,x,-9.35
lsel,r,loc,y,5
lesize,all,,,4
allsel
lsel,s,loc,x,9.35
lsel,r,loc,y,-5
lesize,all,,,4
allsel          !垫层0.5
mshape,0
mshkey,1
vsel,s,loc,z,-0.51,0.01
vatt,1,1,1
vmesh,all
allsel          !2.7
mshape,0
mshkey,1
vsel,s,loc,z,-0.01,2.71
vatt,2,2,1
vmesh,all
allsel          !2.7
mshape,0
mshkey,1
vsel,s,loc,z,2.69,5.41
vatt,3,3,1
vmesh,all
allsel          !2.6
mshape,0
mshkey,1
vsel,s,loc,z,5.39,8.01
vatt,4,4,1
vmesh,all
!输入热生成率函数
*SET,_FNCNAME,'heats'   
*DIM,%_FNCNAME%,TABLE,6,16,1
! Begin of equation: 2400*0.9462*290.5*0.175*exp(-0.8426*({TIME}))
%_FNCNAME%(0,0,1)= 0.0, -999
%_FNCNAME%(2,0,1)= 0.0
%_FNCNAME%(3,0,1)= 0.0
%_FNCNAME%(4,0,1)= 0.0
%_FNCNAME%(5,0,1)= 0.0
%_FNCNAME%(6,0,1)= 0.0
%_FNCNAME%(0,1,1)= 1.0, -1, 0, 2400, 0, 0, 0
%_FNCNAME%(0,2,1)= 0.0, -2, 0, 0.9462, 0, 0, -1
%_FNCNAME%(0,3,1)=   0, -3, 0, 1, -1, 3, -2
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 290.5, 0, 0, -3
%_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -3, 3, -1
%_FNCNAME%(0,6,1)= 0.0, -1, 0, 0.175, 0, 0, -2
%_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,8,1)= 0.0, -1, 0, 0, 0, 0, 0
%_FNCNAME%(0,9,1)= 0.0, -2, 0, 1, 0, 0, -1
%_FNCNAME%(0,10,1)= 0.0, -4, 0, 1, -1, 2, -2
%_FNCNAME%(0,11,1)= 0.0, -1, 0, 0.8426, 0, 0, -4
%_FNCNAME%(0,12,1)= 0.0, -2, 0, 1, -4, 3, -1
%_FNCNAME%(0,13,1)= 0.0, -1, 0, 1, -2, 3, 1
%_FNCNAME%(0,14,1)= 0.0, -1, 7, 1, -1, 0, 0
%_FNCNAME%(0,15,1)= 0.0, -2, 0, 1, -3, 3, -1
%_FNCNAME%(0,16,1)= 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 2400*0.9462*290.5*0.175*exp(-0.8426*({TIME}))
/solu
ANTYPE,4                    !选择瞬态求解                 
TRNOPT,FULL
LUMPM,0
alls
esel,s,mat,,2,4
nsle,s
ic,all,temp,tt     !给3层混凝土施加初始温度
nropt,full
esel,s,mat,,3,4   !杀死2/3层混凝土单元
ekill,all
esel,s,live      
esel,s,mat,,1     !选择活单元
nsle,s            
ic,all,temp,tt2     !给基础施加初始温度
alls
asel,s,,,4
asel,a,,,9
asel,a,,,29,33,4
asel,a,,,66
asel,a,,,76,78
asel,a,,,81
asel,a,,,91,92
asel,a,,,94,96
asel,a,,,98,100
asel,a,,,120
asel,a,,,165,167,2
nsla,s
sfa,all,1,conv,cc,te  !砼\空气对流系数,4侧面 环境温度
alls
asel,s,loc,z,2.699,2.701
nsla,s
sfa,all,1,conv,cf,te     !砼\空气对流系数,顶面 环境温度
alls
vsel,s,loc,z,-0.01,2.701
bfv,all,hgen,%HEATS%
alls
ALLSE                         !选择所有
TRNOPT,FULL  
*DO,DAY,1,3,1        !设置循环计算
TIME,DAY                  !设置荷载步时间
AUTOTS,0                 !关闭自动时间步长
KBC,0                    !设置阶梯荷载步
TSRES,ERASE              !Defines an array of keytimes at which the time-stepping strategy changes
OUTRES,ALL,ALL           !输出所有项
SOLVE   
*ENDDO                   !终止循环计算

提示没有找到 .rth 文件，说明你的热分析有问题，没有得到结果。
    你可以验证一下：在热分析后，立即到后处理查看是否有温度计算结果。如果确实没有温度计算结果，可以到 Utility Menu > List > Files > err File，打开错误文件看一下有什么问题。
    顺便说一句，你的循环语句中，循环变量 DAY，可能你认为它的单位是 "天"，实际上不是，根据你的材料等数据判断，你的时间单位应该是 秒。

是这样，热分析后有结果，1-3的温度分布图都有
但保存后，再从ANSYS读入db文件，提示前几步的数据就不存在，就是没保存，只有最后一步的
请问这是什么情况呢

试了一下，不存在你说的问题。求解结束时有三个载荷步的结果。
存盘，退出，再进入 ANSYS，恢复数据库 s.db，到后处理查看，还是有三个载荷步的结果。

哦，这样啊，你用的是哪个版本的ANSYS呢
难道是我电脑有问题