3. 转子侧变流器数学模型
做双馈变流器控制,绕不开数学模型。
说实话,我刚入行那会儿,看到一堆矩阵和坐标变换就头疼。后来在项目里吃过亏——有一次调试,转子电流波形怎么都调不好,最后发现是模型里漏了一项交叉耦合项。从那以后,我老老实实把数学模型从头捋了一遍。
这一章,我们就来聊聊DFIG的数学模型。从最直观的三相静止坐标系开始,再到工程上最常用的两相旋转坐标系。你想想看,搞懂了这些,后面的控制策略才能站得住脚。
3.1 DFIG在三相静止坐标系下的数学模型
三相静止坐标系,也叫ABC坐标系。说白了,就是直接拿A、B、C三相的物理量来描述电机。
DFIG在这个坐标系下的模型,由电压方程、磁链方程和转矩方程三部分组成。我习惯先把磁链方程写出来,因为电压方程本质上就是磁链对时间的导数加上电阻压降。
3.1.1 电压方程
定子和转子的电压方程,写成矩阵形式是这样的:
┌ u_sa ┐ ┌ R_s 0 0 0 0 0 ┐ ┌ i_sa ┐ ┌ ψ_sa ┐
│ u_sb │ │ 0 R_s 0 0 0 0 │ │ i_sb │ d │ ψ_sb │
│ u_sc │ = │ 0 0 R_s 0 0 0 │ │ i_sc │ + ── │ ψ_sc │
│ u_ra │ │ 0 0 0 R_r 0 0 │ │ i_ra │ dt │ ψ_ra │
│ u_rb │ │ 0 0 0 0 R_r 0 │ │ i_rb │ │ ψ_rb │
└ u_rc ┘ └ 0 0 0 0 0 R_r ┘ └ i_rc ┘ └ ψ_rc ┘
嗯,这里要注意:R_s是定子电阻,R_r是转子电阻。电阻矩阵是对角阵,没什么好说的。关键是磁链那一项,它包含了自感和互感,耦合关系比较复杂。
3.1.2 磁链方程
磁链方程才是重头戏。定子和转子之间的互感,是转子位置角θ_r的函数。这就意味着,随着转子旋转,互感一直在变。
┌ ψ_sabc ┐ ┌ L_s M_sr(θ_r) ┐ ┌ i_sabc ┐
│ │ = │ │ │ │
└ ψ_rabc ┘ └ M_rs(θ_r) L_r ┘ └ i_rabc ┘
其中:
L_s:定子自感矩阵(常数)L_r:转子自感矩阵(常数)M_sr(θ_r):定转子互感矩阵(随θ_r变化)
关键点:互感矩阵M_sr(θ_r)包含cos(θ_r)、cos(θ_r+120°)、cos(θ_r-120°)等项。这就是为什么三相模型这么复杂——时变系数让分析和控制变得非常麻烦。
我在项目里遇到过一个问题:直接用三相模型做仿真,步长稍微大一点,结果就发散。后来才意识到,时变互感导致系统刚性很强,必须用小步长或者改用dq模型。
3.1.3 转矩方程
电磁转矩的表达式:
T_e = -n_p · L_m · [ (i_sa·i_ra + i_sb·i_rb + i_sc·i_rc) · sin(θ_r) + ... ]
(完整表达式很长,这里只给核心形式)
说白了,转矩也是转子位置角的函数。这给转矩控制带来了麻烦——你没法直接解耦控制。
3.2 两相旋转坐标系(dq)下的数学模型
为什么要转到dq坐标系?
原因很简单:把时变系数变成常数。你想想看,一个时变系统和一个定常系统,哪个好控制?当然是后者。
dq变换的核心思想:把三相静止的物理量,投影到两个随转子旋转的轴上。d轴和q轴互相垂直,而且以同步转速旋转。
我的经验:dq变换的物理意义,可以理解为「把交流量变成直流量」。这样一来,PI控制器就能实现无静差跟踪。我在调试网侧变流器时,深刻体会到了这一点——没有dq变换,电流环根本稳不住。
3.2.1 变换矩阵
从ABC到dq的变换矩阵(等幅值变换):
┌ x_d ┐ 2 ┌ cos(θ) cos(θ-120°) cos(θ+120°) ┐ ┌ x_a ┐
│ │ = ── │ │ │ │
│ x_q │ 3 └ -sin(θ) -sin(θ-120°) -sin(θ+120°) ┘ │ x_b │
└ x_c ┘
其中θ是d轴与A相轴线的夹角。对于转子侧,θ = θ_s - θ_r(转差角)。
3.2.2 dq坐标系下的电压方程
经过变换后,电压方程变得清爽多了:
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 - ω_slip · ψ_rq
u_rq = R_r · i_rq + d(ψ_rq)/dt + ω_slip · ψ_rd
这里:
ω_s:同步角速度ω_slip = ω_s - ω_r:转差角速度
注意:方程中出现了-ω_s·ψ_sq和+ω_s·ψ_sd这样的交叉耦合项。这就是dq模型的特点——d轴和q轴不是完全独立的,它们之间存在耦合。控制时需要进行解耦补偿。
我记得第一次做转子侧控制时,没加解耦补偿,结果d轴电流一变,q轴电流也跟着抖。后来加上了前馈解耦项,波形才漂亮起来。
3.2.3 dq坐标系下的磁链方程
磁链方程在dq坐标系下变成了常数矩阵:
ψ_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
其中:
L_s = L_ls + L_m:定子自感L_r = L_lr + L_m:转子自感L_m:励磁互感
你看,所有电感都变成了常数。这就是dq变换的魔力——把时变系统变成了定常系统。
3.2.4 dq坐标系下的转矩方程
转矩方程也简化了:
T_e = 1.5 · n_p · L_m · (i_sq · i_rd - i_sd · i_rq)
或者用磁链表示:
T_e = 1.5 · n_p · (ψ_sd · i_sq - ψ_sq · i_sd)
这个形式简单多了。而且你会发现,转矩只和d、q轴分量有关,和转子位置角无关了。这就是为什么矢量控制能实现转矩和磁链的解耦控制。
3.3 磁链与电压方程的关系
把磁链方程代入电压方程,可以得到更实用的形式。以转子侧为例:
u_rd = R_r · i_rd + L_r · di_rd/dt + L_m · di_sd/dt - ω_slip · (L_r · i_rq + L_m · i_sq)
u_rq = R_r · i_rq + L_r · di_rq/dt + L_m · di_sq/dt + ω_slip · (L_r · i_rd + L_m · i_sd)
这个形式在控制器设计时非常有用。你可以看到:
- 第一项是电阻压降
- 第二、三项是电感上的暂态压降
- 第四项是旋转电动势(反电势)
避坑指南:我曾经在调试时忽略了反电势项,结果电流环带宽一高就振荡。后来才明白,反电势是前馈补偿的关键。不加前馈,控制器相当于在跟一个不断变化的扰动做斗争,很难调好。
3.4 知识体系总览
下面这张图,是我梳理的本章知识结构。你可以把它当作一个导航图:
从这张图可以看得很清楚:三相模型是基础,但工程实现必须转到dq坐标系。中间的桥梁就是坐标变换。
3.5 小结
这一章我们聊了DFIG的数学模型,从ABC到dq,从时变到定常。说白了,数学模型的建立,就是为了给控制器设计铺路。
我个人觉得,理解磁链方程比电压方程更重要。因为磁链是状态变量,电压只是输入。你抓住了磁链的变化规律,就抓住了电机的本质。
下一章,我们会基于这些数学模型,推导转子侧变流器的控制策略。到时候你会发现,今天打下的基础,全都能用上。
公众号:蓝海资料掘金营,微信deep3321