3、DFIG数学模型(上):三相静止坐标系下的电压方程、磁链方程、转矩方程

各位工程师朋友,咱们今天来啃一块硬骨头——双馈电机的数学模型。

说实话,我刚入行那会儿,看到这些方程就头疼。密密麻麻的变量,绕来绕去的耦合关系,简直像一团乱麻。但后来我发现,只要把坐标系理清楚,这团乱麻其实有迹可循。

今天我们先聊三相静止坐标系下的模型。这是最原始、最直观的模型,也是后面所有坐标变换的基础。

3.1 为什么非要建这个模型?

你想想看,我们要控制一个电机,总得知道它内部发生了什么吧?电压加进去,电流怎么变?磁链怎么建立?转矩怎么产生?

数学模型就是干这个的。它用数学语言,把电机的物理规律描述出来。

我个人习惯,在搭建模型前,先做几个假设,把问题简化一下:

  • 忽略磁路饱和:说白了,就是认为磁路是线性的,电感是常数。实际中当然有饱和,但建模时先不考虑,否则方程就太复杂了。
  • 忽略铁耗:涡流和磁滞损耗,我们先放一边。
  • 三相绕组对称:A、B、C三相,参数完全一样,空间上相差120度。
  • 气隙磁场正弦分布:只考虑基波,忽略高次谐波。

嗯,这些假设在工程上是可接受的。我在做兆瓦级风机仿真时,用的就是这套简化模型,结果和实测数据吻合得还不错。

3.2 电压方程——电机的“欧姆定律”

电压方程,说白了就是基尔霍夫电压定律在电机上的应用。绕组的端电压,等于电阻压降加上感应电动势。

双馈电机有两套绕组:定子绕组转子绕组。我们先看定子侧:

u_A = R_s * i_A + dψ_A / dt
u_B = R_s * i_B + dψ_B / dt
u_C = R_s * i_C + dψ_C / dt

这里:

  • u_A, u_B, u_C 是定子三相电压
  • i_A, i_B, i_C 是定子三相电流
  • ψ_A, ψ_B, ψ_C 是定子三相磁链
  • R_s 是定子每相电阻

转子侧也类似,但要注意:转子绕组是短路的(通过变流器接电网),而且转子在旋转,所以方程里多了一个旋转项。

u_a = R_r * i_a + dψ_a / dt
u_b = R_r * i_b + dψ_b / dt
u_c = R_r * i_c + dψ_c / dt

其中 R_r 是转子每相电阻。u_a, u_b, u_c 是转子三相电压,注意这里是小写字母,表示转子侧的量。

⚠️ 避坑指南: 我曾经在仿真时,把定子和转子的电阻搞反了。结果算出来的电流波形怎么看怎么不对。后来查了半天,才发现是电阻值标错了。各位切记,定转子电阻不是一个数量级,定子电阻通常比转子电阻大不少。

3.3 磁链方程——耦合关系的核心

磁链方程是DFIG模型里最绕的部分。为什么?因为定子和转子之间有互感,而且这个互感还随着转子位置变化。

先看定子磁链:

ψ_A = L_s * i_A + L_m * [i_a * cos(θ_r) + i_b * cos(θ_r + 120°) + i_c * cos(θ_r - 120°)]
ψ_B = L_s * i_B + L_m * [i_a * cos(θ_r - 120°) + i_b * cos(θ_r) + i_c * cos(θ_r + 120°)]
ψ_C = L_s * i_C + L_m * [i_a * cos(θ_r + 120°) + i_b * cos(θ_r - 120°) + i_c * cos(θ_r)]

这里:

  • L_s 是定子每相自感(包括漏感和主电感)
  • L_m 是定转子之间的互感幅值
  • θ_r 是转子位置角(电角度)

转子磁链也类似,只是把角色互换:

ψ_a = L_r * i_a + L_m * [i_A * cos(θ_r) + i_B * cos(θ_r - 120°) + i_C * cos(θ_r + 120°)]
ψ_b = L_r * i_b + L_m * [i_A * cos(θ_r + 120°) + i_B * cos(θ_r) + i_C * cos(θ_r - 120°)]
ψ_c = L_r * i_c + L_m * [i_A * cos(θ_r - 120°) + i_B * cos(θ_r + 120°) + i_C * cos(θ_r)]

你发现没有?这些方程里全是cos函数。这意味着磁链是时变的,而且和转子位置强相关。这就是为什么三相静止坐标系下的模型这么复杂——耦合太严重了

💡 核心要点: 三相静止坐标系下,定转子之间的互感是转子位置角θ_r的函数。这个时变耦合,是后面所有坐标变换的出发点。

3.4 转矩方程——能量转换的桥梁

转矩方程,说白了就是电机怎么把电能变成机械能。对于DFIG,电磁转矩可以表示为:

T_e = -n_p * L_m * [ (i_A * i_a + i_B * i_b + i_C * i_c) * sin(θ_r) 
                     + (i_A * i_b + i_B * i_c + i_C * i_a) * sin(θ_r + 120°)
                     + (i_A * i_c + i_B * i_a + i_C * i_b) * sin(θ_r - 120°) ]

其中 n_p 是极对数。

这个方程看着就让人头大,对吧?但它的物理意义其实很清晰:转矩是定转子电流相互作用的结果。而且,转矩也和转子位置有关。

我记得有一次调试,发现转矩波形有奇怪的振荡。后来一查,是电流采样有延迟,导致转矩计算里的角度不对。嗯,这种问题在工程现场很常见。

3.5 知识体系总览

为了帮你理清思路,我画了一张图,把今天讲的内容串起来:

DFIG三相静止坐标系数学模型 电压方程 u = R·i + dψ/dt 定子:u_A, u_B, u_C 转子:u_a, u_b, u_c 磁链方程 ψ = L·i + L_m·i·cos(θ_r) 自感 + 互感(时变) 耦合核心 转矩方程 T_e = f(i_s, i_r, θ_r) 能量转换桥梁 机电耦合 核心特征 时变耦合 | 非线性 | 强耦合 | 三相不对称(转子旋转) 下一站:坐标变换 → 简化模型

3.6 小结与个人体会

好了,今天的内容就到这里。总结一下:

  • 电压方程:就是绕组的KVL,定转子各一套。
  • 磁链方程:自感+互感,互感随转子位置变化,是耦合的根源。
  • 转矩方程:定转子电流相互作用的结果,也是时变的。

说实话,三相静止坐标系下的模型,在实际控制中几乎不用。为什么?因为太复杂了,全是时变系数,没法设计控制器。但它是基础,是理解后面所有变换的起点。

我个人建议,你可以在仿真软件里搭一下这个模型,看看磁链和转矩的波形。亲眼看到那些时变耦合的效果,比看公式印象深得多。

🔧 实用技巧: 如果你用Matlab/Simulink搭建这个模型,建议用S-function或者自定义函数模块来实现磁链方程。直接用基本模块搭,连线会多到你怀疑人生。

下一节,我们会引入坐标变换,把这个时变耦合模型变成常系数模型。那才是真正能用的控制模型。


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