*HEADING
 BARREL VAULT - COARSE IRREGULAR MESH S4R     
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE,NSET=CORNER
1, 300.,0.,0., 1.,0.,0.
9, 300.,0.,150., 1.,0.,.0
25, 300.,0.,300., 1.,0.,0.
1601, 300.,40.,0., 1.,0.,0.
1609, 300.,20.,150., 1.,0.,0.
1625, 300.,20.,300., 1.,0.,0.
2409, 300.,40.,0., 1.,0.,0.
2425, 300.,40.,300., 1.,0.,0.
*NGEN,NSET=APEX1
1,9
*NGEN,NSET=APEX2
9,25
*NSET,NSET=APEX
APEX1,APEX2
*NGEN,NSET=WALL
1,1601,100
*NGEN,NSET=FREE
2409,2425
*NGEN,NSET=SYMM1
25,1625,100
*NGEN,NSET=SYMM2
1625,2425,100
*NSET,NSET=SYMM
SYMM1,SYMM2
*NGEN,NSET=NODES1
9,1609,100
*NGEN,NSET=NODES2
1609,1625
*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=S4R
 1, 1,5,405,401
*ELGEN,ELSET=ROOF
 1, 6,4,1,4,400,100
*ELEMENT,TYPE=S4R
 1001, 1609,1613,2013,2009
*ELGEN,ELSET=ROOF
 1001, 4,4,1,2,400,100
*ELSET,ELSET=ROOF1
 1, 
*SHELL SECTION,MATERIAL=STEEL,ELSET=ROOF
3.0,3
*MATERIAL,NAME=STEEL
*ELASTIC
3.E6, 
*BOUNDARY
 APEX,YSYMM
 SYMM,ZSYMM
 WALL,1,2
*MPC
 BEAM,2009,1605
 BEAM,2409,1601
*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