*HEADING
 BARREL VAULT - COARSE IRREGULAR MESH SC8R     
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE,NSET=CORNER
      1, 298.5, 0.,  0., 1.,0.,0.
      9, 298.5, 0.,150., 1.,0.,0.
     25, 298.5, 0.,300., 1.,0.,0.
   1601, 298.5,40.,  0., 1.,0.,0.
   1609, 298.5,20.,150., 1.,0.,0.
   1625, 298.5,20.,300., 1.,0.,0.
   2409, 298.5,40.,  0., 1.,0.,0.
   2425, 298.5,40.,300., 1.,0.,0.
 100001, 301.5, 0.,  0., 1.,0.,0.
 100009, 301.5, 0.,150., 1.,0.,0.
 100025, 301.5, 0.,300., 1.,0.,0.
 101601, 301.5,40.,  0., 1.,0.,0.
 101609, 301.5,20.,150., 1.,0.,0.
 101625, 301.5,20.,300., 1.,0.,0.
 102409, 301.5,40.,  0., 1.,0.,0.
 102425, 301.5,40.,300., 1.,0.,0.
*NGEN,NSET=APEX1
 1,9
 100001,100009
*NGEN,NSET=APEX2
 9,25
 100009,100025
*NSET,NSET=APEX
 APEX1,APEX2
*NGEN,NSET=WALL
 1,1601,100
 100001,101601,100
*NGEN,NSET=FREE
 2409,2425
102409,102425
*NGEN,NSET=SYMM1
 25,1625,100
 100025,101625,100
*NGEN,NSET=SYMM2
 1625,2425,100
 101625,102425,100
*NSET,NSET=SYMM
 SYMM1,SYMM2
*NGEN,NSET=NODES1
 9,1609,100
 100009,101609,100
*NGEN,NSET=NODES2
 1609,1625
 101609,101625
*NFILL,NSET=ALL
 NODES1,SYMM1,16
 NODES2,FREE,8,100
 WALL,NODES1,8
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
 0.,0.,0., 0.,0.,300.
 300.,0.,0.

*ELEMENT,TYPE=SC8R
 1,100001,100005,100405,100401, 1, 5, 405, 401
*ELGEN,ELSET=ROOF
 1, 6,4,1,4,400,100
*ELEMENT,TYPE=SC8R
 1001, 101609,101613,102013,102009, 1609, 1613, 2013, 2009
*ELGEN,ELSET=ROOF
 1001, 4,4,1,2,400,100
*ELSET,ELSET=ROOF1
 1, 
*SHELL SECTION,MATERIAL=STEEL,ELSET=ROOF,POISSON=ELASTIC
 3.0,3
*MATERIAL,NAME=STEEL
*ELASTIC
 3.E6, 
*BOUNDARY
 APEX,YSYMM
 SYMM,ZSYMM
 WALL,1,2
*MPC
 BEAM,2009,1605
 BEAM,2409,1601
 BEAM,102009,101605
 BEAM,102409,101601
*STEP
 SELF-WEIGHT
*STATIC
*DLOAD
 ROOF,BX,-0.20833
*EL PRINT,FREQUENCY=0
*EL FILE,ELSET=ROOF1
 S, 
*NODE FILE,NSET=CORNER
 U, 
 RF, 
*END STEP