*HEADING
Steady state dynamic analysis of the cantilever beam with dashpot - using SIM based SSD,Subspace
and multiple load cases -  Units: mm, N, MPa
*RESTART, WRITE
*NODE
 1,   0.000,  0.000                  
 3,  40.000,  0.000                  
 5,  80.000,  0.000                  
 7,  120.00,  0.000                  
 9,  160.00,  0.000                  
11,  200.00,  0.000                  
12,  200.00,  -10.000                  
*ELEMENT, TYPE=B21, ELSET=BEAMS
 1,   1,  3
 2,   3,  5
 3,   5,  7
 4,   7,  9
 5,   9, 11
*ELEMENT, TYPE=DASHPOTA, ELSET=DASH
 6,  11, 12
**
** Section and Material data for beams with damping
**
*BEAM SECTION, SECTION=RECT, ELSET=BEAMS, MATERIAL=MATEA
50., 5.
*MATERIAL, NAME=MATEA
*ELASTIC
2.E5,.3
*DENSITY
2.3E-6
*DAMPING, ALPHA=0.05, BETA=0.0
**
** Data for Dashpot element (Frequency dependent)
**
*DASHPOT, ELSET=DASH,DEPENDENCIES=1

0.01,0
0.001,100
0.0005,200
*nset,nset=rf
1,12
*BOUNDARY
1,1,2
1,6,6
12,1,2
*NSET, NSET=TIP
11
**
** FREQUENCY step to extract eigenvalues with AMS
**
*STEP
*FREQUENCY, EIGENSOLVER=AMS, NORMALIZATION=MASS, SIM
0,0,200
*END STEP
**
** STEADY STATE DYNAMICS with multiple loadcases
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE=ALL
1,100,1
*NODE PRINT, NSET=rf,freq=100
RF,
*OUTPUT, HISTORY, FREQUENCY=1
*NODE OUTPUT, NSET=TIP
U,
*NODE OUTPUT, NSET=rf
RF,
*LOAD CASE, NAME="REAL LOAD"
*CLOAD, LOAD CASE=1
11, 2, -1200. 
*END LOADCASE
*LOAD CASE, NAME="IMAG LOAD"
*CLOAD, LOAD CASE=2
11, 2, -1200. 
*END LOADCASE
*END STEP
**
** STEADY STATE DYNAMICS STEP with a single load case for verification
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE=ALL
1,100,1
*OUTPUT, HISTORY, FREQUENCY=1
*NODE OUTPUT, NSET=TIP
U,
*CLOAD, LOAD CASE=1
11, 2, -1200. 
*END STEP