基于模态追踪算法的波导泄漏波频散曲线高效计算方法研究

《Computer Methods in Applied Mechanics and Engineering》：Computing leaky waves in semi-analytical waveguide models by exponential residual relaxation

本文针对多层波导与无界域耦合的复杂非线性特征值问题，提出一种基于常微分方程（ODE）的模态追踪算法。该方法通过将频散关系转化为初值问题，有效解决了传统多参数特征值方法计算量大的瓶颈，在保证精度的同时显著提升计算效率，为超声导波无损检测中的模态识别提供了新工具。

　　

在超声无损检测和结构健康监测领域，导波传播特性的精确建模始终是核心挑战。当弹性波在板状结构中传播时，会与周围介质发生能量交换，形成具有复数波数的泄漏模式（leaky modes）。传统多参数特征值方法虽能精确求解此类问题，但其计算成本随模型复杂度急剧增长，特别是对于多层各向异性材料或复杂耦合边界条件的情况，计算效率成为实际应用的瓶颈。

为解决这一难题，研究团队在《Computer Methods in Applied Mechanics and Engineering》发表论文，提出了一种创新的模态追踪算法。该方法的核心突破在于将非线性特征值问题重新表述为常微分方程初值问题：通过引入衰减参数χ，将残差最小化转化为动力学系统，利用ODE数值解法（如MATLAB的ode15s）沿频率轴连续追踪模态演化。关键技术环节包括：1）采用高精度谱有限元离散波导截面控制方程；2）构建包含辐射边界条件的耦合矩阵系统；3）设计自适应步长控制策略处理特征曲线奇点；4）开发双版本算法（Version i/ii）分别处理可微与非可微区间。

数值验证与性能分析

通过五个典型算例系统验证算法有效性：算例I（单层黄铜板水浸）显示算法能准确捕捉准Scholte模的急剧衰减转折；算例II（固-特氟龙半空间耦合）验证了高频区模态分离的追踪能力；算例III（钛-黄铜-钛三层结构油耦合）成功复现多模态耦合效应；算例IV（钛层与固相半空间强耦合）证明算法对高衰减模态（η>30,000 dB/m）的稳定性；算例V（T800-913复合材料[0/90]_2s铺层）则展现了各向异性材料的复杂频散特性重构能力。与多参数特征值法对比，新方法在保持精度的同时将计算时间降低1-2个数量级（如算例IV从288秒降至0.5秒）。

算法鲁棒性优化

研究特别解决了特征曲线交叉区域的数值稳定性问题：通过引入正则化参数δ=10^-12（对应材料阻尼），平滑处理特征根导数不连续点；采用均值近似法（mean-value approximation）估算初值，结合衰减参数自适应调整（χ₁=100, χ₂=10），确保在特征曲线快速变化区域仍保持L²误差低于10^-3。对于难以追踪的准Scholte模，额外引入近似解（k≈ω/c_j^f）作为补充初始条件。

本研究发展的模态追踪框架突破了传统频散曲线计算方法的效率限制，其核心价值体现在三个维度：计算效率层面，ODE转化策略将非线性求解转化为序列化计算，显著降低内存需求；应用广度层面，算法兼容各向异性材料、多层结构和复杂边界条件，可直接应用于复合材料检测、地下勘探等实际场景；方法论层面，提出的初值构造策略和正则化技术为同类非线性特征值问题提供了普适性解决思路。该工作为高频超声导波定量化应用奠定了坚实基础，未来可进一步扩展至热-力-电多场耦合波导建模等前沿领域。