*HEADING
 CONTACT PROBLEM - DISPLACEMENT CONTROL, COARSE MESH
 TOP SURFACE CONSTRAINED WITH KINEMATIC COUPLING
*RESTART,WRITE,FREQUENCY=999
*NODE
*NODE,NSET=FIX
*NODE
*NODE,SYSTEM=C
*NGEN,LINE=C
*NGEN
*NSET,NSET=NY
*ELEMENT,TYPE=CPE8
*ELGEN
*ELSET,ELSET=TOP
*ELSET,ELSET=BOT
*ELSET,ELSET=CYLINDER
TOP,BOT
*ELEMENT,TYPE=R2D2
*ELGEN,ELSET=RIGID
*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
*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