*HEADING
STEADY STATE VIBRATION OF 2-DOF SYSTEM (VISCOELASTIC DAMPING)
*NODE,NSET=FILE
 1,10.
 2,15.
 3,20.
*ELEMENT,TYPE=DASHPOT2,ELSET=D
1,2,1
11,3,2
*DASHPOT,ELSET=D,DEPENDENCIES=1
1,1
*INCLUDE,INPUT=vibration_dampdata2.inp
**File contains frequency-dependent dashpot coefficients
**c(omega) = spring_stiffness*(2.3508E-2)*(omega/2*pi)**(-1.38366)
*ELEMENT,TYPE=SPRING2,ELSET=S
2,2,1
12,3,2
*SPRING,ELSET=S
1,1
*INCLUDE,INPUT=vibration_springdata.inp
**File contains frequency-dependent spring stiffness
**k(omega)=spring_stiffness*
**                 [(1-omega*(6.5001E-2)*(omega/2*pi)**(-1.38366)]
*ELEMENT,TYPE=MASS,ELSET=M
3,2
13,3
*MASS  ,ELSET=M
0.02588, 
*NSET,NSET=NALL
1,2,3
*BOUNDARY
 1,1,3
NALL,2,3
*STEP
FORCED OSCILLATION WITH VISCOELASTIC DAMPING 
*STEADY STATE DYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
0.77, 14., 378, 1
*CLOAD
2,1,1.0
*NODE PRINT,FREQ=40
U,PU
V,
A,
PRF,
**NODE FILE,F=1
**U,PU
*NODE FILE,F=40
V,
A,
PRF,
U,
PU,
*EL PRINT,FREQ=40
S,PHS
E,PHE
*EL FILE,FREQ=40
S,PHS
E,PHE
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=NALL
U,
*OUTPUT,HISTORY,FREQUENCY=40
*NODE OUTPUT, NSET=NALL
V,
A,
*ELEMENT OUTPUT,ELSET=D
S,
E,
*END STEP