Research Article | Open Access

Volume 2019 |Article ID 4157930 | https://doi.org/10.1155/2019/4157930

Xue-Qin Li, Wei Zhang, Xiao-Dong Yang, Lu-Kai Song, "A Unified Approach of Free Vibration Analysis for Stiffened Cylindrical Shell with General Boundary Conditions", Mathematical Problems in Engineering, vol. 2019, Article ID 4157930, 14 pages, 2019. https://doi.org/10.1155/2019/4157930

# A Unified Approach of Free Vibration Analysis for Stiffened Cylindrical Shell with General Boundary Conditions

Revised16 Jun 2019
Accepted17 Jun 2019
Published10 Jul 2019

#### Abstract

A unified approach of free vibration analysis for stiffened cylindrical shell with general boundary conditions is presented in this paper. The vibration of stiffened cylindrical shell is modeled mathematically involving the first-order shear deformation shell theory. The improved Fourier series is selected as the admissible displacement function while the arbitrary boundary conditions are simulated by adjusting the equivalent spring stiffness. The natural frequencies and modal shapes of the stiffened shell are obtained by solving the dynamic model with the Rayleigh-Ritz procedure. Various numerical results of free vibration analysis for stiffened cylindrical shell are obtained, including natural frequencies and modes under simply supported, free, and clamped boundary conditions. Moreover, the effects of stiffener on natural frequencies are discussed. Compared with several state-of-the-art methods, the feasibility and validity of the proposed method are verified.

#### 1. Introduction

As one critical load-bearing component in structural engineering, stiffened cylindrical shells are widely used in various industrial equipment such as aerospace, rocket, and infrastructure [1ā3]. On the premise of guaranteeing the durability and service life of whole structures, the stiffened structure is conducive to save the consumption of materials used, reduce weight, and improve the mechanical properties of the whole structures. However, compared with simple cylindrical shells, the application of stiffeners also brings in high-complexity in describing its dynamic behaviors and performing its performance analysis. Therefore, it is urgently desired to develop an efficient approach to analyze the dynamic characteristics of stiffened cylindrical shells. As an important analysis technique, free vibration analysis is frequently required and has always been one study focus since it can provide valuable insight into dynamic behaviors and vibration control of stiffened cylindrical shells [4ā7]. In past several decades, the increasing application of stiffened cylindrical shells has motivated great efforts in developing more accurate mathematical models and corresponding approaches for analyzing their dynamic behaviors. In this case, a large variety of modern theories and advanced approaches for vibration analysis have been achieved and the dynamic response characteristics of stiffened cylindrical shells have been well studied [8ā10]. For instance, Al-Najafi [8] chose each ring rib as a discrete element and utilized the finite element theory to study the free vibration of stiffened cylindrical shells with simply supported boundary at both ends. Jafari [9] selected power series as Ritz functions and investigated the effects of initial circumferential stress, centrifugal force, and rotation on the natural frequencies of ring-stiffened cylindrical shells. Chen [10] solved theoretical model and achieves vibration analysis of the ring-stiffened cylindrical shell by using wave method. The free vibration characteristics of ring-stiffened cylindrical shells under several classical boundary conditions were obtained, in which ring-stiffened and frame-stiffened plates were used. However, it appears that most of the previous studies on stiffened cylindrical shells are confined to classic boundary conditions. In fact, as we know, a variety of possible boundary conditions encountered in engineering may not always be classical in nature [11]. Besides, the existing solution procedures are often customized for a specific set of different boundary conditions; the constant model modifications would lead to massive solution procedure operations, which result in very tedious calculations because even only the classic cases still need a total of hundreds of different combinations [12, 13]. Therefore, to simplify the calculation tasks and obtain the accurate vibration characteristics, a unified method is urgently desired to tackle with the general boundary conditions problems in free vibration analysis of stiffened cylindrical shells. In this paper, to address general boundary conditions issues and build unified free vibration analysis method, two obstacles should be overcome: establishment of a high-accuracy dynamic model for stiffened cylindrical shells; development of an efficient solution procedure to resolve the aforementioned dynamic model.

To achieve unified vibration analysis and obtain accurate solution results, as we know, it should be the first to construct high-accuracy mathematical models. The current shell deformation theories on model formulation are mainly classified into three categories: Classical Shell Deformation Theory (CSDT), First-order Shear Deformation Theory (FSDT), and Higher-order Shear Deformation Theory (HSDT). In view of the Kirchhoff-Love assumptions, CSDT neglects neglect of transverse shear strains and possesses the advantages of simple calculation tasks [14]. To date, various subcategory CSDTs were derived and employed in vibration analysis of thin cylindrical shells [15, 16]. Although CSDT can achieve accurate vibration results for thin shells, it is unsuitable to describe the vibration characteristics of the stiffened cylindrical shell with larger thickness value. As one important alternative, the managerial theory-based FSDT were developed, which abandons the hypothesis in CSDT that the transverse normal is still perpendicular to the mid-surface after structural deformation, thus overcoming the defect that the neglect of transverse shear deformation in CSDT [17]. In light of the increasing options of admissible displacement function, the FSDT have been widely used in modeling and analysis of medium thick plates, laminated plates, and moderately thick shells [18ā20]. Moreover, to further enhance the analysis precision, HSDT were developed, in which the transverse normal is no longer perpendicular to the middle after deformation [21, 22]. Although the HSDTs are capable of solving the global dynamic problem of shells more accurately, it also brings sophisticated formulations and boundary terms, which extremely increases the computational complexity and tasks [23, 24]. Therefore, HSDT is deemed to be unpractical in vibration analysis of complex shell structures like stiffened cylindrical shells. In the present work, the FSDT with high computing accuracy and little computing burden was chosen to formulate the theoretical model of stiffened cylindrical shells.

Apart from the aforementioned models with shear deformation theories, to accomplish high-accuracy and high-efficiency unified vibration analysis under general boundary conditions, the corresponding efficient and accurate solution approaches also should be developed to determine the vibration characteristics of stiffened cylindrical shells. A number of computational methods are available for vibration analysis of cylindrical shells, such as polynomial Ritz method (PRM) [25, 26], wave propagation method (WPM) [27ā29], transfer matrix method (TMM) [30, 31], Galletly Ritz method (GRM) [32, 33], Galerkin method [34ā36], finite element method (FEM) [37ā41], the continuous element method (CEM) [42], and Rayleigh-Ritz method [43ā45]. By less computational effort to achieve more accurate results and smooth derivative terms, Rayleigh-Ritz method possesses the potential of alleviating computational burdens and improving solution accuracy in the unified vibration analysis [46]. Nevertheless, although lots of virtues the Rayleigh-Ritz procedure held, it is still insufficient to deal with general boundary conditions owing to the used polynomial displacement function, which can hardly form a complete solution set or obtain stable results due to computer rounding errors [47, 48]. Unfortunately, to the best knowledge of the authors, the corresponding topic on displacement functions is very limited. Recently, Li [49] derived a complementary function and applied it to Fourier cosine series to express the dynamic behavior of beams, which solved the transverse vibration problem with elastic constraints. Zhang [50] proposed an improved Fourier cosine method in which the displacement function was assumed as the superposition of Fourier cosine series and supplementary functions. Dai [51] combined displacement functions of plateās transverse vibration and in-plane vibration to adapt to the vibration analysis of thin cylindrical shells and solved the problems of cylindrical shellās vibration by treating the both kinematic equations and boundary condition equations simultaneously. In this case, in view of the major benefits of only needs the kinetic energy and potential energy calculations rather than the specific stress variation calculations, the improved Fourier cosine series and its supplementary functions are deemed suitable as admissible displacement function in Rayleigh-Ritz procedure [52ā54].

The objective of this paper is to develop a unified vibration analysis method for stiffened cylindrical shells with general boundary conditions. The first shear deformation theory is employed to formulate the theoretical model, the improved Fourier cosine series and two supplementary functions are adopted as admissible displacement function, and the corresponding Rayleigh-Ritz procedure is performed to obtain the exact vibration behaviors of stiffened cylindrical shells. Without remodeling and corresponding procedures, the unified vibration analysis can be efficiently performed and the corresponding accurately vibration behaviors are obtained by only assigning the spring stiffness values of general boundary conditions. The convergence and accuracy of the developed method are validated by numerous examples in respect of frequency comparisons and modal comparisons.

The remainder of this paper is structured as follows. Based on the first-order shear deformation theory and the Rayleigh-Ritz procedure, the essential methodology of the unified vibration analysis method for stiffened cylindrical shells with general boundary conditions is investigated in Section 2. In Section 3, the proposed method is validated by comparing its free vibration analysis results with that of finite element method and several state-of-the-art methods. Some conclusions on this study are discussed and summarized in Section 4.

#### 2. Theoretical Formulations

##### 2.1. Description of Stiffened Cylindrical Shell

The schematic diagram of stiffened cylindrical shell is shown in Figure 1. The basic parameters of stiffened cylindrical shell involve length , radius , thickness , width of two adjacent stiffeners , and thickness of stiffeners . To analyze the vibration characteristics of stiffened cylindrical shell, the cylindrical coordinate system (, , ) is built on the left-edge middle surface; herein, represents the axis direction of the cylindrical shell; and indicate the circumferential and radial directions of cylindrical shell. Moreover, to simulate various boundary conditions, the equivalent linear springs , , and and torsion springs and are set at both edges of stiffened cylindrical shell.

##### 2.2. Displacement Field and Geometric Equation

To describe and quantify the vibration behaviors of moderately thick shell, we construct the FSDT dynamic model by considering transverse shear deformation. Assuming that the transverse normal would not remain vertical and mid-plane after shell deformation, the displacement field is established aswhere , , and indicate middle surface displacements along the -axis, -axis, and -axis, respectively; and indicate the rotation angle along the -axis and -axis, respectively; indicates the time variable.

In view of the von Karman geometric nonlinearity rule, the geometric equation is built aswhere

##### 2.3. Constitutive Relation and Internal Force Relation

In light of the general Hookās law, to further reveal the stress behavior in terms of displacement, the constitutive relation of stiffened cylindrical shell is indicated aswhere elastic constants is denotes as

Assuming that the stiffeners are of the same material with shell structures, the stress-strain relationship of stiffeners can be expressed aswhere is Youngās modulus of stiffener.

By integrating the in-plane stress in thickness and the stress in cross section, the relationship between force and moment, strain, and curvature on the middle plane is indicated aswhere , , and represent the normal force in -direction, -direction, and shear force resultants, respectively; , , and represent the bending moment in -direction, -direction, and twisting moment resultants, respectively; and represent the membrane stress and moment respect to stiffens, respectively; represents the shear correction factor used to correct the discrepancy between the constant stress state and real stress state, which is supposed as 5/6 in this paper. Moreover, the coefficients , , , , , and can be further expressed as

It should be noted that, for a symmetric cylindrical shell, =0, the bending and stretching vibration of the shell are uncoupled. Moreover, for the stiffeners distributed in shell, Bsā 0, the calculation burdens have increased dramatically. Thus, the general boundary conditions of the elastically restrained shell can be described as follows:āāAt left edge (=0):āāAt right edge (=L):

To this extent, all the needed parts of the first-order shear deformation theory (FSDT) are present, and they may be combined to obtain the desired form of energy expressions.

##### 2.4. Energy Expressions

In view of the efficiency and reliability of the results in modeling and solution procedure, the energy-oriented Rayleigh-Ritz method is employed in the present work. To define the energy expressions of the stiffened cylindrical shell, we assume that the spring retains a certain stiffness without considering mass; the total kinetic energy of the stiffened cylindrical shell is consisting of the kinetic energy of cylindrical shell and the kinetic energy of stiffener . In light of kinetic energy theory, the total kinetic energy of stiffened cylindrical shell can be described aswhere the and represent the density of cylindrical shell and stiffener, respectively. The inertia moment (= 0, 1, 2) can be further represented as

Similarly, by integrating the strain potential energy induced by deformation and the elastic potential energy induced by springs, the potential energy of stiffened cylindrical shell is expressed aswhere the strain potential energy is

in whichBesides, the elastic potential energy at the boundary is obtained as

##### 2.5. Admissible Displacement Functions and Deriving Solutions

To derive the discrete motion equation in the Rayleigh-Ritz procedure, the displacement and rotation components should be expanded with one admissible displacement function. For shell problem, the most commonly used form of admissible function is orthogonal polynomial and Fourier series. However, the low-order polynomials cannot form a complete solution set while the high-order polynomials tend to be numerically unstable due to computer rounding errors, which may lead to very tedious calculations and may be inundated with various boundary conditions. Although the Fourier series constitutes a complete solution set and holds good numerical stability by the Fourier series expansion, the convergence difficulties and low solution accuracy problems are still emerged in vibration analysis with general boundary condition. In this paper, to overcome the difficulty and satisfy the general boundary conditions, by using two supplementary functions, the displacement and rotation components of the middle surface of stiffened cylindrical shells are expanded into an improved form of Fourier cosine serieswhere and indicate the axial wave number and circumferential wave number; indicates the nondimensional parameter related to the axial wavenumber whose value is mĻ/L; , , , , and indicate the Fourier expansion coefficients; , , , , and the auxiliary function coefficients. The supplementary function can be indicated asThe uniformly convergent series expansions can be obtained by differentiating operation of series expansion. In fact, the solution holds arbitrary precision since the series expression has to be truncated numerically. In general, to obtain acceptably accurate solution, the infinite series expression is truncated to and in actual calculations.

Once the admissible displacement functions and energy functions of the stiffened cylindrical shell are established, the following tasks are to determine the coefficients in the admissible functions. In light of (12) and (14), the Lagrange energy function can be indicated in terms of kinetic energy and potential energy of the stiffened cylindrical shell as

In view of energy-oriented Rayleigh-Ritz procedure, the total expression of the Lagrangian energy function is minimized with respect to the undetermined coefficients

Assuming that the vector of unknown coefficients are , the corresponding coefficient equations can be retrieved aswhere

The natural frequencies and modes can be determined readily by solving the characteristic equation (22). Each column of the eigenvector matrix is the set of coefficients for the corresponding modes. By substituting the solution from (23) into (18), the corresponding mode shapes can be retrieved accordingly.

#### 3. Numerical Study

To verify the feasibility and validity of the presented approach, we compare the analysis results with several state-of-the-art methods from multiple perspectives, including the convergence analysis performed in Section 3.1, the frequency comparisons illustrated in Section 3.2, and modal comparisons depicted in Section 3.3. Besides, we summarize the detailed discussions in Section 3.4.

##### 3.1. Convergence Analysis

To verify the convergence of this method, the natural frequencies under elastic boundary conditions are calculated and compared with the finite element method. Axial direction is set as elastic and the stiffness coefficients are set as 105 while others are set as infinite large. The geometrical and material parameters of cylindrical shells are set as =39.45Ć10ā2 m, =4.9759Ć10ā2 m, =0.1651Ć10ā2 m, =2762 kg/m3, =0.3, =68.95 GPa, =0.3175Ć10ā2m, =0.5334Ć10ā2 m, =1.9725Ć10ā2 m, and =5. The calculation results are shown in Table 1; the data reveal good convergence and the results between present method and FEM fit well. It is clarified that all computations are performed on an Intel (R) Pentium (R) Desktop Computer (CPU G3260 3.30 GHz and 4.00 GB RAM). It is obvious that the present solution results hold excellent convergence and possesses high computing efficiency. For instance, the modes can be obtained within 0.5 s even the series are truncated as =15.

 Truncated Modes Time consuming, s n=1 n=2 n=3 n=4 n=5 n=6 n=7 n=8 M=5 1290.34 1661.19 2298.11 3612.72 3693.16 4455.88 4472.75 4756.00 0.26357 M=7 1261.44 1645.61 2252.24 3578.23 3587.75 4454.88 4460.87 4715.33 0.34305 M=9 1244.95 1635.48 2226.08 3533.22 3560.78 4453.71 4453.94 4692.21 0.37890 M=11 1235.15 1628.70 2209.30 3500.31 3550.77 4448.73 4453.12 4676.70 0.39036 M=13 1225.68 1621.08 2189.77 3465.19 3540.46 4443.05 4451.44 4658.80 0.41890 M=15 1225.44 1620.93 2189.49 3465.14 3540.21 4442.22 4451.90 4657.96 0.45036
##### 3.2. Frequency Comparisons

In this subsection, to validate the feasibility and effectiveness of the presented method, the free vibration analysis of stiffened cylindrical shell with general boundary conditions is performed by the presented method, and the obtained natural frequencies are compared with several state-of-the-art methods firstly. By gained confidence in the present method, we further investigate the effects of stiffeners or not on the natural frequencies. After that we the contribution of first-order shear deformation theory (FSDT) on moderate-thickness stiffened cylindrical shells in contrast with classic shear deformation theory (CSDT).

Classic boundary support can be considered as special cases of elastic boundary conditions. Different types of classic boundary support can be readily generated by adjusting the stiffness values of the translational springs and rotational springs which are uniformly distributed along the ends of cylindrical shell. Although we can gain the exact solutions for stiffened cylindrical shell with general boundary conditions, to facilitate the method comparisons, in this paper we select three typical boundary conditions that are frequently encountered in engineering practices: free-free (F-F), simple-simple (S-S), and clamped-clamped (C-C). Taking edge =0 as an example, the corresponding spring stiffness of three classical boundary conditions are introduced as follows:

At first, regarding the geometrical and material parameters as =39.45Ć10ā2 m, =4.9759Ć10ā2 m, =0.1651Ć10ā2 m, =2762 kg/m3, =0.3, =68.95 GPa, =15, N=5, =0.3175Ć10ā2m, =0.5334Ć10ā2 m, =1.9725Ć10ā2 m, axial wave number =1, 2, 3, and circumference wave number =1, 2, 3, 4, 5, the natural frequencies of cylindrical shell with S-S boundary conditions are studied. The frequency analysis results of polynomial Ritz method (PRM) [25], wave propagation method (WPM) [27], finite element method (FEM), and the present method are shown in Table 2. Then, with the F-F boundary condition and the geometrical and material parameters =47.09Ć10ā2 m, =10.37Ć10ā2 m, =0.119Ć10ā2 m, =7700 kg/m3, =0.3, =2.06 GPa, =15, =5, =0.218Ć10ā2 m, and =3.14Ć10ā2 m, through introducing dimensionless parameters , the corresponding analysis results are obtained by the present method are compared with that of Galletly Ritz method (GRM) [32], PRM [26], and FEM, which are illustrated in Table 3. Finally, to validate the generality of the presented method in elastic boundary conditions, we compare the presented method with a state-of-the-art Galerkin method (GM) [34]; the comparison results are shown in Tables 4ā6. Herein, the elastic boundary conditions are simulated by F-F, S-S and C-C boundary conditions, and the geometrical and material parameters are chosen as L=0.502 m, R=0.0635 m, h=0.00163 m, Ļ0=7800 kg/m3, Ī¼=0.28, E=2.1e+11 N/m2, M=15, N=5, ds=0.01 m, hs=0.003 m, and bs=L/7. From the above tables, it is obvious that the current solutions are consistent with that of referential methods, which validates the feasibility and effectiveness of the present method.

 n m=1 m=2 m=3 PRM WPM FEM present PRM WPM FEM present PRM WPM FEM present 1 1199 1216 1244 1223.49 3493 3536 3634 3435 5839 5907 6080 5943 2 1564 1635 1642 1596.01 2113 2176 2202 21535 3378 3430 3190 3555 3 4387 4578 4015 4356.80 4400 4573 4109 4439 4595 4788 4428 4527 4 8377 8781 7409 8105.72 8392 8731 7221 8468 8449 8728 8229 8538 5 13490 14172 13143 12979 13508 14119 13356 1387 13555 14069 13589 14389
 n GRM PRM FEM present 1 0.1997 0.1983 0.2065 0.2064 2 0.0851 0.0851 0.0895 0.0895 3 0.0689 0.0704 0.0742 0.0744 4 0.1092 0.1147 0.1123 0.1157 5 0.1728 0.1836 0.1569 0.1656
 =15, =5, =1015 FEM GM present 1 983.34 988.32 980.25 2 1718.24 1717.66 1715.78 3 1736.55 1738.45 1734.99 4 1947.46 1949.45 1942.42
 =15, =5, =0 FEM GM present 0 598.17 595.34 599.84 1 625.00 623.45 624.24 2 662.26 661.42 663.69 3 1098.20 1096.37 1099.48
 =1 =2 FEM GM present FEM GM present 1 1147.0 1143.5 1145.4 1257.7 1253.3 1255.5 2 1021.0 1019.5 1020.2 1681.8 1678.6 1679.2 3 1960.5 1958.3 1959.3 2214.5 2213.6 2212.1 4 3440.5 3435.6 3438.9 3795.2 3792.9 3794.3

Having gained confidence in the present method, considering the geometrical and material parameters as L=0.502 m, R=0.0635 m, h=0.00163 m, Ļ0=7800 kg/m3, Ī¼=0.28, E=2.1e+11 N/m2, M=15, N=5, ds=0.01 m, hs=0.003 m, and bs=L/7, we analyze the effects of stiffeners on the natural frequencies of cylindrical shells. Therein, by setting stiffener thickness hs as 0 in (13), the stiffened dynamics model is degraded as an unstiffened dynamics model, and the natural frequencies of unstiffened cylindrical shell are obtained. The frequency differences between cylinders with and without stiffening are shown in Figure 2. From the table, we can see that the natural frequencies have a distinct increase by adding stiffeners. To further investigate the stiffener thickness variation effects on natural frequencies, we further investigate the vibration characteristics of stiffened cylindrical shell under three stiffener thickness values hs= 0.291 m, 0.582 m, and 0.873 m. The solution results of WPM [27], modified variational method (MVM) [55], FEM, and the presented method have been presented in Table 7. In this case, the geometrical and material parameters of cylindrical shells are =47.09Ć10ā2 m, =10.37Ć10ā2 m, =0.119Ć10ā2 m, =7700 kg/m3, =0.3, =2.06 GPa, =15, =5, =0.218Ć10ā2 m, and =3.14Ć10ā2 m. It is obvious that the present solutions have excellent agreements with the results of several literatures. Besides, the natural frequencies of stiffened cylindrical shell are increasing along with the heightening of thickness value and vice versa.

 =0.291 =0.582 =0.873 MVM WPM FEM present MVM WPM FEM present MVM WPM FEM present 2 4550 4470 4559 4558.44 4580 4450 4660 4649.62 5040 4885 5432 5098.44 3 3870 3655 3897 3790.55 6710 6235 6778 6370.77 10330 9500 9564 9842.43 4 6550 5950 6435 6443.51 12830 11790 10098 11345.66 20200 18010 19797 18902.30 5 10000 9510 9667 9559.39 20120 19020 18890 19554.33 31800 25570 29905 30222.30

It is noted that all the numerical cases above are based on FSDT, the analysis results which can hardly show the significant difference with CSDT, since thin stiffened shells are considered in those examples. In this case, to show the contribution of FSDT, we present the examples of moderately thick stiffened cylindrical shells. The natural frequencies with larger parameters are presented in Tables 8 and 9. The geometrical and material parameters of cylindrical shells are L=0.502 m, R=0.0635 m, Ļ0=7800 kg/m3, Ī¼=0.28, E=2.1e+11 N/m2, M=15, N=5, ds=0.01 m, hs=0.003 m, and bs=L/7. When the thickness increases, the solution based on FSDT is more accurate than CSDT. It is obvious to see that the obtained natural frequency based on FSDT, especially the higher-order frequency, is more accurate than that of CSDT.

 CSDT FEM present 1 1095.4 1099.4 1098.8 2 1104.1 1105.9 1106.4 3 1792.3 1793.8 1794.3 4 2076.3 2079.9 2078.4 5 2434.7 2436.3 2437.1
 CSDT FEM present 1 1299.6 1301.6 1302.3 2 1314.2 1316.4 1317.5 3 1820.1 1822.5 1823.1 4 2254.9 2256.2 2257.0 5 2263.2 2266.3 2267.4

It is noted that the partial data in Tables 2ā6, 8, and 9 are reproduced from the reference [42] (under the Creative Commons Attribution License/public domain).

##### 3.3. Modal Comparisons

To verify the validity of the developed method on mode characteristics in free vibration analysis, we compare the mode shapes obtained by the presented method with that of FEM. Fourth-order mode shapes with three types of boundary conditions are presented. The comparisons of two approaches with the F-F, S-S, and C-C boundary conditions are depicted in Figures 3ā5, respectively. It is noted that the right columns are the results obtained by our methods and the left columns are by FEM method in all subfigures. From these figures, the vibration behaviors of the stiffened cylindrical shell can be inspected straightforwardly and an excellent consistency of FEM and the presented method can be concluded.

##### 3.4. Discussions

From the above convergence analysis, frequency comparisons, and modal comparisons, the feasibility and effectiveness of the proposed unified approach have been validated in the vibration analysis of stiffened cylindrical shell with general boundary conditions. Some discussions and findings from the vibration analysis are summarized as follows.

As shown in Table 1, it is obvious that the presented approach solution holds a high calculation efficiency and shows a great agreement with FEM. The excellent convergence of the presented method is attributed to the following key factors: (i) unlike the high-order shear deformation theory (HSDT) in theoretical formulation, the first-order shear deformation theory (FSDT) we adopted in this paper avoids redundant boundary condition terms and complex high-order integral calculations, which is conducive to simplify the computational tasks in theoretical modeling process; (ii) a unified approach combining the modified Fourier series and the Rayleigh-Ritz method only calculates the kinetic energy and potential energy in the motion process rather than the specific stress variation calculations; the infinite degree-of-freedom system is transformed into a generalized multidegree-of-freedom problem by the basis function, thus reducing the computational workloads. Therefore, the proposed method combining the FSDT with Rayleigh-Ritz procedure can efficiently improve the computational efficiency and saves calculation time in free vibration analysis of stiffened cylindrical shells.

As shown in Tables 2ā9 and Figures 2ā5, the solutions obtained by proposed method are consistent with that of several state-of-the-art methods. The accuracy virtues of the presented method might be attributed to the following: (i) by considering the effect of transverse shear deformation with the Reissner-Mindlin displacement hypothesis, the FSDT can build highly accurate theoretical model for stiffened cylindrical shell; (ii) the superposition of Fourier cosine series and supplementary function are chosen as the admissible displacement function in Rayleigh-Ritz procedure, which ensures the convergence of the second derivative and smooths the admissible displacement function adequately throughout the entire solution domain. Therefore, the proposed method can perform highly accurate vibration analysis for stiffened cylindrical shells.

In summary, the presented method greatly improves computational efficiency while maintaining computational accuracy and hereby is feasible and effective approach in the free vibration analysis of stiffened cylindrical shells.

#### 4. Conclusions

The objective of this paper is to propose a unified approach for free vibration analysis of stiffened cylindrical shell with general boundary condition. In view of the first-order shear deformation shell theory, the vibration analysis of stiffened cylindrical shell is modeled mathematically. Without requiring any special remodeling or procedures, the arbitrary boundary conditions including free, simple, and clamped supports are readily simulated by simply adjusting the equivalent spring stiffness. By choosing the superposition function of Fourier series and complementary function as admissible displacement function, the natural frequencies and modes are obtained by Rayleigh-Ritz procedure based on the energy expression of stiffened cylindrical shell. In contrast to several state-of-the-art methods, the feasibility and effectiveness of the presented method are examined and the corresponding vibration characteristics of stiffened cylindrical shells are required. Some conclusions have been summarized as follows:

In view of the fast computing time and high calculation efficiency, the presented approach solution shows an excellent convergence in free vibration analysis.

For stiffened cylindrical shells subject to general boundary conditions, the presented method possesses steady calculation results and high solution accuracy in free vibration analysis.

It is found that the natural frequencies have a distinct increase with adding several stiffeners. Besides, with the higher stiffener thickness, the natural frequencies show an increasing tendency.

For stiffened cylindrical shells with moderately thick, the first-order shear deformation theory which considering the transverse shear deformation holds higher computing accuracy than classic shear deformation theory.

The modal shapes acquired by the presented method show good consistent with that of FEM, which further validate the effectiveness of the method presented from modal shapes perspective.

The current efforts provide a unified accurate alternative to other analytical techniques and shed lights on the subsequent nonlinear vibration analysis.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

#### Acknowledgments

This paper is supported by National Natural Science Foundation of China (NNSFC) through Grants nos. 11672007, 11832002, and 11427801, the Academic Excellence Foundation of BUAA for PhD Students, and the Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality (PHRIHLB).

#### References

1. M. V. Chernobryvko, K. V. Avramov, V. N. Romanenko, T. J. Batutina, and U. S. Suleimenov, āDynamic instability of ring-stiffened conical thin-walled rocket fairing in supersonic gas stream,ā Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 230, no. 1, pp. 55ā68, 2016. View at: Publisher Site | Google Scholar
2. K. Liang, Y. Zhang, Q. Sun, and M. Ruess, āA new robust design for imperfection sensitive stiffened cylinders used in aerospace engineering,ā Science China Technological Sciences, vol. 58, no. 5, pp. 796ā802, 2015. View at: Publisher Site | Google Scholar
3. A. Barut, E. Madenci, A. Tessler, and J. H. Starnes Jr., āNew stiffened shell element for geometrically nonlinear analysis of composite laminates,ā Computers & Structures, vol. 77, no. 1, pp. 11ā40, 2000. View at: Publisher Site | Google Scholar
4. Y. Qu, Y. Chen, X. Long, H. Hua, and G. Meng, āA modified variational approach for vibration analysis of ring-stiffened conical-cylindrical shell combinations,ā European Journal of Mechanics - A/Solids, vol. 37, pp. 200ā215, 2013. View at: Publisher Site | Google Scholar | MathSciNet
5. V. Meyer, L. Maxit, Y. Renou, and C. Audoly, āExperimental investigation of the influence of internal frames on the vibroacoustic behavior of a stiffened cylindrical shell using wavenumber analysis,ā Mechanical Systems and Signal Processing, vol. 93, pp. 104ā117, 2017. View at: Publisher Site | Google Scholar
6. D. Bianco, F. P. Adamo, M. Barbarino et al., āIntegrated aero-vibroacoustics: the design verification process of Vega-C launcher,ā Applied Sciences (Switzerland), vol. 8, no. 88, 2018. View at: Google Scholar
7. R. Citarella and L. Federico, āAdvances in vibroacoustics and aeroacustics of aerospace and automotive systems,ā Applied Sciences (Switzerland), vol. 8, no. 366, 2018. View at: Google Scholar
8. A. M. J. Al-Najafi and G. B. Warburton, āFree vibration of ring-stiffened cylindrical shells,ā Journal of Sound and Vibration, vol. 13, no. 1, pp. 9ā25, 1970. View at: Publisher Site | Google Scholar
9. A. A. Jafari and M. Bagheri, āFree vibration of non-uniformly ring stiffened cylindrical shells using analytical, experimental and numerical methods,ā Thin-Walled Structures, vol. 44, no. 1, pp. 82ā90, 2006. View at: Publisher Site | Google Scholar
10. J. Wei, M. Chen, G. X. K. Hou, and N. Q. Deng, āWave based method for free vibration analysis of cylindrical shells with nonuniform stiffener distribution,ā Journal of Vibration and Acoustics, vol. 135, no. 6, Article ID 061011, 13 pages, 2013. View at: Google Scholar
11. N. D. Duc, āNonlinear dynamic response of imperfect eccentrically stiffened FGM double curved shallow shells on elastic foundation,ā Composite Structures, vol. 99, pp. 88ā96, 2013. View at: Publisher Site | Google Scholar
12. Y. Qu, H. Hua, and G. Meng, āA domain decomposition approach for vibration analysis of isotropic and composite cylindrical shells with arbitrary boundaries,ā Composite Structures, vol. 95, pp. 307ā321, 2013. View at: Publisher Site | Google Scholar
13. Z. Qin, F. Chu, and J. Zu, āFree vibrations of cylindrical shells with arbitrary boundary conditions: a comparison study,ā International Journal of Mechanical Sciences, vol. 133, pp. 91ā99, 2017. View at: Publisher Site | Google Scholar
14. A. R. Ghasemi and M. Mohandes, āFree vibration analysis of rotating fiber-metal laminate circular cylindrical shells,ā Journal of Sandwich Structures & Materials, vol. 21, no. 3, pp. 1009ā1031, 2019. View at: Publisher Site | Google Scholar
15. J. Zhao, K. Choe, Y. K. Zhang, A. L. Wang, C. H. Lin, and Q. S. Wang, āA closed form solution for free vibration of orthotropic circular general boundary conditions,ā Composites Part B-Engineering, vol. 159, pp. 447ā460, 2019. View at: Google Scholar
16. R. Poultangari and M. Nikkhah-Bahrami, āApplication of vectorial wave method in free vibration analysis of cylindrical shells,ā Advances in Applied Mathematics and Mechanics, vol. 9, no. 5, pp. 1145ā1161, 2017. View at: Publisher Site | Google Scholar | MathSciNet
17. K. Daneshjou, A. Nouri, and R. Talebitooti, āAnalytical model of sound transmission through laminated composite cylindrical shells considering transverse shear deformation,ā Applied Mathematics and Mechanics-English Edition, vol. 29, no. 9, pp. 1165ā1177, 2008. View at: Publisher Site | Google Scholar
18. Q. Wang, D. Shao, and B. Qin, āA simple first-order shear deformation shell theory for vibration analysis of composite laminated open cylindrical shells with general boundary conditions,ā Composite Structures, vol. 184, pp. 211ā232, 2018. View at: Publisher Site | Google Scholar
19. G. Y. Jin, T. Ye, X. L. Ma, Y. Chen, Z. Su, and X. Xie, āA unified approach for the vibration analysis of moderately thick composite laminated cylindrical shells with arbitrary boundary conditions,ā International Journal of Mechanical Sciences, vol. 75, pp. 357ā376, 2013. View at: Publisher Site | Google Scholar
20. L. W. Zhang, Z. X. Lei, and K. M. Liew, āVibration characteristic of moderately thick functionally graded carbon nanotube reinforced composite skew plates,ā Composite Structures, vol. 122, pp. 172ā183, 2015. View at: Publisher Site | Google Scholar
21. D. Liu, S. Kitipornchai, W. Chen, and J. Yang, āThree-dimensional buckling and free vibration analyses of initially stressed functionally graded graphene reinforced composite cylindrical shell,ā Composite Structures, vol. 189, pp. 560ā569, 2018. View at: Publisher Site | Google Scholar
22. S. Javed, āFree vibration characteristic of laminated conical shells based on higher-order shear deformation theory,ā Composite Structures, vol. 204, pp. 80ā87, 2018. View at: Publisher Site | Google Scholar
23. K. M. Liew, L. X. Peng, and S. Kitipornchai, āNonlinear analysis of corrugated plates using a FSDT and a meshfree method,ā Computer Methods Applied Mechanics and Engineering, vol. 196, no. 21-24, pp. 2358ā2376, 2007. View at: Publisher Site | Google Scholar
24. A. Bhar, S. S. Phoenix, and S. K. Satsangi, āFinite element analysis of laminated composite stiffened plates using FSDT and HSDT: a comparative perspective,ā Composite Structures, vol. 92, no. 2, pp. 312ā321, 2010. View at: Publisher Site | Google Scholar
25. A. A. Jafari and M. Bagheri, āFree vibration of rotating ring stiffened cylindrical shells with non-uniform stiffener distribution,ā Journal of Sound and Vibration, vol. 296, no. 1-2, pp. 353ā367, 2006. View at: Publisher Site | Google Scholar
26. S. A. Moeini, M. Rahaeifard, M. T. Ahmadian, and M. R. Movahhedy, āFree vibration analysis of functionally graded cylindrical shells stiffened by uniformly and non-uniformly distributed ring stiffeners,ā in Proceedings of the ASME International Mechanical Engineering Congress and Exposition (IMECE '09), pp. 367ā375, Florida, Fla, USA, November 2009. View at: Google Scholar
27. L. Gan, X. B. Li, and Z. Zhang, āFree vibration analysis of ring-stiffened cylindrical shells using wave propagation approach,ā Journal of Sound and Vibration, vol. 326, no. 3ā5, pp. 633ā646, 2009. View at: Publisher Site | Google Scholar
28. M. Chen, J. Wei, K. Xie, N. Deng, and G. Hou, āWave based method for free vibration analysis of ring stiffened cylindrical shell with intermediate large frame ribs,ā Shock and Vibration, vol. 20, no. 3, pp. 459ā479, 2013. View at: Publisher Site | Google Scholar
29. C.-C. Liu, F.-M. Li, Z.-B. Chen, and H.-H. Yue, āTransient wave propagation in the ring stiffened laminated composite cylindrical shells using the method of reverberation ray matrix,ā The Journal of the Acoustical Society of America, vol. 133, no. 2, pp. 770ā780, 2013. View at: Publisher Site | Google Scholar
30. M. Liu, J. Liu, and Y. Cheng, āFree vibration of a fluid loaded ring-stiffened conical shell with variable thickness,ā Journal of Vibration and Acoustics, vol. 136, no. 5, Article ID 051003, 2014. View at: Google Scholar
31. X. Wang, W. Wu, and X. Yao, āStructural and acoustic response of a finite stiffened conical shell,ā Acta Mechanica Solida Sinica, vol. 28, no. 2, pp. 200ā209, 2015. View at: Publisher Site | Google Scholar
32. G. D. Galletly, āOn the in-vacuo vibrations of simply supported, ring-stiffened cylindrical shells,ā in Proceedings of the 2nd National Congress of Applied Mechanics Processing, pp. 1ā9, Ann Arbor, Mich, USA, March 1958. View at: Google Scholar
33. Y.-S. Lee and Y.-W. Kim, āVibration analysis of rotating composite cylindrical shells with orthogonal stiffeners,ā Computers & Structures, vol. 69, no. 2, pp. 271ā281, 1998. View at: Publisher Site | Google Scholar
34. W. Zhang, Z. Fang, X.-D. Yang, and F. Liang, āA series solution for free vibration of moderately thick cylindrical shell with general boundary conditions,ā Engineering Structures, vol. 165, pp. 422ā440, 2018. View at: Publisher Site | Google Scholar
35. A. Y. Tamijani and R. K. Kapania, āVibration analysis of curvilinearly-stiffened functionally graded plate using element free galerkin method,ā Mechanics of Advanced Materials and Structures, vol. 19, no. 1-3, pp. 100ā108, 2012. View at: Publisher Site | Google Scholar
36. G. G. Sheng and X. Wang, āThe dynamic stability and nonlinear vibration analysis of stiffened functionally graded cylindrical shells,ā Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 56, pp. 389ā403, 2018. View at: Publisher Site | Google Scholar | MathSciNet
37. P. Hao, B. Wang, G. Li, Z. Meng, K. Tian, and X. H. Tang, āHybrid optimization of hierarchical stiffened shells based on smeared stiffener method and finite element method,ā Thin-Walled Structures, vol. 82, pp. 46ā54, 2014. View at: Publisher Site | Google Scholar
38. B. G. Prusty and S. K. Satsangi, āFinite element transient dynamic analysis of laminated stiffened shells,ā Journal of Sound and Vibration, vol. 248, no. 2, pp. 215ā233, 2001. View at: Publisher Site | Google Scholar
39. C. Liu, Q. Tian, and H. Hu, āNew spatial curved beam and cylindrical shell elements of gradient-deficient absolute nodal coordinate formulation,ā Nonlinear Dynamics, vol. 70, no. 3, pp. 1903ā1918, 2012. View at: Publisher Site | Google Scholar | MathSciNet
40. C. Liu, Q. Tian, and H. Hu, āDynamics of a large scale rigid-flexible multibody system composed of composite laminated plates,ā Multibody System Dynamics, vol. 26, no. 3, pp. 283ā305, 2011. View at: Publisher Site | Google Scholar
41. M. Vaziri Sereshk and M. Salimi, āComparison of finite element method based on nodal displacement and absolute nodal coordinate formulation (ANCF) in thin shell analysis,ā International Journal for Numerical Methods in Biomedical Engineering, vol. 27, no. 8, pp. 1185ā1198, 2011. View at: Publisher Site | Google Scholar
42. L. T. B. Nam, N. M. Cuong, T. I. Thinh, T. T. Dat, and V. D. Trung, āContinuous element formulations for composite ring-stiffened cylindrical shells,ā Vietnam Journal of Science and Technology, vol. 56, no. 4, pp. 515ā530, 2018. View at: Publisher Site | Google Scholar
43. O. D. de Matos Junior, M. V. Donadon, and S. G. P. Castro, āAeroelastic behavior of stiffened composite laminated panel with embedded SMA wire using the hierarchical RayleighāRitz method,ā Composite Structures, vol. 181, pp. 26ā45, 2017. View at: Publisher Site | Google Scholar
44. G. Eslami and M. Z. Kabir, āMulti-objective optimization of orthogonally stiffened cylindrical shells using optimality criteria method,ā Scientia Iranica, vol. 22, no. 3, pp. 717ā727, 2015. View at: Google Scholar
45. L. Liu, D. Q. Cao, and S. P. Sun, āVibration analysis for rotating ring-stiffened cylindrical shells with arbitrary boundary conditions,ā Journal of Vibration and Acoustics-Transactions of the ASME, Article ID 061010, 10 pages, 2013. View at: Google Scholar
46. Y.-W. Kim and Y.-S. Lee, āTransient analysis of ring-stiffened composite cylindrical shells with both edges clamped,ā Journal of Sound and Vibration, vol. 252, no. 1, pp. 1ā17, 2002. View at: Publisher Site | Google Scholar
47. Y. Y. Chai, Z. G. Song, and F. M. Li, āInvestigations on the aerothermoelastic properties of composite laminated cylindrical shells with elastic boundaries in supersonic airflow based on the Rayleigh-Ritz method,ā Aerospace Science and Technology, vol. 82, pp. 534ā544, 2018. View at: Google Scholar
48. S. Sun, D. Cao, and Q. Han, āVibration studies of rotating cylindrical shells with arbitrary edges using characteristic orthogonal polynomials in the Rayleigh-Ritz method,ā International Journal of Mechanical Sciences, vol. 68, pp. 180ā189, 2013. View at: Publisher Site | Google Scholar
49. W. L. Li, āVibration analysis of rectangular plates with general elastic boundary supports,ā Journal of Sound and Vibration, vol. 273, no. 3, pp. 619ā635, 2004. View at: Publisher Site | Google Scholar
50. Y. Zhang, J. Du, T. Yang, and Z. Liu, āA series solution for the in-plane vibration analysis of orthotropic rectangular plates with elastically restrained edges,ā International Journal of Mechanical Sciences, vol. 79, pp. 15ā24, 2014. View at: Publisher Site | Google Scholar
51. L. Dai, T. Yang, J. Du, W. L. Li, and M. J. Brennan, āAn exact series solution for the vibration analysis of cylindrical shells with arbitrary boundary conditions,ā Applied Acoustics, vol. 74, no. 3, pp. 440ā449, 2013. View at: Publisher Site | Google Scholar
52. L.-K. Song, C.-W. Fei, G.-C. Bai, and L.-C. Yu, āDynamic neural network method-based improved PSO and BR algorithms for transient probabilistic analysis of flexible mechanism,ā Advanced Engineering Informatics, vol. 33, pp. 144ā153, 2017. View at: Publisher Site | Google Scholar
53. L. K. Song, G. C. Bai, C. W. Fei, and R. P. Liem, āTransient probabilistic design of flexible multibody system using a dynamic fuzzy neural network method with distributed collaborative strategy,ā Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, pp. 1ā14, 2018. View at: Google Scholar
54. B. Rong, X. Rui, G. Wang, and F. Yang, āNew efficient method for dynamic modeling and simulation of flexible multibody systems moving in plane,ā Multibody System Dynamics, vol. 24, no. 2, pp. 181ā200, 2010. View at: Publisher Site | Google Scholar | MathSciNet
55. N. L. Basdekam and M. Chi, āResponse of oddly stiffened circular cylindrical shells,ā Journal of Sound and Vibration, vol. 17, no. 2, pp. 187ā206, 1971. View at: Publisher Site | Google Scholar

Copyright © 2019 Xue-Qin Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.