*heading
 MATERIAL COMBINATIONS:
   MAT1: EOS-USUP+EOS COMPACTION
   MAT2: EOS-USUP+EOS COMPACTION+EOS SHEAR (VISCOUS);
   MAT3: EOS-USUP+EOS COMPACTION+EOS SHEAR (ELASTIC);
   MAT4: EOS-USUP+EOS COMPACTION+EOS SHEAR+MISES;
   MAT5: EOS-USUP+EOS COMPACTION+EOS SHEAR+JOHNSON COOK.
   MAT1_TAB: EOS-TABULAR+EOS COMPACTION
   MAT2_TAB: EOS-TABULAR+EOS COMPACTION+EOS SHEAR (VISCOUS);
   MAT3_TAB: EOS-TABULAR+EOS COMPACTION+EOS SHEAR (ELASTIC);
   MAT4_TAB: EOS-TABULAR+EOS COMPACTION+EOS SHEAR+MISES;
   MAT5_TAB: EOS-TABULAR+EOS COMPACTION+EOS SHEAR+JOHNSON COOK.
 LOADING: SIMPLE SHEAR. Initial conditions in PEEQ.
*preprint, model=yes
*part, name=part-1
*node,nset=a1
**
*node,nset=a2
**
*node,nset=a3
**
*node,nset=a4
**
*node,nset=a5
**
**
**  3d solid element
*element, type=c3d8r, elset=c3d1
*element, type=c3d8r, elset=c3d2
*element, type=c3d8r, elset=c3d3
*element, type=c3d8r, elset=c3d4
*element, type=c3d8r, elset=c3d5
*element, type=c3d8r, elset=c3d1001
*element, type=c3d8r, elset=c3d1002
*element, type=c3d8r, elset=c3d1003
*element, type=c3d8r, elset=c3d1004
*element, type=c3d8r, elset=c3d1005
**
**  plane strain element
*element, type=cpe4r, elset=pe2d1
*element, type=cpe4r, elset=pe2d2
*element, type=cpe4r, elset=pe2d3
*element, type=cpe4r, elset=pe2d4
*element, type=cpe4r, elset=pe2d5
*element, type=cpe4r, elset=pe2d1001
*element, type=cpe4r, elset=pe2d1002
*element, type=cpe4r, elset=pe2d1003
*element, type=cpe4r, elset=pe2d1004
*element, type=cpe4r, elset=pe2d1005
**
*solid section,elset=c3d1,material=mat1
*solid section,elset=pe2d1,material=mat1
**
*solid section,elset=c3d2,material=mat2
*solid section,elset=pe2d2,material=mat2
**
*solid section,elset=c3d3,material=mat3
*solid section,elset=pe2d3,material=mat3
**
*solid section,elset=c3d4,material=mat4
*solid section,elset=pe2d4,material=mat4
**
*solid section,elset=c3d5,material=mat5
*solid section,elset=pe2d5,material=mat5
**
*solid section,elset=c3d1001,material=mat1_tab
*solid section,elset=pe2d1001,material=mat1_tab
**
*solid section,elset=c3d1002,material=mat2_tab
*solid section,elset=pe2d1002,material=mat2_tab
**
*solid section,elset=c3d1003,material=mat3_tab
*solid section,elset=pe2d1003,material=mat3_tab
**
*solid section,elset=c3d1004,material=mat4_tab
*solid section,elset=pe2d1004,material=mat4_tab
**
*solid section,elset=c3d1005,material=mat5_tab
*solid section,elset=pe2d1005,material=mat5_tab
**
*nset,nset=bot
**
*nset,nset=left
**
*nset,nset=top
**
*nset,nset=face,generate
**
*elset,elset=eall,generate
**
*elset,elset=qa_test
 eall,
**
*elset, elset=eplastic
 c3d4, pe2d4, c3d5, pe2d5
 c3d1004, pe2d1004, c3d1005, pe2d1005
*end part
**
*assembly, name=assembly-1
*instance, name=part-1-1, part=part-1
*end instance
*end assembly
**
*parameter
 c0 = 1480
 s = 1.93
 Gamma0 = 0.880
 rho0 = 2070.0
 f2 = Gamma0
 evol1 = 0.0
 evol2 = -0.1
 evol3 = -0.2
 evol4 = -0.3
 evol5 = -0.4
 evol6 = -0.5
 eta1 = 1.0 - exp(evol1)
 eta2 = 1.0 - exp(evol2)
 eta3 = 1.0 - exp(evol3)
 eta4 = 1.0 - exp(evol4)
 eta5 = 1.0 - exp(evol5)
 eta6 = 1.0 - exp(evol6)
 f11 = rho0*c0*c0*eta1*(1.0-Gamma0*eta1/2)/((1.0-s*eta1)*(1.0-s*eta1))
 f12 = rho0*c0*c0*eta2*(1.0-Gamma0*eta2/2)/((1.0-s*eta2)*(1.0-s*eta2))
 f13 = rho0*c0*c0*eta3*(1.0-Gamma0*eta3/2)/((1.0-s*eta3)*(1.0-s*eta3))
 f14 = rho0*c0*c0*eta4*(1.0-Gamma0*eta4/2)/((1.0-s*eta4)*(1.0-s*eta4))
 f15 = rho0*c0*c0*eta5*(1.0-Gamma0*eta5/2)/((1.0-s*eta5)*(1.0-s*eta5))
 f16 = rho0*c0*c0*eta6*(1.0-Gamma0*eta6/2)/((1.0-s*eta6)*(1.0-s*eta6))
**
*material,name=mat1
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
**
*material,name=mat2
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
*eos shear,type=viscous
**
*material,name=mat3
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
*eos shear,type=elastic
**
*material,name=mat4
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
*eos shear,type=elastic
*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=mat5
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
*eos shear,type=elastic
*plastic,hardening=johnson cook
 1.97e6, 4.28e6, 0.32, 0.55, 1538., 25.
**
*material,name=mat1_tab
*density
 <rho0>
*eos, type=tabular
 <f11>,<f2>,<evol1>
 <f12>,<f2>,<evol2>
 <f13>,<f2>,<evol3>
 <f14>,<f2>,<evol4>
 <f15>,<f2>,<evol5>
 <f16>,<f2>,<evol6>
*eos compaction
**
*material,name=mat2_tab
*density
 <rho0>
*eos, type=tabular
 <f11>,<f2>,<evol1>
 <f12>,<f2>,<evol2>
 <f13>,<f2>,<evol3>
 <f14>,<f2>,<evol4>
 <f15>,<f2>,<evol5>
 <f16>,<f2>,<evol6>
*eos compaction
*eos shear,type=viscous
**
*material,name=mat3_tab
*density
 <rho0>
*eos, type=tabular
 <f11>,<f2>,<evol1>
 <f12>,<f2>,<evol2>
 <f13>,<f2>,<evol3>
 <f14>,<f2>,<evol4>
 <f15>,<f2>,<evol5>
 <f16>,<f2>,<evol6>
*eos compaction
*eos shear,type=elastic
**
*material,name=mat4_tab
*density
 <rho0>
*eos, type=tabular
 <f11>,<f2>,<evol1>
 <f12>,<f2>,<evol2>
 <f13>,<f2>,<evol3>
 <f14>,<f2>,<evol4>
 <f15>,<f2>,<evol5>
 <f16>,<f2>,<evol6>
*eos compaction
*eos shear,type=elastic
*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=mat5_tab
*density
 <rho0>
*eos, type=tabular
 <f11>,<f2>,<evol1>
 <f12>,<f2>,<evol2>
 <f13>,<f2>,<evol3>
 <f14>,<f2>,<evol4>
 <f15>,<f2>,<evol5>
 <f16>,<f2>,<evol6>
*eos compaction
*eos shear,type=elastic
*plastic,hardening=johnson cook
 1.97e6, 4.28e6, 0.32, 0.55, 1538., 25.
**
*boundary
 assembly-1.part-1-1.bot,  1, 3
 assembly-1.part-1-1.top,  2, 3
**
*restart,timemarks=no,write,number=10
*initial conditions, type=hardening
 assembly-1.part-1-1.eplastic, 0.02
**
*step
*dynamic,explicit
 ,1.0
** shear
*boundary, type=velocity
 assembly-1.part-1-1.top, 1, 1, 0.1
**
*output,field
*element output,elset=assembly-1.part-1-1.qa_test
 s, peeq, palph, palphmin
*node output
 u
*output,history
*element output,elset=assembly-1.part-1-1.qa_test
 s12, peeq, palph, palphmin
*energy output, var=preselect
**
*end step