*Heading                                                                                            
 Multistep reference solution test for s4 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=eall, generate                                                                        
        1,       9,       1,
*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=s4, elset=shellSecL1                                                                 
        1,       1,       2,       6,       5,
        4,       5,       6,      10,       9,
        7,       9,      10,      14,      13,
*Shell Section, orient=orientL1, material=mat1, elset=shellSecL1, offset=-0.50000E+00, nodal thickness             
  1.,
*Orientation, name=orientL1                                                                         
   0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00
  3,0.
**------------------------------------------------------------------------------------------        
*Element, type=s4, elset=shellSecL2                                                                 
        2,       2,       3,       7,       6,
        5,       6,       7,      11,      10,
        8,      10,      11,      15,      14,
*Shell Section, orient=orientL2, material=mat1, elset=shellSecL2, offset=-0.50000E+00               
   0.01,
*Orientation, name=orientL2                                                                         
  0.0000000E+00,   0.1000000E+01,   0.0000000E+00, -0.1000000E+01,   0.0000000E+00,   0.0000000E+00
  3,0.
**------------------------------------------------------------------------------------------        
*Element, type=s4, elset=shellSecL3                                                                 
        3,       3,       4,       8,       7,
        6,       7,       8,      12,      11,
        9,      11,      12,      16,      15,
*Shell Section, orient=orientL3, material=mat1, elset=shellSecL3, offset=-0.50000E+00, nodal thickness                
   1.,
*Orientation, name=orientL3             
   -0.1000000E+01,   0.0000000E+00,   0.0000000E+00, 0.0000000E+00,   -0.1000000E+01,   0.0000000E+00
  3,0.
*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.1
*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
**