伺服电机仿真（3）：PMSM在三相静止坐标系下的数学模型
3.1 引言：从物理实体到数学描述
建立精确的数学模型是进行伺服电机仿真的第一步。在三相静止坐标系（ABC坐标系）下建立的PMSM模型，直接反映了电机的物理本质，描述了定子三相绕组与旋转转子永磁体之间的完整电磁相互作用关系。这个模型虽然复杂，但它是所有后续分析、坐标变换和控制设计的起点。
3.2 建模前提与假设
为建立可行的数学模型，我们对实际电机作以下合理简化：
建模假设前提 ┌─────────────────────────────────────────────────────────────┐ │ 1. 磁路线性：忽略磁饱和，电感为常数 │ │ 2. 忽略铁损：不计涡流和磁滞损耗 │ │ 3. 三相对称：绕组在空间互差120°，参数完全相同 │ │ 4. 正弦分布：气隙磁场按正弦规律分布 │ │ 5. 无阻尼绕组：转子为实心永磁体 │ │ 6. 电动机惯例：电流流入为正，转矩驱动为正 │ └─────────────────────────────────────────────────────────────┘3.3 坐标系与变量定义
3.3.1 三相静止坐标系
以定子A相绕组轴线为参考轴（0°位置）
B相轴线滞后A相120°电角度
C相轴线滞后A相240°电角度
转子d轴（永磁体N极轴线）与A相轴线的夹角为θ_e
3.3.2 关键物理量符号
电气量： u_A, u_B, u_C: 三相端电压 i_A, i_B, i_C: 三相相电流 ψ_A, ψ_B, ψ_C: 三相总磁链 ψ_f: 永磁体磁链幅值 R_s: 定子相电阻 L_{AA}, L_{BB}, L_{CC}: 自感 L_{AB}, L_{BC}, L_{CA}: 互感 机械量： T_e: 电磁转矩 T_L: 负载转矩 J: 转动惯量 B: 粘性摩擦系数 θ_e: 转子电角度 ω_e: 转子电角速度 p: 极对数3.4 电压方程
基于基尔霍夫电压定律，三相电压方程描述了端电压、电阻压降和感应电动势的关系：
电压方程结构框图 ┌─────────────────────────────────────────────────────┐ │ A相电压平衡 │ │ u_A ───────────────────────────────────────┐ │ │ │ │ │ │ ▼ R_s·i_A (电阻压降) │ │ │ ┌──────┐ │ │ │ │ R_s │ │ │ │ └──────┘ │ │ │ │ │ │ │ ▼ dψ_A/dt (感应电动势) │ │ │ ┌──────┐ │ │ │ │ d/dt│ │ │ │ └──────┘ │ │ │ │ │ │ │ └────────────────────────────────────┘ │ │ ∑ = 0 (KVL) │ │ │ │ B相、C相同理可得 │ └─────────────────────────────────────────────────────┘数学表达式：
⎩⎨⎧uA=RsiA+dtdψAuB=RsiB+dtdψBuC=RsiC+dtdψC
矩阵形式：
uAuBuC=RsiAiBiC+dtdψAψBψC
3.5 磁链方程
磁链方程揭示了系统耦合与时变的物理本质。每相磁链由两部分组成：电流产生的磁链（自感+互感）和永磁体产生的磁链。
3.5.1 磁链组成结构
三相磁链生成机制 ┌─────────────────────────────────────────────────────────┐ │ A相总磁链ψ_A │ │ ┌──────────────────────────────────────────────────┐ │ │ │ 1. 电流产生的磁链 │ │ │ │ ┌─────────┐ ┌─────────┐ ┌─────────┐ │ │ │ │ │ L_AA·i_A│ + │ L_AB·i_B│ + │ L_AC·i_C│ │ │ │ │ └─────────┘ └─────────┘ └─────────┘ │ │ │ │ 自感磁链 互感磁链 互感磁链 │ │ │ └──────────────────────────────────────────────────┘ │ │ ┌──────────────────────────────────────────────────┐ │ │ │ 2. 永磁体产生的磁链 │ │ │ │ ψ_f·cos(θ_e) │ │ │ └──────────────────────────────────────────────────┘ │ │ │ │ B相：ψ_B = L_BA·i_A + L_BB·i_B + L_BC·i_C + ψ_f·cos(θ_e-120°) │ │ C相：ψ_C = L_CA·i_A + L_CB·i_B + L_CC·i_C + ψ_f·cos(θ_e+120°) │ └─────────────────────────────────────────────────────────┘3.5.2 电感矩阵的时变性
由于转子旋转，电感参数是转子位置θ_e的函数：
电感参数随转子位置变化 ┌─────────────────────────────────────────────────────┐ │ 自感变化规律：L_AA(θ_e) = L_s0 + L_s2·cos(2θ_e) │ │ 互感变化规律：L_AB(θ_e) = M_s0 + L_s2·cos(2θ_e-120°)│ │ │ │ 表贴式电机(SPM)：L_s2 ≈ 0，电感为常数 │ │ 内置式电机(IPM)：L_s2 ≠ 0，电感随位置变化 │ └─────────────────────────────────────────────────────┘完整磁链方程：
ψAψBψC=LAA(θe)LBA(θe)LCA(θe)LAB(θe)LBB(θe)LCB(θe)LAC(θe)LBC(θe)LCC(θe)iAiBiC+ψfcos(θe)cos(θe−2π/3)cos(θe+2π/3)
3.6 电磁转矩方程
电磁转矩是机电能量转换的枢纽，由磁场能量对机械角位移的偏导数求得：
电磁转矩生成原理 ┌─────────────────────────────────────────────────────────┐ │ 电磁转矩T_e = 磁阻转矩 + 永磁转矩 │ │ ┌──────────────────────────────────────────────────┐ │ │ │ 磁阻转矩：T_rel = ½·p·i^T·(∂L/∂θ_m)·i │ │ │ │ - 源于电感变化(凸极效应) │ │ │ │ - 内置式电机(IPM)存在此项 │ │ │ └──────────────────────────────────────────────────┘ │ │ ┌──────────────────────────────────────────────────┐ │ │ │ 永磁转矩：T_pm = p·i^T·(∂ψ_pm/∂θ_m) │ │ │ │ - 定子电流与永磁磁场相互作用 │ │ │ │ - 所有永磁电机都有此项 │ │ │ └──────────────────────────────────────────────────┘ │ └─────────────────────────────────────────────────────────┘数学表达式：
Te=21p⋅[iAiBiC]∂θm∂L(θe)iAiBiC+p⋅[iAiBiC]∂θm∂ψpm(θe)
3.7 机械运动方程
描述转子机械运动的动力学方程：
旋转系统动力学 ┌─────────────────────────────────────────────────────┐ │ 转矩平衡关系 │ │ │ │ 电磁转矩T_e ──┐ │ │ ├─►[∑]─────►[J·dω_m/dt]───►ω_m │ │ 负载转矩T_L ──┘ │ │ │ │ [B·ω_m]◄─────────┘ │ │ │ │ 运动学关系：dθ_m/dt = ω_m │ │ 电-机械关系：θ_e = p·θ_m, ω_e = p·ω_m │ └─────────────────────────────────────────────────────┘数学表达式：
Jdtdωm+Bωm=Te−TL
dtdθm=ωm,θe=p⋅θm
3.8 完整数学模型结构框图
综合以上方程，得到PMSM在三相静止坐标系下的完整数学模型框图：
PMSM三相静止坐标系完整数学模型 ┌─────────────────────────────────────────────────────────────────────┐ │ 系统输入 │ ├─────────────────────────────────────────────────────────────────────┤ │ 三相电压：u_A, u_B, u_C 负载转矩：T_L │ ├─────────────────────────────────────────────────────────────────────┤ │ 核心数学模型 │ ├─────────────────────────────────────────────────────────────────────┤ │ │ │ ┌─────────────────────────────────────────────────────────────┐ │ │ │ 电压方程模块 │ │ │ │ 输入：ψ_A,ψ_B,ψ_C, i_A,i_B,i_C │ │ │ │ 输出：dψ_A/dt, dψ_B/dt, dψ_C/dt │ │ │ │ (实现：u_ABC = R_s·i_ABC + dψ_ABC/dt) │ │ │ └──────────────┬──────────────────────────────────────────────┘ │ │ │ │ │ ┌──────────────▼──────────────────────────────────────────────┐ │ │ │ 磁链方程模块 │ │ │ │ ┌──────────────────────────────────────────────────────┐ │ │ │ │ │ 电感计算子模块 │ │ │ │ │ │ 输入：θ_e │ │ │ │ │ │ 输出：L_AA(θ_e), L_AB(θ_e), ... L_CC(θ_e) │ │ │ │ │ └──────────────────────────────────────────────────────┘ │ │ │ │ ┌──────────────────────────────────────────────────────┐ │ │ │ │ │ 永磁磁链子模块 │ │ │ │ │ │ 输入：θ_e │ │ │ │ │ │ 输出：ψ_f·cos(θ_e), ψ_f·cos(θ_e-120°), ... │ │ │ │ │ └──────────────────────────────────────────────────────┘ │ │ │ │ 输入：i_A,i_B,i_C, θ_e │ │ │ │ 输出：ψ_A,ψ_B,ψ_C │ │ │ │ (实现：ψ_ABC = L_ABC(θ_e)·i_ABC + ψ_pm(θ_e)) │ │ │ └──────────────┬──────────────────────────────────────────────┘ │ │ │ │ │ ┌──────────────▼──────────────────────────────────────────────┐ │ │ │ 电流求解模块 │ │ │ │ (从u_ABC和ψ_ABC求解i_ABC，需处理代数环) │ │ │ │ 输入：u_A,u_B,u_C, θ_e │ │ │ │ 输出：i_A,i_B,i_C │ │ │ └──────────────┬──────────────────────────────────────────────┘ │ │ │ │ │ ┌──────────────▼──────────────────────────────────────────────┐ │ │ │ 转矩计算模块 │ │ │ │ ┌──────────────────────────────────────────────────────┐ │ │ │ │ │ 磁阻转矩计算 │ │ │ │ │ │ T_rel = ½·p·i^T·(∂L/∂θ_m)·i │ │ │ │ │ └──────────────────────────────────────────────────────┘ │ │ │ │ ┌──────────────────────────────────────────────────────┐ │ │ │ │ │ 永磁转矩计算 │ │ │ │ │ │ T_pm = p·i^T·(∂ψ_pm/∂θ_m) │ │ │ │ │ └──────────────────────────────────────────────────────┘ │ │ │ │ 输入：i_A,i_B,i_C, θ_e │ │ │ │ 输出：T_e = T_rel + T_pm │ │ │ └──────────────┬──────────────────────────────────────────────┘ │ │ │ │ │ ┌──────────────▼──────────────────────────────────────────────┐ │ │ │ 机械运动模块 │ │ │ │ 输入：T_e, T_L │ │ │ │ 输出：ω_m, θ_m, ω_e, θ_e │ │ │ │ (实现：J·dω_m/dt = T_e - T_L - B·ω_m) │ │ │ └──────────────┬──────────────────────────────────────────────┘ │ │ │ │ │ └──────────────────────────────────────────────┐ │ │ │ │ │ ┌─────────────────────────────────────────────────────────────┘ │ │ │ 反馈回路 │ │ │ │ θ_e, ω_e 反馈到电感计算、永磁磁链、转矩计算模块 │ │ │ │ i_A,i_B,i_C 反馈到电压方程、磁链方程模块 │ │ │ └─────────────────────────────────────────────────────────────┘ │ ├─────────────────────────────────────────────────────────────────────┤ │ 系统输出 │ ├─────────────────────────────────────────────────────────────────────┤ │ 三相电流：i_A, i_B, i_C 电磁转矩：T_e │ │ 转子位置：θ_m, θ_e 转子速度：ω_m, ω_e │ │ 各相磁链：ψ_A, ψ_B, ψ_C │ └─────────────────────────────────────────────────────────────────────┘3.9 模型特性分析
3.9.1 复杂性来源
PMSM三相模型复杂性分析 ┌─────────────────────────────────────────────────────────┐ │ 1. 强耦合性 │ │ - 磁链方程中，每相磁链受三相电流共同影响 │ │ - 电感矩阵为满阵，存在相互耦合 │ │ │ │ 2. 时变性 │ │ - 电感L_AA, L_AB, ... 是θ_e的函数 │ │ - 永磁磁链ψ_pm也是θ_e的函数 │ │ - 参数随时间周期性变化 │ │ │ │ 3. 非线性 │ │ - 转矩方程包含电流的二次项(i^T·L·i) │ │ - 运动方程与电磁方程通过转矩耦合 │ └─────────────────────────────────────────────────────────┘3.9.2 仿真挑战
代数环问题：电压方程需要dψ/dt，而ψ是i的函数，i又需要通过电压方程求解
数值稳定性：时变参数可能导致数值计算困难
计算复杂度：实时计算3×3时变电感矩阵及其导数
3.10 总结
本部分建立了PMSM在三相静止坐标系下的完整数学模型，包括：
电压方程：描述电路特性
磁链方程：揭示电磁耦合本质
转矩方程：体现机电能量转换
运动方程：描述机械动力学
该模型的价值：
作为物理理解的基准，揭示电机内部真实物理过程
作为仿真验证的"黄金标准"，用于检验简化模型的正确性
在故障分析、不对称运行等特殊场合仍需使用
模型的局限性：
过于复杂，难以直接用于控制器设计
时变参数使系统分析困难
强耦合性阻碍了独立控制回路的建立
正是由于这些复杂性，在下一部分我们将通过坐标变换，将这个复杂的模型转换到dq旋转坐标系，从而得到一个参数恒定、方程解耦的简化模型，为高性能的矢量控制奠定数学基础。