Mathematical Problems in Engineering
Volume 2019 (2019), Article ID 4157930, 14 pages
https://doi.org/10.1155/2019/4157930

Research Article

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

Xue-Qin Li,¹ Wei Zhang,¹ Xiao-Dong Yang,¹ and Lu-Kai Song²

¹College of Mechanical Engineering, Beijing University of Technology, Beijing 100124, China
²School of Energy and Power Engineering, Beihang University, Beijing 100191, China

Correspondence should be addressed to Xiao-Dong Yang; jxdyang@163.com

Received 12 April 2019; Revised 16 June 2019; Accepted 17 June 2019; Published 10 July 2019

Academic Editor: Roberto G. Citarella

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.

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.

Figure 1: Schematic diagram 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 as

where , , 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 as

where

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 as

where 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 as

where 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 as

where , , 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 as

where 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 as

where the strain potential energy is

in which

Besides, 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 series

where 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 as

The 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 as

where

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 10⁵ while others are set as infinite large. The geometrical and material parameters of cylindrical shells are set as =39.45×10⁻² m, =4.9759×10⁻² m, =0.1651×10⁻² m, =2762 kg/m³, =0.3, =68.95 GPa, =0.3175×10⁻²m, =0.5334×10⁻² m, =1.9725×10⁻² 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.

Table 1: Comparison of natural frequencies of stiffened cylindrical shells with free boundaries.

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⁻² m, =4.9759×10⁻² m, =0.1651×10⁻² m, =2762 kg/m³, =0.3, =68.95 GPa, =15, N=5, =0.3175×10⁻²m, =0.5334×10⁻² m, =1.9725×10⁻² 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⁻² m, =10.37×10⁻² m, =0.119×10⁻² m, =7700 kg/m³, =0.3, =2.06 GPa, =15, =5, =0.218×10⁻² m, and =3.14×10⁻² 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, ρ₀=7800 kg/m³, μ=0.28, E=2.1e+11 N/m², 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.

Table 2: Comparison of natural frequencies of stiffened cylindrical shells with simple support boundaries.

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

Table 3: Comparison of natural frequencies of stiffened cylindrical shells with free boundaries with m=1.

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

Table 4: Comparison of natural frequencies of cylindrical shells with clamped boundaries.

=15, =5, =10¹⁵

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

Table 5: Comparison of natural frequencies of cylindrical shells with free boundaries.

=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

Table 6: Comparison of natural frequencies (HZ) of cylindrical shells with simply supported boundaries.

=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, ρ₀=7800 kg/m³, μ=0.28, E=2.1e+11 N/m², 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⁻² m, =10.37×10⁻² m, =0.119×10⁻² m, =7700 kg/m³, =0.3, =2.06 GPa, =15, =5, =0.218×10⁻² m, and =3.14×10⁻² 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.

Table 7: Comparison of natural frequencies by different methods.

=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

(a) F-F

(b) C-C

(c) S-S

Figure 2: Frequency differences of cylinders with and without stiffening.

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, ρ₀=7800 kg/m³, μ=0.28, E=2.1e+11 N/m², 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.

Table 8: Comparison of natural frequencies of stiffened cylindrical shells with clamped boundaries = 0.00326 m.

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

Table 9: Comparison of natural frequencies of stiffened cylindrical shells with clamped boundaries =0.005 m.

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.

(a) First-order

(b) Second-order

(c) Third-order

(d) Fourth-order

Figure 3: First four modes with free boundary.

(a) First-order

(b) Second-order

(c) Third-order

(d) Fourth-order

Figure 4: First four modes with clamped boundary.

(a) First-order

(b) Second-order

(c) Third-order

(d) Fourth-order

Figure 5: First four modes with simply supported boundary.

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

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.

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.

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.

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.

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.

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.

R. Citarella and L. Federico, “Advances in vibroacoustics and aeroacustics of aerospace and automotive systems,” Applied Sciences (Switzerland), vol. 8, no. 366, 2018.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

S. Javed, “Free vibration characteristic of laminated conical shells based on higher-order shear deformation theory,” Composite Structures, vol. 204, pp. 80–87, 2018.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.