接触问题分析

中升肖慧丽

遇到接触分析的问题：玻璃与圆柱橡胶接触，命令流是书上的，但就是运行不出来，哪位大侠帮忙看一下问题在哪里？
FINISH
/CLEAR
/FILNAME,pt_to_surf_contact,0   
/TITLE,example of point to surface contact analysis

/PREP7   
ET,1,hyper56  
MP,EX,1,2.8e6  
MP,PRXY,1,0.498  
TB,MOON,1,1,9,  
TBTEMP,0
TBDATA,,0.293e6,0.177e6,,,,
TBDATA,,,,,,,   
!创建几何模型  
PCIRC,0.2,,0,90
RECTNG,,0.4,0.2,0.21,  
!
LESIZE,2, , ,16,5, , , ,1
LESIZE,3, , ,16,5, , , ,1  
MSHAPE,0,2D
MSHKEY,1
AMESH,1   

!创建接触对
/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,2
MAT,1   
R,3
REAL,3  
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,2   
R,3,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
LSEL,S,,,4  
CM,_TARGET,LINE
TYPE,2  
LATT,-1,3,2,-1  
TYPE,2  
LMESH,ALL   
! Generate the contact surface  
LSEL,S,,,1  
CM,_CONTACT,LINE
TYPE,3  
NSLL,S,1
ESLN,S,0
ESURF   
*SET,_REALID,3  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
LSEL,S,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  
LSEL,S,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  
/COM, CONTACT PAIR CREATION - END   
/MREP,EPLOT

/SOL
ANTYPE,0
NLGEOM,1
OUTRES,ERASE
OUTRES,ALL,1

DL,3, ,SYMM
DL,2, ,UY,0.1  
NSUBST,20,10000,10
TIME,100  
SOLVE   
FINISH

koalabear

应该不会有问题，我用11版本执行，一直到求解没有问题，估计是你的版本问题，或者计算机硬件问题。