*HEADING
 CONTACT PROBLEM  -  RIGID SURFACE  - STABILIZE
*RESTART,WRITE,FREQUENCY=9999
*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
*rigid body,analytical surface=BSURF,REF NODE=999
*surface,TYPE=SEGMENTS,NAME=BSURF
START,0.,-11.
LINE,4.,-11.
*surface,NAME=ASURF
BOT,S2
*CONTACT PAIR,INTERACTION=SMOOTH,SMALL SLIDING
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*SURFACE BEHAVIOR, AUGMENTED LAGRANGE
*EQUATION
*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
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=A1
1.,
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,0.3
*STEP,INC=20
*STATIC
 0.5,1.
*CLOAD
1,2,-4.E6
*BOUNDARY
999,1,6
*CONTACT PRINT,SLAVE=ASURF
*CONTACT CONTROLS,STABILIZE,TANGENT FRACTION=0.,MASTER=BSURF,SLAVE=ASURF
*OUTPUT,HISTORY
*ENERGY OUTPUT
*END STEP