谱有限元模拟正交各向异性黏弹性材料中导波传播特性

研究了导波在正交各向异性黏弹性复合板中传播的色散特性、波结构及功率流密度。基于二维平面运动方程,采用谱有限元方法得到了导波色散的特征方程。分析了正交各向异性黏弹性板中各向异性和黏性对能量速度和波结构的影响,以及基底对导波功率流密度的影响。数值研究结果表明:导波沿纤维方向传播的能量速度大,材料的黏性对速度影响较小,但会减小波结构的幅度;在高频时,基底的存在使两个基本模态的功率流密度分别集中到波导的上下表面,形成弱色散、高衰减及无色散、零衰减的表面波。数值模拟结果为导波用于复合材料定量无损检测和性能评价提供理论依据。      

Circumferential thermoelastic waves in orthotropic cylindrical curved plates without energy dissipation

In this article, the propagation of guided thermoelastic waves in the circumferential direction of orthotropic cylindrical curved plates subjected to stress-free, isothermal boundary conditions is investigated in the context of the Green-Naghdi (GN) generalized thermoelastic theory (without energy dissipation). The coupled wave equations and heat conduction equation are solved by the Legendre orthogonal polynomial series expansion approach. The convergency of the method is discussed through a numerical example. The dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Dispersion curves of the corresponding pure elastic cylindrical plate are also shown to analyze the influence of the thermoelasticity on elastic modes. The displacement, temperature and stress distributions are shown to discuss the differences between the elastic modes and thermal modes. A thermoelastic cylindrical plate with a different ratio of radius to thickness is considered to discuss the influence of the ratio on the characteristics of circumferential thermoelastic waves.     

Formulas for the H/V ratio (ellipticity) of Rayleigh waves in orthotropic elastic half-spaces

ABSTRACT

This paper concerns the propagation of Rayleigh waves in compressible and incompressible orthotropic elastic half-spaces. The main aim of the paper is to derive formulas for the H/V ratio (ellipticity) of Rayleigh waves. First, the expression of the H/V ratio in terms of the Rayleigh wave velocity for compressible orthotropic elastic half-spaces is derived using the surface impedance matrix. From this expression and the secular equation of Rayleigh waves, the equation for the H/V ratio is obtained. It is a cubic equation. Employing the theory of cubic equations, the formula for the H/V ratio for compressible orthotropic elastic half-spaces has been derived. The H/V ratio formula for incompressible orthotropic elastic half-spaces has been obtained using the incompressible limit method. Since the H/V ratio of Rayleigh waves is a good tool for nondestructively evaluating mechanical properties of structures before and during loading, the obtained formulas will be very significant in practical applications.     

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.     

Guided circumferential waves in orthotropic cylindrical curved plate and the mode conversion by the end-reflection

In ultrasonic nondestructive inspection of large-diameter pipes and curved plate, longitudinal cracks are detected more efficiently by using guided circumferential waves. In the present, the study of guided circumferential waves and their application in detecting longitudinal defect were relative adequate when pipe material is isotropic. Based on linear three-dimensional elasticity, an orthogonal polynomial series expansions approach is used for determining the guided circumferential waves dispersion curves in homogeneous infinitely long orthotropic hollow cylinders. Results are compared with those published earlier and with the finite element simulation to check up the accuracy and range of applicability of this polynomial approach. Through the analysis of the displacements distributions and finite element simulation, the mode conversion of guided circumferential waves by end-reflection in cylindrical curved plate is discovered.     

Research on the propagation of guided waves in multi-layered plates

The dispersion characteristics and excitation mechanisms of the guided waves in multilayered plates are studied in this paper. Firstly, the dispersion equation is obtained by the propagator matrix method. Then, the bisection technique is employed to find all the roots of the dispersion equation. The dispersion characteristics of the guided waves are investigated and analyzed. For the multilayered plates in which the S-wave velocity increases or decreases from top to bottom, it is found that the velocity limits in high frequency of the first and high modes are equal to the Rayleigh wave and S-wave velocities of the low-velocity layer, respectively. It is also found that the velocity limits in the high frequency of all modes are equal to the S-wave velocity of the low-velocity layer for the plate with a low-velocity middle layer. The normal displacement spectrum of all the modes excited by the normal force source with a definite width on the surface of the plate is also investigated. It is proved that the dominant mode is the first mode when the S-wave velocity increases from top to bottom layer and the dominant mode is different in different frequency range for the plate with a low-velocity middle layer.     

Power flow model of flexural waves in finite orthotropic plates

In this paper, an approximate power flow model is developed for the analysis of the flexural waves in finite orthotropic plate transversely vibrating in the medium-to-high frequency ranges. The derived energy equation for the model is expressed with time- and locally space-averaged farfield wave energy density. It could be the more general form than that of the conventional power flow analysis model seen in the isotropic plate. With the derived model, dynamic characteristics varying with the direction can be expressed. To verify the validity and accuracy of the model, numerical analyses are performed for the case where a finite rectangular plate is excited by a harmonic point force, and the calculated results expressed with the energy levels are compared with classical modal solutions by changing the frequency and the damping loss factor of the plate. The dominant power transmission paths in the plate are also predicted from the distribution of the approximate intensity fields.     

On a technique for deriving the explicit secular equation of Rayleigh waves in an orthotropic half-space coated by an orthotropic layer

The secular equation of Rayleigh propagating in an orthotropic half-space coated by an orthotropic layer has been obtained by Sotiropolous [Sotiropolous, D. A. (1999), The e®ect of anisotropy on guided elastic waves in a layered half-space, Mechanics of Materials 31, 215–233] and by Sotiropolous & Tougelidis [Sotiropolous, D. A. and Tougelidis, G. (1998), Guided elastic waves in orthotropic surface layer, Ultrasonics 36, 371–374]. However, it is not totally explicit and some misprints have occurred in this secular equation in both papers. This secular equation was derived by expanding directly a six-order determinant originated from the traction-free conditions at the top surface of the layer and the continuity of displacements and stresses through the interface between the layer and the half-space. Since the expansion of this six-order determinant was not shown in both two papers, it has been difficult to readers to recognize these misprints. This paper presents a technique that provides a totally explicit secular equation of the wave. The technique makes clear the way from the traction-free and continuity conditions to the secular equation and enables us to recognize the misprints appearing in the reported secular equation. The technique can be employed to obtain explicit secular equations of Rayleigh waves for many other cases. Moreover, the paper introduces a transfer matrix in explicit form for an orthotropic layer that is much simpler in form than the one obtained previously.     

Study of Rayleigh type surface wave in the initial stressed irregular bottom of ocean due to point source

Rayleigh type surface wave propagation in the irregular bottom of ocean model which is the interface of homogeneous liquid layer over laying an irregular boundary of homogeneous orthotropic half space under initial stresses has been discussed in this paper. Three different dispersion equations are obtained in the form of simple equation using and not using Perturbation technique. Some special cases have been considered. The effect of irregularity, initial stressed, point source, and depth of liquid layer on the propagation of Rayleigh waves has been analyzed and results of numerical discussion have been presented graphically for three different dispersion equations. Mainly the graphs are shown the variation of phase velocity with wave number in different cases.     

Normal mode expansion method for laser-generated ultrasonic Lamb waves in orthotropic thin plates

A quantitative theory for modeling the laser-generated transient ultrasonic Lamb waves, which propagates along arbitrary directions in orthotropic plates, is presented by employing an expansion method of generalized Lamb wave modes. The displacement field is expressed by a summation of the symmetric and antisymmetric modes in the surface stress-free orthotropic plate, and therefore the theory is particularly appropriate for waveform analyses of Lamb waves in thin plates because one needs only to evaluate several lower modes. The transient waveforms excited by the thermoelastic expansion and the oil-coating evaporation are analyzed for a transversely isotropic thin plate. The results show that the theory provides a quantitative analysis to characterize anisotropic elastic stiffness properties of orthotropic plates by laser-generated Lamb wave detection.     

Flexural edge waves in semi-infinite elastic plates

This paper presents a detailed analysis of the dispersion for flexural edge waves in semi-infinite isotropic elastic plates. A solution to the dynamic equations of motion is constructed by the superposition of two partial solutions, each providing zero shear stresses at the plate faces. A dispersion equation is expressed via the determinant of an infinite system of linear algebraic equations. The system is reduced to a finite one by taking into account the asymptotic behaviour of unknown coefficients. The accuracy of the solution is confirmed by a good agreement with the available experimental data and by a proper satisfaction of the prescribed boundary conditions.A detailed analysis of dispersion properties for the edge wave and corresponding displacements at various frequencies is carried out. In addition to the well-known results it is shown that the plate height does not influence the existence of the edge wave at high frequencies and, as the frequency increases, the phase velocity of the edge wave in a semi-infinite plate asymptotically approaches the velocity of an edge wave in a right-angled wedge. The performed analysis allows evaluating the plate theories such as the Kirchhoff theory or other refined plate theories developed for modeling edge waves in semi-infinite elastic plates at low frequencies.     

Generalized thermoelastic extensional and flexural wave motions in homogenous isotropic plate by using asymptotic method

In this paper the asymptotic method has been applied to investigate propagation of generalized thermoelastic waves in an infinite homogenous isotropic plate. The governing equations for the extensional, transversal and flexural motions are derived from the system of three-dimensional dynamical equations of linear theories of generalized thermoelasticity. The asymptotic operator plate model for extensional and flexural free vibrations in a homogenous thermoelastic plate leads to sixth and fifth degree polynomial secular equations, respectively. These secular equations govern frequency and phase velocity of various possible modes of wave propagation at all wavelengths. The velocity dispersion equations for extensional and flexural wave motion are deduced from the three-dimensional analog of Rayleigh-Lamb frequency equation for thermoelastic plate. The approximation for long and short waves along with expression for group velocity has also been obtained. The Rayleigh-Lamb frequency equations for the considered plate are expanded in power series in order to obtain polynomial frequency and velocity dispersion relations and its equivalence established with that of asymptotic method. The numeric values for phase velocity, group velocity and attenuation coefficients has also been obtained using MATHCAD software and are shown graphically for extensional and flexural waves in generalized theories of thermoelastic plate for solid helium material.     

正交异性钢桥面板声发射波三维谱元法模拟及损伤定位*

针对大跨度桥梁正交异性钢桥面板的疲劳损伤评估与结构健康监测需求,开展基于声发射波场谱元法模拟的大型复杂板类结构损伤定位研究。采用Legendre高阶插值三维时域谱元法模拟声发射波在正交异性钢桥面板中的传播过程,验证了其内部显著的反射、衍射和频散现象,并代替人工预断铅实测试验获得大量声发射数据。然后,利用赤池信息准则判定声发射波到达各传感器的时间,通过高斯过程回归建立到达时差与声发射源位置的关系模型,用于未知损伤的定位监测。数值模型实验结果表明,赤池信息准则和高斯过程回归改进的时差图法在正交异性钢桥面板中的平均定位误差为37.3 mm(25 dB信噪比工况),平板的定位精度高于U肋。谱元法模拟有望代替繁琐的预断铅实测试验,提升声发射时差图系列损伤定位方法的实用性。     

Localized bending waves in a transversely isotropic plate

The problem of bending waves localized near the free edge of a transversely isotropic plate is investigated using the Ambartsumian higher-order plate theory which takes account of the transverse shears generated by flexural deformation. Unlike the first-order Reissner-Mindlin theory, which also takes account of transverse shears, Ambartsumian's analysis does not demand that plane normal cross-sections remain plane during bending. Within this analysis the existence of localized bending waves in transversely isotropic plates is established, and solutions of the dispersion equation obtained for different values of the elastic parameters.The analysis of frequencies of localized bending waves shows that for thick plates the effect of anisotropy can be considerable. For the particular case of vibrations of a narrow plate, from the long wave approximation a new beam vibration equation of the Timoshenko type is obtained for a transversally isotropic plate.     

铝板中Lamb波检测的实验研究

Lamb波在传播过程中具有频散及多模态特性,若相关参数选择不当,会导致在实际应用中信号相互叠加而无法识别。本文基于Lamb波的频散曲线是其频散方程实数解分布的特点,采用二分法绘制了铝板中Lamb波的频散曲线、波结构曲线和入射角曲线。根据曲线选择S0模态的Lamb波对1mm厚铝板中不同类型的缺陷进行检测。实验结果表明,S0模态的Lamb波对裂纹型缺陷和贯穿型缺陷十分敏感,但对于裂纹型缺陷,其幅值变化并不与缺陷大小成线性关系,并且S0模态Lamb波的声场指向性十分集中,在偏离声束轴线时无法检测到缺陷。     

FREE VIBRATION OF ORTHOTROPIC SQUARE PLATES WITH A SQUARE HOLE

An approximate method is extended for analyzing the free vibration problem of orthotropic square plate with a square hole. In this paper, a square plate with a square hole is transformed into an equivalent square plate with non-uniform thickness by considering the hole as an extremely thin part of the equivalent plate. Therefore, the dynamic characteristics of a plate with a hole can be obtained by analyzing the equivalent plate. The Green function, which is the discrete solution for the deflection of the equivalent plate, is used to obtain the characteristic equation of the free vibration. The effects of the side to thickness ratio, hole side to plate side ratio and the variation of the thickness on the frequencies are considered. Some numerical analyses are carried out for the simply supported orthotropic square plate with a square hole. The efficiency and accuracy of the numerical solutions by the present method are investigated.     

Measurements of dispersion curves of leaky Lamb waves using a lensless line-focus transducer.

This paper presents a new method of measuring dispersion curves of leaky Lamb waves propagating in a thin plate. The measurement system is based on a lens-less line-focus transducer and its defocusing measurement technique. A new waveform processing method known as V(f, z) method is introduced for analyzing the measured waveforms and for extracting the dispersion relation. For two thin-plate samples, namely a 150 microm thick alumina plate and a 320 microm thick stainless steel plate, the dispersion curves of several lower order modes of leaky Lamb waves have been experimentally determined in the frequency range of few to 25 MHz. The experimental dispersion curves have then been compared with the theoretical ones. Excellent agreement is observed. It is then verify that this method indeed provides a convenient and accurate way for experimentally measuring dispersion curves of leaky Lamb waves of thin-plate samples.     

Localized bending vibrations of piezoelectric plates

This paper investigates the existence and propagation of electro-elastic bending waves localized at the free edge of a piezoelectric plate. The problem is considered within the framework of the high-order refined plate theory introduced by Ambartsumian. The condition for existence of a localized bending wave is obtained, and the dispersion equation solved with respect to a dimensionless frequency. It is shown that the piezoelectric effect can increase the attenuation coefficient for a localized wave by up to 70% compared with that for a purely elastic plate, thus significantly decreasing the depth of penetration. The problem is also solved within the classical Kirchhoff theory. A comparison of results is carried out between two theories.     

Effects of dispersion on the inference of metal texture from S0 plate mode measurements. Part I. Evaluation of dispersion correction methods.

Ultrasonic S0 waves (fundamental symmetric Lamb modes) are being considered in several laboratories for the nondestructive characterization of the texture (preferred grain orientation) and formability of metal sheets and plates. In a typical experimental setup, the velocities of the S0 waves are measured as a function of wave propagation angle with respect to the rolling direction of the plate. However, the S0 waves are known to be dispersive, and that dispersion must be considered in order to isolate the small, texture-induced shifts in the S0 wave velocity. Currently, there are two approximate dispersion correction methods, one proposed by Thompson et al. [Met. Trans. A 20, 2431 (1989)] and the other introduced by Hirao and Fukuoka [J. Acoust. Soc. Am. 85, 2311 (1989)]. In this paper, these two methods will be evaluated using an exact theory for wave propagation in orthotropic plates. Through the evaluation, the limits of the current texture measurement techniques are established. It is found that when plate thickness to wavelength ratio is less than 0.15, both Thompson's and Hirao's methods work satisfactorily. When the thickness to wavelength ratio exceeds 0.3, neither Thompson's nor Hirao's dispersion correction method provides adequate corrections for the current texture measurement techniques. Within the range of 0.15-0.3, Thompson's method is recommended for weakly anisotropic sheets and plates and Hirao's method may be more appropriate for some strongly anisotropic cases.     

Flexural and transversal wave motion in homogeneous isotropic thermoelastic plates by using asymptotic method

The present investigation is concerned with the flexural and transversal wave motion in an infinite, transversely isotropic, thermoelastic plate by asymptotic method. The governing equations for the flexural and transversal motions have been derived from the system of three-dimensional dynamical equations of linear theory of coupled thermoelasticity. The asymptotic operator plate model for free vibrations; both flexural and transversal, in a homogenous thermoelastic plate leads to fifth degree and cubic polynomial secular equations, respectively, that governs frequency and phase velocity of various possible modes of wave propagation at all wavelengths. All the coefficients of differential operator have been expressed as explicit functions of the material parameters. The velocity dispersion equations for the flexural and transversal wave motion have been deduced from the three-dimensional analog of Rayleigh-Lamb frequency equation for thermoelastic plate waves. The approximations for long and short waves and expression for group velocity have also been derived. The thermoelastic Rayleigh-Lamb frequency equations for the considered plate are expanded in power series in order to obtain polynomial frequency and velocity dispersion relations whose equivalence is established with that of asymptotic method. The dispersion curves for phase velocity, group velocity and attenuation coefficient of various flexural and transversal wave modes are shown graphically for aluminum-epoxy material elastic and thermoelastic plates.