【PMSLM建模与控制】PMSLM模型

永磁同步直线电机数学模型

1.1 坐标变换：从三相静止到两相旋转坐标系

永磁同步直线电机（PMSLM）是一个复杂的三相电磁系统，其定子绕组（初级）中的电压、电流和磁链方程都是时变交流量，直接分析与控制非常困难。坐标变换的核心思想在于，通过数学上的坐标映射，将三相静止坐标系（ABC）中的变量变换到两相坐标系中，从而简化数学模型，并实现对电机转矩和磁场的独立控制。本节将详细推导最关键的两种变换：Clarke变换（3s/2s）和Park变换（2s/2r）。

1.1.1 Clarke变换（3s/2s变换）

1. 物理思想
Clarke变换的目标是在保持合成磁动势（或空间矢量）不变的前提下，将三相静止坐标系（A, B, C）中的变量转换到两相静止坐标系（α, β）中。它将互差120°的三个标量变量转换为一个在α-β平面上旋转的矢量，该矢量的两个正交分量即为变换后的变量。

2. 变换前提与约束
假设系统为三相对称且满足三相瞬时值之和为零（即无零序分量），即：
$$i_a(t)+i_b(t)+i_c(t)=0\text{\hspace{1em}(1)}$$
对于PMSLM这类三相三线制系统，此条件通常成立。

3. 推导过程
我们的目标是找到变换矩阵 $C_{3s/2s}$，使得：
$$\left[\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right]=C_{3s/2s}\left[\begin{array}{c}{i}_a\\ i_b\\ i_c\end{array}\right]\text{\hspace{1em}(2)}$$

步骤一：确定α轴与A轴重合
为简化计算，我们令α轴在空间上的方向与A相绕组的轴线重合。β轴则超前α轴90°电角度。

步骤二：磁动势等效原则
根据磁动势等效原则，变换前后合成的基波磁动势应相等。

• 三相绕组（A, B, C）通以电流 $i_a$, $i_b$, $i_c$ 产生的合成磁动势为 $F_{abc}$。
• 两相绕组（α, β）通以电流 $i_{\alpha }$, $i_{\beta }$ 产生的合成磁动势为 $F_{\alpha \beta }$。
  令 $F_{abc}=F_{\alpha \beta }$。

步骤三：投影计算
将三相磁动势 $F_{abc}$ 投影到α轴和β轴上，应与两相磁动势的分量相等。

• 在α轴上的投影：A相磁动势全部在α轴上，B相和C相磁动势在α轴上的投影需要乘以它们与α轴夹角的余弦（$\cos 120°=-1/2$, $\cos 240°=-1/2$）。

$$N_2{i}_{\alpha }=N_3\left[i_a\cdot \cos 0°+i_b\cdot \cos 120°+i_c\cdot \cos 240°\right]=N_3\left(i_a-\frac{1}{2}{i}_b-\frac{1}{2}{i}_c\right)\text{\hspace{1em}(3)}$$

• 在β轴上的投影：A相磁动势在β轴投影为0，B相和C相磁动势在β轴上的投影需要乘以它们与β轴夹角的正弦（$\sin 120°=\sqrt{3}/2$, $\sin 240°=-\sqrt{3}/2$）。

$$N_2{i}_{\beta }=N_3\left[i_a\cdot \sin 0°+i_b\cdot \sin 120°+i_c\cdot \sin 240°\right]=N_3\left(0+\frac{\sqrt{3}}{2}{i}_b-\frac{\sqrt{3}}{2}{i}_c\right)\text{\hspace{1em}(4)}$$

其中，$N_3$ 和 $N_2$ 分别是三相和两相绕组的有效匝数。

步骤四：等幅值变换与化简
为了使变换前后的矢量幅值相等（即 $|i_{\alpha \beta }|=|i_{abc}|$），我们令匝数比 $N_3/N_2=\sqrt{2/3}$。同时，利用 $i_a+i_b+i_c=0$，可得 $i_c=-i_a-i_b$，将其代入α轴投影方程进行化简：
$$\begin{array}{rl}{i}_{\alpha }& =\frac{N_3}{N_2}\left(i_a-\frac{1}{2}{i}_b-\frac{1}{2}{i}_c\right)\\ & =\sqrt{\frac{2}{3}}\left(i_a-\frac{1}{2}{i}_b-\frac{1}{2}(-i_a-i_b)\right)\\ & =\sqrt{\frac{2}{3}}\left(i_a-\frac{1}{2}{i}_b+\frac{1}{2}{i}_a+\frac{1}{2}{i}_b\right)\\ & =\sqrt{\frac{2}{3}}\left(\frac{3}{2}{i}_a\right)\\ & =\sqrt{\frac{3}{2}}\cdot i_a\end{array}\text{\hspace{1em}(5)}$$
实际上，我们更希望保留变换与 $i_b$ 和 $i_c$ 的关系。通常直接使用未代入的表达式并乘以系数。

综合以上，得到等幅值Clarke变换公式：
$$\begin{array}{rl}{i}_{\alpha }& =\sqrt{\frac{2}{3}}\left(i_a-\frac{1}{2}{i}_b-\frac{1}{2}{i}_c\right)\\ i_{\beta }& =\sqrt{\frac{2}{3}}\left(0+\frac{\sqrt{3}}{2}{i}_b-\frac{\sqrt{3}}{2}{i}_c\right)=\sqrt{\frac{2}{3}}\cdot \frac{\sqrt{3}}{2}(i_b-i_c)=\frac{1}{\sqrt{2}}(i_b-i_c)\end{array}\text{\hspace{1em}(6)}$$
写成矩阵形式：
$$\left[\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right]=\sqrt{\frac{2}{3}}\left[\begin{array}{ccc}1& -\frac{1}{2}& -\frac{1}{2}\\ 0& \frac{\sqrt{3}}{2}& -\frac{\sqrt{3}}{2}\end{array}\right]\left[\begin{array}{c}{i}_a\\ i_b\\ i_c\end{array}\right]\text{\hspace{1em}(7)}$$

在实际数字控制中，为减少计算量，常采用忽略零序分量的简化形式（$i_a+i_b+i_c=0$）：
$$\begin{array}{rl}{i}_{\alpha }& =i_a\\ i_{\beta }& =\frac{1}{\sqrt{3}}(i_a+2i_b)\text{或}\frac{1}{\sqrt{3}}(i_b-i_c)\end{array}\text{\hspace{1em}(8)}$$
此变换非等幅值，但计算简单，且最终控制效果经PI调节后并无差异，因而被广泛采用。若要使用等幅值变换，需要使 $N_3/N_2=\sqrt{2}$。

1.1.2 Park变换（2s/2r变换）

1. 物理思想
Park变换的目标是将两相静止坐标系（α, β）中的变量转换到与转子磁场同步旋转的两相旋转坐标系（d, q）上。经过此变换，所有在ABC坐标系中正弦变化的交流量（如电流、电压、磁链）在DQ坐标系下都表现为直流量。这使得控制器可以使用简单的PI调节器来实现无静差控制。

2. 推导过程
d轴（直轴）通常与转子永磁体产生的磁场方向重合，q轴（交轴）超前d轴90°电角度。设转子d轴与静止α轴之间的电角度为 $\theta$（对于旋转电机是旋转角度，对于直线电机是位置 $x$ 除以极距 $\tau$ 再乘以 $\pi$，即 $\theta =\pi x/\tau$）。

变换的任务是：将静止的 $i_{\alpha }$ 和 $i_{\beta }$ 投影到旋转的d轴和q轴上。

步骤一：投影计算
将 $i_{\alpha }$ 和 $i_{\beta }$ 合成的电流空间矢量 $i_s$ 分别向d轴和q轴进行投影。

• d轴电流 $i_d$：是 $i_s$ 在d轴上的投影，等于 $i_{\alpha }$ 和 $i_{\beta }$ 在d轴上投影之和。

$$i_d=i_{\alpha }\cos \theta +i_{\beta }\sin \theta \text{\hspace{1em}(9)}$$

• q轴电流 $i_q$：是 $i_s$ 在q轴上的投影，等于 $i_{\alpha }$ 和 $i_{\beta }$ 在q轴上投影之和。q轴超前d轴90°，故其方向单位向量为 $(-\sin \theta ,\cos \theta )$。

$$i_q=-i_{\alpha }\sin \theta +i_{\beta }\cos \theta \text{\hspace{1em}(10)}$$

步骤二：矩阵形式
将上式写成矩阵形式，即得到Park变换方程：
$$\left[\begin{array}{c}{i}_d\\ i_q\end{array}\right]=\left[\begin{array}{cc}\cos \theta & \sin \theta \\ -\sin \theta & \cos \theta \end{array}\right]\left[\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right]\text{\hspace{1em}(11)}$$
相应的逆Park变换（用于将控制电压从DQ系转换回αβ系）为：
$$\left[\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right]=\left[\begin{array}{cc}\cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{array}\right]\left[\begin{array}{c}{i}_d\\ i_q\end{array}\right]\text{\hspace{1em}(12)}$$

1.1.3 复合变换：从ABC到DQ

结合Clarke变换和Park变换，可以直接得到从三相静止ABC坐标系到两相旋转DQ坐标系的复合变换。
$$\left[\begin{array}{c}{i}_d\\ i_q\end{array}\right]=\left[\begin{array}{cc}\cos \theta & \sin \theta \\ -\sin \theta & \cos \theta \end{array}\right]\cdot \sqrt{\frac{2}{3}}\left[\begin{array}{ccc}1& -\frac{1}{2}& -\frac{1}{2}\\ 0& \frac{\sqrt{3}}{2}& -\frac{\sqrt{3}}{2}\end{array}\right]\left[\begin{array}{c}{i}_a\\ i_b\\ i_c\end{array}\right]\text{\hspace{1em}(13)}$$

将其展开后，得到最常用的形式：
$$\begin{array}{rl}{i}_d& =\sqrt{\frac{2}{3}}\left[i_a\cos \theta +i_b\cos (\theta -\frac{2\pi }{3})+i_c\cos (\theta +\frac{2\pi }{3})\right]\\ i_q& =\sqrt{\frac{2}{3}}\left[-i_a\sin \theta -i_b\sin (\theta -\frac{2\pi }{3})-i_c\sin (\theta +\frac{2\pi }{3})\right]\end{array}\text{\hspace{1em}(14)}$$

1.1.4 小结

坐标变换是矢量控制的灵魂。通过Clarke变换，我们将复杂的三相变量简化为两相正交变量；通过Park变换，我们将时变的交流量转化为恒定的直流量，极大地简化了PMSLM的分析与控制设计，为实现高性能的转矩和推力控制奠定了坚实的数学基础。在后续的章节中，我们将基于DQ坐标系建立PMSLM的数学模型。