大家讨论一下这个接触问题
浏览：212415 回答：16
我的问题是：两个立方体之间夹有一个圆柱体。下立方体的下表面施加固定约束，上立方体的上表面施加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
/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
继续讨论。如果是整个轴承又该如何做？