*HEADING
 MULTISTEP VERIFICATION TEST FOR CPE4R ELEMENTS USING DISTRIBUTIONS
*NODE
*NSET, NSET=SIDE1, GENERATE
*NSET, NSET=SIDE3, GENERATE
*ELSET, ELSET=EALL, GENERATE
**------------------------------------------------------------------------------------------
*MATERIAL, NAME=MAT1
*ELASTIC, TYPE=ENGINEERING CONSTANTS
  0.300E+07,  0.400E+07,  0.500E+07,  0.100E+00,  0.200E+00,  0.300E+00,  0.300E+06,  0.400E+06
  0.500E+06,
*DENSITY
  0.750E-03,
**------------------------------------------------------------------------------------------
*ELEMENT, TYPE=CPE4R, ELSET=EALL
*DISTRIBUTION TABLE, NAME=ORIENTATIONTABLE
 COORD3D, COORD3D
*DISTRIBUTION, NAME=DISTORIENT, LOCATION=ELEMENT, TABLE=ORIENTATIONTABLE
*DISTRIBUTION TABLE, NAME=ANGLETABLE
 ANGLE
*DISTRIBUTION, NAME=DISTANGLE, TABLE=ANGLETABLE, LOCATION=ELEMENT
*SOLID SECTION, MATERIAL=MAT1, ELSET=EALL, ORIENTATION=ORI
1.,
*ORIENTATION, NAME=ORI
 DISTORIENT
 3, DISTANGLE
*BOUNDARY
 SIDE1, 1, 2
*PSD-DEFINITION, NAME=WHITE, TYPE=FORCE
*SPECTRUM, NAME=SPEC,TYPE=ACCELERATION
**
** --------------------------------------------------
**
*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
 SIDE3, 2, 100.
*SELECT EIGENMODES, GENERATE
*MODAL DAMPING, RAYLEIGH
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-3
 MODAL DYNAMIC
*MODAL DYNAMIC
 0.1,1.
*CLOAD
 SIDE3, 2, 100.
*SELECT EIGENMODES, GENERATE
*MODAL DAMPING, RAYLEIGH
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-4
 SSD DIRECT
*STEADY STATE DYNAMICS, DIRECT, INTERVAL=EIGENFREQUENCY
 60.,
*CLOAD
 SIDE3, 2, 100.
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-5
 SSD SUBSPACE
*STEADY STATE DYNAMICS, SUBSPACE PROJECTION, INTERVAL=EIGENFREQUENCY
 60.,
*CLOAD
 SIDE3, 2, 100.
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-6
 RANDOM RESPONSE
*RANDOM RESPONSE
*SELECT EIGENMODES, GENERATE
*CLOAD, LOAD CASE=1
 SIDE3, 2, 100.
*CLOAD, LOAD CASE=2
 SIDE3, 2, 100.
*CORRELATION, PSD=WHITE, TYPE=UNCORRELATED, COMPLEX=NO
*END STEP
**
** --------------------------------------------------
**
*STEP, NAME=STEP-7
 RESPONSE SPECTRUM
*RESPONSE SPECTRUM, SUM=ABS, COMP=ALG
 SPEC,0.,1.,0.,1.
*SELECT EIGENMODES, GENERATE
*END STEP
**
** --------------------------------------------------
**
*STEP, NLGEOM=YES, NAME=STEP-8
 STATIC
*STATIC
 1.0, 1.0
*BOUNDARY
 SIDE3, 1, 1, 0.025
 SIDE3, 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