*HEADING FULLY COUPLED ACOUSTIC ANALYSIS OF A MUFFLER, EXPLICIT ** ** This input file needs the following files: ** **muffler_intair_nodes_xpl.inp **muffler_intair_elem_xpl.inp **muffler_shell_nodes.inp **muffler_shell_elem.inp **muffler_extair_nodes_xpl.inp **muffler_extair_elem_xpl.inp ** ** It is recommended to run the analysis ** in DOUBLE PRECISION !!! ** ** PART INSTANCE: INTERIORAIR-1 ** *Parameter f=180.0 T=1.0/f omega=2*3.14159*f *NODE,INPUT=muffler_intair_nodes_xpl.inp *ELEMENT,TYPE=AC3D4,INPUT=muffler_intair_elem_xpl.inp ** ------------------------------------------- *ELSET,ELSET=INTERIORAIR-1_I1,GENERATE 549, 1457, 1 *SOLID SECTION,ELSET=INTERIORAIR-1_I1,MATERIAL=AIRABSORB 1., ** ------------------------------------------- *ELSET,ELSET=INTERIORAIR-1_I2,GENERATE 1, 548, 1 *SOLID SECTION,ELSET=INTERIORAIR-1_I2,MATERIAL=AIR 1., ** ------------------------------------------- ** ** PART INSTANCE: MUFFLERSHELL-1 ** *NODE,INPUT=muffler_shell_nodes.inp *ELEMENT,TYPE=S4R,INPUT=muffler_shell_elem.inp ** ------------------------------------------- *ELSET,ELSET=MUFFLERSHELL-1_I1,GENERATE 1458, 1500, 1 *SHELL SECTION,ELSET=MUFFLERSHELL-1_I1,MATERIAL=STEEL 0.00075, 5 ** ------------------------------------------- ** ** PART INSTANCE: OUTERAIR-1 ** *NODE,INPUT=muffler_extair_nodes_xpl.inp *ELEMENT,TYPE=AC3D8R,INPUT=muffler_extair_elem_xpl.inp ** ------------------------------------------- *ELSET,ELSET=OUTERAIR-1_I1,GENERATE 1501, 2734, 1 *SOLID SECTION,ELSET=OUTERAIR-1_I1,MATERIAL=AIR 1., ** ------------------------------------------- *NSET,NSET=SHELLSYM,UNSORTED 3301,3376,3375,3299,3373,3297,3340,3339,3338, 3337,3336,3335,3334,3333,3332,3331,3330,3329, 3328,3327,3326,3325,3324,3323,3322,3321,3320, 3319,3318,3317,3316,3315,3314,3313,3312,3311, 3310,3309,3295,3307,3292,3304,3303,3291 ** ------------------------------------------- *NSET,NSET=SHELLCONSTRAINT,UNSORTED 3302,3377,3378,3300,3374,3298,3341,3342,3343, 3344,3345,3346,3347,3348,3349,3350,3351,3352, 3353,3354,3355,3356,3357,3358,3359,3360,3361, 3362,3363,3364,3365,3366,3367,3368,3369,3370, 3371,3372,3296,3308,3293,3305,3306,3294 ** ------------------------------------------- *NSET,NSET=SHELLDISP 3304,3325,3375 ** ------------------------------------------- *NSET,NSET=CENTERLINE 10, 11, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 831, 832, 833, 834, 835, 836, 837, 838, 839, 840, 841, 842, 843, 844, 845, 846, 847, 848, 849, 850, 851, 852, 853, 854, 855, 856, 857, 858, 859, 860, 861, 862, 863, 864, 865, 866, 867, 868, 869, 870, 871, 872, 873, 874, 875, 876, 877, 878, 879, 880, 881, 882, 883, 884, 885, 886, 887, 888, 889, 890 ** ------------------------------------------- *NSET,NSET=INLETPRESSURE 11, 12, 13, 177, 186, 618, 891, 892, 893, 913, 914, 915 ** ------------------------------------------- *ELSET,ELSET=INNERAIROUTLET_S4 288, 299, 367 *SURFACE,NAME=INNERAIROUTLET INNERAIROUTLET_S4, S4 ** ------------------------------------------- *ELSET,ELSET=INNERAIR2SHELL_S1 367, 374, 384, 385, 386, 387, 388, 392, 393, 894, 895, 896, 1192, 1195, 1196, 1197, 1198, 1201, 1202, 1203, 1210, 1213, 1214, 1216, 1237, 1238, 1239, 1240, 1243, 1248, 1250, 1330, 1331, 1334, 1337, 1341, 1359 *ELSET,ELSET=INNERAIR2SHELL_S4 370, 373, 375, 380, 383, 391, 559, 634, 894, 895, 1193, 1194, 1200, 1204, 1206, 1208, 1211, 1212, 1215, 1217, 1235, 1242, 1244, 1246, 1247, 1249, 1251, 1276, 1298, 1320, 1321, 1322, 1323, 1326, 1329, 1335, 1336, 1338, 1340, 1343, 1347, 1349 *ELSET,ELSET=INNERAIR2SHELL_S2 379, 592, 635, 858, 917, 1199, 1205, 1209, 1236, 1241, 1252, 1253, 1254, 1325, 1339, 1344, 1345, 1346, 1348, 1350, 1351, 1352, 1353, 1354, 1355, 1356, 1357, 1358, 1360, 1361, 1362 *ELSET,ELSET=INNERAIR2SHELL_S3 377, 378, 389, 390, 557, 560, 561, 568, 683, 888, 897, 1207, 1245, 1277, 1324, 1327, 1328, 1332, 1333, 1342, 1387, 1388 *SURFACE,NAME=INNERAIR2SHELL INNERAIR2SHELL_S1, S1 INNERAIR2SHELL_S4, S4 INNERAIR2SHELL_S2, S2 INNERAIR2SHELL_S3, S3 ** ------------------------------------------- *ELSET,ELSET=SHELL2INNERAIR_SNEG,GENERATE 1458, 1500, 1 *SURFACE,NAME=SHELL2INNERAIR SHELL2INNERAIR_SNEG, SNEG *ELSET,ELSET=SHELL2OUTERAIR_SPOS,GENERATE 1458, 1500, 1 *SURFACE,NAME=SHELL2OUTERAIR SHELL2OUTERAIR_SPOS, SPOS ** ------------------------------------------- *ELSET,ELSET=OUTERAIR2SHELL_S6 2100,2116,2420 *ELSET,ELSET=OUTERAIR2SHELL_S3 2100,2420,2421,2422 *ELSET,ELSET=OUTERAIR2SHELL_S4,GENERATE 1519, 1975, 19 *SURFACE,NAME=OUTERAIR2SHELL OUTERAIR2SHELL_S6, S6 OUTERAIR2SHELL_S3, S3 OUTERAIR2SHELL_S4, S4 ** ------------------------------------------- *ELSET,ELSET=SPHERE1_S3,GENERATE 2340, 2355, 1 *ELSET,ELSET=SPHERE1_S4,GENERATE 2664, 2734, 5 *SURFACE,NAME=SPHERE1 SPHERE1_S3, S3 SPHERE1_S4, S4 ** ------------------------------------------- *ELSET,ELSET=SPHERE2_S3,GENERATE 2036, 2051, 1 *ELSET,ELSET=SPHERE2_S4,GENERATE 1979, 2035, 4 *SURFACE,NAME=SPHERE2 SPHERE2_S3, S3 SPHERE2_S4, S4 ** ------------------------------------------- *ELSET,ELSET=CYLINDER_S6,GENERATE 1501, 1957, 19 *SURFACE,NAME=CYLINDER CYLINDER_S6, S6 ** ------------------------------------------- *ELSET,ELSET=INNERAIRABSORB_S3 382, 394, 403, 404, 405, 408, 411, 415, 417, 418, 419, 423, 429, 430, 431, 437, 438, 439, 448, 451, 452, 456, 457, 458, 459, 461, 462, 463, 464, 467, 468, 473, 477, 478, 482, 483, 484, 498, 499, 502, 503, 505, 506, 507 *ELSET,ELSET=INNERAIRABSORB_S2 407, 428, 436, 443, 471, 476, 504 *ELSET,ELSET=INNERAIRABSORB_S1 402, 410, 414, 416, 422, 426, 434, 435, 440, 441, 444, 450, 453, 455, 460, 466, 472, 475, 479, 481, 485, 486, 487, 488, 490, 491, 492, 495, 496, 497, 500 *ELSET,ELSET=INNERAIRABSORB_S4 427, 432, 433, 442, 445, 446, 447, 449, 454, 465, 469, 470, 474, 480, 489, 493, 494, 501 *SURFACE,NAME=INNERAIRABSORB INNERAIRABSORB_S3, S3 INNERAIRABSORB_S2, S2 INNERAIRABSORB_S4, S4 INNERAIRABSORB_S1, S1 ** ------------------------------------------- ** ** MATERIALS ** *MATERIAL,NAME=AIR *ACOUSTICMEDIUM,BULKMODULUS 1.42E5, *DENSITY 1.2, *MATERIAL,NAME=AIRABSORB *ACOUSTICMEDIUM,BULKMODULUS 1.42E5, *ACOUSTICMEDIUM,VOLUMETRICDRAG 1.2, *DENSITY 1.2, *MATERIAL,NAME=STEEL *DENSITY 7920., *ELASTIC 1.9E+11, 0.3 ** ------------------------------------------- *TIE,NAME=ACOST1 INNERAIR2SHELL, SHELL2INNERAIR OUTERAIR2SHELL, SHELL2OUTERAIR ** ** BOUNDARY CONDITIONS AND CONSTRAINTS ** ** The transform and boundary replaces the equation ** used to create ZSYMM condition. ** *TRANSFORM,NSET=SHELLCONSTRAINT 0.0344683,0.6,0.00607769,-0.00607769 ,0.6, 0.0344683 *BOUNDARY SHELLCONSTRAINT,ZSYMM **** *BOUNDARY SHELLSYM,ZSYMM ***EQUATION ** 2 **SHELLCONSTRAINT,1,0.176326981,SHELLCONSTRAINT,3,-1. ** 2 **SHELLCONSTRAINT,2,1.,SHELLSYM,2,-1. **** ------------------------------------------- *AMPLITUDE,NAME=SINE,DEF=PERIODIC 1,,0.0,0.0 0.0,1.0,0,0 ** ** STEP: STEP-1 ** *NSET,NSET=QA_TEST CENTERLINE *STEP *DYNAMIC,EXPLICIT , *BOUNDARY,AMP=SINE INLETPRESSURE,8,8,1.0 *Simpedance,non=circular cylinder,0.7 *Simpedance,non=spherical sphere1,0.7 *Simpedance,non=spherical sphere2,0.7 *SIMPEDANCE INNERAIROUTLET, ** *OUTPUT,FIELD,NUM=1 *NODE OUTPUT U,POR *NODE OUTPUT,NSET=QA_TEST POR, ***OUTPUT,HISTORY,FRE=1 ***NODE OUTPUT,NSET=CENTERLINE ** POR, ***NODE OUTPUT,NSET=SHELLDISP ** U, **** *END STEP