*HEADING
 CONTACT PROBLEM - DISPLACEMENT CONTROL, COARSE MESH  
*RESTART,WRITE,FREQUENCY=999
*NODE,NSET=REM
 300, 0.     ,-11.
*NODE,NSET=FIX
 301, 0.     ,-10.
 302, 0.26177,-10.
 303, 0.52336,-10.
 304, 0.78459,-10.
 305, 1.04528,-10.
 306, 1.30526,-10.
 307, 1.56434,-10.
 308, 1.82236,-10.
 309, 2.07912,-10.
*NODE
   1, 0.,0.
   5, 2.,0.
   9, 5.5,0.
  13, 10.,0.
  41, 0.,-3.
  45, 2.,-3.
  49, 5.2,-3.
  81, 0.,-5.5
  85, 2.,-5.5
  89, 4.,-5.5
 121, 0.,-7.5
 129, 2.,-7.5
 133, 3.5,-7.5
 101, 0.,-6.5
 103, 1.,-6.5
 105, 2.,-6.5
 107, 3.75,-6.5
 201, 0.,-10.
*NODE,SYSTEM=C
 33,10.,-15.
 53,10.,-30.
 73,10.,-37.5
 93,10.,-45.
113,10.,-55.
213,10.,-65.
209,10.,-78.
*NGEN,LINE=C
 209,213
 201,209
*NGEN
 1,5,2
 5,9,2
 9,13,2
 41,45,2
 45,49,2
 49,53,2
 81,85
 85,89,2
 89,93,2
 121,129
 129,133
 133,213,40
 131,211,40
 1,41,20
 5,45,20
 9,49,20
 41,81,20
 45,85,20
 49,89,20
 81,121,20
 121,201,20
 123,203,20
 125,205,20
 127,207,20
 129,209,20
 161,163
 163,165
 165,167
 167,169
*NSET,NSET=NY
 1,21,41,61,81,101,121,141,161,181,201
*ELEMENT,TYPE=CPE8
 1, 1,41,45,5,21,43,25,3
 7, 81,121,125,83,101,123,103,82
 8, 83,125,129,85,103,127,105,84
 9, 85,129,133,89,105,131,107,87
 10, 89,133,213,93,107,173,113,91
 11, 121,161,163,123,141,162,143,122
 19, 129,209,211,131,169,210,171,130
 20, 131,211,213,133,171,212,173,132
*ELGEN
 1,3,4,1
 1,2,40,3
 4,3,4,1
 11,4,2,1
 11,2,40,4
 15,4,2,1
*ELSET,ELSET=TOP
 1,2,3,4,5,6,7,8,9,10,11,12,13,14,19,20
*ELSET,ELSET=BOT
 15,16,17,18
*ELSET,ELSET=CYLINDER
TOP,BOT
*ELEMENT,TYPE=R2D2
101,301,302
*ELGEN,ELSET=RIGID
101,8,1,1
*RIGID BODY,ELSET=RIGID,REF NODE=300
*surface,NAME=BOTSUR,TRIM=YES
BOT,
*surface,NAME=RIGSUR
RIGID,SPOS
*CONTACT PAIR,INTERACTION=NO_FRI,SMALL SLIDING
BOTSUR,RIGSUR
*SURFACE INTERACTION,NAME=NO_FRI
1.0,
*SURFACE BEHAVIOR, AUGMENTED LAGRANGE
**
**
*EQUATION
 2, 
 3,2,-1.,1,2,1.
 2, 
 5,2,-1.,1,2,1.
 2, 
 7,2,-1.,1,2,1.
 2, 
 9,2,-1.,1,2,1.
 2, 
11,2,-1.,1,2,1.
 2, 
13,2,-1.,1,2,1.
*MPC
QUADRATIC,82,81,83,85
QUADRATIC,84,81,83,85
QUADRATIC,122,121,123,125
QUADRATIC,124,121,123,125
QUADRATIC,126,125,127,129
QUADRATIC,128,125,127,129
QUADRATIC,130,129,131,133
QUADRATIC,132,129,131,133
QUADRATIC,149,129,169,209
QUADRATIC,189,129,169,209
*BOUNDARY
 NY,1
 300,1,6
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=A1
1.0, 
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,0.3
**
*STEP
*STATIC,DIRECT
 1., 
*BOUNDARY
 1,2,2,-.4
*CONTACT PRINT ,SLAVE=BOTSUR,MASTER=RIGSUR
*OUTPUT,FIELD, OP=NEW
*CONTACT OUTPUT, VARIABLE=ALL,SLAVE=BOTSUR,MASTER=RIGSUR
*ELEMENT OUTPUT,VARIABLE=PRESELECT
*NODE OUTPUT
U,
*END STEP