CN113612420B - Integration saturation resistant non-inductive permanent magnet synchronous motor rotor position detection method - Google Patents
Integration saturation resistant non-inductive permanent magnet synchronous motor rotor position detection method Download PDFInfo
- Publication number
- CN113612420B CN113612420B CN202110989711.6A CN202110989711A CN113612420B CN 113612420 B CN113612420 B CN 113612420B CN 202110989711 A CN202110989711 A CN 202110989711A CN 113612420 B CN113612420 B CN 113612420B
- Authority
- CN
- China
- Prior art keywords
- integral
- value
- operator
- offset
- rotor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/14—Electronic commutators
- H02P6/16—Circuit arrangements for detecting position
- H02P6/18—Circuit arrangements for detecting position without separate position detecting elements
- H02P6/182—Circuit arrangements for detecting position without separate position detecting elements using back-emf in windings
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/18—Estimation of position or speed
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2203/00—Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
- H02P2203/03—Determination of the rotor position, e.g. initial rotor position, during standstill or low speed operation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
The invention discloses an integral saturation resistant noninductive permanent magnet synchronous motor rotor position detection method, which comprises the following steps: obtaining a voltage VαAnd voltage VβAnd current IαAnd current Iβ(ii) a Constructing a projected magnetic flux Ψ on an axis alpha with respect to a motor magnetic flux ΨαAnd the projected magnetic flux Ψ on the β axisβThe integral operator is used for operating an integral term in the mathematical model; adjusting the bias value in real time according to the integral value of the integral arithmetic unit; for the magnetic flux ΨαAnd magnetic flux ΨβPerforming high-pass filtering and performing arc tangent operation to obtain the original angle position theta of the rotorraw(ii) a For the original angular position theta of the rotorrawCarrying out differential operation to obtain the original angular velocity omega of the rotorraw(ii) a Rotor original angular velocity omega along with pairrawPerforming a low-pass filter to obtain a rotor angular velocity omega; and obtaining the rotor angular position theta according to the rotor angular speed omega.
Description
Technical Field
The invention relates to the technical field of motor control, in particular to an integral saturation resistant non-inductive permanent magnet synchronous motor rotor position detection method and system.
Background
To maximize efficiency in motor control, it is necessary to know the real-time position information of the rotor. The position information can be obtained from a sensor or estimated by a certain algorithm. Among them, the former is called a sensible scheme, and the latter is called a non-sensible scheme. Compared with an inductive scheme, the non-inductive scheme has the advantages of low cost, adaptation to a high-temperature and high-pressure environment in which the sensor may fail and the like.
In the non-inductive solution, the algorithm module that detects the rotor position is called observer. According to the principle of the algorithm, observers can be divided into two categories: the first type is a saliency-based observer. The algorithm extracts position information by a high-frequency injection method by means of the characteristic that the inductance of the motor is not isotropic, and has the advantages of being suitable for low-speed scenes and insensitive to motor parameters. The disadvantage is that when the motor saliency is weak, the position estimation accuracy is greatly degraded, and high frequency injection generates noise. The second type is an observer based on a back electromotive force model. This type of algorithm extracts the position information contained therein by estimating the back emf generated as the motor rotates. Common extraction algorithms include a kalman filter algorithm, an extended flux linkage algorithm, a direct integration algorithm based on a stator voltage model, and the like. In which the direct integration algorithm is simple to calculate, but the dc offset present in the system saturates the integrator, thereby greatly reducing the position estimation accuracy. The Kalman filtering algorithm, the extended flux linkage algorithm and the like do not have the phenomenon of integrator saturation, but the operation is too complex, and the calculation cost paid in engineering is too large. And for different motors, values of some parameters introduced in the algorithms are different, and how to select proper parameters causes great challenges for users.
Based on the principle and the defect analysis of the existing position observer, the invention provides a noninductive permanent magnet synchronous motor position detection method capable of resisting integrator saturation based on a back electromotive force model. The method solves the problem of integrator saturation in a direct integration algorithm based on a stator voltage model, has small operand, saves logic resources and is easy to realize in engineering. Meanwhile, the algorithm has better parameter robustness, does not need a user to adjust the algorithm parameters, and is easy to popularize.
Disclosure of Invention
The invention aims to solve the technical problem that an integrator is saturated in a direct integration algorithm for detecting the position of a motor rotor based on a stator voltage model, has small operand, saves logic resources and is easy to realize in engineering.
In order to solve the technical problem, the invention provides an integral saturation resistant rotor position detection method for a non-inductive permanent magnet synchronous motor, which comprises the following steps:
acquiring applied voltage V on alpha and beta axes under motor stator coordinate systemαAnd voltage VβAnd current I on alpha and beta axesαAnd current Iβ;
Constructing a projected magnetic flux Ψ on an axis alpha with respect to a motor magnetic flux ΨαAnd the projected magnetic flux Ψ on the β axisβThe mathematical model of (a), said mathematical model being represented as:
in the mathematical model, R is the phase resistance of the motor, L is the phase inductance of the motor, psi is the magnetic flux of the motor, and theta is the angular position of the electronic rotor;in order to be the first integral term,as a second integral term, OffsetαIs the first Offset value, OffsetβIs a second offset value;
using a first integral operator to operate the first integral term, and using a second integral operator to operate the second integral term; and real-time adjusting the first Offset value Offset according to the integral values of the first integral operator and the second integral operatorαAnd a second Offset value OffsetβAdjusting;
magnetic flux ΨαAnd magnetic flux ΨβInputting the magnetic flux Ψ filtered by the high-pass filterα_hpfAnd magnetic flux Ψβ_hpf(ii) a And to the magnetic flux Ψα_hpfAnd a magnetic flux Ψβ_hpfPerforming arc tangent operation to obtain the original angle position theta of the rotorraw;
For the original angular position theta of the rotorrawCarrying out differential operation to obtain the original angular velocity omega of the rotorraw;
Then the original angular speed omega of the rotor is measuredrawInputting the high-frequency noise to a low-pass filter to eliminate high-frequency noise introduced by a difference arithmetic unit so as to obtain a rotor angular speed omega;
obtaining a rotor angle position compensation value delta theta according to the rotor angular velocity omega, and comparing the rotor angle position compensation value delta theta with the original rotor angle position thetarawThe rotor angle position θ is obtained by addition.
In one embodiment, the first Offset value Offset is based on the operation values of the first and second integral operatorsαAnd a second Offset value OffsetβThe adjustment is specifically as follows:
let the maximum allowable value of the first integral arithmetic unit and the second integral arithmetic unit be TmaxMinimum tolerance value of Tmin;
When the integral value of the first integral operator is inToIn between, the first Offset value OffsetαSet to 0; when the integral value of the first integral operator is inToIn between, the first Offset value OffsetαA is set as a, and a is a positive number; when the integral value of the first integral operator is inToIn between, the first Offset value OffsetαSet to-a;
when the integral value of the second integral operator is inToIn between, the second Offset value OffsetβSet to 0; when the integral value of the second integral operator is inToIn between, the second Offset value OffsetβA is set as a, and a is a positive number; when the integral value of the second integral operator is inToIn between, the second Offset value OffsetβSet to-a.
In one embodiment, the first Offset value OffsetαAnd a second Offset value OffsetβThe value of (a) is 1.
In one embodiment, the magnetic flux ΨαAnd magnetic flux ΨβFiltered by the first Offset value Offset after being input into the high-pass filterαAnd a second Offset value OffsetβThe resulting dc offset will be filtered out.
In one embodiment, the angular position compensation value delta theta at the rotor angular velocity omega is obtained according to the rotor angular velocity omega and by searching the angular velocity and compensation value numerical relation table,
according to another aspect of the present invention, there is also provided an anti-integral saturation non-inductive permanent magnet synchronous motor rotor position detection system, which includes a first gain operator 1, a first addition operator 2, a first integrator operator 3, a second gain operator 4, and a second addition operator 5; when the applied voltage on the alpha axis under the motor stator coordinate system is electricityPressure VαAnd the applied current on the alpha axis is the current IαWhen the current is over; current IαPasses through a first gain operator 1 and then is input to a first addition operator 2, and the voltage VαAnd the first Offset value OffsetαSimultaneously input to a first adder 2, the output value of the first adder 2 is input to a first integration operator 3, the integration result of the first integration operator 3 is input to a second adder 5, and a current I is simultaneously inputαThe magnetic flux passes through the second gain operator 4 and is input to the second adder operator 5, and the projection amount Ψ of the motor magnetic flux Ψ on the a axis is obtained by the calculation of the second adder operator 5α;
The system further comprises a third gain operator 6, a third adder operator 7, a second integrator operator 8, a fourth gain operator 9 and a fourth adder operator 10; when the applied voltage on the beta axis under the motor stator coordinate system is a voltage VβAnd the current applied to the beta axis is the current IβWhen the current is over; current IβPasses through a third gain operator 6 and is input to a third addition operator 7, the voltage VαAnd a second Offset value OffsetβSimultaneously input to a third addition operator 7, the output value of the third addition operator 7 is input to a second integral operator 8, the integral operation result of the second integral operator 8 is input to a fourth addition operator 10, and simultaneously the current IβThe magnetic flux is input to a fourth addition operator 10 after passing through a fourth gain operator 9, and a projection amount Ψ of the motor magnetic flux Ψ on the β axis is obtained through calculation of the fourth addition operator 10β
The system further comprises a high-pass filter 11 and a second high-pass filter 12, the magnetic flux ΨαThe input to the first high pass filter 11 results in a filtered magnetic flux Ψα_hpf(ii) a Magnetic flux ΨβInput to a second high pass filter 12 resulting in a filtered magnetic flux Ψβ_hpf;
The system further comprises an arctangent function operator 13, a difference operator 14, a low-pass filter 15 and a fifth addition operator 16; magnetic flux Ψ to be high-pass filteredα_hpfAnd a magnetic flux Ψβ_hpfSending the angle to an arctangent function operator 13 to obtain the original angle position of the rotorSet thetaraw(ii) a The original angular position thetarawThe input difference arithmetic unit 14 can obtain the original angular velocity ω of the rotorrawThe original angular velocity ωrawInputting the angular speed to a low-pass filter 15 to obtain the rotor angular speed omega; obtaining an angle position compensation value delta theta according to the rotor angular speed omega; the original angular velocity ωrawAnd the angle position compensation value delta theta are sent to a fifth arithmetic unit 16, and the angle position theta of the electronic rotor is obtained.
In one embodiment, the integration operation in the first integrator operator 3 includes the Offset of the first Offset valueαIntegral operation of (1); the integration operation in the second integration operator 8 includes a second OffsetβIntegral operation of (1); let the maximum allowable value of the first integral operator 3 and the second integral operator 8 be TmaxMinimum tolerance value of Tmin(ii) a Then:
when the integral value of the first integral operator 3 is atToIn between, the first Offset value OffsetαSet to 0; when the integral value of the first integral operator 3 is atToIn between, the first Offset value OffsetαA is set as a, and a is a positive number; when the integral value of the first integral operator 3 is atToIn between, the first Offset value OffsetαSet to-a;
when the second integral arithmetic unitThe integral value of 8 isToIn between, the second Offset value OffsetβSet to 0; when the integral value of the second integral operator 8 is atToIn between, the second Offset value OffsetβA is set as a, and a is a positive number; when the integral value of the second integral operator 8 is atToIn between, the second Offset value OffsetβSet to-a.
In one embodiment, the first Offset value OffsetαAnd a second Offset value OffsetβThe value of (a) is 1.
One or more embodiments of the present invention may have the following advantages over the prior art:
1. in the invention, a dynamic bias value is introduced into the integral operation, thereby realizing the effective integral saturation resistance effect, and the dynamic bias value can be filtered and eliminated by a high-pass filter without influencing the final calculation result;
2. compared with the traditional Kalman filtering algorithm, the calculation amount of the method is greatly reduced by the extended flux linkage algorithm, and the integral saturation resistant effect is kept equal.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 is a schematic structural diagram of a rotor position detection system of an anti-integral saturation non-inductive PMSM according to an embodiment of the present invention;
fig. 2 is a schematic flow chart of a rotor position detection method of an anti-integral saturation non-inductive permanent magnet synchronous motor according to an example of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to fig. 1-2.
Fig. 1 is a schematic structural diagram of an anti-integral saturation sensorless permanent magnet synchronous motor rotor position detection system according to an embodiment of the invention. The present invention will be described with reference to fig. 1.
In the field of motor technology, in a stator coordinate system that uses a motor stator as a reference object, that is, an α - β coordinate system, a voltage model expression of a motor is as follows:
wherein, the voltage VαAnd voltage VβIs applied voltage on alpha, beta axis, current IαAnd current IβThe current on the alpha and beta axes acquired by the ADC circuit is shown, R is the phase resistance of the motor, L is the phase inductance of the motor, psi is the magnetic flux of the motor, and theta is the angular position of the electronic rotor.
As shown in FIG. 1, current IαPasses through a first gain operator 1 and then is input to a first adderLaw operator 2, voltage VαAnd the first Offset value OffsetαSimultaneously input to a first adder 2, the output value of the first adder 2 is input to a first integration operator 3, the integration result of the first integration operator 3 is input to a second adder 5, and a current I is simultaneously inputαThe magnetic flux passes through the second gain operator 4 and is input to the second adder operator 5, and the projection amount Ψ of the motor magnetic flux Ψ on the a axis is obtained by the calculation of the second adder operator 5α. The mathematical expression of the above process is:
from the above equation, the first gain operator 1 is used to complete the term First addition operator 2 for completing the termThe first integral operator 3 is used to complete the termOf the second gain operator 4, the termOf the second adder 5, the second adder completing the termAnd (4) performing the operation of (1).
Likewise, current IβPasses through a third gain operator 6 and is input to a third addition operator 7, the voltage VαAnd a second Offset value OffsetβSimultaneously input to a third addition operator 7, the output value of the third addition operator 7 is input to a second integral operator 8, the integral operation result of the second integral operator 8 is input to a fourth addition operator 10, and simultaneously the current IβBy a fourth gainThe arithmetic unit 9 is input to a fourth addition arithmetic unit 10, and the projection amount Ψ of the motor magnetic flux Ψ on the β axis is obtained through the arithmetic operation of the fourth addition arithmetic unit 10β. The mathematical expression of the above process is:
from the above equation, the third gain operator 6 is used to complete the termIs used to complete the term, the third addition operator 7 is used to complete the termA second integral operator 8 for performing the termOf the fourth gain operator 9Is performed by a fourth addition operator 10And (4) performing the operation of (1).
Note that the maximum allowable value of the first integration operator 3 and the second integration operator 8 is TmaxMinimum tolerance value of TminThen, the present embodiment provides:
when the integral value of the first integral operator 3 is atToIn between, the first Offset value OffsetαIs set to 0, i.e. the first Offset value Offset at this timeαNo offset is superimposed on the first integrating operator 3.
When the integral value of the first integral operator 3 is atToIn between, the first Offset value OffsetαA is set to a, a is positive, i.e. the first Offset value Offset at this timeαThe first integral arithmetic unit 3 is superposed with forward bias to ensure that the integral result is continuously increased along with the integral operation, thereby avoiding the operation result of the first integral arithmetic unit 3 from reaching the minimum accommodating value Tmin。
When the integral value of the first integral operator 3 is atToIn between, the first Offset value OffsetαIs set to-a, i.e. the first Offset value Offset at this timeαThe first integral arithmetic unit 3 is superposed with negative bias to ensure that the integral result is continuously reduced along with the integral operation, thereby avoiding the operation result of the first integral arithmetic unit 3 from reaching the maximum accommodating value Tmax。
Similarly, when the integral value of the second integral operator 8 is atToIn between, the second Offset value OffsetβIs set to 0, i.e. the second Offset value Offset at this timeβNo offset is superimposed on the second integrator operator 8.
When the integral value of the second integral operator 8 is atToIn between, the second Offset value OffsetβA is set to a, a is positive, i.e. the second Offset value Offset at this timeβThe second integral arithmetic unit 8 is superposed with forward bias to ensure that the integral result is continuously increased along with the integral operation, thereby avoiding the operation result of the second integral arithmetic unit 8 from reaching the minimum accommodating value Tmin。
When the integral value of the second integral operator 8 is atToIn between, the second Offset value OffsetβIs set to-a, i.e. the second Offset value Offset at this timeβThe second integral arithmetic unit 8 is superposed with negative bias to ensure that the integral result is continuously reduced along with the integral operation, thereby avoiding the operation result of the second integral arithmetic unit 8 from reaching the maximum accommodating value Tmax。
Projection amount psi on alpha axis of obtained motor magnetic flux psiαAnd the projection amount Ψ on the β axisβThen, the magnetic flux ΨαThe input to the first high pass filter 11 results in a filtered magnetic flux Ψα_hpf(ii) a Magnetic flux ΨβInput to a second high pass filter 12 resulting in a filtered magnetic flux Ψβ_hpf. Due to magnetic flux ΨαWith the magnetic flux ΨβIncluding a DC Offset, i.e. the first Offset value Offset introduced in the previous integration operationαAnd a second Offset value OffsetβTherefore, the influence of the dc offset on the operation result needs to be eliminated in the subsequent calculation. And the first Offset value OffsetαAnd a second Offset value OffsetβDuring the integration, the frequency changes among three values of-a, 0 and a along with the change of the integrated value, and the change belongs to a slowly-changing signal and is extremely low. Thus in the magnetic flux ΨαWith the magnetic flux ΨβAfter passing through the high frequency filter, is determined by the first Offset value OffsetαAnd a second Offset value OffsetβThe resulting dc offset will be filtered out.
When the first Offset value OffsetαAnd a second Offset value OffsetβThe larger the value of (a), the more obvious the bias effect is, i.e. the larger the value of (b), the faster the integral value can be brought away from the integral saturation state, but the larger the change of the signal is. An Offset value Offset is added to the input term of the integrator for dynamically adjusting the range of the integration value to avoid the integration value approaching the upper or lower bound of the integrator, i.e. to avoid saturation.
The offset value acts as an artifact, and the influence of the offset value on the calculated value needs to be eliminated by a high-pass filter of a later stage. According to the high-pass filter characteristic, in order to filter the bias value influence, it should have a "slow-varying" characteristic. According to the principle of dynamically adjusting the offset value, the smaller the amplitude of the offset value, the slower the whole integrator results to jump between the upper and lower limits, and the correspondingly lower the adjustment frequency of the offset value, the more easily it introduces effects that can be eliminated by the high-pass filter. Therefore, the bias value set should have three elements: the minimum positive value that the system can represent, the minimum negative value that the system can represent, and 0. For a fixed-point digital system, 1 is the minimum non-zero positive value that it can represent, so the value of a in this embodiment is + 1.
The high-pass filtered magnetic flux Ψ is thenα_hpfAnd a magnetic flux Ψβ_hpfIs fed to an arctangent function operator 13 to obtain the original angular position theta of the rotorrawNamely:
the original angular position thetarawThe input difference operator 14 can obtain the original angular velocity ω of the rotorrawNamely:
wherein T is the updating time interval of the observer of the non-inductive motor control system, and when the system observes and estimates the rotor position every fixed time interval T, thetaraw[n]Is the original angular position estimated the nth time.
The original angular velocity ω is then measuredrawAnd inputting the signal into a low-pass filter 15, and eliminating high-frequency noise introduced by a difference arithmetic unit to obtain the angular velocity omega of the motor rotor.
And obtaining an angular position compensation value delta theta under the rotor angular speed omega according to a rotor angular speed omega and rotor angular position compensation relation curve known in the art.
The original angular velocity ωrawAnd the angle position compensation value delta theta are sent to a fifth arithmetic unit 16, and the angle position theta of the electronic rotor is obtained.
Thus, the position detection of the motor rotor is completed after the angular speed omega and the angular position theta of the motor rotor are obtained.
The above description is only an embodiment of the present invention, and the protection scope of the present invention is not limited thereto, and any person skilled in the art should modify or replace the present invention within the technical specification of the present invention.
Claims (10)
1. An integral saturation resistant sensorless permanent magnet synchronous motor rotor position detection method is characterized by comprising the following steps:
acquiring applied voltage V on alpha and beta axes under motor stator coordinate systemαAnd voltage VβAnd current I on alpha and beta axesαAnd current Iβ;
Constructing a projected magnetic flux Ψ on an axis alpha with respect to a motor magnetic flux ΨαAnd the projected magnetic flux Ψ on the β axisβThe mathematical model of (a), said mathematical model being represented as:
in the above mathematical model, RIs the phase resistance of the motor, L is the phase inductance of the motor, psi is the magnetic flux of the motor;in order to be the first integral term,as a second integral term, OffsetαIs the first Offset value, OffsetβIs a second offset value;
using a first integral operator to operate the first integral term, and using a second integral operator to operate the second integral term; introducing direct current offset into the first offset value and the second offset value in integral operation, and adjusting the first offset value in real time according to the integral value of the first integral operator; adjusting the second offset value in real time according to the integral value of the second integral arithmetic unit; the first bias value and the second bias value change among three values of-a, 0 and a along with the change of the integral value in the integration process;
magnetic flux ΨαAnd magnetic flux ΨβThe magnetic flux psi is obtained after being filtered by an input high-pass filterα_hpfAnd magnetic flux Ψβ_hpf(ii) a And to the magnetic flux Ψα_hpfAnd a magnetic flux Ψβ_hpfPerforming arc tangent operation to obtain the original angle position theta of the rotorraw;
For the original angular position theta of the rotorrawCarrying out differential operation to obtain the original angular velocity omega of the rotorraw;
Then the original angular speed omega of the rotor is measuredrawInputting the high-frequency noise to a low-pass filter to eliminate high-frequency noise introduced by a difference arithmetic unit so as to obtain a rotor angular speed omega;
obtaining a rotor angle position compensation value delta theta according to the rotor angular velocity omega, and comparing the rotor angle position compensation value delta theta with the original rotor angle position thetarawThe rotor angle position θ is obtained by addition.
2. The sensorless permanent magnet synchronous motor rotor position detecting method according to claim 1, wherein the first offset value is adjusted in real time according to an integrated value of the first integration operator; the adjusting the second offset value in real time according to the integral value of the second integral arithmetic unit specifically comprises:
let the maximum allowable value of the first integral arithmetic unit and the second integral arithmetic unit be TmaxMinimum tolerance value of Tmin;
When the integral value of the first integral operator is inToIn between, the first Offset value OffsetαSet to 0; when the integral value of the first integral operator is inToIn between, the first Offset value OffsetαA is set as a, and a is a positive number; when the integral value of the first integral operator is inToIn between, the first Offset value OffsetαSet to-a;
when the integral value of the second integral operator is inToIn between, the second Offset value OffsetβSet to 0; when the integral value of the second integral operator is inToIn between, the second Offset value OffsetβA is set as a, and a is a positive number; when the integral value of the second integral operator is inToIn between, the second Offset value OffsetβSet to-a.
3. The sensorless permanent magnet synchronous motor rotor position detection method according to claim 2, characterized in that the first Offset value OffsetαAnd a second Offset value OffsetβThe value of (a) is 1.
4. The sensorless permanent magnet synchronous motor rotor position detection method according to claim 1, characterized in that magnetic flux Ψ is configured toαAnd magnetic flux ΨβFiltered by the first Offset value Offset after being input into the high-pass filterαAnd a second Offset value OffsetβThe resulting dc offset will be filtered out.
5. The method for detecting the rotor position of the sensorless permanent magnet synchronous motor according to claim 1, wherein the angular position compensation value Δ θ at the rotor angular velocity ω is obtained according to the rotor angular velocity ω and by searching a relationship table of the angular velocity and the compensation value.
6. The system is characterized by comprising a first gain operator (1), a first addition operator (2) and a first integral operationA device (3), a second gain operator (4) and a second adder (5); when the applied voltage on the alpha axis under the motor stator coordinate system is voltage VαAnd the applied current on the alpha axis is the current IαWhen the current is over; current IαThe voltage V is input to a first addition operator (2) after passing through a first gain operator (1)αAnd the first Offset value OffsetαSimultaneously input to a first addition operator (2), the output value of the first addition operator (2) is input to a first integration operator (3), the integration result of the first integration operator (3) is input to a second addition operator (5), and simultaneously current IαThe magnetic flux is input to a second adder (5) after passing through a second gain operator (4), and a projection amount Ψ of the motor magnetic flux Ψ on an axis alpha is obtained by the operation of the second adder (5)α;
The system also comprises a third gain operator (6), a third addition operator (7), a second integral operator (8), a fourth gain operator (9) and a fourth addition operator (10); when the applied voltage on the beta axis under the motor stator coordinate system is a voltage VβAnd the current applied to the beta axis is the current IβWhen the current is over; current IβPasses through a third gain operator (6) and is input to a third addition operator (7) to obtain a voltage VαAnd a second Offset value OffsetβSimultaneously input to a third addition operator (7), the output value of the third addition operator (7) is input to a second integral operator (8), the integral operation result of the second integral operator (8) is input to a fourth addition operator (10), and simultaneously current IβThe magnetic flux is input into a fourth addition arithmetic unit (10) after passing through a fourth gain arithmetic unit (9), and the projection amount psi of the motor magnetic flux psi on the beta axis is obtained through the operation of the fourth addition arithmetic unit (10)β;
The integration operation in the first integration operator (3) comprises an integration operation on a first bias value; the integration operation in the second integration operator (8) includes an integration operation of a second offset value; introducing direct current offset into the first offset value and the second offset value in integral operation, and adjusting the first offset value in real time according to the integral value of the first integral operator; adjusting the second offset value in real time according to the integral value of the second integral arithmetic unit; the first bias value and the second bias value change among three values of-a, 0 and a along with the change of the integral value in the integration process;
the system further comprises a first high-pass filter (11) and a second high-pass filter (12), the magnetic flux ΨαIs input to a first high-pass filter (11) to obtain a filtered magnetic flux Ψα_hpf(ii) a Magnetic flux ΨβIs input to a second high-pass filter (12) to obtain a filtered magnetic flux Ψβ_hpf;
The system also comprises an arctangent function arithmetic unit (13), a difference arithmetic unit (14), a low-pass filter (15) and a fifth addition arithmetic unit (16); magnetic flux Ψ to be high-pass filteredα_hpfAnd a magnetic flux Ψβ_hpfFeeding the rotor into an arctangent function operator (13) to obtain the original angular position theta of the rotorraw(ii) a The original angular position thetarawThe original angular velocity omega of the rotor can be obtained by inputting a difference arithmetic unit (14)rawThe original angular velocity ωrawInputting the angular speed to a low-pass filter (15) to obtain the angular speed omega of the rotor; obtaining an angle position compensation value delta theta according to the rotor angular speed omega; the original angular position thetarawAnd the angle position compensation value delta theta are sent to a fifth arithmetic unit (16), and then the angle position theta of the motor rotor is obtained.
7. The sensorless permanent magnet synchronous motor rotor position detection system according to claim 6, characterized in that the maximum allowable value of the first integration operator (3) and the second integration operator (8) is TmaxMinimum tolerance value of Tmin(ii) a And has the following components:
when the integral value of the first integral operator (3) is inToIn between, the first Offset value OffsetαSet to 0; when the integral value of the first integral operator (3) is inToIn between, the first Offset value OffsetαA is set as a, and a is a positive number; when the integral value of the first integral operator (3) is inToIn between, the first Offset value OffsetαSet to-a;
when the integral value of the second integral operator (8) is inToIn between, the second Offset value OffsetβSet to 0; when the integral value of the second integral operator (8) is inToIn between, the second Offset value OffsetβA is set as a, and a is a positive number; when the integral value of the second integral operator (8) is inToIn between, the second Offset value OffsetβSet to-a.
8. The sensorless permanent magnet synchronous motor rotor position detection system according to claim 6, characterized in that the first Offset value OffsetαAnd a second Offset value OffsetβThe value of (a) is 1.
9. A control method of a sensorless permanent magnet synchronous motor, characterized in that the motor rotor position detection method in the control method of the sensorless permanent magnet synchronous motor uses the rotor position detection method of the sensorless permanent magnet synchronous motor according to any one of claims 1 to 5.
10. A control system of a sensorless permanent magnet synchronous motor, characterized in that the sensorless permanent magnet synchronous motor rotor position detection system according to any one of claims 6 to 8 is used in the control system for the sensorless permanent magnet synchronous motor to detect the motor rotor position.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110989711.6A CN113612420B (en) | 2021-08-26 | 2021-08-26 | Integration saturation resistant non-inductive permanent magnet synchronous motor rotor position detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110989711.6A CN113612420B (en) | 2021-08-26 | 2021-08-26 | Integration saturation resistant non-inductive permanent magnet synchronous motor rotor position detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113612420A CN113612420A (en) | 2021-11-05 |
CN113612420B true CN113612420B (en) | 2022-02-18 |
Family
ID=78342128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110989711.6A Active CN113612420B (en) | 2021-08-26 | 2021-08-26 | Integration saturation resistant non-inductive permanent magnet synchronous motor rotor position detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113612420B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114244195A (en) * | 2021-12-30 | 2022-03-25 | 中原内配(上海)电子科技有限公司 | PID-based brushless direct current motor control method, system, device and storage medium |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1339871A (en) * | 2000-08-18 | 2002-03-13 | Lg产电株式会社 | Vector control system without sensor for inductive motor and method |
CN101286725A (en) * | 2008-04-17 | 2008-10-15 | 戴政 | Deriving method of motor rotative velocity and location of rotor in synchronous electric machine vector control system |
EP2180386A1 (en) * | 2008-10-24 | 2010-04-28 | Abb Research Ltd. | Control device and method for operating a control device |
CN101719753A (en) * | 2009-12-30 | 2010-06-02 | 国网电力科学研究院 | Method for detecting initial position of rotor of static frequency-conversion starting motor of pumped storage power station |
CN107171610A (en) * | 2017-06-29 | 2017-09-15 | 广东美芝制冷设备有限公司 | Rotor position estimate method, rotor position estimate device and motor |
CN107809191A (en) * | 2017-10-24 | 2018-03-16 | 厦门金龙汽车新能源科技有限公司 | A kind of angular observation method of synchronous motor zero-velocity sensor |
CN208782746U (en) * | 2018-03-09 | 2019-04-23 | 核工业理化工程研究院 | The control system of permanent-magnet synchronous motor rotor position angle |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7265507B1 (en) * | 2006-02-20 | 2007-09-04 | Hamilton Sundstrand Corporation | Angular position and velocity estimation for synchronous machines based on extended rotor flux |
-
2021
- 2021-08-26 CN CN202110989711.6A patent/CN113612420B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1339871A (en) * | 2000-08-18 | 2002-03-13 | Lg产电株式会社 | Vector control system without sensor for inductive motor and method |
CN101286725A (en) * | 2008-04-17 | 2008-10-15 | 戴政 | Deriving method of motor rotative velocity and location of rotor in synchronous electric machine vector control system |
EP2180386A1 (en) * | 2008-10-24 | 2010-04-28 | Abb Research Ltd. | Control device and method for operating a control device |
CN101719753A (en) * | 2009-12-30 | 2010-06-02 | 国网电力科学研究院 | Method for detecting initial position of rotor of static frequency-conversion starting motor of pumped storage power station |
CN107171610A (en) * | 2017-06-29 | 2017-09-15 | 广东美芝制冷设备有限公司 | Rotor position estimate method, rotor position estimate device and motor |
CN107809191A (en) * | 2017-10-24 | 2018-03-16 | 厦门金龙汽车新能源科技有限公司 | A kind of angular observation method of synchronous motor zero-velocity sensor |
CN208782746U (en) * | 2018-03-09 | 2019-04-23 | 核工业理化工程研究院 | The control system of permanent-magnet synchronous motor rotor position angle |
Non-Patent Citations (2)
Title |
---|
无刷交流同步电机起动阶段无位置传感器控制技术研究;杨袁钰;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20190215(第2期);第21-23页 * |
高速永磁无刷直流电机换相校正技术探究;施晓青;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20200215(第2期);第11-15页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113612420A (en) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
An et al. | Frequency-adaptive complex-coefficient filter-based enhanced sliding mode observer for sensorless control of permanent magnet synchronous motor drives | |
Raca et al. | Robust magnet polarity estimation for initialization of PM synchronous machines with near-zero saliency | |
CN105227010B (en) | A kind of permagnetic synchronous motor position-sensor-free position detection error harmonic pulse removing method | |
JP2004135494A (en) | Sensor-less position control algorithm for ac machine | |
CN108900127B (en) | IPMSM low-speed section position-sensor-free control method considering cross coupling effect | |
Chen et al. | Adaptive sliding mode observer-based sensorless control for SPMSM employing a dual-PLL | |
KR20070107539A (en) | Apparatus and method for controlling high speed driving of sensorless permanent magnet synchoronous mortor | |
CN108347207B (en) | Position and speed estimation method of permanent magnet synchronous motor based on complex PI controller | |
CN111106767B (en) | Sensorless starting control method of permanent magnet synchronous motor | |
CN102668361B (en) | Motor control apparatus and magnetic-pole position detection method therefor | |
CN109600082B (en) | Permanent magnet synchronous motor rotor position full-order sliding mode observation device and method | |
CN113612420B (en) | Integration saturation resistant non-inductive permanent magnet synchronous motor rotor position detection method | |
CN108574444A (en) | One kind being used for initial position detection method for permanent magnet synchronous electric motor rotor | |
WO2020186656A1 (en) | Method for estimating position of linear time-variant rotor of low-speed permanent magnet synchronous motor | |
Brandstetter et al. | Sensorless control of permanent magnet synchronous motor using voltage signal injection | |
CN110620533A (en) | Surface-mounted permanent magnet synchronous motor sensorless control method | |
CN110034709B (en) | Rotor position estimation method suitable for surface-mounted permanent magnet synchronous motor | |
US9966884B2 (en) | Method and device for determining the rotor position and speed of a rotating field machine | |
CN114598206A (en) | Design method of permanent magnet synchronous motor wide-speed-domain rotor position observer | |
CN108512475A (en) | A kind of rotor position detecting method based on pulsating high frequency signal injection method | |
CN113364375B (en) | Sensorless control method for PMSM (permanent magnet synchronous motor) driving system of variable-structure current regulator | |
CN109150044A (en) | A kind of control method of the permanent magnet synchronous motor based on high frequency electrocardiography and immune PID | |
CN101132163B (en) | DTC magnetic linkage estimation filtering method with compensation ability | |
CN109194226A (en) | Parameter of electric machine on-line identification method, system, electric machine controller and storage medium | |
CN115940721A (en) | Permanent magnet synchronous motor observation and filtering method based on supercoiled sliding-mode observer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |