*HEADING
 BOUNCING BALL - DYNAMIC ANALYSIS     
*NODE
 300,0.     , -11.
*NODE,NSET=FIX
 301,0.     , -10.00706
 302,0.26177, -10.00706
 303,0.52336, -10.00706
 304,0.78459, -10.00706
 305,1.04528, -10.00706
 306,1.30526, -10.00706
 307,1.56434, -10.00706
 308,1.82236, -10.00706
 309,2.07912, -10.00706
*NODE,NSET=NALL
 401,0.,5.
 405,2.,4.89898
 409,5.5,4.175825
 501,0.,10.
 503 , 1.00509 , 9.949361
 505,2.,9.79796
 507,4.6619,8.84684
 509,5.5,8.35165
 513,8.8034,4.7434
*NODE,NSET=CPEGRN
 9999, 0., 0.
*NSET,NSET=NY
 401,501
*ELEMENT,TYPE=CPEG8 ,ELSET=CYLINDER
 21,1,5,505,501,3,405,503,401
 22,5,9,509,505,7,409,507,405
 23,9,9,13,509,9,11,513,409
*NODE,NSET=NALL
   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,NSET=NALL
 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,NSET=NALL
 209,213
 201,209
*NGEN,NSET=NALL
 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=CPEG8
 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
BOT,S2
*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,REFNODE=CPEGRN
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,.3
*DENSITY
 7.5E-4, 
*NSET,NSET=TOPNODES ,GENERATE
 1,13
*INITIAL CONDITIONS,TYPE=VELOCITY
NALL,2,-352.
*BOUNDARY
CPEGRN,3,5
*STEP,INC=100
  DYNAMIC HERTZ CONTACT
*DYNAMIC,HAFTOL=1.E10,APPLICATION=TRANSIENT FIDELITY
 2.E-5,1.E-3,,2.E-5
*ENERGY PRINT,FREQUENCY=2
*CONTACT PRINT ,SECONDARY=BOTSUR,MAIN=RIGSUR
*PRINT,contact=yes
*NSET,NSET=OUTFIL1,GENERATE
 201,209
*contact controls, stiff=10.0
*END STEP