声弹性理论下铝板中初始剪应力对超声导波传播特性的影响：模式耦合与色散曲线分裂机制研究

《Ultrasonics》：Effect of initial shear stress on acoustoelastic guided waves in plates

本文研究了一个在无损检测和结构健康监测中至关重要但此前鲜少被系统探讨的问题：初始剪应力如何影响板状波导中超声导波的传播。研究团队采用基于声弹性理论的半解析有限元方法，系统分析了铝板在不同平面内三种剪应力工况下的导波色散曲线、波结构与模式耦合行为。结果表明，与轴向应力相比，剪应力对导散曲线的影响较弱，且主要呈现与应力大小平方相关的二次性行为。其机理可通过将剪应力分解为拉压应力分量来理解。更重要的是，剪应力会引入不同导波模式（如Lamb波与SH波）之间的耦合，这种耦合在色散曲线原本相交的频率附近尤为显著，导致曲线发生分离和转向，表现出高敏感性。该研究不仅深化了对声弹性效应下导波特性的理解，还基于模式耦合的幅值比，提出了一种新颖的剪切应力测量新原理。

　　

在工程结构的安全评估与健康监测领域，超声导波扮演着“侦察兵”的重要角色。通过分析其在材料中传播速度、波形等特征的变化，可以推断出结构内部的“健康状况”，比如是否存在裂纹或应力集中。然而，这个“侦察兵”本身也会受到环境的影响——结构所承受的初始应力（例如拉伸、压缩或剪切力）会改变超声波的传播行为，这就是所谓的“声弹性效应”。如果忽视这种效应，就可能导致对监测信号的误判；反过来，如果善加利用，则可以通过超声波的变化来精确反推结构所受的应力大小，实现“一举两得”。

目前，学术界对于轴向（拉压）应力如何影响导波已经有了相当深入的研究，但对于另一种同样常见且重要的载荷形式——剪切应力，其对导波传播特性的系统性影响却鲜有问津。特别是在剪应力作用下，不同类型的导波模式之间会发生怎样的相互作用？这种相互作用能否被用来开发新的测量技术？这些问题都亟待解答。为此，来自巴西里约热内卢天主教大学电气工程系的研究团队Alan C. Kubrusly、Luis Paulo Brasil de Souza和Fotis Kopsaftopoulos在《Ultrasonics》期刊上发表了一项开创性研究，旨在全面揭示初始剪应力对板中超声导波传播的影响，重点关注其波结构和模式耦合行为随应力大小的变化规律。

为了深入探究这一问题，研究人员采用了一种强大的数值工具——基于声弹性理论的半解析有限元方法。该方法的核心思想是：在波导的横截面上进行有限元离散，同时假设波沿传播方向是谐波形式的解析解。这使得研究者能够有效求解在初始应力存在下的波动方程特征值问题，从而获得导波模式的色散关系（即波速随频率变化的曲线）以及波结构（即位移场在横截面上的分布）。研究针对一块铝板，设定了三种典型的剪应力加载情况：板水平面内的剪应力(s₁₂ⁱ)、平行于波传播方向的剪应力(s₁₃ⁱ)和垂直于波传播方向的剪应力(s₂₃ⁱ)。应力大小从0变化到150兆帕，以系统分析其影响。

研究取得了多项重要发现，具体体现在以下几个方面：

剪应力导致导波模式耦合与色散曲线分裂

研究最引人注目的发现之一是，剪应力会打破各向同性板中导波模式原有的独立性。在无应力状态下， Lamb波（其位移具有对称或反对称性）和SH波（水平剪切波）是相互解耦、独立传播的。然而，施加剪应力后，应力诱导的“表观各向异性”使得这些原本“井水不犯河水”的模式之间发生了耦合。这种耦合在色散曲线原本相交的频率点附近表现得尤为剧烈，导致曲线不再相交，而是发生“转向”或“分裂”——即两条曲线相互靠近后交换路径。例如，在s₁₂ⁱ剪应力下，S0模式与SH0模式在约3300 kHz·mm的频率-厚度积附近发生耦合与分裂；在s₂₃ⁱ剪应力下，S0模式与SH1模式在约2100 kHz·mm附近发生类似现象。这种分裂区域的波速对应力变化极为敏感，为高灵敏度应力测量提供了潜在窗口。

波结构演变揭示耦合物理机制

仅观察波速变化不足以理解耦合的本质。SAFE方法天然能提供每个导波模式的“波结构”，即位移场在板厚度方向上的分布。分析表明，在剪应力作用下，原本“纯净”的模式会“混入”其他模式的位移分量。例如，无应力时，S0模式只有纵向(U₁)和垂向(U₃)位移，水平横向(U₂)位移为零；SH0模式则只有U₂位移。但在s₁₂ⁱ剪应力下，S0模式出现了非零的U₂分量，而SH0模式则出现了原本属于S0模式的U₁和U₃分量。这种波结构的“杂交”直观证实了模式间的能量交换与耦合。

剪应力影响以二次行为为主导且影响较弱

总体而言，与轴向应力相比，剪应力引起的导波相速度变化要小得多（大约低一个数量级）。更重要的是，其变化规律并非与应力大小呈简单的线性关系。通过将相速度变化(Δc)分别除以应力(σ)和应力的平方(σ²)进行分析，研究人员发现Δc/σ²在不同应力水平下几乎完全重叠，而Δc/σ则不然。这清晰地表明，剪应力引起的波速变化主要呈现二次方依赖性。这一现象可以通过应力张量分解来理解：任何一种剪应力状态都可以分解为两个大小相等、方向相互垂直的拉应力和压应力的叠加。由于拉压应力对波速的影响近似线性且方向相反，它们的线性效应相互抵消，留下的主要就是与应力平方相关的非线性（二次）项贡献。

不同剪应力方向导致不同的耦合模式

三种剪应力工况的影响各有特点。s₁₂ⁱ（水平面内剪应力）主要引发Lamb波与SH波之间的耦合（如S0-SH0， S1-SH2等）。s₁₃ⁱ（沿传播方向的剪应力）则主要引发同类型但不同对称性的模式之间的耦合，例如反对称A2模式与对称S2模式之间的耦合，而对SH波影响甚微。s₂₃ⁱ（垂直于传播方向的剪应力）的影响与s₁₂ⁱ类似，也会引起Lamb波与SH波的耦合，但具体的耦合对象不同（如S0与SH1耦合）。

提出基于模式耦合比的剪切应力测量新原理

基于对模式耦合行为的深入理解，本研究提出了一种创新的应力测量思路。在模式强烈耦合的频率点，导波的波结构是多个位移分量的混合体。例如，在s₁₂ⁱ作用下，原本的S0模式中会“混入”SH波的特征位移分量(U₂)，反之亦然。研究人员指出，可以通过测量这些耦合位移分量的幅值比来反演剪应力的大小。由于该比值与应力大小直接相关，且可能在灵敏度上优于传统的基于波速变化的测量方法。论文明确给出了为实现这种测量原理的理想工作点：对于s₁₂ⁱ，建议频率-厚度积约为3300 kHz·mm；对于s₂₃ⁱ，建议约为2100 kHz·mm。

综上所述，这项研究首次系统、全面地揭示了初始剪应力对板中超声导波传播的多维度影响。它不仅证实了剪应力会导致显著的导波模式耦合和色散曲线分裂，阐明了其以二次方行为为主导的影响规律，并创新性地提出了利用耦合位移分量比值进行应力测量的新原理。该工作极大地深化了声弹性理论在复杂应力场下的内涵，为发展下一代基于超声导波的高精度、多参数结构应力监测与无损评估技术奠定了重要的理论基础，具有显著的学术价值与应用前景。