PERMANENT MAGNET SYNCHRONOUS MOTOR PARAMETER MEASUREMENT METHOD

Information

  • Patent Application
  • 20240007036
  • Publication Number
    20240007036
  • Date Filed
    December 18, 2020
    3 years ago
  • Date Published
    January 04, 2024
    4 months ago
Abstract
A permanent magnet synchronous motor parameter measurement method is provided. The method comprises: performing a maximum torque per ampere and a vector control on a permanent magnet synchronous motor, causing the permanent magnet synchronous motor to run stably and normally; when the permanent magnet synchronous motor is stably and normally running, an α-axis current and a β-axis current are obtained by means of three phase current sampling of the permanent magnet synchronous motor undergoing three phase to two phase transformation, and a rotor position angle and an electrical rotational speed of the permanent magnet synchronous motor are measured and obtained by means of reading a sensor on the permanent magnet synchronous motor; the six physical quantities mentioned-above are taken and a recursive least squares method is used to simultaneously obtain four parametric results for stator resistance, d-axis inductance, q-axis inductance, and flux linkage.
Description
BACKGROUND
Technical Field

The disclosure relates to a permanent magnet synchronous motor parameter measurement method in the field of permanent magnet synchronous motor measurement. In particular, the disclosure provides a parameter measurement method based on voltage equation and recursive least squares (RLS) method under the αβ coordinate system of permanent magnet synchronous motor.


DESCRIPTION OF RELATED ART

Online motor parameter measurement plays an important role in real-time motor control and fault diagnosis. Generally, the motor parameters that need to be identified online in permanent magnet synchronous motors include stator resistance Rs, d-axis inductance Ld, q-axis inductance Lq, and flux linkage ψf. At present, most of the literature on motor parameter measurement focuses on the motor voltage equation under the dq coordinate system. However, the voltage and current under the dq coordinate system are basically constant. Typically, only one parameter may be processed and obtained by using a d-axis or q-axis voltage equation, so only two motor parameters may be processed and obtained simultaneously in the steady state under the dq coordinate system.


Therefore, there is a problem of lack of rank in the parameter measurement of the permanent magnet synchronous motor under the dq coordinate system, and it is impossible to measure all parameters simultaneously. Some studies have shown that the method of signal input may increase the equation to measure all parameters simultaneously, but this approach will affect the normal operation of the motor.


SUMMARY

In order to solve the problems in the related art, the purpose of the present disclosure is to provide a parameter measurement method based on the recursive least squares method under the αβ coordinate system. Different from the dq coordinate system, two parameters may be obtained simultaneously by a sinusoidal AC phasor equation tinder the αβ coordinate system (on the condition that the phase difference between the two phasors in the equation is not 0 and 180 degrees). Therefore, four motor parameters may be obtained in the alp coordinate system. Compared with the full-parameter measurement under the dq coordinate system, the full-parameter measurement tinder the αβ coordinate system allows for faster convergence speed and shorter calculation time.


As shown in FIG. 1, the technical solution of the present disclosure is as follows:


(1) Maximum torque per ampere and vector control are performed on the permanent magnet synchronous motor, causing the permanent magnet synchronous motor to run stably and normally.


(2) When the permanent magnet synchronous motor is stably and normally running, the α-axis current iα and the β-axis current if are obtained through the three-phase current sampling of the permanent magnet synchronous motor undergoing three-phase to two-phase transformation. The rotor position angle θe, and electric rotational speed ωe of the permanent magnet synchronous motor are obtained by reading the sensor measurement on the permanent magnet synchronous motor.


(3) According to the α-axis current iα, β-axis current iβ, θe and ωe, obtained in step (2) and according to the preset input α-axis voltage vα and β-axis voltage vβ, the voltage equation under αβ coordinate system of permanent magnet synchronous motor is established. The six described physical quantities are taken and a recursive least squares method is used to simultaneously obtain four estimated values custom-character, custom-character, custom-character, and custom-character for four parameters including a stator resistance, a d-axis inductance, a q-axis inductance, and a flux linkage.


The specific process of the step (3) is as follows: the voltage equation under the αβ coordinate system is constructed in the following form:






v
α
−Ŵ
α2
X
α2

α1
X
α1






v
B
−Ŵ
α1
X
β2

α2
X
β1






Ŵ
α1=[custom-charactercustom-character]






Ŵ
α2=[custom-charactercustom-character]






X
α1
=[i
αpiβ−ωeiα)sin 2θepiα+(½piαeiβ)cos 2θe]T






X
α2=[½piα−(½piβ−ωeiα)sin 2θe−(½piαeiβ)cos 2θe−ωe sin θe]T






X
β1=[½piβ+(½piβ−ωeiα)cos 2θe−(½piαeiβ)sin 2θeωe cos θe]T






X
β2
=[i
βpiαeiβ)sin 2θepiβ−(½piβ−ωeiα)cos 2θe]T


In the equation: custom-character, custom-character, custom-character, and custom-character are respectively the estimated values of stator resistance, d-axis inductance, q-axis inductance and flux linkage; p represents the differential operator; Xα1 represents the first electrical input of the α-axis, Xα2 represents the second electrical input of the α-axis, Xβ1 represents the first electrical input of the β-axis, Xβ2 represents the second electrical input of the β-axis: Ŵα1 represents the estimated value between the stator resistance and the d-axis inductance, Ŵα2 represents the estimated value between the q-axis inductance Lq and the flux linkage ψf; and T represents the matrix transpose.


Then the recursive least squares method is used to solve the above equations simultaneously to obtain the estimated values of the four parameters, namely stator resistance Rs, d-axis inductance Ld, q-axis inductance Lq and flux linkage ψf.


An encoder is disposed on the rotor of the permanent magnet synchronous motor.


In the normal operation state of the permanent magnet synchronous motor, the voltage equation in the Δβ coordinate system of the permanent magnet synchronous motor and the recursive least square method are used for parameter measurement. The implementation steps of the method are as follows:


(1) Maximum torque per ampere (MTPA) and vector control are performed on the permanent magnet synchronous motor, causing the permanent magnet synchronous motor to run stably and normally.


(2) When the permanent magnet synchronous motor is stably and normally running, the α-axis current iα and the β-axis current iβ are obtained through the three-phase current sampling of the permanent magnet synchronous motor undergoing three-phase to two-phase transformation (abc/αβ transformation). The rotor position angle θe and electric rotational speed ωe of the permanent magnet synchronous motor are obtained by reading the sensor measurement on the permanent magnet synchronous motor.


An encoder is disposed on the rotor of the permanent magnet synchronous motor. The encoder is an absolute encoder.


(3) According to the α-axis current iα, β-axis current iβ, θe and ωe obtained in step (2) and according to the preset input α-axis voltage vα and β-axis voltage vβ, the voltage equation under αβ coordinate system of permanent magnet synchronous motor is established. The six described physical quantities are taken and a recursive least squares method is used to simultaneously obtain four parameters respectively, namely stator resistance, d-axis inductance, q-axis inductance, and flux linkage. In this way, simultaneous measurement is achieved.


The α-axis voltage vα and β-axis voltage vβ are given by the space vector pulse width modulation SVPWM input in the permanent magnet synchronous motor control.


The specific process of the full parameter measurement based on the recursive least squares method in the step (3) is as follows: the voltage equation under the αβ coordinate system is constructed in the following form:






v
α
−Ŵ
α2
X
α2

α1
X
α1






v
B
−Ŵ
α1
X
β2

α2
X
β1






Ŵ
α1=[custom-charactercustom-character]






Ŵ
α2=[custom-charactercustom-character]






X
α1
=[i
αpiβ−ωeiα)sin 2θepiα+(½piαeiβ)cos 2θe]T






X
α2=[½piα−(½piβ−ωeiα)sin 2θe−(½piαeiβ)cos 2θe−ωe sin θe]T






X
β1=[½piβ+(½piβ−ωeiα)cos 2θe−(½piαeiβ)sin 2θeωe cos θe]T






X
β2
=[i
βpiαeiβ)sin 2θepiβ−(½piβ−ωeiα)cos 2θe]T


In the equation: custom-character, custom-character, custom-character and custom-character are respectively the estimated values of stator resistance, d-axis inductance, q-axis inductance and flux linkage; p represents the differential operator; Xα1 represents the first electrical input of the α-axis, Xα2 represents the second electrical input of the α-axis, Xβ1, represents the first electrical input of the β-axis, Xβ2 represents the second electrical input of the β-axis; Ŵα1 represents the estimated value between the stator resistance and the d-axis inductance, Ŵα2 represents the estimated value between the q-axis inductance Lq and the flux linkage ψf; and T represents the matrix transpose.


Then the recursive least squares method is used to solve the above equations simultaneously to obtain the estimated values of the four parameters, namely stator resistance Rs, d-axis inductance Ld, q-axis inductance Lq and flux linkage ψf.


The measurement method of the present disclosure is different from the dq coordinate system. Under the αβ coordinate system, the voltage and current quantities are all sinusoidal AC phasors. When using linear regression strategies such as the recursive least square method for parameter measurement, two parameters may be obtained simultaneously by a sinusoidal AC phasor equation (provided that the phase difference between the two phasors in the equation is not 0 and 180 degrees). Therefore, four motor parameters may be obtained in the alp coordinate system. Compared with the full-parameter measurement under the dq coordinate system, the full-parameter measurement under the αβ coordinate system allows for faster convergence speed and shorter calculation time.


Advantageous effects of the present disclosure are as follows:


The disclosure combines the voltage equation set under the αβ coordinate system and the recursive least square method, and is able to directly perform full-parameter measurement on the permanent magnet synchronous motor. Compared with the full-parameter measurement under the dq coordinate system, the full-parameter measurement under the alp coordinate system of the present disclosure allows for faster convergence speed and shorter calculation time, and achieves higher instantaneity.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is an overall control block diagram of realizing a motor of the present disclosure.



FIG. 2 is a block diagram realizing the recursive least square method algorithm of the present disclosure.



FIG. 3 is a simulation verification diagram for realizing the algorithm of the present disclosure.





DESCRIPTION OF THE EMBODIMENTS

The present disclosure will be further elaborated below in conjunction with the accompanying drawings and embodiments.


According to the content of the present disclosure, the practical embodiment and the implementation of the complete embodiment is as follows.


In order to verify the reliability of this method, related experiments were carried out. The parameters of the interior permanent magnet motor (IPMSM) used as an example in the experiment are shown in Table I below.









TABLE 1





Motor Parameters


















Motor type
IPMSM











Stator resistance
0.065
Ω



Flux linkage
0.02
Vs



d-axis inductance
37.3
μH



q-axis inductance
48.8
μH



DC bus voltage
60
V



Rated torque
5
Nm



Rated speed
1500
rpm










The process flow is as follows.


(1) Maximum torque per ampere (MTPA) and vector control are performed on the permanent magnet synchronous motor, causing the permanent magnet synchronous motor to run stably and normally.


(2) An absolute encoder is provided on the rotor of the permanent magnet synchronous motor. When the permanent magnet synchronous motor is stably and normally running, the α-axis current iα and the β-axis current iβ are obtained through the three-phase current sampling of the permanent magnet synchronous motor undergoing three-phase to two-phase transformation (abc/αβ transformation). The rotor position angle θe and electric rotational speed ωe of the permanent magnet synchronous motor are obtained by reading the sensor measurement on the permanent magnet synchronous motor.


(3) According to the α-axis current iα, β-axis current iβ, θe and ωe obtained in step (2) and according to the preset input α-axis voltage vα and β-axis voltage vβ, the voltage equation under αβ coordinate system of permanent magnet synchronous motor is established. The six described physical quantities are taken and a recursive least squares method is used to simultaneously obtain four parameters respectively for stator resistance Rs, d-axis inductance Ld, q-axis inductance Lq, and flux linkage ψf. The voltage equation under the αβ coordinate system is constructed in the following form.






v
α
−Ŵ
α2
X
α2

α1
X
α1






v
B
−Ŵ
α1
X
β2

α2
X
β1






Ŵ
α1=[custom-charactercustom-character]






Ŵ
α2=[custom-charactercustom-character]






X
α1
=[i
αpiβ−ωeiα)sin 2θepiα+(½piαeiβ)cos 2θe]T






X
α2=[½piα−(½piβ−ωeiα)sin 2θe−(½piαeiβ)cos 2θe−ωe sin θe]T






X
β1=[½piβ+(½piβ−ωeiα)cos 2θe−(½piαeiβ)sin 2θeωe cos θe]T






X
β2
=[i
βpiαeiβ)sin 2θepiβ−(½piβ−ωeiα)cos 2θe]T


In the equation: custom-character, custom-character, custom-character, and custom-character are respectively the estimated values of stator resistance, d-axis inductance, q-axis inductance and flux linkage; p represents the differential operator: Xα1 represents the first electrical input of the α-axis, Xα2 represents the second electrical input of the α-axis, Xβ1 represents the first electrical input of the β-axis, Xβ2 represents the second electrical input of the β-axis; Ŵα1 represents the estimated value between the stator resistance and the d-axis inductance, Ŵα2 represents the estimated value between the q-axis inductance Lq and the flux linkage ψf; and T represents the matrix transpose.


Then the recursive least squares method is used to solve the above equations simultaneously to obtain the four parameters, namely stator resistance Rs, d-axis inductance Ld, q-axis inductance Lq and flux linkage ψf.


The specific block diagram of the full-parameter measurement method is shown in FIG. 2, and the solution process is as follows.


First, iα, iβ, θe, and ωe obtained by measurement and calculation were used to calculate Xα1, Xα2, Xβ1, and Xβ2. According to vα and Xα2 as well as Ŵα2β1, vα−Ŵα2Xα2 was calculated, and vα−Ŵα2Xα2 and Xα1 were used as the input of the recursive least squares method to obtain custom-character and custom-character; according to vβ and Xβ2 as well as Ŵβ2α1, vβ−Ŵβ2Xβ2 were calculated. vβ−Ŵβ2Xβ2 and Xβ1 were taken as the input of the recursive least squares method to obtain custom-character and custom-character.


Then the α-axis voltage equation is used to process the stator resistance custom-character and d-axis inductance custom-character and the β-axis voltage equation is used to process the q-axis inductance custom-character and flux linkage custom-character. Ŵα1=[custom-charactercustom-character] obtained by processing the α-axis voltage equation is used to calculate Ŵβ2 on the left side of the β-axis voltage equation. ωβ1=[custom-charactercustom-character] obtained by processing the β-axis voltage equation is used to calculate Ŵα2 on the left side of the α-axis voltage equation.


In this way, all parameters may be processed and obtained simultaneously.



FIG. 3 is the simulation result of the full-parameter measurement by the method of the present disclosure when the motor speed is 1500 rpm and the torque is 5 Nm. Based on the results, it can be seen that the method of the present disclosure is able to quickly obtain four motor parameters, and the error is very small while the convergence time is very short.

Claims
  • 1. A permanent magnet synchronous motor parameter measurement method, comprising: (1) performing a maximum torque per ampere and a vector control on a permanent magnet synchronous motor, causing the permanent magnet synchronous motor to run stably and normally;(2) when the permanent magnet synchronous motor is stably and normally running, obtaining an α-axis current iα and a β-axis current iβ through three-phase current sampling of the permanent magnet synchronous motor undergoing a three-phase to two-phase transformation, and obtaining a rotor position angle θe and an electric rotational speed ωe of the permanent magnet synchronous motor by reading a sensor measurement on the permanent magnet synchronous motor;(3) establishing a voltage equation under a αβ coordinate system of the permanent magnet synchronous motor according to the α-axis current iα, the β-axis current iβ, the θe and the ωe obtained in the step (2) and according to a preset input α-axis voltage vα and a preset input β-axis voltage vβ, wherein six physical quantities mentioned above are taken and a recursive least squares method is used to simultaneously obtain estimated values , , , for four parameters, wherein the four parameters are a stator resistance Rs, a d-axis inductance Ld, a q-axis inductance Lq, and a flux linkage ψf.
  • 2. The permanent magnet synchronous motor parameter measurement method according to claim 1, wherein a process of the step (3) is as follows: constructing the voltage equation under the αβ coordinate system in the following form: vα−Ŵα2Xα2=Ŵα1Xα1 vB−Ŵα1Xβ2=Ŵα2Xβ1 Ŵα1=[]Ŵα2=[]Xα1=[iα(½piβ−ωeiα)sin 2θe+½piα+(½piα+ωeiβ)cos 2θe]T Xα2=[½piα−(½piβ−ωeiα)sin 2θe−(½piα+ωeiβ)cos 2θe−ωe sin θe]T Xβ1=[½piβ+(½piβ−ωeiα)cos 2θe−(½piα+ωeiβ)sin 2θeωe cos θe]T Xβ2=[iβ(½piα+ωeiβ)sin 2θe+½piβ−(½piβ−ωeiα)cos 2θe]T wherein , , and are the estimated values of the stator resistance Rs, the d-axis inductance Ld, the q-axis inductance Lq, and the flux linkage ψf, respectively; p represents a differential operator; Xα1 represents a first electrical input of an α-axis, Xα2 represents a second electrical input of the α-axis, Xβ1 represents a first electrical input of a β-axis, Xβ2 represents a second electrical input of the β-axis; Ŵα1 represents an estimated value between the stator resistance and the d-axis inductance, Ŵα2 represents an estimated value between the q-axis inductance Lq and the flux linkage ψf; and T represents a matrix transpose;using the recursive least squares method to simultaneously solve above-mentioned equations to obtain the estimated values for the four parameters, wherein the four parameters are the stator resistance Rs, the d-axis inductance Ld, the q-axis inductance Lq, and the flux linkage ψf.
  • 3. The permanent magnet synchronous motor parameter measurement method according to claim 1, wherein an encoder is disposed on a rotor of the permanent magnet synchronous motor.
Priority Claims (1)
Number Date Country Kind
202011268499.6 Nov 2020 CN national
PCT Information
Filing Document Filing Date Country Kind
PCT/CN2020/137420 12/18/2020 WO