4、DFIG数学模型(下):坐标变换理论(Clark变换、Park变换)、两相旋转坐标系(dq)下的数学模型

好,咱们接着聊。上一节我们把DFIG在三相静止坐标系下的数学模型给捋了一遍。说实话,那个模型看着就头疼,对不对?一堆互感、时变系数,耦合得跟麻花似的。你要是拿那个模型去设计控制器,我估计你调参调到怀疑人生。

所以这一节,我们要干一件大事——坐标变换。说白了,就是换个角度看问题。把交流量变成直流量,把时变系统变成定常系统。嗯,这才是矢量控制的精髓所在。

4.1 为什么要做坐标变换?

我个人习惯,在讲技术之前,先想清楚“为什么”。

你看DFIG的定转子绕组,它们在空间上差一个角度,转子还在转。这就导致互感系数是转子位置角θr的函数。你想想看,一个控制系统里,参数随着时间变来变去,这控制器得多难设计?

坐标变换的核心思想就一句话:把三相静止坐标系下的交流量,等效变换到两相旋转坐标系下的直流量

我在项目中遇到过这种情况:刚开始用三相模型做仿真,电流波形看着还行,但一上硬件就崩。后来发现,就是因为模型太复杂,控制器根本跟不上那个时变的速度。换成dq模型之后,问题迎刃而解。

4.2 Clark变换:从三相到两相静止

Clark变换,也叫3s/2s变换。就是把三相静止坐标系(abc)下的量,变换到两相静止坐标系(αβ)下。

为什么要先变到αβ?因为三相系统其实有冗余。三相对称系统里,ia + ib + ic = 0,所以三个变量里只有两个是独立的。我们把它压缩成两个正交的变量,方便后续处理。

变换公式长这样(等幅值变换):

[ i_α ]   [ 1    -1/2    -1/2 ] [ i_a ]
[ i_β ] = [ 0    √3/2   -√3/2 ] [ i_b ]
                                [ i_c ]

注意,这里有个系数2/3。为什么?因为我们要保证变换前后幅值不变。如果你做等功率变换,系数就是√(2/3)。我个人习惯用等幅值变换,因为调试的时候看波形更直观。

我的小技巧: 实际编程的时候,别每次都算三角函数。把变换矩阵存成常量,直接做矩阵乘法,效率高很多。

4.3 Park变换:从两相静止到两相旋转

Clark变换做完,我们得到了αβ坐标系下的量。但αβ坐标系还是静止的,里面的电流电压依然是交流量。怎么办?再转一次。

Park变换,也叫2s/2r变换。就是把两相静止坐标系(αβ)下的量,变换到两相旋转坐标系(dq)下。

这个变换的核心,就是让坐标系跟着转子磁场一起转。这样,原本旋转的交流量,在旋转坐标系下看就是静止的直流量了。

公式如下:

[ i_d ]   [ cosθ   sinθ ] [ i_α ]
[ i_q ] = [ -sinθ  cosθ ] [ i_β ]

这里的θ,是d轴与α轴之间的夹角。在DFIG矢量控制中,这个θ通常取定子磁链的位置角。

我曾经犯过一个低级错误:把θ算反了。结果电机转起来之后,转矩方向完全反了,差点把联轴器扭断。嗯,从那以后,我每次做Park变换都要先检查θ的符号。

4.4 两相旋转坐标系(dq)下的DFIG数学模型

好了,重头戏来了。经过Clark和Park变换,我们终于可以在dq坐标系下写出DFIG的数学模型。这个模型,才是矢量控制真正用的模型。

先看电压方程:

定子电压:
u_sd = R_s * i_sd + d(ψ_sd)/dt - ω_s * ψ_sq
u_sq = R_s * i_sq + d(ψ_sq)/dt + ω_s * ψ_sd

转子电压:
u_rd = R_r * i_rd + d(ψ_rd)/dt - (ω_s - ω_r) * ψ_rq
u_rq = R_r * i_rq + d(ψ_rq)/dt + (ω_s - ω_r) * ψ_rd

你发现没有?这里多出了两项:-ωψ和+ωψ。这就是旋转坐标系带来的“反电动势”项,也叫交叉耦合项。说白了,d轴和q轴之间不是独立的,它们会互相影响。

再看磁链方程:

ψ_sd = L_s * i_sd + L_m * i_rd
ψ_sq = L_s * i_sq + L_m * i_rq

ψ_rd = L_r * i_rd + L_m * i_sd
ψ_rq = L_r * i_rq + L_m * i_sq

这里Ls是定子自感,Lr是转子自感,Lm是互感。注意,在dq坐标系下,这些电感都是常数!不再随转子位置变化了。这就是坐标变换最大的好处。

最后,转矩方程:

T_e = 1.5 * n_p * (ψ_sd * i_sq - ψ_sq * i_sd)

如果采用定子磁链定向(即ψsq = 0),转矩公式可以简化为:

T_e = 1.5 * n_p * ψ_sd * i_sq

你看,转矩和q轴电流成正比!这就是矢量控制“解耦”的本质——通过坐标变换,把转矩和磁链的控制分开了。

核心结论: 在dq坐标系下,DFIG的数学模型从时变系统变成了定常系统。转矩由q轴电流控制,磁链由d轴电流控制。这就是双馈风机矢量控制的理论基础。

4.5 知识体系总览

为了让你更直观地理解整个坐标变换的脉络,我画了一张图:

DFIG坐标变换知识体系 三相静止坐标系 (abc) 时变、强耦合 Clark变换 两相静止坐标系 (αβ) 正交、仍交流 Park变换 两相旋转坐标系 (dq) 直流、解耦 dq坐标系下DFIG数学模型 电压方程:u = Ri + dψ/dt + ωψ(交叉耦合项) 磁链方程:ψ = Li(电感为常数) 转矩方程:Te ∝ ψsd · isq(解耦控制)

4.6 避坑指南

坐标变换看着简单,但实际做起来坑不少。我把自己踩过的坑分享给你:

  • 角度对齐问题: 我曾经在项目里把转子位置角θr和定子磁链角θs搞混了。结果变换出来的dq分量完全不对,电机根本转不起来。记住:Park变换用的θ,是d轴与α轴的夹角,不是转子位置角。
  • 变换系数选择: 等幅值变换和等功率变换,系数不一样。如果你从论文里抄公式,一定要看清楚人家用的是哪种。我建议你从头到尾统一用一种,别混着用。
  • 初始相位: 仿真刚开始的时候,所有状态变量都是0。但Park变换里的θ怎么取?我习惯让θ从0开始,然后根据磁链观测器实时更新。
重要提醒: 坐标变换本身不改变系统的物理本质。它只是换了一个数学描述方式。所以,变换前后的功率、转矩等物理量必须守恒。如果你发现变换前后功率对不上,赶紧检查变换矩阵的系数。

4.7 小结

这一节我们干了三件事:

  1. 用Clark变换把三相变两相静止
  2. 用Park变换把两相静止变两相旋转
  3. 在dq坐标系下写出了DFIG的数学模型

现在你手里有了一个定常、解耦的数学模型。下一节,我们就可以基于这个模型,开始设计矢量控制器了。嗯,那才是真正有意思的部分。

我的建议: 如果你刚开始学,别急着写代码。先在纸上把Clark和Park变换的手算一遍,用一组简单的三相电流(比如三相对称正弦波)算一遍,看看变换后的dq分量是不是直流量。这一步走通了,后面就顺了。

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