*Heading                                                                                            
 Multistep verification 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=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,
**------------------------------------------------------------------------------------------        
*Element, type=s4, 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, composite, elset=midEl, offset=-0.50000E+00, orient=ori, shell thickness=distThick             
 1.0, 3, mat1, 90, ply-1
 1.0, 3, mat1, 0, ply-2
*shell section, composite, elset=sideEl, nodal thickness, offset=-0.5, orientation=ori
 1.0, 3, mat1, 90, ply-1
 1.0, 3, mat1, 0, ply-2
*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
**
** --------------------------------------------------
**