*HEADING
Double-sided surface, general contact 
Compression of cylindrical shells in rigid box
**
*PREPRINT, CONTACT=YES, MODEL=YES, HISTORY=YES, ECHO=YES
**
*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.
600.E6,1.
*NSET,NSET=RBREF
10001,10002
*BOUNDARY
10001,1,1
10001,3,6
10002,1,6
**
*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  1,Default
**
*STEP
*DYNAMIC,EXPLICIT
0,0.060
*BOUNDARY,TYPE=VELOCITY
10001,2,2,-25.0
**
*CONTACT
*CONTACT INCLUSIONS, ALL ELEMENT BASED
**
*MONITOR,NOD=10001,DOF=2
**
** Output requests for postprocessing
**
*OUTPUT,FIELD,VAR=PRESELECT,NUM=2
*OUTPUT,HISTORY,TIME INTERVAL=2.E-4
*ENERGY OUTPUT, VAR=PRESELECT
**
** Output requests for qa testing
**
*ELSET,ELSET=QA_TESTEL
75,76,95,96,
*NSET,NSET=QA_TESTN
100,121,142,
*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,
*FILE OUTPUT,NUMBER=4,TIMEMARKS=YES
*ENERGY FILE
*NODE FILE
U,RF
*END STEP