*HEADING
PINCHED CYLINDER WITH DIAPHRAGMS AT ENDS - S4R - 20 X 20 MESH
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*RESTART,TIMEMARKS=YES,WRITE,NUMBER INTERVAL=1
*NODE,NSET=CORNER
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,2,102,101
*ELGEN,ELSET=EALL
1,20,1,1,20,100,100
*ELSET,ELSET=PRINT
 1, 
*ELSET,ELSET=ELOUT,GENERATE
1401,1406,1
1501,1506,1
1601,1606,1
1701,1706,1
1801,1806,1
1901,1906,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.5
*CLOAD
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