*HEADING
  BARREL VAULT 4 X 4 MESH S9R5            
*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=S9R5,ELSET=ROOF1
 1,1,21,23,3,11,22,13,2,12
*ELGEN,ELSET=ROOF
 1,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