*HEADING
 DATA DECK TO VERIFY JOHNSON COOK PLASTICITY MODEL, damage
 ELEMENTS TESTED: C3D8R, CPE4R, CPS4R, CAX4R
                  S4R, M3D4R
 LOADING: UNIAXIAL TENSION
*NODE,NSET=A1
*NODE,NSET=A2
**
**  3D SOLID ELEMENT
*ELEMENT, TYPE=C3D8R, ELSET=S3D1
*ELEMENT, TYPE=C3D8R, ELSET=S3D2
**
**  PLANE STRAIN ELEMENT
*ELEMENT, TYPE=CPE4R, ELSET=PE2D1
*ELEMENT, TYPE=CPE4R, ELSET=PE2D2
**
**  PLANE STRESS ELEMENT
*ELEMENT, TYPE=CPS4R, ELSET=PS2D1
*ELEMENT, TYPE=CPS4R, ELSET=PS2D2
**
**  AXISYMMETRIC SOLID ELEMENT
*ELEMENT, TYPE=CAX4R, ELSET=CAX1
*ELEMENT, TYPE=CAX4R, ELSET=CAX2
**
**  SHELL ELEMENT
*ELEMENT, TYPE=S4R, ELSET=SHE1
*ELEMENT, TYPE=S4R, ELSET=SHE2
**
**  MEMBRANE
*ELEMENT, TYPE=M3D4R, ELSET=MEMB1
*ELEMENT, TYPE=M3D4R, ELSET=MEMB2
*SHELL SECTION, SECTION INT=GAUSS,ELSET=SHE1,MATERIAL=A1
 0.1, 3
*SOLID SECTION,ELSET=S3D1,MATERIAL=A1
*SOLID SECTION,ELSET=PE2D1,MATERIAL=A1
*SOLID SECTION,ELSET=PS2D1,MATERIAL=A1
*SOLID SECTION,ELSET=CAX1,MATERIAL=A1
*MEMBRANE SECTION,ELSET=MEMB1,MATERIAL=A1
 1.0,
*SHELL SECTION, SECTION INT=GAUSS,ELSET=SHE2,MATERIAL=A2
 0.1, 3
*SOLID SECTION,ELSET=S3D2,MATERIAL=A2
*SOLID SECTION,ELSET=PE2D2,MATERIAL=A2
*SOLID SECTION,ELSET=PS2D2,MATERIAL=A2
*SOLID SECTION,ELSET=CAX2,MATERIAL=A2
*MEMBRANE SECTION,ELSET=MEMB2,MATERIAL=A2
 1.0,
**
** Mises Material Properties (OFHC copper)
*MATERIAL,NAME=A1
*ELASTIC
 124.E9, 0.34
** Young's modulus unit: Pa
*PLASTIC
  90.0000E6,   0.0000
 142.7256E6,   0.0040
 155.3643E6,   0.0080
 164.1189E6,   0.0120
 171.0327E6,   0.0160
 176.8365E6,   0.0200
 187.4375E6,   0.0290
 195.9536E6,   0.0380
** Yield stress unit: Pa
*Damage Initiation, Criterion = Ductile
20.82,  	-3.33,		0.001
20.79706,	-0.333,		0.001
11.87949,	-0.267,		0.001
6.78571,	-0.2,		0.001
3.87612,	-0.133,		0.001
2.21419,	-0.0667,	0.001
1.26496,	0,		0.001
.7229,		0.0667,		0.001
*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL
 1E8
*DENSITY
 8.96E3,
** unit: kg/m^3
**
** JOHNSON-COOK PLASTICITY (OFHC copper)
*MATERIAL,NAME=A2
*ELASTIC
 124.E9, 0.34
** Young's modulus unit: Pa
*PLASTIC,HARDENING=JOHNSON COOK
 90.E6, 292.E6, 0.31, 1.09, 1083., 25.
** A & B unit: Pa
*Damage Initiation, Criterion = Ductile
20.82,  	-3.33,		0.001
20.79706,	-0.333,		0.001
11.87949,	-0.267,		0.001
6.78571,	-0.2,		0.001
3.87612,	-0.133,		0.001
2.21419,	-0.0667,	0.001
1.26496,	0,		0.001
.7229,		0.0667,		0.001
*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL
 1E8
*DENSITY
 8.96E3,
** unit: kg/m^3
**
*NSET, NSET=BOTTOM
*NSET, NSET=TOP
*NSET,NSET=LEFT
*NSET,NSET=FACE
*nset, nset=rot
*ELSET, ELSET=EALL
*ELSET, ELSET=QA_TEST
EALL,
**
*BOUNDARY
 LEFT, 1, 1
 BOTTOM, 2, 2
 FACE, 3, 3
 rot, 4,6
*AMPLITUDE,NAME=LOAD
 0., 0., 1., 1.
**
*STEP,nlgeom=yes,inc=1000
*static
0.01,1.,,0.01
** Uniaxial tension
*BOUNDARY,AMPLITUDE=LOAD
 TOP, 2, 2, 0.1
********************************************
*OUTPUT,FIELD,freq=10
*ELEMENT OUTPUT,ELSET=QA_TEST
PEEQ,TRIAX,DMICRT,SDEG
*OUTPUT,HISTORY
*ENERGY OUTPUT, VAR=PRESELECT
**********************
*END STEP