*Heading
*Node
*Elset, elset=cube1
*Elset, elset=cube2
*Elset, elset=cube3
*Elset, elset=cube4
*Elset, elset=cube5
*Elset, elset=cube6
*Element, type=C3D4, elset=cube
*Solid Section, elset=cube, material=Material-1
 1.,
*Distribution, name=distOri, type=orientation, location=element
 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.
*Element Properties
 orientation, distOri
*Nset, nset=Fix
*Nset, nset=Load
*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
*psd-definition, name=white, type=force
*spectrum, name=spec,type=acceleration
*Boundary
 Fix, 1,3
**
** --------------------------------------------------
**
*step, perturbation, name=step-1
 perturbation step for residual modes
*static
*Output, field
*Node output, Nset=Load
 U,
*load case, name=step2
*cload
 load, 2, 1.e5
*end load case
*load case, name=step4
*cload
 load, 3, 1.e5
*end load case
*load case, name=step5to8
*cload
 load, 2, 100.
*end load case
*load case, name=step10bx
*dload
 eall, bx, 1.e3
*end load case
*load case, name=step10grav
*dload
 eall, grav, 1.e6, 0.,1.,0.
*end load case
*load case, name=step10cent
*dload
 eall, cent, 100.,0.,0.,0.,1.,0.
*end load case
*load case, name=step10centrif
*dload
 eall, centrif, 1.e6,0.,0.,0.,1.,0.
*end load case
*end step
**
** --------------------------------------------------
**
*step, name=step-2
 frequency analysis with residual modes
*frequency, residual modes
 5,
*El print
 s,
*Node print
 u,
*Output, field
*Node output
 u,
*Element output
 s,
*end step
**
** --------------------------------------------------
**
*step, name=step-3
 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
*modal damping, rayleigh
*end step
**
** --------------------------------------------------
**
*step, name=step-4
 modal dynamic
*modal dynamic
 0.1,1.
*Output, field
*Node output, Nset=Load
 U,
*cload
 load, 2, 1.e5
 load, 3, 1.e5
*select eigenmodes, generate
*modal damping, rayleigh
*end step
**
** --------------------------------------------------
**
*step, name=step-5
 SSD direct
*steady state dynamics, direct, interval=eigenfrequency
 60.,
*Output, field
*Node output, Nset=Load
 U,
*cload
 load, 2, 100.
*end step
**
** --------------------------------------------------
**
*step, name=step-6
 SSD subspace
*steady state dynamics, subspace projection, interval=eigenfrequency
 60.,
*Output, field
*Node output, Nset=Load
 U,
*cload
 load, 2, 100.
*end step
**
** --------------------------------------------------
**
*step, name=step-7
 random response
*random response
*select eigenmodes, generate
*Output, field
*Node output, Nset=Load
 U,
*cload, load case=1
 load, 2, 100.
*cload, load case=2
 load, 2, 100.
*correlation, psd=white, type=uncorrelated, complex=no
*end step
**
** --------------------------------------------------
**
*step, name=step-8
 response spectrum
*response spectrum, sum=abs, comp=alg
 spec,0.,1.,0.,1.
*select eigenmodes, generate
*Output, field
*Node output, Nset=Load
 U,
*end step
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-9
 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-10
 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