Enabling topography-resolving structural dynamic contact simulation

这是一份关于论文《Enabling topography-resolving structural dynamic contact simulation》（实现拓扑解析的结构动态接触仿真）的详细技术总结。

1. 研究背景与问题 (Problem)

• 核心挑战：机械结构（如螺栓连接、铆接）的阻尼主要来源于连接处的摩擦耗散。准确预测这种阻尼是工程界的长期难题。
• 表面形貌的影响：实际制造表面的微观形貌（粗糙度、波纹度）对接触界面的初始压力分布、间隙以及动态载荷下的摩擦 - 单向接触行为有决定性影响。忽略这些形貌会导致预测偏差。
• 现有方法的局限性：
  • 计算成本：要解析真实的表面形貌，需要极细的网格（微米级），而结构整体尺寸较大（米级）。传统的有限元（FE）全模型仿真在处理这种多尺度问题时，计算量巨大，甚至不可行，尤其是对于长时程瞬态动力学仿真。
  • 数值误差：在细网格上应用集合值接触律（如库仑 - 西格诺里尼条件）时，传统 FE 方法容易出现数值不收敛、数值阻尼（Numerical Damping）和虚假高频振荡。数值阻尼会严重干扰对物理摩擦阻尼的准确评估。
  • 准静态方法的不足：现有的准静态模态分析（QSMA）虽然有效，但假设模态振型不变且法向载荷恒定，无法处理强非线性、振型显著变化或振动引起的“沉降”（Settling）现象。

2. 方法论 (Methodology)

作者提出了一种扩展的多尺度方法（Multi-Scale Method, MSM），将之前仅用于准静态分析的方法扩展到了动态分析领域。

• 混合建模策略 (FE-BE Coupling)：

  • 接触区域：使用基于**半空间理论（Half-space theory）的边界元（BE）**模型。在精细的规则网格上求解接触应力，接触柔度矩阵可以解析表达。
  • 主体结构：使用相对**粗网格的有限元（FE）**模型描述结构的整体动力学。
  • 耦合：通过满足远场的位移协调和力平衡条件将两者耦合。

• 动力学扩展的关键技术：

  • 无质量边界模态综合（Massless-boundary Component Mode Synthesis）：采用改进的 Craig-Bampton 方法，将 FE 模型缩减。关键在于移除接触边界节点的质量（生成奇异质量矩阵）。这消除了接触闭合时的有限冲量，从而避免了由接触引起的虚假高频振荡，使得无需引入人为的数值阻尼即可稳定仿真。
  • 半显式时间积分方案（Semi-explicit Time Integration）：
    • 采用蛙跳格式（Leapfrog scheme）。
    • 显式步：对内部动力学（模态坐标）进行显式积分，计算成本低。
    • 隐式步：对边界接触问题（准静态接触平衡）进行隐式求解，利用增广拉格朗日法（Augmented Lagrangian）处理接触约束。
    • 优势：允许使用较大的时间步长（仅受限于最高保留模态频率，而非网格尺寸），且时间步长与网格无关。
  • 谐波平衡法（Harmonic Balance, HB）：作为对比验证工具，直接计算周期稳态响应，用于验证时间积分方案是否存在数值阻尼。

3. 关键贡献 (Key Contributions)

1. 方法创新：首次将基于半空间理论的 FE-BE 多尺度方法成功扩展到动态时域积分和谐波平衡分析中，解决了拓扑解析动态接触仿真的计算瓶颈。
2. 数值稳定性：通过“无质量边界”模态综合和半显式积分方案，实现了无数值阻尼的仿真。这使得在长时程仿真中能够准确捕捉微小的物理摩擦阻尼，避免了传统方法中数值误差掩盖物理现象的问题。
3. 计算效率：相比全有限元（Full-FE）分析，该方法将计算效率提高了2-3 个数量级（约 400 倍于隐式 Full-FE，80 倍于显式 Full-FE），且允许使用更大的时间步长。
4. 新现象发现：利用该方法，首次在数值上揭示了**振动诱导的沉降（Vibration-induced Settling）**现象，即结构在经历冲击后，会收敛到一个与预紧状态不同的静态平衡位置，并伴随独特的残余接触应力场。

4. 研究结果 (Results)

• 基准验证 (S4 Beam)：使用著名的 S4 梁（具有实测表面形貌数据）作为基准系统。
  • 一致性：多尺度方法（MSM）的结果与全有限元（Full-FE）分析在频率和阻尼比上高度一致。
  • 数值阻尼验证：通过对比谐波平衡法（无数值阻尼）和时间积分法的结果，证实了 MSM 的时间积分方案没有明显的数值阻尼。相比之下，显式 Full-FE 分析显示出显著的数值阻尼，导致阻尼比被高估。
  • 效率：模拟 1ms 的衰减过程，MSM 仅需 20 秒，而显式 Full-FE 需 27.5 分钟，隐式 Full-FE 需 2 小时 23 分钟。
• 沉降现象分析：
  • 研究发现，经过冲击载荷后，系统会“沉降”到一个新的平衡点。
  • 残余应力：预紧后的残余切向应力几乎为零，但经过振动/冲击后，接触面会产生非零的**残余牵引力（Residual Traction）**分布。
  • 历史依赖性：系统的最终平衡状态取决于载荷历史（冲击次数和强度）。重复冲击会导致系统逐渐稳定，而大冲击（导致宏观滑移）会彻底重置接触状态。
  • QSMA 的局限：传统的 QSMA 基于预紧后的初始状态，无法捕捉这种由振动引起的平衡点偏移，导致其预测的固有频率与瞬态衰减（Ring-down）分析结果存在约 2% 的偏差。

5. 意义与结论 (Significance)

• 工程应用价值：该方法为复杂连接结构（如航空航天、汽车中的螺栓连接）的阻尼预测提供了高效、准确的工具，使得在考虑真实表面形貌的情况下进行长时程动力学仿真成为可能。
• 科学突破：
  • 解决了多尺度接触动力学中“计算可行性”与“数值精度”之间的矛盾。
  • 揭示了振动诱导沉降这一长期被忽视的物理现象，解释了实验测试中结果不可重复性以及准静态分析与动态分析结果差异的部分原因。
  • 证明了在弱非线性区域，接触界面的局部应力场细节虽然重要，但在投影到全局模态参数时，多尺度方法的近似是足够精确的。
• 未来展望：该方法为研究振动诱导的模态混合、不同尺度形貌对动力学的影响以及连接结构的长期可靠性评估开辟了新的途径。

总结：这篇论文提出了一种高效、稳健的多尺度动态仿真框架，成功克服了传统方法在解析表面形貌时的计算障碍和数值误差，不仅验证了其在预测结构阻尼方面的准确性，还通过数值实验揭示了振动导致结构平衡状态发生不可逆改变的物理机制。