*HEADING
DRUCKER PRAGER PLASTICITY WITH DIFFERENT PARAMETERS OF
BETA, PSI AND K
*NODE, NSET=QA_TEST_ALLNODES
101,0.0,0.,0.
102,1.0,0.,0.
103,1.0,1.,0.
104,0.0,1.,0.
105,0.0,0.,1.
106,1.0,0.,1.
107,1.0,1.,1.
108,0.0,1.,1.
**
*NODE,NSET=FIX3
201,1.5,0.,0.
202,2.5,0.,0.
203,2.5,1.,0.
204,1.5,1.,0.
205,1.5,0.,1.
206,2.5,0.,1.
207,2.5,1.,1.
208,1.5,1.,1.
*NODE, NSET=QA_TEST_ALLNODES
301,3.0,0.
302,4.0,0.
303,4.0,1.
304,3.0,1.
**
401,4.5,0.
402,5.5,0.
403,5.5,1.
404,4.5,1.
**
501,0.0,1.5,0.
502,1.0,1.5,0.
503,1.0,2.5,0.
504,0.0,2.5,0.
505,0.0,1.5,1.
506,1.0,1.5,1.
507,1.0,2.5,1.
508,0.0,2.5,1.
**
*NODE,NSET=FIX3
601,1.5,1.5,0.
602,2.5,1.5,0.
603,2.5,2.5,0.
604,1.5,2.5,0.
605,1.5,1.5,1.
606,2.5,1.5,1.
607,2.5,2.5,1.
608,1.5,2.5,1.
*NODE, NSET=QA_TEST_ALLNODES
701,3.0,1.5
702,4.0,1.5
703,4.0,2.5
704,3.0,2.5
**
801,4.5,1.5
802,5.5,1.5
803,5.5,2.5
804,4.5,2.5
*NSET,NSET=BOTL
101,102,105,106, 201,202,205,206, 301,302,401,402
*NSET,NSET=BOTU
501,502,505,506, 601,602,605,606, 701,702,801,802
*NSET,NSET=TOPL
103,104,107,108, 203,204,207,208, 303,304,403,404
*NSET,NSET=TOPU
503,504,507,508, 603,604,607,608, 703,704,803,804
*ELEMENT,TYPE=C3D8R,ELSET=C3D_1
1, 101,102,103,104,105,106,107,108
*ELEMENT,TYPE=C3D8R,ELSET=C3D_2
2, 201,202,203,204,205,206,207,208
*ELEMENT,TYPE=C3D8R,ELSET=C3D_3
5, 501,502,503,504,505,506,507,508
*ELEMENT,TYPE=C3D8R,ELSET=C3D_4
6, 601,602,603,604,605,606,607,608
*ELEMENT,TYPE=CPE4R,ELSET=GPE_1
3, 301,302,303,304
*ELEMENT,TYPE=CPE4R,ELSET=GPE_2
7, 701,702,703,704
*ELEMENT,TYPE=CPS4R,ELSET=PS_1
4, 401,402,403,404
*ELEMENT,TYPE=CPS4R,ELSET=PS_2
8, 801,802,803,804
*NODE, NSET=QA_TEST_ALLNODES
1101,0.0,-3.,0.
1102,1.0,-3.,0.
1103,1.0,-2.,0.
1104,0.0,-2.,0.
1105,0.0,-3.,1.
1106,1.0,-3.,1.
1107,1.0,-2.,1.
1108,0.0,-2.,1.
**
*NODE,NSET=FIX3
1201,1.5,-3.,0.
1202,2.5,-3.,0.
1203,2.5,-2.,0.
1204,1.5,-2.,0.
1205,1.5,-3.,1.
1206,2.5,-3,1.
1207,2.5,-2.,1.
1208,1.5,-2.,1.
*NODE, NSET=QA_TEST_ALLNODES
1301,3.0,-3.
1302,4.0,-3.
1303,4.0,-2.
1304,3.0,-2.
**
1401,4.5,-3.
1402,5.5,-3.
1403,5.5,-2.
1404,4.5,-2.
**
1501,0.0,-1.5,0.
1502,1.0,-1.5,0.
1503,1.0,-0.5,0.
1504,0.0,-0.5,0.
1505,0.0,-1.5,1.
1506,1.0,-1.5,1.
1507,1.0,-0.5,1.
1508,0.0,-0.5,1.
**
*NODE,NSET=FIX3
1601,1.5,-1.5,0.
1602,2.5,-1.5,0.
1603,2.5,-0.5,0.
1604,1.5,-0.5,0.
1605,1.5,-1.5,1.
1606,2.5,-1.5,1.
1607,2.5,-0.5,1.
1608,1.5,-0.5,1.
*NODE, NSET=QA_TEST_ALLNODES
1701,3.0,-1.5
1702,4.0,-1.5
1703,4.0,-0.5
1704,3.0,-0.5
**
1801,4.5,-1.5
1802,5.5,-1.5
1803,5.5,-0.5
1804,4.5,-0.5
*NSET,NSET=BOTL
1101,1102,1105,1106, 1201,1202,1205,1206, 1301,1302,1401,1402
*NSET,NSET=BOTU
1501,1502,1505,1506, 1601,1602,1605,1606,1701,1702,1801,1802
*NSET,NSET=TOPL
1103,1104,1107,1108, 1203,1204,1207,1208, 1303,1304,1403,1404
*NSET,NSET=TOPU
1503,1504,1507,1508, 1603,1604,1607,1608, 1703,1704,1803,1804
*ELEMENT,TYPE=C3D8R,ELSET=C3D_11
11, 1101,1102,1103,1104,1105,1106,1107,1108
*ELEMENT,TYPE=C3D8R,ELSET=C3D_12
12, 1201,1202,1203,1204,1205,1206,1207,1208
*ELEMENT,TYPE=C3D8R,ELSET=C3D_13
15, 1501,1502,1503,1504,1505,1506,1507,1508
*ELEMENT,TYPE=C3D8R,ELSET=C3D_14
16, 1601,1602,1603,1604,1605,1606,1607,1608
*ELEMENT,TYPE=CPE4R,ELSET=GPE_11
13, 1301,1302,1303,1304
*ELEMENT,TYPE=CPE4R,ELSET=GPE_12
17, 1701,1702,1703,1704
*ELEMENT,TYPE=CPS4R,ELSET=PS_11
14, 1401,1402,1403,1404
*ELEMENT,TYPE=CPS4R,ELSET=PS_12
18, 1801,1802,1803,1804
*SOLID SECTION,ELSET=C3D_1,MATERIAL=M1
*SOLID SECTION,ELSET=C3D_2,MATERIAL=M1
*SOLID SECTION,ELSET=C3D_3,MATERIAL=M2
*SOLID SECTION,ELSET=C3D_4,MATERIAL=M2
*SOLID SECTION,ELSET=GPE_1,MATERIAL=M1
*SOLID SECTION,ELSET=GPE_2,MATERIAL=M2
*SOLID SECTION,ELSET=PS_1,MATERIAL=M1
*SOLID SECTION,ELSET=PS_2,MATERIAL=M2
*SOLID SECTION,ELSET=C3D_11,MATERIAL=M11
*SOLID SECTION,ELSET=C3D_12,MATERIAL=M11
*SOLID SECTION,ELSET=C3D_13,MATERIAL=M12
*SOLID SECTION,ELSET=C3D_14,MATERIAL=M12
*SOLID SECTION,ELSET=GPE_11,MATERIAL=M11
*SOLID SECTION,ELSET=GPE_12,MATERIAL=M12
*SOLID SECTION,ELSET=PS_11,MATERIAL=M11
*SOLID SECTION,ELSET=PS_12,MATERIAL=M12
*MATERIAL,NAME=M1
*DENSITY
1000., 
*ELASTIC
2.E7,0.3
*DRUCKER PRAGER
40.,0.8,20.
*DRUCKER PRAGER HARDENING
4.E4,0.0
*MATERIAL,NAME=M11
*DENSITY
1000., 
*ELASTIC
2.E7,0.3
*DRUCKER PRAGER
40.,1.0,20.
*DRUCKER PRAGER HARDENING
4.E4,0.0
*MATERIAL,NAME=M2
*DENSITY
1000., 
*ELASTIC
2.E7,0.3
*DRUCKER PRAGER
40.,0.8,40.
*DRUCKER PRAGER HARDENING
4.E4,0.0
*MATERIAL,NAME=M12
*DENSITY
1000., 
*ELASTIC
2.E7,0.3
*DRUCKER PRAGER
40.,1.0,40.
*DRUCKER PRAGER HARDENING
4.E4,0.0
*BOUNDARY
FIX3,3,3
**
**
*ELSET,ELSET=QA_TEST_ALLELEMS
C3D_1, C3D_2, C3D_3, C3D_4
GPE_1, GPE_2, PS_1, PS_2
C3D_11, C3D_12, C3D_13, C3D_14
GPE_11, GPE_12, PS_11, PS_12
**
*NSET, NSET=QA_TEST_ALLNODES
QA_TEST_ALLNODES, FIX3, FIX3, BOTL
BOTU, TOPL, TOPU, FIX3
FIX3, BOTL, BOTU, TOPL
TOPU, 
**
*STEP
*DYNAMIC,EXPLICIT,improved DT=no
,0.5
*BOUNDARY,TYPE=VELOCITY
 TOPU,2,, -0.08
 BOTU,2,,0.08
 TOPL,1,1,0.4
 TOPL,2,2,0.0
 BOTL,1,2,0.0
**
*FILE OUTPUT,NUMBER INTERVAL=2,TIME MARKS=YES
*EL FILE
 MISES,PEEQ,
 PE,
 ELEN,
*NODE FILE
 U,
**
*OUTPUT, FIELD, NUMBER INTERVAL=2, TIMEMARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_ALLELEMS
 MISES,PEEQ,
 PE,
 ELEN,
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
 U,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.25
***OUTPUT,HISTORY,VARIABLE=PRESELECT
*END STEP