*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,MARLOW,MODULI=LONG TERM
**,POISSON=0.49
*UNIAXIAL TEST DATA, SMOOTH
  0.0,   0.0
  53.7160030052593, 0.1        
  97.7407407407408, 0.2     
  135.118798361402, 0.3       
  167.790087463557, 0.4     
  197.037037037037, 0.5    
  223.734375000000, 0.6    
  248.495013230206, 0.7     
  271.758573388203, 0.8     
  293.846916460126, 0.9     
  315.000000000000, 1.0     
  335.399632869021, 1.1     
  355.185574755823, 1.2     
  374.466672145969, 1.3     
  393.328703703704, 1.4     
  411.840000000000, 1.5     
  430.055530268548, 1.6     
  448.019915663263, 1.7     
  465.769679300292, 1.8     
  483.334946082251, 1.9     
  500.740740740741, 2.0     
  518.007988989964, 2.1     
  535.154296875000, 2.2     
  552.194562706960, 2.3     
  569.141461428862, 2.4     
  586.005830903791, 2.5     
  602.796982167353, 2.6     
  619.522950269481, 2.7     
  636.190698352530, 2.8     
  652.806284664273, 2.9     
  669.375000000000, 3.0     
*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,MARLOW,MODULI=LONG TERM
**,POISSON=0.49
*UNIAXIAL TEST DATA, SMOOTH
  0.0,   0.0
  53.7160030052593, 0.1        
  97.7407407407408, 0.2     
  135.118798361402, 0.3       
  167.790087463557, 0.4     
  197.037037037037, 0.5    
  223.734375000000, 0.6    
  248.495013230206, 0.7     
  271.758573388203, 0.8     
  293.846916460126, 0.9     
  315.000000000000, 1.0     
  335.399632869021, 1.1     
  355.185574755823, 1.2     
  374.466672145969, 1.3     
  393.328703703704, 1.4     
  411.840000000000, 1.5     
  430.055530268548, 1.6     
  448.019915663263, 1.7     
  465.769679300292, 1.8     
  483.334946082251, 1.9     
  500.740740740741, 2.0     
  518.007988989964, 2.1     
  535.154296875000, 2.2     
  552.194562706960, 2.3     
  569.141461428862, 2.4     
  586.005830903791, 2.5     
  602.796982167353, 2.6     
  619.522950269481, 2.7     
  636.190698352530, 2.8     
  652.806284664273, 2.9     
  669.375000000000, 3.0        
*VISCOELASTIC,FREQUENCY=TABULAR,PRELOAD=UNIAXIAL
6.773503,	602.1475,	1,	0	   
1.03557,	586.4045,	100,	0	   
5.46328,	391.8666,	1,	0.25	   
0.835256,	379.1688,	100,	0.25	   
4.423274,	244.4492,	1,	0.75	   
0.676254,	234.1686,	100,	0.75	   
3.948341,	184.0831,	1,	2	   
0.603644,	174.9064,	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