*HEADING
 NAFEMS R0031/3 THREE-LAYER SANDWICH SHELL UNDER NORMAL PRESSURE LOADING [SC6R]
 STACKED
*NODE
   1, 0., 0., -0.403
   9, 5., 0., -0.403
  81, 0., 5., -0.403
  89, 5., 5., -0.403
1001, 0., 0., -0.375
1009, 5., 0., -0.375
1081, 0., 5., -0.375
1089, 5., 5., -0.375
2001, 0., 0., 0.375
2009, 5., 0., 0.375
2081, 0., 5., 0.375
2089, 5., 5., 0.375
3001, 0., 0., 0.403
3009, 5., 0., 0.403
3081, 0., 5., 0.403
3089, 5., 5., 0.403
*NGEN,NSET=BOT
   1,   9,1
1001,1009,1
2001,2009,1
3001,3009,1
*NGEN,NSET=TOP
  81,  89,1
1081,1089,1
2081,2089,1
3081,3089,1
*NSET,NSET=LEFT,GEN
   1,  81,10
1001,1081,10
2001,2081,10
3001,3081,10
*NSET,NSET=RIGHT,GEN
   9,  89,10
1009,1089,10
2009,2089,10
3009,3089,10
*NFILL,NSET=NALL
BOT,TOP,8,10
*ELEMENT,TYPE=SC6R,ELSET=ELALL
   1,1,2,11,1001,1002,1011
   2,2,12,11,1002,1012,1011
1001,1001,1002,1011,2001,2002,2011
1002,1002,1012,1011,2002,2012,2011
2001,2001,2002,2011,3001,3002,3011
2002,2002,2012,2011,3002,3012,3011
*ELGEN,ELSET=ELFACES
1,8,1,2,8,10,16
2,8,1,2,8,10,16
2001,8,1,2,8,10,16
2002,8,1,2,8,10,16
*ELGEN,ELSET=ELCORE
1001,8,1,2,8,10,16
1002,8,1,2,8,10,16
*ELSET,ELSET=ELALL
ELCORE,
ELFACES,
*ELSET,ELSET=ELBOTFACE,GEN
1,178,1
*ELSET,ELSET=OUTPUT
 127,
 128,
1127,
1128,
2127,
2128,
*NSET,NSET=NOUT
89,
1089,
2089,
3089, 
*SHELL SECTION,ELSET=ELFACES, MATERIAL=FACE
.028,3
*SHELL SECTION,ELSET=ELCORE, MATERIAL=CORE
.75,3
*MATERIAL,NAME=FACE
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 10.E6,4.E6,4.E6,.3,.3,.3,1.875E6,1.875E6,
 1.875E6,
*MATERIAL,NAME=CORE
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 10.,10.,10.,0.,0.,0.,10.,3.E4,
 1.2E4,
*BOUNDARY
BOT,1
BOT,3
TOP,2
LEFT,2,3
RIGHT,1
**
*STEP
*STATIC
*DLOAD
ELBOTFACE,P1,100
*OUTPUT, FIELD
*ELEMENT OUTPUT, POSITION=NODES
1,2,3
 S,
*NODE OUTPUT
 U,
*EL FILE,ELSET=OUTPUT,POSITION=AVERAGED AT NODES
1,2,3
S, 
*END STEP