2. 双馈风机数学模型:从三相静止到两相旋转

做双馈风机控制,第一步就是搞懂它的数学模型。说实话,我刚入行那会儿,看着一堆微分方程和坐标变换矩阵,头都大了。后来在项目里调试过几次谐波问题,才真正体会到——模型建得准不准,直接决定了控制策略好不好使。

这一节,咱们就聊聊DFIG的数学模型。我会从最直观的三相静止坐标系开始,然后带你做坐标变换,最后落到两相旋转坐标系上。嗯,这是后面所有控制算法的基础。

2.1 三相静止坐标系下的DFIG模型

双馈风机本质上是个绕线式异步电机。定子直接挂电网,转子通过背靠背变流器接电网。我习惯把它的模型分成两部分看:电压方程和磁链方程。

2.1.1 电压方程

先看定子侧。三相定子绕组的电压方程长这样:

u_sa = R_s * i_sa + d(ψ_sa)/dt
u_sb = R_s * i_sb + d(ψ_sb)/dt
u_sc = R_s * i_sc + d(ψ_sc)/dt

转子侧也类似,只不过转子是旋转的:

u_ra = R_r * i_ra + d(ψ_ra)/dt
u_rb = R_r * i_rb + d(ψ_rb)/dt
u_rc = R_r * i_rc + d(ψ_rc)/dt

这里R_s和R_r分别是定转子电阻。ψ代表磁链。你可能会问,为什么转子方程里没有反电动势项?其实它藏在磁链对时间的导数里了。

我的经验: 实际项目中,定子电阻R_s通常很小,但千万别忽略它。我曾经在弱电网工况下吃过亏,R_s上的压降对控制精度影响挺大的。

2.1.2 磁链方程

磁链方程就复杂一些了。定转子之间有互感,而且互感会随着转子位置变化。写成矩阵形式:

[ψ_s]   [L_s  M_sr] [i_s]
[ψ_r] = [M_rs  L_r] [i_r]

其中L_s是定子自感,L_r是转子自感,M_sr是定转子互感矩阵。这个互感矩阵里全是转子位置角θ_r的正余弦函数。说白了,就是转子转到哪里,互感就变到哪里。

注意: 三相静止坐标系下的模型,互感矩阵是时变的。这意味着你直接用它做控制,每个控制周期都要重新计算一大堆三角函数。计算量很大,实时性很难保证。

2.2 坐标变换:Clark和Park

为了解决时变问题,我们需要坐标变换。说白了,就是把三相交流量变成两相直流量。这样控制起来就简单多了。

2.2.1 Clark变换(3s/2s)

Clark变换把三相静止坐标系(abc)变到两相静止坐标系(αβ)。我习惯用等幅值变换:

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

为什么要做这个变换?因为三相系统里有很多冗余信息。你想想看,三相对称系统里,i_a + i_b + i_c = 0,实际上只有两个自由度。Clark变换就是把这个冗余去掉。

关键点: Clark变换后,αβ坐标系下的模型仍然是交流量。对于工频50Hz的电网,αβ轴上的电压电流还是50Hz的正弦波。要变成直流量,还得再做一次Park变换。

2.2.2 Park变换(2s/2r)

Park变换把两相静止坐标系变到两相旋转坐标系(dq)。旋转角速度等于电网角频率ω_e:

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

θ_e是电网电压的相位角。这个角度通常通过锁相环(PLL)获得。我在项目里调试PLL时踩过不少坑,后面会专门讲。

做完Park变换,基波分量就变成了直流量。d轴和q轴分别对应有功和无功的控制通道。这才是我们想要的。

2.3 两相旋转坐标系下的DFIG模型

好了,现在我们把Clark和Park变换应用到DFIG模型上。结果会清爽很多。

2.3.1 定子电压方程(dq坐标系)

u_sd = R_s * i_sd + d(ψ_sd)/dt - ω_e * ψ_sq
u_sq = R_s * i_sq + d(ψ_sq)/dt + ω_e * ψ_sd

看到没?多了两项-ω_e * ψ_sq和+ω_e * ψ_sd。这就是旋转坐标系带来的交叉耦合项。我刚开始学的时候,总觉得这两项是多余的。直到有一次做电流环参数整定,发现d轴和q轴互相影响,才明白耦合项必须解掉。

2.3.2 转子电压方程(dq坐标系)

u_rd = R_r * i_rd + d(ψ_rd)/dt - (ω_e - ω_r) * ψ_rq
u_rq = R_r * i_rq + d(ψ_rq)/dt + (ω_e - ω_r) * ψ_rd

这里ω_r是转子转速。注意看,转子方程里的耦合项是(ω_e - ω_r),也就是转差频率。这个物理意义很明确:转子绕组切割磁场的速度就是转差速度。

避坑指南: 我曾经在双馈风机并网调试时,发现转子电流一直有低频振荡。查了两天才发现,是转子侧PLL的带宽设置太高,把转差频率的波动引入了电流环。后来把PLL带宽降到10Hz以下,问题就解决了。

2.3.3 磁链方程(dq坐标系)

在dq坐标系下,磁链方程变得非常简单:

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

L_m是激磁电感。所有互感都变成了常数,不再随转子位置变化。这就是坐标变换最大的好处——模型变成了线性时不变系统。

2.4 本章知识体系

为了让你更直观地理解整个模型推导的脉络,我画了一张图:

DFIG数学模型推导脉络 三相静止坐标系 电压方程 磁链方程(时变互感) 模型复杂,计算量大 Clark变换 两相静止坐标系 αβ轴分量 去除冗余信息 仍是交流量 Park变换 两相旋转坐标系 dq轴分量 直流量 线性时不变系统 应用:矢量控制策略 d轴控制有功功率 / 转子励磁 q轴控制无功功率 / 电磁转矩 从三相静止到两相旋转,模型从时变系统变为时不变系统 为后续的PI控制器设计和谐波抑制策略奠定基础

2.5 模型参数表

为了方便你查用,我把DFIG模型的主要参数整理成了一张表:

符号 含义 典型值范围 备注
R_s 定子电阻 0.01~0.05 pu 大功率机组偏小
R_r 转子电阻 0.01~0.05 pu 与定子相当
L_s 定子自感 2.0~4.0 pu 含漏感和激磁电感
L_r 转子自感 2.0~4.0 pu 含漏感和激磁电感
L_m 激磁电感 1.8~3.8 pu 决定无功能力
ω_e 电网角频率 314 rad/s (50Hz) 由PLL锁定
ω_r 转子电角速度 0.7~1.3 pu 变速运行范围
总结一下: 三相静止坐标系下的DFIG模型,互感矩阵是时变的,不适合直接做控制。通过Clark和Park变换,我们得到了dq坐标系下的模型,所有参数都变成了常数。这个模型是矢量控制的基础。后面讲谐波抑制时,我们还会在这个模型上做文章。

好了,这一节的内容就到这儿。模型是基础,但光有模型还不够。下一节我们会聊怎么用这个模型做控制——也就是矢量控制策略。到时候我会结合项目里的实际案例,讲讲参数怎么整定、耦合怎么解。

个人建议: 如果你刚开始接触DFIG,建议在仿真软件里搭一下这个模型。把Clark和Park变换的每一步都算一遍,看看波形从交流变直流的过程。我当年就是这么学的,效果比死记公式好得多。

专注资料整理