Detection of defects in cylindrical structures using a time reverse method and a finite-difference approach

The detection and characterization of defects in structures is an important issue in non-destructive testing. To avoid the scanning of large samples, guided elastic waves, which propagate along the structure, are excited. These waves interact with a defect, which results in a scattered wave field. In an experiment, the displacements of these scattered waves are recorded over time for a fixed axial coordinate at a number of circumferential positions of a circular cylindrical tube. Since in complex structures it is difficult to determine the axial and particularly circumferential position of the defect directly from the time signals, a time reversed numerical simulation is performed. There the measured displacement histories are reversed in time and used as displacement excitations in a simulation of the tested structure. A three-dimensional code in cylindrical coordinates, based on a velocity-stress finite-difference method, is used to simulate the wave propagation. As long as the geometric and material parameters are chosen equivalent to the performed experiment, the scattered waves travel back through the simulated structure and interfere, even if no defect is present in the numerical model. The result is an increase of the amplitudes of the stress and displacement components at the location where the defect was in the tested sample.     

Guided waves in a stratified half-space

The dispersion and excitation mechanisms and the energy distribution of guided waves in a stratified half-space are studied. All possible guided waves excited by a symmetric point source in two or three-layer medium models and their relation to the medium parameters are analyzed in detail. The excitation and propagation characteristics, as well as the energy distribution along the depth direction, of all modes of the surface waves and trapped waves are numerically investigated and analyzed thoroughly not only in the case when the shear wave velocity increases from up to down layers but also when a low-velocity layer is contained in halfspace, especially when the shear wave velocity decreases from up to down layers. It is found that there exist many guided wave modes in the case where the shear wave velocity of each layer increases from up to down layers. However, there is less than one guided wave mode in the case where the shear wave velocity of each layer decreases from up to down layers. The trapped waves exist and propagate along the low-velocity structure in the stratified half-space. It is also found that the characteristic of a mode is related to the source frequency. It is possible that a surface wave at one value of frequency is like a trapped wave at another value of frequency. Finally, the relation of the characteristics of all guided waves (surface waves and trapped waves) to the parameters of media is studied.     

Damage detection in composite panels based on mode-converted Lamb waves sensed using 3D laser scanning vibrometer

This paper introduces damage identification approach based on guided ultrasonic waves and 3D laser Doppler vibrometry. The method is based on the fact that the symmetric and antisymmetric Lamb wave modes differ in amplitude of the in-plane and out-of-plane vibrations. Moreover, the modes differ also in group velocities and normally they are well separated in time. For a given time window both modes can occur simultaneously only close to the wave source or to a defect that leads to mode conversion. By making the comparison between the in-plane and out-of-plane wave vector components the detection of mode conversion is possible, allowing for superior and reliable damage detection. Experimental verification of the proposed damage identification procedure is performed on fuel tank elements of Reusable Launch Vehicles designed for space exploration. Lamb waves are excited using low-profile, surface-bonded piezoceramic transducers and 3D scanning laser Doppler vibrometer is used to characterize the Lamb wave propagation field. The paper presents theoretical background of the proposed damage identification technique as well as experimental arrangements and results.     

Influence of the slow wave on the relation between the angular resonances and the leaky guided modes properties for a poroelastic plate embedded in water

A numerical study of the guided modes in a water-saturated poroelastic plate that obeys the Biot theory is presented. In the first part, we study the leaky guided modes and the angular resonances when the slow wave does not propagate. Two types of guided modes exist. The first ones occur from coupling of the fast longitudinal wave with the shear wave; most of them propagate whatever the frequency is, provided that it is not close to their cut-off frequencies. The leaky guided modes of the second type occur from coupling of the two longitudinal waves and the shear wave. These modes do not propagate (they are highly damped) as long as the slow wave remains diffusive. We show that the characteristics of the angular resonances can be linked to the leaky guided waves of the first type in the same way as for an elastic plate. The guided modes of the second type may not be associated to angular resonances. In the second part, we consider a thinner plate in a higher frequency range so that the slow wave can propagate. Once again its influence is studied both on the leaky guided modes and on the angular resonances.     

Detection of cracks at rivet holes using guided waves

Guided Lamb waves can be used for a fast inspection of large areas, e.g. the detection of cracks at rivet holes in the fuselage of airplanes. When the guided wave hits a discontinuity like a hole, a typical scattered displacement field is obtained. A change of the scattered field indicates the development of a crack. In the experiments, the first anti-symmetric mode A0 of Lamb waves in plates is excited selectively by means of a piezoelectric transducer. The used frequency range is below the cut-off frequencies of higher wave modes. The scattered field around undamaged and damaged holes is measured on a grid around the hole with a heterodyne laser interferometer. Two types of damage are introduced: a notch cut with a very fine saw blade, and a fatigue grown crack. A significant change in the scattered field due to the defect is seen. Good agreement of the experimental results with theoretical calculations is obtained. The wave propagation is studied using Mindlin's theory of plates. The scattered field is calculated analytically and using finite difference methods (FDMs).     

A new multichannel time reversal focusing method for circumferential Lamb waves and its applications for defect detection in thick-walled pipe with large-diameter

This paper proposes a new multichannel time reversal focusing (MTRF) method for circumferential Lamb waves which is based on modified time reversal algorithm and applies this method for detecting different kinds of defects in thick-walled pipe with large-diameter. The principle of time reversal of circumferential Lamb waves in pipe is presented along with the influence from multiple guided wave modes and propagation paths. Experimental study is carried out in a thick-walled and large-diameter pipe with three artificial defects, namely two axial notches on its inner and outer surface respectively, and a corrosion-like defect on its outer surface. By using the proposed MTRF method, the multichannel signals focus at the defects, leading to the amplitude improvement of the defect scattered signal. Besides, another energy focus arises in the direct signal due to the partial compensation of dispersion and multimode of circumferential Lamb waves, alongside the multichannel focusing, during MTRF process. By taking the direct focus as a time base, accurate defect localization is implemented. Secondly, a new phenomenon is exhibited in this paper that defect scattered wave packet appears just before the right boundary of truncation window after time reversal, and to which two feasible explanations are given. Moreover, this phenomenon can be used as the theoretical basis in the determination of defect scattered waves in time reversal response signal. At last, in order to detect defects without prior knowing their exact position, a large-range truncation window is used in the proposed method. As a result, the experimental operation of MTRF method is simplified and defect detection and localization are well accomplished.     

Hidden corrosion detection in aircraft aluminum structures using laser ultrasonics and wavelet transform signal analysis

Preliminary results of hidden corrosion detection in aircraft aluminum structures using a noncontact laser based ultrasonic technique are presented. A short laser pulse focused to a line spot is used as a broadband source of ultrasonic guided waves in an aluminum 2024 sample cut from an aircraft structure and prepared with artificially corroded circular areas on its back surface. The out of plane surface displacements produced by the propagating ultrasonic waves were detected with a heterodyne Mach-Zehnder interferometer. Time-frequency analysis of the signals using a continuous wavelet transform allowed the identification of the generated Lamb modes by comparison with the calculated dispersion curves. The presence of back surface corrosion was detected by noting the loss of the S(1) mode near its cutoff frequency. This method is applicable to fast scanning inspection techniques and it is particularly suited for early corrosion detection.     

Three-Dimensional Scattering of an Incident Plane Shear Horizontal Guided Wave by a Partly through-Thickness Hole in a Plate

We investigate the three-dimensional(3D) scattering problem of an incident plane shear horizontal wave by a partly through-thickness hole in an isotropic plate,in which the Lamb wave modes are also included due to the mode conversions by the scattering obstacle in the 3D problem.An analytical model is presented such that the wave fields are expanded in all of propagating and evanescent SH modes and Lamb modes,and the scattered far-fields of three fundamental guided wave modes are analyzed numerically for different sizes of the holes and frequencies.The numerical results are verified by comparing with those obtained by using the approximate Poisson/Mindlin plate model for small hole radius and low frequency.It is also found that the scattering patterns are different from those of the SO wave incidence.Our work is useful for quantitative evaluation of the plate-like structure by ultrasonic guided waves.     

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.     

TM plane wave scattering and diffraction from a randomly rough half-plane: (part II) An evaluation of the diffraction kernel

In a previous paper (part I), it has been shown that a random wavefield from a randomly rough half-plane for a TM plane wave incidence is written in terms of a Wiener-Hermite expansion with three types of Fourier integrals. This paper studies a concrete representation of the random wavefield by an approximate evaluation of such Fourier integrals, and statistical properties of scattering and diffraction. For a Gaussian roughness spectrum, intensities of the coherent wavefield and the first-order incoherent wavefield are calculated and shown in figures. It is then found that the coherent scattering intensity decreases in the illumination side, but is almost invariant in the shadow side. The incoherent scattering intensity spreads widely in the illumination side, and have ripples at near the grazing angle. Moreover, a major peak at near the antispecular direction, and associated ripples appear in the shadow side. The incoherent scattering intensity increases rapidly at near the random half-plane. These new phenomena for the incoherent scattering are caused by couplings between TM guided waves supported by a slightly random surface and edge diffracted waves excited by a plane wave incidence or by free guided waves on a flat plane without any roughness.     

Guided waves in anisotropic and quasi-isotropic aerospace composites: Three-dimensional simulation and experiment

Three-dimensional (3D) elastic wave simulations can be used to investigate and optimize nondestructive evaluation (NDE) and structural health monitoring (SHM) ultrasonic damage detection techniques for aerospace materials. 3D anisotropic elastodynamic finite integration technique (EFIT) has been implemented for ultrasonic waves in carbon fiber reinforced polymer (CFRP) composite laminates. This paper describes 3D EFIT simulations of guided wave propagation in undamaged and damaged anisotropic and quasi-isotropic composite plates. Comparisons are made between simulations of guided waves in undamaged anisotropic composite plates and both experimental laser Doppler vibrometer (LDV) wavefield data and dispersion curves. Time domain and wavenumber domain comparisons are described. Wave interaction with complex geometry delamination damage is then simulated to investigate how simulation tools incorporating realistic damage geometries can aid in the understanding of wave interaction with CFRP damage. In order to move beyond simplistic assumptions of damage geometry, volumetric delamination data acquired via X-ray microfocus computed tomography is directly incorporated into the simulation. Simulated guided wave interaction with the complex geometry delamination is compared to experimental LDV time domain data and 3D wave interaction with the volumetric damage is discussed.     

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.     

Non-destructive testing of tubes using a time reverse numerical simulation (TRNS) method

A method for the detection of defects in cylindrical structures and the determination of their positions and orientations is presented in this paper. The scattered field, which is generated by the interaction of excited guided waves with a defect, is evaluated with an approach named time reverse numerical simulation method (TRNS). Since the excited waves and the scattered field propagate along the sample, the time-consuming scanning of the whole tube can be eliminated. The scattered displacement field is measured in three dimensions over time with a laser vibrometer at different locations distributed equally around the circumference at a fixed axial coordinate far away from the defect. Instead of analyzing the complicated time signals directly, they are played back in time. If the recorded displacement histories of the scattered field are reversed in time and played back in an identical structure, the waves travel back the same path and interfere to a maximum at their origin. The result is an amplitude increase at the position of the defect where the scattered field was generated. Instead of playing back the recorded time signals in an experiment, this step is replaced by a numerical simulation. Only this enables the visualization and detection of the amplitude increase. As long as the simulation is of high accuracy, the position of the maximum interference corresponds exactly to the location of the defect in the experiment, although no defect is implemented in the simulation.     

Three-dimensional finite element modeling of guided ultrasound wave propagation in intact and healing long bones

The use of guided waves has recently drawn significant interest in the ultrasonic characterization of bone aiming at supplementing the information provided by traditional velocity measurements. This work presents a three-dimensional finite element study of guided wave propagation in intact and healing bones. A model of the fracture callus was constructed and the healing course was simulated as a three-stage process. The dispersion of guided modes generated by a broadband 1-MHz excitation was represented in the time-frequency domain. Wave propagation in the intact bone model was first investigated and comparisons were then made with a simplified geometry using analytical dispersion curves of the tube modes. Then, the effect of callus consolidation on the propagation characteristics was examined. It was shown that the dispersion of guided waves was significantly influenced by the irregularity and anisotropy of the bone. Also, guided waves were sensitive to material and geometrical changes that take place during healing. Conversely, when the first-arriving signal at the receiver corresponded to a nondispersive lateral wave, its propagation velocity was almost unaffected by the elastic symmetry and geometry of the bone and also could not characterize the callus tissue throughout its thickness. In conclusion, guided waves can enhance the capabilities of ultrasonic evaluation.     

Diffuse wavefields in cylindrical coordinates

A diffuse wavefield is normally defined in terms of plane waves--to quote one textbook definition "plane waves are incident from all directions with equal probability and random phase." In some vibro-acoustic problems the response of a two-dimensional component such as a plate is more conveniently expressed in terms of cylindrical waves, and it is not immediately obvious what properties should be assigned to the cylindrical waves to constitute a diffuse field. It is shown here that a diffuse wavefield can be modeled as a summation of statistically independent cylindrical waves, apart from the fact that each outgoing wave of a particular order is fully correlated to an incoming wave of the same order. A simple relationship is derived between the energy flow P in each wave component and the energy density e of the wavefield: P = ec(g)/k, where c(g) is the group velocity and k is the wavenumber. This result is shown to hold true for both bending waves and in-plane waves (longitudinal and shear) in a plate. The work has application to the calculation of coupling loss factors in statistical energy analysis.     

Guided Waves in a Multi-Layered Cylindrical Elastic Solid Medium

We investigate the guided waves in a multi-layered cylindrical elastic solid medium. The dispersion function of guided waves is usually complex and the dispersion curves of all modes are not conveniently obtained. Here we present an effective method to obtain the dispersion curves of all modes. First, the dispersion function of the guided waves is transformed into a real function. The dispersion curves are then calculated for all the modes of the guided waves by the bisection method. The modes with the orders n = 0, 1, and 2 are analysed in two- and three-layer media. The existence condition of Stoneley wave is discussed. The modes of the guided waves are also investigated in a two-layer medium, in which the velocity of shear wave in the outer layer is less than that in the inner layer.     

Evolution of modes in double-clad Raman fiber amplifier

Stimulated Raman scattering in a double cladding optical fiber is studied with a continuous wave laser used as a pump source. Under various launch conditions, pump modes are differently excited. Considering the mode coupling effect among the pump modes, the evolution of the power in the Stokes modes is studied. The results show that the scattered waves （the Stokes waves） in the fiber core with 9%tm diameter and 0.14 NA could propagate predominantly in the fundamental mode of the fiber by carefully adjusting the pump light launching conditions.     

High-frequency guided ultrasonic waves for hidden defect detection in multi-layered aircraft structures

Aerospace structures often contain multi-layered metallic components where hidden defects such as fatigue cracks and localized disbonds can develop, necessitating non-destructive testing. Employing standard wedge transducers, high frequency guided ultrasonic waves that penetrate through the complete thickness were generated in a model structure consisting of two adhesively bonded aluminium plates. Interference occurs between the wave modes during propagation along the structure, resulting in a frequency dependent variation of the energy through the thickness with distance. The wave propagation along the specimen was measured experimentally using a laser interferometer. Good agreement with theoretical predictions and two-dimensional finite element simulations was found. Significant propagation distance with a strong, non-dispersive main wave pulse was achieved. The interaction of the high frequency guided ultrasonic waves with small notches in the aluminium layer facing the sealant and on the bottom surface of the multilayer structure was investigated. Standard pulse-echo measurements were conducted to verify the detection sensitivity and the influence of the stand-off distance predicted from the finite element simulations. The results demonstrated the potential of high frequency guided waves for hidden defect detection at critical and difficult to access locations in aerospace structures from a stand-off distance.     

Extremely asymmetrical scattering of slab modes in periodic Bragg arrays

The steady-state extremely asymmetrical scattering of electromagnetic modes in a slab with a periodically corrugated boundary is analyzed theoretically. A new approach, based on allowance for the diffractional divergence of a scattered wave, is used with the approximation of slowly varying amplitudes and a Fourier analysis. The structure of the incident and scattered waves inside and outside the array is determined. The amplitudes of the scattered waves are found to be much larger than the amplitude of the incident wave. The typical time of relaxation to steady-state scattering is found to depend on the distance from the array boundary through which the incident wave enters the array. Conditions of applicability of the results obtained are also presented.     

分层固体板中导波的激发与频散特性

针对无限大弹性分层固体板,研究了结构中导波的频散和激发特性。首先使用传递矩阵法推导分层板模型中导波的频散方程,然后用二分法求取导波各模式的频散曲线,进而分析结构中导波的频散特性。结果表明:在速度递增或递减的分层板中,基阶模式和高阶模式的高频极限分别等于低速层的瑞利波速和横波波速。对于含低速夹层的分层板,所有模式的高频极限都等于低速层的横波速度。在导波激发特性方面,研究了在具有一定宽度的法向力源作用下的分层板中导波各模式在结构中的法向位移谱。发现在速度递增的分层板结构中基阶模式是主导模式,而对于速度递减和含低速夹层模型,主导模式在不同的频段范围内对应不同的导波模式。