*HEADING                                                          
 :  TABULAR FREQUENCY DOMAIN VISCOELASTICITY, HYPERELASTIC, 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., 
*HYPERELASTIC,N=1,MODULI=LONG TERM
33.333333E9,,12.E-12
*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
40.,100.,5
*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,FREQ=0
*NODE 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
40.,100.,5
*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
40.,100.,5
*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 5
*STEP 
*STEADY STATE DYNAMIC,SUBSPACE PROJECTION=EIGENFREQUENCY,INTERVAL=EIGENFREQUENCY
40.,100.,5
*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
*****************************
** 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
40.,100.,5
*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 
40.,100.,5
*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