## Friday, September 8, 2017

### MTPA, MTPF, and Speed Range, Part 2

Previously on Motor Diddlers, we learned that interior PM motors spend most of their time under field weakening; namely, the maximum attainable torque is limited by inductance and available voltage, not necessarily by maximum allowable current. This post will cover the basic case of computing the MPTF trajectory in the absence of saturation and higher-order back EMF components.

Starting from the motor equations $$\begin{array}{lcl} \tau=\frac{3}{2}n_p(\lambda I_q+(L_d-L_q)I_d I_q)\\ V_d=R_s I_d-\omega L_q I_q\\ V_q=R_s I_q+\omega L_d I_d+\omega\lambda\\ V_s=\sqrt{V_d^2+V_q^2}\\ I_s=\sqrt{I_d^2+I_q^2} \end{array}$$ we note that in the voltage limited operating regime, maximum torque must be achieved when the voltage vector is on the boundary of the allowable area, in which case the optimization problem becomes an equality: $$\begin{cases} \frac{3}{2}n_p(\lambda I_q+(L_d-L_q)I_d I_q) = \tau_0\\ (R_s I_d-\omega L_q I_q)^2 + (R_s I_q+\omega L_d I_d+\omega\lambda)^2 = V_0^2 \end{cases}$$ If we ignore saturation, this system is the intersection of a constant-torque hyperbola (which is independent of speed) and a series of shrinking ellipses.

This system is polynomial (and in fact only a quartic) and can be solved in many ways. For example, one element of the reduced Groebner basis of the ideal generated by the two equations is the very long $$(L_d-L_q)^2(R_s^2+L_q^2w^2)I_q^4 + (L_q^2\lambda^2\omega^2+R_s^2\lambda^2+(L_d-L_q)^2(2R_s\tau_0\omega-V_0^2))I_q^2- 2\lambda\tau_0(R_s^2+L_d L_q\omega^2)+ R_s^2\tau_0^2+L_d^2\tau_0^2\omega^2$$ where we have included the factor of $$3/2n_p$$ in $$\tau_0$$ for brevity's sake. The roots of this polynomial are easily found with a variety of numeric or analytic methods; in the case of multiple real solutions the correct one is the $$(I_d,I_q)$$ vector with the shortest length.
Unfortunately, the system is no longer polynomial if the motor saturates (that is, $$L_d$$ and $$L_q$$ are functions $$l_d(I_d,I_q)$$, $$l_q(I_d,I_q)$$ of the axis currents). In order to solve the system in this case, we will have to turn to more complicated (and less reliable) numeric methods, the nature of which will be the subject of the next post in this series.