你搜“飞轮储能建模”看到的往往是一堆精致的结构图、复杂的公式推导然后就没了。自己动手在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因为表贴式结构适合高速旋转转子结构对称没有磁阻转矩配合id0控制非常简单高效。如果转速要求特别高还会用内置式IPMSM配合弱磁控制。我的建议是初学建模就从表贴式PMSM开始先跑通id0控制下的充放电流程后面再去碰弱磁和内置式的复杂控制。这就像学跑步先走直线别一上来就搞弯道冲刺。1.3 建模时要把系统拆成哪四层飞轮储能系统在仿真里不是一个孤立电机而是一个完整的“机械-电磁-控制-电网”耦合系统。我的习惯是把它拆成四个层次来建模第一层机械层。核心是飞轮转子本体、转动惯量J、轴承摩擦、风阻损耗这一层决定了飞轮的转动特性。第二层电磁层。即PMSM本体包括电压方程、磁链方程、转矩方程以及从abc三相到dq两相旋转坐标系的坐标变换。第三层控制层。这是整个仿真的灵魂包括速度外环、电流内环、PI调节器、SVPWM调制、id0矢量控制以及充放电模式切换逻辑。第四层电网交互层。包括逆变器、直流母线电容、电网侧变流器以及母线电压稳定控制。这样分层的好处是每一层都能单独验证。机械层可以单独测空载转速跌落曲线电磁层可以用恒转矩负载验证电磁转矩控制层可以暂时抛开电网直接给定参考转速最后再整体联调。如果你把三层揉在一起出了问题根本分不清是控制算法的问题还是电网侧设置的问题排查起来非常痛苦。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]对于表贴式PMSMLd约等于Lq所以转矩方程简化为Te 1.5·p·ψf·iq这就是为什么id0控制那么重要——id等于零时磁链完全由永磁体提供转矩只与iq成正比转矩线性、可控、无磁阻分量控制最简单。有几个特别容易被新手忽略的点第一机械角速度和电角速度不是一回事。电角速度ωe等于机械角速度ωm乘以极对数p即ωe p·ωm。转矩公式里的速度项全部要用电角速度而机械运动方程里的转速是指机械角速度。很多人在仿真里把这两个量搞混直接导致转矩计算差了好几倍。第二ud和uq方程里的交叉耦合项即ωe·Lq·iq和ωe·Ld·id在转速较高时非常显著。如果想提升仿真精度控制器里应该加入前馈解耦项把耦合项抵消掉。否则高速时你会发现设定iq不变实际的转矩响应却随着转速漂移这就是耦合项在捣乱。2.3 仿真参数怎么取值才靠谱参数是仿真能否反映真实系统的关键乱填参数等于拿着橡皮子弹练枪法。下面是一组我实际使用的是飞轮储能仿真参数可供参考参数符号取值确定方法直流母线电压Vdc700 V按变流器容量和电机反电动势确定飞轮转动惯量J0.05 kg·m²飞轮体电机转子估算额定转速nN20000 r/min按飞轮线速度极限确定定子电阻Rs0.5 Ω电机绕组手册或实测d轴电感Ld5 mH电机手册/有限元计算q轴电感Lq5 mH表贴式电机Ld≈Lq永磁磁链ψf0.1 Wb由反电动势常数推导极对数p4电机设计参数额定功率P5 kW按应用场景确定磁链ψ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.4Kic 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.3Kis一般取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 常见问题排查速查表现象可能原因排查步骤解决办法转速发散或为负损耗模型符号处理错误检查运动方程中损耗项符号空载减速测试加sign(ω)判断损耗阻止反向转动电流快速振荡发散PI参数不匹配、离散化步长过大先单独调电流环观察dq电流跟踪降低Ki或改用连续平均模型验证充电和放电方向反了坐标变换符号约定或转矩符号错误核对dq轴正方向和Park变换修正变换矩阵或转矩公式正负号转速稳态误差大速度环积分增益不足或摩擦项干扰加大Ki观察静态误差是否减小适当提高Kis必要增加前馈摩擦补偿母线电压放电时飙升回馈功率大于负载吸收能力观察功率流动和母线电压趋势加卸荷电阻或强化电网侧母线电压外环仿真速度极慢使用了开关模型且步长过小观察占耗时模块评估是否必须开关级精度换平均模型或只在验证谐波时才用开关模型最后再分享一个小细节仿真跑通不等于参数合理一定要做能量守恒校验。把充放电全过程的能量变化和电机输出功积分对比误差超过5%就要回头看损耗模型或数值积分选项。储能系统的核心是能量能量对不上账一切结论都不可靠。我自己调试时能量闭环校验几乎每次都能揪出至少一个隐藏的参数问题。