*HEADING DEEP DRAWING OF A SQUARE BOX (1/4 SYMMETRY MODEL) ** ** MESH FOR THE PUNCH ** *PREPRINT, CONTACT=YES, MODEL=YES, HISTORY=YES *RESTART,WRITE,NUMBER INTERVAL=2,TIMEMARKS=NO ***PREPRINT,ECHO=NO,MODEL=NO,HISTORY=NO ** *NODE 1, 0., 0., 0.00041 13, 0.04000, 0., 0.00041 25, 0.05000, 0., 0.01041 37, 0.05000, 0., 0.07000 99991, 0.04000, 0., 0.01041 445, 0., 0.04000, 0.00041 457, 0.04000, 0.04000, 0.00041 469, 0.05000, 0.04000, 0.01041 481, 0.05000, 0.04000, 0.07000 99992, 0.04000, 0.04000, 0.01041 1057, 0., 0.04000, 0.00041 1069, 0., 0.05000, 0.01041 1081, 0., 0.05000, 0.07000 99993, 0., 0.04000, 0.01041 99901, 0., 0., 0.00041 *NGEN 1, 13, 1 25, 37, 1 445, 457, 1 469, 481, 1 1069, 1081, 1 *NGEN,LINE=C 13, 25, 1, 99991 457, 469, 1, 99992 1057, 1069, 1, 99993 *NSET,NSET=P1,GENERATE 1,37,1 *NSET,NSET=P2,GENERATE 445,481,1 *NFILL P1,P2,12,37 ** *NSET,NSET=P2C,GENERATE 457,481,1 *NCOPY,CHANGE NUMBER=25,OLD SET=P2C,SHIFT,MULTIPLE=12 0.,0.,0. 0.040,0.040,0.,0.040,0.040,1.,7.5 *NSET,NSET=P3,GENERATE 757,781,1 *NSET,NSET=P4,GENERATE 1057,1081,1 *NFILL P3,P4,12,25 ** *NSET,NSET=PXZ P1, *NSET,NSET=PYZ,GENERATE 1,445,37 1058,1081,1 ** ** *ELEMENT,TYPE=R3D4,ELSET=PUNCH 1, 1, 2, 39, 38 433, 457, 457, 458, 483 434, 457, 457, 483, 508 435, 457, 457, 508, 533 436, 457, 457, 533, 558 437, 457, 457, 558, 583 438, 457, 457, 583, 608 439, 457, 457, 608, 633 440, 457, 457, 633, 658 441, 457, 457, 658, 683 442, 457, 457, 683, 708 443, 457, 457, 708, 733 444, 457, 457, 733, 758 445, 456, 457, 758, 783 446, 455, 456, 783, 808 447, 454, 455, 808, 833 448, 453, 454, 833, 858 449, 452, 453, 858, 883 450, 451, 452, 883, 908 451, 450, 451, 908, 933 452, 449, 450, 933, 958 453, 448, 449, 958, 983 454, 447, 448, 983, 1008 455, 446, 447, 1008, 1033 456, 445, 446, 1033, 1058 457, 483, 458, 459, 484 *ELGEN,ELSET=PUNCH 1, 36, 1, 1, 12, 37, 36 457, 24, 25, 1, 23, 1, 24 *RIGID BODY,ELSET=PUNCH,REF NODE=99901 ** ** MESH FOR THE HOLDER ** *NODE 1101, 0.06000, 0., 0.01000 1103, 0.06000, 0., 0.00041 1115, 0.13000, 0., 0.00041 1117, 0.13000, 0., 0.01000 1305, 0.06000, 0.05000, 0.01000 1307, 0.06000, 0.05000, 0.00041 1319, 0.13000, 0.05000, 0.00041 1321, 0.13000, 0.05000, 0.01000 1713, 0., 0.06000, 0.01000 1715, 0., 0.06000, 0.00041 1727, 0., 0.13000, 0.00041 1729, 0., 0.13000, 0.01000 99902, 0.06000, 0., 0.01000 ** *NGEN,NSET=H1 1101,1103,1 1103,1115,1 1115,1117,1 *NGEN,NSET=H2 1305,1307,1 1307,1319,1 1319,1321,1 *NFILL H1,H2,12,17 *NCOPY,CHANGE NUMBER=17,OLD SET=H2,SHIFT,MULTIPLE=12 0.,0.,0. 0.050,0.050,0.,0.050,0.050,1.,7.5 *NSET,NSET=H3,GENERATE 1509,1525,1 *NGEN,NSET=H4 1713,1715,1 1715,1727,1 1727,1729,1 *NFILL H3,H4,12,17 ** *ELEMENT,TYPE=R3D4 1010,1101,1118,1119,1102 *ELGEN,ELSET=HOLDER 1010,16,1,1,36,17,16 *RIGID BODY,ELSET=HOLDER,REF NODE=99902 ** *ELEMENT,TYPE=MASS,ELSET=EMASS 1009,99902 *MASS,ELSET=EMASS 0.6396, ** ** MESH FOR THE DIE ** *NODE 1801, 0.05125, 0., -0.06000 1813, 0.05125, 0., -0.01041 1825, 0.06125, 0., -0.00041 1837, 0.13000, 0., -0.00041 99994, 0.06125, 0., -0.01041 2245, 0.05125, 0.04125, -0.06000 2257, 0.05125, 0.04125, -0.01041 2269, 0.06125, 0.04125, -0.00041 2281, 0.13000, 0.04125, -0.00041 99995, 0.06125, 0.04125, -0.01041 3133, 0., 0.05125, -0.06000 3145, 0., 0.05125, -0.01041 3157, 0., 0.06125, -0.00041 3169, 0., 0.13000, -0.00041 99996, 0., 0.06125, -0.01041 99903, 0.05125, 0., -0.06000 *NGEN 1801, 1813, 1 1825, 1837, 1 2245, 2257, 1 2269, 2281, 1 3133, 3145, 1 3157, 3169, 1 *NGEN,LINE=C 1813, 1825, 1, 99994 2257, 2269, 1, 99995 3145, 3157, 1, 99996 *NSET,NSET=D1,GENERATE 1801,1837,1 *NSET,NSET=D2,GENERATE 2245,2281,1 *NFILL D1,D2,12,37 *NCOPY,CHANGE NUMBER=37,OLD SET=D2,SHIFT,MULTIPLE=12 0.,0.,0. 0.04125,0.04125,0.,0.04125,0.04125,1.,7.5 *NSET,NSET=D3,GENERATE 2689,2725,1 *NSET,NSET=D4,GENERATE 3133,3169,1 *NFILL D3,D4,12,37 ** *ELEMENT,TYPE=R3D4 1586,1801,1802,1839,1838 *ELGEN,ELSET=DIE 1586,36,1,1,36,37,36 *RIGID BODY,ELSET=DIE,REF NODE=99903 ** ** MESH FOR THE BLANK (70 X 70) ** *NODE 4001, 0., 0., 0. 4071, 0.10000, 0., 0. 8971, 0., 0.10000, 0. 9041, 0.10000, 0.10000, 0. *NGEN,NSET=B1 4001,4071,1 *NGEN,NSET=B2 8971,9041,1 *NFILL B1,B2,70,71 *NSET,NSET=B4,GENERATE 4001,8971,71 ** *ELEMENT,TYPE=S4R 3001,4001,4002,4073,4072 *ELGEN,ELSET=BLANK 3001,70,1,1,70,71,70 *SHELL SECTION,MATERIAL=STEEL,ELSET=BLANK,SECTION INTEGRA=GAUSS, CONTROLS=HGC .00082,5 *SECTION CONTROLS,NAME=HGC,HOURGLASS=ENHANCED ** ** Define material properties for STEEL ** *MATERIAL,NAME=STEEL *DENSITY 7800., *ELASTIC 2.1E11,0.3 *PLASTIC 0.91294E+08, 0.00000E+00 0.10129E+09, 0.21052E-03 0.11129E+09, 0.52686E-03 0.12129E+09, 0.97685E-03 0.13129E+09, 0.15923E-02 0.14129E+09, 0.24090E-02 0.15129E+09, 0.34674E-02 0.16129E+09, 0.48120E-02 0.17129E+09, 0.64921E-02 0.18129E+09, 0.85618E-02 0.19129E+09, 0.11080E-01 0.20129E+09, 0.14110E-01 0.21129E+09, 0.17723E-01 0.22129E+09, 0.21991E-01 0.23129E+09, 0.26994E-01 0.24129E+09, 0.32819E-01 0.25129E+09, 0.39556E-01 0.26129E+09, 0.47301E-01 0.27129E+09, 0.56159E-01 0.28129E+09, 0.66236E-01 0.29129E+09, 0.77648E-01 0.30129E+09, 0.90516E-01 0.31129E+09, 0.10497E+00 0.32129E+09, 0.12114E+00 0.33129E+09, 0.13916E+00 0.34129E+09, 0.15919E+00 0.35129E+09, 0.18138E+00 0.36129E+09, 0.20588E+00 0.37129E+09, 0.23287E+00 0.38129E+09, 0.26252E+00 0.39129E+09, 0.29502E+00 0.40129E+09, 0.33054E+00 0.41129E+09, 0.36929E+00 0.42129E+09, 0.41147E+00 0.43129E+09, 0.45729E+00 0.44129E+09, 0.50696E+00 0.45129E+09, 0.56073E+00 0.46129E+09, 0.61881E+00 0.47129E+09, 0.68145E+00 0.48129E+09, 0.74890E+00 0.49129E+09, 0.82142E+00 0.50129E+09, 0.89928E+00 0.51129E+09, 0.98274E+00 0.52129E+09, 0.10721E+01 ** ** Define surfaces ** *surface, NAME=PUNCH PUNCH,SNEG *surface, NAME=DIE DIE,SPOS *surface, NAME=BLANK BLANK, ** ** Apply symmetry boundary conditions ** *BOUNDARY D1,YSYMM PXZ,YSYMM H1,YSYMM D4,XSYMM PYZ,XSYMM H4,XSYMM ** B1,YSYMM B4,XSYMM ** 99901,1,2 99901,4,6 99902,1,2 99902,4,6 99903,1,6 ** *ELSET,ELSET=ELHIST 5273, *ELSET, ELSET=QA_TEST ELHIST, 5277, 5844, 3001, 3070, 7831 *NSET, NSET=QA_TEST 6522, 6912, 4042, 4071, 4001, 8993 ** ** Simulate the deep drawing operation ** ************************************************************* **General contact definition *CONTACT *CONTACT INCLUSIONS, ALL ELEMENT BASED *CONTACT PROPERTY ASSIGNMENT PUNCH, BLANK, PUNCH_T DIE, BLANK, DIE_B *SURFACE PROPERTY ASSIGNMENT, PROPERTY=THICKNESS BLANK,THINNING *CONTACT FORMULATION, TYPE=PURE MASTER-SLAVE PUNCH, BLANK, MASTER DIE, BLANK, MASTER *SURFACE INTERACTION, NAME = PUNCH_T *FRICTION 0.25, *SURFACE INTERACTION, NAME = DIE_B *FRICTION 0.125, ************************************************************* *NSET,NSET=NPUNCH 99901, *NSET,NSET=HOLDER 99902, ************************************************************* *STEP *DYNAMIC,EXPLICIT ,.0036 *AMPLITUDE,DEFINITION=SMOOTH STEP,NAME=SMOOTH 0.,0., 0.0036, 1. *BOUNDARY, AMPLITUDE=SMOOTH 99901,3,3,-0.036 *CLOAD 99902,3,-2.287E4 ***************************** **** OUTPUT REQUESTS **** ***************************** *OUTPUT,FIELD,VAR=PRESELECT, NUMBER INTERVAL=2 *ELEMENT OUTPUT STH,PEEQ *OUTPUT,HISTORY,TIME INTERVAL=1.8E-4 *ELEMENT OUTPUT,ELSET=ELHIST STH, PEEQ, *NODE OUTPUT,NSET=NPUNCH U3,RF3 *NODE OUTPUT,NSET=HOLDER U3, ************************************** **** OUTPUT FOR QA PURPOSES **** ************************************** *OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1 *ELEMENT OUTPUT, ELSET=QA_TEST PEEQ, STH *NODE OUTPUT, NSET=QA_TEST U, *OUTPUT, HIST, TIME INTERVAL=1.8E-5 *ENERGY OUTPUT, VAR=PRESELECT *END STEP