*HEADING
 BARREL VAULT - FINE IRREGULAR MESH STRI65  
*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=STRI65
    1,     1,     5,   405,     3,   205,   203
50001,     1,   405,   401,   203,   403,   201
*ELGEN,ELSET=ROOF
    1,     6,     4,     1,     4,   400,   100
50001,     6,     4,     1,     4,   400,   100
*ELEMENT,TYPE=STRI65
 1001,  1609,  1613,  2013,  1611,  1813,  1811
51001,  1609,  2013,  2009,  1811,  2011,  1809
*ELGEN,ELSET=ROOF
 1001,     4,     4,     1,     2,   400,   100
51001,     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,1809,1607
 BEAM,2009,1605
 BEAM,2209,1603
 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