飞轮储能这个方向这几年在电网调频、UPS、轨道交通能量回收这些场景里越来越热。做仿真模型的时候最常见也最核心的一套方案就是用永磁同步电机当能量转换核心通过控制电机来实现飞轮的充放电。我见过不少人卡在这个控制策略的搭建上尤其是充电外环转速、内环电流这套双闭环逻辑看起来简单真要把参数调稳、把充放电切换做顺还是有不少细节坑。这篇把我自己的建模和调试经验完整拆开讲讲。1. 项目核心思路与整体方案拆解1.1 为什么飞轮储能系统选择永磁同步电机飞轮储能系统的本质是用一个高速旋转的飞轮转子把电能转化成机械能存起来需要的时候再把机械能转化回电能。这套系统里电机既是电动机也是发电机充电的时候电动运行把飞轮加速放电的时候发电运行把飞轮的动能变成电能送出去。整个能量转换链路里电机和它的控制器决定了系统的效率上限、动态响应速度和稳定边界。选永磁同步电机而不是异步电机核心原因有三点第一是效率永磁电机的转子没有铜耗和铝耗励磁损耗几乎为零高效率区间宽这对储能系统来说很关键飞轮储能本来拼的就是循环效率第二是功率密度同样的体积和重量永磁电机能输出的功率和转矩更大飞轮系统为了减小风阻损耗通常会把电机和飞轮封装在真空腔体里体积约束很紧第三是控制精度永磁同步电机的转子磁场位置可以通过编码器精确获得矢量控制下电流环和转速环都能做到很高的控制带宽这对频繁充放电切换的动态性能是刚需。1.2 飞轮储能充放电控制的双闭环结构选型这套模型采用的控制架构是很标准的转速-电流双闭环结构。充电过程外环控制转速内环控制电流这是大部分飞轮储能充电控制的基本骨架。外环转速环根据目标转速和实际转速的偏差输出一个转矩电流指令q轴电流参考值内环电流环跟踪这个指令通过电压源逆变器输出电压矢量驱动电机加速。放电的时候控制目标从转速维持变成功率输出或母线电压维持控制逻辑需要切换。整体系统框图如下转速指令 → [转速环PI] → iq_ref → [电流环PI] → dq电压指令 → SVPWM → 逆变器 → PMSM ↑ ↓ 转速反馈 ← 编码器 电流采样 → iq/id反馈这套结构的好处是动静分离转速环负责宏观的能量管理电流环负责微观的转矩快速响应。电流环带宽通常整定到几百赫兹转速环带宽在几十赫兹的量级两者差一个数量级以上才能保证不互相耦合振荡。放电模式的切换逻辑也在这套模型里做了处理当系统收到放电指令后外环从转速环切换为直流母线电压环控制目标从飞轮转速变为母线电压恒定内环仍然控制电流。这个时候电机会自然进入发电工况转速下降动能转化为电能输出。1.3 这套控制方案解决了什么问题飞轮储能最怕的就是充放电切换瞬间的大电流冲击和转矩突变直接表现就是母线电压跌落、电流尖峰、甚至机械轴系扭振。双闭环控制通过两个层面的协调把这些问题控制在合理范围转速环让转速按设定好的斜坡或S型曲线平滑升降避免了直接给全转矩导致飞轮转子承受过大的剪切应力。飞轮的储能量和转速的平方成正比转速控制得稳能量管理就稳。电流环则把电磁转矩限制在设定范围内防止大电流冲击损坏功率器件和永磁体退磁。另外这套模型里面把电机模型、逆变器模型、控制算法全部放进了Simulink环境可以方便地验证不同工况下的控制性能。我在搭这个模型之前也参考过不少现成的飞轮储能仿真例子但很多模型要么控制策略太简单没有完整的充放电切换逻辑要么电机模型用理想模型替换控制效果失真。这套模型好就好在完整保留了永磁同步电机的电磁动态和逆变器的开关特性控制算法也按实际工程的做法来设计仿真结果有说服力。2. 核心控制参数设计与关键细节2.1 永磁同步电机在Simulink中的建模方式选择Simulink里做永磁同步电机的建模有三条路用电气专业的永磁同步电机模块、用Simscape Electrical里的物理模型、或者自己搭状态方程。我在这套模型里用的是电气专业Simscape Electrical老版本叫SimPowerSystems或者Specialized Power Systems的永磁同步电机模块配合三相两电平逆变器和直流母线构成主电路。这个选择的原因很简单电机模块内部已经考虑了凸极效应、永磁体磁场谐波、阻尼绕组等实际因素比自己搭状态方程准确得多。而且逆变器的IGBT开关动作可以通过PWM信号驱动的MOSFET/IGBT桥臂来建模比用理想的受控电压源替代要真实得多能观察到电枢电流的谐波成分和开关纹波。用基于物理的建模有个问题是仿真速度慢尤其是开关频率高的时候。我试验过Continuous和Discrete两种模式连续模式下仿真精度高但速度较慢离散模式下需要设定合理的采样时间推荐用1e-6秒的步长和离散求解器速度提升明显而且便于后续部署到实时仿真平台。2.2 转速环和电流环PI参数的计算过程双闭环控制的核心就是两个PI调节器。我直接给出我这里用的参数整定过程整个过程全部可以套公式算不需要靠感觉调。电流环采用PI控制其控制对象是永磁同步电机的定子绕组简化为一阶惯性环节。电流环的目标是快、准、稳因此采用典型I型系统整定。以d轴电流环为例给定电流环PI参数为Kp_i和Ki_i电流环控制对象传递函数近似为G_i(s) (Kp_i * s Ki_i) / s * (1 / (Ls R))其中L是d轴电感表贴式PMSM中LdLqR是定子电阻。按典型I型系统整定取Kp_i / Ki_i L / R即零极点对消后电流环开环增益为Kp_i/L要求闭环带宽为ω_ci。最终参数计算公式Kp_i ω_ci * LKi_i ω_ci * R我这里电机的参数是L 0.000835H8.35mHR 0.0325Ω等效到标幺值后按p.u.系统整定更方便。实际在模型里我用的是标幺值系统基值电流取峰值电流基值电压取额定相电压峰值。电流环带宽ω_ci取500Hz即3140rad/sKp_i 3140 * 0.000835 2.62标幺约等于0.52 Ki_i 3140 * 0.0325 102转速环的整定思路类似但需要把电机运动方程考虑进去控制对象变成电流环闭环加上机械转动惯量。转速环带宽ω_cn一般取电流环带宽的1/5到1/10我这里取50Hz即314rad/s。转速环PI参数计算公式Kp_n J * ω_cn / (1.5 * p * Ψf) Ki_n Kp_n * ω_cn / 5其中J是飞轮加电机的总转动惯量我这里飞轮转子和电机的合成转动惯量约为0.052kg·m²p是极对数这里取4Ψf是永磁体磁链约0.192Wb。代入Kp_n 0.052 * 314 / (1.5 * 4 * 0.192) 7.09 Ki_n 7.09 * 314 / 5 445.7实际上在Simulink里做标幺化之后PI参数要除以对应的基值折算系数但比例关系不变。初版模型里就按这个计算值先填进去出来的波形基本稳然后微调。提示电流环带宽不能取太高受限于逆变器开关频率一般不超过开关频率的1/10。如果开关频率是10kHz电流环带宽取500Hz-1kHz是合理的。超过这个范围数字控制延迟会带来相位裕度不足系统就会振荡。2.3 SVPWM调制与坐标变换的工程细节模型里电流环输出的dq轴电压指令需要经过反Park变换和SVPWM调制才能生成逆变器的开关信号。SVPWM的实现逻辑我直接用了Simulink里的SVPWM模块但在参数设置上有几个点值得注意直流母线电压的利用率SVPWM的线性调制范围比SPWM高15%最大输出相电压基波幅值可以达到Udc/√3。这意味着同样的直流母线电压SVPWM能让电机跑得更高转速这对飞轮储能来说很有意义因为飞轮转速通常很高。矢量切换的零矢量分配SVPWM采用七段式还是五段式调制对开关损耗和电流谐波影响很大。七段式每个开关周期内零矢量对称分布电流谐波抑制效果更好但开关损耗大五段式通过减少开关次数降低损耗但谐波略大。我的模型里用的是七段式追求控制精度优先。过调制当电机运行在弱磁区域时电压指令超过线性调制范围如果不过处理直接限幅输出波形会畸变。模型里做了过调制处理把电压矢量限幅到SVPWM的最大线性范围即Udc/√3对应的圆内同时保留一个弱磁控制器来调整d轴电流确保高速区也能稳定工作。坐标变换环节看起来简单但转子位置角对齐关系到整个控制能否建立。我习惯把电机的零位对齐到A相绕组轴线然后通过编码器信号模型里用增量编码器模块模拟获得转子位置再换算成电角度。有一点容易被忽略如果电机的极对数是4机械角度和电角度的换算系数是4在反馈环节漏乘这个系数电流环直接发散这个坑我踩过。2.4 飞轮转动惯量与转速环动态特性的匹配飞轮储能系统的机械时间常数通常很大因为飞轮的转动惯量不小转速响应相对电流会很慢。这个特性既是好事也是需要处理的问题。好处在于转速环对被控对象参数变化不敏感转动惯量的测量误差对转速环稳定性的影响不大PI参数有较大的裕量。坏处在于如果按转速环带宽50Hz去整定会出现转速响应过快、电流限幅导致加速转矩不够的情况。所以我在模型里对转速环做了两个处理。第一转速指令发生器采用斜坡/梯形波限制转速变化率让转速环的给定变化平缓。这个限速率的做法在工程上等同于给飞轮系统一个最大加速度限制保护机械轴系。第二转速环PI输出端串联了一个限幅器把iq指令幅值限制在逆变器和电机的最大允许电流范围内。限幅值按电机峰值转矩电流设置一般取额定电流的1.2到1.5倍。这两个处理合在一起相当于在原有双闭环外面又加了一层动态特性和机械约束协调层。很多初学者只盯着PI参数调忽略了限幅和给定整形结果仿真出来转速环不稳定电流波形全是限幅饱和的方波这个现象我后面在常见问题里详细讲。3. 仿真模型搭建与充放电实现3.1 主电路与电机参数的设置模型的主电路结构是直流电压源模拟母线或是前级变换器输出→ 三相两电平逆变器 → 永磁同步电机 → 飞轮负载转动惯量摩擦阻尼。直流母线电压我这里设的是700V这是一个比较典型的中小型飞轮储能系统母线电压等级。逆变器的IGBT导通压降设为1.8V开关频率10kHz死区时间2μs。这些参数都按照实际工程器件水平来设不能设成理想开关否则仿真结果偏乐观后续做硬件在环会有落差。永磁同步电机参数按一台额定功率约55kW、额定转速12000rpm的高速电机来设置具体参数如下参数数值说明定子电阻R0.0325Ω常温下的绕组电阻d轴电感Ld0.000835H表贴式结构LdLqq轴电感Lq0.000835H与Ld相等永磁体磁链0.192Wb以转子位置0度时A相反电动势峰值为准极对数4高速电机常用4-8极转动惯量0.052kg·m²电机转子飞轮折算后的合成惯量这里我需要提醒一下不同文献里永磁同步电机的参数基准可能不一样有的用线电压有效值有的用相电压峰值。从电机厂家拿到的参数表需要先统一化到Simulink模块使用的坐标系中。最简单的验证办法空载反电动势仿真结果和理论值对比设一个500rpm的转速测A相反电动势波形幅值和0.19242π*500/60对比一下误差在3%以内说明参数输入正确。3.2 充电过程的转速外环电流内环控制实现充电过程的控制目标是把飞轮从静止加速到额定转速并在达到目标转速后维持转速不变此时系统处于待机状态飞轮以额定转速空转仅克服风阻摩擦损耗维持转速。转速外环的输入是转速指令与实际转速的偏差输出是q轴电流指令。具体到Simulink里的实现转速指令发生器用阶跃斜坡限制模块设定加速度限制为500rpm/s从0加速到12000rpm需要24秒。这个加速度限制是根据电机最大转矩电流和飞轮转动惯量计算出来的最大加速转矩大约140Nm折算到转速变化率就是约600rpm/s留10%裕量后取500rpm/s。转速环PI含积分限幅和输出限幅iq指令上限设为200A对应约280Nm的峰值转矩下限-200A。电流内环d轴电流指令恒为0表贴式PMSM最大转矩电流比控制退化为id0控制q轴电流跟踪转速环输出。两个电流环各自使用PI控制器输出为dq轴电压。解耦补偿在电流环输出处加入交叉解耦项补偿dq轴之间的耦合电动势。解耦项的计算方法Vd_comp -ω_e * Lq * iqVq_comp ω_e * (Ld * id Ψf)。这一项不加的话高速时dq轴电流控制会互相干扰动态响应会明显变差。整套逻辑跑下来充电过程三相电流波形呈逐步增大的正弦波包络转速线性上升转矩电流维持在设定值附近直流母线功率从零爬升到稳定值再下降整体过程清晰稳定。3.3 放电过程的控制切换与能量回馈实现放电过程是整个模型里最容易出问题的地方。控制目标从转速恒定切换到母线电压恒定本质上是从调速模式切到电压控制模式两者需要平滑衔接。我的实现方案是把转速环的PI控制器旁路电机进入发电模式外环改为直流母线电压环。电压环的输入是母线电压给定与实际母线电压的偏差输出是q轴电流指令发电工况下iq为负值内环电流环仍然不变跟踪这个负的电流指令电机转矩反向转速下降动能转化为电能回馈到直流母线。切换的关键细节是切换瞬间的电流指令缓冲。当发电控制使能时电压环PI的输出初始值不能从零开始而应该从当前iq值开始否则控制器积分器充电过程会造成数十毫秒的电流断档母线电压会瞬间跌落。我在模型里用了一个带外部初值设定的积分器切换信号到来时把积分初值设置为切换前的iq值这样就实现了无扰切换。放电工况的功率大小由负载决定。模型中用了一个可变功率负载来模拟功率吸收需求母线电压闭环会自动增大放电电流来维持电压。如果飞轮转速下降到最低允许转速模型设定为6000rpm触发欠速保护停止放电飞轮系统进入待机状态。这个最低转速与最大放电深度相关放电深度越深能量释放越多但转速过低会导致母线电压无法维持需要提前切出。3.4 转速-功率协调与S型速度曲线规划基础的转速斜坡指令可以做但充放电过程更柔顺的做法是用S型速度曲线。原因在于飞轮系统的机械轴承和轴系在加速度突变时会有冲击长期运行会造成疲劳损伤。我在模型的速度指令发生器里做了一个简单的S型曲线规划加速过程中先以线性增加的加速度启动中间恒加速度末端加速度线性减小到零。对应到转速曲线上就是一段S型过渡。这个曲线可以通过Simulink里的信号生成模块实现或者直接用Matlab函数写一段数值积分。S型曲线的好处不仅是机械冲击小对电网侧的能量冲击也小。充电开始时如果直接给满转矩直流母线相当于短时功率冲击对前端整流器和电网都不友好。而S型曲线让功率缓慢爬升对系统友好得多。我做了一个简化处理设定加速时间为24秒但把前3秒作为加速度递增阶段后3秒作为加速度递减阶段。仿真波形上可以看到转速指令在起始和结束阶段呈现出平滑的弧度而不是折线拐角。这个细节虽然不影响控制稳定性但对后续的并网功率质量很有意义。4. 常见问题与排查技巧4.1 电流环发散或振荡怎么办最典型的现象是仿真开始后iq或者id电流振荡发散直接飞掉。排查顺序是这样的先看坐标变换方向尤其是转子位置角的正负号。把d轴方向和转子磁极方向对齐如果符号反了电流反馈就变成正反馈系统必发散。验证方法给一个常数iq指令看电机是否向预期方向加速同时看id是否保持在零附近。再看电流环PI参数如果计算出来的Kp_i和Ki_i是按照连续域算的而Simulink里用的是离散PI模块需要先做离散化转换。直接用连续域的增益放到离散控制器里采样频率不够高时系统会振荡。最后看解耦项有没有加对。解耦项的符号错了会引入额外的正反馈。常见错误是交叉项漏乘转速电角速度导致低速时正常、高速时振荡。4.2 转速环响应慢且超调转速环的PI参数按带宽50Hz计算但实际效果发现转速跟踪慢且到达目标转速时超调量达到3%-5%。这个问题多半是电流环限幅提前饱和了。加速过程中转速误差大转速环PI输出想要一个很大的iq指令但限幅器把它截在200A这时候转速环相当于开环控制等转速接近目标时误差变小PI才从饱和状态退出来。退出过程中积分器还在持续累积所以会产生超调。解决办法有两个一是给转速环PI的积分器做抗积分饱和设计当输出达到限幅值时冻结积分器二是在转速接近目标时提前切换控制策略比如在目标转速的95%开始减小加速度。我在模型里用的是积分限幅条件冻结切换后超调量压到1%以内。4.3 充电到放电切换瞬间母线电压跌落这个问题在不少做飞轮储能仿真的朋友那里反复出现。现象是母线电压在切换后的前50ms内跌落了20V以上然后才慢慢恢复。原因是切换时刻电压环PI的初始输出和当前工况不匹配。调试方法把放电开始的时刻记录下来观察切换前瞬间的iq值、母线电压和直流负载功率。然后把电压环PI的积分初始值设置成切换前iq对应的值或者切换前的实际占空比对应的等效值切换就能做到基本无扰。另一种做法是切换后短时间内保持切换前的电流指令用软启动的方式让电压环慢慢接管也能改善。4.4 SVPWM电压指令饱和高速运行时反电动势接近直流母线电压能输出的极限电流环输出的dq电压指令会超出SVPWM的线性调制范围导致实际输出电压不够电流跟踪失败。这个现象通常发生在额定转速以上。处理办法是给电流环输出电压加限幅保证电压矢量幅值在最大线性范围内。同时需要考虑是否需要弱磁。如果飞轮的工作转速不超过额定转速不加弱磁也能跑如果要超速运行到120%额定转速就必须引入弱磁控制。我的模型里加了一个简单的弱磁控制器当电压指令超过限幅值时减少d轴电流指令负值通过d轴去磁效应来降低反电动势。弱磁控制器在调试时要格外小心负的d轴电流很容易把电流限幅值占满反而挤占q轴电流空间。一般设一个优先级在总电流限幅的约束下q轴电流优先弱磁电流在剩余空间里调整。4.5 仿真速度过慢的优化经验高速PMSM仿真尤其是带开关器件和PWM调制的时候仿真速度很容易让人崩溃。一次完整的充放电仿真要模拟几十秒的物理时间如果步长太小几个小时都跑不完。我优化之后的配置是功率电路和解算器用离散模式采样时间设为1μs对应10kHz开关频率的PWM分辨率和足够的电流环带宽机械部分转速、位置用较慢的采样率模块单独处理。Simulink里可以设置不同的采样时间这样做之后仿真速度提升了3-5倍波形精度损失很小。另外就是避免在循环里面加连续积分器用离散积分器替代数值稳定性也更好。如果模型确实太大还可以考虑把充放电过程拆分成两个独立的模型来跑一个是充电过程模型一个是放电过程模型中间用初始状态继承来衔接。这样每个模型的计算量小调试更方便最后再把两个模型合并验证。4.6 快速排查表我把做这个仿真过程中最容易踩的坑整理成一个速查表按优先级排序排查现象可能原因排查手段仿真初瞬即发散坐标变换符号错误、电流反馈极性反检查Park变换角度符号断开转速环只测电流环电流振荡电流环PI增益过高、解耦项符号错误降低Kp_i关闭解耦项验证转速超调大转速环积分饱和、限幅过低给PI加抗饱和提高限幅到1.5倍额定电流母线电压跌落模式切换无初始值匹配切换前记录iq预置电压环PI初值高速时转矩不足电压饱和、未加弱磁检查SVPWM线性范围启用弱磁控制仿真极慢连续求解器步长过小改用离散模式合理设置采样时间5. 充电阶段深度解析与功率转矩特性5.1 充电过程的三段式控制策略把充电过程进一步拆开看完整的充电控制其实是三个阶段每个阶段对控制器有不同的要求第一阶段是最大转矩加速阶段。转速从静止到某一预设定转速的过程中转速环一直处于饱和状态iq指令稳定在限幅值电机输出恒定最大转矩转速线性上升。这个阶段电流环的跟踪性能和直流母线的供电能力是主要约束。如果直流母线电压不够高大电流会造成母线电压跌落反而限制了最大转矩输出。第二阶段是恒功率加速阶段。随着转速升高反电动势增大电压余量逐渐减少电机从恒转矩区过渡到恒功率区。这时候q轴电流会自然回落因为同样的电压极限下高转速能输出的最大转矩变小了。我在模型里观察到的现象是转速达到约8000rpm后iq指令从200A开始回落转矩功率乘积基本保持恒定。第三阶段是转速维持阶段。达到目标转速后转速环从饱和退出iq减小到维持摩擦力矩和风阻损耗的数值。对于高速飞轮系统这个维持电流并不为零因为高速下的风阻损耗虽然真空环境可以降低但不完全为零轴承摩擦损耗同样存在。充电策略里还有一个放电深度管理的问题需要考虑也就是不能让飞轮一直充到过高的转速要给最高转速留一个保护余量通常控制在额定转速的105%以内。5.2 充电功率特性与直流侧需求充电过程的功率特性对前端变换器选型有直接影响。最大转矩加速阶段机械功率Pm T * ω转速上升的同时功率上升直流侧输入功率从零平滑爬升到最大值。这个特性对前级AC/DC变换器比较友好因为不存在阶跃式的功率冲击。我算过一组数据验证功率匹配关系最大转矩140Nm在8000rpm时机械功率Pm 140 * (2π * 8000/60) ≈ 117kW。考虑电机和逆变器损耗后直流侧需要输入的功率大约在120kW左右。这就是为什么这台55kW额定功率的电机在峰值转矩状态下短时输入功率会明显超过额定值这是飞轮快速充能的特性需求也是电机和功率器件需要短时过载能力的原因。模型里把直流母线的电容和电压源内阻也做了建模目的是验证大电流充电时母线电压的波动幅度。实际仿真下来700V母线在峰值充电电流条件下电压波动控制在3%以内如果有超出需要增加母线电容或前级响应速度。5.3 充电结束的转速保持策略飞轮储能系统充满电之后不是立刻放电而是要维持在额定转速待命。这个阶段有个参数值得关注转速保持的精度和待机损耗。转速保持仍然用转速环但这时候转速指令保持不变干扰来自风阻和轴承摩擦的缓慢变化。转速环的积分作用可以把稳态误差消除到接近零但代价是持续的电流输出。我做了一个估算待机时所需的维持转矩大概是额定转矩的1%-2%对应的iq维持电流约2-4A直流侧输入功率约1-2kW。这个待机损耗决定了飞轮储能的自放电时间常数。为了降低待机损耗工程上会在真空腔体内维持高真空度把风阻降到忽略不计。仿真模型里我给了一个风速相关的阻尼系数并且在模型里留了接口可以在不同真空等级下切换。实际调试中发现如果真空度足够好气压低于1Pa风阻转矩可以忽略待机时只需要补偿轴承损耗这样储存的能量可以保持较长时间这也是飞轮储能系统适合做短时高频次储能的根本原因。6. 放电阶段控制精讲与母线支撑逻辑6.1 母线电压闭环的实现逻辑放电时飞轮储能系统的角色从电压接受方变成电压支撑方。在独立直流母线系统中飞轮系统通过控制逆变器输出功率来维持母线电压稳定。母线电压环的传递链是Udc_ref与Udc_feedback的误差经过PI控制器输出iq_ref电流环跟踪iq电机发电产生负转矩转速下降机械能转化为电能输入直流母线。这里有一个因果关系值得注意母线电压环控制的本质是功率平衡并不是直接控制电压。当外部负载功率增大时母线电压会先下降一个小的偏差电压环感知到偏差后增大电流输出补充功率缺口电压恢复。所以这个闭环的物理延时取决于母线电容的大小和电压环带宽。放电时电流环的dq轴控制和充电时是一样的只是q轴电流参考值是负的。为保护电机和逆变器放电电流同样需要限幅。因为发电工况下如果突然甩负荷机械功率无法瞬间降为零会产生过电压。所以模型中给放电工况做了直流母线过压保护当母线电压超过740V时立即封锁PWM并切到能耗制动模式防止飞轮系统失控。6.2 放电深度与转速下限保护放电过程飞轮转速从12000rpm下降存储动能E 0.5 * J * (ω²_max - ω²_min)。这个公式说明放电深度对释放能量有平方关系影响转速从12000rpm放到6000rpm可以释放75%的动能而放到9000rpm只能释放43.75%的动能。但转速下限不能为零原因有三条第一永磁同步电机反电动势与转速成正比转速过低时反电动势太小无法维持母线电压发电能力丧失。第二转子永磁体在高速旋转下如果突然降到很低转速轴承润滑油膜和动压轴承的支承可能失效机械安全性存疑。第三飞轮系统的控制策略通常在深度放电末端需要切换回电机模式来维持飞轮基础转速避免停转后需要重新从零加速。模型里我把转速下限设为6000rpm对应最大放电深度50%。这个值对应放电结束后的母线电压维持能力足够反电动势还有额定值的一半而且再次充电的加速时间相对较短。实际工程中这个下限要根据具体系统和工况来设定有的飞轮储能系统为了追求更大的放电能量会把下限压到额定转速的40%。6.3 放电工况的功率平滑策略模拟真实负载工况时负载往往是波动的母线电压闭环需要应对这种波动。为了验证模型的动态响应能力我设置了一个可编程负载序列模拟功率在20kW到100kW之间阶跃变化。实测母线电压的恢复时间大约在10-20ms之间超调量小于5%。这个动态响应指标对于飞轮储能系统来说是可以接受的。如果把母线电压环带宽再调高理论上可以获得更快的恢复速度但过高的带宽会和电流环带宽接近引发控制环路之间的耦合振荡。在功率平滑策略上我建议使用一阶低通滤波对负载功率参考值做平滑处理限制功率变化率。这样做的好处是可以减小储能单元的冲击电流延长功率器件寿命代价是动态响应慢一些。对于电网调频这类需要快速响应的应用场景平滑时间常数要取得很小或者干脆不用滤波直接让电压环快调。7. 系统性能评估与待机损耗分析7.1 充放电效率的仿真评估方法储能系统的效率是最核心的性能指标。仿真中评估效率的方法是记录整个充放电过程中的能量值变化然后对比放电输出能量和充电输入能量的比值。我在模型中加了若干个积分器来记录能量。充电过程中直流侧的输入能量通过积分Udc * Idc计算放电过程中负载消耗的能量通过积分负载功率得到。仿真结束后用能量比值就能算出系统效率。以我这组参数为例充放电深度50%的循环效率仿真结果大约是88%-90%。这个效率的量级比较合理铜损、铁损、开关损耗、摩擦损耗都算进去了。如果把中间待机时间拉长待机损耗会让循环效率明显下降这也解释了为什么飞轮储能适合高频次短时储能而不是长时间储能。7.2 待机损耗的来源与参数敏感性待机损耗主要有三块风阻损耗、轴承摩擦损耗、电机铁损和控制系统的自耗电。仿真模型中风阻损耗用转速相关的阻尼系数模拟轴承摩擦损耗用恒定摩擦转矩模拟铁损通过电机模块内部参数计算。我特意做了参数敏感性分析看待机损耗对哪些参数最敏感。结论是风阻损耗对转速最敏感转速降10%风阻损耗下降约30%。所以如果应用场景允许降低待机转速比维持满转速待机要划算得多哪怕需要快速升速时再加速回来。这也是有些飞轮系统采用低转速待机、高速充能策略的原因。7.3 模型的扩展方向这套模型搭建完成之后可以往几个方向做扩展。第一个方向是接入上层调度。把充放电控制封装成一个可接受功率指令的储能单元模型然后放到微电网系统仿真里验证飞轮储能在削峰填谷、调频辅助服务等场景中的作用。第二个方向是并联扩容。实际工程中单台飞轮机组功率不够时会把多台飞轮并联接入同一个直流母线。并联系统的控制难点在于多台飞轮的荷电状态不一致需要做功率分配策略。基于这套模型可以复制出多台机组在功率分配算法上做研究。第三个方向是做硬件在环。用实时仿真机运行电机和主电路模型控制器跑在实际的控制板上验证实际控制器代码的正确性。硬件在环的重点就是模型离散化程度和I/O配置基于这套模型做离散化改造是可行的。7.4 模型验证与数据校准的实操心得仿真模型做得再精细终究是简化模型要保证仿真结果对工程有参考价值必须做验证。我的做法是分三步验证第一步是空载验证。把飞轮负载去掉只保留电机空载运行对比仿真反电动势波形和电机参数表中给出的数值是否一致。这一步如果对不上先查坐标变换和参数换算。第二步是阶跃响应验证。给一个小的转速阶跃记录iq响应波形和转速响应波形检查超调量、上升时间是否和理论计算一致。这一步验证PI参数整定是否正确。第三步是能量闭环验证。做完整的充放电循环检查初始转速和最终转速是否一致母线电压是否维持稳定能量积分结果是否符合效率预期。这一步验证系统级的能量管理是否正确。三次验证都通过之后仿真模型才算是可信度较高的工具后面调控制策略、做工况分析才能有依据。盲目相信仿真结果的后果我在做硬件在环时遇到过控制器参数在仿真里表现很好上到实际系统就振荡原因就是仿真模型忽略了数字控制的离散延迟。所以做仿真的时候要刻意把离散采样、PWM死区这些非理想因素加进去。最后再说一个容易被忽略的细节飞轮储能系统的仿真最好能加入低转速下的数值稳定处理。有些求解器在转速接近零时电机反电动势接近零状态变量可能出现数值奇异。解决的办法是给机械方程式加一个很小的阻尼项或者在仿真初始化时给出一个初始转速值比如100rpm避免从完全静止开始仿真时的数值抖动。这个小技巧能让仿真的启动过程顺畅很多尤其对后续要扩展到硬件在环的朋友来说这个处理几乎是必须的。
