*HEADING
HILL PLASTICITY, UNIAXIAL TENSION TEST DIRECTION 1, C3D8
WITH LINEAR PERTURBATIONS AND MULTIPLE LOAD CASES
*NODE,NSET=ALLN
1,0.,0.,0.
2,1.,0.,0.
3,1.,1.,0.
4,0.,1.,0.
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
*ELEMENT,TYPE=C3D8,ELSET=ALLE
1,1,2,3,4,5,6,7,8
**local orientation matching global system 
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 1, 0.0
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE,ORIENT=RECT
*MATERIAL,NAME=ALLE
*ELASTIC
200.E3,.3
*PLASTIC
200.,0.
220.,.0009
220.,.0029
*POTENTIAL
1.5,1.2,1.,1.,1.,1.
*STEP,INC=10
STEP 1 - NONLINEAR
*STATIC,DIRECT
1.,10.
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,1,,.002
2,1,,.002
6,1,,.002
3,1,,.002
*EL PRINT
S, 
SINV, 
E, 
PE, 
EE, 
*NODE PRINT
U,RF
*END STEP
*STEP,PERTURBATION,INC=1
STEP 2 - LINEAR PERTURBATION
*STATIC,DIRECT
.1,.1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,1,,.0002
2,1,,.0002
6,1,,.0002
3,1,,.0002
*EL PRINT
S, 
SINV, 
E, 
*NODE PRINT
U,RF
*END STEP
*STEP,INC=10
STEP 3 - REST OF NONLINEAR
*STATIC,DIRECT
1.,10.
*BOUNDARY
7,1,,.004
2,1,,.004
6,1,,.004
3,1,,.004
*END STEP

*STEP,PERTURBATION,INC=1
STEP 4 - MULTIPLE LOAD CASES
*STATIC,DIRECT
.1,.1
*LOAD CASE,NAME=STRETCH-1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,1,,.0004
2,1,,.0004
6,1,,.0004
3,1,,.0004
*END LOAD CASE
*LOAD CASE,NAME=COMPRESS-1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,1,,-.0004
2,1,,-.0004
6,1,,-.0004
3,1,,-.0004
*END LOAD CASE
*LOAD CASE,NAME=STRETCH-2
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,1,,.0004
2,1,,.0004
6,1,,.0004
3,1,,.0004
*END LOAD CASE
*LOAD CASE,NAME=COMPRESS-2
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,1,,-.0004
2,1,,-.0004
6,1,,-.0004
3,1,,-.0004
*END LOAD CASE
*EL PRINT
S, 
SINV, 
E, 
*NODE PRINT
U,RF
*END STEP