*HEADING
 BARREL VAULT - FINE IRREGULAR MESH S4 
*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.

*NSET,NSET=MID
2425
*NSET,NSET=QA_TEST_MID
2425
*ELEMENT,TYPE=S4
 1, 1,3,203,201
*ELGEN,ELSET=ROOF
 1, 12,2,1,8,200,100
*ELEMENT,TYPE=S4
 1001, 1609,1611,1811,1809
*ELGEN,ELSET=ROOF
 1001, 8,2,1,4,200,100
*ELSET,ELSET=ROOF1
 1, 
*SHELL SECTION,MATERIAL=STEEL,ELSET=ROOF
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,HISTORY,FREQ=5
*NODE OUTPUT,NSET=MID
 U,
*ENERGY OUTPUT
ALLAE,ALLCD,ALLKE,ALLIE,ALLWK,ETOTAL,ALLSE
*OUTPUT, FIELD, NUMBER INTERVAL=5
*NODE OUTPUT,NSET=QA_TEST_MID
 U,
*FILE OUTPUT,NUMBER INTERVAL=5
*NODE FILE,NSET=QA_TEST_MID
U
*END STEP