*Heading
 Verification test for s3r 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=s3r, elset=eall                                                                
        1,       1,       2,       6
        2,       2,       3,       7
        3,       3,       4,       8
        4,       5,       6,      10
        5,       6,       7,      11
        6,       7,       8,      12
        7,       9,      10,      14
        8,      10,      11,      15
        9,      11,      12,      16
       11,       6,       5,       1
       12,       7,       6,       2
       13,       8,       7,       3
       14,      10,       9,       5
       15,      11,      10,       6 
       16,      12,      11,       7
       17,      14,      13,       9
       18,      15,      14,      10 
       19,      16,      15,      11
**
*Elset, Elset=shell1
 1, 11
*Elset, Elset=shell2
 2, 12
*Elset, Elset=shell3
 3, 13
*Elset, Elset=shell4
 4, 14
*Elset, Elset=shell5
 5, 15
*Elset, Elset=shell6
 6, 16
*Elset, Elset=shell7
 7, 17
*Elset, Elset=shell8
 8, 18
*Elset, Elset=shell9
 9, 19
**
*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, offset=-0.5, 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, offset=-0.5, 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
**
*nonstructural mass, elset=shellSecL3, units=mass per area
 1.0
**
*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
** -----------------------------------------------------------------------------
**