⚡ electrical engineering

S-Shaped Bifurcation Cascades and Loop Proliferation in an Active Magnetic Bearing Rotor with Cubic Shaft Stiffness and Nonlinear Aerodynamic Drag

本文利用谐波二次型法（Harmonic Quadrature Method）和弗洛凯稳定性分析（Floquet stability analysis），揭示了比例增益、微分增益、三次轴刚度和非线性气动阻力如何相互作用，从而驱动主动磁轴承支撑柔性转子中复杂的S型分岔级联、包含超过36个鞍节点分岔的广泛环路增殖以及独特的阻尼机制。

原作者： Farouk Thaljaoui, Thabet Guesmi, Thaljaoui Wathek

发布于 2026-08-05

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Farouk Thaljaoui, Thabet Guesmi, Thaljaoui Wathek

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一个沉重的旋转机器完全不接触地面的世界。它们不是在轮子上滚动或在油液中滑动，而是悬浮在半空中，由无形的磁力之手支撑着。这就是主动磁轴承（AMB）的魔力，这项技术被广泛应用于从高速喷气发动机到能量存储飞轮的各种设备中。但问题在于：磁铁非常难以捉摸。它们产生的拉力并不是线性增长的；随着距离缩短，力量会变得异常强大，就像一根被拉得越长、回弹力就越强的橡皮筋。当这种磁性的“回弹”与旋转轴的自然晃动以及高速运动物体的空气阻力结合在一起时，系统可能会表现得像一个混乱的舞伴。它可能会突然发生剧烈的、危险的震颤，或者陷入一种不知道该往哪边走的循环状态。工程师需要理解这些跳跃和循环，才能防止他们的机器因剧烈震动而解体。

这篇论文深入研究了这场混乱的舞蹈，利用强大的计算机模拟，精确绘制出当调整控制参数时，悬浮转子是如何表现的。研究人员将转子视为一根中间带有重型圆盘的柔性杆，并由使用“PD控制器”（一种基于位置和速度来调整磁力强度的智能大脑）的磁轴承支撑。他们还加入了两个额外的成分：轴在弯曲时会变硬（就像一个越压越硬的弹簧）以及随转子振动速度变化的空气阻力。通过运行数千次模拟，他们发现机器的“大脑”有两个旋钮：一个是比例控制旋钮（P），控制磁力的大小；另一个是微分控制旋钮（D），控制对运动的阻尼。

研究团队发现，调大比例旋钮（P）会让系统变得越来越狂野。在低设置下，转子的表现很温顺，像一个平静的单摆。但随着他们调高 P 值，转子运动的平滑曲线开始自我折叠，形成了“S形”。这就像是一段突然回环自折的过山车轨道。起初只有一个环（1S形），接着有两个（2S形），最后变成了三个（3S形）。在这些 S 形区域内，转子对于以相同频率振动的选择是多种多样的。它既可以轻微振动，也可以剧烈摇晃，而一个微小的扰动就能让它瞬间从安全状态跳跃到危险状态。研究人员计算出，当 P 从 2 增加到 4 时，最大振幅从 0.2328 增长到了 0.3839，并且这些“跳跃点”（称为鞍节点分岔）的数量也成倍增加，创造了一个由于机器行为完全取决于其历史路径而导致的混乱地带。

然而，微分旋钮（D）就像是一个相反的超级英雄。虽然 P 会制造混乱，但 D 会消灭混乱。论文显示，增加 D 就像是在火上浇水，能抚平那些狂野的 S 形环路。当 D 较低时，系统是一团乱麻；但随着研究人员调高 D，这些环路逐渐缩小。在 D = 1.23 的特定设置下，混乱完全消失了，系统回归到一种平静、可预测的状态。这就像是 P 旋钮构建了一个迷宫，而 D 旋钮则是开启通往直线路径之门的钥匙。

最令人惊叹的发现发生在他们开启“轴向硬化”（三次刚度）时。当他们将磁性的“软化”与轴的“硬化”结合在一起时，发生了一些奇特的事情：S 形环路不仅发生了折叠，还扩散成了一个巨大的、连通的闭合环路。想象一个缠绕在一起且最终回到自身的毛线球。在这个特定区域，研究人员发现了一个解的分支，它扭曲、转弯，拥有超过 36 个“折点”（路径折返之处）。当轴的刚度增强到足以与磁力抗衡时，这种情况就会发生，达到一个轴力约为磁力 8.6% 的比例。在这个区域内，转子拥有数十种不同的稳定振动方式，使得其行为在没有完美地图的情况下极难预测。

最后，团队研究了空气阻力。他们发现，空气对转子的作用力取决于其振动程度的不同。当转子剧烈震动（主共振）时，空气表现得像粘稠的泥浆（非线性阻尼），这有助于抑制大幅度的跳跃。但当转子进行微小的、高频的颤动（超谐波）时，空气则表现得像稀薄的水浆（线性阻尼）。如果仅用其中一种模型来模拟空气，结果将会是错误的。论文证明，你需要两种模型才能得到完整的图景，因为“泥浆”主导了大振幅，而“水浆”则处理了微小振动。

简而言之，这篇论文为磁轴承那危险而迷人的领域提供了一份详细的地图。它告诉工程师哪里存在“S 形”陷阱，如何利用微分旋钮来逃脱，并警告他们在轴刚度与磁力相互博弈时可能出现的隐藏“环路增殖”。这是一本防止高速机器变成不可预测的震动噩梦的指南。

技术摘要：有源磁轴承转子中的 S 型分岔级联与环路增殖

问题陈述
有源磁轴承（AMB）是高速旋转机械的关键部件，通过比例-微分（PD）反馈控制实现无接触悬浮。然而，AMB 的电磁力本质上是非线性的（与气隙平方成反比），这可能诱发复杂的动力学行为，如多值频率响应、振幅跳跃和滞后现象。现有文献大多孤立地或简单组合地研究电磁非线性、三次轴刚度和非线性气动阻尼。本研究旨在填补理解这些竞争性非线性机制——特别是磁力引起的软化特性与三次轴刚度引起的硬化特性，以及随振幅变化的气动阻尼之间的相互作用——在单一柔性转子模型中的空白。主要目标是表征由此产生的分岔拓扑结构，并识别产生复杂现象（如具有众多鞍节点分岔的闭环分支）的条件。

方法论
作者将受 AMB 支持的柔性转子建模为一个四自由度系统，包含：

非线性电磁力： 由 PD 反馈控制驱动，其中电流是转子位移和速度的函数。
三次轴刚度： 代表中等或大挠度下的几何非线性。
非线性气动阻尼： 建模为线性粘性阻尼（低速/小振幅）与二次湍流阻尼（高速/大振幅）的结合。

利用拉格朗日形式推导运动方程并进行无量纲化处理。为了求解该系统，研究采用了谐波正交法（HQM），这是一种时域离散化技术，在保持高阶精度的同时避免了代价高昂的 Galerkin 积分。该方法结合伪弧长连续法来追踪参数空间内稳定和不稳定周期解的分支，并结合 Floquet 稳定性分析来验证稳定性并检测鞍节点分岔。该数值框架通过关于谐波阶数和弧长步长的收敛性研究进行了验证。

核心贡献与结果
研究系统地调查了四个关键参数的影响：比例增益 ( $P$ )、微分增益 ( $D$ )、非线性刚度系数 ( $\eta_{NL}$ ) 以及气动阻尼系数。

比例增益 ( $P$ ) 与 S 型级联：
增加 $P$ 会驱动一个渐进的分岔级联，经历四个不同的阶段：
- 线性阶段 ( $P < 2.68$ )： 单值、稳定的响应。
- 1S 型阶段 ( $2.69 \le P \le 2.75$ )： 出现单个 S 型折叠，具有两个鞍节点分岔。
- 2S 型阶段 ( $2.76 \le P \le 3.0$ )： 两个 S 型分支，增加了多稳态性。
- 3S 型阶段 ( $3.1 \le P \le 4.1$ )： 三个 S 型分支，具有八个折叠点。
  随着 $P$ 从 2 增加到 4，主共振振幅从 0.2328 上升到 0.3839，且多稳态频带显著变宽。
微分增益 ( $D$ ) 与逆级联：
微分增益起到了稳定性的对抗措施作用。增加 $D$ 会产生与 $P$ 完全相反的逆级联过程：
- 在低 $D$ 时 ( $0.5 \le D \le 1.16$ )，系统处于 3S 型阶段。
- 在中等 $D$ 时 ( $1.16 \le D \le 1.23$ )，系统过渡到 1S 型阶段。
- 当 $D \ge 1.23$ 时，多稳态性被破坏，恢复为准线性行为。
  这揭示了一种对称的、互补的关系，即 $P$ 诱发复杂性，而 $D$ 系统地抑制复杂性。
非线性刚度 ( $\eta_{NL}$ ) 与环路增殖：
当系统固定在 3S 型阶段并增加 $\eta_{NL}$ 时，硬化机械非线性与软化磁非线性展开竞争。在一个临界阈值（ $\eta_{NL} \approx 0.55$ ）处，系统在主响应曲线之上生成了一个连通的闭环分支。
- 该环路包含超过 36 个经 Floquet 验证的鞍节点分岔点，分布在频率区间 $\Omega \in [0.2466, 0.2922]$ 内。
- 力对比证实这是一个物理现象：在环路出现时，三次项力与磁力的比值从 0.002% 增加到 8.6%，表明硬化刚度已足够强大，足以扭曲软化磁响应。
气动阻尼阶段转换：
研究确定了阻尼主导地位随振幅变化的转换：
- 非线性湍流阻尼支配主共振，有效地抑制了大振幅 S 型环路。
- 线性粘性阻尼控制超谐波峰值，此时振幅较小。
  在 $\Omega \approx 0.27$ 附近存在一个交叉点。仅使用一种阻尼机制的模型会系统性地误判响应振幅（例如，仅线性模型预测的主共振振幅比非线性模型高出约 57%）。

意义与主张
本文声称提供了在同时存在软化电磁力、硬化三次刚度和随振幅变化阻尼的情况下，对 AMB 转子动力学进行全面表征的首次尝试。作者断言：

所识别的 S 型分岔级联 为调节 PD 控制器提供了定量设计图表，定义了保证单值运行的具体阈值（例如， $P < 2.68$ ）。
通过 $D$ 实现的逆级联建立了消除多稳态性的定量稳定性边界（ $D \ge 1.23$ ），为在高增益阶段运行的系统提供了实用的控制策略。
包含超过 36 个折叠点的连通闭环分支的发现，代表了一种不同于简单 S 型折叠的复杂拓扑变换，突显了在特定参数范围内由于历史依赖性导致的转子响应的不可预测性。
双重阻尼公式被认为是高转速 AMB 转子建模的必要手段，以准确捕捉主共振和超谐波行为。

研究结论认为，这些发现表征了 AMB 系统的非线性分岔拓扑，为 PD 控制器的安全调优和在复合非线性存在下的转子动力学预测提供了直接指导。未来的工作建议探索被动稳定（挤压膜阻尼器）和减轻上述环路增殖现象的活性控制策略。

类似论文