*Heading                                                                                            
 Multistep verification test for s4r elements with nodal 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,
*Elset, Elset=shellSecL1, generate
        1,       7,       3,
*Elset, Elset=shellSecL2, generate
        2,       8,       3,
*Elset, Elset=shellSecL3, generate
        3,       9,       3,
*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                                                                
        1,       1,       2,       6,       5,
        2,       2,       3,       7,       6,
        3,       3,       4,       8,       7,
        4,       5,       6,      10,       9,
        5,       6,       7,      11,      10,
        6,       7,       8,      12,      11,
        7,       9,      10,      14,      13,
        8,      10,      11,      15,      14,
        9,      11,      12,      16,      15,
*Distribution Table, name=orientationTable
 coord3d, coord3d
*Distribution, name=distOrient, location=element, table=orientationTable
 , 1., 0., 0., 0., 1.0, 0.0       
 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
 3, 0.0
*Distribution Table, name=ThicknessTable
 length
*Distribution, name=distThick, location=element, table=thicknessTable
  , 1.0     
 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                                                                                
 3,                                                                                              
 tshr,            
*Output, hist, filter=antialiasing, Time Int=0.001
*Element Output, Elset=QA_TEST
 s11,                                                                                                 
 sf1,
 ssavg1,
 se1,
 sk1,
 sth,
*Element output, Elset=QA_TEST                                                                  
 3,                                                                                              
 tshr13,            
*Node Output, Nset=QA_TEST
RF3,
*Energy output, VAR=PRESELECT
*end step
**