飞轮储能这个方向最近两年在我接触的储能同行里聊得越来越多。它不像锂电池那样依赖电化学反应而是把能量以动能形式“存”在一个高速旋转的飞轮里靠电机实现电能和机械能的双向转换。做这个仿真项目的初衷是想把飞轮储能系统里最核心的“电机—飞轮—变流器”这套电气-机械耦合关系摸清楚。我选择永磁同步电机作为飞轮驱动电机在Simulink里搭了一套可复现、可扩展的仿真模型这篇文章就是把从零搭建到调试验证的全过程做一个复盘包括坐标变换原理、SVPWM逆变器设计、转动惯量匹配、PI参数整定这些关键细节以及我在调参过程中踩过的坑。这套内容适合刚接触飞轮储能的研究生、做电机驱动应用开发的工程师也适合对储能系统建模仿真感兴趣但还没找到切入点的入门者。我在文中会尽量把“为什么要这么做”讲透而不是只给一堆模型截图毕竟Simulink模型的核心价值不在模块连线而在你清楚知道每根线背后对应的物理方程和控制逻辑。1. 项目背景与整体设计思路1.1 飞轮储能的能量转换链路本质飞轮储能系统听起来“高大上”但物理本质其实一句话就能说清楚把一个质量较大的转子加速到很高的转速把电能变成动能“存起来”需要的时候再让转子带动电机发电把动能变回电能。整个系统的核心矛盾在于“转速范围”和“能量转换效率”之间的博弈。转速越高单位质量储存的能量越大但机械损耗、轴承损耗、电机高频损耗也会显著上升。从能量角度来说飞轮储存的动能公式是 E 0.5 × J × ω²其中 J 是飞轮转动惯量ω 是角速度。由于转速是平方项转速提升对储能容量的贡献是决定性的。但高转速对电机和轴承的材料强度要求也更苛刻。这就解释了为什么飞轮储能系统的设计从来不是单点优化而是一个涉及机械结构、电机设计、电力电子变换器的多学科耦合问题。在Simulink里仿真飞轮储能本质上是把这套多学科耦合系统拆解成几个子模块飞轮本身的机械运动方程、电机的电磁转矩输出、逆变器/整流器的电能变换环节、以及连接电网或负载的系统级接口。我选择永磁同步电机作为飞轮驱动电机是因为它具备高功率密度、高效率和无励磁损耗的天然优势非常适合需要在充电、放电两个方向频繁切换且转速范围变化大的飞轮储能场景。1.2 为什么选永磁同步电机而不是异步电机或开关磁阻电机在飞轮储能电机选型上我见过有人用异步电机也有人用开关磁阻电机但我个人更推荐永磁同步电机。异步电机的转子存在转差损耗在高转速工况下转子铜损和铁损都比较明显效率会打折扣。开关磁阻电机虽然结构简单、耐高温但它的转矩脉动较大对飞轮这种要求稳定转矩输出的场景不太友好而且噪声问题突出。永磁同步电机的好处在于转子侧不需要励磁电流所以没有转子铜耗配合矢量控制之后电磁转矩的输出非常平滑。飞轮储能系统恰恰需要电机在充电时提供稳定加速转矩、在放电时提供稳定制动转矩这种高动态响应的要求和永磁同步电机的控制性能非常契合。在Simulink里PMSM模块自带dq坐标系下的电压方程和磁链方程建模精度足够覆盖飞轮储能系统的研究需求不需要自己写电机内部方程的S-function能省不少时间。另外一点永磁同步电机在弱磁控制下可以拓宽转速范围这一点对飞轮储能来说非常实用。飞轮从低转速冲到高转速如果全程采用恒转矩控制转速到了一定程度之后电压就会达到逆变器上限这时候需要弱磁控制来继续提升转速。用PMSM的好处就是弱磁控制的仿真实现相对成熟模型里可以直接沿用电动车主驱的弱磁控制思路调试起来不至于两眼一摸黑。1.3 Simulink建模仿真的总体架构规划我搭的这套仿真模型整体划分为四层机械层、电气层、控制层和系统调度层。机械层就是飞轮的本体运动方程包括转动惯量、阻尼系数电气层是永磁同步电机本身和三相逆变器控制层是核心包含转速外环、电流内环的级联控制结构以及SVPWM调制算法系统调度层则是模拟飞轮储能系统的“充电/保持/放电”三种工作模式切换逻辑这层决定了整个系统在仿真中按怎样的时序运行。这样分层的理由很直接调试的时候能够分层排查问题。模型发散的时候先看是机械层参数错乱还是电气层出了问题波形不对的时候先确认控制层的给定和反馈是否正确。如果所有模块混在一起搭成大杂烩模型后期排查几乎无从下手。我在Simulink里还养成了一个习惯——每个子模块旁边都加一个写有核心公式的注释框这能帮助自己在三个月后回头看模型的时候还能秒懂当初的设计意图。2. 永磁同步电机建模与坐标变换核心原理2.1 三相静止坐标系到dq旋转坐标系的数学基础要理解永磁同步电机在Simulink里的模型行为坐标系变换是绕不开的第一道坎。很多人一开始接触PMSM控制时被Clark变换和Park变换的公式吓到其实从物理意义上理解并不复杂。三相静止坐标系下的电机方程里互感系数是随转子位置变化的直接处理的话模型是时变的控制设计会非常困难。Clark变换先把三相变为两相静止坐标系αβ坐标系相当于把三相绕组等效成互相垂直的两套绕组这解决了一部分耦合问题但αβ坐标系下反电动势仍然和转子位置相关。Park变换的核心价值在于把αβ静止坐标系旋转到与转子同步旋转的dq坐标系让基波分量在这个坐标系下变成直流量。这样电感参数不再随位置变化电流环就可以设计成线性时不变系统用常规的PI控制器来处理。在Simulink里PMSM内部模块已经封装好了计算过程但从控制角度你需要明白你输给逆变器的电压 uq、ud 并不是直接的物理电压而是dq坐标系下的合成电压控制量PWM模块负责把这个控制量反变换回三相电压。我在仿真中设置转子初始角度为0同时在控制端做了转子角度补偿对齐这样d轴电流和q轴电流的物理含义就是对应的励磁分量和转矩分量。由于飞轮储能电机基本运行在额定转速以下我使用的是Id0的转子磁场定向控制让全部电流都用来产生转矩这既简化了控制又保证了效率。2.2 PMSM模块在Simulink中的配置要点解析Simulink里的PMSM模块有几种实现形式最常用的是Simscape Electrical中的永磁同步电机模块它包含了电气端口和机械端口可以直接和电力电子电路连接。我建议不要贪图默认参数直接开跑哪怕只是仿真研究也要把电机参数设置成与实际飞轮系统量级匹配的值。比如飞轮储能电机的极对数通常是2到4对极额定转速在8000到20000rpm之间而定子电阻和d、q轴电感则决定了电流环的带宽范围。我的做法是先根据目标储能容量反推出电机功率等级再依据现有商用高速电机产品的参数表选择一组参考值。这个“贴着现实建模”的习惯很重要如果随意填一组参数仿真波形虽然也能跑出来但控制器的PI参数、逆变器直流母线电压等设计全都失去了工程参考价值。仿真可以用来验证概念但目标是让仿真结果能够指导实际系统设计而不仅仅是画几个漂亮的波形图。在Simulink里我特别关注PMSM模块的“反电动势常数”和“永磁体磁链”这两个参数的统一。这两个参数本质上是同一个物理量在不同单位制下的表达设置冲突会直接导致转矩计算错误。我实际踩过一次坑——磁链用的单位是Wb反电动势常数却按V/krpm填结果转矩常数差了近10倍仿真出来的转速响应曲线完全不对。2.3 逆变器与SVPWM调制策略的建模选择逆变器部分我使用了Simscape Electrical的三相两电平桥臂直流侧用理想直流源模拟母线电容简化了电容动态。调制策略上没有采用简单的SPWM正弦脉宽调制而是选用了SVPWM空间矢量脉宽调制这是因为SVPWM的直流母线电压利用率比SPWM高约15%并且能够更自然地处理PMSM控制输出的dq电压合成矢量。SVPWM的Simulink实现路径比较标准先根据dq轴电压和转子角度计算αβ轴电压再计算参考电压矢量幅值和角度然后判断所在扇区计算相邻矢量的作用时间最终生成逆变器开关管的PWM脉冲。这一套逻辑在Matlab里可以用函数模块或Stateflow来实现但我推荐使用Matlab Function块代码清晰且便于修改扇区判断阈值或死区补偿参数。我在实际搭建中发现SVPWM模块最容易被忽视的是载波频率与采样周期的匹配。如果控制器的采样周期是1e-4秒也就是10kHz的控制频率那么PWM载波频率也不能设置得太低否则电流谐波会明显增大。我最终将开关频率设为10kHz和电流采样、控制计算周期保持一致这既方便同步采样也减少了电流纹波对转速控制精度的影响。3. 飞轮本体建模与系统参数匹配计算3.1 转动惯量选取与储能容量换算飞轮的转动惯量是连接机械系统和电气系统最关键的一个参数。转动惯量确定之后电机可以在多长时间内把飞轮从静止加速到额定转速也就意味着充电功率在这个时间段内的大小。设计的时候不能只追求“储得多”还要考虑电机能否提供足够的加速转矩否则配置出来的系统就是一个“大飞轮配小电机”充电过程会非常漫长。我在设计参数时按照这样的顺序来推导先定义目标储能容量比如设计一个100Wh级别的实验样机模型然后明确转速范围比如从最低转速10000rpm到最高转速30000rpm再根据能量公式反推出需要的转动惯量 J 2 × E /ω_max² − ω_min²。算下来需要的转动惯量非常小因为转速高平方效应明显。这个计算过程提示了一个工程上的关键判断飞轮储能虽然能量密度低但功率密度高它更适合短时高频次的充放电调节场景而不是长时间的能量备用。在做Simulink仿真时如果以电池储能的容量思维配置飞轮参数模型里的充电时间会非常长仿真也会很无聊。合理的做法是在仿真里凸显飞轮的“功率型”特性比如用作电网调频或者不间断电源的功率缓冲。3.2 阻尼系数与摩擦损耗对仿真模型的影响机械阻尼系数在飞轮仿真里对应的是轴承摩擦和风阻损耗这个参数的实际量级非常小但如果你在模型里把这个参数设置为0系统的单位阶跃响应就会永远不衰减仿真中的转速会一直漂移。实际飞轮储能系统对真空腔体的依赖非常强因为空气阻力与转速的平方成正比在高转速下风阻损耗会急剧放大所以实际的飞轮储能装置基本都运行在真空环境下。在Simulink机械建模时我把阻尼系数设置为一个与转速相关的变量而不是一个常数。具体做法是用一个查表模块根据转速范围分段对应不同的阻尼系数低转速阶段以轴承摩擦为主阻尼系数较小高转速阶段风阻贡献显著增加阻尼系数变大。虽然这会增加一点模型复杂度但仿真得到的自放电时间常数更接近真实系统对研究“待机保持模式”下的能量损失速率非常有帮助。3.3 充放电工作模式的数学描述飞轮储能系统有两种基本工作模式充电模式下电机工作在电动状态从电网吸收电能输出加速转矩飞轮转速上升放电模式下电机工作在发电状态飞轮减速电机把动能转化为电能回馈到直流母线或电网。在Simulink里的实现方式有两种一种是通过切换电机的负载转矩方向来实现自动切换另一种是通过改变转速环的给定值来实现主动控制。我采用的是后者——给定转速求导的方式。在充电阶段给定转速是一条斜率明确的斜坡曲线转速环不断跟踪上升的参考值在放电阶段给定转速下降到目标下限电机自然进入发电状态。这个方式的优势在于模式切换不需要修改硬件拓扑或控制结构只需要修改参考转速曲线的逻辑。仿真的时候可以把三段逻辑封装在一个调度模块里一点都不复杂而且不会因为硬切换导致电流冲击。4. 双闭环控制策略设计与PI参数整定4.1 转速环与电流环的带宽分配原则飞轮储能电机控制采用典型的转速外环、电流内环级联结构。两个环路的带宽选择是有原则的电流内环需要比转速外环快得多一般差5到10倍。这个差别的逻辑在于电流环控制的是电磁转矩转矩的快速响应是转速平稳的前提如果两个环路带宽接近转速波动会导致电流指令剧烈变动最终可能引发系统振荡。我在仿真里把电流内环的带宽设置在500Hz左右转速外环带宽在50Hz左右。这个搭配让电流响应在1ms内跟上指令而转速的调整过程在几十毫秒量级。调试时先孤立调内环把转速环断开、给定固定电流观察电流跟随情况然后再把外环接上调整转速PI参数。这种“先内后外”的调试顺序是整个控制调参过程中最省时间的做法。4.2 基于工程设计法的PI参数计算过程很多初学者在设置PI参数时习惯“试凑”但试凑法是低效且容易让人沮丧的。我采用的是经典的工程设计法对电流环把电机定子电阻和电感作为一阶惯性环节来处理要求闭环带宽为 ωc那么PI参数可以按如下方式近似计算电流环在d轴和q轴的模型是一致的传递函数近似为 G(s) 1 /Ls R。闭环期望带宽为 ωc 时PI控制器的比例系数 Kp L × ωc积分系数 Ki R × ωc。可能有人觉得这个公式太简单但它确实能给出一个非常接近收敛的初始值剩下只需要微调。以我用的电机参数为例定子电阻R0.05Ωd、q轴电感L0.1mH。设定电流环带宽500Hz也就是约3141.6rad/s。计算得到 Kp 0.0001 × 3141.6 0.31416Ki 0.05 × 3141.6 157.08。把这两组值输入到PI控制器里电流环就能比较快地跟随给定。这个计算过程在Simulink的离散控制器里依然适用。需要注意的是工程设计法的前提是忽略反电动势耦合和延迟效应。实际系统里dq轴之间存在交叉耦合电压项当转速升高时这项耦合的影响会变大。解决的办法是在电流环输出端加上解耦补偿项把 ωLdiq 和 ωLqid 这样的交叉项前馈补偿掉这样在高速区域电流环仍然能保持稳定表现。4.3 转速环参数整定与运行效果转速环的设计以飞轮机械运动方程为基础J × dω/dt Te − TL − Bω其中 Te 是电机电磁转矩TL 是负载转矩B 是阻尼系数。把电流环等效成一个快速响应环节之后转速环的对象近似为一个积分环节加一个惯性环节。转速环控制器的参数直接从期望的闭环系统时间常数来定比例增益决定了转速环的响应刚度积分项则用来消除稳态速度误差。在仿真测试中我把给定转速从0渐变到20000rpm转速环跟踪过程中最大超调量控制在5%以内调整时间在0.6秒左右。这个响应速度对于飞轮储能系统的充电阶段是合理的。太快了会导致充电功率冲击过大太慢了开机体验会很难受。后面做放电模式的时候转速从20000rpm降到5000rpm电流环顺利进入负转矩区域发电功率平稳反馈到直流母线整个切换过程没有出现电流突变。如果你在调试中遇到转速超调很大的情况大概率是转速环比例增益偏大。这时候不要单纯降低Kp另一个思路是给转速参考信号加一个斜坡限幅让飞轮的加速过程不是线性阶跃而是一个平滑爬坡过程。这既保护了机械系统不被过大的加速度冲击也让充电电流保持在一个相对平稳的范围从系统层面看是更合理的控制方式。5. 充放电模式切换与系统级仿真实战5.1 双闭环控制在充电、保持、放电三种模式下的行为对比我在系统调度层面设置了完整的“三段式”运行曲线从0到0.5秒系统处于充电模式给定转速从0加速到20000rpm0.5秒到1.0秒进入保持模式转速给定保持不变考察系统的转速跌落和功率损耗情况1.0秒到1.5秒进入放电模式转速从20000rpm降到5000rpm能量回馈到直流母线。充电阶段电机电磁转矩保持正值dq轴的id分量维持在0附近iq分量稳定输出正转矩。保持阶段转速达到给定值后转矩迅速衰减到接近零只剩下用来克服阻尼损耗的微小正转矩维持转速。放电阶段转速环输出反转电磁转矩变为负值能量从飞轮侧流向电气侧直流母线电压由DC/DC变换器维持稳定。实际仿真波形验证下来转速曲线、电磁转矩曲线、dq轴电流曲线和行为预期完全吻合。特别值得说的是放电阶段电流的平滑性由于转速偏差是逐步累积的参考转速是斜线下降电流指令不会突然跳变生成的PWM占空比变化也比较平滑这在模拟真实系统运行时是很有参考价值的。5.2 储能系统进出功率计算与效率分析飞轮储能系统的瞬时充电功率可以用直流母线侧的功率来度量P Udc × Idc。在充电阶段直流侧输入功率为正值一部分转化为飞轮的动能增量一部分损耗在电机铜耗、铁耗、逆变器开关损耗和机械摩擦损耗上。放电阶段直流侧功率变为负值回馈功率小于飞轮释放的机械功率两者之差就是往返损耗。我在仿真模型里增加了几个方便数据分析的测量点电机输入电功率、输出机械功率、飞轮动能变化率、直流母线功率。把这些信号送到Simulink的Scope和Data Inspector里对比能清楚地看到能量流的方向和大小。实际跑下来系统从电网到飞轮动能的单向转换效率大概在85%左右这个数值和文献报道的高速永磁电机飞轮储能系统效率范围基本吻合说明模型参数设置是合理的。有个容易被忽略的细节需要提一下放电电效率通常低于充电效率因为放电阶段的电流方向反转后逆变器的导通损耗和开关损耗分布发生了改变。在做系统级效率评估时不能简单地把往返总效率除以2当作单向效率最好在模型里把充电和放电两个阶段的损耗分别统计这样才能准确评估系统的热负荷和电气设计余量。6. 仿真调试中常见问题与排查套路6.1 模型发散“仿真跑着跑着就去NaN了”仿真发散这个问题我太有经验了。Simulink模型发散的原因通常有四个参数设置差几个数量级、控制器的采样时间设置不当、代数环问题、以及积分步长过大。在我这个模型里最常见的发散原因是PWM开关频率和控制采样频率不匹配导致电流反馈信号出现混叠效应电流环误判方向最后正反馈失控。排查发散问题的标准流程是先看哪个信号最先变成Inf或NaN。把Scope里所有关键信号都拉出来找一个时间点最早的异常往回倒推是哪一步计算出了问题。如果电流最先异常大概率是电感值或PI增益计算错误如果转速最先异常多半是转动惯量或阻尼参数数量级不对如果是PWM模块输出异常要检查载波频率与采样步长的匹配。我强烈建议在模型里加几个Saturation饱和模块作为保护转速输出饱和、电流指令饱和、占空比饱和。这些保护模块不会影响正常运行但能在异常情况下阻止模型彻底崩掉相当于给仿真模型上了个保险丝。调试的时候有饱和限制和没饱和限制体验完全不一样。6.2 电流环跟踪滞后严重甚至振荡电流环跟踪滞后通常表现为给定iq指令上升时实际电流响应慢且出现明显的振荡超调。这时要回到电流环PI参数上找问题。首先检查采样时间设置Simulink中的离散控制器采样时间要和实际控制系统一致我用的采样时间是1e-4秒对应10kHz控制频率。其次检查dq轴是否实现了正确的解耦。在PMSM矢量化控制中当转速超过一定值后dq轴交叉耦合电压ωLdiq和ωLqid的量级会接近甚至超过PI控制器的输出电压能力。如果不做前馈解耦电流环带宽就会受限动态响应变差。很多跑高速电机控制的朋友遇到电流环不稳定第一时间怀疑PI增益但漏掉了解耦补偿的问题结果怎么调都调不稳。6.3 飞轮转速保持不住保持模式下的转速跌落保持模式下转速跌落的核心原因是阻尼损耗和系统损耗的存在。理论上如果阻尼系数为0飞轮转速可以在保持阶段纹丝不动但真实系统里阻尼损耗始终存在所以转速会缓慢下降。在仿真模型里如果发现保持阶段转速掉得比预期快需要回头查一下阻尼系数的设置是否偏大或者机械端口是否意外连接了一个额外的负载转矩。解决这个问题有两个思路一个是提高阻尼模型的精度把真空环境下的风阻损耗系数设在比较小的量级另一个是在控制策略上做补偿保持模式下让电流环输出一个小幅度的正转矩来平衡阻尼损耗维持转速稳定。第二种思路更接近实际工程的做法因为飞轮系统无论如何都会有损耗完全不做补偿的话转速必然会持续衰减。6.4 常见问题速查表现象可能原因排查步骤与解决策略仿真直接发散或NaN参数数量级错误、采样时间不匹配、代数环检查最先异常的信号逐级回溯加入饱和保护模块电流响应振荡明显电流环PI增益过大、dq耦合未解耦按工程设计法重算PI参数增加前馈解耦项转速跟踪滞后转速环带宽过低、电流环带宽不够先提升电流环带宽再调整转速环Kp放电时直流侧电流冲击大转速参考曲线变化过剧烈对转速参考加斜坡限幅提高电流环响应速度保持阶段转速持续跌落阻尼系数偏大或缺少转速维持控制改进阻尼模型增加零转速差转矩补偿SVPWM输出缺相扇区判断逻辑错误、死区时间配置异常对比三相PWM波形检查SVPWM的扇区计算代码7. 仿真加速、参数辨识与代码生成方向扩展7.1 固定步长与变步长求解器的选用逻辑做电力电子仿真时求解器的选择直接影响仿真速度和稳定性。我刚开始跑这个模型时用了默认的变步长求解器跑是能跑但每次到PWM开关瞬间步长都会自动缩小很多仿真时间长到让人怀疑人生。后来把模型改为固定步长离散求解器步长设为1e-6秒仿真速度明显提升。在Simulink里跑PWM逆变器这类强非线性系统推荐使用离散固定步长求解器如ode3或离散模型并且把PWM采样逻辑改为“基于时间的均匀采样”。这样既避免了变步长求解器在开关边沿反复缩小步长导致的效率问题也让模型的数值行为和实际数字控制系统更接近。需要注意的是固定步长太小仿真时间太长步长太大会导致PWM波形失真需要根据开关频率做权衡。7.2 外部模式实时仿真与代码生成扩展如果后续想从纯仿真走向硬件在环测试Simulink的外部模式和代码生成功能可以无缝衔接这个飞轮储能模型。利用Simulink Coder可以把控制算法部分转速环、电流环、SVPWM生成嵌入式的C代码部署到单片机或DSP上运行电机和逆变器部分仍然在PC端仿真或硬件在环设备上运行。这个“分拆部署”的思路是很多电机控制开发项目从仿真走向实物的标准路径。我在实际项目中用过一次这种流程遇到最多的问题是在代码生成阶段处理SVPWM扇区判断的浮点运算。单片机尤其是低成本MCU的浮点运算能力有限工程师往往需要在代码生成之后把浮点运算改成定点运算。如果在Simulink模型里一开始就使用定点数据类型代码生成之后的资源占用和运行效率会好得多。做仿真研究的时候可以不考虑这些细节但如果是做产品化的嵌入式控制代码生成这是无法回避的现实问题。7.3 基于仿真模型的电机参数辨识与自适应控制最后聊一个近期热度很高的方向——永磁同步电机参数辨识。在飞轮储能系统中永磁体磁链会随温度变化定子电阻随温度升高而增大导致控制器参数偏离设计值。既然我们已经搭好了PMSM的仿真模型完全可以利用Simulink的Parameter Estimation工具箱把模型输出和参考数据做对比反推出定子电阻、电感、磁链等参数的真实变化情况。我在模型里预留了参数辨识的接口把三相电流和转子角度记录下来离线运行参数估计工具就可以实时获取电机参数变化曲线。这种离线辨识的结果可以更新到控制器的PI参数计算中形成“参数辨识—PI参数在线更新—控制性能优化”的闭环逻辑。虽然这个闭环在Simulink里完全实现需要一些编程工作但相比直接做硬件实验试错成本低得多也安全得多。从我的实际使用体验来看Simulink的价值不只是画几条漂亮的波形而是通过分层建模、系数推导和问题排查把整个飞轮储能系统的物理逻辑和控制逻辑彻底打通。最初我做这个仿真只是为了完成一个课程设计但后来发现这套模型变成了一个可以反复扩展的基础平台——换一组电机参数就是另一种容量的飞轮系统改一段调度逻辑就是另一种应用场景。我个人体会最深的一点是建模仿真最大的收获不是模型本身跑通了而是你在建模、调试、改参数过程中逐步建立起来的对系统整体行为的直觉判断力。这种直觉是看多少篇论文都替代不了的。