This application is a 371 of international application of PCT application serial no. PCT/CN2021/082326, filed on Mar. 23, 2021, which claims the priority benefit of China application no. 202110195977.3, filed on Feb. 22, 2021. The entirety of each of the above mentioned patent applications is hereby incorporated by reference herein and made a part of this specification.
The present invention relates to the field of control of a bearingless motor, and in particular, to compensation control for vibration of a bearingless permanent magnet synchronous motor by using dead-time compensation control and rotor eccentricity control technologies for the bearingless permanent magnet synchronous motor.
Bearingless permanent magnet synchronous motors are a new type of special motors that have high speed and high precision and require no lubrication. Their application prospects in aerospace, chemical manufacturing, semiconductor industry, and other fields in need of special environments increasingly grow. The bearingless permanent magnet synchronous motor is used as a rotary drive motor. Due to problems such as material unevenness, processing errors, and assembly errors, a certain degree of rotor mass eccentricity inevitably exists and a centrifugal excitation force at the same frequency as a rotation speed is produced when the motor is in rotation. Meanwhile, in a control process of the bearingless permanent magnet synchronous motor, a dead time must be set to avoid a short circuit between upper and lower bridge arms of an inverter. The introduction of the dead time causes an increase in current harmonics, and the amplitude of an unbalanced force is further increased, which results in unbalanced vibration of a rotor and affects suspension control precision of the rotor.
Regarding control for unbalanced vibration of the rotor in the bearingless permanent magnet synchronous motor, compensation control is mainly performed against unbalanced vibration caused by rotor mass eccentricity in the prior art, while unbalanced vibration caused by a dead-time effect is rarely concerned. Chinese Patent Publication No. CN104659990A discloses a method for extracting an unbalanced vibration displacement of a bearingless motor through adaptive filtering, which paves the way for the primary condition of vibration compensation control of a bearingless motor. Chinese Patent Publication No. CN105048913A discloses a current compensation-based unbalanced vibration control system for a bearingless asynchronous motor, wherein compensation control for suspension vibration is realized by adjusting a compensation current. However, in these solutions, the vibration compensation control of the bearingless motor mainly focuses on detection and compensation for vibration caused by eccentricity and does not concern vibration caused by the dead-time effect. To improve the control precision of the unbalanced vibration displacement of the bearingless permanent magnet synchronous motor, compensation not only needs to be made for a rotor eccentricity displacement caused by rotor mass eccentricity, but also needs to be made for unbalanced vibration of the rotor caused by the dead-time effect, which is critical to the implementation of high-precision control of the bearingless permanent magnet synchronous motor.
An objective of the present invention is to provide a vibration compensation controller with neural network band-pass filters for a bearingless permanent magnet synchronous motor. The controller performs vibration compensation to suppress vibration of the bearingless permanent magnet synchronous motor, thereby solving the problem in the prior art that compensation is only made for vibration caused by rotor mass eccentricity and vibration caused by a dead-time effect is ignored in the vibration compensation control of the bearingless permanent magnet synchronous motor. Therefore, stable suspension and efficient operation of a rotor of the motor are realized, the control precision of the motor is improved, and better application in an electric drive system is achieved.
The vibration compensation controller with neural network band-pass filters for a bearingless permanent magnet synchronous motor provided by the present invention adopts the following technical solution. The controller comprises a displacement controller and a rotating speed controller. The displacement controller includes a vibration force compensation control module and a dead-time vibration compensation module.
The vibration force compensation control module receives, as input, actual displacements x, y in x and y directions and a rotor mechanical angle θm and outputs corresponding vibration compensation forces Fxh and Fyh. The vibration force compensation control module comprises a first neural network band-pass filter, a second neural network band-pass filter, a third proportional-integral-derivative (PID) controller, and a fourth PID controller. The first neural network band-pass filter receives, as input, the actual displacement x in the x direction and the rotor mechanical angle θm and outputs a vibration displacement {circumflex over (x)}. A difference between a specified value 0 and the vibration displacement {circumflex over (x)} is input to the third PID controller, and the third PID controller outputs the vibration compensation force Fxh. The second neural network band-pass filter receives, as input, the actual displacement y in the y direction and the rotor mechanical angle θm and outputs a vibration displacement ŷ. A difference between the specified value 0 and the vibration displacement ŷ is input to the fourth PID controller, and the fourth PID controller outputs the vibration compensation force Fyh. A sum of the vibration compensation force Fxh and a specified force value Fx of a suspension winding in the x direction is input to a force/current conversion module, a sum of the vibration compensation force Fyh and a specified force value Fy of the suspension winding in the y direction is input to the force/current conversion module, and the current conversion module obtains a specified quadrature axis current value i*Bq and a specified direct axis current value i*Bd.
The dead-time vibration compensation module receives, as input, a rotor electrical angle θe, and an actual quadrature axis current iBq, and an actual direct axis current iBd and outputs a quadrature axis compensation voltage uBqh and a direct axis compensation voltage uBdh. The dead-time vibration compensation module comprises a third neural network band-pass filter in a direct axis direction, a fourth neural network band-pass filter in a quadrature axis direction, a sixth proportional-integral (PI) controller, and a seventh PI controller. The third neural network band-pass filter receives, as input, the actual current iBd in the direct axis direction and 6 times of the rotor electrical angle θe and obtains a harmonic current îBd in the direct axis direction. A difference between the specified value 0 and the harmonic current îBd is input to the sixth PI controller, and the sixth PI controller obtains the direct axis compensation voltage uBdh. A sum of a control voltage uBd in the direct axis direction and the direct axis compensation voltage uBdh serves as a direct axis command voltage u*Bd. The fourth neural network band-pass filter receives, as input, the actual current iBq in the quadrature axis direction and 6 times of the rotor electrical angle θe and obtains a harmonic current ÎBq in the direct axis direction. A difference between the specified value 0 and the harmonic current ÎBq is input to the seventh PI controller, and the seventh PI controller obtains the quadrature axis compensation voltage uBqh. A sum of a control voltage uBq in the quadrature axis direction and the quadrature axis compensation voltage uBqh serves as a quadrature axis command voltage u*Bq.
The present invention has the following beneficial effects:
1) By adopting dead-time vibration compensation control, the present invention not only compensates for the dead time, but also effectively suppresses vibration during the operation of the bearingless permanent magnet synchronous motor, thereby improving the suspension control precision.
2) The neural network band-pass filters adopted by the present invention have simple working principles and concise calculation processes and can obtain required signals according to real-time speeds of the motor.
3) The present invention adopts the PI controllers to regulate vibration. The controllers have simple principles, their coefficients can be adjusted conveniently, and they have strong robustness.
4) In the vibration compensation control of the bearingless permanent magnet synchronous motor, generally only the vibration caused by eccentricity is considered and compensation control is implemented, while the vibration caused by the dead-time effect is not concerned, which affects the entire suspension control precision. To achieve higher suspension control precision of the bearingless permanent magnet synchronous motor, the present invention not only performs analysis and compensation for vibration caused by eccentricity, but also performs compensation control for vibration caused by the dead-time effect, thereby effectively improving the suspension control precision.
In order to make the content of the present invention more obvious and understandable, the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
The specific ideas and implementation steps of the present invention are illustrated below.
Referring to
As shown in
wherein Ts is an interrupt cycle of the rotating speed controller 2 and Le is the number of lines of the coder.
A difference between the calculated actual speed n and a specified speed value n* serves as a speed error, and the error is input to the first PI controller 21. The first PI controller 21 makes adjustment to obtain a specified quadrature axis current value i*Mq of a torque winding. Meanwhile, a current sensor collects torque currents i2A and i2C of the two-phase torque winding of the bearingless permanent magnet synchronous motor 3, and inputs the torque currents i2A and i2C to the second coordinate transformation module 25. The second coordinate transformation module 25 is configured for performing Clarke transform and Park transform. The second coordinate transformation module 25 transforms i2A and i2C into an actual quadrature axis current value iMq of the torque winding and an actual direct axis current value iMd of the torque winding in a rotating reference frame. An error between the specified quadrature axis current value i*Mq the torque winding and the actual quadrature axis current value iMq of the torque winding is input to the second PI controller 22 to obtain a specified quadrature axis voltage value u*Mq of the torque winding. When a specified direct axis current value of the torque winding is i*Md=0, an error between i*Md and the actual direct axis current value iMd of the torque winding is input to the third PI controller 23 to obtain a specified direct axis voltage value u*Md of the torque winding. Output ends of the second PI controller 22 and the third PI controller 23 are both connected to an input end of the first coordinate transformation module 24. The first coordinate transformation module 24 is configured for performing inverse Park transform, through which the specified quadrature axis voltage value u*Mq of the torque winding and the specified direct axis voltage value u*Md of the torque winding can be transformed into voltages uMα and uMβ of the torque winding in a stationary reference frame. An output end of the first coordinate transformation module 24 is sequentially connected in series with the first SVPWM inverter 26 and the bearingless permanent magnet synchronous motor 3. The first coordinate transformation module 24 inputs the voltages uMα and uMβ to the first SVPWM inverter 26. An output of the first SVPWM inverter 26 is connected to an input of the bearingless permanent magnet synchronous motor 3. The first SVPWM inverter 26 obtains three-phase input voltages u2A, u2B, and u2C of the bearingless permanent magnet synchronous motor 3.
As shown in
The output end of the coder 27 is further connected to the angle calculation module 17. A pulse signal output by the coder 27 is input to the angle calculation module 17 to obtain a rotor mechanical angle θm. The rotor mechanical angle at a moment k is calculated as follows:
wherein ΔP is the accumulation result of pulses output by the coder 27.
Output ends of the angle calculation module 17 and the displacement calculation module 92 are both connected to an input end of the vibration force compensation control module 5. The rotor mechanical angle θm output by the angle calculation module 17 and the actual rotor displacements x, y output by the displacement calculation module 92 are input to the vibration force compensation control module 5 to obtain compensation forces Fxh and Fyh.
As shown in
{circumflex over (x)}(k)=ωx_1(k)·cos θm(k)+ωx_2(k)·sin θm(k) (3).
The weights ωx_1 and ωx_2 are calculated by the following formulas:
wherein ex is a component in the x direction after harmonics are filtered out; ωx_1 and ωx_2 are updated weights in the x direction; μ1 is a step factor.
Therefore, the vibration displacement {circumflex over (x)} in the x direction is obtained. As shown in
The second neural network band-pass filter 53 is identical to the first neural network band-pass filter 51 in structure and principle. Likewise, the displacement in the y direction and the rotor mechanical angle θm are input to the second neural network band-pass filter 53.
ŷ(k)=ωy_1(k)·cos θm(k)+ωy_2(k)·sin θm(k) (5).
The weights ωy_1 and ωy_2 are calculated by the following formulas:
wherein ey is a component in the y direction after harmonics are filtered out; ωy_1 and ωy_2 are updated weights in the y direction; μ1 is the step factor.
Therefore, the vibration displacement signal ŷ in the y direction is obtained. As shown in
A sum of the force Fx in the x direction output by the first PID controller 11 and the vibration compensation force Fxh in the x direction output by the vibration force compensation module 5 and a sum of the force Fy in the y direction output by the second PID controller 12 and the vibration compensation force Fyh in the y direction output by the vibration force compensation module 5 are input to the force/current conversion module 13 to obtain a specified quadrature axis current value i*Bq and a specified direct axis current value i*Bd of the suspension winding.
Differences between the obtained a specified quadrature axis current value i*Bq, and a specified direct axis current value i*Bd and an actual quadrature axis current value iBq, an actual direct axis current value iBd of the suspension winding are obtained respectively. The current sensor collects currents i1A and i1C of the two-phase suspension winding of the bearingless permanent magnet synchronous motor 3 and inputs the collected currents to the fourth coordinate transformation module 91. The fourth coordinate transformation module 91 is configured for performing Clarke transform and Park transform. The fourth coordinate transformation module 91 processes i1A and i1C to obtain the actual quadrature axis current iBq and the actual direct axis current iBd of the suspension winding. The difference between i*Bq and iBq is input to the fifth PI controller 15 to obtain a quadrature axis control voltage uBq of the suspension winding. The difference between i*Bd and iBd is input to the fourth PI controller 14 to obtain a direct axis control voltage uBd of the suspension winding.
A rotor electrical angle θe, the actual quadrature axis current value iBq of the suspension winding, and the actual direct axis current value iBd of the suspension winding are input to the dead-time vibration compensation module 6 to obtain compensation voltages uBqh and uBdh. The angle calculation module 17 processes the pulse signal, collected by the coder 27, of the bearingless permanent magnet synchronous motor 3 to obtain the rotor electrical angle θe, which is calculated as follows:
θe(k)=PMθm(k) (7)
wherein θm(k) is the rotor mechanical angle at the moment k according to the formula (2) and PM is the number of pole-pairs of the torque winding.
The obtained rotor electrical angle θe and the actual quadrature axis current value iBq and the actual direct axis current value iBd are input to the dead-time vibration compensation module 6. The dead-time vibration compensation module 6 comprises a third neural network band-pass filter 61 in the direct axis direction, a fourth neural network band-pass filter 63 in the quadrature axis direction, a sixth PI controller 62, and a seventh PI controller 64. In the dead-time vibration compensation module 6, compensations in the direct axis direction and the quadrature axis direction are shown in
îBd(k)=ωd6_1(k)·cos 6θe(k)+ωd6_2(k)·sin 6θe(k) (8).
The weights ωd6_1 and ωd6_2 are calculated by the following formulas:
wherein eBd is a component in the direct axis direction after harmonics are filtered out; ωd6_1 and ωd6_2 are updated sixth-harmonic weights in the direct axis direction; μ2 is a step factor.
Therefore, the harmonic current signal îBd in the direct axis direction is obtained. As shown in
The fourth neural network band-pass filter 63 in the quadrature axis direction is identical to the third neural network band-pass filter 61 in structure. Likewise, the current iBq in the quadrature axis direction and 6 times of the rotor electrical angle θe are input to the fourth neural network band-pass filter 63 in the quadrature axis direction to obtain a harmonic current signal îBq in the direct axis direction.
îBq(k)=ωq6_1(k)·cos 6θe(k)+ωq6_2(k)·sin 6θe(k) (10).
The weights ωd6_1 and ωd6_2 are calculated by the following formulas:
wherein eBq is a component in the quadrature axis direction after harmonics are filtered out; ωg6_1 and ωq6_2 are updated sixth-harmonic weights in the quadrature axis direction; μ2 is the step factor.
Therefore, the harmonic current signal îBq in the direct axis direction is obtained. As shown in
A sum of the direct axis voltage uBd output by the fourth PI controller 14 and the direct axis compensation voltage uBdh output by the dead-time vibration compensation module serves as a direct axis command voltage u*Bd. A sum of the quadrature axis voltage uBq output by the fifth PI controller 15 and the quadrature axis compensation voltage uBqh output by the dead-time vibration compensation module serves as a quadrature axis command voltage u*Bq. The obtained u*Bd and u*Bq are input to the third coordinate transformation module 16. The third coordinate transformation module 16 is configured for performing inverse Park transform. The third coordinate transformation module 16 processes u*Bd and u*Bq to obtain voltages uBα and uBβ of the suspension winding in the stationary reference frame.
The voltages uBα and uBβ of the suspension winding are input to the second SVPWM inverter 90. An output of the second SVPWM inverter 90 is connected to the input of the bearingless permanent magnet synchronous motor 3. The second SVPWM inverter 90 obtains three-phase input voltages u1A, u1B, u1C of the bearingless permanent magnet synchronous motor 3.
The present invention can be implemented based on the above descriptions. Other changes and modifications made by persons skilled in the art without departing from the spirit and protection scope of the present invention still fall within the protection scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
202110195977.3 | Feb 2021 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2021/082326 | 3/23/2021 | WO |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2022/174488 | 8/25/2022 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
20090116136 | Zhang | May 2009 | A1 |
Number | Date | Country |
---|---|---|
103684179 | Mar 2014 | CN |
104579042 | Apr 2015 | CN |
104659990 | May 2015 | CN |
105048913 | Nov 2015 | CN |
110380658 | Oct 2019 | CN |
111245318 | Jun 2020 | CN |
2005229717 | Aug 2005 | JP |
5488043 | May 2014 | JP |
Entry |
---|
“International Search Report (Form PCT/ISA/210) of PCT/CN2021/082326,” dated Nov. 24, 2021, pp. 1-5. |
Number | Date | Country | |
---|---|---|---|
20230008153 A1 | Jan 2023 | US |