*HEADING
  - HYPERELASTICITY (POLYNOMIAL), FRQ DOMAIN VISCOELASTICITY
*PREPRINT,MODEL=YES
*NODE,NSET=ALL
1,
2,1.
3,1.,1.,
4,0.,1.,
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
**
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
1,2,3,4,11,12,13,14
*NSET,NSET=FACE2
5,6,7,8,15,16,17,18
*NSET,NSET=FACE3
1,2,5,6,11,12,15,16
*NSET,NSET=FACE4
2,3,6,7,12,13,16,17
*NSET,NSET=FACE4-1
2,3,6,7
*NSET,NSET=FACE4-2
12,13,16,17
*NODE, NSET=OUT1
1000
*NODE, NSET=OUT2
1001
*EQUATION
2
FACE4-1,1,1.0,OUT1,1,-1.0
2
FACE4-2,1,1.0,OUT2,1,-1.0
*NSET,NSET=FACE5
3,4,7,8,13,14,17,18
*NSET,NSET=FACE6
4,1,8,5,14,11,18,15
*ELEMENT,TYPE=C3D8RH,ELSET=REF
1,1,2,3,4,5,6,7,8
*ELEMENT,TYPE=C3D8RH,ELSET=TEST
10,11,12,13,14,15,16,17,18
*SOLID SECTION,ELSET=REF,MATERIAL=POLY1
*MATERIAL,NAME=POLY1
*DENSITY
8.E-7, 
*HYPERELASTIC,N=3,MODULI=LONG TERM
80.,20.,10.,15.,20.,12.,18.,16.,
20.
**,0.00001
*VISCOELASTIC,FREQUENCY=TABULAR
1.161E-2,-3.21E-2,,,1.
1.775E-3,-5.116E-3,,,100.0
*SOLID SECTION,ELSET=TEST,MATERIAL=POLY2
*MATERIAL,NAME=POLY2
*DENSITY
8.E-7, 
*HYPERELASTIC,N=3,MODULI=LONG TERM
80.,20.,10.,15.,20.,12.,18.,16.,
20.
**,0.00001
*VISCOELASTIC,FREQUENCY=TABULAR,PRELOAD=UNIAXIAL
6.966, 619.26, 1., 0.
1.065, 603.0696,100.,0.
8.264366,625.5973,1.,0.25
1.263501,606.3892,100.,0.25
41.48156,3178.671,1.,0.75
6.341927,3082.259,100.,0.75
505.441,38144.53,1.,2.
77.27458,36969.78,100.,2.
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
1.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
100.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,0.25
OUT2,1,1,0.25
*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
1.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
100.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,0.75
OUT2,1,1,0.75
*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
1.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*EL PRINT
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
100.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,2.
OUT2,1,1,2.
*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
1.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
100.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
OUT1,1,1,1.
OUT2,1,1,1.
*EL PRINT
S,
E,
*NODE PRINT, NSET=OUT1
 RF1,U1
*NODE PRINT, NSET=OUT2
 RF1,U1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,NE
*END STEP