CN113098350B - Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer - Google Patents

Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer Download PDF

Info

Publication number
CN113098350B
CN113098350B CN202110563565.0A CN202110563565A CN113098350B CN 113098350 B CN113098350 B CN 113098350B CN 202110563565 A CN202110563565 A CN 202110563565A CN 113098350 B CN113098350 B CN 113098350B
Authority
CN
China
Prior art keywords
value
alpha
electromotive force
observation
beta axis
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
Application number
CN202110563565.0A
Other languages
Chinese (zh)
Other versions
CN113098350A (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.)
Wuxi Weifu High Technology Group Co Ltd
Original Assignee
Wuxi Weifu High Technology Group Co Ltd
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 Wuxi Weifu High Technology Group Co Ltd filed Critical Wuxi Weifu High Technology Group Co Ltd
Priority to CN202110563565.0A priority Critical patent/CN113098350B/en
Publication of CN113098350A publication Critical patent/CN113098350A/en
Application granted granted Critical
Publication of CN113098350B publication Critical patent/CN113098350B/en
Active 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
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • 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
    • H02P2203/00Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/03Determination of the rotor position, e.g. initial rotor position, during standstill or low speed operation
    • 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
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/72Electric energy management in electromobility

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 provides a method and a device for calculating the rotor position of a surface-mounted permanent magnet synchronous motor based on a sliding-mode observer, wherein a cosine saturation piecewise function of a self-adaptive boundary layer thickness is introduced into a sliding-mode control rate, and a quasi-per-unit algorithm is combined, so that high-frequency buffeting is effectively reduced, and the calculation accuracy of the rotor position and the rotating speed is improved; firstly, calculating by a sliding mode observer according to a stator voltage reference value, a stator current actual value and a rotor electrical angular velocity observed value to obtain a back electromotive force sliding mode observation standard per unit value; calculating by a second-order generalized integrator according to the back electromotive force sliding mode observation quasi-per-unit value and the rotor electrical angular velocity observation value to obtain a back electromotive force observation quasi-per-unit value; finally, calculating through a phase-locked loop according to the back electromotive force observation quasi-per-unit value to obtain a rotor electrical angular velocity observation value and a rotor electrical angle observation value; the rotor position calculation method provided by the invention has high reliability and strong robustness, and meets the requirements of the driving field of permanent magnet synchronous motors such as high-speed hydrogen pumps, air compressors and the like on the reliability and efficiency of the system.

Description

Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer
Technical Field
The invention relates to the technical field of motor control, in particular to a sliding-mode observer-based permanent magnet synchronous motor rotor position calculation method.
Background
The permanent magnet synchronous motor has the characteristics of high power density, high efficiency, compact structure, high reliability and the like, and is widely applied to the fields of aircrafts, wind power generation, household appliances, electric automobiles and the like. In order to realize high-performance control of the permanent magnet synchronous motor and meet the requirement of accurate control of rotating speed and current, a motor driving system must accurately acquire the position of a rotor in real time.
The traditional mechanical sensor for direct position detection mainly comprises a rotary transformer, an encoder, a Hall sensor and the like, so that the cost of the motor is increased, the size and the weight of the motor are increased, the motor is easily influenced by the working environment, and the reliability of the system is reduced.
The current technology for estimating the position of the rotor without a position sensor comprises two main categories: based on non-ideal characteristics of the motor and based on a fundamental wave mathematical model. The rotor position can be accurately estimated under the conditions of zero speed and extremely low speed of the motor based on the motor nonideal characteristic technology, the rotor position mainly comprises a high-frequency signal injection method, a low-frequency injection method and the like, but the signal injection can cause extra harmonic waves and noise to the motor, the calculation amount is large, and the voltage utilization rate is low. The method is suitable for medium-high speed operation working conditions and mainly comprises a direct formula calculation method, a model reference self-adaption method, an extended Kalman filtering method, a sliding-mode observer method and the like.
The sliding-mode observer is insensitive to system parameters, strong in robustness to external parameter disturbance and internal parameter perturbation, simple in structural algorithm, high in engineering control reliability and easy to achieve. However, the control of the sliding-mode observer is discontinuous control, and the high-frequency buffeting exists in the system due to the non-ideality of a switch and the inertia of the system, so that the estimation accuracy of the position of the rotor is greatly reduced. In addition, the change of the rotating speed can influence the amplitude of the back electromotive force, thereby influencing the bandwidth of a system, reducing the anti-rotating speed disturbance capability of fixed parameters and influencing the observation precision. Therefore, the sliding mode observer sensorless algorithm capable of realizing high-frequency buffeting suppression and improvement of the anti-rotation speed disturbance capacity has wide application prospect.
Disclosure of Invention
In order to solve the technical problem, the invention provides a surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on a sliding-mode observer; according to the method, a cosine saturation piecewise function of the adaptive boundary layer thickness is introduced into the sliding mode control rate, and a quasi-unitary algorithm is combined, so that high-frequency buffeting is effectively reduced, the calculation accuracy of the position and the rotating speed of the rotor is improved, and the reliability and the robustness of the method for calculating the position of the rotor by the sliding mode observer without sensing are enhanced.
According to the method, a cosine saturation piecewise function of the self-adaptive boundary layer thickness is introduced into the sliding mode control rate, and a quasi-per-unit algorithm is combined, so that high-frequency buffeting is effectively reduced, and the calculation accuracy of the position and the rotating speed of the rotor is improved; firstly, calculating by a sliding mode observer according to a stator voltage reference value, a stator current actual value and a rotor electrical angular velocity observed value to obtain a back electromotive force sliding mode observation standard per unit value; calculating by a second-order generalized integrator according to the back electromotive force sliding mode observation quasi-per-unit value and the rotor electrical angular velocity observation value to obtain a back electromotive force observation quasi-per-unit value; finally, calculating through a phase-locked loop according to the back electromotive force observation quasi-per-unit value to obtain a rotor electrical angular velocity observation value and a rotor electrical angle observation value; the rotor position calculation method provided by the invention has high reliability and strong robustness, and meets the requirements of the driving field of permanent magnet synchronous motors such as high-speed hydrogen pumps, air compressors and the like on the reliability and efficiency of the system.
As a first aspect of the present invention, an embodiment of the present invention provides a surface-mount permanent magnet synchronous motor rotor position calculation method based on a sliding-mode observer, including the following steps:
step S1, according to the reference value of the stator voltage of the alpha-beta axis
Figure BDA0003079900570000021
Actual value i of stator current of alpha and beta axisα/iβObserved value of rotor electrical angular velocity
Figure BDA0003079900570000022
Calculating by a sliding mode observer to obtain an alpha-beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure BDA0003079900570000023
Step S2, observing quasi-per-unit value according to alpha and beta axis back electromotive force sliding mode
Figure BDA0003079900570000024
And rotor electrical angular velocity observations
Figure BDA0003079900570000025
Calculating to obtain alpha beta axis back electromotive force observation quasi-per-unit value through a second-order generalized integrator
Figure BDA0003079900570000026
Step S3, observing the standard mark according to the alpha beta axis counter electromotive forcePer value
Figure BDA0003079900570000027
Obtaining the observed value of the rotor electrical angular velocity by phase-locked loop calculation
Figure BDA0003079900570000028
And rotor electrical angle observation value
Figure BDA0003079900570000029
Further, in step S1, an α β axis back electromotive force sliding mode observation quasi-per-unit value is calculated
Figure BDA00030799005700000210
The method comprises the following steps:
step S1.1, according to the feedback obtained alpha beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure BDA00030799005700000211
Alpha beta axis stator voltage reference value
Figure BDA00030799005700000212
And the actual value i of the stator current of the alpha and beta axisα/iβObtaining an alpha beta axis back electromotive force sliding mode observed value through reverse quasi-unitary calculation
Figure BDA00030799005700000213
As shown in formula (1) - (2):
Figure BDA00030799005700000214
wherein,
Figure BDA00030799005700000215
where λ is the adaptive sliding mode gain, τ is a non-zero smaller positive tuning parameter, eαeq/eβeqRespectively alpha and beta axes of counter currentCalculated value of the kinetic equation, | eeqI is the calculated amplitude of the back electromotive force formula, RsIs stator resistance, LsIs a stator inductance;
step S1.2, according to the reference value of the stator voltage of the alpha and beta axis
Figure BDA00030799005700000216
And the alpha beta axis back electromotive force sliding mode observed value obtained by feedback
Figure BDA00030799005700000217
Calculating to obtain an alpha beta axis stator current observed value through a formula (3)
Figure BDA00030799005700000218
The following:
Figure BDA00030799005700000219
step S1.3, according to the stator current observed value of the alpha-beta axis
Figure BDA00030799005700000220
And the actual value i of the stator current of the alpha and beta axisα/iβThe sliding mode surface s is defined by equation (4) as follows:
Figure BDA0003079900570000031
in the formula,
Figure BDA0003079900570000032
respectively alpha and beta axis stator current observation error values;
step S1.4, observing error value according to alpha and beta axis stator current
Figure BDA0003079900570000033
And rotor electrical angular velocity observations
Figure BDA0003079900570000034
By self-runningCalculating a cosine saturation piecewise function adaptive to the thickness of the boundary layer to obtain an alpha beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure BDA0003079900570000035
As shown in equations (5) - (7):
Figure BDA0003079900570000036
wherein,
Figure BDA0003079900570000037
then the process of the first step is carried out,
Figure BDA0003079900570000038
where ε is the first boundary layer, δ is the second boundary layer, δ0Is the second boundary layer at rated rotor electrical angular velocity, ω0For nominal rotor electrical angular velocity, u(s) is a cosine saturation piecewise function.
Further, in step S2, an α β axis back electromotive force observation quasi-per-unit value is calculated
Figure BDA0003079900570000039
The method of (2) is shown in formula (8) - (10):
Figure BDA00030799005700000310
Figure BDA00030799005700000311
Figure BDA00030799005700000312
in the formula,kzIn order to adjust the parameters for the filtering,
Figure BDA00030799005700000313
respectively are primary filtering observation quasi-per-unit values of alpha and beta axis counter electromotive force,
Figure BDA00030799005700000314
respectively, alpha beta axis back electromotive force secondary filtering observation quasi-per-unit values.
Further, in step S3, an observed value of the rotor electrical angular velocity is calculated
Figure BDA00030799005700000315
And rotor electrical angle observation value
Figure BDA00030799005700000316
The method comprises the following steps:
step S3.1, according to the alpha beta axis counter electromotive force observed value
Figure BDA00030799005700000317
And feeding back the obtained rotor electrical angle observed value
Figure BDA00030799005700000318
Calculating to obtain an observation deviation value delta theta of the rotor electrical angle through a formula (11)eThe following were used:
Figure BDA0003079900570000041
in the formula, thetaeThe actual value of the rotor electrical angle is obtained;
s3.2, observing the deviation value delta theta according to the rotor electrical angleeObtaining the observed value of the rotor electrical angular velocity through the calculation of the formula (12) - (13)
Figure BDA0003079900570000042
And rotor electrical angle observation value
Figure BDA0003079900570000043
The following were used:
Figure BDA0003079900570000044
Figure BDA0003079900570000045
in the formula, KpFor proportional adjustment of parameters for PI, KiThe parameters are adjusted for PI integration.
As a second aspect of the present invention, an embodiment of the present invention provides a surface-mount permanent magnet synchronous motor rotor position calculation apparatus based on a sliding-mode observer, including:
a memory storing a computer program;
a processor for executing the computer program, the computer program when executed performing the steps of the method as described above.
The invention has the following advantages:
(1) the cosine saturation piecewise function of the self-adaptive boundary layer thickness is introduced into the sliding mode control rate, so that the high-frequency buffeting of the system caused by the discontinuity of a sliding mode surface change-over switch and the unstable transition near a boundary layer is weakened while the response speed and the robustness of the original system are not influenced;
(2) the back electromotive force observation value is subjected to quasi-per-unit processing, the influence of system bandwidth change caused by back electromotive force amplitude change is eliminated, the anti-rotation speed disturbance capability of fixed adjustment parameters is enhanced, and the rotor position and rotation speed calculation accuracy is improved.
Drawings
Fig. 1 is a block diagram of a control system in an embodiment of the present invention.
Fig. 2 is a schematic block diagram of a sliding-mode observer in the embodiment of the present invention.
Fig. 3 is a schematic diagram of a cosine saturation piecewise function in the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Firstly, the embodiment of the invention provides a surface-mounted permanent magnet synchronous motor rotor position calculation method based on a sliding-mode observer, which comprises the following steps:
step S1, according to the reference value of the stator voltage of the alpha-beta axis
Figure BDA0003079900570000046
Actual value i of stator current of alpha-beta axisα/iβObserved value of rotor electrical angular velocity
Figure BDA0003079900570000047
Calculating by a sliding mode observer to obtain an alpha-beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure BDA0003079900570000048
Calculating alpha beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure BDA0003079900570000049
The method comprises the following steps:
step S1.1, according to the feedback obtained alpha beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure BDA00030799005700000410
Alpha beta axis stator voltage reference value
Figure BDA00030799005700000411
And the actual value i of the stator current of the alpha and beta axisα/iβObtaining an alpha beta axis back electromotive force sliding mode observed value through reverse quasi-unitary calculation
Figure BDA0003079900570000051
As shown in formula (1) - (2):
Figure BDA0003079900570000052
wherein,
Figure BDA0003079900570000053
where λ is the adaptive sliding mode gain, τ is a non-zero smaller positive tuning parameter, eαeq/eβeqRespectively, alpha beta axis back electromotive force formula calculation value, | eeqI is the calculated amplitude of the back electromotive force formula, RsIs stator resistance, LsIs a stator inductance;
step S1.2, according to the reference value of the stator voltage of the alpha and beta axis
Figure BDA0003079900570000054
And the alpha beta axis back electromotive force sliding mode observed value obtained by feedback
Figure BDA0003079900570000055
Calculating to obtain an alpha beta axis stator current observed value through a formula (3)
Figure BDA0003079900570000056
The following were used:
Figure BDA0003079900570000057
step S1.3, according to the stator current observed value of the alpha-beta axis
Figure BDA0003079900570000058
And the actual value i of the stator current of the alpha and beta axisα/iβThe sliding mode surface s is defined by equation (4) as follows:
Figure BDA0003079900570000059
in the formula,
Figure BDA00030799005700000510
respectively alpha and beta axis stator current observation error values;
step S1.4, observing error value according to alpha and beta axis stator current
Figure BDA00030799005700000511
And rotor electrical angular velocity observations
Figure BDA00030799005700000512
Calculating to obtain an alpha beta axis back electromotive force sliding mode observation quasi-per unit value through a cosine saturation piecewise function of the self-adaptive boundary layer thickness
Figure BDA00030799005700000513
As shown in equations (5) - (7):
Figure BDA00030799005700000514
wherein,
Figure BDA00030799005700000515
then the process of the first step is carried out,
Figure BDA00030799005700000516
where ε is the first boundary layer, δ is the second boundary layer, δ0Is the second boundary layer at rated rotor electrical angular velocity, ω0For nominal rotor electrical angular velocity, u(s) is a cosine saturation piecewise function;
step S2, observing quasi-per-unit value according to alpha and beta axis back electromotive force sliding mode
Figure BDA0003079900570000061
And rotor electrical angular velocity observations
Figure BDA0003079900570000062
Alpha-beta axis inverse is calculated by a second-order generalized integratorStandard per unit value for electromotive force observation
Figure BDA0003079900570000063
Calculating alpha beta axis back electromotive force observation quasi-per-unit value
Figure BDA0003079900570000064
The method of (2) is shown in formula (8) - (10):
Figure BDA0003079900570000065
Figure BDA0003079900570000066
Figure BDA0003079900570000067
in the formula, kzIn order to adjust the parameters for the filtering,
Figure BDA0003079900570000068
respectively alpha beta axis back electromotive force primary filtering observation quasi-per-unit values,
Figure BDA0003079900570000069
respectively are alpha beta axis counter electromotive force secondary filtering observation quasi-per-unit values;
step S3, observing the quasi-per unit value according to the alpha beta axis back electromotive force
Figure BDA00030799005700000610
Obtaining the observed value of the rotor electrical angular velocity through phase-locked loop calculation
Figure BDA00030799005700000611
And rotor electrical angle observation value
Figure BDA00030799005700000612
Calculating observed value of rotor electrical angular velocity
Figure BDA00030799005700000613
And rotor electrical angle observation value
Figure BDA00030799005700000614
The method comprises the following steps:
step S3.1, according to the alpha beta axis counter electromotive force observed value
Figure BDA00030799005700000615
And feeding back the obtained rotor electrical angle observed value
Figure BDA00030799005700000616
Calculating to obtain an observation deviation value delta theta of the rotor electrical angle through a formula (11)eThe following were used:
Figure BDA00030799005700000617
in the formula, thetaeThe actual value of the rotor electrical angle is obtained;
s3.2, observing the deviation value delta theta according to the rotor electrical angleeObtaining the observed value of the rotor electrical angular velocity through the calculation of the formula (12) - (13)
Figure BDA00030799005700000618
And rotor electrical angle observation value
Figure BDA00030799005700000619
The following were used:
Figure BDA00030799005700000620
Figure BDA00030799005700000621
in the formula, KpFor proportional adjustment of parameters for PI, KiThe parameters are adjusted for PI integration.
Secondly, an embodiment of the present invention further provides a sliding-mode observer-based surface-mounted permanent magnet synchronous motor rotor position calculation apparatus, including: a processor and a memory; the processor and the memory communicate with each other, for example, by being connected to and communicating with each other via a communication bus; the memory has stored therein a computer program; the processor is adapted to run the computer program, which when run performs the steps of the method as described above.
Finally, it should be noted that the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, and although the present invention has been described in detail with reference to examples, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, which should be covered by the claims of the present invention.

Claims (4)

1. A surface-mounted permanent magnet synchronous motor rotor position calculation method based on a sliding-mode observer is characterized by comprising the following steps:
step S1, according to the reference value of the stator voltage of the alpha-beta axis
Figure FDA0003543330940000011
Actual value i of stator current of alpha and beta axisα/iβObserved value of rotor electrical angular velocity
Figure FDA0003543330940000012
Calculating by a sliding mode observer to obtain an alpha-beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure FDA0003543330940000013
Step S2, observing quasi-per-unit value according to alpha and beta axis back electromotive force sliding mode
Figure FDA0003543330940000014
And rotor electrical angular velocity observations
Figure FDA0003543330940000015
Calculating to obtain alpha beta axis back electromotive force observation quasi-per-unit value through a second-order generalized integrator
Figure FDA0003543330940000016
Step S3, observing quasi-per-unit value according to alpha and beta axis back electromotive force
Figure FDA0003543330940000017
Obtaining the observed value of the rotor electrical angular velocity through phase-locked loop calculation
Figure FDA0003543330940000018
And rotor electrical angle observation value
Figure FDA0003543330940000019
In step S1, an α β axis back electromotive force sliding mode observation quasi-per-unit value is calculated
Figure FDA00035433309400000110
The method comprises the following steps:
step S1.1, according to the feedback obtained alpha beta axis back electromotive force sliding mode observation quasi-per-unit value
Figure FDA00035433309400000111
Alpha beta axis stator voltage reference value
Figure FDA00035433309400000112
And the actual value i of the stator current of the alpha and beta axisα/iβObtaining an alpha beta axis back electromotive force sliding mode observed value through reverse quasi-unitary calculation
Figure FDA00035433309400000113
As shown in formula (1) - (2):
Figure FDA00035433309400000114
wherein,
Figure FDA00035433309400000115
where λ is the adaptive sliding mode gain, τ is a non-zero smaller positive tuning parameter, eαeq/eβeqRespectively, alpha beta axis back electromotive force formula calculation value, | eeqI is the calculated amplitude of the back electromotive force formula, RsIs stator resistance, LsA stator inductor;
step S1.2, according to the reference value of the alpha and beta axis stator voltage
Figure FDA00035433309400000116
And the alpha beta axis back electromotive force sliding mode observed value obtained by feedback
Figure FDA00035433309400000117
Calculating to obtain an alpha beta axis stator current observed value through a formula (3)
Figure FDA00035433309400000118
The following were used:
Figure FDA00035433309400000119
step S1.3, according to the stator current observed value of the alpha-beta axis
Figure FDA00035433309400000120
And the actual value i of the stator current of the alpha and beta axisα/iβThe sliding mode surface s is defined by equation (4) as follows:
Figure FDA00035433309400000121
in the formula,
Figure FDA0003543330940000021
respectively alpha and beta axis stator current observation error values;
step S1.4, observing error value according to alpha and beta axis stator current
Figure FDA0003543330940000022
And rotor electrical angular velocity observations
Figure FDA0003543330940000023
Calculating to obtain alpha beta axis back electromotive force sliding mode observation quasi-per-unit value through cosine saturation piecewise function of self-adaptive boundary layer thickness
Figure FDA0003543330940000024
As shown in equations (5) - (7):
Figure FDA0003543330940000025
wherein,
Figure FDA0003543330940000026
then the user can use the device to make a visual display,
Figure FDA0003543330940000027
where ε is the first boundary layer, δ is the second boundary layer, δ0Is the second boundary layer at rated rotor electrical angular velocity, ω0For nominal rotor electrical angular velocity, u(s) is a cosine saturation piecewise function.
2. The method for calculating the rotor position of the surface-mounted permanent magnet synchronous motor based on the sliding-mode observer according to claim 1,
in step S2, an α β axis back electromotive force observation quasi-per-unit value is calculated
Figure FDA0003543330940000028
The method (2) is shown in the formula (8) - (10):
Figure FDA0003543330940000029
Figure FDA00035433309400000210
Figure FDA00035433309400000211
in the formula, kzIn order to adjust the parameters for the filtering,
Figure FDA00035433309400000212
respectively alpha beta axis back electromotive force primary filtering observation quasi-per-unit values,
Figure FDA00035433309400000213
respectively, alpha beta axis back electromotive force secondary filtering observation quasi-per-unit values.
3. The method for calculating the rotor position of the surface-mounted permanent magnet synchronous motor based on the sliding-mode observer according to claim 2,
in step S3, an observed value of the rotor electrical angular velocity is calculated
Figure FDA00035433309400000214
And rotor electrical angle observation value
Figure FDA00035433309400000215
The method comprises the following steps:
step S3.1, according to the alpha beta axis counter electromotive force observed value
Figure FDA0003543330940000031
And feeding back the obtained rotor electrical angle observed value
Figure FDA0003543330940000032
Calculating to obtain an observation deviation value delta theta of the rotor electrical angle through a formula (11)eThe following were used:
Figure FDA0003543330940000033
in the formula, thetaeThe actual value of the rotor electrical angle is obtained;
s3.2, observing the deviation value delta theta according to the rotor electrical angleeObtaining the observed value of the rotor electrical angular velocity through the calculation of the formula (12) - (13)
Figure FDA0003543330940000034
And rotor electrical angle observation value
Figure FDA0003543330940000035
The following:
Figure FDA0003543330940000036
Figure FDA0003543330940000037
in the formula, KpFor proportional adjustment of parameters for PI, KiThe parameters are adjusted for PI integration.
4. A surface-mounted permanent magnet synchronous motor rotor position calculating device based on a sliding-mode observer is characterized by comprising:
a memory storing a computer program;
a processor for running the computer program, the computer program when running performing the steps of the method of any one of claims 1 to 3.
CN202110563565.0A 2021-05-24 2021-05-24 Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer Active CN113098350B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110563565.0A CN113098350B (en) 2021-05-24 2021-05-24 Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110563565.0A CN113098350B (en) 2021-05-24 2021-05-24 Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer

Publications (2)

Publication Number Publication Date
CN113098350A CN113098350A (en) 2021-07-09
CN113098350B true CN113098350B (en) 2022-05-31

Family

ID=76666038

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110563565.0A Active CN113098350B (en) 2021-05-24 2021-05-24 Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer

Country Status (1)

Country Link
CN (1) CN113098350B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103199779A (en) * 2013-04-22 2013-07-10 哈尔滨工业大学 Position observation device and method for rotor of built-in permanent magnetic synchronous motor based on adaptive filtering
CN105227010A (en) * 2015-10-23 2016-01-06 哈尔滨工业大学 A kind of permagnetic synchronous motor position-sensor-free position detection error harmonic pulse removing method
CN106487304A (en) * 2016-10-27 2017-03-08 江苏大学 A kind of permagnetic synchronous motor method for estimating state based on sliding formwork back-EMF observer device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100237817A1 (en) * 2009-03-23 2010-09-23 Jingbo Liu Method and Apparatus for Estimating Rotor Position in a Sensorless Synchronous Motor

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103199779A (en) * 2013-04-22 2013-07-10 哈尔滨工业大学 Position observation device and method for rotor of built-in permanent magnetic synchronous motor based on adaptive filtering
CN105227010A (en) * 2015-10-23 2016-01-06 哈尔滨工业大学 A kind of permagnetic synchronous motor position-sensor-free position detection error harmonic pulse removing method
CN106487304A (en) * 2016-10-27 2017-03-08 江苏大学 A kind of permagnetic synchronous motor method for estimating state based on sliding formwork back-EMF observer device

Also Published As

Publication number Publication date
CN113098350A (en) 2021-07-09

Similar Documents

Publication Publication Date Title
CN110198150B (en) Permanent magnet synchronous motor multi-parameter online identification method
CN103560738B (en) A kind of permanent-magnet synchronous motor rotor position observation system and observation procedure thereof suppressing position microseismic observation error
CN110311608B (en) High-frequency square wave voltage injection permanent magnet synchronous motor position-sensorless control method with optimal injection angle
WO2021174582A1 (en) Position sensorless permanent magnet synchronous motor control method and automobile power system
CN110022106B (en) Permanent magnet synchronous motor position sensorless control method based on high-frequency signal injection
CN111193448B (en) Permanent magnet synchronous motor load torque observation method based on extended Kalman filter
CN105227010A (en) A kind of permagnetic synchronous motor position-sensor-free position detection error harmonic pulse removing method
CN111181458A (en) Surface-mounted permanent magnet synchronous motor rotor flux linkage observation method based on extended Kalman filter
CN108306570A (en) Direct torque control method for permanent magnetic synchronous electric machine and system
CN108551285A (en) Direct Torque Control System for Permanent Magnet Synchronous Motor and method based on double synovial membrane structures
CN114598206A (en) Design method of permanent magnet synchronous motor wide-speed-domain rotor position observer
CN111769779A (en) PMSM direct torque control method based on improved Luenberger observer
CN109510539A (en) One kind predicting magnetic linkage control system and method based on novel gain matrix norm type
Wang et al. Rotor position estimation method for permanent magnet synchronous motor based on super-twisting sliding mode observer
CN108599661A (en) A kind of permanent magnet synchronous motor is without sensor composite control method
CN110768600A (en) PMSM speed sensorless rotor detection method
CN118157526A (en) Permanent magnet synchronous motor control method based on improved linear supercoiled
CN113098350B (en) Surface-mounted permanent magnet synchronous motor rotor position calculation method and device based on sliding-mode observer
CN116317747B (en) Full-rotation-speed range tracking method for ultra-high-speed permanent magnet synchronous motor
CN114157193B (en) Optimization interpolation type synchronous motor torque pulsation suppression control method and system
CN114024382B (en) Permanent magnet synchronous motor single-rotating-speed-loop speed regulating device based on load torque feedback
Yang et al. Research about the Sensorless Vector Control of Permanent Magnet Synchronous Motor Based on Two-stage Filter Sliding Mode Observer
CN113179068B (en) Permanent magnet synchronous motor position sensorless control method based on improved square wave injection
CN116032169B (en) Torque ripple suppression method for self-adaptive anti-interference control high-speed permanent magnet synchronous motor
Junfeng et al. Research on PMSM Sensorless Control Based on Improved Synovial Observer Algorithm

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