*HEADING
  - INCOMPRESSIBLE HYPERELASTICITY (OGDEN), 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=MOONEY
*MATERIAL,NAME=MOONEY
*Hyperelastic, yeoh
 1.326,   -0.326,   0.1319, 0.000725,       0.,       0.
*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
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,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,
*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