*heading
 Verification tests for user defined hardening (VUHARD)
*preprint,model=yes
*node,nset=nall
**
**
**
*nset, nset=nout
*nset,nset=face1
*nset,nset=face2
*nset,nset=face3
*nset,nset=face4
*nset,nset=face5
*nset,nset=face6
*element,type=c3d8r,elset=Mises
*element,type=c3d8r,elset=Hill
*element,type=c3d8r,elset=OrthoElasPlas
*element,type=c3d8r,elset=EosUsupMises
*element,type=c3d8r,elset=EosJwlMises
*element,type=c3d8r,elset=EosPalphaMises
*element,type=c3d8r,elset=EosGasMises
*element,type=c3d8r,elset=FeFp
*element,type=c3d8r,elset=CombHard
*element,type=c3d8r,elset=Gurson
*element,type=c3d8r,elset=Mises_ref
*element,type=c3d8r,elset=Hill_ref
*element,type=c3d8r,elset=OrthoElasPlas_ref
*element,type=c3d8r,elset=EosUsupMises_ref
*element,type=c3d8r,elset=EosJwlMises_ref
*element,type=c3d8r,elset=EosPalphaMises_ref
*element,type=c3d8r,elset=EosGasMises_ref
*element,type=c3d8r,elset=FeFp_ref
*element,type=c3d8r,elset=CombHard_ref
*element,type=c3d8r,elset=Gurson_ref
**
*elset,elset=eall
 Mises,Hill,OrthoElasPlas,
 EosUsupMises,EosJwlMises,EosPalphaMises,EosGasMises
 FeFp,CombHard,Gurson
 Mises_ref,Hill_ref,OrthoElasPlas_ref,
 EosUsupMises_ref,EosJwlMises_ref,EosPalphaMises_ref,EosGasMises_ref
 FeFp_ref,CombHard_ref,Gurson_ref
**
*orientation,name=Ori-1
**
*solid section,elset=mises,material=Mises
 1.0
*solid section,elset=hill,material=Hill,orientation=Ori-1
 1.0
*solid section,elset=orthoelasplas,material=OrthoElasPlas,orientation=Ori-1
 1.0
*solid section,elset=EosUsupMises,material=EosUsupMises
 1.0
*solid section,elset=EosJwlMises,material=EosJwlMises
 1.0
*solid section,elset=EosPalphaMises,material=EosPalphaMises
 1.0
*solid section,elset=EosGasMises,material=EosGasMises
 1.0
*solid section,elset=FeFp,material=FeFp
 1.0
*solid section,elset=CombHard,material=CombHard
 1.0
*solid section,elset=Gurson,material=Gurson
 1.0
*solid section,elset=mises_ref,material=Mises_ref
 1.0
*solid section,elset=hill_ref,material=Hill_ref,orientation=Ori-1
 1.0
*solid section,elset=orthoelasplas_ref,material=OrthoElasPlas_ref,orientation=Ori-1
 1.0
*solid section,elset=EosUsupMises_ref,material=EosUsupMises_ref
 1.0
*solid section,elset=EosJwlMises_ref,material=EosJwlMises_ref
 1.0
*solid section,elset=EosPalphaMises_ref,material=EosPalphaMises_ref
 1.0
*solid section,elset=EosGasMises_ref,material=EosGasMises_ref
 1.0
*solid section,elset=FeFp_ref,material=FeFp_ref
 1.0
*solid section,elset=CombHard_ref,material=CombHard_ref
 1.0
*solid section,elset=Gurson_ref,material=Gurson_ref
 1.0
**
*material,name=Mises
*density
 2.070e3,
*elastic
 120e6, 0.3
*plastic,hardening=user,properties=208
 1.970e+06, 0.000
 2.701e+06, 0.004
 2.882e+06, 0.008
 3.009e+06, 0.012
 3.109e+06, 0.016
 3.193e+06, 0.020
 3.348e+06, 0.029
 3.472e+06, 0.038
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=Mises_ref
*density
 2.070e3,
*elastic
 120e6, 0.3
*plastic,hardening=isotropic
 1.970e+06, 0.000
 2.701e+06, 0.004
 2.882e+06, 0.008
 3.009e+06, 0.012
 3.109e+06, 0.016
 3.193e+06, 0.020
 3.348e+06, 0.029
 3.472e+06, 0.038
**
*material,name=Hill
*density
 2.070e3,
*elastic
 120e6, 0.3
*plastic,hardening=user,properties=208
 1.970e+06, 0.000
 2.701e+06, 0.004
 2.882e+06, 0.008
 3.009e+06, 0.012
 3.109e+06, 0.016
 3.193e+06, 0.020
 3.348e+06, 0.029
 3.472e+06, 0.038
*potential
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=Hill_ref
*density
 2.070e3,
*elastic
 120e6, 0.3
*plastic,hardening=isotropic
 1.970e+06, 0.000
 2.701e+06, 0.004
 2.882e+06, 0.008
 3.009e+06, 0.012
 3.109e+06, 0.016
 3.193e+06, 0.020
 3.348e+06, 0.029
 3.472e+06, 0.038
*potential
**
*material,name=OrthoElasPlas
*density
 500.,
*elastic, type=engineering constants
 3.0e6,3.0e6,3.0e6,0.3,0.3,0.3,1153846.154,1153846.154,
 1153846.154
*plastic,hardening=user,properties=16
 0.15e5, 0.0
 0.6e5,  1.0
*potential
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
**
*material,name=OrthoElasPlas_ref
*density
 500.,
*elastic, type=engineering constants
 3.0e6,3.0e6,3.0e6,0.3,0.3,0.3,1153846.154,1153846.154,
 1153846.154
*plastic,hardening=isotropic
 0.15e5, 0.0
 0.6e5,  1.0
*potential
**
*material,name=EosUsupMises
*density
 1.9e3,
*eos,type=usup
*elastic,type=shear
 80.232558e9,
*plastic,hardening=user,properties=8
 175.e6, 0.0
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=EosUsupMises_ref
*density
 1.9e3,
*eos,type=usup
*elastic,type=shear
 80.232558e9,
*plastic,hardening=isotropic
 175.e6, 0.0
**
*material,name=EosJwlMises
*density
 1.9e3,
*eos,type=jwl
*detonation point
*elastic,type=shear
 80.232558e9,
*plastic,hardening=user,properties=8
 175.e6
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=EosJwlMises_ref
*density
 1.9e3,
*eos,type=jwl
*detonation point
*elastic,type=shear
 80.232558e9,
*plastic,hardening=isotropic
 175.e6
**
*material,name=eospalphamises
*density
 2.070e3,
*eos, type=usup
*eos compaction
*elastic,type=shear
 124e6, 0.3
*plastic,hardening=user,properties=208
 1.970E+06, 0.000
 2.701E+06, 0.004
 2.882E+06, 0.008
 3.009E+06, 0.012
 3.109E+06, 0.016
 3.193E+06, 0.020
 3.348E+06, 0.029
 3.472E+06, 0.038
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=eospalphamises_ref
*density
 2.070e3,
*eos, type=usup
*eos compaction
*elastic,type=shear
 124e6, 0.3
*plastic,hardening=isotropic
 1.970E+06, 0.000
 2.701E+06, 0.004
 2.882E+06, 0.008
 3.009E+06, 0.012
 3.109E+06, 0.016
 3.193E+06, 0.020
 3.348E+06, 0.029
 3.472E+06, 0.038
**
*material,name=EosGasMises
*eos,type=ideal gas
*elastic,type=shear
 80.232558e9,
*specific heat
 0.3,
*density
7800
*plastic,hardening=user,properties=8
175.e6
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=EosGasMises_ref
*eos,type=ideal gas
*elastic,type=shear
 80.232558e9,
*specific heat
 0.3,
*density
7800
*plastic,hardening=isotropic
175.e6
**
*material,name=CombHard
*density
 1.0
*elastic
 200.e3,.3
*plastic,hardening=combined,data type=parameters
 200.,0.,0.
*cyclic hardening,user,properties=8
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=CombHard_ref
*density
 1.0
*elastic
 200.e3,.3
*plastic,hardening=combined,data type=parameters
 200.,0.,0.
*cyclic hardening
**
*material,name=Gurson
*density
 2.070e3,
*elastic
 120e6, 0.3
*plastic,hardening=user,properties=208
 1.970e+06, 0.000
 2.701e+06, 0.004
 2.882e+06, 0.008
 3.009e+06, 0.012
 3.109e+06, 0.016
 3.193e+06, 0.020
 3.348e+06, 0.029
 3.472e+06, 0.038
*porous metal plasticity,relative density=1.
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=Gurson_ref
*density
 2.070e3,
*elastic
 120e6, 0.3
*plastic,hardening=isotropic
 1.970e+06, 0.000
 2.701e+06, 0.004
 2.882e+06, 0.008
 3.009e+06, 0.012
 3.109e+06, 0.016
 3.193e+06, 0.020
 3.348e+06, 0.029
 3.472e+06, 0.038
*porous metal plasticity,relative density=1.
**
*material,name=FeFp
*density
 1.0e-3
*hyperelastic,mooney-rivlin
 0.749527319, 0.496199379, 1.616263681E-02
*mullins effect
*plastic,hardening=user,properties=8
 1.5
*depvar,delete=3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=FeFp_ref
*density
 1.0e-3
*hyperelastic,mooney-rivlin
 0.749527319, 0.496199379, 1.616263681E-02
*mullins effect
*plastic,hardening=isotropic
 1.5
**
*amplitude,name=bc
 0.0,0.0,0.5,0.75
*step
uniaxial tension
*dynamic, explicit
 ,0.125
*boundary
 face1,3
 face2,3
 face3,2
 face6,1
*boundary,amp=bc
 face4,1,1,1.0
*output,history,freq=1
*element output, elset=eall
 s11,pe11,ne11,peeq,mises,edt,sdv
*node output,nset=nout
 u1,rf1
*output,field, time marks=yes
*element output
 s,ne,pe,peeq,sdv
*node output
 u,rf
*end step