*HEADING                                                          
 TABULAR FREQUENCY DOMAIN VISCOELASTICITY, ELASTIC, CPS4
 STEADY STATE DYNAMIC ANALYSIS OF CANTILEVER WITH DISTRIBUTED VERTICAL LOAD
 UNITS: SI
*NODE,NSET=A
1, 
6,,1.
11,,2.
16,,3.
*NODE,NSET=B
5,10.
10,10.,1.
15,10.,2.
20,10.,3.
*NFILL,NSET=SIDE
A,B,4
*NSET,NSET=END,GENERATE
1,16,5
101,116,5
*ELEMENT,TYPE=CPS4
1,1,2,7,6
*ELGEN,ELSET=BEAM
1,4,,,3,5,4
*ELSET,ELSET=TOP,GENERATE
9,12
*SOLID SECTION,ELSET=BEAM,MATERIAL=VISCO
*MATERIAL,NAME=VISCO
*DENSITY
8000., 
*ELASTIC,MODULI=LONG TERM
200.E9,0.3
*VISCOELASTIC,FREQUENCY=TABULAR
1.161E-2,-3.21E-2,,,1.
7.849E-3,-2.222E-2,,,15.8
5.354E-3,-1.533E-2,,,25.1
3.639E-3,-1.062E-2,,,39.8
2.543E-3,-7.382E-3,,,63.1
1.775E-3,-5.116E-3,,,100.
*BOUNDARY
END,1,2
**
** STEP 1 - Direct SSD
** 
*STEP
 SSD, DIRECT
**
*STEADY STATE DYNAMIC, DIRECT, INTERVAL=RANGE, FREQUENCY SCALE=LINEAR
1.,100.,50
*DLOAD
TOP,P3,1.E9
*NODE PRINT,NSET=B, FREQ=5
U, 
PU,  
*OUTPUT,FIELD
*NODE OUTPUT,NSET=B
U, V, A,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=B
U, V, A,
*END STEP
***
** STEP 2 - EXTRACT ALL EIGENMODES AT 50 Hz
***
*STEP
 FREQUENCY EXTRACTION at 50 Hz USING Lanczos
**
*FREQUENCY, eigensolver=LANCZOS, NORM=MASS, PROPERTY EVALUATION=50
16, ,1000.,
*OUTPUT,FIELD
*NODE OUTPUT
U,
*END STEP
***
** STEP 3 - SSD USING EIGENMODES BY Lanczos
** 
*STEP
 SSD with stiffness properties evaluated at 50 Hz
**
*STEADY STATE DYNAMIC, SUBSPACE PROJECTION=CONSTANT, INTERVAL=RANGE, FREQUENCY SCALE=LINEAR
1.,100.,50
*DLOAD
TOP,P3,1.E9
*NODE PRINT,NSET=B, FREQ=5
U,  
PU,  
*OUTPUT,FIELD
*NODE OUTPUT,NSET=B
U, V, A,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=B
U, V, A,
*END STEP