*HEADING
: NAFEMS TEST LE2, Cylindrical Shell Bending Patch Test  [S4R5]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*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=S4R5,ELSET=ALLELS
1,  1,3,101,15
2,  3,5,7,101
3,  101,7,9,11
4,  15,101,11,13
*MATERIAL,NAME=MAT
*ELASTIC,TYPE=ISOTROPIC
2.1E5,0.3
*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.
*BOUNDARY
EDGEAB,1,6
EDGEAD,ZSYMM
EDGEBC,ZSYMM
*STEP,PERTURBATION
*STATIC
*CLOAD
 1,6,0.125E3
 3,6,0.25E3
 5,6,0.125E3
*NODE PRINT
*EL FILE,POSITION=AVERAGED AT NODES
 S,
 E,
*EL PRINT,POSITION=AVERAGED AT NODES
 S,
 E,
*END STEP