*HEADING
**
** ABAQUS INPUT FILE UPGRADE: This input file was upgraded from
** version 5.8 or later to version 6.3 with 0 warning(s) and 0 
** error(s).
** 
**
: RANDOM RESPONSE PROCEDURE FOR BASE MOTION
**
**  THIS INPUT FILE TESTS
**    1) CORRELATED AND UNCORRELATED LOADINGS FOR 
**       BASE MOTION
**    2) LOAD CASE NUMBERING
**  
**  PROBLEM: CANTILEVER BEAM 
**    1) TWO B21 ELEMENTS OF UNIT LENGTH ALIGNED VERTICALLY
**    2) BASE MOTION APPLIED TO GLOBAL DOFS 1 AND 6
**    3) WHITE NOISE POWER SPECTRAL DENSITY INPUT
**    4) TWO MODE SHAPES USED IN RANDOM RESPONSE
**
**  RESULTS: AGREE WITH HAND CALCULATIONS
**    1) RESULTS FROM STEP 2 AND STEP 4 SHOULD BE THE SAME
**       (TESTS LOAD CASE NUMBERING)
**    2) RESULTS FROM STEP 3 AND STEP 5 SHOULD BE THE SAME
**       (TESTS LOAD CASE NUMBERING)
**
*NODE,NSET=ALL
10,0., 0.,0.
20,0.,10.,0.
30,0.,20.,0.
*transform, nset=all
1., 1., 0., -1., 1., 0.
*ELEMENT,ELSET=ALL,TYPE=B21
10,10,20
20,20,30
*BEAM SECTION,SECTION=RECT,ELSET=ALL,MATERIAL=STEEL
1.,1.
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,0.3
*DENSITY
0.283,
*BOUNDARY
10,ENCASTRE
*PSD-DEFINITION,NAME=WHITE,TYPE=BASE
1.0E+2,0.,0.
1.0E+2,0.,3000.
*STEP
***
   STEP 1: FREQUENCY STEP
***
*FREQUENCY, eigensolver=lanczos
2, , 
*NODE PRINT,FREQ=1
U, 
*EL PRINT,FREQ=0
*NODE FILE,FREQ=0
*EL FILE,FREQ=0
*MODAL FILE
*ENDSTEP
*STEP
***
   STEP 2: CORRELATED CLOADS
***
*RANDOM RESPONSE
1,25,10,5,1
*modal damping, modal=direct
1,2,0.01
*select eigenmodes, generate
1, 2
*CORRELATION,PSD=WHITE,COMPLEX=YES,TYPE=CORRELATED
1,1.,1.
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=1,DOF=1
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=1,DOF=6
**
*EL PRINT,FREQ=10
S,
SF,
*EL FILE,ELSET=ALL
S,
SF,
*NODE PRINT,NSET=ALL,FREQ=10
U,
TU,
RTU,
*NODE PRINT,NSET=ALL,FREQ=10,GLOBAL=YES
U,
TU,
RTU,
*NODE FILE,NSET=ALL
U,
TU,
RTU
*OUTPUT,FIELD,VAR=PRESELECT
*ENDSTEP
*STEP
***
   STEP 3: UNCORRELATED
***
*RANDOM RESPONSE
1,25,10,5,1
*modal damping, modal=direct
1,2,0.01
*select eigenmodes, generate
1, 2
*CORRELATION,PSD=WHITE,COMPLEX=YES,TYPE=UNCORRELATED
2,1.,1.
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=2,DOF=1
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=2,DOF=6
**
*EL PRINT,FREQ=10
S,
SF,
*EL FILE,ELSET=ALL
S,
SF,
*NODE PRINT,NSET=ALL,FREQ=10
U,
TU,
RTU,
*NODE PRINT,NSET=ALL,FREQ=10,GLOBAL=YES
U,
TU,
RTU,
*NODE FILE,NSET=ALL
U,
TU,
RTU
*ENDSTEP
*STEP
***
   STEP 4: CORRELATED WITH DIFFERENT LOAD CASE NUMBER
***
*RANDOM RESPONSE
1,25,10,5,1
*modal damping, modal=direct
1,2,0.01
*select eigenmodes, generate
1, 2
*CORRELATION,PSD=WHITE,COMPLEX=YES,TYPE=CORRELATED
1000,1.,1.
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=1000,DOF=1
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=1000,DOF=6
**
*EL PRINT,FREQ=10
S,
SF,
*EL FILE,ELSET=ALL
S,
SF,
*NODE PRINT,NSET=ALL,FREQ=10
U,
TU,
RTU,
*NODE PRINT,NSET=ALL,FREQ=10,GLOBAL=YES
U,
TU,
RTU,
*NODE FILE,NSET=ALL
U,
TU,
RTU
*ENDSTEP
*STEP
***
   STEP 5: UNCORRELATED WITH DIFFERENT LOAD CASE NUMBER
***
*RANDOM RESPONSE
1,25,10,5,1
*modal damping, modal=direct
1,2,0.01
*select eigenmodes, generate
1, 2
*CORRELATION,PSD=WHITE,COMPLEX=YES,TYPE=UNCORRELATED
2000,1.,1.
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=2000,DOF=1
*BASE MOTION,TYPE=ACCELERATION,LOAD CASE=2000,DOF=6
**
*EL PRINT,FREQ=10
S,
SF,
*EL FILE,ELSET=ALL
S,
SF,
*NODE PRINT,NSET=ALL,FREQ=10
U,
TU,
RTU,
*NODE PRINT,NSET=ALL,FREQ=10,GLOBAL=YES
U,
TU,
RTU,
*NODE FILE,NSET=ALL
U,
TU,
RTU
*ENDSTEP