*HEADING
 COMPARISON BETWEEN EOS+EOS SHEAR+DRUCKER PRAGER and ELASTIC+DRUCKER PRAGER;
 LOADING: UNIAXIAL TENSION
*preprint, model=yes
*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.
*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.
*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.
*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.
*NODE,NSET=A6
 501, 10., 0., 1.
 502, 11., 0., 1.
 503, 11., 0., 0.
 504, 10., 0., 0.
 505, 10., 1., 1.
 506, 11., 1., 1.
 507, 11., 1., 0.
 508, 10., 1., 0.
 511, 10., 2.
 512, 11., 2.
 513, 11., 3.
 514, 10., 3.
**
**  3D SOLID ELEMENT
*ELEMENT, TYPE=C3D8R, ELSET=S3D1
 1, 1, 2, 3, 4, 5, 6, 7, 8
*ELEMENT, TYPE=C3D8R, ELSET=S3D2
 101, 101, 102, 103, 104, 105, 106, 107, 108
*ELEMENT, TYPE=C3D8R, ELSET=S3D3
 201, 201, 202, 203, 204, 205, 206, 207, 208
*ELEMENT, TYPE=C3D8R, ELSET=S3D4
 301, 301, 302, 303, 304, 305, 306, 307, 308
*ELEMENT, TYPE=C3D8R, ELSET=S3D5
 401, 401, 402, 403, 404, 405, 406, 407, 408
*ELEMENT, TYPE=C3D8R, ELSET=S3D6
 501, 501, 502, 503, 504, 505, 506, 507, 508
**
**  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=PE2D6
 502, 511, 512, 513, 514
**
*SOLID SECTION,ELSET=S3D1,MATERIAL=A1
*SOLID SECTION,ELSET=PE2D1,MATERIAL=A1
*MATERIAL,NAME=A1
** G=0.5*E/(1+anu)=80.232558E9
** K=0.33333*E/(1-2*anu)=164.28571E9
*ELASTIC
 207.E9, 0.29
** Young's modulus unit: Pa
*DRUCKER PRAGER, SHEAR CRITERION=LINEAR
 40.0,
*DRUCKER PRAGER HARDENING, TYPE=COMPRESSION
 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
 317.8417E6,   0.0470
 326.0789E6,   0.0560
 333.4586E6,   0.0650
 340.1725E6,   0.0740
 346.3518E6,   0.0830
 352.0907E6,   0.0920
 357.4595E6,   0.1010
 362.5121E6,   0.1100
 367.2909E6,   0.1190
 371.8298E6,   0.1280
 376.1565E6,   0.1370
 380.2941E6,   0.1460
 384.2617E6,   0.1550
 388.0756E6,   0.1640
 391.7496E6,   0.1730
 395.2959E6,   0.1820
 398.7249E6,   0.1910
 402.0457E6,   0.2000
*DENSITY
 7.890E3,
** unit: kg/m^3
**
*SOLID SECTION,ELSET=S3D2,MATERIAL=A2
*SOLID SECTION,ELSET=PE2D2,MATERIAL=A2
*MATERIAL,NAME=A2
** G=0.5*E/(1+anu)=80.232558E9
** K=0.33333*E/(1-2*anu)=164.28571E9
*EOS,TYPE=USUP
 4.5631148E3, 0., 0.
*ELASTIC, TYPE=SHEAR
 80.232558E9,
** Shear modulus unit: Pa
*DRUCKER PRAGER, SHEAR CRITERION=LINEAR
 40.0,
*DRUCKER PRAGER HARDENING, TYPE=COMPRESSION
 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
 317.8417E6,   0.0470
 326.0789E6,   0.0560
 333.4586E6,   0.0650
 340.1725E6,   0.0740
 346.3518E6,   0.0830
 352.0907E6,   0.0920
 357.4595E6,   0.1010
 362.5121E6,   0.1100
 367.2909E6,   0.1190
 371.8298E6,   0.1280
 376.1565E6,   0.1370
 380.2941E6,   0.1460
 384.2617E6,   0.1550
 388.0756E6,   0.1640
 391.7496E6,   0.1730
 395.2959E6,   0.1820
 398.7249E6,   0.1910
 402.0457E6,   0.2000
*DENSITY
 7.890E3,
** unit: kg/m^3
**
*SOLID SECTION,ELSET=S3D3,MATERIAL=A3
*SOLID SECTION,ELSET=PE2D3,MATERIAL=A3
*MATERIAL,NAME=A3
** G=0.5*E/(1+anu)=80.232558E9
** K=0.33333*E/(1-2*anu)=164.28571E9
*ELASTIC
 207.E9, 0.29
** Young's modulus unit: Pa
*DRUCKER PRAGER, SHEAR CRITERION=EXPONENT 
 1.19175,1.0,
*DRUCKER PRAGER HARDENING, TYPE=COMPRESSION
 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
 317.8417E6,   0.0470
 326.0789E6,   0.0560
 333.4586E6,   0.0650
 340.1725E6,   0.0740
 346.3518E6,   0.0830
 352.0907E6,   0.0920
 357.4595E6,   0.1010
 362.5121E6,   0.1100
 367.2909E6,   0.1190
 371.8298E6,   0.1280
 376.1565E6,   0.1370
 380.2941E6,   0.1460
 384.2617E6,   0.1550
 388.0756E6,   0.1640
 391.7496E6,   0.1730
 395.2959E6,   0.1820
 398.7249E6,   0.1910
 402.0457E6,   0.2000
*DENSITY
 7.890E3,
** unit: kg/m^3
**
*SOLID SECTION,ELSET=S3D4,MATERIAL=A4
*SOLID SECTION,ELSET=PE2D4,MATERIAL=A4
*MATERIAL,NAME=A4
** G=0.5*E/(1+anu)=80.232558E9
** K=0.33333*E/(1-2*anu)=164.28571E9
*EOS,TYPE=USUP
 4.5631148E3, 0., 0.
*ELASTIC, TYPE=SHEAR
 80.232558E9,
** Shear modulus unit: Pa
*DRUCKER PRAGER, SHEAR CRITERION=EXPONENT 
 1.19175,1.0
*DRUCKER PRAGER HARDENING, TYPE=COMPRESSION
 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
 317.8417E6,   0.0470
 326.0789E6,   0.0560
 333.4586E6,   0.0650
 340.1725E6,   0.0740
 346.3518E6,   0.0830
 352.0907E6,   0.0920
 357.4595E6,   0.1010
 362.5121E6,   0.1100
 367.2909E6,   0.1190
 371.8298E6,   0.1280
 376.1565E6,   0.1370
 380.2941E6,   0.1460
 384.2617E6,   0.1550
 388.0756E6,   0.1640
 391.7496E6,   0.1730
 395.2959E6,   0.1820
 398.7249E6,   0.1910
 402.0457E6,   0.2000
*DENSITY
 7.890E3,
** unit: kg/m^3
**
*SOLID SECTION,ELSET=S3D5,MATERIAL=A5
*SOLID SECTION,ELSET=PE2D5,MATERIAL=A5
*MATERIAL,NAME=A5
** G=0.5*E/(1+anu)=80.232558E9
** K=0.33333*E/(1-2*anu)=164.28571E9
*ELASTIC
 207.E9, 0.29
*DRUCKER PRAGER, SHEAR CRITERION=HYPERBOLIC
 20.0,100.0e6
*DRUCKER PRAGER HARDENING, TYPE=COMPRESSION
 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
 317.8417E6,   0.0470
 326.0789E6,   0.0560
 333.4586E6,   0.0650
 340.1725E6,   0.0740
 346.3518E6,   0.0830
 352.0907E6,   0.0920
 357.4595E6,   0.1010
 362.5121E6,   0.1100
 367.2909E6,   0.1190
 371.8298E6,   0.1280
 376.1565E6,   0.1370
 380.2941E6,   0.1460
 384.2617E6,   0.1550
 388.0756E6,   0.1640
 391.7496E6,   0.1730
 395.2959E6,   0.1820
 398.7249E6,   0.1910
 402.0457E6,   0.2000
*DENSITY
 7.890E3,
** unit: kg/m^3
**
*SOLID SECTION,ELSET=S3D6,MATERIAL=A6
*SOLID SECTION,ELSET=PE2D6,MATERIAL=A6
*MATERIAL,NAME=A6
** G=0.5*E/(1+anu)=80.232558E9
** K=0.33333*E/(1-2*anu)=164.28571E9
*EOS,TYPE=USUP
 4.5631148E3, 0., 0.
*ELASTIC, TYPE=SHEAR
 80.232558E9,
** Shear modulus unit: Pa
*DRUCKER PRAGER, SHEAR CRITERION=HYPERBOLIC
 20.0,100.0e6
*DRUCKER PRAGER HARDENING, TYPE=COMPRESSION
 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
 317.8417E6,   0.0470
 326.0789E6,   0.0560
 333.4586E6,   0.0650
 340.1725E6,   0.0740
 346.3518E6,   0.0830
 352.0907E6,   0.0920
 357.4595E6,   0.1010
 362.5121E6,   0.1100
 367.2909E6,   0.1190
 371.8298E6,   0.1280
 376.1565E6,   0.1370
 380.2941E6,   0.1460
 384.2617E6,   0.1550
 388.0756E6,   0.1640
 391.7496E6,   0.1730
 395.2959E6,   0.1820
 398.7249E6,   0.1910
 402.0457E6,   0.2000
*DENSITY
 7.890E3,
** unit: kg/m^3
**
*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
 501,502,503,504,511,512
*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
 501,504,505,508,511,514
*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
 505,506,507,508,513,514
*NSET,NSET=FACE,GENERATE
 1,501,100
 2,502,100
 5,505,100
 6,506,100
*ELSET,ELSET=EALL,GENERATE
 1,501,100
 2,502,100
*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.
**
*STEP
*DYNAMIC,EXPLICIT,ELEMENT BY ELEMENT
 ,1.
**
** Uniaxial 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=4,TIMEMARKS=YES
*NODE OUTPUT
U
*ELEMENT OUTPUT,ELSET=QA_TEST
PEEQ,S
*OUTPUT,HISTORY,TIME INTERVAL=0.25
*ENERGY OUTPUT, VAR=PRESELECT
**********************
*OUTPUT,HISTORY,FREQUENCY=1
*ELEMENT OUTPUT,ELSET=EALL
PEEQ,S,MISES,PRESS,PE
**
*END STEP