*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,
**------------------------------------------------------------------------------------------        
*Element, type=s4, elset=shellSecL1                                                                 
        1,       1,       2,       6,       5,
        4,       5,       6,      10,       9,
        7,       9,      10,      14,      13,
*Shell Section, orient=orientL1, composite, elset=shellSecL1, offset=-0.50000E+00, nodal thickness             
 1.0, 3, mat1, 90, ply-1
 1.0, 3, mat1, 0, ply-2
*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, composite, elset=shellSecL2, offset=-0.50000E+00               
 0.005, 3, mat1, 90, ply-1
 0.005, 3, mat1, 0, ply-2
*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, composite, elset=shellSecL3, offset=-0.50000E+00, nodal thickness                
 1.0, 3, mat1, 90, ply-1
 1.0, 3, mat1, 0, ply-2
*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            
** --------------------------------------------------
**
*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                                                                                     
 3,                                                                                              
 tshr,            
*Output, history                                                                                     
*Element output,elset=shellSecL1                                                                                   
 s11,                                                                                                 
 sf1,
 ssavg1,
 se1,
 sk1,
 sth,
*Element output,elset=shellSecL1                                                                    
 3,                                                                                              
 tshr13,            
*end step
**
** --------------------------------------------------
**