Nonlinear Incompressible Shear Wave Models in Hyperelasticity and Viscoelasticity Frameworks, with Applications to Love Waves

本文提出了适用于任意应变能函数的不可压缩超弹性与粘弹性材料中的非线性剪切波方程，将其应用于 Love 波传播分析，并通过数值模拟揭示了界面波与表面波在完全非线性情形下的传播特性及渐近行为。

要点

这篇文章就像是在给地球和橡胶做了一次深度的"CT 扫描”，只不过这次我们不看骨头，而是看地震波和剪切波在复杂的材料中是如何“跳舞”的。

想象一下，你手里有一块橡皮泥（或者一块果冻），当你用手指快速划过它的表面时，会产生一种波浪。在传统的物理课本里，我们通常假设这块橡皮泥是“听话”的：你推它一下，它就动一下，力越大动得越快，而且这种关系是直线的（线性）。

但现实世界中的材料（比如橡胶、生物组织、甚至地壳深处的岩石）并不这么“听话”。当你用力过猛时，它们会表现出非线性：推得越狠，反应可能越奇怪，甚至会出现像水波破碎那样的现象。

这篇文章就是为了解决这个问题：当这些材料发生大变形，并且具有粘弹性（像蜂蜜一样有粘性）时，那些沿着界面传播的“剪切波”（Love 波）到底是怎么跑的？

以下是用通俗语言和比喻对文章核心内容的解读：

1. 核心角色：Love 波（拉夫波）

想象一下，地球表面覆盖着一层较软的“皮”（比如地壳），下面是一层较硬的“肉”（地幔）。

• 线性世界（旧理论）： 如果这层皮和肉是完美的弹簧，那么一种叫"Love 波”的波浪只能在特定的速度范围内存在。就像你只能在特定的速度下骑滑板才不会摔跟头。如果速度太快或太慢，波就传不过去，或者会散开。
• 非线性世界（新发现）： 作者发现，当材料被剧烈挤压或拉伸（大变形）时，这个“速度范围”的规则依然有效，但波的行为变得更复杂、更有趣。

2. 数学模型：从“直线”到“曲线”

• 旧模型（胡克定律）： 就像拉一根完美的弹簧，拉得越长，弹力越大，比例是固定的。
• 新模型（超弹性 + 粘弹性）： 作者建立了一个更复杂的数学公式。
  • 超弹性（Hyperelasticity）： 想象拉橡皮筋，刚开始很容易，拉到最后变得非常硬。这个模型能描述这种“越拉越硬”的非线性特性。
  • 粘弹性（Viscoelasticity）： 想象拉一块口香糖。它不仅会弹回来，还会因为内部的摩擦（粘性）而消耗能量，慢慢停下来。
  • 结果： 他们推导出的方程里充满了“立方”和“五次方”项。简单说，就是波的速度不再仅仅取决于材料本身，还取决于波有多强。波越强，跑得越快或越慢，甚至波形会扭曲。

3. 计算机模拟：在数字世界里“爆炸”

为了验证这些复杂的公式，作者没有去挖地球，而是在电脑里建了一个虚拟的实验室。

• 实验设置： 他们模拟了一个两层结构（上层软，下层硬），然后在某个点制造了一个“高斯爆炸”（想象在果冻里突然点了一下，产生一个圆形的扰动）。
• 观察现象：
  • 波的分裂： 波在传播过程中，一部分沿着界面跑，一部分沿着表面跑。
  • 速度的归宿： 无论一开始波跑得有多乱，随着时间推移，界面波和表面波的速度最终都会趋向于两层材料中较快的那一层的速度。就像两辆车在高速公路上，慢车最终会被快车甩开，或者被快车的尾流带动，最终大家都以快车的速度跑。
  • 粘性的作用： 当加入“粘性”（像蜂蜜一样的阻力）后，波的能量会迅速消散，波形变得更平滑，不再那么剧烈地震荡。

4. 对称性与精确解：寻找“完美舞步”

在文章的后半部分，作者尝试用一种叫“李群对称性”的高级数学工具，去寻找这些复杂方程的“完美舞步”（精确解）。

• 虽然他们找到了一些数学上完美的解，但这些解在物理上往往是“无限大”的（比如波高无限高），这在现实中是不可能的。
• 不过，通过计算机模拟，他们发现真实的波虽然不像数学解那样完美，但形状和走势与这些理论解非常相似。这就像虽然你不能画出完美的圆，但你的手绘圆看起来确实像个圆。

5. 为什么这很重要？（现实意义）

• 地震预警： 地震波在地球内部传播时，往往伴随着巨大的能量和变形。理解非线性效应，能让我们更准确地预测地震波到达地面的时间和强度，从而更好地保护建筑物。
• 医学与材料： 这种模型不仅适用于地球，也适用于人体组织（如肌肉、肌腱）或新型橡胶材料。医生可以用它来通过超声波更精准地诊断组织病变，工程师可以设计更抗震的材料。

总结

这就好比以前我们以为波浪在沙滩上跑是像火车一样按固定轨道走（线性）；现在作者告诉我们，波浪其实更像是一群在拥挤人群中奔跑的人（非线性），他们互相推挤、速度会变、还会因为摩擦停下来。

这篇文章通过建立更聪明的数学公式和强大的电脑模拟，揭示了这些“调皮”的波浪在复杂材料中的真实舞步。虽然公式看起来很吓人（充满了立方和五次方），但它们最终告诉我们：无论材料多复杂，波最终还是会找到它该走的路，只是过程比我们要想的更精彩、更曲折。

技术摘要

这是一份关于《不可压缩超弹性和粘弹性框架下的非线性剪切波模型及其在勒夫波（Love Waves）中的应用》论文的详细技术总结。

1. 研究问题 (Problem)

该研究旨在解决在有限变形（Finite Deformations）条件下，不可压缩超弹性及超粘弹性材料中剪切波传播的数学建模问题。具体关注点包括：

• 非线性效应：传统的线性弹性理论（胡克定律）仅适用于小变形，无法准确描述聚合物、橡胶、生物组织及地质层在大变形下的行为。地震波（特别是表面波）在震源附近或大振幅传播时表现出显著的非线性特征。
• 界面波传播：研究两种具有不同机械性质材料的界面（如地壳与地幔的莫霍面）上非线性勒夫波（Love waves，即水平偏振的剪切波）的传播特性。
• 耗散机制：在长时间尺度下，材料（如地质构造、生物组织）表现出粘弹性，需要考虑能量耗散对波形的影响。
• 现有模型的局限：现有的线性勒夫波模型假设相位速度恒定且与深度无关，但在大振幅和非线性介质中，波速会随时间和空间变化，且线性模型无法捕捉波破碎（wave breaking）或色散平衡等复杂现象。

2. 方法论 (Methodology)

论文采用连续介质力学框架，结合理论推导与数值模拟，主要步骤如下：

• 理论推导框架：

  • 超弹性模型：基于拉格朗日（材料）坐标系，利用应变能密度函数 $W^H$ 描述材料本构关系。针对不可压缩材料，推导了运动方程。
  • 本构模型选择：重点研究了广义 Neo-Hookean 材料（$W^H$ 仅依赖于第一不变量 $I_1$），并具体应用了三次 Yeoh 模型（Cubic Yeoh model），其应变能密度包含 $I_1$ 的线性和三次项，能够捕捉非线性应力 - 应变曲线的立方项和五次项效应。
  • 超粘弹性扩展：引入伪应变能密度函数 $W^V$（依赖于变形率张量 $\dot{C}$），构建了超粘弹性本构模型，以模拟能量耗散。
  • 控制方程推导：推导了描述 $Y$ 方向位移 $v$ 在 $(X, Z)$ 平面传播的非线性波动方程。对于三次 Yeoh 模型，方程包含立方和五次微分多项式项。

• 数值模拟：

  • 方法：采用直线法（Method of Lines），将偏微分方程（PDE）在空间上离散化为常微分方程组（ODE），然后使用 MATLAB 的 ode23 求解器进行时间积分。
  • 网格策略：使用非均匀网格，在波源附近和界面处加密网格点，以捕捉大振幅和非线性效应，同时节省计算资源。
  • 初始条件：模拟了“高斯爆炸”（Gaussian explosion）作为初始扰动，研究波在双层介质（上层覆盖层和下层半空间）中的传播。
  • 对称性分析：针对一维简化模型，利用Lie 群对称性方法推导精确的不变解（Invariant solutions），并与数值解进行对比验证。

• 边界条件：

  • 上表面为自由表面（Neumann 边界条件，应力为零）。
  • 界面处位移和应力连续。
  • 远场边界设为刚性壁（在足够大的计算域内模拟无限远）。

3. 关键贡献 (Key Contributions)

1. 广义非线性剪切波方程的推导：

   • 提出了适用于任意应变能密度函数的不可压缩超弹性材料的一般剪切波方程。
   • 针对三次 Yeoh 模型，推导出了包含立方和五次非线性项的具体控制方程，并明确了粘性项（混合导数 $u_{xxt}$ 等）的形式。
   • 澄清了广义 Neo-Hookean 材料假设的必要性，以避免在 $W^H$ 依赖第二不变量 $I_2$ 时出现的静水压力定义不唯一的问题（通过附录中的勘误说明）。

2. 非线性勒夫波的存在性与传播特性：

   • 通过数值模拟证明，即使在完全非线性情况下，界面波和表面波的变量波速 $v$ 在长时间演化后，仍然满足线性勒夫波的存在条件：$c_1 < |v| < c_2$（其中 $c_1, c_2$ 分别为上下层材料的剪切波速）。
   • 发现波速随时间演化，最终趋向于两层材料中较大的波速（$\max(c_1, c_2)$）。

3. 粘弹性效应的量化：

   • 展示了粘性项如何显著改变波的传播行为，导致波包快速色散和振幅衰减，防止了纯弹性模型中可能出现的数值波破碎现象。
   • 揭示了粘性系数 $\eta$ 对波形的平滑作用：随着 $\eta$ 增加，波形变得更加可预测和稳定。

4. 精确解与数值解的对比：

   • 利用 Lie 对称性分析获得了一维模型的一族精确解（尽管部分解无界）。
   • 数值模拟显示，有限振幅的波在演化过程中，其包络形状与有界版本的精确解高度吻合，验证了模型的物理合理性。

4. 主要结果 (Results)

• 线性 vs. 非线性：
  • 在 $c_1 < c_2$ 的情况下，非线性模型产生的波形与线性模型定性相似，但振幅衰减更快（由于非线性色散和粘性）。
  • 在 $c_2 < c_1$（违反线性勒夫波存在条件）的情况下，线性模型显示能量主要被困在下层，上层波迅速衰减；非线性模型表现出类似的能量分布特征，但波速演化过程更为复杂。
• 波速演化：
  • 无论是界面波还是表面波，初始时刻波速可能接近较小值或为零（取决于波源位置），但随着时间推移，反射波和折射波相互作用，波速最终稳定在 $\max(c_1, c_2)$。
  • 数值结果证实了非线性勒夫波的速度范围依然受限于线性理论给出的界限。
• 粘性影响：
  • 引入粘性后，波的能量耗散明显，波峰不再像纯弹性模型那样保持尖锐，而是变得平滑。
  • 在 $c_2 < c_1$ 且波源位于下层时，粘性导致位移 $v$ 在更多区域呈现负值，表明波的反射和干涉模式发生了改变。
• 一维模型行为：
  • 一维数值解表现出稳定的行波特征，振幅随粘性衰减。
  • 附录中发现，在特定径向对称条件下，数值解可能出现“波破碎”（梯度 $u_r$ 出现不连续），这通常发生在高振幅且缺乏足够耗散的情况下。

5. 意义与影响 (Significance)

• 地震学与地质勘探：该模型为理解地震表面波（特别是勒夫波）在复杂地质结构（如地壳 - 地幔界面）中的传播提供了更精确的理论工具。它解释了为何在大振幅或近震源区域，线性模型可能失效，并提供了修正方案。
• 生物力学与医学：模型适用于软组织和生物材料的剪切波成像（如弹性成像），有助于更准确地反演组织的非线性粘弹性参数。
• 材料科学：为橡胶、聚合物等超弹性材料在大变形下的动力学行为提供了新的数学描述，特别是引入了五次非线性项和粘性耦合项。
• 理论价值：
  • 证明了非线性波动方程在特定本构关系下，其波速仍受线性理论界限的约束，这是一个重要的理论发现。
  • 展示了 Lie 对称性方法在处理非线性粘弹性波动方程中的有效性，为寻找精确解提供了新途径。
  • 通过数值模拟揭示了非线性与粘性在平衡色散和防止波破碎方面的竞争机制。

综上所述，该论文成功地将超弹性和超粘弹性理论应用于非线性勒夫波的研究，建立了包含高阶非线性项和粘性耗散项的数学模型，并通过数值和解析手段验证了其在描述复杂波传播现象中的有效性，填补了从线性理论到完全非线性有限变形理论之间的空白。