*HEADING
  - INCOMPRESSIBLE HYPERELASTICITY (YEOH), BIAXIAL TENSION
*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.
*NSET,NSET=FACE1
1,2,3,4
*NSET,NSET=FACE2
5,6,7,8
*NSET,NSET=FACE3
1,2,5,6
*NSET,NSET=FACE4
2
*NSET,NSET=FACE42
3,6,7
*NSET,NSET=FACE5
3
*NSET,NSET=FACE52
4,7,8
*NSET,NSET=FACE6
4,1,8,5
*EQUATION
** Since the S11 output is Cauchy or true stress, we need to
** determine the nominal stress for post-processing.
** Nodes 3,6,7 are tied to node 2 in dof 1 so that:
** Nominal stress (dof 1) = RF1 (@ node 2) / Original area
** (w/c is 1 x 1 = 1)
2
FACE42,1,1, 2,1,-1
*EQUATION
2
FACE52,2,1,3,2,-1
*ELEMENT,TYPE=C3D8RH,ELSET=ONE
1,1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ONE,MATERIAL=YEOH
*MATERIAL,NAME=YEOH
*HYPERELASTIC,YEOH,TEST DATA
*UNIAXIAL TEST DATA
   140.704    ,  0.102558    
   218.612    ,  0.214180    
   283.232    ,  0.323707    
   345.659    ,  0.431249    
   410.154    ,  0.535531    
   476.085    ,  0.635031    
   545.034    ,  0.732017    
   616.713    ,  0.826472    
   691.090    ,  0.918368    
   768.445    ,   1.00796    
   848.693    ,   1.09535    
   931.569    ,   1.18236    
   1015.91    ,   1.26772    
   1101.28    ,   1.35419    
   1188.59    ,   1.44020    
   1276.28    ,   1.52645    
   1363.51    ,   1.61423    
   1449.85    ,   1.70203    
   1535.32    ,   1.79182    
   1618.30    ,   1.88415    
*MULLINS EFFECT,TEST DATA
*INCLUDE,INPUT=mullins_calibrate_testdata.inp
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.099
*EL PRINT,F=10
S, 
E, 
*NODE PRINT,F=10
U,RF
*EL FILE,F=10
S,E
*NODE FILE,F=10
U,RF
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=FACE4
U,RF
*ELEMENT OUTPUT,ELSET=ONE
S,
DMENER,
SENER,
*END STEP
*STEP,NLGEOM
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.0
*END STEP
*STEP,NLGEOM
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.26
*END STEP
*STEP,NLGEOM
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.0
*END STEP
*STEP,NLGEOM
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.51
*END STEP
*STEP,NLGEOM
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.0
*END STEP