ANSYS模拟压力产生的变形对温度场的影响

wxwieandy

采用热-结构直接耦合法来模拟压力产生的变形对温度场的影响，但是温度场的分布明显不正确。

命令流如下：
/PREP7
K,1
K,2,4
K,3,4,2
K,4,0,2
K,5,0,-2.5
K,6,4,-2.5
K,7,4,-0.5
K,8,0,-0.5
K,9,1.33
K,10,2.67
K,11,2.33,-0.5
K,12,1.67,-0.5
A,4,1,9,12,11,10,2,3
A,5,6,7,8
ET,1,13,4,,2     !定义单元类型设置自由度
LSEL,S,LENGTH,,2
LESIZE,ALL,,,30
LSEL,S,LENGTH,,4
LESIZE,ALL,,,60
LSEL,S,LENGTH,,1.33
LESIZE,ALL,,,20
LSEL,S,LINE,,3,5,1
LESIZE,ALL,,,15


！定义材料特性
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.11e11
MPDATA,PRXY,1,,0.33
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,1,REFT,,,  
MPDATA,ALPX,1,,5e-5
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,1,,66.6  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,1.9e11
MPDATA,PRXY,2,,0.33
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,2,REFT,,,  
MPDATA,ALPX,2,,9e-5
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,2,,230

！划分网格
TYPE,   1   
MAT,       1
AMESH,1
ALLS
TYPE,1
MAT,2
AMESH,2


!创建接触对
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
MP,MU,1,
MAT,1   
MP,EMIS,1,  
R,3
REAL,3  
ET,2,169
ET,3,172
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,10000,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,0   
KEYOPT,3,11,0   
KEYOPT,3,12,5   
KEYOPT,3,2,0
KEYOPT,3,1,1
LSEL,S,,,11
CM,_TARGET,LINE
TYPE,2  
NSLL,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
LSEL,S,,,4  
CM,_CONTACT,LINE
TYPE,3  
NSLL,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
CMDEL,_TARGET   
CMDEL,_CONTACT  

！设置边界条件
ALLS
/SOLU
ANTYPE,0
NLGEOM,ON
DL,9,,ALL
SFL,8,PRES,1000
BFL,8,TEMP,800
BFL,9,TEMP,293
ALLS
SOLVE