最近在做一个飞轮储能系统的仿真项目电机端用的是永磁同步电机PMSM任务是把从物理储能、电机控制到Simulink仿真的整条链路打通。飞轮储能这个概念并不新但真正能把模型搭到能用的状态、并且波形和能量流都对得上还是有不少门道。这篇文章我把整个建模思路、参数计算、控制策略和踩过的坑完整捋一遍适合正在做飞轮储能、PMSM变频驱动、或者想用Simulink做机电储能仿真的同行参考。先说结论飞轮储能仿真的核心难点不在电机本体而在三个地方——转动惯量与能量区间怎么匹配、充放电双工况下的控制结构如何平滑切换、以及仿真步长与模型精度之间的平衡。把这三点想清楚模型基本就成了一半。1. 飞轮储能系统拓扑与PMSM的定位逻辑1.1 从电网到飞轮的功率链路飞轮储能系统的整体结构并不复杂电网侧经过双向AC/DC变流器接到直流母线母线上再挂电机侧逆变器逆变器驱动永磁同步电机电机转轴与飞轮刚性连接飞轮封装在真空腔体里。常规的容量等级大概在几十千瓦到几兆瓦之间转速从一万转到几万转不等。在这条链路上电机不是普通意义上的拖动电机它承担了电能与机械能的双向转换。充电时电机工作在电动状态把直流母线的电能变成转子动能放电时电机变成发电机飞轮的惯性带动转子旋转把动能变回电能通过逆变器回馈到直流母线。换句话说电机就是飞轮储能系统的能量转换器而且是双向的。为什么选永磁同步电机而不是异步电机核心原因有两个。第一PMSM转子无绕组没有转子铜耗高速工况下效率明显占优。飞轮储能最怕的就是损耗一转起来就是几万转异步机的转子铜耗会直接拉低往返效率。第二PMSM功率密度高、响应快磁场定向控制下转矩可以精确到毫牛米级别这对于充放电功率的精确控制非常关键。而且PMSM可以和飞轮共用一根轴做成一体化结构省去齿轮箱之类的中间环节。1.2 充电、悬浮保持、放电三工况的能量变化关系飞轮储能的一个显著特点是储能多少和转速平方成正比。动能公式是E 1/2 × J × ω²其中J是转动惯量kg·m²ω是飞轮角速度rad/s。因为转速是平方项所以存在一个非常实用的工作区间设计原则实际工程中不会让飞轮把能量全部放光。通常最低转速取最高转速的一半左右这样可释放的能量比例是多少用公式算一下E_release / E_max ω_max² - ω_min²/ ω_max²如果ω_min 0.5 × ω_max那么可释放能量占75%。剩下25%的动能属于保底用来维持飞轮旋转和低速时的反电动势支撑。这三个工况对应的控制本质完全不同。充电是电动加速电磁转矩为正转速上升悬浮保持是让飞轮维持额定转速转矩只需克服摩擦和风阻损耗数值很小放电是发电制动电磁转矩为负转速下降。在Simulink里做仿真时这三个工况不是三个独立模型而是同一套控制系统下转矩指令的符号切换和幅值调节这一点后面第四章细说。2. 飞轮本体建模与储能参数计算先把惯量和能量算准2.1 为什么储能计算要盯住转动惯量而不是质量很多人一上来就问飞轮多重其实飞轮设计里真正决定储能的是转动惯量J而不是单纯的质量。同样质量的飞轮材料分布越靠外转动惯量越大储能越多。这就是为什么飞轮通常做成轮缘厚、中心薄的形状。对于实心圆柱飞轮转动惯量是J 1/2 × m × R²其中m是质量R是圆柱半径。如果是轮缘式飞轮内径R1、外径R2公式是J 1/2 × m ×R2² R1²从结构上看半径的贡献是平方级的所以同等重量下把材料放在外缘远比增加厚度有效。这也解释了为什么高速飞轮大多采用细长轴加厚轮缘的构型。2.2 一个50kW级飞轮的具体参数计算全程我做仿真时按一个50kW、3kWh左右的飞轮系统来设定参数。飞轮材料用钢材密度取7850kg/m³设计成实心圆柱半径 R 0.25 m轴向长度 L 0.3 m质量 m ρ × π × R² × L 7850 × 3.1416 × 0.0625 × 0.3 ≈ 462 kg转动惯量 J 1/2 × m × R² 1/2 × 462 × 0.0625 ≈ 14.4 kg·m²最高转速取12000 rpm换算成角速度ω_max 12000 × 2π / 60 ≈ 1256.6 rad/s最高转速下的总储能为E_max 1/2 × J × ω_max² 1/2 × 14.4 × 1256.6² ≈ 11.4 MJ ≈ 3.17 kWh这个储能配50kW功率按75%放电深度算能支撑大约3分钟左右的满功率放电符合储能调频的典型应用场景。仿真时可以直接把这个惯量值填进PMSM模块的机械参数里。2.3 风阻摩擦损耗怎么在仿真中体现飞轮在空气中的风阻损耗和转速近似成立方关系所以真机都会把飞轮放在真空腔体中。但在实际仿真中如果完全忽略机械损耗结果会偏理想化——尤其是悬浮保持阶段的电流波形会和实际差很多。建议在电机轴系上增加两部分损耗一部分是常数摩擦转矩T_f模拟轴承摩擦另一部分是与转速平方或三次方相关的风阻转矩T_w k_w × ω²。在Simulink里可以把它建模成一个受转速控制的力矩源作用方向与旋转方向相反。真空环境下k_w取一个很小的值即可但如果仿真目标是做真空腔体设计或者长时间保持工况分析这一步就不能省。3. 永磁同步电机的dq模型与Simulink参数配置3.1 PMSM的电压方程与表贴式电机的关键简化PMSM在三相静止坐标系下的方程是非线性的交流量不方便做PI控制。所以仿真的第一步是把它变换到dq旋转坐标系。在dq坐标系下电压方程为ud Rs × id Ld × d(id)/dt - ωe × Lq × iquq Rs × iq Lq × d(iq)/dt ωe ×Ld × id ψf电磁转矩方程为Te 1.5 × p × [ψf × iq Ld - Lq× id × iq]这里面p是极对数ψf是永磁体磁链Rs是定子电阻Ld和Lq分别是d轴和q轴电感。我做仿真时选的是表贴式PMSM它的特点是Ld和Lq基本相等。这样转矩方程可以大幅简化为Te 1.5 × p × ψf × iq这个简化的物理意义非常直接电磁转矩只和q轴电流成正比和d轴电流无关。这就是id0控制的根本原因——把电流全部用在产生转矩的q轴上d轴电流保持为零既不产生有效转矩还增加铜耗。飞轮储能这种低速大转矩、高速低转矩的应用场景用id0控制既简单又高效完全够用。3.2 Simulink中PMSM模块的核心参数设置Simulink的电气库里有现成的永磁同步电机模型参数配置界面大致这几项参数数值我用的说明定子电阻 Rs0.05 Ω影响铜耗和PI整定d轴电感 Ld0.35 mH表贴式q轴电感 Lq0.35 mH与Ld相同即为表贴式永磁体磁链 ψf0.12 Wb决定反电动势和转矩常数极对数 p2高速机常用2对极转动惯量 J14.45 kg·m²必须包含飞轮惯量初始转速628.3 rad/s对应6000rpm直接从中速起步这里最容易出错的是转动惯量。电机模块默认的J通常只有0.0几kg·m²如果忘了改飞轮一个阶跃指令就能从零飙到几万转仿真结果完全失真。记住J用的是电机转子加飞轮的总惯量。我习惯把电机自身转子惯量忽略相对飞轮太小直接把14.4 kg·m²填进去。3.3 坐标变换与SVPWM在仿真链路中的真实位置从控制结构上看矢量控制链路是采样三相电流 → Clark变换abc→αβ → Park变换αβ→dq → dq轴PI调节器 → 反Park变换 → SVPWM → 逆变器 → 电机。在Simulink里搭建时这一串模块每个都能在库里找到现成的。重要的是理解每个环节的作用Clark变换和Park变换是为了把交流量变成可做PI控制的直流量PI调节器输出的是dq电压指令反Park变换把它变成αβ坐标系下的电压矢量SVPWM再把电压矢量转成三相桥臂的开关信号。如果你用的是仿真库自带的高精度PMSM模型其实可以直接给定子侧加三相电压源代替逆变器那SVPWM环节就可以省略。但这样无法反映母线电压波动和开关纹波的影响无法观察充放电切换时母线电压的跌落特性。所以我建议第一版保留逆变器SVPWM等整条控制链路调试稳定了再做平均值模型的简化版本。4. 充放电控制策略与PI参数整定可以抄作业的部分4.1 外环是电压环还是功率环飞轮储能系统的控制需求取决于应用场景。如果是做平滑新能源功率波动外环是功率指令如果是做调频或母线电压支撑外环是母线电压。我这次做的是直流母线电压支撑场景外环采集直流母线电压与额定值750V比较经过PI调节器输出有功功率指令再折算成q轴电流指令。充电时q轴电流为正电机加速放电时q轴电流为负电机发电。从控制结构上看充放电切换根本不涉及切换控制律只是指令符号的自然翻转。这里要注意一个工程细节转矩指令不能从正到负直接阶跃否则机械轴系冲击极大仿真波形也会出现剧烈的电压过冲。我习惯在转矩给定通道上加一阶惯性限制斜坡时间取0.1~0.5s模拟实际变流器的功率变化率限制。4.2 电流内环PI参数的一个直接给法电流环的被控对象近似是一个RL惯性环节时间常数τ L/R。对于表贴式PMSMd轴和q轴对称所以两个电流环的PI参数相同。工程整定方法很简单设期望电流环带宽为αrad/s则Kp α × LKi α × R这里有个经验法则α不要超过开关频率的1/10。比如开关频率取10kHzα取1000rad/s。用我前面的参数Kp 1000 × 0.00035 0.35Ki 1000 × 0.05 50这个参数可以直接作为初值。注意α选得太大会导致离散化后相位裕量不足电流环高频振荡选得太小则电流跟踪慢储能响应速度跟不上。电流环之外还要设两个限幅输出限幅对应最大相电流和积分限幅。积分限幅很关键——启动时如果不限制积分项q轴电流会严重超调转速爬升过程中出现明显的转矩尖峰。我的经验是积分限幅取额定转矩对应iq值的1.2倍。4.3 转速保护与最低转速约束飞轮储能和普通电机驱动还有一个不同点转速不仅是被控量也是储能水平的直接反映。仿真中必须加上转速保护逻辑当转速达到1.05倍额定转速时封锁PWM或强制进入悬浮保持当转速低于0.5倍额定转速时禁止继续放电防止反电动势过低导致母线电压塌陷。这层逻辑在Simulink里用几个比较器和一个RS触发器就能搭出来。别小看它没有这层保护你做长时间充放电循环仿真时转速一冲到上限电流环就饱和了波形就会乱掉。5. Simulink模型的顶层搭建骨架与求解器选择5.1 模型顶层结构和推荐模块把模型分成四个部分来搭主电路、电机、控制、观测结构清晰也方便挑错。主电路用直流电压源直流电容通用三相桥IGBT模块组成。母线电容参数对电压波动影响很大我这个50kW系统取了C20mF后面会专门分析电容对母线电压波动的定量关系。电机直接用电气库的PMSM模型机械端口接到飞轮惯量模块。控制部分用离散PID控制器采样时间取1e-4s与10kHz开关频率对应。观测部分做一个SOC荷电状态计算模块。飞轮SOC不建议用转速线性表示而应该用动能比例SOC ω² - ω_min²/ω_max² - ω_min²× 100%这个公式直接对应还能放出多少能量的物理含义是飞轮储能区别于电池储能的特殊指标。5.2 求解器选择连续还是离散步长多大这是新手最容易卡死的环节。PMSM模型加逆变器是典型的非线性快速系统如果用连续变步长求解器默认最大步长往往过大导致仿真结果剧烈震荡如果强行调小最大步长仿真速度又慢到无法接受。我的做法是第一版直接用离散固定步长求解器步长取1e-5s开关周期10kHz的十分之一。这样能保证开关纹波被正确采样同时避免连续求解器反复试探步长带来的额外开销。等到第一版调通再根据需求做平均值模型提速。第二个实用技巧是仿真时长的控制。一次完整的充放电循环按前面算的从6000rpm加速到12000rpm需要约250秒如果直接用开关级模型跑完仿真步数达到2500万步耗时完全不可接受。解决思路是分两步走先用小步长比如0.5~2s的物理时间验证控制环动态再用平均值模型或者直接给转速一个斜坡设定值快速跑完整个储能循环观察SOC和能量流。5.3 效率计算模块里容易踩的坑效率计算模块看起来简单实际上很多人算错铜损。总铜损的正确公式是Pcu 1.5 × Rs ×id² iq²为什么是1.5因为dq坐标系下的电流幅值是相电流幅值相电流有效值是幅值的√2分之一三相总损耗是3 × Rs × I_rms²I_rms² id² iq²)/2所以系数是3/2。有人直接套3×Rs×(id²iq²)算出来的损耗直接翻倍后面效率分析全歪了。6. 仿真波形解读与关键参数影响量化6.1 一次完整充放电循环的波形观察点第一版验证用短时仿真设定初始转速628 rad/s充电50s悬浮保持20s放电30s。重点观察四个信号转速、q轴电流、直流母线电压、SOC。充电段的预期行为是转速从628 rad/s线性升到约800 rad/s50s内q轴电流在限幅范围内保持恒定母线电压稳定在750V附近。如果q轴电流出现锯齿状波动多半是PI参数或坐标变换角度有问题。悬浮保持段q轴电流降到一个很小值只补偿损耗SOC几乎不变。放电段q轴电流反向转速下降母线电压会出现小幅跌落——跌落的幅度由母线电容和放电功率共同决定。6.2 母线电容对电压波动的定量影响母线电压波动这个问题需要量化地看。电容上的电荷变化满足C × ΔV I × Δt直流侧电流I约等于P/Vdc即50kW/750V ≈ 67A。如果要求在10ms内完成从空载到满载的功率阶跃且允许电压波动不超过33V约4.4%C I × Δt / ΔV 67 × 0.01 / 33 ≈ 0.02 F这就是我取20mF电容的来源。如果电容选小了比如只有5mF同样条件下电压波动会达到130V以上母线电压直接突破工作范围。在仿真里做参数扫描时把C从5mF扫到50mF看母线电压的峰值跌落实在直观。6.3 转动惯量变化对动态特性的影响飞轮系统的功率响应速度和惯量直接相关。同样的充电电流下J从14.4 kg·m²减半到7.2 kg·m²转速上升速度翻倍储能容量减半放电时的母线电压跌落则更深更快。这个特性的工程意义是如果想提升飞轮的快速响应能力减小惯量是手段之一但代价是储能容量下降。仿真中做两组参数扫描对比可以很清楚地看到这个权衡关系。这也是为什么实际飞轮设计必须在强度、惯量和运行转速三个变量之间反复折中。7. 仿真踩坑实录从发散到效率对不上的解决链路7.1 启动抖动转子初始位置与坐标变换角度不一致有一次模型一启动转速没有按斜坡上升反而先在零附近来回震荡了两三秒q轴电流出现了很大的正负交替尖峰。排查后发现PMSM模块的初始转子位置设成了0但坐标变换模块的Park变换初始角度也是0两者理论上应该一致问题出在初始速度——电机模块默认初始转速为0我从零开始加速时低速阶段反电动势很小电流环在强行推一个几乎没有反电动势的电机导致电流尖峰。解决办法仿真初始转速直接设为628rad/s对应SOC75%而不是从零开始。这样一开始就处于正常的工作区间启动瞬间的控制量是平滑的。这也符合飞轮储能的实际运行场景——飞轮正常情况下不会停机始终维持在最低工作转速以上。7.2 照抄的PI参数发散了我前面给出的Kp0.35、Ki50是基于10kHz开关频率。如果把开关频率改成5kHz同样的带宽α1000rad/s就会有问题——因为离散化后的控制延迟增大系统相位裕量不足电流环出现高频振荡。我也踩过这个坑。当时把开关频率从10kHz改到8kHz带宽没变仿真一跑id和iq直接发散。解决办法不是调PI增益而是把带宽α降到开关频率的1/20也就是400rad/s左右。记住这个规律α要与开关频率成比例降低而不是独立整定。7.3 母线电压在充放电切换瞬间过冲第一次做完整的充转放切换时母线电压在切换瞬间冲到了800V以上超过了额定值的6%。原因不是PI问题而是转矩指令阶跃变化时电流环响应延迟造成短暂的能量不平衡——飞轮还在高速转动往母线回馈能量而负载还没来得及吸收。解决方法是给转矩指令加斜坡限制。我把斜坡时间设为0.2s后母线电压过冲降到了40V以内。这个场景也直接对应实际变流器的功率斜坡控制不是仿真的bug而是真实的工程约束。7.4 储能效率算出来超过100%调试效率模块时出现过一次非常离谱的结果放出的能量比充进去的还多。查了一圈发现是两个模型叠加了机械损耗但放电时损耗的符号写反了——制动工况下机械损耗的方向与旋转方向相同不对机械损耗力矩永远与旋转方向相反在放电工况下它消耗的是机械能所以损耗项在能量平衡方程中应该作为从机械能里扣掉的项目。符号写反之后损耗变成了额外注入的能量效率自然就爆表了。这个问题的排查过程提醒我效率模块最好用积分能量来验证而不是瞬时功率。将充电段P_in对时间积分、放电段P_out对时间积分、再扣除全部损耗项能量守恒才对得上。这比盯着瞬时波形更可靠。工作区间的保护逻辑在长时间仿真里尤其重要。有一次我把最低转速保护关掉让系统持续放电到接近零转速结果转速过零后模型出现了明显的数值不稳定母线电压也开始剧烈振荡。因为PMSM的反电动势与转速成正比接近零速时电流环已经没有能力控制母线电压了。实际工程中不会让飞轮放到这个程度仿真里更不该这样做。总结一下我这轮项目下来的整体感受飞轮储能仿真真正花时间的不是搭模型本身而是参数之间的匹配——惯量、电容、带宽、限幅值每个参数都和你想要的工作区间绑定。建议新手先照着我这个流程把一套完整参数跑通再根据自己项目的功率和储能需求去调整不要一上来就拿具体真机数据去怼仿真模型。模型调顺以后再逐步加入风阻、铁损、逆变器损耗这些非理想项也不迟。