求助：瞬态热应力问题？

hanxiao0413

高手帮忙：瞬态热应力问题？

/PREP7  
ET,1,SOLID70
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,1,,43.27
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,2,,0.81  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,3,,1.687
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7850
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,C,1,,450
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,1830
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,C,2,,860
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,3,,2390
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,C,3,,1186
*SET,r2,0.0889        
*SET,t1,0.00919      
*SET,r1,r2-t1         
*SET,t2,0.03            
*SET,r3,r2+t2           
  
CYLIND,r2,r1,-5,5,0,90,
CYLIND,r3,r2,-5,5,0,90,
BLOCK,,10,,10,-5,5,
FLST,2,3,6,ORDE,2   
FITEM,2,1   
FITEM,2,-3  
VPTN,P51X   
VDELE,       4, , ,1
FLST,2,3,6,ORDE,3   
FITEM,2,1   
FITEM,2,5   
FITEM,2,-6  
VGLUE,P51X  
NUMCMP,VOLU

CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y  
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
CM,_Y,VOLU  
VSEL, , , ,       2
CM,_Y1,VOLU
CMSEL,S,_Y  
CMSEL,S,_Y1
VATT,       2, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
CM,_Y,VOLU  
VSEL, , , ,       3
CM,_Y1,VOLU
CMSEL,S,_Y  
CMSEL,S,_Y1
VATT,       3, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  

SMRT,6                     
FLST,5,3,6,ORDE,2   
FITEM,5,1   
FITEM,5,-3  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
/UI,MESH,OFF

FINISH  
/SOL

ALLSEL,ALL  
ANTYPE,4
TRNOPT,FULL
LUMPM,0

TUNIF,50,               
CSYS,1  
ASEL,S,LOC,X,r1
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
SFA,P51X,1,CONV,10e3,350     

ALLSEL,ALL  
CSYS,0  
OUTRES,BASIC,ALL,   
TIME,3600*24*20           
AUTOTS,1
NSUBST,80, , ,1
KBC,0   
TSRES,ERASE
SOLVE  


FINISH  
/POST1
FINISH  
/PREP7  


ETCHG,TTS   
MPTEMP,,,,,,,,  
MPTEMP,1,50
MPTEMP,2,150
MPTEMP,3,250
MPTEMP,4,350
MPDATA,EX,1,,195.7e9         
MPDATA,EX,1,,169.7e9
MPDATA,EX,1,,143.7e9
MPDATA,EX,1,,117.7e9
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,1,REFT,,,  
MPDATA,ALPX,1,,12.1e-6         
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,2,REFT,,,  
MPDATA,ALPX,2,,10.3e-6  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,3,REFT,,,  
MPDATA,ALPX,3,,10.3e-6  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,15e9            
MPDATA,PRXY,2,,0.25
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,3,,10e9              
MPDATA,PRXY,3,,0.25  
FINISH  
/SOL

ALLSEL,ALL  
CSYS,1  
ASEL,S,LOC,X,r1
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
SFA,P51X,1,PRES,12e6           

ALLSEL,ALL  
CSYS,0  
ASEL,S,LOC,Z,5  
FLST,2,3,5,ORDE,3   
FITEM,2,2   
FITEM,2,24  
FITEM,2,30  
/GO
DA,P51X,UZ,

ALLSEL,ALL  
ASEL,S,LOC,Z,-5
FLST,2,3,5,ORDE,3   
FITEM,2,1   
FITEM,2,23  
FITEM,2,29  
/GO
DA,P51X,UZ,

ALLSEL,ALL  
ASEL,S,LOC,X,10
FLST,2,1,5,ORDE,1   
FITEM,2,18  
/GO
SFA,P51X,1,PRES,24.26e6     

ALLSEL,ALL  
ASEL,S,LOC,Y,10
FLST,2,1,5,ORDE,1   
FITEM,2,16  
/GO
SFA,P51X,1,PRES,17.423e6     

ALLSEL,ALL  
ASEL,S,LOC,Y,0  
FLST,2,3,5,ORDE,3   
FITEM,2,5   
FITEM,2,26  
FITEM,2,28  
DA,P51X,SYMM      

ALLSEL,ALL  
ASEL,S,LOC,X,0  
FLST,2,3,5,ORDE,3   
FITEM,2,6   
FITEM,2,25  
FITEM,2,27  
DA,P51X,SYMM

ALLSEL,ALL  
TREF,50,         
LDREAD,TEMP,,, , ,'ZHS','rth',' '   
SOLVE
  
FINISH  
/POST1  


各位高手：模型中间为一钢管，管内流体温度350度，流动时间20天，现在要求每一天所对应的钢管温度、应力。（这里建的是一个四分之一模型）

用上面的命令流进行模拟时，存在下面的问题：

1、流体流动20天时（3600*24*20）的中间时刻（10天时）对应的管子的温度值，与直接流体通过10天时（3600*24*10）算出来的温度结果不相等？而实际情况下应该是相等的。

2、用上面的命令流计算出来的管子内壁的应力结果感觉也不合理。


    哪位高手对瞬态热应力求解比较熟悉的请不吝赐教啊！！！
   

    也请各位简要介绍一下模型瞬态温度场（包括不同的载荷步时）、应力场的计算步骤。