Open access peer-reviewed chapter

Vibration Characteristics of Fluid-Filled Functionally Graded Cylindrical Material with Ring Supports

Written By

Muzamal Hussain, Aamir Shahzad, Muhammad Nawaz Naeem and Maogang He

Submitted: 07 June 2017 Reviewed: 02 November 2017 Published: 14 February 2018

DOI: 10.5772/intechopen.72172

From the Edited Volume

Computational Fluid Dynamics - Basic Instruments and Applications in Science

Edited by Adela Ionescu

Chapter metrics overview

971 Chapter Downloads

View Full Metrics

Abstract

Vibration analysis of fluid-filled functionally graded material (FGM) cylindrical shells (CSs) is investigated with ring supports. The shell problem is formulated by deriving strain and kinetic energies of a vibrating cylindrical shell (CS). The method of variations of Hamiltonian principle is utilized to change the shell integral problem into the differential equation (DE) expression. Three differential equations (DE) in three unknown for displacement functions form a system of partial differential equations (PDEs). The shells are restricted along the thickness direction by ring supports. The polynomial functions describe the influence of the ring supports and have the degree equal to the number of ring supports. Fluid loaded terms (FLT) are affixed with the shell motion equations. The acoustic wave equation states the fluid pressure designated by the Bessel functions of first kind. Axial modal deformation functions are specified by characteristic beam functions which meet end conditions imposed on two ends of the shell. The Galerkin method is employed to get the shell frequency equation. Natural frequency of FGM cylindrical shell is investigated by placing the ring support at different position with fluid for a number of physical parameters. For validity and accuracy, results are obtained and compared with the data in open literature. A good agreement is achieved between two sets of numerical results.

Keywords

  • functionally graded material (FGM)
  • ring supports
  • cylindrical material
  • Galerkin technique (GT)
  • Hamiltonian principle

1. Introduction

All over the word, applications of fluid-filled cylindrical shells have grown in engineering and science. Amendments in shell physical quantities are inducted to enhance strength and stability of cylindrical shells (CSs) [1, 2]. Additional burden due to fluid factor on a physical system may cause damage to it. In the recent times, this feature has appealed to scientists doing research on dynamical properties of materials to explore more about specific strength, stiffness and super corrosion resistance [3, 4, 5, 6, 7]. It has been acquired by highly developed complex materials. Study of vibratory response of cylindrical shells (CSs) containing fluids is very beneficial to study dynamic behavior for their applications. This presents a direct contact between a solid composition and a liquid material [8]. In a lot of fields of engineering and technology (mechanical, civil, aeronautics), its useful implications can be seen. Thin-walled cylindrical shells (CSs) have extensive applications in engineering and industry. They are found in chimney design, pipe flow, nuclear reactors and submarines.

For theoretical point of view, study of cylindrical shell vibrations is done to investigate analytical results and their closeness with experimental ones. Here the Galerkin procedure is employed to solve the shell governing equations. For the present cylindrical shells, functionally graded materials are utilized for their structure construction. In the radial thickness direction, material distribution is handled by the exponential volume fraction law. Due to this law, special types of integrals are evolved and are approximated numerically or analytically to evaluate material stiffness modulus. These integrals involve the material parameters of thickness variable by assuming the Poisson ratios of functionally graded constituent materials. For simplifying the integrals, these are presumed to be nearly equal to each other. This assumption simplifies the material stiffness integrals. Shell dynamical equations are framed by applying the Hamilton’s variational principle to the Lagrangian functional that is obtained from the shell strain and kinetic energy expressions. These equations govern the shell vibration behavior. For this problem, a suitable and effective method is employed to achieve the shell frequency equation in the eigenvalue problem expression. Normally energy variational approaches are used to solve cylindrical shells problem. They consist of the Raleigh - Ritz method and the Galerkin method. The axial deformation functions are estimated by characteristic beam functions. They are achieved from the solutions of beam differential equation.

Pioneering research work on vibrations of cylindrical shells (CSs) has performed by Arnold and Warburton [8]. The consequence of end conditions on vibration characteristics of a circular cylindrical shell (CS) was considered by Fosberg in 1964 by using shell equation.

Najafizadeh and Isvandzibaei [5] analyzed vibration characteristics of functionally graded cylindrical shells with ring supports. They based their analysis higher order shear deformation theory of shells. It was perceived that the influence of ring supports and fluid terms was very significant on shell frequencies. Vibration characteristics of cylindrical shells containing fluid were studied experimentally and theoretically by Chung et al. [9]. Goncalves and Batista [10] presented a theoretical vibration study of cylindrical shells (CSs) partially filled and submerged in a fluid. Simply supported end conditions were imposed on both edges. Goncalves et al. [11] investigated the transient stability of empty and fluid-filled cylindrical shells and used to study the non-linear dynamic behavior of shallow cylindrical shells under axial loading. Gasser [12] studied the frequency spectra of bi-layered cylindrical shells by taking different materials in both layers such as isotropic as well as functionally graded material (FGM) and by taking two different FGM at the inner and outer layers of the CSs respectively. Sharma and Johns [13] explored the vibrations of CSs with clamped-free and clamped-ring-stiffeners conditions by applying the Raleigh-Ritz method and estimated the axial displacement deformation with the help of beam functions. Xi et al. [14] have studied vibrations of cross-ply plastic-coated circular fluid-filled CSs by applying a semi-analytical method based on Reissner-Mindline theory. Zhang et al. [4] studied vibrations of CSs and applied the wave propagation approach (WPA) to solve shell dynamical equations. This method depended on the eigenvalue of characteristic beam functions. Axial wave number was designated to a boundary condition (BC) by a simple formula. They compared the results determined by this method to ones found by a FEM to check the efficiency, robustness and accuracy of the procedure. It was seen by making these comparisons that their approach is more victorious and exact for shell vibration difficulties. It was concluded that the proposed approach could be useful for a problem with compound end states and also for fluid-filled cylindrical shells.

Zhang et al. [15] examined vibrations of CSs containing fluid by applying WPA. After that a similarity evaluation was conducted between uncoupled frequencies with the numerical outcome obtained in the literature. They also put side by side the coupled frequencies estimated by the WPA with those attained by FEM. Xiang et al. [16] accessed exact solution for the vibration characteristics of CS placed at intermediate position and used the domain decomposition technique for the sake of ordering in the segment of the shells. Zhao et al. [17] investigated the effects of vibration with ring stiffeners and stringer for the laminated cross-ply rotating CS and two methods: variational method and averaging approach are used for these effects. They determined that averaging method produced the inexact values and was sensitive whereas the fast and better results were deduced with variational method. Xiang et al. [18] accessed the vibration characteristics of CS placed at intermediate position with axially dense ring supports and used Flügge shell theory and the Timoshenko thin shell theory to analyze the buckling shells as composite materials. Due to FGM their composition vary constantly and smoothly through thickness.

Vibration characteristic of FGM shell with ring supports has investigated by Isvandzibaei and Awasare [19] and they used third order deformation shear theory and Hamilton’s principle for free-free end. Lee and Chang [20] gave a numerical study of coupled problems of fluid conveying dual walled carbon annotates and examined the effects of characteristic ratio and Van der Waals forces on basic frequencies. Silva et al. [21] investigated the nonlinear vibrant behavior and instabilities of partially fluid-filled CS constrained to axial load and resulting in a distinct low-dimensional model for the analysis of the vibrations to observe the shell vibration. Shah et al. [22] gave a vibration analysis of a functionally graded CS containing a fluid. The shells were rested on elastic foundations. They analyzed effects of Winkler and Pasternak moduli on shell vibration characteristics. Xiong et al. [23] investigated the free vibration analysis of fluid-filled elliptical cylindrical shells and explained the sensitivity of frequency parameters to the elliptical parameter with length of CS. The cylinder is filled with a compressible non-viscous fluid and may be subjected to arbitrary time-harmonic on-surface mechanical drives is investigated by Hasheminejad and Alaei-Varnosfadrani [24]. The free vibration of fluid-filled CS covered partially in elastic foundation is investigated by Kim [25] and the elastic foundation of partial axial and angular dimensions is represented by the Pasternak model. The variation of the frequency parameters with respect to the layer thickness, the length-to-radius ratio, the length-to-thickness ratio, and circumferential node number are analyzed by Izyan et al. [26]. Soutis et al. [27] investigated influence of ring supports on free vibration of FGM which is placed on the middle layer and Study is carried out for placing ring support in different position of FGM’s, to find the natural frequencies by Rayleigh–Ritz approach.

In the present paper, vibration frequency characteristics of fluid-filled CSs are investigated. The shells are constrained in the radial direction by ring supports. The present problem is formulated in integral form and is converted into a system of three partial differential equations (PDEs) with the unknown displacement functions. Modal forms for the three unknown functions are assumed such that the special and temporal variables are separated. Energy variation approach is used to solve PDEs so that an eigenvalue problem is cropped up. Axial modal dependence is roughly estimated by trigonometric functions for a simply supported CSs. For other end conditions, characteristic beam functions are taken. The radial constraints are presumed by the polynomial functions having degree equal to number of ring supports. Fluid pressure is stated by the acoustic wave equation and Bessel’s functions of first kind. Axial modal displacement deformations are measured by beam characteristic functions which ensure to meet boundary conditions. The Galerkin technique is implemented to form the shell frequency equation which is solved by using MATLAB coding. The radial deflection is restrained by ring supports. This factor is expressed by the polynomial functions which carry the degree equal to the number of ring supports.

Advertisement

2. Formation of shell problem

2.1. Functionally graded shells

In practice a CS is constructed from a FGM which consists of two constituent materials. Two constituent materials having material parameters: E1, E2, ν1, ν2, and ρ1, ρ2. Then the effective material quantities: Efgm, vfgm and ρfgm are given as:

Efgm=E1E2zh+12p+E2,νfgm=ν1ν2zh+12p+ν2,ρfgm=ρ1ρ2zh+12p+ρ2E1

The value of z lies as 0<z< in the radial direction. The volume fraction Vr for a functionally graded constituent material can be defined by the following function:

Vr=zh+12pE2

where p is the power law exponent which indicates the material variation profile through the shell thickness and keeps its real between zero and infinity.

2.2. Theoretical investigation

Consider Figure 1 in which a geometrical sketch of a CS is given. L, R, and h are termed as shell geometrical parameters. Other shell basic quantities are material parameters and are designate by E, v and ρ.

Figure 1.

Geometrical sketch of a FG cylindrical shell with ring support.

where the forces and moments are designated by N and M have the directions along the longitudinal, tangential and shear directions correspondingly.

NxNΨN=h2h2σxσΨσdz,MxMΨM=h2h2σxσΨσzdzE3

where σx, σΨ are the linear stresses along x and Ψ-directions respectively and σ represents the shear stress along xΨ-direction. For a cylindrical shell, the stresses defined in Eq. (3) are defined by the two dimensional Hook’s law.

σxσΨσ=Q11Q120Q12Q22000Q66exeΨeE4

where the strains along x and Ψ directions are labeled by ex and eΨ respectively and the shear strain is denoted by e in the xΨ - direction. The first thin shell theory was developed by Love [7] which is based on Kirchhoff’s perception for plates. Various thin shell theories have been deduced from this theory by modifying the geometrical and physical parameters. The members of the strain vector e in Eq. (4) have been defined as linear functions of thickness coordinate z which taken from Love’s [32] theory are stated as:

ex=e1+zκ1,eΨ=e2+zκ2,e=γ+2E5

here e1, e2 and γ denote strains with regard to the shell middle reference surface. κ1, κ2 and τ stand for the surface curvatures. The expressions for strain and curvature displacement relationship are written as:

e1e2γ=ux1RvΨ+wvx+1RuΨ
κ1κ2τ=2wx21R22wΨ2vΨ1R2wxΨvxE6

By substituting Eqs. (5) and (6) into Eq. (4) and then substituting the resulting equation into Eq. (3). The force and moment results can be written as:

NxNΨNMxMΨM=A11A120B11B120A12A220B12B22000A6600B66B11B120D11D120B12B220D12D22000B6600D66e1e2γκ1κ22τE7

where Aij, Bij and Dij (i, j = 1, 2 and 6) are the extentional, coupling and bending stiffness defined respectively, as:

AijBijDij=h2h2Qij1zz2dz,E8

where coupling stiffness, Bij’s vanish for a CS structured from isotropic materials where they exist for heterogeneous and an isotropic materials such as laminated FGM. For an isotropic material, Qij (i, j = 1, 2 and 6) are expressed as

Q11=Q22=E1v2,Q12=vE1v2,Q66=E21+vE9

where E and ν are the Young’s modulus and Poisson’s ratio for the shell’s material. In this study, the cylindrical shell is considered thin and valid with thickness-to- radius ratio is less than 0.05. For vibrating thin cylindrical shell, the strain energy expressed as

S=R20L02π[A11e12+A22e22+2A12e1e2+A66γ2+2B11e1k1+2B12e1k2+2B12e2k1+2B22e2k2+4B66γτ+D11k12+D22k22+2D12k1k2+4D66τ2]dΨdxE10

Substituting the expression for the surface strains and the curvatures from the relationships (6)

S=R20L02πA11ux+A22R2vΨ+w2+2A12RuxvΨ+w+A66vx+1RuΨ22B11ux2wx22B12Rux2wxΨvx2B12RvΨ+w2wx22B22R3vΨ+w2wΨ2vΨ8B66Rvx+1RuΨ2wxΨvx+D112wx22+D22R42wΨ2vΨ2+2D12R22wx22wΨ2vΨ+4D66R22wxΨvx2dΨdxE11

Also for a cylindrical shell, its kinetic energy expression represented by T, is given as:

T=R20L02πρTut2+vt2+wt2dΨdxE12

where ρT is expressed as:

ρT=h2h2ρdzE13

Now the shell problem is framed by the Lagrangian energy functional which is the difference between the shell kinetic and strain energies and is given as:

=TSE14

where denotes the Lagrangian functional. Substituting the expressions for strains and kinetic energies of the shell from Eqs. (11) and (12) respectively into Eq. (14). Calculus of variations process is applied to the integral terms to derive the Euler-Lagrange equations. Hamilton’s variational principle is a process in which the variations in the variables are assumed to be zero. Implementation this principle to the Lagrangian functional which is an integral expression. The subsequent dynamical equations are obtained in the following system of PDEs:

A112ux2+A66R22uΨ2+A12+A66R+B12+2B66R22vxΨ+A12RwxB113wx3B12+2B66R23wxΨ2=ρT2ut2
A12+A66R+B12+B66R22uxΨ+A66+3B66R+3D66R22vx2+A22R2+2B22R3+D22R42vΨ2B12+2B66R+D12+2D66R23wx2Ψ+A22R2+B22R3wΨB22R3+D22R43w3Ψ=ρT2vt2
B113ux3A12Rux+B12+2B66R23uxΨ2+B12+2B66R+D12+4D66R23vx2Ψ+B22R3+D22R43vΨ3A22R2+B22R3vΨD114wx42D12+2D66R24wx2Ψ2D22R44wΨ4+2B12R2wx2+2B22R32wΨ2A22R2w=ρT2wt2E15
Advertisement

3. Application of Galerkin technique

Two energy variational techniques namely viz., the Rayleigh-Ritz method and the Galerkin technique (GT) are exploited to solve shell motion equations because these methods yield results fast with enough accuracy. The present CS is analyzed by applying the Galerkin technique (GT) for their vibrations. This approach is very expedient, simple and convenient to use to find vibration frequencies and has been widely engaged by numerous mathematicians [25, 28, 29]. Modal displacement forms are designated by x, Ψ and t. The following modal deformation displacement functions for u, v and w are adopted as:

uxΨt=pmdxsincosωt
vxΨt=qmϕxcoscosωt
wxΨt=rmϕxi=1kxaizisincosωtE16

Here the parameters pm, qm and rm present vibration amplitudes in the x, Ψ and z directions respectively. The position of ith ring support with the circular direction of the shell is denoted by ai and zi has value, 1 when a ring support exists and is 0, when no rings hold up. For this purpose the modal displacement forms for u,v and w given in the relation (12) respectively and their corresponding partial derivatives are substituted into Eq. (15) by taking, zi=1 for a single ring support, the resulting equations are integrated with respect to x from 0 to L, the following equations are got:

A11I1n2A66R2I2pmnA12+A66R+B12+2B66R2I2qm+A12RI6+I7B11I8+3I9+n2B12+2B66R2I6+I7rm=ω2ρTI2pm
nA12+A66R+B12+2B66R2I3pm+A66+3B66R+4D66R2I3n2A22R2+2B22R3+D22R4I4qm+nA22R2+B22R3I10+n3B22R3+D22R4I10nB12+2B66R+D12+4D66R2I11+2I12rm=ω2ρTI4pm
A12RI18+B11I19n2B12+2B66R2I18pm+nA22R2+B22R3I10+n3B22R3+D22R4I10nB12+2B66R+D12+4D66R2I18qm+A22R2I3+2B12RI14+2I152n2B22R3I13D11I16+4I17+2n2D12+2D66R2I14+2I15D22R4I13rm=ω2ρTI13rmE17

where the integral terms are listed in Appendix-I. Terms in Eq. (17) are arranged to form the homogeneous algebraic linear equations (HALEs) in pm, qm and rm. This leads to the formation of the shell frequency equation the eigenvalue shape as:

d11pm+d12qm+d13rm=ω2ρTpmI2d21pm+d22qm+d23rm=ω2ρTqmI4d31pm+d32qm+d33rm=ω2ρTrmI13E18

So the above equations are written in the eigenvalue problem notation as:

d11d12d13d21d22d23d31d32d33pmqmrm=ω2ρTI2000I4000I13pmqmrmE19

3.1. Annexation of fluid terms

The acoustic pressure generated by a fluid is described by the wave equation in cylindrical coordinate system xΨr and is presented as:

1rrrqr+1r22qΨ2+2qx2=1c22qt2E20

where t, q, c represents respectively the time, acoustic pressure and speed of sound in the fluid. The acoustic pressure produced in the fluid meets the equation of motion Eq. (20) and is articulated by the following modal function expression:

q=qmsinJnkrrΨxcosωtE21

Here Jnkrr represents the Bessel’s function of first kind with order n. It is the same number as the circumferential wave number. ω denotes the natural frequency for the CS, kr stands for the radial wave number and axial wave number km is designated by for a number of boundary conditions that has been indicates in Ref. [10]. There exists a relation between km and kr is written as:

krR2=Ω2cLcf2kmR2E22

where cL and cf denote the speeds of the sound in the empty and fluid-filled cylindrical shells respectively. Ω is the non-dimensional frequency parameter. An assumption is made that the fluid and the shell radial displacement must be the same at the borderline of the interior wall of the shell and it is made sure that the fluid is kept with the interaction between the fluid and the shell wall. A coupling condition relation of a fluid with the shell wall exists and is defined by the following expression:

1iωρfqr=wtE23

By applying this condition at r=R, the fluid loading term (FLT) owing to the existence of the fluid pressure is given by

FLT=ω2ρfJnkrRkrJnkrRrmE24

JnKrr denotes the differentiation of the Bessel’s function with respect to the argument (krR). The fluid loaded term (FLT), which represents the fluid pressure, is annexed with the frequency Eq. (18) for an empty CS. Ultimately the shell frequency equation for fluid-filled functionally graded CS is articulated in the following forms:

d11d12d13d21d22d23d31d32d33pmqmrm=ω2ρTI2000I4000I13+ρfkrJnkrRJnkrRpmqmrmE25

The expressions for the terms dij's,I2,I4 and I13 are given in Appendix-I. This is an eigenvalue problem involving the shell frequency. Presently, the MATLAB is generally used to calculate the physical problems in engineering and science. In our case, MATLAB computer software has been used to compute the shell frequencies through eigenvalues and eigenvectors. A single command ‘eig’ furnishes shell frequencies and mode shapes by calculating eigenvalues and eigenvectors respectively. It is significant that MATLAB can compute integrals that fairly easy to apply the Galerkin method with polynomial basis functions. The exact solution of integral equation is known and presents a sample of MATLAB code to illustrate the success of the method. In this way the Galerkin method is implemented to form the shell frequency equation that is solved by employing MATLAB software.

3.2. Effective material

Materials of cylindrical shells have a paramount role in analyzing shell vibrations. They impress their stability. In practice isotropic, laminated and functionally graded materials are used to manufacture them. Here functionally graded materials are benefitted to form the shells. These materials are advanced and useful in a highly environs. Their material properties are temperature - dependents. One of their material properties C is stated by the following relation:

C=C0C1T1+C1T+C2T2+C3T3E26

where C1,C0,C1,C2 and C3 designate the temperature dependent constants. T is calculated in the Kelvin scale. These constants differ from material to material. Formula (26) is due to Toulokian [38]. A functionally graded cylindrical shell comprising of two constituent materials can be classified into two categories. This depends upon the arrangement of the two materials forming the shell. It is known that the stainless steel and nickel are used for structuring such types of shells. Stainless steel and nickel are used in its external and internal surfaces respectively for Category- I (C-I) CS structure, while for Category-II (C-II) cylindrical shell, stainless steel and nickel are taken for constituting its internal and external surfaces respectively. At temperature 300 K, the material properties for stainless steel and nickel of functionally graded cylindrical shell are: E, v, ρ for nickel are 2.05098×1011N/m2, 0.31, 8900 kg/m3 and stainless steel are 2.07788×1011N/m2, 0.317756 and 81,666N/m3. These values have been taken from Ref. [3].

Advertisement

4. Numerical results

In this section, to check the validity and accuracy, for the determination of the natural frequencies with present methodology, empty and fluid-filled cylindrical shells with ring supports are analyzed and the results are compared with experimental and other numerical values found in literature. The stainless steel and nickel are used for structuring such types of shells. Stainless steel and nickel are used in its external and internal surfaces respectively for Category- I (C-I) CS structure, while for Category-II (C-II) cylindrical shell, stainless steel and nickel are taken for constituting its internal and external surfaces respectively. At temperature 300 K, the material properties for stainless steel and nickel of functionally graded cylindrical shell are: E, v, ρ for nickel are 2.05098×1011N/m2, 0.31, 8900 kg/m3 and stainless steel are 2.07788×1011N/m2, 0.317756 and 81,666N/m3.

A few comparisons of analytical frequencies for isotropic cylindrical shells are exhibited to validate the present Galerkin technique. As a first example, the lowest dimensionless frequency Ω=ωR1ν2ρ/E of a simply supported empty cylinder are compared with the solution derived using generalized differential quadrature method (DQM) by Loy et al. [2] as shown in Table 1. They are varied with circumferential wave mode, n for axial wave number m = 1. There is an excellent agreement is seen between two sets of frequency parameters results as the percentage difference is negligible.

nLoy et al. [2]PresentDifference %
10.0161010.0161020.006
20.0093820.0093830.010
30.0221050.0221060.004
40.0420950.0420970.004
50.0680080.0680090.001
60.0997300.0997320.002
70.1372390.1372410.001
80.1805270.1805280.000
90.2295940.2295960.000
100.2844350.2844360.000

Table 1.

Comparison of frequency parameter Ω=ωR1v2ρ/E for a simply supported CS (m=1, L/R=20, h/R=0.01, ν=0.3).

As another example, the natural frequency of a simply supported empty cylinder is compared with Loy et al. [3], Ansari et al. [30] and Warburton [35] as shown in Figure 2(a). They are varied with axial wave number, m for circumferential wave mode, n = 2, 3. For the same shell, the present results for a fluid-filled shell are compared with those numerical results obtained by Gonçalves and Batista [10], Gonçalves et al. [11] and obtained experimentally by Gasser [12] in Figure 2(b). They are varied with circumferential wave mode, n for axial wave number, m = 1. There is an excellent agreement is seen between two sets of frequency parameters results.

Figure 2.

(a) Comparison of natural frequencies (Hz) for a simply supported CS (L = 8 in, h = 0.1 in, R = 2 in, E = 30 × 106 lbf in−2, ν = 0.3, ρ = 7.35 × 10−4 lbf s2 in−4). (b) Comparison of natural frequencies (Hz) for a fluid-filled CS (m = 1, L = 0.41 m, h = 0.001 m, R = 0.3015 m, E = 2.1 × 1011 N/m2, v = 0.3, ρ = 7850 kg/m3).

When n = 2, 3, m = 1, the present values are 2046.4, 2199.0 and the values of Loy et al. [3] are 2050.7, 2204.0 which is good agreement between each other. Figure 2(b) exhibits a good coincidence between respective counters of the frequencies. The values at n = 8, 9 are less significant and on enhancing the value of n, the present numerical values are smaller than those of Gonçalves and Batista [10] and Gonçalves et al. [11] and Gasser [12]. This difference is the result of two separate analytical techniques.

4.1. Frequency analysis of fluid-filled cylindrical shell

Figure 3(a) reveals variations of natural frequencies (Hz) for a CS including fluid against L/R for circumferential wave modes n = 1, 2, 3, 4, 5 and the longitudinal wave mode, m = 1. The value of length-to-radius ratio is composed as 1L/R20, the values of L/R=1, at CWN: n = 1, 2, 3, 4, 5 are 325.625, 252.019, 151.816, 132.818, 99.128 and L/R=20 are 6.2054, 2.3505, 2.8798, 6.3196, 10.0948. For the wave mode, n = 1, 2, 3, 4, 5 as L/R is made to grow i.e., the shell becomes longer and longer, the frequency reduces for each tangential wave mode.

Figure 3.

Variations of natural frequency (a) for fluid-filled CS against L/R. (m = 1, h = 0.004 m, R = 1 m, E = 2.1 × 1011 N/m2, v = 0.3, ρ = 7850 kg/m3, ρf = 1000 kg/m3) (b) for fluid-filled CS with h/R (m = 1, L = 8 m, R = 1 m, E = 2.1 × 1011 N/m2, v = 0.3, ρ = 7850 kg/m3, ρf = 1000 kg/m3).

Now the values of length-to-radius ratio 1L/R20, at CWN: n = 1, the frequencies are 325.675, 220.745, 99.005, 33.693, 16.348, 10.770, 6.2054 and at CWN: n = 5 are 99.128, 33.017, 13.756, 10.3603, 10.1506, 10.1148, 10.0948 respectively. These results show that, on increasing the values of L/R the frequencies reduce to each CWN mode [36]. Figure 3(b) demonstrates variations of natural frequencies (Hz) for a CS enclosing fluid versus the h/R for the tangential wave modes, n = 1, 2, 3, 4, 5 and the longitudinal wave mode, m = 1. The value of height-to-radius ratio is composed as 0.002h/R0.05, the values of h/R = 0.002, at CWN: n = 1, 2, 3, 4, 5 are 3.36929, 4.1911, 5.0052, 4.5812, 5.4979 and h/R=0.05, are 5.6977, 19.6873, 35.4274, 79.3621, 126.739 respectively. As h/R is made to increase i.e., the shell gets thicker and thicker, the frequency increases. Now the values of height-to-radius ratio 0.002h/R0.05, at CWN: n = 1, the frequencies are 3.36929, 3.7929, 4.6937, 4.9046, 5.6977 and at CWN: n = 5 are 5.4979, 10.3603, 50.7402, 76.0671, 126.7139. As h/R is made to increase i.e., these results shows that the frequency increases very slowly for n = 1 but for higher values of n, there is seen an appreciable increments in frequency values with each tangential wave mode. Previous study reveals that when the value of h/R increases then frequencies also increases [37].

4.2. Frequency analysis of fluid-filled cylindrical shells with ring support

In this section, frequency analysis for an isotropic CS enclosing fluid is performed by attaching some ring supports. These rings are located at some distance from one end of a cylindrical shell as shown in Figure 4.

Figure 4.

FGM shell with two ring supports (Ansari et al. [30]).

In Figure 5(a) variations of natural frequencies (Hz) for a CS containing fluid with ring supports versus L/R are demonstrated for three positions of the ring supports at a = 0.3L, 0.5L, L. The value of L/R is taken as 1L/R20, when L/R at ring support a = 0.3L, 0.5L, L without fluid, the frequencies are 790.31, 863.378, 745.085 and L/R = 20 the values are 420.639, 516.775, 208.819 and now with fluid at L/R = 1 the values are 352.975, 362.764, 314.385 and L/R = 20, the values are 187.549, 217.132, 86.109 respectively. There is a substantial decrease in frequency values as the shell length gets higher [36]. Influence of fluid on vibration frequency is seen significantly visible that reduced them approximately to (40~50%). In Figure 5(b) natural frequencies (Hz) for a CS containing fluid without and with ring supports are exhibited with h/R. The value of h/R is taken as: 0.002h/R0.05, when h/R = 0.002 at ring support a = 0.3L, 0.5L, L without fluid are 413.2180, 548.4688, 215.4291 and h/R = 0.05 the values are 413.3080, 548.4888, 215.5019 respectively. Now with fluid at h/R = 0.002 the values are 273.6543, 365.6459, 144.5833 and h/R = 0.05, the values are 273.7139, 365.6592, 144.6321 respectively. The ring supports are located at a = 0.3L, 0.5L, L. It is observed that frequencies enhance minutely for the thicker CSs [37]. Again the frequency has been considerably reduced when the fluid is added. From Figure 5(a) and (b), it is observed that in both the cases, without fluid and with fluid for L/R and h/R, the value of ring support a = 0.3L which is sandwich between a = 0.5L, L. It is seen that the influences of ring supports and fluid terms are converse to each other. The ring supports increase the frequencies whereas the fluid loaded terms lower them [32, 33, 34].

Figure 5.

Variations of natural frequencies (Hz) for a WOF and WF cylindrical shell with ring support against (a) L/R (b) h/R (n = 1, m = 1, L = 8 m, h = 0.004 m, R = 1 m, E = 2.1 × 1011 N/m2, v = 0.3, ρ = 7850 kg/m3, ρf = 1000 kg/m3).

4.3. Frequency analysis of empty and fluid-filled cylindrical shells with ring supports

It is noted that the natural frequencies of the shells have varied by the location of ring support in the shell and this change also varies WOF and WF as shown in Figure 6(a)(d). It is evident from these figures that when the values of a 0<a<0.5 increases then the natural frequencies also increases and a=0.5L, it reaches its peak value but for 0.5L<a<0.1L, on increasing the value of a, it begins to decrease and rust itself as bell shape symmetric curve. It is clear seen that when the ring support is positioned at the center of the fluid-filled cylindrical shell, the natural frequency attains its extreme value; it is observed that the natural frequency decreases as the ring supports shifts from center toward right/left side of the fluid-filled cylindrical shell. Thus, when different exponent is adopted, the frequency curve is also symmetrical about the center of fluid-filled cylindrical shell for simply supported symmetric end condition imposed on the both shell ends.

Figure 6.

Variation of natural frequencies (Hz) (a) against the positions of the ring supports for WOF and WF cylindrical shell. (n = 1, m = 1, L = 8 m, h = 0.004 m, R = 1 m, E = 2.1 × 1011 N/m2, v = 0.3, ρ = 7850 kg/m3, ρf = 1000 kg/m3) (b) versus the positions of ring supports for WOF and WF cylindrical shell with ring support: Category-I (n = 1, m = 1, L = 8 m, h = 0.004 m, R = 1 m) (c) Category-II (d) Category-I and II.

In Figure 6(a), variations of natural frequencies (Hz) with locations of the ring supports are listed with ring supports for WOF and WF. The ring support is composed as 0aL and for the values at a=0 for WOF and WF are 215.429, 140.665 and at a=L the values are 215.429, 146.095 respectively. At 0.4L<a<0.6L means a=0.5L for WOF and WF have a peak value 548.468, 548.358 but for a=0.4L the values are 518.609, 496.419 and a=0.6L are 518.609, 504.140 very closed to each other. In each case the shell frequencies goes up with the position of the ring supports to the highest values at the mid of the shell and then start to lower down to gain their initial values. It is noticed that the fluid addition has made frequency to decrease. From the previous data, the shell frequencies are affected highly as the fluid quantities [11, 31, 32, 33, 34, 35] and ring supports are appended [2, 5, 27, 36]. A beam type vibration of CSs crops up when the addition of ring supports are made.

Figure 6(b) and (c) display the comparisons of variations of natural frequencies (Hz) with ring supports for WOF and WF for C-I and C-II respectively. The ring support is composed as 0aL and for volume fraction law, the exponent p = 0.5, 1, 15 is adopted, the values at a=0 for WOF and WF are 207.058, 135.199 (p = 0.5), 205.346, 134.082 (p = 1), 201.012, 131.252 (p = 15) and at a=L the values are 207.058, 140.46 (p = 0.5), 205.346, 139.255 (p = 1), 201.012, 136.316 (p = 15) for C-I and for C-II: 204.048, 133.234 (p = 0.5), 205.729, 134.332 (p = 1), 210.311, 137.324 (p = 15) and at a=L the values are 204.048, 138.375 (p = 0.5), 205.729, 139.515 (p = 1), 210.311, 142.622 (p = 15) respectively. For C-I, at 0.4L<a<0.6L means a=0.5L have a extreme value of frequency 523.924, 523.819 (p = 0.5), 519.591, 519.487 (p = 1), 508.626, 508.524) (p = 15), for both WOF and WF but for a=0.4L the values of frequencies are 496.188, 474.957 (p = 0.5), 492.085, 471.030 (p = 1), 481.699, 461.088 (p = 15) and a=0.6L are 496.188, 482.345 (p = 0.5), 492.085, 478.356 (p = 1), 481.699, 468.260 (p = 15) very closed to each other for both WOF and WF. For C-II: at 0.4L<a<0.6L means a=0.5L have a maximum frequency value 513.823, 513.720 (p = 0.5), 518.059, 517.955 (p = 1), 529.595, 488.154 (p = 15) for both WOF and WF but for a=0.4L the values are 487.212, 466.366, 491.226, 470.208 (p = 1), 502.164, 480.678 (p = 15) and a=0.6L are 487.212, 473.619, 491.226, 477.521 (p = 1) 502.164, 488.154 (p = 15) very closed to each other for both WOF and WF. In both cases, the frequencies first rise, attain their maximum at the mid position and then fall down to original values at the other shell end. There is observed a reduction in frequency values owing to fluid term induction [11, 32, 33, 34]. It is concluded that the maximum point and other relative value of C-II are bit smaller than that of corresponding value of C-I. This conclusion is also shown in Figure 6(d) where C-I and C-II are illustrated in one graph. In Figure 6(d) variations of frequencies for two categories with ring supports for fluid-filled for C-I and C-II are cataloged with the locations of the ring supports. The ring support is composed for fluid-filled as 0aL, for volume fraction law, the exponent p = 0.5, 1, 15 is adopted, the values at a=0 are 135.99, 133.234 (p = 0.5), 134.082, 132.332 (p = 1) 131.252, 130.324 (p = 15) and at a=L for category-I and for category-II, the values are 140.416, 138.375 (p = 0.5), 139.255, 139.515 (p = 1), 136.316, 133.622 (p = 15). For C-I, C-II at 0.4L<a<0.6L means a=0.5 have a maximum value 513.819, 513.720 (p = 0.5), 519.487, 517.955 (p = 1), 508.524, 505.489 (p = 15) but for a=0.4L the values are 474.957, 466.366 (p = 0.5), 471.030, 470.208 (p = 1), 461.088, 450.678 (p = 15) and a=0.6 are 482.342, 473.342 (p = 0.5), 478.356, 477.521 (p = 1), 468.260, 466.154 (p = 15) respectively. Here for C-I cylindrical shells, frequencies are some bit higher for those of C-II ones [32, 33, 34]. The effects of fluid-filled cylindrical shell with ring supports are perceived to be very conspicuous in the variation of natural frequencies.

Advertisement

5. Summary

Theoretical vibration analysis of cylindrical shells (CSs) has a significant importance in applied mathematics and mechanics in view of their practical uses. Here the shell problem has been associated with investigation vibrations of cylindrical shells with ring supports. They shells have been supposed to be constructed from functionally graded materials. They materials are advanced and smart for their physical properties. Moreover these shells have been assumed to contain fluid. Here the Galerkin technique is employed to obtain the shell frequency equation. Effect of ring supports on shell vibration is inducted by a polynomial function which has degree equal to the number of ring supports. From the study of results, it is observed that vibration frequencies of cylindrical shells decrease significantly when the fluid loaded terms are appended. However their variations are alike to that behavior which is noticed for the cylindrical shells not containing fluid and attachment of ring supports boost much the vibration frequencies. For functionally graded material cylindrical shells, variations of frequencies with the circumferential wave modes for length-to-radius ratio, height-to-radius ratio for fluid-filled cylindrical shell and also for fluid-filled cylindrical shell with ring supports has been analyzed. It is seen that the influences of ring supports and fluid terms are converse to each other. The ring supports increase the frequencies whereas the fluid loaded terms lower them. However, the increments and decrements in the shell frequency depend upon the order of functionally graded material constituents forming a cylindrical shell. The induction of fluid-filled with ring support on the cylindrical shell has a prominent effect on the natural frequency as compared to the shell frequency without fluid attached with ring support. Different position of ring supports with and without fluid with various exponent law for category-I, II or both are analyzed. It is concluded that in both cases, the frequencies first rise, attain their maximum values at the mid position and then fall down to original values at the other shell end. Moreover, there is observed a reduction in frequency values owing to fluid term induction. It is also concluded that the maximum point and other relative value of category-II are bit smaller than that of corresponding value of category-I An extension of this analysis procedure can be performed to investigate vibrations of rotating fluid isotropic and functionally graded cylindrical shells with ring supports.

Advertisement

d11=A11I1n2A66R2I2
d12=nA12+A66R+B12+2B66R2I2
d13=A12RI6+I7B11I8+3I9+n2B12+2B66R2I6+I7
d21=nA12+A66R+B12+2B66R2I3
d22=A66+3B66R+4D66R2I3n2A22R2+2B22R3+D22R4I4
d23=nA22R2+B22R3I10+n3B22R3+D22R4I10nB12+2B66R+D12+4D66R2I11+2I12
d31=A12RI18+B11I19n2B12+2B66R2I18
d23=nA22R2+B22R3I10+n3B22R3+D22R4I10nB12+2B66R+D12+4D66R2I18
d33=A22R2I3+2B12RI14+2I152n2B22R3I13D11I16+4I17+2n2D12+2D66R2I14+2I15D22R4I13

where

I1=0Ld3ϕdx3dxdxI2=0LdxdxdxI3=0Ld2ϕdx2ϕxdx
I4=0LϕxϕxdxI5=0Ld4ϕdx4ϕxdxI6=0Lxadxdxdx
I7=0LϕxdxdxI8=0Lxad3ϕdx3dxdxI9=0Ld2ϕdx2dxdx
I10=0LxaϕxϕxdxI11=0Ld2ϕdx2xaϕxdxI12=I7
I13=0Lϕ2xxa2dxI14=0Ld2ϕdx2ϕxxa2dxI15=0Ldxϕxxadx
I16=0Ld4ϕdx4ϕxxa2dxI17=0Ld3ϕdx3ϕxxadxI18=0Ld2ϕdx2ϕxxadx
I19=0Ld4ϕdx4ϕxxadx

Advertisement

AijExtensional stiffness
BijCoupling stiffness
DijBending stiffness
EYoung’s modulus
hShell thickness
LShell length
vPoisson’s ratio
ωNatural angular frequency
dijij=123Shell parameters
qAcoustic pressure
cSpeed of sound in the fluid
JnkrrBessel’s function of first kind
cfSpeeds of the sound in fluid – filled cylindrical shell
ΩNon-dimensional frequency parameter
E1, E2Young’s moduli
ν1, ν2Poisson’s ratios.
ρ1, ρ2Mass densities
Efgm, vfgm, ρfgmEffective material quantities
TKelvin scale
CrMaterial properties
VrVolume fractions
kNumber of ring supports
R0, R1Outer and inner radii of the shell
h/RThickness to radius ratio
L/RLength to radius ratio
nCircumferential wave number
mHalf- axial wave number
CMaterial property
C0, C−1, C1, C2, C3Coefficients of temperature
pVolume fraction law or Power law exponent
QijReduced stiffness
RShell radius
T(K)Temperature function of Kelvin
xAxial coordinate
ΨCircumferential coordinate
zRadial coordinate or thickness variable
u(x, Ψ, t)Deformation displacement functions in x direction
v(x, Ψ, t)Deformation displacement functions in Ψ direction
w(x, Ψ, t)Deformation displacement functions in z direction
VfiVolume fraction
ρMass density of shell material
ρtMass density per unit length
S11,S22,S12Surface curvatures
krRadial wave number
kmAxial wave number
cLSpeeds of the sound in the empty cylindrical shell
ρfIndicates the density of the fluid density
pm, qm and rmVibration amplitudes in the x, Ψ and z directions
CSsCylindrical shells
FGMFunctionally graded material
WFWith fluid
FLTFluid loading term
C-ICategory-I
WPAWave propagation approach
ODEsOrdinary differential equations
PDEsPartial differential equations
FFCSFluid-filled cylindrical shell
WOFWithout fluid
CWNCircumferential wave number
C-IICategory-II
FEMFinite element method

References

  1. 1. Greif R, Chung H. Vibrations of constrained cylindrical shells. AIAA Journal. 1975;13(9):1190
  2. 2. Loy CT, Lam KY, Shu C. Analysis of cylindrical shells using generalized differential quadrature method. Shock and Vibrations. 1997;4:193-198
  3. 3. Loy CT, Lam KY, Reddy JN. Vibration of functionally graded cylindrical shells. International Journal of Mechanical Sciences. 1999;41:309-324
  4. 4. Zhang XM, Liu GR, Lam KY. Vibration analysis of thin cylindrical shells using wave propagation approach. Journal of Sound and Vibration. 2001a;239(3):397-403
  5. 5. Najafizadeh MM, Isvandzibaei MR. Vibration of functionally graded cylindrical shells based on higher order shear deformation plate theory with ring support. Acta Mechanica. 2007;191(1–2):75-91
  6. 6. Rayliegh JWS. Theory of Sound. London: Macmillan; 1884
  7. 7. Love AEH. On the small free vibrations and deformation of thin elastic shell. Philosophical Transactions. Royal Society of London. 1888;A179:491-549
  8. 8. Arnold RN, Warburton GB. The flexural vibrations of thin cylinders. Proceedings of the Institution of Mechanical Engineers A. 1953;167:62-80
  9. 9. Chung H, Turula P, Mulcahy TM, Jendrzejczyk JA. Analysis of cylindrical shell vibrating in a cylindrical fluid region. Nuclear Engineering and Design. 1981;63(1):109-1012
  10. 10. Goncalves PB, Batista RC. Frequency response of cylindrical shells partially submerged or filled with liquid. Journal of Sound and Vibration. 1987;113(1):59-70
  11. 11. Goncalves PB, da Silva Frederico MA, del Prado Zenon JGN. Transient stability of empty and fluid-filled cylindrical shells. J. of the Braz. Soc. of Mech. Sci. & Eng. 2006;XXVIII(3/331):1-8
  12. 12. Gasser LFF. Free vibrations on thin cylindrical shells containing liquid [M.S. thesis]. Rio de Janeiro: Federal Univ. of Rio de Janeiro, peccoppe-ufrj; 1987 (in Portuguese)
  13. 13. Sharma CB, Johns DJ. Vibrations characteristics of clamped-free and clamped-ring-stiffened circular cylindrical shells. Journal of Sound and Vibration. 1971;14:459-474
  14. 14. Xi ZC, Yam LH, Leung TP. Free vibration of a partially fluid-filled cross-ply laminated composite circular cylindrical shells. The Journal of the Acoustical Society of America. 1997;101(2):909-917
  15. 15. Zhang XM, Liu GR, Lam KY. Coupled vibration of fluid-filled cylindrical shells using the wave propagation approach. Applied Acoustics. 2001b;62:229-243
  16. 16. Xiang Y, Ma YF, Kitipornchal S, Lim CW, Lau CWH. Exact solutions for vibration of cylindrical shells with intermediate ring supports. International Journal of Mechanical Sciences. 2002;44:1907-1924
  17. 17. Zhao X, Liew KM, Ng TY. Vibrations of rotating cross-ply laminated circular cylindrical shells with strings and ring stiffeners. International Journal of Solids and Structures. 2002;39:529-545
  18. 18. Xiang Y, Wang CM, Lim CW, Kitipornchai S. Buckling of intermediate ring supported cylindrical shells under axial compression. Thin-Walled Structures. 2005;43(3):427-443
  19. 19. Isvandzibaei MR, Awasare PJ. Effects edge and free-free boundary conditions for analysis, free vibration of functionally graded cylindrical shell with ring based on third order shear deformation theory using Hamilton’s principle. World Academy of Science, Engineering and Technology. 2010;61:219-225
  20. 20. Lee HL, Chang WJ. Vibration analysis of fluid-conveying double-walled carbon nanotubes based on nonlocal elastic theory. Journal of Physics: Condensed Matter. 2009;21:115302, 1-5
  21. 21. Silva FMA, Goncalves PB, Del Prado ZJGN. Nonlinear vibrations of axially loaded cylindrical shells partially filled with fluid. In: Proceedings of DINCON’10 9th Brazilian Conference on Dynamics, Control and their Applications, June 07–11; 2010. pp. 643-663
  22. 22. Shah AG, Mahmood T, Naeem MN, Arshad SH. Characteristics of fluid-filled cylindrical shells based on elastic foundations. Acta Mechanica. 2011;216:17-28
  23. 23. Xiong L, Li TY, Zhu X, Zhu XM. Free vibration analysis of fluid-filled elliptical cylindrical shells. Applied Mechanics and Materials. 2013;437:102-109
  24. 24. Hasheminejad SM, Alaei-Varnosfadrani M. Active radiation control of a submerged piezocomposite hollow sphere. Journal of Intelligent Material System and Structure. 2015;26:2073-2091
  25. 25. Kim Y. Effect of partial elastic foundation on free vibration of fluid-filled functionally graded cylindrical shells. Acta Mechanica Sinica. 2015;31:920
  26. 26. Izyan N, Viswaanathan MD, Aziz KK, Prabaker K. Free vibration of layered cylindrical shell filled with fluid. Applied Mathematics and Mechanics (English Edition). 2016;37:803
  27. 27. Soutis C, Shahid HA, Naeem MN. Influence of ring support on free vibration of sandwich functionally graded cylindrical shell with middle layer of isotropic material. Journal of Engineering Research. 2016;4(1):159-186
  28. 28. Amabili M, Pellicano F, Paidoussis MP. Nonlinear vibrations of simply supported, circular cylindrical shells, coupled to quiescent fluid. Journal of Fluids and Structure. 1998;12(7):883-918
  29. 29. Goncalves PB, Batista. Non-linear vibration analysis of fluid-filled cylindrical shells. Journal Sound and Vibration. 1988;127(1):133-143
  30. 30. Ansari R, Rahimi GH, Hemmatnezhada M. Vibration of functionally graded cylindrical shells with ring support. Scientia Iranica B. 2011;18(6):1313-1320
  31. 31. Gonacalves PB, Del Prado ZJGN. Low-dimensional Galerkin models for nonlinear and instability analysis of cylindrical shell. Nonlinear Dynamics. 2005;41(1):129-145
  32. 32. Weiqiu C, Ding HJ, Guo YM, Yang QD. Free vibrations of fluid-filled orthotropic cylin-25 drical shells. Journal of Engineering Mechanics. 1997;123(11):1130
  33. 33. Shah SA. Axially symmetric vibrations of fluid-filled poroelastic circular cylindrical shells. Journal of Sound and Vibration. 2008;318:389-405
  34. 34. Pellicano F, Amabili M. Stability and vibration of empty and fluid-filled cylindrical shells under static and period axial loads. International Journal of Solids and Structures. 2003;40(13?14):3229-3251
  35. 35. Warburton GB. Vibration of thin cylindrical shells. Journal of Mechanical Engineering Science. 1965;7:399-407
  36. 36. Zhang L, Xiang Y. Vibration of open circular cylindrical shells with intermediate ring supports. International Journal of Solids and Structures. 2006;43:3705-3722
  37. 37. Liu Y, Chu F. Nonlinear vibrations of rotating thin circular cylindrical shell. Nonlinear Dynamics. 2012;67:1467-1479
  38. 38. Toulokian YS. Thermophysical Properties of High Temperature Solid Materials. New York: Macmillan; 1967

Written By

Muzamal Hussain, Aamir Shahzad, Muhammad Nawaz Naeem and Maogang He

Submitted: 07 June 2017 Reviewed: 02 November 2017 Published: 14 February 2018