ANSYS关于接触分析的一个命令流

lixinfu

关于接触分析的一个命令流，感兴趣的可以自己拿去练习一下，模型建立也许本身不难，只是单元实常数定义之类的比较恼火。
命令流如下：
pressure=33.1
/PREP7   
ET,1,PLANE42                                                
!*                                             ！定义材料属性  
MPTEMP,1,0  
MPDATA,EX,1,,127000
MPDATA,PRXY,1,,0.3  
MPDATA,MU,1,,0                                         
BLC5,0,-0.49,19,0.98
CYL4,0,50.8,50.8,180,0,360
K, ,-4.75,0,,   
K, ,4.75,0,,
K, ,4.75,9.5,,  
K, ,-4.75,9.5,,
LSTR,       8,      11  
LSTR,      10,      11  
LSTR,       9,      10  
lsel,s,loc,y,0.1,9.5
asbl,2,all
asel,s,loc,y,10,50.8
adele,all,,,1
allsel,all
numcmp,area
numcmp,line
LESIZE,2,0.1, , , , , , ,1
LESIZE,1,0.2, , , , , , ,1
mshkey,1                                   !设定为用映射网格划分
amesh,1
lsel,s,loc,y,0.5,9.4
LESIZE,all,0.2, , , , , , ,1
allsel,all
lsel,s,loc,y,9.5,9.5
LESIZE,all,0.2, , , , , , ,1
allsel,all
numcmp,kp
amap,2,5,6,7,8                               !一下是log文件里的命令创造接触对
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.5
MAT,1   
MP,EMIS,1,7.88860905221e-031
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,0,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,3
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,1   
KEYOPT,3,11,0   
KEYOPT,3,12,0   
KEYOPT,3,2,0
! Generate the target surface   
LSEL,S,,,8  
CM,_TARGET,LINE
TYPE,2  
NSLL,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
LSEL,S,,,3  
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  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
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   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT      
FINISH  
/SOL                                                        !定义边界条件和加力
lsel,s,loc,y,0.5,9.4
nsll,s,1
d,all,,,,,,ux,,,,,
allsel,all
lsel,s,loc,x,-9.5,-9.5  
nsll,s,1
d,all,,,,,,ux,uy,,,,
allsel,all
lsel,s,loc,y,-0.98,-0.98
nsll,s,1
dl,1,,symm
allsel,all
SFL,5,PRES,pressure,
slove
。。。