综述：一种用于近似功能梯度板中兰姆波传播的谱配置法

《Scientific African》：A Spectral Collocation Method for Approximating Lamb Wave Propagation in Functionally Graded Plates

　　

研究人员开发了一种谱配置法（SCM）框架，用于功能梯度材料（FGM）板中的兰姆波传播。连续的厚度方向梯度由均匀层表示，其材料常数从Voigt混合定律中采样，而每层内的Chebyshev配置将控制方程简化为稀疏广义特征值问题。该方法通过与二分法、色散计算器（DC）软件以及已发表的各向同性和梯度板FGM基准进行验证。收敛性和性能研究表明，层细化比增加配置点数量更能控制精度：对于验证案例，层数NL=20就足够了（与DC相比最大相对误差低于0.62%），而NL=30-50在参数研究中提供了稳健的精度-成本折衷。通过使用MATLAB的移位-反演eigs求解器代替直接QZ（广义Schur）分解，计算保持了稀疏性，并实现了超过12,000倍的加速，同时均方根误差（RMSE）因子接近1。该框架为无损检测（NDT）和结构健康监测（SHM）应用中FGM板的快速色散曲线评估提供了实用指导。

**理论框架**
**问题描述**：研究人员研究了功能梯度材料（FGM）板中兰姆波的传播，板厚为h，由连续陶瓷-铬（Cer-Cr）混合物组成。顶面为纯陶瓷，底面为纯铬，两面均为自由表面。坐标系中x₁为厚度方向，x₃为波传播方向。
**功能梯度材料**：厚度方向材料变化遵循Voigt混合定律，材料属性（密度ρ和弹性刚度常数C₁₁、C₁₂、C₄₄）由陶瓷和铬的体积分数加权平均得到。梯度指数n∈[0,∞)控制材料组成：n→0为陶瓷主导，n→∞为铬主导，研究中考虑n=0.1至10以及均匀极限n=∞。
**波动方程控制方程**：FGM板被分解为多个均匀层，每层厚度Δh。在每层l中，兰姆波传播由无体力的弹性动力学方程控制，简化为(x₁, x₃)平面内的二维公式。对于沿x₃方向传播的兰姆波，位移分量假设为简谐形式，包含波数k和角频率ω。
**SCM的数学公式**：谱配置法（SCM）用于求解波动方程，将位移分量表达为特征值问题。每层厚度离散为N个Chebyshev配置点，配置点坐标在[-1,1]上定义，并映射到层厚度坐标。通过Chebyshev近似，得到每层的刚度子块矩阵，并组装成全局矩阵。对矩阵进行归一化处理以保证数值稳定性。上下表面施加无牵引条件，通过牵引算子施加边界条件。层间假设完美粘接，位移和应力连续，通过替换全局矩阵的相应行实施连续性条件，最终形成修正特征值问题。
**结果与讨论**
**各向同性板验证**：使用SCM计算1 mm厚各向同性铬板的相位速度色散曲线，与二分法和色散计算器（DC）软件结果高度一致。12个兰姆波模式的最大相对误差为0.62%（A1模式），平均值低于0.26%，验证了方法精度。
**FGM板中的波传播：方法验证**：将SCM应用于FGM板（陶瓷-Cr混合物），线性梯度情况（n=1）与文献参考数据[28]对比，相位速度色散曲线在无量纲波数（kh）域内几乎完全重合，验证了多层公式的有效性。
**收敛分析与计算性能**：影响SCM精度的参数为配置点数量N和层数NL。固定N=5，增加NL从5到100，均方根误差（RMSE）下降，NL≥20时误差低于0.6%。增加NL比增加N更有效地提高精度。对于计算性能，比较了直接求解器（eig，QZ算法）和迭代求解器（eigs，移位-反演Arnoldi/Lanczos）。迭代求解器保持稀疏性，计算时间随NL增长近似线性，而直接求解器时间随问题规模立方增长。在NL=100、N=20时，迭代求解器加速比超过12,000倍，且RMSE因子接近1，即精度相同。
**相位速度对材料梯度的敏感性**：梯度指数n对相位速度有显著影响。对于S0和A0模式，n增加使相位速度系统降低，尤其S0模式敏感。高频率时，色散曲线趋于渐近线，受局部材料属性主导。高阶模式（A1、S1）的截止频率随n增加而降低，改变模态频谱。
**群速度对梯度分布的依赖性**：群速度方面，A0模式在不同梯度下相对稳定，S0模式则表现出显著敏感性，最小群速度（MGV）随n增加而降低并向低频移动。A1和S1模式也呈现类似差异，可用于无损检测中的材料表征。
**梯度对波数和截止频率的影响**：梯度增加导致一阶和二阶兰姆模式（A1、S1、A2、S2）的截止频率系统降低，超过30%相对于均匀情况（n=0）。这丰富了模态内容，便于在较低频率利用高阶模式进行检测。
**厚度方向位移和应力场**：在1 MHz频率下，分析A0和S0模式在不同梯度下的位移（U₁、U₃）和应力（σ₁₁、σ₁₃）场。A0模式场分布相对稳定，表面能量集中，适合表面缺陷检测。S0模式场分布对梯度敏感，位移和应力模式随n变化显著，可用于材料梯度剖面重建。两者互补：S0用于梯度表征，A0用于缺陷检测。