【讨论】多层的梯度材料在室温下还存在应力？

hangch

最近在做多层梯度材料的温度场和应力场，取一种最极端的情况，我把两层梯度材料放在室温条件下计算，温度没有变化，竟然还存在应力场，不知是何缘故？
采用间接法计算
温度场的
/CLEAR,NOSTART                                                !清零
/FILNAM,PLASMASPRAY                                           !制定工作文件
/TITLE,RESIDUAL STRESS ANALYSIS !制定工作标题
/UNITS,SI  
/CONFIG,NRES,10000
/PREP7  
ET,1,SOLID70  
ET,2,SOLID70  
  
MP,DENS,1,5250  
MP,C,1,800   
MP,KXX,1,2.5  
  
MP,DENS,2,7800  
MP,C,2,448  
MP,KXX,2,8  
/VIEW,1,1,1,1  
BLOCK,0,0.018,0,0.016,0,0.002      BLOCK,0,0.018,0,0.016,0.002,0.0023   
vglue,all  
SAVE  
LESIZE,1,,,16
LESIZE,4,,,18  
LESIZE,9,,,2
MSHKEY,1  
VATT,1,1,1
VMESH,1
LESIZE,8,,,16
LESIZE,5,,,18
LESIZE,25,,,1
MSHKEY,1  
VATT,2,1,2
VMESH,3  
/PNUM,DEFA  
EPLOT  
FINISH  
  
/SOLU
sfa,1,,conv,5,25, !w/m**2*c,下底面
sfa,8,,conv,10,25, !上表面
sfa,3,,conv,10,25, !材料1的四个侧面，3-6
sfa,4,,conv,10,25,
sfa,5,,conv,10,25,
sfa,6,,conv,10,25,
sfa,13,,conv,10,25, !材料2的四个侧面，13-16
sfa,14,,conv,10,25,
sfa,15,,conv,10,25,
sfa,16,,conv,10,25,
!开始计算   
  ANTYPE,TRANS  
  TUNIF,25                         !初始温度25度  
  KBC,1  
   
  *SET,DELT,1  
   TIMINT,ON  
   AUTOTS,OFF  
   
*DO,TM,1E-6,10.45,0.05  
  /VIEW,1,1,1,1  
  /ANG,1  
  /REP,FAST  
   
  !FLST,2,1,5,ORDE,1  
  !FITEM,2,8  
   
  /GO  
  !*  
  TIME,TM  
  AUTOTS,0  
  DELTIM, , , ,1  
  KBC,0  
  !*  
  TSRES,ERASE  
  !FLST,2,1,5,ORDE,1  
  !FITEM,2,8  
  /GO  
  !*  
  !*  
  !SFA,P51X,1,HFLUX,%FLUX2%  
  /STATUS,SOLU  
  SOLVE  
  OUTRES,ALL,ALL,  
  /PSF,HFLUX,,2  
   
  /REPLOT  
      
*ENDDO  
OUTRES,ALL,ALL,  
save
FINISH  
  
应力场的
/SOLU
LSCLEAR,SOLID
  
/PREP7  
ETCHG,TTS
KEYOPT,1,2,0
KEYOPT,1,4,0
KEYOPT,1,7,0
KEYOPT,1,8,0
  
KEYOPT,2,2,0
KEYOPT,2,4,0
KEYOPT,2,7,0
KEYOPT,2,8,0
  
MP,EX,1,3.66E11  
MP,ALPX,1,7.962E-6  
MP,NUXY,1,0.25  
  
MP,EX,2,1.2E11  
MP,ALPX,2,11.5E-6  
MP,PRXY,2,0.3  
  
ALLSEL,ALL
FINISH  
  
/SOLU
    DK,1,UX,0,,UY,UZ      !固定底面左上角的顶点
   *DO,tm,1E-6,5,0.05 !设置时间步长和区间  
   TIME,tm !当前时间为tm  
  LDREAD,TEMP,,,TIME,0,'PLASMASPRAY','RTH' !读入时间tm时的温度分布  
  DELTIM,,,,20 !初始步长  
  SOLVE  
  *ENDDO  
  FINISH  
  
温度场当然没有变化，100步时x向的应力场如下图
请高手指教原因，是不是没有设参考温度所致。

hangch

温度场如下：