*HEADING
R313SHL: NAFEMS TEST 3, three-layer sandwich shell under
normal pressure loading [S4R]
*NODE
1, 
9,5.
81,0.,5.
89,5.,5.
*NGEN,NSET=BOT
1,9,1
*NGEN,NSET=TOP
81,89,1
*NSET,NSET=LEFT,GEN
1,81,10
*NSET,NSET=RIGHT,GEN
9,89,10
*NFILL,NSET=NALL
BOT,TOP,8,10
*ELEMENT,TYPE=S4R,ELSET=ELALL
1,1,2,12,11
*ELGEN,ELSET=ELALL
1,8,1,1,8,10,8
*ELSET,ELSET=OUTPUT
64, 
*NSET,NSET=NOUT
89, 
*SHELL SECTION,ELSET=ELALL,COMPOSITE
.028,3,FACE
.75,3,CORE
.028,3,FACE
*MATERIAL,NAME=FACE
*ELASTIC,TYPE=LAMINA
10.E6,4.E6,.3,1.875E6,1.875E6,1.875E6
*DENSITY
1.e-4, 
*MATERIAL,NAME=CORE
*ELASTIC,TYPE=LAMINA
10.,10.,0.,10.,3.E4,1.2E4
*DENSITY
1.e-4, 
*BOUNDARY
BOT,1
BOT,3
BOT,5
TOP,2
TOP,4
LEFT,2,4
RIGHT,1
RIGHT,5
**
*AMPLITUDE, NAME=RAMP
0.,0.,0.1,1.0
*NSET, NSET=CHECK
45, 89 
*NSET,NSET=QA_TEST_CHECK 
CHECK ,
**
*STEP
*DYNAMIC, EXPLICIT
,0.2
*DLOAD, AMP=RAMP
ELALL,P,100
*FILE OUTPUT, NUMBER=1, TIME MARKS=YES
*NODE FILE, NSET=CHECK 
U, 
*ENERGY FILE
***OUTPUT,HISTORY,TIME INTERVAL=0.001
***ENERGY OUTPUT
**ALLIE,ALLWK,ALLKE,ETOTAL,ALLVD,ALLPD,DT
***ELEMENT OUTPUT,ELSET=OUTPUT
**S,
***NODE OUTPUT,NSET=NOUT
**U,
*RESTART, WRITE, NUM=10
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_CHECK 
U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.2
*END STEP