# Dynamic Response of Pipe Conveying Fluid with Lateral Moving Supports

Authors: Baohui Li; Zhengzhong Wang; Lina Jing
Journal: Shock and Vibration
 (2018)

---

## Abstract

The Galerkin method is proposed to reveal the dynamic response of pipe conveying fluid (PCF), with lateral moving supports on both ends of the pipe. Firstly, the dynamic equation is derived by the Newtonian method after calculating the acceleration of the fluid element via the dynamics approach. Secondly, the discrete form of the dynamic equation is formulated by the Galerkin method. Thirdly, the numerical analysis of the system is carried out through the fourth-order Runge–Kutta method, and the effectiveness of the proposed method is validated by comparison with the analytical results obtained by the mode superposition method. In the example analysis, the responses of the lateral deflection and bending moment are investigated for the pinned-pinned, clamped-pinned, and clamped-clamped PCF. The effects of fluid velocity and the moving frequencies of supports are discussed. Especially, the deflection responses are analyzed under extreme condition; i.e., the moving frequency of a support is identical to the natural frequency of PCF.


---

## Body

## 1. Introduction

Transportation pipes are widely used in various industry fields, such as oil exploitation, aviation fuel piping, and nuclear plant cooling systems. The flow-induced vibration, however, seriously threatens the safe operation of the systems. Each year, the vibration-induced leakage and failure cause huge losses in the economy worldwide. Due to this reason, studies on pipe dynamics have become popular in recent years, and a new dynamical model gradually formed [1–4].Usually, the dynamic equation of PCF can be established via the Newtonian method or generalized Hamilton’s principle for an open system [5, 6]. Several dynamic models have been developed during researching PCF. To date, the researchers have already reached a consensus on the classic dynamic model of pipe, such as the linear free vibration model for PCF, which can be described by the following equation [2]:(1)EI∂4w∂x4+mfV2∂2w∂x2+2mfV∂2w∂x∂t+mf+mp∂2w∂t2=0.At present, three aspects, i.e., the numerical solution techniques, pipe dynamic characters, and the extended dynamic model for more specific problems, attract much attention. For instance, Sreejith et al. [7] proposed the finite-element method for solving the dynamic responses of the pipeline system in nuclear reactors. Olson and Jamison [8] used the nonlinear finite-element method to predict the response of elastic PCF. Lee and Oh [9] developed a spectral element model of the pipeline. The method provides more accuracy than the classical finite-element model. Wang and Bloom [10] developed a spatial finite-difference scheme for the stability analysis of submerged and inclined concentric PCF. Naguleswaran and Williams [11] adopted Rayleigh–Ritz and Fourier series solutions to study the lateral vibration of several common constrained PCF with consideration of axial tension and internal pressure. The dynamic properties of PCF are also investigated considering different conditions. Tang et al. [12] applied the multiple scales method to study the fractional dynamics of fluid-conveying pipes made of polymer-like materials. Mnassri and Baroudi [13] developed a theoretical method for the three-dimensional modal analysis of compressible fluid within pipes. The method can be used for defect detection. Zhou et al. [14] studied axially functional cantilevered PCF and analyzed the effects of the elastic modulus gradient and density gradient on the critical flow velocity for flutter instability. Oke and Khulief [15] discussed the vibration behavior of a composite PCF having internal surface damage.The mechanisms discovered in PCF have been gradually applied to other fields. By studying the lift effects of fluid force of a cantilevered elastic pipe for conveying fluid, Texier and Dorbolo [16] tried to introduce motion in soft robotics. The Coriolis effect in a Coriolis mass flow meter under lateral vibration was simulated by the finite-element method based on Timoshenko beam theory [17]. Recently, the extended dynamic models were applied to study the axially moving pipe [18–20] and helical tube conveying fluid [21]. Applications of the theory of PCF are not only limited to traditional pipes but also limited to micropipes and nanotubes conveying fluid [22–27] and even to the novel pipe composed of two different materials [28].In practice, a pipeline usually works in the vibration environment. For instance, the aviation fuel pipeline lays on an aircraft engine, airborne pipes, and other fluid machinery equipment. In such a situation, the pipe supports are in the moving state. It may lead to constraints loosening, pipe fatigue failure, and/or leakage under such circumstances. In this study, the characteristics of dynamic response of PCF are investigated with consideration of both coupling effect between fluid and pipe and the external moving supports. Firstly, the dynamic equation for the pipe constrained by moving supports is deduced by the Newtonian method after finding the acceleration of fluid elements in Section1. Secondly, by virtue of the Galerkin method, in Section 2, the discrete dynamic model is built and the fourth-order Runge–Kutta method is employed for numerical analysis. Finally, the influences of fluid velocity and the movement frequency of supports on the dynamic response of pipe are investigated in Section 3. The pipes with common boundaries, such as pinned-pinned, clamped-pinned, and clamped-clamped pipes, are discussed, respectively. In Section 4, some conclusions are drawn.


## 2. Dynamic Equation of PCF with Lateral Moving Supports

The diagram of PCF with lateral moving supports is shown in Figure1. The global inertial reference frame is XOY with unit vectors I and J. The moving reference frame xoy is attached to the supports, with unit vectors i and j, lateral velocity v0, and lateral acceleration a0. The lateral deflection w is measured in the local frame xoy.Figure 1
PCF with lateral moving supports.
(a)
(b)The acceleration of the fluid element, shown in Figure2, is essential for derivation of the dynamic equation of PCF by the Newtonian method. Here, we introduce another moving reference frame x′o′y′ with unit vectors i′andj′, which is attached to the pipe element. Parameters φ and ρ denote the slope angle of the pipe cross section and the curvature radius of the deformed pipe element, respectively. Without considering variation along radial direction, the acceleration of the fluid element can be represented by the acceleration of the fluid particle at the point O′.Figure 2
Fluid and pipe element.Using the composition theorem in dynamics [29], the acceleration of the fluid element can be written as follows:(2)af=ar+ae+2ω×V,where ar is the relative acceleration, i.e., the fluid acceleration relative to the reference frame x′o′y′, and convected acceleration ae is the acceleration of the point O′ fixed on the reference frame x′o′y′ relative to the inertial reference frame XOY. ω denotes the angular velocity of the pipe cross section, that is ω=dφ/dtk′, in which k′ is the unit vector perpendicular to plane XOY outwardly. V is the velocity of fluid in the pipe, that is, V=Vi′.The terms in Equation (2) can be expressed as follows [1]:(3)ar=dVdti′−V2ρj′,ae=∂2δ∂t2i+a0+∂2w∂t2j≃a0+∂2w∂t2j,2ω×V=2Vdφdtk′×i′=2V∂2w∂s∂tj′,where δ is the longitudinal displacement that is one order smaller than the lateral deflection w and can be neglected.Substituting Equation (3) and 1/ρ=−∂2w/∂s2/1+∂w/∂s23/2≈−∂2w/∂s2 into Equation (2), we have(4)af=dVdti′+2V∂2w∂s∂t+V2∂2w∂s2j′+a0+∂2w∂t2j.Equation (4) can be further simplified once the transformation relationship between unit vectors i′,j′ and i,j is given. As the pipe has small deformation, the transformation can be expressed as follows:(5)i′j′=cosφsinφ−sinφcosφij≈1∂w∂s−∂w∂s1ij.Substituting Equation (5) into Equation (4) and neglecting the second/higher order of infinitesimal, the expression of fluid acceleration yields(6)af=dVdti+dVdt∂w∂s+2V∂2w∂s∂t+V2∂2w∂s2+∂2w∂t2+a0j.If the pipe contains steady fluid and is nonextendable along axial direction, i.e.,s ≈ x and dV/dt = 0, Equation (6) can be simplified further as follows:(7)af=2V∂2w∂x∂t+V2∂2w∂x2+∂2w∂t2+a0j.After obtaining the fluid acceleration, the dynamic equation can be deduced by the Newtonian method. The free body diagram of a pipe element is shown in Figure3, where Q and M are the transverse shear force and bending moment. The pipe contains fluid with pressure p¯ and flow area A. mf and mp are the masses of fluid and pipe per length. The inertia force on the fluid and pipe elements is labeled as FIf and FIp, respectively. Hence, the lateral vibration equation of PCF can be formulated as follows:(8)∂Q∂x−p¯A∂2w∂x2=FIf+FIp=mfaf+mpa0+∂2w∂t2.Figure 3
Free body diagram of a pipe element.Substituting Equation (7) into Equation (8) and using the relationship between the shear force and bending moment, i.e., ∂Q/∂x=−∂2M/∂x2=−∂2/∂x2EI∂2w/∂x2, the lateral vibration equation is finally expressed as follows:(9)EI∂4w∂x4+mfV2+p¯A∂2w∂x2+2mfV∂2w∂x∂t+mf+mp∂2w∂t2=−mf+mpa0.When the supports vibrate harmonically, e.g.,A0sinω0t, the acceleration reads a0=−A0ω02sinω0t, and the dynamic equation becomes(10)EI∂4w∂x4+mfV2+p¯A∂2w∂x2+2mfV∂2w∂x∂t+mf+mp∂2w∂t2=mf+mpA0ω02sinω0t.For convenience, the following dimensionless parameters are introduced to obtain the dimensionless form of Equation (10):(11)η=wL,ξ=xL,u=MEI1/2LV,P=p¯AEIL2,β=mfmf+mp1/2,τ=EImf+mp1/2tL2,ω=mf+mpEI1/2ω0L2,γ=A0L,where L is the length of the pipe.Finally, the dimensionless form of Equation (10) becomes(12)∂4η∂ξ4+u2+P∂2η∂ξ2+2β1/2u∂2η∂ξ∂τ+∂2η∂τ2=γω2sinωτ.


## 3. Galerkin Discretization

The Galerkin method is applied to deal with the partial differential equation, and the dimensionless lateral deflection is represented as follows:(13)ηξ,τ≃∑i=1Nϕiξqiτ,where ϕiξ are the comparison functions and satisfy the boundary condition of the pipe and qiτ are the generalized coordinates of the discretized pipe structure.Substituting Equation (13) back into Equation (12) and minimizing the residual value, one will obtain the following discrete form of dynamic equation:(14)Mq¨+Cq˙+Kq=F,where M, C, K, and F are the matrices of mass, damping, stiffness, and external load, respectively. Meanwhile, q˙ and q¨ are the generalized velocity and acceleration.The specific expressions of the elements of above matrices are as follows:(15)Mij=∫01ϕiξϕjξdξ,Cij=2β1/2u∫01ϕi′ξϕjξdξ,Kij=∫01ϕi″″ξ+u2+Pϕi″ξϕjξdξ,Fj=∫01γω2sinωτϕjξdξ,where ′=∂/∂x.The fourth-order Runge–Kutta method is employed to solve Equation (14) by changing the equation into the following form:(16)Y˙τ=AYτ+Pτ,where Yτ=qτq˙τ,A=0N×NIN×N−M−1K−M−1C,andPτ=0N×1M−1F. IN×N is the N-ordered identity matrix.In numerical analysis, the dimensionless time step is set asΔτ=0.002 and the initial conditions are qτ=0andq˙τ=0. According to Paidoussis and Semler [30], Wang et al. [31], and Ni et al. [20], the convergence can be guaranteed when the order of the Galerkin truncation error of Equation (13) satisfies N≥4. In this study, we choose N=8 for the following analysis.


## 4. Examples

In the examples, the following parameters are chosen for PCF: elastic modulus of the pipe materialE = 68 GPa, the Poisson ratio μ = 0.33, and the density ρp = 2700 kg·m−3. The outer and inner diameters of the pipe are D = 46 mm and d = 40 mm, respectively. The pipe length is L = 2 m. The fluid density is ρf = 870 kg·m−3, and the mean inner pressure p¯=2MPa. The supports are moved harmonically, i.e., A0 sin(ω0t), with the amplitude of A0 = 0.005 m and angular frequency of ω0=300πrads−1.


### 4.1. The Pinned-Pinned PCF

In beam theory, the comparison function for the pinned-pinned pipe, shown in Figure4, is ϕiξ=siniπξ. To show the effectiveness of the proposed method, we compare the results obtained by both the analytical mode superposition method for classical beam problems and the proposed method. For a solid beam, it can be considered as a PCF with zero velocity of fluid (u = 0) and without inner pressure (p¯=0). Comparison of the results will be carried out when u = 0 and p¯=0. It is easy to find out Equation (12) draws the state of a pipe with stationary supports and subjected to uniform load pξ,τ=γω2sinωτ. In this case, based on the mode superposition method and beam theory [32], the lateral deflection of the pipe can be expressed as follows:(17)ηξ,τ=∑iϕiξqiτ.Figure 4
Pinned-pinned PCF.Considering the uniform external loadpξ,τ=γω2sinωτ and initial condition, i.e., ηξ,0=η˙ξ,0=0, the generalized coordinate can be expressed analytically [32] as follows:(18)qiτ=1−−1i1ωi2γω2iπ∫0τsinωζsinωiτ−ζdζ=2γω2iπ1−−1iωi2−ω2sinωτ−ωωisinωiτ,where ωi is the ith natural frequency of the pinned-pinned pipe with stationary fluid and zero inner pressure.Now, the analytical solution for the pipe deflection can be obtained by substituting Equation (18) into Equation (17). We find that the results obtained by the mode superposition method are stable when i ≥ 3. Hence, we set i = 4 for the comparison. The results calculated by the two methods (Figure 5) match each other very well. Hence, the validation of the present method is verified.Figure 5
Comparison of deflection of the pinned-pinned pipe by two methods.We choose three different truncations of Equation (13), with N = 4, 6, and 8, in finding the responses of the pinned-pinned pipe. As demonstrated in Figure 6, the three results are almost the same and converge to the analytical solutions as long as N≥4.Figure 6
Comparison of deflections of the pinned-pinned pipe with differentN.The lateral deflections of the pinned-pinned pipe with three different fluid velocities (u = 0, u = 1, and u = 2) at three different moments (τ = 10, τ = 20, and τ = 30) are involved in discussion in this paper. The detailed results are illustrated in Figure 7.Figure 7
The deflections of the pinned-pinned pipe with different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)The three modes of deflection curves in Figure7 are different in each case. Meanwhile, the deflection curves in Figures 7(b) and 7(c) are asymmetric due to flow of fluid within the beam. It is easy to find out in Figure 7 that the deflections gradually increase with fluid velocity. This implies that the rigidity of the pipe becomes lower when the fluid speed is higher. The effect also appears on the natural frequencies of the pipe with different fluid velocities. In this paper, the first four orders of natural frequency are calculated through the wave approach which was introduced by Li et al. [33, 34]. The results are listed in Table 1.Table 1
The dimensionless natural frequencies of the pinned-pinned PCF.
u
ω
ω
1
ω
2
ω
3
ω
4
0
9.05
38.68
88.04
157.13
1
8.44
38.22
87.61
156.71
2
6.37
36.81
86.31
155.46In practice, when the external vibration frequency is in coincidence with the natural frequencies of structures, then structural resonance will happen, and large deformation will be generated. To evaluate such an effect on pipe structures, we investigate the deflections of PCF with supports moving in the same frequencies as its natural frequencies.For simplicity, the deflections of the pinned-pinned pipe with fluid velocityu = 1 are analyzed, and three cases with respect to the moving frequency of the supports, i.e., ω = 8.44, 38.22, and 87.61, are considered. Figure 8 shows the resonance responses in the three cases.Figure 8
The resonance responses of the pinned-pinned PCF withu = 1 under different supports moving frequencies: (a) ω = 8.44; (b) ω = 38.22; (c) ω = 87.61.
(a)
(b)
(c)Compared with the results in Figure7, the amplitudes of deflections in Figure 8 become much larger, especially when the supports moving frequency reaches the third natural frequency of PCF. At the times τ = 20 and τ = 30 in Figure 8(c), the largest amplitude of deflection appears near the right end of the pipe rather than at the center. The same phenomenon can be found in Figure 7(c).We also obtain the history curves of pipe deflection at different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) under different fluid velocities (u = 0, u = 1, and u = 2) (Figure 9).Figure 9
History curves of deflection at different locations on the pinned-pinned pipe having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)Figure9 suggests that the amplitude of deflection increases with the fluid velocity. The maximum deflection occurs at the center of the pipe. Deflection of a point becomes smaller when the point is closer to a pipe end. With the increase of fluid velocity, the response curves are gradually getting sparse. It implies that the response period becomes larger; i.e., the natural frequency is getting smaller.Bending moment is another important dynamic response for the beam. Once it is obtained, the normal stress on the pipe cross section will be determined. Furthermore, the fatigue analysis can be conducted. From Equation (13), the expression of dimensionless bending moment becomes(19)M¯=MLEI=∂2η∂ξ2.Figure10 illustrates the history curves of bending moment on the pinned-pinned pipe at three different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) with respect to three different fluid velocities (u = 0, u = 1, and u = 2).Figure 10
History curves of bending moment on the pinned-pinned pipe at different locations and having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)In Figure10(a), when u = 0, the bending moment at ξ = 1/2 is larger than those at ξ = 1/4 and ξ = 1/8, but the bending moments at ξ = 1/4 and ξ = 1/8 are different slightly. As the fluid velocity increases, the bending moments at the three locations increase, simultaneously. However, at ξ = 1/4 and ξ = 1/8, they tend to have a faster increase than that at location ξ = 1/2.


### 4.2. The Clamped-Pinned PCF

Another type of fluid-conveying pipe is the clamped-pinned pipe in Figure11. Its comparison function [35] becomes(20)ϕiξ=sinβiξ−sinhβiξ+αicoshβiξ−cosβiξ,where αi=sinβil−sinhβil/cosβil−coshβil and βil≈i+1/4π when i ≥ 1. The values of βil are listed in Table 2 [35]. At three different times (τ = 10, τ = 20, and τ = 30), the deflection responses of the pipe having three different fluid velocities (u = 0, u = 1, and u = 2) are illustrated in Figure 12.Figure 11
The model for a clamped-pinned PCF.Table 2
The values ofβil for clamped-pinned PCF.
β
1
l
β
2
l
β
3
l
β
4
l
3.9266
7.0686
10.2102
13.3518Figure 12
The deflection responses of the clamped-pinned pipe under different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)There is no symmetric mode in any case of Figure12 due to the asymmetric supports. Deflection of the pipe does not vary monotonously with the fluid velocity. At time τ = 10, with the increase of fluid velocity, the deflection first increases and then turns to decrease. However, at time τ = 20, the deflection always decreases. At time τ = 30, the deflection firstly decreases and then increases.The first four natural frequencies of the clamped-pinned PCF are calculated by the wave approach [33, 34] (Table 3). The resonance analysis with respect to fluid velocity u = 1 is also carried out when the supports moving frequencies ωi are equal to 14.36, 48.89, and 103.15, respectively.Table 3
The first four orders of dimensionless natural frequencies of the clamped-pinned PCF.
u
ω
ω
1
ω
2
ω
3
ω
4
0
14.82
49.29
103.54
177.55
1
14.36
48.89
103.15
177.16
2
12.94
47.69
101.98
176.00The deflection curves in Figures13(a) and 13(c) are much larger than that in Figure 13(b), and the largest deflection appears in Figure 13(c). In other words, the resonance is very strong once the supports moving frequency equals the first or the third natural frequency of the pipe. The largest deflection occurs at the right side of the pipe center in Figure 13(a), whilst at the left side of the pipe center in Figures 13(b) and 13(c).Figure 13
The resonance responses of the clamped-pinned pipe withu = 1 and different supports moving frequencies: (a) ω = 14.36; (b) ω = 48.89; (c) ω = 103.15.
(a)
(b)
(c)The history curves of the deflection of the clamped-pinned pipe conveying fluid at different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) and having different fluid velocities (u = 0, u = 1, and u = 2) are illustrated in Figure 14.Figure 14
History curves of the deflection of the clamped-pinned pipe at different locations and having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)It can also be observed that the deflections increase with the fluid velocity. Compared with the pinned-pinned pipe, however, the clamped-pined pipe has lower variation of deflection due to larger support stiffness at the clamped end of the clamped-pinned pipe, and the same fluid velocity produces smaller effect on the pipe deflection.The history curves of bending moment on the clamped-pinned pipe in different cases are given in Figure15. Figure 15 indicates that, in each case, the largest bending moment on the clamped-pinned pipe appears at ξ = 1/2, and this is the same as that of the pinned-pinned pipe. The bending moment at ξ = 1/8, however, is larger than that at ξ = 1/4, and this is different from that of the pinned-pinned pipe. And the magnitude of bending moment on the clamped-pinned pipe is much lower than that of the pinned-pinned pipe.Figure 15
History curves of bending moment response for the clamped-pinned pipe at different locations and under different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)


### 4.3. The Clamped-Clamped PCF

The dynamic model for the clamped-clamped pipe, shown in Figure16, is different from that for either the pinned-pinned pipe or the clamped-pinned pipe because the inextensible assumption is not suitable anymore. It is easy to find from Figure 16 that while vibration happens, the axial extension is inevitably permitted as the beam has transverse vibration. If so, an additional tensile force, which equals T=2μp¯Af, should be induced on the cross section of the pipe [1]. Hence, the dynamic equation of the clamped-clamped PCF needs to be modified as follows:(21)EI∂4w∂x4+mfV2+p¯A1−2μ∂2w∂x2+2mfV∂2w∂x∂t+mf+mp∂2w∂t2=mf+mpA0ω02sinω0t.Figure 16
Schematic of a clamped-clamped PCF.Correspondingly, the dimensionless form of Equation (21) can be written as follows:(22)∂4η∂ξ4+u2+P1−2μ∂2η∂ξ2+2β1/2u∂2η∂ξ∂τ+∂2η∂τ2=γω2sinωτ.In the discrete form of Equation (14), only the stiffness matrix needs to be changed into the following form:(23)Kij=∫01ϕi″″ξ+u2+P1−2μϕi″ξϕjξdξ.The comparison function of the clamped-clamped pipe is expressed as follows [35]:(24)ϕiξ=sinhβiξ−sinβiξ+λicoshβiξ−cosβiξ,where λi=sinhβil−sinβil/cosβil−coshβil and βil≈i+1/2π, when i ≥ 2. The values of βil are listed in Table 4. With different fluid velocities (u=0,u=1,andu=2), the deflections of the clamped-clamped pipe at different moments (τ = 10, τ = 20, and τ = 30) are shown in Figure 17.Table 4
The values ofβil for the clamped-clamped PCF.
β
1
l
β
2
l
β
3
l
β
4
l
4.7300
7.8532
10.9956
14.1372Figure 17
The deflections of the clamped-clamped pipe having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)The deflections of the clamped-clamped pipe are the smallest among the three types of pipes due to two reasons: one is the clamped-clamped pipe having the largest stiffness and the other is the existence of the axial tensile force on the cross section of the pipe. In Figures17(a) and 17(b), the largest amplitudes are different slightly because of counteracting influences of fluid and the axial tensile force. In Figures 17(b) and 17(c), the maximum deflection occurs at the left side of the pipe center.The first four orders of natural frequencies of the clamped-clamped PCF are listed in Table5. In resonance analysis, the resonance responses with fluid velocity u = 2 and supports movement frequencies ω = 20.79, 60.10, and 119.27 are considered. The detailed reflection curves are illustrated in Figure 18.Table 5
The first four dimensionless natural frequencies of the clamped-clamped PCF.
u
ω
ω
1
ω
2
ω
3
ω
4
0
22.23
61.47
120.69
199.63
1
21.87
61.13
120.33
199.27
2
20.79
60.10
119.27
198.20Figure 18
The resonance responses of the clamped-clamped pipe withu = 2 under different supports movement frequencies: (a) ω = 20.79; (b) ω = 60.10; (c) ω = 119.27.
(a)
(b)
(c)The strongest resonance appears atω = 119.27 and the weakest at ω = 60.10. The deflection curves in Figure 18(a) are almost symmetric, which are different from that of either the pinned-pinned pipe or the clamped-pinned pipe. In Figures 18(b) and 18(c), the maximum deflection appears at the right side of the pipe center.Figure19 gives the history curves of deflection of the clamped-clamped fluid-conveying pipe at different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) when having different fluid velocities (u = 0, u = 1, and u = 2).Figure 19
History curves of the deflection of the clamped-pinned pipe at different locations when having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)From Figure19, one can find that the flexural rigidity of the pipe is weakened by flowing fluid. The influence of fluid on the pipe is lower than those on the other two pipes.History curves of bending moment on the clamped-clamped fluid-conveying pipe are shown in Figure20.Figure 20
History curves of bending moment on the clamped-clamped pipe with different fluid velocities at different locations: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)When the pipe is with stationary fluid, i.e.,u = 0 (Figure 20(a)), the maximum bending moment appears at ξ = 1/2 and the minimum at ξ = 1/4. When the fluid velocity becomes u = 1 (Figure 20(b)), the bending moment at ξ = 1/4 tends to be identical to that at ξ = 1/8. When the fluid velocity reaches u = 2 (Figure 20(c)), the bending moment at ξ = 1/4 is higher than that at ξ = 1/8 and is almost the same as that at ξ = 1/2.


## 4.1. The Pinned-Pinned PCF

In beam theory, the comparison function for the pinned-pinned pipe, shown in Figure4, is ϕiξ=siniπξ. To show the effectiveness of the proposed method, we compare the results obtained by both the analytical mode superposition method for classical beam problems and the proposed method. For a solid beam, it can be considered as a PCF with zero velocity of fluid (u = 0) and without inner pressure (p¯=0). Comparison of the results will be carried out when u = 0 and p¯=0. It is easy to find out Equation (12) draws the state of a pipe with stationary supports and subjected to uniform load pξ,τ=γω2sinωτ. In this case, based on the mode superposition method and beam theory [32], the lateral deflection of the pipe can be expressed as follows:(17)ηξ,τ=∑iϕiξqiτ.Figure 4
Pinned-pinned PCF.Considering the uniform external loadpξ,τ=γω2sinωτ and initial condition, i.e., ηξ,0=η˙ξ,0=0, the generalized coordinate can be expressed analytically [32] as follows:(18)qiτ=1−−1i1ωi2γω2iπ∫0τsinωζsinωiτ−ζdζ=2γω2iπ1−−1iωi2−ω2sinωτ−ωωisinωiτ,where ωi is the ith natural frequency of the pinned-pinned pipe with stationary fluid and zero inner pressure.Now, the analytical solution for the pipe deflection can be obtained by substituting Equation (18) into Equation (17). We find that the results obtained by the mode superposition method are stable when i ≥ 3. Hence, we set i = 4 for the comparison. The results calculated by the two methods (Figure 5) match each other very well. Hence, the validation of the present method is verified.Figure 5
Comparison of deflection of the pinned-pinned pipe by two methods.We choose three different truncations of Equation (13), with N = 4, 6, and 8, in finding the responses of the pinned-pinned pipe. As demonstrated in Figure 6, the three results are almost the same and converge to the analytical solutions as long as N≥4.Figure 6
Comparison of deflections of the pinned-pinned pipe with differentN.The lateral deflections of the pinned-pinned pipe with three different fluid velocities (u = 0, u = 1, and u = 2) at three different moments (τ = 10, τ = 20, and τ = 30) are involved in discussion in this paper. The detailed results are illustrated in Figure 7.Figure 7
The deflections of the pinned-pinned pipe with different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)The three modes of deflection curves in Figure7 are different in each case. Meanwhile, the deflection curves in Figures 7(b) and 7(c) are asymmetric due to flow of fluid within the beam. It is easy to find out in Figure 7 that the deflections gradually increase with fluid velocity. This implies that the rigidity of the pipe becomes lower when the fluid speed is higher. The effect also appears on the natural frequencies of the pipe with different fluid velocities. In this paper, the first four orders of natural frequency are calculated through the wave approach which was introduced by Li et al. [33, 34]. The results are listed in Table 1.Table 1
The dimensionless natural frequencies of the pinned-pinned PCF.
u
ω
ω
1
ω
2
ω
3
ω
4
0
9.05
38.68
88.04
157.13
1
8.44
38.22
87.61
156.71
2
6.37
36.81
86.31
155.46In practice, when the external vibration frequency is in coincidence with the natural frequencies of structures, then structural resonance will happen, and large deformation will be generated. To evaluate such an effect on pipe structures, we investigate the deflections of PCF with supports moving in the same frequencies as its natural frequencies.For simplicity, the deflections of the pinned-pinned pipe with fluid velocityu = 1 are analyzed, and three cases with respect to the moving frequency of the supports, i.e., ω = 8.44, 38.22, and 87.61, are considered. Figure 8 shows the resonance responses in the three cases.Figure 8
The resonance responses of the pinned-pinned PCF withu = 1 under different supports moving frequencies: (a) ω = 8.44; (b) ω = 38.22; (c) ω = 87.61.
(a)
(b)
(c)Compared with the results in Figure7, the amplitudes of deflections in Figure 8 become much larger, especially when the supports moving frequency reaches the third natural frequency of PCF. At the times τ = 20 and τ = 30 in Figure 8(c), the largest amplitude of deflection appears near the right end of the pipe rather than at the center. The same phenomenon can be found in Figure 7(c).We also obtain the history curves of pipe deflection at different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) under different fluid velocities (u = 0, u = 1, and u = 2) (Figure 9).Figure 9
History curves of deflection at different locations on the pinned-pinned pipe having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)Figure9 suggests that the amplitude of deflection increases with the fluid velocity. The maximum deflection occurs at the center of the pipe. Deflection of a point becomes smaller when the point is closer to a pipe end. With the increase of fluid velocity, the response curves are gradually getting sparse. It implies that the response period becomes larger; i.e., the natural frequency is getting smaller.Bending moment is another important dynamic response for the beam. Once it is obtained, the normal stress on the pipe cross section will be determined. Furthermore, the fatigue analysis can be conducted. From Equation (13), the expression of dimensionless bending moment becomes(19)M¯=MLEI=∂2η∂ξ2.Figure10 illustrates the history curves of bending moment on the pinned-pinned pipe at three different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) with respect to three different fluid velocities (u = 0, u = 1, and u = 2).Figure 10
History curves of bending moment on the pinned-pinned pipe at different locations and having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)In Figure10(a), when u = 0, the bending moment at ξ = 1/2 is larger than those at ξ = 1/4 and ξ = 1/8, but the bending moments at ξ = 1/4 and ξ = 1/8 are different slightly. As the fluid velocity increases, the bending moments at the three locations increase, simultaneously. However, at ξ = 1/4 and ξ = 1/8, they tend to have a faster increase than that at location ξ = 1/2.


## 4.2. The Clamped-Pinned PCF

Another type of fluid-conveying pipe is the clamped-pinned pipe in Figure11. Its comparison function [35] becomes(20)ϕiξ=sinβiξ−sinhβiξ+αicoshβiξ−cosβiξ,where αi=sinβil−sinhβil/cosβil−coshβil and βil≈i+1/4π when i ≥ 1. The values of βil are listed in Table 2 [35]. At three different times (τ = 10, τ = 20, and τ = 30), the deflection responses of the pipe having three different fluid velocities (u = 0, u = 1, and u = 2) are illustrated in Figure 12.Figure 11
The model for a clamped-pinned PCF.Table 2
The values ofβil for clamped-pinned PCF.
β
1
l
β
2
l
β
3
l
β
4
l
3.9266
7.0686
10.2102
13.3518Figure 12
The deflection responses of the clamped-pinned pipe under different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)There is no symmetric mode in any case of Figure12 due to the asymmetric supports. Deflection of the pipe does not vary monotonously with the fluid velocity. At time τ = 10, with the increase of fluid velocity, the deflection first increases and then turns to decrease. However, at time τ = 20, the deflection always decreases. At time τ = 30, the deflection firstly decreases and then increases.The first four natural frequencies of the clamped-pinned PCF are calculated by the wave approach [33, 34] (Table 3). The resonance analysis with respect to fluid velocity u = 1 is also carried out when the supports moving frequencies ωi are equal to 14.36, 48.89, and 103.15, respectively.Table 3
The first four orders of dimensionless natural frequencies of the clamped-pinned PCF.
u
ω
ω
1
ω
2
ω
3
ω
4
0
14.82
49.29
103.54
177.55
1
14.36
48.89
103.15
177.16
2
12.94
47.69
101.98
176.00The deflection curves in Figures13(a) and 13(c) are much larger than that in Figure 13(b), and the largest deflection appears in Figure 13(c). In other words, the resonance is very strong once the supports moving frequency equals the first or the third natural frequency of the pipe. The largest deflection occurs at the right side of the pipe center in Figure 13(a), whilst at the left side of the pipe center in Figures 13(b) and 13(c).Figure 13
The resonance responses of the clamped-pinned pipe withu = 1 and different supports moving frequencies: (a) ω = 14.36; (b) ω = 48.89; (c) ω = 103.15.
(a)
(b)
(c)The history curves of the deflection of the clamped-pinned pipe conveying fluid at different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) and having different fluid velocities (u = 0, u = 1, and u = 2) are illustrated in Figure 14.Figure 14
History curves of the deflection of the clamped-pinned pipe at different locations and having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)It can also be observed that the deflections increase with the fluid velocity. Compared with the pinned-pinned pipe, however, the clamped-pined pipe has lower variation of deflection due to larger support stiffness at the clamped end of the clamped-pinned pipe, and the same fluid velocity produces smaller effect on the pipe deflection.The history curves of bending moment on the clamped-pinned pipe in different cases are given in Figure15. Figure 15 indicates that, in each case, the largest bending moment on the clamped-pinned pipe appears at ξ = 1/2, and this is the same as that of the pinned-pinned pipe. The bending moment at ξ = 1/8, however, is larger than that at ξ = 1/4, and this is different from that of the pinned-pinned pipe. And the magnitude of bending moment on the clamped-pinned pipe is much lower than that of the pinned-pinned pipe.Figure 15
History curves of bending moment response for the clamped-pinned pipe at different locations and under different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)


## 4.3. The Clamped-Clamped PCF

The dynamic model for the clamped-clamped pipe, shown in Figure16, is different from that for either the pinned-pinned pipe or the clamped-pinned pipe because the inextensible assumption is not suitable anymore. It is easy to find from Figure 16 that while vibration happens, the axial extension is inevitably permitted as the beam has transverse vibration. If so, an additional tensile force, which equals T=2μp¯Af, should be induced on the cross section of the pipe [1]. Hence, the dynamic equation of the clamped-clamped PCF needs to be modified as follows:(21)EI∂4w∂x4+mfV2+p¯A1−2μ∂2w∂x2+2mfV∂2w∂x∂t+mf+mp∂2w∂t2=mf+mpA0ω02sinω0t.Figure 16
Schematic of a clamped-clamped PCF.Correspondingly, the dimensionless form of Equation (21) can be written as follows:(22)∂4η∂ξ4+u2+P1−2μ∂2η∂ξ2+2β1/2u∂2η∂ξ∂τ+∂2η∂τ2=γω2sinωτ.In the discrete form of Equation (14), only the stiffness matrix needs to be changed into the following form:(23)Kij=∫01ϕi″″ξ+u2+P1−2μϕi″ξϕjξdξ.The comparison function of the clamped-clamped pipe is expressed as follows [35]:(24)ϕiξ=sinhβiξ−sinβiξ+λicoshβiξ−cosβiξ,where λi=sinhβil−sinβil/cosβil−coshβil and βil≈i+1/2π, when i ≥ 2. The values of βil are listed in Table 4. With different fluid velocities (u=0,u=1,andu=2), the deflections of the clamped-clamped pipe at different moments (τ = 10, τ = 20, and τ = 30) are shown in Figure 17.Table 4
The values ofβil for the clamped-clamped PCF.
β
1
l
β
2
l
β
3
l
β
4
l
4.7300
7.8532
10.9956
14.1372Figure 17
The deflections of the clamped-clamped pipe having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)The deflections of the clamped-clamped pipe are the smallest among the three types of pipes due to two reasons: one is the clamped-clamped pipe having the largest stiffness and the other is the existence of the axial tensile force on the cross section of the pipe. In Figures17(a) and 17(b), the largest amplitudes are different slightly because of counteracting influences of fluid and the axial tensile force. In Figures 17(b) and 17(c), the maximum deflection occurs at the left side of the pipe center.The first four orders of natural frequencies of the clamped-clamped PCF are listed in Table5. In resonance analysis, the resonance responses with fluid velocity u = 2 and supports movement frequencies ω = 20.79, 60.10, and 119.27 are considered. The detailed reflection curves are illustrated in Figure 18.Table 5
The first four dimensionless natural frequencies of the clamped-clamped PCF.
u
ω
ω
1
ω
2
ω
3
ω
4
0
22.23
61.47
120.69
199.63
1
21.87
61.13
120.33
199.27
2
20.79
60.10
119.27
198.20Figure 18
The resonance responses of the clamped-clamped pipe withu = 2 under different supports movement frequencies: (a) ω = 20.79; (b) ω = 60.10; (c) ω = 119.27.
(a)
(b)
(c)The strongest resonance appears atω = 119.27 and the weakest at ω = 60.10. The deflection curves in Figure 18(a) are almost symmetric, which are different from that of either the pinned-pinned pipe or the clamped-pinned pipe. In Figures 18(b) and 18(c), the maximum deflection appears at the right side of the pipe center.Figure19 gives the history curves of deflection of the clamped-clamped fluid-conveying pipe at different locations (ξ = 1/2, ξ = 1/4, and ξ = 1/8) when having different fluid velocities (u = 0, u = 1, and u = 2).Figure 19
History curves of the deflection of the clamped-pinned pipe at different locations when having different fluid velocities: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)From Figure19, one can find that the flexural rigidity of the pipe is weakened by flowing fluid. The influence of fluid on the pipe is lower than those on the other two pipes.History curves of bending moment on the clamped-clamped fluid-conveying pipe are shown in Figure20.Figure 20
History curves of bending moment on the clamped-clamped pipe with different fluid velocities at different locations: (a)u = 0; (b) u = 1; (c) u = 2.
(a)
(b)
(c)When the pipe is with stationary fluid, i.e.,u = 0 (Figure 20(a)), the maximum bending moment appears at ξ = 1/2 and the minimum at ξ = 1/4. When the fluid velocity becomes u = 1 (Figure 20(b)), the bending moment at ξ = 1/4 tends to be identical to that at ξ = 1/8. When the fluid velocity reaches u = 2 (Figure 20(c)), the bending moment at ξ = 1/4 is higher than that at ξ = 1/8 and is almost the same as that at ξ = 1/2.


## 5. Concluding Remarks

The dynamic responses of three types of fluid-conveying pipes with lateral moving supports are investigated in this paper. The dynamic equations are deduced after obtaining the acceleration of fluid. The Galerkin method and fourth-order Runge–Kutta method are adopted in the numerical analysis. The responses with respect to deflection and bending moment are obtained considering different fluid velocities. Some conclusions are drawn according to the numerical analysis:(1)
Both the deflection and bending moment responses increase with the fluid velocity.(2)
The largest bending moment on the pinned-pinned pipe appears atξ = 1/2. The difference between the bending moments at ξ = 1/4 and ξ = 1/8 is slight, which is independent of the fluid velocity.(3)
The magnitudes of bending moment on the clamped-pinned and clamped-clamped pipes satisfyM¯ξ=1/2>M¯ξ=1/8>M¯ξ=1/4 when the fluid is in the stationary state. However, the bending moment at ξ=1/4 increases sharply with the fluid velocity and tends to be identical to that at ξ=1/8.(4)
Strong resonance will happen when the supports movement frequency equals the odd-ordered natural frequencies of the PCF.



---

*Source: 3295787-2018-11-04.xml*