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

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

The group velocity variation of Lamb wave in fiber reinforced composite plate

Experimentally measured Lamb wave group velocities in composite materials with anisotropic characteristics are not the same as the theoretical group velocities which is calculated with the Lamb wave dispersion equation. This discrepancy arises from the fact that the angle between the group velocity direction and the phase velocity direction in anisotropic materials exists. Wave propagation in a composite material with anisotropic characteristics should be considered with respect to magnitude correction in addition to direction correction. In this study, S0 mode phase velocity dispersion curves are depicted with the variation of degree with respect to the fiber direction using a Lamb wave dispersion relation in the unidirectional, bidirectional, and quasi-isotropic composite plates. Slowness surface is sketched by the reciprocal value of the phase velocity curves. The magnitude and direction of the group velocity could be calculated from the slowness surface. The recalculated group velocities with consideration of the magnitude and direction from the slowness surface are compared with experimentally measured group velocities. The proposed method shows good agreements with theoretical and experimental results.     

Lamb wave extraction of dispersion curves in micro/nano-plates using couple stress theories

In this paper, Lamb wave propagation in a homogeneous and isotropic non-classical micro/nano-plates is investigated. To consider the effect of material microstructure on the wave propagation, three size-dependent models namely indeterminate-, modified- and consistent couple stress theories are used to extract the dispersion equations. In the mentioned theories, a parameter called ‘characteristic length’ is used to consider the size of material microstructure in the governing equations. To generalize the parametric studies and examine the effect of thickness, propagation wavelength, and characteristic length on the behavior of miniature plate structures, the governing equations are nondimensionalized by defining appropriate dimensionless parameters. Then the dispersion curves for phase and group velocities are plotted in terms of a wide frequency-thickness range to study the lamb waves propagation considering microstructure effects in very high frequencies. According to the illustrated results, it was observed that the couple stress theories in the Cosserat type material predict more rigidity than the classical theory; so that in a plate with constant thickness, by increasing the thickness to characteristic length ratio, the results approach to the classical theory, and by reducing this ratio, wave propagation speed in the plate is significantly increased. In addition, it is demonstrated that for high-frequency Lamb waves, it converges to dispersive Rayleigh wave velocity.     

Wave propagation in an anisotropic metamaterial with single-sheeted hyperboloid dispersion relation

We investigate the wave propagation in an anisotropic metamaterial with single-sheeted hyperboloid dispersion relation. Based on boundary conditions and dispersion relations, we find that the opposite amphoteric refraction, such that E (or H)-polarized waves are positively refracted whereas H (or E)-polarized waves are negatively refracted, can occur at the interface associated with the anisotropic metamaterial. Under a certain condition, both E- and H-polarized waves can exhibit the same single-sheeted hyperboloid or straight line dispersion relation, while the two polarized waves exhibit different propagation characteristics. We expect some potential device applications can be derived based on the unique amphoteric refraction properties. PACS 78.20.Ci; 41.20.Jb; 42.25.Gy     

An alternative method for plotting dispersion curves

Solving the frequency equation and plotting the dispersion curves in problems of wave propagation in cylinders and plates, particularly when the material is anisotropic, are complicated tasks. The traditional numerical methods are usually based on determination of the zeros of the frequency equation by using an iterative find-root algorithm. In this paper, an alternative method is proposed which extracts the solution of the frequency equation in the form of dispersion curves from the three-dimensional illustration of the frequency equation. For this purpose, a three-dimensional representation of the real roots of the frequency equation is first plotted. The dispersion curves, which are the numerical solutions of the frequency equation, are then obtained by a suitable cut in the velocity-frequency plane. The advantages of this method include simplicity, high speed, low possibility of numerical error, and presentation of the results in a graphical form that promotes ease of interpretation. This method is not directly applicable to problems which incorporate high damping or leaky waves. However, if the damping is not very high, it could be a good estimate of the true dispersion curves.     

两侧有固体层负载时板中Lamb波的传播

本文研究了薄板二面有固体导负载时板中Ｌａｍｂ波的传播，从弹性波理论出发并结合应的边界条件，导出板中Ｌａｍｂ波的色散方程，数值计算表示，不管作为自由状态时板中Ｌａｍｂ波相速（板厚取定时）是大于或小于外层固体的声表面波波速，板中对称及反对称模式的Ｌａｍｂ波相速都随着外层固体层厚度增加而变化并且渐近于外层固体的声表面波波速，数值计算又表明，对很薄的板，板中对称及反对称模式的相速均随负载板的厚度呈线性变化     

The propagation of coupled Lamb waves in multilayered arbitrary anisotropic composite laminates

Based on linear three-dimensional elasticity theory, the wave equations of coupled Lamb waves in multilayered arbitrary anisotropic composite laminates are derived using a Legendre orthogonal polynomial approach. The elastodynamic solution for the propagation of coupled Lamb waves in composite plates is also presented to determine the characteristics of coupled Lamb waves. To verify the applicability and validity of the method, two cases of bi-layered plates formed with isotropic components and anisotropic components, respectively, are primarily manipulated for comparison with earlier known results. Next, the dispersion curves, displacements and stress distributions of Lamb waves in multilayered anisotropic laminates are calculated. The effects of coupling and fiber orientation on the characteristics of the Lamb waves are illustrated. The potential usefulness of the fundamental modes of the coupled Lamb waves is discussed in detail.     

The propagation of Lamb waves in a plate bordered with layers of a liquid.

The influence of liquid layers on the propagation of Lamb waves in a plate of finite thickness is studied theoretically. The dispersion equations of Lamb waves in a plate bordered with layers of liquids are derived. Numerical solutions of the equations show that the phase velocity of Lamb waves changes with the thickness of the liquid layers. For the lowest antisymmetrical mode of very thin plates, the numerical results calculated from the dispersion equations are compared with those derived from the bending wave acoustic impedance approach. The limitation of the latter is discussed. Applications of Lamb waves pertinent to biosensing are also presented.     

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

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

Hybridization of acoustic waves in piezoelectric plates

The interaction (hybridization) of different types of acoustic waves of zero and higher orders propagating in lithium niobate piezoelectric plates is theoretically investigated. Different crystallographic orientations of the plates and different directions of wave propagation in them are considered. It is shown that, for an electrically free plate with the propagation direction along any of the crystallographic axes, the dispersion curves have intersection points and hybridization is absent. However, when the propagation direction slightly changes or when one of the plate surfaces is short-circuited, the dispersion curves separate and the waves become coupled. A quantitative coefficient characterizing the degree of wave hybridization with allowance for both mechanical and electric coupling is introduced. It is shown that the dependence of this coefficient on the product of the plate thickness by the wave frequency determines the extent of separation of the dispersion curves of interacting waves. The phenomenon under study is of both fundamental and practical interest, for example, in connection with the problem of an efficient excitation of nonpiezoactive acoustic waves in piezoelectric plates.     

Analysis of the dispersion spectrum of fluid-loaded anisotropic plates: leaky-wave branches

The supersonic complex-velocity versus real-frequency dispersion spectrum of the leaky waves in fluid-loaded anisotropic plates is discussed. Utilizing the sextic plate formalism provides approximate solutions for leaky-wave velocity in a form that reveals their basic features, such as the unique correspondence of the signs of its imaginary part and of the free-plate group velocity, the relation between the leakage and the rate of frequency dispersion, and the principal trends at low, high and near-cutoff frequencies in arbitrary anisotropic plates. A particular thrust of the study is the derivation of closed-form asymptotics for the fundamental leaky-wave velocity branch(es) at low frequency and for the continuum of leaky-wave branches near the fluid-coupled and fluid-uncoupled thickness resonances. Conditions for the asymptotics accuracy are analysed, and a comparison between an analytical approximation and exact numerical curves is presented for various cases.     

Development of an advanced multimode automatic ultrasonic texture measurement system for laboratory and production line application

We present work on the development of an ultrasonic texture measurement system for sheet metals using non-contact transducers, suitable for use both in the laboratory and on the production line. Variation of the velocity of the zero-order symmetric (S0) Lamb wave is used to determine the crystallographic texture of polycrystalline metal sheets ranging in thickness from 0.1 to 3 mm. This system features improvements on previous state-of-the-art ultrasonic technology in that it probes velocity over a continuous range of angles using only two electromagnetic acoustic transducers (EMATs). This is demonstrated to offer a significant improvement in accuracy and allows the detection and investigation of asymmetric anisotropies in the sheets. Another advantage of the system is its potential for combining several different measurements using a single pair of transducers. The capability is demonstrated for through-thickness shear wave measurements as well as the zero-order symmetric Lamb wave measurements which are the primary means of determining the texture. The change between generating Lamb and through-thickness bulk waves can be made entirely by changing the electrical circuit connected to the EMATs without modifying the transducer assembly in any way. Measurement of all of the above waves can provide information on the sheet thickness and other physical properties of the sheet in addition to texture. Certain texture parameters can be calculated from both Lamb and shear wave velocities, allowing self-calibration of the system.     

Time-frequency analysis of the dispersion of Lamb modes.

Accurate knowledge of the velocity dispersion of Lamb modes is important for ultrasonic nondestructive evaluation methods used in detecting and locating flaws in thin plates and in determining their elastic stiffness coefficients. Lamb mode dispersion is also important in the acoustic emission technique for accurately triangulating the location of emissions in thin plates. In this research, the ability to characterize Lamb mode dispersion through a time-frequency analysis (the pseudo-Wigner-Ville distribution) was demonstrated. A major advantage of time-frequency methods is the ability to analyze acoustic signals containing multiple propagation modes, which overlap and superimpose in the time domain signal. By combining time-frequency analysis with a broadband acoustic excitation source, the dispersion of multiple Lamb modes over a wide frequency range can be determined from as little as a single measurement. In addition, the technique provides a direct measurement of the group velocity dispersion. The technique was first demonstrated in the analysis of a simulated waveform in an aluminum plate in which the Lamb mode dispersion was well known. Portions of the dispersion curves of the A0, A1, S0, and S2 Lamb modes were obtained from this one waveform. The technique was also applied for the analysis of experimental waveforms from a unidirectional graphite/epoxy composite plate. Measurements were made both along and perpendicular to the fiber direction. In this case, the signals contained only the lowest order symmetric and antisymmetric modes. A least squares fit of the results from several source to detector distances was used. Theoretical dispersion curves were calculated and are shown to be in good agreement with experimental 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.     

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

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

小样本字典学习的兰姆波模态识别方法*

针对兰姆波多模态识别问题,提出了基于小样本字典学习的模态识别方法。将多层复合板的频散特性看作一个线性时不变系统,首先,根据频散知识模拟各个模态传播特定距离后的信号,提取走时和能量特征创建字典;其次,获取待测信号的走时特征,通过查询字典来识别兰姆波模态;最后,根据能量参数估计结果,实现待测信号中各模态信号的分离和重构。通过对三层粘接的AAA板(铝板-亚克力板-铝板,每层厚度为2mm)中传播距离为0.3m、0.5m的直达波和反射波的实验验证,结果显示该方法对A0、S0模态的有效识别和各个波包信号的准确重建。     

非线性兰姆波在厚度缓慢变化和衰减下的特性分析

该文运用解析的方式推导了考虑声波衰减时兰姆波二次谐波的累积和传播规律,并用半解析方式将该理论推广到缓慢变厚度板的情况。由于色散特性,兰姆波二次谐波和基频波相速度不匹配,传播通常会产生拍频效应,使得二次谐波的振幅沿着传播距离周期性的归零。当考虑声波衰减或板的厚度缓慢变化的情况时,拍频效应将不再严格地被满足。二次谐波的振幅依然会沿着传播距离而振荡,但不会归零。该研究可以用于分析如何高效地激发和接收兰姆波的二次谐波,表征和评估不同厚度变化的结构中的微观结构损伤。     

Size-dependent dispersion characteristics in piezoelectric nanoplates with surface effects

In this paper, surface effects on the dispersion characteristics of elastic waves propagating in an infinite piezoelectric nanoplate are investigated by using the surface piezoelectricity model. Based on the surface piezoelectric constitutive theory, the presence of surface stresses and surface electric displacements exerting on the boundary conditions of the piezoelectric nanoplate is taken into account in the modified mechanical and electric equilibrium relations. The partial wave technique is employed to obtain the general solutions of governing equations, and the dispersion relations with surface effects are expressed in an explicit closed form. The impacts of surface piezoelectricity, residual surface stress and plate thickness on the propagation properties of elastic waves are analyzed in detail. Numerical results show that the dispersion behaviors in piezoelectric nanoplates are size-dependent, and there exists a critical plate thickness above which the surface effects may vanish.     

Shear-horizontal waves in a rotated Y-cut quartz plate in contact with a viscous fluid

We study shear-horizontal (SH) waves in a crystal plate of rotated Y-cut quartz in contact with a semi-infinite viscous fluid. The crystal plate and the fluid are governed by the equations of anisotropic elasticity and the theory of Newtonian fluids. A transcendental equation that determines the dispersion relations of the waves is obtained. Approximate analytical solutions to the equation are presented for the case of low viscosity fluids and the case of long waves whose wavelength is much larger than the plate thickness. The effects of the fluid viscosity and density on the dispersion relations of the waves are examined. The results obtained are fundamental and useful to the understanding and design of acoustic wave fluid sensors for measuring fluid viscosity or density.