*Heading
 Test for pressure and triaxiality dependent hardening using *USER DEFINED FIELD, TYPE=SPECIFIELD
*Preprint, echo=NO, model=YES, history=YES, contact=NO
*Node, nset=nall
     11,          -5.,           0.,          10.
     12,          -5.,           5.,          10.
     13,          -5.,           0.,           0.
     14,          -5.,           5.,           0.
     15,           0.,           0.,          10.
     16,           0.,           5.,          10.
     17,           0.,           0.,           0.
     18,           0.,           5.,           0.
**
     21,          -5.,           0.,          30.
     22,          -5.,           5.,          30.
     23,          -5.,           0.,          20.
     24,          -5.,           5.,          20.
     25,           0.,           0.,          30.
     26,           0.,           5.,          30.
     27,           0.,           0.,          20.
     28,           0.,           5.,          20.
**
     31,          -5.,           0.,          50.
     32,          -5.,           5.,          50.
     33,          -5.,           0.,          40.
     34,          -5.,           5.,          40.
     35,           0.,           0.,          50.
     36,           0.,           5.,          50.
     37,           0.,           0.,          40.
     38,           0.,           5.,          40.
**
     41,          -5.,           0.,          70.
     42,          -5.,           5.,          70.
     43,          -5.,           0.,          60.
     44,          -5.,           5.,          60.
     45,           0.,           0.,          70.
     46,           0.,           5.,          70.
     47,           0.,           0.,          60.
     48,           0.,           5.,          60.
*Element, type=C3D8R, elset=e1
 11, 15, 16, 18, 17, 11, 12, 14, 13
*Element, type=C3D8R, elset=e2
 21, 25, 26, 28, 27, 21, 22, 24, 23
*Element, type=C3D8R, elset=e3
 31, 35, 36, 38, 37, 31, 32, 34, 33
*Element, type=C3D8R, elset=e4
 41, 45, 46, 48, 47, 41, 42, 44, 43
*Solid Section, elset=e1, material=MAT1
*Solid Section, elset=e2, material=MAT2
*Solid Section, elset=e3, material=MAT3
*Solid Section, elset=e4, material=MAT4
**  
*Nset, nset=NXSYMM, generate
 11,  14,  1
 21,  24,  1
 31,  34,  1
 41,  44,  1
*Nset, nset=NYSYMM, generate
 11,  17,  2
 21,  27,  2
 31,  37,  2
 41,  47,  2
*Nset, nset=NZSYMM
 13, 14, 17, 18
 23, 24, 27, 28
 33, 34, 37, 38
 43, 44, 47, 48
*Elset, elset=eall
 e1, e2, e3, e4
*Surface, type=ELEMENT, name=Surf
 eall, S3
*Amplitude, name=LOAD, definition=SMOOTH STEP
 0., 0.,  1., 1.
*Amplitude, name=UNLOAD, definition=SMOOTH STEP
 0., 1.,  1., 0.
** 
*Material, name=MAT1
*Density
 7.8e-06,
*Elastic
 200000., 0.3
*Plastic, dependencies=1
 100., 0., ,-1.
 150., 1., ,-1.
 350., 3., ,-1.
 200., 0., , 1.
 300., 1., , 1.
 500., 3., , 1.
*User defined field,type=specified
 1, PRESS
**
*Material, name=MAT2
*Density
 7.8e-06,
*Elastic
 200000., 0.3
*Plastic, dependencies=2
 100., , ,0.,-1.
 150., , ,1.,-1.
 350., , ,3.,-1.
 200., , ,0., 1.
 300., , ,1., 1.
 500., , ,3., 1.
*User defined field,type=specified
 1, PEEQ
 2, PRESS
** 
*Material, name=MAT3
*Density
 7.8e-06,
*Elastic
 200000., 0.3
*Plastic, dependencies=1
 200., 0., ,-0.3
 300., 1., ,-0.3
 500., 3., ,-0.3
 100., 0., , 0.3
 150., 1., , 0.3
 350., 3., , 0.3
*User defined field,type=specified
 1, TRIAX
**
*Material, name=MAT4
*Density
 7.8e-06,
*Elastic
 200000., 0.3
*Plastic, dependencies=2
 200., , ,0.,-0.3
 300., , ,1.,-0.3
 500., , ,3.,-0.3
 100., , ,0., 0.3
 150., , ,1., 0.3
 350., , ,3., 0.3
*User defined field,type=specified
 1, PEEQ
 2, TRIAX
** 
*Boundary
 NXSYMM , XSYMM
*Boundary
 NYSYMM, YSYMM 
*Boundary
 NZSYMM, ZSYMM
** 
*Step, name=Step-1
*Dynamic, Explicit
 , 1.
*Bulk Viscosity
 0.06, 1.2
**
*Dsload, amplitude=LOAD, constant resultant=yes
 Surf, P, 250.
** 
*Output, field
*Node Output
 RF, U
*Element Output, directions=YES
 FV, LE, P, PE, PEEQ, S
*Contact Output
 CSTRESS, 
**
*Output, history
*Energy Output, Variable=All
** 
*Node Output, nset=NALL
 UT, 
** 
*Element Output, elset=eall
 FV, MISES, PEEQ, S33
 PRESS, TRIAX
*End Step
** 
*Step, name=Step-2
*Dynamic, Explicit
 , 1.
*Bulk Viscosity
 0.06, 1.2
** 
*Dsload, amplitude=UNLOAD, constant resultant=yes
 Surf, P, 250.
** 
*Output, field
*Node Output
 RF, U
*Element Output, directions=YES
 FV, LE, P, PE, PEEQ, S
*Contact Output
 CSTRESS, 
*Output, history
*Energy Output, Variable=Preselect
*Node Output, nset=NALL
 UT, 
*Element Output, elset=eall
 FV, MISES, PEEQ, S33
 PRESS, TRIAX
*End Step
** 
*Step, name=Step-3
*Dynamic, Explicit
 , 1.
*Bulk Viscosity
 0.06, 1.2
**
*Dsload, amplitude=LOAD,constant resultant=yes
 Surf, P, -150.
** 
*Output, field
*Node Output
 RF, U
*Element Output, directions=YES
 FV, LE, P, PE, PEEQ, S
*Contact Output
 CSTRESS, 
** 
*Output, history
*Energy Output, Variable=Preselect
*Node Output, nset=NALL
 UT, 
*Element Output, elset=eall
 FV, MISES, PEEQ, S33
 PRESS, TRIAX
*End Step