分层半空间表面非线性瑞利波的激发*

针对非线性瑞利波在均匀分层半空间结构中的激发和传播规律进行研究。根据摄动理论和模态分解将分层半空间结构中瑞利波的二次谐波声场表示为二倍频瑞利波模式的线性组合,经由互易关系得到各模式的展开系数表达式。对不同分层半空间结构中瑞利波二次谐波的激发和传播特性进行讨论,结果表明相速度匹配的瑞利波模式其二次谐波分量随传播距离线性增长,非匹配模式的二次谐波分量则沿传播方向周期震荡传播。此外,文中定义非线性参数表征瑞利波模式产生的非线性程度,这有利于选择出具有明显非线性效应的匹配点,为后续检测工作提供理论依据,具有指导意义。     

Evaluation of the acoustic nonlinearity parameter of materials with Rayleigh waves excited directly

An experimental technique for the determination of the relative acoustic nonnnearity parameter of materials with Rayleigh waves excited directly is presented. Rayleigh surface waves were directly generated in materials by the specific piezoelectric transducer fixed at the specimen＇s edge, and were measured with a laser interferometer system. After the Rayleigh wave signals were processed with FFT method, the relative acoustic nonlinearity parameter of materials was then determined from the absolute magnitudes of the fundamental and second harmonics of Rayleigh surface waves. This procedure was used to determine the acoustic non- linearity parameters of aluminum alloys 2024 and 6061. It is shown that the results comply well with those available in the literature; this method can thus be used to evaluate the acous- tic nonlinearity parameter of materials effectively. This technique can provide a practical way in the nondestructive characterization of degradation of materials and structures in the early fatigue life.     

Second Harmonic Generation of Lamb Wave in Numerical Perspective

The influences of phase and group velocity matching on cumulative second harmonic generation of Lamb waves are investigated in numerical perspective.Finite element simulations of nonlinear Lamb wave propagation are performed for Lamb wave mode pairs with exact and approximate phase velocity matching,with and without group velocity matching,respectively.The evolution of time-domain second harmonic Lamb waves is analyzed with the propagation distance.The amplitudes of primary and second harmonic waves are calculated to characterize the acoustic nonlinearity.The results verify that phase velocity matching is necessary for generation of the cumulative second harmonic Lamb wave in numerical perspective,while group velocity matching is demonstrated to not be a necessary condition.     

Estimation methodology for amplitude of second-order harmonic in nonlinear surface acoustic wave measured by optical diffraction method

An estimation methodology for amplitude from the fundamental to fourth-order harmonic composed of nonlinear surface acoustic waves on a LiNbO₃ substrate with a 128° Y-X propagation axis was demonstrated. First, equations for estimation of amplitudes from diffracted light intensity were shown and then numerical simulations of main factors were carried out. Normalized diffraction light intensities due to components of nonlinear surface acoustic waves were measured and their amplitudes were obtained from measured values. Finally, the amplitude of the second harmonic propagating on a gold thin film of about 50 nm thickness was estimated. It was found that the amplitude of the second harmonic increased continuously from 0.8 to 2.0 Å with propagation distance in the thin film.     

固体板中SH板波非线性效应的实验观察

采用微扰近似和导波的模式展开分析方法,从理论上简要分析了SH板波的二次谐波发生效应;尽管在无限大固体介质中单个切变波的二次谐波发生效应非常微弱,但在一定条件下由两个切变波构成的SH板波可具有强烈的非线性效应;本文的主要工作就是对此结论加以实验验证。试制了激发SH板波的切变波斜劈换能器和接收二次谐波信号的液体斜劈换能器,建立了非线性SH板波的实验研究系统;通过详细的理论分析和对比实验研究,阐明了在一定条件下实验观察到的显著二次谐波信号来源于SH板波传播过程中的强烈非线性效应。此外,针对不同的SH板波传播距离,在远场条件下分别测量了相应的二次谐波幅频曲线;在基频SH板波与二倍频对称兰姆波相速度相等所对应的频率值附近,分析了二次谐波的振幅随传播距离的变化关系,结果证明在一定条件下SH板波的二次谐波振幅可随传播距离积累增长,即SH板波可具有强烈的非线性效应。      

Testing for inhomogeneity of the nonlinear parameter in an acoustic medium by ultrasonic phase conjugation

The principle of applying a selective phase conjugation of the second harmonic of a focused ultrasonic wave to diagnosing inhomogeneity of the nonlinear parameter in an acoustic medium is considered. A solution to the three-dimensional problem of harmonic generation by phase-conjugated waves in a nonlinear medium with a localized isoechogenous inclusion is obtained. The signal amplitudes detected by a transmitting-receiving transducer at the second and forth harmonics of a probe wave are calculated for varying position of the inclusion relative to the focus.     

非线性电磁超声对铝合金拉伸变形评价研究

非线性系数是描述材料中微纳尺度损伤的特征参量,非线性系数常通过接触式压电超声进行检测,但耦合剂引起的非线性一般是未知的,针对这一问题,提出了一种非接触式电磁超声非线性纵波检测方法。该方法基于洛伦兹力机理在试件表面产生的振动弹性波,利用不同拉伸载荷下所制备的损伤试件,分别利用压电超声、电磁超声进行非线性超声系数测量。实验结果表明:利用两种非线性超声检测的相对非线性系数与铝合金的拉伸形变呈单调关系,同时也论证了电磁超声纵波基于非线性理论对塑性变形评估的可行性。      

Nonlinear lamb waves in a metal plate with defects

Experimental results on the propagation of finite-amplitude Lamb waves in a solid plate made of polycrystalline aluminum alloy with defects are presented. The Lamb waves are recorded and visualized using a scanning laser vibrometer. The dependences of the higher harmonic amplitudes, both averaged over the plate surface and measured at a point far from the defect and at the site of the defect, on the fundamental frequency amplitude are studied. A threshold character of the higher harmonic generation and a power-law behavior of their amplitudes are revealed, the latter feature being unconventional for the nonlinearity associated with the anharmonicity of the crystal lattice of the material. The possibility of locating the distribution of individual defects from the measured spatial distribution of structural nonlinearity in the material under study is experimentally demonstrated. The results of the experiments are explained in terms of the bilinear medium model.     

Effects of harmonic phase on nonlinear surface acoustic waves in the (111) plane of cubic crystals

Spectral evolution equations are used to perform numerical studies of nonlinear surface acoustic waves in the (111) plane of several nonpiezoelectric cubic crystals. Nonlinearity matrix elements which describe the coupling of harmonic interactions are used to characterize velocity waveform distortion. In contrast to isotropic solids and the (001) plane of cubic crystals, the nonlinearity matrix elements usually cannot be written in a real-valued form. As a result, the harmonic components are not necessarily in phase, and dramatic variations in waveforms and propagation curves can be observed. Simulations are performed for initially monofrequency surface waves. In some directions the waveforms distort in a manner similar to nonlinear Rayleigh waves, while in other directions the velocity waveforms distort asymmetrically and the formation of shocks and cusped peaks is less distinct. In some cases, oscillations occur near the shocks and peaks because of phase differences between harmonics. A mathematical transformation based on the phase of the matrix elements is shown to provide a reasonable approximation of asymmetric waveform distortion in cases where the matrix elements have similar phase.     

Baseline-free estimation of residual fatigue life using a third order acoustic nonlinear parameter

Prediction of crack growth and fatigue life estimation of metals using linear/nonlinear acousto-ultrasound methods is an ongoing issue. It is known that by measuring nonlinear parameters, the relative accumulated fatigue damage can be evaluated. However, there is still a need to measure two crack propagation states to assess the absolute residual fatigue life. A procedure based on the measurement of a third-order acoustic nonlinear parameter is presented to assess the residual fatigue life of a metallic component without the need of a baseline. The analytical evaluation of how the cubic nonlinear-parameter evolves during crack propagation is presented by combining the Paris law to the Nazarov-Sutin crack equation. Unlike other developed models, the proposed model assumes a crack surface topology with variable geometrical parameters. Measurements of the cubic nonlinearity parameter on AA2024-T351 specimens demonstrated high sensitivity to crack propagation and excellent agreement with the predicted theoretical behavior. The advantages of using the cubic nonlinearity parameter for fatigue cracks on metals are discussed by comparing the relevant results of a quadratic nonlinear parameter. Then the methodology to estimate crack size and residual fatigue life without the need of a baseline is presented, and advantages and limitations are discussed.     

非线性声波方程二次谐波高阶近似解及实验验证

非线性声波方程的传统二次谐波二阶近似解不够精确,对测量材料的非线性系数造成较大误差。利用摄动法展开非线性声波方程得到一系列非齐次偏微分方程,根据低阶解的性质拟定高阶特解的形式,通过符号计算求出高次谐波特解,对所有二次谐波成分求和最终获得二次谐波的高阶近似解。在水中开展非线性声学实验验证高阶近似解,结果表明:二次谐波相对幅度随传播距离和输入信号均呈现先增加后减小的趋势,高阶近似解即使在传播距离不太短和输入信号不太弱的宽泛条件下也能更好地吻合实验结果。得到的高阶近似解不仅弥补了现有二次谐波幅度理论公式的不足,进行非线性声学实验时还可扩展测量范围,充分利用实验数据也能提高非线性系数的测量精度。      

Nonlinear surface waves and solitons

First a general introduction on the notion of surface waves on solids (types of different waves), a reminder on the simplest familiar nonlinear dispersive model equations, and another on the basic equations of nonlinear elasticity are given. Then attention is focused on the linear surface wave problem. The main properties of nonlinear surface waves in the absence of dispersion are studied next by use of several asymptotic techniques. The additional effects of dispersion are then considered and combined with those of nonlinearity with an emphasis on the case of so-called shear-horizontal surface waves and solitary-wave solutions for envelope signals. Finally, typical nonlocality is introduced for nonlinear Rayleigh surface waves, and general comments on more general two-dimensional (in propagation space) nonlinear strain waves on structures are evoked by way of conclusion.     

Two-dimensional modeling of wave propagation in materials with hysteretic nonlinearity

A multiscale model for the two-dimensional nonlinear wave propagation in a locally microdamaged medium is presented, and numerical simulations are analyzed in view of nondestructive testing applications. The multiscale model uses a statistical distribution of hysterons and upscales their microscopic stress-strain relations to a mesoscopic level. Macroscopic observations are then predicted by finite integration techniques. The influence of a small region with hysteretic nonlinearity on the generation of harmonics is investigated, and numerical results for different amplitudes of the input signal and different analysis techniques of the response signal are presented. Second, a study is conducted on the interaction of a Rayleigh wave with a microdamaged zone with hysteretic nonlinearity at the surface of an otherwise linear body, and the influence of the microdamaged zone on the surface wave velocity and on the generation of harmonics is examined. It is found that the effect of hysteresis on the Rayleigh wave propagation can be barely seen in the surface wave velocity measurement, but shows up nicely in the wave spectrum. The potential of a nonlinearity based depth profiling technique is explored by evaluating the nonlinear responses at different frequencies for a vertically stratified medium with spatially varying hysteresis properties.     

一维非线性声波传播特性

针对一维非线性声波的传播问题进行了有限元仿真和实验研究. 首先推导了一维非线性声波方程的有限元形式, 含有高阶矩阵的非线性项导致声波具有波形畸变、谐波滋生、基频信号能量向高次谐波传递等非线性特性. 编制有限元程序对一维非线性声波进行了计算并对仿真得到的畸变非线性声波信号进行处理, 分析其传播性质和物理意义. 为验证有限元计算结果, 开展了水中的非线性声波传播的实验研究, 得到了不同输入信号幅度激励下和不同传播距离的畸变非线性声波信号. 然后对基波和二次谐波的传播性质进行详细讨论, 分析了二次谐波幅度与传播距离和输入信号幅度的变化关系及其意义, 拟合出二次谐波幅度随传播距离变化的方程并阐述了拟合方程的物理意义. 结果表明, 数值仿真信号及其频谱均与实验结果有较好的一致性, 证实计算方法和结果的正确性, 并提出了具有一定物理意义的二次谐波随传播距离变化的简单数学关系. 最后还对固体中的非线性声波传播性质进行了初步探讨. 本研究工作可为流体介质中的非线性声传播问题提供理论和实验依据.     

Nonlinear acoustic response through minute surface cracks: FEM simulation and experimentation

The second harmonic of a Rayleigh wave passing through a minute surface crack has been numerically analyzed by semi-explicit FEM including special elements which account for a nonlinear stress-strain relation at crack surfaces. Minute cracks perpendicular to a free, flat surface close under compressive stress when width of the crack opening is less than the longitudinal amplitude of the Rayleigh wave. Thereafter, compressive and shear stresses are partially transmitted through the closed cracks, whereas tensile and shear stresses are not transmitted through cracks that remain open. This leads to marked nonlinear ultrasonic response. Calculation was performed for an aluminum block having a surface crack. The transverse component of the Rayleigh wave propagating through the cracks shows distorted waveforms, making the second harmonic amplitude clearly noticeable. In an experiment, the second harmonic component of the leaky Rayleigh wave was detected for a simple crack model consisting of two aluminum blocks, by use of a PVDF line-focused transducer. The results of the experiment show that the second harmonic amplitude is a second-order function of the fundamental wave amplitude, and is more pronounced for low compressive stress applied to close the crack surfaces.     

Second-harmonic generation in a layer with variable susceptibility

The weak nonlinearity limit of the second-harmonic generation in a medium with variable second-order susceptibility, changing along the direction of propagation of interacting waves according to the hyperbolic secant law, is considered. It is shown that the variation of the second harmonic’s normalized intensity with the spatial coordinate is given in terms of the solution to an auxiliary linear problem, which is known as the Rosen-Zener quantum two-level model. The final intensity of the second harmonic at the exit from the medium is calculated and analyzed. It is shown that, for the particular Rosen-Zener profile under consideration, because of its inherent properties, the variation range of the peak susceptibility of the layer for which the medium can be strictly considered as weakly nonlinear is narrower than the corresponding parameter range for a homogeneous medium.     

Propagation of acoustic wave in viscoelastic medium permeated with air bubbles

Based on the modification of the radial pulsation equation of an individual bubble, an effective medium method (EMM) is presented for studying propagation of linear and nonlinear longitudinal acoustic waves in viscoelastic medium permeated with air bubbles. A classical theory developed previously by Gaunaurd (Gaunaurd GC and \"{U}berall H, {\em J. Acoust. Soc. Am}., 1978; 63: 1699--1711) is employed to verify the EMM under linear approximation by comparing the dynamic (i.e. frequency-dependent) effective parameters, and an excellent agreement is obtained. The propagation of longitudinal waves is hereby studied in detail. The results illustrate that the nonlinear pulsation of bubbles serves as the source of second harmonic wave and the sound energy has the tendency to be transferred to second harmonic wave. Therefore the sound attenuation and acoustic nonlinearity of the viscoelastic matrix are remarkably enhanced due to the system's resonance induced by the existence of bubbles.     

Simulation of laser-generated ultrasonic wave propagation in solid media and air with application to NDE

Ultrasonic methods are well known as powerful and reliable tool for defect detection. In the previous decades focus and interest have been directed to non-contact sensors and methods, showing many advantages over contact techniques where inspection depends on contact conditions (pressure, coupling medium, contact area). The non-contact hybrid ultrasonic method described here is of interest for many applications, requiring periodic inspection in service or after manufacturing. Despite the potential impact of laser-generated ultrasound in many areas of industry, robust tools for studying the phenomenon are lacking and thus limit the design and optimization of non-destructive testing and evaluation techniques. Here a specific numerical method is presented to efficiently and accurately solve ultrasound wave propagation problems with frequencies in the MHz range traveling in relatively large bodies and through air. This work improves a previous numerical model where propagation of the acoustic waves through air had not been considered, allowing us to simulate the presence of a non-contact transducer in reception in order to simulate numerically the complete experimental setup. It is very important to limit the amount of air to be considered in the FE analyses; otherwise the computational cost would often exceed the resources available. A way to solve the problem is to implement non-reflecting boundary conditions. A non-reflecting boundary condition allows all outgoing waves to exit the domain at the boundary where they have been imposed without reflection; thus, it is possible to model only the portion of air between the non-contact transducer and the solid under testing. Several numerical and experimental analyses were conducted on a 136 lb AREMA rail; here we study in detail two fully non-contact testing configurations for the rail head and web. The information that can be acquired is very valuable for choosing the right setup and configuration when performing non-contact hybrid ultrasonic inspection.     

Harmonic generation of an obliquely incident ultrasonic wave in solid-solid contact interfaces

The conventional acoustic nonlinear technique to evaluate the contact acoustic nonlinearity (CAN) at solid-solid contact interfaces (e.g., closed cracks), which uses the through-transmission of normally incident bulk waves, is limited in that access to both the inner and outer surfaces of structures for attaching pulsing and receiving transducers is difficult. The angle beam incidence and reflection technique, where both the pulsing and receiving transducers are located on the same side of the target, may allow the above problem to be overcome. However, in the angle incidence technique, mode-conversion at the contact interfaces as well as the normal and tangential interface stiffness should be taken into account. Based on the linear and nonlinear contact stiffness, we propose a theoretical model for the reflection of an ultrasonic wave angularly incident on contact interfaces. In addition, the magnitude of the CAN-induced second harmonic wave in the reflected ultrasonic wave is predicted. Experimental results obtained for the contact interfaces of A16061-T6 alloy specimens at various loading pressures showed good agreement with theoretical predictions. Such agreement proves the validity of the suggested oblique incidence model.     

矩形表面波探头声场的高斯声束叠加法

利用矩形压电晶片和有机玻璃楔块折射可激励出超声表面波,广泛用于固体近表面缺陷检测和材料特性测量.由于描述表面波三维声场的理论方法还鲜有报道,因而主要采用简化的表面波二维声场模型来定量分析这类问题.高斯声束模型近些年被广泛应用于解决超声体波传播的各种复杂问题,然而,目前还没有将其扩展应用到超声表面波的声场的计算中.通过结合表面波格林方程和矩形换能器的高斯声束模型,推导出基于高斯声束叠加的表面波三维声场解析解.进一步,将该方法与点源叠加的数值解进行了分析比较,计算结果表明表面波声场的高斯声束叠加方法在具有较好计算精度的同时,还具有更快的计算效率.