*HEADING
 Procedure tests for *DISTRIBUTION and *ELEMENT PROPERTIES element SAX1
*NODE,NSET=NALL
 1, 0., 0.
 2, 1., 0.
 3, 2., 0.
 4, 3., 0.
 5, 4., 0.
*Elset, elset=shellSec1
 1,
**---------------------------------------------------------------------
*Element, type=sax1,elset=eall
 1, 1, 2
 2, 2, 3
 3, 3, 4
 4, 4, 5
**---------------------------------------------------------------------
*SHELL SECTION,MATERIAL=MAT1, ELSET=eall, offset=distOffset, shell thickness=distThick
 .1,
**---------------------------------------------------------------------
*Distribution Table, name=thicknessTable
 length
*Distribution, name=distThick, location=element, table=thicknessTable                                        
         ,   0.1,
        2,   0.2,
        3,   0.3,
        4,   0.4,
*Distribution Table, name=offsetTable
 ratio
*Distribution, name=distOffset, location=element, table=offsetTable
         ,  -0.05,
        2,  -0.05,
        3,  -0.025,
        4,  -0.025,
**---------------------------------------------------------------------
*NSET,NSET=EDGE
 5,
*NSET,NSET=CENTER
 1,
*NSET,NSET=MID
 3,
*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-03,
*Damping, Alpha=0.5
**---------------------------------------------------------------------
*Boundary
CENTER,1,1
EDGE,1,2
*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.e7, 300., 0.
 1.e7, 1000., 0.
**
** --------------------------------------------------
**
*step, name=step-1
 frequency analysis
*frequency
 5,
*El print
 s,                                                                                                 
 sf,
 ssavg,
 se,
*Node print
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
 sf,
 ssavg,
 sth,
*end step
**
** --------------------------------------------------
**
*step, name=step-2
 SSD mode-based
*steady state dynamics, interval=eigenfrequency
 60,
*El print, freq=999
 s,                                                                                                 
 sf,
 ssavg,
 se,
 sth,      
*Node print, freq=999
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
 sf,
 ssavg,
 sth,
*cload
 center,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
 center,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
 center,2,0.5
*end step
**
** --------------------------------------------------
**
*step, name=step-5
 SSD subspace
*steady state dynamics, subspace projection, interval=eigenfrequency
 60.,
*cload
 center,2,0.5
*end step
**
** --------------------------------------------------
**
*step, name=step-6
 random response
*random response
 60.,
*select eigenmodes, generate
 1, 5
*cload, load case=1
 center,2,0.5
*cload, load case=2
 edge,6,0.2
*correlation, psd=white, type=uncorrelated, complex=no
 1,1
 2,1
*El print, freq=999                                                                                 
 sf,
 se,
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 sf,
*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
 edge,6,6,0.1
*El print, freq=999                                                                                 
 s,                                                                                                 
 sf,
 ssavg,
 se,
 sth,      
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
 sf,
 ssavg,
 sth,
*Output, history                                                                                     
*Element output,elset=shellSec1                                                                                   
 s11,                                                                                                 
 sf1,
 ssavg1,
 se1,
 sk1,
 sth,
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-9
 static 
*static
*El print, freq=999                                                                                 
 s,                                                                                                 
 sf,
 ssavg,
 se,
 sth,      
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
 sf,
 ssavg,
 sth,
*load case, name=load1
*dload
 eall, br, 1.e                                                                        
*end load case
*load case, name=load2
*dload
 eall, grav, 1.e6, 0.,1.,0.
*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