飞轮储能系统一直是被低估的储能方案它不依赖化学反应靠的是“电-动能-电”的转换瞬时功率大、循环寿命长、响应快这几年在数据中心的应急供电、电网一次调频、轨道交通再生制动回收这些场景里越来越常见。我在做某飞轮储能项目时驱动电机选的是永磁同步电机PMSM而不是更常见的地面异步电机因为飞轮储能的工况实在太适合PMSM了转速高、损耗小、控制准确。这篇文章就把这套飞轮储能系统建模与Simulink仿真的完整思路、电机参数、控制结构和实操步骤全部拆开讲适合正在做飞轮储能、高速电机驱动或者储能变流器仿真的学生和工程师参考也适合刚接触PMSM矢量控制的读者入门。1. 飞轮储能系统整体方案与建模思路1.1 系统组成不只是“一个电机带一个轮子”飞轮储能听起来结构简单实际拆开来看一个完整的电化学储能之外的机械储能系统至少要包括五大环节飞轮转子、电机发电/电动一体、双向变流器、控制器、辅助支撑系统轴承、真空腔、冷却等。在仿真里我一般把重点放在前四部分辅助系统用等效摩擦系数和阻尼折算到机械侧就可以。飞轮转子负责储存动能它的核心参数是转动惯量J和转速ω电机负责电能和动能的双向转换充电时它是电动机把电网或母线电能变成飞轮转速放电时它是发电机把飞轮动能变成电能送回母线。双向变流器负责电机侧和电网/直流母线之间的电压电流变换控制器则承担模式判断、转速闭环、电流闭环和PWM脉冲生成。这个结构对应到Simulink仿真里就是五个子系统电机模型、逆变器模型、坐标变换、双闭环控制器、PWM调制。讲到后面我会专门给出搭建顺序这里先记住一个大原则飞轮储能建模的本质是“机械侧转动惯量”和“电磁侧转矩电流”之间的闭环匹配。机械侧决定转速变化多快电磁侧决定转矩给多大两者通过电机的电磁转矩方程耦合在一起。1.2 为什么选永磁同步电机做飞轮驱动电机很多储能飞轮项目选异步电机或开关磁阻电机各有各的理由但我在这个项目里坚定选PMSM。核心原因有三点。第一是高速性能。飞轮储能密度和转速平方成正比要让同样重量的飞轮存更多能量最直接的办法就是把转速拉高达到每分钟几万转都不是稀罕事。PMSM的转子没有励磁绕组结构坚固很适合高速旋转。所以高速电机结构坚固升速容易这是飞轮储能最需要的。第二是效率。PMSM转子不需要励磁电流因此没有转子铜耗只有定子铜耗和铁耗在充放电反复循环中效率优势非常明显。飞轮储能系统本身的损耗主要来自轴承摩擦、风阻和电机损耗电机效率哪怕只提升两三个百分点对长时间待机运行的工况来说节能效果也很可观。第三是控制简单。表贴式PMSM配合id0控制策略电磁转矩和定子q轴电流严格成正比完全没有异步电机那样的转差控制问题也不像开关磁阻电机那样需要精确处理非线性转矩模型。这就意味着从仿真到工程样机控制代码的迁移路径非常清晰。实际工程中高速PMSM也有一些需要留意的地方比如弱磁控制、永磁体高温退磁风险、高速下反电动势过高等但那是后续优化问题不影响它在飞轮储能中的主流地位。1.3 能量关系与转速SOC用转速来“管”电量飞轮储存的动能表达式是E0.5Jω²。设计一套飞轮储能系统不能从最高转速一直放到零转速因为电机反电动势和变流器输出能力会限制低转速段的功率输出所以要定义一个工作转速区间。以我做仿真时用的参数为例飞轮转动惯量取J5kg·m²工作转速区间2000~6000r/min。最高角速度ω_max628.3rad/s最低角速度ω_min209.4rad/s那么可用能量是ΔE0.5×J×ω_max²-ω_min²0.5×5×628.3²-209.4²≈877kJ换算成度电约0.244kWh。也就是说一台5kW的电机可以在约3分钟内把这部分能量放完或者用5kW功率充约3分钟充满。这个能量等级适合给小型UPS、线路调频这类应用做概念验证不适合做长时间能量支撑。正因为储能和转速的平方相关飞轮的荷电状态SOC就不能像电池那样按安时累积必须用转速平方归一化。我习惯用下面这个公式SOC(ω²-ω_min²)/(ω_max²-ω_min²)×100%。这个定义的好处是直观SOC100%对应最高转速6000r/minSOC0对应最低转速2000r/min控制器只要测量转速就能知道还能放多少能量。很多仿真初学者直接把SOC按转速线性处理低速段会明显低估剩余能量这也是我在实际调试中踩过的一个坑。2. 永磁同步电机建模与磁场定向控制参数2.1 dq坐标系下的电机方程与id0策略PMSM在三相静止坐标系下的方程是时变的直接控制三相交流量很麻烦。实际建模几乎都走坐标变换路线先把三相ia、ib、ic通过Clarke变换变成两相静止αβ轴电流再利用转子位置角θ做Park变换得到同步旋转dq坐标系下的id、iq。这样交流量变成直流量控制器的设计难度大幅下降。dq坐标系下电压方程是udRs×idLd×d(id)/dt-ωe×Lq×iq uqRs×iqLq×d(iq)/dtωe×(Ld×idψf)其中ωe是电角速度ωep×ωp是极对数ψf是永磁磁链。电磁转矩方程为Te1.5×p×[ψf×iq(Ld-Lq)×id×iq]如果电机是表贴式结构LdLq转矩表达式简化为Te1.5×p×ψf×iq转矩完全由iq控制这就是id0控制的物理基础。即使内置式PMSM存在磁阻转矩飞轮储能这种恒功率区间明确的负载id0策略依然最简单可靠。我在仿真中设置的一组典型PMSM参数如下参数数值说明额定功率5kW飞轮充放电功率等级额定转速6000r/min对应飞轮最高工作转速极对数p2高速电机极对数通常少定子电阻Rs0.05Ω绕组铜阻d轴电感Ld0.5mH表贴式时LdLqq轴电感Lq0.5mH同上永磁磁链ψf0.2Wb决定反电动势大小最大电流Imax40A过载能力约3倍额定用转矩方程校核一下额定5kW在6000r/min时负载转矩约7.96N·m对应q轴电流7.96/(1.5×2×0.2)≈13.3A在这个参数下电流裕量很充足。最大电流40A换算成最大转矩约24N·m足够做快速加减速和短时过载功率输出。2.2 双闭环控制结构与PI参数整定飞轮电机控制我采用转速外环、电流内环的双闭环结构。充电时转速环给定目标转速比如6000r/min放电时根据母线功率需求生成转速偏差转速环输出q轴电流参考值电流环再控制实际iq跟踪参考。id的参考值恒为0。电流环的对象可以看作一阶惯性环节G(s)1/(LsR)用零极点对消法整定PI参数。设定期望电流环时间常数τ_i1ms则KpL/τ_i0.5mH/1ms0.5 KiR/τ_i0.05Ω/1ms50这里Kp单位是V/AKi单位是V/(A·s)在离散化后要注意把Ki换算成增量形式。电流环带宽约1000rad/s量级开关频率10kHz时这个带宽是合理的。转速环整定要慢得多我一般把电流环带宽设为转速环的5到10倍以上。转速环对象近似为1/(Js)考虑转矩常数kt1.5×p×ψf0.6N·m/A按期望转速环时间常数τ_n0.05s粗略计算Kp≈J/(kt×τ_n)5/(0.6×0.05)≈167Ki按10倍带宽关系再调。实际仿真必须在限幅条件下重新微调不要迷信计算值。这里还要强调一件事dq轴之间存在交叉耦合项-ωe×Lq×iq和ωe×(Ld×idψf)。低速时影响小高速飞轮转速达到每分钟几千转时这两项会明显影响电流环动态。我的做法是在电流环PI输出后叠加前馈补偿项ud_ff-ωe×Lq×iq uq_ffωe×(Ld×idψf)然后把前馈值加到PI输出上。这样高速段iq电流跟踪能力会明显改善转速上升时电流不容易发散。3. 从机械能到电能的转换飞轮-变流器-母线功率平衡3.1 充电、待机、放电三种运行模式飞轮储能不像电池它的模式切换必须考虑转速约束。我在仿真里用有限状态机实现三种运行模式。模式触发条件控制目标电机工作状态充电SOC100%有功需求为正按给定功率把飞轮加速到目标转速电动机电磁转矩为正待机转速达到目标无能量交换需求维持转速恒定抵消摩擦损耗电动机输出小转矩补偿损耗放电有功需求为负或母线电压偏低按给定功率释放动能飞轮减速发电机电磁转矩为负实际仿真中很多人忽略待机模式导致充完电后电机转矩直接掉到零飞轮转速会因摩擦缓慢下降等下一次需要放电时转速已经不够。仿真里待机模式其实就是让转速环继续工作参考值保持当前转速PI控制器自动输出一个很小的正转矩来平衡摩擦阻力。放电模式有两种实现思路。第一种是把转速参考值直接设低一些比如目标从6000r/min降到2000r/min转速环自然输出负转矩电机进入发电状态。第二种更接近实际工程用直流母线电压外环替代转速外环母线电压偏低时电压环自动给出负的q轴电流参考电机回馈能量。我在初始做概念仿真时优先用第一种逻辑清晰、容易调试等验证完电机控制后再切换到第二种研究母线电压动态。3.2 直流母线功率平衡与电容选型飞轮的机械功率表达式是PeTe×ω。充电时Te0电能变成动能放电时Te0动能变回电能。高速飞轮在发电状态下机械功率很容易超过kW级为了实现功率平衡仿真中直流母线就不能简单当成理想电压源还要考虑母线电容的动态。母线电容的电压方程为C×d(Udc)/dtI_in-I_out其中I_in是变流器从飞轮侧回馈到母线的电流I_out是负载或电网吸收的电流。放电瞬间如果飞轮回馈功率大于负载功率多余的能量会灌入母线电容导致Udc升高相反如果回馈功率不足Udc跌落。所以母线电容起到“缓冲池”作用。我做母线稳压放电仿真时电容取值2200μF直流母线电压等级540V。用5kW放电功率估算如果没有外部负载母线电压从540V升高到安全上限的暂态时间大约几十毫秒量级这个时间必须大于控制环的响应时间否则过压保护会触发。这也是为什么实际飞轮储能并网必须用跟网型变流器把能量稳定送入电网而不能只对母线放电。仿真早期为了简化我会先接一个理想直流电压源验证电机电动/发电能力后面再换成电容加负载模型建议你也按这个顺序做。3.3 模式切换平滑处理避免转速突变模式切换最容易出现的问题是切换瞬间转矩、电流跳变严重时会导致整个模型发散。我从仿真里总结出的经验是三条。一是参考值一律经过斜坡或一阶惯性环节后再进入转速环不要让转速参考值在切换瞬间从6000r/min直接跳到2000r/min。斜坡速率按电机最大转矩能力设定考虑J5kg·m²最大转矩24N·m最大转速变化率约为24/54.8rad/s²对应约45.8r/min/s。如果工况要求5秒内从6000r/min降到2000r/min平均降速约800r/min/s已经超过这个能力必须限制放电功率否则转速环会饱和。二是切换瞬间把积分项做初始化或复位处理。双闭环PI控制器在不同模式下的输出初值不一致直接切换会产生积分器跳变。我在状态机切换动作里加了一步切换时把转速环积分器清零让输出从零开始按新的参考值重新建立代价是切换过程会有几百毫秒的动态调整但换来了稳定。三是如果切换过于频繁可以考虑不用硬切换而是始终用“目标功率”作为统一外环参考充电时功率为正放电时功率为负。这样状态机退化成功率调度器只需要把功率参考值映射成q轴电流给定切换点处自然连续不会出现模式切换震荡。4. Simulink仿真模型搭建步骤与关键模块配置4.1 五个子系统的搭建顺序我使用Simulink里Simscape Electrical的电气库搭建这个模型分成五个子系统依次搭建。第一步搭PMSM与逆变器。从库中拖出Permanent Magnet Synchronous Machine模块按前面表格设置电机参数输出端口接机械负载。机械输入侧直接设置转动惯量J5kg·m²、摩擦系数B取一个很小的值例如0.001N·m·s/rad这样飞轮的转动惯量特性就集成在电机模型里了。逆变器用Universal Bridge模块桥臂数量选3开关器件选IGBT/Diode直流侧接电压源或母线电容。第二步搭坐标变换。用两个Subsystem实现Clarke和Park变换输入是三相电流和转子位置角θ。电机模块通常输出转子机械角度需要乘以极对数p得到电角度。第三步搭电流环。把id的参考值设为0iq_ref由外环输出接入电流误差分别进两个PI控制器PI输出叠加前馈解耦项ud_ff、uq_ff得到dq轴的电压给定。第四步搭外环。从电机模块取转速wm和目标转速做差经转速PI后输出iq_ref同时加限幅到±40A这是整个控制里面最关键的限幅环节。第五步搭PWM与脉冲。把dq轴电压给定做Park反变换得到三相调制波与10kHz三角载波比较生成六路PWM脉冲送给Universal Bridge的门极。为节约仿真时间我第一版用SPWM正弦调制调制度超过0.9以后才考虑SVPWM。4.2 求解器设置与仿真加速技巧飞轮储能系统横跨电磁和机械两个时间尺度。电流环响应是毫秒级PWM开关是10kHz即0.1ms周期而飞轮转速爬升需要几十秒甚至上百秒如果老老实实把整个物理过程仿完仿真步长选1e-5秒仿真时长100秒那是1000万步效率低得难以接受。我用的加速技巧有三个。第一仿真测试阶段把转动惯量J从5kg·m²临时改成0.5kg·m²这样转速爬升时间缩短到原来的1/10控制策略验证通过后再改回真实值做精细仿真。第二求解器选择离散定步长步长取2e-5秒和10kHz开关周期匹配既能观察电流纹波又不浪费算力如果只看宏观能量变化可以进一步用5e-5秒配合均值模型。第三尽量用速度源或角度源代替部分反馈回路时要确认反电动势影响仍然存在不要为了加速把物理过程简化掉。4.3 监控量与SOC估算模块仿真调试过程中我最常看的波形是实际转速wm、三相电流ia/ib/ic、dq轴电流id/iq、电磁转矩Te、直流母线电压Udc。另外一定要加一路SOC估算用转速经Function模块计算(ω²-ω_min²)/(ω_max²-ω_min²)并与模式状态信号一起送入Scope实时观察。一个很实用的小技巧是把目标转速、实际转速、SOC放在同一个Scope里显示能直接判断充电/待机/放电各阶段是否有逻辑漏洞。比如充电结束SOC应该等于100%如果SOC还在缓慢下降说明待机模式没有正确补偿摩擦损耗。这种问题单看转速可能不能立刻察觉但SOC波形一对比就很明显。5. 典型工况仿真与问题排查经验5.1 充电-待机-放电循环仿真我用一组典型工况展示流程。初始飞轮转速2000r/minSOC0%充电指令下达后目标转速以斜坡从2000r/min升至6000r/min参考斜坡斜率取40r/min/s约140秒充满。实际调试时因为转动惯量缩小到0.5kg·m²充满时间约14秒便于观察。充电阶段转速环输出正q轴电流电磁转矩约10N·m三相电流幅值随转速变化基本稳定。待机阶段保持6000r/min约2秒此时q轴电流降到接近摩擦补偿的小值电磁转矩几乎为0母线电流也很小。然后下达放电指令目标转速斜坡降到3000r/min转速环输出负转矩电机进入发电状态三相电流相位与反电动势相反机械能持续转化为电能送回直流母线。如果母线侧接的是耗能负载放电期间母线电压能稳定在540V附近如果只接电容母线电压会明显升高这正好验证了功率平衡设计的正确性。5.2 问题排查速查表以下是我在这个飞轮储能仿真项目中实际遇到的问题和排查结果整理成表格供参考。故障现象最可能原因处理办法启动瞬间电流超限初始转子位置角设置不当dq电流理解错误初始化时把转子角置0并给电流环加限幅充电时转速爬升太慢转动惯量太大、转矩限幅太小或摩擦系数偏大检查J和B参数适当放宽电流限幅到40A高速区q轴电流跟踪不上dq耦合前馈缺失调制比饱和加入前馈解耦项必要时改SPWM为SVPWM放电瞬间母线电压飙升母线电容太小回馈功率大于负载吸收能力加大电容限制放电功率或接耗能负载模式切换时转速突变转速参考阶跃变化积分器初值不一致参考值加斜坡切换时转速环积分器清零三相电流波形正弦度差PWM载波频率太低或步长太大载波10kHz定步长降到2e-5秒以下SOC计算为负或超过100%转速低于ω_min或高于ω_max时公式溢出SOC估算模块加饱和限幅到0~100%5.3 几点实操心得这套模型我从零开始搭前前后后改了四版才稳定。最大的经验教训是调试顺序不能乱先确保PMSM模型单开环能转动再闭环调电流环最后才合上转速环。如果一开始就是双闭环全套电机参数和PI参数同时出错波形乱成一片根本不知道从哪下手。电机参数我建议除了Rs、Ld、Lq、ψf这四项还要把摩擦系数B记录清楚。很多仿真用默认值0.001实际电机轴承系统在高速段阻尼不是恒定的至少要设成转速的弱函数否则高速段待机功率会失真。仿真模型里可以先用常数近似工程样机测试后再修正。另外PI限幅一定要在所有PI输出端口加上。仿真新手最容易犯的错是把限幅只加在转速环输出上而忽视电流环输出的电压限幅。一旦调制度超过1逆变器输出已经失真电流环PI还在拼命积分最后必然是失控发散。电压限幅我设为±0.45倍母线电压对应SPWM的线性调制区。最后再分享一个小技巧把整台飞轮系统的参数做成一个.m脚本用Matlab的变量区统一管理。每个子系统里的PI参数、限幅值、电容值都引用工作区变量调整参数只用改脚本文件不用拆模型连线。我后来跑参数扫描时靠这个习惯至少节省了三分之一的调试时间。下一步如果要接电网侧做调频仿真这套模型只需把直流电压源替换成网侧变流器加入锁相环和功率外环就能扩展底子已经打得比较稳了。
飞轮储能系统一直是被低估的储能方案它不依赖化学反应靠的是“电-动能-电”的转换瞬时功率大、循环寿命长、响应快这几年在数据中心的应急供电、电网一次调频、轨道交通再生制动回收这些场景里越来越常见。我在做某飞轮储能项目时驱动电机选的是永磁同步电机PMSM而不是更常见的地面异步电机因为飞轮储能的工况实在太适合PMSM了转速高、损耗小、控制准确。这篇文章就把这套飞轮储能系统建模与Simulink仿真的完整思路、电机参数、控制结构和实操步骤全部拆开讲适合正在做飞轮储能、高速电机驱动或者储能变流器仿真的学生和工程师参考也适合刚接触PMSM矢量控制的读者入门。1. 飞轮储能系统整体方案与建模思路1.1 系统组成不只是“一个电机带一个轮子”飞轮储能听起来结构简单实际拆开来看一个完整的电化学储能之外的机械储能系统至少要包括五大环节飞轮转子、电机发电/电动一体、双向变流器、控制器、辅助支撑系统轴承、真空腔、冷却等。在仿真里我一般把重点放在前四部分辅助系统用等效摩擦系数和阻尼折算到机械侧就可以。飞轮转子负责储存动能它的核心参数是转动惯量J和转速ω电机负责电能和动能的双向转换充电时它是电动机把电网或母线电能变成飞轮转速放电时它是发电机把飞轮动能变成电能送回母线。双向变流器负责电机侧和电网/直流母线之间的电压电流变换控制器则承担模式判断、转速闭环、电流闭环和PWM脉冲生成。这个结构对应到Simulink仿真里就是五个子系统电机模型、逆变器模型、坐标变换、双闭环控制器、PWM调制。讲到后面我会专门给出搭建顺序这里先记住一个大原则飞轮储能建模的本质是“机械侧转动惯量”和“电磁侧转矩电流”之间的闭环匹配。机械侧决定转速变化多快电磁侧决定转矩给多大两者通过电机的电磁转矩方程耦合在一起。1.2 为什么选永磁同步电机做飞轮驱动电机很多储能飞轮项目选异步电机或开关磁阻电机各有各的理由但我在这个项目里坚定选PMSM。核心原因有三点。第一是高速性能。飞轮储能密度和转速平方成正比要让同样重量的飞轮存更多能量最直接的办法就是把转速拉高达到每分钟几万转都不是稀罕事。PMSM的转子没有励磁绕组结构坚固很适合高速旋转。所以高速电机结构坚固升速容易这是飞轮储能最需要的。第二是效率。PMSM转子不需要励磁电流因此没有转子铜耗只有定子铜耗和铁耗在充放电反复循环中效率优势非常明显。飞轮储能系统本身的损耗主要来自轴承摩擦、风阻和电机损耗电机效率哪怕只提升两三个百分点对长时间待机运行的工况来说节能效果也很可观。第三是控制简单。表贴式PMSM配合id0控制策略电磁转矩和定子q轴电流严格成正比完全没有异步电机那样的转差控制问题也不像开关磁阻电机那样需要精确处理非线性转矩模型。这就意味着从仿真到工程样机控制代码的迁移路径非常清晰。实际工程中高速PMSM也有一些需要留意的地方比如弱磁控制、永磁体高温退磁风险、高速下反电动势过高等但那是后续优化问题不影响它在飞轮储能中的主流地位。1.3 能量关系与转速SOC用转速来“管”电量飞轮储存的动能表达式是E0.5Jω²。设计一套飞轮储能系统不能从最高转速一直放到零转速因为电机反电动势和变流器输出能力会限制低转速段的功率输出所以要定义一个工作转速区间。以我做仿真时用的参数为例飞轮转动惯量取J5kg·m²工作转速区间2000~6000r/min。最高角速度ω_max628.3rad/s最低角速度ω_min209.4rad/s那么可用能量是ΔE0.5×J×ω_max²-ω_min²0.5×5×628.3²-209.4²≈877kJ换算成度电约0.244kWh。也就是说一台5kW的电机可以在约3分钟内把这部分能量放完或者用5kW功率充约3分钟充满。这个能量等级适合给小型UPS、线路调频这类应用做概念验证不适合做长时间能量支撑。正因为储能和转速的平方相关飞轮的荷电状态SOC就不能像电池那样按安时累积必须用转速平方归一化。我习惯用下面这个公式SOC(ω²-ω_min²)/(ω_max²-ω_min²)×100%。这个定义的好处是直观SOC100%对应最高转速6000r/minSOC0对应最低转速2000r/min控制器只要测量转速就能知道还能放多少能量。很多仿真初学者直接把SOC按转速线性处理低速段会明显低估剩余能量这也是我在实际调试中踩过的一个坑。2. 永磁同步电机建模与磁场定向控制参数2.1 dq坐标系下的电机方程与id0策略PMSM在三相静止坐标系下的方程是时变的直接控制三相交流量很麻烦。实际建模几乎都走坐标变换路线先把三相ia、ib、ic通过Clarke变换变成两相静止αβ轴电流再利用转子位置角θ做Park变换得到同步旋转dq坐标系下的id、iq。这样交流量变成直流量控制器的设计难度大幅下降。dq坐标系下电压方程是udRs×idLd×d(id)/dt-ωe×Lq×iq uqRs×iqLq×d(iq)/dtωe×(Ld×idψf)其中ωe是电角速度ωep×ωp是极对数ψf是永磁磁链。电磁转矩方程为Te1.5×p×[ψf×iq(Ld-Lq)×id×iq]如果电机是表贴式结构LdLq转矩表达式简化为Te1.5×p×ψf×iq转矩完全由iq控制这就是id0控制的物理基础。即使内置式PMSM存在磁阻转矩飞轮储能这种恒功率区间明确的负载id0策略依然最简单可靠。我在仿真中设置的一组典型PMSM参数如下参数数值说明额定功率5kW飞轮充放电功率等级额定转速6000r/min对应飞轮最高工作转速极对数p2高速电机极对数通常少定子电阻Rs0.05Ω绕组铜阻d轴电感Ld0.5mH表贴式时LdLqq轴电感Lq0.5mH同上永磁磁链ψf0.2Wb决定反电动势大小最大电流Imax40A过载能力约3倍额定用转矩方程校核一下额定5kW在6000r/min时负载转矩约7.96N·m对应q轴电流7.96/(1.5×2×0.2)≈13.3A在这个参数下电流裕量很充足。最大电流40A换算成最大转矩约24N·m足够做快速加减速和短时过载功率输出。2.2 双闭环控制结构与PI参数整定飞轮电机控制我采用转速外环、电流内环的双闭环结构。充电时转速环给定目标转速比如6000r/min放电时根据母线功率需求生成转速偏差转速环输出q轴电流参考值电流环再控制实际iq跟踪参考。id的参考值恒为0。电流环的对象可以看作一阶惯性环节G(s)1/(LsR)用零极点对消法整定PI参数。设定期望电流环时间常数τ_i1ms则KpL/τ_i0.5mH/1ms0.5 KiR/τ_i0.05Ω/1ms50这里Kp单位是V/AKi单位是V/(A·s)在离散化后要注意把Ki换算成增量形式。电流环带宽约1000rad/s量级开关频率10kHz时这个带宽是合理的。转速环整定要慢得多我一般把电流环带宽设为转速环的5到10倍以上。转速环对象近似为1/(Js)考虑转矩常数kt1.5×p×ψf0.6N·m/A按期望转速环时间常数τ_n0.05s粗略计算Kp≈J/(kt×τ_n)5/(0.6×0.05)≈167Ki按10倍带宽关系再调。实际仿真必须在限幅条件下重新微调不要迷信计算值。这里还要强调一件事dq轴之间存在交叉耦合项-ωe×Lq×iq和ωe×(Ld×idψf)。低速时影响小高速飞轮转速达到每分钟几千转时这两项会明显影响电流环动态。我的做法是在电流环PI输出后叠加前馈补偿项ud_ff-ωe×Lq×iq uq_ffωe×(Ld×idψf)然后把前馈值加到PI输出上。这样高速段iq电流跟踪能力会明显改善转速上升时电流不容易发散。3. 从机械能到电能的转换飞轮-变流器-母线功率平衡3.1 充电、待机、放电三种运行模式飞轮储能不像电池它的模式切换必须考虑转速约束。我在仿真里用有限状态机实现三种运行模式。模式触发条件控制目标电机工作状态充电SOC100%有功需求为正按给定功率把飞轮加速到目标转速电动机电磁转矩为正待机转速达到目标无能量交换需求维持转速恒定抵消摩擦损耗电动机输出小转矩补偿损耗放电有功需求为负或母线电压偏低按给定功率释放动能飞轮减速发电机电磁转矩为负实际仿真中很多人忽略待机模式导致充完电后电机转矩直接掉到零飞轮转速会因摩擦缓慢下降等下一次需要放电时转速已经不够。仿真里待机模式其实就是让转速环继续工作参考值保持当前转速PI控制器自动输出一个很小的正转矩来平衡摩擦阻力。放电模式有两种实现思路。第一种是把转速参考值直接设低一些比如目标从6000r/min降到2000r/min转速环自然输出负转矩电机进入发电状态。第二种更接近实际工程用直流母线电压外环替代转速外环母线电压偏低时电压环自动给出负的q轴电流参考电机回馈能量。我在初始做概念仿真时优先用第一种逻辑清晰、容易调试等验证完电机控制后再切换到第二种研究母线电压动态。3.2 直流母线功率平衡与电容选型飞轮的机械功率表达式是PeTe×ω。充电时Te0电能变成动能放电时Te0动能变回电能。高速飞轮在发电状态下机械功率很容易超过kW级为了实现功率平衡仿真中直流母线就不能简单当成理想电压源还要考虑母线电容的动态。母线电容的电压方程为C×d(Udc)/dtI_in-I_out其中I_in是变流器从飞轮侧回馈到母线的电流I_out是负载或电网吸收的电流。放电瞬间如果飞轮回馈功率大于负载功率多余的能量会灌入母线电容导致Udc升高相反如果回馈功率不足Udc跌落。所以母线电容起到“缓冲池”作用。我做母线稳压放电仿真时电容取值2200μF直流母线电压等级540V。用5kW放电功率估算如果没有外部负载母线电压从540V升高到安全上限的暂态时间大约几十毫秒量级这个时间必须大于控制环的响应时间否则过压保护会触发。这也是为什么实际飞轮储能并网必须用跟网型变流器把能量稳定送入电网而不能只对母线放电。仿真早期为了简化我会先接一个理想直流电压源验证电机电动/发电能力后面再换成电容加负载模型建议你也按这个顺序做。3.3 模式切换平滑处理避免转速突变模式切换最容易出现的问题是切换瞬间转矩、电流跳变严重时会导致整个模型发散。我从仿真里总结出的经验是三条。一是参考值一律经过斜坡或一阶惯性环节后再进入转速环不要让转速参考值在切换瞬间从6000r/min直接跳到2000r/min。斜坡速率按电机最大转矩能力设定考虑J5kg·m²最大转矩24N·m最大转速变化率约为24/54.8rad/s²对应约45.8r/min/s。如果工况要求5秒内从6000r/min降到2000r/min平均降速约800r/min/s已经超过这个能力必须限制放电功率否则转速环会饱和。二是切换瞬间把积分项做初始化或复位处理。双闭环PI控制器在不同模式下的输出初值不一致直接切换会产生积分器跳变。我在状态机切换动作里加了一步切换时把转速环积分器清零让输出从零开始按新的参考值重新建立代价是切换过程会有几百毫秒的动态调整但换来了稳定。三是如果切换过于频繁可以考虑不用硬切换而是始终用“目标功率”作为统一外环参考充电时功率为正放电时功率为负。这样状态机退化成功率调度器只需要把功率参考值映射成q轴电流给定切换点处自然连续不会出现模式切换震荡。4. Simulink仿真模型搭建步骤与关键模块配置4.1 五个子系统的搭建顺序我使用Simulink里Simscape Electrical的电气库搭建这个模型分成五个子系统依次搭建。第一步搭PMSM与逆变器。从库中拖出Permanent Magnet Synchronous Machine模块按前面表格设置电机参数输出端口接机械负载。机械输入侧直接设置转动惯量J5kg·m²、摩擦系数B取一个很小的值例如0.001N·m·s/rad这样飞轮的转动惯量特性就集成在电机模型里了。逆变器用Universal Bridge模块桥臂数量选3开关器件选IGBT/Diode直流侧接电压源或母线电容。第二步搭坐标变换。用两个Subsystem实现Clarke和Park变换输入是三相电流和转子位置角θ。电机模块通常输出转子机械角度需要乘以极对数p得到电角度。第三步搭电流环。把id的参考值设为0iq_ref由外环输出接入电流误差分别进两个PI控制器PI输出叠加前馈解耦项ud_ff、uq_ff得到dq轴的电压给定。第四步搭外环。从电机模块取转速wm和目标转速做差经转速PI后输出iq_ref同时加限幅到±40A这是整个控制里面最关键的限幅环节。第五步搭PWM与脉冲。把dq轴电压给定做Park反变换得到三相调制波与10kHz三角载波比较生成六路PWM脉冲送给Universal Bridge的门极。为节约仿真时间我第一版用SPWM正弦调制调制度超过0.9以后才考虑SVPWM。4.2 求解器设置与仿真加速技巧飞轮储能系统横跨电磁和机械两个时间尺度。电流环响应是毫秒级PWM开关是10kHz即0.1ms周期而飞轮转速爬升需要几十秒甚至上百秒如果老老实实把整个物理过程仿完仿真步长选1e-5秒仿真时长100秒那是1000万步效率低得难以接受。我用的加速技巧有三个。第一仿真测试阶段把转动惯量J从5kg·m²临时改成0.5kg·m²这样转速爬升时间缩短到原来的1/10控制策略验证通过后再改回真实值做精细仿真。第二求解器选择离散定步长步长取2e-5秒和10kHz开关周期匹配既能观察电流纹波又不浪费算力如果只看宏观能量变化可以进一步用5e-5秒配合均值模型。第三尽量用速度源或角度源代替部分反馈回路时要确认反电动势影响仍然存在不要为了加速把物理过程简化掉。4.3 监控量与SOC估算模块仿真调试过程中我最常看的波形是实际转速wm、三相电流ia/ib/ic、dq轴电流id/iq、电磁转矩Te、直流母线电压Udc。另外一定要加一路SOC估算用转速经Function模块计算(ω²-ω_min²)/(ω_max²-ω_min²)并与模式状态信号一起送入Scope实时观察。一个很实用的小技巧是把目标转速、实际转速、SOC放在同一个Scope里显示能直接判断充电/待机/放电各阶段是否有逻辑漏洞。比如充电结束SOC应该等于100%如果SOC还在缓慢下降说明待机模式没有正确补偿摩擦损耗。这种问题单看转速可能不能立刻察觉但SOC波形一对比就很明显。5. 典型工况仿真与问题排查经验5.1 充电-待机-放电循环仿真我用一组典型工况展示流程。初始飞轮转速2000r/minSOC0%充电指令下达后目标转速以斜坡从2000r/min升至6000r/min参考斜坡斜率取40r/min/s约140秒充满。实际调试时因为转动惯量缩小到0.5kg·m²充满时间约14秒便于观察。充电阶段转速环输出正q轴电流电磁转矩约10N·m三相电流幅值随转速变化基本稳定。待机阶段保持6000r/min约2秒此时q轴电流降到接近摩擦补偿的小值电磁转矩几乎为0母线电流也很小。然后下达放电指令目标转速斜坡降到3000r/min转速环输出负转矩电机进入发电状态三相电流相位与反电动势相反机械能持续转化为电能送回直流母线。如果母线侧接的是耗能负载放电期间母线电压能稳定在540V附近如果只接电容母线电压会明显升高这正好验证了功率平衡设计的正确性。5.2 问题排查速查表以下是我在这个飞轮储能仿真项目中实际遇到的问题和排查结果整理成表格供参考。故障现象最可能原因处理办法启动瞬间电流超限初始转子位置角设置不当dq电流理解错误初始化时把转子角置0并给电流环加限幅充电时转速爬升太慢转动惯量太大、转矩限幅太小或摩擦系数偏大检查J和B参数适当放宽电流限幅到40A高速区q轴电流跟踪不上dq耦合前馈缺失调制比饱和加入前馈解耦项必要时改SPWM为SVPWM放电瞬间母线电压飙升母线电容太小回馈功率大于负载吸收能力加大电容限制放电功率或接耗能负载模式切换时转速突变转速参考阶跃变化积分器初值不一致参考值加斜坡切换时转速环积分器清零三相电流波形正弦度差PWM载波频率太低或步长太大载波10kHz定步长降到2e-5秒以下SOC计算为负或超过100%转速低于ω_min或高于ω_max时公式溢出SOC估算模块加饱和限幅到0~100%5.3 几点实操心得这套模型我从零开始搭前前后后改了四版才稳定。最大的经验教训是调试顺序不能乱先确保PMSM模型单开环能转动再闭环调电流环最后才合上转速环。如果一开始就是双闭环全套电机参数和PI参数同时出错波形乱成一片根本不知道从哪下手。电机参数我建议除了Rs、Ld、Lq、ψf这四项还要把摩擦系数B记录清楚。很多仿真用默认值0.001实际电机轴承系统在高速段阻尼不是恒定的至少要设成转速的弱函数否则高速段待机功率会失真。仿真模型里可以先用常数近似工程样机测试后再修正。另外PI限幅一定要在所有PI输出端口加上。仿真新手最容易犯的错是把限幅只加在转速环输出上而忽视电流环输出的电压限幅。一旦调制度超过1逆变器输出已经失真电流环PI还在拼命积分最后必然是失控发散。电压限幅我设为±0.45倍母线电压对应SPWM的线性调制区。最后再分享一个小技巧把整台飞轮系统的参数做成一个.m脚本用Matlab的变量区统一管理。每个子系统里的PI参数、限幅值、电容值都引用工作区变量调整参数只用改脚本文件不用拆模型连线。我后来跑参数扫描时靠这个习惯至少节省了三分之一的调试时间。下一步如果要接电网侧做调频仿真这套模型只需把直流电压源替换成网侧变流器加入锁相环和功率外环就能扩展底子已经打得比较稳了。
平方复杂度无法被常规优化消除｜TaoToken 统一 Key 通道实测)
