*HEADING
SOFTENED CONTACT - EXPONENTIAL - (ISL)
INPUT DECK: .INP
********************
*SYSTEM
.02197,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=ISL21A,ELSET=CONTACT
701,1401,1403
702,1403,1405
703,1405,1407
704,1407,1409
705,1409,1411
706,1411,1413
707,1413,1415
708,1415,1417
709,1417,1419
710,1419,1421
711,1421,1423
712,1423,1425
713,1425,1225
714,1225,1025
715,1025,825
716,825,625
717,625,425
718,425,225
719,225,25
*SOLID SECTION, ELSET=PREFORM, MATERIAL=RUBBER
*MATERIAL, NAME=RUBBER
*HYPERELASTIC,POLYNOMIAL,N=1
80.,20.
********************
*NODE
1901,-.02197, 0.03175
1905,-.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
2504,2804,100
2505,2805,100
*NSET,NSET=DIETOP,GEN
1905,3005,100
*NSET,NSET=DIESIDE,GEN
1901,1905,1
*ELEMENT,TYPE=CAX4,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.E6,0.3
********************
*SLIDE LINE,ELSET=CONTACT,TYPE=LINEAR
2900,2901,2902,2903,2904,2804,2704,2604,2504,2404,2403,2402,2401,
2301,2201,2101,2001,1901
*INTERFACE, ELSET=CONTACT
*FRICTION
.1, 
*SURFACE BEHAVIOR, PRESSURE-OVERCLOSURE=EXPONENTIAL
0.001,12.5
*BOUNDARY
DIETOP, 1
DIESIDE, 1
AXIS, 1
BOTTOM, 2
********************
*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
*NODE PRINT
U, RF
*END STEP