前阵子整理模型库翻到我最早一版飞轮储能Simulink仿真里面还留着不少“凑合能跑”的痕迹。飞轮储能最近在电网调频、应急电源、轨道交通制动能量回收这类场景里被反复提起原理一句话就能说明白用一台可以双向运行的电机把电网或直流母线的电能变成飞轮高速旋转的机械能存起来需要时再让能量反向流动。其中电机选择上永磁同步电机(PMSM)凭借高效率和高转矩密度成了主流。但这套系统的仿真建模很多人卡住的地方不是电机本身而是怎么把“储能”这个整机概念落进模型结构。本文会从我实际搭建模型的顺序出发把PMSM驱动的飞轮储能建模要解决的事讲透先搞清楚飞轮储能数学模型里哪些量起决定性作用再谈PMSM在飞轮场景下和普通调速系统的差别接着是Simulink里的分层搭建过程最后是PI整定、步长、仿真效率以及我踩过最多的坑。适合正在做飞轮储能仿真课题的研究生以及想快速建立系统级模型的工程师。1. 飞轮储能系统的模型边界能量存储靠惯量和转速不是靠“把转盘画出来”1.1 能量存哪、可放多少先算清楚飞轮储能的物理本质很朴素一个旋转刚体在转速升高时把能量加速存入转速降低时再把能量放出来。存储在转动部分的动能是E (1/2) J ω²其中J是旋转体总转动惯量kg·m²ω是旋转角速度rad/s。这里有个很容易被忽略的工程含义储能容量和转速的平方成正比所以飞轮系统一般愿意把转速推得比较高而不是单纯把飞轮做得更重。实际工程项目很少从零转速一直充到最高转速寿命和轴承/电机约束会让飞轮在最低工作转速ω_min和最高工作转速ω_max之间运行。那么可用储能增量是ΔE (1/2) J (ω_max² - ω_min²)这个公式是做系统方案时最先要算的量。举个例子假如飞轮转动惯量为10 kg·m²工作转速范围取1500 rpm到3000 rpm换算成角速度就是157.08 rad/s到314.16 rad/s那么ΔE 5 × (314.16² - 157.08²) ≈ 370 kJ如果希望这个飞轮在10秒内完成一次充电或放电平均功率大约37 kW。别看能量数值不大放到短时高频调频或制动能量回收场景里这种“小容量、大功率”特征正好匹配需求。做Simulink模型之前先把这个数算出来后面很多参数设置就不会离谱。1.2 一个完整飞轮储能单元有哪些部分很多仿真初学者把“飞轮储能”理解成“飞轮加一台电机”但落到Simulink模型时需要更完整的系统眼光。工程上一个飞轮储能单元通常包含飞轮转子本体真正的能量存储元件通常做成大惯量盘体或细长圆柱体需要考虑材料强度和惯量分布永磁同步电机同一台电机既当电动机又当发电机电机转子与飞轮同轴连接双向变流器连接电机侧交流电和直流母线负责把电能形态和方向管理好轴承系统低损耗轴承或磁悬浮轴承减少静止时和低速时的摩擦损耗真空/冷却辅助系统高速旋转时降低风阻和散热这部分在详细热仿真里才需要纳入。如果目标是验证能量吞吐控制策略模型里可以只保留“飞轮惯量 PMSM 机侧变流器 直流母线”这一条主链。如果目标是研究并网波动平抑那还要加上电网侧变流器、滤波器和上层功率分配策略。模型边界不同仿真复杂度差一个数量级第一步先明确这个。1.3 第一版模型可以省去什么对于刚起步的仿真项目我强烈建议第一版不要追求“面面俱到”而是把注意力放在机电能量转换主回路上。可以暂时省略电机铁损、磁滞损耗和涡流损耗认为反电动势正弦且磁链恒定飞轮转子的柔性模态假设它是刚性轴系轴承的复杂摩擦力矩先等效成粘性阻尼系数B和库仑摩擦项磁悬浮系统的动态直接用低阻尼系数代替温度对永磁体磁链的影响。这样机电侧的核心方程就只剩下一组J_total (dω/dt) Te - TL - Bω其中J_total是电机转子惯量、飞轮惯量和联轴器惯量之和。Te是PMSM输出的电磁转矩TL是外部负载转矩Bω是等效粘性阻尼。飞轮储能时TL可以取零或很小如果仿真想更真实加一个与转速正相关的风阻损耗项也不会破坏模型结构。这组方程看着简单却是后面整个Simulink模型是否可信的基石。飞轮储能的“充”与“放”在机械侧就是电磁转矩Te符号的变化能量到底增加还是减少看的是Te和ω的乘积方向而不是单纯看转速值。2. PMSM在飞轮工况下和普通变频驱动有本质区别2.1 dq坐标系下的PMSM电压与转矩方程永磁同步电机的三相电压和磁链本身是交流量控制上通常把它变换到随转子旋转的dq坐标系。静止的A/B/C三相量先经过Clarke变换变成αβ静止坐标再经过Park变换变成dq坐标。在dq旋转坐标系下表贴式或内置式PMSM的定子电压方程为u_d R_s i_d L_d (di_d/dt) - ω_e L_q i_qu_q R_s i_q L_q (di_q/dt) ω_e (L_d i_d ψ_f)电磁转矩方程为T_e 1.5 p [ ψ_f i_q (L_d - L_q) i_d i_q ]其中R_s是定子电阻L_d和L_q是d/q轴电感ψ_f是永磁磁链p是极对数ω_e是电角速度。对于表贴式PMSML_d和L_q近似相等转矩就主要由q轴电流i_q决定这就是为什么常见FOC控制采用id0的策略。做飞轮储能仿真时牵涉这些恒等式不是为了考试而是为了看懂模块参数表。Simulink里的PMSM模块内部已经替你算好这些方程但你需要把电阻、电感、磁链、极对数填对填错一个值仿真结果就完全不可信。另外电角速度和机械转速之间是ω_e p ω_m极对数设置错误会让控制器输出频率和实际电机频率对不上。2.2 飞轮运行跨越四象限双向功率流才是常态普通调速系统的PMSM长期工作在电动状态偶尔需要制动减速工作点一般稳定在第一象限或者只在附近波动。飞轮储能系统却要求电机频繁在电动和发电两个方向切换而且切换往往发生在几分钟甚至几十秒的持续过程中。充电时电机吸收直流母线能量转子带动飞轮加速电磁转矩方向和旋转方向一致电功率从直流侧流向机械侧。放电时飞轮靠惯量维持旋转电机变成发电机电磁转矩方向和旋转方向相反机械能转为电能回馈到直流侧。在Simulink里体现最直接的就是转矩符号和电流方向会反向如果不理解这一点很容易出现“模型能转但放不出电”的问题。你可以先这样调试给一个手动转矩阶跃正转矩时观察转速和直流母线电流方向再给负转矩观察转速下降时电流是否反向。把这组逻辑理清了比盲目加复杂控制环更有效。2.3 为什么“额定转速”在飞轮系统里特别敏感普通变频电机通常在额定转速附近一个较小范围内运行使用范围可能只有0.8到1.2倍额定转速。飞轮储能系统不是这样它的释放能量能力依赖于转速范围转速范围直接决定可用储能容量。如果只在额定转速以下运行控制方式可以用简单的id0矢量控制逻辑清晰电流和转矩近线性对应。但实际飞轮要尽量提高能量密度最高转速往往会超出“基速”此时反电动势已经接近直流母线电压上限必须进入弱磁区域也就是给定负的i_d通过削弱气隙磁场来限制反电动势。弱磁区的控制参数整定、PI限幅和电压保护都要重新考虑。做第一版Simulink模型时可以先限定只仿真基速以下工况把主能量循环验证通再扩展弱磁段。这个顺序能少踩很多坑。2.4 飞轮储能为什么常选PMSM而不是异步电机飞轮储能电机需要低损耗、高效率和宽调速范围。异步电机在转子侧会有铜耗空载和轻载时损耗明显电励磁同步电机需要额外的励磁回路和滑环在高速真空系统中维护麻烦开关磁阻电机噪声和转矩脉动问题比较突出。PMSM因为转子是永磁体结构没有转子铜耗定子电流大部分都用于产生有效转矩能量转换效率高功率密度也高因此在飞轮储能这个“电能和动能反复转换”的场景里非常合适。不过PMSM也有代价高速端需要弱磁控制永磁体可能失磁而且电机造价较高。仿真中这些差异不会直接暴露但它们会影响你在模型里选择“磁链是否恒定”“是否包含损耗电阻”等边界条件。3. Simulink模型分层搭建电机、飞轮、变流器和FOC是怎么接起来的3.1 主电路层、控制层和参数层的结构设计拿过一个新的Simulink仿真任务我习惯于先画清楚三层结构而不是把所有模块堆在一个窗口里。主电路层包含直流电压源或直流母线电容、三相桥式逆变器、PMSM、飞轮转动惯量和机械参考控制层包含转速采样、电流采样、Clarke/Park变换、PI调节器、SVPWM或PWM调制参数层把电机参数、控制参数、飞轮惯量、阻尼、转速限制都集中到MATLAB工作区变量或模型回调中避免在模块里写死数字。分层的直接好处是调一种控制策略时不用动主电路验证不同参数时也不用翻几十个窗口。很多项目后来需要做参数扫描、批量仿真、硬件在环时这种分层结构的价值会成倍放大。3.2 机械侧到底怎么接扭矩信号接入和惯量模块接入是两码事这是新手最容易糊涂的地方。Simulink/Simscape里PMSM模块的机械端口有几种处理习惯取决于你用哪个库里的模块。如果用Simscape Electrical库里的电机模块机械侧通常要接到转动域。做法是在PMSM的机械旋转端口上依次连接转动惯量(Inertia)模块和机械参考(Mechanical Rotational Reference)再加上Solver Configuration。飞轮的惯量不是输入一个转速补偿值而是要表现为真实旋转体模块。你可以在机械转动域里把PMSM转子和飞轮看成一个共同旋转的刚体分别用两个Inertia模块串在同一根轴上也可以直接把总惯量合进一个模块。建议单独保留飞轮惯量模块这样建模时能清楚地区分“电机固有惯量”和“飞轮额外加上的惯量”。如果用老版本SimPowerSystems里的Permanent Magnet Synchronous Motor模块机械输入往往是一个Torque端子电机转速由内部运动方程算出来没有外部的机械轴域。这种情况下飞轮惯量就要通过修改PMSM模块内部的转子惯量参数来表示。两种方法在原理上等价但在模块连接方式上完全不兼容动手前先看清楚你用的是哪一套库。我个人的习惯是新版模型优先用Simscape Electrical带物理端口的方式结构更直观将来也方便扩展成多质量轴系。3.3 电机侧主回路接线和FOC控制链路这里给出一个适合做储能验证的模型接线思路。主电路从直流母线开始直流母线电压源/电容 → Universal Bridge
飞轮储能PMSM建模：Simulink系统仿真从原理到实现
前阵子整理模型库翻到我最早一版飞轮储能Simulink仿真里面还留着不少“凑合能跑”的痕迹。飞轮储能最近在电网调频、应急电源、轨道交通制动能量回收这类场景里被反复提起原理一句话就能说明白用一台可以双向运行的电机把电网或直流母线的电能变成飞轮高速旋转的机械能存起来需要时再让能量反向流动。其中电机选择上永磁同步电机(PMSM)凭借高效率和高转矩密度成了主流。但这套系统的仿真建模很多人卡住的地方不是电机本身而是怎么把“储能”这个整机概念落进模型结构。本文会从我实际搭建模型的顺序出发把PMSM驱动的飞轮储能建模要解决的事讲透先搞清楚飞轮储能数学模型里哪些量起决定性作用再谈PMSM在飞轮场景下和普通调速系统的差别接着是Simulink里的分层搭建过程最后是PI整定、步长、仿真效率以及我踩过最多的坑。适合正在做飞轮储能仿真课题的研究生以及想快速建立系统级模型的工程师。1. 飞轮储能系统的模型边界能量存储靠惯量和转速不是靠“把转盘画出来”1.1 能量存哪、可放多少先算清楚飞轮储能的物理本质很朴素一个旋转刚体在转速升高时把能量加速存入转速降低时再把能量放出来。存储在转动部分的动能是E (1/2) J ω²其中J是旋转体总转动惯量kg·m²ω是旋转角速度rad/s。这里有个很容易被忽略的工程含义储能容量和转速的平方成正比所以飞轮系统一般愿意把转速推得比较高而不是单纯把飞轮做得更重。实际工程项目很少从零转速一直充到最高转速寿命和轴承/电机约束会让飞轮在最低工作转速ω_min和最高工作转速ω_max之间运行。那么可用储能增量是ΔE (1/2) J (ω_max² - ω_min²)这个公式是做系统方案时最先要算的量。举个例子假如飞轮转动惯量为10 kg·m²工作转速范围取1500 rpm到3000 rpm换算成角速度就是157.08 rad/s到314.16 rad/s那么ΔE 5 × (314.16² - 157.08²) ≈ 370 kJ如果希望这个飞轮在10秒内完成一次充电或放电平均功率大约37 kW。别看能量数值不大放到短时高频调频或制动能量回收场景里这种“小容量、大功率”特征正好匹配需求。做Simulink模型之前先把这个数算出来后面很多参数设置就不会离谱。1.2 一个完整飞轮储能单元有哪些部分很多仿真初学者把“飞轮储能”理解成“飞轮加一台电机”但落到Simulink模型时需要更完整的系统眼光。工程上一个飞轮储能单元通常包含飞轮转子本体真正的能量存储元件通常做成大惯量盘体或细长圆柱体需要考虑材料强度和惯量分布永磁同步电机同一台电机既当电动机又当发电机电机转子与飞轮同轴连接双向变流器连接电机侧交流电和直流母线负责把电能形态和方向管理好轴承系统低损耗轴承或磁悬浮轴承减少静止时和低速时的摩擦损耗真空/冷却辅助系统高速旋转时降低风阻和散热这部分在详细热仿真里才需要纳入。如果目标是验证能量吞吐控制策略模型里可以只保留“飞轮惯量 PMSM 机侧变流器 直流母线”这一条主链。如果目标是研究并网波动平抑那还要加上电网侧变流器、滤波器和上层功率分配策略。模型边界不同仿真复杂度差一个数量级第一步先明确这个。1.3 第一版模型可以省去什么对于刚起步的仿真项目我强烈建议第一版不要追求“面面俱到”而是把注意力放在机电能量转换主回路上。可以暂时省略电机铁损、磁滞损耗和涡流损耗认为反电动势正弦且磁链恒定飞轮转子的柔性模态假设它是刚性轴系轴承的复杂摩擦力矩先等效成粘性阻尼系数B和库仑摩擦项磁悬浮系统的动态直接用低阻尼系数代替温度对永磁体磁链的影响。这样机电侧的核心方程就只剩下一组J_total (dω/dt) Te - TL - Bω其中J_total是电机转子惯量、飞轮惯量和联轴器惯量之和。Te是PMSM输出的电磁转矩TL是外部负载转矩Bω是等效粘性阻尼。飞轮储能时TL可以取零或很小如果仿真想更真实加一个与转速正相关的风阻损耗项也不会破坏模型结构。这组方程看着简单却是后面整个Simulink模型是否可信的基石。飞轮储能的“充”与“放”在机械侧就是电磁转矩Te符号的变化能量到底增加还是减少看的是Te和ω的乘积方向而不是单纯看转速值。2. PMSM在飞轮工况下和普通变频驱动有本质区别2.1 dq坐标系下的PMSM电压与转矩方程永磁同步电机的三相电压和磁链本身是交流量控制上通常把它变换到随转子旋转的dq坐标系。静止的A/B/C三相量先经过Clarke变换变成αβ静止坐标再经过Park变换变成dq坐标。在dq旋转坐标系下表贴式或内置式PMSM的定子电压方程为u_d R_s i_d L_d (di_d/dt) - ω_e L_q i_qu_q R_s i_q L_q (di_q/dt) ω_e (L_d i_d ψ_f)电磁转矩方程为T_e 1.5 p [ ψ_f i_q (L_d - L_q) i_d i_q ]其中R_s是定子电阻L_d和L_q是d/q轴电感ψ_f是永磁磁链p是极对数ω_e是电角速度。对于表贴式PMSML_d和L_q近似相等转矩就主要由q轴电流i_q决定这就是为什么常见FOC控制采用id0的策略。做飞轮储能仿真时牵涉这些恒等式不是为了考试而是为了看懂模块参数表。Simulink里的PMSM模块内部已经替你算好这些方程但你需要把电阻、电感、磁链、极对数填对填错一个值仿真结果就完全不可信。另外电角速度和机械转速之间是ω_e p ω_m极对数设置错误会让控制器输出频率和实际电机频率对不上。2.2 飞轮运行跨越四象限双向功率流才是常态普通调速系统的PMSM长期工作在电动状态偶尔需要制动减速工作点一般稳定在第一象限或者只在附近波动。飞轮储能系统却要求电机频繁在电动和发电两个方向切换而且切换往往发生在几分钟甚至几十秒的持续过程中。充电时电机吸收直流母线能量转子带动飞轮加速电磁转矩方向和旋转方向一致电功率从直流侧流向机械侧。放电时飞轮靠惯量维持旋转电机变成发电机电磁转矩方向和旋转方向相反机械能转为电能回馈到直流侧。在Simulink里体现最直接的就是转矩符号和电流方向会反向如果不理解这一点很容易出现“模型能转但放不出电”的问题。你可以先这样调试给一个手动转矩阶跃正转矩时观察转速和直流母线电流方向再给负转矩观察转速下降时电流是否反向。把这组逻辑理清了比盲目加复杂控制环更有效。2.3 为什么“额定转速”在飞轮系统里特别敏感普通变频电机通常在额定转速附近一个较小范围内运行使用范围可能只有0.8到1.2倍额定转速。飞轮储能系统不是这样它的释放能量能力依赖于转速范围转速范围直接决定可用储能容量。如果只在额定转速以下运行控制方式可以用简单的id0矢量控制逻辑清晰电流和转矩近线性对应。但实际飞轮要尽量提高能量密度最高转速往往会超出“基速”此时反电动势已经接近直流母线电压上限必须进入弱磁区域也就是给定负的i_d通过削弱气隙磁场来限制反电动势。弱磁区的控制参数整定、PI限幅和电压保护都要重新考虑。做第一版Simulink模型时可以先限定只仿真基速以下工况把主能量循环验证通再扩展弱磁段。这个顺序能少踩很多坑。2.4 飞轮储能为什么常选PMSM而不是异步电机飞轮储能电机需要低损耗、高效率和宽调速范围。异步电机在转子侧会有铜耗空载和轻载时损耗明显电励磁同步电机需要额外的励磁回路和滑环在高速真空系统中维护麻烦开关磁阻电机噪声和转矩脉动问题比较突出。PMSM因为转子是永磁体结构没有转子铜耗定子电流大部分都用于产生有效转矩能量转换效率高功率密度也高因此在飞轮储能这个“电能和动能反复转换”的场景里非常合适。不过PMSM也有代价高速端需要弱磁控制永磁体可能失磁而且电机造价较高。仿真中这些差异不会直接暴露但它们会影响你在模型里选择“磁链是否恒定”“是否包含损耗电阻”等边界条件。3. Simulink模型分层搭建电机、飞轮、变流器和FOC是怎么接起来的3.1 主电路层、控制层和参数层的结构设计拿过一个新的Simulink仿真任务我习惯于先画清楚三层结构而不是把所有模块堆在一个窗口里。主电路层包含直流电压源或直流母线电容、三相桥式逆变器、PMSM、飞轮转动惯量和机械参考控制层包含转速采样、电流采样、Clarke/Park变换、PI调节器、SVPWM或PWM调制参数层把电机参数、控制参数、飞轮惯量、阻尼、转速限制都集中到MATLAB工作区变量或模型回调中避免在模块里写死数字。分层的直接好处是调一种控制策略时不用动主电路验证不同参数时也不用翻几十个窗口。很多项目后来需要做参数扫描、批量仿真、硬件在环时这种分层结构的价值会成倍放大。3.2 机械侧到底怎么接扭矩信号接入和惯量模块接入是两码事这是新手最容易糊涂的地方。Simulink/Simscape里PMSM模块的机械端口有几种处理习惯取决于你用哪个库里的模块。如果用Simscape Electrical库里的电机模块机械侧通常要接到转动域。做法是在PMSM的机械旋转端口上依次连接转动惯量(Inertia)模块和机械参考(Mechanical Rotational Reference)再加上Solver Configuration。飞轮的惯量不是输入一个转速补偿值而是要表现为真实旋转体模块。你可以在机械转动域里把PMSM转子和飞轮看成一个共同旋转的刚体分别用两个Inertia模块串在同一根轴上也可以直接把总惯量合进一个模块。建议单独保留飞轮惯量模块这样建模时能清楚地区分“电机固有惯量”和“飞轮额外加上的惯量”。如果用老版本SimPowerSystems里的Permanent Magnet Synchronous Motor模块机械输入往往是一个Torque端子电机转速由内部运动方程算出来没有外部的机械轴域。这种情况下飞轮惯量就要通过修改PMSM模块内部的转子惯量参数来表示。两种方法在原理上等价但在模块连接方式上完全不兼容动手前先看清楚你用的是哪一套库。我个人的习惯是新版模型优先用Simscape Electrical带物理端口的方式结构更直观将来也方便扩展成多质量轴系。3.3 电机侧主回路接线和FOC控制链路这里给出一个适合做储能验证的模型接线思路。主电路从直流母线开始直流母线电压源/电容 → Universal Bridge
