The effect of boundary conditions on guided wave propagation in two-dimensional models of healing bone

Guided wave propagation has recently drawn significant interest in the ultrasonic characterization of bone. In this work, we present a two-dimensional computational study of ultrasound propagation in healing bones aiming at monitoring the fracture healing process. In particular, we address the effect of fluid loading boundary conditions on the characteristics of guided wave propagation, using both time and time-frequency (t-f) signal analysis techniques, for three study cases. In the first case, the bone was assumed immersed in blood which occupied the semi-infinite spaces of the upper and lower surfaces of the plate. In the second case, the bone model was assumed to have the upper surface loaded by a 2mm thick layer of blood and the lower surface loaded by a semi-infinite fluid with properties close to those of bone marrow. The third case, involves a three-layer model in which the upper surface of the plate was again loaded by a layer of blood, whereas the lower surface was loaded by a 2mm layer of a fluid which simulated bone marrow. The callus tissue was modeled as an inhomogeneous material and fracture healing was simulated as a three-stage process. The results clearly indicate that the application of realistic boundary conditions has a significant effect on the dispersion of guided waves when compared to simplified models in which the bone's surfaces are assumed free.     

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.     

无限流体中孔隙介质圆柱周向导波的传播特性

为研究无限大流体约束的孔隙圆柱中周向导波的传播规律,分析孔隙参数对导波传播特性的影响,建立了无限流体中孔隙介质圆柱的理论模型,利用孔隙介质弹性波动理论,建立了周向导波频散方程,通过数值模拟计算得到无限流体中孔隙介质圆柱的频散曲线,探讨了圆柱半径和孔隙参数对导波传播特性的影响,并对导波的衰减特性进行了分析;通过数值计算,得到了周向导波的时域波形,讨论了孔隙参数对波形的影响.结果表明,孔隙介质圆柱半径的改变影响圆柱尺度,孔隙度的改变影响孔隙介质中体声波的波速,都对周向导波频散曲线产生一定的影响,所得到的频散曲线特征及衰减曲线与时域波形吻合.研究结果对开展无限流体中孔隙介质圆柱的超声无损评价提供了一定的理论参考.     

A numerical study on the propagation of Rayleigh and guided waves in cortical bone according to Mindlin's Form II gradient elastic theory

Cortical bone is a multiscale heterogeneous natural material characterized by microstructural effects. Thus guided waves propagating in cortical bone undergo dispersion due to both material microstructure and bone geometry. However, above 0.8 MHz, ultrasound propagates rather as a dispersive surface Rayleigh wave than a dispersive guided wave because at those frequencies, the corresponding wavelengths are smaller than the thickness of cortical bone. Classical elasticity, although it has been largely used for wave propagation modeling in bones, is not able to support dispersion in bulk and Rayleigh waves. This is possible with the use of Mindlin's Form-II gradient elastic theory, which introduces in its equation of motion intrinsic parameters that correlate microstructure with the macrostructure. In this work, the boundary element method in conjunction with the reassigned smoothed pseudo Wigner-Ville transform are employed for the numerical determination of time-frequency diagrams corresponding to the dispersion curves of Rayleigh and guided waves propagating in a cortical bone. A composite material model for the determination of the internal length scale parameters imposed by Mindlin's elastic theory is exploited. The obtained results demonstrate the dispersive nature of Rayleigh wave propagating along the complex structure of bone as well as how microstructure affects guided waves.     

轨道板中超声兰姆波传播的有限元计算和实验*

该文针对我国高速铁路轨道板缺陷的非接触动态检测问题,研究了空气耦合超声兰姆波在轨道板中的传播规律。首先,给出了轨道板中超声兰姆波的相速度和群速度频散曲线,结果表明:随着频厚积的增加,频散现象越明显,并且A0相速度收敛于Rayleigh波的波速。然后,建立轨道板中波传播的有限元模型,计算得到兰姆波传播的群速度为2220 m/s,且二维傅里叶变换系数的较大值沿Rayleigh波的频散曲线分布。最后,在沪杭高铁嘉兴南站进行了现场测试,以8.8°倾斜角向轨道板激励产生超声兰姆波,激发产生的兰姆波模态群速度为2325 m/s,且二维傅里叶变换分析其系数的较大值沿Rayleigh波的频散曲线分布。有限元计算结果和实验结果均与理论计算结果一致。该研究为后续轨道板缺陷的非接触动态检测提供了理论依据和实验方法。     

Wave propagation along transversely periodic structures

The dispersion curves for guided waves have been of constant interest in the last decades, because they constitute the starting point for NDE ultrasonic applications. This paper presents an evolution of the semianalytical finite element method, and gives examples that illustrate new improvements and their importance for studying the propagation of waves along periodic structures of infinite width. Periodic boundary conditions are in fact used to model the infinite periodicity of the geometry in the direction normal to the direction of propagation. This method allows a complete investigation of the dispersion curves and of displacement/stress fields for guided modes in anisotropic and absorbing periodic structures. Among other examples, that of a grooved aluminum plate is theoretically and experimentally investigated, indicating the presence of specific and original guided modes.     

长骨中振动声激发超声导波的方法

为了实现一定频段内任意低频下在长骨中激励导波信号,本文提出一种采用聚焦高频(5 MHz)超声换能器在长骨皮质骨中激发低频(150 kHz)超声导波的振动声方法.首先介绍了板状超声导波理论和双声束共聚焦法与单声束调幅法激发振动声的基本原理;进而采用三维有限元仿真方法分析振动声激发低频超声导波的基本现象,然后结合牛胫骨板离体实验,验证振动声激发低频超声导波的可行性.结果均表明,双声束共焦与单声束振动超声均可在骨板中激发低频超声导波.相关研究方法有助于提高空间域长骨中超声导波测量精度,以及在一定频段内实现任意频率激励等,对发展低频超声导波在体测量长骨皮质骨的新技术具有一定的指导意义.     

Rayleigh wave and detection of low-velocity layers in a stratified half-space

The excitation and propagation of the guided waves in a stratified half-space and a Rayleigh wave exploration method in shallow engineering seismic exploration are studied in this paper. All the modes of the guided waves are calculated by the bisection method in the case where the low velocity layers are contained in a stratified half-space. Cases when the formation shear wave velocity gradually decreases from the top to the bottom layers are also studied. The dispersion curves obtained in actual Rayleigh wave exploration are usually noncontinual zigzag curves, but the dispersion curves given by the elastic theory for given modes of the guided waves are smooth and continual curves. In this paper, the mechanism of zigzag dispersion curves in Rayleigh wave exploration is investigated and analyzed thoroughly. The zigzag dispersion curves can give not only the possible positions of the low-velocity layers but also the other information on the formation structure (fractures, oil, gas, etc.). It is found that the zigzag dispersion curves of the Rayleigh wave are the result of the leap of the modes and the existence of low velocity layers in a stratified half-space. The effects of the compressional wave velocity, shear wave velocity, and density of each layer on zigzag dispersion curves and the relationship of the low velocity layers to zigzag dispersion curves are also investigated in detail. Finally, the exploration depth of the Rayleigh wave is discussed. The exploration depth of the Rayleigh wave is equal to the wavelength multiplied by a coefficient that is variable and usually given by the work experience and the formation properties of the local work area.     

Excitation mechanisms and dispersion characteristics of guided waves in multilayered cylindrical solid media

This paper first reviews a method of simulating the propagation characteristics of guided waves in multilayered coaxial cylindrical elastic solid media. Secondly, this method is used to investigate the properties of the guided waves for the ultrasonic long-range non-destructive evaluation techniques for rockbolts. To do so, the special case of non-leaky guided modes in open waveguides is considered. The method explains how the complex dispersion function is converted into a real function: hence the bisection technique can be employed to search for all the real roots. The model is used to (i) characterize the low dispersion range and anomalous dispersion of normal and Stoneley modes and (ii) analyze the excitation mechanisms of guided waves from axisymmetric and non-axisymmetric acoustic sources. The results are used to select suitable excitation frequency ranges associated with dominant modes with large amplitudes, low dispersion, and distinguishable propagation velocities to reduce signal distortion. The results suggest the lowest order flexural mode, excited by a radial force source, has potential to be used in practice. Also, the highly dispersive Stoneley mode propagating along a cylindrical interface is defined and distinguished from the normal mode using two properties, velocity high-frequency asymptotes and amplitude distributions along the radial direction.     

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.     

谱有限元数值分析各向异性复合板中导波色散特性

提出谱有限元方法研究层状各向异性复合板中导波的色散特性和波结构。基于三维弹性动力学方程,用有限元方法离散波导截面,波传播方向的位移用简谐波表示,得到了导波色散的特征方程。分析了单层和双层复合板中导波沿不同方向传播的色散特性和波结构,讨论了双层复合板中层厚比对相速度的影响。数值研究结果表明:导波的对称模态沿纤维方向传播时在较宽的频率范围内保持弱色散状态。双层复合板中导波基本模态的相速度在低频时受层厚比的影响较明显,随着频率的增加趋向于相速度较低的材料。数值模拟结果为导波用于复合材料定量无损检测和性能评价提供理论依据。      

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.     

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.     

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.     

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.     

Propagation characteristics of guided waves in a rod surrounded by an infinite solid medium

The dispersion and excitation characteristics of the guided waves in a rod surrounded by an infinite solid medium (cladding) are investigated. First, the bisection technique is employed to find all the roots of the dispersion function on the basis of theoretical analysis and to obtain the complex phase and group velocity dispersion curves of the guided modes. Second, according to their different dispersion characteristics, the guided modes are divided into two categories: normal modes and Stoneley modes. And it is concluded that the normal modes merely exist in the “hard cladding” model in which the cladding’s shear velocity is larger than the rod’s; while the Stoneley modes in cylindrical interface are highly dispersive and merely exist in the model whose acoustical parameters satisfied the existence condition of the Stoneley waves. Third, the seldom discussed issue, the excitation mechanisms of the guided waves, excited by three source models: symmetric point source, axial and radial force sources, are simulated respectively. Attention is paid on the dominant mode which has better excitation sensitivity and the suitable excitation frequency range. Moreover, the propagation characteristics of the Stoneley modes, ignored in previous references, are analyzed and compared with those of the normal modes.     

Quantum and geometric effects on the symmetric and anti-symmetric modes of the surface plasma wave

The propagation of the surface quantum plasma waves is investigated in a thin quantum plasma slab. The symmetric and anti-symmetric dispersion modes of the quantum surface wave are obtained by the plasma dielectric function with the kinetic dispersion model for the slab geometry. The quantum mechanical and slab geometric effects on the symmetric and anti-symmetric modes are also discussed.     

长骨裂纹宽度和角度对低阶超声导波的影响

超声导波检测长骨骨折和监测骨折愈合已成为一种极具前景的课题,但高频厚积下导波多模式混叠带来的模式识别和分离问题一直是一个难点。为避免上述问题,在低频下仅激励两个低阶兰姆波SO和AO模式,采用二维时域有限差分(2DFDTD)法定量分析裂纹宽度和骨折角度对SO和AO幅度的影响。结果表明,SO与AO模式的幅度均随裂纹宽度的增加而下降;AO的幅度随骨折角度增大而上升,SO幅度随骨折角度增大先下降后上升,转折点约为45°。SO与AO的幅度比值在不同骨折角度下均能较好地表征裂纹宽度的变化,可为横断型和斜切型长骨骨折状况的超声评价及骨折愈合监测提供依据。      

时频方法分析长骨中的超声导波及皮质骨厚度

超声导波在长骨中传播时,接收信号中含有相互叠加的多个导波模式。本文提出将希尔伯特-黄变换(HHT)用于分析叠加的多模式导波信号,将其分解成许多单个独立的模式,然后对分解出的模式求得其对应频率下的群速度,并与短时傅里叶变换所得的结果进行比较。通过与理论计算结果比较,可得到长骨皮质骨的厚度。研究结果表明,实验得到的厚度与实际厚度一致。说明通过测量导波模式的速度可以评价皮质骨的厚度,也说明HHT方法是一种识别叠加多模式导波信号的有效方法。      

Borehole guided waves in a non-Newtonian (Maxwell) fluid-saturated porous medium

<正>The property of acoustic guided waves generated in a fluid-filled borehole surrounded by a non-Newtonian (Maxwell) fluid-saturated porous formation with a permeable wall is investigated.The influence of non-Newtonian effects on acoustic guided waves such as Stoneley waves,pseudo-Rayleigh waves,flexural waves,and screw waves propagations in a fluid-filled borehole is demonstrated based on the generalized Biot-Tsiklauri model by calculating their velocity dispersion and attenuation coefficients.The corresponding acoustic waveforms illustrate their properties in time domain.The results are also compared with those based on generalized Biot's theory.The results show that the influence of non-Newtonian effect on acoustic guided wave,especially on the attenuation coefficient of guided wave propagation in borehole is noticeable.