*HEADING
 Cylindrical bending Pagano problem using SC8R elements
 Three layered plate
 6 elements through the thickness (2 per composite layer)
 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,  501,  505, 10505, 10501
  901, 901, 905, 10905, 10901, 1301, 1305, 11305, 11301
 1701,1701,1705, 11705, 11701, 2101, 2105, 12105, 12101
*ELGEN, ELSET=LAYER1
  101, 10, 4, 1, 2, 400, 400
*ELGEN, ELSET=LAYER2
  901, 10, 4, 1, 2, 400, 400
*ELGEN, ELSET=LAYER3
 1701, 10, 4, 1, 2, 400, 400
*ELSET, ELSET=BOTTOM, GEN
 101, 140, 1
*SHELL SECTION, ELSET=LAYER1, MATERIAL=MAT1, ORIENT=ZERO, THICKNESS MODULUS=1.E6, POISSON=0.0
 0.8333333,3
*SHELL SECTION, ELSET=LAYER2, MATERIAL=MAT2, ORIENT=ZERO, THICKNESS MODULUS=1.E6, POISSON=0.0
 0.8333333,3
*SHELL SECTION, ELSET=LAYER3, MATERIAL=MAT1, ORIENT=ZERO, THICKNESS MODULUS=1.E6, POISSON=0.0
 0.8333333,3
*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
 SF,
*EL PRINT, ELSET=LEDGE
 COORD,
*EL FILE, ELSET=LEDGE
 SF,
*EL PRINT, ELSET=LEDGE
 COORD,
*EL PRINT, ELSET=LEDGE
 1,2,3
 CTSHR13,
 SSAVG4,
*EL FILE, ELSET=LEDGE
 COORD,
*OUTPUT, FIELD
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 SF,
*ELEMENT OUTPUT, ELSET=LEDGE
 1,2,3
 CTSHR,
*END STEP