*HEADING
 Cylindrical bending Pagano problem using SC8R elements
 Three layered plate
 1 elements through the thickness (composite section)
 10 elements along the length
*PREPRINT, MODEL=YES
*NODE
   101, 0., 0.
   141,10., 0.
  2501, 0., 2.5
  2541,10., 2.5
 10101, 0., 0. , -1.0
 10141,10., 0. , -1.0
 12501, 0., 2.5, -1.0
 12541,10., 2.5, -1.0
*NGEN, NSET=LEFT
   101,  2501, 100
 10101, 12501, 100
*NGEN, NSET=RIGHT
   141,  2541, 100
 10141, 12541, 100
*NSET, NSET=NCENTER
 121,
*NFILL, NSET=NALL
 LEFT, RIGHT, 40, 1
*ELSET, ELSET=LEDGE, GEN
 101, 2501, 100
*ELEMENT,TYPE=SC8R
  101, 101, 105, 10105, 10101, 2501, 2505, 12505, 12501
*ELGEN, ELSET=LAYER1
  101, 10, 4, 1
*ELSET, ELSET=BOTTOM, GEN
 101, 140, 1
*SHELL SECTION, ELSET=LAYER1, composite, THICKNESS MODULUS=1.D6, POISSON=0.0
 0.8333333,11,mat1
 0.8333333,11,mat2
 0.8333333,11,mat1
*MATERIAL, NAME=MAT1
*ELASTIC,TYPE=LAMINA
** E1,  E2, n12, G12 , G13 , G23
25.E6,1.E6,0.25,0.5E6,0.5E6, 0.2E6
*MATERIAL, NAME=MAT2
*ELASTIC,TYPE=LAMINA
** E1,  E2 , n12, G12 , G13 , G23
 1.E6,25.E6,0.01,0.5E6,0.2E6, 0.5E6
*ORIENTATION, NAME=ZERO
 1.,0.,0., 0.,0.,-1.
 3,0. 
*BOUNDARY
  NALL, 3,3
  LEFT, 2,2
 RIGHT, 2,2
 NCENTER, 1,1
*STEP,PERTURBATION
*STATIC
*DLOAD
 BOTTOM, P1NU
*EL PRINT, POSITION=AVERAGED AT NODES
 S,
*EL FILE, POSITION=AVERAGED AT NODES
 S,
*EL PRINT, ELSET=LEDGE
 SF,
*EL PRINT, ELSET=LEDGE
 1,2,3,4,5,6,7,8,9,10,11
 TSHR,
*EL PRINT, ELSET=LEDGE
 12,13,14,15,16,17,18,19,20,21,22
 TSHR,
*EL PRINT, ELSET=LEDGE
 23,24,25,26,27,28,29,30,31,32,33
 TSHR,
*EL FILE, ELSET=LEDGE
 1,2,3,4,5,6,7,8,9,10,11
 TSHR,
*EL FILE, ELSET=LEDGE
 12,13,14,15,16,17,18,19,20,21,22
 TSHR,
*EL FILE, ELSET=LEDGE
 23,24,25,26,27,28,29,30,31,32,33
 TSHR,
*EL PRINT, ELSET=LEDGE
 SF,
 COORD,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 SF,
*ELEMENT OUTPUT
 1,2,3,4,5,6,7,8,9,10,11
 TSHR,
*ELEMENT OUTPUT
 12,13,14,15,16,17,18,19,20,21,22
 TSHR,
*ELEMENT OUTPUT
 23,24,25,26,27,28,29,30,31,32,33
 TSHR,
*END STEP