Acessibilidade / Reportar erro

Effect of Open Crack on Vibration Behavior of a Fluid-Conveying Pipe Embedded in a Visco-Elastic Medium

Abstract

In this paper vibration behavior of a fluid-conveying cracked pipe surrounded by a visco-elastic medium has been considered. During this work, the effect of an open crack parameters and flow velocity profile shape inside the pipe on natural frequency and critical flow velocity of the system has been analytically investigated. An explicit function for the local flexibility of the cracked pipe has been offered using principle of the fracture mechanics. Comparison between the results of the present study and the experimental data reported in the literature reveals success and high accuracy of the implemented method. It is demonstrated that the existence of the crack in the pipe, decreases the natural frequency and the critical flow velocity so that the system instability onsets at a lower flow velocity in comparison with the intact pipe. Results indicate that the flow velocity profile shape inside the pipe caused by the viscosity of real fluids, significantly affects the critical flow velocity of both intact and fluid-conveying cracked pipe. For instance, as the flow-profile-modification factor decreases from 1.33 to 1.015, the dimensionless critical flow velocity of intact clamped-clamped pipe increases from 5.45 to 6.24.

Keywords:
Cracked Fluid-conveying Pipe; Vibration Analysis; Natural Frequency; Velocity Profile; Critical Flow Velocity

NOMENCLATURE

1 INTRODUCTION

Dynamic behavior of the fluid-conveying pipes has significant importance due to its extensive application in different industries such as petrol and gas transportation systems (Païdoussis 1998Païdoussis, M. P. (1998). Fluid-Structure Interactions. Fluid-Structure Interactions, Academic Press., Ibrahim 2010Ibrahim, R. A. (2010). "Overview of Mechanics of Pipes Conveying Fluids-Part I:Fundamental Studies" Journal of Pressure Vessel Technology 132: 034001-034031., Ghaitani and Esmaeili 2013Ghaitani, M. M. and H. A. Esmaeili (2013). "Vibration response of oil pipes embedded in viscous fluid using DQM".2nd International Conference on Emerging Trends in Engineering & Technology. College of Engineering, Teerthanker Mahaveer University.). These pipes often lie on a foundation and would be cracked as a result of various internal and external loads while conveying fluid. So, analyzing the dynamic behavior of the cracked pipes would be practically important.

Inspired by the work of Tada et al. (1973)Tada, H., P. C. Paris and G. R. Irwin (1973). Stress Analysis of Cracks Handbook. New York, American Society of Mechanical Engineers., many researchers have performed the vibration analysis of the cracked structures such as thin walled tubes during the last five decades. Recently, some reviews on vibration of cracked structures were reported by Dimarogonas (1996)Dimarogonas, A. D. (1996). "Vibration of cracked structures: a state of the art review." Engineering fracture mechanics 55(5): 831-857. and Kumar et al. (2010)Kumar, C. and V. Rastogi (2010). "A brief review on dynamics of a cracked rotor" International Journal of Rotating Machinery 2009. . In some researches the crack has been simulated by a local reduced section or stiffness of the structure around the crack (Sinha, Friswell et al. 2002Sinha, J. K., M. I. Friswell and S. Edwards (2002). "Simplified Models for the Location of Cracks in Beam Structures Using Measured Vibration Data." Journal of Sound and Vibration 251(1): 13-38.). The local flexibility of the structure is reduced due to the crack and this concept has been used to model the vibration behavior of the cracked pipe in some studies (Xie 1998Xie, Y. (1998). "A theory on cracked pipe." International journal of pressure vessels and piping 75(12): 865-869., Mahjoob and Shahsavari 2007Mahjoob, M. J. and a. A. Shahsavari (2007). A Vibration-Based Damage Detection Method for Pipes Conveying Fluid. 48th AIAA/ ASME/ ASCE/ AHS/ ASC Structures, Structural Dynamics, and Materials Conference. Honolulu, Hawaii., Hu et al. 2011Hu, J.-S., W.-Y. Sun and J. Zhou (2011). "Vibration analysis and crack identification of a cantilever pipe with a circumferential part-through crack." Journal of Vibration and Shock 4: 008., Yumin et al. 2011Yumin, H., Y. Junjie and H. Zhengjia (2011). "Research on Quantitative Identification of Pipe Crack Based on Vibration." Measuring Technology and Mechatronics Automation (ICMTMA), 2011 Third International Conference on, IEEE., liu et al. 2012Liu, D., T. Zhu, J. Zhou And X.-G. Zhao (2012). "Study on Vibration Characteristics of a Pipe with a Crack" China Offshore Platform 2: 007., Bai et al. 2013Bai, Q., X. Shang and L. Yin (2013). "Free Vibration Analysis of Elastic Pipe with Crack Defects." Journal of Applied Sciences 13: 5440-5445.). Literature review indicates that available studies about the vibration behavior of the cracked pipes can be divided into three groups;

  • Having no fluid inside,

  • Filled with fluid without flow (zero flow velocity)

  • Having fluid flow inside (fluid-conveying pipe).

There are fairly large amount of researches dedicated to the first two groups; Liu et al. (2003)Liu, D., H. Gurgenci and M. Veidt (2003). "Crack detection in hollow section structures through coupled response measurements." Journal of sound and vibration 261(1): 17-29. examined the stability of thin-walled pipe due to circumferential cracks analytically and experimentally based on the coupled-response measurements. Ye et al. (2010)Ye, J., Y. He, X. Chen, Z. Zhai, Y. Wang and Z. He (2010). "Pipe crack identification based on finite element method of second generation wavelets." Mechanical Systems and Signal Processing 24(2): 379-393. investigated the vibration behavior of the cracked pipe in the absence of the fluid flow using the finite element method. Their results show that the presence of the crack in the pipe reduces the natural frequency, and its effect on the upper natural frequencies is more than the lower ones. Murigendrappa et al. (2004aMurigendrappa, S. M., S. K. Maiti and H. R. Srirangarajan (2004). "Experimental and theoretical study on crack detection in pipes filled with fluid." Journal of Sound and Vibration 270(4-5): 1013-1032.; 2004bMurigendrappa, S. M., S. K. Maiti and H. R. Srirangarajan (2004). "Frequency-based experimental and theoretical identification of multiple cracks in straight pipes filled with fluid." NDT & E International 37(6): 431-438.) studied the vibration behavior of the cracked pipe theoretically and experimentally. In their works, the crack has been considered as a torsional spring and the transfer matrix method with Rayleigh's quotient has been used for vibration analysis. ) Naniwadekar et al. (2008)Naniwadekar, M. R., S. S. Naik and S. K. Maiti (2008). "On prediction of crack in different orientations in pipe using frequency based approach." Mechanical Systems and Signal Processing 22(3): 693-708. developed an experimental method for crack detection based on changing the natural frequency and presented a numerical method to analyze the vibration behavior of a straight horizontal steel hollow pipe filled with an incompressible fluid (without flow).

In contrast to the first two groups, there are fairly limited amount of investigation about the vibration behavior of the cracked fluid-conveying pipes (third group). Yoon and Soo are two researchers that have contributed in most papers published about vibration behavior of the fluid-conveying cracked pipes. They investigated the effects of the open crack parameters and the moving mass on the dynamic behavior of simply supported pipe conveying fluid using Timoshenko beam model. Numerical methods based on the transfer matrix have been preferred in their works (Son, Lee et al.Son, I.-S., S.-P. Lee, J.-K. Lee and D.-S. Bae "Dynamic Behavior of Forced Vibration of Elastically Restrained Pipe Conveying Fluid with Crack and Concentrated Mass." , Yoon and Son 2004Yoon, H.-I. and I.-S. Son (2004). "Dynamic behavior of cracked pipe conveying fluid with moving mass based on Timoshenko beam theory." KSME international journal 18(12): 2216-2224., Yoon, Son et al. 2004Yoon, H. I., I. S. Son and J. T. Jin (2004). "A study on dynamic behavior of simply supported fluid flow pipe with crack and moving mass.", Yoon and Son 2005Yoon, H.-I. and I.-S. Son (2005). "Influence of tip mass on dynamic behavior of cracked cantilever pipe conveying fluid with moving mass." Journal of mechanical science and technology 19(9): 1731-1741., Yoon and Son 2006Yoon, H.-I. and I.-S. Son (2006). "Dynamic behavior of cracked simply supported pipe conveying fluid with moving mass." Journal of Sound and Vibration 292(3-5): 941-953., Son, Cho et al. 2007Son, I., J. Cho and H. Yoon (2007). "Effects of Attached Mass on Stability of Pipe Conveying Fluid with Crack." Transactions of KSNVE 17(10): 1002-1009., Yoon and Son 2007Yoon, H.-I. and I.-S. Son (2007). "Dynamic response of rotating flexible cantilever pipe conveying fluid with tip mass." International journal of mechanical sciences 49(7): 878-887., Yoon, Son et al. 2007Yoon, H., I. Son and S. Ahn (2007). "Free Vibration Analysis of Euler-Bernoulli Beam with Double Cracks." Journal of Mechanical Science and Technology 21: 476-485., Son and Yoon 2008Son, I. and H. Yoon (2008). "Dynamic Stability of Elastically Restrained Cantilever Pipe Conveying Fluid with Crack." Transactions of KSNVE 18(2): 177-184., Son, Yoon et al. 2010Son, I. S., H. I. Yoon, S. P. Lee and D. J. Kim (2010). "Effects of Tip Mass on Stability of Rotating Cantilever Pipe Conveying Fluid with Crack." International Journal of Modern Physics B 24(15n16): 2609-2614.). So, available research in literature did not consider the effects of the crack parameters on the critical flow velocity which is our main objective in this paper.

In present paper the vibration behavior of the cracked fluid-conveying pipe is examined via Euler-Bernoulli beam model focusing on natural frequency and critical flow velocity. Also, in order to have a more realistic model, the surrounding elastic medium is modeled via a visco-elastic foundation. Our study is carried out through an analytical method. The crack is considered using a massless torsional spring. An explicit function for obtaining the local flexibility of the cracked pipe is presented using the theory of fracture mechanics. Also, we demonstrate that how the fluid flow velocity profile would influence the critical flow velocity of the system.

2 LOCAL FLEXIBILITY OF THE CRACKED PIPES

Dimarogonas and Papadopoulos (1983)Dimarogonas, A. S. and C. A. Papadopoulos (1983). "Vibration of cracked shafts in bending." Journal of Sound and Vibration 91(4): 583-593. proposed the massless torsional spring model of the cracked beam and shaft, and derived the additional flexibility due to the crack using the principle of fracture mechanics. Following this procedure, a small vertical sectional strip is considered in the cracked section of the pipe. The strain energy of the strip is obtained using the theory of fracture mechanics. By integrating of the strain energy over the cracked section, for the local flexibility, an explicit function is obtained as a function of the crack depth.

The additional strain energy due to the crack can be considered in the form of a flexibility coefficient expressed in terms of the stress intensity factor, which can be derived by Castigliano's theorem in the linear elastic range. Therefore, the local flexibility in the presence of the crack is defined by (Zou et al. 2003Zou, J., J. Chen, J. C. Niu and Z. M. Geng (2003). "Discussion on the local flexibility due to the crack in a cracked rotor system." Journal of Sound and Vibration 262(2): 365-369.):

Where E is Young module, v is Poisson's ratio, a and b are the crack dimensions as shown in Figure 1. KI is the stress intensity factor corresponding to the first fracture mode due to the bending moment M.

Figure 1:
Cross section of the cracked pipe.

It is seen from Eq. (1) that the local flexibility depends on the stress intensity factor. In order to calculate the stress intensity factor for the cracked pipes, consider a small vertical rectangular strip, having a small cross-sectional width dx and height hx with an arbitrary offset distance x, as shown in Figure 1. In Ref. (Tada, Paris et al. 2000) the stress intensity factor for a rectangular section is presented as:

where, a is the crack depth, i and h are the second moment of area and height of the rectangular cross-section, respectively and F(a/h) is a mathematical function of crack dimension. Employing Eq. (2) for the rectangular strip shown in Figure 1 leads:

in which, ξ is local coordinate which varies from ξ = 0 to ξ = ax swiping the crack depth in the given strip, and F/hx) is defined as (Tada et al. 1973Tada, H., P. C. Paris and G. R. Irwin (1973). Stress Analysis of Cracks Handbook. New York, American Society of Mechanical Engineers.):

Substituting Eq. (3) into Eq. (1) leads to:

The equivalent rotational spring stiffness, Kt, is defined as (Dimarogonas and Papadopoulos, 1983Dimarogonas, A. S. and C. A. Papadopoulos (1983). "Vibration of cracked shafts in bending." Journal of Sound and Vibration 91(4): 583-593.):

Taking into account the coordinate system illustrated in Figure 1, hx and b calculated as:

and the upper boundary of the crack, ax , can be expressed as:

By defining the dimensionless variables; = x/Ro, = hx/Ro, = b/Ro, = ax/ Ro and αc= a/ Ro, the dimensionless form of Eqs. (7)-(9) becomes:

Substitution of Eqs. (10)-(12) into Eq. (5) leads:

Where = ξ/Ro. For the flexibility coefficient, the dimensionless form introduced as:

it can be derived from Eq. (13):

Eq. (15) represents an integral form for the local flexibility of the cracked pipes. This integral expression is only a function of the dimensionless crack depth, αc= a/ Ro. By numerical integration of Eq. (15) and fitting an appropriate curve to the numerical values of and then using the least square method, an explicit expression for the dimensionless flexibility coefficient of the cracked pipe is obtained as:

It should be noted that according to the main assumption in the fracture mechanics, Eq. (16) is satisfied up to middle of the pipe wall thickness i.e.a = (Ro - Ri)/2.

This explicit function may be applied for any pipe with specified mechanical properties and geometrical dimensions. Considering Eqs. (16) and (14), the equivalent stiffness of the cracked pipe, Kt is calculated straightforwardly from Eq. (6). In order to verify the accuracy of this model, the predicted values of Kt and the experimental results reported by Murigendrappa et al. (2004)Murigendrappa, S. M., S. K. Maiti and H. R. Srirangarajan (2004). "Frequency-based experimental and theoretical identification of multiple cracks in straight pipes filled with fluid." NDT & E International 37(6): 431-438. for a cracked pipe with a parameters Ro = 16.5mm, Ri = 10mm, E = 65GPa and v = 0.33 have been compared in Table 1. The comparison shows that there is a good agreement between the equivalent stiffness obtained from the proposed method and the experimental data. So, the maximum difference between the predicted results and the experimental data of Murigendrappa et al. (2004)Murigendrappa, S. M., S. K. Maiti and H. R. Srirangarajan (2004). "Frequency-based experimental and theoretical identification of multiple cracks in straight pipes filled with fluid." NDT & E International 37(6): 431-438. is less than 5%.

Table 1:
The local stiffness due to crack, Comparison between the calculated results and the experimental data of Murigendrappa et al. (2004)Murigendrappa, S. M., S. K. Maiti and H. R. Srirangarajan (2004). "Frequency-based experimental and theoretical identification of multiple cracks in straight pipes filled with fluid." NDT & E International 37(6): 431-438..

The effect of the crack depth on the dimensionless local flexibility coefficient is shown in Figure 2. The figure indicates that an increase in the crack depth enhances the flexibility coefficient of the cracked pipe because of the reduction in equivalent stiffness (Eq. 6).

Figure 2:
Effect of the crack depth on the dimensionless local flexibility coefficient.

3 MATHEMATICAL MODEL

Figures 3a and 3b illustrate schematically a fluid-conveying cracked pipe resting on a visco-elastic foundation and its corresponding mathematical model, respectively. Using the Euler-Bernoulli theory, the modified governing equation of motion for the uniform fluid-conveying pipe which is subjected to an external force, Fext, can be written as (Guoa et al. 2010Guoa,C.Q, C.H. Zhanga and M.P. Paidoussis (2010). "Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profile." Journal of Fluids and Structures 26: 793-803.):

where mp and mf are the mass per unit length of pipe and fluid, respectively. t stands for time, X is axial coordinate, W(X ,t) and EI are the transverse deflection and the flexural stiffness of the pipe, respectively. α denotes a dimensionless parameter named as the flow-profile-modification factor dependent on fluid velocity profile inside pipe. Considering the power-law profile for the time-average flow velocity profile (Streeter et al. 1998Streeter, V. L., E. B. Wylie and K. W. Bedford (1998). Fluid mechanics, WCB, McGraw-Hill.) as:

α can be determined by (Guoa et al. 2010Guoa,C.Q, C.H. Zhanga and M.P. Paidoussis (2010). "Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profile." Journal of Fluids and Structures 26: 793-803.):

Figure 3:
(a) Schematic view of the cracked clamped-clamped pipe, (b) Mathematical model.

Where n is dependent on the Reynolds number. The well-known one-seven power-law velocity profile introduced by Prandtl can be achieved for n=7 that is valid for the smooth-wall circular pipes in the range 3×103≤Re≤3×105 (Re denotes the Reynolds number). Also, n=4 and n=5 provide the velocity profile of a rough-wall circular pipe based on the experiments of Nikuradse. For n=7, α =1.020; and for n=4 and n=5, α≈1.037-1.055. We will discuss more about the impact of α on the vibration behavior of the intact and cracked fluid-conveying pipes in subsection 6.2.

In the present work, the external force related to the visco-elastic foundation modeled using two-parameter Kelvin-Voigt model. Based on this model Fext can be written as:

where, Km and C stand for the Winkler stiffness constant and the damping coefficient of the visco-elastic foundation, respectively. Ks is an additional parameter defining the foundation, usually termed as the shear constant of the foundation (Kargarnovin et al. 2005Kargarnovin, M. H., D. Younesian, D. J. Thompson and C. J. C. Jones (2005). "Response of beams on nonlinear viscoelastic foundations to harmonic moving loads." Computers & Structures 83(23-24): 1865-1877., Nguyen and Duhamel 2008Nguyen, V.-H. and D. Duhamel (2008). "Finite element procedures for nonlinear structures in moving coordinates.Part II: Infinite beam under moving harmonic loads." Computers & Structures 86(21-22): 2056-2063., Ansari et al. 2011Ansari, M., E. Esmailzadeh and D. Younesian (2011). "Frequency analysis of finite beams on nonlinear Kelvin-Voight foundation under moving loads." Journal of Sound and Vibration 330(7): 1455-1471., Basu and Kameswara 2013Basu, D. and N. S. V. Kameswara Rao (2013). "Analytical solutions for Euler-Bernoulli beam on visco-elastic foundation subjected to moving load." International Journal for Numerical and Analytical Methods in Geomechanics 37(8): 945-960.). By substituting of Eq. (18) into Eq. (17), the governing equation of the fluid-conveying pipe resting on a visco-elastic foundation is obtained as:

It is assumed that the solution of Eq. (19) is:

where (X) and ω indicate the amplitude and natural frequency, respectively. The dimensionless variables are defined as:

Substituting Eq. (20) into Eq. (19) and using the above dimensionless variables, Eq. (19) can be rewritten as follows:

The solution of the Eq. (22) is:

where A and λ are constants. Substituting this expression in Eq. (22) leads to the following characteristic equation:

The Eq. (24) is a quartic polynomial equation, which has generally four complex roots (consider a typical real root as a complex root with zero imaginary part). The explicit algebraic formulas for a general solutions of Eq. (24) are expressed as (Persidis 2007Persidis, S. (2007). Mathematical Handbook. ESPI, ATHENS, ESPI PUBLISHING.):

where the parameters p1, p2 and p3 which depend on the dimensionless natural frequency β, have been defined in appendix. The general solution of Eq. (22) is

In which Ai, i = 1,2,3,4 are unknown constants and will be obtained by applying the boundary conditions.

4. FREQUENCY EQUATION

Consider the corresponding mathematical model of the cracked pipe conveying fluid resting on a visco-elastic medium as shown in Figure 3b. For analyzing the free vibration of the cracked pipe, the entire pipe is divided from the cracked section into the two pipe segments. The crack is represented by a torsional mass-less spring that its stiffness was calculated previously in section 2. The two pipe segments can be treated separately. The equations of motion for two mentioned intact segments are obtained from Eq. (26) as follows:

where wL(ζ) and wR(ζ) are the vibration mode of each segment in the left and the right side of the crack, respectively. These vibration modes have eight unknown constants Ai, i = 1,2,...,8, which can be obtained by applying both boundary conditions and compatibility conditions at the crack location. The compatibility of displacement, moment and shear force of both segments at the dimensionless crack location (i.e. ζ = ζc), are respectively given by:

and the discontinuity in the slope of the pipe at the cracked location may be implemented as

where is the dimensionless parameter indicating the effect of crack depth. General boundary conditions can be applied for the present model. In this paper, clamped-clamped supported pipe is considered. Thus the four boundary conditions may be written as:

Substituting Eqs. (27) and (28) into the boundary conditions Eq. (31) and into the compatibility conditions Eqs. (29) and (30), gives a set of eight homogeneous linear algebraic equations for the eight unknown constants Ai, i = 1,2,...,8 as follows:

where

The system of equations (32) has the nontrivial solution if and only if the determinant of the coefficients becomes zero. So the characteristic equation of the system will be:

Now, by solving the Eq. (34), the dimensionless natural frequency, β, can be obtained.

5 VALIDATION

To the best knowledge of the authors, there is no analytical or experimental result in the literature for the vibration analysis of a problem that is completely the same of one stated here. So in order to validate the aforementioned theoretical model and analytical procedure, we use available data for a similar problem. For instance, comparison made between the present work results and the experimental data of Mahjoob and Shahsavari (2007)Mahjoob, M. J. and a. A. Shahsavari (2007). A Vibration-Based Damage Detection Method for Pipes Conveying Fluid. 48th AIAA/ ASME/ ASCE/ AHS/ ASC Structures, Structural Dynamics, and Materials Conference. Honolulu, Hawaii. for the case of a clamped-clamped cracked pipe conveying fluid without foundation (i.e. km = ks = c = 0). The implemented mechanical and geometrical properties required to calculate the results are listed in Table 2. The maximum difference between the experimental data of Mahjoob and Shahsavari (2007)Mahjoob, M. J. and a. A. Shahsavari (2007). A Vibration-Based Damage Detection Method for Pipes Conveying Fluid. 48th AIAA/ ASME/ ASCE/ AHS/ ASC Structures, Structural Dynamics, and Materials Conference. Honolulu, Hawaii. and the present model results are 0.29%, 0.47% and 0.45% for the first, second and third natural frequencies, respectively (Table 3). Therefore, this comparison indicates an excellent agreement between present study and the experimental data.

Table 2:
Mechanical properties and geometrical dimensions of the investigated pipe, fluid and foundation.

Table 3:
Natural frequencies of the clamped-clamped fluid-conveying cracked pipe for ζc = 0.4. Comparison between results of present analytical method and the experimental data of Mahjoob and Shahsavari (2007)Mahjoob, M. J. and a. A. Shahsavari (2007). A Vibration-Based Damage Detection Method for Pipes Conveying Fluid. 48th AIAA/ ASME/ ASCE/ AHS/ ASC Structures, Structural Dynamics, and Materials Conference. Honolulu, Hawaii. for various flow velocity.

6 RESULTS AND DISCUSSION

Having insured the accuracy of the present model, it is used to examine the vibration behavior of the cracked pipe conveying fluid resting on a visco-elastic foundation. The elastic, shear and damping constants used for the foundation as a Kelvin-Voigt model are respectively Km=3.08 MPa, Ks=12550 kN and C=3110 Pa.s as given in (Gerolymos and Gazetas 2006Gerolymos, N. and G. Gazetas (2006). "Winkler model for lateral response of rigid caisson foundations in linear soil." Soil Dynamics and Earthquake Engineering 26(5): 347-361., Younesian et al. 2012Younesian, D., S. R. Marjani and E. Esmailzadeh (2012). "Nonlinear vibration analysis of harmonically excited cracked beams on viscoelastic foundations." Nonlinear Dynamics 71(1-2): 109-120.). In all numerical calculations, the flow-profile-modification factor, α, has been taken 1.33 except where it has been clearly specified.

6.1 Effect of the Crack Parameters

Figure 4, illustrates first natural frequency of the cracked pipe conveying fluid as a function of the dimensionless crack location, for three different values of αc =0.025, αc =0.25 and αc =0.4, where the foundation parameters were ks = km = 0.01 and c = 0. It is found that the crack parameters affect the vibration behavior of the cracked pipe. It is obvious that the frequencies of a system generally decrease due to the crack. Results show that when the crack location approaches to the fixed ends of the pipe, the effect of the crack on the frequency reduction is increased. This tendency is intensified with an increase in the crack depth. Also, as seen in Figure 4, for the crack located at ζc = 0.225 and ζc = 0.775, there is no reduction in the first natural frequency. Indeed, the mentioned locations are the inflection points for the first vibration mode of a clamped-clamped pipe. In these locations, the second derivation of the deflection function is equal to zero (d2W (X)/dX2 = 0) causing the bending moment to be completely vanished. Furthermore, since the major factor that decreases the natural frequencies caused by the crack is the bending moment, thus for the first vibration mode, the crack located at ζc = 0.225 and ζc = 0.775 has no effect on the first natural frequency. The same behavior is observed for the alteration of the second natural frequency shown in Figure 5.

Figure 4:
First dimensionless natural frequency of the fluid-conveying cracked pipe as a function of the dimensionless crack location.

Figure 5:
Second dimensionless natural frequency of the fluid-conveying cracked pipe as a function of the dimensionless crack location.

Figure 6 presents the natural frequency of the cracked pipe for the first mode as a function of the dimensionless crack depth for different values of ζc. When the crack depth,ac, increases, the local flexibility of the pipe reduces according to Eq. (7) causing reduction in the natural frequency.

Figure 6:
First natural frequency of the fluid-conveying cracked pipe.

6.2 Effect of the Flow Velocity

As seen from frequency equation (Eq. 34), the natural frequency of the cracked pipe conveying fluid in visco-elastic medium depends on the fluid flow velocity. Figure 7 shows the natural frequency of the intact pipe versus the dimensionless flow velocity for the first and second vibration modes. For the clamped-clamped intact pipe, the dimensionless critical flow velocity is = 2π (Doare and Langer 2002Doare, O. and E. D. Langer (2002). "Local and global instability of fluid-conveying pipes on elastic foundation." Journal of Fluids and Structures 16(1): 1-14.). As it is obvious from the Figure 7, when the flow velocity increases, the pipe natural frequency decreases smoothly until it completely vanishes for u=2π. This equals to the lowest critical flow velocity that causes instability in the system. It is worth to note that u=2π corresponding to the mean flow velocity of U=498.5 m/s. This high velocity for a typical compressible fluid corresponds to Mach numbers greater than unit (Ma>1). Since for Ma>0.3, the compressibility should be taken into account, the present results for the critical flow velocity are not valid for the compressible fluids.

Figure 7:
Variation of the natural frequency of the clamped-clamped fluid-conveying cracked pipe against the dimensionless flow velocity.

Experiments of Nikuradse revealed that fluid velocity profile in the circular pipes changes based on the quality of wall roughness. In order to take into account the effect of the mentioned change in the velocity profile, on vibration behavior of the fluid-conveying pipes, Guoa et al. (2010)Guoa,C.Q, C.H. Zhanga and M.P. Paidoussis (2010). "Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profile." Journal of Fluids and Structures 26: 793-803. derived a modified equation of motion for the pipe. The expression, a, appears in the centrifugal term of the modified equation (Eq. 17). As depicted in Figure 8, the flow-profile-modification factor would influence the critical flow velocity of both intact and cracked fluid-conveying pipe. As Guoa et al. (2010)Guoa,C.Q, C.H. Zhanga and M.P. Paidoussis (2010). "Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profile." Journal of Fluids and Structures 26: 793-803. indicated, when a decreases from 1.33 to 1.015, the velocity profile in the pipe cross section tends to a typical turbulent flow profile. This trendency corresponds to increase in critical flow velocity according to Figure 8. Thus it can be deduced that the instability in a typical fluid-conveying pipe would be delayed in turbulence regime in comparison with the laminar one. For an instance, the dimensionless critical flow velocity of intact pipe increases from 5.45 to 6.24 as the flow-profile-modification factor decreases from 1.33 to 1.015.

Figure 8:
Flow-profile-modification factor influences the impact of crack depth, φ, on the critical flow velocity, ucr for ζc = 0.01 (left vertical axis). Also, it affects the critical flow velocity of an intact clamped-clamped pipe (right vertical axis).

In the case of a cracked fluid-conveying pipe, the crack depth influences significantly the impact of flow-profile-modification factor on critical flow velocity. An interesting behavior observed in Figure 8 is about the effect of crack and its depth on the system instability. As it is seen, existence of a crack on the pipe causes the instability in a lower fluid flow velocity and this effect is intensified proportionally with crack depth for a given crack location (ζc = 0.001 in Figure 8). This behavior shows that the crack and the instability would practically reinforce each other in a pipe system causing its destruction.

6.3 Effect of the Visco-Elastic Foundation

Figure 9 shows the effect of the foundation parameters on the critical velocity as a function of the dimensionless damping factor, c, for the different foundation parameters. The comparison is made for the crack parameters αc = 0.2 and ζc = 0.1. As seen in Figure 9, divergence instability for the first mode occurs when ucr≃1.01 but by increasing the dimensionless damping parameter of the visco-elastic foundation, c and/or decreasing the dimensionless stiffness parameter, ks, the instability of the pipe takes place at a lower flow velocity.

Figure 9:
Variation of the dimensionless critical flow velocity of the clamped-clamped cracked pipe against the dimensionless damping factor for different values of the dimensionless stiffness, ks.

Figure 10 shows the first three dimensionless frequencies of the cracked pipe as a function of the dimensionless parameters km where c = 0.01, u = 0. The results indicate that any augmentation in the foundation stiffness constant is followed by increasing the natural frequencies. This means that increasing the elastic constant makes the pipe stiffer. Also, effect of the foundation stiffness is more pronounced in the lower natural frequencies than the upper ones.

Figure 10:
Effect of the foundation stiffness, km, on the first three natural frequencies of the fluid-conveying cracked pipe.

7 CONCLUSION

In the current study, the vibration behavior of a cracked pipe conveying fluid resting on a visco-elastic foundation was investigated analytically. A new approach for calculating the local flexibility caused by a crack was developed. The mentioned approach is applicable for various pipes with different mechanical properties. Comparison between the analytical results and available experimental data in literature shows good agreement for a wide range of the crack parameters.

The results indicate that increasing the crack depth improves the flexibility, and therefore the local stiffness reduces at the crack location. This leads to reduction in the both natural frequency and critical flow velocity of the cracked pipe system. This feature can be utilized for crack detection based on the vibration analysis. Also, when the crack location approaches to the fixed ends of the pipe, reduction in the natural frequency becomes more considerable. Furthermore, the results show that a decrease in the natural frequency reduces when the crack gets closer to the inflection points until it completely vanishes in the inflection points.

The fluid flow inside the pipe and its profile shape influences the natural frequency of the cracked pipe significantly. Moreover, the crack decreases the critical flow velocity so that the system instability onsets at a lower flow velocity.

Acknowledgement

The authors would like to appreciate Ahar Branch, Islamic Azad University for the financial support of this research, which is based on a research project contract.

References

  • Ansari, M., E. Esmailzadeh and D. Younesian (2011). "Frequency analysis of finite beams on nonlinear Kelvin-Voight foundation under moving loads." Journal of Sound and Vibration 330(7): 1455-1471.
  • Bai, Q., X. Shang and L. Yin (2013). "Free Vibration Analysis of Elastic Pipe with Crack Defects." Journal of Applied Sciences 13: 5440-5445.
  • Basu, D. and N. S. V. Kameswara Rao (2013). "Analytical solutions for Euler-Bernoulli beam on visco-elastic foundation subjected to moving load." International Journal for Numerical and Analytical Methods in Geomechanics 37(8): 945-960.
  • Guoa,C.Q, C.H. Zhanga and M.P. Paidoussis (2010). "Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profile." Journal of Fluids and Structures 26: 793-803.
  • Dimarogonas, A. D. (1996). "Vibration of cracked structures: a state of the art review." Engineering fracture mechanics 55(5): 831-857.
  • Dimarogonas, A. S. and C. A. Papadopoulos (1983). "Vibration of cracked shafts in bending." Journal of Sound and Vibration 91(4): 583-593.
  • Doare, O. and E. D. Langer (2002). "Local and global instability of fluid-conveying pipes on elastic foundation." Journal of Fluids and Structures 16(1): 1-14.
  • Ettefagh,Mir.M., Hossein Akbari, Keivan Asadi and Farshid Abbasi (2015). "New structural damage-identification method using modal updating and model reduction." Proc IMechE Part C: J Mechanical Engineering Science 229(6): 1041-1059.
  • Gerolymos, N. and G. Gazetas (2006). "Winkler model for lateral response of rigid caisson foundations in linear soil." Soil Dynamics and Earthquake Engineering 26(5): 347-361.
  • Ghaitani, M. M. and H. A. Esmaeili (2013). "Vibration response of oil pipes embedded in viscous fluid using DQM".2nd International Conference on Emerging Trends in Engineering & Technology. College of Engineering, Teerthanker Mahaveer University.
  • Hu, J.-S., W.-Y. Sun and J. Zhou (2011). "Vibration analysis and crack identification of a cantilever pipe with a circumferential part-through crack." Journal of Vibration and Shock 4: 008.
  • Ibrahim, R. A. (2010). "Overview of Mechanics of Pipes Conveying Fluids-Part I:Fundamental Studies" Journal of Pressure Vessel Technology 132: 034001-034031.
  • Kargarnovin, M. H., D. Younesian, D. J. Thompson and C. J. C. Jones (2005). "Response of beams on nonlinear viscoelastic foundations to harmonic moving loads." Computers & Structures 83(23-24): 1865-1877.
  • Kumar, C. and V. Rastogi (2010). "A brief review on dynamics of a cracked rotor" International Journal of Rotating Machinery 2009.
  • Liu, D., H. Gurgenci and M. Veidt (2003). "Crack detection in hollow section structures through coupled response measurements." Journal of sound and vibration 261(1): 17-29.
  • Liu, D., T. Zhu, J. Zhou And X.-G. Zhao (2012). "Study on Vibration Characteristics of a Pipe with a Crack" China Offshore Platform 2: 007.
  • Mahjoob, M. J. and a. A. Shahsavari (2007). A Vibration-Based Damage Detection Method for Pipes Conveying Fluid. 48th AIAA/ ASME/ ASCE/ AHS/ ASC Structures, Structural Dynamics, and Materials Conference. Honolulu, Hawaii.
  • Murigendrappa, S. M., S. K. Maiti and H. R. Srirangarajan (2004). "Experimental and theoretical study on crack detection in pipes filled with fluid." Journal of Sound and Vibration 270(4-5): 1013-1032.
  • Murigendrappa, S. M., S. K. Maiti and H. R. Srirangarajan (2004). "Frequency-based experimental and theoretical identification of multiple cracks in straight pipes filled with fluid." NDT & E International 37(6): 431-438.
  • Naniwadekar, M. R., S. S. Naik and S. K. Maiti (2008). "On prediction of crack in different orientations in pipe using frequency based approach." Mechanical Systems and Signal Processing 22(3): 693-708.
  • Nguyen, V.-H. and D. Duhamel (2008). "Finite element procedures for nonlinear structures in moving coordinates.Part II: Infinite beam under moving harmonic loads." Computers & Structures 86(21-22): 2056-2063.
  • Païdoussis, M. P. (1998). Fluid-Structure Interactions. Fluid-Structure Interactions, Academic Press.
  • Persidis, S. (2007). Mathematical Handbook. ESPI, ATHENS, ESPI PUBLISHING.
  • Sinha, J. K., M. I. Friswell and S. Edwards (2002). "Simplified Models for the Location of Cracks in Beam Structures Using Measured Vibration Data." Journal of Sound and Vibration 251(1): 13-38.
  • Son, I.-S., S.-P. Lee, J.-K. Lee and D.-S. Bae "Dynamic Behavior of Forced Vibration of Elastically Restrained Pipe Conveying Fluid with Crack and Concentrated Mass."
  • Son, I., J. Cho and H. Yoon (2007). "Effects of Attached Mass on Stability of Pipe Conveying Fluid with Crack." Transactions of KSNVE 17(10): 1002-1009.
  • Son, I. and H. Yoon (2008). "Dynamic Stability of Elastically Restrained Cantilever Pipe Conveying Fluid with Crack." Transactions of KSNVE 18(2): 177-184.
  • Son, I. S., H. I. Yoon, S. P. Lee and D. J. Kim (2010). "Effects of Tip Mass on Stability of Rotating Cantilever Pipe Conveying Fluid with Crack." International Journal of Modern Physics B 24(15n16): 2609-2614.
  • Streeter, V. L., E. B. Wylie and K. W. Bedford (1998). Fluid mechanics, WCB, McGraw-Hill.
  • Tada, H., P. C. Paris and G. R. Irwin (1973). Stress Analysis of Cracks Handbook. New York, American Society of Mechanical Engineers.
  • Xie, Y. (1998). "A theory on cracked pipe." International journal of pressure vessels and piping 75(12): 865-869.
  • Ye, J., Y. He, X. Chen, Z. Zhai, Y. Wang and Z. He (2010). "Pipe crack identification based on finite element method of second generation wavelets." Mechanical Systems and Signal Processing 24(2): 379-393.
  • Yoon, H.-I. and I.-S. Son (2004). "Dynamic behavior of cracked pipe conveying fluid with moving mass based on Timoshenko beam theory." KSME international journal 18(12): 2216-2224.
  • Yoon, H.-I. and I.-S. Son (2005). "Influence of tip mass on dynamic behavior of cracked cantilever pipe conveying fluid with moving mass." Journal of mechanical science and technology 19(9): 1731-1741.
  • Yoon, H.-I. and I.-S. Son (2006). "Dynamic behavior of cracked simply supported pipe conveying fluid with moving mass." Journal of Sound and Vibration 292(3-5): 941-953.
  • Yoon, H.-I. and I.-S. Son (2007). "Dynamic response of rotating flexible cantilever pipe conveying fluid with tip mass." International journal of mechanical sciences 49(7): 878-887.
  • Yoon, H., I. Son and S. Ahn (2007). "Free Vibration Analysis of Euler-Bernoulli Beam with Double Cracks." Journal of Mechanical Science and Technology 21: 476-485.
  • Yoon, H. I., I. S. Son and J. T. Jin (2004). "A study on dynamic behavior of simply supported fluid flow pipe with crack and moving mass."
  • Younesian, D., S. R. Marjani and E. Esmailzadeh (2012). "Nonlinear vibration analysis of harmonically excited cracked beams on viscoelastic foundations." Nonlinear Dynamics 71(1-2): 109-120.
  • Yumin, H., Y. Junjie and H. Zhengjia (2011). "Research on Quantitative Identification of Pipe Crack Based on Vibration." Measuring Technology and Mechatronics Automation (ICMTMA), 2011 Third International Conference on, IEEE.
  • Zou, J., J. Chen, J. C. Niu and Z. M. Geng (2003). "Discussion on the local flexibility due to the crack in a cracked rotor system." Journal of Sound and Vibration 262(2): 365-369.

Publication Dates

  • Publication in this collection
    Jan 2016

History

  • Received
    15 Mar 2015
  • Reviewed
    29 Sept 2015
  • Accepted
    30 Sept 2015
Individual owner www.lajss.org - São Paulo - SP - Brazil
E-mail: lajsssecretary@gmsie.usp.br