你选择的什么接触单元?
我刚按照书上做过一个类似的例子(ANSYS机械工程应用25例,电子工业出版社——第13例 平行圆柱体承受法向载荷时的接触应力分析),书上选用的contac49,用于节点对面的接触。
我换成面面接触,在分析时提示有刚性位移。然后加了一定的初始渗透,可以分析出结果。
书上提供的命令流如下,可以参考:
/CLEAR
/PREP7
ET,1,PLANE42
ET,2,SOLID45
ET,3,CONTAC49
MP,EX,1,2E11
MP,PRXY,1,0.3
CYL4,0,0,0.1,0,0,90
CYL4,0,0.15,0.05,-90,0,0
WPOFF,0.01
WPROT,0,0,90
ASBW,ALL
AATT,1,1,1
ESIZE,0.0075
MSHAPE,0
MSHKEY,1
AMESH,ALL
KREFINE,2,4,2,2
TYPE,2
EXTOPT,ESIZE,10
EXTOPT,ACLEAR,1
VEXT,ALL,,,,,0.1
FINISH
/SOLU
DA,16,UY
DA,8,UY
DA,15,UX
DA,20,UX
DA,1,UZ
DA,9,UZ
DA,13,UZ
DA,17,UZ
ASEL,S,,,11,19,8
NSLA,S,1
SFA,ALL,1,PRES,2E5
CP,1,UY,ALL
ALLSEL,ALL
FINISH
/PREP7
R,1,2E10
ASEL,S,,,14
NSLA,S,1
CM,SNAPINS,NODE
ASEL,S,,,18
NSLA,S,1
CM,SNAPPRG,NODE
ASEL,ALL
NSEL,ALL
TYPE,3
REAL,1
GCGEN,SNAPINS,SNAPPRG
GCGEN,SNAPPRG,SNAPINS
FINISH
/SOLU
PRED,ON
AUTOTS,ON
!NLGEOM,ON
NSUBST,6,100,4
OUTRES,ALL,ALL
SOLVE
FINISH
/POST1
SET,LAST
ESEL,S,TYPE,,3
ETABLE,FN,SMISC,1
ETABLE,ST,NMISC,1
ETABLE,GAP,NMISC,3
ESORT,ETAB,FN,1
PRETAB,FN,ST,GAP
FINISH |