飞轮储能听起来有点像“把一个大转盘装进电机里转着存能量”但真做起来远不是“转个圈”那么简单。如果你正在用 Simulink 搭飞轮储能系统仿真还选了永磁同步电机作为驱动电机那你大概率已经感受到真正麻烦的不是让它转起来而是让充电、保持、放电整个过程中的能量边界在模型里也能转得清清楚楚。这篇文章不是从软件自带的 Demos 里抄一份模型给你而是把整套系统按“机械层-电机层-控制层-场景层”拆开把飞轮储能系统建模与 Simulink 仿真里最容易被忽略的物理量、参数关系和调试方法摊开说透。无论你是研究生做储能课题还是在做电网调频、UPS 或再生电能回收的工程师沿着这套思路都能搭出一版能跑、能看结果、能写进报告里的仿真模型。1. 项目概述与需求拆解飞轮储能的核心思想非常简单把电能变成“旋转动能”存起来需要的时候再让电机反过来当发电机用把旋转动能重新变成电能。这种储能方式的特点是功率密度高、响应快、循环寿命不受化学衰减影响特别适合短时高频次的功率吞吐场景。你可能会问储能介质那么多锂电池、超级电容、超导线圈为什么要用飞轮在 Simulink 里做仿真时把这几种储能方式放在一起看就非常直观锂电池的能量密度高但充放电循环次数和功率响应是硬约束超级电容功率大但能量密度偏低自放电也快飞轮刚好落在两者中间它的“容量”由转动惯量和最高转速决定“功率”由电机和变流器决定两者可以独立设计。这个特性在系统级仿真里特别宝贵因为你可以只改转动惯量来调容量只改电流环限幅来调功率不用像电池那样反复改化学模型。那为什么飞轮储能系统的驱动电机一定要选永磁同步电机从仿真模型的角度看PMSM 有几个非常诱人的特点。首先它采用永磁体励磁没有转子励磁绕组电机模型里不需要额外算励磁电流这会让仿真速度更快也不会出现感应电机那样因为转差率带来的滑差损耗项。其次PMSM 的转矩与 q 轴电流近似成正比转矩控制非常直接特别适合飞轮储能这种需要明明在“充电”又要随时准备“放电”的工况。第三PMSM 在低速段也能输出稳定转矩飞轮从零转速启动时不会像异步电机那样出现启动转矩不足的问题。有些资料喜欢把飞轮电机叫做“电动/发电一体机”这个说法在 PMSM 的模型里其实是同一个东西。电机工作在电动状态时电磁转矩方向和转速方向一致能量从直流母线流向飞轮工作在发电状态时电磁转矩方向和转速方向相反能量从飞轮回馈到直流母线。换向的本质不是改变电机的结构而是改变 q 轴电流的符号。仿真里只要控制量给对模型自动会呈现这两种状态不需要像传统双轴系统那样切换动力路径。所以这本项目的需求可以拆成三层第一层是飞轮本体的能量特性和损耗特性建模第二层是 PMSM 的电气与机械动态方程建模第三层是控制策略与仿真场景设计。三层之外还要解决一个隐性需求怎么判断仿真结果是合理的而不是各个模块“各转各的”但整体能量对不上。这些都是下文要重点展开的内容。2. 系统构成与数学建模2.1 飞轮本体的能量模型飞轮储存的能量由转动惯量和转速共同决定E 0.5 × J × ω²其中 J 是飞轮与电机转子折算到同一根轴上的等效转动惯量单位是 kg·m²ω 是机械角速度单位是 rad/s。这个式子看着简单但里面藏着一个设计原则因为能量与转速的平方成正比提升最高转速对储能容量的贡献远比增大转动惯量更明显。仿真时你会观察到同样的现象转速从 6000 rpm 提升到 8000 rpm转速只增加了 33%储能量却增加了约 78%。机械层面上飞轮在运行时会有两类损耗。一类是近似与转速成正比的阻尼损耗比如空气阻力、轴承摩擦模型里可以写成 B × ω另一类是近似恒定的摩擦损耗比如密封装置和轴承的静摩擦模型里可以当成一个常量转矩 T_f。在 Simulink 里更贴近实际的做法是把机械方程写成(J_fw J_motor) × dω/dt T_e - T_load - B × ω - T_f其中 T_e 是电机的电磁转矩T_load 是外部负载转矩充电工况下 T_load 可以设为 0 或负值放电工况下则用电机转矩直接承担负载需求。仿真时如果完全不计损耗储能曲线会好看但太理想化了建议至少保留一个线性阻尼项否则后续做充电保持工况时转速会一直不变体现不出工程意义。飞轮的角速度是系统的状态变量在 Simulink 里就是那个积分器的输出。初始条件设为 0 表示从静止开始充电如果模拟飞轮已经转起来、开始放电就把初始转速设成一个非零值。这点很容易被忽略但直接决定了仿真是从“满电”开始还是从“空电”开始。2.2 永磁同步电机的数学模型PMSM 在 abc 三相坐标系下的方程是一组时变系数的微分方程直接放进 Simulink 里很难处理通常都会先做 Clark 变换和 Park 变换把模型转换到 d-q 旋转坐标系下。旋转坐标系与转子磁场同步旋转原来随时间变化的电感和反电动势项就变成常量系数了仿真会稳定很多运算量也更小。在 d-q 旋转坐标系下PMSM 的电压方程是V_d R_s × i_d L_d × di_d/dt - ω_e × L_q × i_q V_q R_s × i_q L_q × di_q/dt ω_e × (L_d × i_d λ_pm)其中 R_s 是定子电阻L_d 和 L_q 是 d、q 轴电感λ_pm 是永磁磁链ω_e 是电角速度与机械角速度的关系是 ω_e p × ω_mp 是极对数。电磁转矩方程为T_e 1.5 × p × [λ_pm × i_q (L_d - L_q) × i_d × i_q]。表面式永磁同步电机的 L_d 和 L_q 基本相等公式后半项为零转矩就简化为 T_e 1.5 × p × λ_pm × i_q。这意味着 q 轴电流和转矩之间就是纯比例关系控制起来非常清爽。这也是飞轮储能系统用表贴式 PMSM 比较多的原因之一——不需要考虑磁阻转矩的复杂分配只要把 i_d 控制在 0 附近就行。不过在 Simulink 里建模型时没必要把电流方程用积分器全部手搭一遍。如果你用的是 Simscape 电气模块库直接拖一个永磁同步电机模型设置好参数和极对数就行如果你想完全使用基本模块搭建自己可见的模型也可以把上面的微分方程拆成积分器、增益、乘法器构建一个 d-q 轴电流计算子系统再通过逆 Park 变换输出三相电压。两种路线我都试过。用 Simscape 模块的优势是搭得快适合系统级验证用基本模块手搭的优势是每个增益都能看到方便写论文中的数学模型推导也方便做故障注入。文章后面以“基本模块加控制子系统”的思路展开因为这样你能看到每一步信号是怎么流的排查问题也更直接。2.3 四象限运行与充放电模式飞轮储能系统需要支持电动和发电两种状态对应电机的四象限运行。这里的核心逻辑用一句话说就是系统根据功率需求决定 q 轴电流参考值的方向和大小。充电模式下外环速度控制器输出正的 q 轴电流参考值 i_q_ref电机输出正向电磁转矩推动飞轮加速转速上升飞轮动能增加。放电模式下系统给出负的 i_q_ref电机输出反向电磁转矩飞轮减速机械能转化为电能回馈到直流母线。保持模式下i_q_ref 只需要克服摩擦损耗和阻尼损耗维持转速基本不变这时候能量只是在弥补系统自身的损耗。在模型里充电还是放电并不需要刻意切换电机模型只需要改变 i_q_ref 的符号。这就是协调控制层的任务。实际项目里常见做法是给出一个功率参考值 P_ref然后根据当前转速折算成转矩参考值T_ref P_ref / ω_m再把 T_ref 除以转矩系数得到 i_q_ref。这种功率转转矩、转矩转电流的链路非常直观飞轮低速时同样的功率需要更大的转矩所以控制器会自动把电流限制住。有些仿真会直接把 i_q_ref 当成输入那当然更简单但工程意义弱很多。我建议至少做成功率参考输入这样你面对的就不再是“电流有多大”而是“系统需要充多少功率”这个真实工程问题。2.4 直流母线侧的能量平衡研究飞轮储能系统光看电机本身是不够的必须把直流母线也纳入仿真视野。飞轮充电时母线上其他电源需要向电机侧输送能量飞轮放电时电机侧反馈的能量会抬高母线电压。如果母线电容不够大又没有放电负载或缺电支路电压可能升到保护值这在仿真里也能看得很清楚。Simulink 里常见做法是把直流母线简化成一个带电容的电压源或者直接用受控电流源把电机功率折算到母线上。前者适合做完整的“电源-变流器-电机-飞轮”一体化仿真后者适合只看能量流向不考虑电压波动的快速验证。飞轮储能的系统级仿真里母线电压波动本身就是重要观察指标建议至少保留一个电容模型否则放电时电压变化看不出来。3. Simulink 建模流程与核心模块配置3.1 顶层架构设计一个完整的飞轮储能系统仿真模型顶层可以分成四个子系统飞轮机械子系统、PMSM 电气子系统、PWM 逆变器与直流母线子系统、矢量控制子系统。飞轮机械子系统负责把电磁转矩转化为转速核心是一个积分器输入是转矩差输出是角速度。PMSM 电气子系统负责把 d-q 轴电压转换成三相电流再算回电磁转矩。逆变器与母线子系统负责把控制输出的电压指令转换为 d-q 轴电压同时体现直流母线电容的电压变化。矢量控制子系统负责产生 PWM 电压指令并由上层的速度环/功率环给出电流参考。顶层不需要一开始就把 PWM 开关频率级细节全部建出来。先跑一个连续域的平均模型把能量流和控制逻辑调试正确再去加 SVPWM 和开关器件这样定位问题会快得多。很多仿真最后“结果不对”不是模型错而是从一开始就把开关噪声和个人调试问题混在一起看不出哪一部分的影响。推荐分两步走第一步是不含开关的平均模型控制周期和电气周期统一观察能量曲线和电流包络第二步再加入载波调制、死区和采样保持观察高频纹波。文章里主要讲第一步的系统框架因为这是飞轮储能仿真的主干。3.2 求解器与步长设置PMSM 模型里有一条容易踩坑的线电角速度 ω_e 等于机械转速乘以极对数。如果你的仿真额定转速是 8000 rpm约 838 rad/s极对数为 2那么电角速度就达到 1676 rad/s对应频率约 267 Hz。电流环的带宽必须覆盖这个频率仿真步长也必须足够小否则电流波形会出现明显的阶梯失真。对于连续模型我习惯先用 ode23t 或 ode15s 这类带刚性解的求解器因为 PMSM 电流回路的电感很小、时间常数很短普通 ode45 很容易步长反复收缩跑得非常慢。对于离散模型逆变器开关频率如果是 10 kHz仿真步长至少取 1e-5 秒到 1e-6 秒量级如果只做平均模型步长取 1e-4 秒也可以接受但前提是控制周期也要在这个尺度上。还有一点要注意Simulink 默认的“自动求解器”不适合这种电动机储能系统。自动求解器为了保证稳定会把步长压得很小仿真时间一长就会出现半天跑不完一秒钟工况的情况。手动指定固定步长或者合适的变步长求解器配合你预期的控制周期来设置是必须做的一步。3.3 矢量控制策略实现飞轮储能 PMSM 的常用控制策略是 id0 的转子磁场定向控制。它的含义是让定子电流的 d 轴分量始终为零把全部电流都用来产生 q 轴转矩。对于表面式 PMSM 来说id0 就是最优策略模型里只需要一个很小的 PI 控制器把 i_d 拉回零即可。整体控制链路是转速参考 ω_ref 与实测转速 ω_m 做差经过速度环 PI 得到 q 轴电流参考 i_q_refi_q_ref 再经过限幅与实测 i_q 做差经过电流环 PI 得到 q 轴电压参考 v_q_refd 轴电流参考设为 0与实测 i_d 做差经过电流环 PI 得到 v_d_ref。最后经过 Park 逆变换得到 α-β 轴电压送入 SVPWM 或直接用于平均模型。这里的转速参考值在充电工况是一个斜坡上升的目标函数比如从 0 线性升到 8000 rpm时长 30 秒保持工况是保持在 8000 rpm放电工况可以是一个功率参考 P_ref通过 T_ref P_ref / ω_m 折算出 i_q_ref需要注意低速时功率折算出来的转矩可能超限必须先做限幅。在一个仿真环境里同时展示充电、保持、放电三种模式建议做一个“模式选择器”输入可以是 1、0、-1分别代表充电、保持、放电。这样做的好处是可以写测试用例比如“0-30 s 充电30-31 s 切换31-50 s 放电”一次仿真跑完三段工况直接看到飞轮能量上升又下降的完整曲线。3.4 状态观测与数据记录飞轮储能仿真里最该看的指标包括飞轮转速 ω_m、电磁转矩 T_e、d-q 轴电流、直流母线电压、飞轮动能变化量、输入/输出功率。飞轮动能不要直接从转速用“0.5×J×ω×ω”反复计算最好是用积分器算能量即对功率做时间积分。因为功率本身就是电压和电流的乘积积分器输出直接就是“已充入能量”或“已放出能量”和理论公式做对比时更有说服力。用 Simulink 的 Data Inspector 记录数据时建议勾选“保存所有输出”但只把关键信号标成高亮否则模型里的中间变量太多会生成海量数据导致仿真后期卡顿。我习惯把转速、转矩、直流母线电压、动能这四路信号单独拎出来其余信号只在需要检查时才打开。4. 典型参数与仿真实例验证4.1 一个可复现的标称参数体系下面这套参数是我在仿真项目里常用的标称值是一套逻辑自洽但不需要对应任何具体真实产品的设计值。你可以原样跑也可以按比例改。参数数值说明飞轮等效转动惯量 J1.0 kg·m²包含电机转子折算惯量额定机械转速8000 rpm对应约 838 rad/sPMSM 极对数 p2决定电角速度倍数定子电阻 R_s0.02 Ω影响铜耗和电流响应d 轴电感 L_d0.3 mH表贴式可令 L_dL_qq 轴电感 L_q0.3 mH同上永磁磁链 λ_pm0.12 Wb决定转矩系数直流母线电压750 V逆变器输入侧开关频率10 kHzPWM 载波频率阻尼系数 B0.005 N·m·s/rad风阻、轴承损耗等效恒摩擦转矩 T_f0.5 N·m密封件、轴承静摩擦用这组参数算一下最大储能E_max 0.5 × 1.0 × 838²大约 351 kJ。如果平均功率做到 50 kW理论上能充满约 7 秒的容量。这个水平放在飞轮储能里属于中小型的装置但对于学术研究、系统级仿真和功能验证完全够用。转矩系数 K_t 1.5 × p × λ_pm 1.5 × 2 × 0.12 0.36 N·m/A。也就是说q 轴电流每增加 1 A电磁转矩增加 0.36 N·m。如果希望额定转矩是 100 N·m额定电流就是约 278 A对模型里的电流限幅和功率电平等参数都有指导意义。这些参数并不需要一次性全部确定。在 Simulink 里我建议通过 MATLAB 脚本定义参数变量模型里直接用变量名这样调参数时只需要改脚本不需要改模型。模型设计成“参数脚本-子系统-场景输出”三层结构比把所有数字都写死在模块里高效得多。4.2 充电工况实测体会充电工况的标准设置是仿真时长 30 s速度参考从 0 斜坡上升到 8000 rpm斜坡时间 25 s之后保持稳定。我实际跑下来的结果是启动阶段电流保护限幅立刻触发电流环把 i_q 顶在限幅值上电磁转矩保持在最大值飞轮加速最快随着转速升高反电动势逐渐增大电流环可用的电压余量减少充电功率逐渐逼近最大值。如果你看到的转速曲线是一条非常完美的直线没有出现任何加速渐缓的趋势那多半是模型里没有考虑反电动势对电流环的影响或者电流环增益太大直接忽略了实际电压约束。飞轮充电的物理特征是越到高速越难继续加速因为反电动势和功率都在上升电流却受限于逆变器容量。仿真里出现这个趋势才能说明模型真实。充电过程中最容易出错的地方是速度环 PI 参数。速度环输出是 i_q_ref如果比例增益 P 设得太大转速还没有跟上参考i_q_ref 就已经被限幅饱和后续转速到达目标值时积分器还没有退饱和就会出现明显的超调。此时加一个抗饱和机制或者把速度环 P 调小、I 调慢都比直接加限幅要好。4.3 放电工况与能量回馈验证放电工况我常用的做法是前 30 秒完成充电第 30 秒发出一个 30 kW 的放电功率指令持续 10 秒观察母线电压和飞轮转速变化。现实问题来了仿真里如果没有给放电能量安排“去处”母线电容就会被充高电压可能直接超过保护值。这时有两种选择一是把飞轮电机模型连接到一个小电阻负载让回馈电能消耗在电阻上观察母线电压受控上升二是连接回网变流器模型把能量回馈到交流电网侧但这就必须额外搭建一个网侧变流器复杂度会上升。做第一级验证时我建议先用电阻负载代替电网。母线电压从 750 V 略微抬升飞轮转速从 8000 rpm 开始下降动能曲线持续下降同时电阻上检测到放电功率。这样你就验证了“机电能量转换-母线支撑-负载消耗”这条完整链路后面再替换成真实的网侧变流器模型逻辑也不会有大变化。放电时还要特别注意电流环正负方向的对称性。充电时输出的是正 i_q放电时输出负 i_q如果限幅只做单边就会在切换瞬间出现跳变。我通常会把限幅写成对称形式i_q_ref 的限幅为 [-Imax, Imax]避免出现只能正不能负的尴尬。4.4 三次仿真结果如何确认可信判断仿真结果是否可信不能只看曲线形状。我总结了三个自检方法每次跑完都过一遍。第一能量守恒。把输入到系统的电功率积分值和飞轮动能的增加量、各种损耗的总和做对比误差应该小于 1%。如果误差很大大概率是模型中某个增益、系数或单位折算写错了。飞轮储能仿真里能量不守恒是最常见也最致命的错误。第二电流和转矩的关系。确认 T_e 与 i_q 之间满足 T_e 1.5 × p × λ_pm × i_q。用 Data Inspector 把一个点的 i_q 和 T_e 调出来手算一遍能对上说明电机模块接线正确。第三母线电压变化的方向。充电时如果母线电压应该下降或经电源稳压放电时母线上应该出现能量回馈抬升。如果方向反了说明功率方向或者电流参考符号反了这种错误经常藏在不显眼的增益符号里。5. 调试实录与常见坑5.1 代数环和初始条件错误Simulink 里把飞轮机械子系统做成“转矩输入、转速输出”的积分环节时输入端口和输出端口之间没有直接代数依赖通常不会形成代数环。但如果把摩损公式写成 T_e - B × ω 而又把同一个 ω 同时用于计算和反馈就可能出现代数环。代数环会导致仿真报错或运行极慢。解决方法是给阻尼转矩反馈路径加一个 memory 块或单位延迟阻断瞬时的代数依赖。也可以把机械方程改写成 (T_e - T_load) × 1/(J × s B) 的形式用传递函数模块代替积分器这样阻尼项就并入了传递函数分母自然避免了代数环。初始条件的问题更隐蔽。飞轮从零转速启动没问题但如果模拟的是“转速已经 5000 rpm 继续充电”或“已经 8000 rpm 开始放电”就必须在积分器模块里设置初始转速。若疏漏积分器会从默认 0 开始积分最后曲线全都对不上。排查时我会先看转速初值是否与预期工况一致再往下查电流和功率。5.2 电流环与速度环带宽没有拉开PMSM 双闭环控制有一个基本原理内环电流环带宽要远高于外环速度环带宽。一般电流环带宽设计在数百到上千赫兹速度环带宽只有它的十分之一甚至更低。如果两个 PI 控制器的响应速度接近就会出现“电流环还没跟上速度环又给出新的参考”整个环路产生耦合振荡。调参实际表现是转速曲线出现高频抖动或者 i_q_ref 一直在剧烈波动但转速变化很小。这时候第一反应不是继续调 PI 数值而是检查两个环路的采样时间电流环采样时间应该比速度环小一个数量级例如电流环 0.1 ms速度环 1 ms。把采样时间拉开再重新整定 PI症状通常会消失大半。5.3 数值振动与死区影响加入 PWM 开关和死区之后电流会出现明显的高频纹波这是正常现象。但如果纹波大到让电流环 PI 输出震荡就要检查电流采样是否跟开关周期不同步以及仿真步长是否足够小。在 Simulink 里若采用变步长求解器建议把最大步长设置为开关周期的 1/10 左右例如 10 kHz 开关频率对应最大步长 1e-5 秒。如果仍出现毛刺可以先在电流反馈路径加一阶低通滤波器截止频率放在开关频率的 1/3 左右。这种方法虽然会带来一点相位延迟但能让 PI 调控器的输出平滑很多。仿真阶段这样做没问题硬件实现时再去考虑采样同步和延迟补偿。5.4 参数脚本管理混乱刚开始我习惯在 Simulink 模型的模块里直接输入常量数字改一次参数翻遍所有模块效率很低而且容易漏改一个导致结果莫名其妙。后来我改成统一使用 MATLAB 脚本设定变量在模型参数框里填变量名。例如飞轮惯量填 J_fw极对数填 p磁链填 lambda_pm。这样每次调整都只需要在脚本里改一行并且可以用循环参数跑参数扫描。完整脚本不建议放在博文里贴满但结构可以参照% 飞轮与电机参数 J_fw 1.0; % 飞轮电机等效转动惯量 [kg·m^2] R_s 0.02; % 定子电阻 [Ω] L_d 0.3e-3; % d轴电感 [H] L_q 0.3e-3; % q轴电感 [H] lambda_pm 0.12; % 永磁磁链 [Wb] p 2; % 极对数 Vdc 750; % 直流母线电压 [V] B_damp 0.005; % 阻尼系数 [N·m·s/rad] T_friction 0.5; % 恒摩擦转矩 [N·m] I_max 300; % 电流限幅 [A]这类脚本配合 Simulink 的sim命令批量执行后就可以高效地对比不同 J 或不同 I_max 下的充放电特性曲线。很多参数的规律比如惯量越大充电时间越长、电流限幅越大功率越高都会在这种参数扫描里得到非常直观的验证。5.5 常见问题速查表现象可能原因处理办法充电时转速不上升i_q_ref 一直为 0检查速度环 PI 输出限幅与模式选择信号转速上升但能量不守恒阻尼损耗或摩擦损耗未计入在机械方程中加入 B×ω 和 T_f 项母线电压在放电时冲到过高没有为回馈能量设置负载接电阻负载或网侧回馈变流器模型电流波形高频振荡电流环带宽过高或采样时间过大降低电流环 PI、减小仿真步长仿真速度极慢使用了不适配的求解器换 ode23t 并设置最大步长转速在目标值附近反复穿越速度环 PI 与电流环采样混叠把速度环采样时间降到电流环 1/10 以下我在实际仿真里踩得最深的一个坑是一开始为了图方便把飞轮机械模型里的阻尼系数 B 设成了零结果在“保持工况”测试时转速纹丝不动能量曲线也没有任何下降。看起来像是很完美的运行实际却掩盖了飞轮储能“维持转速也需要消耗少量功率”这个重要事实。后来把 B 加回去之后才发现保持工况下逆变器其实要有一定的电流输出而这个电流在真实系统中直接影响系统效率和热损耗。飞轮储能系统的 Simulink 仿真做到最后最考验人的往往不是模型里某个模块怎么连而是你能不能把“能量从哪儿来、到哪儿去、中途损耗多少”这条线索在每一帧数据里追踪清楚。只要这条链路清晰了齿轮式地加进网侧变流器、SOC 估算、充放电策略切换都会自然地长在这套框架上。关于“PMSM 参数怎么辨识”“SVPWM 死区补偿怎么做”这类更细的模块后面我再单独展开分享。