关于热--应力耦合结果问题的求助

曾大兵

我在做一个热--应力耦合的分析,现在得出的结果(见附件云图)很奇怪,大家帮我分析看看,为什么进行热应力耦合计算后,结构的变形方向和没考虑温度是完全相反,而且应力大小相差很大!!!

下面是我进行热应力耦合分析时的命令流文件,大家帮我看看是不是那里有问题:
*SET,a1,0.55
*SET,b1,0.37
*SET,r1,2.4
*SET,r2,1.1
*SET,dens,7.8e3
/PREP7  
ET,1,SOLID87
MP,KXX,1,49.3   
CYLIND,3.0,0.5,0,a1,0,360,  
CYL4,0,2.627,0.13, , , ,a1  
CSYS,1  
FLST,3,1,6,ORDE,1   
FITEM,3,2   
VGEN,12,P51X, , , ,30, , ,0
FLST,3,12,6,ORDE,2  
FITEM,3,2   
FITEM,3,-13
VSBV,       1,P51X  
CYLIND,r1,r2,0,b1,0,20,
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN, ,P51X, , , ,5, , , ,1
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN,12,P51X, , , ,30, , ,0
FLST,3,12,6,ORDE,2  
FITEM,3,1   
FITEM,3,-12
VSBV,      14,P51X  
CYLIND,0.15,0,0,b1,0,360,   
CSYS,0  
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN, ,P51X, , , ,0.85, , , ,1  
CSYS,1  
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN, ,P51X, , , ,15, , , ,1
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN,12,P51X, , , ,30, , ,0
FLST,3,12,6,ORDE,2  
FITEM,3,1   
FITEM,3,-12
VSBV,      13,P51X  
CYLIND,1.85,0,a1,a1+0.37,0,360,
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,14  
VGLUE,P51X  
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL  
/VIEW,1,,1  
/ANG,1  
/REP,FAST   
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
ESIZE,0.2,0,
MSHAPE,1,3D
MSHKEY,0
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
FINISH  
/SOL
FLST,2,1,5,ORDE,1   
FITEM,2,131
/GO
SFA,P51X,1,CONV,7,27
ASEL,S,LOC,Z,0.910009,0.920001  
APLOT   
FLST,2,1,5,ORDE,1   
FITEM,2,50  
/GO
SFA,P51X,1,CONV,300,160
ALLSEL,ALL  
/REPLOT
SAVE
ANTYPE,0
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  

/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0
KEYW,PR_CFD,0   
/GO

/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Structural   

FINISH  
/PREP7  
LSCLEAR,SOLID   
ETCHG,TTS   

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.06e11
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,1,REFT,,,27
MPDATA,ALPX,1,,1.62e-5  
  
FINISH  
/SOL
FLST,2,24,5,ORDE,24
FITEM,2,8   
FITEM,2,-9  
FITEM,2,12  
FITEM,2,-13
FITEM,2,16  
FITEM,2,-17
FITEM,2,20  
FITEM,2,-21
FITEM,2,24  
FITEM,2,-25
FITEM,2,28  
FITEM,2,-29
FITEM,2,32  
FITEM,2,-33
FITEM,2,36  
FITEM,2,-37
FITEM,2,40  
FITEM,2,-41
FITEM,2,44  
FITEM,2,-45
FITEM,2,48  
FITEM,2,-49
FITEM,2,52  
FITEM,2,-53
  
/GO
DA,P51X,ALL,
LDREAD,TEMP,,, , ,'file','rth',' '  
ASEL,S,LOC,Z,0.919999,0.920001  
APLOT   
FLST,2,1,5,ORDE,1   
FITEM,2,50  
/GO
SFA,P51X,1,PRES,-3.56e6  
ALLSEL,ALL  
/REPLOT

曾大兵

再补发一个应变图 1# 曾大兵