*HEADING
Double-sided surface 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
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
*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
*elset, elset=testel
 QA_TESTEL,
*NSET,NSET=QA_TESTN,ELSET=QA_TESTEL
*NSET,NSET=TESTRB
RBREF,
*nset, nset=test_10001
 10001,
*nset, nset=test_10002
 10002,
*OUTPUT,FIELD,NUM=3
*ELEMENT OUTPUT,ELSET=TESTEL
PEEQ,
*NODE OUTPUT,NSET=QA_TESTN
U,
*NODE OUTPUT,NSET=TESTRB
U,RF
*Output, History, time interval=5e-5
*Node Output, nset=test_10001
 U2,RF2
*Node Output, nset=test_10002
 RF2,
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP