求助：加载高斯热源后怎么分析热内应力

skyprincee

   /SOLU
  LSCLEAR,SOLID  清除数据库中的载荷和载荷步选项

/PREP7
ETCHG,TTS  变热单元为结构单元


MP,EX,1,2.8E4
MP,ALPX,1,7.962E-6
MP,NUXY,1,0.25
FINISH

/SOLU
   ANTYPE,TRANS                                                
   DK,1,UX,0,,UY,UZ                                                !固定底面左上角的顶点
  *DO,TM,1E-6,TT,TINC                                         !设置时间步长和区间
                                 
  TIME,TM                                                               !指定载荷步结束的时间
  LDREAD,TEMP,,,TM,0,'stress','RTH'                           !读入时间tm时的温度分布 原来是TIME，

TM处
                                                                                    应该需要指定参考温度
  DELTIM,,,,1                                                           !初始步长
  SOLVE
  *ENDDO
  FINISH
哪里不对，请高手指正