*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.
**
*shell section, composite, elset=shell1, nodal thickness, offset=-0.5, 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, offset=-0.5, 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, offset=-0.5, 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, offset=-0.5, 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, offset=-0.5, 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, offset=-0.5, 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
**
*nonstructural mass, elset=shell7, units=mass per area
 1.0
*nonstructural mass, elset=shell8, units=mass per area
 1.0
*nonstructural mass, elset=shell9, 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
** -----------------------------------------------------------------------------
**