为什么没有接触反应呢？

lhxing

我的问题是：两个立方体之间夹有一个圆柱体。下立方体的下表面施加固定约束，上立方体的上表面施加100N的压力，两个侧面固定约束，圆柱体与两个立方体之间建立接触对，立方体为目标面。求解后只有一个立方体上有应力显示，而圆柱体和下立方体无应力显示。不知问题出在哪儿，真心请教斑竹和高手。以下是我建立的模型和结果图，谢谢！
/BATCH  
/COM,ANSYS RELEASE  9.0    UP20041104       10:40:36    07/14/2006
/PREP7  
!*  
ET,1,SOLID45
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.1e5  
MPDATA,PRXY,1,,0.3  
BLOCK,0,30,0,15,0,25,   
BLOCK,0,30,25,40,0,25,  
CYL4,15,20,5, , , ,25   
ESIZE,0,20,
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
VPLOT   
ESIZE,0,10,
CM,_Y,VOLU  
VSEL, , , ,       3
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
FLST,2,1,5,ORDE,1   
FITEM,2,3   
!*  
/GO
DA,P51X,ALL,0   
FLST,2,2,5,ORDE,2   
FITEM,2,11  
FITEM,2,-12
!*  
/GO
DA,P51X,ALL,0   
FLST,2,1,5,ORDE,1   
FITEM,2,10  
/GO
!*  
SFA,P51X,1,PRES,100
!*  
!*  
/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.1
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3  
ET,2,170
ET,3,174
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,4,0
KEYOPT,3,5,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,2   
KEYOPT,3,11,0   
KEYOPT,3,12,0   
KEYOPT,3,2,0
KEYOPT,2,5,0
! Generate the target surface   
ASEL,S,,,4  
CM,_TARGET,AREA
TYPE,2  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,16
CM,_CONTACT,AREA
TYPE,3  
NSLA,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  
/COM, CONTACT PAIR CREATION - END   
!*  
!*  
/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.1
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,4
REAL,4  
ET,4,170
ET,5,174
R,4,,,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,5,4,0
KEYOPT,5,5,0
KEYOPT,5,7,0
KEYOPT,5,8,0
KEYOPT,5,9,0
KEYOPT,5,10,2   
KEYOPT,5,11,0   
KEYOPT,5,12,0   
KEYOPT,5,2,0
KEYOPT,4,5,0
! Generate the target surface   
ASEL,S,,,9  
CM,_TARGET,AREA
TYPE,4  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,15
CM,_CONTACT,AREA
TYPE,5  
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,4  
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
SAVE
FINISH  
/SOL
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
!*  
/EFACET,1   
PLNSOL, S,EQV, 0,1.0

lhxing

结果图

zlg

调整一下圆柱与下立方体的接触方向

sysh320

肯定是你的接触单元或接触目标单元有问题