管道弯头对低频纵向导波传播特性影响分析

管道弯头显著改变了导波传播特性,影响了对检测信号的解读,研究弯头对导波传播特性的影响是实现复杂管道系统导波检测的基础。采用半解析有限元法计算弯管导波频散曲线,分析了弯管导波频散曲线所呈现的不同特征,并基于弯管导波频散曲线,以低频L(0,1)模态导波为研究对象,实验研究了低频L(0,1)模态导波通过管道弯头时的模态变换特征。研究结果发现,当L(0,1)模态导波通过管道弯头时,不仅会发生L(0,1)到F(1,1)的模态变换,还会模态变换出反向L(0,1)模态导波,即弯头反射现象,且随着激励频率的降低和弯头弯曲半径的减小,弯头反射现象愈发明显。研究结果将深化对弯管导波传播特性的认识,推动导波检测技术在复杂管道系统检测中的应用。      

杆中导波声弹性效应分析与数值模拟研究*

大型钢结构在役应力检测意义重大,基于超声导波声弹性效应进行应力检测具有潜在的优势。本文基于等效弹性常数法研究了杆中超声导波的声弹性效应。通过计算典型钢结构构件钢杆的频散曲线,确定了检测频率范围,对不同工作应力状态下L（0,1）、F（1,1）与T（0,1）模态的群速度值进行了理论分析与数值模拟。结果表明： L（0,1）模态较适合于钢杆轴力检测, F（1,1）、T（0,1）模态声弹性效应较弱,不适合于应力检测; L（0,1）模态的声弹性效应随着频率的增加而减弱;理论分析与有限元计算对声弹性效应的验证都与波结构轴向位移对声弹性效应的判断得到了较好的吻合。     

基于周向SH0导波的管道轴向缺陷定量表征*

随着管道服役时间的增加,其损伤逐渐累积,最终导致泄漏或爆炸事故。引入导波在线监测管道损伤是确保其安全运行的重要保障。周向零阶水平剪切(CSH0)波具有不易频散、对管道轴向缺陷敏感等优秀特性,非常适用于管道轴向缺陷的定量表征。该文 基于周向SH0导波在管道的传播特性,通过有限元模拟和试验对管道周向SH0模态导波与轴向缺陷的定量关系展开研究。结果表明：管道曲率半径越大,周向SH0模态的传播特性越接近板中SH0模态的传播特性;导波信号的反射系数和透射系数对轴向缺陷的长度和深度均呈现单调变化。依此建立了缺陷定量表征的关系云图,可用于评估缺陷尺寸。     

无限液体介质内管道轴对称纵向导波激发与传播特性研究

利用声-结构耦合有限元法,分别对轴对称分布径的向及轴向外力作用在无限液体介质内未充液及充液管道内壁所激发导波进行了模拟,并进一步利用短时傅里叶变换技术对瞬态波形作时频分析,在此基础上,结合色散及衰减曲线,探讨了外力分布特征对导波激发的影响。研究表明,无限液体介质内未充液管中L(0,2)模式及充液管中L(0,3)~L(0,4)模式皆具有高群速度、弱色散且弱衰减频带,适宜用于缺陷探测,而通过控制轴对称分布外力的频率,并令其沿轴向作用于管内壁或外壁,可实现上述导波模式的高效激发。      

周向超声导波对管道纵向缺陷检测的研究

在弹性动力学理论的基础上,研究了管道中周向超声导波的传播及其频散特性,并数值计算了周向导波的频散曲线。通过所建立的超声导波实验系统,在理论上与实验上研究了斜探头的入射角、频率及周向导波模态的关系,并在外径为88.8 mm、内径为80.8 mm的管道中激励出特定的单一模态周向导波;同时,利用单一模态周向导波对该管道表面上长25 mm、宽1 mm、深0.7 mm的人工缺陷进行了检测。结果表明,同在板中激励单一Lamb波模态类似,选用特定角度的斜探头与特定频率可在管道中激励出单一周向导波模态并有望应用于管道缺陷的检测。     

带粘弹性包覆层充液管道中的超声导波纵向模态

理论分析和实验研究了超声导波纵向模态在带粘弹性包覆层充液管道中的传播特性。得到了纵向模态的频散曲线。以此确定了适合带粘弹性包覆层充液管道缺陷检测的一定频带的纵向模态。经分析认为,频散小,衰减低的频带0～50 kHz的L(0,1)模态和未受干扰的L(0,2)模态分支部分,如频带170～210 kHz的L(0,4)模态,适合检测外直径25 mm,壁厚1．2 mm,外壁涂覆0．35 mm厚环氧树脂的充水钢管中的缺陷。而频散和衰减大,能量主要在水或环氧树脂粘弹性层中传播的纵向模态则不适合检测该类管道中的缺陷。     

液态金属快堆主管道中超声导波传播特性研究

讨论了钠冷快堆(Sodium-cooled Fast Reactor,SFR)主管道的整体温度和内部液态金属钠流动速度的变化对管道导波传播特性的影响。推导了充液管道中导波频散方程的一般形式,并给出了管道内液态金属钠处于流动状态下的导波频散方程。采用数值计算方法获得了管内液态金属钠处于不同温度和不同流速时的导波纵向模式频散曲线和导波时域波形。结果表明,温度变化对基阶纵向模式的影响较小,但对高阶纵向模式的影响较大;液态钠流速增大会使导波频散曲线向高频轻微移动,但在实际检测中可以忽路管内液体流动速度的影响。通过对时域接收波形的模拟计算,进一步考察了液态金属钠的温度及流动速度变化对导波传播的影响,并通过对比不同模态的激发特点和不同频段的导波时域波形特点,结合导波频散曲线,给出了适用于SFR管道超声无损检测的导波模态和声源激发频段选择方案。      

充粘液管材中超声纵向轴对称导波的频散特性分析

模拟分析充粘液管内径-壁厚比和超声纵向轴对称导波波包宽度的关系。考虑了衰减的情况下,在各传播距离、频厚积和激发脉冲周数上,采用脉冲回波法,分析了导波的频散特性。结果表明:对于特定的材料,当频厚积一定时,充粘液管材中传播的导波特性只与内径-壁厚比有关;各高阶模式的截止频厚积随内径-壁厚比的增大而减小。用L(0,1)模式检测内径-壁厚比在4附近的充粘液管时灵敏度较低;而用L(0,2)模式检测内径-壁厚比在9.2附近的管时灵敏度较低。对于不同内径-壁厚比的充粘液管,检测中应当用不同的纵向轴对称导波模式,所用的激发脉冲周数和频厚积也应当不同。     

小波时-频变换的高强钢丝弹性波传播模态分析

超声导波是近年来桥梁拉索无损检测研究的重要方法之一。针对弹性波在高强钢丝介质中传播的多模态频散问题,采用单点时域波形的小波时频变换进行混叠信号的模态识别分离。通过数值求解Pochhammer频率超越方程,计算得到0~1.5 MHz范围内纵向导波模态理论频散曲线;采用有限元模拟半波正弦脉冲激励导波在钢丝中传播过程,由小波时-频变换得到导波模态分布,并进行了不同腐蚀程度钢丝实验对比分析。结果表明,经小波时-频变换得到的第1、2、3阶纵向导波模态与理论值对应吻合,单点时域波形的小波时-频变换结果能够有效识别高强钢丝中的导波模态;钢丝在无腐蚀状态下,一阶纵向导波模态能量占比达57.74%,随腐蚀程度增加,能量更为集中到一阶纵波模态,二阶模态能量逐渐减小。     

充水粘弹性管道的频散曲线计算分析*

针对谱方法分析计算充水粘弹性管道的广义特征值问题,根据Chebyshev多项式及微分矩阵、位移和应力连续条件,将波动方程离散为相应的线性方程。利用MATLAB数值编程计算充水弹性和粘弹性管道对应频率下的轴对称纵向导波频散曲线和衰减曲线。分析表明,波传播在粘弹性管道中不仅具有衰减特性,而且由于水和粘弹性壳体交叉耦合作用,在一定频率范围内产生两种截断模态。     

Circumferential and longitudinal defect detection using T(0,1) mode excited by thickness shear mode piezoelectric elements

Different kinds of defects, such as corrosions, notches and cracks etc, exist in pipes. Mode choice is important since unfortunately not all ultrasonic guided wave modes are suitable for these kinds of defect detection. T(0,1) mode which is non-dispersive is the lowest and fastest torsional mode and most suitable for defect detection in pipes. Two completely different artificial defects including longitudinal and circumferential defects are processed successively in a 4-m-long, 60-mm-OD, 3.5-mm-wall steel pipe. T(0,1) mode at 45 kHz is excited to detect these defects using thickness shear mode piezoelectric elements. Experimental results show that two kinds of defects are detectable using T(0,1) mode. Comparing with longitudinal modes, torsional modes are dominant in pipe inspection for their sensitivities to different kinds of defects.     

The effect of pipe bend on low-frequency longitudinal mode guided wave propagation

The pipe bend significantly changes the propagation characteristics of guided wave,and makes the interpretation of the received signals difficult.Therefore,better understanding of guided wave propagating in bended pipe is essential for the inspection of pipeline comprising bends.First of all,the different features of dispersion curves derived with the semi-analytical finite element method for guided wave in bended pipes are summarized.Secondly,based on the dispersion curves for guided wave in bended pipes,experiments are performed to investigate the mode conversions of L(0,1) mode guided wave traveling through pipe bends.It is found that,except for the mode conversion from L(0,1) to F(1,1),the L(0,1) reflections of bends are also observed in some cases,which are proven to be the mode converted negative L(0,1)mode guided wave,and the negative L(0,1) mode guided wave becomes more obvious with the decrease of excitation frequency and bending radius.The findings of this paper will provide some insight for guided wave behavior in bended pipe,and generalize the application of guided wave inspection in practical pipelines.     

Ultrasonic guided wave parameters for detection of axial cracks in feeder pipes of PHWR nuclear power plants

The dispersion curves for the feeder pipes in PHWR nuclear power plants were determined. The wave modes used for the detection of notches in the feeder pipe were confirmed as F(m,2) and/or L(0,1) by an analysis of short time Fourier transformation (STFT). The axial notches in the straight pipe were not detectable, but an axial notch in a bent pipe was detected with the mode at the frequency of 500 kHz. Initial F(m,2) and/or L(0,1) modes contains a circumferential displacement and might be converted to certain complicated modes in the bent region, which is sensitive to the axial notch. The circumferential guided wave technique was also applied for quantitative evaluation of the axial notches. The waves generated by a rocking motion of the transducer along the circumferential direction were estimated as the circumferential guided waves after a review of the acquired data and the dispersion curves.     

Natural beam focusing of non-axisymmetric guided waves in large-diameter pipes

The propagation of non-axisymmetric guided waves in larger diameter pipes is studied in this paper by treating the guided waves as corresponding Lamb waves in an unwrapped plate. This approximation leads to a simpler method for calculating the phase velocities of hollow cylinder guided waves, which reveals a beam focusing nature of non-axisymmetric guided waves generated by a partial source loading. The acoustic fields in a pipe generated by a partial-loading source includes axisymmetric longitudinal modes as well as non-axisymmetric flexural modes. The circumferential distribution of the total acoustic field, also referred as an angular profile, diverges circumferentially while guided waves propagate with dependence on such factors as mode, frequency, cylinder size, propagation distance, etc. Exact prediction of the angular profile of the total field can only be realized by numerical calculations. In particular cases, however, when the wall thickness is far less than the cylinder diameter and the wavelength is smaller than or comparable to the pipe wall thickness, the acoustic field can be analyzed based on the characteristics of Lamb waves that travel along a periodic unwrapped plate. Based on this assumption, a simplified model is derived to calculate the phase velocities of non-axisymmetric flexural mode guided waves. The model is then applied to discussions on some particular characteristics of guided-wave angular profiles generated by a source loading. Some features of flexural modes, such as cutoff frequency values are predicted with the simpler model. The relationship between the angular profiles and other factors such as frequency, propagation distance, and cylinder size is obtained and presented in simple equations. The angular profile rate of change with respect to propagation distance is investigated. In particular, our simplified model for non-axisymmetric guided waves predicts that the wave beam will converge to its original circumferential shape after the wave propagates for a certain distance. A concept of "natural focal point" is introduced and a simple equation is derived to compute the 1st natural focal distance of non-axisymmetric guided waves. The applicable range of the simplified equation is provided. Industrial pipes meet the requirement of wall thickness being far less than the pipe diameter. The approximate analytical algorithms presented in this paper provides a convenient method enabling quick acoustic field analysis on large-diameter industrial pipes for NDE applications.     

Edge resonance in semi-infinite thick pipe: numerical predictions and measurements

This paper presents theoretical and experimental studies of axisymmetric longitudinal guided wave L(0,2) interaction with the free edge of the pipe. A numerical method based on normal mode superposition is applied to predict the edge resonance by an analysis of dispersion relations of separate modes. In parallel, the finite element analysis and experimental measurements prove the existence of edge resonance in the pipe in case of L(0,2) wave incidence. It is shown that the edge resonance is mainly caused by the first pair of complex modes. Additionally the behavior of edge resonance phenomenon as a function of the curvature of the pipe is studied. The displacement amplitudes measured at the edge demonstrate that the edge resonance is affected by the frequency and thickness to midradius ratio of the pipe, and it is losing its strength in thicker pipes, as the growing difference between the outer and inner radii destroys symmetry. The reflected energy amplitudes show that at the resonance frequencies the incident wave is strongly converted to L(0,1) and L(0,3) modes, depending also on the curvature parameter of the pipe.     

The reflection of the fundamental torsional mode from cracks and notches in pipes

A quantitative study of the reflection of the T(0,1) mode from defects in pipes in the frequency range 10-300 kHz has been carried out, finite element predictions being validated by experiments on selected cases. Both cracklike defects with zero axial extent and notches with varying axial extents have been considered. The results show that the reflection coefficient from axisymmetric cracks increases monotonically with depth at all frequencies and increases with frequency at any given depth. In the frequency range of interest there is no mode conversion at axisymmetric defects. With nonaxisymmetric cracks, the reflection coefficient is a roughly linear function of the circumferential extent of the defect at relatively high frequencies, the reflection coefficient at low circumferential extents falling below the linear prediction at lower frequencies. With nonaxisymmetric defects, mode conversion to the F(1,2) mode is generally seen, and at lower frequencies the F(1,3) mode is also produced. The depth and circumferential extent are the parameters controlling the reflection from cracks; when notches having finite axial extent, rather than cracks, are considered, interference between the reflections from the start and the end of the notch causes a periodic variation of the reflection coefficient as a function of the axial extent of the notch. The results have been explained in terms of the wave-number-defect size product, ka. Low frequency scattering behavior is seen when ka < 0.1, high frequency scattering characteristics being seen when ka > 1.     

Guided wave propagation in single and double layer hollow cylinders embedded in infinite media

Millions of miles of pipes are being used for the transportation, distribution, and local use of petroleum products, gas, water, and chemicals. Most of the pipes are buried in soil, leading to the significance of the study on the subject of guided wave propagation in pipes with soil influence. Previous investigations of ultrasonic guided wave propagation in an elastic hollow cylinder and in an elastic hollow cylinder coated with a viscoelastic material have led to the development of inspection techniques for bare and coated pipes. However, the lack of investigation on guided wave propagation in hollow cylinders embedded in infinite media like soil has hindered the development of pipe inspection methods. Therefore the influence of infinite media on wave propagation is explored in this paper. Dispersion curves and wave structures of both axisymmetric and nonaxisymmetric wave modes are developed. Due to the importance of the convergence of numerical calculations, the requirements of thickness and element number of the finite soil layer between hollow cylinder and infinite element layer are discussed, and an optimal combination is obtained in this paper. Wave structures are used for the mode identification in the non-monotonic region caused by the viscoelastic properties of coating and infinite media.     

Guided wave propagation and mode differentiation in hollow cylinders with viscoelastic coatings

Guided wave propagation theories have been widely explored for about one century. Earlier theories on single-layer elastic hollow cylinders have been very beneficial for practical nondestructive testing on piping and tubing systems. Guided wave flexural (nonaxisymmetric) modes in cylinders can be generated by a partial source loading or any nonaxisymmetric discontinuity. They are especially important for guided wave mode control and defect analysis. Previous investigations on guided wave propagation in multilayered hollow cylindrical structures mostly concentrate on the axisymmetric wave mode characteristics. In this paper, the problem of guided wave propagation in free hollow cylinders with viscoelastic coatings is solved by a semianalytical finite element (SAFE) method. Guided wave dispersion curves and attenuation characteristics for both axisymmetric and flexural modes are presented. Due to the fact that dispersion curve modes obtained from SAFE calculations are difficult to differentiate from each other, a mode sorting method is established to distinguish modes by their orthogonality. Theoretical proof of the orthogonality between guided wave modes in a viscoelastic coated hollow cylinder is provided. Wave structures are also calculated and discussed in view of wave mechanics in multilayered cylindrical structures containing viscoelastic materials.     

圆管结构中周向导波非线性效应的模式展开分析

在二阶微扰近似条件下, 采用导波模式展开分析方法研究了圆管结构中周向导波的非线性效应. 伴随基频周向导波传播所发生的二次谐波, 可视为由一系列二倍频周向导波模式叠加而成. 从动量定理出发, 结合柱坐标系下非线性应力张量及其散度的数学表达式, 针对圆管中某一基频周向导波模式, 推导出相应的二倍频应力张量及二倍频彻体驱动力的数学表达式, 建立了确定二倍频周向导波模式展开系数的控制方程, 得到了伴随基频周向导波传播所发生的二次谐波声场的形式解. 理论分析和数值计算表明, 当构成二次谐波声场的某一二倍频周向导波模式与基频周向导波的相速度匹配时, 该二倍频周向导波模式的位移振幅表现出随传播周向角积累增长的性质; 当两者的相速度失配时, 二倍频周向导波的振幅随传播周向角表现出“拍”效应.     

Interaction of low-frequency axisymmetric ultrasonic guided waves with bends in pipes of arbitrary bend angle and general bend radius

The use of ultrasonic guided waves for the inspection of pipes with elbow and U-type bends has received much attention in recent years, but studies for more general bend angles which may also occur commonly, for example in cross-country pipes, are limited. Here, we address this topic considering a general bend angle φ, a more general mean bend radius R in terms of the wavelength of the mode studied and pipe thickness b. We use 3D Finite Element (FE) simulation to understand the propagation of fundamental axisymmetric L(0, 2) mode across bends of different angles φ. The effect of the ratio of the mean bend radius to the wavelength of the mode studied, on the transmission and reflection of incident wave is also considered. The studies show that as the bend angle is reduced, a progressively larger extent of mode-conversion affects the transmission and velocity characteristics of the L(0, 2) mode. However the overall message on the potential of guided waves for inspection and monitoring of bent pipes remains positive, as bends seem to impact mode transmission only to the extent of 20% even at low bend angles. The conclusions seem to be valid for different typical pipe thicknesses b and bend radii. The modeling approach is validated by experiments and discussed in light of physics of guided waves.