飞轮储能PMSM建模与Simulink仿真：从原理到参数调试全解析

你搜“飞轮储能建模”，看到的往往是一堆精致的结构图、复杂的公式推导，然后就没了。自己动手在MATLAB/Simulink里搭的时候，各种问题就冒出来了：转速负值为什么会有、初始仿真步长怎么设、PI参数怎么给初始值、永磁同步电机的dq轴方程到底该用哪套符号约定，每一个细节都能卡住你大半天。

我最近把飞轮储能系统的完整仿真模型从头到尾搭了一遍，电机选的就是永磁同步电机（PMSM），从飞轮本体的机械方程、PMSM的电磁模型，到双闭环控制器、SVPWM、充放电逻辑，最后在Simulink里跑通并调出了理想波形。这篇文章就把整个建模思路、关键公式、参数取值、仿真配置，以及调试过程中踩过的坑全部摊开讲。内容面向电气工程、储能方向的学生，以及刚接触飞轮储能的工程师，目标是让你看完以后，能独立把一套可复现的“PMSM飞轮储能仿真模型”搭出来，而不是停留在看图说话的层面。

1. 先想清楚再建模：飞轮储能系统的整体设计思路

1.1 飞轮储能到底在存什么

飞轮储能的核心原理就是一句话：把电能转成旋转体的动能，需要时再把这部分动能变回电能。很多人听到“储能”第一反应是电池，但飞轮完全是另一条技术路线，它用的是机械能。

能量公式是 E = ½Jω²，J是飞轮转动惯量，ω是飞轮旋转的角速度。从这个公式能看出两个关键点：第一，转动惯量越大，储能越多；第二，转速的提升对储能的贡献是平方级的，所以现代飞轮储能都在拼命追求高转速，而不是单纯把飞轮做大。

我举个例子，假设飞轮转动惯量取 J = 0.05 kg·m²，额定转速设为 20000 r/min，换算成角速度是：

ω = 2π × 20000 / 60 ≈ 2094.4 rad/s

那么额定转速下存储的能量就是：

E = 0.5 × 0.05 × 2094.4² ≈ 109.7 kJ

如果系统充放电功率是5kW，理论上从静止加速到额定转速需要约22秒，从额定转速放电到停机大约也是22秒。这就是飞轮储能的时间尺度，它适合短时、高频、大功率的功率型应用，比如电网一次调频、UPS应急供电、轨道交通再生制动能量回收、新能源波动平抑。它的优势和劣势都很鲜明：功率密度高、响应快、循环寿命长，但能量密度远低于电池，持续放电时间短。

1.2 为什么飞轮电机大多是永磁同步电机

飞轮储能的关键就是电机既要当电动机又要当发电机，而且要在高速下长时间稳定运行。电机选型是整个建模的前提。

对比几个常见方案：

• 异步电机结构简单，高速运行稳定性好，但转子有铜耗和铁耗，效率偏低，功率因数控制复杂，作为发电机时还需要励磁电流，控制难度上升。
• 电励磁同步电机励磁灵活，可控性强，但转子绕组需要励磁电流，高速运行时滑环或无刷励磁系统可靠性是问题。
• 开关磁阻电机结构最坚固，适合超高速，但转矩脉动和噪声比较明显，控制模型是非线性的，建模和仿真难度都更高。
• 永磁同步电机（PMSM）转子本身就是永磁体，不需要励磁电流，没有转子铜耗，效率在整个负载范围内都很高；而且它的功率密度大、转矩响应快、控制精度高，矢量控制技术非常成熟。

所以目前的飞轮储能系统大量采用PMSM，尤其是表贴式永磁同步电机（SPMSM），因为表贴式结构适合高速旋转，转子结构对称，没有磁阻转矩，配合id=0控制非常简单高效。如果转速要求特别高，还会用内置式（IPMSM）配合弱磁控制。

我的建议是：初学建模就从表贴式PMSM开始，先跑通id=0控制下的充放电流程，后面再去碰弱磁和内置式的复杂控制。这就像学跑步先走直线，别一上来就搞弯道冲刺。

1.3 建模时要把系统拆成哪四层

飞轮储能系统在仿真里不是一个孤立电机，而是一个完整的“机械-电磁-控制-电网”耦合系统。我的习惯是把它拆成四个层次来建模：

第一层，机械层。核心是飞轮转子本体、转动惯量J、轴承摩擦、风阻损耗，这一层决定了飞轮的转动特性。

第二层，电磁层。即PMSM本体，包括电压方程、磁链方程、转矩方程，以及从abc三相到dq两相旋转坐标系的坐标变换。

第三层，控制层。这是整个仿真的灵魂，包括速度外环、电流内环、PI调节器、SVPWM调制、id=0矢量控制，以及充放电模式切换逻辑。

第四层，电网交互层。包括逆变器、直流母线电容、电网侧变流器，以及母线电压稳定控制。

这样分层的好处是：每一层都能单独验证。机械层可以单独测空载转速跌落曲线；电磁层可以用恒转矩负载验证电磁转矩；控制层可以暂时抛开电网，直接给定参考转速；最后再整体联调。如果你把三层揉在一起，出了问题根本分不清是控制算法的问题还是电网侧设置的问题，排查起来非常痛苦。

2. 机械和电磁这两套数学模型，一个都不能错

2.1 飞轮本体的机械运动方程

飞轮转子的运动可以用一个方程描述：

J·dω/dt = Te - TL - Tloss

其中Te是电机的电磁转矩，TL是负载转矩，Tloss代表所有损耗转矩。

在飞轮储能系统里，没有外部机械负载，所以TL直接取0，剩下的是损耗项。损耗这部分，很多论文里只写一个阻尼系数B，把它和转速成正比，写成Bω。但我实际建模时发现，真实飞轮损耗更复杂，至少要分三部分：

• 恒定摩擦转矩：主要来自轴承摩擦，近似为一个恒定值，方向始终与旋转方向相反。
• 粘滞阻尼：和转速成正比，写为Bω。
• 风阻损耗：高速飞轮转速可达每分钟两万转以上，风扇效应明显，损耗近似与转速的平方或三次方成正比，写为Cω²。

所以完整的机械模型建议写成：

J·dω/dt = Te - Tf·sign(ω) - Bω - Cω²·sign(ω)

注意符号处理。飞轮在充电时，Te为正，机械能增加，转速上升；放电时，Te为负，机械能减小，转速下降。损耗项的方向始终和转速方向相反，所以在离散仿真里一定要加符号判断，否则当转速降到接近零时模型还会继续耗能，结果转速变成负的，那就是明显错误。

转动惯量J不是随便给的。如果是“飞轮+电机转子”一体化设计，J就是两者之和。估算方法可以用圆柱体转动惯量公式：

J = ½·m·r²

m是飞轮质量，r是质心到轴线的等效半径。要注意，飞轮材料选择直接决定安全转速上限，高强度钢材可以做到每分钟上万转，碳纤维复合材料可以做到更高，但成本也随之上升。在仿真里，J越大，系统时间常数越大，转速响应越慢，启动需要的时间越长；J太小，则转速变化太快，对控制器带宽要求更高。这是一个需要在项目初期就定下来的核心参数。

2.2 永磁同步电机的dq坐标系数学模型

PMSM建模的核心是把三相静止坐标系下的交流量，变换成两相旋转坐标系下的直流量。这个变换过程叫Park变换和Clark变换。变换之后，电机方程就从时变的交流方程变成了定常的直流方程，控制精度和仿真效率都大幅提升。

标准PMSM在dq坐标系下的电压方程为：

ud = Rs·id + Ld·did/dt - ωe·Lq·iq uq = Rs·iq + Lq·diq/dt + ωe·Ld·id + ωe·ψf

磁链方程：

ψd = Ld·id + ψf ψq = Lq·iq

电磁转矩方程：

Te = 1.5·p·[ψf·iq + (Ld - Lq)·id·iq]

对于表贴式PMSM，Ld约等于Lq，所以转矩方程简化为：

Te = 1.5·p·ψf·iq

这就是为什么id=0控制那么重要——id等于零时，磁链完全由永磁体提供，转矩只与iq成正比，转矩线性、可控、无磁阻分量，控制最简单。

有几个特别容易被新手忽略的点：

第一，机械角速度和电角速度不是一回事。电角速度ωe等于机械角速度ωm乘以极对数p，即ωe = p·ωm。转矩公式里的速度项全部要用电角速度，而机械运动方程里的转速是指机械角速度。很多人在仿真里把这两个量搞混，直接导致转矩计算差了好几倍。

第二，ud和uq方程里的交叉耦合项，即ωe·Lq·iq和ωe·Ld·id，在转速较高时非常显著。如果想提升仿真精度，控制器里应该加入前馈解耦项，把耦合项抵消掉。否则高速时你会发现，设定iq不变，实际的转矩响应却随着转速漂移，这就是耦合项在捣乱。

2.3 仿真参数怎么取值才靠谱

参数是仿真能否反映真实系统的关键，乱填参数等于拿着橡皮子弹练枪法。下面是一组我实际使用的是飞轮储能仿真参数，可供参考：

磁链ψf的工程取值方法是：先看电机额定反电动势有效值E0，然后 ψf = E0·√2/ωe。注意这里的ωe是电角速度，不是机械角速度，别又把极对数漏了。

还有一个反直觉的细节：飞轮储能的PMSM额定转速比普通伺服电机高得多，普通工业伺服电机往往额定3000r/min左右，而飞轮储能电机动辄上万转。在这个转速区间，电机设计要考虑高速离心应力、转子结构强化和风摩耗散热问题。所以在仿真里把转速上限压得很保守，可能反而掩盖了高速段控制交互的问题。

3. 控制策略和充放电切换逻辑的设计

3.1 双闭环PI控制的带宽分配

飞轮储能的PMSM控制，最常见的是“速度外环、电流内环”的双闭环结构。之所以要双闭环，是因为电流环可以快速抑制扰动，限制电流冲击，而速度环负责跟踪转速指令，两者配合才能兼顾快速性和稳定性。

电流内环把PI调节器的输出作为dq轴电压给定，再进行Park逆变换得到三相电压参考，经过SVPWM驱动逆变器。速度外环则把转速误差转化为iq给定，也就是转矩电流指令。

带宽分配的原则是：内环带宽要远高于外环，一般内环截止频率是外环的5到10倍。这样做能保证时间尺度分离，外环变化时内环能迅速跟随，不会出现内外环互相纠缠的振荡。

给个从零开始的设计步骤：

第一步，确定电流环截止频率。假设开关频率是10kHz，电流环截止频率通常取开关频率的1/10到1/20，这里取fcc = 1kHz，即ωcc ≈ 6283 rad/s。

第二步，计算电流环PI参数。工程上常用零极点对消思想，把PI调节器的零点抵消掉定子电感的极点，所以：

Kpc = ωcc·L = 6283 × 0.005 ≈ 31.4

Kic = Kpc·Rs/L = 31.4 × 0.5 / 0.005 ≈ 3140

但这是连续时间域的算例，实际在Simulink离散仿真中，离散PI的Ki数值含义不同，往往需要缩小若干倍再微调。这也是为什么你不能把教材公式算出的参数直接塞进离散模型就指望它完美工作。

第三步，确定转速外环参数。转速环截止频率取电流环的1/10左右，比如ωcn ≈ 628 rad/s。把电磁转矩系数视为常数Kt = 1.5·p·ψf = 0.6，于是：

Kps ≈ J·ωcn / Kt = 0.05 × 628 / 0.6 ≈ 52.3

Kis一般取Kps的1/5到1/3，这里取10左右。

反正记住：这是一组可以起调的初始值，在仿真里要根据波形继续修正。调参的秩序永远是先电流环再速度环，如果电流环都振荡，速度环参数再准也没用。

3.2 充电、放电、待机三种工况怎么切换

飞轮储能系统一共有三种工况：充电、放电、待机。充电时，变流器从电源取电，电机工作在电动状态，电磁转矩方向与飞行方向一致，飞轮加速，动能增加。放电时，电机进入发电状态，电磁转矩方向与飞行方向相反，飞轮减速，动能通过电机和变流器回馈到直流母线，再回馈给负载或电网。待机时，飞轮维持在额定转速附近，电磁转矩几乎为零，只用来补偿摩擦和风阻损耗。

在Simulink里，最直接的模式切换方式是：通过一个转速指令生成器来定义整个运行过程。比如：

• 0~20秒：转速指令从0斜坡上升到20000r/min，对应充电过程。
• 20~40秒：转速指令保持20000r/min不变，对应待机，母线侧维持电压，功率接近零。
• 40~55秒：转速指令从20000r/min斜坡下降到12000r/min，对应放电过程，储能部分释放。

这样做的优点是：不需要频繁切换控制器结构，只需要改变转速指令和负载方向即可，实现简单，现场经验也比较好对应。需要注意的是，转速指令绝不能给阶跃突变，而必须是斜坡或积分限速。阶跃指令会让惯性极大的飞轮瞬间产生巨大的电流冲击，结果就是母线电压跌穿、电流保护动作、系统停机。

关于功率方向，我用一个简单方法判断：充电时iq给定为正，放电时iq给定为负。因为表贴式PMSM的转矩只和iq成正比，正负iq就是正负转矩。你只要在控制器里把转速误差倒过来，放电工况自然就能出来。

3.3 和电网侧怎么协同

完整飞轮储能系统不是电机直接接电网，而是通过背靠背变流器接入。电网侧变流器负责维持直流母线电压稳定，电机侧变流器负责控制PMSM的转矩和转速。两者通过直流母线电容解耦，母线电压就在两者之间搭桥。

电网侧的控制核心是电压外环、电流内环的双闭环，目标是让直流母线电压保持在给定值附近。当电机侧放电回馈能量时，母线电压会上升，电网侧控制就反向增大输出功率，把能量送到交流母线；当电机侧充电时，电网侧从电网取能量，母线电压下降，然后被电网侧控制拉回额定值。

在仿真里，如果你只关心飞轮本体的控制策略，可以先不接电网侧，用一台直流电压源代替直流母线，并在母线电容上并联一个可编程负载或泄放电阻。这样模型更简单，也不影响验证电机侧的充放电逻辑。只有当你想研究功率双向流动或电网频率支撑时，才需要把电网侧变流器完整建模进去。

4. MATLAB/Simulink仿真搭建与参数配置实例

4.1 仿真模型用哪种结构

我个人强烈建议把模型搭成子系统分层的结构，而不是把所有模块平铺在主模型里，否则后面排查问题无异于大海捞针。

我的Simulink模型分为以下几个子系统：

• 飞轮机械子系统：输入为电磁转矩，输出为机械转速，内部实现运动方程和损耗模型。
• PMSM模型子系统：输入为机械转速和定子电压，输出为dq电流和电磁转矩。
• 控制器子系统：包括速度PI、电流PI、坐标变换、SVPWM或平均调制。
• 直流母线子系统：包括母线电容、负载、电网侧接口。

关于电机本体，MATLAB的电气仿真模块库里就有PMSM模型，参数配置直观，适合快速出结果。但如果你想把底层细节全部掌握，可以自己搭基于方程的PMSM模型，即上一节那些dq方程。我建议初学者先用库里的模型，跑通系统；有条件再自建电机方程模块，两者对比验证。自建模型有个好处是你能完全控制每个中间变量，排错能力会大幅提升。

关于逆变器模型，推荐两种方式：

• 平均模型：逆变器被简化成按调制波输出的理想电压源，仿真速度快，适合验证控制策略和系统级动态。
• 开关模型：使用实际开关管和SVPWM调制器，可以看到谐波和开关纹波，但仿真步长要小、速度慢，模型也更复杂。

初学就从平均模型起步，先把控制系统调好，再升级到开关模型。

4.2 核心模块配置与仿真步长选择

仿真步长的选择是一个经常影响成败的细节。如果用的是平均模型，离散定步长取1e-5秒（10微秒）通常足够，既能看到机电动态，又不会太慢；如果用了开关模型，想要看到PWM纹波，步长要压到1e-6秒或更小。

其实这里要区分两件事：开关频率导致的电磁纹波时间尺度，和机械转速变化的机电时间尺度。飞轮转动惯量大，机电时间常数在几秒甚至几十秒；但如果步长取太大，电机电流和控制器离散化误差会积累，仿真结果会出现虚假振荡。所以推荐的策略是：离散定步长10微秒，先跑通，再根据结果调节。

信号观测方面，我建议使用Simulink的Data Inspector，或者直接在示波器里记录这些信号：机械转速、电磁转矩、dq轴电流、三相电流、母线电压、储能能量，以及PI调节器的输出限幅状态。记录PI输出限幅状态非常关键，因为很多PID振荡和发散问题，根源是输出达到了限幅值而你不知道。

我在模型里还加了一个“转速上限保护”模块：一旦转速超过允许上限，立即停止功率输出或切换为能耗制动。这不是多此一举，飞轮超速在物理世界是灾难级的故障，仿真阶段养成这个习惯，后面做硬件不会吃亏。

4.3 仿真结果怎么看

我用前面那组参数，跑完充电-待机-放电全流程后，几个关键波形的预期是这样的。

充电阶段（0~20秒），转速指令按斜坡从0升到20000r/min，机械转速曲线线性上升，电磁转矩基本恒定在一个正数附近，用来克服惯性和损耗。直流母线持续向电机侧输出功率，母线电压略降但电网侧控制会快速拉回。

待机阶段，转速保持在20000r/min附近，电磁转矩趋近于摩擦和损耗项的补偿值，功率明显减小，母线电压平稳。

放电阶段（40~55秒），转速指令斜坡下降，机械转速跟随下降。电磁转矩变为负值，对应发电制动。此时母线电压会有短暂上升，如果接了电网侧模型，功率会回馈到交流侧；如果只接泄放电阻，则能量消耗在电阻上。

判断仿真是否合理的几个标准：转速指令与实际转速误差是否在预期范围；dq电流是否快速收敛且无高频振荡；母线电压是否维持在5%以内波动；充放电过程的能量变化是否大致符合E=½Jω²的理论曲线。

5. 调试中踩过的坑和排查方法论

5.1 转速指令突变导致系统失控

我第一次用阶跃转速指令测试时，模型直接崩了，转速曲线像过山车一样上下跳，母线电压跌穿到接近零。原因其实很简单：飞轮转动惯量太大，转速不可能突变，但控制器为了追阶跃参考，瞬间输出最大转矩，电流达到限幅值，母线侧提供的功率跟不上，最终系统失去平衡。

解决方法是给转速指令加斜坡生成器，限制最大加/减速度。具体斜坡斜率怎么定？根据额定功率和储能容量：

充电从0到20000r/min按5kW计算需要约22秒，所以上升斜率约为909 r/min·s⁻¹。你就用这个数量级设定斜坡，留出控制器调节余量。

这个坑给我的经验是：仿真模型里的指令约束，本质上就是现实系统里的“机组运行规程”。在仿真里乱下阶跃指令，和现实中一上来就满功率冲击设备是一回事，都不合理。

5.2 双闭环PI参数怎么调都振荡

我遇到的第二个大坑是：用了理论计算的PI初始参数，模型一跑，电流环直接发散。排查了很久才发现问题不在PI增益本身，而在于离散化步长和参数匹配。

电流环时间常数相当于L/R = 0.005/0.5 = 0.01秒，也就是10毫秒；我这个系统的电流动态在毫秒量级，而仿真步长10微秒足够覆盖。但理论参数在连续域的Kic算出来是3140，在离散仿真里，这个Ki按离散积分离散化后效果相当于很大的积分增益，很容易振荡。应对办法是：把Ki降低约一半，或者先只加比例项，让电流环稳定到跟随给定，再逐步增加积分项，直到静态误差消失。

调试顺序再强调一遍：先把电流环的iq给定固定成一个常数，用示波器看dq电流是否跟踪给定，调稳定了再放开速度环。速度环也类似，先把pi设为纯比例，调到转速振荡临界，再按比例退一半，然后加积分消除稳态误差。这套方法比盲目复制公式参数再乱调要靠谱得多。

5.3 仿真波形和理论预期对不上

出现过几次“理论计算说的是A，仿真结果是B”的情况。排查后基本都出在这几类：坐标变换方向反了（dq轴顺序或正负号约定不同）、机械角速度没有乘以极对数、损耗模型取了过大的阻尼系数。

PMSM模型库和自建模型，对dq轴的正方向约定不完全一致，有的模型id为负表示弱磁，有的则相反。你在看模块文档时要特别注意转矩表达式里的符号，否则电磁转矩方向一反转，充放电就全反了。

排查这类问题的方法很简单但很有效：做开环验证。断开控制器，直接给电机加一个恒定的电压，看电机能否正常起动、转速方向是否正确、空载电流是否符合预期。开环验证过了，基本可以判断电机模型和参数没错，接下来再回到闭环结构排查控制器。

5.4 母线电压飞升和变流器饱和问题

放电回馈时，我发现母线电压快速上升，超出电容允许范围。这是因为放电回馈功率大于负载吸收功率，直流侧能量无处可去，全堆在电容里，电压当然飙升。这在真实系统里会直接打坏电容和功率管。

解决方案有两条路：一是在直流母线加卸荷电路，当母线电压超限时接入泄放电阻消耗多余能量，仿真里用一个电压比较器加受控开关切换电阻就行；二是把电网侧变流器的功率指令改为“母线电压外环”，让电网侧自动根据母线电压高低调整回馈功率。实际工程中两者常同时使用，卸荷电路作为最后保护措施。

另外注意，仿真里变流器其实有容量限制，如果不给电流限幅，PI输出会任意增大。我在控制器里给iq限幅设为额定电流的1.5倍，给id限幅设为额定电流，这样充放电功率不会超出可实现的电磁转矩范围，仿真结果才贴合真实系统。

5.5 常见问题排查速查表

最后再分享一个小细节：仿真跑通不等于参数合理，一定要做能量守恒校验。把充放电全过程的能量变化和电机输出功积分对比，误差超过5%就要回头看损耗模型或数值积分选项。储能系统的核心是能量，能量对不上账，一切结论都不可靠。我自己调试时，能量闭环校验几乎每次都能揪出至少一个隐藏的参数问题。