*heading
 MATERIAL COMBINATIONS:
   MAT1: EOS-USUP+EOS COMPACTION
   MAT2: EOS-USUP+EOS COMPACTION+VISCOSITY;
   MAT3: EOS-USUP+EOS COMPACTION+ELASTIC,TYPE=SHEAR;
   MAT4: EOS-USUP+EOS COMPACTION+ELASTIC,TYPE=SHEAR+MISES;
   MAT5: EOS-USUP+EOS COMPACTION+ELASTIC,TYPE=SHEAR+JOHNSON COOK.
   MAT1_TAB: EOS-TABULAR+EOS COMPACTION
   MAT2_TAB: EOS-TABULAR+EOS COMPACTION+VISCOSITY;
   MAT3_TAB: EOS-TABULAR+EOS COMPACTION+ELASTIC,TYPE=SHEAR;
   MAT4_TAB: EOS-TABULAR+EOS COMPACTION+ELASTIC,TYPE=SHEAR+MISES;
   MAT5_TAB: EOS-TABULAR+EOS COMPACTION+ELASTIC,TYPE=SHEAR+JOHNSON COOK.
 LOADING: HYDROSTATIC PRESSURE / CYCLIC LOADING
*preprint, model=yes
*part, name=part-1
*node,nset=a1
    1,  0., 0., 1.
    2,  1., 0., 1.
    3,  1., 0., 0.
    4,  0., 0., 0.
    5,  0., 1., 1.
    6,  1., 1., 1.
    7,  1., 1., 0.
    8,  0., 1., 0.
**
 1001,  0., 4., 1.
 1002,  1., 4., 1.
 1003,  1., 4., 0.
 1004,  0., 4., 0.
 1005,  0., 5., 1.
 1006,  1., 5., 1.
 1007,  1., 5., 0.
 1008,  0., 5., 0.
*node,nset=a2
  101,  2., 0., 1.
  102,  3., 0., 1.
  103,  3., 0., 0.
  104,  2., 0., 0.
  105,  2., 1., 1.
  106,  3., 1., 1.
  107,  3., 1., 0.
  108,  2., 1., 0.
**
 1101,  2., 4., 1.
 1102,  3., 4., 1.
 1103,  3., 4., 0.
 1104,  2., 4., 0.
 1105,  2., 5., 1.
 1106,  3., 5., 1.
 1107,  3., 5., 0.
 1108,  2., 5., 0.
*node,nset=a3
  201,  4., 0., 1.
  202,  5., 0., 1.
  203,  5., 0., 0.
  204,  4., 0., 0.
  205,  4., 1., 1.
  206,  5., 1., 1.
  207,  5., 1., 0.
  208,  4., 1., 0.
**
 1201,  4., 4., 1.
 1202,  5., 4., 1.
 1203,  5., 4., 0.
 1204,  4., 4., 0.
 1205,  4., 5., 1.
 1206,  5., 5., 1.
 1207,  5., 5., 0.
 1208,  4., 5., 0.
*node,nset=a4
  301,  6., 0., 1.
  302,  7., 0., 1.
  303,  7., 0., 0.
  304,  6., 0., 0.
  305,  6., 1., 1.
  306,  7., 1., 1.
  307,  7., 1., 0.
  308,  6., 1., 0.
**
 1301,  6., 4., 1.
 1302,  7., 4., 1.
 1303,  7., 4., 0.
 1304,  6., 4., 0.
 1305,  6., 5., 1.
 1306,  7., 5., 1.
 1307,  7., 5., 0.
 1308,  6., 5., 0.
*node,nset=a5
  401,  8., 0., 1.
  402,  9., 0., 1.
  403,  9., 0., 0.
  404,  8., 0., 0.
  405,  8., 1., 1.
  406,  9., 1., 1.
  407,  9., 1., 0.
  408,  8., 1., 0.
**
 1401,  8., 4., 1.
 1402,  9., 4., 1.
 1403,  9., 4., 0.
 1404,  8., 4., 0.
 1405,  8., 5., 1.
 1406,  9., 5., 1.
 1407,  9., 5., 0.
 1408,  8., 5., 0.
**
**  3d solid element
*element, type=c3d8r, elset=c3d1
 1, 1, 2, 3, 4, 5, 6, 7, 8
*element, type=c3d8r, elset=c3d2
 101, 101, 102, 103, 104, 105, 106, 107, 108
*element, type=c3d8r, elset=c3d3
 201, 201, 202, 203, 204, 205, 206, 207, 208
*element, type=c3d8r, elset=c3d4
 301, 301, 302, 303, 304, 305, 306, 307, 308
*element, type=c3d8r, elset=c3d5
 401, 401, 402, 403, 404, 405, 406, 407, 408
*element, type=c3d8r, elset=c3d1001
 1001, 1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008
*element, type=c3d8r, elset=c3d1002
 1101, 1101, 1102, 1103, 1104, 1105, 1106, 1107, 1108
*element, type=c3d8r, elset=c3d1003
 1201, 1201, 1202, 1203, 1204, 1205, 1206, 1207, 1208
*element, type=c3d8r, elset=c3d1004
 1301, 1301, 1302, 1303, 1304, 1305, 1306, 1307, 1308
*element, type=c3d8r, elset=c3d1005
 1401, 1401, 1402, 1403, 1404, 1405, 1406, 1407, 1408
**
*solid section,elset=c3d1,material=mat1
*solid section,elset=c3d2,material=mat2
*solid section,elset=c3d3,material=mat3
*solid section,elset=c3d4,material=mat4
*solid section,elset=c3d5,material=mat5
*solid section,elset=c3d1001,material=mat1_tab
*solid section,elset=c3d1002,material=mat2_tab
*solid section,elset=c3d1003,material=mat3_tab
*solid section,elset=c3d1004,material=mat4_tab
*solid section,elset=c3d1005,material=mat5_tab
**
*nset,nset=bot
 1,2,3,4
 101,102,103,104
 201,202,203,204
 301,302,303,304
 401,402,403,404
**
 1001,1002,1003,1004
 1101,1102,1103,1104
 1201,1202,1203,1204
 1301,1302,1303,1304
 1401,1402,1403,1404
*nset,nset=left
 1,4,5,8
 101,104,105,108
 201,204,205,208
 301,304,305,308
 401,404,405,408
**
 1001,1004,1005,1008
 1101,1104,1105,1108
 1201,1204,1205,1208
 1301,1304,1305,1308
 1401,1404,1405,1408
*nset,nset=right
 2,3,6,7
 102,103,106,107
 202,203,206,207
 302,303,306,307
 402,403,406,407
*nset,nset=top
 5,6,7,8
 105,106,107,108
 205,206,207,208
 305,306,307,308
 405,406,407,408
**
 1005,1006,1007,1008
 1105,1106,1107,1108
 1205,1206,1207,1208
 1305,1306,1307,1308
 1405,1406,1407,1408
*nset,nset=face,generate
 1,401,100
 2,402,100
 5,405,100
 6,406,100
**
 1001,1401,100
 1002,1402,100
 1005,1405,100
 1006,1406,100
*nset,nset=back,generate
 3,403,100
 4,404,100
 7,407,100
 8,408,100
**
 1003,1403,100
 1004,1404,100
 1007,1407,100
 1008,1408,100
*elset,elset=eall,generate
 1,401,100
 2,402,100
**
 1001,1401,100
 1002,1402,100
*elset,elset=qa_test
 eall,
*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
 600, 0.049758, 0.0, 6.5e6
**
*material,name=mat2
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
 600, 0.049758, 0.0, 6.5e6
*viscosity
 1.0e4
**
*material,name=mat3
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
 600, 0.049758, 0.0, 6.5e6
*elastic, type=shear
 124e6
**
*material,name=mat4
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*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=mat5
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
 600, 0.049758, 0.0, 6.5e6
*elastic, type=shear
 124e6
*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
 600, 0.049758, 0.0, 6.5e6
**
*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
 600, 0.049758, 0.0, 6.5e6
*viscosity
 1.0e4
**
*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
 600, 0.049758, 0.0, 6.5e6
*elastic, type=shear
 124e6
**
*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
 600, 0.049758, 0.0, 6.5e6
*elastic, type=shear
 124e6
*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=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
 600, 0.049758, 0.0, 6.5e6
*elastic, type=shear
 124e6
*plastic,hardening=johnson cook
 1.97e6, 4.28e6, 0.32, 0.55, 1538., 25.
**
*boundary
 assembly-1.part-1-1.left, 1, 1
 assembly-1.part-1-1.bot,  2, 2
 assembly-1.part-1-1.back, 3, 3
**
*restart,timemarks=no,write,number=10
**
*amplitude,name=load,definition=smooth step
 0., 0., 0.1, 1.
**
*amplitude,name=unload,definition=smooth step
 0., 1., 0.1, 0.
**
*initial conditions, type=specific energy
 assembly-1.part-1-1.eall, 0.0
*initial conditions, type=porosity
 assembly-1.part-1-1.eall, 0.049758
*step
*dynamic,explicit
 ,0.1
** hydrostatic loading
*dload, amplitude = load
 assembly-1.part-1-1.eall, P2, 2.0e6
 assembly-1.part-1-1.eall, P3, 2.0e6
 assembly-1.part-1-1.eall, P4, 2.0e6
**
*output,field
*element output,elset=assembly-1.part-1-1.qa_test
 peeq, palph, palphmin
*element output,elset=assembly-1.part-1-1.eall
 s
*node output
 u
*output,history
*element output,elset=assembly-1.part-1-1.qa_test
 peeq, palph, palphmin
*element output,elset=assembly-1.part-1-1.eall
 s22
*energy output, var=preselect
**
*end step
**
*step
*dynamic,explicit
 ,0.1
** unloading
*dload, amplitude = unload
 assembly-1.part-1-1.eall, P2, 2.0e6
 assembly-1.part-1-1.eall, P3, 2.0e6
 assembly-1.part-1-1.eall, P4, 2.0e6
*end step
**
*step
*dynamic,explicit
 ,0.1
** loading
*dload, amplitude = load
 assembly-1.part-1-1.eall, P2, 4.0e6
 assembly-1.part-1-1.eall, P3, 4.0e6
 assembly-1.part-1-1.eall, P4, 4.0e6
*end step
**
*step
*dynamic,explicit
 ,0.1
** unloading
*dload, amplitude = unload
 assembly-1.part-1-1.eall, P2, 4.0e6
 assembly-1.part-1-1.eall, P3, 4.0e6
 assembly-1.part-1-1.eall, P4, 4.0e6
*end step
**
*step
*dynamic,explicit
 ,0.1
** loading
*dload, amplitude = load
 assembly-1.part-1-1.eall, P2, 7.0e6
 assembly-1.part-1-1.eall, P3, 7.0e6
 assembly-1.part-1-1.eall, P4, 7.0e6
*end step
**