*HEADING
 BARREL VAULT  ( 4 X 4 MESH , SC8R )          
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE,SYSTEM=C,NSET=CORNER
  1,298.5,0.,0.,
  9,298.5,,300.
 81,298.5,40.,
 89,298.5,40.,300.
 1001,301.5,0.,0.,
 1009,301.5,,300.
 1081,301.5,40.,
 1089,301.5,40.,300.
*NGEN,LINE=C,NSET=WALL
 1,81,20
 1001,1081,20
*NGEN,LINE=C,NSET=SYMM
 9,89,20,,,,300.
 1009,1089,20,,,,300.
*NFILL
 WALL,SYMM,4,2
*NGEN,NSET=APEX
 1,9,2
 1001,1009,2
*ELEMENT,TYPE=SC8R, ELSET=ROOF1
 1, 1001, 1003, 1023, 1021, 1,3,23,21
*ELGEN,ELSET=ROOF
 1,4,2,1,4,20,4
*ELSET,ELSET=FEW
 1,2,3,4
*MATERIAL,NAME=STEEL
*ELASTIC
 3.E6, 
*BOUNDARY
 APEX,YSYMM
 SYMM,ZSYMM
 WALL,1,2
*SHELL SECTION,ELSET=ROOF ,MATERIAL=STEEL,POISSON=ELASTIC
 3.0,3
*STEP
 SELF-WEIGHT
*STATIC
*DLOAD
 ROOF,BX,-0.20833
*EL FILE,ELSET=ROOF1
S, 
*NODE FILE,NSET=CORNER
U, 
RF, 
*EL PRINT,ELSET=ROOF
*END STEP