2、PMSM数学模型：三相静止坐标系下的电压方程、磁链方程、转矩方程与坐标变换基础

永磁同步电机（PMSM）的矢量控制之所以能够实现高性能的转矩与转速调节，其根本前提是建立精确且可解的数学模型。在三相静止坐标系（ABC坐标系）下，PMSM的电压、磁链和转矩方程具有清晰的物理意义，但方程中包含随转子位置变化的时变电感系数，直接求解十分困难。坐标变换（Clark变换与Park变换）正是为了消除这种时变性，将三相交流量等效变换到两相旋转坐标系（dq坐标系）下，使磁链与电流解耦，从而为双闭环控制系统的设计奠定数学基础。本章将系统梳理三相静止坐标系下的基本方程，并深入解析坐标变换的数学原理与工程实现要点。

2.1 三相静止坐标系下的PMSM基本方程

在建立数学模型之前，通常假设电机为理想电机：三相绕组完全对称、气隙磁场正弦分布、忽略磁路饱和、铁耗和涡流损耗不计。基于这些假设，PMSM在三相静止坐标系下的电压方程可表示为矩阵形式：

| uA |   | Rs  0   0  | | iA |   | dψA/dt |
| uB | = | 0   Rs  0  | | iB | + | dψB/dt |
| uC |   | 0   0   Rs | | iC |   | dψC/dt |

其中，Rs为定子相电阻，ψA、ψB、ψC为三相绕组的总磁链。磁链方程则包含两部分：定子电流产生的电枢反应磁链和永磁体产生的励磁磁链。以A相为例，其磁链表达式为：

ψA = LAA·iA + MAB·iB + MAC·iC + ψf·cos(θe)

式中，LAA为A相自感，MAB、MAC为A相与B、C相的互感，ψf为永磁体磁链幅值，θe为转子电角度。由于转子旋转，LAA、MAB等电感系数均随θe呈周期性变化，这使得方程成为一组变系数微分方程，直接求解极为复杂。

转矩方程则基于磁场能量或虚位移法推导，三相静止坐标系下的瞬时电磁转矩为：

Te = p·ψf·[iA·sin(θe) + iB·sin(θe - 2π/3) + iC·sin(θe + 2π/3)]

其中p为极对数。该方程表明，转矩是三相电流与转子位置的复杂函数，控制时需要同时协调三相电流的幅值和相位，这在实际工程中难以直接实现。

2.2 Clark变换：从三相到两相静止坐标系

Clark变换（又称3s/2s变换）的核心思想是：在功率不变或幅值不变的前提下，将三相绕组等效为两相正交绕组（αβ轴）。由于三相绕组在空间上互差120°，其合成磁动势与两相正交绕组在任意时刻产生的磁动势可以完全等效。变换矩阵如下（等幅值形式）：

| iα |   | 1    -1/2   -1/2  | | iA |
| iβ | = | 0    √3/2  -√3/2 | | iB |
                  | iC |

等幅值变换的系数为2/3，其物理含义是变换前后电流矢量的幅值保持不变，便于与电机铭牌参数直接对应。若采用等功率变换，则系数变为√(2/3)，此时变换前后功率计算式保持一致，但电流幅值会变为原来的√(3/2)倍。工程中，TI、ST等主流控制芯片的库函数默认采用等幅值变换，而MATLAB/Simulink的某些模块则采用等功率形式，使用时必须明确标注，否则会导致PI参数计算错误。

Clark变换的物理意义在于：三相电流中存在一个冗余自由度（iA + iB + iC = 0），变换后仅需两个变量iα、iβ即可完整描述电枢磁动势。这为后续Park变换提供了前提，同时也简化了电流采样——只需采样两相电流即可重构第三相。

2.3 Park变换：从静止到旋转坐标系

Park变换（又称2s/2r变换）将αβ轴上的交流分量变换为dq轴上的直流分量。关键在于选择一个以电角速度ωe旋转的坐标系，并使d轴与转子永磁体磁链方向对齐。变换矩阵为：

| id |   | cosθe   sinθe | | iα |
| iq | = | -sinθe  cosθe | | iβ |

经过Park变换后，在稳态工况下id、iq均为直流量，PI调节器可以无静差跟踪。更重要的是，变换后电压方程变为：

ud = Rs·id + Ld·(did/dt) - ωe·Lq·iq
uq = Rs·iq + Lq·(diq/dt) + ωe·(Ld·id + ψf)

式中出现了耦合项-ωe·Lq·iq和ωe·Ld·id，以及反电动势项ωe·ψf。对于表贴式PMSM（Ld = Lq），转矩方程简化为Te = 1.5·p·ψf·iq，此时id = 0控制策略即可实现转矩与电流的线性关系。对于内置式PMSM（Ld ≠ Lq），则需利用磁阻转矩，采用最大转矩电流比（MTPA）控制。

2.4 坐标变换的工程实现与案例拆解

在数字控制系统中，坐标变换以周期性中断方式执行，典型执行频率为10kHz~20kHz。每次中断需完成：ADC采样三相电流→Clark变换→Park变换→PI调节→反Park变换→SVPWM调制。其中反Park变换矩阵为Park变换矩阵的转置（正交矩阵性质），计算量极小。

案例：某3kW伺服驱动器电流环调试中的坐标变换问题

某工业伺服项目，额定电流7.5A，采用2500线增量式编码器（四倍频后10000线）。调试中发现电机空载运行时电流波形存在明显6次谐波，且转矩波动超过额定值的5%。初步排查：

• 电流采样正常，三相电流幅值平衡，无偏置误差；
• Clark变换系数采用等幅值形式，与PI参数计算一致；
• 进一步检查Park变换角度，发现编码器零位与电机电气零位偏差约3°电角度。

分析判断：3°的角度偏差导致id、iq轴未完全解耦，d轴电流分量耦合到q轴，产生周期性转矩脉动。实际操作中，通过自整定程序获取编码器零位偏移量：给定id = 0、iq = 小电流，缓慢旋转转子，记录q轴电流最小时的编码器角度作为零位补偿值。修正后，6次谐波幅值降低约80%，转矩波动降至2%以内。

该案例的代价是额外增加了约10分钟的出厂校准时间，但换来了转矩平稳性的显著提升。复盘经验：坐标变换的精度不仅取决于数学公式，更取决于角度反馈的准确性。在批量生产中，应将该校准步骤纳入出厂测试流程。

2.5 常见误区与总结复盘