Abstract

This paper presents elasticity solutions for the vibration analysis of isotropic and orthotropic open shells and plates with arbitrary boundary conditions, including spherical and cylindrical shells and rectangular plates. Vibration characteristics of the shells and plates have been obtained via a unified three-dimensional displacement-based energy formulation represented in the general shell coordinates, in which the displacement in each direction is expanded as a triplicate product of the cosine Fourier series with the addition of certain supplementary terms introduced to eliminate any possible jumps with the original displacement function and its relevant derivatives at the boundaries. All the expansion coefficients are then treated equally as independent generalized coordinates and determined by the Rayleigh-Ritz procedure. To validate the accuracy of the present method and the corresponding theoretical formulations, numerical cases have been compared against the results in the literature and those of 3D FE analysis, with excellent agreements obtained. The effects of boundary conditions, material parameters, and geometric dimensions on the frequencies are discussed as well. Finally, several 3D vibration results of isotropic and orthotropic open spherical and cylindrical shells and plates with different geometry dimensions are presented for various boundary conditions, which may be served as benchmark solutions for future researchers as well as structure designers in this field.

1. Introduction

In many applications and fields of modern technology, for example, civil aviation, space industry, and deep-ocean exploitation, the engineering structures often work in complex environment conditions and can be subjected to various forms of external loads, which may result in violent vibrations and lead the structures to failure. Therefore, it is of particular importance to understand the structural vibrations and reduce them through proper design to ensure a reliable, safe, and lasting structural performance. An important step in vibration design of an engineering structure is the evaluation of its vibration modal characteristics. This modal information plays a key role in the structure design and vibration suppression when subjected to dynamics excitations.

Plate and shells, such as rectangular plates and cylindrical, conical, and spherical shells, are basic structural elements of most engineering structures. A thorough understanding of their vibration characteristics is of great significance for engineers to predict the vibrations of the whole structures. Plates can be viewed as shells with zero curvature (infinite radii of curvature); thus, for simplicity, the plates are treated as a special case of shells in subsequent introductions. Recently, the vibration analysis of shells has received much attention. The literature on this subject is vast. In the literature, many of the studies are based on the thin-shell theories (such as Donnell’s, Love’s, and Sanders’ shell theories), which are developed on the basis of the Kirchhoff-Love’s simplifying assumptions [110]. Much more information on thin-shell theories is available in an influential monograph published by Leissa [11] in 1973. The thin-shell theories are very efficient in calculating vibration results of thin shells; however, they are highly inadequate for the analysis of even slightly thick shells due to the fact that the transverse stress and strain components are omitted. Therefore, the vibration of thick shells has conventionally been solved using the first-order shear deformation theories (FSDTs) and the higher-order shear deformation theories (HSDTs) [1219]. However, shear correction factors have to be incorporated in the FSDTs to adjust the transverse shear stiffness due to the fact that the shear strains along the shell thickness are assumed to be constant. A remedy to this drawback is to use the HSDTs, which are based on the assumption of nonlinear strain variation through the shell thickness. The HSDTs, however, are still inadequate for the analysis of thick shells since the transverse normal stress and strain components are ignored in these theories. To analyze thick shells, the three-dimensional (3D) theory of elasticity which accounts for all the transverse stress and strain components may be the best choice.

In recent years, vibration analysis of shells relying on the 3D theory of elasticity has attracted the attention of many researchers. The 3D analysis not only provides realistic results but also allows solving the whole spectrum of thick shells. Some research works in this research topic include those of Loy and Lam [20] using a layer-wise approach to study the vibration of thick cylindrical shells with simply supported-simply supported and clamped-clamped boundary conditions. Liew and his coauthors [2123] had developed a -Ritz method to calculate vibration results of solid cylinders, spherical panels, open cylindrical shells, and curved panels. In this method, the admissible displacements of a structure are expressed in terms of sets of one- and two-dimensional orthogonal polynomial functions and solved by the Ritz method. Leissa and his coworkers [2426] carried out a 3D method for determining the free vibration frequencies and mode shapes of thick shells with variable thickness. In this approach, displacement components in the meridional and normal directions are taken to be algebraic polynomials, respectively. Strain and kinetic energies of the shells are formulated, and upper bound values of the frequencies are obtained by minimizing the energy functional with respect to the expanded coefficients. Also, Zhou et al. [27, 28] investigated the vibration characteristics of cylindrical shells and doubly curved panels by the Chebyshev-Ritz method, in which the displacement in each direction is taken as a triplicate product of the Chebyshev polynomials in shell dimensions, multiplied by a boundary function along with a set of generalized coefficients determined by Ritz procedure. Other contributions can be seen in [2934]. The development of researches on this topic can be seen in several monographs by Leissa and Qatu [35], Qatu [36], Reddy [37], and Saada [38] and reviews [3941].

Compared with the above researches, it is observed that there exist a few investigations concerning 3D vibration analysis of closed shells and shallow shells with particular geometrical configurations, such as cylindrical, conical, and spherical shells. However, to the best of our knowledge, a unified 3D elasticity solution for isotropic and orthotropic open shells with arbitrary boundary conditions and variable circumferential dimensions does not seem to exist. In addition, most of the previous researches were limited to isotropic shells with a specified set of boundary conditions (such as free-fee, simply supported-simply supported, and clamped-clamped), and thus they typically require constant modifications of the trial functions and corresponding solution procedures to adapt to different boundary cases. Therefore, the use of the existing solution procedures can be easily inundated with various boundary conditions in practical applications due to the fact that there are hundreds of sets of boundary conditions for an open shell. Moreover, the open shells in the engineering applications may be of recognizable thickness-length/thickness-radius ratios and have large curvature (small radius of curvature). Therefore, a unified, reliable, and efficient method for predicting the 3D vibration characteristics of thick deep and shallow open shells with arbitrary boundary conditions would be highly desirable.

In view of these voids, this paper presents an endeavor to investigate the vibrations of isotropic and orthotropic open shells with arbitrary boundary conditions and variable circumferential dimensions in the framework of 3D shell theory.

In the present paper, the 3D shell theory and the energy based Rayleigh-Ritz procedure are combined to develop the elasticity solutions. Since there are many types of open shells, and it is impossible to undertake an all-encompassing survey of every case of them, therefore, the present work is focused on the open cylindrical and spherical shells, which are most frequently used in the engineering practices. Under the current framework, the displacement in each shell direction is expanded as a triplicate product of the cosine Fourier series with the addition of certain supplementary terms introduced to eliminate any possible jumps with the original displacement function and its relevant derivatives at the boundaries. All expansion coefficients are then determined by the Rayleigh-Ritz procedure. Based on the proposed method, numerical cases are presented and compared against results in the literature to validate the accuracy of this method. Effects of boundary conditions, geometric dimensions, and material properties on the natural frequencies of open shells and plates are discussed as well. Finally, several 3D vibration results of isotropic and orthotropic open cylindrical and spherical shells and plates with different geometry dimensions are presented for various boundary conditions, which may serve as benchmark solution for future researches to evaluate the new 2D shell theories and to compare results obtained by approximate numerical methods.

2. Theoretical Formulations

2.1. The Model

As shown in Figure 1(a), an open shell of meridional length , circumferential width , and uniform thickness are selected as the analysis model. An orthogonal shell coordinate system (, , and ) located in the middle surface of the shell is considered to present the geometry dimensions and its deformations, in which the coordinates , , and are taken in the meridional, circumferential, and normal directions, respectively. and are the mean radii of curvature in the and directions. The displacement variations in the meridional, circumferential, and normal directions of an arbitrary point in the shell space are represented by , , and , respectively.

2.2. Kinematic Relations and Stress Resultants

According to the three-dimensional shell theory of elasticity, the general strain-displacement relations for a shell can be defined as follows [33, 38]:where , , , , , and are the normal and shear strains. , , and indicate the Lamé parameters; they can be derived from the rectangular coordinate system (, , and ) by the following relationships: It should be stressed that the Lamé parameters given in (2) are determined by the adopted orthogonal curvilinear coordinate system. Consequently, the strains given in (1) depend on the adopted coordinate system as well. Since the present work is focused on the open cylindrical and spherical shells and plates, which are most frequently used in the engineering practices, therefore, under the current framework, the rectangular, cylindrical, and spherical coordinate systems are chosen to develop the theoretical formulations for the plates and open cylindrical and spherical shells, respectively. According to Figures 1(b)1(d), the coordinates (, , and ), geometry dimensions (, , , and ), and Lamé parameters (, , and ) of the plates and shells under consideration are [38](a)plates: , , , , , and ,(b)cylindrical shell: , , , , , , , , and ,(c)spherical shell: , , , , , , , , and ,where is the included angle in the circumferential direction of open cylindrical and spherical shells. is determined by the ends of a spherical shell (i.e., ). The corresponding strain-displacement relations for the plates and open cylindrical and spherical shells are given in Appendix A.

For open shells and plates, the three-dimensional constitutive relations are given by Hooke’s law aswhere , , , , , and are the normal and shear stresses. And the elastic stiffness coefficients () for the orthotropic materials are given aswhere , , and represent Young’s moduli in the , , and directions, respectively. (, ) are the corresponding Poisson’s ratios. , , and indicate the shear moduli. It should be noted that there are only nine independent material parameters for orthotropic materials, that is, , , , , , , , , and . The other Poisson’s ratios are defined as . By letting , and ), the present formulation can be readily applied to open shells and plates made from isotropic materials.

2.3. Energy Expressions

3D vibration analysis of structures has long been a goal of researchers in the field. Unfortunately, in most cases, it is very difficult to obtain the exact 3D analytical solutions of a structure. Therefore, this compels researchers to seek for approximate solutions by modern techniques such as the Rayleigh-Ritz method, DQ method, RBF method, FEM, and DSC approach. Among them, the Rayleigh-Ritz is widely adopted due to the reliability of its results and efficiency in practical calculation. Since the purpose of this work is to develop a unified, reliable, and efficient method for predicting 3D vibration characteristics of thick isotropic and orthotropic open shells and plates with arbitrary boundary conditions and variable geometry dimensions, the Rayleigh-Ritz procedure is applied in this paper. The deformation strain energy of a shell can be expressed in terms of the strains and stresses given previously asand the kinetic energy can be written aswhere is the density of the materials. In the present work, the general restrained conditions of an open shell are realized by introducing three groups of artificial linear springs (, , and ) at each boundary of the open shell to separately simulate the boundary forces and displacements. These boundary springs are assumed to be continuously distributed through the entire boundary domain. Therefore, arbitrary boundary conditions of the shell can be easily simulated by assigning these springs at proper stiffness. For example, the clamped boundary conditions can be realized by setting all the boundary springs to be infinitely rigid (represented by a very large number , where is the flexural stiffness , defined as ). For the purpose of convenience, symbols , , and (, , , and ) are used to indicate the stiffness (per unit area) of the boundary springs at the boundaries , , , and , respectively. Therefore, the deformation strain energy about the boundary springs () can be defined asIn particular, for open cylindrical shells, the energy functions are defined asand, for open spherical shells, the energy functions are defined aswhere and .

2.4. Admissible Displacement Functions

The Raleigh-Ritz method is one of the most common approximate methods used in the vibration analysis of continuous systems. The selection of suitable admissible displacement functions is of particular importance in this method due to the fact that its convergence and accuracy highly depend on the selected admissible displacement functions. It has been always of great interest for researchers to develop a unified, accurate, and feasible admissible displacement function which can be used to determine the vibration characteristics of a structure with arbitrary boundary conditions. In the present work, an improved Fourier series proposed by Li for the vibration analysis of 1D [42] and 2D [43] structures with general boundary conditions is extended to investigate the 3D vibration characteristics of open shells and plates. The displacements of an open shell/plate can be written asin which is the eigenfrequency and . , , and are the 3D improved Fourier series expansions constructed aswhere , , and . , , and represent the truncation numbers with respect to variables , , and . The magnitude of , , and will be determined by the convergence studies given in Section 3. The good accuracy and fast convergence behavior of the improved Fourier series solution will be shown later. , , and are the triple cosine Fourier series expansion coefficients for the displacement components , , and , respectively. , , , , , , , , and are the corresponding expansion coefficients of the auxiliary functions. Each coefficient is treated equally as an independent generalized coordinate and will be determined by performing the Rayleigh-Ritz procedure in future. , , and () are six auxiliary functions. The objective of introducing these auxiliary functions in each displacement is to eliminate any possible jumps with the original displacement function and its relevant derivatives at the boundaries when the displacement is expanded as a conventional triple cosine series. Thus, the function sets given in (13) are capable of representing any 3D motion of an open shell/plate with increasing accuracy as the truncated numbers are increased. It should be stressed that these jumps are not inherently related to the displacement function over the solution domain; instead, they are the artifact resulting from the Fourier series representation of the displacement solution. The detailed illustration is given in [43]. More information about the improved Fourier series can be seen in [42]. The auxiliary functions are given asIt can be verified thatSimilar conditions can be found in the auxiliary functions related to variables and . More detailed information about the auxiliary functions can be seen in [43, 44].

2.5. Solution Procedure

Since the improved Fourier series expression given in (9) is defined on the interval of , thus, linear transformations for coordinates from and (spherical shells) to and should be introduced for the practical programming and computing. With the energy expressions and admissible displacement functions given in the previous subsections, the remaining work is to obtain the actual displacement functions by determining the coefficients associated with the displacement functions. The Lagrangian energy functional comprising the strain and kinetic energies of the open shell/plate can be defined asSubstituting the energy expressions and the admissible displacement functions into the Lagrangian energy functional and performing the Rayleigh-Ritz operation with respect to the following undetermined coefficients:the vibration characteristic equation can be obtained and summed up as the following matrix form:where and are, respectively, the stiffness and mass matrices. Both of them are symmetric matrices and can be written asThe detailed expressions of submatrices , , , , , , , , and are given in Appendix B, respectively. is the column matrix which contains, in an appropriate order, the unknown expansion coefficients that appear in the series expansions (13); namely,where

The natural frequencies of open shells and plates under consideration can now be solved based on a standard numerical procedure. These frequencies are upper bounds on the exact values. With the truncation numbers in (13) being increased, the current results will approach the exact solution as closely as desired. The mode shape corresponding to each frequency can be easily determined by substituting the eigenvector back into the improved Fourier series expansions and solving for the ratios of coefficients. In addition, the proposed method can be readily applied to force vibration analysis of isotropic and orthotropic open shells and plates with arbitrary boundary conditions by summing the loading vector on the right side of (18) and solving the standing characteristic equation.

3. Numerical Examples and Discussions

According to the theoretical formulations presented in Section 2, a unified computer code based on the scientific software MATLAB has been developed, which is capable of calculating vibration results of isotropic and orthotropic plates and open shells with arbitrary boundary conditions and variable circumferential dimensions. With this code, in this section, several numerical examples are presented to test the convergence, accuracy, and reliability of the proposed improved Fourier series method. In the following calculations, three types of frequently encountered boundary conditions, that is, complete free boundary (), simply supported support (), and full clamped restraint (), are considered. Taking edge , for example, these types of boundary conditions can be written in terms of boundary spring rigidities as(1)free edge (): ,(2)simply supported edge (, ): ,(3)clamped edge (): ,where is the bending rigidity, defined as .

The appropriateness of defining these boundary conditions in terms of boundary spring rigidities will be verified by the numerical applications given later. The present formulation can be applied to shells subjected to uniform and nonuniform elastic boundary conditions, mixed boundary conditions, and their combinations as well. For simplicity, a short letter string is adopted to describe the boundary conditions of a plate/open shell, such as FCSF which denotes that the boundaries , , , and of the structure are completely free, clamped, simply supported, and completely free, respectively (counterclockwise).

3.1. Convergence Study

Theoretically, there are infinite terms in the improved Fourier expansions given in (13); thus, the expansions should be truncated and only finite terms are considered in the actual computing. As the first example, in Table 1, convergence studies of the first six frequencies (Hz) of completely free thick open cylindrical and spherical shells with different thickness-radius ratios are performed to check the convergence behavior of the improved Fourier solutions. Two thickness-radius ratios, that is, and 0.4, are considered in the study. The two shells are assumed to be made from isotropic materials with the following material properties:  GPa, , and  kg/m3. The geometric dimensions are  m,  m, and for the cylindrical shell and , ,  m, and for the spherical one. It can be seen from Table 1 that the frequencies have converged monotonically to four significant figures. For all modes, the differences between results forms “13 × 13 × 9” and “14 × 14 × 11” are very small, less than 0.05%. With the truncated numbers being increased, the current results will approach the exact solution as closely as desired; thus, the proposed improved Fourier method should be better understood as a method with arbitrary precision. In order to check the accuracy of the current method, results obtained by the FEM analysis (ANALYS, SOLID45, element size: 0.02 m) are also listed in the table. It can be found that the two results are in good agreement. Unless otherwise stated, the improved Fourier expansions are truncated as “14 × 14 × 11” in all the following examples.

3.2. Plates with Arbitrary Boundary Conditions

In this subsection, we consider vibrations of isotropic and orthotropic plates. Table 2 compares the lowest seven frequency parameters of isotropic plates with FFFF and SSSS boundary conditions. Two different thickness-length ratios, that is, and 0.5, corresponding to moderately thick and thick plates are considered in the comparison. The geometry parameters and material constants used in this example are ,  GPa, , and  kg/m3. The solutions given by Jin et al. [44] by the Ritz method with modified Fourier series and those of Liew et al. [45, 46] by the 3D Ritz method with general orthogonal polynomials using the Gram–Schmidt process are provided for a direct comparison. It is obvious that the present solutions agree very well with the referential data. The maximum differences between the present results and those of Jin et al. [44] and Liew et al. [45, 46] are less than 0.02% and 0.06%, respectively.

We next give the results for orthotropic plates with various boundary conditions and thickness-length ratios in Table 3 and Figure 2. In Table 3, the lowest five frequency parameters of an orthotropic rectangular plate with different combinations of boundary conditions are presented for three different thickness-length ratios; that is, , 0.2, and 0.3. The other geometrical parameters used for the analysis are . The plate is made of orthotropic material with the following material parameters: ,  GPa, , ,  GPa, and  kg/m3. From Table 3, we can observe that the frequency parameters of the plate decrease with the thickness-radius ratio increases. Figure 2 shows the first five mode shapes of the plate with FCFC boundary conditions. The figure reveals that the mode shapes of rectangular plates varied with thickness-length ratio.

3.3. Open Cylindrical Shells with Arbitrary Boundary Conditions

In this subsection, the accuracy and reliability of the current method are validated by comparing the results with those published in the open literature and the FEM analysis firstly. Then, numerous new results of open isotropic and orthotropic cylindrical shells with different combinations of boundary conditions are presented for various thickness-radius ratios, included angles, and length-radius ratios. Also, the effects of boundary conditions, geometric dimensions, and material properties on the natural frequencies of the cylindrical shells are studied.

As the first example, comparison of the first ten frequency parameters of a thick open cylindrical shell with CFFF, SSSS, and CCCC boundary conditions is presented in Table 4. The shell with length  m, width  m, thickness  m, and included angle is assumed to be made from isotropic materials with the following material properties:  GPa, , and  kg/m3. The comparison is performed between the current results and those provided by Liew et al. [23] by using the continuum and discrete approaches and the FE simulation (MSC/NASTRAN) based on the 3D shell theory. The symbol “” represents frequency parameters that were not considered in the referential work. It is observed that the present results match well with the approximate solutions and FEM results provided by Liew et al. [23]. For all the three cases, the differences between these three results are less than 0.67%. From this comparison, we can see that the presented method is accurate in predicting vibration results of open cylindrical shells with arbitrary boundary conditions. The table also reveals that it is appropriate to define the boundary conditions in terms of boundary springs.

Having validated the accuracy of the improved Fourier series solution, some further vibration results of isotropic open cylindrical shells with different combinations of boundary conditions for various geometric parameters are presented in Tables 57 and Figures 3-4. In Table 5, the first five frequency parameters of an isotropic deep open cylindrical shell subjected to as many as eleven combinations of boundary conditions are presented for various thickness-radius ratios. Five different thickness-radius ratios, that is, , 0.2, 0.3, 0.4, and 0.5, corresponding to moderately thick to very thick open cylindrical shells, are considered in the calculation. The geometrical and material parameters used for the analysis are  m, , ,  GPa, , and  kg/m3. From the table, we can see that the frequency parameters of the shell increase with the thickness-radius ratio increases. The 3D mode shapes for the shell with FCFC boundary conditions are given in Figure 3 as well. Each of these mode shapes is determined by substituting corresponding eigenvector back into the improved Fourier series expansions and solving for the ratios of coefficients then constructed by the SURF function based on MATLAB. The shear deformations in the thickness dimensions can be seen obviously, which may not be obtained by using the 2D shell theories. Unlike many existing approaches limited to cylindrical shells with small curvature (shallow shells), the presented method can be used to deal with deep open shells with variable circumferential dimensions. Therefore, the first five frequency parameters of an isotropic deep open cylindrical shell with eleven combinations of boundary conditions and five different included angles are presented in Table 6. The geometrical and material properties of the shell are the same as the previous one except that the thickness-radius ratio is and the circumferential dimensions under consideration are , 2π/3, π, 4π/3, and 5π/3. Frequency parameters of the shell reduce with increases in general. In order to enhance our understanding on the vibration deformation of open cylindrical shells, some selected mode shapes of the shell with CFCF boundary conditions are given in Figure 4. It is obvious that the mode shapes of the shell varied with circumferential dimensions. For the sake of completeness, the first five frequency parameters of a deep open isotropic cylindrical shell ( m, , ,  GPa, , and  kg/m3) with various boundary conditions and length-radius ratios are listed in Table 7. Three different length-radius ratios performed in the calculation are , 3, and 5. As seen from the table, the increase of the length-radius ratio results in the decreases of the frequency parameters.

To the best of the authors’ knowledge, there is no literature for 3D vibration analysis of orthotropic open cylindrical shells with arbitrary boundary conditions and variable circumferential dimensions. Therefore, some vibration results for the titled problem, which may be used as referential data for future researchers, are presented to fill this void. Unless otherwise stated, in all the following examples, open cylindrical shells under consideration are made from orthotropic materials with the following parameters:  GPa,  GPa, , , and  GPa,  kg/m3. In order to validate the accuracy of the current approach for the vibration analysis of orthotropic open cylindrical shells, Table 8 shows a comparison made for a thick deep open cylindrical shell ( m, , and ) with different sets of boundary conditions. The comparison is performed between the current 3D results and those obtained by ANSYS (SOLID 45, element size: 0.02 m). The lowest five frequencies (Hz) are considered in the comparison. From Table 8, it is obvious that the two results are in good agreement with each other. The proposed three-dimensional improved Fourier series method is applicable for the analysis of orthotropic open cylindrical shells with arbitrary boundary conditions.

Based on the above verification, free vibration results of open cylindrical shells with various boundary conditions and geometric parameters are provided in Tables 9-10 and Figures 5-6. Also, these results can be used to validate new 2D theories and new computational techniques in the future. In Table 9, the first five frequencies are listed for a thick shallow cylindrical panel ( m, , and ) with FFFF, FFFC, CFFF, FSFS, FCFC, SFSF, SSSS, SCSC, CFCF, CSCS, and CCCC boundary conditions and different thickness-radius ratios; that is, , 0.2, 0.3, 0.4, and 0.5. It can be seen from the table that frequencies of the shell are considerably increased by increasing the thickness-radius ratio. This is due to the fact that, with the increase of the thickness-radius ratio, the shell stiffness increases. Table 10 shows the first five frequencies (Hz) for a rare thick open cylindrical shell subjected to various boundary conditions and circumferential dimensions (, π/2, 3π/4, π, and 5π/4). The shell is made from the aforementioned orthotropic materials and with the following geometric parameters:  m, , and . From Table 10, we can see that the frequencies of the open cylindrical shell reduce by increasing the circumferential dimensions. This is due to the fact that, with the increase of the included angle, the shell stiffness reduces. Also, the lowest five mode shapes for the shell with full clamped boundary conditions are presented in Figure 5 to further enrich the vibration results of orthotropic cylindrical panels. As the last example, variation of the first ten frequencies (Hz) against material parameter (,  GPa, , , and  GPa,  kg/m3) for a deep open cylindrical shell ( m, , , and ) with FFFF, FFFC, CFFF, CCFF, SSSS, and CCCC is presented in Figure 6. Results of Figure 6 reveal that the effects of the material parameter on the vibration characteristics of open cylindrical shells varied with mode sequences and boundary conditions.

3.4. Open Spherical Shells with Arbitrary Boundary Conditions

In this subsection, the current method is applied to study the vibrations of isotropic and orthotropic open spherical shells with arbitrary boundary conditions. As the first example, accuracy of the presented method is validated by comparing the results with solutions published by other researchers as well as those obtained by the FEM analysis. In Table 11, comparison of frequency parameters is presented for a FFFF open spherical shell with different thickness-radius ratios. Two different thickness-radius ratios, that is, h/R = 0.05 and 0.25, are considered in the comparison. The comparison is carried out between the solutions of Reddy [12], Liew and Lim [13], and Reddy and Liu [14] and the 3D results of Liew et al. [21]. The shell parameters used in the comparison are  m, ,  GPa, and ,  kg/m3. From the table, we can see that the present results are in good agreement with the 3D solutions provided by Liew et al. [21]. However, results from different 2D theories deviate significantly as the thickness-radius ratio increases. The comparisons given in Tables 1 and 11 for deep and shallow open spherical shells with FFFF boundary conditions indicate that the current formulation is accurate on the vibration analysis of isotropic open spherical shells with different circumferential dimensions.

Having validated the accuracy of the present method, some further vibration results for isotropic open spherical shells with various boundary conditions, meridional dimensions, thickness-radius ratios, and included angles are presented in Tables 1214 and Figures 79. In Table 12, the first five frequency parameters of an isotropic deep open spherical shell with thickness-radius ratio , 0.2, 0.3, 0.4, and 0.5 are presented. The spherical panel is made from isotropic material ( GPa, , and  kg/m3) with top angle , end angle , included angle , and mean radius  m. Table 12 shows that any increase in the thickness-radius ratio , from 0.1 to 0.5, always leads to corresponding increases in the frequency parameters. Table 13 contains the frequency parameters of an open spherical shell with different boundary conditions and included angles. The shell is made from steel ( GPa, , and  kg/m3) with the following geometric parameters: , ,  m, and . From Table 13, it can be seen that frequency parameters of the shell reduce as the included angle increased. For the sake of completeness, an isotropic deep open spherical shell with different meridional dimensions (, ) is investigated, and the first five frequency parameters of the shell are given in Table 14. The material and geometric properties used in the investigation are ,  m, ,  GPa, and . In order to enhance our understanding of deep open spherical shell with different meridional dimensions, some selected mode shapes for the shell with CFCF boundary conditions are presented in Figure 7. Because of the curvature of the middle surface, the mode shapes for the doubly curved spherical shells are much more complicated than the singly curved cylindrical ones.

Then, free vibration analyses of deep orthotropic spherical shells with arbitrary boundary conditions are considered in the following discussions, for which very limited amount of results is available in the literature. In Table 15, the lowest five frequencies (Hz) of an orthotropic deep open spherical shell with eleven combinations of boundary conditions are presented for five different thickness-radius ratios; that is, , 0.2, 0.3, 0.4, and 0.5. The other geometrical parameters used for the analysis are , ,  m, and . In order to verify the accuracy of the present work, corresponding frequencies obtained by the FE analysis (ANSYS, SOLID 45, element size: 0.02 m) of the shell with thickness-radius ratio are included in the table as well. From Table 15, one can observe that the present results are in good agreement with the 3D elasticity solutions obtained by the FE analysis. The table also reveals that the frequency parameters of the shell increase with the thickness-radius ratio increases. Then, the first five frequencies for an open orthotropic spherical shell (, ,  m, and ) involving different boundary conditions and circumferential dimensions are given in Table 16. As expected, frequencies of the spherical shell vary from the maximum value for the shell to the minimum value for the one. In Table 17, frequencies of an open orthotropic spherical shell with various boundary conditions and meridional dimensions are presented. The shell with top angle , π/2, π/3, or π/6, end angle , included angle , mean radius  m, and thickness-radius ratio is considered in the analysis. Unlike cylindrical shells, with increasing the meridional dimensions of a spherical shell, the stiffness of the shell may not reduce. The vibration characteristics of thick open orthotropic spherical shells can be further elaborated by Figures 8-9 by considering the mode shapes of open spherical shells with FCFC boundary conditions and CCCC boundary conditions as their frequencies presented in Tables 15 and 16, respectively. Finally, variation of the first ten frequencies (Hz) against material parameter (,  GPa, , ,  GPa, and  kg/m3) of a deep open spherical shell (, ,  m, and , ) with FFFF, FFFC, CFFF, CCFF, SSSS, and CCCC is presented in Figure 10. By comparing Figures 10 and 6, we can see that the effects of the material parameter on the vibration characteristics of open shells varied with mode sequences, boundary conditions, and shell types.

4. Conclusions

This paper presents a unified method for the free vibration analysis of isotropic and orthotropic plates and shells with arbitrary boundary conditions and geometry dimensions in the framework of three-dimensional displacement-based energy formulation, including rectangular plates, open cylindrical shells, and the spherical ones. Regardless of boundary conditions, displacement in each structure direction is expanded as a triplicate product of the cosine Fourier series with the addition of certain supplementary terms introduced to eliminate any possible jumps with the original displacement function and its relevant derivatives at the boundaries. The boundary conditions of the shells and plates are accounted for by using the artificial spring technique and any classical boundary conditions can be easily realized by assigning the springs with corresponding values, without the need of making any modifications to the trial functions and solution program to adapt to different boundary cases. A unified computer code based on the scientific software MATLAB is developed to obtain the accurate 3D vibration results by using the Rayleigh-Ritz procedure. The validation of the current method, theoretical formulations, and program code is proved by comparing results with solutions in the open literature as well as data by the 3D FE analysis. It is seen that the current results are very close to the referential ones. The effects of boundary conditions, material parameters, and geometric dimensions on the frequencies are discussed as well. Finally, several 3D vibration results of isotropic and orthotropic plates and open cylindrical and spherical shells with different geometry dimensions are presented for various boundary conditions, which may serve as benchmark solutions for future researches as well as structure designers in this field.

Appendices

A. Strain-Displacement Relations for the Plates and Cylindrical and Spherical Shells

Based on the rectangular coordinate system given in Figure 1(b), the strain-displacement relations for plates areand, from Figure 1(c), the strain-displacement relations for open cylindrical shells are defined as Similarly, according to Figure 1(d), the strain-displacement relations for open spherical shells are given as

B. Stiffness and Mass Matrices for Cylindrical and Spherical Shells

The displacement components can be rewritten in the vector form aswhereThus, the stiffness and mass matrixes of an open cylindrical shell can be written as where , , and . Similarly, the stiffness and mass matrices of an open spherical shell are where , , and .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (nos. 51175098 and 51279035).