*HEADING
le2_s4rs_f: NAFEMS TEST LE2, Cylindrical Shell Bending Patch Test  [S4RS]
*NODE
 1,    1.,  10.,  0.0,   1.,  0.,  0.
 5,    1.,  10.,  0.5,   1.,  0.,  0.
 9,    1.,   0.,  0.5,   1.,  0.,  0.
 13,   1.,   0.,  0.0,   1.,  0.,  0.
 17,   1.,  10.,  0.0,   1.,  0.,  0.
 101,  1.,  6.66667,  0.3,   1.,  0.,  0.
*NGEN,NSET=EDGECD
 1,5
*NGEN,NSET=EDGEBC
 5,9
*NGEN,NSET=EDGEAB
 9,13
*NGEN,NSET=EDGEAD1
 13,17
*NSET,NSET=EDGEAD
 1,13,14,15,16
*NSET,NSET=OUTER
 EDGECD,EDGEBC,EDGEAB,EDGEAD
*NGEN,NSET=INNER
 1,101,50
 3,101,49
 5,101,48
 7,101,47
 9,101,46
 11,101,45
 13,101,44
 15,101,43
*NSET,NSET=ALL
 INNER,OUTER
*NMAP,NSET=ALL,TYPE=CYLINDRICAL
0.,0.,0.,0.,0.,1.
1.,0.,0.
*ELEMENT,TYPE=S4RS,ELSET=ALLELS
1,  1,3,101,15
2,  3,5,7,101
3,  101,7,9,11
4,  15,101,11,13
*ELSET,ELSET=HIST
2, 
*MATERIAL,NAME=MAT
*ELASTIC,TYPE=ISOTROPIC
210.E9,0.3
*DENSITY
7800., 
*ORIENTATION, NAME=ORI, SYSTEM=CYLINDRICAL
 0.,0.,0., 0.,0.,1.
 1,0.
*SHELL SECTION, MATERIAL=MAT, ELSET=ALLELS, ORIENTATION=ORI
0.01, 
*TRANSFORM,TYPE=C,NSET=ALL
0.,0.,0.,0.,0.,1.
*AMPLITUDE,NAME=RAMP
0.0,0.0,.25,1.0
*BOUNDARY
EDGEAB,1,6
EDGEAD,ZSYMM
EDGEBC,ZSYMM
*RESTART,TIMEMARKS=NO,WRITE,NUM=10
*NSET, NSET=CHECK
11,7,3,15, 101
*NSET,NSET=QA_TEST_CHECK
CHECK,
**
*STEP
*DYNAMIC,EXPLICIT
,.25
*CLOAD,AMP=RAMP
 1,6,0.125E3
 3,6,0.25E3
 5,6,0.125E3
***OUTPUT,HISTORY,TIME INTERVAL=.005
***ELEMENT OUTPUT,ELSET=HIST,SECTIONPOINT=5
**S,
*FILE OUTPUT,TIMEMARKS=YES, NUM=1
*NODE FILE, NSET=CHECK
U, 
*ENERGY FILE
*OUTPUT, FIELD, TIMEMARKS=YES, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_CHECK
U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.25
*END STEP