*HEADING
 : TABULAR FREQUENCY DOMAIN VISCOELASTICITY, ELASTIC, C3D8
 STEADY STATE DYNAMIC ANALYSIS OF CANTILEVER WITH DISTRIBUTED VERTICAL LOAD
 UNITS: SI
*RESTART,WRITE
*NODE,NSET=A
*NODE,NSET=B
*NFILL,NSET=SIDE
A,B,4
*NCOPY,CHANGE NUMBER=100,OLD SET=SIDE,SHIFT
*NSET,NSET=END,GENERATE
*ELEMENT,TYPE=C3D8
*ELGEN,ELSET=BEAM
*ELSET,ELSET=TOP,GENERATE
*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.0
*BOUNDARY
END,1,2
1,3,3
16,3,3
**
*******************************************
** reference solution
** STEP1
*STEP,INC=500
*STEADY STATE DYNAMIC,DIRECT
*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
*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
*****************************
**check for INTERVAL=EIGENFREQUENCY
**STEP 4
*STEP,INC=500
*STEADY STATE DYNAMIC,DIRECT,INTERVAL=EIGENFREQUENCY
*DLOAD
TOP,P4,1.E9
*NODE PRINT,NSET=B
U,
RF,
PU,
PRF,
*EL PRINT
S,
E,
*ENERGY PRINT
*END STEP
*****************************
**check for INTERVAL=EIGENFREQUENCY
**STEP 5
*STEP
*STEADY STATE DYNAMIC,SUBSPACE PROJECTION=EIGENFREQUENCY,INTERVAL=EIGENFREQUENCY
*DLOAD
TOP,P4,1.E9
*NODE PRINT,NSET=B
U,
RF,
PU,
PRF,
*EL PRINT
S,
E,
*ENERGY PRINT
*END STEP
*****************************
** extract fewer modes with stiffness propoerties evalutated at 50 Hz
** STEP 6
*STEP
*FREQUENCY,PROPERTY EVALUATION=50, eigensolver=lanczos, SIM=NO
16, ,
*NODE PRINT
*EL PRINT
*OUTPUT,FIELD
*NODE OUTPUT
*OUTPUT,FIELD
*ELEMENT OUTPUT
*END STEP
*****************************
** subspace solution is not exact but only very close
**STEP 7
*STEP
*STEADY STATE DYNAMIC,SUBSPACE PROJECTION=EIGENFREQUENCY,INTERVAL=RANGE
*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
*****************************
**check for SUBSPACE PROJECTION=PROPERTY CHANGE
**STEP 8
*STEP
*STEADY STATE DYNAMIC,SUBSPACE PROJECTION=PROPERTY CHANGE,INTERVAL=RANGE,
 STIFFNESS CHANGE=0.2, DAMPING CHANGE=0.1
*DLOAD
TOP,P4,1.E9
*NODE PRINT,NSET=B
U,
RF,
PU,
PRF,
*EL PRINT
S,
E,
*ENERGY PRINT
*END STEP