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分别是定转子电阻。ψ代表磁链。你可能会问,为什么转子方程里没有反电动势项?其实它藏在磁链对时间的导数里了。
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变换就是把这个冗余去掉。
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),也就是转差频率。这个物理意义很明确:转子绕组切割磁场的速度就是转差速度。
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 本章知识体系
为了让你更直观地理解整个模型推导的脉络,我画了一张图:
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 | 变速运行范围 |
好了,这一节的内容就到这儿。模型是基础,但光有模型还不够。下一节我们会聊怎么用这个模型做控制——也就是矢量控制策略。到时候我会结合项目里的实际案例,讲讲参数怎么整定、耦合怎么解。