*HEADING 
CONCENTRIC RINGS WITH INITIAL OVERCLOSURE [CPS4R]
TEST OF 2D FINITE SLIDING PENALTY CONTACT WITH TABULAR SOFTENED BEHAVIOR
**
** The density has been modified from 8.0100E-06 to 8.0100E-03,
** such that it is more stable and faster. 
** 
*NODE, NSET=ODB_NSET
999,0.,0.
1,3.5,0.0
17,3.5,0.0
18,4.7,0.0
34,4.7,0.0
35,4.6,0.0
51,4.6,0.0
52,6.1,0.0
68,6.1,0.0
*NGEN,LINE=C,NSET=RAD_1
1,17,1,999,,,,0.0,0.0,1.0
*NGEN,LINE=C,NSET=RAD_2
18,34,1,999,,,,0.0,0.0,1.0
*NGEN,LINE=C,NSET=RAD_3
35,51,1,999,,,,0.0,0.0,1.0
*NGEN,LINE=C,NSET=RAD_4
52,68,1,999,,,,0.0,0.0,1.0
*NSET,NSET=INNER,GEN
1,17,1
*ELEMENT,TYPE=CPS4R
1,1,18,19,2
16,16,33,18,1
17,35,52,53,36
32,50,67,52,35
*ELGEN,ELSET=INELE
1,15,1,1,
*ELGEN,ELSET=OUTELE
17,15,1,1
*ELSET,ELSET=INELE,GEN
1,16,1
*ELSET,ELSET=OUTELE,GEN
17,32,1
*SOLID SECTION, ELSET=INELE, MATERIAL=MAT_1
1.0, 
*SOLID SECTION, ELSET=OUTELE, MATERIAL=STEEL
1.0, 
*MATERIAL, NAME=STEEL
*DENSITY
8.0100E-03, 
*ELASTIC, TYPE=ISO
206800.,.32
*MATERIAL, NAME=MAT_1
*DENSITY
8.0100E-03, 
*ELASTIC,TYPE=ISOTROPIC
206800.,.32
*ELSET,ELSET=ODB_ESET
INELE, OUTELE
*NSET, NSET=ODB_NSET
ODB_NSET, RAD_1, RAD_2, RAD_3
RAD_4, INNER
*NSET,NSET=NODERES,ELSET=INELE
*NSET,NSET=NODERES,ELSET=OUTELE
*SURFACE,NAME=IN_SURF
INELE,S2
*SURFACE,NAME=OUT_SURF
OUTELE,S4
**
**
*STEP
*DYNAMIC,EXPLICIT
,1.0E-2
*BULK VISCOSITY
0.02, 0.4
*CONTACT PAIR,INTERACTION=NOFRIC,MECHANICAL CONSTRAINT=PENALTY
OUT_SURF,IN_SURF
*SURFACE INTERACTION,NAME=NOFRIC,PAD THICKNESS=-9.80785E-02
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=TABULAR
0.0,-9.80785E-02
1.0E3,0.0
*BOUNDARY
INNER,1,2
52,2,2
60,2,2
**
** Output requests for postprocessing
**
*OUTPUT,FIELD,VAR=PRESELECT,NUM=2
**
** Output requests for qa testing
**
*NSET,NSET=QA_TESTN
7,10,1
58,61,1
24,27,1
41,44,1
*OUTPUT,FIELD,NUMBER=1
*NODE OUTPUT,NSET=QA_TESTN
U,
*OUTPUT,HISTORY,TIME INTERVAL=0.005
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP