*HEADING
** General contact and kinematic contact
**  Rigid sphere with radius of .01 m.
*NODE,NSET=BALL,INPUT=sphere_n.inp
*NODE,NSET=RIG_SPH
*ELEMENT,TYPE=R3D4,ELSET=RIG_SPH,INPUT=sphere_e.inp
*ELEMENT,TYPE=MASS,ELSET=SPH_MASS
*MASS,ELSET=SPH_MASS
1.e-2,
**
**  Rigid plate (see analytical rigid surface defn)
*NODE,NSET=RIG_PL
*ELEMENT, TYPE=R3D4, ELSET=RIG_PL
**
**  Define mesh for deformable plate
*NODE
*NGEN,NSET=B1
*NGEN,NSET=B2
*NFILL,NSET=DEF_PL
B1,B2,35,36
*ELEMENT,TYPE=S4R,ELSET=DEF_PL
*ELGEN,ELSET=DEF_PL
*SHELL SECTION,MATERIAL=STEEL,ELSET=DEF_PL,SECTION INTEGRA=GAUSS
.00082,5
**
**  Define material properties for deformable plate
**
*MATERIAL,NAME=STEEL
*DENSITY
7800.,
*ELASTIC
2.1E11,0.3
*PLASTIC
  0.91294E+08, 0.00000E+00
  0.10129E+09, 0.21052E-03
  0.11129E+09, 0.52686E-03
  0.12129E+09, 0.97685E-03
  0.13129E+09, 0.15923E-02
  0.14129E+09, 0.24090E-02
  0.15129E+09, 0.34674E-02
  0.16129E+09, 0.48120E-02
**
*INITIAL CONDITIONS,TYPE=VELOCITY
DEF_PL,3, 100.
**
*SURFACE,TYPE=ELEMENT,NAME=DEF_PL
DEF_PL,
*SURFACE,TYPE=ELEMENT,NAME=RIG_SPH
RIG_SPH,SPOS
*SURFACE, TYPE=ELEMENT, NAME=RIG_PL
RIG_PL, SNEG
*RIGID BODY,ELSET=RIG_PL, REF NODE=1001
*RIGID BODY,ELSET=RIG_SPH,REF NODE=599991
*BOUNDARY
1001, 1,6
599991, 4,6
*STEP
*DYNAMIC,EXPLICIT,element by element
,3.E-4
*CONTACT
*CONTACT INCLUSIONS, ALL ELEMENT BASED
*CONTACT PAIR
DEF_PL,RIG_PL
**
** Output requests for postprocessing
**
*NSET,NSET=ODB_RBH
RIG_PL,RIG_SPH
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,HISTORY,TIME INTERVAL=0.5E-6
*INCREMENTATION OUTPUT
DT,
*NODE OUTPUT,NSET=ODB_RBH
U,V,RF
**
** Output requests for qa testing
**
*ELSET,ELSET=QA_TESTEL,GEN
*NSET,NSET=QA_TESTN,ELSET=QA_TESTEL
*NSET,NSET=QA_TESTRB
RIG_PL,RIG_SPH
*****************************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TESTN
U,
*NODE OUTPUT,NSET=QA_TESTRB
U
*OUTPUT,HISTORY,TIME INTERVAL=1.0E-4
*ENERGY OUTPUT,VAR=PRESELECT
*INCREMENTATION OUTPUT
DT,
******************
*END STEP