伺服电机仿真（4）：PMSM在d-q旋转坐标系下的状态方程与等效电路
4.1 引言：从复杂到时不变的革命性转变
在第三部分，我们建立了PMSM在三相静止坐标系下的完整数学模型。这个模型虽然精确描述了电机的物理本质，但其强耦合、非线性、时变的特性使得控制器设计变得异常困难。通过Clark变换和Park变换，我们将这个复杂系统转换到了与转子同步旋转的d-q坐标系，实现了从时变交流系统到时不变直流系统的根本转变。
本部分将系统性地推导PMSM在d-q旋转坐标系下的状态方程，建立等效电路模型，并分析其数学特性和物理意义。这个模型是现代矢量控制（FOC）的理论基石。
4.2 d-q旋转坐标系下的电压方程
经过坐标变换，PMSM在d-q旋转坐标系下的电压方程变得简洁而富有物理意义：
4.2.1 基本电压方程
d-q旋转坐标系电压方程结构 ┌─────────────────────────────────────────────────────────────┐ │ d轴电压方程：u_d = R_s·i_d + dψ_d/dt - ω_e·ψ_q │ │ q轴电压方程：u_q = R_s·i_q + dψ_q/dt + ω_e·ψ_d │ │ │ │ 各分量物理意义： │ │ 1. R_s·i_d, R_s·i_q：电阻压降 │ │ 2. dψ_d/dt, dψ_q/dt：变压器电动势（感应电动势） │ │ 3. -ω_e·ψ_q, +ω_e·ψ_d：旋转电动势（速度电动势） │ └─────────────────────────────────────────────────────────────┘4.2.2 方程详解
d轴电压方程：
ud=Rsid+dtdψd−ωeψq
q轴电压方程：
uq=Rsiq+dtdψq+ωeψd
旋转电动势项的物理意义：
当转子旋转时，d轴磁链ψ_d在q轴绕组中切割产生电动势：+ω_e·ψ_d
当转子旋转时，q轴磁链ψ_q在d轴绕组中切割产生电动势：-ω_e·ψ_q
这些项体现了d轴和q轴之间的耦合，是矢量控制中需要前馈补偿的关键项
4.3 d-q旋转坐标系下的磁链方程
在d-q坐标系下，磁链方程得到了极大的简化，实现了完全解耦：
4.3.1 基本磁链方程
d-q坐标系磁链关系 ┌─────────────────────────────────────────────────────────────┐ │ d轴磁链：ψ_d = L_d·i_d + ψ_f │ │ q轴磁链：ψ_q = L_q·i_q │ │ │ │ 参数说明： │ │ L_d：直轴电感（d轴方向磁阻最小） │ │ L_q：交轴电感（q轴方向磁阻较大） │ │ ψ_f：永磁体磁链（恒定值，方向在d轴上） │ │ │ │ 重要特性： │ │ 1. 完全解耦：d轴磁链只与i_d有关，q轴磁链只与i_q有关 │ │ 2. 参数恒定：L_d、L_q、ψ_f均为常数，不随时间变化 │ └─────────────────────────────────────────────────────────────┘4.3.2 电机类型对电感的影响
不同类型PMSM的电感特性 ┌───────────────────┬─────────────────┬─────────────────────┐ │ 电机类型 │ L_d │ L_q │ ├───────────────────┼─────────────────┼─────────────────────┤ │ 表贴式(SPM) │ L_d = L_q │ L_q = L_d │ │ (圆柱形转子) │ (气隙均匀) │ (无凸极效应) │ ├───────────────────┼─────────────────┼─────────────────────┤ │ 内置式(IPM) │ L_d < L_q │ L_q > L_d │ │ (凸极转子) │ (d轴磁路短) │ (q轴磁路经过磁桥) │ └───────────────────┴─────────────────┴─────────────────────┘4.4 完整的电压-电流状态方程
将磁链方程代入电压方程，得到用电流表示的状态方程：
4.4.1 标准形式状态方程
d-q轴电压-电流动态方程 ┌─────────────────────────────────────────────────────────────┐ │ d轴：u_d = R_s·i_d + L_d·di_d/dt - ω_e·L_q·i_q │ │ q轴：u_q = R_s·i_q + L_q·di_q/dt + ω_e·(L_d·i_d + ψ_f) │ │ │ │ 状态变量：i_d, i_q │ │ 输入变量：u_d, u_q │ │ 扰动变量：ω_e (可作为已知量或状态) │ └─────────────────────────────────────────────────────────────┘4.4.2 矩阵形式状态方程
[uduq]=[RsωeLd−ωeLqRs][idiq]+[Ld00Lq]dtd[idiq]+[0ωeψf]
状态空间表达：
dtd[idiq]=[−LdRs−LqωeLdLdωeLq−LqRs][idiq]+[Ld100Lq1][uduq]−[0Lqωeψf]
4.5 电磁转矩方程
在d-q坐标系下，电磁转矩方程变得异常简洁直观：
4.5.1 通用转矩方程
Te=23p[ψfiq+(Ld−Lq)idiq]
4.5.2 转矩组成分析
电磁转矩的两种成分 ┌─────────────────────────────────────────────────────────────┐ │ 电磁转矩T_e = 永磁转矩 + 磁阻转矩 │ │ │ │ 1. 永磁转矩：T_pm = (3/2)·p·ψ_f·i_q │ │ - 由永磁磁场与q轴电流相互作用产生 │ │ - 与i_q成正比，是表贴式电机的主要转矩 │ │ │ │ 2. 磁阻转矩：T_rel = (3/2)·p·(L_d - L_q)·i_d·i_q │ │ - 由转子凸极效应产生（L_d ≠ L_q） │ │ - 对于内置式电机(IPM)，L_d < L_q，此项为负 │ │ - 通过合理控制i_d、i_q可优化转矩输出 │ └─────────────────────────────────────────────────────────────┘4.5.3 特殊情况的转矩方程
表贴式电机(SPM)：Ld=Lq，转矩方程简化为：
Te=23pψfiq
内置式电机(IPM)：Ld<Lq，存在磁阻转矩，总转矩为：
Te=23p[ψfiq+(Ld−Lq)idiq]
由于 (Ld−Lq)<0，为产生正的磁阻转矩，需要 id和 iq异号。
4.6 机械运动方程
机械运动方程与坐标系无关，但在d-q模型框架下，与电气方程通过转矩耦合：
4.6.1 基本运动方程
Jdtdωm+Bωm=Te−TL
dtdθm=ωm,θe=p⋅θm,ωe=p⋅ωm
4.6.2 机电系统耦合结构
PMSM机电系统耦合框图 ┌─────────────────────────────────────────────────────────────┐ │ 电气子系统 机械子系统 │ │ ┌─────────────────┐ ┌─────────────────┐ │ │ │ 输入：u_d, u_q │ │ 输入：T_e │ │ │ │ 状态：i_d, i_q │ │ 状态：ω_m, θ_m │ │ │ │ 输出：T_e ├───────►│ 输出：ω_e, θ_e │ │ │ └─────────────────┘ └────────┬────────┘ │ │ │ │ │ ┌─────▼─────┐ │ │ │ 坐标变换 │ │ │ │ 反馈回路 │ │ │ └─────┬─────┘ │ │ │ │ │ ┌─────▼─────┐ │ │ │ ω_e, θ_e │ │ │ └─────┬─────┘ │ │ │ │ │ ┌─────▼─────┐ │ │ │ 前馈补偿 │ │ │ │ 与解耦 │ │ │ └───────────┘ │ └─────────────────────────────────────────────────────────────┘4.7 d轴和q轴等效电路模型
基于d-q轴电压方程，可以建立直观的等效电路模型，这有助于理解各物理量的作用和设计控制器。
4.7.1 d轴等效电路
d轴等效电路模型 ┌─────────────────────────────────────────────────────────────┐ │ u_d ────┬──────[R_s]──────[L_d]──────┤├───┐ │ │ │ │ │ │ │ ├────( -ω_e·L_q·i_q )──────┐ │ │ │ │ │ (旋转电动势) │ │ │ │ │ │ │ │ │ │ │ └─────( 反电动势耦合 )─────┘ │ │ │ │ │ │ │ │ === ψ_f·d/dt? │ │ │ │ │ │ └────┘ │ │ │ │ 物理量说明： │ │ 1. R_s, L_d：d轴绕组的电阻和电感 │ │ 2. -ω_e·L_q·i_q：q轴电流在d轴产生的旋转电动势 │ │ 3. 永磁磁链ψ_f不直接出现在d轴电压方程，但通过 │ │ ω_e·ψ_f项影响q轴方程 │ └─────────────────────────────────────────────────────────────┘4.7.2 q轴等效电路
q轴等效电路模型 ┌─────────────────────────────────────────────────────────────┐ │ u_q ────┬──────[R_s]──────[L_q]──────┤├───┐ │ │ │ │ │ │ │ ├────( +ω_e·L_d·i_d )──────┐ │ │ │ │ │ (旋转电动势1) │ │ │ │ │ │ │ │ │ │ │ ├────( +ω_e·ψ_f )────────┐ │ │ │ │ │ │ (旋转电动势2) │ │ │ │ │ │ │ │ │ │ │ │ │ └─────( 反电动势耦合 )───┘ │ │ │ │ │ │ │ │ │ │ │ === === ??? │ │ │ │ │ │ │ │ └─┴─┘ │ │ │ │ │ │ └────┘ │ │ │ │ 物理量说明： │ │ 1. R_s, L_q：q轴绕组的电阻和电感 │ │ 2. +ω_e·L_d·i_d：d轴电流在q轴产生的旋转电动势 │ │ 3. +ω_e·ψ_f：永磁体旋转在q轴产生的反电动势 │ │ (这是永磁同步电机最主要的反电动势源) │ └─────────────────────────────────────────────────────────────┘4.8 完整的状态空间模型
结合电气和机械方程，得到PMSM在d-q坐标系下的完整状态空间模型：
4.8.1 状态变量选择
系统有四个主要状态变量：
电气状态：d轴电流 id，q轴电流 iq
机械状态：转子角速度 ωm，转子位置 θm
4.8.2 状态方程汇总
⎩⎨⎧dtdid=−LdRsid+LdωeLqiq+Ld1uddtdiq=−LqωeLdid−LqRsiq+Lq1uq−Lqωeψfdtdωm=J1(23p[ψfiq+(Ld−Lq)idiq]−Bωm−TL)dtdθm=ωmωe=p⋅ωm
4.8.3 模型框图表示
PMSM d-q坐标系完整状态模型 ┌─────────────────────────────────────────────────────────────────────┐ │ 系统输入 │ ├─────────────────────────────────────────────────────────────────────┤ │ 控制输入：u_d, u_q 扰动输入：T_L │ ├─────────────────────────────────────────────────────────────────────┤ │ 状态空间模型 │ ├─────────────────────────────────────────────────────────────────────┤ │ ┌─────────────────────────────────────────────────────────────┐ │ │ │ 电气子系统状态方程 │ │ │ │ ┌─────────────────────────────────────────────────────┐ │ │ │ │ │ di_d/dt = f1(i_d, i_q, ω_e, u_d) │ │ │ │ │ │ di_q/dt = f2(i_d, i_q, ω_e, u_q) │ │ │ │ │ └──────────────┬──────────────────────┬──────────────┘ │ │ │ │ │ │ │ │ │ │ i_d, i_q │ ω_e │ │ │ │ │ │ │ │ │ │ ┌──────────────▼──────────────────────▼──────────────┐ │ │ │ │ │ 转矩计算：T_e = f3(i_d, i_q) │ │ │ │ │ └──────────────────────────┬─────────────────────────┘ │ │ │ │ │ T_e │ │ │ └─────────────────────────────┼─────────────────────────────┘ │ │ │ │ │ ┌─────────────────────────────▼─────────────────────────────┐ │ │ │ 机械子系统状态方程 │ │ │ │ ┌─────────────────────────────────────────────────────┐ │ │ │ │ │ dω_m/dt = (T_e - B·ω_m - T_L)/J │ │ │ │ │ │ dθ_m/dt = ω_m │ │ │ │ │ │ ω_e = p·ω_m │ │ │ │ │ └──────┬──────────────────────┬──────────────────────┘ │ │ │ │ │ │ │ │ │ │ ω_m θ_m, ω_e │ │ │ │ │ │ │ │ │ └─────────┼──────────────────────┼─────────────────────────┘ │ │ │ │ │ │ └──────────────────────┘ │ │ 反馈到电气子系统 │ ├─────────────────────────────────────────────────────────────────────┤ │ 系统输出 │ ├─────────────────────────────────────────────────────────────────────┤ │ 电气输出：i_d, i_q, T_e 机械输出：ω_m, θ_m, ω_e │ └─────────────────────────────────────────────────────────────────────┘4.9 模型特性与简化
4.9.1 d-q模型的优势
d-q旋转坐标系模型的革命性优势 ┌─────────────────────────────────────────────────────────────┐ │ 1. 时不变性：所有参数(L_d, L_q, R_s, ψ_f)均为常数 │ │ 2. 解耦性：d轴和q轴方程基本独立，便于设计独立控制器 │ │ 3. 直流化：稳态时i_d, i_q, u_d, u_q为直流量 │ │ 4. 物理直观：d轴控制磁场，q轴控制转矩，概念清晰 │ │ 5. 控制简化：可直接应用成熟的直流电机控制理论 │ └─────────────────────────────────────────────────────────────┘4.9.2 模型简化情况
忽略次要因素时的简化模型：
忽略粘性摩擦：B=0
稳态分析：did/dt=0, diq/dt=0
小信号线性化：在工作点附近线性化
表贴式电机的特殊简化：
对于SPMSM（Ld=Lq=Ls），模型简化为：
⎩⎨⎧ud=Rsid+Lsdtdid−ωeLsiquq=Rsiq+Lsdtdiq+ωe(Lsid+ψf)Te=23pψfiq
id=0控制策略下的简化：
许多应用中采用 id=0控制，此时模型进一步简化为：
⎩⎨⎧ud=−ωeLqiq(仅为耦合项)uq=Rsiq+Lqdtdiq+ωeψfTe=23pψfiq
4.10 总结
本部分系统性地建立了PMSM在d-q旋转坐标系下的完整数学模型：
4.10.1 核心成果
电压方程：包含电阻压降、感应电动势和旋转电动势
磁链方程：实现d轴和q轴的完全解耦
转矩方程：明确分离了永磁转矩和磁阻转矩
状态空间模型：统一了电气和机械动态
等效电路：提供了直观的物理理解工具
4.10.2 模型的重要意义
d-q模型在伺服控制中的核心地位 ┌─────────────────────────────────────────────────────────────┐ │ 理论层面 │ 工程层面 │ ├─────────────────────────────────────────────────────────────┤ │ 1. 建立了交流电机与直流电机的 │ 1. 为矢量控制(FOC)提供 │ │ 数学等价关系 │ 直接理论依据 │ │ │ │ │ 2. 揭示了磁场与转矩的解耦控制 │ 2. 使PI调节器设计成为可能 │ │ 原理 │ │ │ │ │ │ 3. 提供了系统分析和设计的统一 │ 3. 简化了控制器参数整定 │ │ 框架 │ 和优化过程 │ │ │ │ │ 4. 为高性能控制算法（如MPC、 │ 4. 便于实现弱磁控制、 │ │ 自适应控制）奠定基础 │ MTPA等先进策略 │ └─────────────────────────────────────────────────────────────┘4.10.3 后续应用
我们将基于这个d-q模型，深入探讨永磁同步电机矢量控制（FOC）原理与设计，包括电流环、速度环、位置环的设计方法，以及SVPWM调制技术的实现，最终完成从数学模型到控制系统的完整构建。
关键洞见：坐标变换不仅仅是一个数学技巧，它从根本上改变了我们观察和控制交流电机的方式，使得高性能的伺服控制成为可能。d-q模型的美妙之处在于，它将复杂的时变交流系统转化为我们熟悉的时不变直流系统，为控制工程师打开了通往高性能伺服驱动的大门。