*HEADING
 BOUNCING BALL - DYNAMIC ANALYSIS
*NODE
*NODE,NSET=FIX
*NODE,NSET=NALL
*NODE,NSET=CPEGRN
*NSET,NSET=NY
*ELEMENT,TYPE=CPEG8 ,ELSET=CYLINDER
*NODE,NSET=NALL
*NODE,SYSTEM=C,NSET=NALL
*NGEN,LINE=C,NSET=NALL
*NGEN,NSET=NALL
*NSET,NSET=NY
*ELEMENT,TYPE=CPEG8
*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
BOT,S2
*surface,NAME=RIGSUR
RIGID,SPOS
*CONTACT PAIR,INTERACTION=NO_FRI,SMALL SLIDING
BOTSUR,RIGSUR
*SURFACE INTERACTION,NAME=NO_FRI
1.0,
**
**
**
*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
 300,1,6
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=A1,REFNODE=CPEGRN
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,.3
*DENSITY
 7.5E-4,
*NSET,NSET=TOPNODES ,GENERATE
*INITIAL CONDITIONS,TYPE=VELOCITY
NALL,2,-352.
*BOUNDARY
CPEGRN,3,5
*RESTART,WRITE,FREQUENCY=25
*STEP,INC=100
  DYNAMIC HERTZ CONTACT
*DYNAMIC,HAFTOL=1.E10
 2.E-5,1.E-3,,2.E-5
*ENERGY PRINT,FREQUENCY=2
*CONTACT PRINT ,SLAVE=BOTSUR,MASTER=RIGSUR
*CONTACT FILE,SLAVE=BOTSUR,MASTER=RIGSUR
*NODE PRINT,FREQUENCY=20
U,
RF,
*PRINT,RESIDUAL=NO
*NSET,NSET=OUTFIL1,GENERATE
*NODE FILE,NSET=OUTFIL1,FREQUENCY=10
U,
*NSET,NSET=OUTFIL2,GENERATE
*NODE FILE,NSET=OUTFIL2,FREQUENCY=10
V,A
*NSET,NSET=OUTFIL3,GENERATE
*NODE FILE,NSET=OUTFIL3,FREQUENCY=10
RF,
*CONTROLS,PARAMETERS=FIELD,FIELD=DISPLACEMENT
 , , 0.1
*END STEP