*Heading
*Node
      1,          0.5,          0.5,           0.
      2,           0.,          0.5,          0.5
      3,          0.5,          0.5,          0.5
      4,           0.,         -0.5,          0.5
      5,          0.5,         -0.5,           0.
      6,          0.5,         -0.5,          0.5
      7,         -0.5,         -0.5,           1.
      8,         -0.5,         -0.5,          0.5
      9,         -0.5,          0.5,           1.
     10,         -0.5,          0.5,          0.5
     11,          0.5,          0.5,           1.
     12,           0.,          0.5,           1.
     13,          0.5,         -0.5,           1.
     14,           0.,         -0.5,           1.
     15,         -0.5,          0.5,           0.
     16,           0.,          0.5,           0.
     17,         -0.5,         -0.5,           0.
     18,           0.,         -0.5,           0.
*Element, type=C3D6, elset=cube
 1,  3,  2,  1,  6,  4,  5
 2,  8,  7,  4, 10,  9,  2
 3, 12,  2, 11, 14,  4, 13
 4, 16,  2, 15, 18,  4, 17
 5, 14,  4,  7, 12,  2,  9
 6, 10, 15,  2,  8, 17,  4
 7, 18,  4,  5, 16,  2,  1
 8,  6,  4, 13,  3,  2, 11
*Distribution Table, name=OrientationTable
 coord3d, coord3d
*Distribution, name=distOri, table=orientationTable, location=element
  , 1., 0., 0.,  0., 1., 0.
 1,  1., 1., 1., -1., 1.,-1.
 2,  1., 1.,-1., -1., 1., 1.
 3, -1., 1.,-1.,  1., 1., 1.
 4, -1.,-1., 1.,  1.,-1.,-1.
 5, -1., 1., 1.,  1., 1.,-1.
 6, -1.,-1.,-1.,  1.,-1., 1.
 7,  1.,-1.,-1., -1.,-1., 1.
 8,  1.,-1., 1., -1.,-1.,-1.
*Distribution Table, name=AngleTable
 angle
*Distribution, name=distAngle, table=AngleTable, location=element
  ,0
 1, 0.0
 2, 1.0
 3, 2.0
 4, 3.0
 5, 4.0
 6, -1.0
 7, -2.0
 8, -3.0
*Orientation, name=ori
 distOri
 3, distAngle
*Solid Section, elset=cube, material=Material-1, orientation=ori
 1.,
*Nset, nset=Fix
  4,  5,  6,  7,  8, 11, 13, 14, 17, 18
*Nset, nset=Load
  1,  2,  3,  9, 10, 12, 15, 16
*Material, name=Material-1
*Elastic, type=ENGINEERING CONSTANTS
  3e+06,  4e+06,  5e+06,    0.1,    0.2,    0.3,300000.,400000.
 500000.,
*Density
 0.750E-3
*Elset, elset=eall, generate
 1, 8, 1
*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.  
*Boundary
 Fix, 1,3          
**
** --------------------------------------------------
**
*step, name=step-1
 frequency analysis
*frequency
 5,
*El print                                                                               
 s,                                                                                                 
*Node print                                                                           
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
*end step
**
** --------------------------------------------------
**
*step, name=step-2
 SSD mode-based
*steady state dynamics, interval=eigenfrequency
 60,
*El print, freq=999                                                                                 
 s,                 
*Node print, freq=0                                                                                
*Output, field, op=new                                                                                 
*Element output                                                                                     
 s,           
*cload
 load, 2, 1.e5
 load, 3, 1.e5
*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
 load, 2, 1.e5
 load, 3, 1.e5
*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
 load, 2, 100.
*end step
**
** --------------------------------------------------
**
*step, name=step-5
 SSD subspace
*steady state dynamics, subspace projection, interval=eigenfrequency
 60.,
*cload
 load, 2, 100.
*end step
**
** --------------------------------------------------
**
*step, name=step-6
 random response
*random response
 60.,
*select eigenmodes, generate
 1, 5
*cload, load case=1
 load, 2, 100.
*cload, load case=2
 load, 2, 100.
*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.,1.,0.,1.
*select eigenmodes, generate
 1, 5
*end step
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-8
 static
*static
 1.0, 1.0
*Boundary
 load, 1, 1, 0.025
 load, 2, 2, 0.050
*El print, freq=999                                                                                 
 s,                                                                                                 
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field, op=new                                                                                            
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-9
 static 
*static
*El print, freq=999                                                                                 
 s,      
*Node print, freq=0                                                                          
*Output, field, op=new                                                                                            
*Element output                                                                                     
 s,                                                                                                 
*load case, name=load1
*dload
 eall, bx, 1.e3                                                                             
*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