请教高手我这个例子为什么不能生成rth文件

redbee123

我是初学者，看一本书上有个表面残余应力的例子，是先用plane67求温度场，再转换成plan42单元，计算热应力。但是我只能得到rst文件，无论如何得不到rth文件。我用搜索引擎查了很多问答，也没搞太明白。有人说rth是单纯的热计算，有应力计算就会产生rst文件，我看了很久也没查出问题，请高手帮忙看看，谢谢。

/FILEname,exercise6
/TITLE,THERMALS STRESS ANALYSIS OF A MOLTEN NI PARTICLE IN THERMAL SPRAYING
!定义工作标题
/PREP7
et,1,67,,,1           !定义单元
mp,ex,1,0.204         !输入材料1弹性模量
mp,nuxy,1,0.28         !输入材料1材料泊松比
mp,alpx,1,13.36e-6     !输入材料1线膨胀系数
mp,kxx,1,87.86e-6      !输入材料1导热系数
mptemp,1,0,25,355,360,365,1454   !建立温度参数表
mpdata,enth,1,1,0,8.48e-5,1.45e-3,1.50e-3,1.55e-3,6.65e-3 !输入不同温度下材料1的焓值
mp,ex,2,0.215         !输入材料2弹性模量
mp,nuxy,2,0.26         !输入材料2材料泊松比
mp,alpx,2,11.36e-6     !输入材料2线膨胀系数
mp,dens,2,7.9e-9        !输入材料2密度
mp,kxx,2,7.1e-5      !输入材料2导热系数
mp,c,2,460.24        !输入材料2的比热

k,1,,22.5,0           !设置关键点
k,2,10,22
k,3,20,19.3
k,4,30,16.5
k,5,40,13.8
k,6,50,12
k,7,60,9.5
k,8,70,7.5
k,9,80,5.5
k,10,90,4.5
k,11,100,2.8
k,12,110,2.8
k,13,200,0
k,14,0,0
*do,i,1,13
lstr,i,i+1
*enddo
lstr,14,1                 !连接关键点14和1，生成线
al,all                    !选择所有线生成面
RECTING,0,2000,-1000,0,    !生成矩形面
aglue,all                  !将面1与面Glue在一起
numcmp,area                 !压缩面编号
asel,s,area,,2              !选择面2
aatt,2,                     !将面2赋予材料属性2
asel,all                    !选择所有面
amesh,all                   !划分网格
NUMMEG,NODE                 !合并重合节点
SAVE,EXERCISE61,db
FINISH
/solu                       !进入求解器
antype,4                    !指定分析类型为瞬态分析
solcontrol,on               !激活优化的非线性求解控制
ESEL,S,MAT,,1               !选择材料1对应的所有单元
NSLE,S                      !选择所选单元对应的所有节点
ic,all,temp,1454              !在熔滴表面施加初始边界条件温度1454摄氏度
ESEL,S,MAT,,2                !选择材料2对应的所有单元
NSLE,S                      !选择所选单元对应的所有节点
ic,all,temp,25              !在基底施加初始边界条件温度25摄氏度
allsel,all                  !选择全部实体
time,100                     !设置计算时间为100微秒
autots,on                    !打开自动时间步长
deltim,1                     !指定时间步长1
kbc,0                        !按rampped方式加载
outres,all,all               !将每一步计算结果都输入到结果文件
solve                        !开始求解
SAVE,EXERCISE62,db           !将当前计算结果保存到文件
FINISH
/post1
ESEL,S,MAT,,1                 !选择材料1对应所有单元
PLNSOL,TEMP                    !绘制100微秒时熔滴的温度场分布
FINISH
/post26
nsol,2,1,temp,,point-A          !定义变量point-A
nsol,3,42,temp,,point-B          !定义变量point-B
plvar,2,3                        !绘制A,B两点的温度随时间的变化曲线
allsel,all
finish
/prep7                           !进入前处理器，进行应力场分析
etchg,tts                        !转变单元类型
keyopt,1,3,1                     !为单元设置轴对称关键选项
NSEL,S,LOC,X,0                   !选择X=0处所有节点
D,all,0,,,,UX                    !定义UX=0约束
tref,625                         !设置参考温度为625
allsel,all
LDREAD,TEMP,,,,,'rth',''         !从文件exercise6.rth读取温度荷载
/solu
solve
save,exercise63,db
finish

redbee123

在这一步做不下去了
LDREAD,TEMP,,,,,'rth',''         !从文件exercise6.rth读取温度荷载

woomii

单元搞错了吧。67是热电耦合单元，出来的结果就是.rst后缀的，也可以查看温度结果。但这里面不需要电分析啊。
应该用77，纯热分析单元，结果文件才是.rth。