*HEADING
*PREPRINT,MODEL=YES
*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.
*ELEMENT, TYPE=C3D8, ELSET=CUBE
 1, 5, 6, 7, 8, 1, 2, 3, 4
 2,  6, 11, 12,  7,  2,  9, 10,  3
 3,  6, 13, 14,  7, 11, 15, 16, 12
 4, 21,  5, 19, 20, 13,  6, 17, 18
 5,  5,  6, 13, 21,  8,  7, 14, 22
 6, 11,  6, 13, 15, 23, 17, 18, 24
 7, 11,  6, 17, 23,  9,  2, 25, 26
 8,  6,  5, 19, 17,  2,  1, 27, 25
*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
 1, DISTANGLE
*SOLID SECTION, ELSET=CUBE, MATERIAL=MATERIAL-1, ORIENTATION=ORI
 1.,
*NSET, NSET=FIX
 4, 17, 18, 19, 20, 23, 24, 25, 26, 27
*NSET, NSET=LOAD
  3,  7,  8, 10, 12, 14, 16, 22
*ELSET, ELSET=EALL, GENERATE
 1, 8, 1
** 
** MATERIALS
** 
*DISTRIBUTION TABLE,NAME=DISTTABLEORTHO
 MODULUS,MODULUS,MODULUS,MODULUS,MODULUS,MODULUS,MODULUS,
 MODULUS,MODULUS
*DISTRIBUTION, NAME=DISTORTHO, LOCATION=ELEMENT, TABLE=DISTTABLEORTHO
      ,4.03846153846154E6,1.73076923076923E6,4.03846153846154E6,1.73076923076923E6,1.73076923076923E6,4.03846153846154E6,1.15384615384615E6
      1.15384615384615E6,1.15384615384615E6
  CUBE,3.38095238095238E6,0.88888888888889E6,4.74074074074074E6,1.46031746031746E6,2.07407407407407E6,6.26455026455027E6,3.E5,
  4.E5,5.E5
*MATERIAL, NAME=MATERIAL-1
*ELASTIC, TYPE=ORTHOTROPIC
 DISTORTHO
*DENSITY                                                                                            
  0.750E-03,
*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