【永磁同步电机数学模型全程推导】【7 转矩方程】
目录
- 等功率变换
- 转矩方程推导
在物理学中,转矩乘转速等于功率,而转速是已知的,所以可以从转矩相关的功率入手进行推导。
电路中输入给电机的总功率与dq轴的功率之和( udid+uqiqu_di_d+u_qi_qudid+uqiq)有关,注意此处是输入的总功率,总功率不是全部都用于生成转矩的。
等功率变换
在永磁同步电机数学模型的推导过程中,我们总是希望三个坐标系(abc相、αβ轴、dq轴)的尺度(幅值)是一致的,不仅是电压电流幅值一致,功率的幅值也要一致,这些一致是人为乘以一个系数造成的,完全是为了公式好看、推导方便、代码实现方便,关于等幅值和等功率变换的讲解,可以前往查看【无刷电机FOC进阶基础准备】【04 clark变换、park变换、等幅值变换、等功率变换】。αβ轴和dq轴只差一个不会引起尺度变化的park变换,因此只需要αβ轴和abc相的尺度一致即可。
本系列教程的dq轴电压电流都是等幅值变换得到的,
udid+uqiq=uαiα+uβiβ=[uα,uβ][iαiβ]=(23MclarkT∗[uaubuc])T∗23Mclark[iaibic]=49(MclarkT∗[uaubuc])T∗Mclark[iaibic]\begin{aligned}u_di_d+u_qi_q=u_{\alpha}i_{\alpha}+u_{\beta}i_{\beta}=[u_{\alpha},u_{\beta}]\begin{bmatrix} i_{\alpha} \\ i_{\beta} \end{bmatrix}&=(\frac{2}{3}\mathrm{M}_{clark}^T*\begin{bmatrix} u_a \\ u_b \\ u_c \end{bmatrix})^T*\frac{2}{3}\mathrm{M}_{clark}\begin{bmatrix} i_a \\ i_b \\ i_c \end{bmatrix}\\&=\frac{4}{9}(\mathrm{M}_{clark}^T*\begin{bmatrix} u_a \\ u_b \\ u_c \end{bmatrix})^T*\mathrm{M}_{clark}\begin{bmatrix} i_a \\ i_b \\ i_c \end{bmatrix}\end{aligned} udid+uqiq=uαiα+uβiβ=[uα,uβ][iαiβ]=(32MclarkT∗uaubuc)T∗32Mclarkiaibic=94(MclarkT∗uaubuc)T∗Mclarkiaibic
而等功率变换的情况是:
(23MclarkT∗[uaubuc])T∗23Mclark[iaibic]=23(MclarkT∗[uaubuc])T∗Mclark[iaibic](\sqrt{\frac{2}{3}}\mathrm{M}_{clark}^T*\begin{bmatrix} u_a \\ u_b \\ u_c \end{bmatrix})^T*\sqrt{\frac{2}{3}}\mathrm{M}_{clark}\begin{bmatrix} i_a \\ i_b \\ i_c \end{bmatrix}=\frac{2}{3}(\mathrm{M}_{clark}^T*\begin{bmatrix} u_a \\ u_b \\ u_c \end{bmatrix})^T*\mathrm{M}_{clark}\begin{bmatrix} i_a \\ i_b \\ i_c \end{bmatrix} (32MclarkT∗uaubuc)T∗32Mclarkiaibic=32(MclarkT∗uaubuc)T∗Mclarkiaibic
可以看到,两者相差了23\frac{2}{3}32,因此需要人为乘以32\frac{3}{2}23,使得等幅值变成等功率:
P=32(udid+uqiq)P=\frac{3}{2}(u_di_d+u_qi_q)P=23(udid+uqiq)
转矩方程推导
udu_dud和uqu_quq的公式在之前已经推导得到了,在这里代入即可:
P=32(udid+uqiq)=32[(Rsid+Lddiddt−ωLqiq)id+(Rsiq+Lqdiqdt+ωLdid+ωψf)iq]=32Rs(id2+iq2)⏟Pcu铜耗+32(Ldiddiddt+Lqiqdiqdt)⏟磁场速率相关+32(ωLdidiq−ωLqiqid+ωψfiq)⏟Pem电磁功率\begin{aligned}P&=\frac{3}{2}(u_di_d+u_qi_q)=\frac{3}{2}[(R_si_d+L_d\frac{di_d}{dt}-\omega L_qi_q)i_d+(R_si_q+L_q\frac{di_q}{dt}+\omega L_di_d+\omega\psi_f)i_q]\\&=\underbrace{\frac{3}{2}R_s(i^2_d+i^2_q)}_{P_{cu}铜耗}+\underbrace{\frac{3}{2}(L_di_d\frac{di_d}{dt}+L_qi_q\frac{di_q}{dt})}_{磁场速率相关}+\underbrace{\frac{3}{2}(\omega L_di_di_q-\omega L_qi_qi_d+\omega\psi_fi_q)}_{P_{em}电磁功率}\end{aligned}P=23(udid+uqiq)=23[(Rsid+Lddtdid−ωLqiq)id+(Rsiq+Lqdtdiq+ωLdid+ωψf)iq]=Pcu铜耗23Rs(id2+iq2)+磁场速率相关23(Ldiddtdid+Lqiqdtdiq)+Pem电磁功率23(ωLdidiq−ωLqiqid+ωψfiq)
可以看到,输入功率主要由三部分组成。铜耗PcuP_{cu}Pcu纯是发热的,不产生转矩;磁场速率相关的也不产生转矩,因为当堵转时dq轴电流变化率为0;只有电磁功率PemP_{em}Pem才产生转矩。
将电磁功率PemP_{em}Pem摘出来,并化简一下,由于此处的ω\omegaω是电角度,所以写作ωe\omega_eωe:
Pem=32ωe(ψfiq+(Ld−Lq)idiq)P_{em}=\frac{3}{2}\omega_e(\psi_fi_q+(L_d-L_q)i_di_q) Pem=23ωe(ψfiq+(Ld−Lq)idiq)
在物理中,电机转矩相关的功率可以等于转矩乘转速,注意此处是物理转速:
Pem=T∗ωphyP_{em}=T*\omega_{phy}Pem=T∗ωphy
两个PemP_{em}Pem相等,于是转矩方程就得出了:
T∗ωphy=32ωe(ψfiq+(Ld−Lq)idiq)T=32ωeωphy(ψfiq+(Ld−Lq)idiq)=32p(ψfiq+(Ld−Lq)idiq)\begin{aligned}T*\omega_{phy}&=\frac{3}{2}\omega_e(\psi_fi_q+(L_d-L_q)i_di_q)\\T&=\frac{3}{2}\frac{\omega_e}{\omega_{phy}}(\psi_fi_q+(L_d-L_q)i_di_q)\\&=\frac{3}{2}p(\psi_fi_q+(L_d-L_q)i_di_q)\end{aligned} T∗ωphyT=23ωe(ψfiq+(Ld−Lq)idiq)=23ωphyωe(ψfiq+(Ld−Lq)idiq)=23p(ψfiq+(Ld−Lq)idiq)
其中ppp是极对数。