*HEADING
Double-sided surface contact
Compression of cylindrical shells in rigid box, normals
reversed.
**
*NODE,INPUT=shell_compact_ef1.inp
**
**
*ELEMENT, TYPE=S4R, ELSET=RINGS,INPUT=shell_compact_ef3.inp
*ELEMENT,TYPE=R3D4,ELSET=RIGID
**
** lid
**
*NSET, NSET=LID
**
** lid
**
*ELSET, ELSET=LID
*NSET, NSET=BOX,GEN
*ELSET, ELSET=BOX,GEN
*NSET, NSET=RING1,GEN
*ELSET, ELSET=RING1,GEN
*NSET, NSET=RING2,GEN
*ELSET, ELSET=RING2,GEN
*NSET, NSET=RING3,GEN
*ELSET, ELSET=RING3,GEN
*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
*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
*RIGID BODY,ELS=BOX,REFNOD=10002,POSITION=INPUT,DENSITY=1e4
*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
*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