*HEADING
*NODE
      1,          0.5,           0.,           1.
      2,          0.5,           0.,          0.5
      3,          0.5,          0.5,          0.5
      4,          0.5,          0.5,           1.
      5,           0.,           0.,           1.
      6,           0.,           0.,          0.5
      7,           0.,          0.5,          0.5
      8,           0.,          0.5,           1.
      9,          0.5,           0.,           0.
     10,          0.5,          0.5,           0.
     11,           0.,           0.,           0.
     12,           0.,          0.5,           0.
     13,         -0.5,           0.,          0.5
     14,         -0.5,          0.5,          0.5
     15,         -0.5,           0.,           0.
     16,         -0.5,          0.5,           0.
     17,           0.,         -0.5,          0.5
     18,         -0.5,         -0.5,          0.5
     19,           0.,         -0.5,           1.
     20,         -0.5,         -0.5,           1.
     21,         -0.5,           0.,           1.
     22,         -0.5,          0.5,           1.
     23,           0.,         -0.5,           0.
     24,         -0.5,         -0.5,           0.
     25,          0.5,         -0.5,          0.5
     26,          0.5,         -0.5,           0.
     27,          0.5,         -0.5,           1.
     28,           0.,         0.25,           1.
     29,           0.,          0.5,         0.75
     30,           0.,         0.25,          0.5
     31,           0.,           0.,         0.75
     32,          0.5,           0.,         0.75
     33,          0.5,         0.25,          0.5
     34,          0.5,          0.5,         0.75
     35,          0.5,         0.25,           1.
     36,         0.25,           0.,          0.5
     37,         0.25,           0.,           1.
     38,         0.25,          0.5,          0.5
     39,         0.25,          0.5,           1.
     40,           0.,          0.5,         0.25
     41,           0.,         0.25,           0.
     42,           0.,           0.,         0.25
     43,          0.5,           0.,         0.25
     44,          0.5,         0.25,           0.
     45,          0.5,          0.5,         0.25
     46,         0.25,           0.,           0.
     47,         0.25,          0.5,           0.
     48,        -0.25,          0.5,          0.5
     49,         -0.5,         0.25,          0.5
     50,        -0.25,           0.,          0.5
     51,        -0.25,           0.,           0.
     52,         -0.5,         0.25,           0.
     53,        -0.25,          0.5,           0.
     54,         -0.5,           0.,         0.25
     55,         -0.5,          0.5,         0.25
     56,         -0.5,        -0.25,           1.
     57,        -0.25,         -0.5,           1.
     58,           0.,        -0.25,           1.
     59,        -0.25,           0.,           1.
     60,           0.,        -0.25,          0.5
     61,        -0.25,         -0.5,          0.5
     62,         -0.5,        -0.25,          0.5
     63,         -0.5,           0.,         0.75
     64,           0.,         -0.5,         0.75
     65,         -0.5,         -0.5,         0.75
     66,         -0.5,          0.5,         0.75
     67,        -0.25,          0.5,           1.
     68,         -0.5,         0.25,           1.
     69,           0.,         -0.5,         0.25
     70,         -0.5,         -0.5,         0.25
     71,        -0.25,         -0.5,           0.
     72,           0.,        -0.25,           0.
     73,         -0.5,        -0.25,           0.
     74,          0.5,        -0.25,          0.5
     75,          0.5,         -0.5,         0.25
     76,          0.5,        -0.25,           0.
     77,         0.25,         -0.5,          0.5
     78,         0.25,         -0.5,           0.
     79,          0.5,        -0.25,           1.
     80,          0.5,         -0.5,         0.75
     81,         0.25,         -0.5,           1.
*ELEMENT, TYPE=C3D20R, ELSET=CUBE
 1,  5,  6,  7,  8,  1,  2,  3,  4, 31, 30, 29, 28, 32, 33, 34,
    35, 37, 36, 38, 39
 2,  6, 11, 12,  7,  2,  9, 10,  3, 42, 41, 40, 30, 43, 44, 45,
    33, 36, 46, 47, 38
 3,  6, 13, 14,  7, 11, 15, 16, 12, 50, 49, 48, 30, 51, 52, 53,
    41, 42, 54, 55, 40
 4, 21,  5, 19, 20, 13,  6, 17, 18, 59, 58, 57, 56, 50, 60, 61,
    62, 63, 31, 64, 65
 5,  5,  6, 13, 21,  8,  7, 14, 22, 31, 50, 63, 59, 29, 48, 66,
    67, 28, 30, 49, 68
 6, 11,  6, 13, 15, 23, 17, 18, 24, 42, 50, 54, 51, 69, 61, 70,
    71, 72, 60, 62, 73
 7, 11,  6, 17, 23,  9,  2, 25, 26, 42, 60, 69, 72, 43, 74, 75,
    76, 46, 36, 77, 78
 8,  6,  5, 19, 17,  2,  1, 27, 25, 31, 58, 64, 60, 32, 79, 80,
    74, 36, 37, 81, 77
*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
  ,1.0
 1, 1.0
 2, 2.0
 3, 3.0
 4, 4.0
 5, 5.0
 6, 6.0
 7, 7.0
 8, 8.0
*ORIENTATION, NAME=ORI
 DISTORI
 3, DISTANGLE
*SOLID SECTION, ELSET=CUBE, MATERIAL=MATERIAL-1, ORIENTATION=ORI
 1.,
*NSET, NSET=FIX
  4, 17, 18, 19, 20, 23, 24, 25, 26, 27, 57, 61, 64, 65, 69, 70
 71, 75, 77, 78, 80, 81
*NSET, NSET=LOAD
  3,  7,  8, 10, 12, 14, 16, 22, 29, 34, 38, 39, 40, 45, 47
 48, 53, 55, 66, 67
*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,SIM,normalization=mass
 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