*HEADING
 FULLY COUPLED ACOUSTIC ANALYSIS OF A MUFFLER
 M, KG, S
** -------------------------------------------
**
** PART INSTANCE: INTERIORAIR-1
**
*NODE, INPUT=muffler_intair_nodes.inp
*ELEMENT, TYPE=AC3D10,
INPUT=muffler_intair_elem.inp
** -------------------------------------------
*ELSET, ELSET=INTERIORAIR-1_I1, GENERATE
*SOLID SECTION, ELSET=INTERIORAIR-1_I1,
MATERIAL=AIRABSORB
 1.,
** -------------------------------------------
*ELSET, ELSET=INTERIORAIR-1_I2, GENERATE
*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
*SHELL SECTION, ELSET=MUFFLERSHELL-1_I1,
MATERIAL=STEEL
 0.00075, 5
** -------------------------------------------
**
** PART INSTANCE: OUTERAIR-1
**
*NODE,INPUT=muffler_extair_nodes.inp
*ELEMENT, TYPE=AC3D20,
INPUT=muffler_extair_elem.inp
** -------------------------------------------
*ELSET, ELSET=OUTERAIR-1_I1, GENERATE
*SOLID SECTION, ELSET=OUTERAIR-1_I1,
MATERIAL=AIR
 1.,
** -------------------------------------------
*NSET, NSET=SHELLSYM,UNSORTED
** -------------------------------------------
*NSET, NSET=SHELLCONSTRAINT,UNSORTED
** -------------------------------------------
*NSET, NSET=SHELLDISP
** -------------------------------------------
*NSET, NSET=CENTERLINE
** -------------------------------------------
*NSET, NSET=INLETPRESSURE
** -------------------------------------------
*ELSET, ELSET=INNERAIROUTLET_S4
*SURFACE, NAME=INNERAIROUTLET
INNERAIROUTLET_S4, S4
** -------------------------------------------
*ELSET, ELSET=INNERAIR2SHELL_S1
*ELSET, ELSET=INNERAIR2SHELL_S4
*ELSET, ELSET=INNERAIR2SHELL_S2
*ELSET, ELSET=INNERAIR2SHELL_S3
*SURFACE, NAME=INNERAIR2SHELL
INNERAIR2SHELL_S1, S1
INNERAIR2SHELL_S4, S4
INNERAIR2SHELL_S2, S2
INNERAIR2SHELL_S3, S3
** -------------------------------------------
*ELSET, ELSET=SHELL2INNERAIR_SNEG, GENERATE
*SURFACE, NAME=SHELL2INNERAIR
SHELL2INNERAIR_SNEG, SNEG
*ELSET, ELSET=SHELL2OUTERAIR_SPOS, GENERATE
*SURFACE, NAME=SHELL2OUTERAIR
SHELL2OUTERAIR_SPOS, SPOS
** -------------------------------------------
*ELSET, ELSET=OUTERAIR2SHELL_S6
*ELSET, ELSET=OUTERAIR2SHELL_S3
*ELSET, ELSET=OUTERAIR2SHELL_S4, GENERATE
*SURFACE, NAME=OUTERAIR2SHELL
OUTERAIR2SHELL_S6, S6
OUTERAIR2SHELL_S3, S3
OUTERAIR2SHELL_S4, S4
** -------------------------------------------
*ELSET, ELSET=SPHERE1_S3, GENERATE
*ELSET, ELSET=SPHERE1_S4, GENERATE
*SURFACE, NAME=SPHERE1
SPHERE1_S3, S3
SPHERE1_S4, S4
** -------------------------------------------
*ELSET, ELSET=SPHERE2_S3, GENERATE
*ELSET, ELSET=SPHERE2_S4, GENERATE
*SURFACE, NAME=SPHERE2
SPHERE2_S3, S3
SPHERE2_S4, S4
** -------------------------------------------
*ELSET, ELSET=CYLINDER_S6, GENERATE
*SURFACE, NAME=CYLINDER
CYLINDER_S6, S6
** -------------------------------------------
*ELSET, ELSET=INNERAIRABSORB_S3
*ELSET, ELSET=INNERAIRABSORB_S2
*ELSET, ELSET=INNERAIRABSORB_S1
*ELSET, ELSET=INNERAIRABSORB_S4
*SURFACE, NAME=INNERAIRABSORB
INNERAIRABSORB_S3, S3
INNERAIRABSORB_S2, S2
INNERAIRABSORB_S4, S4
INNERAIRABSORB_S1, S1
** -------------------------------------------
**
** MATERIALS
**
*MATERIAL, NAME=AIR
*ACOUSTICMEDIUM,BULKMODULUS
*DENSITY
 1.2,
*MATERIAL, NAME=AIRABSORB
*ACOUSTICMEDIUM,BULKMODULUS
*ACOUSTICMEDIUM,VOLUMETRICDRAG
*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
**
*BOUNDARY
 SHELLSYM,ZSYMM
*EQUATION
 SHELLCONSTRAINT,1,0.176326981,
SHELLCONSTRAINT,3,-1.
 SHELLCONSTRAINT,2,1.,  SHELLSYM,2,-1.
** -------------------------------------------
**
** STEP: STEP-1
**
*STEP
 STEP-1: FREQUENCY RESPONSE
*STEADY STATE DYNAMICS, DIRECT,
FREQUENCY SCALE=LINEAR
 179.0,181.0,12
*BOUNDARY
 INLETPRESSURE,8,8,1.0
*Simpedance,non=circ
 cylinder,0.7
*Simpedance,non=spher
 sphere1,0.7
 sphere2,0.7
*SIMPEDANCE,non=plan
 INNERAIROUTLET,
*OUTPUT,FIELD,OP=NEW,FREQUENCY=6
*NODE OUTPUT
 U,POR
*OUTPUT,HISTORY,OP=NEW,FREQUENCY=1
*NODE OUTPUT,NSET=CENTERLINE
 POR,
*NODE OUTPUT,NSET=SHELLDISP
 U,
*NODE FILE,NSET=CENTERLINE,FREQUENCY=1
 COORD,POR,PPOR
*NODE FILE,NSET=SHELLDISP,FREQUENCY=1
 U,PU
*END STEP