Optimization of Sound Transmission Loss through a Thin Functionally Graded Material Cylindrical Shell

Themaximizing of sound transmission loss (TL) across a functionally graded material (FGM) cylindrical shell has been conducted using a genetic algorithm (GA). To prevent the softening effect from occurring due to optimization, the objective function is modified based on the first resonant frequency. Optimization is performed over the frequency range 1000–4000Hz, where the ear is the most sensitive. The weighting constants are chosen here to correspond to an A-weighting scale. Since the weight of the shell structure is an important concern inmost applications, the weight of the optimized structure is constrained. Several traditional materials are used and the result shows that optimized shells with aluminum-nickel and aluminum-steel FGMare themost effective at maximizing TL at both stiffness and mass control region, while they have minimum weight.


Introduction
Functionally graded material (FGM) is a multiphase material consisting of different material components, such as ceramics and metals that have various mixture ratios and microstructures.It was first produced in Japan in the mid-1980s by Yamanouchi et al. [1] and Koizumi [2].FGM is a class of composite materials with continuous designed variation of properties throughout the volume and thus can alleviate the stress concentration found in laminated composites and can significantly enhance the thermal and mechanical properties of FGM.Furthermore, FGM can be designed to fulfill particular requirements, such as enhancing stiffness, toughness, and resistance to corrosion, wear, and high temperature by utilizing materials or material systems with various properties.Using powder metallurgy, a mixture of ceramic and metal with continuously varying volume fractions can be easily manufactured.Yamaoka et al. [3] and Zhu et al. [4] described various techniques used for the fabrication of FGMs.Consequently, in the last two decades, FGM has been developed for general use as structural components in extremely high temperature environments such as rocket engine components, military armour, space plan body, engine components, thermal barrier coating for turbine blades, and other applications [5].
Thin circular cylindrical shells are widely used in many engineering fields ranging from civil, naval, nuclear, mechanical, chemical, and aerospace applications [6].Vibroacoustic behavior of FGM cylindrical shell is considered as an important subject for researchers.Williamson et al. [7] carried out finite element analysis and studied the residual stresses and strains in ceramic-metal joints.Moya [8] concluded that FGMs are promising candidates for future composites in aerospace, fast computers, environmental sensors, and so forth.Using a micromechanics model, Gasik and Lilius [9] calculated mechanical and thermal properties of W-Cu FGM.Sasaki and Hirai [10] carried out experiments to estimate the fatigue resistance of a Si-graphite FGM.Other researchers [11][12][13][14][15] carried out thermal stress analyses of FGMs under various thermal loading conditions.
A cylindrical shell structure made up of a FGM is one of the basic structural configurations where important studies have been done.Loy et al. [16,17] carried out analysis with strain-displacement relations from Love's shell theory and the eigenvalue governing equation is obtained using the Rayleigh-Ritz method.The results demonstrated that the frequency characteristics are similar to those observed for homogeneous isotropic cylindrical shells and the frequencies are affected by the constituent volume fractions and the configurations of the constituent materials.Noise transmission through a composite cylindrical shell was studied by many researchers.Tang et al. [18] considered an infinite cylindrical sandwich shell excited by an oblique plane sound wave with two independent incident angles.Daneshjou et al. [19] presented an analytical solution for acoustic transmission through relatively thick FGM cylindrical shells using third order shear deformation theory.
Optimization studies of transmission loss (TL) in cylindrical shell structures with respect to their material and geometric properties have been an important topic of research for a few decades.Most of these studies have been related to composite materials.Makris et al. [20] studied transmission loss optimization in acoustic sandwich panels.The results indicated that the optimum TL is typically found at the value of the core density corresponding to the lowest of the permissible range of core densities established at the start of the optimization search.Wang et al. [21] used a genetic algorithm to optimize a sandwich panel with balanced acoustic and mechanical properties at a minimal weight.The optimization of sound transmission through FGM cylindrical shells that were subjected to plane sound wave has not been studied yet.In the present work, an attempt is made to address this problem.Therefore, it is very important to develop an accurate, reliable analysis towards the understanding of the optimization sound transmission through FGM cylindrical shells.Therefore, this paper presents a study on the optimization of sound transmission loss across FGM cylindrical shell using a genetic algorithm.A parametric study of sound transmission through a FGM shell excited by an incident oblique plane sound wave was developed and considered as objective function to determine which combinations of design variables could result in a FGM shell which maximizes the sound insulation.The result shows that an optimized shell with aluminum-nickel or aluminum-steel FG has the largest TL at both stiffness and mass control region while having minimum weight.

Model Specification
The specific problem is an oblique plane wave impinging on a flexible FGM thin cylindrical shell of infinite length and includes the reflection and scattering of the incident wave and the effect of an external airflow in the x-direction (Figure 1).The angle of incidence of the wave is  (measured from the axial axis of cylinder), and the incidence wave approaches from the  =  direction.The inside cavity was assumed to be anechoic and only inward-traveling waves exist.In the analysis, all waves will be assumed to have the same dependence on the axial coordinate  and the cylinder will be assumed infinitely long [22].The fluid media in the external and the internal space are defined by the density and the speed of sound: { 1 ,  1 } and { 3 ,  3 }, respectively.

Functionally Graded Material.
In order to model FG properties for a cylindrical shell composed of  different materials with a uniform thickness h and a middle reference surface, the volume fraction relation   changes continuously through the thickness direction, considered as [16] where  is the power-law exponent (0 ≤  ≤ ∞).The general formulation for a FGM composed of two components can be set as [17]  = ( where  can be substituted with , , or ] in (2) to obtain the elasticity modulus, mass density, and poisson's ratio for the FGM cylindrical shell, respectively.From these equations, when  = ℎ/2,  =  2 ,  =  2 , ] = ] 2 , and when  = −ℎ/2,  =  1 ,  =  1 , ] = ] 1 .The material properties vary continuously from material 2 at the inner surface of the cylindrical shell to material 1 at the outer surface of the cylindrical shell.

Fluid Equations.
Due to the existence of airflow in the external fluid medium, the external pressure, which is the summation of the incident wave   1 and the reflected wave   1 , satisfies the following wave equation [23]: As a result of considering Helmholtz equations in internal and external space and also the equality of particle velocities of the acoustic media and the shell, the boundary conditions of the model can be defined as 2.3.Formulation.Consider a cylindrical shell with radius  and thickness ℎ; see Figure 1.The deformations defined with reference to a coordinate system, taken at the middle surface, are , V, and  in the axial, circumferential, and radial directions, respectively.For a thin cylindrical shell, plane stress condition is assumed and the constitutive relation is given by The stress vector, strain vector, and the reduced stiffness matrix are defined as For isotropic materials, the reduced stiffness   (,  = 1, 2 and 6) is defined as From Love's shell theory, the components in the strain vector are defined as [23,24] where  1 ,  2 , and  are the reference surface strains and  1 ,  2 , and  are the surface curvatures.These surface strains and curvatures are defined as For a thin cylindrical shell, the force and moment resultants are defined as Substituting (7), with substitution from (10), into (12), the constitutive equation is obtained as where   ,   , and   (,  = 1, 2 and 6) are the extensional, coupling, and bending stiffness defined as According to [23,24], a thin cylindrical shell's equations of motion in cylindrical coordinate are as follows: where   ,   , and   are loads in axial, circumferential, and radial direction, respectively.

Vibroacoustic Solution.
The harmonic plane wave   in cylindrical coordinate incident from outside of the shell to the direction shown in Figure 1 can be expressed as [25]   (, , , ) where  1 is the wave number in the external medium of shell and   is the Bessel function of the first kind of integer order ,   is the Neumann factor,  0 is the amplitude of the incident wave,  = −1,  = 0, 1, 2, 3, . .., and  is the angular frequency.The waves radiated from the shell to the outside and into the cavity,   1 and   3 , can be represented as where  2  ,  1  are the Hankel functions of the first and second kind of integer order , respectively.Three components of the shell displacements can be expressed as [24]  (, , ) = Because the traveling waves in the acoustic media and inside the shell are driven by the incident traveling wave, the wave numbers (or trace velocities) in the  direction should match throughout the system; therefore,  3 =  1 .
Equations ( 16)-( 18) are substituted into the three equations of motion and two boundary conditions (( 5)-( 6)).These five equations involve six variables: the amplitudes of the outgoing and incoming waves in the exterior cavity, the transmitted wave in the interior cavity, and three displacements of the shell structure.Therefore, the solutions can be obtained as the ratios to the one of the variables, in this case, the pressure amplitude of the incoming wave [18].

Sound Transmission Loss. Sound transmission loss (STL)
can be defined as the ratio of the incoming and transmitted sound powers per unit length of the cylinder.Consider STL = 10 log 10 where   and   are the transmitted and incident power flow per unit length of the shell.The equations for   and   are developed by Kim and Lee [25].However, because of the incident angle in the present studies, the equation of Kim and Lee's work should be changed to Therefore, the final quantity of transmission loss can be defined as where

The Optimization Problem
To optimally design a structure so that its STL must be minimal, we have to deal with a multivariable constrained nonlinear optimization problem.The optimization problem can be stated as where the   are  equality constraints and the   are (pm) inequality constraints.() is the objective function to be minimized with respect to the vector  of the design variables.In the present paper, we study the optimization problem of the shell with respect to the geometric and material properties.The design variables for this problem are chosen as where , ℎ, and  are radius, thickness (geometric design variable), and power-law exponent, respectively. lb and  ub are vectors which define the lower and upper bounds of the design variables.Since the weight of the shell structure is an important concern in most applications, a constraint is introduced as follows: where  0 is a constant and refers to the mass per unit length of shell.
A large number of optimization techniques exist for different problems.In this study, the genetic algorithm is used as optimization method.The genetic algorithm is a method for solving both constrained and unconstrained optimization problems that is based on natural selection, the process that drives biological evolution.It repeatedly modifies a population of individual solutions.At each step, the genetic algorithm selects individuals at random from the current population to be parents and uses them to produce the children for the next generation.Over successive generations, the population "evolves" toward an optimal solution.The genetic algorithm uses three main types of rules at each step to create the next generation from the current population: selection rules select the individuals, called parents that contribute to the population at the next generation.Crossover rules combine two parents to form children for the next generation.Mutation rules apply random changes to individual parents to form children.

Result and Discussion
Numerical simulations have been carried out to study the optimization problem of a FGM cylindrical shell.With proper qualitative interpretations, the theoretical model developed can be used very effectively in the basic design stage of cylinder-shape vibroacoustic systems.As a demonstration of such applications, design parameter studies are conducted.The basic shell dimensions and simulation conditions used in the study are listed in Table 1.
The TL of a FGM Brass-Al shell is illustrated in Figure 2. As shown, three important frequencies are recognized: ring frequency   (where the wavelength of a longitudinal wave in the shell is equal to the circumference), critical pseudocoincidence  pc (spatial coincidence in radial direction between the wave vector projection of excitation and the shell circumferential wave number), and coincidence frequency   (where the trace velocity of the acoustic wave is equal to the bending wave velocity in the shell wall).The various zones controlled by the stiffness (0 −   ), the mass (  − pc ), and the coincidence ( >   ) and the resonant modes zone enclosed between  pc and   will be noted [26].
Figure 3 shows the effect of the shell materials on the STL.Materials chosen for the comparison are aluminum, steel, nickel, zirconia, and brass as shown in Table 1.The figure shows that the brass is the most effective in the high frequency range and zirconia is the most effective in the low frequency range.This is as expected because the density of the brass and the stiffness of the zirconia are the largest, which make them most effective in the mass and stiffness controlled regions, respectively.The figure also shows that the aluminum, which has the lowest stiffness, is the least effective in the low frequency range, which is again as expected because the low frequency range is controlled by the stiffness.Steel and nickel have the same behavior in the low and high frequency range and are more effective in both mass and stiffness controlled regions.
Although steel and nickel are more effective, the structural weight of a cylindrical shell that is made of them is too high.Since the weight of the shell structure is an important concern in most applications, it is necessary to compromise between sound transmission loss of cylindrical shell and its weight.The results, however, lead to a structure with a lower stiffness than the base shell.To prevent the softening  effect from occurring due to optimization, we introduce an additional term in the objective function: where () is the first resonant frequency of the optimized shell,  1 = 7Hz is the first resonant frequency of the base shell, and  is a weighting constant.Note that the first resonant frequency of a structure is a common measure of the stiffness of the structure.This modified objective function will penalize the reduction of the stiffness of the structure.The optimization is then performed at octave band frequency under 1000 Hz.The results evaluate an equal weighting of  = 0, 0.4, 0.5, 0.6 given to both the sound transmission loss and the stiffness of the structure.The optimization results for other values of  are similar and are not presented.Table 2 shows the results of the optimized shell at the frequency.In this example, stainless steel and aluminum are considered as the constituent materials of the FGM cylindrical shell.
The results indicate that the optimum values are typically found at the value of  = 0.5, so it is used for other analyses in this study.The results shown in Table 2 are the shell optimized at discrete frequency.It is also common to optimize a frequency-weighted average of the transmission loss.For practical purposes hearing encompasses a range of frequencies from about 16 Hz to somewhat less than 20000 Hz.However, low frequencies do not affect humans very strongly, while those from 500-5000 Hz (where the ear is most sensitive) are very important.
This study uses the frequency range 1-4 KHz.However, instead of the 21 frequencies (150 Hz steps) used then, only the seven standardized frequencies from the STC rating in the 1000 to 4000 Hz range will be incorporated.Namely, the following frequencies are considered: 1000, 1250, 1600, 2000, 2500, 3150, and 4000 Hz.The smaller number of frequencies does not lessen the accuracy in the calculations of TL because a sophisticated integration algorithm was used to find the field-incidence average TL.Also, the use of these frequencies parallels the portion of the spectrum emphasized in the recent ASTM standard on isolation between neighboring rooms [25].The lower limit of the range (1000 Hz) was chosen where   represent weighting constants and   is the number of frequencies [29].  are normalized so that the sum of all the coefficients is unity.The weighting constants   are commonly chosen to correspond to an A-weighting.It should be mentioned here that with A-weighting, the sound-level meter is relatively less sensitive to sounds of frequency below 1000 and above 4000 Hz.The weighting and discrete frequencies that were chosen in the bound 1-4 kHz are identical to those used to defined single insulation measure in an ASTM standard [27,28].The normalized weights   are given in Table 3.
In this study, based on the materials listed in Table 1, different functionally graded cylindrical shells (Table 4) should be optimized based on (27) and ( 24)- (25).
As the weight of structure is an important parameter in aerospace vehicles, aluminum is used as inner surface of constituent materials because its density is lowest among the metals in Table 1.The optimized results of the all cases are tabulated in Table 5.
The result shows that case 2 is the most effective for the STL and case 3 is worse than the other cases.Figure 4 shows the variation of the STL for the all cases with optimized design parameters of Table 5.
As indicated, the results are different for the optimized shells with different kinds of FGM, although STL improvement is achieved in all cases.As obviously illustrated in Figure 4, TLs characteristics have been significantly improved    over frequency of 1000 Hz, which have been weighted according to Table 3. Figure 5 compares the TLs of FGM shells made from different materials with optimized design parameters found by the GA.Since the aluminum-brass FGM shell was the heaviest, its optimized curve is the highest on the mass control region but its behavior in the stiffness control is the worst.The figure shows that case 1 and case 2 are the most effective in the low and high frequencies range and so they should be used as the first selection material tabulated in Table 4 with optimized value by (26) and (24).Figure 6 shows the variations of the volume fractions in the thickness direction  for all cases of FG cylindrical shell.
In Figure 6, the material properties on the inner surface of the all optimized cylindrical shell are those of aluminum and on the outer surface are those of as in Table 3.As  increases  = −0.5, the amount of aluminum constituent for the optimized case 1 and case 3 is the highest and the lowest compared to other optimized cases, respectively.This means that the weight constraint of optimized shell at case 1 is better than other cases, while maximum STL is achieved.

Concluding Remarks
The optimization study of sound transmission across a FGM cylindrical shell has been presented in this paper.The main objective of the optimization is to maximize the sound transmission loss across the structure subjected to weight and the thickness constraint.In order to prevent the softening effect occurrence due to optimization, an additional term is added to the objective function based on the first resonant frequency.Several traditional materials have been chosen for investigation.The power-law exponent, shell radius, thickness of shell, and first resonant frequency have been selected as design variables.Optimization of transmission loss is performed over the frequency range 1000-4000 Hz, where the ear is the most sensitive.The results show that the optimized shell made of aluminum-nickel and aluminum-steel is more effective on the in a wide frequency range, while aluminum-brass and aluminum-zircon are only effective in mass control region and stiffness control region, respectively.Based on the result of optimized shell geometries and its volume fraction, minimum weight design with best STL is found for aluminum-nickel and aluminum-steel functionally graded materials.

Figure 4 :
Figure 4: TL as a function of frequency  for the FGM shells listed in Table4.

Figure 5 :Figure 6 :
Figure 5: Comparison between TL-f curves for the FGM shells listed in Table4.

Table 1 :
Material and environmental properties.

Table 2 :
Effect of  values on the design parameters of the optimized shell.

Table 3 :
Normalized weights using the A-weighting.

Table 5 :
Optimized value of FGM combinations.