ansys接触问题小案例深度解析与技巧分享

BATCH
/COM,ANSYS RELEASE 5.6 UP19991022 13:25:41 08/08/2002
!*
/PREP7
ET,1,SOLID95
!*
!*
UIMP,1,EX, , ,212000,
UIMP,1,NUXY, , ,0.3,
UIMP,1,ALPX, , , ,
UIMP,1,REFT, , , ,
UIMP,1,MU, , ,0.2 ,
UIMP,1,DAMP, , , ,
UIMP,1,DENS, , , ,
UIMP,1,VISC, , , ,
UIMP,1,SONC, , , ,
UIMP,1,PERX, , , ,
*SET,u,0.5
!u是一个参数，这段程序是为了确定u为何值时，接触应力最小，上面语句中u是可变的
!定义材料
K, ,0,14,,
K, ,12,14-u,,
K, ,13.5,14-u,,
K, ,13.5,7.5,,
K, ,13,7,,
K, ,-13,7,,
K, ,-13.5,7.5,,
K, ,-13.5,14-u,,
K, ,-12,14-u,,
LARC, 9, 2, 1
LSTR, 2, 3
LSTR, 3, 4
LSTR, 4, 5
LSTR, 5, 6
LSTR, 6, 7
LSTR, 7, 8
LSTR, 8, 9
!*
LFILLT,2,3,0.3, ,
!*
LFILLT,8,7,0.3, ,
FLST,2,10,4
FITEM,2,1
FITEM,2,2
FITEM,2,9
FITEM,2,3
FITEM,2,4
FITEM,2,5
FITEM,2,6
FITEM,2,7
FITEM,2,10
FITEM,2,8
AL,P51X
K, ,-10,,,
K, ,10,,,
FLST,2,1,5,ORDE,1
FITEM,2,1
FLST,8,2,3
FITEM,8,14
FITEM,8,15
VROTAT,P51X, , , , , ,P51X, ,180, ,
WPSTYLE,,,,,,,,1
wpstyle,0.05,1,-10,10,0.003,0,0,,5
wpro,,,90.000000
wpoff,0,0,12
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VSBW,P51X
wpoff,0,0,1
FLST,2,2,6,ORDE,2
FITEM,2,3
FITEM,2,-4
VSBW,P51X
wpoff,0,0,-25
FLST,2,2,6,ORDE,2
FITEM,2,5
FITEM,2,-6
VSBW,P51X
wpoff,0,0,-1
FLST,2,2,6,ORDE,2
FITEM,2,9
FITEM,2,-10
VSBW,P51X
WPCSYS,-1,0
FLST,5,14,4,ORDE,6
FITEM,5,41
FITEM,5,-50
FITEM,5,57
FITEM,5,65
FITEM,5,74
FITEM,5,81
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,12,3, , , ,1
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,14,4,ORDE,6
FITEM,5,21
FITEM,5,-30
FITEM,5,36
FITEM,5,53
FITEM,5,71
FITEM,5,78
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,12,1/3, , , ,1
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,4,4,ORDE,4
FITEM,5,1
FITEM,5,12
FITEM,5,31
FITEM,5,73
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,20, , , , ,1
CMDELE,_Y
CMDELE,_Y1
!*
ESIZE,1,0,
FLST,5,10,6,ORDE,4
FITEM,5,1
FITEM,5,-8
FITEM,5,11
FITEM,5,-12
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
!*
VCLEAR,_Y1
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!滚轮建模及分网结束
wpro,,,90.000000
CYLIND,7, ,-19.65,19.65,0,360,
wpro,,,90.000000
FLST,2,2,5,ORDE,2
FITEM,2,46
FITEM,2,-47
ASBW,P51X
wpro,,,-90.000000
wpoff,0,0,14
FLST,2,2,5,ORDE,2
FITEM,2,49
FITEM,2,59
ASBW,P51X
wpoff,0,0,-28
FLST,2,2,5,ORDE,2
FITEM,2,46
FITEM,2,-47
ASBW,P51X
WPCSYS,-1,0
!滚轮销切面结束
ESIZE,1.5,0,
CM,_Y,VOLU
VSEL, , , , 9
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!滚轮销划网结束
BLOCK,-20,20,-15,15,14,30,
FLST,5,2,4,ORDE,2
FITEM,5,114
FITEM,5,-115
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,12,4, , , ,1
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,2,4,ORDE,2
FITEM,5,113
FITEM,5,116
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,12,1/4, , , ,1
CMDELE,_Y
CMDELE,_Y1
!*
ESIZE,2,0,
CM,_Y,VOLU
VSEL, , , , 10
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!立方体建模及划网结束
FLST,5,4,5,ORDE,3
FITEM,5,49
FITEM,5,59
FITEM,5,-61
ASEL,S, , ,P51X
NSLA,S,1
NPLOT
FLST,2,170,1,ORDE,9
FITEM,2,5522
FITEM,2,-5528
FITEM,2,5536
FITEM,2,5545
FITEM,2,-5589
FITEM,2,5797
FITEM,2,-5801
FITEM,2,6036
FITEM,2,-6147
!*
/GO
F,P51X,FZ,46.18
ALLSEL,ALL
!加力为7850N
FLST,2,2,3,ORDE,2
FITEM,2,26
FITEM,2,-27
!*
/GO
DK,P51X, , , ,0,UX,UY, , , , ,
FLST,2,1,5,ORDE,1
FITEM,2,47
!*
/GO
DA,P51X,ALL,
！定义的约束
FINISH
/PREP7
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_LINECM,LINE
CM,_AREACM,AREA
/GSAV,cwz,gsav,,temp
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,2,170
ET,3,174
RMODIF,3,1,,,0.1,0.1,,
RMODIF,3,7,,,1.0e20,0.0,1.0
KEYOPT,3,2,0
KEYOPT,3,6,0
KEYOPT,3,7,0
KEYOPT,3,8,
KEYOPT,3,9,0
KEYOPT,3,11,
KEYOPT,3,12,0
! Generate the target surface
ASEL,S,,,57
ASEL,A,,,58
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,R,0
ESURF,ALL
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,27
ASEL,A,,,44
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,R,0
ESURF,ALL
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,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
!*
!*
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_LINECM,LINE
CM,_AREACM,AREA
/GSAV,cwz,gsav,,temp
MP,MU,1,0.2
MAT,1
R,4
REAL,4
ET,4,170
ET,5,174
RMODIF,4,1,,,0.1,0.1,,
RMODIF,4,7,,,1.0e20,0.0,1.0
KEYOPT,5,2,0
KEYOPT,5,6,0
KEYOPT,5,7,0
KEYOPT,5,8,
KEYOPT,5,9,0
KEYOPT,5,11,
KEYOPT,5,12,0
! Generate the target surface
ASEL,S,,,46
CM,_TARGET,AREA
TYPE,4
NSLA,S,1
ESLN,R,0
ESURF,ALL
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,2
ASEL,A,,,13
CM,_CONTACT,AREA
TYPE,5
NSLA,S,1
ESLN,R,0
ESURF,ALL
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,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
！定义接触
FINISH
/SOLU
EQSLV,PCG,1E-8
NEQIT,50
!定义求解器