Abstract

In this paper, the Ritz method is adopted to investigate the vibration characteristics of isotropic moderately thick annular spherical shell with general boundary conditions. The energy expressions of the annular spherical shell were established based on the first-order shear deformation theory (FSDT). The spring stiffness method is introduced to guarantee continuity and simulate various boundary conditions on the basis of the domain decomposition method. Under the current framework, the displacement admissible function along axial direction and circumferential direction of the shell structure are, respectively, expanded as the unified Jacobi polynomials and Fourier series. The final solutions can be obtained according to the Ritz method. The validity of the proposed method is proved by comparing the results of the same condition with those obtained by the finite element method (FEM) and published literatures. The results show that the current method has fast convergence and delightful accuracy through the comparative study. On this basis, the vibration characteristics of isotropic moderately thick annular spherical shell are further studied by a series of numerical examples.

1. Introduction

The isotropic moderately thick annular spherical shell structure, as a basic component, has been widely used in many fields of engineering, such as drive shaft, cooling towers, water tanks, pressure vessels, and rotor system. In practical engineering, the dynamical behavior and stability are very important for the structural design because of the structure exposed to various complex environments and bear various dynamic loads. Based on this background, it is very important and necessary to investigate the free vibration characteristics of isotropic moderately thick annular spherical shell under general boundary conditions.

For vibration analysis of thin spherical shell, Bryan [1] developed a new method to investigate the differential equations of motion about a thin spherical shell, in which the Legendre polynomials are utilized to obtain the analytical solutions for the spatial differential equation. Xie et al. [2] proposed a unified approach to study the vibration characteristics of spherical of revolution stiffened by rings with T cross-section. Based on Flugge’s thin shell theory, Wang et al. [3] presented the vibration characteristics of coupled doubly curved shell structures which subject to arbitrary boundaries. Pang et al. [4] applied the Rayleigh–Ritz method to analyze the free vibration of doubly curved shells. Li and Yuan [5, 6] applied quasi-Green’s function method to analyze the free vibration of shallow spherical shell based on the fundamental solution and boundary equation of the problem. By using a hybrid finite element method, Menaa and Lakis [7] presented the free vibration characteristics of the spherical shell by uniting the thin shell theory and classical finite element method. By using the separation of variable method, Zaera et al. [8] obtained the free vibration solutions of closed thin spherical. Polyakov et al. [9] derived the formulation of a closed spherical shell to analyze the free vibrations along the thickness direction of structure. According to von Karman's thin shallow shells theory, Thomas et al. [10, 11] derived the nonlinear vibration equations of shallow spherical shell under large amplitude displacement.

For vibration analysis related to moderately thick spherical shell structures, Tornabene et al. [1216] applied the Generalized Differential Quadrature (GDQ) method to investigate the free vibration characteristics of the rotating shell structure based on FSDT. Buchanan and Rich [17] figured out the frequency parameters of thick isotropic spherical shells with simple supports by formulating Lagrangian finite element in spherical coordinates. Based on FSDT, Wang et al. [1821] carried out a great deal of research studies on spherical shell structures under different boundary conditions, and the analytical solutions were obtained by the Ritz method. Jin et al. [2225] also conducted a large number of studies to analyze free vibration characteristics of spherical shell structures under general boundary conditions based on the Rayleigh–Ritz method, which each admissible functions is invariantly expanded as a modified Fourier series according to the three-dimensional shell theory of elasticity. Hosseini-Hashemi and Fadaee [26] applied the separation of variables method for free vibration analysis of moderately thick spherical shell according to the first-order shear deformation theory. Zenkour [27] expressed the dynamic responses of anisotropic spherical shells under a uniformly distributed transverse load in the frame work of mixed first-order shear deformation theory. Based on higher order shear deformation theory (HSDT), Fazzolari [28] presented the free vibration characteristics of doubly curved laminated composite spherical shallow shells by using the dynamic stiffness method (DSM). Panda et al. [2932] investigated the nonlinear free vibration characteristics of laminated composite shallow spherical shell using the variational method, and the analytical model was established on the basis of Green–Lagrange nonlinear kinematics. For space reasons, other papers on rotating cylindrical and spherical shells can be found in references [3335].

It can be concluded that lots of published literatures focused on the free vibration response of thin spherical shell which subject to general edge restraints, and the free vibration characteristics of moderately thick spherical shell structures have been investigated in some research studies. However, to the authors’ knowledge, the aforementioned review for vibration solutions of moderately thick annular spherical shell structure with general boundary conditions based on the FSDT is infrequent. Considering this background, it is necessary to analyze the free vibration characteristics of moderately thick annular spherical shell structure with general boundary conditions. In addition, this paper presents a series of unpublished numerical results of the effects on the boundary conditions, shell segments, and structural parameters, which can be used as basic data to provide comparative data for future researchers in this field.

2. Theoretical Formulations

2.1. The Mathematical Model of Moderately Thick Annular Spherical Shell

An isotropic moderately thick annular spherical shell with radius R is considered, as shown in Figure 1. The deformation displacements of the annular spherical shell with respect to the coordinate system can be, respectively, defined by u, , and in the φ, θ, and δ directions, and the system (φ, θ, δ) represent axial, circumferential, and normal directions, accordingly. Cs and h, respectively, represent the geometry center and thickness of the annular spherical shell. φ0 and φ1, respectively, denote the center angle correspond to the top and bottom of annular spherical shell. The formulas of this paper are derived on the basis of the domain decomposition method [3641], the spherical shell is divided into H segments along axial direction, and the artificial springs are arranged at both ends of each segment; the continuity condition and various boundary conditions can be simulated by assigning the stiffness values of springs.

2.2. Energy Expressions of Moderately Thick Annular Spherical Shell

According to FSDT [4244], the displacements of ith segment of moderately thick annular spherical shell can be written as follows:

The strains of moderately thick annular spherical shell in this study can be written as follows:where , , , , , , , and can, respectively, be expressed as below:

In this paper, the symbols A and B are Lamé parameters, and they can be expressed as follows [45, 46]:

The stresses of the structure can be written as follows:where symbols of and , respectively, represent normal and shear stresses. The Qij(δ) are defined as follows:where E and μ are elastic modulus and Poisson’s ratio, respectively. The force and moment resultants can be expressed as below:where represents the shear correction factor, and it is set as 5/6 in this study. The symbols of , , and can be written as follows:

The strain energy of the select segment can be expressed as follows:

The strain energy expression can be rewritten as , and the detailed expressions are shown in equations (10)–(12):

The kinetic energy of the select segment can be written as follows:where the dot on the symbols represents differentiation about time, and the related symbols can be expressed as follows:

As mentioned above, the application of artificial springs can ensure the convergence of result. Each edge of the structure is restrained by three linear springs (ku, , ) and two rotational springs (kφ, kθ). Then, the boundary potential energy for moderately thick annular spherical shell is expressed as follows:where and , respectively, denote the spring stiffness value of the moderately thick annular spherical shell.

The potential energy stored in two adjacent segments can be shown as follows:

Thus, the total potential energy can be written as follows:

2.3. Admissible Displacement Functions

In this paper, the authors try to introduce the unified Jacobi polynomials to improve the selection of displacement functions based on the domain decomposition method, and the value of φ are defined in the range of [−1, 1]. The related recurrence formulas are shown as below:where and

Therefore, the displacement functions can be generalized according to Jacobi polynomial:where , , , , and are unknown coefficients; m and n, respectively, represent the semiwave number in axial and circumferential direction, and M is highest degrees of m.

2.4. Solution Procedure

The Lagrangian energy functions can be expressed as follows:

Differentiate with regard to undetermined coefficients, and the formula can be obtained:

Substituting the formulas mentioned above into equation (21), the matrix form can be shown as follows:where K, M, and Q of equation (22) represent stiffness matrix, mass matrix, and unknown coefficients matrix.

3. Numerical Results and Discussion

According to the above research, a unified analytical model for moderately thick annular spherical shell is established. The main purpose of this section is to discuss the numerical results. According to the experience of published literatures, the various boundary conditions can be represented by the first letter of a word. For example, the free and elastic edge conditions are denoted by the letter F and Ei (i = 1, 2, 3).The material properties are chosen as , , and . In addition, other related parameters appeared in this paper are defined as M = 8, α = β = 0.5, H = 6; R = 1 m, φ0 = π/6, φ1 = π/2, and h = 0.06 m. The results of this paper are handle by .

3.1. Convergence Study

In order to investigate the convergence of different spring stiffness values, Figure 2 illustrates the frequency parameters of moderately thick annular spherical shell with different boundary parameters. No matter for boundary spring and connective spring, it is obvious that the edge condition changes from free to clamped case with the spring stiffness values in range of 10−11E to 1010E. For example, the boundary condition is clamped when spring stiffness values are defined in 103E to 1010E. Thus, the various boundary conditions can be easily obtained through the analysis. All edge conditions used in this study are shown in Table 1.

Figure 3 shows results of moderately thick annular spherical shell with different segment numbers, and the boundary condition is set as clamped-clamped case. The example shows clearly that the current method converges quickly with the increase of H, that is to say, it is not necessary for very high value of H. Considering the solution precision of moderately thick annular spherical shell, the value of H is set as 6 in this paper.

The results of moderately thick annular spherical shell with different truncation numbers are shown in Figure 4. The example shows clearly that the current method converges quickly with the increase of M. To obtain the accurate solution of moderately thick annular spherical shell, the value of M is defined as 8 in this study.

Figure 5 displays the results of moderately thick annular spherical shell with different Jacobi parameters α and β. The ordinate in this example is the relative percentage value, and the results of α=β= 0.5 are selected as the comparative object. It is apparently that the different Jacobi parameters have little effect on free vibration of moderately thick annular spherical shell. That is to say, all the Jacobi polynomials are capable to construct the displacement functions, and it can be regarded as one of the advantages of the current method.

3.2. Free Vibration Behavior of Moderately Thick Annular Spherical Shell Structure

Table 2 shows the natural frequency of moderately thick annular spherical shell with classical boundary conditions, and all the results are compared with those obtained in commercial FEM software ABAQUS. The geometrical parameters of model is defined as φ0 = π/6, φ1 = π/2, and h = 0.06 m, and the meshes size are set as 0.01 m. Table 3 exhibits the results of moderately thick annular spherical shell with classical boundary conditions, and the results are defined as . The geometrical parameters are defined as φ0 = π/8, φ1 = π/2, h/R = 0.05, and R = 1 m, and all the results are compared with Qu et al. [37]. From the above research, it is clear that the proposed method is very accurate in solving the free vibration of moderately thick annular spherical shell.

Tables 4 and 5 show the results of moderately thick annular spherical shell, respectively, which subject to classical and elastic edge restraints. Figure 6 displays the frequency parameters Ω under different boundary conditions which varies with the increasing of mode of n. From the examples, it is obviously that the different boundary conditions have a clear effect on the results of moderately thick annular spherical shell. To improve the understanding on free vibration characteristics of moderately thick annular spherical shell, some mode shapes under different boundary conditions are given in Figure 7.

Tables 6 and 7 and Figure 8 display the results with different h/R ratios in moderately thick annular spherical shell. From the examples, it is obvious that that the results have a tendency to increase with h/R increasing. Not surprisingly, the frequency parameters Ω increased obviously with the increase of mode of m. Tables 8 and 9 and Figure 9 exhibit the results with different φ0 in moderately thick annular spherical shell. In addition, the results with different φ1 in moderately thick annular spherical shell are shown in Tables 10 and 11, in which φ0 = π/8. From the examples, it is remarkable that the frequency parameters Ω tend to increase with the increase of φ0, and the results usually decrease monotonously as the φ1 increases, that is to say higher frequency parameters can be obtained as the structural scale decreases. Beyond that we can also see that the different boundary conditions have apparent effects on the frequency parameters Ω.

4. Conclusions

This paper introduced a unified computational Ritz method to investigate free vibration characteristics of isotropic moderately thick annular spherical shell with general boundary conditions. The energy expressions of the annular spherical shell were established based on FSDT. The spring stiffness method is introduced to guarantee continuity and simulate various boundary conditions on the basic of the domain decomposition method. Beyond that the unified Jacobi polynomials and Fourier series are, respectively, applied to represent the displacement admissible function along axial direction and circumferential direction of the structure. The final solutions can be obtained according to the Ritz method. In order to prove the validity of the proposed method, the results of the same condition are compared with those obtained by the finite element method and published literatures. The results show that the current method has fast convergence and delightful accuracy through the comparative study. The most discovery of the current method is that the selection of admissible displacement functions is generalized by applying Jacobi polynomial. In addition, the effects on the boundary conditions, shell segments, and structural parameters are presented through a large number of numerical examples, which can be used as basic data to provide reference for future staff of such research in this field.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

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

Acknowledgments

This study was funded by National Key Research and Development Program (2016YFC0303406), Key Research and Development Program of Shandong Province (2019JZZY010125), National Key Laboratory on Ship Vibration & Noise Fund of China (6142204190207), National Natural Science Foundation of China (51709063), and Ph.D. Student Research and Innovation Fund of the Fundamental Research Funds for the Central Universities (HEUGIP201902).