*HEADING                                                          
 : TABULAR FREQUENCY DOMAIN VISCOELASTICITY, ELASTIC, C3D8
 STEADY STATE DYNAMIC ANALYSIS OF CANTILEVER WITH DISTRIBUTED VERTICAL LOAD
 UNITS: SI
*RESTART,WRITE
*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
*NCOPY,CHANGE NUMBER=100,OLD SET=SIDE,SHIFT
,,-1.

*NSET,NSET=END,GENERATE
1,16,5
101,116,5
*ELEMENT,TYPE=C3D8
1,1,6,7,2,101,106,107,102
*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.775E-3,-5.116E-3,,,100.0
*USER MATERIAL,CONSTANTS=4
200.E9,0.3,1.775E-3,-5.116E-3
*BOUNDARY
END,1,2
1,3,3
16,3,3
**
*******************************************
** reference solution
** STEP1 
*STEP,INC=500
*STEADY STATE DYNAMIC,DIRECT
100.
*DLOAD
TOP,P4,1.E9
*NODE PRINT,NSET=B
U, 
RF, 
PU, 
PRF, 
*EL PRINT
S, 
E, 
*ENERGY PRINT
*OUTPUT,FIELD
*NODE OUTPUT,NSET=B
U, 
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=B
U, 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=END
RF, 
*OUTPUT,HISTORY
*NODE OUTPUT, NSET=END
RF, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S, 
E, 
*OUTPUT,HISTORY
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP
*****************************
**STEP 2 - extract all modes
*STEP
*FREQUENCY, eigensolver=lanczos, SIM=NO
102, , 
*NODE PRINT
*EL PRINT
*OUTPUT,FIELD
*NODE OUTPUT
*OUTPUT,FIELD
*ELEMENT OUTPUT
*END STEP
*****************************
** the subspace solution uses the entire eigenspace
** solution should be identical to the SSD,direct solution in the 
** previous step
** STEP 3
*STEP
*STEADY STATE DYNAMIC,SUBSPACE PROJECTION=ALL,INTERVAL=RANGE
100.
*DLOAD
TOP,P4,1.E9
*NODE PRINT,NSET=B
U, 
RF, 
PU, 
PRF, 
*EL PRINT
S, 
E, 
*ENERGY PRINT
*OUTPUT,FIELD
*NODE OUTPUT,NSET=B
U, 
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=B
U, 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=END
RF, 
*OUTPUT,HISTORY
*NODE OUTPUT, NSET=END
RF, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S, 
E, 
*OUTPUT,HISTORY
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP