Research Article  Open Access
Free Vibration of Layered Circular Cylindrical Shells of Variable Thickness Using Spline Function Approximation
Abstract
Free vibration of layered circular cylindrical shells of variable thickness is studied using spline function approximation by applying a point collocation method. The shell is made up of uniform layers of isotropic or specially orthotropic materials. The equations of motions in longitudinal, circumferential and transverse displacement components, are derived using extension of Love's first approximation theory. The coupled differential equations are solved using Bickleytype splines of suitable order, which are cubic and quintic, by applying the point collocation method. This results in the generalized eigenvalue problem by combining the suitable boundary conditions. The effect of frequency parameters and the corresponding mode shapes of vibration are studied with different thickness variation coefficients, and other parameters. The thickness variations are assumed to be linear, exponential, and sinusoidal along the axial direction. The results are given graphically and comparisons are made with those results obtained using finite element method.
1. Introduction
Circular cylindrical shells are used in various fields like aviation, missiles, ship buildings, and chemical industries. Shells made of composite materials with variable thickness are used increasingly, since composite structures are having high specific stiffness, better damping, and shock absorbing characters over the homogeneous ones. The study of vibrational behavior of such shells is very important. The effect of variation of thickness on frequency parameter of the shell, which is made up of different layered materials, has been studied by very few researchers. Baker and Herrmann [1] analysed three layered (Sandwich) shells, including the effects of shear deformation, rotary inertia, and initial stress. Sivadas and Ganesan [2] studied the vibration of circular cylindrical shells having the thickness variations of linear and quadratic along the axial direction. A series of studies has been made on vibration of cylindrical shells by Tonin and Bies [3], Takahashi et al. [4], Suzuki et al. [5] and Sivadas and Ganesan [6]. Hinton et al. [7] presented free vibration analysis of variable thickness of plates and curved shells using a finite strip formulation. The fundamental frequencies of laminated anisotropic circular cylindrical shells are presented by Sun et al. [8] using finite element method (FEM). Zhang [9] used a propagation approach to analyse the crossply laminated composite cylindrical shells. Hufenbach et al. [10] presented a study on vibration and damping behaviour of multilayered composite cylindrical shells using analytic calculation method. Sakiyama et al. [11] and Tsuiji and Sueoka [12] analysed the vibration of cylindrical panel using the RaleighRitz method. Recently, Tizzi [13] applied the Ritz procedure for optimization of cylindrical shell profile under a frequency constraint, and Toorani and Lakis [14] studied the vibrations of nonuniform composite shells applying the combination of hybrid finite element analysis and shearable shell theory.
Mizusawa and Kito [15] applied the spline strip method to study the vibration of crossply laminated cylindrical panels. This method involved expressing displacement functions in a strip element as the product of basic function series in the axial direction and Bspline functions in the circumferential direction. However, there seems to be no work carried out so far on vibration of symmetric angleply layered cylindrical shells with variable thickness using Bickley spline function, which is done in the present study.
The present work analyses the flexural free vibration of layered circular cylindrical shells of variable thickness. The equations of motion are derived using Love’s first approximation theory for homogeneous shells. The layers are considered to be thin, elastic, specially orthotropic, or isotropic and assumed to be perfectly bonded together and move without interface slip. Three different thickness variations (linear, exponential, and sinusoidal) are considered along the axial direction of the cylinder. The governing coupled differential equations are obtained in terms of the reference surface displacements which are in longitudinal, circumferential, and transverse directions. Assuming the displacement functions in a separable form, they reduce to a system of ordinary differential equations on a set of displacement functions which are functions of meridional coordinate only. Two sets of boundary conditions are imposed and two types of materials are used to analyse the problem. In general, the equations have no closed form solution, so that the numerical solution techniques have to be resorted to.
The spline function technique is adopted to solve the coupled differential equations which are in three displacement functions. Bickley [16] successfully tested the spline collocation method over a twopoint boundary value problem with cubic spline. Viswanathan and Navaneethakrishnan [17] and Viswanathan and Kim [18] have also demonstrated this, along with its attractive features of elegance in handling and convergence. Recently, Viswanathan et al. [19] studied the vibration of crossply cylindrical shell walls including shear deformation theory using the spline function techniques. The advantage of this method is that a chain of lowerorder approximations than the global higher order approximation.
The three displacement functions are approximated using cubic and quintic splines. Collocation with these splines yields a set of field equations which, along with the equations of boundary conditions, reduce to a system of homogeneous simultaneous algebraic equations on assumed spline coefficients which results in a generalized eigenvalue problem. This eigenvalue problem is solved using eigensolution technique to obtain as many frequencies as required, starting from the least. From the eigenvectors, the spline coefficients can be found to construct the mode shapes.
2. Formulation of the Problem
The system of differential equations in terms of longitudinal, circumferential, and transverse displacements components is derived, which characterise the vibration of a thin shell of revolution. The general line of procedure of Ambartsumyan [20] for the classical theory of thin shell is adopted. The development is based on the Love’s first approximation theory in which the rotatory inertia and transverse shear deformation are neglected. Such an approach results in an analytically simpler procedure, by way of less number of equations of motion and avoidance of nonlinear terms, thereby conveniencing the application of spline function method. The coordinate system and the geometric parameters of the laminated cylindrical shells of constant thickness are shown in Figure 1.
In general, the thickness of the th layer of the shell is assumed in the form where is a constant thickness and Here is the length of the cylinder, , and are the coefficients of linear, exponential, and sinusoidal variations, respectively.
The thickness of the layers is not completely independent. Their dependence is given by where is the mass density of the th layer and is the distance of the outer boundary of the th layer from the reference surface. This may be interpreted as determining one of the in terms of the rest of . If the shell wall has only two layers, in particular, one can obtain Here . It may be noted that for linear and sinusoidal variation of thickness .
The stress resultants and moment resultants are expressed in terms of the longitudinal, circumferential, and transverse displacements , , and of the reference surface. The displacements are assumed in a separable form given by where and are the longitudinal and rotational coordinates, is the time, is the angular frequency of vibration and is the circumferential node number. Using (2.5) in the constitutive equations and the resulting expressions for the stress resultants and the moment resultants in the equilibrium equations, the governing differential equations of motion are obtained in the form The operators are defined in Appendix A.
3. Method of Solution
The differential equations on the displacement functions of (2.6) contain derivatives of third order in , second order in , and fourth order in . Therefore, the present form is not suitable to the solution procedure we propose to adopt. Hence, the equations are combined within themselves and a modified set of equations is derived. The modified equations now become as 2nd order in , 2nd order in , and 4th order in , and is given by The new operators and are given in Appendix B.
The parameters are nondimensionalised as Here is the radius of the cylinder and is the total thickness of the shell. Also we define and , since we consider only two layers.
The displacement functions , and are approximated by the cubic and quintic spline functions , and as stated below The boundary conditions are used as follows: (i) both the edges are clamped (C–C) and (ii) both the edges are hinged (H–H). The resulting field and boundary conditions give rise to the generalized eigenvalue problem of the form where and are matrices of order is a matrix of order , and is the number of knots of the splines on axial direction. The parameter is the eigenparameter and the eigenvector whose elements are the spline coefficients. Only twolayer shells are considered with ratio of thickness of the first mentioned layer to the total thickness, at the one edge of the cylinder.
4. Results and Discussion
Convergence study is made for the frequency parameter value to fix the number of knots of the spline function. As mentioned earlier, only twolayered shells are considered and it is tested for High strength graphite (HSG) and Sglass epoxy (SGE) material combinations. Table 1 shows the material properties of High strength graphite (HSG) and Sglass epoxy (SGE).

After a number of trials, it is found that the number of knots could be taken as 14, since for the next value of the percent change in values of are very low, the maximum being 0.35%. The results are not furnished here due to space constraints.
Comparative studies are next made for homogeneous cylindrical shells of exponential variation in thickness. Table 2 presents the natural frequencies obtained for the shell with C–C boundary conditions and compared with the results obtained by Takahashi et al. [4] and Sivadas and Ganesan [2]. The present frequency parameter converted into the parameter that has been given by Sivadas and Ganesan [2], where . The percentage changes between the present results and Sivadas and Ganesan [2] for = 0.4, 0.6, 0.8, and 1.0 are 9.8%, 8.4%, 6.3%, and 4%, respectively. It shows that, when the value of increases the difference in percentage decreases. This may be due to the method we adopted. In this paper the value of is taken as 1 in all the cases since the difference in percentage is minimum for = 1 when compared with results obtained by Sivadas and Ganesan [2]. This indicates that the correctness of the analysis and accuracy of the results by using spline function techniques.
In this work, asymmetric free vibration of layered circular cylindrical shells of variable thickness is studied. Only twolayered materials with HSG and SGE combinations are used in this analysis, and the first three meridional modes are considered in all the analyses that follow.
In Figure 2, the variation of frequency parameter with respect to the increase of the relative thickness ratio for the layered cylindrical shells under linear variation in thickness is displayed. In this case, it can be assumed that , where is the taper ratio . The values of the ratio of the shell’s constant thickness to radius and the ratio of the shell length to the radius are fixed as 0.02 and 1.5, respectively. The value of taper ratio is fixed as . The two layers of the shell are arranged in the order of HSG and SGE materials. Thus, when , the inner layer disappears, and the shell is homogeneous, which is made of SGE material. When , the outer layer disappears, again the shell is homogeneous, made of HSG material. Figures 2(a) and 2(b) correspond to the node number , with C–C and H–H boundary conditions, respectively. It is clearly seen that as increase, decreases for a for all values of and decrease for 0.8 and a small increase for 0.8. If the order of the materials are inverted (SGEHSG) then the frequency parameter increases as increase for all .
(a)
(b)
The results are not shown here for want of space. For the extreme values of , equal to 0 or 1, the shell becomes homogeneous, with the material of either of the two layers. It is seen that it is possible to attain a desired frequency, between these two extreme values by suitably choosing the value of . This is interesting from a design point of view.
Figure 3 presents the variation of frequency parameters with relative thickness of layer for HSGSGE materials by fixing and under C–C and H–H boundary conditions. The thickness varies exponentially and coefficient is fixed as 0.2. Figures 3(a) and 3(b) correspond to the node number , with C–C and H–H boundary conditions, respectively. In Figure 4, the nature of the frequency parameter for sinusoidal variation in thickness is depicted. The other parameters and are fixed with . The effect of frequency parameters presented in Figures 3 and 4 almost has the similar pattern as discussed in Figure 2 (linear variation).
(a)
(b)
(a)
(b)
Figures 5(a)–5(f) show the manner of variation of the frequency parameter with reference to the circumferential node number . The range of is considered between 0 and 10. A shell of HSGSGE lamination under C–C and H–H boundary conditions is considered with , and . Figures 5(a)–5(c) show the effect of on for C–C conditions and Figures 5(d)–5(f) show the effect of on for H–H conditions. All the three types of variation in thickness of layers are considered, as indicated in the diagrams. It is seen that all the frequency parameter values decrease upto and then increase. The curvature at the turning points seems to be greater for lower modes. The absolute and relative differences between the maximum and minimum values of , caused in the range of values of considered, is more in the case of C–C boundary conditions than with that of H–H boundary conditions. The kind of thickness variation in layers does not seem to greatly affect the nature of variation of with .
(a)
(b)
(c)
(d)
(e)
(f)
The frequency parameter is explicitly a function of the length of the cylinder. Hence, when studying the influence of the length of the cylinder on its vibrational behaviour, the actual frequency , and not , is considered. Figures 6(a)–6(c) depict the manner of variation of the actual frequency (in 10^{3} Hz) with respect to the length parameter for HSGSGE layered cylinders under C–C boundary conditions with and . All the three types of variation in thickness of layers are considered. As increases, is observed to decrease, in general. The decrease is fast for very short shells (for herein), the rate of decrease increasing with higher modes and then the decrease is very low for . The percent changes in over the range of for three cases (linear, exponential and sinusoidal variations) depicted are, respectively: (a) 452.034%, (b) 450.476%, and (c) 453.573% for C–C conditions. Similar phenomenon is observed in the case of layered cylindrical shells of variable thickness under H–H boundary conditions shown in Figure 7. The percent changes in over the range of considered, for three different variations, are (a) 272.268%, (b) 270.533%, and (c) 273.858% for H–H conditions.
(a)
(b)
(c)
(a)
(b)
(c)
In Figure 8, the influence of the nature of variation of thickness of the layers of the shell on its vibrational behaviour is studied. A HSGSGE shell held under C–C boundary conditions with the three types of variation in thickness of layers is considered, with , and . Figure 8(a) relates to linear variation in thickness of layers. The thickness is constant when the taper ratio . Variation of with respect to for is studied. It is seen that is almost constant for all the values of . The effect of exponential variation in thickness of layers is analysed in Figure 8(b). When , the thickness is uniform. The thickness at the end of the cylinder is higher or lower than the thickness at the other end according to . The effect of sinusoidal variation in thickness of layers on frequency parameters is studied in Figure 8(c). These effects are almost similar to those due to the exponential variation discussed above. Here, the coefficient of thickness variation is considered over the range . The thickness of the shell is the same at 0 and ; the surface of the shell is convex or concave for .
(a)
(b)
(c)
In Figure 9, the influence of the taper ratio , the coefficient of exponential variation of thickness and the coefficient of sinusoidal variation on are depicted, along with the effect of the H–H boundary conditions. The effect of is almost same for all the cases of linear and exponential variation, as described in Figure 8. In this variation, the C–C boundary conditions contribute slightly higher values to the influence of the coefficients of thickness variation on frequencies than the H–H conditions contributing values to the influence of the coefficients of thickness variation on frequencies.
(a)
(b)
(c)
5. Conclusion
The influence of the natural frequencies of the vibration of layered cylindrical shells of variable thickness has been analysed. The materials of the layers, length of the shell, and coefficients of variable thickness affect the frequency. A desired frequency of vibration may be obtained by a proper choice of the relative thickness of the layers, length parameter, and the coefficient of thickness variations. The clampedclamped (C–C) boundary conditions give rise to higher frequencies in comparison with hingedhinged (H–H) boundary conditions. The nature of variation in thickness of layers considerably affects the natural frequencies. When the circumferential node number is increased, the frequencies initially decrease and then increase. The effect of increasing the length of the cylinder is a decrease in frequencies, for all kinds of variation in thickness of layers. This study also shows the elegance and usefulness of the spline functions with application of the collocation method for boundary value problems.
Appendices
A. The operators are:
The differential operators appearing in (2.6) are where and , and are the elastic coefficients of constant thickness, which are extensional rigidities, the bendingstretching coupling rigidities and the bending rigidities, respectively.
B. The operators are:
The differential operators and appearing in (3.1) are
Acknowledgment
This work was supported by Inha University Research grant.
References
 E. H. Baker and G. Herrmann, “Vibrations of orthotropic cylindrical sandwich shells under initial stress,” AIAA Journal, vol. 29, pp. 963–975, 1966. View at: Google Scholar  Zentralblatt MATH
 K. R. Sivadas and N. Ganesan, “Free vibration of circular cylindrical shells with axially varying thickness,” Journal of Sound and Vibration, vol. 147, no. 1, pp. 73–85, 1991. View at: Publisher Site  Google Scholar
 R. F. Tonin and D. A. Bies, “Free vibration of circular cylinders of variable thickness,” Journal of Sound and Vibration, vol. 62, no. 2, pp. 165–180, 1979. View at: Publisher Site  Google Scholar
 S. Takahashi, K. Suzuki, T. Kosawada, and E. Anzai, “Vibration of cylindrical shells with variable thickness,” Bulletin of Japan Society of Mechanical Engineers, vol. 24, pp. 1826–1836, 1981. View at: Google Scholar
 K. Suzuki, E. Anzai, and S. Takahashi, “Vibration of cylindrical shells with varying thickness,” Bulletin of Japan Society of Mechanical Engineers, vol. 25, pp. 1108–1119, 1982. View at: Google Scholar
 K. R. Sivadas and N. Ganesan, “Axisymmetric vibration analysis of thick cylindrical shell with variable thickness,” Journal of Sound and Vibration, vol. 160, no. 3, pp. 387–400, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Hinton, M. Özakça, and N. V. R. Rao, “Free vibration analysis and shape optimization of variable thickness plates, prismatic folded plates and curved shells. Part 1: finite strip formulation,” Journal of Sound and Vibration, vol. 181, no. 4, pp. 553–566, 1995. View at: Publisher Site  Google Scholar
 G. Sun, P. N. Bennett, and F. W. Williams, “An investigation on fundamental frequencies of laminated circular cylinders given by shear deformable finite elements,” Journal of Sound and Vibration, vol. 205, no. 3, pp. 265–273, 1997. View at: Publisher Site  Google Scholar
 X. M. Zhang, “Vibration analysis of crossply laminated composite cylindrical shells using the wave propagation approach,” Applied Acoustics, vol. 62, no. 11, pp. 1221–1228, 2001. View at: Publisher Site  Google Scholar
 W. Hufenbach, C. Holste, and L. Kroll, “Vibration and damping behaviour of multilayered composite cylindrical shells,” Composite Structures, vol. 58, no. 1, pp. 165–174, 2002. View at: Publisher Site  Google Scholar
 T. Sakiyama, X. X. Hu, H. Matsuda, and C. Morita, “Vibration of twisted and curved cylindrical panels with variable thickness,” Journal of Sound and Vibration, vol. 254, no. 3, pp. 481–502, 2002. View at: Publisher Site  Google Scholar
 T. Tsuiji and T. Sueoka, “Free vibrations of twisted thin cylindrical panels (numerical analysis by using the RayleighRitz method),” Transactions of the Japan Society of Mechanical Engineers, Part C, vol. 55, no. 514, pp. 1325–1329, 1989. View at: Google Scholar
 S. Tizzi, “A Ritz procedure for optimisation of cylindrical shells, formed by a nearly symmetric and balanced angleply composite laminate, with fixed minimum frequency,” Computers & Structures, vol. 84, no. 3132, pp. 2159–2173, 2006. View at: Publisher Site  Google Scholar
 M. H. Toorani and A. A. Lakis, “Free vibrations of nonuniform composite cylindrical shells,” Nuclear Engineering and Design, vol. 236, no. 17, pp. 1748–1758, 2006. View at: Publisher Site  Google Scholar
 T. Mizusawa and H. Kito, “Vibration of crossply laminated cylindrical panels by the spline strip method,” Computers and Structures, vol. 57, no. 2, pp. 253–265, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 W. G. Bickley, “Piecewise cubic interpolation and twopoint boundary problems,” Computer Journal, vol. 11, pp. 206–208, 1968. View at: Google Scholar  Zentralblatt MATH
 K. K. Viswanathan and P. V. Navaneethakrishnan, “Free vibration study of layered cylindrical shells by collocation with splines,” Journal of Sound and Vibration, vol. 260, no. 5, pp. 807–827, 2003. View at: Publisher Site  Google Scholar
 K. K. Viswanathan and K. S. Kim, “Free vibration of antisymmetric angleplylaminated plates including transverse shear deformation: spline method,” International Journal of Mechanical Sciences, vol. 50, no. 1011, pp. 1476–1485, 2008. View at: Publisher Site  Google Scholar
 K. K. Viswanathan, K. S. Kim, J. H. Lee, H. S. Koh, and J. B. Lee, “Free vibration of multilayered circular cylindrical shell with crossply walls, including shear deformation by using spline function method,” Journal of Mechanical Science and Technology, vol. 22, no. 11, pp. 2062–2075, 2008. View at: Publisher Site  Google Scholar
 S. A. Ambartsumyan, “Theory of anisotropic shells,” Tech. Rep. NASA TTF118, 1964. View at: Google Scholar
Copyright
Copyright © 2010 K. K. Viswanathan 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.