1. 从飞轮储能到Simulink建模为什么这个仿真值得认真做这两年做微电网、新能源并网、电梯节能这些项目绕不开一个词——储能。电池储能大家最熟但有一类机械储能经常被低估飞轮储能。它本质上是把一个很大惯量的转子带到高速让能量以旋转动能的形式存起来需要放能时再让电机反过来当发电机用把动能送回电网或负载。我在前面某次某实验室项目里做的就是飞轮储能系统的建模与Simulink仿真驱动电机选用永磁同步电机这句话看着挺长拆开就是三个关键部分储能本体飞轮、机电能量转换接口PMSM电机、以及控制系统也就是让电机既能电动又能发电的那套逻辑。很多刚接触的人会问既然有电池为什么还要飞轮答案在于功率密度和循环寿命。电池适合大容量、慢充慢放的场景但像电网调频、电梯回馈制动、列车再生制动、UPS瞬时功率补偿这类功率冲击大、充放电频繁、时间短几秒到几十秒的场景电池反而吃力频繁大倍率充放电会严重缩短寿命。飞轮对循环次数几乎不敏感几万次甚至几十万次充放都没问题而且响应速度快毫秒级就能从充电切换成放电。永磁同步电机在飞轮储能里几乎是标配因为效率高、功率密度高、响应快而且可以四象限运行——正转电动把飞轮加速储能正转发电让飞轮减速放能还能在转速方向上灵活切换。这篇内容我就按实际做仿真时的思路来写先讲飞轮储能系统从机械到电气的模型怎么搭再讲Simulink里模块怎么选、怎么连线、参数怎么设接着讲控制策略里面最关键的两个环路怎么设计最后把我在调试中踩过的坑和验证手段一起摆出来。适合正在做电力电子、电机控制、储能系统仿真课题的同学以及想快速建立飞轮储能仿真平台的项目工程师。我默认你用过Simulink基础模块但不会假设你熟悉电机建模所以中间必要的理论基础我会尽量用大白话讲清楚。2. 飞轮储能系统的本体建模机械侧先搞清楚电机侧才不会乱2.1 飞轮储能的能量关系与转速范围设计飞轮存储的能量公式很简单E1/2Jω²。J是飞轮转动惯量单位kg·m²ω是角速度单位rad/s。从这里就能看出两个设计方向一是加大转动惯量J堆质量、堆半径二是提高转速ω因为能量和转速的平方成正比所以把转速从10000rpm提到20000rpm储能直接变四倍这也是现代飞轮都往高速走的根本原因。但转速不能无限抬高受两个因素限制材料强度离心力会把飞轮拉裂和轴承/电机的损耗。实际仿真里我们一般会设定一个最低转速ω_min和最高转速ω_max。为什么不放到零转速因为低速时电机反电动势太低母线电压很难稳住而且飞轮在极低速时输出功率能力很差。通常放电的下限转速是最高转速的一半左右这样能量的可用深度是1-(1/2)²75%即最多能放出总储存能量的75%。仿真模型里飞轮本体就是一个积分环节加一个一阶惯性环节。转动惯量J直接决定了转速变化的快慢dω/dt (T_m - T_load - B*ω) / JT_m是电机电磁转矩T_load是负载或扰动转矩在纯仿真里可以设成0或一个阶跃值B是粘滞摩擦系数。这里有个容易忽略的点Bω是线性摩擦模型但在高速飞轮里风阻损耗其实和转速的平方更接近风速高时湍流主导精确模型至少要加一项Cω²。我在仿真里做的是把摩擦损耗分为线性项Bω和非线性项Cω²这样高速区间的转速跌落曲线更接近实际。2.2 永磁同步电机的dq轴模型从三相世界到旋转坐标系的降维PMSM在仿真里有两种常见做法直接用Simulink自带的永磁同步电机模块或者自己用数学方程搭一个。对于飞轮储能这种需要深度控制策略研究的项目我建议一开始就用自带的PMSM模块把精力放在外环控制和参数匹配上。但如果想在论文里写清楚原理自己推导dq模型是基本功。dq轴模型的核心思想把三相静止坐标系下的交流量变换到与转子磁场同步旋转的dq坐标系下这样本来随时间变化的正弦量就变成了直流分量控制起来就像控制直流电机一样简单。变换过程涉及Park变换和Clarke变换具体公式这里不展开但最终得到的电压方程就是U_d R_sI_d L_ddI_d/dt - ω_eL_qI_qU_q R_sI_q L_qdI_q/dt ω_eL_dI_d ω_e*ψ_f其中R_s是定子电阻L_d、L_q是d轴和q轴电感ψ_f是永磁体磁链ω_e是电角速度。电磁转矩的表达式也简洁T_e 1.5p[ψ_f*I_q (L_d - L_q)I_dI_q]对于表贴式永磁同步电机L_d约等于L_q所以转矩主要靠q轴电流I_q控制I_d设成0就最简单。对于内置式永磁同步电机L_d 通常小于L_q还要利用磁阻转矩控制逻辑要加最大转矩电流比MTPA策略。飞轮储能电机一般高速运行内置式更常见但在做基础仿真验证时可以先用表贴式模型把框架跑通。2.3 机械轴模型与电机模型的接口转动惯量、阻尼和初始速度把飞轮本体和电机连接起来的关键参数有三个总转动惯量飞轮加电机转子、阻尼系数B、初始转速。总转动惯量不是拍脑袋定的要满足储能容量和功率需求。比如目标容量1kWh3.6MJ工作转速范围15000rpm到30000rpm那么J值从能量差公式能反推出来ΔE 0.5J(ω_max² - ω_min²)代入ω_max 30000rpm ≈ 3141.6 rad/sω_min 15000rpm ≈ 1570.8 rad/s可得J ≈ 3.6e6 / (0.5*(3141.6² - 1570.8²)) ≈ 1.02 kg·m²。这个转动惯量不算大但对高速轴来说已经是很实际的量级。在Simulink里设置PMSM模块参数时机械输入端口选Torque转矩输入把飞轮的动力学方程写在外面用积分器算转速再把转速反馈给电机模块这样自由度最大。如果直接用模块内部的机械参数也没问题但要额外在外部加一个转动惯量附加量——因为模块内部默认惯量只算电机转子不含飞轮。3. Simulink仿真的架构搭建主电路、控制电路和信号流怎么组织3.1 主电路拓扑三相逆变器、直流母线、负载电阻飞轮储能系统主电路的具体方案我先说明一下常见的就是交-直-交结构。并网运行时电网交流电经整流器变成直流电直流母线再经过三相逆变器给电机供电放电时电机发电能量通过逆变器流回直流母线再由逆变器或DC-DC变换器回馈电网。孤立微电网项目里直流母线会接电池或超级电容、光伏发电单元、直流负载飞轮就在这里承担瞬时功率缓冲角色。为了让仿真聚焦于飞轮本体和电机控制主电路简化为直流电压源模拟母线、三相两电平逆变器通用桥臂模块、PMSM。充电模式时直流母线电压不变逆变器工作在逆变状态把功率送给电机放电模式时PMSM发电回馈到直流母线母线再接一个可调的电阻负载或回馈负载消耗能量。直流母线电压选多少这要和电机反电动势匹配。电机额定转速30000rpm时的线反电动势峰值一般不能超过母线电压否则逆变器无法控制电流。以某模拟电机的参数为例反电动势系数Ke0.02V·s/rad那么30000rpm对应反电动势峰值≈0.02*3141.6≈62.8V线反电动势还需乘以√3约108.8V。所以直流母线电压取200V留足裕量。实际项目里母线电压越高电流越小损耗越低但IGBT和MOSFET耐压成本也上去了。3.2 永磁同步电机与逆变器在Simulink中的具体连接方式Simulink里要用到Simscape Electrical库或者Specialized Power Systems库。我的建议是用Specialized Power Systems现在被整合到Simscape Electrical里因为它的电机模块和电力电子模块配合更顺手控制信号可以直接用常规Simulink信号线。具体连接是这样直流电源正极接逆变器桥臂的()端负极接(-)端。逆变器的三相输出A、B、C接PMSM模块的定子三相端口。PMSM模块的机械端口通常是m端口输出转速和转角控制部分从这里采样转速。逆变器的门极信号来自PWM发生器一般是6路脉冲控制三个桥臂。逆变器模块选Three-Level Bridge或Two-Level Bridge都可以但飞轮储能高速运行开关频率一般在10kHz以上仿真时如果步长太大谐波会很明显。建议PWM开关频率设10kHz仿真步长1e-6秒左右不要超过2e-6秒否则电流波形失真严重。3.3 控制回路信号流从转速/电压环到电流环再到PWM控制部分我采用经典的双闭环结构外环可以是转速环充电工况或者直流母线电压环放电工况和母线稳压工况内环是dq轴电流环。之所以把转速环和电压环分开是因为飞轮系统在充电时监控的是转速防止超速放电时监控的是直流母线电压防止母线过压两者不能同时做被控量。这里有个小技巧外环切换的时候电流内环的计算逻辑不用变只需切换给定来源——一个多路选择器就能搞定。外环输出的是q轴电流参考值Iq_ref给定d轴电流参考Id_ref 0表贴式PMSM。电流环PI控制器的输出是dq轴电压参考Ud_ref和Uq_ref再做坐标反变换得到三相电压参考送入PWM发生器。电角度从PMSM模块的m端口输出得到后续经过反馈滤波和微分得到电角速度。电流环PI参数设计有标准方法根据电机电感和电阻以及期望的电流环带宽来确定。比如L_d0.5mHR_s0.05Ω电流环带宽取2000rad/s约318Hz那么Kp L_d × ω_bw ≈ 0.0005 × 2000 1Ki R_s × ω_bw ≈ 0.05 × 2000 100。先按这个初值放进模块再微调即可。转速环带宽一般取电流环带宽的1/5到1/10这样内外环时间尺度分离不容易震荡。4. 充电、待机和放电三种工况的仿真策略与模型参数示例4.1 充电工况恒转速控制下的加速过程飞轮充电的本质是电机电动运行把电能转化为飞轮动能。充电策略一般分两段初始阶段用最大转矩恒加速快速储能转速接近设定上限时切换成恒转速控制电流逐渐减小直到稳定。仿真里我习惯用阶跃给定来验证功率需求从0阶跃到15kW电机加速转矩基本恒定转速线性上升达到设定转速后转矩自动下降。转速给定上限建议加一个斜坡限制器不要直接给阶跃否则转速环PI输出瞬间饱和电流冲击很大。斜坡时间选0.5秒就够平滑了。参数示例某模拟电机配置额定功率20kW额定转速30000rpm额定转矩6.4N·m极对数3定子电阻0.05Ωd轴电感0.5mHq轴电感0.5mH永磁体磁链0.02Wb转动惯量含飞轮1.02kg·m²粘滞摩擦系数B0.0001N·m·s/rad用这个参数跑充电仿真30000rpm转速稳定时间大约在数秒量级具体要看转矩限幅和转动惯量。转矩限幅设6.4N·m额定值充电功率就稳定在约20kW速度从静止加速到30000rpm的理论时间 J×Δω/T_e ≈ 1.02×3141.6/6.4 ≈ 500秒。不过这是理想无摩擦时间实际会更长。看到这个数字有人会惊讶但飞轮储能确实更合适短时功率冲击而非长时间大能量吞吐——它的优势是功率密度极高以20kW功率持续几分钟能放出的能量是有限的所以储能容量算出来也就那么多。4.2 待机工况母线电压稳定与损耗补偿飞轮待机时转速维持在一个较高水平此时控制器仍需工作一方面通过电机施加微小转矩补偿摩擦损耗保持转速基本不变另一方面让变流器给直流母线提供稳定电压。待机工况的仿真重点在于观察损耗模型是否合理以及母线电压纹波是否在可接受范围。我在待机仿真里把摩擦模型拆成线性项风阻项设置B0.0001和C0.00000002高速时风阻损耗会显著增加。比如30000rpm时风阻损耗 C×ω³ ≈ 2e-8×3.1e10 ≈ 620W线性摩擦损耗 B×ω² ≈ 0.0001×9.87e6 ≈ 987W合计约1.6kW。这意味着待机时变流器需要持续输入约1.6kW来维持转速这部分就是系统自放电损耗。仿真结果可以清楚看到待机模式下逆变器输出功率稳定在1.5-1.7kW之间波形有微小脉动但平均功率非常稳定。如果只是做数学验证把B和C都设成0也行但那样就无法评估系统的长期待机能力和热设计需求。我建议至少保留线性摩擦项。提示摩擦参数B和C的取值不能凭感觉最好参考电机样本和实际测试数据。很多飞轮储能厂家的公开数据里会有空载转速衰减曲线由衰减曲线上可拟合出等效B和C。4.3 放电工况电机发电模式下的转速下降与功率输出放电是飞轮的高光时刻。电机从电动切换到发电逆变器把飞轮动能转换成电能送回直流母线母线电压由电压环稳住。放电过程中转速持续下降为了维持恒功率输出q轴电流需要随着转速下降而增大——因为功率PT_e×ω转速低了要维持同样功率就必须更大转矩也就是更大电流。仿真里放电场景我会设置成初始转速30000rpm负载功率需求阶跃为10kW持续放电10秒。那么放电结束时转速是多少用能量方程算ΔE 10kW × 10s 100kJΔE 0.5J(ω_start² - ω_end²)100000 0.51.02(3141.6² - ω_end²)解出来ω_end ≈ √(3141.6² - 200000/1.02) ≈ √(9.87e6 - 1.96e5) ≈ √9.67e6 ≈ 3110 rad/s折合约29650rpm。也就是说10kW放电10秒才让转速从30000rpm掉到29650rpm掉速很少。这就是飞轮储能的特性短时间功率冲击没问题但长时间大能量输出就力不从心了。如果把功率需求提高到100kW并且持续10秒那结束转速是ω_end √(3141.6² - 2e6/1.02) ≈ √(9.87e6 - 1.96e6) ≈ √7.91e6 ≈ 2813 rad/s ≈ 26860rpm掉到26860rpm已经发了100kJ的能量换算约0.278kWh。从中能直观看出飞轮跑长时恒功率转速和端电压都会快速下跌系统设计要根据这个特点规划可用能量窗口。4.4 放电过程中母线电压控制的难点为什么只用电压环不够放电时直流母线不能直接靠电压环控制逆变器就稳稳当当因为逆变器是电压源母线电压主要由外部源和电容支撑。当电机回馈电流很大时母线电压会泵升。要控制住母线电压通常要加制动斩波器crowbar或让回馈能量流入电池组否则母线电容会过压烧毁。我在仿真模型里加了一个母线电压PI控制器它以给定母线电压为参考如200V输出PMSM的Iq_ref。当母线电压因回馈功率升高时Iq_ref增大发电功率增大从而把电压拉回参考点。母线电容大小很关键电容越小动态电压波动越大电压环PI参数越难调。经验做法是电容值让时间常数R_C×C在0.05秒左右R_C是母线负载等效电阻。比如母线电压200V负载功率10kW等效电阻R_CU²/P4Ω电容C≈0.05/4≈0.0125F12.5mF。这个数值对仿真来说不算大但如果实际硬件里这么小纹波会很大所以仿真里经常会把电容设得比理论值大一些比如20mF让母线电压稳一些。5. 电流环参数整定与坐标变换的仿真实现细节5.1 电流内环PI参数的工程整定方法前面给的Kp和Ki初值公式是常规做法但实际整定还要考虑离散化和延迟。Simulink仿真里PWM采用离散化控制系统用BFSBackward Euler或Tustin离散化采样周期取10微秒100kHz比较合适。此时电流环PI参数可以比连续域算出来的值稍做调整因为离散化引入了半拍延迟。调试技巧先把电流环打开外环断开给一个阶跃Iq_ref10A观察Iq响应。好阶跃响应应该是无超调或小超调10%调节时间在2-3ms内。如果振荡按顺序做先减小Kp保持Ki不变看到波形稳定后再逐步加大Ki直到出现振荡再退10%。这样操作比直接套公式要快得多尤其模型参数与铭牌参数有偏差时。5.2 Park变换和Clarke变换在Simulink里的实现方式坐标变换我推荐直接用Simulink自带的abc_to_dq0 Transformation模块输入三相电流和电角度theta输出Id、Iq。该模块的变换公式可以有不同变体等幅变换/等功率变换仿真里默认是等幅变换和PMSM模块内部定义一致即可。需要特别谨慎的是电角度采样。PMSM模块的输出端通常给的是转子位置机械角度需要乘以极对数p得到电角度theta_e。如果把机械角度直接送进Park变换Id和Iq就不是真正的直流分量环PI会失效。这一点是新手最容易忽略的坑我见过不止一次因为少了这个p乘积导致电流环完全发散。5.3 离散求解器选择对仿真结果的影响飞轮储能模型是典型的高频开关电气模型加中频机电模型混合系统。Simulink求解器我推荐用离散求解器如ode4或ode45固定步长。步长选择逻辑开关频率10kHz对应周期100微秒步长至少小于周期的1/10即10微秒但为了PWM脉冲边缘的分辨率我通常取2微秒到5微秒。如果电机功率等级小、可接受慢一点仿真速度取5微秒即可功率大、谐波多建议2微秒。还有一个大坑连续求解器变步长在Simulink里碰到电力电子开关时虽然能自动减小步长但会带来大量BDF/NDF积分器开销仿真速度反而很慢且开关噪声会让步长不断震荡导致仿真卡死。所以对于逆变器加PMSM这种系统固定步长离散求解器是稳定和速度兼顾的选择。6. 仿真实验设计与结果分析怎么验证模型正确性6.1 从能量守恒角度验证模型的正确性仿真模型不是搭完能跑就行必须验证物理正确性。我习惯用能量守恒来做第一道检查。把整个系统框起来统计电气输入能量、机械动能变化、摩擦损耗三者应满足输入电能 飞轮动能变化 摩擦损耗 输出电能在Simulink里可以用积分模块分别对功率做时间积分然后对比。如果偏差超过2%-3%说明某个环节的连接或参数有误。最常见问题PMSM模块的电气功率符号定义和机械功率方向不一致导致能量路径计算错误。6.2 电压、电流、转速、转矩波形的观测与判读波形观察要注意分量和归一化。通常在示波器Scope里同时放转速rpm或rad/s和时间曲线充电阶段应近似线性上升放电阶段应缓慢下降。d轴和q轴电流稳态时Id≈0Iq为正充电或负值放电。电磁转矩与Iq成比例充电时为正驱动加速放电时为负制动减速。三相电流开关纹波叠加在基波上频率应接近电频率转速乘以极对数。直流母线电压应有较小纹波不应有明显的电压跌落或过冲。比如放电工况三相电流频率 转速(rpm)/60×极对数p。初始30000rpm、p3时三相电流基波频率30000/60×31500Hz。如果你看到的三相电流频率不对往往说明坐标变换角度或极对数设置错了。6.3 转矩限幅、电流限幅和超速保护逻辑的仿真飞轮储能系统安全性是重中之重仿真里也应该体现保护逻辑。我给模型加了三重保护电流限幅Iq_ref经饱和模块限幅到±20A防止过流。转矩限幅电磁转矩限制在±6.4N·m防止机械过载。超速保护转速超过31000rpm高于额定上限1000rpm时触发PWM封锁逆变器六路脉冲全部置零电机进入自由停车damping由摩擦承担。在仿真里验证超速保护逻辑可以把转速给定突然提高或者把转矩限幅调大人为制造超速环境观察PWM是否被及时封锁。这些保护逻辑虽然简单但要从一开始就放在模型里否则后期再加容易遗漏关键路径。7. 经验教训与易出错环节调试飞轮储能仿真时那些让人抓狂的问题7.1 逆变器与电机接口的电气基础问题Simulink里如果遇到Simulink cannot determine block output...之类的报错先检查逆变器门极信号类型是Simulink逻辑信号0/1还是电气信号。通常PWM发生器输出的是Simulink信号可以直接接入逆变器模块的Gate端口。但如果用了Specialized Power Systems的变流器模块门极输入可能需要0/1的脉冲信号注意信号类型不匹配会导致仿真直接报错。另外直流电源和逆变器之间建议加一个大的母线电容如20mF和一个小电感如0.1mH模拟实际母线的分布电感和限流作用。不加母线电容的话逆变器开关瞬间的电流冲击直接加载到直流源上容易导致电压跌落和数值振荡。7.2 转子初始角度对启动仿真的影响PMSM模型默认初始位置一般为0度但如果从静止0rpm启动且施加的是Iq0的控制策略则转矩为0转子不会自己转起来。永磁同步电机不像异步电机有自启动能力这一点必须注意。我把启动阶段控制策略改成速度开环切换给定一个固定频率的正弦电压或直接给定一个很小的Iq让电机转动起来后再切换为闭环控制。这种开环启动方式简单说就是在仿真里先让电机转起来再切入闭环避免静止状态下的控制死锁。如果你从静止直接给Iq_ref一个阶跃也会正常转动但初始状态如果Id和Iq都是0且初始位置在某个特定角度可能因为无刷直流式的换相问题导致启动失败。最稳妥方式是在仿真开始前给一个很小的初始转速比如1rpm避免数值上的奇异。7.3 电流环和转速环带宽配合不当导致的震荡这个坑非常典型。如果你把转速环带宽设得太高比如和电流环带宽同一量级两个环互相耦合系统容易出现持续振荡。我从调试经验里提炼一个准则转速环的穿越频率建议控制在电流环的1/10以下保持两环时间尺度分离否则跟父子吵架一样谁都无法稳定。用波特图核实电流环闭环带宽约2000rad/s时转速环带宽设在200rad/s以内。在Simulink里可以对线性化模型用Control System Toolbox做bode分析但大多数时候用阶跃响应试凑也行。如果出现振荡优先降转速环Kp而不是降Ki——Ki主要影响稳态误差和抗扰Kp影响相位裕度振荡多半是Kp太大。7.4 仿真速度太慢时的策略模型简化与并行加速飞轮储能系统模型里最耗时的部分是PWM高频开关和逆变器电力电子模块。如果只是研究控制策略而非开关细节可以改用平均模型用可控电压源替代逆变器其输出电压由PWM占空比和母线电压决定不再有开关纹波。平均模型可以把仿真速度提升几十倍且不影响外环动态特性的准确性。我在做长时间几百秒充放电工况仿真时几乎都用平均模型只有做低电平谐波分析时才切换回详细开关模型。另一个技巧是分层仿真先做电气快速动态步长几微秒监测波形和响应再做外环慢动态平均模型验证整体能量管理逻辑。两个模型的参数保持一致结论可以互相对照。8. 从仿真到实物哪些仿真结论靠得住哪些要小心仿真和实物之间永远有一条鸿沟飞轮储能尤其明显。仿真里最可信的是控制逻辑、能量关系、PI参数初值最不可信的是损耗模型、温升以及机械振动。比如我在仿真里设置的风阻系数C换算到实际飞轮系统要经过CFD或实测才能标定否则待机损耗的绝对数值一定会偏。飞轮储能的精确损耗数据通常来源于真空腔体的真空度、转子表面粗糙度、以及轴承摩擦系数这些很难从Simulink里凭空得到。我曾经参与过一个某模拟飞轮项目仿真阶段算出来的待机损耗是1.2kW实际台架测试是1.8kW差了50%。原因就是高速轴承的交叉耦合损耗没法用简单BC模型描述。所以仿真结果的价值更多在于系统设计边界验证和控制策略可行性而不在绝对效率预测。提示如果你在毕设或项目中需要给出可靠的效率数据请使用厂商提供的实测损耗曲线或者实测空载转速衰减曲线反推摩擦模型不要直接把仿真里的B和C当真实值。9. 后续可以扩展的方向飞轮储能系统仿真平台做到这一步已经是一个相对完整的闭环。在此基础上还能扩展几个方向一是加入磁悬浮轴承或主动磁轴承模型研究转子悬浮控制对损耗和稳定性的影响二是加电网侧变流器模型把系统连接到三相电网或微电网研究有功无功控制、频率支撑三是加入与电池的混合储能协调用飞轮扛高频分量电池扛低频分量四是做损耗敏感度分析和参数辨识找到影响系统效率的关键参数。对于准备深入做此课题的同学建议先把基础模型做扎实能量守恒验证、波形一致性验证、保护逻辑验证都通过了再往混合储能或电网互动方向走。基础模型的所有参数都应该记录下来形成一份模型参数说明表否则后面改一处就乱一处。我个人的体会是飞轮储能仿真最大的价值不是跑通而是通过建模强迫你把机械、电气、控制三个学科的知识串起来并且在Simulink里建立一条完整的能量流转认知。你看到母线电压、转速、电流这些物理量在几个仿真场景里如何联动理解会比读十篇文献都深。最后再分享一个小技巧把每个典型工况的截图和关键数据存成一个固定的实验记录表可以帮你快速向导师或团队成员解释系统行为也方便自己回查不同参数版本的差异。