Abstract
Vibrational characteristics of functionally graded cylindrical shells filled with fluid and placed on Winkler and Pasternak elastic foundations are investigated. Love's thin-shell theory is utilized for strain-displacement and curvature-displacement relationships. Shell dynamical equations are solved by using wave propagation approach. Natural frequencies for both empty and fluid-filled functionally graded cylindrical shells based on elastic foundations are determined for simply-supported boundary condition and compared to validate the present technique. Results obtained are in good agreement with the previous studies. It is seen that the frequencies of the cylindrical shells are affected much when the shells are filled with fluid, placed on elastic foundations, and structured with functionally graded materials. The influence of Pasternak foundation is more pronounced than that of Winkler modulus.
1. Introduction
Cylindrical shells have vast applications in many engineering fields such as aerospace, civil, chemical, mechanical, naval, and nuclear. First thin-shell theory was developed by Love [1]. He based his theory on Kirchhoff's approximations for plate and beam problems. Arnold and Warburton [2] did some pioneering work on the vibration characteristics of cylindrical shells. Best reviews on the shell vibration characteristics have been given by Leissa [3], Blevins [4], and Markus [5]. Different types of materials like isotropic, laminated and composite like functionally graded materials are used in fabricating shells.
Loy et al. [6] did vibration frequency analysis of functionally graded cylindrical shells. They have employed the Raleigh-Ritz numerical technique to evaluate natural frequencies for simply-supported cylindrical shells. The axial modal dependence was approximated by trigonometric functions. They studied their characteristics with regard to circumferential wave numbers and geometrical parameters of length-to- radius and thickness-to- radius ratios. Pradhan et al. [7] did work on vibration characteristics of functionally graded cylindrical shells made of stainless steel and zirconia for a number of boundary conditions. They also used the Rayleigh-Ritz method to solve the shell governing equations. Characteristic beam functions were utilized for the axial modal displacement deformations. In 2002, Naeem [8] worked on vibration analysis of nonrotating and rotating FGM circular cylindrical shells using Rayleigh-Ritz method and Galerkin technique, respectively. Ritz polynomial functions for nonrotating shells and characteristic beam functions for rotating ones approximate axial modal dependences. Najafizadeh and Isvandzibaei [9] have studied the vibration characteristics of functionally graded material cylindrical shells with ring support. They based their analysis on higher-order plate theory. Arshad et al. [10] studied the natural frequencies of FGM cylindrical shells by exponential volume fraction law with base . Shah et al. [11] further modified and extended this volume fraction law to a general base (). Arshad et al. [12] also extended their work of reference [10] to study the natural frequencies of FGM cylindrical shells under the influence of various boundary conditions with exponential volume fraction law. Arshad et al. [13, 14] also evaluated natural frequencies of bi-layered cylindrical shells by inserting different materials in the layers of the shells. Rayleigh-Ritz numerical technique has been employed to solve the governing equations of motions in the references [9β14].
Fluid-filled cylindrical shells have paramount importance in engineering and industries to design pressure vessels, fluid tanks and underground pipe fitting. Junger [15, 16] maidenly studied the structural response of thin-shell in an acoustic fluid media by studying the free and forced vibrations of cylinders submerged in an acoustic medium. Chung et al. [17] gave an analytical and experimental study of the vibrational characteristics of thin circular cylindrical shells filled with fluid. Goncalves and Batista [18] gave a theoretical study for the free vibration of vertical shells with simply-supported end condition, partially filled with or submerged in a fluid. They used Sandersβ shell theory for acquiring shell frequency equation coupling with shell dynamical equation and the acoustic medium of fluid equations. Amabili [19] has studied the natural frequencies and mode shapes of a simply-supported circular cylindrical shell placed horizontally partially filled with liquid. He [20] further extended his study to investigate the free vibrations of circular cylindrical shells and tubes entirely filled with a dense and partially immersed in various fluids under different end conditions. Sharma et al. [21] have determined the natural frequencies of vertical cantilever fluid-filled composite shells by utilizing the Fourier series of trigonometric functions to approximate the axial modal dependence. Zhang et al. [22] solved the shell frequency equations with wave propagation approach to study the vibrational response of cylindrical shells containing fluid. Natsuki et al. [23] used this technique for the vibrational analysis of carbon nanotubes having fluid inside. Iqbal et al. [24] investigated the vibration characteristics of FGM circular cylindrical shells using wave propagation approach. Shah et al. [25] studied the vibration characteristics of functionally graded cylindrical shells resting on Winkler and Pasternak elastic foundations by using wave propagation approach. This approach seem to be a simple and easily applicable to study vibration characteristics of shells involving fluids. In the wave propagation method, the axial modal dependence is approximated by simple wave number characteristic formulas. Paliwal et al. [26] and Paliwal and Pandey [27] have studied the free vibrations of circular cylindrical shells resting on elastic foundations.
In the present study, the vibration characteristics of fluid-filled cylindrical shells fabricated from functionally graded materials resting on the elastic foundations are investigated. The acoustic wave equation of fluid is coupled with the shell dynamical equations. This gives the mathematical relationship between solid and fluid dynamics. The influence of FGM configurations on shell vibrations is analyzed for various physical and geometrical parameters. For strain-displacement and curvature displacement relationships, Love's thin-shell theory is used and wave propagation approach is employed to solve the shell dynamical equations involving the Winkler and Pasternak elastic foundation and fluid loading terms. A good agreement is found between the present and the literature's results on comparison.
2. Theoretical Formulation
A circular cylindrical shell is considered for its vibration characteristics with length , thickness , and mean radius as shown in Figure 1. The material properties are Young's modulus , the Poisson ratio , and the mass density . An orthogonal coordinate system () is assumed to be established at the shell mid surface where , , and represent the axial, circumferential, and radial coordinates, respectively. , and are axial, circumferential and radial deformation displacements, respectively, of the shell mid surface.
The equations of motion for a cylindrical shell derived from the Loveβs thin-shell theory can be written as where , , and are force resultants and , , moment resultants along the normal and shear directions, respectively. They are defined as where , , and are the reference surface strains and , and are the surface curvatures. , , and ( and 6) are the membrane, coupling, and flexural stiffnesses defined as
The reduced stiffnesses ( and 6) for isotropic materials are defined as
The cylindrical shell considered in this study is assumed to be of thin with thickness to radius ratio less than 1/20. For strain displacements and the curvature displacements Love's thin-shell theory is utilized and these quantities are defined as
The above-mentioned expressions for the surface strains and curvatures are substituted into (2) to get the forces and moments of the shell and then putting them into shell dynamical (1)β(3) and introducing the fluid acoustic pressure inside the cylinder and the Winkler and Pasternak foundations terms () in the -direction outside the cylinder. The equations of motion for a modified cylindrical shell can be written in the differential operator form as: where (, 2, 3) are the differential operators with respect to and given in Appendix A. represents the shear modulus of the material used for the elastic foundation and for the Winkler foundation modulus, and the expression for the differential operator is
The Winkler model is a special case of the Pasternak model when . The equation of motion of the fluid can be written in the cylindrical coordinate system (, , ) as: where is the time, is the acoustic pressure, and is the speed of sound in the fluid.
Wave propagation approach has been employed by many researchers (Zhang et al. [22, 29, 30], Natsuki et al. [23], Li et al. [31], Liu et al. [32], and Li [28]) to evaluate the approximate solution of the shell dynamic equations. This numerical technique is used in this study as it is very simple and easily applicable to determine the shell frequencies. For separating the spatial and temporal variables, the following shapes of modal displacement functions are assumed as
The coefficients denote the wave amplitudes in the , , and directions, respectively. denotes the circumferential waves number, is the natural angular frequency for the cylindrical shell, and is the axial wave number that has been specified in the reference [18] for a number of boundary conditions. The acoustic pressure field in the cylindrical shell containing fluid, which satisfies the acoustic wave (8), can be expressed as where is the Bessel function of first kind with order , and is radial wave number. Coupling condition of the fluid with the shell wall is given by at and where stands for the density of the contained fluid and the prime on the represents the differentiation of the Bessel function with respect to the argument . By substituting the expressions for , , and from (9) into (6) and simplifying the algebraic expressions and rearranging the terms, the shell frequency equation leads to the following eigenvalue form: where are some matrix coefficients depending on shell parameters and the type of boundary conditions specified at the ends of an FG cylindrical shell and are given in Appendix B. Equation (13) is solved for shell frequencies and mode shapes by using MATLAB software. The three frequencies are obtained corresponding to the axial, circumferential, and radial displacements.
3. Functionally Graded Materials
Functionally graded materials are advanced materials and are used in engineering and technology applications due to their mechanical properties. Their best use is found in the thermal environment for their superb properties. They are fabricated from two or more materials. The material properties of their constituents are functions of the temperature and the volume fractions. If represents a material property of the ith constituent material of an FGM consisting of constituent materials, then its effective material is written as where is the volume fraction of the ith constituent material. Also the sum of volume fractions of the constituent materials is equal to one, that is, The volume fraction depends on the thickness variable and is defined as for a cylindrical shell. and denote inner and outer radii of the shell and is the thickness variable in the radial direction. is known as power law exponent and is a nonnegative real number and lies between zero and infinity. When a cylindrical shell is considered to be consisting of two constituent functionally graded materials and , the volume fraction of the outer shell surface is obtained from (16) as where is the shell uniform thickness. The effective Young's modulus , Poisson's ratio , and the mass density are given by where , are Young's moduli, , Poisson's ratios, and , the mass densities of the constituent materials and , respectively. From the relations (18), the following things are noted, that is, at the inner surface the FGM properties are those of the constituent material and at the outer surface, they are those of material . Thus the FGM properties change continuously from the material at the inner surface to the material at the outer surface.
4. Results and Discussions
In Table 1, a comparison of frequency parameter is made against circumferential wave number with Loy et al. [6] for a simply-supported isotropic cylindrical shell with shell parameters , , , and . Table 2 shows the comparison of the variation of frequency parameter of isotropic cylindrical shells with circumferential wave numberand ratio to that of Loy et al. [33] for simply-supported end condition having shell parameters (, ). In Table 3, comparison of natural frequencies (Hz) for simply-supported type I and type II FGM cylindrical shell is made with Loy et al. [6] against circumferential wave number for various values of power law exponent , 2, and 15 with shell geometrical and material parameters (, , , and ). It is observed that results agree very well with one another. This shows accuracy and validity of the present approach.
In Tables 4, 5, and 6, variations of natural frequencies (Hz) against circumferential wave number are given for types I & II FG cylindrical shells on elastic foundation. Shells parameters are given in the tables. In Table 4, Pasternak foundation is zero and Winkler foundation is taken equal to . Values of exponential volume fraction law are , 0.7, 1, 2, 5, 15, and 30 in all cases. It is observed that frequencies of both shells increase vertically with increasing values of circumferential wave number , whereas it slowly decreases for shell type I but increases for shell type II with increasing value of exponential volume fraction law. Frequency of the pure stainless steel is more than frequency of pure nickel. Frequency of FG cylindrical shell remains between the frequency of stainless steel and nickel when FG cylindrical shell is on elastic foundation in Tables 5 and 6; all parameters of shell are same as given in Table 4 but value of and . In this case, frequency increases rapidly with increasing values of circumferential wave number as compared with frequency when , . It means that influence of is more than that of on the shell frequencies. Here again frequencies decrease and increase with increasing values of power law exponent, for shell types I & II, respectively, but remain between the frequencies of pure stainless steel and pure nickel. In Table 6, variation of natural frequencies (Hz) against circumferential wave number is given when both and are nonzero. Natural frequencies increase rapidly with increasing value of circumferential wave number but it decreases and increases with increasing values of and it always remains between the frequencies of pure stainless steel and nickel.
In Tables 7, 8, and 9, variations of natural frequencies (Hz) against circumferential wave number of type I & II functionally graded fluid-filled cylindrical shells on elastic foundation are given. Parameters are given in the tables. In Table 7, and . Value of and in Table 8, and , in Table 9. In Table 7, it is observed that frequencies increase with circumferential wave number . It decreases and increases very slowly with increasing value of . It is observed that their frequencies remain within the frequencies of pure stainless steel and pure nickel fluid-filled cylindrical shells on elastic foundation, and it is also observed that frequency of fluid-filled FGM cylindrical shell is much less than that of empty FG cylindrical shell on the elastic foundation. Tables 8 and 9 also show that the influence of is more than that of when fluid-filled FG cylindrical shell is on elastic foundation. Frequency increases with increasing values of and . These tables show that influence of Pasternak foundation is more than that of Winkler foundation. Also it is clear that frequency of fluid-filled FG cylindrical shell is much less than that of empty FG cylindrical shell on elastic foundation.
Figure 2 shows a similar and opposite behaviour of the natural frequencies (Hz) of both the configurations of fluid-filled FGM type I & type II cylindrical shells for the same values of Winkler and Pasternak elastic foundations and shell parameters as shown in this figure. It is clear from this figure that natural frequencies of shell type I decrease while that of shell type II increase with the increasing values of power law exponent . Rate of decrement of shell type I and increment of shell type II is the same for the ascending values of and frequencies of both types of the shells seem to intersect at .
In Figures 3(a)β3(d), variation of natural frequencies (Hz) of fluid-filled isotropic as well as FGM type I cylindrical shells resting on elastic foundations is sketched against ratio with shell parameters , , and , by keeping Winkler elastic modulus () constant and Pasternak elastic modulus varied from to with class interval, respectively. Nickel and stainless steel are used in isotropic cylindrical shells, whereas these isotropic materials are also used in the fabrication of FGM cylindrical shells as constituents. Influence of is studied in both isotropic and three FGM cylindrical shells corresponding to power law exponents , 1, 5, respectively.
(a)
(b)
(c)
(d)
It is observed from these figures that natural frequencies (Hz) of these cylindrical shells first decrease rapidly from to , for these frequencies decline very slowly by increasing values of ratio. Natural frequencies of FG cylindrical shells always lie between the natural frequencies of isotropic cylindrical shells to which FG cylindrical shells are fabricated. By increasing the values of , natural frequencies of the all types of cylindrical shells are also considerable increased. It is further noted that for lower values of ratio say , natural frequencies of the cylindrical shells seem to be converged but for higher values of ratio these frequencies show diverging behaviour. Natural frequencies of isotropic cylindrical shell made with nickel are low as compared to that of stainless steel isotropic cylindrical shells for small value of say , but for higher values of say this order of natural frequencies is changed. It is further observed that for small values of say , natural frequencies of type I FGM cylindrical shells with lower values of power law exponent are higher than those with higher values of , but for higher values of say , the order of frequencies of the FGM shells is reversed.
In Figures 4(a)β4(d), variation of natural frequencies (Hz) of fluid-filled isotropic as well as FGM type I cylindrical shells resting on elastic foundations is drawn with ratio having shell parameters , , and , by keeping Pasternak elastic modulus () while Winkler elastic modulus is varied from to with class interval, respectively. Influence of is studied on these cylindrical shells.
(a)
(b)
(c)
(d)
These figures have all the properties of Figures 3(a)β3(d) except that Winkler modulus does not affect the order of natural frequencies of isotropic as well as FGM cylindrical shells.
In Figures 5(a)β5(d), variation of natural frequencies (Hz) of fluid-filled isotropic as well as FGM type I cylindrical shells placed on elastic foundations is sketched against ratio with shell parameters , , and , by taking Winkler elastic modulus () constant and Pasternak elastic modulus varied from to with class interval . Influence of is studied on the natural frequencies of these cylindrical shells.
(a)
(b)
(c)
(d)
It is observed from these figures that natural frequencies (Hz) of these cylindrical shells decrease simultaneously from to , afterward it decreases sharply from to 0.02, 0.025, 0.03, and 0.035 in figures 5(a)β5(d) respectively and afterwards these frequencies increase rapidly with ratio. It is further noted that natural frequencies of FG cylindrical shells always lie between the natural frequencies of isotropic cylindrical shells to which FG cylindrical shells are fabricated. For lower values of , natural frequencies of cylindrical shells with small ratio are low as compared to high ratio but for higher values of , this order of frequencies is reversed with ratio. For lower values of ratio, natural frequencies of the cylindrical shells seem to be converged, but for higher values of ratio, these frequencies show diverging behaviour. Natural frequencies (Hz) of isotropic cylindrical shell made with nickel are low as compared to that of stainless steel isotropic cylindrical shells for small value of say ; for higher values of say , this order of natural frequencies remains the same. Similarly, natural frequencies of FG cylindrical shells with low value of power law exponent say are higher in type I and lower in type II FGM cylindrical shells, whereas for higher value of power law exponent say , the natural frequencies are lower for type I and higher in type II FGM cylindrical shells. It is clear from these figures that variation of does not affect the order of natural frequencies of these cylindrical shells. It is further noted that for increasing values of , natural frequencies of the cylindrical shells also increased; this rate of increment is higher for low values of ratio than those of higher values of ratio. It is further observed that for smaller values of say , natural frequencies of type I & II FGM cylindrical shells with lower values of ratio are smaller than those with higher values of , but for higher values of say , the order of frequencies of the FGM shells is reversed. It means that Pasternak modulus has prominent effect on the natural frequencies of the cylindrical shells.
In Figures 6(a)β6(d), variation of natural frequencies (Hz) of fluid-filled isotropic as well as FGM type I cylindrical shells placed on elastic foundations is drawn with increasing values of ratio having shell parameters , , , considering Pasternak elastic modulus () constant and Winkler elastic modulus varied from to with class interval . Influence of is studied on the natural frequencies of these cylindrical shells.
(a)
(b)
(c)
(d)
The same behaviour is observed as in Figures 5(a)β5(d) except that Winkler modulus also has no effect on the order of natural frequencies of isotropic as well as FGM cylindrical shells. Natural frequencies of cylindrical shells increased with but very slowly. Rate of increment for lower values of ratio is the same for higher values of ratio with increasing values of . Rate of increment of natural frequencies of the cylindrical shells with Pasternak modulus is very high as compared to Winkler modulus. It meansthat, Pasternak modulus directly influences the natural frequencies of the cylindrical shells under considered. Similar and opposite behaviour of shell frequencies can be observed for FGM shell type II on elastic foundations.
5. Conclusion
In this study, vibration frequency analysis of empty as well as fluid-filled FGM cylindrical shells based on elastic foundations is presented for simply-supported boundary conditions. Wave propagation approach is used to derive the shell frequency equation in the form of the eigenvalue problem including the elastic foundation and fluid loading terms. The influence of elastic moduli is more pronounced on the shell frequencies. The frequency of the fluid-filled FG cylindrical shell is much less than that of empty FG cylindrical shell based on elastic foundations.