*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,ARRUDA-BOYCE,MODULI=LONG TERM
200., 5.
**, 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,ARRUDA-BOYCE,MODULI=LONG TERM
200., 5.
**, 0.00001
*VISCOELASTIC,FREQUENCY=TABULAR,PRELOAD=UNIAXIAL 
7.139762,	634.707,	1,	0	   
1.091566,	618.1127,	100,	0	   
6.00157,	433.3527,	1,	0.25	   
0.917553,	419.4038,	100,	0.25	   
5.341648,	306.3967,	1,	0.75	   
0.81666,	293.9816,	100,	0.75	   
5.550269,	283.8338,	1,	2	   
0.848555,	270.9338,	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