2、电机数学模型:永磁同步电机(PMSM)的数学模型、Clark变换与Park变换、电压方程与磁链方程

好,咱们直接进入正题。做无传感器速度估算,说白了,你得先懂电机本身是怎么“转”起来的。数学模型就是电机的“身份证”,你搞不清楚它,后面所有算法都是空中楼阁。

我个人习惯,拿到一个新项目,第一件事不是调代码,而是先把电机的数学模型在白板上推一遍。为什么?因为后面所有观测器、估算器,都是从这几个方程里“长”出来的。

2.1 永磁同步电机(PMSM)的物理模型

先看个最简单的图。PMSM的定子上有三相绕组,转子上有永磁体。你给三相绕组通上电,就会产生旋转磁场,这个磁场“拽”着转子上的永磁体一起转。

核心要点: 电机的电磁转矩,本质上是定子磁场和转子磁场“较劲”的结果。你控制电流,就是在控制这个“较劲”的力度和角度。

嗯,这里要注意:PMSM的转子是永磁体,所以它没有励磁绕组,也就没有励磁电流。这一点和异步电机完全不同。我在项目中遇到过有人把异步电机的控制思路直接套到PMSM上,结果转矩控制一塌糊涂。

PMSM 数学模型与坐标变换关系图 三相静止坐标系 A-B-C 坐标系 电压方程: v_a, v_b, v_c i_a, i_b, i_c Clark 3s → 2s 两相静止坐标系 α-β 坐标系 电压方程: v_α, v_β i_α, i_β Park 2s → 2r 两相旋转坐标系 d-q 坐标系 电压方程: v_d, v_q i_d, i_q 核心数学模型(d-q 坐标系下) 电压方程: v_d = R_s · i_d + L_d · (di_d/dt) - ω_e · L_q · i_q v_q = R_s · i_q + L_q · (di_q/dt) + ω_e · (L_d · i_d + ψ_f) 磁链方程: ψ_d = L_d · i_d + ψ_f ψ_q = L_q · i_q

2.2 Clark变换:从ABC到αβ

三相系统分析起来太麻烦。你想想看,三个变量互相耦合,解耦都费劲。所以我们要做坐标变换,把三个变量变成两个。

Clark变换,就是把三相静止的ABC坐标系,映射到两相静止的αβ坐标系。说白了,就是“降维”。

数学上很简单,就是下面这个矩阵乘法:

// Clark变换 (等幅值变换)
i_α = i_a
i_β = (i_a + 2*i_b) / √3

// 或者写成矩阵形式
[i_α]   [ 1      -1/2    -1/2   ] [i_a]
[i_β] = [ 0      √3/2   -√3/2  ] [i_b]
                                [i_c]

我的经验: 实际代码里,我一般用等幅值变换,因为这样变换前后电流幅值不变,调试时看波形更直观。等功率变换虽然数学上更“干净”,但幅值会变,容易把人搞晕。

Clark变换之后,我们得到了α轴和β轴上的电流。这两个轴在空间上是垂直的,互相独立。嗯,到这里,问题简化了一半。

2.3 Park变换:从αβ到dq

Clark变换只是第一步。αβ坐标系还是静止的,而电机转子在转。你想想看,如果坐标系跟着转子一起转,那会怎样?

Park变换就是干这个的。它把静止的αβ坐标系,旋转到与转子同步的dq坐标系。d轴对准转子磁极方向,q轴超前90度。

为什么要这么干?因为变换之后,所有交流量都变成了直流量!

// Park变换
i_d =  i_α * cos(θ_e) + i_β * sin(θ_e)
i_q = -i_α * sin(θ_e) + i_β * cos(θ_e)

// 其中 θ_e 是转子电角度

关键理解: Park变换的本质是“解耦”。变换之后,d轴电流控制磁通,q轴电流控制转矩。这两个通道可以独立设计PI控制器,这就是矢量控制(FOC)的数学基础。

我曾经在调试一个高速电机项目时,发现q轴电流响应特别慢。查了半天,发现是Park变换的角度θ_e有延迟。你想想看,角度不准,变换就不准,电流环能好才怪。

2.4 电压方程与磁链方程

好,现在坐标系已经转到dq轴上了。我们来看看电机在这个坐标系下长什么样。

电压方程:

v_d = R_s * i_d + L_d * (di_d/dt) - ω_e * L_q * i_q
v_q = R_s * i_q + L_q * (di_q/dt) + ω_e * (L_d * i_d + ψ_f)

这里面每个项都有物理意义:

  • R_s * i_d/q:电阻压降,这个好理解
  • L_d/q * (di/dt):电感压降,电流变化引起的
  • ω_e * L_q * i_q:交叉耦合项,d轴和q轴之间的“串扰”
  • ω_e * ψ_f:反电动势项,这是速度估算的关键!

避坑指南: 我曾经在低速时忽略交叉耦合项,结果电流环震荡得一塌糊涂。低速时耦合效应不明显,但转速一高,ω_e * L_q * i_q 这一项会变得很大。如果你不做解耦补偿,PI控制器根本扛不住。

磁链方程:

ψ_d = L_d * i_d + ψ_f
ψ_q = L_q * i_q

这里ψ_f是永磁体产生的磁链,是个常数。你看,d轴磁链由两部分组成:电枢反应(L_d * i_d)加上永磁体磁链。q轴磁链纯粹是电枢反应。

有了磁链,电磁转矩就出来了:

T_e = (3/2) * p * (ψ_d * i_q - ψ_q * i_d)
     = (3/2) * p * [ψ_f * i_q + (L_d - L_q) * i_d * i_q]

这个公式很有意思。第一项是永磁转矩,由q轴电流产生。第二项是磁阻转矩,由d轴和q轴电感差异产生。对于表贴式PMSM(L_d = L_q),磁阻转矩为零。对于内嵌式PMSM(L_d < L_q),你可以利用磁阻转矩来提高效率。

2.5 小结:这些方程怎么用?

你可能会问,搞这么多变换和方程,到底有什么用?

我简单总结一下:

  1. Clark + Park变换:把三相交流量变成两相直流量,方便控制
  2. 电压方程:告诉你电压和电流、转速之间的关系。反电动势项ω_e * ψ_f是速度估算的“金矿”
  3. 磁链方程:告诉你磁链怎么算,转矩怎么来

在无传感器速度估算中,我们就是利用电压方程,从已知的电压和电流中,把转速ω_e“反推”出来。具体怎么推?后面的章节会详细讲。

嗯,数学模型就到这里。你把这些方程吃透了,后面讲观测器、估算器的时候,你会觉得“哦,原来就是这么回事”。


专注资料整理