Abstract

In order to improve the spherical thin shells’ vibrations analysis, we introduce a new analytical method. In this method, we take into consideration the terms of the inertial couples in the stress couples’ differential equations of motion. These inertial couples are omitted in the theories provided by Naghdi–Kalnins and Kunieda. The results show that the current method can solve the axisymmetric vibrations’ equations of elastic thin spherical shells. In this paper, we focus on verifying the current method, particularly for free vibrations with free edge and clamped edge boundary conditions. To check the validity and accuracy of the current analytical method, the natural frequencies determined by this method are compared with those available in the literature and those obtained by a finite element calculation.

1. Introduction

With the fast and continuous advancement of technology, it becomes necessary to have more efficient and effective innovative devices and engineering structures. Many publications have focused on modeling lightweight structures [13]. Thin spherical shells are used in civil engineering and a variety of engineering applications such as submarine hulls, aircraft, smart devices, and health monitoring systems [47]. It is well known that the spherical shells can be used in complex environmental conditions and can be exposed to a variety of dynamic excitations inducing excessive vibrations [8]. Undesirable vibrations cause structural failure, cracks in parts of machine elements, malfunction of electronic devices, and many other issues [9]. By measuring structure or system frequencies, designers can select speeds remote from resonance and can prevent many adverse effects of vibration [9]. Determining the frequency spectrum of a structure is an essential step of its design. In addition, developing an accurate and unified analytical solution for axisymmetric thin spherical shells with free edge or clamped edge boundary conditions remains an extremely challenging task. This solution is useful to provide reference data for future research in related fields. Hence, a survey of the literature depicts the existence of many studies on free vibrations analysis [1014]. Besides, the research for curved shells, especially for cylindrical and spherical shells, have attracted much attention recently [7, 1520].

Shell structures' vibrations analysis can be carried out by analytical or numerical methods and can be tackled by different theories and approximations. The linear theories for plates and shells are as follows:(i)Love–Kirchhoff’s theory that does not take account of transverse shear deformation that generalizes Euler–Bernoulli’s theory for beams [2123].(ii)Mindlin–Reissner’s theory for plates and Naghdi–Reissner’s one for shells, taking into account the transverse shear deformation extending Timoshenko’s theory for beams [2430].

They all assume an elastic behavior. They are said to be of the first-order: the displacement fields vary linearly with variable t (Figure 1) in the thickness, while the thickness does not vary.

However, higher-order theories have been developed like the theory HSDTs (Reddy) and many others [31].

Exact (Naghdi–Kalnins, Kunieda, current method) or approximate (Kunieda) analytical resolution methods are based on Legendre functions [3234]. One cites among the numerical methods, finite differences, weighted residuals, variational methods (Rayleigh–Ritz), and finite element method [3538].

Using the finite element method, El Harif studied the transverse shear effect in plates and shells of revolution [39]. Batoz et al. focused on structure modeling [40]. Juber investigated the analytical solution for the dynamic analysis of thin composite cylindrical and spherical shells. He presented the fundamental natural frequencies and mode shapes for simply supported composite cylindrical and spherical shells, by applying the third-order theory of Reddy [41]. Moreover, Sammoura et al. developed an analytical solution for curved piezoelectric micromachined ultrasonic transducers by applying Naghdi–Kalnins’ theory [17].

The present paper concentrate on a new analytical method inspired by the Naghdi–Kalnins’ theory [32] and the model presented in [39, 40]. This method solves elastic thin spherical shells’ axisymmetric vibrations’ differential equations. For this purpose, we introduce the differential equations of motion in terms of stress, where we take into consideration the terms of the inertial couples in the stress couples’ equations. These terms are omitted in the analytical methods proposed by Naghdi–Kalnins and Kunieda [32, 33]. Then, we introduce a new stress function different from Naghdi–Kalnins’ one and derivate new equations depending on the displacement and the stress function. Moreover, by using the Legendre functions, we analytically solve these differential equations in the specific case of the hemispherical shell’s harmonic axisymmetric vibrations, under the free and clamped edges boundary conditions. To check the validity and accuracy of the current method, we compare the current analytical results with Naghdi–Kalnins’ results [32] and Kunieda’s results [33]. It is worth to mention that Naghdi–Kalnins’ theory is one of the most important references used in different recent studies [17, 42, 43]. After that, we compare the results with the finite element method (FEM) ones. The FEM numerical simulation results confirm the correctness and accuracy of the current analytical method.

The remainder of this paper is organized as follows: in Section 2, we introduce the analytical model formulation. In Section 3, we propose the general solution for harmonic axisymmetric vibration. After that, we present and discuss the results in Section 4. The conclusions of this study are described in Section 5.

2. Analytical Model Formulation

2.1. Basic Concept

Inspired by Naghdi–Kalnins’ theory of axisymmetric vibrations of the elastic spherical shell, we introduce here new basic concepts that are used to analyze the axisymmetric vibrations response of the elastic spherical shells as illustrated in Figure 2. The main purpose of the approaches presented below is to derivate the expression of the differential equations by applying the motion equations with their stress couples equations. These stress couples equations take into account the terms of the inertial couples. These terms are omitted by Naghdi–Kalnins. In addition, we propose a new stress function that helps to derivate the new equations, different from Naghdi–Kalnins’ equations. On the other side, as the thickness of the shells is small in comparison with the radius , we use Love’s first-order approximation for the dynamic analysis of the thin spherical shell. We assume that the transverse shear and normal strains are negligible. We recall here that the basic equations characterize the elastic deformation of the spherical shell in the weightlessness and the spherical coordinate system .

2.2. Strain-Displacement Relations in Spherical Coordinates

By applying Love-Kirchhoff’s assumption, the strain can be expressed as a function of the displacement vector components and their derivatives in the spherical coordinate system, where , and are the displacement vector components in , and directions, respectively.

Furthermore, the total strain of a spherical shell can be expressed into two components: membrane strains and flexural strains as mentioned in [44]:where .

As shown in Figure 1, is the radial position measured from the middle surface along the spherical shell thickness, wherewhere M is a point on the middle surface.

The membrane and the flexural strains are given by the following expressions [44]:

2.3. Stress Resultants and Stress Couples in Spherical Coordinate System

The tangential stresses of the spherical shells in spherical coordinates can be expressed aswhere is Young’s modulus, is Poisson’s ratio.

As shown in Figure 3, it is convenient to introduce the stress resultants , transverse shear stress resultants , and the stress couples , where , by integrating the infinitesimal stresses along the thickness of the spherical shell as [44].

Moreover, the stress resultants can be written as

The stress couples can be expressed bywhere is the flexural rigidity defined as [17]

2.4. Stress Differential Equations of Motion

The stress differential equations of motion in the weightlessness can be expressed as [39, 40]where is the mass density of the material and is the received pressure perpendicular to the spherical shell.

2.5. Derivation of Differential Equations for Axisymmetric Vibrations

According to the torsionless axisymmetry with respect to the z-axis, we have

The compatibility relation can be deduced from (4), (5), and (12) by eliminating and replacing the strains by their expressions according to the stress resultants:

The expression of the shear stress resultant as a function of displacements and can be obtained from (15) and (18):

Inspired by [32], we introduce a new stress function through the following new relations:where, as proposed in [32], we have

The tangential displacement can be related to the radial displacement , by elimination of from (15) as

By using (4), (5), (15), (22), (23), and (24), we get the following differential equation:where is expressed as

By substituting the equations (23) and (24) in the compatibility equation (21), another differential equation related also to and can be defined as

The solution of (27) and (29) can be written as , where the coefficients A and B can be set to zero without loss in generality. Then, the first differential equation in and can be defined as

To obtain the second differential equation in and , the shear stress resultant (22) into (17) and using (23) and (24) can be substituted:

3. General Solution for Harmonic Axisymmetric Vibrations

The steady-state problem of the spherical shells can be represented by the system of the two differential equations (30) and (31). Thus, the solution of this system in and can be written in the following form:where and are the magnitudes of the radial displacement and the stress function and is the circular frequency.

The natural frequency can be noted as

By eliminating from (30) and (31) and using the notation of the natural frequency, a differential equation in of the sixth-order can be obtained aswhere

In the particular case of the free vibration analysis, the surface loads are equal to zero. Equation (34) is then written as follows:

Therefore, the general solution of the free vibration differential equation (36) can be expressed in terms of Legendre functions aswhere each radial displacement can be written aswhere and are the Legendre functions of the first and the second kind of order , respectively. can be defined as

With the aid of the identities,

The characteristic equation of the differential free vibrations equation of spherical shells (34) is

4. Results and Discussion

In this section, we test the validity and accuracy of the method exposed previously by determining the natural frequencies of spherical shells. Then, to demonstrate the correctness of the aforementioned formulation, two appropriate boundary conditions are treated: free edge and clamped edge. The current results are compared, on the one hand, with those available in literature, that is, the analytical results obtained by Naghdi–Kalnins [32] and by Kunieda [33]. It should be noted that these two analytical methods do not take into consideration the terms of the inertial couples in their equations of motion. On the other hand, the current results are also compared with the numerical results obtained by two different axisymmetric finite elements: the first one is SHELL61 that does not take into account transverse shear deformation and the second is SHELL209 that takes it into account. Furthermore, these two shell finite elements are also suitable for the modeling of axisymmetric thin shells. SHELL61 has two nodes and four degrees of freedom per node: displacements in the x (radial), y (axial), and z (circular) directions and rotation around the z-axis. The loading may be axisymmetric or nonaxisymmetric [45, 46]. SHELL209 has three nodes and three degrees of freedom per node: displacements in the x, y directions and rotation about the z-axis [45]. More information on the shape functions is available in [45, 46]. The FEM solutions are computed via ANSYS Mechanical APDL. The mesh is fine enough to guarantee the accuracy of the numerical results. In this study, only the thin shells whose thickness is less than the twentieth of the radius of curvature will be considered. The thickness to radius ratio in this work varies in the range 0.0025 to .

4.1. Free Vibration with Free Edge
4.1.1. Frequency Equation

In this section, we focus on solving the problem of axisymmetric vibrations of a thin hemispherical shell with a free edge with a free circular boundary at .

In the case of thin hemispherical shells , the Legendre functions of the second kind can tend towards infinity for . To eliminate the singularity of , the coefficients are set to zero, where is as (38).

In the free edge boundary condition at , the stress resultant , the stress couple , and the transverse shear stress resultant are set to zero :whereand the coefficients of the frequency equationThen, they can be written aswhere the coefficient was obtained by the linear combination between the two equations (42) and (43).

The Legendre function of a complex degree can be calculated by the following integral [28]:

To calculate the lowest natural frequency from equation (46), for and , we specify the value of and and then evaluate the determinant by incrementing . This procedure is repeated until the value of the determinant changes sign. Table 1 and Figure 4 present an example of this procedure, where . As shown in Table 1, the real parts of the natural frequencies are negligible. For this reason, in this case, we consider just the imaginary part. Then, the small box of Figure 4 shows that there are two sign changes between 0.82 and . Therefore, it means that in this interval, we have two natural frequencies.

It is worthy to mention that a determinant is a complex number, as shown in Table 2. In this case, the real parts of the determinants are not negligible. For this reason, we used the determinant’s modulus to determine the lowest natural frequencies, which correspond to the minimum values closest to zero. The variation of the determinant’s modulus as a function of natural frequency for and is also plotted in Figure 5. As seen in the zoomed panel, the results indicate the existence of five minimum values, which correspond to the lowest natural frequencies. Their values are presented in Table 3.

4.1.2. Comparison of Results

In order to verify the accuracy of the current solution, the results are compared with those obtained by Naghdi and Kalnins [32], Kunieda [33], and calculus by the finite elements SHELL61 and SHELL209. The difference estimators are defined as follows:where is the frequency obtained by the Kunieda method, is the frequency obtained by the current method, is the frequency obtained by the Naghdi–Kalnins method, is the frequency obtained by the FE-SHELL61, and is the frequency obtained by the FE-SHELL209. It should be noted that the analytical results obtained by the Kunieda method are taken as a reference.

Table 4 shows a comparison study of the lowest natural frequency for spherical elastic thin shells under free vibration with free edge boundary condition for and , calculated by the current method, the Naghdi–Kalnins method [32], and the Kunieda method [33]. We observe that the current results and Kunieda’s ones are identical. It is also seen that the Naghdi–Kalnins results are very close to these results, with .

A comparative study of the first seventeen natural frequencies between the analytical solutions (the current solution and Kunieda’s [33] solution) and numerical solutions (FE-SHELL61 and FE-SHELL209), is carried out in Table 5 for thin spherical shells under free vibration with the free edge boundary condition at and .

From Table 5, we observe a very good agreement between the different results confirming the accuracy of the current analytical method. On the one hand, we see that the current results are close to the results of Kunieda’s method. The difference between the results of these two analytical methods is very small and does not exceed . This closeness is verified even for higher vibration modes. On the other hand, we see also that the FE-SHELL61 results are very close to the SHELL209 results particularly for the first modes. Indeed, when the ratio is very small (here ) the transverse shear becomes negligible and the two finite elements give almost the same results. We observe also that, in order, the current results are closer, first, to the results of the finite element SHELL61, then to those of Kunieda, and finally to those of the finite element SHELL209.

The results of Table 5 are instructive in more than one way. The current method reveals a new natural frequency (6th mode) which does not appear in Kunieda’s method and finite element calculation while the frequency (13th mode) obtained by the finite element SHELL61 is nonexistent for the two analytical calculations and the finite element SHELL209 calculation.

4.2. Effect of the Thickness to Radius Ratio on the Lowest Natural Frequency

The results in Table 6, which are plotted in Figure 6, show the effect of the ratio on the lowest natural frequency of thin spherical shells under free vibrations with the free edge. We observe that the lowest frequency, calculated by the current method and FE methods, is almost independent of the ratio, while that calculated by the Naghdi–Kalnins method is sensitive to this ratio. This strong variation in frequency, calculated by the Naghdi–Kalnins method, as a function of the ratio is not legitimate.

All the methods give almost the same lowest frequency when the ratio is very small.

We also observe that the two finite element methods give exactly the same lowest frequency even for high values of the ratio . Which leads us to say that taking into account the transverse shear does not have effect on the lowest frequency and therefore that the current method can be applied in this case even for high ratios.

4.3. Free Vibrations with Clamped Edge Boundary Condition
4.3.1. Frequency Equation

For free vibrations with clamped boundary conditions, the edge of the shell cannot translate in the and directions and cannot undergo any rotation around the axis. The boundary conditions, in this case, can be summarized aswhere . Also, in this case of clamped boundary conditions, the Legendre functions of the second kind can tend towards infinity for . To eliminate the singularity of , the coefficients are set to zero as mentioned in the free vibrations with free edge condition section.

Using the boundary conditions (50) to (52) along with the expression for , the following form is obtained:

We note that the displacement can be expressed as a function of (as derived in the appendix).

The coefficients of the frequency equation:and they can be written as

To calculate lowest natural frequency from equation (54)for and , we specify the value of and and then evaluate the determinant by incrementing as done in the section of frequency equation for free vibrations with free edge.

Table 7 presents the lowest natural frequencies for thin spherical shells at and . The results are obtained by the current analytical method under axisymmetric free vibrations with clamped edge boundary condition.

4.3.2. Comparison of Results

In order to verify the accuracy of the current method, in the case of axisymmetric free vibrations for hemispherical thin shells with clamped edge boundary condition at and , the frequencies of the first sixteen modes are compared with those obtained by Kunieda [33] and calculus by the finite elements SHELL61 and SHELL209, as shown in Table 8.

The current solutions are in close agreement with the solutions obtained by Kunieda, FE-SHELL61, and FE-SHELL209. The difference between the results obtained by the current method and Kunieda and between the current method and FE-SHELL61 are closely zero. For the higher modes, a negligible difference appears in the results of FE-SHELL209 compared to the others results, probably due to the fact that it takes into account the transverse shear.

The results of Table 8 are instructive in more than one way. The current method reveals a new natural frequency (5th mode) which does not appear in Kunieda’s method and in the the finite element calculation while the frequency 1.240 obtained by the FE-SHELL61 is nonexistent for the three other methods.

4.3.3. Influence of the Thickness to Radius Ratio

The results in Table 9 show the effect of the ratio on the four (or five) lowest natural frequencies of thin spherical shells under free vibrations with clamped edge.

We observe that, for ratio , varying between and , the natural frequencies obtained by the current method are very close to those obtained by Kunieda’s method and those of the calculation by finite elements.

The effect of transverse shear is nonexistent even for the highest ratio : the element SHELL209 which takes account of the transverse shear gives the same results. Therefore, the current method can be applied in this case even for high ratios and even for high modes, with negligible error as observed above.

The effect of the ratio is weak on the lowest frequency. The variation of this frequency is when the ratio varies from the value 0.0025 to while the influence of the ratio on the frequency of mode 4 is more important. The variation of this frequency is , when the ratio varies in the same interval. Kunieda [33] considers that the effect of is not substantial for thin shells.

This effect is much greater for thin shells with a clamped edge than for those with a free edge. When the ratio varies from the value to 0.02, the percentage variation of the lowest frequency is for clamped edge and only for free edge.

The current method reveals a new natural frequency (1.067, mode 4) which is not obtained by the finite element calculus and which is absent from Kunieda’s results.

5. Conclusions

The present method solves the differential equations of the axisymmetric vibrations of elastic thin spherical shells. These differential equations of motion in terms of stress take into account the terms of inertial couples in the equations of stress couples. These terms are omitted in the analytical methods proposed by Naghdi–Kalnins and by Kunieda. A new stress function is introduced and new equations as a function of the radial displacement and the stress function are derived. The harmonic axisymmetric free vibrations are described by a sixth-order differential equation in terms of the radial displacement , the natural frequency , the thickness to radius ratio of the shell , and Poisson’s ratio . The general solution is expressed by the Legendre function of degree . Taking into account the boundary conditions makes it possible to reach the analytical solution in the specific case of the axisymmetric harmonic vibrations of the hemispherical shell. In this work, we tested two boundary conditions: free and clamped edges.

The results demonstrate the correctness and high accuracy of the current analytical method compared with the Naghdi–Kalnins and Kunieda analytical methods and the finite element calculation.

This method is distinguished by taking into account the inertia of rotation (inertial couple in the stress couples’ differential equation of motion) in addition to the two inertias of translation while the methods of Naghdi–Kalnins and Kunieda only take into account the two inertias of translation, which probably explains the evidence by this method of new frequencies. This method can be applied, with good approximation, for the calculation of the first natural frequencies even for moderately high values of ratios for free edge or clamped edge cases, despite the fact that it does not take into account the effect of the transverse shear.

It can also be extremely useful for the validation of numerical methods.

Appendix

A: Proof of the Tangential Displacement

Substituting the stress resultant expressions (23) and (24) in the expression (26) of the tangential displacement , then, by using the first differential equation in terms of F and as in (30), the tangential displacement can be written as

Using equations (30), (31), and (37), the stress function can be expressed in terms of the radial displacement as

From (A.2) and (38), it can be then found that

After substituting (A.3) in (A.1), the tangential displacement can be expressed as

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors would like to thank Mr. Vishal Gupta for his technical support.