*HEADING
 BARREL VAULT - FINE IRREGULAR MESH S4R  
*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.

*NSET, NSET=N1, GEN
 201, 1601, 200
*NSET, NSET=N2, GEN
 100201, 101601, 200
*EQUATION
 2,
 N1, 1, 1.0, N2, 1, 1.0
 2,
 N1, 2, 1.0, N2, 2, 1.0
*ELEMENT,TYPE=SC8R
 1, 100001,100003,100203,100201, 1,3,203,201
*ELGEN,ELSET=ROOF
 1, 12,2,1,8,200,100
*ELEMENT,TYPE=SC8R
 1001, 101609,101611,101811,101809, 1609,1611,1811,1809
*ELGEN,ELSET=ROOF
 1001, 8,2,1,4,200,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,1809,1607
 BEAM,2009,1605
 BEAM,2209,1603
 BEAM,2409,1601
 BEAM,101809,101607
 BEAM,102009,101605
 BEAM,102209,101603
 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