*HEADING COPPER ROD IMPACT (3D MODEL) ** ** This example simulates the high velocity impact of a copper rod ** onto a rigid wall. Similar tests have been performed to determine ** the material constants for high-pressure equations of state, and ** the test is sometimes described as the Taylor bar experiment. ** Extremely high plastic strains develop at the crushed end of the ** rod, resulting in severe local mesh distortion. ** *RESTART, WRITE, NUMBER INTERVAL=4, TIMEMARKS=NO *PREPRINT,ECHO=NO,MODEL=NO,HISTORY=NO ** ** Generate the nodes at one end ** *NODE, NSET=CORE 1, 0.0000E+00, 0.0000E+00, 0.0000E+00 2, 1.6000E-04, 0.0000E+00, 0.0000E+00 3, 1.5624E-04, 1.5987E-04, 0.0000E+00 4, 0.0000E+00, 1.6000E-04, 0.0000E+00 5, 3.2000E-04, 0.0000E+00, 0.0000E+00 6, 3.1343E-04, 1.5827E-04, 0.0000E+00 7, 3.0810E-04, 3.1461E-04, 0.0000E+00 8, 1.5342E-04, 3.1829E-04, 0.0000E+00 9, 0.0000E+00, 3.2000E-04, 0.0000E+00 10, 4.8000E-04, 0.0000E+00, 0.0000E+00 11, 4.7133E-04, 1.5558E-04, 0.0000E+00 12, 4.6374E-04, 3.0945E-04, 0.0000E+00 13, 4.5698E-04, 4.6205E-04, 0.0000E+00 14, 3.0371E-04, 4.6949E-04, 0.0000E+00 15, 1.5129E-04, 4.7562E-04, 0.0000E+00 16, 0.0000E+00, 4.8000E-04, 0.0000E+00 17, 6.4000E-04, 0.0000E+00, 0.0000E+00 18, 6.2970E-04, 1.5215E-04, 0.0000E+00 19, 6.2002E-04, 3.0331E-04, 0.0000E+00 20, 6.1080E-04, 4.5372E-04, 0.0000E+00 21, 6.0190E-04, 6.0364E-04, 0.0000E+00 22, 4.5076E-04, 6.1379E-04, 0.0000E+00 23, 2.9995E-04, 6.2340E-04, 0.0000E+00 24, 1.4965E-04, 6.3223E-04, 0.0000E+00 25, 0.0000E+00, 6.4000E-04, 0.0000E+00 26, 8.0000E-04, 0.0000E+00, 0.0000E+00 27, 7.8830E-04, 1.4835E-04, 0.0000E+00 28, 7.7660E-04, 2.9667E-04, 0.0000E+00 29, 7.6490E-04, 4.4496E-04, 0.0000E+00 30, 7.5320E-04, 5.9323E-04, 0.0000E+00 31, 7.4150E-04, 7.4150E-04, 0.0000E+00 32, 5.9315E-04, 7.5331E-04, 0.0000E+00 33, 4.4482E-04, 7.6509E-04, 0.0000E+00 34, 2.9651E-04, 7.7682E-04, 0.0000E+00 35, 1.4823E-04, 7.8847E-04, 0.0000E+00 36, 0.0000E+00, 8.0000E-04, 0.0000E+00 37, 9.6000E-04, 0.0000E+00, 0.0000E+00 38, 9.4691E-04, 1.4455E-04, 0.0000E+00 39, 9.3320E-04, 2.9002E-04, 0.0000E+00 40, 9.1901E-04, 4.3619E-04, 0.0000E+00 41, 9.0450E-04, 5.8282E-04, 0.0000E+00 42, 8.8984E-04, 7.2969E-04, 0.0000E+00 43, 8.7518E-04, 8.7657E-04, 0.0000E+00 44, 7.2979E-04, 8.8977E-04, 0.0000E+00 45, 5.8441E-04, 9.0297E-04, 0.0000E+00 46, 4.3889E-04, 9.1638E-04, 0.0000E+00 47, 2.9307E-04, 9.3022E-04, 0.0000E+00 48, 1.4682E-04, 9.4469E-04, 0.0000E+00 49, 0.0000E+00, 9.6000E-04, 0.0000E+00 50, 1.1200E-03, 0.0000E+00, 0.0000E+00 51, 1.1053E-03, 1.4111E-04, 0.0000E+00 52, 1.0895E-03, 2.8386E-04, 0.0000E+00 53, 1.0728E-03, 4.2785E-04, 0.0000E+00 54, 1.0556E-03, 5.7267E-04, 0.0000E+00 55, 1.0382E-03, 7.1791E-04, 0.0000E+00 56, 1.0207E-03, 8.6316E-04, 0.0000E+00 57, 1.0035E-03, 1.0080E-03, 0.0000E+00 58, 8.6066E-04, 1.0232E-03, 0.0000E+00 59, 7.1809E-04, 1.0380E-03, 0.0000E+00 60, 5.7552E-04, 1.0529E-03, 0.0000E+00 61, 4.3269E-04, 1.0681E-03, 0.0000E+00 62, 2.8933E-04, 1.0841E-03, 0.0000E+00 63, 1.4519E-04, 1.1013E-03, 0.0000E+00 64, 0.0000E+00, 1.1200E-03, 0.0000E+00 65, 1.2800E-03, 0.0000E+00, 0.0000E+00 66, 1.2632E-03, 1.3839E-04, 0.0000E+00 67, 1.2451E-03, 2.7868E-04, 0.0000E+00 68, 1.2261E-03, 4.2038E-04, 0.0000E+00 69, 1.2064E-03, 5.6304E-04, 0.0000E+00 70, 1.1865E-03, 7.0618E-04, 0.0000E+00 71, 1.1665E-03, 8.4934E-04, 0.0000E+00 72, 1.1469E-03, 9.9205E-04, 0.0000E+00 73, 1.1278E-03, 1.1338E-03, 0.0000E+00 74, 9.8683E-04, 1.1521E-03, 0.0000E+00 75, 8.4646E-04, 1.1694E-03, 0.0000E+00 76, 7.0639E-04, 1.1863E-03, 0.0000E+00 77, 5.6633E-04, 1.2032E-03, 0.0000E+00 78, 4.2595E-04, 1.2206E-03, 0.0000E+00 79, 2.8498E-04, 1.2389E-03, 0.0000E+00 80, 1.4310E-04, 1.2585E-03, 0.0000E+00 81, 0.0000E+00, 1.2800E-03, 0.0000E+00 82, 1.4400E-03, 0.0000E+00, 0.0000E+00 83, 1.4204E-03, 1.3677E-04, 0.0000E+00 84, 1.3998E-03, 2.7496E-04, 0.0000E+00 85, 1.3785E-03, 4.1422E-04, 0.0000E+00 86, 1.3568E-03, 5.5419E-04, 0.0000E+00 87, 1.3348E-03, 6.9453E-04, 0.0000E+00 88, 1.3128E-03, 8.3489E-04, 0.0000E+00 89, 1.2910E-03, 9.7490E-04, 0.0000E+00 90, 1.2697E-03, 1.1142E-03, 0.0000E+00 91, 1.2491E-03, 1.2525E-03, 0.0000E+00 92, 1.1097E-03, 1.2742E-03, 0.0000E+00 93, 9.7097E-04, 1.2950E-03, 0.0000E+00 94, 8.3272E-04, 1.3150E-03, 0.0000E+00 95, 6.9469E-04, 1.3347E-03, 0.0000E+00 96, 5.5667E-04, 1.3543E-03, 0.0000E+00 97, 4.1842E-04, 1.3744E-03, 0.0000E+00 98, 2.7971E-04, 1.3951E-03, 0.0000E+00 99, 1.4031E-04, 1.4169E-03, 0.0000E+00 100, 0.0000E+00, 1.4400E-03, 0.0000E+00 101, 1.6000E-03, 0.0000E+00, 0.0000E+00 102, 1.5766E-03, 1.3660E-04, 0.0000E+00 103, 1.5532E-03, 2.7320E-04, 0.0000E+00 104, 1.5298E-03, 4.0980E-04, 0.0000E+00 105, 1.5064E-03, 5.4640E-04, 0.0000E+00 106, 1.4830E-03, 6.8300E-04, 0.0000E+00 107, 1.4596E-03, 8.1960E-04, 0.0000E+00 108, 1.4362E-03, 9.5620E-04, 0.0000E+00 109, 1.4128E-03, 1.0928E-03, 0.0000E+00 110, 1.3894E-03, 1.2294E-03, 0.0000E+00 111, 1.3660E-03, 1.3660E-03, 0.0000E+00 112, 1.2294E-03, 1.3894E-03, 0.0000E+00 113, 1.0928E-03, 1.4128E-03, 0.0000E+00 114, 9.5620E-04, 1.4362E-03, 0.0000E+00 115, 8.1960E-04, 1.4596E-03, 0.0000E+00 116, 6.8300E-04, 1.4830E-03, 0.0000E+00 117, 5.4640E-04, 1.5064E-03, 0.0000E+00 118, 4.0980E-04, 1.5298E-03, 0.0000E+00 119, 2.7320E-04, 1.5532E-03, 0.0000E+00 120, 1.3660E-04, 1.5766E-03, 0.0000E+00 121, 0.0000E+00, 1.6000E-03, 0.0000E+00 *NSET, NSET=INNER, GENERATE 101, 121, 1 *NODE, NSET=ENDPTS 311, 0.3200E-02, 0.0000E+00, 0.0000E+00 331, 0.0000E+00, 0.3200E-02, 0.0000E+00 *NGEN, NSET=OUTER, LINE=C 311, 331, 1, 1 *NFILL INNER, OUTER, 10, 21 *NSET, NSET=HEAD, GENERATE 1, 331, 1 ** ** Extrude a three-dimensional node mesh ** *NCOPY, SHIFT, OLD SET=HEAD, NEW SET=TAIL, CHANGE NUMBER=23832 0., 0., .0324 0., 0., -1., 0., 0., 1., 0. *NFILL, NSET=ROD HEAD, TAIL, 72, 331 *NSET, NSET=SPLANE1, GENERATE 1, 23833, 331 2, 23834, 331 5, 23837, 331 10, 23842, 331 17, 23849, 331 26, 23858, 331 37, 23869, 331 50, 23882, 331 65, 23897, 331 82, 23914, 331 101, 23933, 331 ** 122, 23954, 331 143, 23975, 331 164, 23996, 331 185, 24017, 331 206, 24038, 331 227, 24059, 331 248, 24080, 331 269, 24101, 331 290, 24122, 331 311, 24143, 331 *NSET, NSET=SPLANE2, GENERATE 1, 23833, 331 4, 23836, 331 9, 23841, 331 16, 23848, 331 25, 23857, 331 36, 23868, 331 49, 23881, 331 64, 23896, 331 81, 23913, 331 100, 23932, 331 121, 23953, 331 ** 142, 23974, 331 163, 23995, 331 184, 24016, 331 205, 24037, 331 226, 24058, 331 247, 24079, 331 268, 24100, 331 289, 24121, 331 310, 24142, 331 331, 24163, 331 *NSET, NSET=BODY, GENERATE 332, 24163, 1 ** ** Fill-in the elements ** *ELEMENT, TYPE=C3D8R 1, 1, 2, 3, 4, 332, 333, 334, 335 2, 2, 5, 6, 3, 333, 336, 337, 334 3, 5, 10, 11, 6, 336, 341, 342, 337 4, 10, 17, 18, 11, 341, 348, 349, 342 5, 17, 26, 27, 18, 348, 357, 358, 349 6, 26, 37, 38, 27, 357, 368, 369, 358 7, 37, 50, 51, 38, 368, 381, 382, 369 8, 50, 65, 66, 51, 381, 396, 397, 382 9, 65, 82, 83, 66, 396, 413, 414, 397 10, 82, 101, 102, 83, 413, 432, 433, 414 11, 4, 3, 8, 9, 335, 334, 339, 340 12, 3, 6, 7, 8, 334, 337, 338, 339 13, 6, 11, 12, 7, 337, 342, 343, 338 14, 11, 18, 19, 12, 342, 349, 350, 343 15, 18, 27, 28, 19, 349, 358, 359, 350 16, 27, 38, 39, 28, 358, 369, 370, 359 17, 38, 51, 52, 39, 369, 382, 383, 370 18, 51, 66, 67, 52, 382, 397, 398, 383 19, 66, 83, 84, 67, 397, 414, 415, 398 20, 83, 102, 103, 84, 414, 433, 434, 415 21, 9, 8, 15, 16, 340, 339, 346, 347 22, 8, 7, 14, 15, 339, 338, 345, 346 23, 7, 12, 13, 14, 338, 343, 344, 345 24, 12, 19, 20, 13, 343, 350, 351, 344 25, 19, 28, 29, 20, 350, 359, 360, 351 26, 28, 39, 40, 29, 359, 370, 371, 360 27, 39, 52, 53, 40, 370, 383, 384, 371 28, 52, 67, 68, 53, 383, 398, 399, 384 29, 67, 84, 85, 68, 398, 415, 416, 399 30, 84, 103, 104, 85, 415, 434, 435, 416 31, 16, 15, 24, 25, 347, 346, 355, 356 32, 15, 14, 23, 24, 346, 345, 354, 355 33, 14, 13, 22, 23, 345, 344, 353, 354 34, 13, 20, 21, 22, 344, 351, 352, 353 35, 20, 29, 30, 21, 351, 360, 361, 352 36, 29, 40, 41, 30, 360, 371, 372, 361 37, 40, 53, 54, 41, 371, 384, 385, 372 38, 53, 68, 69, 54, 384, 399, 400, 385 39, 68, 85, 86, 69, 399, 416, 417, 400 40, 85, 104, 105, 86, 416, 435, 436, 417 41, 25, 24, 35, 36, 356, 355, 366, 367 42, 24, 23, 34, 35, 355, 354, 365, 366 43, 23, 22, 33, 34, 354, 353, 364, 365 44, 22, 21, 32, 33, 353, 352, 363, 364 45, 21, 30, 31, 32, 352, 361, 362, 363 46, 30, 41, 42, 31, 361, 372, 373, 362 47, 41, 54, 55, 42, 372, 385, 386, 373 48, 54, 69, 70, 55, 385, 400, 401, 386 49, 69, 86, 87, 70, 400, 417, 418, 401 50, 86, 105, 106, 87, 417, 436, 437, 418 51, 36, 35, 48, 49, 367, 366, 379, 380 52, 35, 34, 47, 48, 366, 365, 378, 379 53, 34, 33, 46, 47, 365, 364, 377, 378 54, 33, 32, 45, 46, 364, 363, 376, 377 55, 32, 31, 44, 45, 363, 362, 375, 376 56, 31, 42, 43, 44, 362, 373, 374, 375 57, 42, 55, 56, 43, 373, 386, 387, 374 58, 55, 70, 71, 56, 386, 401, 402, 387 59, 70, 87, 88, 71, 401, 418, 419, 402 60, 87, 106, 107, 88, 418, 437, 438, 419 61, 49, 48, 63, 64, 380, 379, 394, 395 62, 48, 47, 62, 63, 379, 378, 393, 394 63, 47, 46, 61, 62, 378, 377, 392, 393 64, 46, 45, 60, 61, 377, 376, 391, 392 65, 45, 44, 59, 60, 376, 375, 390, 391 66, 44, 43, 58, 59, 375, 374, 389, 390 67, 43, 56, 57, 58, 374, 387, 388, 389 68, 56, 71, 72, 57, 387, 402, 403, 388 69, 71, 88, 89, 72, 402, 419, 420, 403 70, 88, 107, 108, 89, 419, 438, 439, 420 71, 64, 63, 80, 81, 395, 394, 411, 412 72, 63, 62, 79, 80, 394, 393, 410, 411 73, 62, 61, 78, 79, 393, 392, 409, 410 74, 61, 60, 77, 78, 392, 391, 408, 409 75, 60, 59, 76, 77, 391, 390, 407, 408 76, 59, 58, 75, 76, 390, 389, 406, 407 77, 58, 57, 74, 75, 389, 388, 405, 406 78, 57, 72, 73, 74, 388, 403, 404, 405 79, 72, 89, 90, 73, 403, 420, 421, 404 80, 89, 108, 109, 90, 420, 439, 440, 421 81, 81, 80, 99, 100, 412, 411, 430, 431 82, 80, 79, 98, 99, 411, 410, 429, 430 83, 79, 78, 97, 98, 410, 409, 428, 429 84, 78, 77, 96, 97, 409, 408, 427, 428 85, 77, 76, 95, 96, 408, 407, 426, 427 86, 76, 75, 94, 95, 407, 406, 425, 426 87, 75, 74, 93, 94, 406, 405, 424, 425 88, 74, 73, 92, 93, 405, 404, 423, 424 89, 73, 90, 91, 92, 404, 421, 422, 423 90, 90, 109, 110, 91, 421, 440, 441, 422 91, 100, 99, 120, 121, 431, 430, 451, 452 92, 99, 98, 119, 120, 430, 429, 450, 451 93, 98, 97, 118, 119, 429, 428, 449, 450 94, 97, 96, 117, 118, 428, 427, 448, 449 95, 96, 95, 116, 117, 427, 426, 447, 448 96, 95, 94, 115, 116, 426, 425, 446, 447 97, 94, 93, 114, 115, 425, 424, 445, 446 98, 93, 92, 113, 114, 424, 423, 444, 445 99, 92, 91, 112, 113, 423, 422, 443, 444 100, 91, 110, 111, 112, 422, 441, 442, 443 *ELGEN 1, 72, 331, 300 2, 72, 331, 300 3, 72, 331, 300 4, 72, 331, 300 5, 72, 331, 300 6, 72, 331, 300 7, 72, 331, 300 8, 72, 331, 300 9, 72, 331, 300 10, 72, 331, 300 11, 72, 331, 300 12, 72, 331, 300 13, 72, 331, 300 14, 72, 331, 300 15, 72, 331, 300 16, 72, 331, 300 17, 72, 331, 300 18, 72, 331, 300 19, 72, 331, 300 20, 72, 331, 300 21, 72, 331, 300 22, 72, 331, 300 23, 72, 331, 300 24, 72, 331, 300 25, 72, 331, 300 26, 72, 331, 300 27, 72, 331, 300 28, 72, 331, 300 29, 72, 331, 300 30, 72, 331, 300 31, 72, 331, 300 32, 72, 331, 300 33, 72, 331, 300 34, 72, 331, 300 35, 72, 331, 300 36, 72, 331, 300 37, 72, 331, 300 38, 72, 331, 300 39, 72, 331, 300 40, 72, 331, 300 41, 72, 331, 300 42, 72, 331, 300 43, 72, 331, 300 44, 72, 331, 300 45, 72, 331, 300 46, 72, 331, 300 47, 72, 331, 300 48, 72, 331, 300 49, 72, 331, 300 50, 72, 331, 300 51, 72, 331, 300 52, 72, 331, 300 53, 72, 331, 300 54, 72, 331, 300 55, 72, 331, 300 56, 72, 331, 300 57, 72, 331, 300 58, 72, 331, 300 59, 72, 331, 300 60, 72, 331, 300 61, 72, 331, 300 62, 72, 331, 300 63, 72, 331, 300 64, 72, 331, 300 65, 72, 331, 300 66, 72, 331, 300 67, 72, 331, 300 68, 72, 331, 300 69, 72, 331, 300 70, 72, 331, 300 71, 72, 331, 300 72, 72, 331, 300 73, 72, 331, 300 74, 72, 331, 300 75, 72, 331, 300 76, 72, 331, 300 77, 72, 331, 300 78, 72, 331, 300 79, 72, 331, 300 80, 72, 331, 300 81, 72, 331, 300 82, 72, 331, 300 83, 72, 331, 300 84, 72, 331, 300 85, 72, 331, 300 86, 72, 331, 300 87, 72, 331, 300 88, 72, 331, 300 89, 72, 331, 300 90, 72, 331, 300 91, 72, 331, 300 92, 72, 331, 300 93, 72, 331, 300 94, 72, 331, 300 95, 72, 331, 300 96, 72, 331, 300 97, 72, 331, 300 98, 72, 331, 300 99, 72, 331, 300 100, 72, 331, 300 *ELEMENT, TYPE=C3D8R 101, 101, 122, 123, 102, 432, 453, 454, 433 *ELGEN 101, 10, 21, 1, 20, 1, 10, 72, 331, 300 *ELSET, ELSET=ROD, GENERATE 1, 21600, 1 ** ** Representative material properties for copper ** *SOLID SECTION, ELSET=ROD, MATERIAL=COPPER *MATERIAL, NAME=COPPER *ELASTIC 110.E9, .3 *PLASTIC 314.E6, *DENSITY 8.97E3, ** *ELSET, ELSET=ELOUT 1,110,300 *NSET, NSET=NOUT 1,311,331 ** ** Zero displacement b.c.'s are prescribed on the symmetry planes. ** To simulate the impact of the rod on a (frictionless) rigid wall, ** zero axial displacements are prescribed at one end of the rod, ** while all other nodes are subjected to an initial velocity field. ** *BOUNDARY SPLANE1, 2 SPLANE2, 1 HEAD, 3 *INITIAL CONDITIONS, TYPE=VELOCITY BODY, 3, -227. ** ** Run the problem for 80 microseconds ** *STEP *DYNAMIC,EXPLICIT ,80.E-6 ** ***OUTPUT,HISTORY,TIME INTERVAL=80.E-9 ***NODE OUTPUT,NSET=NOUT **U,V ***ELEMENT OUTPUT,ELSET=ELOUT **PEEQ,MISES ***ENERGY OUTPUT **ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLCD,ALLPD,ALLFD,ALLWK,ETOTAL,DT *OUTPUT,FIELD,VAR=PRESELECT,NUMBER INTERVAL=5 ** ** Output requests for qa testing ** *NSET,NSET=QA_TEST NOUT, *ELSET,ELSET=QA_TEST ELOUT, *OUTPUT,FIELD,NUMBER INTERVAL=1 *NODE OUTPUT,NSET=QA_TEST U, *ELEMENT OUTPUT,ELSET=QA_TEST PEEQ, *OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=2.E-5 *FILE OUTPUT, TIMEMARKS=YES, NUM=1 *EL FILE,ELSET=ELOUT PEEQ, *NODE FILE,NSET=NOUT U, *ENERGY FILE *END STEP