A New Method Based on Laplace Transform and Its Application to Stability of Pipe Conveying Fluid

Wen, H. B.; Yang, Y. R.; Li, P.; Li, Y. D.; Huang, Y.

doi:https://doi.org/10.1155/2017/1472601

Shock and Vibration

On this page

Abstract Introduction Discussion Conclusions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2017 | Article ID 1472601 | https://doi.org/10.1155/2017/1472601

A New Method Based on Laplace Transform and Its Application to Stability of Pipe Conveying Fluid

H. B. Wen,^1,2Y. R. Yang,²P. Li,²Y. D. Li,^1,3and Y. Huang^2,4

Academic Editor: Lorenzo Dozio

Received31 Oct 2016

Accepted26 Dec 2016

Published24 Jan 2017

Abstract

A new differential transformation method is developed in this paper and is applied for free vibration problem of pipes conveying fluid. The natural frequencies, critical flow velocities, and vibration mode functions of such pipes with several typical boundary conditions are obtained and compared with the results predicted by Galerkin method and finite element method (FEM) and with other results archived. The results show that the present method is of high precision and can serve as an analytical method for the vibration of pipes conveying fluid.

1. Introduction

The dynamics of pipes conveying fluid receives considerable attention always because it contains abundant mechanical phenomena. The equations of motion of the fluid-pipe system are high-order partial differential equations, so the analysis is difficult. Many analytical methods have been proposed to study this dynamic model since Ashley and Haviland [1] began to study the dynamics of straight tubes conveying fluid. These methods include the finite element method (FEM), analytical method, multiple scales method, Galerkin method, transfer matrix method, differential quadrature method, and generalized differential quadrature rule.

Gregory and Paidoussis [2] studied unstable oscillation of tubular cantilevers conveying fluid by an analytical method and Bubnov-Galerkin method. Both parameter and combination resonances for a continuously flexible pipe were studied by Paidoussis and Sundararajan [3], where Bolotin’s method and a numerical Floquet analysis were presented. Hannoyer and Paidoussis [4] studied the dynamics and stability of cylindrical tubular beams conveying fluid and simultaneously subjected to axial external flow, where the solutions were obtained by means of a method similar to Galerkin method, with eigenfunctions approximated by Fourier series. Misra et al. [5, 6] have studied the dynamics and stability of the cantilevered pipes conveying fluid by finite element method; the conclusions were the same as those of Gregory and Paidoussis [2]. Olson and Jamison [7] introduced a general-purpose finite element program to simulate the dynamics of elastic pipes conveying fluid with different boundary conditions. Kuiper and Metrikine [8] employed a D-decomposition method for stability analysis of a riser conveying oil from seabed to a floating platform. Wang and Ni [9] established the nonlinear vibrations of a cantilevered pipe conveying fluid with motion-limiting nonlinear constraints by extending the differential quadrature method (DQM). Utilizing the differential quadrature method again, Qian et al. [10] analyzed the instability of the pinned-pinned pipe conveying fluid under thermal loads. Tornabene et al. [11] extracted the critical flow velocity of straight pipes conveying fluids for various fluid densities by the generalized differential quadrature (GDQ) method. Ni et al. [12] analyzed the free vibration problem of pipes conveying fluid with several typical boundary conditions by a relatively new semianalytical method called differential transformation method (DTM). Sinir and Demir [13] used the method of multiple scales to investigate the vibrations of highly tensioned pipe subject to external force. Quite recently, Li et al. [14] have developed the steady-state Green’s functions for forced vibrations of Timoshenko beam with damping effects. Subsequently, Li and Yang [15] employed the same method to study the forced vibration of pipes conveying fluid. Three examples of the pipe conveying fluid were analyzed, the results of which show the validity of the Green function method.

In present paper, a new differential transformation method is developed to analyze the free vibration of pipes conveying fluid. Starting from the basic equations, the variable separation technique and Laplace transforms are employed to derive vibration modes of the pipe conveying fluid, where all the constants involved are determined by the boundary conditions. The vibration mode is a closed and explicit form. Furthermore, a linear system with four equations is obtained by the application of boundary conditions. For the free vibration problem of pipes conveying fluid, the determinant of coefficient matrix of the equations is zero. So the natural frequency equations are derived, which is an implicit form in terms of . The natural frequency equation is solved by a computer program. In general, the frequency of the fluid-pipe system will be complex and the dynamic system will be stable or unstable according to the real part or the imaginary part of the frequency being positive or negative.

2. The Equation of Small Motions

The equation of small motion of pipes conveying fluid is [16]where is flexural rigidity; is the fluid flow velocity; and are the mass-per-unit length of fluid and the pipe, respectively; is the structure damping effect; is the transverse deflection of the pipe; is the coordinate along the center line of the pipe; and is the pipe length.

For the sake of illustration, the following dimensionless quantities are introduced:

Substitution of these terms into equation (1) yields the dimensionless governing equation of motionThe corresponding dimensionless boundary conditions are as follows:(i)At a pinned end(ii)At a clamped end(iii)At a free end

3. Vibration Mode Equations and Frequency Equations of the Pipe Conveying Fluid

Let us consider the motions of the pipe for which the dimensionless lateral deflection is given bywhere is a dimensionless circular frequency defined by being the circular frequency of motion. In general, will be a complex number and the system will be stable or unstable depending on the imaginary component of , , being positive or negative; in the case of neutral stability, is wholly real and is the mode function.

Substituting equation (7) into equation (3) yieldswhere

To derive the mode function of the pipe conveying fluid, the Laplace transform method with respect to the variable is applied in (9):where the parameter in the transformed domain is a complex variable; , , , and are constants which can be determined by the boundary conditions of the pipes.

To obtain the inverse transform of , let us assume thatThus, the following results can be obtained by the method in [14, 17]:where are functions specified byThe inverse transform of (11), which is the mode function of the pipe conveying fluid, can be obtained as follows:Inserting (13a)–(13d) into (15) yieldswhere are defined byIf the right of equation (3) is the Dirac delta function, the Green function of equation (3) will be derived by the same method as above. Green’s function has high precision and computational efficiency in the forced vibration analysis of beam [14, 15].

To derive the frequency equation of the pipe and to determine the constants , , , and , it is necessary to calculate the various order derivatives of . From (16), we can obtainwhere the various order derivatives of are as follows:

From (18a) to (18d), the intrinsic relations between the physical quantities at the boundary and can be obtained

From (4) to (7), we can determine , , , and at and as follows:(i)At a pinned end(ii)At a clamped end(iii)At a free end

For simplicity, the boundary conditions as (21) are considered. Inserting (21) into (20) yieldsIn order that a nontrivial solution of these equations exists, the determinant of the coefficient matrix of (24) must be zero. Thus, we can obtain the frequency equation of a pinned-pinned pipe conveying fluid as follows:Similarly, the frequency equations of a clamped-clamped pipe conveying fluid and a cantilevered pipe conveying fluid can be obtained as follows:From (24), we can obtain the ratio of and as follows:Utilizing (16) and (26), one can obtain the mode function equation of a pinned-pinned pipe conveying fluid as follows:Similarly, the mode function equations of a clamped-clamped pipe conveying fluid and a cantilevered pipe conveying fluid can be obtained as follows:

4. Numerical Results and Discussion

4.1. The Stability of Pipe

The natural frequencies of pipe conveying fluid can be obtained from (25a) to (25c). It is impossible to obtain solutions of (25a)–(25c) by direct methods because the roots of the problem cannot be expressed in simple explicit form in terms of . So a trial and error procedure is adopted to seek the roots of (25a)–(25c). In this case, keeping and fixed for each set of calculations, the calculation of the real and imaginary part of the complex circular frequency is required for a series of values of . For most of the calculation, a computer program is used along the following lines [18]:(i)Starting with and trial values of and corresponding to , a minimizing procedure finds the appropriate values of and which meet (25a)–(25c).(ii)Taking an increment of in and using the values of and found above as first approximations, the minimizing procedure finds the complex frequency for .

Thus the computer works continuously and the values of and are obtained for increasing values of , for any mode function of the system and any given and . Some results are shown in Figures 1–3, in which the natural frequencies correspond to the three cases of support conditions specified in (21)–(23), respectively.

Figure 1 shows the real and imaginary components of the dimensionless frequencies of the lowest three modes of a pinned-pinned pipe as functions of the dimensionless flow velocity. Based on the calculated results of Galerkin method and the present method, the first mode frequency vanishes altogether at , indicating the onset of first mode divergence; beyond the critical , the first mode frequency becomes wholly imaginary, indicating the system buckles in the first mode; if the flow velocity is increased further, the second mode frequency vanishes altogether at , indicating the onset of second mode divergence; beyond the critical of , the frequencies of the first and second modes coalesce and frequencies become complex conjugates, indicating the system loses stability by couple-mode flutter. Based on the calculated results of FEM, the frequency as function of the flow velocity is the same to the calculated results of other two methods except and ; the real component of the first frequency is 0.0246 at ; the real component of the second frequency is 0.6326 at . So the FEM has a certain error at some flow velocities, such as and , in comparison with Galerkin method and the present method.

Figure 2 shows that the real and imaginary components of the dimensionless frequencies of the lowest three modes of a clamped-clamped pipe as functions of the dimensionless flow velocity. Based on the calculated results of the present method, the first mode frequency vanishes altogether at , indicating the onset of first mode divergence; beyond the critical , the first mode frequency becomes wholly imaginary, indicating the system buckles in the first mode; if the flow velocity is increased further, the system loses stability by couple-mode flutter in the first mode and second mode at . The calculated results of the three methods are different at ; the real components of the first frequency are 0.9994 and 0.2111, respectively, with Galerkin method truncating four modes and FEM. According to Paidoussis and Sundararajan [3], the first mode frequency vanishes altogether. Therefore, Galerkin method and FEM have a certain error at some flow velocities, such as . As the flow velocity increases, the difference is more and more, between the calculated results of Galerkin method and the calculated results of other methods. The reason is that the number of the modes truncated is insufficient, according to the literature [19]. In this paper, the lower four modes are used in Galerkin method to solve the equation of small motion of pipes conveying fluid. The first and second modes of the clamped-clamped pipe conveying fluid represent the fifth and the higher modes of a clamped-clamped beam, when dimensionless flow velocity is relatively larger. So, the calculated results of Galerkin method are the difference from those of the other two methods.

Figure 3 shows the real and imaginary components of the dimensionless frequencies of the lowest three modes of a clamped-free pipe as functions of the dimensionless flow velocity. Based on the calculated results, beyond , the imaginary part of the second frequency becomes negative, and the real part of the second frequency is still positive, indicating the system flutters in the second mode through a Hopf bifurcation. Galerkin method has a certain error in comparison with FEM and the present, since the number of the modes truncated is insufficient.

4.2. The Mode Function of Pipe

From (27a) to (27c), one can find that the mode function of the pipe system is the function of , , , and . These variables are not independent because they have an implicit relationship as (25a)–(25c). So, in generally, it is impossible to obtain the mode functions.

When flow velocity is zero, (3) can be degenerated into the free vibration equation of tubular beam. For the free vibration of tubular beam, the natural frequencies are wholly real, and the vibration mode functions can be obtained in analytical form. The vibration mode functions of tubular beam with several typical boundary conditions are obtained at as follows:(i)A pinned-pinned tubular beam(ii)A clamped-clamped tubular beam(iii)A clamped-free tubular beamwhere , corresponding to the various mode functions of the system. These results agree with literatures [20, 21].

The modal displacement patterns of pipe conveying fluid can be obtained with (27a)–(27c), and the modal displacement patterns of a clamped-clamped pipe and a clamped-free pipe can be seen in Figures 4 and 5, respectively. When flow velocity is zero, the modal displacement patterns of pipe are classical normal modes, as seen in Figures 4(a) and 5(a). When is not zero, the equation of motion of pipe contains the Coriolis term, which renders classical normal modes impossible, as seen in Figures 4(b) and 5(b). The motions associated with the symmetric modes of a clamped-clamped pipe will produce the Coriolis force, which will induce the motion of antisymmetric modes. The Coriolis force, produced by the motions associated with the modes of a clamped-free pipe, is still in the direction against the motion along the whole pipe. This will tend to the oscillations. In general, the effect of the Coriolis force is more and more as the flow velocity increased.

(a)

(b)

(a)

(b)

5. Conclusions

The paper presented a new differential transformation method to study the free vibration of the pipe conveying fluid with various boundary conditions. The mode function equations and the frequency equations of the pipe conveying fluid were obtained by the present method. The frequency equation of pipe conveying fluid is implicit form in terms of . A trial and error computer program was written for the of the pipe conveying fluid, so the solving procedure is simple and rapid. The vibration mode equation has an explicit form in terms of , so the vibration mode function of the pipe conveying fluid will be calculated when the of pipe conveying fluid is known. By comparing with the calculated results of Galerkin method and FEM, this analytical method is valid for the problem of dynamics and stability of pipes conveying fluid.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 11102170); the project of Science & Technology Department of Sichuan Province, China (Grant no. 2015GZ0037); the Opening Project of Sichuan Province University Key Laboratory of Bridge Non-Destruction Detecting and Engineering Computing, China (Grant no. 2016QZJ03).

References

H. Ashley and G. Haviland, “Bending vibrations of a pipe line containing flowing fluid,” Journal of Applied Mechanics, vol. 17, pp. 229–232, 1950.
View at: Google Scholar | MathSciNet
R. W. Gregory and M. P. Paidoussis, “Unstable oscillation of tubular cantilevers conveying fluid. I. theory,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 293, no. 1435, pp. 512–527, 1966.
View at: Publisher Site | Google Scholar
M. P. Paidoussis and C. Sundararajan, “Parametric and combination resonances of a pipe conveying pulsating fluid,” Journal of Applied Mechanics, Transactions ASME, vol. 42, no. 4, pp. 780–784, 1975.
View at: Publisher Site | Google Scholar
M. J. Hannoyer and M. P. Paidoussis, “Instabilities of tubular beams simultaneously subjected to internal and external axial flows,” Journal of Mechanical Design, Transaction of ASME, vol. 100, no. 2, pp. 328–336, 1978.
View at: Publisher Site | Google Scholar
K. S. Van, Dynamics and stability of curved pipes conveying fluid [M.S. thesis], MoGill University, 1986.
A. K. Misra, M. P. Païdoussis, and K. S. Van, “On the dynamics of curved pipes transporting fluid. Part I: inextensible theory,” Journal of Fluids and Structures, vol. 2, no. 3, pp. 221–244, 1988.
View at: Publisher Site | Google Scholar
L. G. Olson and D. Jamison, “Application of a general purpose finite element method to elastic pipes conveying fluid,” Journal of Fluids and Structures, vol. 11, no. 2, pp. 207–222, 1997.
View at: Publisher Site | Google Scholar
G. L. Kuiper and A. V. Metrikine, “On stability of a clamped-pinned pipe conveying fluid,” Heron, vol. 49, no. 3, pp. 211–232, 2004.
View at: Google Scholar
L. Wang and Q. Ni, “The nonlinear dynamic vibrations of a restrained pipe conveying fluid by Differential Quadrature Method,” Journal of Dynamics and Control, vol. 2, no. 4, pp. 56–61, 2004.
View at: Google Scholar
Q. Qian, L. Wang, and Q. Ni, “Instability of simply supported pipes conveying fluid under thermal loads,” Mechanics Research Communications, vol. 36, no. 3, pp. 413–417, 2009.
View at: Publisher Site | Google Scholar | Zentralblatt MATH
F. Tornabene, A. Marzani, and E. I. Elishakoff, “Critical flow speeds of pipes conveying fluid using the generalized differential quadrature method,” Advances in Theoretical and Applied Mechanics, vol. 3, pp. 121–138, 2010.
View at: Google Scholar
Q. Ni, Z. L. Zhang, and L. Wang, “Application of the differential transformation method to vibration analysis of pipes conveying fluid,” Applied Mathematics and Computation, vol. 217, no. 16, pp. 7028–7038, 2011.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
B. G. Sinir and D. D. Demir, “The analysis of nonlinear vibrations of a pipe conveying an ideal fluid,” European Journal of Mechanics. B. Fluids, vol. 52, pp. 38–44, 2015.
View at: Publisher Site | Google Scholar | MathSciNet
X. Y. Li, X. Zhao, and Y. H. Li, “Green's functions of the forced vibration of Timoshenko beams with damping effect,” Journal of Sound and Vibration, vol. 333, no. 6, pp. 1781–1795, 2014.
View at: Publisher Site | Google Scholar
Y.-D. Li and Y.-R. Yang, “Forced vibration of pipe conveying fluid by the Green function method,” Archive of Applied Mechanics, vol. 84, no. 12, pp. 1811–1823, 2014.
View at: Publisher Site | Google Scholar
F. J. Bourrieres, Sur un Phenomene d'Oscillations Auto-Entrenue en Mecanique des Fluides Reels, vol. 147 of Publications Scientifiques et Techniques du Ministere de l'Air, 1939.
V. A. Ditkin and A. P. Prudnikov, Fromulaire pour le calcul Operationnel, Traduction du Russe par G. Carvallo, Masson & Cie, Paris, France, 1967.
View at: MathSciNet
R. W. Gregory and M. P. Paidoussis, “Unstable oscillation of tubular cantilevers conveying fluid. I. Theory,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 293, no. 1435, pp. 512–527, 1966.
View at: Publisher Site | Google Scholar
H.-H. Qi and J. Xu, “Effect of Galerkin modal truncation on natural frequency analysis of a cantilevered pipe conveying fluid,” Journal of Vibration and Shock, vol. 30, no. 1, pp. 148–151, 2011.
View at: Google Scholar
R. E. Bishop and D. C. Johnson, The Mechanics of Vibration, Cambridge University Press, New York, NY, USA, 1960.
View at: MathSciNet
M. P. Paidoussis, Fluid-Structure Interactions: Slender Structures and Axial Flow, vol. 1, Academic Press, California, Calif, USA, 1998.

Copyright

Copyright © 2017 H. B. Wen 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.

PDF Download Citation

Download other formats

Order printed copies