*HEADING
 Procedure tests for *DISTRIBUTION and *ELEMENT PROPERTIES element SAXA22
*NODE,NSET=NALL
 1, 0., 0.
 2, 0.5, 0.
 3, 1., 0.
 4, 1.5, 0.
 5, 2., 0.
 6, 2.5, 0.
 7, 3., 0.
 8, 3.5, 0.
 9, 4., 0.
*Elset, elset=shellSec1
 1,
**---------------------------------------------------------------------
*Element,TYPE=SAXA22,elset=eall, offset=100
 1, 1, 2, 3
 2, 3, 4, 5
 3, 5, 6, 7
 4, 7, 8, 9
**---------------------------------------------------------------------
*SHELL SECTION,MATERIAL=MAT1, ELSET=eall, orientation=ori, offset=distOffset, shell thickness=distThick                    
 .1,
**---------------------------------------------------------------------
*Distribution Table, name=orientationTable
 coord3d, coord3d
*Distribution, name=distOrient, location=element, table=orientationTable
   , 0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00
  1, 0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00
  2, 0.9993908E+00,   0.3489950E-01,   0.0000000E+00,  -0.3489950E-01,   0.9993908E+00,   0.0000000E+00
  3, 0.9961947E+00,   0.8715574E-01,   0.0000000E+00,  -0.8715574E-01,   0.9961947E+00,   0.0000000E+00
  4, 0.9902681E+00,   0.1391731E+00,   0.0000000E+00,  -0.1391731E+00,   0.9902681E+00,   0.0000000E+00
*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,
*Distribution Table, name=OriAngleTable
 angle
*Distribution, name=distOriAngle, location=element, table=OriAngleTable
         ,   0.0,                                     
        1,  -0.5
        2,  -0.5
        3,  -1.0
        4,  -1.0
*orientation, name=ori
 distOrient
 3, distOriAngle
**---------------------------------------------------------------------
*NSET,NSET=EDGE
 9, 109, 209
*NSET,NSET=CENTER
 1, 101, 201
*NSET,NSET=MID
 5, 105, 205
*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,
*Conductivity
  1000.,
*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, SIM=NO
 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
*Boundary                                                                                           
 CENTER,2,,0.1
*end load case
*load case, name=load2
*Boundary
 EDGE,6,,0.2
*end load case
*load case, name=load3
*Boundary                                                                                           
 MID,6,,-0.1
*end load case
*end step