*Heading
 Verification test for s4r elements with distributed layer angle and thickness
*Node
*Nset, nset=side1, generate
*Nset, nset=side3, generate
*Nset, nset=nodeThick1, generate
*Nset, nset=nodeThick2, generate
*Nodal Thickness
 nodeThick1, 0.01
 nodeThick2, 0.02
**------------------------------------------------------------------------------------------
*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
**
*Elset, Elset=shell1
*Elset, Elset=shell2
*Elset, Elset=shell3
*Elset, Elset=shell4
*Elset, Elset=shell5
*Elset, Elset=shell6
*Elset, Elset=shell7
*Elset, Elset=shell8
*Elset, Elset=shell9
**
*Elset, Elset=shellSecL1
 shell1, shell2, shell3
*Elset, Elset=shellSecL2
 shell4, shell5, shell6
*Elset, Elset=shellSecL3
 shell7, shell8, shell9
**
*orientation, name=ori
**
*Distribution Table, name=thicknessTable
 length
*Distribution, name=distThickPly1, location=element, table=thicknessTable
 shell4,   0.175
 shell6,   0.165
*Distribution, name=distThickPly2, location=element, table=thicknessTable
 shell7,   1.250
 shell8,   1.100
*Distribution, name=distThickPly4, location=element, table=thicknessTable
 shell5,   0.095
 shell6,   0.075
**
*Distribution Table, name=orientAngTable
 angle
*Distribution, name=distAngPly1, location=element, table=orientAngTable
shell2,  40.
shell3,  70.
*Distribution, name=distAngPly2, location=element, table=orientAngTable
shell7,  60.
shell9,  90.
*Distribution, name=distAngPly3, location=element, table=orientAngTable
shell4, -80.
shell5, -70.
*Distribution, name=distAngPly4, location=element, table=orientAngTable
shell7, -60.
shell8, -90.
**
*Shell general Section, composite, elset=shellSecL1, nodal thickness, offset=-0.5, orientation=ori
        1.000, 3, mat1, distAngPly1,   ply-1
        1.000, 3, mat1,  78.,          ply-2
        1.000, 3, mat1,  18.,          ply-3
        1.000, 3, mat1,  28.,          ply-4
**
*Shell general Section, composite, elset=shellSecL2, offset=-0.5, orient=ori
distThickPly1, 3, mat1,  38.,          ply-1
        1.250, 3, mat1,  78.,          ply-2
        1.250, 3, mat1, distAngPly3,   ply-3
distThickPly4, 3, mat1,  58.,          ply-4
**
*Shell general 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
 tshr,
*Output, history
*Element output,elset=shellSecL1
 s11,
 sf1,
 ssavg1,
 se1,
 sk1,
 sth,
*Element output,elset=shellSecL2
 tshr13,
*end step
**
** --------------------------------------------------
**