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轴电流成正比!这就是矢量控制“解耦”的本质——通过坐标变换,把转矩和磁链的控制分开了。
4.5 知识体系总览
为了让你更直观地理解整个坐标变换的脉络,我画了一张图:
4.6 避坑指南
坐标变换看着简单,但实际做起来坑不少。我把自己踩过的坑分享给你:
- 角度对齐问题: 我曾经在项目里把转子位置角θr和定子磁链角θs搞混了。结果变换出来的dq分量完全不对,电机根本转不起来。记住:Park变换用的θ,是d轴与α轴的夹角,不是转子位置角。
- 变换系数选择: 等幅值变换和等功率变换,系数不一样。如果你从论文里抄公式,一定要看清楚人家用的是哪种。我建议你从头到尾统一用一种,别混着用。
- 初始相位: 仿真刚开始的时候,所有状态变量都是0。但Park变换里的θ怎么取?我习惯让θ从0开始,然后根据磁链观测器实时更新。
4.7 小结
这一节我们干了三件事:
- 用Clark变换把三相变两相静止
- 用Park变换把两相静止变两相旋转
- 在dq坐标系下写出了DFIG的数学模型
现在你手里有了一个定常、解耦的数学模型。下一节,我们就可以基于这个模型,开始设计矢量控制器了。嗯,那才是真正有意思的部分。
公众号:蓝海资料掘金营,微信deep3321