*HEADING
  - INCOMPRESSIBLE HYPERELASTICITY (YEOH), UNIAXIAL TENSION
*PREPRINT,MODEL=YES
*NODE,NSET=ALL
1,
2,1.
3,1.,1.,
4,0.,1.,
11,
21,1.
31,1.,1.,
41,0.,1.,
*ELEMENT,TYPE=CPS4R,ELSET=ONE
1,1,2,3,4
*ELEMENT,TYPE=CPS4R,ELSET=TWO
2,11,21,31,41
*NSET,NSET=SIDE1
1,2
*NSET,NSET=SIDE2
2,3
*NSET,NSET=SIDE3
3,4
*NSET,NSET=SIDE4
4,1
*NSET,NSET=SIDE11
11,21
*NSET,NSET=SIDE21
21,31
*NSET,NSET=SIDE31
31,41
*NSET,NSET=SIDE41
41,11
*SOLID SECTION,ELSET=ONE,MATERIAL=YEOH1
*SOLID SECTION,ELSET=TWO,MATERIAL=YEOH2
*MATERIAL,NAME=YEOH1
*Hyperelastic, yeoh
 1.326,   -0.326,   0.1319, 0.,       0.,       0.
*MULLINS EFFECT,USER,PROPERTIES=3
1.1,100.,0.
*MATERIAL,NAME=YEOH2
*Hyperelastic, yeoh
 1.326,   -0.326,   0.1319, 0.,       0.,       0.
*MULLINS EFFECT,USER,PROPERTIES=3
1.1,100.,0.
*DEPVAR
1
*ELSET,ELSET=ALL
1,2
*INITIAL CONDITIONS,TYPE=SOLUTION
TWO,283.01
*STEP,NLGEOM,INC=20
UNIAXIAL TENSION
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,1,1,2.
**
SIDE41,1
SIDE11,2
SIDE21,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=ALL
S,
*END STEP
*STEP
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,1,1,
**
SIDE41,1
SIDE11,2
SIDE21,1,1,
*END STEP
*STEP
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,1,1,3.
**
SIDE41,1
SIDE11,2
SIDE21,1,1,3.
*END STEP
*STEP
*STATIC,DIRECT
1.,20.
*BOUNDARY,OP=NEW
SIDE4,1
SIDE1,2
SIDE2,1,1,
**
SIDE41,1
SIDE11,2
SIDE21,1,1,
*END STEP