ISSN: 2641-3078
Annals of Limnology and Oceanography
Research Article       Open aaacess      Peer-Reviewed

Corrosion Detection in Pipelines Based on Measurement of Natural Frequencies

Zhenghao Yang and Erkan Oterkus*

University of Strathclyde, UK
*Corresponding author: Dr. Erkan Oterkus, 100 Montrose Street Glasgow, University of Strathclyde, G4 0LZ United Kingdom, Tel: +44-771-470-3872; E-mail:
Received: 31 December, 2016 | Accepted: 21 January, 2017 | Published: 24 January, 2017
Keywords: Corrosion detection; Pipeline; Rayleigh’s law; Natural frequency

Cite this as

Yang Z, Oterkus E (2017) Corrosion Detection in Pipelines Based on Measurement of Natural Frequencies. Ann Limnol Oceanogr 2(1): 001-006. DOI: 10.17352/alo.000004

Natural frequency of structure mainly depends on mass and stiffness. Stiffness is bound to change after structural damage. Hence, natural frequency starts to decline.

This study presents a new method to determine the location and degree of the corrosion damage by measuring the natural frequencies of the damaged pipeline. With this method, only measurement of the first and second natural frequencies of damaged pipeline is required. The formulation is based on Rayleigh’s Law to determine a relationship between the degree of damage and damage location. The formulation is validated by comparing against beam and solid finite element models.


Pipeline transport is important in modern industry, and pipelines are widely used in the fields of petroleum, natural gas, coal gas, chemical fluid, water, coal, etc. They are especially common in the petrochemical and natural gas industries. Furthermore, urban water supplies and gas systems, which are an integral element of modern society, rely on even larger pipeline networks. Various factors, including pipeline corrosion, external force and loose joint parts, may lead to pipeline leaks. Such leaks result, not only in resource consumption and economic loss, but also serious environmental pollution. In order to avoid pipeline accidents, there is a distinct requirement for systems by which industrial pipelines can be monitored on a real-time basis and safety assessments and durability predictions can be produced. The two main challenges associated with the online monitoring of pipeline operations are the detection of pipeline defects and the accurate identification of the location of the pipeline damage.

Structural damage results from a variety of different factors including operation overload, impingement, cracks, corrosion, strains, production defects, etc. These deficiencies will often lead to changes in the physical properties (rigidity, mass, damping) of a structure and will be accompanied by changes in its dynamic behaviours. This fact is commonly acknowledged in processes that involve monitoring the health of engineering structures and identifying any damage. Pipeline systems, which are part of large-scale structural engineering applications, are a fundamental element of lifeline engineering; therefore, significant attention has been invested in processes by which any damage to a given pipeline structure can be rapidly and accurately detected. When the operation of a pipeline structure exceeds a certain amount of time, its structure can be measured dynamically with the help of vibration testing. The dynamic behaviours of the data obtained can subsequently be employed to assess the overall health of the pipeline structure. Also, this form of detection, which is based on vibration features, has already been applied to assess various other structures including bridges, architecture, etc. Along with the development of modern sensor technology, microelectronics and computer technology, data collection, transmission, real-time analysis and processing technology have also been used widely. Vibration testing technology represents an automated and miniaturized detection system and, as a result, the detection of structural damage through the assessment of dynamic behaviours has consistently been a hot topic in the international academic community as well as within engineering circles. Coupled with the development of Finite Element Theory, it is anticipated that damage detection based on vibration technology will play a much more prominent role in the future [1-3].

At present, vibration methods are widely used to detect the corrosion of engineering structures. Some of the existing techniques are based on studying changes in vibration mode shape or damping factors. For example, Vitaly and Vladislav [4], developed a method of monitoring the vibration distortions that could be observed in a system’s symmetry as a result of the development of defects. Panteliou [5] and Kyriazoglou [6], investigated the effects of vibration damping factors on crack detection. Another approach to corrosion detection is based on monitoring changes in resonant frequencies. For example, Liang [7] and Barad [9], employed a natural frequency method by starting with a differential equation of free vibration [4-10].

This study describes a new corrosion detection method that involves the measurement of natural frequencies using an approach that is based on Rayleigh’s Law. This method is sampled on a pipeline and employed to assess the relationship between the location of the corrosion and degree of corrosion by adopting sectional integration and integral transformation. The main advantage of this approach is its simplicity while allowing a more accurate detection of the location of damage than the traditional natural frequency method.

Procedures for the development of the corrosion detection method

The use of Rayleigh’s Law to calculate natural frequency is based on the principle of the conservation of energy. When an elastic body is free of vibration without damping, its total mechanical energy (sum of strain energy and kinetic energy) remains constant.

For a distributed mass beam with a uniform section, its displacement function can be given by:

y( x,t )=Y( x )sin(ωt+α)          (1) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaGqacabaaaaaaaaapeGaa8xEamaabmaapaqaa8qacaWF4bGaaiilaiaa=rhaaiaawIcacaGLPaaacqGH9aqpcaWFzbWaaeWaa8aabaWdbiaa=HhaaiaawIcacaGLPaaacaqGZbGaaeyAaiaab6gacaGGOaqedmvETj2BSbacfiGae4xYdCNaa8hDaiabgUcaRiab+f7aHjaacMcacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeikaiaabgdacaqGPaaaaa@5A7F@

where Y( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiaadMfadaqadaqaaiaadIhaaiaawIcacaGLPaaaaaa@3F92@ is the amplitude, ω is the natural frequency, and ω and α MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiabeg7aHbaa@3DCD@ depend on initial condition. Natural frequency is the frequency at which the system oscillates without any external loading. Natural frequency of the structure mainly depends on mass and stiffness. To measure the natural frequency of a pipeline, a pulse analyzer system can be equipped with accelerometer sensors, which are used to conduct the structural dynamic vibration response for monitoring system and data acquisition [11].

The strain energy, U , of the free vibration beam can be expressed as:

U= 1 2 0 L EI [ y(x)'' ] 2 dx= 1 2 si n 2 (ωt+α) 0 L EI [ Y''(x) ] 2 dx           (2) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaGqacabaaaaaaaaapeGaa8xvaiabg2da9maalaaapaqaa8qacaaIXaaapaqaa8qacaaIYaaaamaawahabeWcpaqaa8qacaaIWaaapaqaa8qacaWFmbaan8aabaWdbiabgUIiYdaakiaa=veacaWFjbWaamWaa8aabaWdbiaa=LhacaGGOaGaa8hEaiaacMcacaWFNaGaa83jaaGaay5waiaaw2faa8aadaahaaWcbeqaa8qacaaIYaaaaOGaa8hzaiaa=HhacqGH9aqpdaWcaaWdaeaapeGaaGymaaWdaeaapeGaaGOmaaaacaWFZbGaa8xAaiaa=5gapaWaaWbaaSqabeaapeGaaGOmaaaakiaacIcarmWu51MyVXgaiuGacqGFjpWDcaWF0bGaey4kaSIae4xSdeMaaiykamaawahabeWcpaqaa8qacaaIWaaapaqaa8qacaWFmbaan8aabaWdbiabgUIiYdaakiaa=veacaWFjbWaamWaa8aabaWdbiaa=LfacaWFNaGaa83jaiaacIcacaWF4bGaaiykaaGaay5waiaaw2faa8aadaahaaWcbeqaa8qacaaIYaaaaOGaa8hzaiaa=HhacaWFGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqGYaGaaeykaaaa@77D0@

where 1 is the length of the beam, E is the elastic modulus, andis the moment of inertia.

Moreover, kinetic energy, T can be expressed as:

T= 1 2 0 L ρA [ y(x,t) . ] 2 dx= 1 2 ω 2 co s 2 (ωt+α) 0 L ρA [ Y(x) ] 2 dx          (3) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaGqacabaaaaaaaaapeGaa8hvaiabg2da9maalaaapaqaa8qacaaIXaaapaqaa8qacaaIYaaaamaawahabeWcpaqaa8qacaaIWaaapaqaa8qacaWFGaGaa8hiaiaa=bcacaWFGaGaa8hiaiaa=bcacaWFGaGaa8hiaiaa=bcacaWFGaGaa8hiaiaa=bcacaWFmbaan8aabaWdbiabgUIiYdaarmWu51MyVXgaiuGakiab+f8aYjaa=feadaWadaWdaeaadaWfGaqaa8qacaWF5bGaaiikaiaa=HhacaGGSaGaa8hDaiaacMcaaSWdaeqabaWdbiaac6caaaaakiaawUfacaGLDbaapaWaaWbaaSqabeaapeGaaGOmaaaakiaa=rgacaWF4bGaeyypa0ZaaSaaa8aabaWdbiaaigdaa8aabaWdbiaaikdaaaGae4xYdC3damaaCaaaleqabaWdbiaaikdaaaGccaWFJbGaa83Baiaa=nhapaWaaWbaaSqabeaapeGaaGOmaaaakiaacIcacqGFjpWDcaWF0bGaey4kaSIae4xSdeMaaiykamaawahabeWcpaqaa8qacaaIWaaapaqaa8qacaWFGaGaa8hiaiaa=bcacaWFGaGaa8hiaiaa=bcacaWFGaGaa8hiaiaa=bcacaWFGaGaa8hiaiaa=Xeaa0WdaeaapeGaey4kIipaaOGae4xWdiNaa8xqamaadmaapaqaa8qacaWFzbGaaiikaiaa=HhacaGGPaaacaGLBbGaayzxaaWdamaaCaaaleqabaWdbiaaikdaaaGccaWFKbGaa8hEaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaae4maiaabMcaaaa@8A7B@

The strain energy reaches a peak value when the beam has the maximum amplitude; meanwhile, the kinetic energy is 0. When the amplitude is 0, the kinetic energy has the maximum value; meanwhile, the strain energy is 0. The maximum strain energy, U max MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiaadwfadaWgaaWcbaGaciyBaiaacggacaGG4baabeaaaaa@4008@ , and kinetic energy, T max MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiaadsfadaWgaaWcbaGaciyBaiaacggacaGG4baabeaaaaa@4007@ , can be respectively given as:

U max = 1 2 0 L EI [ Y '' ( x )] 2 dx          (4a) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaGqacabaaaaaaaaapeGaa8xva8aadaWgaaWcbaWdbiaa=1gacaWFHbGaa8hEaaWdaeqaaOWdbiabg2da9maalaaapaqaa8qacaaIXaaapaqaa8qacaaIYaaaamaawahabeWcpaqaa8qacaaIWaaapaqaa8qacaWFmbaan8aabaWdbiabgUIiYdaakiaa=veacaWFjbGaai4waiaa=LfapaWaaWbaaSqabeaapeGaa83jaiaa=DcaaaGcdaqadaWdaeaapeGaa8hEaaGaayjkaiaawMcaaiaac2fapaWaaWbaaSqabeaapeGaaGOmaaaakiaa=rgacaWF4bGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqG0aGaaeyyaiaabMcaaaa@5C47@


T max = 1 2 ω 2 0 L m [Y( x )] 2 dx          (4b) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaGqacabaaaaaaaaapeGaa8hva8aadaWgaaWcbaWdbiaa=1gacaWFHbGaa8hEaaWdaeqaaOWdbiabg2da9maalaaapaqaa8qacaaIXaaapaqaa8qacaaIYaaaaeXatLxBI9gBaGqbciab+L8a39aadaahaaWcbeqaa8qacaaIYaaaaOWaaybCaeqal8aabaWdbiaaicdaa8aabaWdbiaa=Xeaa0WdaeaapeGaey4kIipaaOWdamaaxacabaWdbiaa=1gaaSWdaeqabaWdbiabgkHiTaaakiaacUfacaWFzbWaaeWaa8aabaWdbiaa=HhaaiaawIcacaGLPaaacaGGDbWdamaaCaaaleqabaWdbiaaikdaaaGccaWFKbGaa8hEaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaaeinaiaabkgacaqGPaaaaa@60EF@

where m ¯ MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiqad2gagaqeaaaa@3D38@ is the mass per unit length.

According to the principle of conservation of energy, we have U max = T max MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiaadwfadaWgaaWcbaGaciyBaiaacggacaGG4baabeaakiabg2da9iaadsfadaWgaaWcbaGaciyBaiaacggacaGG4baabeaaaaa@44F1@ , hence:

ω 2 = 0 Λ ΕΙ [ Ψ ( ξ )] 2 δξ 0 Λ μ [Ψ( ξ )] 2 δξ           (5) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaeXatLxBI9gBaGqbcabaaaaaaaaapeGae8xYdC3damaaCaaaleqabaWdbiab=jdaYaaakiab=1da9maalaaapaqaa8qadaqfWaqabSWdaeaapeGae8hmaadapaqaa8qacqWFBoata0WdaeaapeGae83kIipaaOGae8xLduKae8xMdKKae83waSLae8hQdK1damaaCaaaleqabaWdbiab=1GiNiab=1GiNaaakmaabmaapaqaa8qacqWF+oaEaiaawIcacaGLPaaacqWFDbqxpaWaaWbaaSqabeaapeGae8NmaidaaOGae8hTdqMae8NVdGhapaqaa8qadaqfWaqabSWdaeaapeGae8hmaadapaqaa8qacqWFBoata0WdaeaapeGae83kIipaaOWdamaaxacabaWdbiab=X7aTbWcpaqabeaapeGae8NeI0caaOGae83waSLae8hQdK1aaeWaa8aabaWdbiab=57a4bGaayjkaiaawMcaaiab=1faD9aadaahaaWcbeqaa8qacqWFYaGmaaGccqWF0oazcqWF+oaEaaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqG1aGaaeykaaaa@757C@

When considering a corroded pipeline, it can be assumed that the corrosion pit decreases the bending stiffness at the damaged position of the pipe, as shown in Figure 1,

where EI and (EI)’ are the bending stiffnesses of the intact and corroded elements of the pipe, respectively, L is the length of the pipe, x is the distance between the midpoint of the corroded segment to the left end of the pipeline, and D is the length of the corroded segment. If we neglect the mass loss due to the corrosion and apply these parameters into Eq. (5), the natural frequency of the corroded pipe can be obtained as follows:

ω c 2 = 0 xD/2 EI [ Y''(x) ] 2 dx+ xD/2 x+D/2 ( EI ) ' [ Y''(x) ] 2 dx+ x+D/2 L EI [ Y''(x) ] 2 dx 0 L m [ Y(x) ] 2 dx           (6) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaeXatLxBI9gBaGqbcabaaaaaaaaapeGae8xYdC3damaaBaaaleaaieGapeGaa43yaaWdaeqaaOWaaWbaaSqabeaapeGaaGOmaaaakiabg2da9maalaaapaqaa8qadaqfWaqabSWdaeaapeGaaGimaaWdaeaapeGaa4hEaiabgkHiTiaa+reacaGGVaGaaGOmaaqdpaqaa8qacqGHRiI8aaGccaGFfbGaa4xsamaadmaapaqaa8qacaGFzbGaa43jaiaa+DcacaGGOaGaa4hEaiaacMcaaiaawUfacaGLDbaapaWaaWbaaSqabeaapeGaaGOmaaaakiaa+rgacaGF4bGaey4kaSYaaubmaeqal8aabaWdbiaa+HhacqGHsislcaGFebGaai4laiaaikdaa8aabaWdbiaa+HhacqGHRaWkcaGFebGaai4laiaaikdaa0WdaeaapeGaey4kIipaaOWaaeWaa8aabaWdbiaa+veacaGFjbaacaGLOaGaayzkaaWdamaaCaaaleqabaWdbiaa+DcaaaGcdaWadaWdaeaapeGaa4xwaiaa+DcacaGFNaGaaiikaiaa+HhacaGGPaaacaGLBbGaayzxaaWdamaaCaaaleqabaWdbiaaikdaaaGccaGFKbGaa4hEaiabgUcaRmaavadabeWcpaqaa8qacaGF4bGaey4kaSIaa4hraiaac+cacaaIYaaapaqaa8qacaGFmbaan8aabaWdbiabgUIiYdaakiaa+veacaGFjbWaamWaa8aabaWdbiaa+LfacaGFNaGaa43jaiaacIcacaGF4bGaaiykaaGaay5waiaaw2faa8aadaahaaWcbeqaa8qacaaIYaaaaOGaa4hzaiaa+Hhaa8aabaWdbmaavadabeWcpaqaa8qacaaIWaaapaqaa8qacaGFmbaan8aabaWdbiabgUIiYdaak8aadaWfGaqaa8qacaGFTbaal8aabeqaa8qacqGHsislaaGcdaWadaWdaeaapeGaa4xwaiaacIcacaGF4bGaaiykaaGaay5waiaaw2faa8aadaahaaWcbeqaa8qacaaIYaaaaOGaa4hzaiaa+HhaaaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqG2aGaaeykaaaa@9A84@

If we introduce a factor, c, to describe the degree of the corrosion pit, where, c=( EI )'/( EI ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiaadogacqGH9aqpdaqadaqaaiaadweacaWGjbaacaGLOaGaayzkaaGaai4jaiaac+cadaqadaqaaiaadweacaWGjbaacaGLOaGaayzkaaaaaa@45BC@ Eq. (6) can be rewritten as:

ω c 2 = 0 L EI [ Y(x)'' ] 2 dx xD/2 x+D/2 ( 1c )EI [ Y( x )'' ] 2 dx 0 L m [Y(x)] 2 dx            (7) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaeXatLxBI9gBaGqbcabaaaaaaaaapeGae8xYdC3damaaBaaaleaaieGapeGaa43yaaWdaeqaaOWaaWbaaSqabeaapeGaaGOmaaaakiabg2da9maalaaapaqaa8qadaqfWaqabSWdaeaapeGaaGimaaWdaeaapeGaa4htaaqdpaqaa8qacqGHRiI8aaGccaGFfbGaa4xsamaadmaapaqaa8qacaGFzbGaaiikaiaa+HhacaGGPaGaa43jaiaa+DcaaiaawUfacaGLDbaapaWaaWbaaSqabeaapeGaaGOmaaaakiaa+rgacaGF4bGaeyOeI0Yaaubmaeqal8aabaWdbiaa+HhacqGHsislcaGFebGaai4laiaaikdaa8aabaWdbiaa+HhacqGHRaWkcaGFebGaai4laiaaikdaa0WdaeaapeGaey4kIipaaOWaaeWaa8aabaWdbiaaigdacqGHsislcaGFJbaacaGLOaGaayzkaaGaa4xraiaa+LeadaWadaWdaeaapeGaa4xwamaabmaapaqaa8qacaGF4baacaGLOaGaayzkaaGaa43jaiaa+DcaaiaawUfacaGLDbaapaWaaWbaaSqabeaapeGaaGOmaaaakiaa+rgacaGF4baapaqaa8qadaqfWaqabSWdaeaapeGaaGimaaWdaeaapeGaa4htaaqdpaqaa8qacqGHRiI8aaGcpaWaaCbiaeaapeGaa4xBaaWcpaqabeaapeGaeyOeI0caaOGaai4waiaa+LfacaGGOaGaa4hEaiaacMcacaGGDbWdamaaCaaaleqabaWdbiaaikdaaaGccaGFKbGaa4hEaaaacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqG3aGaaeykaaaa@861C@


ω c 2 = ω I 2 ( 1c )EI m xD/2 x+D/2 [Y(x)''] 2 dx 0 L [Y( x )] 2 dx            (8) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaeXatLxBI9gBaGqbcabaaaaaaaaapeGae8xYdC3damaaBaaaleaaieGapeGaa43yaaWdaeqaaOWaaWbaaSqabeaapeGaaGOmaaaakiabg2da9iab=L8a39aadaqhaaWcbaWdbiaa+Leaa8aabaWdbiaaikdaaaGccqGHsisldaWcaaWdaeaapeWaaeWaa8aabaWdbiaaigdacqGHsislcaGFJbaacaGLOaGaayzkaaGaa4xraiaa+Leaa8aabaWaaCbiaeaapeGaa4xBaaWcpaqabeaapeGaeyOeI0caaaaakiabgwSixpaalaaapaqaa8qadaqfWaqabSWdaeaapeGaa4hEaiabgkHiTiaa+reacaGGVaGaaGOmaaWdaeaapeGaa4hEaiabgUcaRiaa+reacaGGVaGaaGOmaaqdpaqaa8qacqGHRiI8aaGccaGGBbGaa4xwaiaacIcacaGF4bGaaiykaiaa+DcacaGFNaGaaiyxa8aadaahaaWcbeqaa8qacaaIYaaaaOGaa4hzaiaa+Hhaa8aabaWdbmaavadabeWcpaqaa8qacaaIWaaapaqaa8qacaGFmbaan8aabaWdbiabgUIiYdaakiaacUfacaGFzbWaaeWaa8aabaWdbiaa+HhaaiaawIcacaGLPaaacaGGDbWdamaaCaaaleqabaWdbiaaikdaaaGccaGFKbGaa4hEaaaacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqG4aGaaeykaaaa@7CE5@

where ω I MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaiabeM8a3naaBaaaleaacaWGjbaabeaaaaa@3EF5@ is the natural frequency of the intact pipe with the same dimensions. Then, the degree of corrosion, c, can be represented by the following expression:

c= (EI)' EI =1 m ( ω In 2 ω Cn 2 ) EI 0 L [ Y n ( x )] 2 dx xD/2 x+D/2 [ Y n (x)''] 2            (9) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaapeGaae4yaiabg2da9maalaaapaqaa8qacaGGOaacbiGaa8xraiaa=LeacaGGPaGaa83jaaWdaeaapeGaa8xraiaa=LeaaaGaeyypa0JaaGymaiabgkHiTmaalaaapaqaamaaxacabaWdbiaa=1gaaSWdaeqabaWdbiabgkHiTaaakiaacIcarmWu51MyVXgaiuGacqGFjpWDpaWaa0baaSqaa8qacaWFjbGaa8NBaaWdaeaapeGaaGOmaaaakiabgkHiTiad0bOFjpWDpaWaa0baaSqaa8qacaWFdbGaa8NBaaWdaeaapeGaaGOmaaaakiaacMcaa8aabaWdbiaa=veacaWFjbaaaiabgwSixpaalaaapaqaa8qadaqfWaqabSWdaeaapeGaaGimaaWdaeaapeGaa8htaaqdpaqaa8qacqGHRiI8aaGccaGGBbGaa8xwa8aadaWgaaWcbaWdbiaa=5gaa8aabeaak8qadaqadaWdaeaapeGaa8hEaaGaayjkaiaawMcaaiaac2fapaWaaWbaaSqabeaapeGaaGOmaaaakiaa=rgacaWF4baapaqaa8qadaqfWaqabSWdaeaapeGaa8hEaiabgkHiTiaa=reacaGGVaGaaGOmaaWdaeaapeGaa8hEaiabgUcaRiaa=reacaGGVaGaaGOmaaqdpaqaa8qacqGHRiI8aaGccaGGBbGaa8xwa8aadaWgaaWcbaWdbiaa=5gaa8aabeaak8qacaGGOaGaa8hEaiaacMcacaWFNaGaa83jaiaac2fapaWaaWbaaSqabeaapeGaaGOmaaaaaaGccaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqG5aGaaeykaaaa@868D@

Here ‘n’ is the nth vibration mode, ωIn and ωCn are the nth natural frequenies of the intact and corroded pipe, respectively, and Y n (x) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaceaaymacbiaeaaaaaaaaa8qacaWFzbWdamaaBaaaleaapeGaa8NBaaWdaeqaaOWdbiaacIcacaWF4bGaaiykaaaa@4106@ is the nth vibration mode function. The research presented in this study is based on examples of a simply supported constrain and the first and second mode of vibration for this boundary condition is shown in Figure 2 [12].

Thus, the vibration mode functions of the first two modes can be expressed as:

Mode 1:            Y 1 ( x )= a 1 sin( πx L )          (10a) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaapeGaaeytaiaab+gacaqGKbGaaeyzaiaabckacaaIXaGaaiOoaiaabckacaqGGcGaaeiOaiaabckacaqGGcGaaeiOaiaabckacaqGGcGaaeiOaiaabckacaqGGcacbiGaa8xwa8aadaWgaaWcbaWdbiaaigdaa8aabeaak8qadaqadaWdaeaapeGaa8hEaaGaayjkaiaawMcaaiabg2da9iaa=fgapaWaaSbaaSqaa8qacaaIXaaapaqabaGcpeGaci4CaiaacMgacaGGUbWaaeWaa8aabaWdbmaalaaapaqaaeXatLxBI9gBaGqbc8qacqGFapaCcaWF4baapaqaa8qacaWFmbaaaaGaayjkaiaawMcaaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaaeymaiaabcdacaqGHbGaaeykaaaa@6BAC@

Mode 2:            Y 2 ( x )= a 2 sin( 2πx L )          (10b) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaapeGaaeytaiaab+gacaqGKbGaaeyzaiaabckacaaIYaGaaiOoaiaabckacaqGGcGaaeiOaiaabckacaqGGcGaaeiOaiaabckacaqGGcGaaeiOaiaabckacaqGGcacbiGaa8xwa8aadaWgaaWcbaWdbiaaikdaa8aabeaak8qadaqadaWdaeaapeGaa8hEaaGaayjkaiaawMcaaiabg2da9iaa=fgapaWaaSbaaSqaa8qacaaIYaaapaqabaGcpeGaae4CaiaabMgacaqGUbGaaiikamaalaaapaqaa8qacaaIYaqedmvETj2BSbacfiGae4hWdaNaa8hEaaWdaeaapeGaa8htaaaacaGGPaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqGXaGaaeimaiaabkgacaqGPaaaaa@6C18@

If corrosive length, D, can be estimated to a reasonable level, by combining Eqs. (9), (10a) and (10b), the extent of corrosion factor, c, and corrosion location, x, can be predicted in consideration of the natural frequency of the intact pipeline ωI, and the corroded pipeline ωc.

For example, a corroded pipeline with simply supported constrains has the following properties: pipe length L = 10 m, Young’s modulus E = 210 GPa, mass density ρ = 7850 kg/m3, 0.25 m of outer radius, and 0.02m of thickness. Its first and second natural frequencies were measured as ωc1 =53.89rad/s and ωc2 =213.69rad/s, and its corresponding natural frequencies of intact pipeline were calculated using Eq. (5) as: ωI1 =55.81Rad/s and ωI2 =221.30Rad/s, and corrosion length was assumed to be D = 0.50 m, then the variation of c with corroded region position x was plotted and is presented in Figure 3. Since physically there is only one unique c value to represent the degree of corrosion at the damaged locations, the potential corrosion region positions are denoted by the intersection points of the curves (Figure 3), and the numerical value of c at the intersection points represents the degree of pipeline corrosion. Due to the symmetry of this pipeline structure, there are two intersection points at about x =3. 30 m and x = 6.70 m and the actual corrosion region position is one of them.

Finite element analyses

In order to investigate and verify the accuracy of the pipeline damage detection through the use of natural frequency measurements, a commercial finite element program, ANSYS, was employed.

A simply supported corroded pipeline was modelled using Beam188 element with the following characteristics: Pipeline length, L=10 m outer radius, r = 0.16 m wall thickness, t = 0.01m mass density ρ=7850 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaeXatLxBI9gBaGqbaabaaaaaaaaapeGae8xWdiNaeyypa0JaaG4naiaaiIdacaaI1aGaaGimaaaa@44A2@ kg/m3, elastic modulus of non-corroded segment, E = 210 GPa elastic modulus of corroded segment E’ = c . E. Here, c represents the degree of corrosion, ranging from 0 to 1. The model was discretized with 20 elements along the axis, as shown in Figure 4. Five simulation cases were analysed with five different corrosion positions and corrosion levels. Their first and second natural frequencies, ωC1 and ωC2 , were computed individually using ANSYS, as shown in Table 1.

Based on the first and second natural frequencies that were computed using ANSYS ( 1) and the corresponding natural frequencies of the intact pipeline, ω I1 =55.81 rad/s,   ω I2 =221.30 rad/s MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaeXatLxBI9gBaGqbcabaaaaaaaaapeGae8xYdC3damaaBaaaleaaieGapeGaa4xsaiaaigdaa8aabeaak8qacqGH9aqpcaaI1aGaaGynaiaac6cacaaI4aGaaGymaiaa+bkacaqGYbGaaeyyaiaabsgacaGGVaGaae4CaiaacYcacaGFGcGaaeiOaiab=L8a39aadaWgaaWcbaWdbiaa+LeacaaIYaaapaqabaGcpeGaeyypa0JaaGOmaiaaikdacaaIXaGaaiOlaiaaiodacaaIWaGaa4hOaiaabkhacaqGHbGaaeizaiaac+cacaqGZbaaaa@5E8D@ , Eq. (9) can be solved to determine the variation of (EI’/EI) versus the location of the corrosion for the first two modes. The results of Cases 1 to 5 are shown in Figure 5, which depict the location of the corrosion and the degree of corrosion at the intersection points. There are two intersection points in each graph and this is due to the symmetry of the structure. The solution and plots were calculated using Mathcad.

As can be seen clearly from Table 2, the calculations in terms of the location of the corrosion and the degree of corrosion were accurately predicted.

It is noteworthy that, from the characteristic Eq. (9), there is an important parameter, the effect of estimated corrosion length ‘D’, which should not be neglected. The FEA described above reveals that, ideally, if an estimated ‘D’ is merely equal to the practical corrosion length, both corrosion location and degree of corrosion can be accurately predicted. However, in practical engineering applications, it is almost impossible to accurately estimate corrosion length. Thus, it is also important to take into consideration the estimated corrosion length, ‘D’, in order to investigate how this affects the prediction.

A further finite element analysis was carried to investigate the influence of ‘D’. The finite element model was largely based on Case 3 of the previous problem. The only changes were that the corrosion segment was increased to 1 m, and the midpoint of the corrosion was positioned at x = 3.50 m, as shown in Figure 6.

The first two natural frequencies of this corroded pipeline were computed using ANSYS, and were ω C1 =51.86 rad/s MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaGqacabaaaaaaaaapeGaa8xYd8aadaWgaaWcbaWdbiaa=neacaaIXaaapaqabaGcpeGaeyypa0JaaGynaiaaigdacaGGUaGaaGioaiaaiAdacaWFGcGaaeOCaiaabggacaqGKbGaai4laiaabohaaaa@49D7@ and ω C2 =209.93 rad/s MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaiKc9yrVq0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabaGaciaacaqabeaadaqaaqaaaOqaaGqacabaaaaaaaaapeGaa8xYd8aadaWgaaWcbaWdbiaa=neacaaIYaaapaqabaGcpeGaeyypa0JaaGOmaiaaicdacaaI5aGaaiOlaiaaiMdacaaIZaGaa8hOaiaabkhacaqGHbGaaeizaiaac+cacaqGZbaaaa@4A95@ . A total of five cases of varying values of ‘D’ were incorporated into the characteristic, Eq. (9). These were D = 0.50 m, 0.75 m, 1 m, 1.25 m and 1.5 m. The results of Cases 1 to 5 were calculated using Mathcad, as shown in Figure 7, and detailed data is presented in Table 3.

As can be clearly seen from Table 3, in all cases, the corrosion position were predicted correctly, even though the predicted corrosion length, D, was smaller or greater than the actual corrosion length. However, over or under estimating D will affect the accuracy of the predicted degree of corrosion.

Finite element analyses of other geometries

In order to expand the investigation of the current corrosion detection method, a 3D corrosive pipeline model was studied using ANSYS, and this model was deemed to be more reflective of reality.

The finite element model was designed with a 20-node isoparametric Solid186 element. The length was 10 m, the inner and outer radius of the pipe were RI = 0.18 m and RO = 0.2 m, the elastic modulus was E=210 GPa, mass density ρ = 7850 kg/m3, and the Poisson ratio ν = 0.3. There were 24 cells along the circumference, 5 cells along the thickness, and 80 along the axis. The corrosion pit was idealized into a rectangular shape, and a typical finite element mesh is shown in Figure 8a.

Four different parameters of,x lc , dp and () were considered in the FEA. Here, x represents the location of the corrosion in terms of the distance from the left end of the pipeline to the midpoint of the dp corrosion pit; represents the length of the corrosion pit; lc dp represents the depth of the corrosion pit, anddescribes the width of the corrosion pit. As shown in Figure 8b,c, all four cases were studied by varying the above parameters as summarized in Table 4.

The constraints of the two ends of the finite element model were set to be simply supported. For all cases, the first and second natural frequency, ωC1and ωC2, were computed with the help of ANSYS, and the corresponding dimensional intact pipeline model was also created by using ANSYS to compute its first and second natural frequency, ωI1and ωI2. The details of FEA results are summarized in Table 5.

The computed natural frequencies are given in Table 5 by using an estimated corrosion length of D = 0.5 m. This information was plugged into Eq. (9) to calculate the variation of corrosion degree (EI)’/(EI) versus the corrosion pit location. Evaluated results are presented in Figure 9. The intersection points represent the corrosion pit location and the corrosion degree. Again, due to the symmetry, there are two potential corrosion locations for Case 1-4, and the actual location is at one of them.

The numerical results for all cases are shown in Table 6. This data clearly indicates that the predictions of the corrosion pit location are pretty accurate and that there is a good agreement between the finite element model and the proposed method.

Moreover, this FEA more closely resembles reality because the corrosion is modelled by creating a pit on the wall of the pipeline rather than decreasing the flexural rigidity of a segment. However, it can be simply estimated based on common sense that the larger and deeper the corrosion pit is, the more serious the degree of corrosion. According to the corrosion pit dimensions of each case, the degree of corrosion sequence from high to low is Case 2 < Case1 < Case 3 < Case 4, and this rank is reflected by the given results.

Discussion and Conclusions

This study presents a new method for detecting corrosion damage within a pipeline structure that is based on measurements of first and second natural frequency. The method is based on Rayleigh’s Law to determine a relationship between c (degree of corrosion) and corrosion location. Since only one unique value of c is physically permissible at a given corrosion location, the intersections of the various values of c at various natural frequencies along the axial direction of the pipeline provided insights into the location of the corrosion damage.

Two different types of finite element models were employed to verify the feasibility of this method: one considering the corrosion segment as a decreased bending stiffness on a uniform section pipeline, and the other using 3D Solid186 element within which the corrosion pit was modelled directly on the wall of the pipeline. More specifically, for the first type of model, the corrosion segment was created by decreasing its elastic modulus E, and the second model type by decreasing the moment of inertia. Both supported the application and accuracy of this method.

Note that an imprecise estimated value of corrosion length D does not affect the accuracy of the prediction of the corrosion location. However, it will lead to an inaccurate prediction of degree of corrosion. In this case, after determining the location of the corrosion, the length of the corrosion region should be measured in the pipeline and by using this actual length, Eq. (9) can provide an accurate degree of corrosion.

  1. Pandy AK, Biswas M, Samman MM (1991) Damage detection from changes in curvature mode shapes. Journal of Sound and Vibration 142: 321–332. Link:
  2. Yam LH, Leunga TP, Lib DB, Xueb KZ (1996) Theoretical and experimental study of model strain analysis. Journal of Sound and vibration 191: 251–260. Link:
  3. Salawa OS (1997) Detection of structural damage through changes in frequencies, A Review. Engineering Structures 19: 718-723. Link:
  4. Vitaly B, Valadislav Y (2013) Vibration methods of damage detection in initially symmetric structures. 569: 1116-1123. Link :
  5. Panteliou SD, Chondros TG, Argyrakis VC, Dimarogonas AD (2001) Damping factor as an indicator of crack severity. Journal of Sound and Vibration 241: 235-245. Link :
  6. Kyriazoglou C, Page BH, Guild FJ (2004) Vibration damping for crack detection in composite laminates. Composites Part A-Applied Science and Manufacturing 35: 945-953. Link:
  7. Liang RY, Choy FK, Hu J, Choy FK, Hu L (1991) Detection of Cracks in Beam structures using measurements of natural frequencies. 328: 505-518. Link:
  8. Hu J, Robert Liang Y (1993) An integrated approach to detection of cracks using vibration characteristics. 330: 841-853. Link:
  9. Barad KH, Sharma DS, Vishal Vyas (2013) Crack detection in cantilever beam by frequency based method. 51:770-775. Link:
  10. Bovsunovsky AP, Surace C (2005) Considerations regarding superharmonic vibrations of cracked beam and the variation in damping caused by the presence of the crack. Journal of Sound and Vibration 288: 865-886. Link:
  11. Mahmoodi R, Shahriari M, Zarghami R (2010) Natural frequencies behavior of pipeline system during LOCA in nuclear power plants, Proceedings of the World Congress on Engineering II, WCE, London, U.K. Link :
  12. Eaton DCG (1997) An Overview of Structural Acoustics and Related High-Frequency-Vibration Activities. Link:
© 2017 Yang Z, et al. This is an open-aaacess article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.