*HEADING
  - COMPRESSIBLE HYPERELASTICITY (MARLOW), UNIAXIAL TENSION
*PREPRINT,MODEL=YES
*RESTART,WRITE
*NODE,NSET=ALL
1,
2,1.
3,1.,1.,
4,0.,1.,
*ELEMENT,TYPE=CPS4R,ELSET=ONE
1,1,2,3,4
*NSET,NSET=SIDE1
1,2
*NSET,NSET=SIDE2
2,3
*NSET,NSET=SIDE3
3,4
*NSET,NSET=SIDE4
4,1
*SOLID SECTION,ELSET=ONE,MATERIAL=MARLOW
*MATERIAL,NAME=MARLOW
*HYPERELASTIC,MARLOW
*UNIAXIAL TEST DATA
  0.000000000000000,       0.000000000000000
  0.231494000000000,       0.030000000000000
  0.448560000000000,       0.060000000000000
  0.747177000000000,       0.105000000000000     
   1.13595000000000,       0.172500000000000     
   1.59760000000000,       0.273750000000000     
   2.08089000000000,       0.423750000000000     
   2.47973000000000,       0.573750000000000     
   3.03767000000000,       0.723750000000000     
   4.07610000000000,       0.873750000000000     
   5.99936000000000,        1.02375000000000     
   9.30106000000000,        1.17375000000000     
   14.5712000000000,        1.32375000000000     
   22.5031000000000,        1.47375000000000     
   33.9013000000000,        1.62375000000000     
   49.6884000000000,        1.77375000000000     
   70.9126000000000,        1.92375000000000     
   98.7549000000000,        2.07375000000000     
   134.537000000000,        2.22375000000000     
   179.726000000000,        2.37375000000000     
   235.948000000000,        2.52375000000000     
   304.986000000000,        2.67375000000000     
   388.797000000000,        2.82375000000000     
   489.510000000000,        2.97375000000000     
   509.043000000000,        3.00000000000000     
*VOLUMETRIC TEST DATA
  0.000000000000000,       1.00000000000000     
   82.7586206896552,       0.970000000000000     
   165.517241379310,       0.940000000000000     
   248.275862068965,       0.910000000000000     
   331.034482758621,       0.880000000000000     
   413.793103448276,       0.850000000000000     
   496.551724137931,       0.820000000000000     
   579.310344827586,       0.790000000000000     
   662.068965517241,       0.760000000000000     
   744.827586206897,       0.730000000000000     
   827.586206896552,       0.700000000000000     
*Mullins effect,dependencies=1
1.01,91.,0.01,0.,0.
1.09,99.,0.09,40.,0.
1.11,101.,0.11,0.,1.
1.19,109.,0.19,40.,1.
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL,20.
*INITIAL CONDITIONS,TYPE=FIELD
ALL,0.5
*EQUATION
 2,
 2,1,1,3,1,-1
*nset,nset=n1
3
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
**SIDE2,1,1,2.
3,1,1,2.
*EL PRINT,F=10
S, 
E, 
*NODE PRINT,F=10
U,RF
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
U,RF
*OUTPUT,HISTORY,FREQ=1
*ELEMENT OUTPUT,ELSET=ONE
S,
*NODE OUTPUT,NSET=n1
U1,RF1
*END STEP
*STEP
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,1,1,
*END STEP
*STEP
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,1,1,3.
*END STEP
*STEP
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,1,1,
*END STEP