*HEADING
 CONTACT PROBLEM - DISPLACEMENT CONTROL, COARSE MESH
 TOP SURFACE CONSTRAINED WITH KINEMATIC COUPLING CONSTRAINT
*RESTART,WRITE,FREQUENCY=999
*NODE
 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, 
** *KINEMATIC COUPLING, REF NODE = 1
** 3,2
** 5,2
** 7,2
** 9,2
** 11,2
** 13,2
*NSET, NSET=NREF
 1,
*ELSET, ELSET=ELEMTOP
 1,2,3
*SURFACE, NAME=SURFTOP
 ELEMTOP,S4
*COUPLING, CONSTRAINT NAME=COUPLING_1, REF NODE=NREF, SURFACE=SURFTOP
*KINEMATIC
 2,
*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
 1,6,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
*CONTACT FILE ,SLAVE=BOTSUR,MASTER=RIGSUR
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=BOTSUR,MASTER=RIGSUR
*END STEP