2. 磁悬浮轴承系统建模:单自由度系统动力学模型、电磁力线性化模型、系统传递函数推导

做磁悬浮轴承控制,第一步就是建模。

很多人一上来就调PID参数,结果系统抖得像筛子。为什么?因为你根本不知道被控对象长什么样。我当年刚入行时也犯过这个错,后来被老工程师骂了一顿,才老老实实回来推公式。

说白了,建模就是搞清楚两件事:转子怎么动,以及电磁力怎么给。搞清楚了这两点,后面的控制器设计才有根有据。

2.1 单自由度系统动力学模型

先看最简单的单自由度系统。一个转子,悬浮在磁场中,只考虑一个方向的运动——比如竖直方向。

根据牛顿第二定律:

m * d²x/dt² = F_mg - mg

其中:

  • m — 转子质量
  • x — 转子在竖直方向的位移
  • F_mg — 电磁合力
  • g — 重力加速度

嗯,这里要注意:我们通常把平衡位置设为原点。也就是说,当转子稳定悬浮时,电磁力刚好抵消重力,此时 x = 0。

所以更常用的写法是:

m * d²x/dt² = F_mg - mg

这个方程看起来简单,但它是整个控制系统的根基。你想想看,如果连这个都搞错了,后面所有的传递函数、控制器设计全是白搭。

核心要点:单自由度模型是磁悬浮轴承控制的基础。实际工程中,我们通常把五自由度系统解耦成五个单自由度系统分别控制。我做过的一个项目里,就是因为忽略了耦合项,导致系统在高速旋转时突然失稳——那次教训太深刻了。

2.2 电磁力线性化模型

接下来是重头戏——电磁力怎么算?

电磁力公式长这样:

F = k * (i / x)²

其中:

  • k — 电磁常数(跟线圈匝数、磁极面积、气隙长度有关)
  • i — 线圈电流
  • x — 气隙长度(转子与磁极之间的距离)

问题来了:这是个非线性关系。i 和 x 都在变,F 跟它们的平方成正比。直接拿这个去设计控制器?别闹了,非线性系统控制可不是闹着玩的。

所以我们需要线性化

怎么做?泰勒展开。在平衡点 (i₀, x₀) 附近展开:

F ≈ F₀ + k_i * Δi + k_x * Δx

其中:

  • F₀ — 平衡点电磁力(等于 mg)
  • k_i — 电流刚度系数(∂F/∂i 在平衡点的值)
  • k_x — 位移刚度系数(∂F/∂x 在平衡点的值)
  • Δi — 电流变化量
  • Δx — 位移变化量

个人经验:我建议你在实际项目中,k_i 和 k_x 不要只靠理论计算。最好用有限元仿真校核一下,有条件的话做实验标定。我曾经在一个精密加工项目中,理论算的 k_i 跟实测差了 15%,导致控制器参数全得重调。

线性化后的电磁力模型,说白了就是一个电流力加上一个位移力。电流力是我们能控制的,位移力是负刚度效应——转子越靠近磁极,吸力越大,这是个不稳定的正反馈。

2.3 系统传递函数推导

有了动力学模型和线性化电磁力模型,接下来就可以推导传递函数了。

把线性化电磁力代入动力学方程:

m * d²Δx/dt² = k_i * Δi + k_x * Δx

整理一下:

m * d²Δx/dt² - k_x * Δx = k_i * Δi

做拉普拉斯变换(初始条件为零):

m * s² * X(s) - k_x * X(s) = k_i * I(s)

所以从电流到位移的传递函数为:

G(s) = X(s) / I(s) = k_i / (m * s² - k_x)

你看,分母是 m*s² - k_x,不是 m*s² + k_x。这意味着什么?

系统有一个极点位于正半平面。

换句话说,开环系统是不稳定的。这就是为什么磁悬浮轴承必须用闭环控制——没有反馈,转子根本悬浮不起来。

避坑指南:我曾经见过有人直接把 G(s) 当成普通二阶系统来设计控制器,结果系统怎么调都不稳定。原因就是没注意到 k_x 前面的负号。记住,磁悬浮系统的传递函数分母是 m*s² - k_x,不是 m*s² + k_x。这个负号决定了系统的本质特性。

2.4 知识体系总览

为了让你更直观地理解本章的知识结构,我画了一张图:

磁悬浮轴承系统建模知识体系 磁悬浮轴承系统建模 单自由度动力学模型 电磁力线性化模型 系统传递函数推导 核心方程 m·d²x/dt² = F_mg - mg 牛顿第二定律应用 确定被控对象运动规律 线性化方法 F ≈ F₀ + k_i·Δi + k_x·Δx 泰勒展开在平衡点 电流刚度 + 位移刚度 传递函数形式 G(s) = k_i / (m·s² - k_x) 拉普拉斯变换 开环不稳定,需闭环控制 核心结论:开环不稳定 → 必须闭环控制 → 控制器设计的基础 三个模型层层递进,缺一不可

2.5 关键参数总结

为了方便你查阅,我把本章涉及的关键参数整理成了一张表:

参数符号 物理含义 单位 典型取值范围
m 转子质量 kg 0.1 ~ 100
k_i 电流刚度系数 N/A 10 ~ 1000
k_x 位移刚度系数 N/m 10⁴ ~ 10⁶
x₀ 平衡气隙 mm 0.1 ~ 1.0
i₀ 偏置电流 A 0.5 ~ 10

一个小建议:拿到一个新项目,我习惯先把 m、k_i、k_x 这三个参数搞清楚。有了它们,系统的开环特性基本就定了。剩下的就是怎么设计控制器去压住那个正半平面的极点。

好了,这一章的内容就到这里。建模是控制的基础,基础打牢了,后面调参才能有的放矢。


公众号:蓝海资料掘金营,微信deep3321