*HEADING
SOFTENED CONTACT - EXPONENTIAL - (GAPUNI)
INPUT DECK: .INP
********************
*SYSTEM
.00667,0.0,0.0
*NODE
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=CAX4H, 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
*ELEMENT,TYPE=GAPUNI,ELSET=CONTACT
701,2001,1401
702,2101,1403
703,2201,1405
704,2301,1407
705,2401,1409
706,2501,1411
707,2601,1413
708,2701,1415
709,2801,1417
710,2901,1419
711,3001,1421
712,3101,1423
713,3201,1425
*SOLID SECTION, ELSET=PREFORM, MATERIAL=RUBBER
*MATERIAL, NAME=RUBBER
*HYPERELASTIC,POLYNOMIAL,N=1
80.,20.
********************
*NODE
1901,-0.00667, 0.03175
1905,-0.00667, 0.0500
3301, 0.08667, 0.03175
3305, 0.08667, 0.0500
*NGEN,NSET=NL
1901,1905
*NGEN,NSET=NR
3301,3305
*NFILL,NSET=NDIE
NL,NR,14,100
*NSET,NSET=DIETOP,GEN
1905,3305,100
*NSET,NSET=DIEBOT,GEN
1901,3301,100
*NSET,NSET=DIESIDE,GEN
1901,1905,1
*ELEMENT,TYPE=CAX4,ELSET=EDIE
4997,1901,2001,2002,1902
*ELGEN,ELSET=EDIE
4997,4,1,1,14,100,4
*ELSET,ELSET=RFACE1,GEN
4997,5049,4
*ELSET,ELSET=RFACE2,GEN
5049,5052,1
*SOLID SECTION,ELSET=EDIE,MATERIAL=DIEMAT
*MATERIAL,NAME=DIEMAT
*ELASTIC
30.E6,0.3
********************
*GAP,ELSET=CONTACT
0.0,0.,-1.,0.
*FRICTION
.1, 
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=EXPONENTIAL
0.001,0.0445
*BOUNDARY
DIETOP, 1
DIESIDE, 1
AXIS, 1
BOTTOM, 2
TOP,3
DIEBOT,3
********************
*STEP, INC=200, NLGEOM
*STATIC
1.E-2, 1.0,  , 0.1
*BOUNDARY
DIETOP,2, ,-0.01223
*PRINT, CONTACT=YES
*NODE FILE, NSET=TOP,FREQUENCY=5
U, RF
*EL FILE,ELSET=CONTACT,FREQUENCY=5
S, 
E, 
*EL PRINT, ELSET=PREFORM
S, MISES
E, PEEQ
*EL PRINT,ELSET=CONTACT
S,
E, 
*NODE PRINT
U, RF
*END STEP