*HEADING
 ELEMENT PLASTICITY TEST ON - UNIAXIAL LOAD [SC8R]
*NODE
   1, 0.,0.,-0.0005
   5, 1.,0.,-0.0005
  21, 0.,1.,-0.0005
  25, 1.,1.,-0.0005
 101, 0.,0.,+0.0005
 105, 1.,0.,+0.0005
 121, 0.,1.,+0.0005
 125, 1.,1.,+0.0005
*NGEN,NSET=NX
 1,5,2
 101,105,2
*NGEN,NSET=NY
 1,21,10
 101,121,10
*NGEN
 21,25,2
 121,125,2
 5,25,10
 105,125,10
 3,23,10
 103,123,10
*NSET,NSET=FIXR
 23,25,15
 123,125,115
*NSET, GEN, NSET=NALL
 1,25,1
 101,125,1
*ELEMENT,TYPE=SC8R,ELSET=EALL
 1, 1,3,13,11, 101,103,113,111
 2, 3,5,15,13, 103,105,115,113
 3, 13,15,25,23, 113,115,125,123
 4, 11,13,23,21, 111,113,123,121
*SHELL SECTION,MATERIAL=A1,ELSET=EALL
 .001, 
*MATERIAL,NAME=A1
*ELASTIC
 1.E7, 0.3
*PLASTIC,HARDENING=KINEMATIC
 1.E4, 
 5.01E6,0.5
*EQUATION
 2, 
 25,1,1.,15,1,-1.
 2, 
 125,1,1.,115,1,-1.
 2, 
 5,1,1.,15,1,-1.
 2, 
 105,1,1.,115,1,-1.
 2, 
 15,1,1.,115,1,-1.
*AMPLITUDE,NAME=PULL
.10 , 10000. ,  .40 , 25000. ,  .60 , 100. ,  .80 , 30000.
 1. , 100.
*BOUNDARY
 NX,2
 NY,1
*RESTART,WRITE,FREQUENCY=2
*STEP,INC=100
  PULL
*STATIC,DIRECT
  .1  ,  1.0
*CLOAD,AMPLITUDE=PULL
  15,1,1.0E-3
*NODE PRINT,SUMMARY=NO
U1,U2,RF1,RF2
*EL PRINT,SUMMARY=NO
S,
SINV, 
E ,IE
PE, 
NFORC, 
*EL FILE
S, 
SINV,
E, IE
PE,PEP
NFORC,
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=EALL
 S,
 E,
 PE, 
*END STEP