*HEADING
PINCHED CYL. WITH DIAPH. AT ENDS - FINE IRREGULAR MESH - S4R,TYPE 1  
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*RESTART,TIMEMARKS=YES,WRITE,NUMBER INTERVAL=1
*NODE,NSET=P1
1,300.,0.,0., 1.,0.,0.
9,300.,63.2616,0., 1.,0.,0.
25,300.,90.,0., 1.,0.,0.
1601,300.,0.,300., 1.,0.,0.
1609,300.,63.2616,140., 1.,0.,0.
1625,300.,90.,140., 1.,0.,0.
2409,300.,0.,300., 1.,0.,0.
2425,300.,90.,300., 1.,0.,0.
*NGEN,NSET=BND1A
1,9
*NGEN,NSET=BND1B
9,25
*NSET,NSET=BND1
BND1A,BND1B
*NGEN,NSET=BND2
1,1601,100
*NSET,NSET=BND22,GEN
1,1501,100
*NGEN,NSET=BND3
2409,2425
*NGEN,NSET=BND4A
25,1625,100
*NGEN,NSET=BND4B
1625,2425,100
*NSET,NSET=BND4
BND4A,BND4B
*NGEN,NSET=NODES1
9,1609,100
*NGEN,NSET=NODES2
1609,1625
*NFILL,NSET=ALL
BND2,NODES1,8
NODES2,BND3,8,100
NODES1,BND4A,16
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
0.,0.,0., 300.,0.,0.
0.,300.,0.

*ELEMENT,TYPE=S4R
1,1,2,102,101
*ELGEN,ELSET=EALL
1,24,1,1,16,100,100
*ELEMENT,TYPE=S4R
 8001,1609,1610,1710,1709
*ELGEN,ELSET=EALL
 8001,16,1,1,8,100,100
*ELSET,ELSET=PRINT
 1, 
*ELSET,ELSET=ELOUT,GENERATE
20,420,100
21,421,100
22,422,100
23,423,100
24,424,100
*NSET,NSET=NOUT
BND1,BND22,BND3,BND4
*NSET,NSET=PRINT
25,
*TRANSFORM,TYPE=C,NSET=ALL
.0,.0,.0,300,.0,.0
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT,SECTION INTEGRATION=GAUSS
3.0, 
*MATERIAL,NAME=MAT
*ELASTIC
3.0E6,0.3
*DENSITY
3.0E-2, 
*BOUNDARY
BND22,YSYMM
BND1,ZSYMM
BND4,YSYMM
BND3,2
BND3,1
*AMPLITUDE,NAME=RAMP
0.0,0.0, 0.8,1.0, 10.0,1.0
*MPC
 TIE,1601,2409
 TIE,1709,1608
 TIE,1809,1607
 TIE,1909,1606
 TIE,2009,1605
 TIE,2109,1604
 TIE,2209,1603
 TIE,2309,1602
*NSET,NSET=QA_TEST_BND1	
BND1	,
**
*STEP,NLGEOM=YES
*DYNAMIC,EXPLICIT
 ,10.0
*DLOAD
EALL,VP,3.4
*CLOAD, AMPLITUDE=RAMP
25,1,-0.25
***OUTPUT,HISTORY,TIME INTERVAL=0.01
***ELEMENT OUTPUT,ELSET=PRINT
**S,LE
***NODE OUTPUT,NSET=PRINT
**U,RF
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLWK,ETOTAL,DT
*FILE OUTPUT,TIMEMARKS=YES, NUMBER INTERVAL=4
*NODE FILE,NSET=BND1	
U,
*ENERGY FILE
*OUTPUT, FIELD, TIMEMARKS=YES, NUMBER INTERVAL=4
*NODE OUTPUT, NSET=QA_TEST_BND1	
U,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=2.5
*END STEP