摸索了好久，热-结构-接触的问题，有命令流，高手指点一下了，

huntfora

FINI
/CLEAR

*AFUN,RAD
/UNITS,SI                     
/PREP7  
ET,1,13
KEYOPT,1,1,4
MP,EX,1,10E9
MP,PRXY,,0.3
MP,KXX,1,339

Et,2,169
ET,3,171
KEYOPT,3,1,1
KEYOPT,3,12,5   
!R,3,0,0,0.1, 0.1,0,0
R,3,,,
RMORE,,,, ,,

!MODELING
AA=1E-7
BB=1E-6
pi=3.14159265
*do,i,0,100,1                 
*SET,x,AA* i                    
*SET,y,BB*sin(i*pi/20)               
k,i+1,x,y                       
*enddo  

allsel
bspline,all           !样条插值，生成一条光滑的曲线
!以上生成一条光滑的正弦曲线

K,102,0,BB
K,103,100*AA,BB
K,104,0,10*BB
K,105,100*AA,10*BB
L,102,103
L,102,104
L,103,105
L,104,105

K,106,0,-10*BB
K,107,100*AA,-10*BB

L,1,106
L,101,107
L,106,107
LPLOT

AL,1,6,7,8
AL,2,3,4,5
APLOT

ESIZE,0.5*BB
MSHKEY,0
AMESH,ALL            !划分网格

LSEL,S,,,2   
CM,_TARGET,AREA  
TYPE,2   
NSLL,S,1
ESLN,S,0
ESURF,ALL   

LSEL,S,,,1  
CM,_CONTACT,AREA
TYPE,3   
NSLL,S,1
ESLN,S,0
ESURF,ALL   

CMDEL,_TARGET   
CMDEL,_CONTACT   
ALLSEL,ALL   
! 以上生成接触对   

ALLSEL
!耦合结合处的温度自由度，从而使得温度连续
NSEL,S,LOC,Y,BB,0.95*BB    !我也不知这样耦合对不对
CP,1,TEMP,ALL

FINISH

/SOLU
ALLSEL
!DL,5, ,UY,0
DA,2,UY,0
DA,2,UX,0
SFL,8,PRES,1e3,          !设为1e9,不收敛
！当调整向上作用的压力时，von mises却没有成比例变化

DL,8, ,TEMP,1000,0  

SFL,5,CONV,30000, ,25,           !上表面施加对流边界条件


TIME,200
NSUBST,100,1000,100
kbc,1
autots,on               !打开自动时间步长
ALLSEL

CNVTOL,F,1,0.005,2, ,        !定义收敛准则
CNVTOL,HEAT,1,0.005,2, ,   

SOLVE

/POST1
SET,LAST
PLNSOL,TEMP

PLNSOL, S,EQV, 0,1.0
当调整向上作用的压力时，von mises却没有成比例变化，理论上应该是成比例变化的啊
而温度场也没有变化。
版主指点指点啊
另外，我发现VM63里只用了一个接触单元CONTA178，没有成对使用
难道不建立接触对也能进行接触分析吗？

以上是网格。

上面是von mises, 由于上传要小于512k,故截掉了一些

上面是温度场

huntfora

自已顶下，哪们有好点的想法啊？