*HEADING
Double-sided surface, penalty contact 
Compression of cylindrical shells in rigid box
**
*NODE,INPUT=shell_compact_ef1.inp
**
**
*ELEMENT, TYPE=S4R, ELSET=RINGS,INPUT=shell_compact_ef2.inp
*ELEMENT,TYPE=R3D4,ELSET=RIGID
     601,     694,     695,     697,     696
     602,     698,     699,     701,     700
     603,     698,     702,     705,     700
     604,     699,     703,     704,     701
     605,     698,     699,     703,     702
     606,     700,     701,     704,     705
**
** lid
**
*NSET, NSET=LID
     694,     695,     696,      697
**
** lid
**
*ELSET, ELSET=LID
     601,
*NSET, NSET=BOX,GEN
698,705
*ELSET, ELSET=BOX,GEN
602,606
*NSET, NSET=RING1,GEN
1,231
*ELSET, ELSET=RING1,GEN
1,200
*NSET, NSET=RING2,GEN
232,462
*ELSET, ELSET=RING2,GEN
201,400
*NSET, NSET=RING3,GEN
463,693
*ELSET, ELSET=RING3,GEN
401,600
*NSET, NSET=RINGS
RING1,RING2,RING3
*SHELL SECTION,ELSET=RING1,MATERIAL=DEFORM
0.075, 
*SHELL SECTION,ELSET=RING2,MATERIAL=DEFORM
0.05, 
*SHELL SECTION,ELSET=RING3,MATERIAL=DEFORM
0.075, 
*MATERIAL,NAME=DEFORM
*DENSITY
7850., 
*ELASTIC
200.E9,.3
*PLASTIC
250.E6,0.
700.E7,1.
*NSET,NSET=RBREF
10001,10002
*BOUNDARY
10001,1,1
10001,3,6
10002,1,6
**
** step  1,Default
**
*SURFACE,TYPE=ELEMENT,NAME=RING1
RING1, 
*SURFACE,TYPE=ELEMENT,NAME=RING2
RING2, 
*SURFACE,TYPE=ELEMENT,NAME=RING3
RING3, 
*SURFACE,TYPE=ELEMENT,NAME=LID
LID, 
*SURFACE,TYPE=ELEMENT,NAME=BOX
BOX, 
*RIGID BODY,ELS=LID,REFNOD=10001,POSITION=CENTEROFMASS,DENSITY=1e4
0.01, 
*RIGID BODY,ELS=BOX,REFNOD=10002,POSITION=INPUT,DENSITY=1e4
0.02, 
*STEP
*DYNAMIC,EXPLICIT
0,0.010
*BOUNDARY,TYPE=VELOCITY
10001,2,2,-130.0
*CONTACT PAIR,MECHANICAL CONSTRAINT=PENALTY
RING1,RING3
RING2,RING3
RING1,BOX
RING2,BOX
RING3,BOX
RING1,LID
RING1,RING2
*MONITOR,NOD=10001,DOF=2
**
** Output requests for postprocessing
**
*OUTPUT,FIELD,VAR=PRESELECT,NUM=2
*OUTPUT,HISTORY,TIME INTERVAL=2.E-4
*ENERGY OUTPUT
ALLIE,ALLKE,ALLPD,ALLWK
**
** Output requests for qa testing
**
*ELSET,ELSET=QA_TESTEL,GEN
301,320,1
*NSET,NSET=QA_TESTN,ELSET=QA_TESTEL
*NSET,NSET=QA_TESTRB
RBREF,
*OUTPUT,FIELD,NUM=3
*ELEMENT OUTPUT,ELSET=QA_TESTEL
PEEQ,
*NODE OUTPUT,NSET=QA_TESTN
U,
*NODE OUTPUT,NSET=QA_TESTRB
U,RF
*OUTPUT,HISTORY,TIME INTERVAL=2.5E-3
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP