*heading
 Verification tests for user defined hardening (VUHARD)
*preprint,model=yes
*node,nset=nall
 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.
 201,0.,10.,0.
 202,1.,10.,0.
 203,1.,11.,0.
 204,0.,11.,0.
 205,0.,10.,1.
 206,1.,10.,1.
 207,1.,11.,1.
 208,0.,11.,1.
**
 11,2.,0.,0.
 12,3.,0.,0.
 13,3.,1.,0.
 14,2.,1.,0.
 15,2.,0.,1.
 16,3.,0.,1.
 17,3.,1.,1.
 18,2.,1.,1.
 211,2.,10.,0.
 212,3.,10.,0.
 213,3.,11.,0.
 214,2.,11.,0.
 215,2.,10.,1.
 216,3.,10.,1.
 217,3.,11.,1.
 218,2.,11.,1.
**
 21,4.,0.,0.
 22,5.,0.,0.
 23,5.,1.,0.
 24,4.,1.,0.
 25,4.,0.,1.
 26,5.,0.,1.
 27,5.,1.,1.
 28,4.,1.,1.
 221,4.,10.,0.
 222,5.,10.,0.
 223,5.,11.,0.
 224,4.,11.,0.
 225,4.,10.,1.
 226,5.,10.,1.
 227,5.,11.,1.
 228,4.,11.,1.
**
 31,6.,0.,0.
 32,7.,0.,0.
 33,7.,1.,0.
 34,6.,1.,0.
 35,6.,0.,1.
 36,7.,0.,1.
 37,7.,1.,1.
 38,6.,1.,1.
 231,6.,10.,0.
 232,7.,10.,0.
 233,7.,11.,0.
 234,6.,11.,0.
 235,6.,10.,1.
 236,7.,10.,1.
 237,7.,11.,1.
 238,6.,11.,1.
**
 41,8.,0.,0.
 42,9.,0.,0.
 43,9.,1.,0.
 44,8.,1.,0.
 45,8.,0.,1.
 46,9.,0.,1.
 47,9.,1.,1.
 48,8.,1.,1.
 241,8.,10.,0.
 242,9.,10.,0.
 243,9.,11.,0.
 244,8.,11.,0.
 245,8.,10.,1.
 246,9.,10.,1.
 247,9.,11.,1.
 248,8.,11.,1.
**
 51,10.,0.,0.
 52,11.,0.,0.
 53,11.,1.,0.
 54,10.,1.,0.
 55,10.,0.,1.
 56,11.,0.,1.
 57,11.,1.,1.
 58,10.,1.,1.
 251,10.,10.,0.
 252,11.,10.,0.
 253,11.,11.,0.
 254,10.,11.,0.
 255,10.,10.,1.
 256,11.,10.,1.
 257,11.,11.,1.
 258,10.,11.,1.
**
 61,12.,0.,0.
 62,13.,0.,0.
 63,13.,1.,0.
 64,12.,1.,0.
 65,12.,0.,1.
 66,13.,0.,1.
 67,13.,1.,1.
 68,12.,1.,1.
 261,12.,10.,0.
 262,13.,10.,0.
 263,13.,11.,0.
 264,12.,11.,0.
 265,12.,10.,1.
 266,13.,10.,1.
 267,13.,11.,1.
 268,12.,11.,1.
**
 71,14.,0.,0.
 72,15.,0.,0.
 73,15.,1.,0.
 74,14.,1.,0.
 75,14.,0.,1.
 76,15.,0.,1.
 77,15.,1.,1.
 78,14.,1.,1.
 271,14.,10.,0.
 272,15.,10.,0.
 273,15.,11.,0.
 274,14.,11.,0.
 275,14.,10.,1.
 276,15.,10.,1.
 277,15.,11.,1.
 278,14.,11.,1.
**
 81,16.,0.,0.
 82,17.,0.,0.
 83,17.,1.,0.
 84,16.,1.,0.
 85,16.,0.,1.
 86,17.,0.,1.
 87,17.,1.,1.
 88,16.,1.,1.
 281,16.,10.,0.
 282,17.,10.,0.
 283,17.,11.,0.
 284,16.,11.,0.
 285,16.,10.,1.
 286,17.,10.,1.
 287,17.,11.,1.
 288,16.,11.,1.
**
 91,18.,0.,0.
 92,19.,0.,0.
 93,19.,1.,0.
 94,18.,1.,0.
 95,18.,0.,1.
 96,19.,0.,1.
 97,19.,1.,1.
 98,18.,1.,1.
 291,18.,10.,0.
 292,19.,10.,0.
 293,19.,11.,0.
 294,18.,11.,0.
 295,18.,10.,1.
 296,19.,10.,1.
 297,19.,11.,1.
 298,18.,11.,1.
*nset, nset=nout
 2,12,22,32,42,52,62,72,82,92
*nset,nset=face1
  1, 2, 3, 4
 11,12,13,14
 21,22,23,24
 31,32,33,34
 41,42,43,44
 51,52,53,54
 61,62,63,64
 71,72,73,74
 81,82,83,84
 91,92,93,94
 201,202,203,204
 211,212,213,214
 221,222,223,224
 231,232,233,234
 241,242,243,244
 251,252,253,254
 261,262,263,264
 271,272,273,274
 281,282,283,284
 291,292,293,294
*nset,nset=face2
  5, 6, 7, 8
 15,16,17,18
 25,26,27,28
 35,36,37,38
 45,46,47,48
 55,56,57,58
 65,66,67,68
 75,76,77,78
 85,86,87,88
 95,96,97,98
 205,206,207,208
 215,216,217,218
 225,226,227,228
 235,236,237,238
 245,246,247,248
 255,256,257,258
 265,266,267,268
 275,276,277,278
 285,286,287,288
 295,296,297,298
*nset,nset=face3
  1, 2, 5, 6
 11,12,15,16
 21,22,25,26
 31,32,35,36
 41,42,45,46
 51,52,55,56
 61,62,65,66
 71,72,75,76
 81,82,85,86
 91,92,95,96
 201,202,205,206
 211,212,215,216
 221,222,225,226
 231,232,235,236
 241,242,245,246
 251,252,255,256
 261,262,265,266
 271,272,275,276
 281,282,285,286
 291,292,295,296
*nset,nset=face4
  2, 3, 6, 7
 12,13,16,17
 22,23,26,27
 32,33,36,37
 42,43,46,47
 52,53,56,57
 62,63,66,67
 72,73,76,77
 82,83,86,87
 92,93,96,97
 202,203,206,207
 212,213,216,217
 222,223,226,227
 232,233,236,237
 242,243,246,247
 252,253,256,257
 262,263,266,267
 272,273,276,277
 282,283,286,287
 292,293,296,297
*nset,nset=face5
  3, 4, 7, 8
 13,14,17,18
 23,24,27,28
 33,34,37,38
 43,44,47,48
 53,54,57,58
 63,64,67,68
 73,74,77,78
 83,84,87,88
 93,94,97,98
 203,204,207,208
 213,214,217,218
 223,224,227,228
 233,234,237,238
 243,244,247,248
 253,254,257,258
 263,264,267,268
 273,274,277,278
 283,284,287,288
 293,294,297,298
*nset,nset=face6
  4, 1, 8, 5
 14,11,18,15
 24,21,28,25
 34,31,38,35
 44,41,48,45
 54,51,58,55
 64,61,68,65
 74,71,78,75
 84,81,88,85
 94,91,98,95
 204,201,208,205
 214,211,218,215
 224,221,228,225
 234,231,238,235
 244,241,248,245
 254,251,258,255
 264,261,268,265
 274,271,278,275
 284,281,288,285
 294,291,298,295
*element,type=c3d8r,elset=Mises
 1,1,2,3,4,5,6,7,8
*element,type=c3d8r,elset=Hill
 11,11,12,13,14,15,16,17,18
*element,type=c3d8r,elset=OrthoElasPlas
 21,21,22,23,24,25,26,27,28
*element,type=c3d8r,elset=EosUsupMises
 31,31,32,33,34,35,36,37,38
*element,type=c3d8r,elset=EosJwlMises
 41,41,42,43,44,45,46,47,48
*element,type=c3d8r,elset=EosPalphaMises
 51,51,52,53,54,55,56,57,58
*element,type=c3d8r,elset=EosGasMises
 61,61,62,63,64,65,66,67,68
*element,type=c3d8r,elset=FeFp
 71,71,72,73,74,75,76,77,78
*element,type=c3d8r,elset=CombHard
 81,81,82,83,84,85,86,87,88
*element,type=c3d8r,elset=Gurson
 91,91,92,93,94,95,96,97,98
*element,type=c3d8r,elset=Mises_ref
 201,201,202,203,204,205,206,207,208
*element,type=c3d8r,elset=Hill_ref
 211,211,212,213,214,215,216,217,218
*element,type=c3d8r,elset=OrthoElasPlas_ref
 221,221,222,223,224,225,226,227,228
*element,type=c3d8r,elset=EosUsupMises_ref
 231,231,232,233,234,235,236,237,238
*element,type=c3d8r,elset=EosJwlMises_ref
 241,241,242,243,244,245,246,247,248
*element,type=c3d8r,elset=EosPalphaMises_ref
 251,251,252,253,254,255,256,257,258
*element,type=c3d8r,elset=EosGasMises_ref
 261,261,262,263,264,265,266,267,268
*element,type=c3d8r,elset=FeFp_ref
 271,271,272,273,274,275,276,277,278
*element,type=c3d8r,elset=CombHard_ref
 281,281,282,283,284,285,286,287,288
*element,type=c3d8r,elset=Gurson_ref
 291,291,292,293,294,295,296,297,298
**
*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
 1.0,0.0,0.0,0.0,1.0,0.0
 3,0.0
**
*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
 3.578e+06, 0.047
 3.671e+06, 0.056
 3.754e+06, 0.065
 3.829e+06, 0.074
 3.899e+06, 0.083
 3.964e+06, 0.092
 4.024e+06, 0.101
 4.081e+06, 0.110
 4.135e+06, 0.119
 4.186e+06, 0.128
 4.234e+06, 0.137
 4.281e+06, 0.146
 4.326e+06, 0.155
 4.369e+06, 0.164
 4.410e+06, 0.173
 4.450e+06, 0.182
 4.489e+06, 0.191
 4.526e+06, 0.200
*depvar,delete=3
 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
 3.578e+06, 0.047
 3.671e+06, 0.056
 3.754e+06, 0.065
 3.829e+06, 0.074
 3.899e+06, 0.083
 3.964e+06, 0.092
 4.024e+06, 0.101
 4.081e+06, 0.110
 4.135e+06, 0.119
 4.186e+06, 0.128
 4.234e+06, 0.137
 4.281e+06, 0.146
 4.326e+06, 0.155
 4.369e+06, 0.164
 4.410e+06, 0.173
 4.450e+06, 0.182
 4.489e+06, 0.191
 4.526e+06, 0.200
**
*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
 3.578e+06, 0.047
 3.671e+06, 0.056
 3.754e+06, 0.065
 3.829e+06, 0.074
 3.899e+06, 0.083
 3.964e+06, 0.092
 4.024e+06, 0.101
 4.081e+06, 0.110
 4.135e+06, 0.119
 4.186e+06, 0.128
 4.234e+06, 0.137
 4.281e+06, 0.146
 4.326e+06, 0.155
 4.369e+06, 0.164
 4.410e+06, 0.173
 4.450e+06, 0.182
 4.489e+06, 0.191
 4.526e+06, 0.200
*potential
 1.0, 1.2, 0.9, 1.02, 1.1, 1.1
*depvar,delete=3
 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
 3.578e+06, 0.047
 3.671e+06, 0.056
 3.754e+06, 0.065
 3.829e+06, 0.074
 3.899e+06, 0.083
 3.964e+06, 0.092
 4.024e+06, 0.101
 4.081e+06, 0.110
 4.135e+06, 0.119
 4.186e+06, 0.128
 4.234e+06, 0.137
 4.281e+06, 0.146
 4.326e+06, 0.155
 4.369e+06, 0.164
 4.410e+06, 0.173
 4.450e+06, 0.182
 4.489e+06, 0.191
 4.526e+06, 0.200
*potential
 1.0, 1.2, 0.9, 1.02, 1.1, 1.1
**
*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
 1.0, 1.2, 0.9, 1.02, 1.1, 1.1
*depvar,delete=3
 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
 1.0, 1.2, 0.9, 1.02, 1.1, 1.1
**
*material,name=EosUsupMises
*density
 1.9e3,
*eos,type=usup
 4.5631148e3, 0., 0.
*elastic,type=shear
 80.232558e9,
*plastic,hardening=user,properties=8
 175.e6, 0.0
*depvar,delete=3
 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
 4.5631148e3, 0., 0.
*elastic,type=shear
 80.232558e9,
*plastic,hardening=isotropic
 175.e6, 0.0
**
*material,name=EosJwlMises
*density
 1.9e3,
*eos,type=jwl
 7.596,5.206,.053,.35,4.1,1.2,3.63e-2
*detonation point
 5.5,0.5,0.5,0.
 7.5,0.5,0.0,0.
 9.5,0.5,0.0,0.
*elastic,type=shear
 80.232558e9,
*plastic,hardening=user,properties=8
 175.e6
*depvar,delete=3
 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
 7.596,5.206,.053,.35,4.1,1.2,3.63e-2
*detonation point
 5.5,0.5,0.5,0.
 7.5,0.5,0.0,0.
 9.5,0.5,0.0,0.
*elastic,type=shear
 80.232558e9,
*plastic,hardening=isotropic
 175.e6
**
*material,name=eospalphamises
*density
 2.070e3,
*eos, type=usup
 1480, 1.93, 0.880
*eos compaction
 600, 0.049758, 0.0, 6.5e6
*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
 3.578E+06, 0.047
 3.671E+06, 0.056
 3.754E+06, 0.065
 3.829E+06, 0.074
 3.899E+06, 0.083
 3.964E+06, 0.092
 4.024E+06, 0.101
 4.081E+06, 0.110
 4.135E+06, 0.119
 4.186E+06, 0.128
 4.234E+06, 0.137
 4.281E+06, 0.146
 4.326E+06, 0.155
 4.369E+06, 0.164
 4.410E+06, 0.173
 4.450E+06, 0.182
 4.489E+06, 0.191
 4.526E+06, 0.200
*depvar,delete=3
 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
 1480, 1.93, 0.880
*eos compaction
 600, 0.049758, 0.0, 6.5e6
*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
 3.578E+06, 0.047
 3.671E+06, 0.056
 3.754E+06, 0.065
 3.829E+06, 0.074
 3.899E+06, 0.083
 3.964E+06, 0.092
 4.024E+06, 0.101
 4.081E+06, 0.110
 4.135E+06, 0.119
 4.186E+06, 0.128
 4.234E+06, 0.137
 4.281E+06, 0.146
 4.326E+06, 0.155
 4.369E+06, 0.164
 4.410E+06, 0.173
 4.450E+06, 0.182
 4.489E+06, 0.191
 4.526E+06, 0.200
**
*material,name=EosGasMises
*eos,type=ideal gas
 0.2,  
*elastic,type=shear
 80.232558e9,
*specific heat
 0.3,
*density
7800
*plastic,hardening=user,properties=8
175.e6
*depvar,delete=3
 3
 1,YIELD,"Yield Stress"
 2,EQPS,"Equivalent Plastic Strain"
 3,MPSTATUS,"Material Point Status"
*material,name=EosGasMises_ref
*eos,type=ideal gas
 0.2,  
*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
 200.
*depvar,delete=3
 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
 200.
**
*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
 3.578e+06, 0.047
 3.671e+06, 0.056
 3.754e+06, 0.065
 3.829e+06, 0.074
 3.899e+06, 0.083
 3.964e+06, 0.092
 4.024e+06, 0.101
 4.081e+06, 0.110
 4.135e+06, 0.119
 4.186e+06, 0.128
 4.234e+06, 0.137
 4.281e+06, 0.146
 4.326e+06, 0.155
 4.369e+06, 0.164
 4.410e+06, 0.173
 4.450e+06, 0.182
 4.489e+06, 0.191
 4.526e+06, 0.200
*porous metal plasticity,relative density=1.
 1.5,1.,2.25
*depvar,delete=3
 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
 3.578e+06, 0.047
 3.671e+06, 0.056
 3.754e+06, 0.065
 3.829e+06, 0.074
 3.899e+06, 0.083
 3.964e+06, 0.092
 4.024e+06, 0.101
 4.081e+06, 0.110
 4.135e+06, 0.119
 4.186e+06, 0.128
 4.234e+06, 0.137
 4.281e+06, 0.146
 4.326e+06, 0.155
 4.369e+06, 0.164
 4.410e+06, 0.173
 4.450e+06, 0.182
 4.489e+06, 0.191
 4.526e+06, 0.200
*porous metal plasticity,relative density=1.
 1.5,1.,2.25
**
*material,name=FeFp
*density
 1.0e-3
*hyperelastic,mooney-rivlin
 0.749527319, 0.496199379, 1.616263681E-02
*mullins effect
 2., 0.0, 0.5
*plastic,hardening=user,properties=8
 1.5
*depvar,delete=3
 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
 2., 0.0, 0.5
*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