一个升温降温的热应力分析模拟，怎么不对呀

飞剑客

！升温与降温，参考温度为10，路面在参考温度基础上升温与降温
！需要得到温度分布图
fini
/cle
/PREP7   
keyw,pr_therm,1

ET,1,PLANE42  
!辅助网格划分之用
ET,2,SOLID5
!…………………………面层模量与泊松比等参数
mp,ex,1,1200e6   
mp,nuxy,1,0.25  
mp,alpx,1,21e-6   
mp,kxx,1,1.2

！…………………………土基层模量与泊松比等参数
mp,ex,2,23000e6   
mp,nuxy,2,0.15  
mp,alpx,2,10e-6   
mp,kxx,2,1.5
！…………………………下基层模量与泊松比等参数
mp,ex,3,2000e6   
mp,nuxy,3,0.25
mp,alpx,3,8e-6   
mp,kxx,3,1
！…………………………土基模量与泊松比等参数
mp,ex,4,200e6   
mp,nuxy,4,0.35
mp,alpx,4,5e-6   
mp,kxx,4,1

k,1,0,0
k,2,8,0
k,3,8,4
k,4,0,4

a,1,2,3,4

lsel,s,,,2,4,2
lesize,all,,,20

lsel,s,,,1,3,2
lesize,all,,,10

MSHKEY,1
amesh,all

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~拉伸下基层  
extopt,esize,4,0
extopt,aclear,1
type,2  
mat,1
vext,1,,,,,0.2
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~拉伸基层  
extopt,esize,4,0
extopt,aclear,1
type,2  
mat,2
vext,2,,,,,0.2

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~拉伸面层  
extopt,esize,4,0
extopt,aclear,1
type,2  
mat,3
vext,7,,,,,0.2

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~拉伸土基  
extopt,esize,4,0
extopt,aclear,1
type,2  
mat,4
vext,1,,,,,-6

extopt,esize,1,0
extopt,aclear,1
type,2  
mat,4
vext,18,,,,-1

extopt,esize,1,0
extopt,aclear,1
type,2  
mat,4
vext,20,,,,1

extopt,esize,1,0
extopt,aclear,1
type,2  
mat,4
vext,19,24,5,1
vext,31,,,1

extopt,esize,1,0
extopt,aclear,1
type,2  
mat,4
vext,21,29,8,-1
vext,26,,,-1

nummrg,all    !合并同位置或等效项目
numcmp,all     ！压缩项目编号
！进入求解器
/solu         

antype,static   
！指定分析类型
asel,s,loc,z,-6  
aplot
DA,ALL,all           
！施加位移约束
asel,s,loc,x,-1  
aplot
DA,ALL,ux        
！施加位移约束
asel,s,loc,x,9  
aplot
DA,ALL,ux         
！施加位移约束
asel,s,loc,y,-1  
aplot
DA,ALL,uy           
！施加位移约束
asel,s,loc,y,5
aplot
DA,ALL,uy
tref,10                  
!参考温度
asel,s,loc,z,0.6
nsll,s,1
bf,all,temp,200         
!施加温度荷载，
allsel
solve


哪里出现问题了，谢谢，指教，我是个初学者

飞剑客

在线等……