*HEADING
RATE DEPENDENT HILL PLASTICITY, UNIAXIAL TENSION DIRECTION 1, C3D8
AND MULTIPLE LOAD CASES TEST FOR HILL'S PLASTICITY
*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
*Distribution Table, name=distTableElastic
 modulus, ratio
*Distribution, name=distElastic, table=distTableElastic, location=element
 , 10.e3, 0.25
 alle, 200.0e3, 0.3
**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
distElastic,
*PLASTIC
200.,0.
220.,.0009
220.,.0029
*POTENTIAL
1.5,1.2,1.,1.,1.,1.
*RATE DEPENDENT,TYPE=YIELD RATIO
1.0 , 0.0
1.190, 0.01
1.302, 0.1000    
1.478,  1.000    
1.758,  10.00    
2.201,  100.0    
2.904,  1000.    
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
*STEP,INC=20
*STATIC,DIRECT
1.E-3,20.E-3
*BOUNDARY
7,1,,.004
2,1,,.004
6,1,,.004
3,1,,.004
*EL PRINT
S, 
SINV, 
E, 
PE, 
EE, 
*NODE PRINT
U,RF
*END STEP
**
*STEP,perturbation
*STATIC
1,1
*LOAD CASE,name=STRETCH-1
*BOUNDARY
7,1,,.004
2,1,,.004
6,1,,.004
3,1,,.004
*END LOAD CASE
*LOAD CASE,name=STRETCH-2
*BOUNDARY
7,1,,.008
2,1,,.008
6,1,,.008
3,1,,.008
*END LOAD CASE
*EL PRINT
S, 
SINV, 
E, 
PE, 
EE, 
*NODE PRINT
U,RF
*END STEP