[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN110112965B - Counter electromotive force observation method for permanent magnet synchronous motor - Google Patents

Counter electromotive force observation method for permanent magnet synchronous motor Download PDF

Info

Publication number
CN110112965B
CN110112965B CN201910376079.0A CN201910376079A CN110112965B CN 110112965 B CN110112965 B CN 110112965B CN 201910376079 A CN201910376079 A CN 201910376079A CN 110112965 B CN110112965 B CN 110112965B
Authority
CN
China
Prior art keywords
current
motor
electromotive force
obtaining
group
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.)
Expired - Fee Related
Application number
CN201910376079.0A
Other languages
Chinese (zh)
Other versions
CN110112965A (en
Inventor
郭磊磊
金楠
秦世耀
王瑞明
代林旺
罗魁
曹玲芝
李琰琰
武洁
吴振军
窦智峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Electric Power Research Institute Co Ltd CEPRI
Zhengzhou University of Light Industry
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
Zhengzhou University of Light Industry
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Electric Power Research Institute Co Ltd CEPRI, Zhengzhou University of Light Industry filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201910376079.0A priority Critical patent/CN110112965B/en
Publication of CN110112965A publication Critical patent/CN110112965A/en
Application granted granted Critical
Publication of CN110112965B publication Critical patent/CN110112965B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14Electronic commutators
    • H02P6/16Circuit arrangements for detecting position
    • H02P6/18Circuit arrangements for detecting position without separate position detecting elements
    • H02P6/182Circuit arrangements for detecting position without separate position detecting elements using back-emf in windings

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

The invention provides a method for observing back electromotive force of a permanent magnet synchronous motor, which comprises the following steps: firstly, the voltage and the current of the stator of the motor are both static by sampling and calculation
Figure DEST_PATH_IMAGE001
A voltage component and a current component in a coordinate system; secondly, calculating according to the voltage component and the current component to obtain estimated stator current and a high-frequency sliding mode signal; secondly, performing low-pass filtering on the high-frequency sliding mode signal twice respectively to obtain an estimated actual rotor position angle and a back electromotive force amplitude of the motor; and finally, calculating to obtain the estimated motor back electromotive force according to the estimated actual rotor position angle and the back electromotive force amplitude of the motor. The invention adopts a double low pass filter method, compensates the problems of back electromotive force phase delay and amplitude attenuation caused by the first low pass filter according to the influence of the second low pass filter on the back electromotive force phase and amplitude, and can realize accurate observation of the back electromotive force.

Description

Counter electromotive force observation method for permanent magnet synchronous motor
Technical Field
The invention relates to the field of power electronics, in particular to a method for observing back electromotive force of a permanent magnet synchronous motor.
Background
In recent years, with the increasing energy crisis, photovoltaic power generation technology, wind power generation technology and new energy electric vehicle technology have been developed vigorously. The permanent magnet synchronous motor has the advantages of high power density, high efficiency and the like, and is widely applied to a wind power generation system and a new energy electric vehicle driving system. However, under severe operating environment and environmental factors such as vibration, humidity and low temperature, the speed sensor often has faults such as disconnection and loss of pulse signals, and further the control system fails to operate. Therefore, in order to improve the operational reliability of the permanent magnet synchronous motor control system, a speed sensorless control technique has been widely studied in recent years. In the traditional permanent magnet synchronous motor non-speed sensor control based on a second-order sliding-mode observer, a calculation method is usually adopted to obtain the back electromotive force of a motor, and the rotating speed information of the motor is required to be used in the calculation process, but in the motor non-speed sensor control, the estimated rotating speed of the motor often has a certain rotating speed error, so that the calculated back electromotive force of the motor inevitably has errors.
At present, a method of a speed sensorless control technology of a permanent magnet synchronous motor exists, for example, application number 201610631269.9, the invention name is a speed sensorless control method based on a sliding-mode observer, a method for observing back electromotive force and compensating a motor phase shift angle by using a low-pass filter is provided, and speed sensorless control of the motor is realized. Therefore, the motor back electromotive force used in the control system has an error. Document [ position sensorless control of an interior permanent magnet synchronous motor for a compressor [ J ]. report of electrical and technical, 2013, 28 (5): 182-187 proposes a method of obtaining back electromotive force information in real time by using a series connection of low pass filters, compensating for the phase shift caused by the low pass filters to obtain an accurate motor angle, thereby implementing a speed sensorless control of the motor. However, the method only carries out phase shift compensation on the angle and does not consider the attenuation of the back electromotive force amplitude of the motor caused by the filter. Document [ wanggulin, yangfeng, in swimming, etc. ] built-in permanent magnet synchronous motor position sensorless control [ J ], chinese motor engineering report, 2010, 30 (30): 93-98 ] a speed sensorless control method of a permanent magnet synchronous motor based on a second-order sliding-mode observer is provided, the method only compensates the phase delay of the back electromotive force, and a compensation algorithm needs to use the estimated rotating speed. In the sensorless control of the motor, the estimated motor rotation speed often has a certain rotation speed error, thereby causing an inevitable large error in the calculated motor back electromotive force.
Disclosure of Invention
Aiming at the technical problem that the error of the counter electromotive force of the motor is larger due to the fact that the amplitude of the uncompensated counter electromotive force exists in the existing counter electromotive force calculation method of the permanent magnet synchronous motor, the invention provides the counter electromotive force observation method of the permanent magnet synchronous motor.
The technical scheme of the invention is realized as follows:
a counter electromotive force observation method of a permanent magnet synchronous motor comprises the following steps:
step one, utilizing a voltage sensor to align a motor statorThe voltage is sampled to obtain the stator voltage u of the motorABAnd uBCCalculating three-phase voltage u of the motorA、uBAnd uCAnd the three-phase voltage u is converted intoA、uBAnd uCObtaining the voltage u under a static α - β coordinate system through coordinate transformationαAnd voltage uβ
Step two, sampling the three-phase current of the motor stator by using a current sensor to obtain the current i of the motor statorA、iBAnd iCAnd obtaining the current i under a stationary α - β coordinate system through coordinate transformationαAnd current iβ
Step three, initializing the current of the motor
Figure BDA0002051709770000021
And current
Figure BDA0002051709770000022
Calculating the current
Figure BDA0002051709770000023
With the current i obtained in step twoαIs calculated from the difference of (1), the current is calculated
Figure BDA0002051709770000024
And the current i obtained in step twoβAnd calculating by a sign function to obtain a high-frequency sliding mode signal sαAnd high frequency sliding mode signal sβ
Step four, utilizing the voltage u obtained in the step oneαSubtracting the high-frequency sliding mode signal s obtained in the step threeαObtain a first set of intermediate variables EUsing the voltage u obtained in step oneβSubtracting the high-frequency sliding mode signal s obtained in the step threeβObtain a first set of intermediate variables E
Step five, obtaining a first group of intermediate variables E according to the step four、EAnd the q-axis inductance of the motor is calculated to obtain a second group of intermediate variables E、E
Step six, estimating according to step threeElectric current
Figure BDA0002051709770000025
Electric current
Figure BDA0002051709770000026
Calculating the resistance of the stator of the motor and the q-axis inductance to obtain a third group of intermediate variables E、E
Step seven, utilizing the second group of intermediate variables E obtained in the step fiveSubtracting the third group intermediate variable E obtained in the step sixObtaining a fourth group of intermediate variables EUsing the second set of intermediate variables E obtained in step fiveSubtracting the third group intermediate variable E obtained in the step sixObtaining a fourth group of intermediate variables E
Step eight, obtaining a fourth group intermediate variable E according to the step seven、EUpdating the current in step three
Figure BDA0002051709770000027
And current
Figure BDA0002051709770000028
Further updating high-frequency sliding mode signal sαAnd high frequency sliding mode signal sβ
Step nine, updating the high-frequency sliding mode signal sαAnd high frequency sliding mode signal sβLow-pass filtering through the first low-pass filter to respectively obtain a fifth group of intermediate variables s、s
Step ten, the fifth group of intermediate variables s obtained in the step nine、sLow-pass filtering through a second low-pass filter to respectively obtain a sixth group of intermediate variables s、s
Step eleven, obtaining a fifth group of intermediate variables s according to the step nine、sCalculating effective back electromotive force q-axis deviation
Figure BDA0002051709770000031
Step twelve, the effective counter electromotive force q-axis deviation obtained in the step eleven
Figure BDA0002051709770000032
Calculating to obtain estimated rotating speed by a proportional integrator
Figure BDA0002051709770000033
Thirteen step, rotating speed obtained in the step twelve
Figure BDA0002051709770000034
The rotor position angle theta after primary phase delay generated by a first low-pass filter is obtained through integral adjustment1
Step fourteen, obtaining a sixth group of intermediate variables s according to the step ten、sCalculating the rotor position angle theta after two phase delays generated by the first low-pass filter and the second low-pass filter2
Step fifteen, obtaining the rotor position angle theta according to the step thirteen1Subtracting the rotor position angle theta obtained in the step fourteen2Obtaining a delay angle delta theta, and then obtaining a rotor position angle theta1Calculating the delay angle delta theta to obtain an estimated actual rotor position angle theta;
sixthly, obtaining a fifth group of intermediate variables s according to the step nine、sAnd a sixth set of intermediate variables s obtained in step ten、sCalculating to obtain the attenuation ratio k of the first low-pass filter and the second low-pass filter to the back electromotive force amplitude of the motor;
seventhly, obtaining a fifth group of intermediate variables s according to the step nine、sAnd calculating the attenuation ratio k obtained in the step sixteen to obtain the back electromotive force amplitude e of the motorm
Eighteen, obtaining the counter electromotive force amplitude e according to the seventeenth stepmAnd calculating the estimated back electromotive force of the motor according to the actual rotor position angle theta obtained in the step fifteen
Figure BDA0002051709770000035
And
Figure BDA0002051709770000036
preferably, the three-phase voltage u in the step oneA、uBAnd uCObtaining the voltage u under a static α - β coordinate system through coordinate transformationαAnd voltage uβThe method comprises the following steps:
Figure BDA0002051709770000037
wherein,
Figure BDA0002051709770000038
uABand uBCIs the motor stator voltage.
Preferably, the stator current i of the motor in the second stepA、iBAnd iCObtaining the current i under a stationary α - β coordinate system through coordinate transformationαAnd current iβThe method comprises the following steps:
Figure BDA0002051709770000039
preferably, the high-frequency sliding mode signal s in step threeαAnd high frequency sliding mode signal sβThe obtaining method comprises the following steps:
Figure BDA0002051709770000041
where M is the sliding mode gain and sgn () is the sign function.
Preferably, said first set of intermediate variables E、EThe obtaining method comprises the following steps:
Figure BDA0002051709770000042
the second set of intermediate variables E、EThe obtaining method comprises the following steps:
Figure BDA0002051709770000043
wherein L isqIs a motor q-axis inductor;
the third set of intermediate variables E、EThe obtaining method comprises the following steps:
Figure BDA0002051709770000044
wherein R issIs a motor stator resistor;
the fourth set of intermediate variables E、EThe obtaining method comprises the following steps:
Figure BDA0002051709770000045
preferably, the current in step eight
Figure BDA0002051709770000046
And current
Figure BDA0002051709770000047
The updating method comprises the following steps:
Figure BDA0002051709770000048
wherein, TsIs the sampling period.
Preferably, said fifth set of intermediate variables s、sThe obtaining method comprises the following steps:
Figure BDA0002051709770000049
wherein, ω iscIs the cut-off frequency of the first low-pass filter and the second low-pass filter, s is the laplacian operator;
the sixth set of intermediate variables s、sThe obtaining method comprises the following steps:
Figure BDA00020517097700000410
preferably, the effective back electromotive force q-axis deviation
Figure BDA00020517097700000411
The obtaining method comprises the following steps:
Figure BDA00020517097700000412
wherein, theta1Is the rotor position angle;
the estimated rotation speed
Figure BDA00020517097700000413
The obtaining method comprises the following steps:
Figure BDA00020517097700000414
kpis a proportionality coefficient, kiIs an integral coefficient;
the rotor position angle theta1The updating method comprises the following steps:
Figure BDA00020517097700000415
the rotor position angle theta2The obtaining method comprises the following steps: theta2=arctan(-s/s);
The method for obtaining the delay angle delta theta comprises the following steps: Δ θ ═ θ12
The method for obtaining the estimated actual rotor position angle theta comprises the following steps: theta is equal to theta1+Δθ。
Preferably, the attenuation ratio k of the low-pass filter to the back electromotive force amplitude of the motor is:
Figure BDA0002051709770000051
the back electromotive force amplitude e of the motormThe obtaining method comprises the following steps:
Figure BDA0002051709770000052
preferably, the estimated motor back electromotive force in the step eighteen
Figure BDA0002051709770000053
And
Figure BDA0002051709770000054
the obtaining method comprises the following steps:
Figure BDA0002051709770000055
the beneficial effect that this technical scheme can produce: the method adopts a second-order sliding mode observer to roll the back electromotive force of the permanent magnet synchronous motor, and designs a method which adopts a double low-pass filter method to compensate the amplitude and the phase of the back electromotive force, thereby overcoming the problem that the phase and the amplitude compensation of the back electromotive force are influenced by the rotating speed estimation error, and improving the observation precision of the back electromotive force.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 shows the current i of the present inventionαAnd (4) an observer block diagram.
FIG. 2 shows the current i of the present inventionβAnd (4) an observer block diagram.
Fig. 3 is an overall block diagram of the present invention.
FIG. 4 is a graph of the relationship between the back EMF and the actual back EMF obtained by the present invention; e.g. of the typeαAnd eβFor the actual back-emf of the motor,
Figure BDA0002051709770000056
and
Figure BDA0002051709770000057
the back electromotive force of the motor obtained by the invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive effort based on the embodiments of the present invention, are within the scope of the present invention.
As shown in fig. 3, the invention provides a method for observing back electromotive force of a permanent magnet synchronous motor, which comprises the steps of firstly, obtaining a voltage component and a current component of a stator voltage and a stator current of the motor under a static alpha-beta coordinate system by sampling calculation; secondly, calculating according to the voltage component and the current component to obtain estimated stator current and a high-frequency sliding mode signal; secondly, performing low-pass filtering on the high-frequency sliding mode signal twice respectively to obtain an estimated actual rotor position angle and a back electromotive force amplitude of the motor; and finally, calculating to obtain the estimated motor back electromotive force according to the estimated actual rotor position angle and the back electromotive force amplitude of the motor. The method comprises the following specific steps:
step one, sampling the motor stator voltage by using a voltage sensor to obtain the motor stator voltage uABAnd uBCCalculating the three-phase voltage u of the motor according to the formula (1)A、uBAnd uC
Figure BDA0002051709770000061
Then according to the formula (2), the three-phase voltage uA、uBAnd uCObtaining the voltage u under a static α - β coordinate system through coordinate transformationαAnd voltage uβ
Figure BDA0002051709770000062
Step two, sampling the three-phase current of the motor stator by using a current sensor to obtain the current i of the motor statorA、iBAnd iCAnd obtaining the current i under a stationary α - β coordinate system through the coordinate transformation of the formula (3)αAnd current iβ
Figure BDA0002051709770000063
Step three, initializing the current of the motor
Figure BDA0002051709770000064
And current
Figure BDA0002051709770000065
As shown in FIGS. 1 and 2, the current was calculated
Figure BDA0002051709770000066
With the current i obtained in step twoαAnd current of
Figure BDA0002051709770000067
With the current i obtained in step twoβThen obtaining a high-frequency sliding mode signal s by symbol function calculation according to a formula (4)αAnd high frequency sliding mode signal sβ
Figure BDA0002051709770000068
Wherein M is sliding mode gain and current
Figure BDA0002051709770000069
And current
Figure BDA00020517097700000610
Are all 0 and sgn () is a sign function.
Step four, utilizing the voltage u obtained in the step oneαSubtracting the high-frequency sliding mode signal s obtained in the step threeαObtain a first set of intermediate variables EUsing the voltage u obtained in step oneβSubtracting the high-frequency sliding mode signal s obtained in the step threeαAnd high frequency sliding mode signal sβObtain a first set of intermediate variables EAs shown in equation (5):
Figure BDA00020517097700000611
step five, obtaining a first group of intermediate variables E according to the step four、EAnd the q-axis inductance of the motor is calculated to obtain a second group of intermediate variables E、EAs inFormula (6):
Figure BDA0002051709770000071
wherein L isqIs the q-axis inductance of the motor.
Step six, estimating the current according to the step three
Figure BDA0002051709770000072
Electric current
Figure BDA0002051709770000073
Calculating the resistance of the stator of the motor and the q-axis inductance to obtain a third group of intermediate variables E、EAs shown in equation (7):
Figure BDA0002051709770000074
wherein R issIs the motor stator resistance, LqIs the q-axis inductance of the motor.
Step seven, utilizing the second group of intermediate variables E obtained in the step fiveSubtracting the third group intermediate variable E obtained in the step sixObtaining a fourth group of intermediate variables EUsing the second set of intermediate variables E obtained in step fiveSubtracting the third group intermediate variable E obtained in the step sixObtaining a fourth group of intermediate variables EAs shown in equation (8):
Figure BDA0002051709770000075
step eight, obtaining a fourth group intermediate variable E according to the step seven、ERefresh current
Figure BDA0002051709770000076
And current
Figure BDA0002051709770000077
As shown in formula (9):
Figure BDA0002051709770000078
wherein, TsFor a sampling period, current
Figure BDA0002051709770000079
And current
Figure BDA00020517097700000710
Is 0.
According to the current i, as shown in FIGS. 1 and 2αObserver and current iβClosed loop of observer using updated current
Figure BDA00020517097700000711
And current
Figure BDA00020517097700000712
Further updating high frequency sliding mode signal sαHigh frequency sliding mode signal sβ
Step nine, updating the high-frequency sliding mode signal sαAnd high frequency sliding mode signal sβLow-pass filtering through the first low-pass filter to respectively obtain a fifth group of intermediate variables s、sAs shown in equation (10):
Figure BDA00020517097700000713
wherein, ω iscS is the laplacian, the cut-off frequency of the first low-pass filter and the second low-pass filter.
Step ten, the fifth group of intermediate variables s obtained in the step nine、sCarrying out second low-pass filtering through a second low-pass filter to respectively obtain a sixth group of intermediate variables s、sAs shown in formula (11):
Figure BDA0002051709770000081
wherein, ω iscS is the laplacian, the cut-off frequency of the first low-pass filter and the second low-pass filter.
Step eleven, obtaining a fifth group of intermediate variables s according to the step nine、sCalculating effective back electromotive force q-axis deviation
Figure BDA0002051709770000082
As shown in equation (12):
Figure BDA0002051709770000083
wherein, theta1As rotor position angle, rotor position angle θ1Is 0.
Step twelve, the effective counter electromotive force q-axis deviation obtained in the step eleven
Figure BDA0002051709770000084
Calculating to obtain estimated rotating speed by a proportional integrator
Figure BDA0002051709770000085
As shown in equation (13):
Figure BDA0002051709770000086
wherein k ispIs a proportionality coefficient, kiIs the integral coefficient and s is the laplacian operator.
Thirteen step, rotating speed obtained in the step twelve
Figure BDA0002051709770000087
The rotor position angle theta after primary phase delay generated by a first low-pass filter is obtained through integral adjustment1I.e. angle theta to rotor position1Updating is performed as shown in equation (14):
Figure BDA0002051709770000088
wherein, TsIs the sampling period.
Step fourteen, obtaining a sixth group of intermediate variables s according to the step ten、sCalculating the rotor position angle theta after two phase delays generated by the first low-pass filter and the second low-pass filter2As shown in equation (15):
θ2=arctan(-s/s) (15)。
fifteenth, according to the formula (16), the rotor position angle theta obtained in the thirteenth step1Subtracting the rotor position angle theta obtained in the step fourteen2Obtaining a retardation angle Δ θ:
Δθ=θ12(16);
and then utilizing the rotor position angle theta obtained in the step thirteen according to the formula (17)1And the delay angle Δ θ to calculate an estimated actual rotor position angle θ:
θ=θ1+Δθ (17)。
sixthly, obtaining a fifth group of intermediate variables s according to the step nine、sAnd a sixth set of intermediate variables s obtained in step ten、sCalculating to obtain the attenuation ratio k of the first low-pass filter and the second low-pass filter to the back electromotive force amplitude of the motor, as shown in the formula (18):
Figure BDA0002051709770000091
seventhly, obtaining a fifth group of intermediate variables s according to the step nine、sAnd calculating the attenuation ratio k obtained in the step sixteen to obtain the back electromotive force amplitude e of the motormAs shown in equation (19):
Figure BDA0002051709770000092
where k is the attenuation ratio.
Eighteen, according to seventeenThe resulting back EMF amplitude emAnd calculating the estimated back electromotive force of the motor according to the actual rotor position angle theta obtained in the step fifteen
Figure BDA0002051709770000093
And
Figure BDA0002051709770000094
as shown in equation (20):
Figure BDA0002051709770000095
in order to verify the effectiveness of the present invention, simulation verification was performed. DC side voltage U of inverter for simulationdc400V, the rated power of the permanent magnet synchronous motor is 6.6kW, the flux linkage is 0.35Wb, the number of pole pairs is 4, the stator inductance is 12mH, and the stator resistance Rs0.5 omega, rated frequency of 50Hz, rated voltage of 190V, and sampling frequency fsIs 10 kHz. FIG. 4 shows the back EMF calculated by the present invention
Figure BDA0002051709770000096
And the actual back emf value eα、eβThe comparison shows that the estimated back electromotive force can be accurately obtained when the back electromotive force changes, and the amplitude and the phase of the back electromotive force are accurate. Meanwhile, the method does not need the estimated rotating speed of the motor, so the method is simpler and has higher precision.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (3)

1. A counter electromotive force observation method of a permanent magnet synchronous motor is characterized by comprising the following steps:
step one, sampling the motor stator voltage by using a voltage sensor to obtain the motor stator voltage uABAnd uBCCalculating the electricityThree-phase voltage u of machineA、uBAnd uCAnd the three-phase voltage u is converted intoA、uBAnd uCObtaining the voltage u under a static α - β coordinate system through coordinate transformationαAnd voltage uβ
Step two, sampling the three-phase current of the motor stator by using a current sensor to obtain the current i of the motor statorA、iBAnd iCAnd obtaining the current i under a stationary α - β coordinate system through coordinate transformationαAnd current iβ
Step three, initializing the current of the motor
Figure FDA0002564304500000011
And current
Figure FDA0002564304500000012
Calculating the current
Figure FDA0002564304500000013
With the current i obtained in step twoαIs calculated from the difference of (1), the current is calculated
Figure FDA0002564304500000014
And the current i obtained in step twoβAnd calculating by a sign function to obtain a high-frequency sliding mode signal sαAnd high frequency sliding mode signal sβ
Figure FDA0002564304500000015
Wherein M is sliding mode gain, sgn () is a sign function;
step four, utilizing the voltage u obtained in the step oneαSubtracting the high-frequency sliding mode signal s obtained in the step threeαObtain a first set of intermediate variables EUsing the voltage u obtained in step oneβSubtracting the high-frequency sliding mode signal s obtained in the step threeβObtain a first set of intermediate variables E
Figure FDA0002564304500000016
Step five, obtaining a first group of intermediate variables E according to the step four、EAnd the q-axis inductance of the motor is calculated to obtain a second group of intermediate variables E、E
Figure FDA0002564304500000017
Wherein L isqIs a motor q-axis inductor;
step six, according to the current in the step three
Figure FDA0002564304500000018
Electric current
Figure FDA0002564304500000019
Calculating the resistance of the stator of the motor and the q-axis inductance to obtain a third group of intermediate variables E、E
Figure FDA00025643045000000110
Wherein R issIs a motor stator resistor;
step seven, utilizing the second group of intermediate variables E obtained in the step fiveSubtracting the third group intermediate variable E obtained in the step sixObtaining a fourth group of intermediate variables EUsing the second set of intermediate variables E obtained in step fiveSubtracting the third group intermediate variable E obtained in the step sixObtaining a fourth group of intermediate variables E
Figure FDA00025643045000000111
Step eight, obtaining a fourth group intermediate variable E according to the step seven、EUpdating the current in step three
Figure FDA00025643045000000112
And current
Figure FDA00025643045000000113
Figure FDA00025643045000000114
Further updating high-frequency sliding mode signal sαAnd high frequency sliding mode signal sβWherein, TsIs a sampling period;
step nine, updating the high-frequency sliding mode signal sαAnd high frequency sliding mode signal sβLow-pass filtering through the first low-pass filter to respectively obtain a fifth group of intermediate variables s、s
Figure FDA0002564304500000021
Wherein, ω iscIs the cut-off frequency of the first low-pass filter and the second low-pass filter, s is the laplacian operator;
step ten, the fifth group of intermediate variables s obtained in the step nine、sLow-pass filtering through a second low-pass filter to respectively obtain a sixth group of intermediate variables s、s
Figure FDA0002564304500000022
Step eleven, obtaining a fifth group of intermediate variables s according to the step nine、sCalculating effective back electromotive force q-axis deviation
Figure FDA0002564304500000023
Figure FDA0002564304500000024
Wherein, theta1Is the rotor position angle;
step twelve, the effective counter electromotive force q-axis deviation obtained in the step eleven
Figure FDA0002564304500000025
Calculating to obtain estimated rotating speed by a proportional integrator
Figure FDA0002564304500000026
Figure FDA0002564304500000027
kpIs a proportionality coefficient, kiIs an integral coefficient;
thirteen step, rotating speed obtained in the step twelve
Figure FDA0002564304500000028
The rotor position angle theta after primary phase delay generated by the first low-pass filter is obtained through integral adjustment1
Figure FDA0002564304500000029
Step fourteen, obtaining a sixth group of intermediate variables s according to the step ten、sCalculating the rotor position angle theta after two phase delays generated by the first low-pass filter and the second low-pass filter2:θ2=arctan(-s/s);
Step fifteen, obtaining the rotor position angle theta according to the step thirteen1Subtracting the rotor position angle theta obtained in the step fourteen2The retardation angle △ theta, △ theta and theta are obtained12According to the rotor position angle theta1And delay angle △ theta to obtain an estimated actual rotor position angle
Figure FDA00025643045000000210
Figure FDA00025643045000000211
Sixthly, obtaining a fifth group of intermediate variables s according to the step nine、sAnd a sixth set of intermediate variables s obtained in step ten、sCalculating to obtain an attenuation ratio k of the first low-pass filter and the second low-pass filter to the back electromotive force amplitude of the motor:
Figure FDA00025643045000000212
seventhly, obtaining a fifth group of intermediate variables s according to the step nine、sAnd calculating the attenuation ratio k obtained in the step sixteen to obtain the back electromotive force amplitude e of the motorm
Figure FDA00025643045000000213
Eighteen, obtaining the counter electromotive force amplitude e according to the seventeenth stepmAnd the actual rotor position angle obtained in step fifteen
Figure FDA0002564304500000037
Calculating to obtain estimated motor back electromotive force
Figure FDA0002564304500000031
And
Figure FDA0002564304500000032
Figure FDA0002564304500000033
2. method for observing back electromotive force of permanent magnet synchronous motor according to claim 1, wherein the three-phase voltage u in the first stepA、uBAnd uCObtaining the voltage u under a static α - β coordinate system through coordinate transformationαAnd voltage uβThe method comprises the following steps:
Figure FDA0002564304500000034
wherein,
Figure FDA0002564304500000035
uABand uBCIs the motor stator voltage.
3. Permanent magnet synchronous machine back-emf according to claim 1The method for observing the dynamic force is characterized in that the motor stator current i in the step twoA、iBAnd iCObtaining the current i under a stationary α - β coordinate system through coordinate transformationαAnd current iβThe method comprises the following steps:
Figure FDA0002564304500000036
CN201910376079.0A 2019-05-07 2019-05-07 Counter electromotive force observation method for permanent magnet synchronous motor Expired - Fee Related CN110112965B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910376079.0A CN110112965B (en) 2019-05-07 2019-05-07 Counter electromotive force observation method for permanent magnet synchronous motor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910376079.0A CN110112965B (en) 2019-05-07 2019-05-07 Counter electromotive force observation method for permanent magnet synchronous motor

Publications (2)

Publication Number Publication Date
CN110112965A CN110112965A (en) 2019-08-09
CN110112965B true CN110112965B (en) 2020-09-04

Family

ID=67488625

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910376079.0A Expired - Fee Related CN110112965B (en) 2019-05-07 2019-05-07 Counter electromotive force observation method for permanent magnet synchronous motor

Country Status (1)

Country Link
CN (1) CN110112965B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110995072B (en) * 2019-12-19 2021-09-10 华中科技大学 Motor rotor position estimation method
CN111969923B (en) * 2020-07-20 2022-04-19 四川虹美智能科技有限公司 Method and device for determining rotating speed of motor and electronic equipment
CN112325442A (en) * 2020-11-05 2021-02-05 广东美的暖通设备有限公司 Motor locked-rotor detection method, air conditioner and readable storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7256564B2 (en) * 2005-09-29 2007-08-14 Agile Systems Inc. System and method for attenuating noise associated with a back electromotive force signal in a motor
CN103904973A (en) * 2014-03-24 2014-07-02 合肥工业大学 Method for controlling salient pole permanent magnet synchronous machine in speed sensorless mode
CN107222133A (en) * 2017-05-18 2017-09-29 哈尔滨工程大学 A kind of method for the suppression brshless DC motor commutation torque ripple that can compensate Hall deviation automatically
CN107342713A (en) * 2017-07-17 2017-11-10 北京航空航天大学 A kind of permanent-magnet synchronous motor rotor position detection means and method based on improvement sliding mode observer

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7429841B2 (en) * 2003-10-02 2008-09-30 The Berquist Torrington Company Method and apparatus for movable element position detection in an electrically commutated machine

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7256564B2 (en) * 2005-09-29 2007-08-14 Agile Systems Inc. System and method for attenuating noise associated with a back electromotive force signal in a motor
CN103904973A (en) * 2014-03-24 2014-07-02 合肥工业大学 Method for controlling salient pole permanent magnet synchronous machine in speed sensorless mode
CN107222133A (en) * 2017-05-18 2017-09-29 哈尔滨工程大学 A kind of method for the suppression brshless DC motor commutation torque ripple that can compensate Hall deviation automatically
CN107342713A (en) * 2017-07-17 2017-11-10 北京航空航天大学 A kind of permanent-magnet synchronous motor rotor position detection means and method based on improvement sliding mode observer

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Sensorless control of PMSM in high speed range with iterative sliding mode observer;Kye-Lyong Kang等;《 Nineteenth Annual IEEE Applied Power Electronics Conference and Exposition, 2004. APEC "04》;20040226(第2期);第1111-1116页 *
改进滑模观测器的永磁同步电机SVM-DTC系统;杨帆等;《电气自动化》;20150131;第37卷(第1期);第4-6页 *

Also Published As

Publication number Publication date
CN110112965A (en) 2019-08-09

Similar Documents

Publication Publication Date Title
CN110429886B (en) Permanent magnet synchronous motor low-speed domain rotor position identification method
CN105915136B (en) Current of electric harmonics restraint control system based on fuzzy neural network and method
CN110112965B (en) Counter electromotive force observation method for permanent magnet synchronous motor
CN109889117B (en) IPMSM position observation method, system and driving system based on rotation high-frequency injection method
CN108631672B (en) Permanent magnet synchronous motor prediction flux linkage control method considering optimal duty ratio modulation
CN110022106B (en) Permanent magnet synchronous motor position sensorless control method based on high-frequency signal injection
CN105827168A (en) PMSM control method and system based on sliding mode observation
CN105915142B (en) A kind of permanent-magnet synchronous motor rotor position and turn count method based on decoupling self-adaptive observer
Wang et al. Improved fast method of initial rotor position estimation for interior permanent magnet synchronous motor by symmetric pulse voltage injection
CN110995102A (en) Direct torque control method and system for permanent magnet synchronous motor
CN111478636A (en) Single current sensor prediction control method for permanent magnet synchronous motor
CN108306569A (en) Based on the permanent magnet synchronous motor of Generalized Integrator without method for control speed and system
CN111211717B (en) IPMSM (intelligent power management system) position-sensorless motor closed-loop control method of nonsingular sliding mode structure
CN105262395A (en) Method and system for controlling permanent magnet synchronous motor based on sliding mode control theory
CN109713971B (en) Disturbance suppression method for permanent magnet synchronous motor
WO2020215588A1 (en) Function-switching-based position-sensorless control method and system for direct-current magnetism-regulating memory machine
CN112422002A (en) Robust permanent magnet synchronous motor single current sensor prediction control method
CN109560741B (en) Permanent magnet synchronous motor system based on measurement error compensator and compensation method
CN115566947A (en) Improved sliding-mode non-inductive six-phase permanent magnet synchronous motor control method based on ESO
Parasiliti et al. Torque ripple compensation in permanent magnet synchronous motors based on Kalman filter
CN106849803B (en) Permanent Magnet Synchronous Motor estimation method is filtered based on edge particles are uniformly distributed
CN110649847A (en) Position-sensorless control method of PMSLM (permanent magnet synchronous Motor) at low-speed stage
CN112928959B (en) Permanent magnet synchronous motor position sensorless control method
De Kock et al. Dynamic control of the permanent magnet-assisted reluctance synchronous machine
CN113965129A (en) Compensation method for current measurement offset error of permanent magnet synchronous motor control system

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200904