Flow Induced Vibration In Pipes Pdf

Posted on  by
Flow Induced Vibration In Pipes Pdf Average ratng: 3,7/5 1492 votes
  • About this Journal ·
  • Abstracting and Indexing ·
  • Aims and Scope ·
  • Article Processing Charges ·
  • Bibliographic Information ·
  • Editorial Board ·
  • Editorial Workflow ·
  • Publication Ethics ·
  • Reviewer Resources ·
  • Submit a Manuscript ·
  • Subscription Information ·

Cross-flow induced vibration is the most important problem in various fields, and is known to have caused many problems in various industrial components. A classification based on the correspondence between the engineering product of concern and the structural model is first presented. For vibration of pipes excited by internal flow the. Design criteria and guidelines to be used for Design & Supporting for lines prone to Acoustic Induced Vibration (AIV) & Flow Induced Vibration (FIV). And guidelines to be used for Design & Supporting for lines prone to Acoustic Induced Vibration (AIV) & Flow Induced Vibration. In Pipe Support Standard. 4.2 Flow Induced Vibration a) Flow. KEY WORDS: flow-induced vibration, pipeline, vibration analysis, vibration measurement, damper, NETPULS, dPIPE INTRODUCTION Vibration of High Energy pipelines very often is the cause of pipe and fitting cracking and other relevant problems.

  • Annual Issues ·
  • Open Special Issues ·
  • Published Special Issues ·

Flow-Induced Vibration Analysis of Supported Pipes Conveying Pulsating Fluid Using Precise Integration Method

1Logistics Engineering College, Shanghai Maritime University, Shanghai 201306, China
2School of Mechanical and Power Engineering, East China University of Science and Technology, Shanghai 200237, China

Received 1 January 2010; Revised 6 May 2010; Accepted 19 August 2010

Academic Editor: Carlo Cattani

Copyright © 2010 Long Liu and Fuzhen Xuan. 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.

Abstract

Dynamic analysis of supported pipes conveying pulsating fluid is investigated in Hamiltonian system using precise integration method (PIM). First, symplectic canonical equations of supported pipes are deduced with state variable vectors composed of displacement and momentum. Then, PIM with linear interpolation formula is proposed to solve these equations. Finally, this approach's precision is testified by several numerical examples of pinned-pinned pipes with different fluid velocities and frequencies. The results show that PIM is an efficient and rapid approach for flow-induced dynamic analysis o f supported pipes.

1. Introduction

As the pipes are widely used in many industrial fields, flow-induced vibration analysis of pipes conveying fluid has been one of the attractive subjects in structural dynamics. It is well known that pipeline systems may undergo divergence and flutter types of instabilities generated by fluid-structure interaction. Over the last sixty years, extensive studies have been carried out on dynamic analysis of pipeline systems subject to different boundary conditions and loadings. Notable contributions in this area include the works of Chen [1] and Paidoussis [2, 3]. At present, most of the research is concentrated on nonlinear dynamic analysis of pipes conveying pulsating fluid. A recent survey on bifurcations for supported pipes can be found in [4]. Folley and Bajaj [5] considered nonlinear spatial dynamic characteristics of cantilever pipes conveying fluid.

In most cases, the corresponding ordinary differential motion equations of fluid-conveyed pipes are deduced using Galerkin’s method in Lagrange system. Then many numerical methods, such as transfer matrix method, finite element method, perturbation method, Runge-Kutta method, and differential quadrature method, are applied to solve these equations.

For example, Jensen [6] analyzed dynamic behaviors of vibrating pipe containing fluid subject to lateral resonant base excitation using the perturbation method of multiple scales. Yang et al. [7] investigated the effect of fluid viscosity and mass ratio on instability regions of a Kelvin-type viscoelastic pipe conveying harmonically pulsating fluid using multiple scales method. Wang et al. [8] studied the nonlinear dynamics of curved fluid conveying pipe with differential quadrature method.

Jeong et al. [9] proposed a finite element model of pipes conveying periodically pulsating fluid and analyzed the influence of fluid velocities on pipe's stability. Stangl et al. [10] solved the extended version of Lagrange nonlinear equations for cantilevered pipes using implicit Runge-Kutta solver HOTINT. Wang [11] explored numerically the effect of the nonlinear motion constraints on dynamics of simply supported pipes conveying pulsating fluid via the fourth-order Runge-Kutta scheme.

Nikolić and Rajković [12] used Lyapunov-Schmidt reduction and singularity theory to investigate the behaviors of extensive fluid-conveying pipe supported at both ends around the neighborhood of the bifurcation points. Furthermore, Modarres-Sadeghi and Païdoussis [13] studied the possible postdivergence flutter instabilities of this complete nonlinear supported pipe’s model with Houbolt’s finite difference method [14] and AUTO Software package. Xu et al. [15] proposed the analytical expression of natural frequencies of fluid-conveying pipes with the help of homotopy perturbation method. Those calculated frequencies were in good agreement with experiment results.

Considering the effect of the internal and external fluids, the three-dimensional nonlinear differential equations of a fluid-conveying pipe undergoing overall motions were derived based on Kane’s equation and the Ritz method [16]. Moreover, the time histories for the displacements were obtained using the incremental harmonic balance method. Based on Timoshenko beam model, Shen et al. [17] studied the band gap properties of the flexural vibration for periodic pipe system conveying fluid using the transfer matrix method. These methods have proved to be effective in analyzing flow-induced vibration of certain pipes.

It is well known that analysis of pipe dynamics could be conducted based on the energy-based approach according to Hamiltonian principle [18, 19]. However, these approximation methods mentioned above are not ideal for Hamiltonian systems [20], because they are not structurally stable, which means that the Hamiltonian system will become dissipative.

Recently, many numerical algorithms, which can inherit the symplectic structure of Hamiltonian system, have been studied. Especially, Zhong and Williams [21] have proposed the precise integration method, which can give the highly precise numerical integration result and approach the full computer precision for these homogeneous equations. Moreover, this approach has been applied to solving complicated inhomogeneous problems with nonlinear time-variant item, for example, Floquet transition matrix, control problems, and so on [22–25].

In this paper, a Hamiltonian model of nonlinear flow-induced dynamics of supported pipes is analyzed numerically using precise integration method. Firstly, nonlinear equations of supported pipes conveying harmonically fluctuating fluid are deduced to two-order ordinary differential equations using the Galerkin’s method. Then the equations are transformed into symplectic canonical equations composed of displacement and momentum. Moreover, PIM with linear interpolation formula is proposed. Finally, several numerical examples of pinned-pinned pipe conveying pulsating fluid are used to testify the precision of this approach. The results are compared with those using traditional Runge-Kutta method. The influence of different fluid parameters on nonlinear behaviors of supported pipes is also discussed.

2. Formulation of Problem in Hamiltonian System

In this section, typical governing equations of supported pipe conveying fluid are deduced in Hamiltonian system.

2.1. Equation of Motion

We consider a straight supported pipe conveying the harmonically pulsating flow (Figure 1). It is assumed that the motion is planar, and the pipe is nominally horizontal. The cross-sectional area of the flow is assumed constant. The effects of gravity and external tension are ignored. Moreover, the pipe behaves like an Euler-Bernoulli beam in transverse vibration and the fluid is assumed to be incompressible.

Figure 1: The simply supported pipe conveying fluid.

The transverse motion equation of the pipe is given by Païdoussis and Issid [26], where is the longitudinal coordinate, the transverse deflection, the moment of flexure of the pipe, the pipe length, the mass of the fluid conveyed per unit length, pipe mass per unit length, the fluid velocity, the fluid pressure, and the cross-sectional area of the flow.

Then the viscoelastic Kelvin-Voigt damping model is introduced, where is the flexural stiffness of the pipe material, and is the coefficient of Kelvin-Voigt viscoelastic damping.

Moreover, the velocity of pulsating fluid is assumed to be harmonically fluctuating, and has the following form: where is the mean flow velocity, the amplitude of the harmonic fluctuation (assumed small), and the fluid pulsating frequency. This fluctuating flow velocity appears as parametric excitation term in the equation of motion and may lead to parametric instabilities.

Substituting (2.2) and (2.3) into (2.1) yields that

Incorporate the following dimensionless quantities:

Then the equation of motion can be nondimensionalized as

The motion equation above is inhomogeneous, as the derivative coefficients of are explicit functions of and .

Then we discretize (2.6) using the Galerkin’s method. Let where () are generalized coordinates of the discretized pipe and are eigenfunctions of the beam with the same boundary conditions.

It has been pointed out that instability boundaries for supported pipes could be determined with adequate precision using the two-mode expansion [2]. So the two-mode expansion of (2.7) is used in the analytical model for simplicity to investigate qualitative behaviors of supported pipes conveying fluid.

Substitute (2.7) with into (2.6). Then according to the orthogonal property of modal modes, the partial differential equation could be transformed into the second-order ordinary differential equation where

In (2.9), and denote the structural damping matrix and stiffness matrix, respectively. These two matrices are associated with systematic parameters, such as dimensionless flow velocity and mass ratio . () are the th eigenvalues of the supported pipe and is the diagonal matrix with elements .

Moreover, , , and are matrices with elements , , and (), respectively. They are defined as

Different value should be taken for those three parameters depending on different boundary conditions of the pipe. For the pinned-pinned pipe, we have

2.2. Symplectic Canonical Equation

According to Hamiltonian principle, the nonlinear equation of supported pipe conveying pulsating fluid can be transformed into symplectic canonical equation with state variable vectors composed of displacement and corresponding momentum , where

So, we can see that is a Hamiltonian symplectic matrix and is a time-variant matrix related to state variable vectors.

3. Precise Integration Method with Linear Interpolation Approximation

In this section, the principle of precise integration method is briefly introduced. For a more detailed explanation, it is suggested that [21, 22] are consulted.

The precise integration method for homogeneous equations with initial value is fundamental, so it is described in the next subsection firstly.

3.1. Integration of Homogeneous Equation

The general solution of homogeneous equation can be expressed as

Suppose that the time step is , and then we have the following recursive steps: where . It is seen that how to compute the exponential matrix is essential for the integration precision.

Then split the time interval into a smaller one. Define and . For example, when . As is small, is an extremely small time interval.

Assume that , execute the cycle where is the identity matrix.

So is no longer a very small matrix. It could be computed by the following function:

The algorithm given above is called precise integration method. It has no serious numerical round-off error and could approach full computer precision [20].

3.2. Integration of Inhomogeneous Equation

In this subsection, PIM with linear interpolation formula would be proposed to solve inhomogeneous equations.

With the solution of homogeneous equation, (2.13) could been written as

Then its solution could be given by the Duhamel’s integration as

Similarly, the duration of structural dynamic response is also divided into small time intervals. The response between can be written as

To solve this inhomogeneous equation, the analytical expression of the time-variant item is required. But it is not always available.

In this study, the linear interpolation formula is used to approximate this nonlinear item. Motion for leave to file demurrer to evidence.

Assume that where

Substituting (3.8) into (3.7) gives the linear interpolation expression

Thus, we have the numerical expression of symplectic canonical equation using PIM with linear interpolation formula.

In the next section, this method would be used to investigate the motion of supported pipes conveying pulsating fluid under different conditions.

4. Numerical Examples

In this section, several numerical examples of pinned-pinned pipes are used to testify the effectiveness of precise integration method.

4.1. Dynamic Response of Pipes Conveying Stable Fluid

In this subsection, this approach is used to analyze the dynamic response of stable fluid-conveying pipes, especially for their computation stability after a long period. In this case, the pipe’s dynamic function is a homogeneous equation. The results are compared with those using traditional forth-order Runge-Kutta method.

Consider that the dimensionless mean flow velocity is 2.0, the mass ratio is 0.32 and the fluid pressure . The initial conditions are chosen to be . Time increases from 0 to 1000 s and the time step is selected as 0.2 s.

Figure 2 illustrates time history of four state variables of pipe’s middle point using the Runge-Kutta method, while Figure 3 shows the results calculated by precise integration method.

Figure 2: Time history of four state variables using Runge-Kutta method.
Figure 3: Time history of four state variables using precise integration method.

It can be found that there are evident differences for four state variables’ amplitudes using two methods. The amplitudes in Figure 2 decrease gradually with time. When the simulation time is long enough, state variables may converge to zero. However, those in Figure 3 still keep constant with time, which are almost unaffected by the time step.

So we can conclude that there is the energy dissipation using traditional Runge-Kutta method, which cannot get the accurate numerical results. However, precise integration method is an energy conservative method and could maintain the stability of the numerical simulation in the long period of time.

Table 1 lists the computation time needed for two methods, as two different time intervals are selected. It can be noted that precise integration method needs much less computing time than Runge-Kutta method.

Table 1: Computation time needed using two methods.
4.2. Dynamic Response of Pipes Conveying Pulsating Fluid

In this subsection, PIM with linear approximation is used to analyze the dynamic response of supported pipes conveying harmonically pulsating fluid. Similarly, the results are compared with those using forth-order Runge-Kutta method.

Consider the dimensionless mean flow velocity is 1, the amplitude , the frequency , the mass ratio , and the fluid pressure . The initial conditions are chosen to be . Time increases from 0 to 100 s, and two time steps involved in this example are selected as s and 0.05 s, respectively. By the way, the relative stable dynamic response after 20 s is considered.

Figure 4 illustrates time history of displacement response and phase planes of pipe’s middle point using the forth-order Runge-Kutta method, while Figure 5 shows the results calculated by PIM, when the time step is 0.01 s. Figures 6 and 7 show the results with time step 0.05 s using two methods, respectively.

Flow Induced Vibration Of Pipeline

Figure 4: Displacement response and phase diagram of pipe’s middle point using Runge-Kutta .
Figure 5: Displacement response and phase diagram of pipe’s middle point using PIM .
Figure 6: Displacement response and phase diagram of pipe’s middle point using Runge-Kutta ().
Figure 7: Displacement response and phase diagram of pipe’s middle point using PIM .

These figures show that phrase planes of two modes would shrink gradually with time and converge to a point if time is long enough. The precise integration method shows nearly the same precision during calculating the dynamic response of supported pipes conveying pulsating fluid.

Furthermore, Table 2 lists the computation time needed for two methods. Similarly, PIM with linear interpolation formula need much less computing time than Runge-Kutta method. This approach is very quick to obtain dynamic response because of running a number of cycles during the computation, which is shown in (3.3). So it is suitable for long-term dynamic analysis of fluid-conveyed pipes.

Table 2: Computation time needed using two methods.
4.3. Stability Analysis of Pipes under Different Fluid Velocities and Frequencies

In this subsection, the influence of different fluid parameters on nonlinear behaviors of pinned-pinned pipes is discussed using PIM with linear interpolation formula.

Energy Institute Flow Induced Vibration

The dimensionless fluid frequency increases from 0 to 70, and three fluid velocities in this example are selected as , 1.5, and 2.0. The fluid frequency step and the time step are selected as and , respectively. Others parameters are the same with the preceding Section 4.2.

The displacement response of the pipe’s middle point under different fluid parameters is calculated. Figure 8 shows the displacement responses of the middle point versus the pulsating fluid frequency. It can be seen that the pipe keeps stable at most frequencies domain. For example, as , the pipe behaves stable on the limit loop condition, which is shown as Figure 9(a). However, when the frequencies lie between and , the pipes are unstable. For example, the pipe will be divergent as Figure 9(b) shows when . So, it is very dangerous for pipes’ operating safety.

Figure 8: Displacement response of pipe’s middle point as the function of the dimensionless pulsating fluid frequency .
Figure 9: Phase diagram of pipe’s middle point as . (a) ; (b) .

Figures 10 and 11 show the displacement response variations of the middle point versus the pulsating fluid frequency as and 2.0, respectively. It can be shown that the instable zone is changing with fluid frequencies. As , the instable zones are at , , and . As , the instable zones lie at , , , and .

Figure 10: Displacement response of pipe’s middle point as the function of the dimensionless pulsating fluid frequency .
Figure 11: Displacement response of pipe’s middle point as the function of the dimensionless pulsating fluid frequency .

It can be seen that with the increasing of fluid velocity, the critical fluid frequency gets smaller and the pipe shows complicated nonlinear vibration.

5. Conclusion

In this study, PIM with linear interpolation formula is presented to analyze nonlinear dynamics of Hamiltonian model of supported straight pipe conveying pulsating fluid. Several numerical examples are used to testify the effectiveness of this approach. The results show this approach could keep stable even for long period of time, and is much more rapid than traditional Runge-Kutta method. Moreover, the pipe’s nonlinear behaviors under the condition of different fluid parameters are discussed.

The work presented here provides an alternative approach for investigating the nonlinear dynamic response of the pipes conveying fluid. However, it should be pointed out that linear interpolation formula is a rough approximation method, and more accurate methods should be studied to analyze nonlinear flow-induced dynamics in Hamiltonian system.

Flow Induced Vibration In Pipes Pdf Free

References

Flow Induced Vibration In Heat Exchanger

  1. S. S. Chen, “Dynamic stability of tube conveying fluid,” Journal of Engineering Mechanical Division, vol. 97, pp. 1469–1485, 1971.View at Google Scholar
  2. M. P. Paidoussis, Fluid-Structure Interactions: Slender Structures and Axial Flow V1, Academic Press, Amsterdam, The Netherlands, 1998.
  3. M. P. Paidoussis, Fluid-Structure Interactions: Slender Structures and Axial Flow V2, Academic Press, Amsterdam, The Netherlands, 2004.
  4. M. Nikolić and M. Rajković, “Bifurcations in nonlinear models of fluid-conveying pipes supported at both ends,” Journal of Fluids and Structures, vol. 22, no. 2, pp. 173–195, 2006.View at Publisher · View at Google Scholar · View at Scopus
  5. C. N. Folley and A. K. Bajaj, “Spatial nonlinear dynamics near principal parametric resonance for a fluid-conveying cantilever pipe,” Journal of Fluids and Structures, vol. 21, no. 5–7, pp. 459–484, 2005.View at Publisher · View at Google Scholar · View at Scopus
  6. J. S. Jensen, “Fluid transport due to nonlinear fluid-structure interaction,” Journal of Fluids and Structures, vol. 11, no. 3, pp. 327–344, 1997.View at Publisher · View at Google Scholar · View at Scopus
  7. X. Yang, T. Yang, and J. Jin, “Dynamic stability of a beam-model viscoelastic pipe for conveying pulsative fluid,” Acta Mechanica Solida Sinica, vol. 20, no. 4, pp. 350–356, 2007.View at Publisher · View at Google Scholar · View at Scopus
  8. L. Wang, Q. Ni, and Y.-Y. Huang, “Hopf bifurcation of a nonlinear restrained curved fluid conveying pipe by differential quadrature method,” Acta Mechanica Solida Sinica, vol. 16, no. 4, pp. 345–352, 2003.View at Publisher · View at Google Scholar · View at Scopus
  9. W. B. Jeong, Y. S. Seo, S. H. Jeong, S. H. Lee, and W. S. Yoo, “Stability analysis of a pipe conveying periodically pulsating fluid using finite element method,” Mechanical Systems Machine Elements and Manufacturing, vol. 49, no. 4, pp. 1116–1122, 2007.View at Publisher · View at Google Scholar · View at Scopus
  10. M. Stangl, J. Gerstmayr, and H. Irschik, “An alternative approach for the analysis of nonlinear vibrations of pipes conveying fluid,” Journal of Sound and Vibration, vol. 310, no. 3, pp. 493–511, 2008.View at Publisher · View at Google Scholar · View at Scopus
  11. L. Wang, “A further study on the non-linear dynamics of simply supported pipes conveying pulsating fluid,” International Journal of Non-Linear Mechanics, vol. 44, no. 1, pp. 115–121, 2009.View at Publisher · View at Google Scholar · View at Scopus
  12. M. Nikolić and M. Rajković, “Bifurcations in nonlinear models of fluid-conveying pipes supported at both ends,” Journal of Fluids and Structures, vol. 22, no. 2, pp. 173–195, 2006.View at Publisher · View at Google Scholar · View at Scopus
  13. Y. Modarres-Sadeghi and M. P. Païdoussis, “Nonlinear dynamics of extensible fluid-conveying pipes, supported at both ends,” Journal of Fluids and Structures, vol. 25, no. 3, pp. 535–543, 2009.View at Publisher · View at Google Scholar · View at Scopus
  14. C. Semler, W. C. Gentleman, and M. P. Païdoussis, “Numerical solutions of second order implicit non-linear ordinary differential equations,” Journal of Sound and Vibration, vol. 195, no. 4, pp. 553–574, 1996.View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  15. M.-R. Xu, S.-P. Xu, and H.-Y. Guo, “Determination of natural frequencies of fluid-conveying pipes using homotopy perturbation method,” Computers and Mathematics with Applications, vol. 60, pp. 520–527, 2010.View at Publisher · View at Google Scholar · View at Scopus
  16. D. Meng, H. Y. Guo, and S. P. Xu, “Nonlinear dynamic model of a fluid-conveying pipe undergoing overall motions,” Applied Mathematical Modelling, vol. 35, no. 2, pp. 781–796, 2011.View at Publisher · View at Google Scholar
  17. H.-J. Shen, J.-H. Wen, D.-L. Yu, and X.-S. Wen, “Flexural vibration property of periodic pipe system conveying fluid based on Timoshenko beam equation,” Acta Physica Sinica, vol. 58, no. 12, pp. 8357–8363, 2009.View at Google Scholar · View at Scopus
  18. M. Stangl, N. A. Beliaev, and A. K. Belyaev, “Applying Lagrange equations and Hamilton’s principle to vibrations of fluid conveying pipes,” in Proceedings of the 33th Summer School on Advanced Problems in Mechanics (APM '05), pp. 269–275, St. Petersburg, Russia, 2005.
  19. M. Stangl and H. Irschik, “Dynamics of an Euler elastic pipe with internal flow of fluid,” Proceedings of Applied Mathematics and Mechanics, vol. 6, pp. 335–336, 2006.View at Publisher · View at Google Scholar
  20. K. Feng and M. Z. Qin, “Hamiltonian algorithms for Hamiltonian dynamical systems,” Progress in Natural Science, vol. 1, no. 2, pp. 105–116, 1991.View at Google Scholar
  21. W. X. Zhong and F. W. Williams, “A precise time step integration method,” Journal of Mechanical Engineering Science, vol. 208, no. 6, pp. 427–430, 1994.View at Publisher · View at Google Scholar · View at Scopus
  22. W. X. Zhong, “On precise integration method,” Journal of Computational and Applied Mathematics, vol. 163, no. 1, pp. 59–78, 2004.View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  23. J. Lin, W. Shen, and F. W. Williams, “A high precision direct integration scheme for structures subjected to transient dynamic loading,” Computers and Structures, vol. 56, no. 1, pp. 113–120, 1995.View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  24. G. Zhou, Y. X. Wang et al., “A homogenized high precise direct integration based on Taylor serials,” Journal of Shanghai Jiaotong University, vol. 35, no. 12, pp. 1916–1919, 2001.View at Google Scholar
  25. T. C. Fung, “Construction of higher-order accurate time-step integration algorithms by equal-order polynomial projection,” Journal of Vibration and Control, vol. 11, no. 1, pp. 19–49, 2005.View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  26. M. P. Païdoussis and N. T. Issid, “Dynamic stability of pipes conveying fluid,” Journal of Sound and Vibration, vol. 33, no. 3, pp. 267–294, 1974.View at Publisher · View at Google Scholar · View at Scopus