MTPA based method for parameterless and position-sensorless control of a permanent magnet synchronous motor

Information

  • Patent Grant
  • 11689132
  • Patent Number
    11,689,132
  • Date Filed
    Thursday, March 3, 2022
    2 years ago
  • Date Issued
    Tuesday, June 27, 2023
    a year ago
Abstract
A parameterless and position-sensorless MTPA control of a permanent magnet synchronous motor including: using three rotating reference frames having different observation angles to parse the current vector; using a target current value and a preset current-rotor angle y that is between the current vector and the qr-axis of the (dr, qr) rotor reference frame to obtain the angles between the current vector, the voltage vector, and the rotor position; obtaining the target voltage value and the target voltage angle by using the obtained angles to obtain the target phase voltage values for regulation. The method is simple in controlling the motor, improves the control efficiency and reliability, and improves the control accuracy.
Description
BACKGROUND

The invention relates to a maximum torque per ampere (MTPA) based method for a position-sensorless control (which is carried out without a position sensor) of a permanent magnet synchronous motor.


Methods for a position-sensorless control of a permanent magnet synchronous motor typically include a constant-torque control mode, a constant-speed control mode, and a constant-air-volume control mode.


For example, U.S. Pat. No. 7,525,269 discloses an apparatus for a position-sensorless control of a permanent magnet synchronous motor drive system, in which a current-torque control is provided for performing a constant torque control.


Chinese patent publication No. 103929109A discloses a constant speed control method for a position-sensorless control of a permanent magnet synchronous motor.



FIG. 1 shows a block diagram of a conventional constant torque control. Referring to FIG. 1, the torque T is related only to the q-axis current, and thus a preset torque T can be converted to a preset q-axis current iq0 according to the formula T=K×iq0. Then a PI controller with the requisite corrective behavior can be used to perform a closed-loop control of the q-axis current, thus achieving the constant torque control.


The conventional field-oriented control (FOC) methods for a position-sensorless control of a permanent magnet synchronous motor are typically based on the rotor frame. The algorithm for deriving the rotor position in the FOC method highly depends on the accuracy of the motor parameters (e.g., the resistance Rs, the q-axis inductance Lq, the d-axis inductance Ld, and the magnetic flux λm), thus results in a large error due to for example the change in the rotor temperature. Also, the vector control algorithm in the FOC method is complicated, which requires time-consuming computation and large sources in the microcontroller unit (MCU). Therefore, a high-demanded microcontroller unit (MCU) is required to operates the FOC method, making the motor control costly. Moreover, because the FOC method is highly dependent on the accuracy of the motor parameters such as the resistance Rs, the q-axis inductance Lq, the d-axis inductance Ld, and the magnetic flux λm, the motor control based on the FOC method has a relative narrow range of applications.


SUMMARY

In view of the above-described problems, the disclosure provides a MTPA based control method for a parameterless and position-sensorless control of a permanent magnet synchronous motor. The method is operated without a position sensor (i.e., position-sensorless control) and does not depend on the motor parameters such as the resistance Rs, the q-axis inductance Lq, the d-axis inductance Ld, and the magnetic flux λm (i.e., parameterless control); its control algorithm is more uncomplicated and time-saving, and is more cost-saving because it does not require a demanded microcontroller unit (MCU).


To achieve the above objectives, in accordance with one embodiment of the invention, there is provided a MTPA based control method for a parameterless and position-sensorless control of a permanent magnet synchronous motor, the method comprising:


1). calculating a target dr-axis current value Idr and a target qr-axis current value Iqr of a target current vector {right arrow over (Idq)} in a (dr, qr) rotor reference frame by using a target current value Idq of the current vector {right arrow over (Idq)} and a preset current-rotor angle γ; in which the preset current-rotor angle y is measured between the current vector {right arrow over (Idq)} and the qr-axis of the (dr, qr) rotor reference frame; the (dr, qr) rotor reference frame is a first (d, q) rotating reference frame that is rotating in synchronism with the rotating magnetic field of the rotor, and the rotor position (i.e., the north polar of the magnetic field of the rotor) is aligned with the dr-axis of the (dr, qr) rotor reference frame; the target current value Idq is input by a user or obtained by using the real-time phase current values Ia, Ib, and Ic of the real-time current vector {right arrow over (Iabc)}; and the preset current-rotor angle γ is input by the user;


2). looking up a target voltage-rotor angle α or a target voltage-current angle β through a MTPA look-up table by referring to the target dr-axis current value Idr and the target qr-axis current value Iqr; in which the target voltage-rotor angle α is an angle between the dv-axis of a (dv, qv) voltage reference frame and the dr-axis of the (dr, qr) rotor reference frame; the (dv, qv) voltage reference frame is a second (d, q) rotating reference frame that is rotating in synchronism with the rotating magnetic field, and a target voltage vector {right arrow over (Vdq)} is aligned with the qv-axis of the (dv, qv) voltage reference frame such that a target dv-axis voltage value Vdv of the target voltage vector {right arrow over (Vdq)} equals to zero and a target qv-axis voltage value Vqv of the target voltage vector {right arrow over (Vdq)} equals to a target voltage value Vdq of the voltage vector {right arrow over (Vdq)} (i.e., Vdv=0 and Vqv=Vdq); the target voltage-current angle β is an angle between the voltage vector {right arrow over (Vdq)} and the current vector {right arrow over (Idq)}; and the MTPA look-up table is obtained in the maximum torque per ampere (MTPA) mode and comprises correspondences between the target dr-axis current value Idr, the target qr-axis current value Iqr, the preset current-rotor angle γ, the target voltage-rotor angle α, and the target voltage-current angle β;


3) calculating a PI error Δ by using the target dr-axis current value Idr, the target qr-axis current value Iqr, the target voltage-rotor angle a, the target voltage-current angle β, and the real-time phase current values Ia, Ib, and Ic; and obtaining the target voltage value Vdq of the voltage vector {right arrow over (Vdq)} by regulating the PI error Δ to be zero through a PI controller via the formula Vdq=PI (Δ);


4) obtaining a target voltage angle θv of the target voltage vector {right arrow over (Vdq)}, in which the target voltage angle θv is an angle between the target voltage vector {right arrow over (Vdq)} and the A-axis of the orthogonal (A, B) stationary reference frame; the A-axis of the orthogonal (A, B) stationary reference frame is aligned with the a-axis of the (a, b, c) stationary reference frame (i.e., the phase-a winding having the real-time phase current Ia); and the target voltage angle θv is obtained by using the target dr-axis current value Idr, the target qr-axis current value Iqr, the target voltage-rotor angle α, and the real-time phase current values Ia, Ib, and Ic, or obtained by using a given rotating speed value Spd of the rotor;


5) calculating the target A-axis voltage value VA and the target B-axis voltage value VB of the target voltage vector {right arrow over (Vdq)} in the orthogonal (A, B) stationary reference frame by using the target voltage value Vdq and the target voltage angle θv through the inverse Park transmission; converting the target A-axis voltage value VA and the target B-axis voltage value VB into target phase voltage values Va, Vb, and Vc, and modulating PWM signals of the inverter by using the phase the target voltage values Va, Vb, and Vc for regulating the real-time phase current values Ia, Ib, and Ic.


In a class of this embodiment, the method is operated under a current-control mode, in which:


in 1), the target current value Idq is input by a user; and the target dr-axis current value Idr and the target qr-axis current value Iqr are calculated via the formulas:

Idr=−Idq×sin(γ), and
Iqr=Idq×cos(γ);

in 3), the PI error Δ equals to the target dv-axis current value Idv of the target current vector {right arrow over (Idq)} minus the dv-axis current value idv_real of the real-time current vector {right arrow over (Iabc)} in the (dv, qv) voltage reference frame (i.e., Δ=Idv−Idr_real); in which the target dv-axis current value Idv and the target qv-axis current value Iqv of the target current vector {right arrow over (Idq)} are calculated via the formulas:

Idv=Idr*cos(α)+Iqr*sin(α), and
Iqv=−Idr*sin(α)+Iqr*cos(α);

the dv-axis current value Idv_real and the qv-axis current value Idv_real of the real-time current vector {right arrow over (Iabc)} are obtained by projecting the real-time current vector {right arrow over (Iabc)} from the orthogonal (A, B) stationary reference frame onto the (dv, qv) voltage reference frame via the Park transmission, by using the formulas:

Idv_real=IA*cos(θdv−qv)+IB*sin(θdv−qv), and


Iqv_real=IB*cos(θdv−qv)−IA*sin(θdv−qv); θdv−qv is the observation angle (i.e., the azimuth angle) for the (dv, qv) voltage reference frame, θdv−qv is an angle between the dv-axis of the (dv, qv) voltage reference frame and the A-axis of the orthogonal (A, B) stationary reference frame, and θvdv−qv+90°;


IA and IB are respectively the A-axis current value and the B-axis current value of the real-time current vector {right arrow over (Iabc)} in the orthogonal (A, B) stationary reference frame, and are converted from the real-time phase current values Ia, Ib, and Ic of the real-time current vector {right arrow over (Iabc)} in the (a, b, c) stationary reference frame via the Clarke transmission, by using the formulas:








I
A

=



2
3



I
a


-


1
3



(


I
b

+

I
c


)




,
and








I
B

=


2

3




(


I
b

-

I
c


)



;





in 4), the target voltage angle θv is obtained by inputting the target qv-axis current value Iqv of the target current vector {right arrow over (Idq)} and the qv-axis current value Iqv_real of the real-time current vector {right arrow over (Iabc)} in the (dv, qv) voltage reference frame into a phase lock loop (PLL) for processing; in which the phase lock loop is carried out via the formulas to decode θv:

Iqv=Iqv_real,
Iqv_real=IB*cos(θdv−qv)−IA*sin(θdv−qv), and
θvdv−qv+90°.


In a class of this embodiment, the method is operated under a speed-control mode, in which:


in 1), the target current value Idq is obtained via:


1a) constructing a (di, qi) current reference frame, in which the (di, qi) current reference frame is a third (d, q) rotating reference frame that is rotating in synchronism with the rotating magnetic field, and the real-time current vector {right arrow over (Iabc)} is aligned with the di-axis of a (di, qi) current reference frame; whereby a current angle θi_real of the real-time current vector {right arrow over (Iabc)} in the orthogonal (A, B) stationary reference frame is the same as the observation angle θdi−qi for the (di, qi) current reference frame (i.e., θi_realdi−qi), and a qi-axis current value Iqi_real of the real-time current vector {right arrow over (Iabc)} is zero; the current angle θi_real is an angle between the real-time current vector {right arrow over (Iabc)} and the A-axis of the orthogonal (A, B) stationary reference frame, and the observation angle θdi−qi is an angle between the di-axis of the (di, qi) current reference frame and the A-axis of the orthogonal (A, B) stationary reference frame;


1b) obtaining the di-axis current value Idi_real and the qi-axis current value Iqi_real of the real-time current vector {right arrow over (Iabc)} in the (di, qi) current reference frame by projecting the real-time current vector {right arrow over (Iabc)} in the orthogonal (A, B) stationary reference frame onto the (di, qi) current reference frame via the Park transmission, by using the formulas:

Idi_real=IA*cos(θi_real)+IB*sin(θi_real), and
Iqi_real=IB*cos(θi_real)−IA*sin(θi_real);


obtaining the current angle θi_real for the real-time current vector {right arrow over (Iabc)} by inputting the qi-axis current value Iqi_real and a zero into a phase lock loop (PLL); wherein the phase lock loop decodes θi_real via the formula

Iqi_real=IB*cos(θi_real)−IA*sin(θi_real)=0; and


1c) calculating the di-axis current value Idi_real by using the current angle θi_real, and using the di-axis current value Idi_real as the target current value Idq of a target current vector {right arrow over (Idq)} through Idq=Idi_real.

  • in 3), the PI error Δ equals to a real-time voltage-current angle θiv minus the target voltage-current angle β (i.e., Δ=θiv−β); wherein the real-time voltage-current angle θiv equals to the target voltage angle θv of the target voltage vector {right arrow over (Vdq)} minus the current angle θi_real of the real-time current vector {right arrow over (Iabc)} (i.e., θivv−θi_real); and
  • in 4), the target voltage angle θv is calculated by using a rotating speed value Spd of the motor that is input by a user through the formula θv=∫spd×(pole_pair×360×Δt÷60)·dt; in which pole pair is the number of the magnetic rotor poles, and Δt is the time variable.


In a class of this embodiment, in 2) the MTPA look-up table is obtained through experiments, theoretical calculations, or computer-aided finite-element-analysis software.


In a class of this embodiment, when the target voltage value Vdq (i.e., qv-axis voltage value Vqv) of the target voltage vector {right arrow over (Vdq)} is larger than or equal to the preset threshold Vmax, the PI controller is operated at a saturated state and the output voltage of the PI controller is limited to be the preset threshold Vmax, the dv-axis current value Idv of the target current vector {right arrow over (Idq)} is not useful for control, and the method is turns into a field-weakening control.


The benefits of the invention include:


1) The control method derives the rotor position without using a magnetic flux observer, thus requiring less CPU time to do the calculation, and being more intuitive and simpler to the position-sensorless control of the motor. The current-control mode and the speed-control mode of the disclosure are operated by two decoupled PI-controllers, achieving better control stability and dynamic response than the multi-stage nested control circuits.


2) In the current-control mode and the speed-control mode of the disclosure, the method regulates the current along the MTPA trajectory that can be calibrated to optimize the motor, and the motor is fully functional and allows starts with a full load. And the method of the invention is carried out over the full operating range from BEMF-free to field-weakening control.


3) The model used in the method of the invention is named a PLSL-MTPA mathematical model. The model does not rely on a single rotor reference frame, it converts the current vector of the motor to a current reference frame and a voltage reference frame, and parses the angles between the voltage vector, the current vector, and the rotor position, for performing the position-sensorless control. The mathematical model and the relative algorithms and calculations are simple in the method, which take less space on the chips, have low requirement for MCU, and make the motor more costly to control.


4) The PLSL-MTPA mathematical model used in the disclosure is an optimized position-sensorless control technology that is not dependent on the motor parameters. The technology solves the conventional bottleneck problem that the position-sensorless control of the motor is highly dependent on the motor parameters such as the resistance Rs, the q-axis inductance Lq, the d-axis inductance Ld, and the magnetic flux λm, such that the mathematical model has a relative wider range of applications.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a block diagram of a conventional position-sensorless FOC method for controlling a permanent magnet synchronous motor;



FIG. 2 is a perspective view of a permanent magnet synchronous motor, in accordance with one embodiment of the invention;



FIG. 3 is a perspective view of a motor controller of a permanent magnet synchronous motor, in accordance with one embodiment of the invention;



FIG. 4 is a sectional view of a permanent magnet synchronous motor, in accordance with one embodiment of the invention;



FIG. 5 is a diagram of a control logic in a motor controller of a permanent magnet synchronous motor, in accordance with one embodiment of the invention;



FIG. 6 is a circuit diagram corresponding to FIG. 5;



FIG. 7 is a schematic diagram of a (a, b, c) stationary reference frame for a three-phase permanent magnet synchronous motor, in accordance with one embodiment of the invention;



FIG. 8 is a schematic diagram of an orthogonal (A, B) stationary reference frame for a three-phase permanent magnet synchronous motor, in accordance with one embodiment of the invention;



FIG. 9 shows a relationship between a (a, b, c) stationary reference frame, an orthogonal (A, B) stationary reference frame, and a (d, q) rotating reference frame for a vector control of a three-phase permanent magnet synchronous motor, in accordance with one embodiment of the invention;



FIG. 10 shows a relationship between an orthogonal (A, B) stationary reference frame and a (d, q) rotating reference frame, in accordance with one embodiment of the invention;



FIG. 11 shows the changes of three angles θv, θi and θr having the same frequency along with the phase-a voltage component Va and the phase-a current la in a time domain, in accordance with one embodiment of the invention;



FIG. 12 shows a voltage vector {right arrow over (Vabc)}, a current vector {right arrow over (Iabc)}, and a rotor position in a (a, b, c) stationary reference frame, in accordance with one embodiment of the invention;



FIG. 13 is a schematic diagram of the control of the voltage vector {right arrow over (Vdq)} and the current vector {right arrow over (Idq)} in the method, in accordance with one embodiment of the invention;



FIG. 14 shows the dr-axis current value Idr and the qr-axis current value Iqr in a (dr, qr) rotor reference frame and the dv-axis current value Idv and the qv-axis current value Iqv in a (dv, qv) voltage reference frame for the current vector {right arrow over (Idq)} in accordance with one embodiment of the invention;



FIG. 15 is a block diagram of the method under the current-control mode, in accordance with one embodiment of the invention;



FIG. 16 is a block diagram of the method under the speed-control mode, in accordance with one embodiment of the invention;



FIG. 17 is a correlation between the voltage-rotor angle α, the dr-axis current value Idr, and the qr-axis current value Iqr, in accordance with one embodiment of the invention;



FIG. 18 is a correlation between the voltage-current angle β, the dr-axis current value Idr, and the qr-axis current value Iqr, in accordance with one embodiment of the invention.





DETAILED DESCRIPTION OF THE EMBODIMENTS

Detailed description of the invention will be given below in conjunction with the drawings.


Referring to FIGS. 2-4, the invention is used for controlling a three-phase permanent magnet synchronous motor that comprises a motor unit 1 and a motor controller 2. The motor unit 1 comprises a stator assembly 12, a rotor assembly 13 and a housing assembly 11. The stator assembly 12 is mounted on the housing assembly 11; the rotor assembly 13 is mounted within or outside the stator assembly 12; the motor controller 2 comprises a controller housing 22, and a control circuit board 21 mounted in the controller housing 22; the control circuit board 21 typically comprises a power supply circuit, a microprocessor, a bus voltage detection circuit, and an inverter; the power supply circuit supplies power to each part of the circuit. The DC bus voltage VDC-bus is applied from the bus voltage detection circuit to the microprocessor; the microprocessor controls the inverter to turn on/off the currents to the phase windings of the stator assembly 12.


Referring to FIGS. 5 and 6, a three-phase permanent magnet brushless DC motor inputs the phase currents Ia, Ib, and Ic through the phase current detection circuit to the microprocessor. An alternating-current input (AC input) flows through a full-wave rectification circuit comprising four diodes D7, D8, D9, and D10, and a DC bus voltage Vdc_bus is output at one end of a capacitor C1, wherein the DC bus voltage Vdc_bus is related to the AC-input voltage. The microprocessor inputs pulse-width modulation (PWM) signals to the inverter comprising six electronic switching tubes Q1, Q2, Q3, Q4, Q5, Q6, and the six control ends of the electronic switching tubes are separately controlled by six PWM signals (P1, P2, P3, P4, P5, P6) transmitted from the microprocessor.


Referring to FIG. 7, the phase currents of phase-a winding, phase-b winding, and the phase-c winding in the motor are labeled as Ia, lb and Ic, respectively, with a phase difference of 120 degrees between the two adjacent windings in the time domain, and a reference frame formed by respectively aligning an axis a, an axis b, and an axis c along the phase-a current Ia, the phase-b current Ib, the phase-c current Ic, is called a (a, b, c) stationary reference frame. Referring to FIG. 8, the three time-domain current quantities Ia, Ib and Ic, can be simplified into two currents IA and IB in an orthogonal (A, B) stationary reference frame. FIG. 9 is a relationship between a (a, b, c) stationary reference frame, an orthogonal (A, B) stationary reference frame, and a (d, q) rotating reference frame for a vector control of the invention. The mathematical relationships between the currents Ia, Ib, Ic in the (a, b, c) stationary reference frame and the currents IA and IB in the orthogonal (A, B) stationary reference frame of the same current vector, are as follows:












I
A

=



2
3



I
a


-


1
3



(


I
b

+

I
c


)




,
and









I
B

=


2

3




(


I
b

-

I
c


)



;





(
1
)







The above mathematical transformation in the formulas (1) is the so-called Clarke Transformation, where IA and IB are two sinusoidal currents that change with time. When IA and IB are observed on a (d, q) rotating reference frame having the same angular frequency as the sinusoidal currents IA and IB, the sinusoidal characteristics are eliminated from the current values and only the phase characteristics are retained in the current values. The sinusoidal current values IA and IB of a current vector in the orthogonal (A, B) stationary reference frame are converted to the direct-current (DC) d-axis and q-axis current values Id and Iq in a (d, q) rotating reference frame by performing the Park Transformation through the formulas:

Id=IA*cos(θd−q)+IB*sin(θd−q), and
Iq=IB*cos(θd−q)−IA*sin(θd−q).   (2)


Referring to FIG. 10, the variable θd−q in the PARK transformation of the formulas (2) is the angle between the d-axis of the (d, q) rotating reference frame and the A-axis of the orthogonal (A, B) stationary reference frame, and is called the azimuth angle or observation angle for the (d, q) rotating reference frame. When observing the current values IA and IB in different (d, q) rotating reference frames having different azimuth angles θd−q, the current values IA and IB are projected onto the different (d, q) rotating reference frames to obtain different current values Id and Iq. When a (d, q) rotating reference frame having an observation angle θd−q is selected for performing the Park transformation, the sinusoidal current values IA and IB are converted to the DC current values Id and Iq. Other sinusoidal A-axis and B-axis variable quantities (e.g., the voltage values of a voltage vector) in the orthogonal (A, B) stationary reference frame can also be converted to direct-current d-axis and q-axis variable quantities in a (d, q) rotating reference frame by using Park transformation through the formulas in (2), is the observation angle θd−q of the (d, q) rotating reference frame is given. In the position-sensorless control of a permanent magnet synchronous motor, the engineers are faced with an important problem of how to select a suitable observation angle θd−q for converting the sinusoidal A-axis and B-axis variable quantities (e.g., the current values of a current vector, and the voltage values of a voltage vector) in the orthogonal (A, B) stationary reference frame to the direct-current d-axis and q-axis variable quantities in a (d, q) rotating reference frame, for the purpose of convenient control of the current and voltage vectors of the motor. In a conventional motor control method, the rotor position that refers to the north polar of the magnetic field of the rotor is originated at 0 degree, and the rotor goes through a full circle is defined as a 360° of rotation. The relationship between the rotor position and the number of the magnetic rotor poles is θr=Pole_pair×θ0, in which θr is the rotor angle representing the angle of the rotor position in the orthogonal (A, B) stationary reference frame with respect to the A-axis of the orthogonal (A, B) stationary reference frame, and is measured between the rotor position and the A-axis of the orthogonal (A, B) stationary reference frame; θ0 is the mechanical angle of the rotor; and pole_pair is the number of the magnetic rotor poles.


In most position-sensorless control of the permanent magnet synchronous motors, the mechanical angle θ0 of the rotor is unknown, the task of the control algorithm is to estimate the angle θr that is the key element to execute the control algorithm of the motor. But the conventional FOC method to estimate θr used in the control is highly dependent on the motor parameters, thus its control algorithm and mathematical model are complicated, which requires time-consuming computation.


The disclosure provides a MTPA based method for parameterless and position-sensorless control of a permanent magnet synchronous motor, which is named a parameterless sensorless MTPA (PLSL-MTPA) method. It is important to note that the sinusoidal quantities are parsed by establishing different (d, q) rotating reference frames having different observation angles. In particular, the current vector and the voltage vectors are projected onto a (dr, qr) rotor reference frame, a (dv, qv) voltage reference frame, and a (di, qi) current reference frame which all are (d, q) rotating reference frames having different observation angles θd−q. In particular, the rotor position is aligned with the dr-axis of the (dr, qr) rotor reference frame such that the rotor angle θr equals to the observation angles θdr−qr of the (dr, qr) rotor reference frame. The voltage vector {right arrow over (Vdq)} is aligned with the qv-axis of the (dv, qv) voltage reference frame, such that the voltage angle θv equals to the observation angles θdv−qv of the (dv, qv) voltage reference frame plus 90° (i.e., θvdv−qv+90°); the dv-axis value Vdv of the voltage vector {right arrow over (Vdq)} is zero (i.e., Vdv=0); and qv-axis value Vqv of the voltage vector {right arrow over (Vdq)} equals to the voltage value Vdq of the voltage vector {right arrow over (Vdq)} (i.e., Vqv=Vdq). The current vector {right arrow over (Idq)} is aligned with the di-axis of the (di, qi) current reference frame, such that the current angle θi is the same as the observation angles θdi−qi of the (di, qi) current reference frame (i.e., θidi−qi); the qi-axis value Iqi of the current vector {right arrow over (Idq)} is zero (i.e., Iqi=0); and di-axis value Idi of the current vector {right arrow over (Idq)} equals to the current value Idq of the current vector {right arrow over (Idq)} (i.e., Idi=Idq).


With reference to FIG. 12, the voltage angle θv represents the angle of the voltage vector {right arrow over (Vdq)} in the orthogonal (A, B) stationary reference frame with respect to the A-axis of the orthogonal (A, B) stationary reference frame, and is measured between the voltage vector {right arrow over (Vdq)} and the A-axis of the orthogonal (A, B) stationary reference frame. The current angle θi represents the angle of the current vector {right arrow over (Idq)} in the orthogonal (A, B) stationary reference frame with respect to the A-axis of the orthogonal (A, B) stationary reference frame, and is measured between the current vector {right arrow over (Idq)} and the A-axis of the orthogonal (A, B) stationary reference frame.


In this invention, the current and voltage vectors are parsed in a (dv, qv) voltage reference frame and a (di, qi) current reference frame as claimed through the Park transformation, and a phase lock loop (PLL) is introduced to decode the observation angle θdv−qv of the (dv, qv) voltage reference frame or the observation angle θdi−qi of the (di, qi) current reference frame so as to arrive at the voltage angle θv or the current angle θi. The working principle of a phase lock loop is opposite to that of the Park Transformation. Particularly, the Park Transformation converts the sinusoidal variable quantities into the DC (direct-current) variable quantities by using the observation angle θd−q of a (d, q) rotating reference frame, while the phase lock loop decodes the observation angle θd−q of a (d, q) rotating reference frame by locking one DC variable quantity on an axis of the (d, q) rotating reference frame.


In the speed control mode of the method, the voltage angle θv is obtained by integrating the open-loop speed Spd for calculating the A-axis and B-axis voltage value VA and VB in the orthogonal (A, B) stationary reference frame and the phase voltage values Va, Vb and Vc in the (a, b, c) stationary reference frame. The phase lock loop outputs the current angle θi_real of the real-time current vector {right arrow over (Iabc)} that is generated by parsing the A-axis and B-axis current value IA and IB of the real-time current vector {right arrow over (Iabc)} in the orthogonal (A, B) stationary reference frame in a (di, qi) current reference frame by letting the qi-axis current value Iqi_real to be 0. In particular, the current angle θi_real is decoded by using the above formulas (2) to fulfill the relationship (Idi_real, 0)=Park transformation (IA, IB) by θi_real. In other words, θi_real is decoded in the PLL by fulfilling the relationship: 0=IB×cos(θi_real)−IA×sin(θi_real). The di-axis current value Idi_real is the calculated by using the formula Idi_real=IA×cos(θi_real)+IB×sin(θi_real) and is used as the target current value Idq of the target current vector {right arrow over (Idq)}.


In the synchronous motors, the angles θv, θi and θr are three time-domain variables with the same frequency. Taking the phase a for example, supposing the phase-a voltage component Va(t) has a phase lead β with respect to the phase-a current component Ia(t), the phase-a voltage component Va(t) and the phase-a current component Ia(t) are expressed as: Va(t)=Vabc×cos(θv)=Vabc×cos(ωt+β) and Ia(t)=Iabc×cos(θi)=Iabc×cos(ωt); in which Vabc is the voltage value of a voltage vector {right arrow over (Vabc)} composing the phase voltage components Va(t), Vb(t), Vc(t) in the (a, b, c) stationary reference frame; Iabc is the current value of a current vector {right arrow over (Iabc)} composing the phase current components Ia(t), Ib(t), Ic(t) in the (a, b, c) stationary reference frame; ω is the angular velocity of the rotor; t is the time; and β is the angle between the voltage component Va (t) and the current component Ia (t) of the phase a.


The relationship between the angles θv, θi and θr is as follows:

β=θv−θi   (3), and
θr=θv−α−90°  (4).



FIG. 11 shows the changes of three angles θv, θi and θr along with the phase-a voltage Va and the phase-a current Ia in a time domain. In one example, the PLSL-MTPA method of the invention control the result of θv−θi to result in a given voltage-current angle β according to the above formula (3), thus forming a speed-control mode. In another example, the PLSL-MTPA method of the invention regulate the projections of the sinusoidal currents IA and IB on a (dv, qv) voltage reference frame to result in a given voltage-rotor angle α according to the above formula (4), thus forming a current-control mode. In the PLSL-MTPA method of the invention, the dv-axis voltage value and the qv-axis voltage value of the voltage vector in the (dv, qv) voltage reference frame are controlled as: Vdv=0 and Vqv=Vdq, and the relationships between the voltage vector, the current vector, and the rotor position are controlled as that in FIG. 13. As shown in FIG. 13, the relation α=β+γ is obtained in the method of the invention.


Referring to FIGS. 14 and 15, in one example of the invention, the method is performed under a current-control mode, a (dv, qv) voltage reference frame in which the voltage vector is aligned with the qv-axis of the frame (i.e., Vdv=0 and Vqv=Vdq), is used to parse the orthogonal sinusoidal currents IA and IB of the real-time current vector {right arrow over (Iabc)}. The qv-axis projection Iqv_real of the real-time current vector {right arrow over (Iabc)} and the qv-axis projection Iqv of the target current vector {right arrow over (Idq)} in the (dv, qv) voltage reference frame are processed with the phase lock loop (PLL) to obtain the voltage angle θv of the target voltage vector {right arrow over (Vdq)}. And the dv-axis projection Idv_real of the real-time current vector {right arrow over (Iabc)} and the dv-axis projection Idv of the target current vector {right arrow over (Idq)} in the (dv, qv) voltage reference frame are processed with the PI controller to obtain the voltage value Vdq of the target voltage vector {right arrow over (Vdq)}. Then, the DC quantities Vdv and Vqv of the target voltage vector {right arrow over (Vdq)} in the (dv, qv) voltage reference frame are converted into the sinusoidal quantities VA and VB in the orthogonal (A, B) stationary reference frame by performing a dq to AB transformation (i.e., a transformation from a (d, q) rotating reference frame to an orthogonal (A, B) stationary reference frame) through the inverse Park transmission by using the voltage value Vdq and the voltage angle θ. Finally, the sinusoidal voltages VA and VB are converted to the phase voltage values Va, Vb, and Vc for modulating PWM signals of the inverter so as to regulate the real-time phase currents Ia, Ib, and Ic.


Referring to FIGS. 14 and 16, in another example of the invention, the method is performed under a speed-control mode, which makes the voltage vector rotate at a given speed Spd to form an open-loop control. The voltage angle θv of the voltage vector {right arrow over (Vdq)} is obtained by: θv=∫spd×(pole_pair×360×Δt÷60)·dt, in which Spd is the speed value, pole_pair is the number of the magnetic rotor poles, and Δt is the time variable. The phase lock loop (PLL) coincides the (real-time or target) current vector {right arrow over (Idq)} with the di-axis of a (di, qi) current reference frame (i.e., Iqi=0 and Idi=Idq) to generate the current angle θi of current vector {right arrow over (Idq)}; then the angle difference θiv is obtained according to θiv=θv−θi and the voltage value Vdq of the target current vector {right arrow over (Vdq)} is obtained according to Vdq=PI (θiv−β); and the DC quantities (Vdv=0, and Vqv=Vdq) in the (dv, qv) voltage reference frame are converted to the sinusoidal quantities voltage values VA and VB in the orthogonal (A, B) stationary reference frame by using the voltage angle θv and the voltage value Vdq.


The PI controller of the two control modes in the method of this invention regulates an angle error or a current error to produce the voltage value Vdq through Vdq=PI (Δ); in the speed-control mode, Δ=θiv−β; and in the current control mode, Δ=Idv−Idv_real. The phase current values Va, Vb, and Vc are derived from the voltage angle θv and the voltage value Vdq, and the voltage angle θv and the voltage value Vdq are calculated by different means in these two control modes.


A static full-load start-up process of the motor with the PLSL-MTPA control of the invention includes: the motor is started with the maximum current value Idq_max and is in a resistive state when the speed is low, the current vector and the voltage vector are in the same phase, and the PI controller tends to be operated at the maximum current value Idq_Max to drive the motor. Then the back electromotive force (BEMF) increases with the increasing rotational speed, causing the phase difference β between the voltage vector {right arrow over (Vdq)} and the current vector {right arrow over (Idq)} to gradually become non-zero, such that the PI controller is operated in a normal condition. The phase difference β between the voltage vector {right arrow over (Vdq)} and the current vector {right arrow over (Idq)} varies with the actual load, and the motor currents change accordingly. The start-up ability of the motor depends on the target speed, the speed increasing slope, current PI gain, and the maximum current limit.


The conventional FOC method is typically carried out through a rotor reference frame, and a position-sensorless algorithm is used to derive the rotor position for performing the FOC method. Compared with the conventional FOC method, the PLSL-MTPA control method of the invention is not dependent on the motor parameters, it employs an angle conversion scheme and regulates phase difference between the current vector and the voltage vector, for realizing the synchronous control of the motor. The PLSL-MTPA scheme of the invention greatly simplifies the processes of a sensorless control of a motor.


The PLSL-MTPA control method of the disclosure employs the following four control schemes:

  • 1) Under synchronous operation conditions, an adjustable angle β between the current vector and the voltage vector is unchanged, ensuring the comprehensive control of the synchronous motor, in which the angle β is also commonly referred to as the power factor angle;
  • 2) Under synchronous operation conditions, the rotor position of a synchronous motor is derivable by using the voltage vector and the adjustable angle α, in which the angle α may lead or lag the voltage vector. For a given angle α, a control of a synchronous motor is fulfilled by applying a voltage vector to the motor;
  • 3) a control under the maximum torque per current (MTPA) mode is carried out by converting the target current into the angles α and β, and the two angles α and β are determined according to the MTPA principle to be used the schemes in 1) and 2);
  • 4) only when the target current is controlled along the MTPA trajectory and is converted to the angles α and β under the MTPA criterion, the rotor position derived through the scheme in 2) is the actual rotor position of the motor, and the torque generated by the motor is called MTPA torque; and
  • 5) When the motor runs at a small current, the ratio of voltage to speed can be approximated as back electromotive force (BEMF). The operating conditions of a synchronous motor can be determined by comparing the ratio with a threshold.


The PLSL-MTPA method under the speed-control mode of the invention is an open-loop speed control that can be carried out under a load and through the field-weakening area. The speed and the position of the rotor only change with control commands, so that the control method is simple and optimizable and has a wide range of applications on electric machinery and electric drives.


The MTPA look-up table, also known as MTPA_Angle_Lookup, is mainly used to look up the angles α and β angle in FIG. 13. The relationship between the angles α and β and the MTPA look-up table is as follows:







(



α




β



)

=

MTPA_Angle

_Lookup







(


I
dr

,

I
qr


)

.






The angles α and β angle are determined by looking up the MTPA look-up table by referring the current Idr and the current Iqr.


The MTPA look-up table is a data set obtained from experiments. For example, a motor of 1/3 HP is used for experiments and its power is measured a dynamometer. A test for determining the included angle γ includes the steps of: the rotational speed of the motor is set at 1450 rpm, and the torque interval of the dynamometer is set at 10, 15, 20, 26, 31 oz-in, respectively. A maximum system efficiency or a maximum MTPA Index is searched under each torque range of the dynamometer, the corresponding current values Idr, Iqr and the angles α, β, and γ are recorded. Next, the test results are subjected to verification, five sets of data are obtained under the rotational speed of 1450 rpm, the set of data that allows for the maximum torque is selected and written to the MTPA look-up table. As shown in Table 1, the set of data including the angle γ, the current Idr, the current Iqr, the angle α, and the angle β that allows for the maximum torque 31.64 are selected and written to the MTPA look-up table as the data correspondence at rotational speed of 1450 rpm. Similarly, the data set of the angle γ, the current Idr, the current Iqr, the angle α, and the angle β that allows for the maximum torque is determined for a different rotational speed of 1400 rpm, 1350 rpm, 1300 rpm, etc., and written to the MTPA look-up table. The data in Table 1 are obtained in the speed control mode.















TABLE 1





Sdq
Torque
γ
Idr
Iqr
α
β





















1450 rpm
10.46
1728
500
1652
3391
13356


1450 rpm
15.98
2640
1240
2344
6592
11308


1450 rpm
20.94
3500
1900
2964
8640
10456


1450 rpm
26.46
4483
2608
3672
10177
9988


1450 rpm
31.64
5498
3356
4388
11328
9604









For the current control mode, the MTPA look-up table is obtained in a similar way from experiments. For example, a motor of 1/3 HP is used for experiments and its power is measured a dynamometer. A test for determining the included angle γ includes the steps of: the current value Idq flowing through the motor is set to 3.2A, and the torque interval of the dynamometer is set at 12, 18, 24, 28, 33 oz-in, respectively. A maximum system efficiency or a maximum MTPA_Index is searched under each torque range of the dynamometer, and the corresponding current values Idr, Iqr and the angles α, β, and γ are recorded. The test results are subjected to verification for the current value Idq that defaults to 3.2 A, five sets of data are obtained under the current value Idq of 3.2 A, the set of data as shown in FIG. 2 that allows for the maximum torque is selected and written to the MTPA look-up table.















TABLE 2





Idq
Torque
γ
Idr
Iqr
α
β





















3.2 A
12
1455
650
2253
2351
13450


3.2 A
18
2228
1440
2698
4323
12011


3.2 A
24
3400
1800
3298
6543
10244


3.2 A
28
4578
2708
3789
8678
951


3.2 A
33
5677
3676
4565
10233
8867









The MTPA look-up table can also be obtained from computer-aided finite element analysis software. As shown in FIGS. 17 and 18, the relevant data is obtained by computer-aided finite element analysis, and, the angles α or β are obtained under various current Idr and current Iqr. The MTPA look-up table can also be obtained from theoretical calculations.


The control method derives the rotor position without using a magnetic flux observer, thus requiring less CPU time to do the calculation, and being more intuitive and simpler to the position-sensorless control of the motor. The current control mode and the speed control mode of the motor are operated by two decoupled PI-controllers, achieving better control stability and dynamic response than the multi-stage nested control circuits. In the current control mode and the speed control mode of the disclosure, the control method regulates the current to flow along the MTPA trajectory that can be calibrated. The motor is fully functional and allows for starting up with a full load, over the full operating range from BEMF-free to field-weakening control. The PLSL-MTPA mathematical model of the disclosure does not reply on a single rotor reference frame, it converts the current vector of the motor to a current reference frame and a voltage reference frame, thus parsing the included angle of the vectors to perform the position-sensorless control. The mathematical models, the relative algorithms, and the calculations in the method are simple, and thus the method requires less space on the chips and low requirement for MCU and makes the motor more costly to control. The PLSL-MTPA mathematical model used in the disclosure is an optimized position-sensorless control technology that is not dependent on the motor parameters. The technology solves the bottleneck problem that the position-sensorless control of the motor is highly dependent on the motor parameters such as the resistance Rs, the d-axis inductance Ld, the q-axis inductance Lq, and the magnetic flux λm, such that the mathematical model has a relative wider range of applications.


Comparison of the field-weakening control method of the disclosure and that of the conventional FOC control theory:


The principle of the field-weakening control method of the conventional FOC control theory is that the control process of a permanent magnet synchronous motors is typically performed in two areas including MTPA area and a field weakening area. As the speed of the motor increases, there may be a maximum torque area or a maximum current area outside the field weakening area, but the two maximum areas are rarely used in practice. The FOC control method of the motor is a method to regulate the current values Id and Iq, that is, a control method with two degrees of freedom. The MTPA control is a theory that controls the two current values Id and Iq to keep a motor in optimized operational condition. When the current enters the field-weakening area, the freedom in the Id direction is locked, and only the current value Iq is proportional to the output torque, so that the motor may no longer be optimized. The relationship between the current value Id and field weakening is: since a permanent magnet is embedded in the rotor, the rotor magnetic field induces a back electromotive force (BEMF) that is offset by the stator voltage when the motor is rotating, the back electromotive force (BEMF) is proportional to the speed. Once the rotational speed exceeds a threshold, the BEMF will also exceed the stator voltage, causing the motor to fail to operate satisfactorily in the electric state. The so-called field-weakening control is a method to increase the current value Id in the negative direction when the back electromotive force (BEMF) exceeds a threshold, thereby generating a magnetic field to specifically weaken the rotor magnetic field. The current value Id depends on the rotational speed and the motor load, but the ultimate goal is to make the back electromotive force (BEMF) less than the maximum stator voltage.


In the classic field-weakening control theory, the following motor electromagnetic equations are used, in which Vd is a d-axis voltage of the voltage vector; Vq is the q-axis voltage; Id is the d-axis current of the current vector; Iq is the q-axis current; r is the motor resistance; O)r is the angular speed of the rotor; λd is the d-axis flux linkage; λq is the q-axis flux linkage; λpm is the permanent magnetic flux linkage; Ld is the d-axis inductance; and Lq is the q-axis inductance; the quantities Vd, Vq, Id, Iq, λd, λq, Ld, and Lq are components in a (d, q) rotating reference frame.

Vd=rId−ωrλq=rId−ωrLqIq;
Vq=rIqrλd=rIqr(LdIdpn).


Because the motor resistance r is very small, the components rId and rIq are negligible in the above electromagnetic equations. And due to the voltage limitation, the current values Id and Iq satisfy the following elliptic equation:








(


V
0



ω
r



L
q



)

2

=


I
q
2

+



(


L
d


L
q


)

2





(


I
d

+


λ
pm


L
d



)

2

.







Under the field-weakening condition with








V
0

=


2


V
DC


π


,





and the single solution of the current value Iq is obtained by the following formula, and the current value Iq can be calculated according to the known current Id







I
q

=





(


V
0



ω
r



L
q



)

2

-



(


L
d


L
q


)

2




(


I
d

+


λ
pm


L
d



)

2




.





It can be seen from the actual test simulated result that, the currents Id and Iq can flow along the MTPA trajectory at low speed, and satisfy the voltage-ellipse limit at high speed. The solution of the current value Id is given by the control logic (such as a PI controller, a look-up table and so on) instead of a given formula. In summary, the key to the field-weakening control is to find the current value Id that ensures the voltage does not exceed the preset threshold, and then the current value Iq is calculated by using the Id, thus achieving the current control of a motor under the FOC theory.


The control method of the disclosure, also known as the PLPS-MTPA control method, has different control means in the field-weakening control compared to the conventional FOC control. The field-weakening control method of the invention directly regulates the voltage instead of finding the current value Id in the conventional method. The voltage value Vdq (Vdq=Vq, Vd=0) is limited below the DC bus voltage Vdc_bus, the resulting current value Id is completely snubbed by the inverter, and the phase angle of the current and voltage vectors are regulated to generate the required torque. The control strategy in the invention completely meets the field-weakening theory, its implementation process is greatly simplified, and it allows for a more reliable control of the motor.


EXAMPLE 1

As shown in FIGS. 14 and 15, the PLSL-MTPA control method of the disclosure is carried out under a current-control mode as follows:


a), a target current value Idq and a preset current-rotor angle γ are input by a user; second, the target dr-axis current value Idr and qr-axis current value Iqr are calculated by projecting target current value Idq onto the dr-axis and the qr-axis of the (dr, qr) rotor reference frame; the target current value Idq and a preset current-rotor angle γ are data conforming to the MTPA mode.


b), a target voltage-rotor angle a is looked up by using a MTPA look-up table according to the target dr-axis current value Idr and qr-axis current value Iqr, in which the angle α is an included angle between the (dr, qr) rotor reference frame and the (dv, qv) voltage reference frame, the MTPA look-up table refers to a set of data obtained in the maximum torque per ampere (MTPA) mode.


c), the target dv-axis current value Idv and the target qv-axis current value Iqv of the target current vector {right arrow over (Idq)} are calculated by using the angle α, the dr-axis current value Idr, and the qr-axis current value Iqr, in which the dv-axis current value idv and the qv-axis current value Iqv are obtained by projecting the target current value Idq onto the (dv, qv) voltage reference frame;


d), the target qv-axis current value Iqv and the feedback qv-axis current value Iqv_real reflecting the operating conditions of the rotor in real time are inputted to a phase lock loop (PLL) to obtain the target voltage angle θv; target dv-axis current value Idv and the feedback dv-axis current value Idv_real are regulated with a PI controller to obtain the target voltage value Vqv. Because Vdv=0 and Vdq=Vqv, the voltage VA and the voltage VB in the orthogonal (A, B) stationary reference frame are calculated according to the angle θv and the target voltage value Vdq, and are used for controlling the phase currents of the motor.


The following formulas are used to calculate the dv-axis current Idv and the qv-axis current Iqv:

Idv=Idr*cos(α)+Iqr*sin(α), and
Iqv=−Idr*sin(α)+Iqr*cos(α).


The following formulas are used to calculate the dr-axis current value Idr and the qr-axis current value Iqr:

Idr=−Idq×sin(γ), and
Iqr=Idq×cos(γ).


In the current control mode, when the target current value Vdq (i.e., the qv-axis voltage value Vqv) is larger than or equal to the preset threshold Vmax, the PI controller is operated at a saturated state, the output voltage of the controller is limited to the preset threshold Vmax, and the dv-axis current value Idv cannot be used for control. Such condition is termed as a “field-weakening control.”


The MTPA look-up table is a data set obtained from experiments, theoretical calculations, or computer-aided finite element analysis softwares.


EXAMPLE 2

As shown in FIGS. 14 and 16, the PLSL-MTPA control method of the disclosure is carried out under a current-control mode as follows:


a), a given speed Spd and a preset angle γ are inputted by a user, and the vector angle of the voltage vector {right arrow over (Vdq)} is rotated at the given speed Spd to determine the angle θv, where the angle θv is an included angle between the voltage vector Vdq and the A-axis the (A, B) stationary coordinate system;


b), a (di, qi) current reference frame is used and the real-time current vector {right arrow over (Iabc)} is aligned with the di-axis of the (di, qi) current reference frame, such that the qi-axis current value θi_real in the frame is 0. A phase lock loop (PLL) is used to decode the current angle θi_real for the real-time current vector {right arrow over (Iabc)}. The di-axis current value Idi_real is calculated and used as the target current value Idq.


c), the target dr-axis current value Idr and the target qi-axis current value Iqr are calculated by using the angle γ and the target current value Idq, in which the target dr-axis current value Idr and the target qr-axis current value Iqr are obtained by converting the target current value Idq to the (dr, qr) rotor reference frame, and the current values Idr and Iqr are data conforming to the MTPA mode.


d), an target voltage-current β is determined by using the MTPA look-up table with reference to the current values Idr and Iqr, in which the angle β is an included angle between the voltage vector {right arrow over (Vdq)} and the current vector {right arrow over (Idq)} in the MTPA mode; and


e) an angle difference θiv between the target voltage angle θv and the current angle θi_real for the feedback real-time current vector Iabc is calculated according to θivv−θi_real. The target angle β and the angle difference θiv are regulated with the PI controller to obtain a target voltage value Vdq. Because Vdv=0 and Vdq=Vqv, the voltage VA and the voltage VB in the orthogonal (A, B) stationary reference frame are calculated according to the angle θv and the target voltage value Vdq, and are used for controlling the phase currents of the motor.


The MTPA look-up table is a data set obtained from experiments, theoretical calculations, or computer-aided finite element analysis softwares.


In the feedback qi-axis current Iqi_real reflecting the operating condition of the motor in real time is inputted to the phase lock loop (PLL), and another input current Iq* defaults to zero, in which the feedback qi-axis current Iqi_real is a projection of the real-time current vector {right arrow over (Iabc)} onto the qi-axis of the (di, qi) current reference frame. Then the phase lock loop (PLL) outputs the angle θi_real that is generated by parsing the currents IA and the IB with the phase lock loop (PLL) by fulfilling Iqi_real=0.


The following formula is used to calculate the angle θv:

θv=∫spd×(pole_pair×360×α60)·dt,

in which Spd is the speed value, pole_pair is the number of the magnetic rotor poles, and Δt is the time variable.


In the speed control mode, when the target current value Vdq (i.e., the qv-axis voltage value Vqv) is larger than or equal to the preset threshold Vmax, the PI controller is operated at a saturated state and the output voltage of the controller is limited to the preset threshold Vmax, and thus the current control mode is automatically converted to the filed-weakening control mode.


While particular embodiments of the invention have been shown and described, it will be obvious to those skilled in the art that changes and modifications may be made without departing from the invention in its broader aspects, and therefore, the aim in the appended claims is to cover all such changes and modifications as fall within the true spirit and scope of the invention.

Claims
  • 1. A maximum-torque-per-ampere (MTPA) based control method for a parameterless and position-sensorless control of a permanent magnet synchronous motor, the method comprising: 1) calculating a target dr-axis current value Idr and a target qr-axis current value Iqr of a target current vector {right arrow over (Idq)} in a (dr, qr) rotor reference frame by using a target current value Idq of the target current vector {right arrow over (Idq)} and a preset current-rotor angle γ; wherein the preset current-rotor angle γ is an angle between the target current vector {right arrow over (Idq)} and the qr-axis of the (dr, qr) rotor reference frame; the (dr, qr) rotor reference frame is a first (d, q) rotating reference frame that is rotating in synchronism with the rotating magnetic field, and the rotor position is aligned with the dr-axis of the (dr, qr) rotor reference frame; the target current value Idq is input by a user or obtained by using the real-time phase current values Ia, Ib, and Ic of a real-time current vector {right arrow over (Iabc)}; and the preset current-rotor angle γ is input by the user;2) looking up a target voltage-rotor angle α or a target voltage-current angle β through a MTPA look-up table by referring to the target dr-axis current value Idr and the target qr-axis current value Iqr; wherein the target voltage-rotor angle α is measured between the dv-axis of α (dv, qv) voltage reference frame and the dr-axis of the (dr, qr) rotor reference frame; the (dv, qv) voltage reference frame is a second (d, q) rotating reference frame that is rotating in synchronism with the rotating magnetic field, and a target voltage vector {right arrow over (Vdq)} is aligned with the qv-axis of the (dv, qv) voltage reference frame; the target voltage-current angle β is measured between the target voltage vector {right arrow over (Vdq)} and the target current vector {right arrow over (Idq)}; and the MTPA look-up table is obtained in the MTPA mode and comprises correspondences between the target dr-axis current value Idr, the target qr-axis current value Iqr, the preset current-rotor angle γ, the target voltage-rotor angle α, and the target voltage-current angle β;3) calculating a PI error Δ by using the target dr-axis current value Idr, the target qr-axis current value Iqr, the target voltage-rotor angle α, the target voltage-current angle β, and the real-time phase current values Ia, Ib, and Ic; and obtaining a target voltage value Vdq of the target voltage vector {right arrow over (Vdq)} by regulating the PI error Δ through a PI controller;4) obtaining a target voltage angle θv of the target voltage vector {right arrow over (Vdq)}, wherein the target voltage angle θv is measured between the target voltage vector {right arrow over (Vdq)} and the A-axis of the orthogonal (A, B) stationary reference frame; the target voltage angle θv is obtained by using the target dr-axis current value Idr, the target qr-axis current value Iqr, the target voltage-rotor angle α, and the real-time phase current values Ia, Ib, and Ic, or obtained by using a given rotating speed value Spd of the rotor; and5) calculating a target A-axis voltage value VA and a target B-axis voltage value VB in the orthogonal (A, B) stationary reference frame by using the target voltage value Vdq and the target voltage angle θv through the inverse Park transmission, converting the target A-axis voltage VA and the target B-axis voltage VB into target phase voltage values Va, Vb, and Vc, and modulating PWM signals of the inverter by using the target voltage values Va, Vb, and Vc, for regulating the real-time phase current values Ia, Ib, and Ic.
  • 2. The method of claim 1, the method is operated in a current-control mode, wherein: in 1), the target current value Idq is input by the user; and the target dr-axis current value Idr and the target qr-axis current value Iqr are calculated via the formulas: Idr=−Idq×sin(γ), and Iqr=Idq×cos(γ).in 3), the PI error Δ equals to a target dv-axis current value Idv of the target current vector {right arrow over (Idq)} minus a dv-axis current value Idv_real of the real-time current vector {right arrow over (Iabc)} in the (dv, qv) voltage reference frame by using Δ=Idv−Idv_real; wherein the target dv-axis current value Idv is calculated via the formula Idv=Idr*cos(α)+Iqr*sin(α); the dv-axis current value Idv_real is obtained via the park transmission by using an A-axis current value IA and a B-axis current value Is of the real-time current vector {right arrow over (Iabc)} in the orthogonal (A, B) stationary reference frame; the A-axis current value IA and the B-axis current value IB are converted by the real-time phase current values Ia, Ib, and Ic; and the dv-axis current value Idv_real is calculated by using the formula Idv_real=IA*cos(θdv−qv)+IB*sin(θdv−qv), θdv−qv is the observation angle of the (dv, qv) voltage reference frame, θdv−qv is measured between the dv-axis of the (dv, qv) voltage reference frame and the A-axis of the orthogonal (A, B) stationary reference frame, and θv=θdv−qv+90°; andin 4), the target voltage angle θv is obtained by inputting a target qv-axis current value Iqv of the target current vector {right arrow over (Idq)} and a qv-axis current value Iqv_real of the real-time current vector {right arrow over (Iabc)} in the (dv, qv) voltage reference frame into a phase lock loop; wherein the target qv-axis current value Iqv is calculated via the formula Iqv=−Idr*sin(α)+Iqr*cos(α); the qv-axis current value Iqv_real is calculated via the formula Iqv_real=IBcos(θdv−qv)−IA*sin(θdv−qv); and the phase lock loop decodes θv by fulfilling the formulas: Iqv=Iqv_real and θv=θdv−qv+90°.
  • 3. The method of claim 1, the method is operated in a speed-control mode, wherein: in 1), the target current value Idq is obtained via: constructing a (di, qi) current reference frame, wherein the (di, qi) current reference frame is a third (d, q) rotating reference frame that is rotating in synchronism with the rotating magnetic field, and the real-time current vector {right arrow over (Iabc)} is aligned with the di-axis of a (di, qi) current reference frame; whereby a current angle θi_real of the real-time current vector {right arrow over (Iabc)} is the same as the observation angle θdi−qi for the (di, qi) current reference frame, and a qi-axis current value Iqi_real of the real-time current vector {right arrow over (Iabc)} is zero; the current angle θi_real is an angle between the real-time current vector {right arrow over (Iabc)} and the A-axis of the orthogonal (A, B) stationary reference frame, and the observation angle θdi−qi is an angle between the di-axis of the (di, qi) current reference frame and the A-axis of the orthogonal (A, B) stationary reference frame;2b) obtaining a qi-axis current value Iqi_real of the real-time current vector {right arrow over (Iabc)} in the (di, qi) current reference frame via the Park transmission by using the formula Iqi_real=IB*cos(θi_real)−IA*sin(θi_real); and obtaining the current angle θi_real by inputting the qi-axis current value Iqi_real and a zero into a phase lock loop, wherein the phase lock loop decodes θi_real via the formula IB*cos(θi_real0−IA*sin(θi_real)=0;1c) obtaining a di-axis current value Idi_real of the real-time current vector {right arrow over (Iabc)} via the Park transmission by using the formula Idi_real=IA*cos(θi_real)+IB*sin(θi_real), and using the di-axis current value Idi_real as the target current value Idq;in 3), the PI error Δ equals to a real-time voltage-current angle θiv minus the target voltage-current angle β; wherein real-time voltage-current angle θiv equals to the target voltage angle θv minus the real-time current angle θi_real; andin 4), the target voltage angle θv is obtained by using the formula θv=∫spd×(pole_pair×360×Δt÷60)·dt; wherein Spd is a given rotating speed value that is input by the user, pole_pair is the number of the magnetic rotor poles, and Δt is the time variable.
  • 4. The method of claim 1, the MTPA look-up table is obtained through experiments, theoretical calculations, or computer-aided finite-element-analysis software.
  • 5. The method of claim 2, when the target voltage value Vdq is larger than or equal to the preset threshold Vmax, the dv-axis current value Idv of the target current vector {right arrow over (Idq)} is not useful for control, and the method is a field-weakening control.
  • 6. The method of claim 3, when the target voltage value Vdq is larger than or equal to the preset threshold Vmax, the dv-axis current value Idv of the target current vector {right arrow over (Idq)} is not useful for control, and the method is a field-weakening control.
Priority Claims (1)
Number Date Country Kind
201911028408.9 Oct 2019 CN national
CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation-in-part of International Patent Application No. PCT/CN2019/114919 with an international filing date of Nov. 1, 2019, designating the United States, now pending, and further claims foreign priority benefits to Chinese Patent Application No. 201911028408.9, filed Oct. 28, 2019. The contents of all of the aforementioned applications, including any intervening amendments thereto, are incorporated herein by reference. Inquiries from the public to applicants or assignees concerning this document or the related applications should be directed to: Matthias Scholl P. C., Attn.: Dr. Matthias Scholl Esq., 245 First Street, 18th Floor, Cambridge, Mass. 02142.

US Referenced Citations (5)
Number Name Date Kind
8339094 Perisic Dec 2012 B2
8618756 Hung Dec 2013 B2
20130169201 Kinoshita Jul 2013 A1
20140333241 Zhao Nov 2014 A1
20160056738 Kim Feb 2016 A1
Foreign Referenced Citations (8)
Number Date Country
101110559 Jan 2008 CN
102594250 Jul 2012 CN
103607155 Feb 2014 CN
103929109 Jul 2014 CN
109889117 Jun 2019 CN
110071674 Jul 2019 CN
2544362 Jan 2013 EP
2003081110 Mar 2003 JP
Non-Patent Literature Citations (1)
Entry
Zhenbo Li, et al., Research of Sensorless Air Compressor MTPA Control, Journal of Taiyuan University of Science and Technology, Aug. 2018, pp. 269-275, vol. 39, No. 4, China Academic Journal Electronic Publishing House, China.
Related Publications (1)
Number Date Country
20220190762 A1 Jun 2022 US
Continuation in Parts (1)
Number Date Country
Parent PCT/CN2019/114919 Nov 2019 US
Child 17686391 US