Sound Transmission Loss of a Honeycomb Sandwich Cylindrical Shell with Functionally Graded Porous Layers

: To examine the acousto-structural behavior of a sandwich cylindrical shell beneﬁting from hexagonal honeycomb structures in its core and functionally graded porous (FGP) layers on its outer and inner surfaces, a comprehensive study based on an analytical model which also considers the effect of an external ﬂow is conducted. A homogenous orthotropic model is used for the honeycomb core while its corresponding material features are found from the modiﬁed Gibson’s equation. The distribution pattern of FGP parts is either even or logarithmic-uneven, and a special rule-of-mixture relation governs their properties. Based on the ﬁrst-order shear deformation theory (FSDT), Hamilton’s principle is exploited to derive the ﬁnal coupled vibro-acoustic equations, which are then solved analytically to allow us to calculate the amount of sound transmission loss (STL) through the whole structure. This acoustic property is further investigated in the frequency domain by changing a set of parameters, i.e., Mach number, wave approach angle, structure’s radius, volume fraction, index of functionally graded material (FGM), and different honeycomb properties. Overall, good agreement is observed between the result of the present study and previous ﬁndings.


Introduction
Sandwich structures are multilayer physical members composed of two stiff face sheets bonded to a core layer.Core layers usually consist of metallic and non-metallic honeycomb core, open and closed cell foams, and balsa wood [1][2][3][4][5].One of the promising ways to reduce material cost and weight is to utilize honeycomb structures of various shapes.Most frequently, such structures are employed instead of or in combination with plates and shells, especially as a lightweight core.There have been several applications of honeycomb materials in various engineering fields, corroborating their numerous benefits including high specific strength/stiffness, incredibly lower weight, energy damping, and proper stability.However, their vibration characteristics, especially in combination with shells and plates, still need further evaluation.Li and Jin [6] investigated the natural frequencies of symmetric rectangular honeycomb panels under clamped and simply supported boundary conditions, applying corrected Gibson formula and third-order shear deformation theory (TSDT).Li et al. [7] analyzed geometrically nonlinear free vibration of symmetric honeycomb sandwich panels considering simply supported boundary condition utilizing homotopy analysis technique and TSDT.Duc et al. [8] analytically studied nonlinear dynamic behavior of sandwich composite cylindrical panels with auxetic honeycomb core layer subjected to blast, damping, and mechanical loads.Duc et al. [9] employed FSDT to calculate the dynamic response of composite double curved shallow shells with negative Poisson's ratios in auxetic honeycombs core layer.Li et al. [10] predicated the natural frequencies of a thin-walled hexagonal honeycomb sandwich cylindrical shell according to Flügge's shell theory, and used finite element method (FEM) and experimental results to validate their models.Eipakchi et al. [11] presented a mathematical process to calculate the critical velocities, dynamic response, and the natural frequencies of composite cylindrical shells with auxetic honeycombs core layer under a moving pressure.Quyen et al. [12] carried out the nonlinear dynamic response of a honeycomb sandwich cylindrical panel on visco-Pasternak substrates under thermal and blast loading.
A shell structure, usually defined as a thin, curved plate, is a generally lightweight element that can be used in different sizes and applications.The main body of nearly all vehicles is built from such shell elements, and these structures have become an indispensable part of other advanced technologies such as reactors and gigantic telescopes.Unsurprisingly, they undergo large loads in different environments and scenarios, especially when placed inside a fluid or when facing aerodynamic loads.The resulting vibrations can specifically create a strong noise in the surroundings, and are thus required to be fully investigated under such circumstances [13][14][15][16][17].In this regard, forced and free vibrations of cylindrical shells and mechanical properties of polymer based nanocomposite have been thoroughly investigated in the past, as in the works of [18][19][20][21][22][23][24].
Functionally graded materials (FGMs) are generally characterized by the gradual changes of composition and structure over volume, thus offering different properties on different sides of the specimen [25][26][27][28][29][30].They have shown promising results in mitigating common issues found in conventional laminated composites such as interfacial debonding and matrix cracking while finding their way into different industries.These materials can be employed in different arrangements, one of which is an FGM sandwich structure formed around a core layer and covered by external/internal faces.This straightforward configuration is easy to manufacture and offers various improvements over traditional FGMs, including better thermal protection and sound absorption.Nevertheless, accurate models are necessary to analyze them in different applications [31][32][33].Loy et al. [34] studied the effect of constituent volume fractions of FGMs on the natural frequencies of a cylindrical shell using Love's shell theory and Rayleigh-Ritz method.Pradhan et al. [35] used volume fraction power law distribution to analyze the free vibration of FGM cylindrical shells under various boundary conditions.Based on the FSDT in conjunction with the wave based method (WBM), Liu et al. [36] investigated the effect of material gradient index on the natural frequencies of FGM cylindrical shells.The influence of the nonlinear temperature variations and porosity volume fraction, and constituent volume on the natural frequencies of FGM cylindrical shells are presented by Wang et al. [37].Baghlani et al. [38] developed and suggested a semi-analytical method to compute the natural frequency of eccentrically stiffened FGM shells considering fluid-structure interaction and higherorder shear deformation theory.Sofiyev [39] investigated the dynamic behavior of the infinitely-long FGM cylindrical shells under moving loads.Baghbadorani and Kiani [40] performed research on the behavior of a cylindrical shell made of a composite material reinforced with graphene platelets employing FSDT and Halpin-Tsai micro-mechanical.Nguyen et al. [41] proposed a semi analytical method to study the nonlinear asymmetric vibration of corrugated sandwich FGM cylindrical shells containing fluid subjected to harmonic radial load.Ni et al. [42] investigated the effects of magneto-electro-thermal loadings and material properties on the free vibration characteristics of FGM cylindrical shells.Cong and Duc [43] used the Galerkin method and FSDT to predict the nonlinear thermo-buckling and post-buckling behavior of eccentrically stiffened FGM double curved shallow auxetic honeycomb sandwich shells.
Another aspect of FGMs worthy of mention is the formation of voids owing to the differences in solidification temperature when sintering them, as discovered by Rodrıguez-Castro [44] through a series scanning electron microscopy studies.The porosity level greatly affects the mechanical behavior of FGMs, and should thus be precisely characterized to be able to predict the dynamic response of FGM-incorporated structures.Wang [45] performed the free vibration of a FGM porous cylindrical shell under different sets of immovable boundary conditions using sinusoidal shear deformation theory.Cuong-Le et al. [46] studied the elastic instability and free vibration of annular plate, cylindrical and conical shell made of FGM porous rock materials using isogeometric analysis (IGA).They assumed that according to a power-law model, the remarkable material properties are associated with the porosity volume fraction and are considered to be constantly changeable through the thickness direction.Shahgholian et al. [47] presented buckling behavior a porous nanocomposite cylindrical shell reinforced with graphene platelets (GPLs) based on the FSDT and Halpin-Tsai micromechanics approach.Ghadiri and SafarPour [48] used modified couple stress theory in conjunction with FSDT to investigate the free vibration analysis of a FGM porous cylindrical microshell under a thermal environment.Keleshteri and Jelovica [49] denoted that both porosity coefficient and porosity distribution have a notable effect on the nonlinear natural frequencies of FG porous cylindrical panels using FSDT and generalized differential quadrature method (GDQM).Li et al. [50] investigated the thermal vibration problem of an FG porous stepped cylindrical shell applying characteristic orthogonal polynomials and FSDT.
Sound transmission through different objects in another interesting area of study for vibro-acoustic engineers [51].More specifically, STL determines the amount of sound (usually in decibels/dB) isolated by a certain structure, and is considered to be an important quantity when acoustic effects are not negligible.Generally, a plethora of analytical equations have been proposed to obtain the sound transmission characteristics of most common structures, such as fuselages and building walls as in the study by Heckl [52].A major hindrance, however, to implement these equations is to find their exact solutions, hence giving rise to various approximate approaches over the course of time.Pellicier and Trompette [53] provided a comprehensive review of such techniques.Ramezani and Talebitooti [54] used the three different distribution models, including power-law, sigmoid, and exponential types, to predict STL across a double-walled FGM porous sandwich cylindrical shell.Daneshjou et al. [55] presented STL of a thick-walled cylindrical shell based on three dimensional (3-D) theory of elasticity.Oliazadeh and Farshidianfar [56] applied Donnell's shell theory to analyze STL through double-and triple-walled cylindrical shells.Ahmadi et al. [57] analyzed the effect of different models of carbon nanotube distribution on the variation of acoustic transmission of FGM carbon nanotube-reinforced composite cylindrical shells using FSDT.Hasheminejad et al. [58] performed STL through a sandwich cylindrical shell with electrorheological fluid core.Fu et al. [59] indicated the effects of gradient index, porosity volume fraction, and porosity distribution type on the variations of STL of a FG porous cylindrical shell under nonlinear thermal loading.Hasheminejad and Jamalpoor [60] improved the diffuse STL across a smart hybrid double concentric sandwich circular cylindrical shell structure under internal and external air gap mean flows using multi-input multi-output (MIMO) sliding mode control (SMC).According to the hyperbolic tangent shear deformation assumption, Li et al. [61] calculated the natural frequencies and STL of FG honeycomb sandwich plates.
To the best of authors' knowledge, no study has considered the sound transmission and wave propagation of FGM-incorporated honeycomb sandwich cylindrical shells until now.To fill this gap, this paper examines the behavior of an incident sound wave passing through a honeycomb sandwich shell filled with air and benefiting from FGP layers.For a more realistic evaluation, the whole structure is assumed to be located in fluid flow, and the structural response is analyzed via the FSDT.A power law model is necessary to describe the varying material properties caused by two porosity patterns (i.e., even distribution and logarithmic-uneven pattern).For the simulation part of this study, the honeycomb core is represented by a homogenous orthotropic material whose different features are explained via a modified Gibson's equation.
The outline of the paper is as follows.Section 2 expresses the acoustic equation of fluids and the governing equations, whereas in Section 3 the STL is defined.Finally, comparisons and parametric results to study the acousto-structural behavior of a honeycomb sandwich cylindrical shell with FGP layers are proposed in Section 4.

Mathematical Model
Figure 1 shows the schematic of the problem at hand.The honeycomb core of the sandwich cylindrical shell has a radius of R and thickness of H.The hexagonal cells are characterized by t c , l c , and θ c .The exterior and the interior of this core are enclosed by porous layers of FG materials, each having a thickness of h.The external FG layer also faces a steady fluid flow, which is air in this case.Finally, the hollow section of the sandwich structure houses air, with the characteristic impedance of (ρ 0 , c 0 ).In addition, a steady flow of air passes over the structure at the velocity V.
the structural response is analyzed via the FSDT.A power law model is necessary to describe the varying material properties caused by two porosity patterns (i.e., even distribution and logarithmic-uneven pattern).For the simulation part of this study, the honeycomb core is represented by a homogenous orthotropic material whose different features are explained via a modified Gibson's equation.
The outline of the paper is as follows.Section 2 expresses the acoustic equation of fluids and the governing equations, whereas in Section 3 the STL is defined.Finally, comparisons and parametric results to study the acousto-structural behavior of a honeycomb sandwich cylindrical shell with FGP layers are proposed in Section 4.

Mathematical Model
Figure 1 shows the schematic of the problem at hand.The honeycomb core of the sandwich cylindrical shell has a radius of  and thickness of .The hexagonal cells are characterized by  ,  , and  .The exterior and the interior of this core are enclosed by porous layers of FG materials, each having a thickness of ℎ.The external FG layer also faces a steady fluid flow, which is air in this case.Finally, the hollow section of the sandwich structure houses air, with the characteristic impedance of ( ,  ).In addition, a steady flow of air passes over the structure at the velocity .

Acoustic Equation of Fluids
The wave equation in the external fluid medium (incident region) considering an airflow with constant velocity () can be expressed as [62]  ∇ ( +  ) + + .∇ ( +  ) = 0, where ∇ =  + + refers to the Laplacian operator in the cylindrical coordinate system.Furthermor  and  are the acoustic pressures associated with the

Acoustic Equation of Fluids
The wave equation in the external fluid medium (incident region) considering an airflow with constant velocity (V) can be expressed as [62] where ∂x 2 refers to the Laplacian operator in the cylindrical coordinate system.Furthermor p I 1 and p R 1 are the acoustic pressures associated with the incident and reflected waves, respectively.The acoustic wave equation in the internal fluid medium (transmitted region) is indicated as where p T 2 decsribes the transmitted wave.For cylindrical coordinate system, the terms of time-harmonic pressure waves are defined as [63] p It should be noted that , p 0 expresses the amplitude of pressure of the incident wave, Ω displays the angular frequency, i = √ −1, J n signifies the cylindrical Bessel function of the first kind and n-th order, and n and H (2) n indicate, respectively, the cylindrical Hankel functions of the first and second kinds.Also, P R 1n , P T 2n , P R 2n , P T 3n are unknown complex coefficients.Furthermore, the radial and axial components of the wavenumbers are stated as in which M = V /c 0 refers to the Mach number of the external flow.

Equations of Motion
Despite the generally acceptable results obtained from FSDT, it poses certain difficulties for some thick laminates in sandwich configurations whose transverse shear modulus is not large enough [64,65].As varying transverse shear strains are not inherently implemented in this approach, it is assumed that virtual transverse shear strains are presented on the surfaces of sandwich structures.To express the displacement fields under such assumptions, one can write [66] u where U and V describe the in-plane deflections of the mid-surface along x and θ directions, respectively.Also, the transverse deflection of the structure is expressed with W. Furthermore, the rotation angles of the middle plane along θ and x directions are signified by φ θ and φ x , respectively.The displacement and rotation terms are stated as [67,68] U, W, in which U, φ x , V, φ θ , W denote the unknown modal factors.However, the strain components in terms of the curvature and mid-surface strain of the sandwich shell are written as follows in which normal and shear strains are signified by (ε xx , ε θθ ) and (γ xθ , γ xz , γ θz ), respec- tively.The stress-strain relation for inner and outer FGP layers can be expressed as where the superscript T and B are the top and bottom layers, respectively.Furthermore, E and ϑ refer to the Young's modulus and Poisson's ratio, respectively.To achieve FG-like properties, the bottom and top surfaces of the sandwich structure are purely made from metal and ceramic, respectively.As already described, two porosity distribution patterns, namely even and logarithmic-uneven, were considered for this study.Accordingly, the corresponding values of elastic modulus, mass density, and Poisson's ratio based on the rule of mixture are described as [69] Even porosity dispersion: Logarithmic uneven porosity dispersion: where m and c signify metal and ceramic phases, respectively.Also, ρ is the mass density.The always-positive gradient index (p) is used in this study to determine the changes of a specific property in the thickness direction.The greater the gradient index, the more metallic the structure.Furthermore, ξ expresses the porosity coefficient.The classical constituent relations between stress and strain tensors for honeycomb core layer can be presented as: where G is the shear modulus.It should be noted that based on the corrected Gibson's formula, the corresponding material properties of honeycomb core layer are defined as [9,10]: The Hamilton's principle is applied to obtain the equilibrium equations of motion in the form [70]: where strain energy, the work done by external forces (the work applied by the incidence sound wave), and the kinetic energy are denoted with X s , X f , and X K , respectively.The variation of kinetic energy of the sandwich honeycomb FGP shell is expressed as
where A is the cross-sectional area, and: The variation of strain energy is presented as: The forces are exerted to the cylinder stem from three sources: (1) the incident sound wave, (2) the reflected external sound pressure, and (3) the transmitted acoustic pressure inside the cylindrical structure.In view of these three forces, one can calculate the work by external loads as: By substituting Equations ( 14), (16), and ( 17) into ( 13) and performing some manipulations, the equations of motion for the sandwich honeycomb FGP shell can be derived as: in which where k s is the shear correction coefficient, which depends on the material and geometric properties, boundary conditions, and forces.This study uses a shear correction coefficient of 5/6 as a reliable estimate, although it should be kept in mind that at small scales, this value turns out to be inaccurate.Lastly, replacing Equation ( 19) (considering relations (8) and ( 11)) into Equation ( 18), the equilibrium equations in terms of displacement and rotations are obtained and provided in Appendix A.

Continuity Conditions of Fluid/Structure
There is no need to consider a boundary condition along the z axis due to the infinite nature of the structure in this specific direction.The structural and acoustic domains are coupled in the following manner in the r direction [63,67]: By replacing Equations (3) and (6) into Equations (A1)-(A5), and Equation (20), after some manipulation, the equilibrium equations in a 7 × 7 matrix format can be obtained as: where parameters Y i,j , f 3 , and f 6 are expressed in Appendix B.

STL
Dividing the transmitted to the incident acoustic powers (i.e., Π tr /Π inc ) yields the transmission coefficient τ.The value of sound transmission loss is then obtained as in which Re refers to the real part of the argument and the superscript "*" signifies the complex conjugate.

Numerical Results and Discussion
Before presenting the chief findings, a series of simplified cases are considered to demonstrate the validity of the developed procedure.

Mode Convergence
To show the effect of the number of modes on the convergence of results, a special case for the incoming wave angle of α = π/4 is considered.According to Figure 2, satisfactory convergence at higher excitation frequencies necessitates a greater mode number.Different assumptions used for this investigation are presented in Table 1.Furthermore, hexagonal cells are assumed to be metal.

Comparison Study
A series of comparisons are made here to show whether the developed procedure is accurate enough.As a first verification study, acoustic effects and honeycomb core layer were assumed to be nonexistent, then the natural frequencies (Hz) of the FGM cylindrical shell were calculated for different values of material gradient index based on the present formulation.The results alongside other numerical findings in Ref. [34] are provided in Table 2, demonstrating the acceptable accuracy of the present formulations.

Comparison Study
A series of comparisons are made here to show whether the developed procedure is accurate enough.As a first verification study, acoustic effects and honeycomb core layer were assumed to be nonexistent, then the natural frequencies (Hz) of the FGM cylindrical shell were calculated for different values of material gradient index based on the present formulation.The results alongside other numerical findings in Ref. [34] are provided in Table 2, demonstrating the acceptable accuracy of the present formulations.As a further verification targeting the structural part of the developed procedure, acoustic and FGP effects were neglected, and the natural frequencies (Hz) of the honeycomb sandwich cylindrical shell structure were extracted for different values of core-to-thickness ratios.As shown in Figure 3, the developed formulation almost matches the previous numerical findings in another study by Li et al. [10] (Flügge shell theory).As a further verification targeting the structural part of the developed procedure, acoustic and FGP effects were neglected, and the natural frequencies (Hz) of the honeycomb sandwich cylindrical shell structure were extracted for different values of core-tothickness ratios.As shown in Figure 3, the developed formulation almost matches the previous numerical findings in another study by Li et al. [10] (Flügge shell theory).Also, as shown in Figure 4a,b, for an elastic cylindrical shell, STL from an incident acoustic wave at  = /4 was found from the derived formulation, then compared with the results of 3D elasticity theory in Ref. [55] and those of an FSDT approach in Ref. [71] in the frequency domain for the properties mentioned in those studies.As another validation, for a simplified case of aluminum shell, a comparison with the present results based on the FSDT the classical shell theory was performed.STL was found for an incident wave angle at  = /4.As observed in Figure 4c, there is practically no difference between the findings of this study and those of a numerical method by Ref. [63] (classical shell theory).Also, as shown in Figure 4a,b, for an elastic cylindrical shell, STL from an incident acoustic wave at α = π/4 was found from the derived formulation, then compared with the results of 3D elasticity theory in Ref. [55] and those of an FSDT approach in Ref. [71] in the frequency domain for the properties mentioned in those studies.As another validation, for a simplified case of aluminum shell, a comparison with the present results based on the FSDT the classical shell theory was performed.STL was found for an incident wave angle at α = π/4.As observed in Figure 4c, there is practically no difference between the findings of this study and those of a numerical method by Ref. [63] (classical shell theory).

Parametric Study
Figure 5 shows the effect of incidence angles on the STL of a honeycomb cylindrical shell with FGP layers over a wide frequency range (1 < f < 10 Hz) w 1.5m , ℎ = 1mm,  = 5mm,  = /6,  = 5mm,  = 0.5mm,  = 1,  = 1Pa,  = STL plots are generally characterized by a set of frequencies that allow us to dif different behaviors observed in related acoustic problems.These frequencies, order of occurrence, are named ring, critical, and coincidence frequency (respec noted by  ,  ,  ).The reason behind the occurrence of such frequencies is in detail in other studies [71,72].However, it should be kept in mind that each formed by them gives rise to distinct phenomena.The frequencies lower than  stiffness-controlled region where the structural stiffness plays the most significa STL.In contrast, the mass-controlled region that starts from  and ends at  affected by the structure's mass.Finally, the coincidence-controlled region lies a The value of  primarily depends on the incidence angle of sound waves.ST and incidence angle are inversely proportional to each other.It is noteworthy tha valleys before the ring frequency are the shell resonances.

Parametric Study
Figure 5 shows the effect of incidence angles on the STL of a honeycomb sandwich cylindrical shell with FGP layers over a wide frequency range (1 < f < 10 4 Hz) when R = 1.5m , h = 1mm, H = 5mm, θ c = π/6, l c = 5mm, t c = 0.5mm, p = 1, p 0 = 1 Pa, M = 0, ξ = 0. STL plots are generally characterized by a set of frequencies that allow us to differentiate different behaviors observed in related acoustic problems.These frequencies, from the order of occurrence, are named ring, critical, and coincidence frequency (respectively denoted by f r , f cr , f co ).The reason behind the occurrence of such frequencies is described in detail in other studies [71,72].However, it should be kept in mind that each interval formed by them gives rise to distinct phenomena.The frequencies lower than f r form the stiffness-controlled region where the structural stiffness plays the most significant role in STL.In contrast, the mass-controlled region that starts from f r and ends at f cr is mostly affected by the structure's mass.Finally, the coincidence-controlled region lies above f co .The value of f co primarily depends on the incidence angle of sound waves.STL values and incidence angle are inversely proportional to each other.It is noteworthy that the STL valleys before the ring frequency are the shell resonances.As shown in Figure 6, the influence of material gradient index on the STL plots was investigated when  = 1.5 m , ℎ = 1 mm,  = 5 mm,  = ,  = 5mm,  = 0.5 mm,  = ,  = 1Pa,  = 0,  = 0. Figure 6 offers a good understanding of how to achieve higher levels of STL in different frequency regions.For example, reducing the FG index, i.e., attaining a stronger ceramic behavior instead of metallic behavior, is one way to increase STL in the stiffness-controlled region.The exact opposite is true for the mass-controlled region where reducing the FG index decreases STL.What is more, the characteristic frequencies can be shifted by changing the parameters of power law.As shown in Figure 6, the influence of material gradient index on the STL plots was investigated when R = 1.5 m , h = 1 mm, H = 5 mm, θ c = π 6 , l c = 5 mm, t c = 0.5 mm, α = π 6 , p 0 = 1Pa, M = 0, ξ = 0. Figure 6 offers a good understanding of how to achieve higher levels of STL in different frequency regions.For example, reducing the FG index, i.e., attaining a stronger ceramic behavior instead of metallic behavior, is one way to increase STL in the stiffness-controlled region.The exact opposite is true for the mass-controlled region where reducing the FG index decreases STL.What is more, the characteristic frequencies can be shifted by changing the parameters of power law.As shown in Figure 6, the influence of material gradient index on the STL plots was investigated when  = 1.5 m , ℎ = 1 mm,  = 5 mm,  = ,  = 5mm,  = 0.5 mm,  = ,  = 1Pa,  = 0,  = 0. Figure 6 offers a good understanding of how to achieve higher levels of STL in different frequency regions.For example, reducing the FG index, i.e., attaining a stronger ceramic behavior instead of metallic behavior, is one way to increase STL in the stiffness-controlled region.The exact opposite is true for the mass-controlled region where reducing the FG index decreases STL.What is more, the characteristic frequencies can be shifted by changing the parameters of power law.It was expected to attain a lower STL for a larger radius before the ring frequency.In a stark contrast, changing the radius has virtually no effect on STL in the coincidence control region.This is due to the fact that at the high frequencies the wavelengths become shorter and therefore radius would incorporate no effect on STL in the coincidence control region.It should be noted that only the ring frequency is shifted as a result of a new radius.The variations of STL under different porosity distributions and porosity coefficient are indicated in Figure 9 when  = 0 , ℎ = 1 mm,  = 5 mm,  = ,  = 5 mm,  = 0.5 mm,  = ,  = 1 Pa,  = 1,  = 1.5.Turning our attention to the porosity coefficient, one observes a degraded STL performance with growing porosity coefficient irrespective of porosity pattern.To justify this behavior, one must pay attention to the decrease in the structural stiffness with growing porosity coefficient.This also explains the higher levels of STL for the case of logarithmic-uneven distribution as it possesses a higher stiffness   It was expected to attain a lower STL for a larger radius before the ring frequency.In a stark contrast, changing the radius has virtually no effect on STL in the coincidence control region.This is due to the fact that at the high frequencies the wavelengths become shorter and therefore radius would incorporate no effect on STL in the coincidence control region.It should be noted that only the ring frequency is shifted as a result of a new radius.The variations of STL under different porosity distributions and porosity coefficient are indicated in Figure 9 when  = 0 , ℎ = 1 mm,  = 5 mm,  = ,  = 5 mm,  = 0.5 mm,  = ,  = 1 Pa,  = 1,  = 1.5.Turning our attention to the porosity coefficient, one observes a degraded STL performance with growing porosity coefficient irrespective of porosity pattern.To justify this behavior, one must pay attention to the decrease in the structural stiffness with growing porosity coefficient.This also explains the higher levels of STL for the case of logarithmic-uneven distribution as it possesses a higher stiffness The variations of STL under different porosity distributions and porosity coefficient are indicated in Figure 9 when M = 0 , h = 1 mm, H = 5 mm, θ c = π 6 , l c = 5 mm, t c = 0.5 mm, α = π 6 , p 0 = 1 Pa, p = 1, R = 1.5.Turning our attention to the porosity coefficient, one observes a degraded STL performance with growing porosity coefficient irrespective of porosity pattern.To justify this behavior, one must pay attention to the decrease in the structural stiffness with growing porosity coefficient.This also explains the higher levels of STL for the case of logarithmic-uneven distribution as it possesses a higher stiffness than even porosity distribution.Nonetheless, for both porosity patterns, a larger porosity coefficient shifts the critical and coincidence frequencies to the right of the frequency spectrum.
than even porosity distribution.Nonetheless, for both porosity patterns, a larger porosity coefficient shifts the critical and coincidence frequencies to the right of the frequency spectrum.
The increase in the wall thickness of each honeycomb cell increases the core stiffness, so that STL grows.Furthermore, improvements in STL in the mass-controlled regions is observable.
Figure 11 demonstrates the variations of STL curves under different values of hexagonal cell length when M = 0 , h = 1 mm, H = 5 mm, θ c = π 6 , t c = 0.5 mm, ξ = 0, α = π 6 , p 0 = 1 Pa, p = 1, R = 1.5.Decreasing the value of l c concurrently increases the elastic modulus, total bending stiffness, and structure's density.In addition, the total bending stiffness of the whole structure is decreased as a result of regulation impact, which increases STL in the stiffness-controlled and mass-controlled regions.

Conclusions
Using the FSDT as a capable means of investigating vibration problems, this article develops a formulation to determine the acousto-structural behavior of a sandwich structure incorporating a hexagonal honeycomb core and face layers made from FGP.A homogenous orthotropic model was used for the honeycomb core while its corresponding material features were found from the modified Gibson's equation.The distribution pattern of FGP parts is either even or logarithmic-uneven, and a special rule-of-mixture rela-

Conclusions
Using the FSDT as a capable means of investigating vibration problems, this article develops a formulation to determine the acousto-structural behavior of a sandwich structure incorporating a hexagonal honeycomb core and face layers made from FGP.A homogenous orthotropic model was used for the honeycomb core while its corresponding material features were found from the modified Gibson's equation.The distribution pattern of FGP parts is either even or logarithmic-uneven, and a special rule-of-mixture rela-

Conclusions
Using the FSDT as a capable means of investigating vibration problems, this article develops a formulation to determine the acousto-structural behavior of a sandwich structure incorporating a hexagonal honeycomb core and face layers made from FGP.A homogenous orthotropic model was used for the honeycomb core while its corresponding material features were found from the modified Gibson's equation.The distribution pattern of FGP parts is either even or logarithmic-uneven, and a special rule-of-mixture relation govern their properties.Then, fluid-structure compatibility relations were implemented followed by the use of Hamilton's principle, leading to the final form of governing equations.A series of simulations and comparisons with previous studies demonstrate the accuracy of the presented method.STL was then calculated over different frequency regions.The results are presented as follows:

•
The increase in the wall thickness of each honeycomb cell increases the core stiffness, so that STL grows.

•
Decreasing the value of l c concurrently increases the elastic modulus, total bending stiffness, and structure's density, which increases STL in the stiffness-controlled and mass-controlled regions.

•
STL reduces with growing porosity coefficient irrespective of porosity pattern.

•
Higher levels of STL were obtained for the case of logarithmic-uneven distribution as it possesses a higher stiffness than even porosity distribution.

•
Reducing the FG index, i.e., attaining a stronger ceramic behavior instead of metallic behavior, is one way to increase STL in the stiffness-controlled region.

•
The radiation damping after f r has the capability of increasing STL with a growing Mach number.

Figure 1 .
Figure 1.The schematic of a honeycomb sandwich cylindrical shell with FGP layers under incidence wave.

Figure 1 .
Figure 1.The schematic of a honeycomb sandwich cylindrical shell with FGP layers under incidence wave.

Figure 3 .
Figure 3.Comparison study of the natural frequencies (Hz) of the honeycomb sandwich cylindrical shell, red solid line: present, black dotted line: Ref [10].

Figure 3 .
Figure 3.Comparison study of the natural frequencies (Hz) of the honeycomb sandwich cylindrical shell, red solid line: present, black dotted line: Ref [10].

Figure 6 .
Figure 6.Variation of STL through the honeycomb FGP sandwich cylindrical shell with material gradient index.To analyze the effect of the external flow Mach on the variations of STL curves, Figure 7 is presented, where  = 1.5 m , ℎ = 1 mm,  = 5 mm,  = ,  = 5mm,  = 0.5mm,  = ,  = 1Pa,  = 1,  = 0. Another noteworthy finding is that the radiation damping after the ring frequency has the capability of increasing STL with a growing Mach number.In other words, a new value of Mach number results in different values of  and  .

Figure 6 .
Figure 6.Variation of STL through the honeycomb FGP sandwich cylindrical shell with material gradient index.To analyze the effect of the external flow Mach on the variations of STL curves, Figure 7 is presented, where  = 1.5 m , ℎ = 1 mm,  = 5 mm,  = ,  = 5mm,  = 0.5mm,  = ,  = 1Pa,  = 1,  = 0. Another noteworthy finding is that the radiation damping after the ring frequency has the capability of increasing STL with a growing Mach number.In other words, a new value of Mach number results in different values of  and  .

Figure 6 .
Figure 6.Variation of STL through the honeycomb FGP sandwich cylindrical shell with material gradient index.To analyze the effect of the external flow Mach on the variations of STL curves, Figure 7 is presented, where R = 1.5 m , h = 1 mm, H = 5 mm, θ c = π 6 , l c = 5 mm, t c = 0.5 mm, α = π 6 , p 0 = 1 Pa, p = 1, ξ = 0. Another noteworthy finding is that the radiation damping after the ring frequency has the capability of increasing STL with a growing Mach number.In other words, a new value of Mach number results in different values of f cr and f co .

Figure 7 .
Figure 7. Effect of the external flow Mach number on the changes of STL.

Figure 8 .
Figure 8. STL plots for different values of the average radius.

Figure 7 .
Figure 7. Effect of the external flow Mach number on the changes of STL.

Figure 8 21 Figure 7 .
Figure8depicts the effect of the average radius on the variations of STL curves when M = 0 , h = 1 mm, H = 5 mm, θ c = π 6 , l c = 5 mm, t c = 0.5 mm, α = π 6 , p 0 = 1 Pa, p = 1, ξ = 0.It was expected to attain a lower STL for a larger radius before the ring frequency.In a stark contrast, changing the radius has virtually no effect on STL in the coincidence control region.This is due to the fact that at the high frequencies the wavelengths become shorter and therefore radius would incorporate no effect on STL in the coincidence control region.It should be noted that only the ring frequency is shifted as a result of a new radius.

Figure 8 .
Figure 8. STL plots for different values of the average radius.

Figure 8 .
Figure 8. STL plots for different values of the average radius.

Figure 9 .
Figure 9.Effect of different porosity distributions and porosity coefficient on the changes of STL.

Figure 10
Figure 10 displays the effect of the wall thickness of honeycomb cell on the variations of STL over frequency range when  = 0 , ℎ = 1 mm,  = 5 mm,  = ,  = 5 mm,  = 0,  = ,  = 1 Pa,  = 1,  = 1.5.The increase in the wall thickness of each honeycomb cell increases the core stiffness, so that STL grows.Furthermore, improvements in STL in the mass-controlled regions is observable.

Figure 9 .
Figure 9.Effect of different porosity distributions and porosity coefficient on the changes of STL.

Figure 10
Figure 10 displays the effect of the wall thickness of honeycomb cell on the variations of STL over frequency range when M= 0 , h = 1 mm, H = 5 mm, θ c = π 6 , l c = 5 mm, ξ = 0, α = π 6 , p 0 = 1 Pa, p = 1, R = 1.5.The increase in the wall thickness of each honeycomb cell increases the core stiffness, so that STL grows.Furthermore, improvements in STL in the mass-controlled regions is observable.Figure11demonstrates the variations of STL curves under different values of hexagonal cell length when M = 0 , h = 1 mm, H = 5 mm, θ c = π 6 , t c = 0.5 mm, ξ = 0, α = π 6 , p 0 = 1 Pa, p = 1, R = 1.5.Decreasing the value of l c concurrently increases the elastic modulus, total bending stiffness, and structure's density.In addition, the total bending stiffness of the whole structure is decreased as a result of regulation impact, which increases STL in the stiffness-controlled and mass-controlled regions.

Figure 10 .
Figure 10.Effect of the wall thickness of honeycomb cell on the changes of STL.

Figure 11
Figure 11 demonstrates the variations of STL curves under different values of hexagonal cell length when  = 0 , ℎ = 1 mm,  = 5 mm,  = ,  = 0.5 mm,  = 0,  = ,  = 1 Pa,  = 1,  = 1.5.Decreasing the value of  concurrently increases the elastic modulus, total bending stiffness, and structure's density.In addition, the total bending stiffness of the whole structure is decreased as a result of regulation impact, which increases STL in the stiffness-controlled and mass-controlled regions.

Figure 11 .
Figure 11.Effect of the hexagonal cell length on the changes of STL.

Figure 10 .
Figure 10.Effect of the wall thickness of honeycomb cell on the changes of STL.

Figure 10 .
Figure 10.Effect of the wall thickness of honeycomb cell on the changes of STL.

Figure 11
Figure 11 demonstrates the variations of STL curves under different values of hexagonal cell length when  = 0 , ℎ = 1 mm,  = 5 mm,  = ,  = 0.5 mm,  = 0,  = ,  = 1 Pa,  = 1,  = 1.5.Decreasing the value of  concurrently increases the elastic modu-lus, total bending stiffness, and structure's density.In addition, the total bending stiffness of the whole structure is decreased as a result of regulation impact, which increases STL in the stiffness-controlled and mass-controlled regions.

Figure 11 .
Figure 11.Effect of the hexagonal cell length on the changes of STL.

Figure 11 .
Figure 11.Effect of the hexagonal cell length on the changes of STL.

Table 2 .
A comparative study of the natural frequencies of an FG cylindrical shell.

Table 2 .
A comparative study of the natural frequencies of an FG cylindrical shell.