*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.
**
*Distribution Table, name=thicknessTable
 length
*Distribution, name=distThickPly1, location=element, table=thicknessTable
       ,   0.150
 shell4,   0.175
 shell6,   0.165
*Distribution, name=distThickPly2, location=element, table=thicknessTable
       ,   1.400
 shell7,   1.250
 shell8,   1.100
*Distribution, name=distThickPly4, location=element, table=thicknessTable
       ,   0.130
 shell5,   0.095
 shell6,   0.075
**
*Distribution Table, name=orientAngTable
 angle
*Distribution, name=distAngPly1, location=element, table=orientAngTable
      ,  10.
shell2,  40.
shell3,  70.
*Distribution, name=distAngPly2, location=element, table=orientAngTable
      ,  30.
shell7,  60.
shell9,  90.
*Distribution, name=distAngPly3, location=element, table=orientAngTable
      , -20.
shell4, -80.
shell5, -70.
*Distribution, name=distAngPly4, location=element, table=orientAngTable
      , -40.
shell7, -60.
shell8, -90.
**
*shell section, composite, elset=shellSecL1, nodal thickness, orientation=ori
        1.000, 3, mat1, distAngPly1,   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=shellSecL2, orient=ori      
distThickPly1, 3, mat1,  38.,          ply-1
        1.250, 3, mat1, ori_ply,       ply-2
        1.250, 3, mat1, distAngPly3,   ply-3
distThickPly4, 3, mat1,  58.,          ply-4
**
*shell section, composite, elset=shellSecL3, nodal thickness, orientation=ori
        1.000, 3, mat1,  68.,          ply-1
distThickPly2, 3, mat1, distAngPly2,   ply-2
        1.000, 3, mat1,  48.,          ply-3
        1.000, 3, mat1, distAngPly4,   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
** -----------------------------------------------------------------------------
**