4、飞轮储能系统建模：从数学到仿真的完整路径

各位工程师朋友，咱们今天聊聊飞轮储能系统的建模。说实话，这个环节是整个飞轮储能参与电网调频的基石。模型建不好，后面所有的控制策略、仿真验证都是空中楼阁。我个人习惯把建模分成三个层次：数学描述、控制实现、仿真验证。咱们一个一个来拆解。

4.1 飞轮储能数学模型：把物理世界翻译成方程

飞轮储能的核心，说白了就是一个旋转的飞轮把电能变成动能存起来。那怎么用数学描述这个过程呢？

最基本的动力学方程：

J * dω/dt = T_m - T_e - B * ω

其中：

• J —— 飞轮的转动惯量 (kg·m²)，这决定了飞轮能存多少能量
• ω —— 角速度 (rad/s)，飞轮转得快慢
• T_m —— 机械转矩 (N·m)，电机给飞轮的驱动力
• T_e —— 电磁转矩 (N·m)，发电时电机从飞轮取走的力
• B —— 阻尼系数，代表轴承摩擦、风阻这些损耗

嗯，这里要注意：B * ω 这一项很多人会忽略。我在项目中遇到过，刚开始建模时觉得阻尼小，直接设成0。结果仿真出来的能量衰减曲线跟实测对不上，差了将近15%。后来才意识到，高速飞轮（比如30000转/分以上）的风阻损耗其实相当可观。

能量存储方程：

E = 0.5 * J * ω²

这个公式大家应该很熟悉。飞轮储能密度跟转速的平方成正比，所以提高转速是提升能量密度的关键。但转速高了，材料强度、轴承寿命都是挑战。你想想看，一个直径1米的钢制飞轮，转速到30000转/分时，边缘线速度超过1500米/秒，比子弹还快。所以材料选择上，我建议优先考虑碳纤维复合材料，强度高、重量轻。

电机模型（永磁同步电机为例）：

d(i_d)/dt = (v_d - R_s * i_d + ω_e * L_q * i_q) / L_d
d(i_q)/dt = (v_q - R_s * i_q - ω_e * L_d * i_d - ω_e * ψ_f) / L_q
T_e = 1.5 * p * (ψ_f * i_q + (L_d - L_q) * i_d * i_q)

这里：

• i_d, i_q —— dq轴电流
• v_d, v_q —— dq轴电压
• R_s —— 定子电阻
• L_d, L_q —— dq轴电感
• ω_e —— 电角速度
• ψ_f —— 永磁体磁链
• p —— 极对数

为什么要用dq轴模型？因为交流电机在三相坐标系下是时变的，控制起来很麻烦。通过坐标变换，把交流量变成直流量，控制就简单多了。这是电机控制领域的基本功。

4.2 飞轮储能控制模型：让飞轮听话地充放电

数学模型建好了，接下来就是怎么控制它。飞轮储能参与电网调频，说白了就是电网频率高了就充电（吸收能量），频率低了就放电（释放能量）。

控制架构通常分三层：

1. 功率外环 —— 根据电网频率偏差计算目标功率
2. 转速内环 —— 根据目标功率计算目标转速
3. 电流内环 —— 根据目标转速控制电机电流

我给大家画个控制框图，这样更直观：

这个框图我用了很多年，每次给新同事讲飞轮控制，都是从这个图开始的。你看，从电网频率到目标功率，再到目标转速，最后到电流控制，一层层往下走，逻辑很清晰。

功率外环的设计要点：

功率外环的输入是频率偏差Δf，输出是目标功率P*。常用的控制策略是下垂控制：

P* = -K_droop * Δf

其中K_droop是下垂系数，单位是MW/Hz。这个系数怎么选？我建议根据飞轮储能系统的额定功率和电网允许的频率偏差范围来定。比如额定功率1MW，频率偏差允许±0.2Hz，那K_droop = 1/0.2 = 5 MW/Hz。

转速内环的设计要点：

转速内环的目标是把飞轮转速控制在合理范围内。飞轮不能一直充到最高速，也不能放到最低速，得留出裕量。我通常把工作范围设定在额定转速的70%~95%。

转速内环的PI参数整定，我推荐用工程整定法。先让系统在开环下跑，测出飞轮的机电时间常数，然后按典型I型系统来设计。具体公式：

K_p = ω_c * J / (K_t * K_fb)
K_i = K_p / T_i

其中ω_c是期望的截止频率，K_t是转矩常数，K_fb是转速反馈系数。嗯，这里要注意，飞轮的转动惯量J很大，所以转速环的响应速度不能太快，否则电流会饱和。

4.3 飞轮储能仿真模型：在电脑里跑一遍

模型建好了，控制策略设计好了，接下来就是仿真验证。我个人习惯用MATLAB/Simulink做仿真，当然也可以用其他工具，但Simulink在电力电子和电机控制领域确实方便。

仿真模型的结构：

一个完整的飞轮储能仿真模型通常包含以下几个模块：

模块名称	功能描述	关键参数
电网模型	模拟电网频率变化，提供调频需求	额定频率50Hz，频率波动范围±0.5Hz
功率外环控制器	根据频率偏差计算目标功率	下垂系数K_droop，死区阈值
转速内环控制器	根据目标功率计算目标转速	PI参数K_p, K_i，转速限幅
电流内环控制器	控制电机dq轴电流	PI参数，解耦系数
电机模型	永磁同步电机的电磁特性	L_d, L_q, R_s, ψ_f
飞轮机械模型	飞轮的转动惯量和阻尼	J, B
变流器模型	模拟IGBT开关和PWM调制	开关频率，直流母线电压

仿真参数设置示例：

% 飞轮参数
J = 100;            % 转动惯量 kg·m²
B = 0.001;          % 阻尼系数 N·m·s/rad
omega_rated = 3000; % 额定转速 rad/s (约28600转/分)

% 电机参数
R_s = 0.01;         % 定子电阻 Ω
L_d = 0.5e-3;       % d轴电感 H
L_q = 0.5e-3;       % q轴电感 H
psi_f = 0.2;        % 永磁体磁链 Wb
p = 4;              % 极对数

% 控制参数
K_droop = 5e6;      % 下垂系数 W/Hz
K_p_speed = 100;    % 转速环比例系数
K_i_speed = 10;     % 转速环积分系数
K_p_current = 50;   % 电流环比例系数
K_i_current = 500;  % 电流环积分系数

% 仿真设置
T_sim = 10;         % 仿真时长 s
T_step = 1e-5;      % 仿真步长 s (PWM开关频率10kHz)

仿真结果分析要点：

仿真跑完了，怎么看结果？我一般关注以下几个指标：

1. 频率响应时间 —— 从频率偏差发生到飞轮输出功率达到目标值的时间，一般要求小于100ms
2. 功率跟踪精度 —— 实际输出功率与目标功率的偏差，稳态误差应小于5%
3. 转速波动范围 —— 飞轮转速是否在安全范围内，不能超速也不能过低
4. 电流谐波含量 —— 电机电流的THD，一般要求小于5%

我记得有一次做仿真，频率响应时间总是达不到要求，调了半天PI参数都没用。后来仔细一看，是电流内环的带宽设得太低了。电流环响应慢，转速环和功率环再快也没用。所以调试的时候，一定要从内环往外环调，先把电流环调好，再调转速环，最后调功率环。

仿真与实测的差异：

仿真毕竟是仿真，跟实际系统肯定有差距。我总结了几点常见的差异来源：

• 损耗模型不准确 —— 仿真里的阻尼系数B是常数，实际中轴承摩擦、风阻都跟转速、温度有关
• 死区时间影响 —— 变流器的死区时间在仿真里往往被忽略，但实际中会导致电流畸变
• 采样延迟 —— 仿真里信号是连续的，实际中AD采样、PWM更新都有延迟
• 温度影响 —— 电机电阻、永磁体磁链都随温度变化，仿真里通常是常数

所以我的建议是：仿真结果只能作为参考，不能完全相信。关键参数一定要留裕量，比如PI参数的稳定裕度至少留30%。

好了，关于飞轮储能系统建模，我就讲这么多。各位在实际项目中遇到什么问题，欢迎随时交流。记住，模型是工具，不是目的。我们的最终目标是让飞轮储能系统在电网调频中发挥实实在在的作用。

---

公众号:蓝海资料掘金营，微信deep3321