*HEADING
PINCHED CYLINDER WITH DIAPHRAGMS AT ENDS - S4R - 10 X 10 MESH
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*RESTART,TIMEMARKS=YES,WRITE,NUMBER INTERVAL=1
*NODE
1, 300.,0.,0., 1.,0.,0.
21, 300.,90.,0., 1.,0.,0.
2001, 300.,0.,300., 1.,0.,0.
2021, 300.,90.,300., 1.,0.,0.
*NGEN,NSET=LOADB
1,21
*NGEN,NSET=BND3
2001,2021
*NGEN,NSET=BND2
1,2001,100
*NGEN,NSET=BND4
21,2021,100
*NFILL,NSET=ALL
BND2,BND4,20
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
0.,0.,0., 300.,0.,0.
0.,300.,0.

*ELEMENT,TYPE=S4R
1,1,3,203,201,2,103,202,101
*ELGEN,ELSET=EALL,ELSET=EALL
1,10,2,1,10,200,100
*ELSET,ELSET=PRINT
 1, 
*ELSET,ELSET=ELOUT,GENERATE
601,604,1
701,704,1
801,804,1
901,904,1
*NSET,NSET=NOUT
BND2,BND3,BND4,LOADB
*NSET,NSET=PRINT
2001,
*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
BND2,YSYMM
BND3,ZSYMM
BND4,YSYMM
LOADB,2
LOADB,1
*AMPLITUDE,NAME=RAMP
0.0,0.0, 0.9,1.0, 10.0,1.0
*NSET,NSET=QA_TEST_BND3
BND3,
**
*STEP,NLGEOM=YES
*DYNAMIC,EXPLICIT
 ,10.0
*DLOAD
EALL,VP,3.4
*CLOAD, AMPLITUDE=RAMP
2001,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=BND3
U,
*ENERGY FILE
*OUTPUT, FIELD, TIMEMARKS=YES, NUMBER INTERVAL=4
*NODE OUTPUT, NSET=QA_TEST_BND3
U,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=2.5
*END STEP