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

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

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

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

北极冰-水耦合系统中声波传播特性

利用固体和流体介质中波传播理论,导出了冰-水两层复合结构中导波频散方程。进一步,利用二分法对频散方程进行了数值求解,得到了ω-k频散曲线(ω与k分别为圆频率和波数),以及相速度和群速度频散曲线。结果表明:冰-水两层复合结构中导波由具有相同厚度水层和冰层中导波耦合而成,但与水层和冰层中导波频散曲线相比,复合结构中导波频散曲线除第1阶模式外,其余高阶模式均发生了很大变化。从原水层第1阶模式的截止频率开始,复合结构第2阶模式的相速度曲线被压低,各高阶(大于2阶)模式的相速度曲线出现一个跃变点,群速度曲线出现一个极大和一个极小值。水层越厚,复合结构各高阶模式的截止频率越低,相同频带内导波模式越丰富。水层厚度保持不变时,复合结构各阶模式的相速度和群速度曲线均随冰层厚度的增加而向低频方向移动。另外,还进一步分析了冰-水复合结构的导波波结构,发现第1阶导波模式的能量主要集中在冰层内和海表面附近,而2阶以上高阶导波模式的振动位移幅度随深度方向呈现周期性特征,并且模式阶数越高,振动越复杂。      

圆管超声导波频散与多模态特性

超声导波检测技术作为一种新兴的无损检测技术广泛应用于圆管类结构。为选择合适于不同缺陷检测的导波模态,推导分析了圆管导波传播的运动方程和频散方程;利用数值计算的方法得到了超声导波在圆管中传播的频散曲线和各模态沿壁厚方向的位移分布图,分析得出各个模态对不同缺陷的敏感程度;以一种特定的圆管为例,建立圆管缺陷有限元模型,对不同类型圆管缺陷对导波传播特性的影响进行仿真计算。结果表明,纵向模态对周向缺陷比较敏感,而扭转模态则对轴向缺陷更敏感,仿真结果与理论分析结果相吻合,为圆管缺陷检测的超声导波模态选择提供了理论依据。     

圆管超声导波频散与多模态特性研究*

超声导波检测技术作为一种新兴的无损检测技术广泛应用于圆管类结构。为选择合适于不同缺陷检测的导波模态,推导分析了圆管导波传播的运动方程和频散方程;利用数值计算的方法得到了超声导波在圆管中传播的频散曲线和各模态沿壁厚方向的位移分布图,分析得出各个模态对不同缺陷的敏感程度;以一种特定的圆管为例,建立圆管缺陷有限元模型,对不同类型圆管缺陷对导波传播特性的影响进行仿真计算。结果表明,纵向模态对周向缺陷比较敏感,而扭转模态则对轴向缺陷更敏感,仿真结果与理论分析结果相吻合,为圆管缺陷检测的超声导波模态选择提供了理论依据。     

管间界面特性对周向超声导波传播特性的影响

采用界面弹簧模型对圆管结构的管间界面特性进行描述, 推导出含弱界面的圆管结构中声波沿周向传播时的位移场及应力场的数学表达式. 在此基础上采用导波的模式展开分析方法, 给出了与管间界面特性及激励源密切相关的周向超声导波模式展开系数的解析表达式. 数值分析了管间界面特性的变化对周向超声导波的频散和声场产生的影响. 理论与数值分析结果表明, 通过选择适当的驱动频率及周向导波模式, 可使周向超声导波的相速度及圆管外表面的位移场随管间界面特性的变化表现出非常敏感且单调的性质. 这一结果有助于采用周向超声导波方法准确定征圆管结构的管间界面特性.     

液态金属快堆主管道中超声导波传播特性研究

讨论了钠冷快堆(Sodium-cooled Fast Reactor,SFR)主管道的整体温度和内部液态金属钠流动速度的变化对管道导波传播特性的影响。推导了充液管道中导波频散方程的一般形式,并给出了管道内液态金属钠处于流动状态下的导波频散方程。采用数值计算方法获得了管内液态金属钠处于不同温度和不同流速时的导波纵向模式频散曲线和导波时域波形。结果表明,温度变化对基阶纵向模式的影响较小,但对高阶纵向模式的影响较大;液态钠流速增大会使导波频散曲线向高频轻微移动,但在实际检测中可以忽路管内液体流动速度的影响。通过对时域接收波形的模拟计算,进一步考察了液态金属钠的温度及流动速度变化对导波传播的影响,并通过对比不同模态的激发特点和不同频段的导波时域波形特点,结合导波频散曲线,给出了适用于SFR管道超声无损检测的导波模态和声源激发频段选择方案。      

基于波传播法的周期复合板振动带隙衰减特性研究

运用波传播法对有限和无限周期对边简支复合板的振动带隙衰减特性进行了研究.在建立相邻板结构边界连续方程的基础上, 分别运用传递矩阵和Bloch定理建立了有限和无限周期复合板的耦合运动方程, 并详细对比分析了有限和无限周期复合板带隙衰减特性的关联关系.研究表明: 周期板结构的振动带隙频率范围与激励方式和激励位置是相关的, 若周期复合板在宽度方向按某阶模态进行线激励, 则该激励下的振动带隙与无限周期复合板在该阶模态下的振动带隙是一致的; 若周期板在点激励作用, 则该点激励下的振动带隙是参与振动的各阶模态振动带隙的交集. 此外, 还进一步研究了结构阻尼对振动衰减带隙的影响. 关键词： 周期复合板 带隙衰减特性 波传播法 结构阻尼     

四边简支含多孔材料双层板隔声特性

有限大含多孔材料和空气层复合板结构隔声特性的研究尚不充分。该文旨在研究四边简支边界条件下复合板结构隔声特性。首先基于流体饱和多孔弹性介质的声传播理论计算了声波在多孔介质中的传播波数；继而采用四边简支边界条件下板结构的模态函数，利用模态叠加法和伽辽金法推导了复合板结构隔声系数理论模型，并数值求解复合板隔声量。将理论模型得到的四边简支复合板隔声量与实验结果对比，验证了理论模型的正确性。最后，详细讨论了边界条件、板结构尺寸和多孔材料厚度等主要参数对隔声特性的影响。结果表明：四边简支复合结构隔声特性曲线上“谷值”出现得较少，并且简支复合板隔声特性的第一个“吻合频率”比固支支撑复合板更靠近低频，当频率超过3000 Hz以后，简支和固支边界条件复合板结构隔声特性趋于一致。     

SH波在中厚板对接焊缝中的传播特性

本文对应用于中厚板对接焊缝缺陷检测的SH波传播特性进行了理论和实验研究。基于半解析有限元法建立标准对接焊缝的有限元模型,获得SH波在结构中的能量分布和衰减特性,并探索不同对接焊缝结构对导波能量集中特性的影响,最后通过实验验证SH波的衰减情况和实际缺陷检测效果。结果表明,SH模态能量集中、衰减小、传播特性稳定,当焊缝结构特征加强时,模态主位移更为集中,实验得出利用SH波进行中厚板长距离对接焊缝的缺陷检测具有一定的研究价值和意义。本研究为提高中厚板对接焊缝缺陷检测效果提供理论依据,也为焊缝特征导波检测的实际应用奠定基础。     

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.     

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 propagating in sandwich structures made of anisotropic,viscoelastic, composite materials

The propagation of Lamb-like waves in sandwich plates made of anisotropic and viscoelastic material layers is studied. A semi-analytical model is described and used for predicting the dispersion curves (phase velocity, energy velocity, and complex wave-number) and the through-thickness distribution fields (displacement, stress, and energy flow). Guided modes propagating along a test-sandwich plate are shown to be quite different than classical Lamb modes, because this structure does not have the mirror symmetry, contrary to most of composite material plates. Moreover, the viscoelastic material properties imply complex roots of the dispersion equation to be found that lead to connections between some of the dispersion curves, meaning that some of the modes get coupled together. Gradual variation from zero to nominal values of the imaginary parts of the viscoelastic moduli shows that the mode coupling depends on the level of material viscoelasticity, except for one particular case where this phenomenon exists whether the medium is viscoelastic or not. The model is used to quantify the sensitivity of both the dispersion curves and the through-thickness mode shapes to the level of material viscoelasticity, and to physically explain the mode-coupling phenomenon. Finite element software is also used to confirm results obtained for the purely elastic structure. Finally, experiments are made using ultrasonic, air-coupled transducers for generating and detecting guided modes in the test-sandwich structure. The mode-coupling phenomenon is then confirmed, and the potential of the air-coupled system for developing single-sided, contactless, NDT applications of such structures is discussed.     

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.     

Natural beam focusing of non-axisymmetric guided waves in large-diameter pipes

The propagation of non-axisymmetric guided waves in larger diameter pipes is studied in this paper by treating the guided waves as corresponding Lamb waves in an unwrapped plate. This approximation leads to a simpler method for calculating the phase velocities of hollow cylinder guided waves, which reveals a beam focusing nature of non-axisymmetric guided waves generated by a partial source loading. The acoustic fields in a pipe generated by a partial-loading source includes axisymmetric longitudinal modes as well as non-axisymmetric flexural modes. The circumferential distribution of the total acoustic field, also referred as an angular profile, diverges circumferentially while guided waves propagate with dependence on such factors as mode, frequency, cylinder size, propagation distance, etc. Exact prediction of the angular profile of the total field can only be realized by numerical calculations. In particular cases, however, when the wall thickness is far less than the cylinder diameter and the wavelength is smaller than or comparable to the pipe wall thickness, the acoustic field can be analyzed based on the characteristics of Lamb waves that travel along a periodic unwrapped plate. Based on this assumption, a simplified model is derived to calculate the phase velocities of non-axisymmetric flexural mode guided waves. The model is then applied to discussions on some particular characteristics of guided-wave angular profiles generated by a source loading. Some features of flexural modes, such as cutoff frequency values are predicted with the simpler model. The relationship between the angular profiles and other factors such as frequency, propagation distance, and cylinder size is obtained and presented in simple equations. The angular profile rate of change with respect to propagation distance is investigated. In particular, our simplified model for non-axisymmetric guided waves predicts that the wave beam will converge to its original circumferential shape after the wave propagates for a certain distance. A concept of "natural focal point" is introduced and a simple equation is derived to compute the 1st natural focal distance of non-axisymmetric guided waves. The applicable range of the simplified equation is provided. Industrial pipes meet the requirement of wall thickness being far less than the pipe diameter. The approximate analytical algorithms presented in this paper provides a convenient method enabling quick acoustic field analysis on large-diameter industrial pipes for NDE applications.     

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.     

Effect of porous surface layer on wave propagation in elastic cylinder immersed in fluid

To study the damage to an elastic cylinder immersed in fluid, a model of an elastic cylinder wrapped with a porous medium immersed in fluid is designed. This structure can both identify the properties of guided waves in a more practical model and address the relationship between the cylinder damage degree and the surface and surrounding medium. The principal motivation is to perform a detailed quantitative analysis of the longitudinal mode and flexural mode in an elastic cylinder wrapped with a porous medium immersed in fluid. The frequency equations for the propagation of waves are derived each for a pervious surface and an impervious surface by employing Biot theory. The influences of the various parameters of the porous medium wrapping layer on the phase velocity and attenuation are discussed. The results show that the influences of porosity on the dispersion curves of guided waves are much more significant than those of thickness, whereas the phase velocity is independent of the static permeability. There is an apparent "mode switching" between the two low-order modes. The characteristics of attenuation are in good agreement with the results from the dispersion curves. This work can support future studies for optimizing the theory on detecting the damage to cylinder or pipeline.