*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