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 知识体系总览
为了让你更直观地理解本章的知识结构,我画了一张图:
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