热分析模型和结构分析节点,单元都是一致的,为什么温度不能读入

andylau436

我热分析和结构分析都是两个四边形,热分析模型两个面是挨着的.但是结构分析的模型,两个面之间有间隙,现在我想把热分析求得的温度场作为体载荷加载到结构模型上,热模型和结构模型节点和单元都对应,但是为什么就是不能将温度加载上去呢?我做热膨胀结构分析,请高手赐教!谢谢了.这是命令流!
/filename,pengzhangfeng
/title,pengzhangfeng
/PREP7  
eq_alpx=1
ET,1,PLANE55
mptemp,1,20,400,800,1200
mp,kxx,1,1.15*1000*60*60*60
mpdata,c,1,1,0.8*60*60e9,0.9*60*60e9,1.02*60*60e9,1.2*60*60e9   
mp,alpx,1,eq_alpx*8.5e-6
mp,dens,1,2.95e-6   
mp,ex,1,6.3e6   
mp,nuxy,1,0.21  

!包底第一永久层!          高铝
mp,kxx,2,0.5*1000*60*60*60           
mpdata,c,2,1,0.61*60*60e9,0.75*60*60e9,1.175*60*60e9,1.32*60*60e9               
mp,alpx,2,eq_alpx*5.8e-6              
mp,dens,2,2.8e-6            
mp,ex,2,5.7e6
mp,nuxy,2,0.21

RECTNG,0,120,0,120,
RECTNG,120,240,0,120,

CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  

LESIZE,ALL, , ,3, ,1, , ,1,
FLST,5,2,5,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,1
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
NUMMRG,NODE
/SOL
!*  
ANTYPE,0
AUTOTS,1
TIME,0.01   
FLST,2,1,4,ORDE,1   
FITEM,2,6   
!*  
/GO
DL,P51X, ,TEMP,1000,0   
FLST,2,1,4,ORDE,1   
FITEM,2,4   
/GO
!*  
SFL,P51X,CONV,12*60*60*60, ,30,
ALLSEL,ALL  
/STATUS,SOLU
SOLVE   
save
/POST1  
!*  
/EFACET,1   
PLNSOL, TEMP,, 0

fini
/clear
/prep7
/title,jiegou
eq_alpx=1
et,1,plane42
KEYOPT,1,3,2 !打开PLANE42平面应变选项
ET,2,TARGE169
ET,3,CONTA171
KEYOPT,3,12,0 !选择接触类型(CONTA171)为STANDARD
R,2,950000,,1,0.01,, !设置实常数,法向接触刚度因子1，最大渗透范围0.001

mptemp,1,20,400,800,1200
mp,alpx,1,eq_alpx*8.5e-6
mp,dens,1,2.95e-6   
mp,ex,1,6.3e6   
mp,nuxy,1,0.21  

!包底第一永久层!          高铝
mp,alpx,2,eq_alpx*5.8e-6              
mp,dens,2,2.8e-6            
mp,ex,2,5.7e6
mp,nuxy,2,0.21

RECTNG,0,120,0,120,
RECTNG,122,242,0,120,  

CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1

LESIZE,ALL, , ,3, ,1, , ,1,
FLST,5,2,5,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,1
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  


LSEL,S, , ,       2
nsll,s,1
TYPE,2 !桩为刚性面TARGET169单元
REAL,2
ESURF,ALL


LSEL,S, , ,       8
nsll,s,1
TYPE,3 !桩为刚性面TARGET169单元
REAL,2
ESURF,ALL

/solu
LDREAD,TEMP,1,1, , ,'pengzhangfeng','rth',' '
ALLSEL,ALL  
LPLOT   
FLST,2,4,4,ORDE,4   
FITEM,2,1   
FITEM,2,3   
FITEM,2,5   
FITEM,2,7   
DL,P51X, ,SYMM  


ALLSEL,ALL !全部选定，避免前面操作的漏选
!**求解参数设置**
ANTYPE,STATIC !分析类型为静力
ACEL,,10e3 !施加重力加速度9.8
NLGEOM,ON !打开大变形选项
NROPT,FULL,,ON !打开牛顿-辛普生选项，加快收敛速度
SOLVE

/post1
plnso,s,1