*HEADING
**
** This test is the same as multpenaltycont1.inp except that one of
** the contact pairs (deformable plate to rigid plate contact) is
** kinematically enforced. This test demonstrates that a rigid
** body (the rigid plate in this case) can be involved in both
** kinematic contact and penalty contact. Results are similar to
** multpenaltycont1 except that the deformable plate does not bounce
** off of the rigid plate.
**
**  Rigid sphere with radius of .01 m.
*NODE,NSET=BALL,INPUT=sphere_n.inp
*NODE,NSET=RIG_SPH
599991,0.,0.,0.01
*ELEMENT,TYPE=R3D4,ELSET=RIG_SPH,INPUT=sphere_e.inp
*ELEMENT,TYPE=MASS,ELSET=SPH_MASS
50000,599991
*MASS,ELSET=SPH_MASS
1.e-2,
**
**  Rigid plate (see analytical rigid surface defn)
*NODE,NSET=RIG_PL
1001,0.,0.,0.025
**
**  Define mesh for deformable plate
*NODE
90001,-.1,-.1,-0.001
90036, .1,-.1,-0.001
91261,-.1, .1,-0.001
91296, .1, .1,-0.001
*NGEN,NSET=B1
90001,90036,1
*NGEN,NSET=B2
91261,91296,1
*NFILL,NSET=DEF_PL
B1,B2,35,36
*ELEMENT,TYPE=S4R,ELSET=DEF_PL
3001, 90001,90002,90038,90037
*ELGEN,ELSET=DEF_PL
3001, 35,1,1, 35,36,35
*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
  0.17129E+09, 0.64921E-02
  0.18129E+09, 0.85618E-02
  0.19129E+09, 0.11080E-01
  0.20129E+09, 0.14110E-01
  0.21129E+09, 0.17723E-01
  0.22129E+09, 0.21991E-01
  0.23129E+09, 0.26994E-01
  0.24129E+09, 0.32819E-01
  0.25129E+09, 0.39556E-01
  0.26129E+09, 0.47301E-01
  0.27129E+09, 0.56159E-01
  0.28129E+09, 0.66236E-01
  0.29129E+09, 0.77648E-01
  0.30129E+09, 0.90516E-01
  0.31129E+09, 0.10497E+00
  0.32129E+09, 0.12114E+00
  0.33129E+09, 0.13916E+00
  0.34129E+09, 0.15919E+00
  0.35129E+09, 0.18138E+00
  0.36129E+09, 0.20588E+00
  0.37129E+09, 0.23287E+00
  0.38129E+09, 0.26252E+00
  0.39129E+09, 0.29502E+00
  0.40129E+09, 0.33054E+00
  0.41129E+09, 0.36929E+00
  0.42129E+09, 0.41147E+00
  0.43129E+09, 0.45729E+00
  0.44129E+09, 0.50696E+00
  0.45129E+09, 0.56073E+00
  0.46129E+09, 0.61881E+00
  0.47129E+09, 0.68145E+00
  0.48129E+09, 0.74890E+00
  0.49129E+09, 0.82142E+00
  0.50129E+09, 0.89928E+00
  0.51129E+09, 0.98274E+00
  0.52129E+09, 0.10721E+01
**
*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=CYLINDER, NAME=RIG_PL
0., 0., 0., 1., 0., 0.
0., 1., 0.
start, 0.12,  0.022
line, -0.12,  0.022
*RIGID BODY,REF NODE=1001, ANALYTICAL SURFACE =RIG_PL
*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 PAIR,MECHANICAL CONSTRAINT=PENALTY
RIG_SPH,DEF_PL
RIG_SPH,RIG_PL
*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
3631,3665,1
*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