*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.775E-3,-5.116E-3,,,100.0
*USER MATERIAL,CONSTANTS=4
*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