Guided wave dispersion curves for a bar with an arbitrary cross-section,a rod and rail example

Theoretical and experimental issues of acquiring dispersion curves for bars of arbitrary cross-section are discussed. Since a guided wave can propagate over long distances in a structure, guided waves have great potential for being applied to the rapid non-destructive evaluation of large structures such as rails in the railroad industry. Such fundamental data as phase velocity, group velocity, and wave structure for each guided wave mode is presented for structures with complicated cross-sectional geometries as rail. Phase velocity and group velocity dispersion curves are obtained for bars with an arbitrary cross-section using a semi-analytical finite element method. Since a large number of propagating modes with close phase velocities exist, dispersion curves consisting of only dominant modes are obtained by calculating the displacement at a received point for each mode. These theoretical dispersion curves agree in characteristic parts with the experimental dispersion curves obtained by a two-dimensional Fourier transform technique.     

Wave structure analysis of guided waves in a bar with an arbitrary cross-section

Both dispersion curves and wave structures, which are displacement distributions on a bar cross-section, are essential for guided wave NDEs. Theoretical dispersion curves and wave structures for a bar with an arbitrary cross-section are derived in this paper using a special modeling technique called the semi-analytical finite element method (SAFEM). The guidelines for guided wave NDEs of bar-like structures are also shown based on wave structure and modal analysis. First, the relationship between the dispersion curves and their corresponding wave structures were obtained for a square rod. Modes with longitudinal vibration have higher group velocities and torsional modes have constant phase and group velocities. Next, the relationship between the dispersion curves and wave structures for a rail are detailed. The rail is used to represent a bar with a complex cross-section. Similar to the square rod results, the rail's longitudinal modes have higher group velocities. However, the rail contains modes with local vibration. Finally, single mode detection and excitation techniques are introduced. A single mode can be obtained by detecting and exciting with a weighted function that corresponds to a specific mode's wave structure.     

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.     

Wideband dispersion removal and mode separation of Lamb waves based on two-component laser interferometer measurement

Ultrasonic Lamb waves are considered as a sensitive and effective tool for nondestructive testing and evaluation of plate-like or pipe-like structures. The nature of multimode and dispersion causes the wave packets to spread, and the modes overlap in both time and frequency domains as they propagate through the structures. By using a two-component laser interferometer technique, in combination with a priori knowledge of the dispersion characteristics and wave structure information of Lamb wave modes, a two-component signal processing technique is presented for implementing dispersion removal and mode separation simultaneously for two modes mixture signals of Lamb waves. The proposed algorithm is first processed and verified using synthetic Lamb wave signals. Then, the two-component displacements test experiment is conducted using different aluminum plate samples. Moreover, we confirm the effectiveness and robustness of this method.     

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.     

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.     

管道结构中的类兰姆波层析成像

理论和实验上研究了层析成像在管道结构导波检测中的适用性。在管的两个圆周带上分别布置发射和接收换能器,形成多点发射-多点接收的几何结构。将管中绕管道沿最短路径的螺旋方向到达接收换能器的模式作为感兴趣的模式,它们看起来像在弯曲的板中传播,称之为类兰姆波模式。给出了类兰姆波各射线路径走时的理论计算方法,并采用联合迭代重建技术(SIRT)对走时数据进行慢度图像重建。数值模拟和实验研究了管中横穿孔缺陷层析成像的输出并对健康管和缺陷管中类兰姆波的走时数据进行比较,发现穿过缺陷区域的类兰姆波因为绕过缺陷传播而呈现出更长的走时。文中也给出了两缺陷管道类兰姆波层析成像的数值模拟和实验结果,相比单缺陷的情况,射线较少穿过区域缺陷的重建图像质量下降。研究结果为采用类兰姆波模式进行管道缺陷无损评价提供了理论依据。      

Propagation properties of elastic waves in semi-infinite phononic crystals and related waveguides

The transmission and reflection coefficients of two-dimensional semi-infinite solid-solid phononic crystal systems and fluid-fluid phononic waveguide structures have been investigated. The numerical results show that the transmission spectra for longitudinally and transversally polarized incident waves are different, and the spectra of the transmission and reflection coefficients of the semi-infinite system agree well with the band structure. The numerical results show that when a guided wave incident, localized modes are excited, and different polarities have different coupling efficiencies with the incident guided wave. At the same time, far from the cutoff frequency, the guided wave couples out of semi-infinite waveguide highly efficiently.     

Formation of guided spin-wave bullets in ferrimagnetic film stripes

The formation of quasi-2D nonlinear spin-wave eigenmodes in longitudinally magnetized stripes of a ferrimagnetic film, the so-called guided spin-wave bullets, was experimentally observed by using time- and space-resolved Brillouin light scattering spectroscopy and confirmed by numerical simulation. They represent stable spin-wave packets propagating along a waveguide structure, for which both transversal instability and interaction with the side edges of the waveguide are important. The experiments and the numerical simulation of the evolution of the spin-wave excitations show that the shape of the formed packets and their behavior are strongly influenced by the confinement conditions. The discovery of these modes demonstrates the existence of quasistable nonlinear solutions in the transition regime between one-dimensional and two-dimensional wave packet propagation.     

三层介质超声谐振模式随材料和界面粘接性能变化的演变规律

在用超声波谐振对粘接材料的粘接强度进行无损评估时, 不同模式对粘接强度的敏感程度受到众多因素和参数的影响, 对检测结果的可靠性至关重要. 基于多层介质中声传播和界面弱粘接边界条件的理论模型, 将一个上下非对称的金属-粘接剂-金属三层结构的平面波反射系数函数中的谐振模式看作是上下铝金属层各自的Lamb波频散模式通过夹心粘接剂层相互耦合后叠加组成. 改变影响结构粘接强度的因素, 即粘接剂的性能参数(声阻抗、密度、厚度)和界面切向劲度系数k_t来分析三层结构谐振模式耦合方式的变化,得出结论: 粘接结构粘接性能的变化基本上不改变与被粘铝层相关的固有部分的Lamb波模式, 而它们的耦合模式则在谐振频率上产生平移并会与固有模式进行交换和替代; 不同参数的变化引起的模式演变有各自的规律, 大多可彼此区分.     

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.     

Efficient omnidirectional couplers achieved by anisotropic photonic crystal waveguides

Photonic crystals (PCs) have many potential applications because of their ability to control light-wave propagation. We have investigated omnidirectional couplers in two-dimensional anisotropic PC structures. The anisotropic PC coupler composed of an anisotropic-dielectric cylinder in air is studied by solving Maxwell's equations using the plane wave expansion method and finite-difference time-domain method. The photonic band structures are found to exhibit absolute bandgaps for the square lattices. Numerical simulations show that the incident light-waves at both transverse electric and transverse magnetic modes have efficient coupling in anisotropic PC coupler with square lattices. The guided modes and coupling length are analyzed by considering various line defect anisotropic PC waveguides and interaction regions of couplers. Such a mechanism of omnidirectional coupling should open up a new application for designing components in photonic integrated circuits.     

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.     

An Efficient Algorithm for Truncating Spatial Domain in Modeling Light Scattering by Finite-Difference Technique

The finite-difference time domain technique is one of the most robust and accurate numerical methods for the solution of light scattering by small particles with arbitrary composition and geometry. In practice, this method requires that the spatial domain for the computation of near-field be truncated. An absorbing boundary condition must be imposed in conjunction with this truncation. The performance of this boundary condition is essential to the stability of numerical computations and the reliability of results. In the present study, a new boundary condition, referred to as the mixed T algorithm, has been developed, which is a generalization of the transmitting boundary condition originally developed by Liao and co-workers. The present algorithm does not require spatial interpolation for wave values at interior grid points. In addition, it produces two minima of spurious reflections at small and large incident angles, allowing efficient absorption of the scattered waves at the boundary for large incident angles. When the third-order mixed T algorithm is used, the reflection coefficient of the boundary is less than 1% for incident angles from 0° to about 70°. We find that the numerical instability associated with the transmitting boundary condition is caused by the location-dependent amplitude of outgoing waves in the vicinity of the boundary. For this reason, the mixed T algorithm is stabilized by consistently introducing diffusive coefficients into the boundary equation. When the stabilized algorithm is applied, the near-field within the truncated domain can be computed by using single-precision arithmetic without overflows for more than 10⁵steps in the time-marching iteration. Finally, the new absorbing boundary condition is validated by carrying out numerical experiments involving the propagation of a TM wave excited by a sinusoidal point source, simultaneous simulation of the wave propagation in small and large domains, and the scattering of a TM wave by an infinite circular cylinder.     

Comment on "Three-dimensional finite element modeling of guided ultrasound wave propagation in intact and healing long bones," [J. Acoust. Soc. Am. 121(6), 3907-3921 (2007)

Protopappas et al. performed finite element (FE) studies on the propagation of guided ultrasound waves in intact and healing long bones, and found that the dispersion of guided modes was significantly influenced by the irregularity and anisotropy of the bone. A time-frequency (t-f) method was applied to the obtained signals and several wave modes were identified. However, this technique was unable to quantify their observations and provide monitoring capabilities. One possible reason of this shortcoming may come from the inherent disadvantage of the t-f method. The objective of this comment is to demonstrate that it is necessary to combine other techniques with FE simulations for the extraction of significant quantitative ultrasonic features. Individual guided modes in an isotropic pipe have been theoretically examined using the normal mode expansion (NME) method, and many modes that are missed by the t-f analysis have been identified. It is concluded that in order to extract quantitative ultrasonic features, FE simulations should be supplemented by other techniques such as the NME.     

Wave Simulation in Frozen Porous Media

We propose a numerical algorithm for simulation of wave propagation in frozen porous media, where the pore space is filled with ice and water. The model, based on a Biot-type three-phase theory, predicts three compressional waves and two shear waves and models the attenuation level observed in rocks. Attenuation is modeled with exponential relaxation functions which allow a differential formulation based on memory variables. The wavefield is obtained using a grid method based on the Fourier differential operator and a Runge–Kutta time-integration algorithm. Since the presence of slow quasistatic modes makes the differential equations stiff, a time-splitting integration algorithm is used to solve the stiff part analytically. The modeling is second-order accurate in the time discretization and has spectral accuracy in the calculation of the spatial derivatives.     

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.     

High order shear horizontal modes for minimum remnant thickness

Thickness mapping in aging structures suffering from corrosion is challenging especially when the structure is only partially accessible. In plates the high order shear horizontal guided wave modes all have a cutoff frequency thickness product below which they cannot propagate. This property is potentially attractive to estimate the minimum remnant thickness between two transducers. When using a source and a sensor array it is possible to control the number of modes being excited and the size of the region interrogated by the technique. Finite element simulations were used to show that by exciting multiple guided wave modes simultaneously and identifying which modes are received by a sensor array it is possible to estimate the minimum remaining thickness along the propagation path. Initial experimental results showed excellent agreement with the finite element simulations when the plate is uniform and with a thickness reduction between the source and the sensor arrays the minimum remnant thickness was underestimated by approximately 20%.     

Dynamics of Spin I=3/2 under Spin-Locking Conditions in an Ordered Environment

We have derived approximate analytic solutions to the master equation describing the evolution of the spin I=3/2 density operator in the presence of a radio-frequency (RF) field and both static and fluctuating quadrupolar interactions. Spectra resulting from Fourier transformation of the evolutions of the on-resonance spin-locked magnetization into the various coherences display two satellite pairs and, in some cases, a central line. The central line is generally trimodal, consisting of a narrow component related to a slowly relaxing mode and two broad components pertaining to two faster relaxing modes. The rates of the fast modes are sensitive to slow molecular motion. Neither the amplitude nor the width of the narrow component is affected by the magnitude of the static coupling, whereas the corresponding features of the broad components depend in a rather complicated manner on the spin-lock field strength and static quadrupolar interaction. Under certain experimental conditions, the dependencies of the amplitudes on the dynamics are seen to vanish and the relaxation rates reduce to relatively simple expressions. One of the promising emerging features is the fact that the evolutions into the selectively detected quadrupolar spin polarization order and the rank-two double-quantum coherence do not exhibit a slowly relaxing mode and are particularly sensitive to slow molecular motion. Furthermore, these coherences can only be excited in the presence of a static coupling and this makes it possible to discern nuclei in anisotropic from those in isotropic environment. The feasibility of the spin-lock pulse sequences with limited RF power and a nonvanishing average electric field gradient has been demonstrated through experiments on sodium in a dense lyotropic DNA liquid crystal.