*HEADING
  - INCOMPRESSIBLE HYPERELASTICITY (YEOH), PLANAR 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
1.326,-0.326,0.1319,0.0,0.0,0.0
*MULLINS EFFECT,TEST DATA INPUT,R=1.25
*UNIAXIAL TEST DATA
84.1596    , 2.         
52.2701    , 1.9        
32.8573    , 1.8        
21.2217    , 1.7        
14.1184    , 1.6        
 9.629     , 1.5        
 6.68741   , 1.4        
 4.70457   , 1.3        
 3.34586   , 1.2        
 2.41153   , 1.1        
 1.77447   , 1.         
 1.34805   , 900.E-03     
 1.06922   , 800.E-03     
 889.632E-03,700.E-03     
 770.943E-03,600.E-03     
 682.284E-03,500.E-03     
 598.896E-03,400.E-03     
 501.378E-03,300.E-03     
 375.081E-03,200.E-03     
 209.883E-03,100.E-03     
 0., 0.
*UNIAXIAL TEST DATA
509.046     ,  3.    
309.315     ,  2.85  
190.291     ,  2.7   
120.689     ,  2.55  
 79.0685    ,  2.4   
 53.0791    ,  2.25  
 36.1182    ,  2.1   
 24.644     ,  1.95  
 16.7047    ,  1.8   
 11.1661    ,  1.65  
  7.32415   ,  1.5   
  4.71004   ,  1.35  
  2.98974   ,  1.2   
  1.91099   ,  1.05  
  1.2786    ,  900.E-03
  934.216E-03, 750.E-03 
  751.632E-03, 600.E-03     
  631.26E-03 , 450.E-03     
  497.431E-03, 300.E-03     
  297.089E-03, 150.E-03     
  0., 0.
*BIAXIAL TEST DATA
483.638     ,  2.     
307.653     ,  1.9    
197.176     ,  1.8    
129.372     ,  1.7    
 87.2224    ,  1.6    
 60.1221    ,  1.5    
 42.0717    ,  1.4    
 29.5985    ,  1.3    
 20.7535    ,  1.2    
 14.3889    ,  1.1    
  9.79593   ,  1.     
  6.51097   ,  900.E-03     
  4.21196   ,  800.E-03     
  2.66104   ,  700.E-03     
  1.67218   ,  600.E-03     
  1.09216   ,  500.E-03     
  788.463E-03, 400.E-03     
  640.653E-03, 300.E-03     
  533.528E-03, 200.E-03     
  354.571E-03, 100.E-03     
  0.,0.
*BIAXIAL TEST DATA
2.52169E+03,  3.      
1.58161E+03,  2.85    
1.00121E+03,  2.7     
650.501    ,  2.55  
435.203    ,  2.4      
298.236    ,  2.25     
207.459    ,  2.1      
145.069    ,  1.95     
100.914    ,  1.8      
 69.3758   ,  1.65     
 46.6858   ,  1.5      
 30.4757   ,  1.35     
 19.1023   ,  1.2      
 11.363    ,  1.05     
  6.33716  ,  900.E-03    
  3.29594  ,  750.E-03     
  1.64788  ,  600.E-03     
  902.022E-03,450.E-03     
  635.266E-03,300.E-03     
  456.6E-03  ,150.E-03     
  0.,0.
*PLANAR TEST DATA
98.9568    ,      2.   
62.4102    ,      1.9  
39.7517    ,      1.8  
25.958     ,      1.7  
17.4286    ,      1.6  
11.9816    ,      1.5  
 8.37799   ,      1.4  
 5.92302   ,      1.3  
 4.21898   ,      1.2  
 3.02851   ,      1.1  
 2.2017    ,      1.   
 1.63739   ,      900.E-03 
 1.26259   ,      800.E-03 
 1.02129   ,      700.E-03 
 868.284E-03,     600.E-03  
 765.72E-03 ,     500.E-03  
 681.104E-03,     400.E-03  
 586.018E-03,     300.E-03  
 455.248E-03,     200.E-03  
 266.316E-03,     100.E-03  
 0.,0.
*PLANAR TEST DATA
561.053    ,  3.   
345.872    ,  2.85 
215.537    ,  2.7  
138.188    ,  2.55 
91.3716    ,  2.4  
61.8989    ,  2.25    
42.5327    ,  2.1     
29.3416    ,  1.95    
20.1378    ,  1.8     
13.6458    ,  1.65    
 9.07634   ,  1.5     
 5.90754   ,  1.35    
 3.77131   ,  1.2     
 2.3926    ,  1.05    
 1.55638   ,  900.E-03  
 1.08902   ,  750.E-03  
 847.405E-03, 600.E-03  
 711.932E-03, 450.E-03  
 581.221E-03, 300.E-03  
 368.084E-03, 150.E-03  
 0.,0.
*STEP,NLGEOM,INC=20
PLANAR TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,2.
FACE5,2,2,0.
*EL PRINT,F=10
S, 
E, 
*NODE PRINT,F=10
U,RF
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=FACE4
U,RF
*ELEMENT OUTPUT,ELSET=ONE
S,
DMENER,
SENER,
*END STEP
*STEP,NLGEOM,INC=20
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.0
FACE5,2,2,0.0
*END STEP
*STEP,NLGEOM,INC=20
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,3.
FACE5,2,2,0.
*END STEP
*STEP,NLGEOM,INC=20
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
FACE1,3
FACE3,2
FACE6,1
FACE4,1,1,0.0
FACE5,2,2,0.0
*END STEP