*Heading                                                                                            
 Multistep verification test for s4 elements using distributions    
*Preprint, model=yes                                
*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,
*Elset, Elset=eall, generate                                                                        
        1,       9,       1,
*Elset, elset=shellSecL1  
        1,
**------------------------------------------------------------------------------------------        
*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,
*Shell General Section, elset=eall, density=.04, orientation=ori                                          
   distStiff,
*Distribution, name=distOrient, location=element, table=distTableOri, input=distrib_s4_orient_in.inp                               
*Distribution, name=distStiff, location=element, table=distTableStiff, input=distrib_s4_stiff_in.inp    
*Distribution Table, name=distTableStiff
             SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF1,SHELLSTIFF2,
 SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF3,
 SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF2,SHELLSTIFF3,SHELLSTIFF3,SHELLSTIFF3
*Distribution Table, name=distTableOri
 COORD3D, COORD3D, 
*Orientation, name=ori
 distOrient
 3, 0.
*Boundary                                                                                           
 side1, 1, 6                                                                                        
*psd-definition, name=white, type=force
 1.0,0.1,0.0
 1.0,0.1,60.0
*spectrum, name=spec,type=acceleration
 0., 0., 0.
 1.e5, 10., 0.
 1.e5, 100., 0.
**
** --------------------------------------------------
**
*step, name=step-1
 frequency analysis
*frequency, SIM=NO
 5,
*El print
 sf,
 se,
*Node print
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 sf,
*end step
**
** --------------------------------------------------
**
*step, name=step-2
 SSD mode-based
*steady state dynamics, interval=eigenfrequency
 60,
*El print, freq=999
 sf,
 se,
*Node print, freq=999
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 sf,
*cload
 side3, 4, 0.1
 side3, 3, 2.0
 side3, 2, 0.5
*select eigenmodes, generate
 1, 5
*modal damping, rayleigh
 1, 5, 0.5 
*end step
**
** --------------------------------------------------
**
*step, name=step-3
 modal dynamic
*modal dynamic
 0.1,1.
*cload
 side3, 4, 0.1
 side3, 3, 0.2
 side3, 2, 0.5
*select eigenmodes, generate
 1, 5
*modal damping, rayleigh
 1, 5, 0.5 
*end step
**
** --------------------------------------------------
**
*step, name=step-4
 SSD direct
*steady state dynamics, direct, interval=eigenfrequency
 60.,
*cload
 side3, 4, 0.1
 side3, 3, 0.2
 side3, 2, 0.5
*end step
**
** --------------------------------------------------
**
*step, name=step-5
 SSD subspace
*steady state dynamics, subspace projection, interval=eigenfrequency
 60.,
*cload
 side3, 4, 0.1
 side3, 3, 0.2
 side3, 2, 0.5
*end step
**
** --------------------------------------------------
**
*step, name=step-6
 random response
*random response
 60.,
*select eigenmodes, generate
 1, 5
*cload, load case=1
 side3, 4, 0.1
 side3, 3, 0.2
 side3, 2, 0.5
*cload, load case=2
 side3, 4, 0.1
 side3, 3, 0.2
 side3, 2, 0.5
*correlation, psd=white, type=uncorrelated, complex=no
 1,1
 2,1
*end step
**
** --------------------------------------------------
**
*step, name=step-7
 response spectrum
*response spectrum, sum=abs, comp=alg
 spec,0.,0.,1.,1.
*select eigenmodes, generate
 1, 5
*end step
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-8
 static
*static
 1.0, 1.0
*Boundary
 side3, 2, 2, 0.025
 side3, 3, 3, 0.050
*El print, freq=999                                                                                 
 sf,
 se,
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 sf,
*Output, history                                                                                     
*Element output,elset=shellSecL1                                                                                   
 sf1,
 se1,
 sk1,
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-9
 static 
*static
*El print, freq=999                                                                                 
 sf,
 se,
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 sf,
*load case, name=load1
*dload
 eall, bz, 1.e3                                                                             
*end load case
*load case, name=load2
*dload
 eall, grav, 1.e6, 0.,0.,1.                                                                                
*end load case
*load case, name=load3
*dload
 eall, cent, 100.,0.,0.,0.,1.,0.                                                                                 
*end load case
*load case, name=load4
*dload
 eall, centrif, 1.e6,0.,0.,0.,1.,0.                                                                                 
*end load case
*end step