*HEADING
*NODE
      1,          0.5,          0.5,           1.
      2,          0.5,         -0.5,           1.
      3,          0.5,         -0.5,           0.
      4,          0.5,          0.5,           0.
      5,         -0.5,         -0.5,           1.
      6,         -0.5,         -0.5,           0.
      7,         -0.5,          0.5,           0.
      8,         -0.5,          0.5,           1.
      9,           0.,           0.,          0.5
*ELSET, ELSET=CUBE1
 1,2
*ELSET, ELSET=CUBE2
 3,4
*ELSET, ELSET=CUBE3
 5,6
*ELSET, ELSET=CUBE4
 7,8
*ELSET, ELSET=CUBE5
 9,10
*ELSET, ELSET=CUBE6
 11,12
*ELEMENT, TYPE=C3D4, ELSET=CUBE
  1, 9, 6, 8, 7
  2, 9, 3, 6, 4
  3, 9, 4, 6, 7
  4, 9, 5, 8, 6
  5, 9, 3, 5, 6
  6, 9, 5, 2, 8
  7, 9, 2, 5, 3
  8, 9, 4, 2, 3
  9, 9, 1, 2, 4
 10, 9, 8, 2, 1
 11, 9, 7, 1, 4
 12, 9, 8, 1, 7
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
 COORD3D, COORD3D
*DISTRIBUTION, NAME=DISTORI, TABLE=ORIENTATIONTABLE, LOCATION=ELEMENT
      , 1., 0., 0.,  0., 1., 0.
 CUBE1,  1., 1., 1., -1., 1.,-1.
 CUBE2,  1., 1.,-1., -1., 1., 1.
 CUBE3, -1., 1.,-1.,  1., 1., 1.
 CUBE4, -1.,-1., 1.,  1.,-1.,-1.
 CUBE5, -1., 1., 1.,  1., 1.,-1.
 CUBE6, -1.,-1.,-1.,  1.,-1., 1.
*DISTRIBUTION TABLE, NAME=ANGLETABLE
 ANGLE
*DISTRIBUTION, NAME=DISTANGLE, TABLE=ANGLETABLE, LOCATION=ELEMENT
  ,0
 CUBE1, 0.0
 CUBE2, 1.0
 CUBE3, 2.0
 CUBE4, 3.0
 CUBE5, 4.0
 CUBE6, -1.0
*ORIENTATION, NAME=ORI
 DISTORI
 3, DISTANGLE
*SOLID SECTION, ELSET=CUBE, MATERIAL=MATERIAL-1, ORIENTATION=ORI
 1.,
*NSET, NSET=FIX
 1, 2, 3, 5, 6
*NSET, NSET=LOAD
 4, 7, 8
*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, 12, 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, PERTURBATION, NAME=STEP-0
 DEFINE RESIDUAL MODES
*STATIC
*LOAD CASE, NAME=LC1
*CLOAD
LOAD, 2, 1.
*END LOAD CASE
*LOAD CASE, NAME=LC2
*CLOAD
LOAD, 3, 1.
*END LOAD CASE
*OUTPUT,FIELD,FREQUENCY=0
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-1
 FREQUENCY ANALYSIS
*FREQUENCY,RESIDUAL MODES
 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