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 知识体系总览

下面这张图,是我梳理的本章知识结构。你可以把它当作一个导航图:

DFIG转子侧数学模型知识体系 DFIG数学模型 三相静止坐标系(ABC) 两相旋转坐标系(dq) 电压方程 磁链方程 转矩方程 电压方程 磁链方程 转矩方程 核心:ABC时变系统 → dq定常系统 变换工具:Park变换 / Clark变换 关键参数:同步角速度ωs、转差角速度ωslip、转子位置角θr 工程应用:矢量控制、直接转矩控制、模型预测控制

从这张图可以看得很清楚:三相模型是基础,但工程实现必须转到dq坐标系。中间的桥梁就是坐标变换。

3.5 小结

这一章我们聊了DFIG的数学模型,从ABC到dq,从时变到定常。说白了,数学模型的建立,就是为了给控制器设计铺路。

我个人觉得,理解磁链方程比电压方程更重要。因为磁链是状态变量,电压只是输入。你抓住了磁链的变化规律,就抓住了电机的本质。

下一章,我们会基于这些数学模型,推导转子侧变流器的控制策略。到时候你会发现,今天打下的基础,全都能用上。


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