*Heading
 Verification test for s4r elements with distributed layer angle and thickness
*Node
        1,  0.0000000E+00,  0.0000000E+00,  0.0000000E+00
        2,  0.3333333E+00,  0.0000000E+00,  0.0000000E+00
        3,  0.6666667E+00,  0.0000000E+00,  0.0000000E+00
        4,  0.1000000E+01,  0.0000000E+00,  0.0000000E+00
        5,  0.0000000E+00,  0.3333333E+00,  0.0000000E+00
        6,  0.3333333E+00,  0.3333333E+00,  0.0000000E+00
        7,  0.6666667E+00,  0.3333333E+00,  0.0000000E+00
        8,  0.1000000E+01,  0.3333333E+00,  0.0000000E+00
        9,  0.0000000E+00,  0.6666667E+00,  0.0000000E+00
       10,  0.3333333E+00,  0.6666667E+00,  0.0000000E+00
       11,  0.6666667E+00,  0.6666667E+00,  0.0000000E+00
       12,  0.1000000E+01,  0.6666667E+00,  0.0000000E+00
       13,  0.0000000E+00,  0.1000000E+01,  0.0000000E+00
       14,  0.3333333E+00,  0.1000000E+01,  0.0000000E+00
       15,  0.6666667E+00,  0.1000000E+01,  0.0000000E+00
       16,  0.1000000E+01,  0.1000000E+01,  0.0000000E+00
*Nset, nset=side1, generate
        1,       4,       1,
*Nset, nset=side3, generate
       13,      16,       1,
*Nset, nset=nodeThick1, generate
        1,      13,       4,
        4,      16,       4,
*Nset, nset=nodeThick2, generate
        2,      14,       4,
        3,      15,       4,
**
*Nodal Thickness
 nodeThick1, 0.5
 nodeThick2, 0.8
** -----------------------------------------------------------------------------
*Material, name=mat1
*Elastic, type=engineering constants
  0.300E+07,  0.400E+07,  0.500E+07,  0.100E+00,  0.200E+00,  0.300E+00,  0.300E+06,  0.400E+06
  0.500E+06,
*Density
  0.750E+01,
** -----------------------------------------------------------------------------
*Element, type=s4r, elset=eall
        1,       1,       2,       6,       5,
        2,       2,       3,       7,       6,
        3,       3,       4,       8,       7,
        4,       5,       6,      10,       9,
        5,       6,       7,      11,      10,
        6,       7,       8,      12,      11,
        7,       9,      10,      14,      13,
        8,      10,      11,      15,      14,
        9,      11,      12,      16,      15,
**
*Elset, Elset=shell1
 1,
*Elset, Elset=shell2
 2,
*Elset, Elset=shell3
 3,
*Elset, Elset=shell4
 4,
*Elset, Elset=shell5
 5,
*Elset, Elset=shell6
 6,
*Elset, Elset=shell7
 7,
*Elset, Elset=shell8
 8,
*Elset, Elset=shell9
 9,
**
*Elset, Elset=shellSecL1
 shell1, shell2, shell3
*Elset, Elset=shellSecL2
 shell4, shell5, shell6
*Elset, Elset=shellSecL3
 shell7, shell8, shell9
** 
*orientation, name=ori
-1.,   0.,   0.,   0.,   -1.,   0.
 3, 15.
*orientation, name=ori_ply
 1.,   1.,   0.,  -1.,    1.,   0.
 3, 78.
**
*shell section, composite, elset=shell1, nodal thickness, orientation=ori
        1.000, 3, mat1,  10.,       ply-1
        1.000, 3, mat1, ori_ply,    ply-2
        1.000, 3, mat1,  18.,       ply-3
        1.000, 3, mat1,  28.,       ply-4
*shell section, composite, elset=shell2, nodal thickness, orientation=ori
        1.000, 3, mat1,  40.,       ply-1
        1.000, 3, mat1, ori_ply,    ply-2
        1.000, 3, mat1,  18.,       ply-3
        1.000, 3, mat1,  28.,       ply-4
*shell section, composite, elset=shell3, nodal thickness, orientation=ori
        1.000, 3, mat1,  70.,       ply-1
        1.000, 3, mat1, ori_ply,    ply-2
        1.000, 3, mat1,  18.,       ply-3
        1.000, 3, mat1,  28.,       ply-4
**
*Shell Section, composite, elset=shell4, orient=ori      
        0.175, 3, mat1,  38.,       ply-1
        1.250, 3, mat1, ori_ply,    ply-2
        1.250, 3, mat1, -80.,       ply-3
        0.130, 3, mat1,  58.,       ply-4
*Shell Section, composite, elset=shell5, orient=ori      
        0.150, 3, mat1,  38.,       ply-1
        1.250, 3, mat1, ori_ply,    ply-2
        1.250, 3, mat1, -70.,       ply-3
        0.095, 3, mat1,  58.,       ply-4
*Shell Section, composite, elset=shell6, orient=ori      
        0.165, 3, mat1,  38.,       ply-1
        1.250, 3, mat1, ori_ply,    ply-2
        1.250, 3, mat1, -20.,       ply-3
        0.075, 3, mat1,  58.,       ply-4
**
*shell section, composite, elset=shell7, nodal thickness, orientation=ori
        1.000, 3, mat1,  68.,       ply-1
        1.250, 3, mat1,  60.,       ply-2
        1.000, 3, mat1,  48.,       ply-3
        1.000, 3, mat1, -60.,       ply-4
*shell section, composite, elset=shell8, nodal thickness, orientation=ori
        1.000, 3, mat1,  68.,       ply-1
        1.100, 3, mat1,  30.,       ply-2
        1.000, 3, mat1,  48.,       ply-3
        1.000, 3, mat1, -90.,       ply-4
*shell section, composite, elset=shell9, nodal thickness, orientation=ori
        1.000, 3, mat1,  68.,       ply-1
        1.400, 3, mat1,  90.,       ply-2
        1.000, 3, mat1,  48.,       ply-3
        1.000, 3, mat1, -40.,       ply-4
**
*Boundary
 side1, 1, 6
**
*amplitude, def=smooth step, name=smooth
0.,  0., 1., 1.0
** -----------------------------------------------------------------------------
*step, nlgeom=yes, name=step-8
 explicit
*dynamic, explicit
 , 1.0
*Boundary, amplitude=smooth
 side3, 3, 3, 0.050
*Output, field
*Node output
 u,
*Element output
 s,
 sf,
 ssavg,
 se,
 sth,
*Element output
 3, 6, 9
 tshr,
*Output, history
*Element output,elset=shellSecL1
 s11,
 sf1,
 ssavg1,
 se1,
 sk1,
 sth,
*Element output,elset=shellSecL1
 3, 6, 9
 tshr13,
*end step
** -----------------------------------------------------------------------------
**