*HEADING
	2-D ELEMENT IN CONTACT WITH ANOTHER DEFORMABLE 2-D ELEMENT 
        REPRESENTED BY A MATRIX. LARGE SLIDING CONTACT.
*NODE
1,
2,1.,
3,1.,1.,
4,,1.
11,-1.,1.01
12,2.,1.01
13,2.,3.
14,-1.,3.
15,5.,1.01
16,5.,3.
*ELEMENT,TYPE=CPE4,ELSET=TOP
11,11,12,13,14
12,12,15,16,13
*SOLID SECTION,ELSET=TOP,MATERIAL=MAT
***ELEMENT,TYPE=CPE4,ELSET=BOT
**1,1,2,3,4
***SOLID SECTION,ELSET=BOT,MATERIAL=MAT
*MATERIAL,NAME=MAT
*ELASTIC
3.E7,0.0
*surface,NAME=TOPS
TOP,S1
*MATRIX INPUT,NAME=MATRIX1,TYPE=SYMMETRIC,INPUT=contact_stiff.inp
*MATRIX ASSEMBLE,STIFFNESS=MATRIX1
*surface,type=node,NAME=SLAVE
3,1.0
4,1.0
*SURFACE INTERACTION,NAME=ISLS
*FRICTION
0.1
*CONTACT PAIR,INTERACTION=ISLS
SLAVE,TOPS
*BOUNDARY
13,1,2
14,1,2
16,1,2
************************
*STEP,NLGEOM,unsymm=yes
*STATIC
0.1,1.0
*BOUNDARY
1,1,1,0.
1,2,2,0.2
2,1,1,0.
2,2,2,0.2
***PRINT,CONTACT=YES
***NODE PRINT
**RF,
**U,
***ELPRINT
**S,
**E
*OUTPUT, FIELD, FREQ=100
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S, 
E,
*OUTPUT,HISTORY
*ENERGY OUTPUT 
*ENDSTEP
************************
*STEP
*STATIC
0.1,1.0
*BOUNDARY
1,1,1,3.0
2,1,1,3.0
***PRINT,CONTACT=YES
***NODE PRINT
**RF,
**U,
***ELPRINT
**S,
**E
*OUTPUT, FIELD, FREQ=100
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S, 
E, 
*OUTPUT,HISTORY
*ENERGY OUTPUT 
*ENDSTEP