*HEADING
EXPLICIT DYNAMICS TEST OF SOFTENED KINEMATIC CONTACT WITH TABULAR
PRESSURE-OVERCLOSURE DATA AND DEFORMABLE SURFACES
********************
*SYSTEM
.02197,0.0,0.0
*NODE, NSET=ODB_NSET
1,     0.0,     0.0
25,    0.08,    0.0
1401,  0.0,     0.03175
1425,  0.08,    0.03175
*NGEN, NSET=BOTTOM
1, 25
*NGEN, NSET=TOP
1401, 1425
*NFILL, NSET=ALL
BOTTOM, TOP, 14, 100
*NSET, NSET=AXIS, GENERATE
1, 1401, 100 
*ELEMENT, TYPE=CAX4R, ELSET=PREFORM
1, 1, 3, 203, 201
*ELGEN, ELSET=PREFORM
1, 12, 2, 1, 7, 200, 100
*ELSET,ELSET=DTOP,GEN
601,612,1
*ELSET,ELSET=DSIDE,GEN
12,612,100
*SOLID SECTION, ELSET=PREFORM, MATERIAL=RUBBER
*MATERIAL, NAME=RUBBER
*HYPERELASTIC,POLYNOMIAL,N=1
80.,20.
*DENSITY
100., 
********************
*NODE, NSET=ODB_NSET
1901,-0.02197, 0.03175
1905,-0.02197, 0.0500
2301, 0.0659, 0.03175
2305, 0.0659, 0.0500
2401, 0.0711, 0.0332
2402, 0.0762, 0.0363
2403, 0.0787, 0.0407
2404, 0.0807, 0.0437
2405, 0.0793, 0.0500
2504, 0.0853, 0.0452
2505, 0.0853, 0.0500
2804, 0.1030, 0.0452
2805, 0.1030, 0.0500
2904, 0.1076, 0.0437
2905, 0.1090, 0.0500
3005, 0.1142, 0.0500
3004, 0.1142, 0.0452
2903, 0.1090, 0.0363
3003, 0.1142, 0.0363
2902, 0.1090, 0.03175
3002, 0.1142, 0.03175
2901, 0.1090, 0.0284
2951, 0.1116, 0.0284
2900, 0.1090, 0.0251
2950, 0.1116, 0.0251
3001, 0.1142, 0.0251
3101, 0.1194, 0.0251
3102, 0.1194, 0.03175
*NGEN,NSET=NL
1901,1905
*NGEN,NSET=NR
2301,2305
*NFILL,NSET=NDIE
NL,NR,4,100
*NGEN, NSET=ODB_NSET
2504,2804,100
2505,2805,100
*NSET,NSET=DIETOP,GEN
1905,3005,100
*NSET,NSET=DIESIDE,GEN
1901,1905,1
*ELEMENT,TYPE=CAX4R,ELSET=EDIE
4997,1901,2001,2002,1902
5017,2404,2504,2505,2405
5023,2903,3003,3004,2904
5024,2902,3002,3003,2903
5025,2901,2951,3002,2902
5026,2900,2950,2951,2901
5027,2950,3001,3002,2951
5028,3001,3101,3102,3002
*ELGEN,ELSET=EDIE
4997,4,1,1,5,100,4
5017,6,100,1
*ELSET,ELSET=RFACE1,GEN
4997,5013,4
*ELSET,ELSET=RFACE2,GEN
5013,5015,1
*ELSET,ELSET=RFACE3,GEN
5017,5021,1
*ELSET,ELSET=RFACE4,GEN
5023,5026,1
*ELSET,ELSET=RFACE5,GEN
5026,5028,1
*SOLID SECTION,ELSET=EDIE,MATERIAL=DIEMAT
*MATERIAL,NAME=DIEMAT
*ELASTIC
30.E+6,0.3
*DENSITY
1000., 
*ELSET,ELSET=ODB_ESET
PREFORM, EDIE
**
**
*NSET, NSET=ODB_NSET
ODB_NSET, BOTTOM, TOP, ALL
AXIS, NL, NR, NDIE
DIETOP, DIESIDE
*SURFACE, NAME=DIE
RFACE1, S1
RFACE2, S2
RFACE3, S1
RFACE4, S4
RFACE5, S1
*SURFACE, NAME=BLANK
DTOP,
DSIDE,
*NSET,NSET=NODERES,ELSET=DTOP
*NSET,NSET=NODERES,ELSET=DSIDE
**
*STEP
*DYNAMIC,EXPLICIT
,1.0E-2
*CONTACT PAIR, INTERACTION=SLIDE,WEIGHT=0.0, CPSET=CP1
 BLANK, DIE
*SURFACE INTERACTION, NAME=SLIDE
*FRICTION
 0.1, 
*SURFACE BEHAVIOR, PRESSURE-OVERCLOSURE=TABULAR
 0.,-.001
 25.,.001
 23750.,.02
*BOUNDARY
 DIETOP, 1
 DIESIDE, 1
 AXIS, 1
 BOTTOM, 2
*AMPLITUDE,NAME=RAMP
 0.0,0.0,1.0E-2,1.0
*BOUNDARY, AMPLITUDE=RAMP
 DIETOP,2, ,-0.01223
**
*OUTPUT,FIELD,NUMBER INTERVAL=2,TIME MARKS=YES
*ELEMENT OUTPUT
 S,
*NODE OUTPUT
 U,
*CONTACT OUTPUT
 CSTRESS,
**
*OUTPUT,HISTORY,TIME INTERVAL=1.E-5
*CONTACT OUTPUT,CPSET=CP1
 CAREA
*NODE OUTPUT,NSET=NODERES
 U,RF
*ENERGY OUTPUT
 ALLKE,ALLSE,ALLIE,ALLAE,ETOTAL
**
*FILE OUTPUT,NUMBER INTERVAL=2,TIME MARKS=YES
*EL FILE
 S,
*NODE FILE
 U,
**
*END STEP