*HEADING
  BARREL VAULT 4 X 4 MESH STRI65                  
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE,SYSTEM=C,NSET=CORNER
 1,300.,0.,0.,
 9,300.,,300.
 81,300.,40.,
 89,300.,40.,300.
*NGEN,LINE=C,NSET=WALL
 1,81,10
*NGEN,LINE=C,NSET=SYMM
 9,89,10,,,,300.
*NFILL
WALL,SYMM,8,1
*NSET,NSET=APEX,GENERATE
1,9,1
*ELEMENT,TYPE=STRI65,ELSET=ROOF1
    1,     1,    21,    23,    11,    22,    12
50001,     1,    23,     3,    12,    13,     2
*ELGEN,ELSET=ROOF
    1,     4,    20,     4,     4,     2,     1
50001,     4,    20,     4,     4,     2,     1
*ELSET,ELSET=FEW
 1,2,3,4
*MATERIAL,NAME=STEEL
*ELASTIC
 3.E6, 
*BOUNDARY
 APEX,YSYMM
 SYMM,ZSYMM
 WALL,1,2
*SHELL SECTION ,MATERIAL=STEEL,ELSET=ROOF
 3.0,3
*STEP
 SELF-WEIGHT
*STATIC
*DLOAD
 ROOF,BX,-0.20833
*ENERGY PRINT
*EL PRINT,ELSET=FEW
 E, S
*EL PRINT
SINV, 
ENER, 
*EL FILE,POSITION=AVERAGED AT NODES
S,E
*NODE FILE,NSET=CORNER
 U, CF
*ENERGY FILE
*END STEP