THERMOELASTIC WAVE PROPAGATION IN PLATES: AN IMPROVED LEGENDRE POLYNOMIAL APPROACH 1)

近年来, 超声导波因其衰减小, 传播距离远和信号覆盖范围广, 成为无损检测领域快速发展的方向之一. 然而, 基于超声导波的高温在线检测和激光超声技术却发展缓慢, 其关键在于热弹耦合波动方程求解难度大、传播与衰减特性研究困难. 作为一种有效的求解方法, 勒让德正交多项式方法已广泛应用于导波传播问题, 但该方法在求解热弹导波传播时存在两个不足, 限制其进一步的发展和应用. 这两个缺陷是: (1)求解过程中大量积分的存在, 致使计算效率低下; (2)仅能处理等热边界条件的热弹导波传播. 针对两项不足之处, 提出一种改进的勒让德正交多项式方法, 以求解分数阶热弹板中的导波传播. 推导求解方法中积分的解析表达式, 以提高计算效率; 引入温度梯度展开式, 发展适合勒让德多项式级数的绝热边界条件处理方法. 与已有文献结果对比表明改进方法的正确性; 与已有方法的计算时间对比说明改进方法的高效性. 最后将改进的方法用于求解分数阶热弹板中的导波传播, 研究分数阶次对频散、衰减曲线和应力、位移、温度分布等的影响.     

功能梯度圆柱板中的周向衰逝导波

与传播波模态不同,衰逝波模态波数为纯虚数或复数,它们随传播距离呈指数或阻尼式衰减。复数根描述的衰逝波对结构缺陷形状和尺寸的导波检测具有重要作用,但其求解却是非常困难的,通常要借助于各种迭代技术。本文提出一种计算求解衰逝波问题的改进的Legendre正交多项式方法,该方法可将复杂的变系数微分方程组计算转换为特征值求解问题,无需迭代便能计算得到包含实波数域、虚波数域和复波数域的完整的频散曲线。通过具体算例验证了提出方法的正确性。应用提出方法计算了不同梯度圆柱板中的衰逝导波,绘制了三维频散曲线,研究了不同径厚比和梯度场对衰逝导波频散的影响,分析了衰逝导波的位移和应力分布,讨论了衰逝导波的传播特性。     

随机结构有限元分析的递推求解方法的改进

将随机结构有限元分析的递推求解方法和伽辽金投影方法相结合,提出了求解随机静力响应的改进的递推求解方法。利用随机收敛的非正交多项式展开表示由于材料、外部荷载或构件几何尺寸的随机性导致的结构随机响应。采用递推求解方法得到响应多项式展开的初始系数,并运用定义的数学算子显式地表达出来。然后,通过定义修正系数,应用伽辽金方法对随机力平衡方程在非正交多项式基上进行投影,得到了和响应展开阶次个数相同的确定的有限元方程,并进行求解得到了修正系数。数值算例表明,通过对递推求解方法中响应表达式系数的修正,以很小的计算代价较大地提高了随机响应的计算精度;与基于正交多项式展开的随机有限元方法相比,在精度相当的前提下新方法耗费的计算时间大大降低。     

饱和多孔地基与矩形板动力相互作用的非轴对称混合边值问题

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.      

高频纵向导波在钢杆中传播特性的研究

计算了钢杆中纵向轴对称导波模态的衰减频散曲线和群 速度频散曲线. 分析了1~3MHz范围内高频纵向轴对称超声导波在钢杆中的传播特 性. 理论分析表明，各高阶纵向模态都存在一个衰减最小值，在此衰减最小值所对应频率下 的高阶模态能传播较远距离，可用于钢杆导波检测. 建立实验系统，采用轴对称同端激 励接收的方法，根据第1次端面回波做出群速度和端面回波幅值随频率变化曲线，实验结 果与理论分析基本吻合. 表明考虑材料衰减特性的钢杆频散曲线可以作为实验指导依据.     

加权残数配点法解正交各向异性板的积分方程

本文推导了一般各向异性板弯曲的积分方程,运用加权残数配点法求解了正交各向异性板弯曲的积发方程,本文将部分配点取在边界上,另一部分配点取在域外,只用关于找度的基本积分方程,而不用关于转角的补充积分方程,简化了方程求解和计算程序,由于正交各向异性板没有争析形式的、实用的基本解,本文提出了两种新的近似基本解;加权双三角级数;广义各向同性板解析形式的基本解和加权双三角级数的叠加,算例表明,本文提出的解法和近似基本解适用于各类边界条件的正交各向异性板,具有简单、可靠、精度高等优点。     

复杂介质中任意阶频率依赖耗散声波的分数阶导数模型

实验现象表明,声波在复杂介质中传播时,其衰减往往呈现频率的任意次幂律依赖现象.鉴于复杂介质的力学和物理性质的记忆性和长程相关性,频率幂律依赖的声波衰减现象难以用经典的声波方程描述,因为经典的阻尼波方程和近似热黏性波方程只能分别描述与频率无关和频率二次方依赖的声衰减.近年来,带有分数阶导数项的声波方程已被成功用于描述这一声衰减现象.基于课题组对声波衰减分数阶导数建模的研究,对已有的分数阶导数声波方程的研究进展及获得的成果做一个系统的综述,重点讨论这些模型的力学本构、统计力学解释等.简述了软物质中声波传播的时间分数阶导数唯象模型和本构模型,空间分数阶导数唯象模型和本构模型,并深入讨论了各种模型之间的联系与区别:介绍了分数阶导数声波模型在多孔介质中的成功应用,该部分内容涉及了均匀和非均匀多孔介质,刚性固体骨架和可变形固体骨架多孔介质等;通过空间分数阶扩散方程与Levy稳定分布之间的联系,给出了频率幂律依赖指数的变化区间为[0,2]的统计力学解释.最后,讨论了声波传播耗散行为的分数阶导数建模领域仍然存在的问题,并对今后的研究方向进行了探讨和展望.     

NON-AXISYMMETRIC MIXED BOUNDARY VALUE PROBLEM FOR DYNAMIC INTERACTION BETWEEN SATURATED POROUS FOUNDATION AND RECTANGULAR PLATE$^{\bf 1)}$

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.     

正交各向异性板中非主轴方向的Lamb波

基于线性三维弹性理论,采用勒让德正交多项式展开法,推导了波沿正交各向异性材料非主轴方向传播时的Lamb波耦合波动方程,并对耦合波动方程进行了数值求解。为验证该方法的适用性和正确性,首先将此方法应用于各向同性材料,并与已知的数据结果进行了比较;然后以单向纤维增强复合材料为例,计算了耦合Lamb波沿不同的非主轴方向传播时的相速度频散曲线,并分别研究了传播方向改变时低阶模态Lamb波和高阶模态Lamb波频散特性的变化。最后,针对潜在用于各向异性复合材料结构健康监测的耦合Lamb波低阶模态,给出了其在不同传播方向时的相速度分布和群速度分布。同时,结合低阶模态Lamb波的位移分布特性和材料的各向异性特点,阐释了S0模态对波的传播方向变化最为敏感的原因。     

热弹耦合圆板非线性振动的研究

对温度场中圆板的非线笥热弹耦合自由振动问题,由非线性振动方程、协调方程及热传导方程出发,动用伽辽金法求解,得出一个关于时间的非线笥常策分方程组。将热弹耦合与非热弹耦合情况进行对比,发现给定初始位移较小时,热弹耦合效应使板的固有频率相对与无热弹耦合情形提高;给定初始位移较大时,热弹耦合２使固有频率降低,该文不还比较了不同热弹参数和边界条件对热弹耦合效应的影响。     

GUIDED THERMOELASTIC WAVE PROPAGATION IN LAYERED PLATES WITHOUT ENERGY DISSIPATION

In this paper, the propagation of guided thermoelastic waves in laminated orthotropic 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 validity of the method is confirmed through a comparison. The dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Dispersion curves of the corresponding pure elastic plate are also shown to analyze the influence of the thermoelasticity on elastic modes. The displacement and temperature distributions are shown to discuss the differences between the elastic modes and thermal modes.     

Investigation of guided waves propagation in orthotropic viscoelastic carbon–epoxy plate by Legendre polynomial method

The effect of viscoelasticity on the guided waves propagation in viscoelastic plate has been investigated according to multi-aspect. To this purpose, an extension of the Legendre polynomial method is proposed to formulate the guided waves equation in orthotropic viscoelastic plate composed of carbon–epoxy. The validity of the proposed Legendre polynomial method is illustrated by comparison with available data. The convergence of the method is discussed through a numerical example. The hysteretic and Kelvin–Voigt viscoelastic models are used to integrate the imaginary part of the complex stiffness matrix associated with the viscoelastic plate in this study. Accordingly, both viscoelastic models do not affect on the dispersion curves results. However, appreciable effects are seen in the attenuation curves. Also, the sensitivity of the guided waves propagation caused by variations of elastic and viscoelastic modulus has been studied in detail. Finally, the advantages of the Legendre polynomial method are described.     

Characteristics of guided waves in graded spherical curved plates

An elastodynamic solution for the stress wave propagation in spherical curved plates composed of functionally graded materials (FGM) is presented. Properties of materials are assumed to vary in the direction of the thickness according to a known radial variation law (gradient field). The formulation is based on the linear three-dimensional elasticity. The Legendre orthogonal polynomial series expansion approach is used for determining the guided waves dispersion curves in functionally graded spherical curved plates. Our results from a homogeneous anisotropic spherical curved plate are compared with those published earlier to confirm the accuracy and range of applicability of this polynomial approach. Guided wave dispersion curves for graded spherical curved plates with different gradient fields are calculated, and the effects of the gradient field on the characteristics of guided waves are illustrated. Finally, dispersion curves for graded spherical curved plates at different ratios of inner radius to thickness are calculated to discover the influences of that ratio on the wave characteristics.     

Transient thermoelastic waves in an anisotropic hollow cylinder due to localized heating

This paper presents a theoretical study of transient ultrasonic guided waves generated by concentrated heating of the outer surface of an infinite anisotropic hollow circular cylinder. Generalized thermoelastic theory proposed by Lord and Shulman is adopted to model the dynamic thermoelastic behavior of the cylinder. The concentrated heat source model used is to represent heating due to a pulsed laser beam, which is focused on the outer surface of the cylinder. A semi-analytical finite element (SAFE) method is employed to evaluate guided wave modes in the cylinder. Using integral transform techniques, the modal wave forms are obtained in frequency and wave number domains. Time histories of the propagating modes are then calculated by applying inverse Fourier transformation in the time domain. Numerical results showing the dispersion curves for the group velocities of the propagating modes and transient radial displacements are presented. For this purpose it is assumed that the cylinder is made of transversely isotropic silicon nitride (Si₃N₄). Attention is focused on the propagation characteristics of longitudinal and flexural modes separately.     

接地网圆钢杆中高频纵向超声导波无损检测方法研究

针对接地网圆钢杆腐蚀检测问题,本文进行了接地网圆钢杆高频纵向超声导波无损检测方法研究。首先对埋地环境下圆钢杆中纵向导波传播特性进行理论研究,分析了不同模态导波的群速度和衰减频散特性。研究发现,高阶纵向模态导波在衰减最小和群速度最大对应的频率处,在圆钢杆中传播能力强,是适合进行地埋圆钢杆导波检测试验的频率范围。在此基础上,进行了埋地圆钢杆高频纵向超声导波无损检测试验研究。结果发现,利用优选的检测参数可以很好地实现埋地圆钢杆中腐蚀缺陷检测。研究工作为接地网运行状态评价提供了很好的技术支撑。     

有限单元法在超声导波检测技术中的应用

超声导波检测技术具有对波导结构中的缺陷进行远距离无损检测的能力,多年来一直是无损检测领域关注的热点之一.有限单元法具有对各种复杂动力学问题进行计算的能力,已成为超声导波检测技术研究的重要工具.本文结合超声导波检测技术研究领域中的热点问题,对相关的有限单元法进行了简要综述.介绍了有限单元法的发展及其在多物理场耦合机制下导波的激励与接收、线弹性和黏弹性结构中导波的传播特性、非线性超声导波等多个方面的应用研究情况. 最后,基于超声导波检测技术研究趋势展望了相关有限单元法的未来研究重点和发展方向.      

一种基于深度学习的超声导波缺陷重构方法

超声导波检测因其传播效率高、耗能少等优势成为了无损检测领域的重要研究方向。目前已有的利用超声导波进行结构缺陷探测和定量化重构的方法主要由相关的导波散射理论推导得出。然而,由于导波散射问题本身的高复杂性,使得在推导上述理论方法时引入一些近似假设,降低了重构结果的质量。另外,有些方法通过优化迭代的方式提高重构精度,又会增加检测的时间成本。有鉴于此,本文探索了一种将卷积神经网络与导波散射理论模型以局部融合的方式实现缺陷定量化重构的新方法。应用样本数据训练后的神经网络实现缺陷定量化重构,弥补缺陷重构过程中的理论模型误差,同时去除在实际检测过程中所存在的环境噪声。本文以利用SH导波重构平板中的减薄缺陷为研究对象,通过数值模拟验证了该方法在缺陷重构时具有高效率和高精度的特点,特别是对矩形缺陷的重构,新方法的结果精度比波数空间域变换法的精度提高了近200%。     

Elasto-thermodiffusive (ETNP) surface waves in semiconductor materials

The present article is devoted to investigate the propagation of elasto-thermodiffusive (ETNP) surface waves in a homogeneous isotropic, thermally conducting semiconductor material of half-space with relaxation of heat and charge carrier fields. The secular equation, a more general functional relation, that governs the propagation of elasto-thermodiffusive (ETNP) surface waves in homogeneous isotropic, thermoelastic semiconductor material halfspace with relaxation of heat and charge carrier fields has been derived by solving a system of coupled partial differential equations. A hybrid numerical technique consisting of Descartes algorithm for solving complex polynomial characteristic equation along with functional iteration scheme has been successfully used to solve the secular equation in order to obtain dispersion curves, attenuation coefficient and specific loss factor of energy dissipation for p-type germanium (Ge) semiconductor. Some particular forms of the general secular equation governing the propagation of elasto-thermodiffusive (ETN/ETP), thermoelastic (ET), elastodiffusive (EP/EN) and thermodiffusive (TP/TN) surface waves have been also deduced and discussed. In order to illustrate the analytical development, the numerical solution of the secular equation and other relevant relations under different situations is also carried out for Ge semiconductor materials to characterize the elasto-thermodiffusive (ETP) and thermodiffusive (TP) surface waves. The computer simulated results have been presented graphically in respect of the dispersion curves, attenuation coefficient and specific loss factor.     

Analysis of damped guided waves using the method of multiple scales

We analytically investigate the influence of damping on Lamb waves, which are a specific type of guided wave in two-dimensional plates. Considering material attenuation, we suppose that Lamé constants are complex numbers. This leads to the associated wavenumbers being complex, with the imaginary part of the wavenumber being associated with effect of attenuation of the guided wave. In this paper, we show how dispersion curves and attenuation coefficients can be obtained using the self-adjointness and the method of multiple scales (MMS), which is a type of perturbation method. Using the self-adjointness and the MMS, we can calculate the frequency- and wavenumber-dependent attenuation coefficients from the integral values and boundary values of a corresponding eigenfunction with respect to each propagation mode. This analytical method can yield not only dispersion curves but also mode-by-mode attenuation coefficients regardless of the numerical initial values, unlike numerical approaches using the Newton method. Thus, the proposed method can more easily calculate the attenuation coefficients with respect to a particular mode than conventional methods. Furthermore, the results obtained by proposed method were in good agreement with those obtained by the semi-analytical finite element (SAFE) method, which validates the proposed method.