Next Article in Journal
Knowledge-Guided Parallel Hybrid Local Search Algorithm for Solving Time-Dependent Agile Satellite Scheduling Problems
Next Article in Special Issue
Seismic Vulnerability Analysis of Concrete-Filled Steel Tube Tied Arch Bridges Using Symmetrically Arranged High-Damping Rubber Bearings
Previous Article in Journal
Energy Efficiency Measurement of Mechanical Crushing Based on Non-Contact Identification Method
Previous Article in Special Issue
Large Deflection Geometrically Nonlinear Bending of Porous Nanocomposite Cylindrical Panels on Elastic Foundation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Study of Four-Layered Cylindrical Shell Comprising Ring Support

by
Madiha Ghamkhar
1,†,
Ahmad N. Al-Kenani
2,† and
Naveed Hussain
1,*,†
1
Department of Mathematics and Statistics, University of Agriculture, Faisalabad 38000, Pakistan
2
Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2024, 16(7), 812; https://doi.org/10.3390/sym16070812
Submission received: 15 April 2024 / Revised: 14 June 2024 / Accepted: 16 June 2024 / Published: 28 June 2024
(This article belongs to the Special Issue Applied Mechanics, Engineering and Modeling - Volume II)

Abstract

:
In this work, the vibration analysis of a layered, cylinder-shaped shell is undertaken. The structure of the shell layers is composed of functionally graded and isotropic materials. The vibrations of four-layered cylindrical shells with a ring support along the axial direction are investigated in this research. The two internal layers are composed of isotropic materials, and the external two layers are composed of functionally graded materials. The outer functionally graded material layers considered are stainless steel, zirconia, and nickel. The inner two isotropic layers considered are aluminum and stainless steel. The shell frequency equation is acquired by employing the Rayleigh–Ritz method under the shell theory of Sanders. The trigonometric volume fraction law is used to sort the functionally graded material composition of the FGM layers. The natural frequencies are attained under two boundary conditions, namely simply supported–simply supported and clamped–clamped.

1. Introduction

In the analysis of structures, a cylindrical shell is a major feature. Various mechanical facets of these shell structures are analyzed in addition to their potential implementation, one of which is shell vibration. In areas of technology such as liquid gases at high pressure, nuclear power stations, networks of pipes, and other marine and aircraft systems, the examination of such shells plays a crucial role. The biography for the past scientific work has been developed based on Love’s theory. Utilizing Sanders’ shell theory and the Ritz composition, many researchers [1,2,3,4,5,6,7,8] have analyzed the frequency of dimensionally stable cylindrical shells under ring support over many end situations. Zang et al. [9] described the study of a hypothesis in which, contrary to Love’s results, all strains disappear from any rigid body. It offered formulations for the stress effects and couples that satisfy the homogeneous equilibrium conditions based on three stress frictions. The effect of the boundary conditions on the free vibrations for a multilayered cylindrical shell was investigated using the concept based on Love’s first theorem and the beam functions utilized for the Ritz as inertial structural features. There were nine boundary conditions taken into account, four of which had the same finishing criterion as their starting criterion. The other five conditions had different end requirements. Zang et al. [10] created heat-resistant structures for space plane fusion reactors and airframes. Zhi et al. [11] compared the frequency parameters with a precise three-dimensional linear elasticity analysis to test the analysis. The vibrational analysis of the functionally graded materials (FGMs) also discussed the behavior of natural frequencies.
Iqbal et al. [12] discussed the vibrational analysis of cylindrical shells and natural frequencies (NFs) under different boundary conditions and investigated new precise programs for the motions of the circular cylindrical shell with intermediate ring support based on the Goldenviezer–Novozhilov shell principle. The analytical method investigated the vibrational behavior of the ring support cylinder, while the state-space technique was used to operate the homogeneous comparative methodology for the shell component, and the domain decomposition technique was investigated for the continuity criteria among shell parts in the sense of the previously described method. Iqbal et al. [13] suggested 3D anisotropic elasticity fundamental equations, from which time visibility with differential equations is constructed in a coherent matrix form, and examined numerous techniques such as power law, sigmoid, and exponential FGMs, and we will focus on the numerical solution that was predominantly derived from the theoretical formulation and calculated using the MARC system. With reference to the effects of the application of P-FGM to S-FGM and E-FGM, Li et al. [14] introduced four sets of in-plane boundary conditions for the strictly assisted FG cylindrical shell for which the free vibration analysis was conducted. Sofiyev et al. [15] investigated physical properties that were thought to be temperature-dependent and eventually modified in the stress distribution of the shell. The results of the temperature increase were studied by specifying an approximate elevated temperature on the cylinder’s outer surface and the average temperature on the cylinder’s inner surface. The properties were temperature-dependent and gradually changed as the intensity of the shell increased.
By evaluating the elevated temperature, Arshad et al. [16] also mentioned the rising temperature effect. The entire study examined the strength of cylindrical shells supported on Winkler–Pasternak frameworks and made of ceramic, FGM, and metal layers that were exposed to an axial load. The revised stability and durability equations of the Donnell form on the Pasternak base were gathered. Galerkin’s process analysis values were obtained for the critical axial load of three-layered cylindrical shells with and without an elastic foundation. Extensive parametric experiments were performed to investigate the effects of FGM sheet thickness differences on the radius-to-thickness ratio. Naeem et al. [17] studied the vibration properties of two-layered cylindrical shells made of isotropic material on a single layer and functionally graded material on another. A bi-layered cylindrical shell’s vibration frequency analysis consisting of two distinct functionally graded layers is proposed by Arshad et al. [18]. The thicknesses of the shell’s layers are considered to be constant and equivalent. It is anticipated that the components of the material of the bi-layered dynamically graded cylindrical shell will differ seamlessly.
Ghamkhar et al. [19] investigated the vibration analysis of a three-layered cylindrical shell, with the anterior and posterior layers made of mechanically graded materials and the buffer part made of isotropic material. Love’s thin shell theorem was applied to the relationships of strain vs. displacement and curvature vs. displacement. To manipulate the intensity formulas, the Rayleigh–Ritz method (RRM) was used to grow the computational Lagrangian to be expressed as an eigenvalue query. The frequency wavelengths for long, simple-assisted end states were calculated by altering the thickness-to-radius and length-to-radius ratios; two non-dimensional geometrical parameters; and small, dense, and thin CSs. Ghamkhar et al. [20] examined a three-layered FGM cylindrical shell’s vibration frequency with an FGM central layer and homogeneous material as the interior and posterior layers. The relationships between the strain and the stress were drawn from the shell theory of Sanders. By utilizing the RRM, the shell frequency result was derived. The effect of natural frequencies (NFs) was evaluated for different middle-layer thicknesses. The beam functions of the features were used to approximate the axial modal dependency. The fundamental natural frequency of a cylindrical shell against the time span and thickness-to-radius percentage in a broad range was reported and investigated through this research.
Liu et al. [21] focused on determining the axial buckling load of stiffened multilayer cylindrical shell panels made of functionally graded graphene-reinforced composites (FG-GPL RCs). They used rings and stringers as stiffening tools. The study utilized Halpin–Tsai relations for elastic properties and implemented the virtual work principle and finite element approach based on the first-order shear deformation theory and Lekhnitskii’s smeared stiffener approach to derive stability equations. It examined four nanofiller dispersion types (FG-X, FG-A, FG-O, and UD) and systematically investigated the effects of geometric and material parameters on buckling loads and mode shapes, providing valuable design insights. The present research investigates the vibrations of four-layered cylindrical shells with ring support along the axial direction. The shell layers consist of a combination of functionally graded and isotropic materials. Specifically, the two internal layers are made of isotropic materials, while the two external layers are composed of functionally graded materials. The outer layers are assumed to be made of stainless steel, zirconia, and nickel, whereas the inner isotropic layers are composed of aluminum and stainless steel. The shell frequency equation is derived using the Rayleigh–Ritz method based on Sander’s shell theory. It utilizes energy principles (strain energy and kinetic energy) to derive the governing equations in an approximate form. The resulting system of algebraic equations is solved as an eigenvalue problem to determine natural frequencies and mode shapes.

2. Materials and Methods

A coordinate mechanism of cylinders x , θ , t is represented as the shell of the center reference surface with x, θ, and z as the axial, angular, and thickness coordinates, respectively, where u 1 x , θ , t , u 2 x , θ , t , and u 3 x , θ , t denote the deformations of the displacement functions in the longitudinal, tangential, and transverse instructions accordingly (see Figure 1).
The strain energy relation is given below
I = 1 2   0 L 0 2 π K S K R d θ d x
where
{ K } = { c 1 , c 2 , l , κ 1 ,   κ 2 ,   2 t }
The matrix [C] is represented as follows:
C = ξ 11 ξ 12 0 ξ 12 ξ 22 0 0 0 ξ 66 ψ 11 ψ 12 0 ψ 12 ψ 22 0 0 0 ψ 66 ψ 11 ψ 12 0 ψ 12 ψ 22 0 0 0 ψ 66 ζ 11 ζ 12 0 ζ 12 ζ 22 0 0 0 ζ 66
where c 1 , c 2 ,   l express the strains that consist of the reference surface, κ 1 , κ 2 , t indicate the curvatures, and ξ i j , ψ i j , and ζ i j (i, j = 1, 2, 6) represent the extensional, coupling, and bending stiffness [5] (Loy et al. 1999).
They are defined by the following formulas:
{ ξ i j , ψ i j , ζ i j } = H 2 H 2 Q i j 1 , z , z 2 d z
The reduced stiffness Q i j for isotropic materials is stated as,
Q 11 = E ( 1 λ 2 ) = Q 22 Q 12 = λ E 1 λ 2
And
Q 66 = E 2 1 λ 2
The matrix ψ i j = 0 is for isotropic circular-shaped CS, and ψ i j   0 is for FG CS; considering its structure and the characteristics of its constituents, the expressions from (2) and (3) in (1), E s t r a i n , are rewritten as:
E s t r a i n = 1 2   0 L 0 2 π { ξ 11 c 1 2 + ξ 22 c 2 2 + 2 ξ 11 c 1 c 2 + ξ 66 l 2 + 2 ψ 11 c 1 K 1 + 2 ψ 12 c 1 K 2 + 2 ψ 12 c 2 K 1 + 2 ψ 22 c 2 K 2 + 4 ψ 66 l t + ζ 11 K 1 2 + ζ 22 K 2 2 + 2 ζ 12 K 1 K 2 + 4 ζ 66 l 2 } R d θ d x
The expressions that follow are extracted from [19,20] and written as follows:
c 1 , c 2 , l = u 1 x , 1 R u 2 θ + u 3 , u 2 x + 1 R u 1 θ
  κ 1 , κ 2 , t = 2 u 3 x 2 , 1 R 2 2 u 3 θ 2 u 2 θ , 2 R 2 u 3 x θ 3 4 u 2 x + 1 4 R u 1 θ
By substituting expressions (7) and (8) into Equation (6), expression (6) becomes as follows:
E s t r a i n = 1 2   0 L 0 2 π [ ξ 11 u 1 x 2 + ξ 22 R 2 u 2 θ + u 3 2 + 2 ξ 12 R u 1 x u 2 θ + u 3 2 + ξ 66 u 2 x + 1 R u 1 θ 2 ψ 11 u 1 x 2 u 3 x 2 2 ψ 12 R 2 u 1 x 2 u 3 θ 2 u 2 θ 2 ψ 12 R u 2 θ + u 3 2 u 3 x 2 2 ψ 22 R 3 u 2 θ + u 3 2 u 3 θ 2 u 2 θ 4 ψ 66 R u 2 x + 1 R u 1 θ 2 u 3 x θ 3 4 u 2 x + 1 4 R u 1 θ + ζ 11 2 u 3 x 2 2 + ζ 22 R 4 2 u 3 θ 2 u 2 θ 2 + 2 d ζ 12 R 2 2 u 3 x 2 2 u 3 θ 2 u 2 θ + 4 ζ 66 R 2 2 u 3 x θ 3 4 u 2 x + 1 4 R u 1 θ 2 ] d x d θ
Shell kinetic energy is stated as
I = 1 2 0 L 0 2 π ρ t u 1 t 2 + u 2 t 2 + u 3 t 2 R d θ d x
where
ρ t = H 2 H 2 ρ   d z
where ρ is the mass density.
The Lagrange energy functional Γ for a CS is described as follows as a function of the strain and kinetic energies:
Γ   =   I E s t r a i n

Numerical Procedure

The RR method controls the normal frequencies of CS. The distortion defects in the longitudinal u 1 , transverse u 2 , and angular directions u 3 are taken into account. The displacement domains are now assumed by the following relationships:
  u 1 x , θ , t = x n D x cos m θ s i n ϖ t , u 2 x , θ , t = y n E x sin m θ c o s ϖ t
And
u 3 x , θ , t = z n F x cos m θ sin ϖ t
where x n , y n , and z n reflect the vibrational amplitudes in the x, θ, and t direction accordingly; the axial and circumferential wave numbers of mode shapes are indicated by n and m; ω symbolizes the angular vibration frequency of the shell wave; and   D ( x ) = d φ ( x ) d x , E ( x ) = φ ( x ) , and F ( x ) = φ ( x ) i = υ ( x d i ) l i , which indicates the axial function that satisfies the end point conditions. In the longitudinal direction at x = d i , i t h ring support is present. Here l i = 1 and l i = 0 represent with and without a ring support, respectively. m is the circumferential wave number and ω is the angular vibration frequency. And n denotes the axial model dependence in the longitudinal, circumferential, and transverse directions, respectively.
  • φ x symbolizes the axial function that meets the geometric edge condition requirement and is defined as
    φ x =   β 1 cosh μ m x + β 2 cos μ m x σ m β 3 sin h μ m x + σ m β 4 sin μ m x
μ m represents the roots of some transcendental equations, σ m parameters depend on the values of μ m , and β i   ( i = 1 , 2 , 3 , 4 ) change with respect to the edge conditions.
The geometric boundary conditions, namely simply supported and clamped boundary conditions, can be expressed mathematically in terms of the characteristic beam function φ x as follows:
Simply supported boundary conditions φ x = φ x = 0 .
Clamped boundary conditions φ x = φ x = 0 .
The following non-dimensional parameters are used to expand this problem.
U = D ( x ) H ,   V = E ( x ) H ,   W = F ( x ) R   a i j = ξ i j H   ,   b i j = ψ i j H 2   ,   d i j = ζ i j H 3 a = R / L , b = H / R , X = x / L , ρ t = ρ t / H
Expression (13) alters the structure as follows:
u 1 ( x , θ , t ) = H x n U cos ( m θ ) sin ( ϖ t ) u 2 ( x , θ , t ) = H y n V sin ( m θ ) cos ( ϖ t ) u 3 ( x , θ , t ) = R z n W cos ( m θ ) sin ( ϖ t )
The Lagrangian function £ m a x is transformed by employing the principle of maximum energy.
Γ max = π H L R 2 R 2 ϖ 2 ρ t 0 1 b 2 x n U 2 + b 2 y n V 2 + b 2 z n W 2 d X 0 1 a 2 b 2 a 11 x n d U d X 2 + a 22 m b y n V + z n W 2 2 a 3 b 2 b 11 x n d U d X z n d 2 W d X 2 + 2 a b a 12 x n d U d X m b y n V + z n W + a 66 a b y n d V d X m b x n U 2 2 a b 2 b 12 x n d U d X m 2 z n W + m b y n V 2 a 2 b b 12 m b y n V + z n W z n 2 d 2 W d X 2 2 b b 22 m b y n V + z n W m 2 z n W + m b y n V + a 4 b 2 d 11 z n 2 d 2 W d X 2 2 4 b b 66 a b y n d V d X + m b x n U n a z n d W d X 3 a b y n 4 d V d X + m b 4 x n U + b 2 d 22 m 2 z n W + m b y n V 2 + 2 a 2 b 2 d 12 z n 2 d 2 W d X 2 m 2 z n W + m b y n V + 4 d 66 m a z n d W d X 3 a b y n 4 d V d X + m b 4 x n U 2 d X
Lagrangian energy functional Γ max is extremized with respect to vibrational amplitudes x n , y n , and z n , as shown below:
Γ max x n = Γ max y n = Γ max z n = 0
The resulting relation is defined in the form of a matrix as
K Ω 2 M X = 0
where
Ω 2 = R 2 ω 2 ρ t ,
K and M are the stiffness and mass matrices corresponding to the cylindrical shell.

3. Classification of Material

In the current study, a CS is taken into consideration, which is made up of four layers: isotropic material is used to construct the internal layers, while nickel, zirconia, and stainless steel are used to construct the outermost layers(see Figure 2).
The two elements using the trigonometric volume fraction law (VFL) are as follows:
V f 1 = sin 2 z H 1 H 2 H 1 N , V f 2 = cos 2 z H 1 H 2 H 1 N , 0 N
This relation fulfills the VFL, i.e., V f 1 + V f 2 = 1 , where h is the shell’s thickness and N indicates the power law exponent. Every layer is assumed to be the same thickness h / 4 . The material parameters are as follows: E 1 , λ 1 , ρ 1 for nickel; E 2 , λ 2 , ρ 2 for zirconia; E 3 , λ 3 , ρ 3 for stainless steel; E 4 , λ 4 , ρ 4 for aluminum; E 5 , λ 5 , ρ 5 for stainless steel and zirconia. Then, the effective material quantities are as follows:
E f g m 1 = E 1 E 2 sin 2 z H 1 H 2 H 1 N + E 2 λ f g m 1 = λ 1 λ 2 sin 2 z H 1 H 2 H 1 N + λ 2 ρ f g m 1 = ρ 1 ρ 2 sin 2 z H 1 H 2 H 1 N + ρ 2 E f g m 2 = E 4 E 5 sin 2 z H 3 H 4 H 3 N + E 5 λ f g m 2 = λ 4 λ 5 sin 2 z H 3 H 4 H 3 N + λ 5 ρ f g m 2 = ρ 4 ρ 5 sin 2 z H 3 H 4 H 3 N + ρ 5
For the expression above at z = H 2 , H 2
E f g m 1 = E 2 , λ f g m 1 = λ 2 , ρ f g m 1 = ρ 2 and
E f g m 2 = E 5 , λ f g m 2 = λ 5 , ρ f g m 2 = ρ 5 , respectively.
The stiffness moduli are altered as follows:
ξ i j = ξ i j F G M + ξ i j i s o + ξ i j i s o + ξ i j F G M , ψ i j = ψ i j F G M + ψ i j i s o + ψ i j i s o + ψ i j F G M , ζ i j = ζ i j F G M + ζ i j i s o + ζ i j i s o + ζ i j F G M
where i . j = 1 , 2 , 6 , and iso signify the two internal isotropic layers, and FGM symbolizes the two outer FGM layers.

4. Results and Discussion

To explain the accuracy and knowledge of the exact outcome, the current outcome is used by FGM and isotropic CS for SS edge conditions. RRM is used to achieve the current result. Table 1 represents the comparison of frequency parameters for simply supported isotropic cylindrical shells (n = 1, L = 20, R = 1, H = 0.01, and λ = 0.3) with those in [4]. Three materials are used to make two types of FG materials known as FGM cylindrical layers: stainless steel, zirconia, and nickel. Two FGM variations for the form of substance are found. Table 2 represents the material distribution.
  • Type 1:
In type 1 FGM CS, inner surfaces are made of nickel and zirconia; on the other hand, outer layers are composed of stainless steel and zirconia.
  • Type 2:
In type 2 FGM CS, inner surfaces are made with zirconia and nickel, and the outer layers are composed of zirconia and stainless steel.
In Table 3, the variation in natural frequencies (NFs) is presented under SS to SS edge conditions observed with (n = 1, R = 1, d = 0.01, H = 0.002, and L = 20). As the value of circumferential wave number m increases from 1 to 10, the value of NFs decreases. For power law exponent ranging from N = 1 to N = 9, the values of NFs also decrease, but slowly. For values ranging from N = 1 to 20, the values change from 366.3000 to 363.0547.
In Table 4, the variation in NFs for length-to-radius ratios is presented against the circumference wave number m for shell I with parameters (n = 2, R = 1, d = 0.01, H = 0.002, and L = 20) under the S-S Edge condition. As the value of m increases from 1 to 10, the values of NFs change, and for N = 1 to 11, the values of NFs change from 362.6197 to 359.2090.
In Table 5, the variation is the same as in Table 3 with parameters (n = 3, R = 1, d = 0.03, H = 0.002, and L = 90).
Table 6 represents the variation in NFs for length-to-radius ratios against m for shell I with parameters (n = 4, R = 1, d = 0.04, H = 0.02, and L = 20) under the S-S Edge condition. NFs increase up to 10. The value of NFs decreases from 4755.773 to 4575.089.
Configurations of shells:
The composition of the four-layered cylindrical shell made by FGM and the isotropic type of material is defined below in Table 2.
Table 7 describes the variation in NFs for length-to-radius ratios against m for shell I with parameters (n = 1, d = 0.4, R = 1, N = 5, and L = 20) under the S-S edge condition, and it was observed that NFs change from 4493.442 to 798.6292 for m.
Table 8 presents the variation in NFs for length-to-radius ratios against m for the shell I with parameters (n = 1, R = 1, N = 5, and L = 20). It was observed and noted that the values of NFs change from 4313.561 to 792.0536 for m = 1 to m = 10. And for exponential law H = 0.001 to H = 0.07, the values change from 4313.561 to 4314.682.
Table 9 shows the variation in NFs for length-to-radius ratios against m for the shell I with parameters (n = 1, R = 1, N = 10, d = 0.06, and L = 20) under the S-S Edge condition, and it is noted that NFs vary from 4280.537 to 4280.562.
Table 10 shows NFs for length-to-radius ratios against m for shell I with parameters (n = 1, R = 1, d = 0.7, N = 20, and L = 20). It is observed that its variation for m = 1 to m = 10 changes from 425.7616 to 425,7638. It has been noticed that a very small change occurs throughout the region.
Table 11 shows NFs for length-to-radius ratios against circumference wave number n for shell I with parameters (n = 1, R = 1, N = 1, H = 0.005, and d = 0.8) under the S-S Edge condition by changing the value of H. The value of NFs changes from 418.758 to 594.9993 as n increases from 1 to 10. As the thickness of the shell varies, NFs remain nearly the same for different thicknesses from 0.001 to 0.07. It is noticed that for small changes in thickness, the values of NFs remain the same. Another noticeable variation for L is that the value of NFs increases when the range is from L = 10 to L = 70, increasing from 418.758 to 439.2311, indicating an increase in the frequency value.
Figure 3 represents the variation in natural frequencies against circumferential wave number m for type 1 shell. The lower line shows the behavior of NFs for N = 1 and the upper line presents the behavior of NFs for N = 2. Natural frequencies increased swiftly from m = 1 to 2 and then increased slowly. Figure 4 represents the same for N = 1 and N = 3 under simply supported end point conditions for type 2 shell.
Figure 5 presents the behavior of NFs against m for type 1 four-layered cylinder-shaped shell. The natural frequencies are obtained under simply supported edge conditions. Natural frequencies rapidly decreased from m = 1 to 2 and then slowly increased for m = 3, after which they behaved smoothly. The results are acquired for thicknesses 0.001 and 0.005 of type 1 shell. In Figure 6, the same behavior of NFs is observed for thicknesses 0.001, 0.005, and 0.007 of type 2 shell.
Figure 7 represents the variation in NFs against m of type 2 SS-SS four-layer CS with n = 2, d = 0.9, N = 1, R = 1, and H = 0.005. It is noticed that NFs are very close for L = 10, 20, and 40. The results increased rapidly from m = 1 to 2 and then increased slowly with increasing m.
Figure 8 shows the vibrations of NFs (Hz) for a four-layered clamped–simply supported circular cylinder-shaped shell against the position of the ring support with different H/R ratios ( m = 10 , n = 1 , L / R = 20 , H / R = 0.001 , 0.01 , 0.05 ) .
Figure 8 describes the vibrations of NFs (Hz) for a four-layered free circular cylinder-shaped shell against the position of the ring support with different H/R ratios ( n = 10 , m = N = 1 , L = 20 , H / R = 0.001 , 0.01 , 0.05 ) . The extreme lower line is obtained for H/R = 0.001, the extreme upper line is obtained for H/R = 0.01, and the middle line represents H/R = 0.05. It is clearly seen that NFs behave symmetrically as they increase from d = 0 to 0.5, attain their maximum value at d = 0.5, and then decrease from d = 0.05 to d = 1.

5. Discussion

The study presented in this work delves into the intricate vibration behavior and natural frequencies of layered cylindrical structures composed of functionally graded and isotropic materials. By investigating the vibrations of four-layered cylindrical shells assisted by a ring over their length, this research aims to uncover valuable insights into the dynamics of such composite structures. One of the distinctive features of this study is the composition of the layered cylindrical shells. The internal two layers are constructed by using isotropic materials, specifically aluminum and stainless steel, while the external two layers consist of functionally graded materials (FGM), namely stainless steel, zirconia, and nickel. This combination offers a diverse range of mechanical properties, making it an intriguing subject for exploration. The research involves a cylindrical shell with each of its four layers having an equal thickness, denoted as h/4.
To effectively analyze the vibration behavior and natural frequencies of these layered cylindrical shells, the Rayleigh–Ritz method (RRM) is employed within the framework of Sander’s shell theory. This method allows for the derivation of the shell frequency equation, providing valuable insights into the system’s dynamic response. Furthermore, the research investigates the natural frequencies of the shells under two distinct boundary conditions: simply supported–simply supported (SS-SS) and clamped–clamped (C-C). By exploring these boundary conditions, the study aims to provide a thorough understanding of how different support configurations impact the vibrational characteristics of the cylindrical shells.
In summary, this research contributes significantly to the understanding of vibration behavior and natural frequencies in layered cylindrical structures composed of functionally graded and isotropic materials. By leveraging advanced analytical techniques and exploring various material compositions and boundary conditions, the study sheds light on the complex dynamics of such composite systems, with potential implications for a wide range of engineering applications.

6. Conclusions

The current study was conducted to check the vibration characteristics of thin circular CS formed from FGM. The FG cylindrical shell constituents are made of stainless steel, nickel, zirconia, and aluminum. These four types were studied by changing the configuration of the layers of FG material. These constituents in type 1 of the cylindrical shells vary across the shell’s thickness constantly, seamlessly, and steadily from the FG layer’s inner to its outermost layer. The basic objectives were to obtain shell frequency equations with the help of these theories and techniques. Sander’s shell theory was applied to solve the shell frequency equations based on Kirchhoff’s assumptions. RRM was used to solve the shell vibration problems by utilizing the Lagrangian energy relation to obtain the shell’s frequency equation in numerical form. Axial modal dependence was approximated by trigonometric functions for simply supported end conditions and characteristic beam functions for number-clamped boundary conditions.
The study found that the minimum frequency occurs at specific wave numbers for all boundary conditions and is influenced directly by increasing the value of N, L, n, d, H, and L/R. The NFs decreased by varying N under different conditions. All the results were very close to each other. This analysis may be extended for the study of different shell problems, such as changing parameters like the thickness or length of the shell and the position of the ring. The frequency curves of the shell with ring support at different positions obtain symmetric shapes because of the same edge conditions. They are not symmetrical around the center because of different end point conditions. The induction of ring support on a cylinder-shaped shell has a significant effect on the NFs compared to the shell frequencies without ring support.
For the presented cylindrical shell, the results were obtained using MATLAB software version 2019 under simply supported–simply supported and clamped–clamped boundary conditions. These results illustrate the natural frequency behavior as a function of the circumferential wave number. Significant findings are as follows: The natural frequency increases as the power law index N″ increases. An increase in the length L″ of the shell leads to an increase in the natural frequency. Increasing the shell thickness h″ results in higher natural frequencies. Varying the ring support position along the longitudinal path shows that the natural frequency increases, reaching its maximum when the ring is at the central position of the shell, and then decreases as the ring position moves beyond the center.
These findings highlight the critical influence of geometric and material parameters on the vibrational behavior of cylindrical shells, providing insights for optimizing design and performance.

Author Contributions

Conceptualization, M.G. and N.H.; methodology, M.G., A.N.A.-K. and N.H.; validation, A.N.A.-K., M.G. and N.H.; writing draft, editing, M.G., A.N.A.-K. and N.H.; review, M.G., A.N.A.-K. and N.H.; funding, A.N.A.-K.; supervision, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Niino, M.; Hirai, T. Functionally gradient materials. In pursuit of super heat-resisting materials for spacecraft. J. Jpn. Soc. Compos. Mater. 1987, 13, 257–264. [Google Scholar]
  2. Lam, K.Y.; Loy, C.T. Effects of boundary conditions on frequencies of a multi-layered cylindrical shell. J. Sound Vib. 1995, 188, 363–384. [Google Scholar] [CrossRef]
  3. Swaddiwudhipong, S.; Tian, J. Vibrations of cylindrical shells with intermediate supports. J. Sound Vib. 1995, 187, 69–93. [Google Scholar] [CrossRef]
  4. Loy, C.T.; Lam, K.Y. Vibration of cylindrical shells with ring support. Int. J. Mech. Sci. 1997, 39, 455–471. [Google Scholar] [CrossRef]
  5. Loy, C.T.; Lam, K.Y. Vibration of functionally graded cylindrical shells. Int. J. Mech. Sci. 1999, 41, 309–324. [Google Scholar] [CrossRef]
  6. Pradhan, S.C.; Loy, C.T. Vibration characteristics of functionally graded cylindrical shells under various boundary conditions. Appl. Acoust. 2000, 61, 111–129. [Google Scholar] [CrossRef]
  7. Xiang, Y.; Ma, Y.F. Exact solutions for vibration of cylindrical shells with intermediate ring supports. Int. J. Mech. Sci. 2002, 44, 1907–1924. [Google Scholar] [CrossRef]
  8. Chen, W.Q.; Bian, Z.G. Three-dimensional vibration analysis of fluid-filled orthotropic FGM cylindrical shells. Int. J. Mech. Sci. 2004, 46, 159–171. [Google Scholar] [CrossRef]
  9. Zhang, L.; Xiang, Y. Vibration of open circular cylindrical shells with intermediate ring supports. Int. J. Solids Struct. 2006, 43, 3705–3722. [Google Scholar] [CrossRef]
  10. Zhang, L.; Xiang, Y. Local adaptive differential quadrature for free vibration analysis of cylindrical shells with various boundary conditions. Int. J. Mech. Sci. 2006, 48, 1126–1138. [Google Scholar] [CrossRef]
  11. Zhi-Yuan, C.; Hua-Ning, W. Free vibration of FGM cylindrical shells with holes under various boundary conditions. J. Sound Vib. 2007, 306, 227–237. [Google Scholar] [CrossRef]
  12. Iqbal, Z.; Naeem, M.N. Vibration characteristics of FGM circular cylindrical shells using wave propagation approach. Acta Mech. 2009, 208, 237–248. [Google Scholar] [CrossRef]
  13. Iqbal, Z.; Naeem, M.N. Vibration characteristics of FGM circular cylindrical shells filled with fluid using wave propagation approach. Appl. Math. Mech. 2009, 30, 1393–1404. [Google Scholar] [CrossRef]
  14. Li, S.R.; Fu, X.H. Free vibration of three-layer circular cylindrical shells with functionally graded middle layer. Mech. Res. Commun. 2010, 37, 577–580. [Google Scholar] [CrossRef]
  15. Sofiyev, A.H.; Avcar, M. The stability of cylindrical shells containing an FGM layer subjected to axial load on the Pasternak foundation. Engineering 2010, 2, 228. [Google Scholar] [CrossRef]
  16. Arshad, S.H.; Naeem, M.N. Vibration of bilayered cylindrical shells with layers of different materials. J. Mech. Sci. Technol. 2010, 24, 805–810. [Google Scholar] [CrossRef]
  17. Naeem, M.N.; Khan, A.G. Vibration of three-layered FGM cylindrical shells with a middle layer of isotropic material for various boundary conditions. World J. Mech. 2014, 4, 315–331. [Google Scholar] [CrossRef]
  18. Arshad, S.H.; Naeem, M.N. Influence of ring support on free vibration of sandwich functionally graded cylindrical shells with a middle layer of isotropic material. J. Eng. Res. 2016, 4, 9. [Google Scholar] [CrossRef]
  19. Ghamkhar, M.; Naeem, M.N. Vibration frequency analysis of three-layered cylinder-shaped shell with the effect of FGM central layer thickness. Sci. Rep. 2019, 9, 1566. [Google Scholar] [CrossRef]
  20. Ghamkhar, M.; Naeem, M.N. Vibration analysis of a three-layered FGM cylindrical shell including the effect of ring support. Open Phys. 2019, 17, 587–600. [Google Scholar] [CrossRef]
  21. Liu, Z.; Tornabene, F.; Dimitri, R.; Babaei, M. Numerical Study of the Buckling Response of Stiffened FG Graphene-Reinforced Multilayer Composite Cylindrical Panels. Processes 2024, 12, 430. [Google Scholar] [CrossRef]
Figure 1. Structure of CS with a ring support that is located along the length.
Figure 1. Structure of CS with a ring support that is located along the length.
Symmetry 16 00812 g001
Figure 2. Component materials of shell layers.
Figure 2. Component materials of shell layers.
Symmetry 16 00812 g002
Figure 3. Variation in NFs against m of type 1 SS-SS four-layer CS (n = 1, d = 0.1, H = 0.002, R = 1, and L = 20).
Figure 3. Variation in NFs against m of type 1 SS-SS four-layer CS (n = 1, d = 0.1, H = 0.002, R = 1, and L = 20).
Symmetry 16 00812 g003
Figure 4. Variation in NFs against m of type 2 SS-SS four-layer CS (n = 2, d = 0.1, H = 0.002, R = 1, L = 20).
Figure 4. Variation in NFs against m of type 2 SS-SS four-layer CS (n = 2, d = 0.1, H = 0.002, R = 1, L = 20).
Symmetry 16 00812 g004
Figure 5. Variation in NFs against m of type 1 SS-SS four-layer CS (n = 1, d = 0.4, N = 1, R = 1, L = 20, and H = 0.001, 0.005).
Figure 5. Variation in NFs against m of type 1 SS-SS four-layer CS (n = 1, d = 0.4, N = 1, R = 1, L = 20, and H = 0.001, 0.005).
Symmetry 16 00812 g005
Figure 6. Variation in NFs against m of type 2 SS-SS four-layer CS (n = 1, d = 0.5, N = 5, R = 1, L = 20, and H = 0.001, 0.005, and 0.007).
Figure 6. Variation in NFs against m of type 2 SS-SS four-layer CS (n = 1, d = 0.5, N = 5, R = 1, L = 20, and H = 0.001, 0.005, and 0.007).
Symmetry 16 00812 g006
Figure 7. Variation in NFs against m of type 2 SS-SS four-layer CS (n = 2, d = 0.9, N = 1, R = 1, and H = 0.005).
Figure 7. Variation in NFs against m of type 2 SS-SS four-layer CS (n = 2, d = 0.9, N = 1, R = 1, and H = 0.005).
Symmetry 16 00812 g007
Figure 8. Vibrations of NFs (Hz) for a three-layered free circular cylinder-shaped shell against the position of the ring support with different H/R ratios ( n = 10 , m = N = 1 , L = 20 , H / R = 0.001 , 0.01 , 0.05 ) .
Figure 8. Vibrations of NFs (Hz) for a three-layered free circular cylinder-shaped shell against the position of the ring support with different H/R ratios ( n = 10 , m = N = 1 , L = 20 , H / R = 0.001 , 0.01 , 0.05 ) .
Symmetry 16 00812 g008
Table 1. Comparison of frequency parameters for simply supported isotropic cylindrical shell (n = 1, L = 20, R = 1, H = 0.01, and λ = 0.3).
Table 1. Comparison of frequency parameters for simply supported isotropic cylindrical shell (n = 1, L = 20, R = 1, H = 0.01, and λ = 0.3).
m[4]PresentDifference
10.0161010.0161000.006%
20.0093820.0093760.06%
30.0221050.0221010.02%
40.0420950.0420920.007%
50.0680080.0680060.009%
Table 2. Configurations of shell types.
Table 2. Configurations of shell types.
Types of
Shell
1st Layer
FGM
2nd Layer
Isotropic
3rd Layer
Isotropic
4th Layer
FGM
Type 1Nickel, ZirconiaStainless SteelAluminumStainless Steel,
Zirconia
Type 2Zirconia, NickelStainless SteelAluminumZirconia,
Stainless Steel
Table 3. Vibration of NFs (Hz) against n for FGM CS under SS-SS boundary conditions when (n = 1, L = 20, H = 0.02, R = 1, and d = 0.1).
Table 3. Vibration of NFs (Hz) against n for FGM CS under SS-SS boundary conditions when (n = 1, L = 20, H = 0.02, R = 1, and d = 0.1).
mN = 1N = 2N = 5N = 7N = 9
1366.3000364.0547363.0547363.0547363.0547
2456.1014454.0393452.3383452.3383452.3383
3477.2445475.0875473.3106473.3106473.3106
4485.2193483.0266481.2239481.2239481.2239
5489.6059487.3932485.5785485.5785485.5785
6493.1564490.9264489.1031489.1031489.1031
7497.0930494.8423493.0085493.0085493.0085
8502.1957499.9165498.0668498.0668498.0668
9509.0963506.7773504.9033504.9033504.9033
10518.3668515.9934514.0839514.0839514.0839
Table 4. Vibration of NFs (Hz) against n for FGM CS under SS-SS boundary conditions when (n = 2, L = 20, H = 0.02, R = 1, and d = 0.1).
Table 4. Vibration of NFs (Hz) against n for FGM CS under SS-SS boundary conditions when (n = 2, L = 20, H = 0.02, R = 1, and d = 0.1).
mN = 1N = 3N = 5N = 7N = 9N = 11
1362.6197359.9578359.2090359.2090359.2090359.2090
2455.8980452.9469452.1372452.1372452.1372452.1372
3477.2085474.1213473.2755473.2755473.2755473.2755
4485.2101482.0728481.2149481.2149481.2149481.2149
5489.6049486.4406485.5777485.5777485.5777485.5777
6493.1596489.9726489.1063489.1063489.1063489.1063
7497.0994493.8851493.0148493.0148493.0148493.0148
8502.2050498.9528498.0760498.0760498.0760498.0760
9509.1086505.8025504.9154504.9154504.9154504.9154
10518.3821515.0016514.0989514.0989514.0989514.0989
Table 5. Behavior of NFs (Hz) against n for FGM CS under SS-SS boundary conditions when (n = 1, L = 20, N = 1, R = 1, and d = 0.4).
Table 5. Behavior of NFs (Hz) against n for FGM CS under SS-SS boundary conditions when (n = 1, L = 20, N = 1, R = 1, and d = 0.4).
mH = 0.001H = 0.005H = 0.007H = 0.01H = 0.05H = 0.07
14493.4424493.4434493.4444493.4454493.464493.469
2785.8470785.7706785.7388785.6991786.0900786.9266
3794.0083793.8696793.8343793.8241798.5644804.2985
4796.3821796.2296796.2626796.4488814.3950833.7276
5797.3917797.3464797.5922798.2960844.7577891.3993
6797.9048798.1816798.8778800.6169897.6645988.5901
7798.1968799.1243800.6216804.1492979.43441131.8350
8798.3810800.4223803.2025809.5400979.43441131.8350
9798.5141802.2857806.9778817.43771242.72681560.6494
10798.6292804.9192812.3117828.50701425.32111842.6299
Table 6. Behavior of NFs (Hz) when (n = 1, L = 20, N = 5, R = 1, and d = 0.5).
Table 6. Behavior of NFs (Hz) when (n = 1, L = 20, N = 5, R = 1, and d = 0.5).
mH = 0.001H = 0.005H = 0.007H = 0.01H = 0.05H = 0.07
14313.5614314.6574314.6584314.6594314.6734314.682
2779.3320838.7948838.7603838.7167839.0724839.9035
3787.4543838.7012838.6600838.6395843.0766848.5493
4789.8100838.6410838.6587838.8150855.6305873.9499
5790.8129838.6853838.8951839.5266883.1276927.3678
6791.3239838.9250839.5478841.1367932.47241019.1568
7791.6158839.4714840.8337844.08301009.73671156.0430
8791.8012840.4554843.0066848.87241009.73671156.0430
9791.9362842.0268846.3550856.06971262.09341571.3251
10792.0536844.3534851.1977866.28191438.92591846.7592
Table 7. Behavior of NFs (Hz) against n for FGM cylindrical shell when (n = 4, H = 0.005, N = 1, R = 1, and d = 1).
Table 7. Behavior of NFs (Hz) against n for FGM cylindrical shell when (n = 4, H = 0.005, N = 1, R = 1, and d = 1).
mL = 10L = 20L = 30L = 40L = 50L = 70
1461.9686329.8816306.7302.7992302.8946304.7281
2420.4381387.3319385.3606385.4567385.6814385.9822
3415.7960406.0455405.6847405.7264405.7787405.8421
4417.1420413.4404413.3143413.3238413.3382413.3561
5419.1609417.3701417.2904417.2846417.2859417.2891
6421.3530420.2377420.1592420.1424420.1365420.1324
7423.9559423.0649422.9713422.9443422.9328422.9233
8427.3108426.4492426.3323426.2948426.2780426.2637
9431.7691430.8393430.6941430.6454430.6232430.6041
10437.6788436.6267436.4497436.3891436.3612436.3372
Table 8. Variation in NFs (Hz) against n for FGM cylindrical shell under the SS-SS boundary condition (n = 1, R = 1, H = 0.002, L = 20, and d = 0.1).
Table 8. Variation in NFs (Hz) against n for FGM cylindrical shell under the SS-SS boundary condition (n = 1, R = 1, H = 0.002, L = 20, and d = 0.1).
mN = 1N = 2N = 3N = 4N = 5N = 10
1356.7561353.1512351.5667350.1086348.8782348.4389
2451.3025447.1279445.3168443.6651442.2836441.7932
3472.5032468.1340466.2384464.5095463.0634462.5501
4480.1950475.7544473.8275472.0697470.5992470.0773
5483.8605479.3859477.4437475.6717474.1890473.6626
6485.9757481.4816479.5303477.7496476.2593475.7301
7487.4275482.9202480.9625479.1754477.6793477.1480
8488.6148484.0970482.1339480.3413478.8402478.3069
9489.7555485.2281483.2599481.4619479.9556479.4204
10490.9919486.4545484.4807482.6771481.1654480.6281
Table 9. Variation in NFs (Hz) against n for FGM cylindrical shell under SS-SS when (n = 2, R = 1, H = 0.002, L = 20, and d = 0.2).
Table 9. Variation in NFs (Hz) against n for FGM cylindrical shell under SS-SS when (n = 2, R = 1, H = 0.002, L = 20, and d = 0.2).
mN = 1N = 2N = 3N = 4N = 5N = 10
1419.8248413.2517410.2073407.3157404.8068403.8956
2542.1721537.1515534.9731532.9863531.3244530.7345
3563.7859558.5723556.3103554.2473552.5216551.9091
4571.2834566.0004563.7079561.6167559.8672559.2462
5574.7881569.4725567.1653565.0602563.2989562.6735
6576.7816571.4475569.1315567.0180565.2491564.6210
7578.1269572.7805570.4583568.3385566.5639565.9337
8579.2045573.8485571.5211569.3958567.6161566.9839
9580.2178574.8530572.5207570.3902568.6054567.9712
10581.2970575.9233573.5858571.4497569.6594569.0231
Table 10. Variation in NFs (Hz) against n for FGM cylindrical shell under the SS-SS boundary condition (N = 1, R = 1, n = 1, H = 0.005, and d = 0.8).
Table 10. Variation in NFs (Hz) against n for FGM cylindrical shell under the SS-SS boundary condition (N = 1, R = 1, n = 1, H = 0.005, and d = 0.8).
mL = 10L = 20L = 40L = 50L = 70
1413.1000419.8248428.9463431.2599434.1054
2540.9617542.1822542.5790542.6293542.6737
3563.6582563.8588563.9190563.9265563.9331
4571.4879571.5395571.5550571.5570571.5587
5575.4332575.4426575.4459575.4463575.4467
6578.1811578.1709578.1688578.1686578.1683
7580.7581580.7346580.7289580.7283580.7277
8583.7200583.6847583.6760583.6750583.6741
9587.4637587.4167587.4051587.4037587.4025
10592.3291592.2701592.2554592.2536592.2521
Table 11. Variation in NFs (Hz) against n for FGM CS under SS-SS boundary condition (N = 10, R = 1, n = 1, H = 0.005, and d = 0.9).
Table 11. Variation in NFs (Hz) against n for FGM CS under SS-SS boundary condition (N = 10, R = 1, n = 1, H = 0.005, and d = 0.9).
mL = 10L = 20L = 40L = 50L = 70
1343.7526348.8778354.1291355.0526355.9365
2441.4266442.2936442.6029442.6428442.6781
3462.9772463.1435463.1979463.2047463.2108
4470.8447470.8886470.9031470.9049470.9065
5474.9309474.9370474.9399474.9403474.9406
6477.8692477.8566477.8541477.8538477.8536
7480.7127480.6867480.6805480.6798480.6791
8484.0531484.0148484.0055484.0043484.0034
9488.3237488.2733488.2608488.2593488.2580
10493.8992493.8362493.8205493.8186493.8170
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ghamkhar, M.; N. Al-Kenani, A.; Hussain, N. Structural Study of Four-Layered Cylindrical Shell Comprising Ring Support. Symmetry 2024, 16, 812. https://doi.org/10.3390/sym16070812

AMA Style

Ghamkhar M, N. Al-Kenani A, Hussain N. Structural Study of Four-Layered Cylindrical Shell Comprising Ring Support. Symmetry. 2024; 16(7):812. https://doi.org/10.3390/sym16070812

Chicago/Turabian Style

Ghamkhar, Madiha, Ahmad N. Al-Kenani, and Naveed Hussain. 2024. "Structural Study of Four-Layered Cylindrical Shell Comprising Ring Support" Symmetry 16, no. 7: 812. https://doi.org/10.3390/sym16070812

APA Style

Ghamkhar, M., N. Al-Kenani, A., & Hussain, N. (2024). Structural Study of Four-Layered Cylindrical Shell Comprising Ring Support. Symmetry, 16(7), 812. https://doi.org/10.3390/sym16070812

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop