1. Introduction
The dispersion relation of a circular cylindrical shell had been a subject of great interest for several decades [1-3]. In some earlier works, numerical procedures were used exclusively in the computation of dispersion relations. Recently, Karczub [4] obtained an analytical expression for the dispersion relations from Flügge shell theory by using a symbolic algebra package, Mathematica. Karczub’s main concern was to check the agreement between his analytical results and results obtained previously from numerical methods, and limited to the harmonic orders
The solutions of the modal patterns for each of the propagating and non-propagating modes are also important and are not discussed in Karczub’s paper. These solutions are crucial to determine the vibration of a finite shell under various admissible boundary conditions and arbitrary external forces. The dispersion relations and the associated eigenvectors are also the means by which to construct transfer matrices for vibroacoustic transmission in cylindrical shell structures or pipe-hose systems [5]. The traditional and standard method of finding eigenvectors was given in detail by Leissa [6]; the eigenvector was normalized in such a way that its radial component of the displacement was unity, while the longitudinal and circumferential components were expressed as displacement ratios relative to the radial displacement. In some cases, these ratios could become exceedingly large and the eigenvector thus could deviate significantly from the eigenvector commonly used in mathematical physics.
In order to improve this problem for finding the normalized eigenvectors with norms equal to unity, an alternative method by using the built-in eigenvalue problem in any numerical package is proposed. A continuous variation of mode patterns for all kind of wave types in the desired frequency range will be rapidly obtained, and the physical meanings of wave propagation can also be shown more clearly and straightforward [8].
2. Methods to determine the propagating mode wavenumbers
2.1. Expression of equation of motion of cylindrical shell
The coordinate system which {
The equations of motion for the free vibration of a cylindrical shell can be written and shown as below
where
where
In order to transfer the expression from time and displacement domain into frequency and wavenumber domain, the following spectral representation can be further to use [7]
The dimensionless spectral listed in Eq.(1) can becomes as following
where Ω=(
2.2. Analytical expressions for dispersion relations
The relationship between s and Ω is so called dispersion relations of cylindrical shell. When a value of dimensionless angular frequency, Ω, is given, non-trivial solutions for the matrix form in the Eq.(4) exist at certain values of s. Finally, the specified roots can be integrated to the polynomial of s2. All the calculating processes and expressions can be rapidly built up by using any kind of symbolic algebra package, such as Matlab or Mathematica. According the physical characteristics of acoustic propagation, two conditions can further be separated and discussed as following.
2.2.1. n=0
When n=0, two independent sets of equations and corresponding characteristic equation can be calculated and shown as below:
For axial and radial displacements:
For the non-trivial solution for {U(s,0,Ω), W(s,0, Ω)}, the determinant of matrix in Eq. (5) should be zero ad obtained as below:
where the coefficients, gi, are the function of n and Ω, and shown as
The solutions of Eq.(5) represent a combination of longitudinally and radially expanding and contracting motions; these are the so-called breathing-modes. The expressions of polynomial s2 is obtained in Eq.(6).
Six characteristic roots are given
where
For the circumferential displacement:
The 2nd equation in Eq. (4) can be separated and shown as:
In the same way, under the non-trivial solution for V(s,0,Ω), the new relation can be rewritten as
Obviously, the solution of Eq.(8) represents a uniformly circumferential motion, which is therefore a torsional-mode. The corresponding characteristic roots may be considered the 4th branch pair and show as following
This torsional waves are known for being non-dispersive, since the wave speed is
where c
2.2.2. n≧1
Under the non-trivial solution of Eq. (4), a 4th polynomial of s2 can be obtained as
These roots yield eight branches of dispersion curves, which are grouped into four branch pairs:
where g
2.3. Methods to determine the propagating mode pattern
As mentioned earlier, the solutions regarding the modal patterns (eigenvectors) for all four root branch pairs are important. The conventional method has been frequently used for a long time. In this method, the mode shapes are expressed as the ratios among U, V and W. By returning to Eq.(4) after each of the roots (±s
A new alternative method, which solves for the eigenvectors by using a commercially available linear algebra package, will be introduced. This method is extremely simple and straightforward. The kernel of processing is by using the built-in formulation of eigenvalue and eigenvector problem in the numerical package.
According to the eigenvalue problem, the eigenvectors can be found by solving the following homogeneous equations, and the simple flow can be introduce by following:
Step 1. For a specific frequency, dimensionless frequency, Ω, can be obtained and rewritten as Ω
where the asterisk in superscript of L
Step 2. After slightly rearranging, Eq.(13) can be further modified as
Step 3. If Ωs2 is assumed to be unknown for the time being, Eq.(14) can be rewritten to be a typical linear eigenvalue problem and shown as:
In most cases in Eq. (15), for any given value of s (the analytical root), there are three distinct eigenvalues in Ω2 by using the built-in eigenvector function in the numerical package, but only one of them will be equal to Ωs2. Consequently, the eigenvector solution of Eq.(13) is one of the three eigenvectors of Eq.(15), which corresponding eigenvalue is equal to Ωs2. This is because, among the three eigenvalues, only the one that is equal to Ωs will yield the specific s given here, the other two will yield different values of s and are not what we want. When s is specified, the corresponding eigenvector can further calculated by using the numerical solution package, such as MatLab.
2.4. Comparsion between alternative and conventional method for modal pattern
The dispersion relations were calculated using the thickness-radius factor β=0.0058 and the Poisson ratio v=0.3. For n=2, the dispersion relation and corresponding modal patterns of Branch pairs 1 are obtained in Fig. 2.(a)and 2(b). The consistency of modal patterns between the alternative and conventional method can also be demonstrated in Fig. 2.(c). The physical characteristics and variations of cylindrical shell at specific circumferential number can be obtained clearly and straightforward.
3. Result and discussion
In order to discuss the characteristics of wave propagation in dispersion relations and modal pattern, the numerical results are calculated by assuming the dimensionless thickness parameter, β, to be equal to 0.0058, and the Poisson ratio equals to 0.3. According to Karczub and Leissa, for the usual range of parameters and n≧1, the roots of Eq.(11) were found to have the four following types:
where b
The dispersion curves for all the roots and the corresponding eigenvectors plotted in the figures shown hereafter are those with positive real parts only. In the case of Type 2 roots which have no real part, only the positive imaginary part will be plotted. The roots with negative real parts are not plotted; their eigenvectors are basically the same as those of positive real roots, except there will be sign changes in the longitudinal components if the signs of the radial and circumferential components are chosen to be the same as those plotted.
3.1. The case when n=0
The dispersion relations and the associated eigenvectors are shown in Fig. 3 and 4, respectively. For simplicity, the eigenvectors [U' V' W'] are plotted by using the absolute values of the displacement components, i.e., [|U| |V| |W|]'. The dispersion relations for the longitudinal and bending waves of a plate, and those for shear waves are also shown for reference. It shows that at frequencies much lower than the ring frequency (Ω<<1), Branch 1 represents a non-dispersive longitudinal wave, from which it gradually deviates as Ω approaches the neighborhood of unity; Branch 1’s wave eventually approaches that of a plate bending wave when Ω>>1. These are illustrated by the mode shapes shown in Fig. 4(a). Below ring frequency, see the mode shape in Fig. 4(b), Branch 2 (Branch pair 2 root with positive real part) represents the Type 4 root, but it becomes a Type 1 root after crossing the ring frequency and represents the longitudinal wave above the ring frequency, say Ω>1. Therefore, when Ω>1, the Branch 2 curve takes over the role of representing the longitudinal wave which was previously the role of Branch 1’s curve. Below the ring frequency, Branch 3 (Branch pair 3 root with positive real part) represents the Type 3 root, but above the ring frequency, it represents the Type 2 root. The modal components of the Branch 3 roots are predominantly radial (see Fig. 4(c)). Branch 4 (Branch pair 4 root with positive real part) represents the characteristic roots of Eq. (13), which therefore represents the non-dispersive torsional wave and has nearly the same dispersion relation as shear waves. The mode shape of torsional mode is not shown; it is simply [U V W]’=[0 1 0]’.
3.2. The case when n=1
Fig. 5 and 6 show the dispersion relations and the associated eigenvectors, respectively, for the n=1 circular harmonics. The equivalent curves for Bernoulli beam (labeled as “Berbeam”), Timoshenko beam (labeled as “Timbeam”), and the longitudinal and shear wave curves of a plate are also shown for reference. Branch 1 has only real roots (Type 1 roots); the dispersion curve representing beam bending (when Ω<0.5) agrees well with that of the Timoshenko beam. The Bernoulli beam, on the other hand, agrees with the shell and Timoshenko beam curves only up to Ω=0.03; its applicable frequency range is therefore limited. At frequencies above Ω=0.5, the originally beam-like bending wave gradually transitions to a plate-like bending wave with a predominantly radial modal displacement.
Below Ω=0.5, the Branch 2 curve is purely imaginary, and represents the Type 2 root. When Ω>0.5, the roots become positive real (Type 1) and gradually approach the shear wave curve, with a predominantly circumferential displacement component when Ω>1.5 (see Fig. 6(b)). Obviously, this branch has a cut off frequency which is about Ω=0.5.
Below the ring frequency, the Branch 3 curve indicates a Type 3 root with a slow and exponentially diminishing wave. When Ω>1.4 it becomes a Type 1 root, and the Branch 3 curve approaches that of longitudinal wave (see Fig. 5 and 6(c)) at further higher frequencies. The propagating waves in Branch 3 thus has a cutoff frequency of about Ω=1.4. The Branch 4 curve indicates a Type 4 root below the ring frequency, but above ring frequency it has only the imaginary part (Type 2) and represents the transversal near field distortion (see Fig. 6(d)).
3.3. The case when n=2
The dispersion relations for the n=2 circular harmonic case are shown in Fig. 7. At very low frequencies, say Ω<0.016, the roots in Branch 1 are the Type 3 roots. At higher frequencies, the roots become Type 1 and represent propagating waves (which thus have a cut off frequency at Ω≒0.016). Fig. 8(a) shows that below ring frequency, Ω<1, the modal patterns of Branch 1 roots have three well-coupled displacement components, although the predominant component is radial. When Ω>1, the modal displacement becomes predominantly radial and the propagating speed asymptotically approaches that of a plate.
At very low frequencies, say Ω<0.016, the roots in Branch 2 are the Type 4 roots. Between say Ω=0.016 and the ring frequency, the roots of this branch are Type 2; they represent the near field distortion of the three coupled components. When Ω >1.15, they become Type 1 propagating waves with the wave speeds that asymptotically approach that of the shear wave, and their modal patterns reflect predominantly circumferential shear motions. The curve of Branch 3 indicates a Type 3 root when 0<Ω<1.15, a Type 2 root when 1.15<Ω<2.2, and a Type 1 root when Ω>2.2. This is a propagating wave with a high cutoff frequency (at Ω≒2.2), and at higher frequencies its wave speed approaches asymptotically to that of the longitudinal wave. Branch 4 represents the Type 4 root when 0<Ω≦1.15, which, however, turns into Type 2 and represents the near field distortion when Ω>1.15.
4. Conclusion
The classic problem of the dispersion of waves in a cylindrical shell had been re-examined with analytical solutions obtained by using a symbolic algebra package. The previous work by Karczub did not include the analytical solutions for the dispersion relations of axisymmetric waves (of circular harmonic order n=0). An axisymmetric wave contains both longitudinal and transverse components; the former are important in the transmission of longitudinal vibration, and the latter are particularly important in acoustics due to their higher acoustic radiation efficiency than the corresponding waves of higher circular harmonic orders (n>0). In this way, a complete set of analytical solutions which based on Flügge shell theory is now available for all orders of circular harmonics, n=0, 1, 2, …, ∞.
A considerable effort has been expended on the solutions for the modal patterns of all propagating and non-propagating modes, because a complete set of properly normalized eigenvectors is crucial to the solution of the vibration problem of a finite shell under various admissible boundary conditions. The eigenvectors obtained by the conventional method are not conveniently normalized as are those commonly used in mathematical physics. A new alternative method to find normalized eigenvectors with norms equal to unity has been proposed and discussed.
Use of analytical solutions has demonstrated the capability of a straightforward continuous tracking on the frequency-dependent changes of the root types and of the corresponding modal patterns for each branch of the dispersion curves associated with a particular analytical root. Therefore, a parallel display of the dispersion curves and of the associated modal patterns used in this paper has provided more insight about the wave phenomena in a cylindrical shell.
Appendix
The coefficients in Eq. (11) are listed as follows:
where
The coefficients in Eq. (12) are given in the following.
where
and
References
- 1.
Forsberg K. 1963 Influence of boundary conditions on the modal characteristics of thin cylindrical shells. Journal of the Acoustical Society of America,119 6 1898 1899 1729-8806 - 2.
Fuller C. R. 1981 The effects of wall discontinuities on the propagation of flexural waves in cylindrical shells. Journal of Sound and Vibration,75 207 228 0002-2460 X. - 3.
Hwang Y. F. 2002 Structural-acoustic wave transmission and reflection in a hose- pipe system, Proceedings of 2002 ASME International Mechanical Engineering Congress and Exposition, IMECE2002-32683,91 96 0-79183-643-6 - 4.
Karczub D. G. 2006 Expressions for direct evaluation of wave number in cylindrical shell vibration studies using the Flügge equations of motion. Journal of the Acoustical Society of America,119 6 3553 3557 0001-4966 - 5.
Lesmez M. W. Wiggert D. C. Hatfield F. J. 1990 Modal analysis of vibrations in liquid-filled piping systems. Journal of Fluids Engineering ASME,112 3 311 318 0098-2202 - 6.
Leissa A. W. 1997 Vibration of Shells, Acoustical Society of America,1-56396-293-4 - 7.
Skelton E. A. James J. H. 1997 Theoretical Acoustics of Underwater Structures, Imperial College Press,1-86094-085-4 - 8.
Liu Y. C. Hwang Y. F. Huang J. H. 2009 Modes of wave propagation and dispersion relations in a cylindrical shell. ASME Journal of Vibration and Acoustic,131 4 1048-9002