*HEADING
 BARREL VAULT - FINE IRREGULAR MESH S3R    
*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=S3R
 1, 1,203,201
 2, 1,3,203
*ELGEN,ELSET=ROOF
 1,12,2,2,8,200,100
 2,12,2,2,8,200,100
*ELEMENT,TYPE=S3R
 1001, 1609,1811,1809
 1002, 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,ELSET=ROOF,MATERIAL=STEEL
 3.0,3
*MATERIAL,NAME=STEEL
*ELASTIC
 3.E6, 
*DENSITY
8.
*BOUNDARY
 APEX,YSYMM
 SYMM,ZSYMM
 WALL,1,2
*MPC
BEAM,1809,1607
BEAM,2009,1605
BEAM,2209,1603
BEAM,2409,1601
*AMPLITUDE,NAME=RAMP,DEFINITION=SMOOTH STEP
 0.0, 0.0,150, 1.0
*STEP,NLGEOM=NO
 SELF-WEIGHT
*DYNAMIC, EXPLICIT
  ,225
*DLOAD,AMP=RAMP
 ROOF,BX,-0.20833
*OUTPUT, FIELD, VAR=PRESELECT, NUMBER INTERVAL=5
**
** OUTPUT REQUESTS FOR QA PURPOSES
**
*NSET,NSET=QA_TESTN
CORNER,
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TESTN
U,
*FILE OUTPUT,NUM=1
*NODE FILE,NSET=CORNER
U, 
*ENERGY FILE
*END STEP