*HEADING
CONTACT PROBLEM
**
** DEFINE NODE MESH
**
*NODE,NSET=NODE1
*NODE,NSET=NODE2
*NODE,NSET=NODE3
**TWO STEP option is need when using quadratic elements
*NFILL,NSET=TOP,BIAS=1.0,TWO STEP
NODE1,NODE2,60,1
*NFILL,NSET=TOP,BIAS=0.8,TWO STEP
NODE2,NODE3,?,?
*NCOPY,OLD SET=TOP,NEW SET=MIDDLE,CHANGE NUMBER=3030,SHIFT
*NCOPY,OLD SET=TOP,NEW SET=BOTTOM,CHANGE NUMBER=6060,SHIFT
*NFILL,NSET=NALL,BIAS=1.0,TWO STEP
TOP,MIDDLE,30,101
*NFILL,NSET=NALL,BIAS=0.8,TWO STEP
MIDDLE,BOTTOM,30,101
*NSET,NSET=SYMM,GENERATE
?,?,?
**
** DEFINE ELEMENTS
**
*ELEMENT,TYPE=CAX8R,ELSET=ELALL
*ELGEN,ELSET=ELALL
*ELSET,ELSET=TOPLINE,GENERATE
** DEFINE RIGID SPHERE WITH D=3MM
**
*NODE,NSET=BALL
*SURFACE,TYPE=SEGMENTS,NAME=BSURF
START,0.0015,0.0015
CIRCL,0.0,0.0,0.0,0.0015
*RIGID BODY,REF NODE=99999, ANALYTICAL SURFACE=BSURF
*SURFACE,TYPE=ELEMENT,NAME=ASURF
?,?
**
** DEFINE CONTACT PAIRS
**
*CONTACT PAIR,INTERACTION=WITHFRIC
?,?
*SURFACE INTERACTION,NAME=WITHFRIC
*FRICTION
0.1
*BOUNDARY
SYMM,1,1,0.0
BOTTOM,2,2,0.0
BALL,1,1,0.0
BALL,6,6,0.0
**
** Define Ideal Plastic Material
**
*SOLID SECTION,ELSET=ELALL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
?,?
*PLASTIC
200.0E+06,0.0
200.0E+06,0.5
**
** Apply Indentation
**
*STEP,NLGEOM=??,INC=1000
*STATIC
0.001,1.0,1E-9,1.0
*BOUNDARY
99999,2,2,-0.0001
*NODE PRINT,NSET=BALL,freq=1000
RF2
*NODE PRINT,NSET=TOP,FREQ=1000
COORD
*EL PRINT,ELSET=TOPLINE,FREQ=1000
S
*CONTACT PRINT,SLAVE=ASURF,FREQ=1000
*OUTPUT, FIELD, VARIABLE=PRESELECT
*OUTPUT, HISTORY, VARIABLE=PRESELECT
*END STEP