*HEADING
 BARREL VAULT  ( 2 X 2 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,20
*NGEN,LINE=C,NSET=SYMM
 9,89,20,,,,300.
*NFILL
WALL,SYMM,4,2
*NGEN,NSET=APEX
 1,9,2
*ELEMENT,TYPE=STRI65,ELSET=ROOF1
    1,     1,     5,    45,     3,    25,    23
50001,     1,    45,    41,    23,    43,    21
*ELGEN,ELSET=ROOF
    1,     2,     4,     1,     2,    40,     2
50001,     2,     4,     1,     2,    40,     2
*ELSET,ELSET=FEW
 1,2
*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
*EL PRINT,ELSET=FEW
*NODE FILE,NSET=CORNER
U, 
RF, 
*EL FILE,ELSET=ROOF1
S,SINV
SF, 
*END STEP