Dynamic stiffness method for free vibration analysis of variable diameter pipe conveying fluid
Liu Yongshou^{1} , Zhang Zijun^{2} , Li Baohui^{3} , Gao Hangshan^{4}
^{1, 2, 3, 4}Department of Engineering Mechanics, Northwestern Polytechnical University, Xi’an, P. R. China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 16, Issue 2, 2014, p. 832845.
Received 7 December 2013; received in revised form 10 January 2014; accepted 24 January 2014; published 31 March 2014
JVE Conferences
The governing equation of flowinduced vibration is deduced in terms of Hamilton’s principle for a variable diameter pipe conveying axial steady flow. Frobenius method is applied to analyze the governing equation. After the recursion formulas of coefficients are obtained, dynamic stiffness method is proposed for free vibration analysis of the variable diameter pipe conveying fluid. In the example, the natural frequencies of uniform pipes conveying fluid are computed and comparisons are made to validate the dynamic stiffness method. Then, the natural frequencies and modal shapes are obtained for the variable diameter pipe conveying fluid with different section variation coefficients and fluid velocities.
Keywords: variable diameter pipe conveying fluid, dynamic stiffness method, natural frequencies, modal shapes.
1. Introduction
Fluidconveying pipe is widely used in nuclear industry, petroleum industry, and many other industries. The failure caused by the interaction between fluid and pipe wall occurs very frequently in practice. So flowinduced vibration of pipe conveying fluid attracts more and more interests of researchers. Ashley and Haviland [1] begun the research of bending vibrations of pipe line containing flowing fluid. Paidoussis [2] has studied the dynamic characteristics of pipe conveying fluid for many years, the great achievement he accomplished in nonlinear vibration of pipes set the foundation of pipe dynamics. Recently, flowinduced vibration of pipe conveying fluid is researched more exhaustively. Doare et al. [3] studied the dissipation effect on the local and global stability of pipe conveying fluid and proposed a numerical method to analyze the stability of finitelength system. Jung, Chung [4] investigated the stability of semicircular pipe conveying harmonically oscillating fluid. They considered the extensibility and nonlinearity of the fluidconveying pipe with different boundary conditions through Hamilton principle. Rinaldi, Prabhakar [5] studied microscale resonators containing internal flow, modeled as microfabricated pipe conveying fluid, and investigated the effects of flow velocity on damping, stability, and frequency shift. Huang, Liu [6] studied the natural frequencies of fluid conveying pipe with different boundary conditions through Garlerkin approach. Ragulskis, Fedaravicius [7] have developed a harmonic balance method in the FEM analysis of the fluid in pipe with taking the interactions between the vibrating pipe and fluid flow into consideration.
Uniform pipe was always adopted in the flowinduced vibration analysis of pipe conveying fluid. However, fluidconveying pipes with variable section are widespread in nuclear engineering, such as pipeexpanding and nozzle etc. Hannoyer and Paidoussis [8, 9] have studied the dynamics of a slender tapered beam with internal and external flow in theory and experiment. Based on the researches of Hannoyer and Paidoussis [8, 9], we studied the dynamic characteristics of variable diameter pipes conveying fluid through stiffness matrix method. Natural frequencies and modal shapes of the variable diameter pipe were calculated in this paper and contrasts were made between the results in this paper and previous researches, which prove that the method employed in this paper is reasonable.
2. Motion equation of variable diameter pipe conveying fluid
Here, we use Hamilton’s principle to build the dynamic model of variable diameter pipe conveying fluid. Because that the fluidconveying pipe is an open system. The Hamilton’s principle for open system is used here, i.e.:
where, $\mathcal{L}$ is the system’s Lagrange function, that is $\mathcal{L}={T}_{k}{T}_{p}$, where ${T}_{k}$, ${T}_{p}$ stand for the kinetic energy and potential energy of the system respectively.$\delta W$ is the work done by nonconservative force, ${m}_{f}$ is the fluid mass per unit length, ${\overrightarrow{r}}_{e}$ is the position vector of pipe exit. ${\overrightarrow{e}}_{t}$ is the unit vector in tangent direction of the deformed pipeline, as shown in Fig. 1.
In the coordinate shown in Fig. 1, the absolute fluid velocity can be written as follow:
In Eq. (2), $V\left(x\right)$ stands for the velocity of the fluid relative to the pipe. ${\overrightarrow{e}}_{1}$ and ${\overrightarrow{e}}_{2}$ are unit vectors in $x$ and $y$directions.$\left(\right)\text{'}$ and $\dot{\left(\right)}$ stand for the partial derivatives with respect to $x$ (coordinate) and $t$ (time). According to the law of conservation of fluid mass, as shown in Eq. (3), the velocity of the fluid is change along the length of the conical pipe. So we use $V\left(x\right)$ and ${V}_{a}\left(x\right)$ here:
where $C$ is a constant.
Fig. 1. The displacement sketch of simply supported pipe
Fig. 2. The sketch of variable diameter pipe
The sketch of a variable diameter pipe is shown in Fig. 2. The fluid flows from left end to right end. The outer and inner diameters of left end of the pipe are ${D}_{i}$ and ${d}_{i}s$ respectively. The nonuniformity of the pipe is characterized by the diameter variation coefficient $\gamma $. The outer and inner diameters of the section at coordinate $x$ can be expressed as $D\left(x\right)=\left(1\gamma \frac{x}{L}\right){D}_{i}$ and $d\left(x\right)=\left(1\gamma \frac{x}{L}\right){d}_{i}$ respectively, where $L$ is length of the pipe. Because that the area of the fluid is ${A}_{f}\left(x\right)=\frac{1}{4}\pi {d}^{2}\left(x\right)$, the velocity of the fluid along the length of the pipe can be expressed as:
Given that the flow velocity at the entrance of the pipe is ${V}_{0}$, then Eq. (4) can be transformed into following form:
Now, we use Hamilton’s principle for open system to deduce the motion equation of variable diameter pipe conveying fluid. The kinetic energy of fluid in pipe can be written in the following form:
Here, the axial inextensible assumption was used, i.e. $\mathrm{}{\left(\frac{\partial u}{\partial x}\right)}^{2}+{\left(\frac{\partial w}{\partial x}\right)}^{2}=1$. And high order infinitesimals have been ignored in Eq. (6). ${\rho}_{f}$ is the density of fluid in pipe, ${A}_{f}\left(x\right)$ is area of the fluid section located at $x$. ${m}_{f}\left(x\right)dx={\rho}_{f}{A}_{f}\left(x\right)dx$ is mass of the fluid element analyzed.
The kinetic energy of the pipe is as follow:
where, ${\rho}_{p}$ is density of the pipe material, ${A}_{p}\left(x\right)$ is area of the pipe cross section at $x$. ${m}_{p}\left(x\right)dx={\rho}_{p}{A}_{p}\left(x\right)dx$ is mass of the pipe element.
The potential energy of the pipe is as follow:
In Eq. (8), ${I}_{x}$ is the inertial moment of the pipe section. Obviously, ${I}_{x}$ also changes along $x$axis. Only considering the work done by fluid pressure, the work done by nonconservative force is zero [2], that is:
Substituting Eq. (6)(9) into Eq. (1), the following equation can be obtained:
Here, we take the simply supported pipe (as shown in Fig. 1) for example to deduce the lateral vibration equation of the variable diameter pipe conveying fluid.
The boundary conditions of simply supported pipe are:
At time ${t}_{1}$ and ${t}_{2}\text{:}$
Combine Eq. (11) with Eq. (12), the following equation is obtained:
The detail manipulation of the third term in Eq. (10) is as follow:
Considering the assumption of axial inextensible, we obtain that ${u}_{e}=\frac{1}{2}\underset{0}{\overset{L}{\int}}{{w}^{\text{'}}}^{2}dx$, where the subscript ‘$e$’ means “exit”, and the following equation can be gained:
Now, we can get the motion equation of variable diameter pipe conveying fluid as follow:
3. Free vibration analysis of variable diameter pipe
3.1. Frobenius method
To harmonic vibration, we have:
Substituting Eq. (17) into Eq. (16), we can get the following ordinary differential equation:
where:
${m}_{p}\left(x\right)={\rho}_{p}{A}_{p}\left(x\right)={\rho}_{f}\pi \frac{{\left(1\frac{\gamma x}{L}\right)}^{2}\left({D}_{i}^{2}{d}_{i}^{2}\right)}{4}=\left(1+{\alpha}_{1}\frac{x}{L}+{\alpha}_{2}\frac{{x}^{2}}{{L}^{2}}\right){\rho}_{p}{A}_{p}\left(0\right),$
and
For the purpose of analyzing Eq. (18) conveniently, the following dimensionless parameters are introduced here:
${v}_{0}=\sqrt{\frac{{\rho}_{f}{A}_{f}\left(0\right)}{EI\left(0\right)}}L{V}_{0},\beta =\frac{{\rho}_{f}{A}_{f}\left(0\right)}{{\rho}_{f}{A}_{f}\left(0\right)+{\rho}_{p}{A}_{p}\left(0\right)}.$
The dimensionless form of Eq. (18) is as follow:
Given that the contraction angle of the conical pipe is very small, that is, ($\gamma \xi $) is sufficiently small so that we can use Tylor series here as:
$\frac{1}{{(1\gamma \xi )}^{2}}\approx 1+2\gamma \xi +3{\left(\gamma \xi \right)}^{2}+{4\left(\gamma \xi \right)}^{3}+{5\left(\gamma \xi \right)}^{4}+{6\left(\gamma \xi \right)}^{5}+7{\left(\gamma \xi \right)}^{6},$
$\frac{1}{{(1\gamma \xi )}^{3}}\approx 1+3\gamma \xi +6{\left(\gamma \xi \right)}^{2}+{10\left(\gamma \xi \right)}^{3}+{15\left(\gamma \xi \right)}^{4}+{21\left(\gamma \xi \right)}^{5}+{28\left(\gamma \xi \right)}^{6},$
$\frac{1}{{(1\gamma \xi )}^{4}}\approx 1+4\gamma \xi +10{\left(\gamma \xi \right)}^{2}+{20\left(\gamma \xi \right)}^{3}+35{\left(\gamma \xi \right)}^{4}+56{\left(\gamma \xi \right)}^{5}+{84\left(\gamma \xi \right)}^{6},$
$\frac{1}{{(1\gamma \xi )}^{5}}\approx 1+5\gamma \xi +15{\left(\gamma \xi \right)}^{2}+{35\left(\gamma \xi \right)}^{3}+{70\left(\gamma \xi \right)}^{4}+{126\left(\gamma \xi \right)}^{5}+{210\left(\gamma \xi \right)}^{6}.$
Here, we apply Frobenius method to analyze Eq. (22). Given that the solution of Eq. (22) takes the following form:
In which ${a}_{n}$ is the coefficient of the polynomial. The derivatives of $\eta \left(\xi \right)$ are:
$\frac{{d}^{3}\eta}{d{\xi}^{3}}=\sum _{n=0}^{\infty}{a}_{n}\left(n+k\right)\left(n+k1\right)\left(n+k2\right){\xi}^{n+k3},$
$\frac{{d}^{4}\eta}{d{\xi}^{4}}=\sum _{n=0}^{\infty}{a}_{n}(n+k)(n+k1)(n+k2)(n+k3){\xi}^{n+k4}.$
When $n=0$, substituting Eq. (24) and Eq. (25) into Eq. (22), because that the lowest power of $\xi $ equals to zero, so the following equation can be obtained:
Because that ${a}_{0}\ne 0$, so we obtain:
Eq. (27) is namely indicial equation and it plays an important role in the analysis [10]. Obviously, Eq. (27) have four roots, i.e. ${k}_{j}=$0, 1, 2, 3 ($j=$1, 2, 3, 4) respectively. To each $k$, only one coefficient ${a}_{n}\left(k\right)$ is related. Substituting Eq. (24) and Eq. (25) into Eq. (22), then all the coefficients can be obtained:
$\frac{{\beta}_{3}\left(k+1\right)k\left(k1\right)+{\alpha}_{1}{\zeta}_{1}k+4{\zeta}_{1}\gamma \left(k1\right)+\left({\zeta}_{1}+{\zeta}_{2}\right)}{\left(k+3\right)\left(k+2\right)\left(k+1\right)}{a}_{0}\left(k\right),$
$\frac{{\beta}_{3}\left(k+2\right)\left(k+1\right)k+{\alpha}_{1}{\zeta}_{1}\left(k+1\right)+4\gamma {\zeta}_{1}k+\left({\zeta}_{1}+{\zeta}_{2}\right)}{\left(k+4\right)\left(k+3\right)\left(k+2\right)}{a}_{1}\left(k\right)\frac{\left[\begin{array}{c}{\beta}_{4}\left(k+2\right)\left(k+1\right)k\left(k1\right)+\left(4{\alpha}_{1}\gamma +10{\gamma}^{2}+{\alpha}_{2}\right){\zeta}_{1}k\left(k1\right)\\ +\left[\left(4{\alpha}_{1}\gamma +2{\alpha}_{2}+5\gamma +{\alpha}_{1}\right){\zeta}_{1}+2\gamma {\zeta}_{2}\right]k+\left({\alpha}_{1}{\zeta}_{3}+{\zeta}_{4}{\zeta}_{5}\right)\end{array}\right]}{(k+4)(k+3)(k+2)(k+1)}{a}_{0}\left(k\right),$
$\frac{\left[\begin{array}{c}{\beta}_{3}\left(k+n+3\right)\left(k+n+2\right)\left(k+n+1\right)\\ +{\alpha}_{1}{\zeta}_{1}\left(k+n+2\right)+4\gamma {\zeta}_{1}\left(k+n+1\right)+\left({\zeta}_{1}+{\zeta}_{2}\right)\end{array}\right]}{\left(k+n+5\right)\left(k+n+4\right)\left(k+n+3\right)}{a}_{n+2}\left(k\right)\frac{\left\{\begin{array}{c}{\beta}_{4}\left(k+n+3\right)\left(k+n+2\right)\left(k+n+1\right)\left(k+n\right)\\ +\left(4{\alpha}_{1}\gamma {\zeta}_{1}+10{\gamma}^{2}{\zeta}_{1}+{\alpha}_{2}\right)\left(k+n+1\right)\left(k+n\right)\\ +\left[\left(4{\alpha}_{1}\gamma +2{\alpha}_{2}+5\gamma +{\alpha}_{1}\right){\zeta}_{1}+2\gamma {\zeta}_{2}\right]\left(k+n+1\right)+\left({\alpha}_{1}{\zeta}_{3}+{\zeta}_{4}{\zeta}_{5}\right)\end{array}\right\}}{\left(k+n+5\right)\left(k+n+4\right)\left(k+n+3\right)\left(k+n+2\right)}{a}_{n+1}\left(k\right)\frac{{\zeta}_{3}\left(2{\alpha}_{2}+2{\alpha}_{1}\gamma \right)+{\zeta}_{4}\left(3\gamma +{\alpha}_{1}\right){\zeta}_{5}{\alpha}_{1}}{\left(k+n+5\right)\left(k+n+4\right)\left(k+n+3\right)\left(k+n+2\right)}{a}_{n}\left(k\right),$
where:
Here, the high orders of $\gamma $ have been neglected. After getting these coefficients, the solution of Eq. (22) can be written in the following form:
where ${A}_{j}$ is constant. The expression of ${f}_{j}\left(\xi \right)$ is as follow:
${f}_{2}{\left(\xi \right)=a}_{0}{(k}_{2}{)\xi +a}_{1}{(k}_{2}{)\xi}^{2}{+a}_{2}{(k}_{2}{)\xi}^{3}{+a}_{3}{(k}_{2}{)\xi}^{4}+\dots ,$
${f}_{3}{\left(\xi \right)=a}_{0}{(k}_{3}{)\xi}^{2}{+a}_{1}{(k}_{3}{)\xi}^{3}{+a}_{2}{(k}_{3}{)\xi}^{4}{+a}_{3}{(k}_{3}{)\xi}^{5}+\dots ,$
${f}_{4}{\left(\xi \right)=a}_{0}{(k}_{4}{)\xi}^{3}{+a}_{1}{(k}_{4}{)\xi}^{4}{+a}_{2}{(k}_{4}{)\xi}^{5}{+a}_{3}{(k}_{4}{)\xi}^{6}+\dots .$
Fig. 3. The sketch of node force
3.2. Dynamic stiffness method
Here, we use dynamic stiffness method to analyze vibration of variable diameter pipe conveying fluid. As shown in Fig. 3, the node displacement of the pipe is:
where:
$P=\left[\begin{array}{cccc}{a}_{0}\left({k}_{1}\right)& 0& 0& 0\\ {a}_{1}\left({k}_{1}\right)& {a}_{0}\left({k}_{2}\right)& 0& 0\\ {f}_{1}\left(1\right)& {f}_{2}\left(1\right)& {f}_{3}\left(1\right)& {f}_{4}\left(1\right)\\ {f}_{1}^{\mathrm{\text{'}}}\left(1\right)& {f}_{2}^{\mathrm{\text{'}}}\left(1\right)& {f}_{3}^{\mathrm{\text{'}}}\left(1\right)& {f}_{4}^{\text{'}}\left(1\right)\end{array}\right].$
As shown in Fig. 3, the shear force and bending moment can be expressed as:
$Q=\frac{\partial M}{\partial x}=E{I}_{x}\frac{{\partial}^{3}w}{\partial {x}^{3}}E\frac{\partial {I}_{x}}{\partial x}\frac{{\partial}^{2}w}{\partial {x}^{2}}.$
Introducing the dimensionless form of shear force and bending moment as follows:
$\stackrel{}{Q}=\frac{Q{L}^{2}}{E{I}_{i}}=\left(1+{\beta}_{1}\xi +{\beta}_{2}{\xi}^{2}+{\beta}_{3}{\xi}^{3}+{\beta}_{4}{\xi}^{4}\right)\frac{{d}^{3}\eta}{d{\xi}^{3}}$
$\left({\beta}_{1}+2{\beta}_{2}\xi +3{\beta}_{3}{\xi}^{2}+4{\beta}_{4}{\xi}^{3}\right)\frac{{d}^{2}\eta}{d{\xi}^{2}}.$
The directions of $M$, $Q$ are shown in Fig. 3, then it is easy to get the vector of node force:
where ${\stackrel{}{Q}}_{1}=\stackrel{}{Q}\left(0\right)$, ${\stackrel{}{M}}_{1}=\stackrel{}{M}\left(0\right)$, ${\stackrel{}{Q}}_{2}=\stackrel{}{Q}\left(1\right)$, ${\stackrel{}{M}}_{2}=\stackrel{}{M}\left(1\right)$.
The relation between vector of node force and vector $A$ can now be expressed as follow:
where:
In which ${\sigma}_{1}=1+{\beta}_{1}+{\beta}_{2}+{\beta}_{3}+{\beta}_{4}$, ${\sigma}_{2}={\beta}_{1}+2{\beta}_{2}+3{\beta}_{3}+4{\beta}_{4}$, and what should be noted here is that, $\left(\right)\text{'}$ stands for the derivative with respect to $\xi $, which is different from that in Section 2.
According to Eq. (32), we can get the coefficient vector as follow:
Now, substituting Eq. (39) into Eq. (37), the following result will be gained:
Here, the matrix $S$ is named as the dynamic stiffness matrix of the system. Once the dynamic stiffness matrix is obtained, the natural frequencies and modal shapes can be easily obtained.
Here, we take the simply supported pipe for an example to illuminate how to compute natural frequencies of the pipe through dynamic stiffness method. The boundary conditions of simply supported pipe can be written as:
$\xi =1:{\eta}_{2}=0,\mathrm{}\mathrm{}\mathrm{}{\stackrel{}{M}}_{2}=0.$
Rewriting Eq. (41) in the matrix form, we have:
The following characteristic equation can be obtained from Eq. (42):
It is easy to find from Eq. (43) that the following equation must be satisfied:
The natural frequencies can be obtained from Eq. (44). In the same way, we can get the characteristic equation for cantilevered pipe and clampedpinned pipe conveying fluid as follows:
Cantilevered pipe (left clamped and right free):
Clampedpinned pipe (left clamped and right pinned):
4. Examples
4.1. Uniform pipe conveying fluid
In order to validate the proposed method, we calculated the natural frequencies of a uniform pipe conveying fluid. The uniform pipe is actually the case of variable diameter pipe with section variation coefficient $\gamma =\text{0}$. The pipe and fluid parameters are chosen the same as previous reference, that is, Yang’s modulus is $E=\text{210}\text{}\text{GPa}$, the outer and inner diameters at origin of the pipe are ${D}_{i}=\text{324}\text{}\text{mm}$ and ${d}_{i}=\text{292}\text{}\text{mm}$ respectively, and the pipe span is $L=\text{32}\text{}\text{m}$, the density of pipe material is ${\rho}_{p}=\text{8200}{\text{}\text{kg/m}}^{\text{3}}$, the density of fluid is ${\rho}_{f}=\text{908.2}{\text{}\text{kg/m}}^{\text{3}}$ [11]. The first five orders of natural frequencies are computed under three different velocities, i.e. $V=\text{0}$, $V=\text{15}\text{}\text{m/s}$, $V=\text{25}\text{}\text{m/s}$ respectively. And comparisons are presented between the results obtained by the proposed method and that published in the research of Housner [11]. The computation results of natural frequencies and error are listed in Table 1.
Table 1. Natural frequencies of uniform simply supported pipe conveying fluid under different fluid velocities (rad/s)
Fluid velocity

Natural frequency


${\omega}_{1}$

${\omega}_{2}$

${\omega}_{3}$

${\omega}_{4}$

${\omega}_{5}$


$V=\text{0}$

This paper

4.373

17.493

39.359

69.971

109.330

Reference

4.3732

17.4928

39.3587

–

–


error (%)

0.004

0.001

0.0007

–

–


$V=\text{15}\text{}\text{m/s}$

This paper

4.287

17.417

39.286

69.900

109.259

Reference

4.2870

17.4171

39.2858

–

–


error (%)

0

0.0006

0.0005

–

–


$V=\text{25}\text{}\text{m/s}$

This paper

4.131

17.282

39.156

69.772

109.133

References

4.1293

17.2816

39.1559

–

–


error (%)

0.04

0.002

0.0003

–

–

It can be found from Table 1 that a well agreement is obtained between the results in this paper and that in reference.
Fig. 4. Curves of $\text{lg}\left(h\right)$frequency under different fluid velocities: a) $V=\text{0}$, b) $V=\text{15}\text{}\text{m/s}$, c) $V=\text{25}\text{}\text{m/s}$
a)
b)
c)
Fig. 4 shows the curve of $\text{lg}h\left(\omega \right)$ versus frequency $\omega $ under different fluid velocities, and the natural frequencies have been marked in the figure. It should be noted here that when the slope of the curve is negative infinite, the real part and the imaginary part of $h\left(\omega \right)$ are both zeros. In another word, the frequencies related to the negative cuspidal points are the natural frequencies of the fluidconveying pipe.
4.2. Variable diameter pipe conveying fluid
To simply supported variable diameter pipes conveying fluid, we choose three different diameter variation coefficients, that is 0.1, 0.2, 0.3, respectively. Given that the Yang’s modulus of the pipe material is $E=\text{70}\text{}\text{GPa}$, the outer and inner diameters of pipe at the origin of the pipe are ${D}_{i}=\text{80}\text{}\text{mm}$ and ${d}_{i}=\text{72}\text{}\text{mm}$ respectively. The pipe length is $L=\text{15}\text{}\text{m}$, the density of pipe material is ${\rho}_{p}=\text{2800}\text{}{\text{kg/m}}^{\text{3}}$, and the density of fluid is ${V}_{0}=\text{25}\text{}\text{m/s}$. The natural frequencies are computed under the following three different fluid velocities: $V=\text{0,}$$V=\text{15}\text{}\text{m/s,}$$V=\text{25}\text{}\text{m/s}$. The results are listed in Table 2.
Table 2. Natural frequencies of variable diameter pipe conveying fluid with different diameter variation coefficients under different fluid velocities
$\gamma $ and fluid velocities

Natural frequencies


${\omega}_{1}$

${\omega}_{2}$

${\omega}_{3}$

${\omega}_{4}$

${\omega}_{5}$


$\gamma =\text{0.1}$

${V}_{0}=\text{0}$

3.54

14.14

31.82

56.56

88.36

${V}_{0}=\text{10}\text{}\text{m/s}$

3.10

13.80

31.52

56.26

88.08


${V}_{0}=\text{25}\text{}\text{m/s}$

–

11.92

29.84

54.68

86.54


$\gamma =\text{0.}2$

${V}_{0}=\text{0}$

3.34

13.48

30.32

53.88

84.18

${V}_{0}=\text{10}\text{}\text{m/s}$

2.88

13.12

29.98

53.56

83.88


${V}_{0}=\text{25}\text{}\text{m/s}$

–

11.10

28.20

51.88

82.84


$\gamma =\text{0.3}$

${V}_{0}=\text{0}$

3.16

12.92

29.04

51.62

80.58

${V}_{0}=\text{10}\text{}\text{m/s}$

2.66

12.52

28.68

51.28

80.42


${V}_{0}=\text{25}\text{}\text{m/s}$

–

10.30

26.76

49.44

77.98

From Table 2, we can find that to the same diameter variation coefficient $\gamma $, the natural frequencies decrease as the fluid velocity increases, and that the larger $\gamma $ we choose, the smaller the natural frequencies we get. The results obtained above agree with that in the researches of Hannoyer and Paidoussis [8].
Additionally, we notice that when the fluid velocity reaches 25 m/s, the first order natural frequency has already vanished. So we can conclude that such velocity has already exceeded the critical velocity of the pipe. Actually, the critical velocities of variable diameter pipe conveying fluid with different diameter variation coefficients are $\gamma =$0.1, ${V}_{cr}=\text{19.96}\text{}\text{m/s;}$$\gamma =$0.2, ${V}_{cr}=\text{18.82m/s}$; $\gamma =\text{}$0.3, ${V}_{cr}=\text{17.77m/s}$, respectively, as shown in Fig. 5.
Fig. 5. Curves of first order natural frequency versus fluid velocity
Fig. 6. Buckling of the variable diameter pipe conveying fluid
When the fluid velocity reaches the critical value, the simply supported variable diameter pipe would yield buckling, as a form of static instability, as shown in Fig. 6. Further, the instability would result in the failure of the pipe.
It is obvious that the larger the diameter variation coefficient is, the smaller the critical velocity we get, which agree with the result in the researches of Hannoyer and Paidoussis [8].
Fig. 7 shows the curves of $\text{lg}h\left(\omega \right)$ versus frequency $\omega $ for variable diameter pipe conveying fluid with different diameter variation coefficients under velocity ${V}_{0}=\text{0}$ and ${V}_{0}=\text{25}\text{}\text{m/s}$.
Fig. 7. Curves of $\text{log}\text{(}h)$frequency with different diameter variation ratio under different fluid velocities: a) $\gamma =\text{0.1}$, ${V}_{0}=\text{0}$; b) $\gamma =\text{0.1}$, ${V}_{0}=\text{25}\text{}\text{m/s}$; c) $\gamma =\text{0.2}$, ${V}_{0}=\text{0}$; d) $\gamma =\text{0.2}$, ${V}_{0}=\text{25}\text{}\text{m/s}$; e) $\gamma =\text{0.3}$, ${V}_{0}=\text{0}$; f) $\gamma =\text{0.3}$, ${V}_{0}=\text{25}\text{}\text{m/s}$
a)
b)
c)
d)
e)
f)
After we obtain the natural frequencies of the variable diameter pipe, the modal shapes can be obtained correspondingly. The modal shapes of variable diameter pipe conveying fluid are shown in Fig. 8 and Table 3.
Table 3. Dimensionless amplitudes of the modal shape under different flow velocities
Fluid velocity

Amplitudes


1st order

2nd order

3rd order

4th order


${V}_{0}=\text{0}$

0.22

0.29

0.21

0.12

${V}_{0}=\text{10}\text{}\text{m/s}$

0.25

0.31

0.23

0.15

${V}_{0}=\text{25}\text{}\text{m/s}$

0.30

0.35

0.26

0.18

Table 3 shows that the amplitude of the modal shape increases with increasing flow velocity. This illuminate that the increasing fluid velocity can weaken the pipe stiffness.
Fig. 8. The first four orders modal shapes of variable diameter pipe conveying fluid with different diameter variation coefficients under different fluid velocities
4.3. Variable diameter pipe with damping
Damping is ignored in above calculations. But in practice, damping cannot be neglected. Here we just consider the viscous damping, then the motion equation governing the vibration of the pipe is:
where the term ($c\dot{w}$) is the viscous damping force, which may be introduced by immersing the pipe in liquid. Suppose that $c=$0.05, then we can calculate the natural frequencies of the variable diameter pipe. A contrast is illustrated in Fig. 9 between damping pipe and nondamping pipe ($c=$0).
Fig. 9. First order natural frequencies of damping pipe and nondamping pipe with variable fluid velocity
From Fig. 9, we conclude that damping would decrease the natural frequencies and critical flow velocity of the system. This has been illustrated by Paidoussis [2] in the research for a uniform fluidconveying pipe. Moreover, damping delays the decreasing speed of natural frequency versus flow velocity. That is, as seen in Fig. 9, there is just a small difference between the two critical flow velocities, but a relatively larger difference between the natural frequencies when flow velocity is zero.
5. Conclusion
In this paper, the dynamic model was built by open system’s Hamilton’s principle on the basis of EulerBernoulli beam model and axial inextensible assumption. The motion equation of variable diameter pipe conveying fluid is more complex than that of uniform pipe. We employed Frobenius method to analyze the motion equation and proposed the dynamic stiffness method for free vibration analysis of variable diameter pipe conveying fluid. By using the presented method, the natural frequencies of a uniform pipe conveying fluid under different fluid velocities were obtained and the comparisons between our results and that in previous researches were made. A well agreement validates our method. The natural frequencies and modal shapes were obtained for variable diameter pipe conveying fluid with different diameter variation coefficients under different fluid velocities. We find that the natural frequencies and critical velocities of variable diameter pipe both decreased with increasing diameter variation coefficient.
Acknowledgements
This work was supported by Northwestern Polytechnical University Basic Research Fund (Grant No. JC201114) and Aeronautical Science Foundation of China (2011ZA53014).
References
 Ashley H., HavilandG. Bending vibrations of a pipe line containing flowing fluid. Journal of Applied MechanicsTransactions of the ASME, Vol. 17, Issue 3, 1950, p. 229232. [Search CrossRef]
 Paidoussis M. P. Fluid structure interactions: slender structures and axial flow. Volume 1. London, Academic Press, 1998. [Search CrossRef]
 Doare O. Dissipation effect on local and global stability of fluidconveying pipes. Journal of Sound and Vibration, Vol. 329, Issue 1, 2010, p. 7283. [Search CrossRef]
 Jung D., Chung J., Mazzoleni A. Dynamic stability of a semicircular pipe conveying harmonically oscillating fluid. Journal of Sound and Vibration, Vol. 315, Issue 12, 2008, p. 100117. [Search CrossRef]
 Rinaldi S., et al. Dynamics of microscale pipes containing internal fluid flow: damping, frequency shift, and stability. Journal of Sound and Vibration, Vol. 329, Issue 8, 2010, p. 10811088. [Search CrossRef]
 Huang Y. M., et al. Natural frequency analysis of ﬂuid conveying pipeline with different boundary conditions. Nuclear Engineering and Design, Vol. 240, Issue 3, 2010, p. 461467. [Search CrossRef]
 Ragulskis M., Fedaravicius A., Ragulskis K. Harmonic balance method for FEM analysis of fluid flow in a vibrating pipe. Communications in Numerical Methods in Engineering, Vol. 22, Issue 4, 2006, p. 347356. [Search CrossRef]
 Hannoyer M. J., Paidoussis M. P. Dynamics of slender tapered beams with internal or external axial flow 1 theory. Journal of Applied MechanicsTransactions of the ASME, Vol. 46, Issue 1, 1979, p. 4551. [Search CrossRef]
 Hannoyer M. J., Paidoussis M. P. Dynamics of slender tapered beams with internal or external axial flow 2 experiments. Journal of Applied MechanicsTransactions of the ASME, Vol. 46, Issue 1, 1979, p. 5257. [Search CrossRef]
 Arfken G. B., Weber H. J. Mathematical methods for physicists. San Diego, CA, Academic Press, 1995. [Search CrossRef]
 Housner G. W. Bending vibrations of a pipe line containing flowing fluid. Journal of Applied MechanicsTransactions of the ASME, Vol. 19, Issue 2, 1952, p. 205208. [Search CrossRef]