*HEADING
 BARREL VAULT - COARSE IRREGULAR MESH SC6R     
*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.
90001, 301.5,0.,0., 1.,0.,0.
90009, 301.5,0.,150., 1.,0.,.0
90025, 301.5,0.,300., 1.,0.,0.
91601, 301.5,40.,0., 1.,0.,0.
91609, 301.5,20.,150., 1.,0.,0.
91625, 301.5,20.,300., 1.,0.,0.
92409, 301.5,40.,0., 1.,0.,0.
92425, 301.5,40.,300., 1.,0.,0.
*NGEN,NSET=APEX1
 1,9
 90001,90009
*NGEN,NSET=APEX2
 9,25
 90009,90025
*NSET,NSET=APEX
 APEX1,APEX2
*NGEN,NSET=WALL
 1,1601,100
 90001,91601,100
*NGEN,NSET=FREE
 2409, 2425
 92409, 92425
*NGEN,NSET=SYMM1
 25,1625,100
 90025,91625,100
*NGEN,NSET=SYMM2
 1625,2425,100
 91625, 92425,100
*NSET,NSET=SYMM
 SYMM1,SYMM2
*NGEN,NSET=NODES1
 9,1609,100
 90009,91609,100
*NGEN,NSET=NODES2
 1609,1625
 91609,91625
*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=SC6R
 1, 90001,90405,90401, 1,405,401
 2, 90001,90005,90405, 1,5,405
*ELGEN,ELSET=ROOF
 1,6,4,2,4,400,100
 2,6,4,2,4,400,100
*ELEMENT,TYPE=SC6R
 1001, 91609,92013,92009, 1609,2013,2009
 1002, 91609,91613,92013, 1609,1613,2013
*ELGEN,ELSET=ROOF
 1001,4,4,2,2,400,100
 1002,4,4,2,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,92009,91605
 BEAM,92409,91601
*STEP
 SELF-WEIGHT
*STATIC
*DLOAD
 ROOF,BX,-0.20833
*EL PRINT,FREQ=0
*EL FILE,ELSET=ROOF1
S, 
*NODE FILE,NSET=CORNER
U, 
RF, 
*END STEP