Frequency Analysis for Functionally Graded Material Cylindrical Shells: A Significant Case Study

Department of Mathematics and Statistics, University of Agriculture, Faisalabad 38000, Pakistan Department of Mathematics, College of Arts and Sciences, University of Jhang, Jhang 35200, Pakistan Department of Mathematics, Lahore College Women University, Lahore 54000, Pakistan Department of Mathematics, Capital University of Science and Technology, Islamabad 44000, Pakistan Department of Mathematics, Art and Science Faculty, Siirt Üniversity, TR-56100 Siirt, Turkey Department of Mathematics, Dr. N.G.P. Institute of Technology, Coimbatore, India


Introduction
A shell is a structure consisting of curved sheets of material. A shell is the best effective way of consuming the material and can be very valuable in case of stowage of liquids and solids. Shells have different structures such as elliptic parboiled, hyperbolic parboiled, and circular cylindrical. Shell structures are generally used structural mechanisms in manufacturing designs. Vibration of cylindrical shells is a general field of research in mechanical dynamics. ese shells are made from various kinds of materials. ere are many types of shells, but CS is most important. e high strength, heat resistance, and corrosion resistance of FGM provide a reliable guarantee for the development of the nextgeneration nuclear industry. Nuclear power generation, nuclear power, and the storage of nuclear weapons are extremely dangerous. As a high-strength, heat-resistant, and shielding material, functionally graded materials show great superiority in the construction materials of nuclear furnaces and inner wall materials of nuclear furnaces, which greatly protect the safety of the nuclear industry.
Nowadays, FGMs are applicable in various fields including automotive, aerospace, mining, energy, and medical fields. Cylindrical shells are applicable in different appliances of mechanics, civil, nuclear, and chemical industries. ey have different layers, e.g., single layer, bilayers, and trilayers. ese layers are made up of "functionally graded" (FG) and "isotropic" materials [1]. Functionally graded materials are mostly utilized in the vast area of engineering and technology. Functionally graded material is a combination of more than one material. ese materials are stainless steel (SS), zirconia (Zr), and nickel (Ni). e scientists from Japan (1884) purposed the group of materials with the conception of FGMs. ese types of materials are manufactured by the combination of different materials employing the powder metallurgy method (powder metallurgy method is an industrial technique for attaining metals, such as powders). ey contain vibration in essential fraction volume that takes to continual transform in the porosity and composition microstructure. ese results are employed as a gradient in thermal and mechanical characteristics. Mostly, FGMs are applicable in the environment of maximum temperature. FGMs are used to construct the barrier of thermal materials. e analysis on FGMs tells that they are massive but mostly restrained to test the deformation and thermal stress [2]. e vibrational characteristics of the CS are composed of FGMs.
ere are multiple uses of FGMs in the industry as they are the most significant materials. ey are auspicious materials and have the ability to continue their presence in a maximum produced temperature condition. An enormous amount of research work has been completed on the functionally graded CS. e vibration analysis depends on four clamped conditions employed on both shell ends. ese conditions were expressed by Swaddiwudhipong et al. [3]. Rayleigh-Ritz method is an energy variation method that is used to find the natural frequency (NF). In 1909, Walther Ritz introduced this method. Rayleigh-Ritz method is also known as Galerkin's method. In other words, Rayleigh-Ritz technique is an integral approach method. It is valuable for resolving the problems of structural mechanics. It is also called the variation approach as given by Preumont [4]. Ng et al. [5] exhibited the construction for the dynamic solidness examination of the FG shell under the periodic axial loading. By the equation of motion, they produced the Mathieu-Hill equations and examined their stability by Bolotin's method. Rahimi et al. [6] examined the behavior of the vibrational characteristic of the CS with the supports of the ring placed at the centre of the shell. ey utilized energy functional to get motion's equation. ey considered the "functionally graded cylindrical shell" composed of the combination of ceramic and metal. Lee and Lu [7] examined the free vibration of cylindrical shells filled with liquid. ey utilized "Donnell's shell theory" to analyze the vibrations of cylindrical shells and employed "Hamilton's principle" to obtain the function of shell dynamics. To modify the equations of differential into algebraic, "Galerkin's method" was used. Tang et al. [8] established the analytical technique for the free vibrational analysis of the cylindrical shell with arbitrary edge condition by reverberation-ray matrix's method. ey utilized the Flugge thin theory to resolve the motion's equation and accomplished the exact solutions of travelling and standing waves from along the direction of axial and circumferential, respectively. ey studied the influence of elastic restraints on the parameters of frequency. Qin et al. [9] scrutinized the free vibrational behavior of CSs with random edge conditions. For the calculation of "elastic strain energy," they utilized Sander's theory. Artificial springs were carried out on the ends of the shells to symbolize the arbitrary edge conditions. A united solution for three specific sorts of enlargement function was developed by using the technique of Rayleigh-Ritz. Wang et al. [10] presented the Jacobi-Ritz formulation to explore the free vibrations of different coupled doubly curved revolt shell structures with subjective edge conditions. ey recognized the multisegment apportioning methodology and Flügge's shell theory to set up the theoretic model. By utilizing the Rayleigh-Ritz strategy, the natural frequencies and mode shapes of the coupled shell structures are chosen. Trabelsi et al. [11] inspected the behavior of thermal buckling of FG cylindrical shells and plates. e governing equations were expanded by considering the theory of first-order shear distortion. To solve the thermal buckling problem, fournode shell elements were assumed. e calculation of the censorious buckling temperature of structures under nonuniform temperature rise depends on Gaussian numerical combination.
e impact of material composition, edge condition, thermal loading, and geometrical parameters of shells on the thermal buckling conduct of FGM structures was also examined. Daud and Viswanathan [12] introduced the vibrational conduct of the symmetric angle-ply layered circular cylindrical shell filled by calm fluid. e motion equation of the cylindrical shell as far as stress and moment resultants was obtained from the first-order shear distortion theory. To attain the differential equations with the rotational function and displacements, the relation of "straindisplacement and stress-strain" was assumed. e investigation of frequency parameters was examined utilizing the parameters the thickness and length ratio, angle-ply, properties of material, and the number of layers with various edge conditions. Sheng and Wang [13] established a theoretic model in the thermal atmosphere to examine the nonlinear vibration and dynamic stability of the stiffened FG cylindrical shell. For the modelling of the stiffened FG cylindrical shell, they utilized different techniques and theories such as Bolotin's technique, smearing stiffener approach, and von Karman nonlinear and first-order shear deformation theory. To achieve the differential equation, modal analysis and Galerkin's method were employed. e influence of the thermal atmosphere and material features on the dynamic constancy and temporary and crucial resonance responses was also observed. Strozzi and Pellicano [14] analyzed the nonlinear vibration of the FGM cylindrical shell. ey applied Sanders-Koiter theory to determine the finite amplitude of vibration and deliberately supported the clamped and free boundary conditions. ey carried out numerical analysis to depict the nonlinear response by considering a diversity of axisymmetric and asymmetric modes. Rawat et al. [15] exhibited the finite element method (FEM) for the vibration analysis of the circular shell. ey employed various shell theories for FEM results to report in the literature. To investigate the computational efficacy, they used block "Lanczos and subspace iteration methods." ey investigated model frequencies for various types of closures. Budak et al. [16] used numerical and experimental approaches to determine the dynamic characteristics of the isotropic cylindrical shell. ey applied the finite element method with Femar software to accomplish numerical calculations. In order to observe the interferential pattern of the vibrational shell, they developed the stroboscopic holographic interferometry technique. Mat Daud and 2 Mathematical Problems in Engineering Viswanathan [12] presented the vibrational behavior of the symmetric angle-ply layered circular shell filled with quiescent fluid. ey derived the equation of motion from the first-order shear deformation theory and stated the wave equation. ey solved the eigenvalue problem for the frequency parameter of spline coefficient and examined the frequency parameter with different edge conditions. For further interests, the reader may read [17][18][19][20][21]. Figure 1 shows a cylinder-shaped shell in which R is the radius, h is the thickness, and L is the length. e orthogonal coordinate system (x, θ, z) is established at the central surface of the CS, where x, θ, and z lie in the axial, circumferential, and radial directions of the shell, and (u 1 , v 1 , w 1 ) are the displacements of the shell in x, θ, and z directions, respectively. e strain energy of a cylindrical shell is stated as

Materials and Methods
where where e 1 , e 2 , κ 1 , κ 2 , and 2τ represent the strains and curvatures' reference surface relations, respectively. A variety of shell theories utilized are [17][18][19][20]. Love's shell theory is the first shell model, and other models were obtained from the Love shell theory by adding few terms. According to Love's shell theory, the strain and curvature displacement relations can be defined as [S] is defined as For isotropic materials, Q ij is the reduced stiffness stated as where E is Young's modulus and ] is Poisson's ratio. For the isotropic cylindrical shell, B ij (coupling stiffness) is zero, and for FGM cylindrical shells, B ij becomes nonzero. With the help of expressions (2) and (5)  Mathematical Problems in Engineering 3 By putting (3) and (4) in (8), M attains the following form: e kinetic energy for a CS can be expressed as where ρ T is the mass density per unit length and is expressed as For a cylindrical shell, Lagrangian energy functional is stated as Axial modal dependence is employed to fulfill the edge conditions. For the approximation of axial modal dependence, different varieties of functions are adopted. e expression for axial, tangential, and radial displacements is supposed as where F, G, and H are vibration amplitudes, n is the circumferential wave number, and ω is the natural frequency (NF) of the shell. e axial function ϕ(x) is chosen to be a beam function as e original formula of the Lagrangian functional is established as where To obtain the shell frequency equation, the energy functional is maximized with regard to the unknowns F, G, and H of the CS. e following expression is attained as e shell frequency is expressed in the eigenform as where MATLAB software is used which is a dominant tool for the quick and converging solution of numerical problems. e MATLAB command, e.g., (K, M), is utilized to extract the eigenvalues and eigenvectors in the procedure of shell natural frequencies and mode shapes, respectively.

Functionally Graded Materials (FGMs).
Functionally graded materials are attained by mixing more than one material. Generally, FGMs are applied in the environment of maximum temperature. e material property P is stated as a function of temperature T, where T is expressed in kelvin [21].
where P 0 , P − 1 , P 1 , P 2 , and P 3 are the coefficients of T. e material property P of FGM and the volume fraction of the composing materials are stated as P � k j�1 P j V fj . (20) e summation of volume fraction of all composing materials must be unity.
e volume fraction for the single-layer CS having thickness h, made up of two composing materials M 1 (stainless steel) and M 2 (nickel), can be written as where z is the thickness variable and N is the power law exponent. E is Young's modulus, ] is Poisson's ratio, and ρ is the mass density, which can be stated for the FGM cylindrical shell consisting of two materials: Volume fraction for the bilayer cylindrical shell is represented by the following three cases.
Case I: When the inner layer is made of FGM and outer is of isotropic, Case II: When the inner layer is of isotropic and outer is of FGM, Case III: When both inner and outer layers are of FGM, Young's modulus, Poisson's ratio, and mass density for the bilayer cylindrical shell can be stated as

Results and Discussion
To justify the exactness and proficiency of the current results, simply supported-simply supported cylindrical shells are matched with [22,23] (see Table 1) and [24] (see Table 2). Table 3 is compared with [2]. ese comparisons show that the current outcomes are accurate, efficient, and valid. Ritz technique is utilized to accomplish the current result. is technique depends on the regulation of minimization of energy.

Mathematical Problems in Engineering 5
In the functionally graded (FG) cylindrical shell, FGM is composed of stainless steel (SS) and nickel (see Figure 2). Two edge conditions are applied on it, which are given as simply supported-simply supported and clamped-free. Table 4 represents the behavior of the NF (Hz) with n for a type I and II functionally graded C-F and SS-SS cylindrical shell. Here, SPs are m � 2, L � 20, N � 1, and R � 1. In both types, natural frequency (NF) rises with an increase in h. e natural frequency increases for different ascending values of h. ese values decrease from n � 1 to 2 and then start to decrease with the increasing values of n. e influence of NFs in types I and II is the same. e bilayer cylindrical shell is composed of two layers. In case I, the inner layer is made up of FGM, and the outer layer is of isotropic material. e FGM layer is made of nickel (Ni) and stainless steel (SS). In case II, both layers are composed   [2] Current [2] Current     of functionally graded material. In FGM, one layer is made up of nickel (Ni) and stainless steel (SS), while the other is of alumina (Al) and stainless steel (SS) (see Figures 3 and 4). Table 5 shows the behavior of the NF (Hz) with n for type-I, II, III, and IV functionally graded and isotropic bilayer cylindrical shells with both edge conditions (C-F and SS-SS). e NF decreases when the value of L rises. e natural frequency (NF) obtained from different values of L increases when the circumferential wave number n rises in all four types. Table 6 represents the behavior of the NF (Hz) with n for a type I and II FG bilayer CS. e NF rises when the value of power exponent N increases. e NF obtained from different values of power exponent N rises when n increases in both types.

Conclusion
In this work, vibration analysis of natural frequencies is made for the bilayered cylindrical shell. e shell layers are constructed from functionally graded materials (FGMs). Different FGMs are used to compose the bilayered cylindrical shell such as SS, Ni, and Al. e present analysis is performed under different boundary conditions. e shell vibration is based on different parameters of geometrical material.
e Lagrangian functional in the Rayleigh-Ritz technique has been used to solve the problem. In the current case, the axial mode is examined with the help of the characteristic beam function. Vibration of NFs is calculated for simply supported-simply supported and clamped-free edge conditions. Love's theory and Rayleigh-Ritz technique are applied to obtain the shell frequency equation. However, natural frequencies are examined with the effect of the volume fraction law. It is noticed that natural frequencies become high with the increase in thickness-to-radius ratios. e converse outcome is obtained with the increase of L/R ratios. Natural frequencies are increased with the increase in power law exponent N.
MATLAB software has been used for the simplification of the equation of eigenvalues. e current study can be prolonged to analyze the cylindrical shell (CS) submerged in the fluid and multilayered FGM cylindrical shell.

CS:
Cylindrical shell FGM: Functionally graded material SS-SS: Simply supported-simply supported C-F: Clamped-free NF: Natural frequency E: Young's modulus N: Power law exponent