有关初始应力，对称边界和接触的实例

cam_1980

以前照书上做的一个实例，现在把log文件整理贴出来。通过这个实例大家可以了解对称边界，接触以及过盈配合产生初始应力的一些知识。
/PREP7
ET,1,SOLID185
MP,EX,1,2.1E5
MP,PRXY,1,0.3
CYL4,0,0,34,0,100,90,25
CYL4,0,0,25,0,35,90,150           !注意这里的模型包含了一个毫米的过盈量
VGEN, ,2, , , , ,-10, , ,1
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.1,0.1,,
NROPT,UNSYM
ASEL,S,,,4
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF,ALL
!* Generate the contact surface
ASEL,S,,,9
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF,ALL
ALLSEL,ALL
EPLOT
FINISH

/SOLU
DA,5,SYMM         !定义面5为对称边界
DA,6,SYMM         !定义面6为对称边界
DA,11,SYMM        !定义面11为对称边界
DA,12,SYMM        !定义面12为对称边界
DA,3,ALL,         !约束圆盘的外圆柱面所有自由度
ANTYPE,0
NLGEOM,1          !开启大变形
AUTOTS,0
TIME,100
SOLVE              !第一个载荷步求解
NSUBST,150,10000,10
OUTRES,ALL,ALL
AUTOTS,1
TIME,250
NSEL,S,LOC,Z,140
D,ALL,UZ,40
ALLSEL,ALL
SOLVE              !第二个载荷步求解        

/EXPAND,4,POLAR,HALF,,90                !将1/4模型还原成整模型
/REPLOT