*HEADING
: NAFEMS TEST LE2, Cylindrical Shell Bending Patch Test  [STRI3]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
 1,    1.,  10.,  0.0
 5,    1.,  10.,  0.5
 9,    1.,   0.,  0.5
 13,   1.,   0.,  0.0
 17,   1.,  10.,  0.0
 101,  1.,  6.6667,  0.3
*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=STRI3,ELSET=ALLELS
1, 101,15,1
2, 101,1,3
3, 101,3,5
4, 101,5,7
5, 101,7,9
6, 101,9,11
7, 101,11,13
8, 101,13,15
*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,FREQ=0
*EL FILE,POSITION=AVERAGED AT NODES
 S,
 E,
*EL PRINT,POSITION=AVERAGED AT NODES
 S,
 E,
*END STEP