*HEADING
 COUPLED ACOUSTIC ANALYSIS OF A MUFFLER, EXPLICIT
 SUBMODEL ANALYSIS--THIS IS THE
 SUBMODEL: INCLUDES THE SHELL STRUCTURE AND 
 THE EXTERIOR AIR
**
** This input file needs the following files:
**  
**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
** -------------------------------------------
** 
** 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=OUTERAIR,ELSET=OUTERAIR-1_I1
** -------------------------------------------
*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
** -------------------------------------------
** 
** 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
OUTERAIR2SHELL,SHELL2OUTERAIR
*** -------------------------------------------
*AMPLITUDE,NAME=SINE,DEF=PERIODIC
1,<omega>,0.0,0.0
0.0,1.0,0,0
** 
** STEP: STEP-1
** 
*NSET,NSET=INTER_P,ELSET=OUTERAIR2SHELL_S4
*NSET,NSET=OUTER_P,ELSET=CYLINDER_S6
*NSET,NSET=TEST_P
INTER_P,OUTER_P
*NSET,NSET=QA_TEST
INTER_P,
*NSET,NSET=shell,ELSET=MUFFLERSHELL-1_I1
*SUBMODEL, ABSOLUTE=1.E-3,GLOBAL ELSET=MUFFLERSHELL-1_I1
Shell,
*STEP
*DYNAMIC,EXPLICIT
,<T>
*Boundary,submodel,step=1
shell,1,3
*Simpedance,nonreflecting=circular
 cylinder,0.7
*Simpedance,nonreflecting=spherical
 sphere1,0.7
*Simpedance,nonreflecting=spherical
 sphere2,0.7
*************
*FILE OUTPUT,NUM=1
*NODE FILE,NSET=INTER_P
POR,
******************
*OUTPUT,FIELD,NUM=4
*NODE OUTPUT
U,POR
*******************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
POR,
*******************
***OUTPUT,HISTORY,FREQ=1
***NODE OUTPUT,NSET=TEST_P
***POR
*END STEP