*Heading
 Multistep verification test for s4r elements with nodal thickness
*Node
*Nset, nset=side1, generate
*Nset, nset=side3, generate
*Nset, nset=nodeThick1, generate
*Nset, nset=nodeThick2, generate
*Elset, Elset=shellSecL1, generate
*Elset, Elset=shellSecL2, generate
*Elset, Elset=shellSecL3, generate
*elset, elset=sideEl
 shellSecL1, shellSecL3
*elset, elset=midEl
 shellSecL2,
*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,
**  0.750E-03,
**------------------------------------------------------------------------------------------
*Element, type=s4r, elset=eall
*Distribution Table, name=orientationTable
 coord3d, coord3d
*Distribution, name=distOrient, location=element, table=orientationTable
 shellSecL1,  0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00
 shellSecL2,  0.0000000E+00,   0.1000000E+01,   0.0000000E+00,  -0.1000000E+01,   0.0000000E+00,   0.0000000E+00
 shellSecL3, -0.1000000E+01,   0.0000000E+00,   0.0000000E+00,  0.0000000E+00,   -0.1000000E+01,   0.0000000E+00
*orientation, name=ori
 distOrient
*Distribution Table, name=ThicknessTable
 length
*Distribution, name=distThick, location=element, table=thicknessTable
 shellSecL2, 0.01
*Shell Section, material=mat1, elset=midEl, offset=-0.50000E+00, orient=ori, shell thickness=distThick
 1.0,
*shell section, material=mat1, elset=sideEl, nodal thickness, offset=-0.5, orientation=ori
 1.,
*Boundary
 side1, 1, 6
**
** --------------------------------------------------
**
**
*Elset, Elset =QA_TEST
eall,
*Nset, Nset =QA_TEST
Side3,
**
*step, nlgeom=yes, name=step-8
 dynamic
*Dynamic, Explicit
, 1.0
*Amplitude, Name=Smooth, Def=Smooth Step
0., 0., 0.5, 1.0, 1., 1.
*Boundary, Amp=Smooth
 side3, 3, 3, 0.050
*Output, field, op=new, Time Int = 0.25, Time Marks =Yes
*Node output, Nset =QA_TEST
 RF,
*Node output
 u,
*Element output, Elset =QA_Test
 s,
 sf,
 ssavg,
 se,
 sth,
*Element output, Elset =QA_Test
 tshr,
*Output, hist, filter=antialiasing, Time Int=0.1
*Element Output, Elset=QA_TEST
 s11,
 sf1,
 ssavg1,
 se1,
 sk1,
 sth,
*Element output, Elset=QA_TEST
 tshr13,
*Node Output, Nset=QA_TEST
RF3,
*Energy output, VAR=PRESELECT
*end step
**