!************对实体划分网格 **************************** !*
LESIZE,17, , ,15, , , , ,1 LESIZE,19, , ,15, , , , ,1 !*
LESIZE,18, , ,2, , , , ,1 LESIZE,20, , ,2, , , , ,1 !*
LESIZE,22, , ,20, , , , ,1 !*
LESIZE,5, , ,10, , , , ,1 LESIZE,7, , ,10, , , , ,1 !*
LESIZE,6, , ,8, , , , ,1 LESIZE,8, , ,8, , , , ,1 !*
LESIZE,10, , ,3, , , , ,1 !*
VSWEEP,ALL
!*************网格划分完毕********************************** ! !
/COM, CONTACT PAIR CREATION - START MP,MU,1,0.2 MAT,1 R,3 REAL,3 ET,2,170 ET,3,174
R,3,,,0.2,0.2,0.9,0 0.2;
! 设置过盈量为 f; RMORE,,,1.0E20,f,1.0,0 RMORE,0.0,0,1.0,,1.0,0 RMORE,0,1.0,1.0,0.0,,1.0 KEYOPT,3,4,0
! Icont 系数设置为 KEYOPT,3,5,0 NROPT,UNSYM KEYOPT,3,7,0 KEYOPT,3,8,0 !KEYOPT(9)=4;
! 不考虑初始几何位置造成的过盈或者间隙， ! 只考虑 CNOF 参数设置的值，即过盈量。 KEYOPT,3,9,4 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 ESLL,U
ESEL,U,ENAME,,188,189 ESURF
CMSEL,S,_ELEMCM ! Generate the contact surface ASEL,S,,,9
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
/COM, CONTACT PAIR CREATION - END FINISH /SOL
FLST,2,4,5,ORDE,4 FITEM,2,5
FITEM,2,-6 FITEM,2,11 FITEM,2,-12
DA,P51X,SYMM FLST,2,1,5,ORDE,1 FITEM,2,3 /GO
DA,P51X,ALL, !*
ANTYPE,0
NLGEOM,1 NSUBST,1,0,0 AUTOTS,0 TIME,100
/STATUS,SOLU SOLVE FINISH /POST1
SET,1,LAST,1,
PLNSOL,s,eqv stress;
esel,s,ename,,174
PLNSOL, CONT,PRES, 0,1.0pressure;
SAVE
! 设置对称约束； ! 必须打开大变形效应； ! 查看单元的Von Mises
在contact 单元上查看contact
! 计算结果：
说明：修改 f 值即可修改过盈量。如果过盈量设置的过大，接触部件可能会产生塑性变形，如果要考虑塑性变形，则应该修改材料属性定义为塑性材料。