*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
    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.
   11,  0., 2.
   12,  1., 2.
   13,  1., 3.
   14,  0., 3.
**
 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.
 1011,  0., 6.
 1012,  1., 6.
 1013,  1., 7.
 1014,  0., 7.
*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.
  111,  2., 2.
  112,  3., 2.
  113,  3., 3.
  114,  2., 3.
**
 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.
 1111,  2., 6.
 1112,  3., 6.
 1113,  3., 7.
 1114,  2., 7.
*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.
  211,  4., 2.
  212,  5., 2.
  213,  5., 3.
  214,  4., 3.
**
 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.
 1211,  4., 6.
 1212,  5., 6.
 1213,  5., 7.
 1214,  4., 7.
*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.
  311,  6., 2.
  312,  7., 2.
  313,  7., 3.
  314,  6., 3.
**
 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.
 1311,  6., 6.
 1312,  7., 6.
 1313,  7., 7.
 1314,  6., 7.
*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.
  411,  8., 2.
  412,  9., 2.
  413,  9., 3.
  414,  8., 3.
**
 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.
 1411,  8., 6.
 1412,  9., 6.
 1413,  9., 7.
 1414,  8., 7.
**
**  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
**
**  plane strain element
*element, type=cpe4r, elset=pe2d1
 2, 11, 12, 13, 14
*element, type=cpe4r, elset=pe2d2
 102, 111, 112, 113, 114
*element, type=cpe4r, elset=pe2d3
 202, 211, 212, 213, 214
*element, type=cpe4r, elset=pe2d4
 302, 311, 312, 313, 314
*element, type=cpe4r, elset=pe2d5
 402, 411, 412, 413, 414
*element, type=cpe4r, elset=pe2d1001
 1002, 1011, 1012, 1013, 1014
*element, type=cpe4r, elset=pe2d1002
 1102, 1111, 1112, 1113, 1114
*element, type=cpe4r, elset=pe2d1003
 1202, 1211, 1212, 1213, 1214
*element, type=cpe4r, elset=pe2d1004
 1302, 1311, 1312, 1313, 1314
*element, type=cpe4r, elset=pe2d1005
 1402, 1411, 1412, 1413, 1414
**
*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
 1,2,3,4,11,12
 101,102,103,104,111,112
 201,202,203,204,211,212
 301,302,303,304,311,312
 401,402,403,404,411,412
**
 1001,1002,1003,1004,1011,1012
 1101,1102,1103,1104,1111,1112
 1201,1202,1203,1204,1211,1212
 1301,1302,1303,1304,1311,1312
 1401,1402,1403,1404,1411,1412
*nset,nset=left
 1,4,5,8,11,14
 101,104,105,108,111,114
 201,204,205,208,211,214
 301,304,305,308,311,314
 401,404,405,408,411,414
**
 1001,1004,1005,1008,1011,1014
 1101,1104,1105,1108,1111,1114
 1201,1204,1205,1208,1211,1214
 1301,1304,1305,1308,1311,1314
 1401,1404,1405,1408,1411,1414
*nset,nset=top
 5,6,7,8,13,14
 105,106,107,108,113,114
 205,206,207,208,213,214
 305,306,307,308,313,314
 405,406,407,408,413,414
**
 1005,1006,1007,1008,1013,1014
 1105,1106,1107,1108,1113,1114
 1205,1206,1207,1208,1213,1214
 1305,1306,1307,1308,1313,1314
 1405,1406,1407,1408,1413,1414
*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
*elset,elset=eall,generate
 1,401,100
 2,402,100
**
 1001,1401,100
 1002,1402,100
*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
 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
*eos shear,type=viscous
 1.0e4
**
*material,name=mat3
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
 600, 0.049758, 0.0, 6.5e6
*eos shear,type=elastic
 124e6, 0.3
**
*material,name=mat4
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*eos compaction
 600, 0.049758, 0.0, 6.5e6
*eos shear,type=elastic
 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
*eos shear,type=elastic
 124e6, 0.3
*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
*eos shear,type=viscous
 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
*eos shear,type=elastic
 124e6, 0.3
**
*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
*eos shear,type=elastic
 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_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
*eos shear,type=elastic
 124e6, 0.3
*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