*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
** -------------------------------------------
*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
** -------------------------------------------
*ELEMENT, TYPE=ACIN3D4, ELSET=INF
101519,4719,4692,3380,3407
101538,4720,4719,3407,3408
101557,4721,4720,3408,3409
101576,4722,4721,3409,3410
101595,4723,4722,3410,3411
101614,4724,4723,3411,3412
101633,4725,4724,3412,3413
101652,4726,4725,3413,3414
101671,4727,4726,3414,3415
101690,4728,4727,3415,3416
101709,4729,4728,3416,3417
101728,4730,4729,3417,3418
101747,4731,4730,3418,3419
101766,4732,4731,3419,3420
101785,4733,4732,3420,3421
101804,4734,4733,3421,3422
101823,4735,4734,3422,3423
101842,4736,4735,3423,3424
101861,4737,4736,3424,3425
101880,4738,4737,3425,3426
101899,4739,4738,3426,3427
101918,4740,4739,3427,3428
101937,4741,4740,3428,3429
101956,4742,4741,3429,3430
101975,4693,4742,3430,3381
102100,4692,4696,3384,3380
102100,4696,4831,3519,3384
102116,4831,4697,3385,3519
102420,4862,4698,3386,3550
102420,4698,4693,3381,3386
102421,4863,4862,3550,3551
102422,4699,4863,3551,3387
102423,3550,4862,4698,3386
102424,3380,4692,4696,3384
*SOLID SECTION, ELSET=INF, MATERIAL=AIR
1.0,0.0,0.0,0.0
*SURFACE, TYPE=ELEMENT, NAME=INFSURF
INF,SPOS
**--------------------------------------------
*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
INFSURF, SHELL2OUTERAIR
*TIE,NAME=ACOST2
INNERAIR2SHELL, SHELL2INNERAIR
** 
** 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,<omega>,0.0,0.0
0.0,1.0,0,0
** 
** STEP: STEP-1
** 
*NSET,NSET=QA_TEST
19,118,119,
*STEP
*DYNAMIC,EXPLICIT
,<T>
*BOUNDARY,AMP=SINE
 INLETPRESSURE,8,8,1.0
*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