*HEADING
le2_s4r_c: NAFEMS TEST LE2, Cylindrical Shell Bending Patch Test  [S4R]
*NODE
 1, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0
 5, 1.0, 30.0, 0.0, 1.0, 0.0, 0.0
 21,1.0, 0.0, 0.5, 1.0, 0.0, 0.0
 25,1.0, 30.0, 0.5, 1.0, 0.0, 0.0
 13,1.0, 20.0, 0.3, 1.0, 0.0, 0.0
*NGEN,NSET=EDGEAB
 1,21,5
*NGEN,NSET=EDGECD
 5,25,5
*NGEN,NSET=EDGEAD
 1,5
*NGEN,NSET=EDGEBC
 21,25,
*NGEN
 11,13
 13,15
 3,13,5
 13,23,5
 6,8
 8,10
 16,18
 18,20
*NSET,NSET=NALL,GEN
 1,25
*NMAP,NSET=NALL,TYPE=CYLINDRICAL
0.,0.,0.,0.,0.,1.
1.,0.,0.
*ELEMENT,TYPE=S4R,ELSET=ALLELS
 1,1,3,13,11, 2,8,12,6
*ELGEN,ELSET=ALLELS
 1,2,2,1,2,10,10
*ELSET,ELSET=HIST
12, 
*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, POISSON=0.3, ORIENTATION=ORI
0.01, 
*TRANSFORM,TYPE=C,NSET=NALL
0.,0.,0.,0.,0.,1.
*AMPLITUDE,NAME=RAMP
0.0,0.0,.25,1.0
*BOUNDARY
 EDGEAD,ZSYMM
 EDGEBC,ZSYMM
 EDGEAB,1,6
*RESTART,TIMEMARKS=NO,WRITE,NUM=10
*NSET, NSET=CHECK
3,11,15,23,13
*NSET,NSET=QA_TEST_CHECK 
CHECK ,
**
*STEP
*DYNAMIC,EXPLICIT
,.25
*CLOAD,AMP=RAMP
5,6,125.0
15,6,250.0,
25,6,125.0
***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