Published: 30 September 2015

# Measurement of stability of a pipe system with flowing fluid

A. Sudintas1
P. Paškevičius2
B. Spruogis3
1, 2Kaunas University of Technology, Kaunas, Lithuania
3, 4Vilnius Gediminas Technical University, Vilnius, Lithuania
Corresponding Author:
A. Sudintas
Views 21

#### Abstract

The two dimensional system of coupled vibrating pipes is investigated. The model of a pipe system consisting from two dimensional beams is employed. The problem of stability because of the reduction of stiffness caused by the flow of fluid is solved and the stability eigenmodes are determined. The places where the deflections of the stability eigenmodes are large are recommended for the location of measurement devices.

## 1. Introduction

The two dimensional system of coupled vibrating pipes is investigated. The model of a pipe system consisting from two dimensional beams is used. The problem of stability because of the reduction of stiffness caused by the flow of fluid is solved and the stability eigenmodes are determined.

Measurement of vibrations of similar systems when the fluid does not flow was investigated in [1, 2]. In  vibrations in the plane of the two dimensional system of pipes were investigated and in  transverse vibrations of the two dimensional system of pipes were investigated. This paper continues the investigations presented in both previous papers.

The model for the analysis of stability of a pipe system is proposed on the basis of the results described in [1, 2] and general relationships presented in [3-5]. Similar problems of dynamics, vibrations and stability are investigated in [6-12].

## 2. In plane stability of a pipe system with flowing fluid

### 2.1. Model for the analysis of stability of a pipe system with flowing fluid

$x$, $y$ and $z$ denote the axes of the system of coordinates. The finite element of a two dimensional beam representing the pipe in the plane $xOz$ has three nodal degrees of freedom: the displacement in the direction of the $x$ axis denoted as $u$, the displacement in the direction of the $z$ axis denoted as $w$ and the rotation about the $y$ axis denoted as ${\theta }_{y}$.

The displacements in the direction of the longitudinal axis of the pipe $s$ and in the direction of the axis perpendicular to $s$ and located in the plane $xOz$ are denoted as $\stackrel{-}{u}$ and $\stackrel{-}{w}$ and they are related with $u$, $w$ as:

1
$\left\{\begin{array}{c}\stackrel{-}{u}\\ \stackrel{-}{w}\\ {\theta }_{y}\end{array}\right\}=\left[T\right]\left\{\begin{array}{c}u\\ w\\ {\theta }_{y}\end{array}\right\}=\left[\begin{array}{ccc}\frac{dx}{d\xi }/\frac{ds}{d\xi }& \frac{dz}{d\xi }/\frac{ds}{d\xi }& 0\\ -\frac{dz}{d\xi }/\frac{ds}{d\xi }& \frac{dx}{d\xi }/\frac{ds}{d\xi }& 0\\ 0& 0& 1\end{array}\right]\left\{\begin{array}{c}u\\ w\\ {\theta }_{y}\end{array}\right\},$

where $\xi$ is the local coordinate of the finite element of the pipe and $ds/d\xi =\sqrt{{\left(dx/d\xi \right)}^{2}+{\left(dz/d\xi \right)}^{2}}.$

The following notation is used:

2
$\left[\begin{array}{l}\left[{\stackrel{-}{N}}_{1}\right]\\ \left[{\stackrel{-}{N}}_{2}\right]\\ \left[{\stackrel{-}{N}}_{3}\right]\end{array}\right]=\left[T\right]\left[\begin{array}{cccc}{N}_{1}& 0& 0& \dots \\ 0& {N}_{1}& 0& \dots \\ 0& 0& {N}_{1}& \dots \end{array}\right],$

where ${N}_{i}$ are the shape functions of the one dimensional finite element of the pipe.

Also it is assumed:

3
$\left[\begin{array}{l}\left[{\stackrel{-}{N}}_{1}^{\text{'}}\right]\\ \left[{\stackrel{-}{N}}_{2}^{\text{'}}\right]\\ \left[{\stackrel{-}{N}}_{2}^{\text{'}}\right]\end{array}\right]=\left[T\right]\left[\begin{array}{cccc}\frac{d{N}_{1}}{ds}& 0& 0& \dots \\ 0& \frac{d{N}_{1}}{ds}& 0& \dots \\ 0& 0& \frac{d{N}_{1}}{ds}& \dots \end{array}\right].$

The stiffness matrix of the pipe has the usual form:

4
$\left[K\right]=\int \left({\left[{\stackrel{-}{N}}_{1}^{\text{'}}\right]}^{T}E\pi \left({R}^{2}-{r}^{2}\right)\left[{\stackrel{-}{N}}_{1}^{\text{'}}\right]+{\left[{\stackrel{-}{N}}_{3}^{\text{'}}\right]}^{T}E\pi \frac{{R}^{4}-{r}^{4}}{4}\left[{\stackrel{-}{N}}_{3}^{\text{'}}\right]$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{\left[B\right]}^{T}\frac{E}{2\left(1+\nu \right)1.2}\pi \left({R}^{2}-{r}^{2}\right)\left[B\right]\right)ds,$

where $E$ is the modulus of elasticity of the pipe, $\nu$ is the Poisson’s ratio of the pipe, $r$ is the internal radius of the pipe, $R$ is the external radius of the pipe and $\left[B\right]=\left[\left[{\stackrel{-}{N}}_{2}^{\text{'}}\right]+\left[{\stackrel{-}{N}}_{3}\right]\right].$

The supplementary stiffness matrix of the model of a pipe because of the motion of the fluid in it has the form:

5
$\left[{K}_{v}\right]=\int {\left[{\stackrel{-}{N}}_{2}^{\text{'}}\right]}^{T}{v}^{2}{\rho }_{f}\pi {r}^{2}\left[{\stackrel{-}{N}}_{2}^{\text{'}}\right]ds,$

where ${\rho }_{f}$ is the density of the fluid filling the pipe and v is the velocity of the fluid inside the pipe.

The total stiffness of the pipe is equal to $\left[K\right]-\left[{K}_{v}\right]$. If the velocity is changed $\sqrt{\lambda }$ times, then the stiffness becomes $\left[K\right]-\lambda \left[{K}_{v}\right]$. This results in the eigenproblem for the determination of critical velocities and stability eigenmodes.

### 2.2. Analysis and measurement of stability of a pipe system with flowing fluid

The structure consists from a lower straight pipe and an upper straight pipe both of them having fixed ends (all of the generalized displacements equal to zero) and they are assumed filled with water. The midpoints of both pipes are connected by an empty pipe perpendicular to them. Length of the structure is 2 m and the distance between the lower and upper pipes is 0.4 m.

The following parameters are assumed: modulus of elasticity of the pipes $E=$ 6·108 Pa, Poisson’s ratio of the pipes $\nu =$ 0.3, density of the fluid filling the two pipes ${\rho }_{f}=$ 998 kg/m3, velocity of the fluid in the lower pipe is assumed to be equal to $v=$ 1 m/s, velocity of the fluid in the upper pipe is assumed to be equal to $v=$ –1 m/s, internal radius of the pipes $r=$ 0.004 m, external radius of the pipes $R=$ 0.006 m.

The first stability eigenmodes are presented in Fig. 1. The fourth and the fifth eigenmodes correspond to the same eigenvalue.

In order to perform precise measurements of stability of a pipe system the measurements of deflections of the pipes are to be performed not at the nodes of the stability eigenmodes and preferably at the places with maximum deflections. In practical applications usually the first eigenmode is important. Thus from the presented images of deflections the locations of measurement devices may be chosen.

Fig. 1The first stability eigenmodes of the pipe system: a) the first eigenmode, b) the second eigenmode, …, j) the tenth eigenmode a) b) c) d) e) f) g) h) i) j)

## 3. Transverse stability of a pipe system with flowing fluid

### 3.1. Model for the analysis of transverse stability of a pipe system with flowing fluid

The finite element of a two dimensional beam representing the pipe in the plane $xOy$ has three nodal degrees of freedom: the displacement in the direction of the $z$ axis denoted as $w$, the rotation about the $x$ axis denoted as ${\theta }_{x}$ and the rotation about the $y$ axis denoted as ${\theta }_{y}$.

The rotations about the longitudinal axis of the pipe $s$ and about the axis $t$ perpendicular to $s$ and located in the plane $xOy$ are denoted as ${\theta }_{s}$ and ${\theta }_{t}$ and they are related with ${\theta }_{x}$, ${\theta }_{y}$ as:

6
$\left\{\begin{array}{c}w\\ {\theta }_{s}\\ {\theta }_{t}\end{array}\right\}=\left[T\right]\left\{\begin{array}{c}w\\ {\theta }_{x}\\ {\theta }_{y}\end{array}\right\}=\left[\begin{array}{lll}1& 0& 0\\ 0& \frac{dx}{d\xi }/\frac{ds}{d\xi }& \frac{dy}{d\xi }/\frac{ds}{d\xi }\\ 0& -\frac{dy}{d\xi }/\frac{ds}{d\xi }& \frac{dx}{d\xi }/\frac{ds}{d\xi }\end{array}\right]\left\{\begin{array}{c}w\\ {\theta }_{x}\\ {\theta }_{y}\end{array}\right\},$

where $ds/d\xi =\sqrt{{\left(dx/d\xi \right)}^{2}+{\left(dy/d\xi \right)}^{2}}.$

The stiffness matrix of the pipe for the analysis of transverse stability has the usual form:

7
$\left[K\right]=\int \left({\left[{\stackrel{-}{N}}_{3}^{\text{'}}\right]}^{T}E\pi \frac{{R}^{4}-{r}^{4}}{4}\left[{\stackrel{-}{N}}_{3}^{\text{'}}\right]+{\left[B\right]}^{T}\frac{E}{2\left(1+\nu \right)1.2}\pi \left({R}^{2}-{r}^{2}\right)\left[B\right]$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{\left[{\stackrel{-}{N}}_{2}^{\text{'}}\right]}^{T}\frac{E}{2\left(1+\nu \right)}\pi \frac{{R}^{4}-{r}^{4}}{2}\left[{\stackrel{-}{N}}_{2}^{\text{'}}\right]\right)ds,$

where $\left[B\right]=\left[\left[{\stackrel{-}{N}}_{1}^{\text{'}}\right]+\left[{\stackrel{-}{N}}_{3}\right]\right].$

The supplementary stiffness matrix of the model of a pipe for the analysis of transverse stability because of the motion of the fluid in it has the form:

8
$\left[{K}_{v}\right]=\int {\left[{\stackrel{-}{N}}_{1}^{\text{'}}\right]}^{T}{v}^{2}{\rho }_{f}\pi {r}^{2}\left[{\stackrel{-}{N}}_{1}^{\text{'}}\right]ds.$

### 3.2. Analysis and measurement of transverse stability of a pipe system with flowing fluid

The displacement in the direction of the $z$ axis is represented in the normal direction to the pipe in the plane $xOy$. The first transverse stability eigenmodes are presented in Fig. 2. The fifth and the sixth eigenmodes correspond to the same eigenvalue.

From the presented images of deflections the locations of measurement devices may be chosen.

Fig. 2The first transverse stability eigenmodes of the pipe system: a) the first eigenmode, b) the second eigenmode, …, j) the tenth eigenmode a) b) c) d) e) f) g) h) i) j)

## 4. Conclusions

The structure consisting from a lower straight pipe and an upper straight pipe both of them having fixed ends (all of the generalized displacements equal to zero) and filled with water is investigated. The midpoints of both pipes are connected by an empty pipe perpendicular to them. It is assumed that the fluid in the two pipes is flowing with given velocities. The first stability eigenmodes are obtained.

In order to perform precise measurements of stability of a pipe system the measurements are to be performed not at the nodes of the stability eigenmodes and preferably at the places with maximum deflections. The acceptable places of measurement are seen from the images of the stability eigenmodes.