Numerical predictions and experiments for optimizing hidden corrosion detection in aircraft structures using Lamb modes

To increase the sensitivity of Lamb waves to hidden corrosion in aircraft structures, a preliminary step is to understand the phenomena governing this interaction. A hybrid model combining a finite element approach and a modal decomposition method is used to investigate the interaction of Lamb modes with corrosion pits. The finite element mesh is used to describe the region surrounding the corrosion pits while the modal decomposition method permits to determine the waves reflected and transmitted by the damaged area. Simulations make easier the interpretation of some parts of the measured waveform corresponding to superposition of waves diffracted by the corroded area. Numerical results permit to extract significant information from the transmitted waveform and thus to optimize the signal processing for the detection of corrosion at an early stage. Now, we are able to detect corrosion pits down to 80-μm depth distributed randomly on a square centimeter of an aluminum plate. Moreover, thickness variations present on aircraft structures can be discriminated from a slightly corroded area. Finally, using this experimental setup, aircraft structures have been tested.     

Modal decomposition method for modeling the interaction of Lamb waves with cracks

The interaction of the low-order antisymmetric (a0) and symmetric (s0) Lamb waves with vertical cracks in aluminum plates is studied. Two types of slots are considered: (a) internal crack symmetrical with respect to the middle plane of the plate and (b) opening crack. The modal decomposition method is used to predict the reflection and transmission coefficients and also the through-thickness displacement fields on both sides of slots of various heights. The model assumes strip plates and cracks, thus considering two-dimensional plane strain conditions. However, mode conversion (a0 into s0 and vice versa) that occurs for single opening cracks is considered. The energy balance is always calculated from the reflection and transmission coefficients, in order to check the validity of the results. These coefficients together with the through-thickness displacement fields are also compared to those predicted using a finite element code widely used in the past for modeling Lamb mode diffraction problems. Experiments are also made for measuring the reflection and transmission coefficients for incident a0 or s0 lamb modes on opening cracks, and compared to the numerical predictions.     

3D analysis of interaction of Lamb waves with defects in loaded steel plates

The objective of the research presented here is the investigation of the interaction of guided waves with welds, defects and other non-uniformities in steel plates loaded by liquid. The investigation has been performed using numerical simulation for 2D and 3D cases by the finite differences method, finite element method and measurement of 3D distributions of acoustic fields. Propagation of the S(0) mode in a steel plate and its interaction with non-uniformities was investigated. It was shown that using the measured leaky wave signals in the water loading of the steel plate and by application of signal processing, the 3D ultrasonic field structure inside and outside of the plate can be reconstructed. The presence of leaky wave signals over the defect caused by the mode conversion of Lamb waves has been proved using the numerical modelling and experimental investigations. The developed signal and data processing enables to visualise dynamics of ultrasonic fields over the plate, and also to estimate spatial positions of defects inside the steel plates.     

The interaction of the S0 Lamb mode with vertical cracks in an aluminium plate

This article presents the interaction of the first symmetric Lamb mode S0 with vertical cracks in an aluminium plate placed in vacuum. The cracks are symmetrical regarding to the median plane of the plate and their heights are increasing from 0% to 100% of the plate thickness, by steps of 25%. The frequency-thickness product is chosen to be lower than the S1 frequency cut-off. A modal decomposition method is used to solve the diffraction problem. The variation with the height of the crack of the reflection and transmission coefficients of modes propagating in the far field is predicted. The displacement fields at both sides of the cracks are also calculated, so that it is possible to quantify the crack-opening displacement. These results are compared to numerical predictions obtained using a finite element software. Measurements are also conducted and compared to the predictions.     

Piezoelectric transducer parameter selection for exciting a single mode from multiple modes of Lamb waves

Excitation and propagation of Lamb waves by using rectangular and circular piezoelectric transducers surface-bonded to an isotropic plate are investigated in this work. Analytical stain wave solutions are derived for the two transducer shapes, giving the responses of these transducers in Lamb wave fields. The analytical study is supported by a numerical simulation using the finite element method. Symmetric and antisymmetric components in the wave propagation responses are inspected in detail with respect to test parameters such as the transducer geometry, the length and the excitation frequency. By placing only one piezoelectric transducer on the top or the bottom surface of the plate and weakening the strength of one mode while enhancing the strength of the other modes to find the centre frequency, with which the peak wave amplitude ratio between the S0 and A0 modes is maximum, a single mode excitation from the multiple modes of the Lamb waves can be achieved approximately. Experimental data are presented to show the validity of the analyses. The results are used to optimize the Lamb wave detection system.     

Lamb waves in phononic crystal slabs: truncated plane parallels to the axis of periodicity

A theoretical study is presented on the propagation properties of Lamb wave modes in phononic crystal slabs consisting of a row or more of parallel square cylinders placed periodically in the host material. The surfaces of the slabs are parallel to the axis of periodicity. The dispersion curves of Lamb wave modes are calculated based on the supercell method. The finite element method is employed to calculate the band structures and the transmission power spectra, which are in good agreement with the results by the supercell method. We also have found that the dispersion curves of Lamb waves are strongly dependent on the crystal termination, which is the position of the cut plane through the square cylinders. There exist complete or incomplete (truncated) layers of square cylinders with the change of the crystal termination. The influence of the crystal termination on the band gaps of Lamb wave modes is analyzed by numerical simulations. The variation of the crystal termination leads to obvious changes in the dispersion curves of the Lamb waves and the widths of the band gaps.     

Reflection and transmission coefficients for guided waves reflected by defects in viscoelastic material plates

The paper presents a Fourier transform-based signal processing procedure for quantifying the reflection and transmission coefficients and mode conversion of guided waves diffracted by defects in plates made of viscoelastic materials. The case of the S(0) Lamb wave mode incident on a notch in a Perspex plate is considered. The procedure is applied to numerical data produced by a finite element code that simulates the propagation of attenuated guided modes and their diffraction by the notch, including mode conversion. Its validity and precision are checked by the way of the energy balance computation and by comparison with results obtained using an orthogonality relation-based processing method.     

On the reflection of coupled Rayleigh-like waves at surface defects in plates

The reflection of coupled Rayleigh-like waves from surface defects in elastic plates is investigated experimentally and analyzed on the basis of an analytical model and finite difference simulations. The propagation of Rayleigh-like waves in plates is characterized by an energy transfer to the opposite plate side and back over a distance called the beat length. Experimental results clearly show this beating effect and its dependency on the frequency-thickness product, and excellent agreement is obtained with existing analytical predictions. The propagation and scattering are modeled separately for the fundamental A(0) and S(0) Lamb modes that constitute the incident Rayleigh-like wave. The reflection coefficients from surface slots are investigated using finite difference simulations and the reflected Rayleigh-like wave is obtained by superposition. The theoretical model reveals strong dependencies of the reflected field on the ratio between excitation distance and beat length and on the cutoff frequencies of specific higher Lamb modes. Standard pulse-echo measurements allow for the detection of small defects from a remote transducer location. Good agreement is obtained between the predicted and measured amplitude spectra of the reflected Rayleigh-like wave. The developed model allows for the evaluation of defect location and damaged plate side using a combination of time-of-flight and frequency measurements.     

Multiple reflections of Lamb waves at a delamination

Wave propagation in laminated plates with delaminations was calculated using the semi-analytical finite element method. The visualization results and deeper numerical analyses revealed the following phenomena on the fundamental Lamb modes at delamination regions of laminated plates. First, Lamb wave propagates toward the delamination, and then splits into two independent waves at the "Entrance" of the delamination with no significant reflections. These two waves reach the "Exit" of the delamination with the different phases and arrival time. Thus reflected and transmitted waves are excited at the Exit. The repetition of such reflections at the Exit causes the multiple reflections at regular intervals corresponding to the delamination size.     

Modeling of lamb waves generated by integrated transducers in composite plates using a coupled finite element-normal modes expansion method

Thin piezoelectric transducers attached to or embedded within composite structures could be used for in situ structural health monitoring. For plate-shaped structures, the useful ultrasonic vibration modes are Lamb waves. Preliminary testing has already demonstrated the suitability and practical feasibility of such integrated transducers, but better control of the generation of Lamb modes seems to be necessary. Therefore, an original modeling approach has been developed, which can be used to design and optimize these "sensitive materials." This modeling technique allows the determination of the amplitude of each Lamb mode excited in a composite plate with surface-bonded or bulk-embedded piezoelectric elements. The method consists of a coupling of the finite element method (FEM) and the normal modes expansion method. The limited finite element mesh of the transducer and its vicinity enables the computation of the mechanical field created by the transducer, which is then introduced as a forcing function into the normal modes equations. The adequacy and accuracy of this modeling method have been numerically and experimentally verified.     

Remote characterization of defects in plates with viscoelastic coatings using guided waves

The remote inspection for defects in large metallic elements such as pipes, tubes and plates is a field where guided waves are being applied with success. There are a number of situations where a surface coating is added for corrosion protection or insulation purposes. Since the coating materials are usually viscoelastic, the guided wave inspection range may be severely reduced unless a proper mode and an adequate frequency range is selected. Previous authors found the existence of low-attenuated modes at certain frequency ranges, which were used to detect and locate defects at reasonably large distances. This paper studies the potential of guided waves for not only locating but also sizing defects in plates with viscoelastic coatings. A hybrid finite element-boundary element method which explicitly includes the attenuating characteristics of the coating is used to determine Lamb and SH mode conversion factors at corrosion defects under the coating. Through parametric studies and analysis of the numerical results, some trends and features are highlighted that could be used for sizing purposes.     

The use of an orthogonality relation for reducing the size of finite element models for 3D guided waves scattering problems

The scattering of guided waves by complex shaped defects in three-dimensional (3D) waveguides is considered. For such problems, analytical solutions do not exist, and modal decomposition techniques based on the establishment of the displacement and stress fields in the vicinity of the scatterer are quite heavy and complicated to perform. On the other hand, finite elements (FE)-based methods constitute a powerful way to obtain solutions, but they are known to be very memory consuming. This paper proposes a post-processing technique, based on a 3D orthogonality relation, to decompose a complex acoustic field produced by a scatterer and predicted by a 3D FE model, into plane waves, the amplitudes of which are quantified in the far field. This technique allows important reductions in the size of the FE models to be made. Two applications are presented to demonstrate the potential of this method. The first one concerns the scattering of the S₀ Lamb wave incident on a flat bottom circular hole. In this example, the amplitude of each mode is calculated via the orthogonality relation-based method, and compared to that obtained by simply monitoring the displacements at appropriate through-thickness positions. In the second application, the incident S₀ Lamb mode is converted into five modes scattered by a defect of complex geometry.     

Laser-ultrasonic generation of Lamb waves in the reaction force range

The Laser-ultrasonic generation of Lamb waves in an elastic plate is investigated theoretically and experimentally for a laser source whose intensity is high enough to create reaction forces (normal tractions) on the illuminated surface of the specimen. The analytical solutions for transient waves are derived using the integral transform method first by considering an arbitrary source shape and time excitation function, and then specifically for circular and line source shapes. The simulation study allows us not only to predict the behavior of individual wave modes but also to construct the overall responses; thus it helps us better understand the wave excitation mechanisms. The dispersive and multi-modal nature of laser-generated Lamb waves is presented by showing the spatiotemporal Fourier transform of displacements obtained by the simulation study. The transform, displayed in the frequency-wave number domain, enunciates the characteristics of the propagating individual Lamb wave modes. The simulation results are then compared with the 2-D Fourier transform of a set of experimental data obtained by scanning an aluminum plate specimen.     

Lamb wave topological imaging for blind hole defects in isotropic plates

The research aims at validating the ability of topological imaging to blind holes in isotropic plates using Lamb waves. Due to the defect is not symmetric around the midplane of the plate, the effect of Lamb mode conversion will have to be taken into account.The imaging method is based on two computations of ultrasonic fields, one forward and one adjoint, performed for the defect-free reference medium. The excited signal and scattered Lamb waves caused by the blind hole, are used as emitting sources to compute the forward problem and the adjoint problem, respectively. With the help of the finite element simulations,the natural refocusing process of the multimode Lamb waves at the defect location is visually demonstrated by the transient acoustic field snapshots at the different moments to strengthen the physical mechanism of the topological imaging method. The numerical results demonstrate that topological imaging has relatively stronger applicability to the blind hole in contrast to classical Delay And Sum(DAS) method and Time Reversal(TR) method. The topological imaging could handle complex Lamb wave signals containing mode conversions without the imaging quality being affected. The proposed imaging method presents a certain developing potential for detecting and imaging asymmetric defects in plate-like configurations using Lamb waves.     

基于模式分离的兰姆波逆时偏移成像

应用频率域逆时偏移方法实现各向同性和各向异性板中缺陷的兰姆波成像.由于缺陷引起的多模态散射信号会在重建图像中形成伪像,根据基本导波模式振动对称性的差别进行了模式分离预处理.基于多元阵列超声技术,开展了铝板和复合板内缺陷频率域逆时偏移超声成像方法的数值仿真研究.首先,建立有限元模型,采用环形传感器数值采集由缺陷引起的兰姆波散射信号,然后,将采集到的多模式散射信号进行模式分离处理,再将模式分离后的兰姆波散射信号经时间反转后并在相应的接收器处重新激励,在频域中运用格林函数反向传播兰姆波散射信号,获取监测区域的声场信息,与正向传播声场进行互相关,重建缺陷图像.首先对铝板中单缺陷以及复合材料板中相邻的两个相同缺陷进行数值仿真,对比有无模式分离处理的缺陷逆时偏移成像效果,体现出模式分离的重要性.在此基础上,采用逆时偏移方法对复合板材内位置邻近、深度不同的双缺陷进行识别.数值结果表明,模式分离预处理后的缺陷重建图像能够有效去除多模式干扰产生的伪像.文中提出的成像方法对各向同性板和各向异性板内缺陷的检测和成像具有很好的发展潜力,可以准确地探测多个缺陷的形状、尺寸和深度.     

Extraction of dispersion curves for waves propagating in free complex waveguides by standard finite element codes

This paper presents a fast and reliable method, for obtaining all the range of dispersion curves for wave propagation usually used in practice, by numerical simulation only, via common commercial finite element codes. Essentially, the method is based on a simple and robust approach, consisting in a few series of modal analyses for a representative part of the inspected structure. In this way, for different wave lengths, one can find the mode shapes and corresponding natural frequencies by solving some real, symmetric and well numerically conditioned eigenvalue problems. The method allows the extraction of propagating modes only and, in spite of not producing continuous dispersion curves, it is not susceptible to aliasing effects, as some similar methods are. Additionally, complete graphical representations of guided waves are possible with some minor calculus effort.     

各向同性板中盲孔缺陷的兰姆波拓扑成像

将拓扑成像用于各向同性板中盲孔缺陷的兰姆波检测。由于缺陷关于板中面是非对称的,兰姆波模式转换的影响不能忽略。成像方法是基于无缺陷参考媒质中两个超声场(直接声场和伴随声场)的计算。以传感器激励信号为声源,计算直接声场;以盲孔引起的兰姆波散射信号为声源,计算伴随声场。拓扑成像通过直接声场与伴随声场的相乘,可消除无缺陷处多模式混叠的干扰,提高成像分辨率。建立有限元模型,通过不同时刻的瞬态声场图可视化地显示多模式兰姆波在缺陷处的聚焦过程,揭示拓扑成像方法的物理机理。数值结果表明,相对于经典的延迟求和以及时间反转方法,拓扑成像对盲孔缺陷有更强的适应能力。即使对于复杂的包含模式转换的兰姆波信号,拓扑图像的质量依然没有受到影响。文中提出的成像方法对板类结构中非对称缺陷的兰姆波检测具有一定的应用潜力。      

Symmetric and Anti-Symmetric Lamb Waves in a Two-Dimensional Phononic Crystal Plate

It is weft known that Lamb waves in a plate with a mirror plane can be separated into two uncoupled sets： symmetric and anti-symmetric modes. Based on this property, we present a revised plane wave expansion method （PWE） to calculate the band structure of a phononie crystal （PC） plate with a mirror plane. The developed PWE method can be used to calculate the band structure of symmetric and anti-symmetric modes separately, by which the depending relationship between the partial acoustic band gap （PABG）, which belongs to the symmetric and anti-symmetric modes alternatively, and the position of the scatterers can be determined. As an example of its application, the band structure of the Lamb modes in a two-dimensional PC plate with two layers of void circular inclusions is investigated. The results show that the band structure for the symmetric and anti-symmetric modes can be changed by the position of the scatterers drastically, and larger PABGs will be opened when the scatterers are inserted into the area of the plate, where the elastic potential energy is concentrated.     

板状结构Lamb波频散补偿与聚焦接收成像检测

对板状结构中Lamb波的频散和多模现象进行了理论分析,采用有限元仿真方法获取了Lamb波与缺陷作用后的回波信号,对信号进行聚焦接收处理,即在一定距离范围内利用导波频散特性对各个模式信号进行频散补偿,提取补偿后信号的幅值,结果表明当所得信号幅值最大时,所对应的补偿距离等于缺陷与换能器之间的实际距离.利用这一结论,提出了以...     

Simple expressions of the reflection and transmission coefficients of fundamental Lamb waves by a rectangular notch

The scattering of Lamb waves by a two-dimensional rectangular notch is investigated for rapid inspection of defects in a structure. To derive the reflection and transmission coefficients of the scattered waves in a simple way, the scattering caused by the notch is analyzed through the composition of individual scattering processes. Linear equations corresponding to the reflection and transmission coefficients are constructed along with scattering graphs. For an illustration of the efficacy of the presented method, the scattering of fundamental symmetric and anti-symmetric modes are inspected according to the depth and width of a notch in a plate. Validity of these expressions is demonstrated by the comparison of the theoretical analysis results with those from the finite element analysis.