*Heading
 Exterior air model for sequentially 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,nset=shell,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
** ----------------------------------------------------------------
*Nset, nset=shellsym,unsorted
** ----------------------------------------------------------------
*Nset, nset=shellconstraint,unsorted
** ----------------------------------------------------------------
*Nset, nset=shelldisp
** ----------------------------------------------------------------
*Nset, nset=centerline
** ----------------------------------------------------------------
*Nset, nset=inletpressure
** ----------------------------------------------------------------
*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=_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
**
** 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
 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=shell,frequency=1
 U,PU
*End step