【讨论】请各位帮我看看这个热应力分析的问题在哪？

anzi888

命令流：
/FILNAME,test
/title,an
/prep7
!指定单元类型
et,1,87,0
et,2,90,0
!指定单元属性
! 1号材料是bean Q235BH
! 2号材料是cloumn Q345BH-16Mn
! 3号材料是weld  Q235BH
mptemp,1,20,200,500,800,1500
mpdata,ex,1,1,2.064e11,1.997e11,1.348e11,8.106e10,1.0e8
mpdata,ex,2,1,2.1e11,1.95e11,1.55e11,1.0e11,5e7  
mp,alpx,1,11e-6
mp,alpx,2,12e-6
mp,kxx,1,68
mp,kxx,2,46
mp,c,1,480
mp,c,2,450
mp,dens,1,7840
mp,dens,2,7815  
mp,prxy,1,0.285
mp,prxy,2,0.29
  
TB,BKIN,1,5         
TBTEMP,20,1         
TBDATA,1,2.45E8,2.064e10
TBTEMP,200,2
TBDATA,1,1.6E8,1.997e10
TBTEMP,600,3
TBDATA,1,0.7E8,1.348e10
TBTEMP,800,4
TBDATA,1,5E7,8.106e9
TBTEMP,1500,5
TBDATA,1,1E6,1.0e7
  
TB,BKIN,2,5         
TBTEMP,20,1         
TBDATA,1,3.9E8,2.1e10
TBTEMP,200,2
TBDATA,1,2.58E8,1.95e10
TBTEMP,600,3
TBDATA,1,1.37E8,1.55e10
TBTEMP,800,4
TBDATA,1,7.5E7,1.0e10
TBTEMP,1500,5
TBDATA,1,1E6,5e6
  
BLOCK,0,130e-3,0,4e-3,0,10e-3  
wprota,,,90
wpoffs,,,70e-3
vsbw,1
wpoffs,,,5e-3
vsbw,3
wpstyl,defa
WPSTYLE,,,,,,,,0
  
mat,2
type,2
esize,1e-3
MSHAPE,1,3D  
MSHKEY,0
vmesh,4
  
mat,1
type,1
esize,3e-3
MSHAPE,1,3D  
MSHKEY,0
vmesh,1
vmesh,2
  
/SOLU
  
ANTYPE,TRANS
TREF,24
TIMINT,on
tintp,0.005,,,1,0.5,0.2
VSEL,S,,,4
NSLV,S,1
D,ALL,TEMP,1450  
NSEL,INVE
D,ALL,TEMP,30
TIME,5
KBC,1  
NSUBST,2,2,2
OUTRES,all,all
ALLSEL,all
SOLVE
  
OUTRES,all,all         
DDELE,ALL,TEMP                           
TIME,40            !phase change
DELTIME,0.2,0.1,2        
AUTOTS,ON            
SOLVE
  
finish
  
/prep7
etchg,tts
finish
/SOLu
ANTYPE,0
  
nsel,s,loc,x,0 ,0
nsel,a,loc,x,130e-3 ,130e-3
d,all,ux,0  
  
LDREAD,TEMP,,,, ,'test','rth',' '   
time,1
nlgeom,on  
nsubst,15
autots,on  
pred,on  
nropt,full,,on      
TREF,24
allsel,all
solve
  
热应力的最大值应该在中间，可这个刚好相反，居然在两边。
请各位帮我看看出了什么问题？？