*HEADING
 BARREL VAULT - FINE IRREGULAR MESH S4R
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE,NSET=CORNER
*NGEN,NSET=APEX1
*NGEN,NSET=APEX2
*NSET,NSET=APEX
 APEX1,APEX2
*NGEN,NSET=WALL
*NGEN,NSET=FREE
*NGEN,NSET=SYMM1
*NGEN,NSET=SYMM2
*NSET,NSET=SYMM
  SYMM1,SYMM2
*NGEN,NSET=NODES1
*NGEN,NSET=NODES2
*NFILL,NSET=ALL
 NODES1,SYMM1,16
 NODES2,FREE,8,100
 WALL,NODES1,8
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
*NSET, NSET=N1, GEN
*NSET, NSET=N2, GEN
*EQUATION
 N1, 1, 1.0, N2, 1, 1.0
 N1, 2, 1.0, N2, 2, 1.0
*ELEMENT,TYPE=SC8R
*ELGEN,ELSET=ROOF
*ELEMENT,TYPE=SC8R
*ELGEN,ELSET=ROOF
*ELSET,ELSET=ROOF1
*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