*HEADING
 BARREL VAULT  (4 X 4 SC6R)               
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE,SYSTEM=C,NSET=CORNER
  110, 298.5,  0.0,  0.0
  150, 298.5,  0.0, 300.
  510, 298.5, 40.0,  0.0
  550, 298.5, 40.0, 300.
 1110, 301.5,  0.0,  0.0
 1150, 301.5,  0.0, 300.
 1510, 301.5, 40.0,  0.0
 1550, 301.5, 40.0, 300.
 2000,,,
*NGEN,LINE=C,NSET=WALL
 110,510,100,2000
 1110,1510,100,2000
*NGEN,LINE=C,NSET=SYMM
  150,550,100,,0.,0.,298.5
 1150,1550,100,,0.,0.,301.5
*NFILL
 WALL,SYMM,4,10
*NGEN,NSET=APEX
 110,150,10
 1110,1150,10
*ELEMENT,TYPE=SC6R,ELSET=ROOF1
 1,1110,1120,1220,110,120,220
 2,1110,1220,1210,110,220,210
*ELGEN,ELSET=ROOF
 1,4,100,2,4,10,20
 2,4,100,2,4,10,20
*MATERIAL,NAME=STEEL
*ELASTIC
 3.E6, 
*SHELL SECTION ,MATERIAL=STEEL,ELSET=ROOF,POISSON=ELASTIC
 3.,3
*BOUNDARY
 APEX,YSYMM
 SYMM,ZSYMM
 WALL,1,2
*STEP
 SELF WEIGHT
*STATIC
*DLOAD
 ROOF,BX,-0.20833
*ENERGY PRINT
*EL PRINT,FREQ=0
*EL FILE,ELSET=ROOF1
 S,SINV
 SF, 
 ENER, 
 E, 
*NODE FILE,NSET=CORNER
 U, 
 RF, 
*ENERGY FILE
*END STEP