*HEADING
 NAFEMS TEST LE2, Cylindrical Shell Bending Patch Test  [SC6R]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
 1,    0.995,      10.,  0.0,   1.,  0.,  0.
 5,    0.995,      10.,  0.5,   1.,  0.,  0.
 9,    0.995,       0.,  0.5,   1.,  0.,  0.
 13,   0.995,       0.,  0.0,   1.,  0.,  0.
 17,   0.995,      10.,  0.0,   1.,  0.,  0.
 101,  0.995,  6.66667,  0.3,   1.,  0.,  0.
 1001, 1.005,      10.,  0.0,   1.,  0.,  0.
 1005, 1.005,      10.,  0.5,   1.,  0.,  0.
 1009, 1.005,       0.,  0.5,   1.,  0.,  0.
 1013, 1.005,       0.,  0.0,   1.,  0.,  0.
 1017, 1.005,      10.,  0.0,   1.,  0.,  0.
 1101, 1.005,  6.66667,  0.3,   1.,  0.,  0.
*NGEN,NSET=EDGECD
 1,5
 1001,1005
*NGEN,NSET=EDGEBC
 5,9
 1005,1009
*NGEN,NSET=EDGEAB
 9,13
 1009,1013
*NGEN,NSET=EDGEAD1
 13,17
 1013,1017
*NSET,NSET=EDGEAD
 1,13,14,15,16
 1001,1013,1014,1015,1016
*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
 1001,1101,50
 1003,1101,49
 1005,1101,48
 1007,1101,47
 1009,1101,46
 1011,1101,45
 1013,1101,44
 1015,1101,43
*NSET,NSET=ALL
 INNER,OUTER
*NMAP,NSET=ALL,TYPE=CYLINDRICAL
0.,0.,0.,0.,0.,1.
1.,0.,0.
*ELEMENT,TYPE=SC6R,ELSET=ALLELS
 1,1,3,101,1001,1003,1101
 1001,1,101,15,1001,1101,1015
 2,3,5,7,1003,1005,1007
 1002,3,7,101,1003,1007,1101
 3,101,7,9,1101,1007,1009
 1003,101,9,11,1101,1009,1011
 4,15,101,11,1015,1101,1011
 1004,15,11,13,1015,1011,1013
*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,3
EDGEAD,ZSYMM
EDGEBC,ZSYMM
*STEP,PERTURB
*STATIC
*CLOAD
 1,2,-0.125E5
 3,2,-0.25E5
 5,2,-0.125E5
 1001,2,+0.125E5
 1003,2,+0.25E5
 1005,2,+0.125E5
*OUTPUT, FIELD
*ELEMENT OUTPUT
 S,
*OUTPUT, FIELD
*ELEMENT OUTPUT, POSITION=NODES
 S,MISESONLY,PRESSONLY
*END STEP