*Heading                                                                                            
 Verification reference for s3r elements with distributed general section stiffness
*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,
**
*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=s3r, elset=eall                                                                
        1,       1,       2,       6
        2,       2,       3,       7
        3,       3,       4,       8
        4,       5,       6,      10
        5,       6,       7,      11
        6,       7,       8,      12
        7,       9,      10,      14
        8,      10,      11,      15
        9,      11,      12,      16
       11,       6,       5,       1
       12,       7,       6,       2
       13,       8,       7,       3
       14,      10,       9,       5
       15,      11,      10,       6 
       16,      12,      11,       7
       17,      14,      13,       9
       18,      15,      14,      10 
       19,      16,      15,      11
**
*Elset, Elset=shell1
 1, 11
*Elset, Elset=shell2
 2, 12
*Elset, Elset=shell3
 3, 13
*Elset, Elset=shell4
 4, 14
*Elset, Elset=shell5
 5, 15
*Elset, Elset=shell6
 6, 16
*Elset, Elset=shell7
 7, 17
*Elset, Elset=shell8
 8, 18
*Elset, Elset=shell9
 9, 19
**
*Elset, Elset=shellSecL1
 shell1, shell2, shell3
*Elset, Elset=shellSecL2
 shell4, shell5, shell6
*Elset, Elset=shellSecL3
 shell7, shell8, shell9
**
*orientation, name=ori
-1.,   0.,   0.,   0.,   -1.,   0.
 3, 15.
**
*shell general section, elset=shell1, density=0.8E+01
   0.32967E+07, 0.98901E+06, 0.32967E+07, 0.00000E+00, 0.00000E+00, 0.11539E+07, 0.00000E+00, 0.00000E+00, 
   0.00000E+00, 0.27473E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.82418E+03, 0.27473E+04, 0.00000E+00, 
   0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.96154E+03
*shell general section, elset=shell2, density=0.8E+01
   0.49451E+07, 0.14835E+07, 0.49451E+07, 0.00000E+00, 0.00000E+00, 0.17308E+07, 0.00000E+00, 0.00000E+00, 
   0.00000E+00, 0.92720E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.27816E+04, 0.92720E+04, 0.00000E+00, 
   0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.32452E+04
*shell general section, elset=shell3, density=0.8E+01
   0.65934E+07, 0.19780E+07, 0.65934E+07, 0.00000E+00, 0.00000E+00, 0.23077E+07, 0.00000E+00, 0.00000E+00, 
   0.00000E+00, 0.21978E+05, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.65934E+04, 0.21978E+05, 0.00000E+00, 
   0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.76923E+04
**
*Shell general section, elset=shellSecL2, density=0.8E+01
   0.32967E+07, 0.98901E+06, 0.32967E+07, 0.00000E+00, 0.00000E+00, 0.11539E+07, 0.00000E+00, 0.00000E+00, 
   0.00000E+00, 0.27473E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.82418E+03, 0.27473E+04, 0.00000E+00, 
   0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.96154E+03
**
*shell general section, elset=shellSecL3, density=0.8E+01
   0.49451E+07, 0.14835E+07, 0.49451E+07, 0.00000E+00, 0.00000E+00, 0.17308E+07, 0.00000E+00, 0.00000E+00, 
   0.00000E+00, 0.92720E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.27816E+04, 0.92720E+04, 0.00000E+00, 
   0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.32452E+04
**
*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    
 1, 3, 4, 6, 7, 9, 10, 12
 s,                                                                                                 
 sf,
 ssavg,
 se,
 sth,
*Element output                                                                                     
 5,                                                                                              
 tshr,            
*Output, history                                                                                     
*Element output,elset=shellSecL1                                                                                   
 s11,                                                                                                 
 sf1,
 ssavg1,
 se1,
 sk1,
 sth,
*Element output,elset=shellSecL1                                                                    
 3,                                                                                              
 tshr13,            
*end step
**
** --------------------------------------------------
**