*HEADING
Double-sided surface contact
Compression of cylindrical shells in rigid box, NOTHICK
**
*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,NOTHICK
RING2, 
*SURFACE,TYPE=ELEMENT,NAME=RING3
RING3, 
*SURFACE,TYPE=ELEMENT,NAME=LID,NOTHICK
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
*CONTACT PAIR,CPSET=QA_TESTCP
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=QA_TESTRB
RBREF,
*nset, nset=testrb
rbref,
*OUTPUT,FIELD,NUM=3
*ELEMENT OUTPUT,ELSET=TESTEL
PEEQ,
*NODE OUTPUT,NSET=QA_TESTN
U,
*NODE OUTPUT,NSET=QA_TESTRB
U,
*output, history, freq=1
*node output, nset=testrb
RF,
*OUTPUT,HISTORY,TIME INTERVAL=2.5E-3
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP