*HEADING
 BARREL VAULT - FINE 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.
10001, 301.5,0.,0., 1.,0.,0.
10009, 301.5,0.,150., 1.,0.,.0
10025, 301.5,0.,300., 1.,0.,0.
11601, 301.5,40.,0., 1.,0.,0.
11609, 301.5,20.,150., 1.,0.,0.
11625, 301.5,20.,300., 1.,0.,0.
12409, 301.5,40.,0., 1.,0.,0.
12425, 301.5,40.,300., 1.,0.,0.
*NGEN,NSET=APEX1
 1,9
 10001,10009
*NGEN,NSET=APEX2
 9,25
 10009,10025
*NSET,NSET=APEX
 APEX1,APEX2
*NGEN,NSET=WALL
 1,1601,100
 10001,11601,100
*NGEN,NSET=FREE
 2409,2425
 12409, 12425
*NGEN,NSET=SYMM1
 25,1625,100
 10025,11625,100
*NGEN,NSET=SYMM2
 1625,2425,100
 11625, 12425,100
*NSET,NSET=SYMM
 SYMM1,SYMM2
*NGEN,NSET=NODES1
 9,1609,100
 10009,11609,100
*NGEN,NSET=NODES2
 1609,1625
 11609,11625
*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, 10001,10203,10201, 1,203,201
 2, 10001,10003,10203, 1,3,203
*ELGEN,ELSET=ROOF
 1,12,2,2,8,200,100
 2,12,2,2,8,200,100
*ELEMENT,TYPE=SC6R
 1001, 11609,11811,11809, 1609,1811,1809
 1002, 11609,11611,11811, 1609,1611,1811
*ELGEN,ELSET=ROOF
 1001,8,2,2,4,200,100
 1002,8,2,2,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,11809,11607
 BEAM,12009,11605
 BEAM,12209,11603
 BEAM,12409,11601
*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