*HEADING
 COMPARISON BETWEEN
   1. EOS-TABULAR+ELASTIC,TYPE=SHEAR and ELASTIC;
   2. EOS-TABULAR+ELASTIC,TYPE=SHEAR+MISES and ELASTIC+MISES; and
   3. EOS-TABULAR+ELASTIC,TYPE=SHEAR+JOHNSON COOK and ELASTIC+JOHNSON COOK.
 LOADING: UNAXIAL TENSION
*parameter
 E = 207.E9                  # Young's modulus (Pa)
 nu = 0.29                   # Poisson's ratio
 G = 0.5*E/(1.0+nu)          # Shear modulus (Pa)
 rho0 = 7.890E3              # Reference density (kg/m^3)
 K = E/(3.0*(1.0-2.0*nu))    # Bulk modulus (Pa)
 f2 = 0.0
 evol_0 = 0.0
 f1_0   = 0.0
 evol_1 = -0.1
 f1_1   = -K * evol_1
*NODE,NSET=A1
*NODE,NSET=A2
*NODE,NSET=A3
*NODE,NSET=A4
*NODE,NSET=A5
*NODE,NSET=A6
**
**  3D SOLID ELEMENT
*ELEMENT, TYPE=C3D8R, ELSET=S3D1
*ELEMENT, TYPE=C3D8R, ELSET=S3D2
*ELEMENT, TYPE=C3D8R, ELSET=S3D3
*ELEMENT, TYPE=C3D8R, ELSET=S3D4
*ELEMENT, TYPE=C3D8R, ELSET=S3D5
*ELEMENT, TYPE=C3D8R, ELSET=S3D6
**
**  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=PE2D6
**
*SOLID SECTION,ELSET=S3D1,MATERIAL=A1
*SOLID SECTION,ELSET=PE2D1,MATERIAL=A1
*MATERIAL,NAME=A1
*ELASTIC
 <E>, <nu>
*DENSITY
 <rho0>
*EXPANSION
 0.005,
**
*SOLID SECTION,ELSET=S3D2,MATERIAL=A2
*SOLID SECTION,ELSET=PE2D2,MATERIAL=A2
*MATERIAL,NAME=A2
*eos, type=tabular
 <f1_0>,<f2>,<evol_0>
 <f1_1>,<f2>,<evol_1>
*ELASTIC,TYPE=SHEAR,DEPENDENCIES=1
 <G>,, 10
 <G>,,100
*DENSITY
 <rho0>
**
*SOLID SECTION,ELSET=S3D3,MATERIAL=A3
*SOLID SECTION,ELSET=PE2D3,MATERIAL=A3
*MATERIAL,NAME=A3
*ELASTIC
 <E>, <nu>
*PLASTIC,HARDENING=ISOTROPIC
 175.0000E6,   0.0000
 239.9296E6,   0.0040
 256.0536E6,   0.0080
 267.2830E6,   0.0120
 276.1817E6,   0.0160
 283.6709E6,   0.0200
 297.3914E6,   0.0290
 308.4487E6,   0.0380
*DENSITY
 <rho0>
**
*SOLID SECTION,ELSET=S3D4,MATERIAL=A4
*SOLID SECTION,ELSET=PE2D4,MATERIAL=A4
*MATERIAL,NAME=A4
*eos, type=tabular
 <f1_0>,<f2>,<evol_0>
 <f1_1>,<f2>,<evol_1>
*ELASTIC,TYPE=SHEAR
 <G>,
*PLASTIC,HARDENING=ISOTROPIC
 175.0000E6,   0.0000
 239.9296E6,   0.0040
 256.0536E6,   0.0080
 267.2830E6,   0.0120
 276.1817E6,   0.0160
 283.6709E6,   0.0200
 297.3914E6,   0.0290
 308.4487E6,   0.0380
*DENSITY
 <rho0>
**
*SOLID SECTION,ELSET=S3D5,MATERIAL=A5
*SOLID SECTION,ELSET=PE2D5,MATERIAL=A5
*MATERIAL,NAME=A5
*ELASTIC
 <E>, <nu>
*PLASTIC,HARDENING=JOHNSON COOK
 175.E6, 380.E6, 0.32, 0.55, 1538., 25.
** A & B unit: Pa
*DENSITY
 <rho0>
**
*SOLID SECTION,ELSET=S3D6,MATERIAL=A6
*SOLID SECTION,ELSET=PE2D6,MATERIAL=A6
*MATERIAL,NAME=A6
*eos, type=tabular
 <f1_0>,<f2>,<evol_0>
 <f1_1>,<f2>,<evol_1>
*ELASTIC,TYPE=SHEAR
 <G>
*PLASTIC,HARDENING=JOHNSON COOK
 175.E6, 380.E6, 0.32, 0.55, 1538., 25.
** A & B: Pa
*DENSITY
 <rho0>
**
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=TOP
*NSET,NSET=FACE,GENERATE
*ELSET,ELSET=EALL,GENERATE
*ELSET,ELSET=QA_TEST
EALL,
**
*BOUNDARY
 LEFT, 1, 1
 BOT,  2, 2
 FACE, 3, 3
*AMPLITUDE,NAME=LOAD,DEFINITION=SMOOTH STEP
 0., 0., 1., 1.
*NSET,NSET=NPE2D2
*INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1
NPE2D2, 50.0
**
*RESTART,TIMEMARKS=NO,WRITE,NUMBER=10
**
*STEP
*DYNAMIC,EXPLICIT,ELEMENT BY ELEMENT
 ,1.
**
** Unaxial tension
*BOUNDARY,AMPLITUDE=LOAD
 TOP, 2, 2, 0.1
**
*FILE OUTPUT,NUMBER=4,TIME MARKS=YES
*EL FILE
 PEEQ,
 MISES,
 PRESS,
********************************************
*OUTPUT,FIELD,NUMBER INTERVAL=20,TIMEMARKS=YES
*ELEMENT OUTPUT
S, PEEQ,
*NODE OUTPUT
U
****QA***
*OUTPUT,FIELD,NUMBER INTERVAL=4,TIMEMARKS=YES
*ELEMENT OUTPUT,ELSET=QA_TEST
PEEQ,
*OUTPUT,HISTORY,TIME INTERVAL=0.05
*ENERGY OUTPUT, VAR=PRESELECT
**********************
**
*END STEP