*HEADING
 Cylindrical bending Pagano problem using SC8R elements
 Two layered plate
 2 elements through the thickness (1 per composite layer)
 10 elements along the length
*PREPRINT, MODEL=YES
*NODE
   101, 0., 0.
   141,10., 0.
  1701, 0., 2.5
  1741,10., 2.5
 10101, 0., 0. , -1.0
 10141,10., 0. , -1.0
 11701, 0., 2.5, -1.0
 11741,10., 2.5, -1.0
*NGEN, NSET=LEFT
   101,  1701, 100
 10101, 11701, 100
*NGEN, NSET=RIGHT
   141,  1741, 100
 10141, 11741, 100
*NSET, NSET=NCENTER
 121,
*NFILL, NSET=NALL
 LEFT, RIGHT, 40, 1
*ELSET, ELSET=LEDGE
 101, 901
*ELEMENT,TYPE=SC8R
  101, 101, 105, 10105, 10101,  901,  905, 10905, 10901
  901, 901, 905, 10905, 10901, 1701, 1705, 11705, 11701
*ELGEN, ELSET=LAYER1
  101, 10, 4, 1
*ELGEN, ELSET=LAYER2
  901, 10, 4, 1
*ELSET, ELSET=BOTTOM, GEN
 101, 140, 1
*SHELL SECTION, ELSET=LAYER1, MATERIAL=MAT1, ORIENT=ZERO, THICKNESS MODULUS=1.E6, POISSON=0.0
 1.25,7
*SHELL SECTION, ELSET=LAYER2, MATERIAL=MAT2, ORIENT=ZERO, THICKNESS MODULUS=1.E6, POISSON=0.0
 1.25,7
*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, ELSET=LEDGE, TOTALS=YES
 SF,
*EL PRINT, ELSET=LEDGE
 1,2,3,4,5,6,7
 CTSHR,
*EL FILE, ELSET=LEDGE
 SF,
*EL FILE, ELSET=LEDGE
 COORD,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 SF,
*ELEMENT OUTPUT, ELSET=LEDGE
 1,4,7
 CTSHR,
*END STEP