*HEADING
  - HYPERELASTICITY (POLYNOMIAL), FRQ DOMAIN VISCOELASTICITY
**PREPRINT,MODEL=YES
*NODE,NSET=ALL
11,10.,0.,0.
12,11.,0.,0.
13,11.,1.,0.
14,10.,1.,0.
15,10.,0.,1.
16,11.,0.,1.
17,11.,1.,1.
18,10.,1.,1.
*NSET,NSET=FACE1
11,12,13,14
*NSET,NSET=FACE2
15,16,17,18
*NSET,NSET=FACE3
11,12,15,16
*NSET,NSET=FACE4
12,13,16,17
*NODE,NSET=OUTNODE
1000
*EQUATION
2
FACE4,1,1.0,1000,1,-1.0
*NSET,NSET=FACE5
13,14,17,18
*NSET,NSET=FACE6
14,11,18,15
*ELEMENT,TYPE=C3D8RH,ELSET=TEST
10,11,12,13,14,15,16,17,18
*SOLID SECTION,ELSET=TEST,MATERIAL=POLY2
*MATERIAL,NAME=POLY2
*DENSITY
8.E-7, 
*HYPERELASTIC,N=1,MODULI=LONG TERM
80.,20.
*VISCOELASTIC,FREQUENCY=TABULAR,PRELOAD=UNIAXIAL
6.,630.,1.,0.
3.,615.,100.,0.
5.,500.,1.,0.24
3.,300.,100.,0.24
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUTNODE,1,1,0.12
*EL PRINT,F=10
S, 
NE, 
*NODE PRINT,F=10
U,RF
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
 50.5
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUTNODE,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUTNODE
RF1,U1
**---------------------------
**  total RF1 on FACE4
**		real part    imaginary part	
**theoretical:   511.25		4.25
**     ABAQUS:   511.2		4.25
**---------------------------
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP