*heading
 MATERIAL COMBINATIONS: 
   MAT1: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=CROSS    
   MAT2: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=USER
*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.
*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.
**
**  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
**
**  plane strain element
*element, type=cpe4r, elset=pe2d1
 2, 11, 12, 13, 14
*element, type=cpe4r, elset=pe2d2
 102, 111, 112, 113, 114
**
*solid section,elset=c3d1,material=mat1
*solid section,elset=pe2d1,material=mat1
**
*solid section,elset=c3d2,material=mat2
*solid section,elset=pe2d2,material=mat2
**
*nset,nset=bot
 1,2,3,4,11,12
 101,102,103,104,111,112
*nset,nset=left
 1,4,5,8,11,14
 101,104,105,108,111,114
*nset,nset=top
 5,6,7,8,13,14
 105,106,107,108,113,114
*nset,nset=face,generate
 1,101,100
 2,102,100
*elset,elset=eall,generate
 1,101,100
 2,102,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
*parameter
 eta_0 = 1.0
 eta_inf = 0.0
 n = 0.392
 lamda = 0.11
**
*material,name=mat1
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*viscosity, definition=Cross
 <eta_0>, <eta_inf>, <lamda>, <n>
**
*material,name=mat2
*density
 <rho0>
*eos, type=usup
 <c0>, <s>, <Gamma0>
*viscosity, definition=User, Properties=3
 <eta_0>, <lamda>, <n>
**
*boundary
 assembly-1.part-1-1.bot,  1, 3
 assembly-1.part-1-1.top,  2, 3
**
*restart,timemarks=no,write,number=10
**
*amplitude, name=ramp
 0.0, 0.0, 1.0, 15.0
**
*step
*dynamic,explicit
 ,1.0
** shear
*boundary, type=velocity, amplitude=ramp
 assembly-1.part-1-1.top, 1, 1, 1.0
**
*output,field
*element output,elset=assembly-1.part-1-1.qa_test
 s
*node output
 u
*output,history
*element output,elset=assembly-1.part-1-1.qa_test
 s12, le12
*energy output, var=preselect
**
*end step