*Heading
 Verification test for s4 elements with distributed general section stiffness
*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=s4, 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=distTableStiff
 SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF2,SHELLSTIFF2,
 SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF3,SHELLSTIFF2,
 SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF3,SHELLSTIFF3
**
*Distribution, name=distStiff, table=distTableStiff, location=element
        shell2, 0.49451E+07, 0.14835E+07, 0.49451E+07, 0.00000E+00, 0.00000E+00, 0.17308E+07, 0.00000E+00,
        shell3, 0.65934E+07, 0.19780E+07, 0.65934E+07, 0.00000E+00, 0.00000E+00, 0.23077E+07, 0.00000E+00,
**
*shell general section, elset=shellSecL1, density=0.8E+01
 distStiff
**
*Shell general section, elset=shellSecL2, density=0.8E+01
   0.32967E+07, 0.98901E+06, 0.32967E+07, 0.00000E+00, 0.00000E+00, 0.11539E+07, 0.00000E+00, 0.00000E+00,
   0.00000E+00, 0.27473E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.82418E+03, 0.27473E+04, 0.00000E+00,
   0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.96154E+03
**
*shell general section, elset=shellSecL3, density=0.8E+01
   0.49451E+07, 0.14835E+07, 0.49451E+07, 0.00000E+00, 0.00000E+00, 0.17308E+07, 0.00000E+00, 0.00000E+00,
   0.00000E+00, 0.92720E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.27816E+04, 0.92720E+04, 0.00000E+00,
   0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.32452E+04
**
*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=shellSecL1
 tshr13,
*end step
**
** --------------------------------------------------
**