Research Article  Open Access
Fuzhen Pang, Haichao Li, Kwangnam Choe, Dongyan Shi, Kwanghun Kim, "Free and Forced Vibration Analysis of Airtight Cylindrical Vessels with Doubly Curved Shells of Revolution by Using JacobiRitz Method", Shock and Vibration, vol. 2017, Article ID 4538540, 20 pages, 2017. https://doi.org/10.1155/2017/4538540
Free and Forced Vibration Analysis of Airtight Cylindrical Vessels with Doubly Curved Shells of Revolution by Using JacobiRitz Method
Abstract
This paper presents free and forced vibration analysis of airtight cylindrical vessels consisting of elliptical, paraboloidal, and cylindrical shells by using JacobiRitz Method. In this research, the theoretical model for vibration analysis is formulated by FlÃ¼ggeâ€™s thin shell theory and the solution is obtained by RayleighRitz method. The vessel structure is divided into shell components (i.e., ellipsoid, parabolic, and cylinder) and their segments, and each displacement field of shell segments is represented by the Jacobi polynomials and the standard Fourier series. The continuous conditions at the interface are modeled by using the spring stiffness technique. The reliability and the accuracy of the present method are verified by comparing the results of the proposed method with the results of the previous literature and the finite element method (FEM). Moreover, some numerical results for free and forced vibration of ellipticalcylindricalelliptical vessel (ECE vessel) and paraboloidalcylindricalelliptical vessel (PCE vessel) are reported.
1. Introduction
The coupled shell consisting of a cylindrical shell and doubly curved shells of revolution are widely used in aerospace and marine structures and civil and mechanical engineering. For example, hermetic capsules, pressure vessels, and storage vessels can be modeled as ellipticalcylindricalelliptical shell structures; rockets and submarines can be viewed as paraboliccylindricalelliptical shell structures. So it is important to analyze the vibration characteristics of the coupled shell mentioned above.
To emphasize the characteristic and purpose of this paper, the existing research related to this paper is reviewed and shown below: Leissa [1] and Qatu [2] analyzed the classical thin shell theories (Loveâ€™s, Reissnerâ€™s, Naghdiâ€™s, Sandersâ€™, and FlÃ¼ggeâ€™s, etc.) and the thick shell theories. Based on these theories, many researchers have proposed free and forced vibration analysis methods for some structures. Liew and Lim [3] proposed a pb2 Ritz method to study the free vibration of the doubly curved shallow shells in form of a rectangular plane. Tan [4] presented an efficient substructuring analysis method to investigate the free vibration characteristics for doubly curved shells of revolution based on the firstorder shear deformation shell theory and the classical thin shell theory. Tang et al. [5] researched the free and forced vibration analysis of multistepped circular cylindrical shell with arbitrary boundary conditions with the employment of the method of reverberationray matrix. Xie et al. [6] presented a semianalytic method to analyze free and forced vibrations of combined conicalcylindricalspherical shells with ring stiffeners and bulkheads. Xie et al. [7] presented an analytic method for free and forced vibration analysis of stepped conical shells with general boundary conditions. Qatu and Asadi [8] used the Ritz method to present vibration frequencies for spherical, cylindrical, and hyperbolic paraboloidal shells with 21 possible boundary conditions by using thin shallow shell theory. Zhang et al. [9] proposed an analytical method to study the free and forced vibration behaviors of a submerged finite elliptic cylindrical shell. Jiang et al. [10] derived a unified solution for the vibration analysis of doubly curved shallow shells with arbitrary elastic supports by using an improved trigonometric series and the RayleighRitz method. Messina [11] studied the free vibrations of multilayered doubly curved shells based on a mixed variational approach and global piecewisesmooth functions. Poultangari and NikkhahBahrami [12] developed the Vectorialwave method (VWM) to study free and forced vibrations of cylindrical shells in the presence of dampers at supports. Chen et al. [13] presented an analytic method to analyze free and forced vibration characteristics of ringstiffened combined conicalcylindrical shells with arbitrary boundary conditions. Tornabene et al. [14â€“16] studied the dynamic behavior of functionally graded or laminated composite doubly curved shells and panels of revolution by using the Generalized Differential Quadrature (GDQ) method. Chen et al. [17] presented wave based method to analyze the free and forced vibration of cylindrical shells with discontinuity in thickness. Ma et al. [18] presented a free and forced vibration analysis of coupled conicalcylindrical shells with arbitrary boundary conditions using a modified FourierRitz method. Dozio [19] deals with the formulation of advanced twodimensional Ritzbased models for accurate prediction of natural frequencies of thin and thick sandwich plates core made of functionally graded material. In [20, 21], research results are presented on dynamic behavior of the composite laminated or the functionally graded material (FGM) doubly curved shells and panels of revolution with arbitrary boundary conditions. In literatures the displacement fields are represented by using several auxiliary functions and a standard Fourier series, the classical boundary conditions and the general elastic restraint are imitated by the coupling spring stiffness technique, and a unified numerical analysis model is presented to solve the free vibration of doubly curved shells and panels by using the RayleighRitz method. Dozio [22] used the statespace levyâ€™s method to solve the free vibration problem of specially orthotropic multilayered cylindrical and spherical panels. Dozio and Carrera [23] described a variablekinematic Ritz formulation to accurately predict free vibration of thick and thin, rectangular, and skew multilayered plates with clamped, free, and simply supported boundary conditions.
Galletly and Mistry [24] investigated the vibration characteristics of the closed cylindrical shells with cones, hemispheres, ellipsoids, and so on by using the finite difference method. Irie et al. [25] represented the vibration equation of the shells as a combined set of firstorder differential equations by using the transfer matrix method, and the research studied vibration characteristics of the coupled conicalcylindrical and annular plateâ€“cylindrical shell. Qu et al. [26] presented a domain decomposition technique for solving vibration problems of uniform and stepped cylindrical shells with arbitrary boundary conditions. Wong and Sze [27] investigated low frequency vibration of the torsion free axisymmetric modes of thin cylindrical shell with hemispherical caps vibrating in vacuo by using membrane approximation. Rougui et al. [28] investigated the nonlinear free and forced vibrations of simply supported thin circular cylindrical shells using Lagrangeâ€™s equations and an improved transverse displacement expansion. Efraim and Eisenberger [29] studied the dynamic behaviors of coupled shells made of conical, cylindrical, and plate segments by using a dynamic stiffness matrix method. Caresta and Kessissoglou [30] presented the vibration characteristics of isotropic coupled conicalcylindrical shells by using a wave solution and a power series expansion method. Lee [31, 32] presented the axisymmetric free vibration analysis of sphere, hemisphere and the joined sphericalcylindrical shell structures by using a pseudospectral method. Qu et al. [33â€“36] proposed a kind of domain decomposition method to study the vibration characteristics of various combination shell structures consisting of conical, cylindrical, spherical shell, and ring stiffeners. Chen et al. [13] presented an analytic method for free and forced vibration analysis of ringstiffened combined conicalcylindrical shells with arbitrary boundary conditions by using FlÃ¼ggeâ€™s theory and the smeared out method.
As discussed above, many of the vibration analyses focus on the simple doubly curved shells, while the study on the vibration analysis for their coupled shell structures is relatively small. Furthermore, most of the mentioned researches studied the vibration problems of the jointed shell structures consisting of cylindrical, conical, and spherical shells. And the authors prefer the free vibration analysis instead of the forced vibration analysis of the structures, especially for the structures of doubly curved shells of revolution. However, the results of these studies provide the feasibility for solving the free and forced vibration problems of airtight cylindrical vessels with doubly curved shells of revolution.
The main purpose of this work is to present a generalized and unified JacobiRitz formulation to investigate the free and forced vibration characteristics of airtight cylindrical vessels consisting of the elliptical, paraboloidal, and cylindrical shells. This paper is characterized by not only simplifying the calculation by applying the RayleighRitz method and the penalty function method, but also generalizing the selection of the permissible displacement function by using the Jacobi polynomial. FlÃ¼ggeâ€™s thin shell theory is employed to do vibration analysis for theoretical model; the continuous conditions at the interface are modeled by using the spring stiffness technique. Based on the domain decomposition method [33â€“36], the vessel structures are divided into ellipsoid, parabola, cylinder, and their segments, and displacement functions are represented by the Jacobi polynomials and the standard Fourier series. The comparisons with the results of the previous paper or the finite element method (FEM) are presented, and some numerical results for free and forced vibration of the ellipticalcylindricalelliptical shell and the paraboloidalcylindricalelliptical shell are reported.
2. Theoretical Formulations
2.1. Geometry of the Airtight Cylindrical Vessel
The geometric model for the airtight cylindrical vessel with doubly curved shell of revolution is shown in Figure 1.
The airtight cylindrical vessel is divided into its components, that is, the elliptical, cylindrical, and paraboloidal shell, depending on the position of the junctions. The orthogonal coordinate systems () for the left and right doubly curved shells are fixed in the middle surface of shells, where are the coordinates along the meridional direction, are the circumferential coordinates, and are normal coordinates. Also the cylindrical shell is described by a cylindrical coordinate system (), where the subscripts , , and represent the left, cylindrical, and right shell, respectively. We need to remember that there are two main curvatures radii (curvature radius of the meridional direction , curvature radius of the circumferential direction ) in geometry of the doubly curved shells [20, 21], and two curvatures radii of the elliptical shell and the paraboloidal shell can be expressed as follows:
Elliptical shell:where and are the lengths of the semimajor and semiminor axes of the elliptic meridian, which are given by and , respectively.
Paraboloidal shell:where is the characteristic parameter of the parabolic meridian.
Specially,where is the coordinate angle along the meridian at the junction between the paraboloidal shell and the cylindrical shell. Also, in this study, the spherical shell is considered to be a special case of the elliptical shell; the extents of the coordinate angle along the meridian for elliptical shell and the paraboloidal shell are, respectively, defined by and . Regardless of each shell component, it is assumed that all shells are composed of homogeneous and isotropic materials; the elastic modulus , mass density , Poissonâ€™s ratio , and thickness have the same values, respectively. The displacement of each shell component is presented by , , and () at coordinate systems. In order to increase the accuracy of calculations, these shell components may be further decomposed into left, cylindrical, and right shell component segments along the meridian or the axial of revolution (see Qu et al. [34]).
2.2. Energy Functions of the Airtight Cylindrical Vessel
In this work, FlÃ¼ggeâ€™s thin shell theory [1] is employed to define the relationships between strain and displacement of the middle surface of the th segment in the shell component. They are written as follows:where , , and are strain components of the meridional, circumferential, and shearing direction in the middle surface of the shell, , , and are the curvature changes of the middle surface in the shell, and , , and are the displacement components in the meridional, circumferential, and normal directions of the shell, respectively. In above equations, the superscript represents the th segment in the shell component and the subscript represents the left, cylindrical, and right shell, respectively. Also, and are the LamÃ© parameters; and indicate the curvature radius in the meridional directions and circumferential directions of middle surface, respectively. The coordinate system and geometric parameters of each shell components in Figure 1 are defined as follows:
(1) Left shell:
(2) Cylindrical shell:
(3) Right shell:
With reference to the FlÃ¼ggeâ€™s thin shell theory, the normal and shear strains of the shell components are expressed as follows:where .
From Hookeâ€™s law on the relationship between stress and strain, the corresponding stresses are expressed as follows:where , , and are the normal and shear stresses and are the constants indicating the relationship between stress and strain, which are defined as follows:where and are Youngâ€™s modulus and Poissonâ€™s ratio, respectively.
From the theory of elasticity, the strain energy equation in the shells can be expressed as
In the case of the thin shells, we can simplify some of the terms in the energy equation, and the energy equation is rearranged through a mathematical process:where
And the kinetic energy of the shell components can be defined as
In order to investigate the forced vibration of the airtight cylindrical vessels under point loads, it is assumed that the external unit force acts on the middle surface of the vessel in the meridional, circumferential, and the normal directions, respectively. The virtual work by the external forces can be written aswhere is the Dirac function and , , and are the external force in the meridional, circumferential, and the normal directions, respectively.
2.3. Continuity Conditions of the Airtight Cylindrical Vessel
In this analysis, the continuity conditions of the airtight cylindrical vessel structures with doubly curved revolution shells are modeled by the penalty parameters, which are characterized by artificial spring stiffness. The penalty parameter permits the flexible selection of the admissible displacement functions of the airtight cylindrical vessels, and the appropriate value of the penalty parameter ensures fast convergence of the accurate solution. Thus, the potential energy stored in the connective springs at the junction between the doubly curved shells and the cylindrical shell can be described aswhere , , , and denote the stiffness of the translational and rotational springs between the shell components, respectively.
If each shell component is further subdivided into shell segments, the potential energy stored in the connective springs of two adjacent shell segments can be described as follows:where , , and represent the potential energy stored in the connective springs of the shell segments at the left, cylindrical, and right shell, respectively. And the superscripts, and , represent the th and th shell segments.
Therefore, the total potential energy reflecting the connective conditions can be expressed aswhere , , and are the number of segments divided in the each shell components.
2.4. Unified Solution and Solution Procedure
The selection of an admissible displacement function is an important factor in guaranteeing a stable convergence and accuracy of the solution. In the handling of the continuity conditions, the use of spring stiffness, which can be seen as penalty parameters, makes the selection of admissible functions very flexible [37â€“47]. Subscript is omitted here for the sake of brevity. In this work, the admissible displacement functions of shell segments are generally extended to the Jacobi orthogonal polynomials and Fourier series, regardless of shape of the shell components. As we know, classical Jacobi polynomials [48] are defined on the interval of and the recurrence formula of the classical Jacobi polynomials of degree is given bywhere and .
The orthogonality condition of classical Jacobi polynomials iswhere and
The Jacoby polynomials are generalized orthonormal polynomials containing some orthonormal polynomials such as Legendre, Chebyshev, and Gegenbauer polynomials. For example, the choice yields the Chebyshev polynomials of the first kind, while choosing gives the Chebyshev polynomials of the second kind. Also the choice yields the Legendre polynomials, while choosing gives the Gegenbauer polynomials. Thus, the selection of an admissible displacement function for the airtight cylindrical vessels is more generalized by the Jacobi polynomials.
The displacement functions of the shell segments can be written in the formswhere , , and are unknown coefficients of the Jacobi polynomials that you want to obtain; is the th order Jacobi polynomial for the displacement function in the meridional direction; is an angular frequency; denotes time. The nonnegative integers and represent the highest degrees taken in the Jacobi polynomials and the Fourier series; the positive integer represents the mode circumference wave number. The subscript means the left, cylindrical, and right shells, respectively. The classical Jacobi polynomials are complete and orthogonal polynomials defined on the interval of . Thus, a linear transformation rule is introduced for coordinate conversion from interval of the th shell segment to interval of the Jacobi polynomials, that is, .
The total Lagrangian energy functions () of the airtight cylindrical vessels with doubly curved shells of revolution can be written in the formswhere represent total potential energy reflecting the connective conditions and represent the virtual work by the external forces.
The total Lagrangian energy function is minimized with respect to the unknown coefficients by using the RayleighRitz method.where K and M are the stiffness and mass matrix of the airtight cylindrical vessel, respectively. And X and F represent the coefficient vector and external force vector. If F is a zero vector, the above equation is the free vibration analysis equation, and if it is not a zero vector, it is the forced vibration analysis equation. A detailed description of (21) is given in Appendix.
3. Numerical Results and Discussions
3.1. Convergence and Comparison Studies
In this section, the reliability and accuracy of this method are verified by comparing the results of the method with those of the previous literature or the finite element method (FEM). From the mathematical point of view, unlike the domain decomposition method which uses the Lagrange multipliers and the weight residual least squares method [35], when the penalty parameters is defined as a very high value the solution of the proposed method may not converge [38â€“40]. And it is necessary to emphasize that the research will focus on the elastic deformation rather than rigid deformation. So some zero eigenvalue obtained in the solution will be ignored in the following research. The convergence characteristics of the nondimensionless frequency parameters for different spring stiffness values in the ECE and the PCE vessels are shown in Figure 2. Where the material parameters and the geometric dimensions of the ellipticalcylindricalelliptical and the paraboloidalcylindricalelliptical vessels are as follows:
(a)
(b)
ECE vessel:â€‰.
PCE vessel:â€‰.
In this analysis, the number of segments of the two structures is taken as , , and , the parameters of the Jacobi polynomials are taken as , and the nondimensionless frequency parameter is defined as .
As can be seen in Figure 2, it can be seen that the connection spring stiffness value ensures stable convergence within the range of = 10^{0}E~10^{6}E. Therefore, unless otherwise noted, in this paper, it is assumed that the connective spring stiffness is the same (i.e., ).
Next, in order to verify the accuracy of the presented method, free vibration of a sphericalcylindricalspherical shell having the same material characteristics and structural dimensions as those of the previous literature is considered. Frequency convergence results for the shell segments () of the sphericalcylindricalspherical shell are shown in Table 1 and the results are compared with those of the previous literature [37].

In Table 1, it can be seen that the frequency quickly converges to correct solution with the number of shell segments increasing. However, for simplicity the calculation, the number of shell segments is limited to eight or less in this paper. The percentage error of the solution for the Jacobi parameters and in the ellipticalcylindricalelliptical shell are pointed out in Figure 3.
(a)
(b)
(c)
As shown in Figure 3, it can be seen that the change of the Jacobi parameters and does not have a great influence on the frequency convergence and the maximum error does not exceed . This means that not only the Chebyshev, Legendre, Bernstein, and Gegenbauer polynomials but also the generalized Jacobi polynomial can be used. In Table 2, the comparison between the results of the ECE and the PCE vessels pointed out above and the finite element method (FEM) are shown. In here, the ECE vessel used the S4R model and the 19433 element, and the PCE vessel used the S4R model and the 19389 element of the ABAQUS.

The validity and accuracy of the presented method were verified by comparisons with the above results. For the sake of understanding, unless otherwise stated, the parameters of the structure and the material properties in the subsequent study are taken as follows:â€‰.
Nondimensional frequency parameter is defined as .
3.2. Free Vibration of the Airtight Cylindrical Vessel
In this section, the study on free vibration analysis of ellipticalcylindricalelliptical vessel and paraboloidalcylindricalelliptical vessel is performed. The change in frequency parameters for the axial wave number and the circumferential wave number in the ECE vessel and the PCE vessel used in Figure 2 is shown in Figure 4.
(a)
(b)
As shown in Figure 4, when , the lowest frequency of the two structures is shown, and, as increases, the frequency parameters increased and the difference between the frequency parameters is decreased. Table 3 shows the frequency parameters according to the ratio between the semimajor axis and the semiminor axis of the right ellipsoid in the ellipticalcylindricalelliptical vessel, and Table 4 shows the frequency parameters according to the ratio between the maximum radius and the length of the paraboloid in the paraboloidalcylindricalelliptical vessel.


From Tables 3 and 4, it can be seen that the frequency parameters of the structures are significantly reduced as two ratios gradually increase. In addition, to facilitate the readerâ€™s understanding, some mode shapes of the ECE vessel and PCE vessel are shown in Figure 5.
(a) ECE vessel
(b) PCE vessel
3.3. Forced Vibration of the Airtight Cylindrical Vessels
Analyzing the forced vibration of a vessel subjected to a point load is an important factor in studying the vibration response characteristics of the vessel. In this study, the forced vibration of the structure is investigated assuming that the axial unit force and the radial unit force act separately on the vessels used in Figure 1. And the structural damping is realized by using a complex modulus instead of the original Youngâ€™s modulus , where is defined as the damping of the vessel structures [33]. Forced vibration analysis is investigated between the frequency domain 20~500â€‰Hz and the calculation frequency step is adopted as 0.2â€‰Hz. The fourpoint coordinate systems D_{1}~D_{4} are set in the ECE vessel and the PCE vessel. In the ECE vessel and the PCE vessel, the radial displacements of point when the radial force acts on point are shown in Figure 5, and this is compared with FEM results (ECE vessel; S4R model, 19433 element, PCE vessel; S4R model, 19389 element).
The comparisons of Figures 5 and 6 show that the results of the presented method are in good agreement with those of ABAQUS. Figure 7 shows the radial displacements of the points , , and when the axial force acts on the point of the ECE and PEC vessels.
(a)
(b)
(a)
(b)
The forced vibrations of the airtight vessels show almost the same resonance peak at the three test points of the airtight vessel. However, as shown in the Figure 7, when the axial force is acted, the displacement is the largest at point , which is far from the point of action of the force, and point is the next. This means that the influence of vibrations at the junction of the shell component is great when the airtight vessel is subjected to external action.
4. Conclusions
In this paper, the unified JacobiRitz Method was presented and the free and forced vibrations of the airtight vessels consisting of elliptical, cylindrical, and paraboloidal shells were analyzed. In this study, FlÃ¼ggeâ€™s thin shell theory was adopted for the calculation model of vessels; the selection of displacement functions was generalized by the Jacobi polynomials. In addition, the continuous conditions at the interface were incorporated into the presented model by using the spring stiffness technique, and the accuracy of solution was improved by the domain decomposition method. The results of the proposed method show that this method can vouch faster convergence, high accuracy, and simplicity of calculation. Furthermore, the results of various free and forced vibrations of the ellipticalcylindricalelliptical and the paraboloidalcylindricalelliptical vessels were presented. We hope that these results can be used as benchmark data for future researchers.
Appendix
The generalized mass and stiffness matrix of the airtight cylindrical vessels used in (21) is given aswhere the subscript represents the left, cylindrical, and right shells and the superscript represents the th segment in the shell components.
Left and right shell:where .where :
Cylindrical shells :