*HEADING
 DATA DECK TO VERIFY JOHNSON COOK PLASTICITY MODEL, Hc damage
 ELEMENTS TESTED: C3D8R, CPE4R, CPS4R, CAX4R  
                  S4R, M3D4R
 Input file name: johnsoncook_dmg_hc_s.inp
 LOADING: UNIAXIAL TENSION
*NODE,NSET=A1
  11,  6.0, 0.0, 1.0
  12,  7.0, 0.0, 1.0
  13,  7.0, 0.0, 0.0
  14,  6.0, 0.0, 0.0
  15,  6.0, 1.0, 1.0
  16,  7.0, 1.0, 1.0
  17,  7.0, 1.0, 0.0
  18,  6.0, 1.0, 0.0
 101,  8.0, 0.0
 102,  9.0, 0.0
 103,  9.0, 1.0
 104,  8.0, 1.0
 201, 10.0, 0.0
 202, 11.0, 0.0
 203, 11.0, 1.0
 204, 10.0, 1.0
 301, 12.0, 0.0
 302, 13.0, 0.0
 303, 13.0, 1.0
 304, 12.0, 1.0
 401, 14.0, 0.0
 402, 15.0, 0.0
 403, 15.0, 1.0
 404, 14.0, 1.0
 501, 16.0, 0.0
 502, 17.0, 0.0
 503, 17.0, 1.0
 504, 16.0, 1.0
*NODE,NSET=A2
 1011,  6.0, 2.0, 1.0
 1012,  7.0, 2.0, 1.0
 1013,  7.0, 2.0, 0.0
 1014,  6.0, 2.0, 0.0
 1015,  6.0, 3.0, 1.0
 1016,  7.0, 3.0, 1.0
 1017,  7.0, 3.0, 0.0
 1018,  6.0, 3.0, 0.0
 1101,  8.0, 2.0
 1102,  9.0, 2.0
 1103,  9.0, 3.0
 1104,  8.0, 3.0
 1201, 10.0, 2.0
 1202, 11.0, 2.0
 1203, 11.0, 3.0
 1204, 10.0, 3.0
 1301, 12.0, 2.0
 1302, 13.0, 2.0
 1303, 13.0, 3.0
 1304, 12.0, 3.0
 1401, 14.0, 2.0
 1402, 15.0, 2.0
 1403, 15.0, 3.0
 1404, 14.0, 3.0
 1501, 16.0, 2.0
 1502, 17.0, 2.0
 1503, 17.0, 3.0
 1504, 16.0, 3.0
**
**  3D SOLID ELEMENT
*ELEMENT, TYPE=C3D8R, ELSET=S3D1
 11, 11, 12, 13, 14, 15, 16, 17, 18
*ELEMENT, TYPE=C3D8R, ELSET=S3D2
 1011,1011,1012,1013,1014,1015,1016,1017,1018
**
**  PLANE STRAIN ELEMENT
*ELEMENT, TYPE=CPE4R, ELSET=PE2D1
 101, 101, 102, 103, 104
*ELEMENT, TYPE=CPE4R, ELSET=PE2D2
 1101, 1101, 1102, 1103, 1104
**
**  PLANE STRESS ELEMENT
*ELEMENT, TYPE=CPS4R, ELSET=PS2D1
 201, 201, 202, 203, 204
*ELEMENT, TYPE=CPS4R, ELSET=PS2D2
 1201, 1201, 1202, 1203, 1204
**
**  AXISYMMETRIC SOLID ELEMENT
*ELEMENT, TYPE=CAX4R, ELSET=CAX1
 301, 301, 302, 303, 304
*ELEMENT, TYPE=CAX4R, ELSET=CAX2
 1301, 1301, 1302, 1303, 1304
**
**  SHELL ELEMENT
*ELEMENT, TYPE=S4R, ELSET=SHE1
 401, 401, 402, 403, 404
*ELEMENT, TYPE=S4R, ELSET=SHE2
 1401, 1401, 1402, 1403, 1404
**
**  MEMBRANE
*ELEMENT, TYPE=M3D4R, ELSET=MEMB1
 501, 501, 502, 503, 504
*ELEMENT, TYPE=M3D4R, ELSET=MEMB2
 1501, 1501, 1502, 1503, 1504
*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
 203.1705E6,   0.0470
 209.4871E6,   0.0560
 215.1371E6,   0.0650
 220.2701E6,   0.0740
 224.9886E6,   0.0830
 229.3661E6,   0.0920
 233.4572E6,   0.1010
 237.3040E6,   0.1100
 240.9393E6,   0.1190
 244.3896E6,   0.1280
 247.6762E6,   0.1370
 250.8171E6,   0.1460
 253.8271E6,   0.1550
 256.7187E6,   0.1640
 259.5029E6,   0.1730
 262.1888E6,   0.1820
 264.7846E6,   0.1910
 267.2973E6,   0.2000
** Yield stress unit: Pa
*Damage Initiation, Criterion = HC 
1.0,0.05,0.01,0.1
*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 = HC  
1.0,0.05,0.01,0.1
*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=EXPONENTIAL
 1E8
*DENSITY
 8.96E3,
** unit: kg/m^3
**
*NSET, NSET=BOTTOM
 11,12,13,14,101,102
 201,202,301,302,401,402,501,502,
 1011,1012,1013,1014,1101,1102
 1201,1202,1301,1302,1401,1402,1501,1502
*NSET, NSET=TOP
 15,16,17,18,103,104
 203,204,303,304,403,404,503,504
 1015,1016,1017,1018,1103,1104
 1203,1204,1303,1304,1403,1404,1503,1504
*NSET,NSET=LEFT
 11,14,15,18,
 101,104,201,204,
 301,304,401,404,
 501,504
 1011,1014,1015,1018,
 1101,1104,1201,1204,
 1301,1304,1401,1404,
 1501,1504
*NSET,NSET=FACE
 11,12,15,16,
 401,402,403,404,501,502,503,504,
 1011,1012,1015,1016,
 1401,1402,1403,1404,1501,1502,1503,1504
*nset, nset=rot
403,1403
*ELSET, ELSET=EALL
 11,101,201,301,401,501
 1011,1101,1201,1301,1401,1501
*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