CN114006557A - Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding-mode observer - Google Patents

Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding-mode observer Download PDF

Info

Publication number
CN114006557A
CN114006557A CN202111165630.0A CN202111165630A CN114006557A CN 114006557 A CN114006557 A CN 114006557A CN 202111165630 A CN202111165630 A CN 202111165630A CN 114006557 A CN114006557 A CN 114006557A
Authority
CN
China
Prior art keywords
equation
permanent magnet
extended
stator
mode observer
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.)
Granted
Application number
CN202111165630.0A
Other languages
Chinese (zh)
Other versions
CN114006557B (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.)
Hunan University of Science and Technology
Original Assignee
Hunan University of Science and Technology
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 Hunan University of Science and Technology filed Critical Hunan University of Science and Technology
Priority to CN202111165630.0A priority Critical patent/CN114006557B/en
Publication of CN114006557A publication Critical patent/CN114006557A/en
Application granted granted Critical
Publication of CN114006557B publication Critical patent/CN114006557B/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/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P21/0007Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using sliding mode control
    • 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/13Observer control, e.g. using Luenberger observers or Kalman filters
    • 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
    • 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/20Estimation of torque
    • 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/34Modelling or simulation for control purposes
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

The invention discloses a permanent magnet synchronous motor mechanical parameter identification method based on an extended sliding-mode observer, which comprises the following steps: establishing a mathematical model of the surface-mounted permanent magnet synchronous motor and adopting vector control through coordinate transformation; designing an extended sliding mode mechanical parameter view and determining related parameters of the observer; and observing system disturbance information by using the proposed extended sliding-mode observer, and extracting a damping coefficient, a rotational inertia and a load torque from the system disturbance information. Firstly, taking a friction coefficient parameter error, a rotational inertia parameter error and a load torque system disturbance as an extended system state, establishing an extended state equation for designing an extended sliding mode observer, wherein the extended sliding mode observer can track the system state in real time and soften a sliding mode variable structure control signal; therefore, the output of the extended sliding mode observer can be directly used for mechanical parameter estimation without generating phase lag, and the self-tuning control of the permanent magnet synchronous motor system is realized by correctly estimating the mechanical parameters in real time, so that the optimal control is realized.

Description

Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding-mode observer
Technical Field
The invention relates to the field of motor control, in particular to a permanent magnet synchronous motor mechanical parameter identification method based on an extended sliding-mode observer.
Background
The permanent magnet synchronous motor has the advantages of simple structure, small volume, high efficiency, high torque current ratio, high power density and the like, and compared with a speed regulating system formed by other motors, the alternating current permanent magnet speed regulating system formed by the permanent magnet synchronous motor has more excellent performance. The vector control technology is a control method commonly used by a permanent magnet synchronous motor system, and can realize the decoupling of the electromagnetic torque and the magnetic flux of the system through vector control, thereby avoiding the influence caused by the characteristics of nonlinearity and strong coupling of an alternating current motor, and enabling the system to have higher tracking capability, higher tracking precision and good robustness. However, for more complex working conditions, parameter changes of the system and external unknown disturbances, the dynamic performance of the system is reduced, and at the moment, the parameters of the controller should be adjusted in time to suppress the system disturbances caused by the parameter changes, so that the good dynamic performance and the steady-state performance of the system are ensured. The adjustment of the parameters of the controller can come from practical experience, but the adjustment method has no real-time property and cannot be applied to the working condition with frequent parameter change. Compared with the prior art, the parameter self-tuning technology of the controller can adjust the parameters on line, has high real-time performance, and the adjusting effect of the method depends on the accuracy of mechanical parameters such as load torque, rotational inertia and the like. Therefore, the system needs to include a mechanism for accurately identifying system parameters, and the controller parameters are adjusted in time according to the identification values, so that the influence of the disturbance on the system can be reduced or even eliminated.
The mechanical parameters of the permanent magnet synchronous motor system mainly refer to load torque, rotational inertia, damping coefficient and the like, wherein the identification method of the rotational inertia can be generally divided into two categories: an off-line identification method and an on-line identification method.
The off-line identification method comprises a power-off deceleration method, an orthogonal integration method, a single-wire torsional oscillation method and the like, but the methods cannot consider the parameter change of the motor in the actual engineering, are generally only used for initial value setting of motor control, and have the main defect of non-real-time property. The online identification eliminates the defect of offline identification by monitoring parameters in real time, and better meets the actual control requirement, and the negative influence is that a large amount of data needs to be stored and calculated in the working process, so that the requirement on a computer is higher and higher, and the identification accuracy is also influenced to a certain extent.
The online identification method comprises a state observer, extended Kalman filtering, model reference adaptive identification, least square identification and the like. The parameter identification method based on the model reference adaptive system is simple and feasible, and is applied to many practical engineering applications, but the method cannot be used for estimating the load torque in real time. The model-referenced adaptive method can estimate parameters by continuously adjusting adaptive law or function cost, but is difficult to be practically applied due to the complexity of practical implementation and the sensitivity of adaptive gain. The extended kalman filter algorithm takes the mechanical parameter as one of the system variables and takes it as a direct output of the extended kalman filter. The recursive least square algorithm iteratively solves the least square estimation when converging, and the dependence of the least square estimation on the initial condition is large. Although the recursive least squares method can be used to estimate the parameters of the permanent magnet synchronous motor control system, a longer estimation time is required in the estimation process. Therefore, this method is limited in scientific practice and practical industrial manufacturing. When load torque identification of the permanent magnet synchronous motor is carried out, the rotational inertia is defined as a known parameter by a disturbance observer, but in actual production practice, the two variables are generally unknown, and a good identification result is difficult to achieve in such a case.
Disclosure of Invention
In order to solve the technical problem, the invention provides a permanent magnet synchronous motor mechanical parameter identification method which can track the system state in real time, can be directly used for mechanical parameter estimation and does not generate phase lag.
The technical scheme for solving the problems is as follows: a permanent magnet synchronous motor mechanical parameter identification method based on an extended sliding-mode observer comprises the following steps:
1) establishing a mathematical model of the surface-mounted permanent magnet synchronous motor, and performing coordinate transformation by adopting vector control: firstly, establishing a mathematical model under an A-B-C three-phase static coordinate system to obtain a phase voltage equation, a flux linkage equation, a torque equation and a mechanical motion equation of each phase, and then obtaining the mathematical model under an alpha-beta two-phase static coordinate system and a d-q axis two-phase rotating coordinate system through Clark conversion and Park conversion in sequence to realize complete decoupling of exciting current and torque current so as to be convenient for control;
2) designing an extended sliding mode observer and determining observer related parameters: according to a mechanical motion equation of the PMSM, taking a damping coefficient parameter error, a rotational inertia parameter error and load torque system disturbance as an extended system state, establishing an extended state equation for designing an extended sliding mode observer, and ensuring the stability of the observer by considering a Lyapunov function;
3) observing system disturbance information by using the proposed extended sliding-mode observer, and extracting a damping coefficient, a rotational inertia and a load torque from the system disturbance information: firstly, the permanent magnet synchronous motor runs at two different steady-state speeds, the damping coefficient is estimated by utilizing the actually measured speed information and the estimated disturbance, then the rotary inertia is estimated by respectively running the permanent magnet synchronous motor at two different constant acceleration or constant deceleration, and finally, the damping coefficient and the rotary inertia are directly used for estimating the load torque.
The permanent magnet synchronous motor mechanical parameter identification method based on the extended sliding mode observer comprises the following specific processes in the step 1):
firstly, establishing a mathematical model under an A-B-C three-phase static coordinate system, wherein a voltage equation of the A-B-C three-phase static coordinate system is as follows:
Figure BDA0003291563990000031
wherein R issAs stator resistance,UA、UB、UCIs a stator three-phase voltage iA、iB、iCFor stator three-phase currents, ΨsA、ΨsB、ΨsCThe stator is a stator three-phase full magnetic linkage, and the vector form is as follows:
Figure BDA0003291563990000041
in the formula psisFor stator flux linkage, USIs the stator voltage; i.e. isIs the stator current;
the flux linkage equation is:
Figure BDA0003291563990000042
in the formula, ΨfA、ΨfB、ΨfCFlux linkage i for three-phase interlinkage of permanent magnet field and statorsA、isB、isCRespectively representing A, B, C three-phase stator currents, LsIn order to realize the synchronous inductance of the stator,
Figure BDA0003291563990000043
LA=LB=LC=L+Lmlwherein L isA、LB、LCFor self-induction of each phase, LmlExciting inductance for each phase; l isIs the leakage inductance of each phase;
the torque equation:
the electromagnetic torque is considered as a result of the interaction between the stator, rotor and armature, and is expressed as:
Figure BDA0003291563990000044
wherein p is the number of pole pairs, TeBeing electromagnetic torque, ΨfMagnetic linkage psi linking the permanent magnet field with the statorsA stator flux linkage;
mechanical equation of motion:
Figure BDA0003291563990000045
ωmas mechanical angular speed of the rotor, TLThe load torque is J, the moment of inertia is J, and the damping coefficient of the motor is B;
and finally obtaining a mathematical model under a d-q axis two-phase rotating coordinate system through Clark transformation and Park transformation to realize complete decoupling of exciting current and torque current so as to perform vector control:
the flux linkage equation:
Figure BDA0003291563990000051
in the formula, Ψd、ΨqIs the quadrature-direct component of the stator flux linkage, Ld、LqComponent of stator inductance in d, q axes, id、iqD and q axis stator currents;
voltage equation:
Figure BDA0003291563990000052
in the formula of Ud、UqThe quadrature-direct component, omega, of the stator voltageeRepresenting the rotor angular velocity;
the torque equation:
Figure BDA0003291563990000053
in the formula, Ld、LqThe inductance values are direct axis and quadrature axis inductance values; therefore, the decoupling of the torque is realized under the d-q axis, and the torque equation is as follows:
Figure BDA0003291563990000054
in the method for identifying mechanical parameters of the permanent magnet synchronous motor based on the extended sliding mode observer, in the step 2), a kinetic equation of the motor is expressed as follows:
Figure BDA0003291563990000055
wherein B is0、J0Is a rough estimation of real parameters, Δ B, Δ J are parameter errors between the real system and its rough estimation, and the moment of inertia J ═ J0+ Δ J, damping coefficient B ═ B0+ Δ B, ω is the rotational speed, d represents disturbances including parameter errors and load disturbances;
considering the disturbance d as an extended system state, in order to perform parameter estimation, the extended sliding mode observer is designed as follows:
Figure BDA0003291563990000061
wherein,
Figure BDA0003291563990000062
is an estimate of the disturbance,
Figure BDA0003291563990000063
is an estimate of the velocity, m is a sliding mode parameter, usmoRepresents the sliding mode observer signal and usmoEta, sgn (S), eta is real number, S is sliding mode surface and negative value, and is designed as
Figure BDA0003291563990000064
The error equation is now obtained as follows:
Figure BDA0003291563990000065
error of the measurement
Figure BDA0003291563990000066
Intermediate volume
Figure BDA0003291563990000067
In order to ensure the occurrence of the sliding mode, the stable condition of the variable structure of the sliding mode must be met; thus, consider the following Lyapunov function: v is 0.5s2V and s represent the speed versus time in a motion system, differentiating V from time t:
Figure BDA0003291563990000068
namely:
Figure BDA0003291563990000069
in order to ensure the stability of the extended sliding-mode observer, a stability condition must be satisfied
Figure BDA00032915639900000610
Eta is<-|e2-B0e1In practical application, the following parameter adaptation law is adopted: eta ═ l | e2-B0e1|,l>1, l is the safety factor of the sliding mode; error e1And derivatives thereof
Figure BDA00032915639900000611
Can converge to zero along a sliding mode occurring within a limited time, i.e.
Figure BDA00032915639900000612
Therefore, the error equation is simplified to
Figure BDA00032915639900000613
Namely, it is
Figure BDA00032915639900000614
Thus, e2The results of (a) are expressed as: e.g. of the type2=e-mt[C+∫r·emtdt]Where e is a mathematical constant and C is a constant, in order to ensure a disturbance estimation error e2Converging to zero, and selecting sliding mode parameters as parameters m of disturbance estimation>0; therefore, the sliding-mode observer parameters are extendedThe choice of number is constrained by η and m;
speed after occurrence of sliding mode
Figure BDA00032915639900000615
Observer equation simplification to
Figure BDA0003291563990000071
Thus, the equation is obtained:
Figure BDA0003291563990000072
the observer is equivalent to a low-pass filter after the sliding mode occurs; therefore, the disturbance observation effect of the extended sliding mode observer is equivalent to the output of the actual system disturbance of the filtering.
In the method for identifying mechanical parameters of the permanent magnet synchronous motor based on the extended sliding mode observer, in the step 3), the permanent magnet synchronous motor is enabled to run at two different stable speeds, so as to estimate the damping coefficient B;
when the motor runs at a first steady-state rotating speed omega (t), based on the extended sliding-mode observer and the disturbance d, the following disturbance estimation is obtained:
Figure BDA0003291563990000073
Figure BDA0003291563990000074
respectively, the estimated values of delta B and delta J;
after a time delay τ of the motor, the disturbance of the second steady speed ω (t + τ) is estimated as follows,
Figure BDA0003291563990000075
when the permanent magnet synchronous motor control system reaches the steady-state rotating speed, the load torque TLRegarded as a constant, the above two are subtracted:
Figure BDA0003291563990000076
Therefore, the temperature of the molten metal is controlled,
Figure BDA0003291563990000077
writing:
Figure BDA0003291563990000078
thus, an estimated parameter is obtained
Figure BDA0003291563990000079
Comprises the following steps:
Figure BDA00032915639900000710
in the above method for identifying mechanical parameters of a permanent magnet synchronous motor based on the extended sliding-mode observer, in step 3), when the motor is in the first constant acceleration state ac1When the method is operated, based on an extended sliding-mode observer and a disturbance equation, the B estimation value obtained in the previous step is substituted, and the following disturbance estimation value can be obtained:
Figure BDA0003291563990000081
then, by ensuring that the two different constant accelerations correspond exactly to two different time periods, respectively, after a time delay τ, the second constant acceleration state ac2The interference estimate of (c) is given by:
Figure BDA0003291563990000082
the estimation parameters can be obtained by subtracting and sorting the two formulas
Figure BDA0003291563990000083
Figure BDA0003291563990000084
In the above method for identifying mechanical parameters of a permanent magnet synchronous motor based on the extended sliding-mode observer, in step 3), when the parameters B and J are accurately estimated, the original estimation B in the state equation is performed0And J0Estimated parameters
Figure BDA0003291563990000085
And
Figure BDA0003291563990000086
instead, the extended sliding-mode observer is then rewritten as:
Figure BDA0003291563990000087
at this time, the updated estimated disturbance
Figure BDA0003291563990000088
Expressed as:
Figure BDA0003291563990000089
it can be seen that the updated extended sliding-mode observer has the capability of estimating the load torque, and therefore, in practical applications, the proposed extended sliding-mode observer is used to estimate B and J, or to estimate the load torque T on-line in real time with B and J being knownL
The invention has the beneficial effects that:
1. according to the problem that the control performance of the permanent magnet synchronous motor is influenced due to variable load conditions in an actual permanent magnet synchronous motor control system, the extended sliding mode observer can observe mechanical parameters of the permanent magnet synchronous motor control system in real time and feed back the mechanical parameters to the controller to adjust the parameters in time, so that system disturbance is restrained, and good dynamic performance and stable performance of the system are guaranteed.
2. The extended sliding-mode observer provided by the invention can be equivalent to a low-pass filter after sliding mode occurs, and the disturbance observation effect of the extended sliding-mode observer is equivalent to the output of filtering actual system disturbance. The cut-off frequency of the low-pass filter is m, and the low-pass filter can be designed arbitrarily according to the buffeting suppression requirement of the observer. Therefore, the output of the observer does not contain sliding mode buffeting caused by a low-pass filter and can be directly used for system control.
3. The invention can reduce the error value of the observer to zero in the effective time, has stronger anti-interference capability to external disturbance, simple structure and low sensitivity to system parameter change.
Drawings
Fig. 1 is a block diagram of estimation of mechanical parameters of a permanent magnet synchronous motor.
Fig. 2 is a schematic diagram of an extended sliding-mode observer.
Fig. 3 is a schematic diagram of an equivalent low-pass filter of an extended sliding-mode observer.
Fig. 4 is a schematic diagram of mechanical parameter estimation.
FIG. 5 is a graph of experimental results of expanding damping coefficients identified by a sliding-mode observer.
FIG. 6 is a graph of experimental results of expanding the rotational inertia identified by the sliding-mode observer.
Fig. 7 is a graph of experimental results of load torque identified by the extended sliding mode observer.
Detailed Description
The invention is further described below with reference to the figures and examples.
As shown in fig. 1, a method for identifying mechanical parameters of a permanent magnet synchronous motor based on an extended sliding-mode observer includes the following steps:
1) establishing a mathematical model of the surface-mounted permanent magnet synchronous motor, and performing coordinate transformation by adopting vector control: firstly, a mathematical model under an A-B-C three-phase static coordinate system is established to obtain an equation of each phase voltage, a flux linkage equation, a torque equation and a mechanical motion equation, and then the mathematical models under an alpha-beta two-phase static coordinate system and a d-q axis two-phase rotating coordinate system are obtained through Clark conversion and Park conversion in sequence, so that complete decoupling of exciting current and torque current is realized to facilitate control.
The specific process of the step 1) is as follows:
firstly, establishing a mathematical model under an A-B-C three-phase static coordinate system, wherein a voltage equation of the A-B-C three-phase static coordinate system is as follows:
Figure BDA0003291563990000101
wherein R issIs stator resistance, UA、UB、UCIs a stator three-phase voltage iA、iB、iCFor stator three-phase currents, ΨsA、ΨsB、ΨsCThe stator is a stator three-phase full magnetic linkage, and the vector form is as follows:
Figure BDA0003291563990000102
in the formula psisFor stator flux linkage, USIs the stator voltage; i.e. isIs the stator current;
the flux linkage equation is:
Figure BDA0003291563990000103
in the formula, ΨfA、ΨfB、ΨfCFlux linkage i for three-phase interlinkage of permanent magnet field and statorsA、isB、isCRespectively representing A, B, C three-phase stator currents, LsIn order to realize the synchronous inductance of the stator,
Figure BDA0003291563990000104
since the distribution of the air gaps is uniform, the rotor does not affect the self-inductance and the mutual inductance of the stator and A, B, C phases, so they are constant. L isA=LB=LC=L+LmlWherein L isA、LB、LCFor self-induction of each phase, LmlExciting inductance for each phase; l isIs the leakage inductance of each phase;
the torque equation:
the electromagnetic torque is considered as a result of the interaction between the stator, rotor and armature, and is expressed as:
Figure BDA0003291563990000105
wherein p is the number of pole pairs, TeBeing electromagnetic torque, ΨfMagnetic linkage psi linking the permanent magnet field with the statorsA stator flux linkage;
mechanical equation of motion:
Figure BDA0003291563990000111
ωmas mechanical angular speed of the rotor, TLThe load torque is J, the moment of inertia is J, and the damping coefficient of the motor is B;
and finally obtaining a mathematical model under a d-q axis two-phase rotating coordinate system through Clark transformation and Park transformation to realize complete decoupling of exciting current and torque current so as to perform vector control:
the flux linkage equation:
Figure BDA0003291563990000112
in the formula, Ψd、ΨqIs the quadrature-direct component of the stator flux linkage, Ld、LqComponent of stator inductance in d, q axes, id、iqD and q axis stator currents;
voltage equation:
Figure BDA0003291563990000113
in the formula of Ud、UqThe quadrature-direct component, omega, of the stator voltageeRepresenting the rotor angular velocity.
The torque equation:
Figure BDA0003291563990000114
in the formula, Ld、LqThe inductance values are direct axis and quadrature axis inductance values; therefore, the decoupling of the torque is realized under the d-q axis, and the torque equation is as follows:
Figure BDA0003291563990000115
2) designing an extended sliding mode observer and determining observer related parameters: according to a mechanical motion equation of the PMSM, a damping coefficient parameter error, a rotational inertia parameter error and a load torque system disturbance are used as an extended system state, an extended state equation is established and used for designing an extended sliding mode observer, and the stability of the observer is ensured by considering a Lyapunov function.
The motor's kinetic equation is expressed as follows:
Figure BDA0003291563990000121
wherein B is0、J0Is a rough estimation of real parameters, Δ B, Δ J are parameter errors between the real system and its rough estimation, and the moment of inertia J ═ J0+ Δ J, damping coefficient B ═ B0+ Δ B, ω is the rotational speed, d represents disturbances including parameter errors and load disturbances;
considering the disturbance d as an extended system state, in order to perform parameter estimation, the extended sliding mode observer is designed as follows:
Figure BDA0003291563990000122
wherein,
Figure BDA0003291563990000123
is an estimate of the disturbance,
Figure BDA0003291563990000124
is an estimate of the velocity, m is a sliding mode parameter, usmoRepresents the sliding mode observer signal and usmoEta, sgn (S), eta is real number, S is sliding mode surface and negative value, and is designed as
Figure BDA0003291563990000125
The error equation is now obtained as follows:
Figure BDA0003291563990000126
error of the measurement
Figure BDA0003291563990000127
In order to ensure the sliding mode, parameters of the observer must be reasonably selected, namely the stable condition of the sliding mode variable structure must be met; thus, consider the following Lyapunov function: v is 0.5s2V and s represent the speed versus time in a motion system, differentiating V from time t:
Figure BDA0003291563990000128
namely:
Figure BDA0003291563990000129
in order to ensure the stability of the extended sliding-mode observer, a stability condition must be satisfied
Figure BDA00032915639900001210
Eta is<-|e2-B0e1In practical application, the following parameter adaptation law is adopted: eta ═ l | e2-B0e1|,l>1, l is the safety factor of the sliding mode; in general, l ═ 2 is sufficient to ensure the stability of the observer.
It can be seen that the error e1And derivatives thereof
Figure BDA0003291563990000131
Can converge to zero along a sliding mode occurring within a limited time, i.e.
Figure BDA0003291563990000132
Therefore, the error equation is simplified to
Figure BDA0003291563990000133
Namely, it is
Figure BDA0003291563990000134
Thus, e2The results of (a) are expressed as: e.g. of the type2=e-mt[C+∫r·emtdt]Where e is a mathematical constant and C is a constant, in order to ensure a disturbance estimation error e2Converging to zero, and selecting sliding mode parameters as parameters m of disturbance estimation>0; therefore, the selection of the parameters of the extended sliding-mode observer is constrained by η and m;
as shown in FIG. 2, the speed after the sliding mode occurs
Figure BDA0003291563990000135
Observer equation simplification to
Figure BDA0003291563990000136
Thus, the equation is obtained:
Figure BDA0003291563990000137
the observer is equivalent to a low-pass filter after the sliding mode occurs, as shown in fig. 3; therefore, the disturbance observation effect of the extended sliding mode observer is equivalent to the output of the actual system disturbance of the filtering. Since the cut-off frequency of the low-pass filter is the same as the value of m, the design can be arbitrarily carried out according to the buffeting suppression requirement of the observer. Therefore, the output of the observer does not contain sliding mode buffeting caused by a low-pass filter and can be directly used for system control.
3) Observing system disturbance information by using the proposed extended sliding-mode observer, and extracting a damping coefficient, a rotational inertia and a load torque from the system disturbance information: firstly, the permanent magnet synchronous motor runs at two different steady-state speeds, the damping coefficient is estimated by utilizing the actually measured speed information and the estimated disturbance, then the rotary inertia is estimated by respectively running the permanent magnet synchronous motor at two different constant acceleration or constant deceleration, and finally, the damping coefficient and the rotary inertia are directly used for estimating the load torque.
As shown in fig. 4, when the motor is operating at a first steady-state rotational speed ω (t), based on the extended sliding-mode observer and the disturbance d, the following disturbance estimate is obtained:
Figure BDA0003291563990000141
Figure BDA0003291563990000142
respectively, the estimated values of delta B and delta J;
after a time delay τ of the motor, the disturbance of the second steady speed ω (t + τ) is estimated as follows,
Figure BDA0003291563990000143
when the permanent magnet synchronous motor control system reaches the steady-state rotating speed, the load torque TLRegarded as a constant, the above two equations are subtracted:
Figure BDA0003291563990000144
therefore, the temperature of the molten metal is controlled,
Figure BDA0003291563990000145
writing:
Figure BDA0003291563990000146
therefore, the temperature of the molten metal is controlled,obtaining an estimated parameter
Figure BDA0003291563990000147
Comprises the following steps:
Figure BDA0003291563990000148
when the motor is in a first constant acceleration state ac1When the method is operated, based on an extended sliding-mode observer and a disturbance equation, the B estimation value obtained in the previous step is substituted, and the following disturbance estimation value can be obtained:
Figure BDA0003291563990000149
then, by ensuring that the two different constant accelerations correspond exactly to two different time periods, respectively, after a time delay τ, the second constant acceleration state ac2The interference estimate of (c) is given by:
Figure BDA00032915639900001410
the estimation parameters can be obtained by subtracting and sorting the two formulas
Figure BDA00032915639900001411
Figure BDA00032915639900001412
Raw estimate B in the equation of state when parameters B and J are accurately estimated0And J0Estimated parameters
Figure BDA00032915639900001413
And
Figure BDA00032915639900001414
instead, the extended sliding-mode observer is then rewrittenComprises the following steps:
Figure BDA0003291563990000151
at this time, the updated estimated disturbance
Figure BDA0003291563990000152
Expressed as:
Figure BDA0003291563990000153
it can be seen that the updated extended sliding-mode observer has the capability of estimating the load torque, and therefore, in practical applications, the proposed extended sliding-mode observer is used to estimate B and J, or to estimate the load torque T on-line in real time with B and J being knownL
The method for identifying the mechanical parameters is simulated by using a Matlab simulation platform, a system template adopts a vector control structure, and i is selecteddThe identification results of the three parameters are respectively shown in fig. 5, 6 and 7, and the comparison of the actual data of the motor can find that the identification results are substantially consistent with the actual values. Therefore, the validity of the mechanical parameter identification algorithm is verified.
In summary, according to the actual control system of the permanent magnet synchronous motor, aiming at the problem that the control performance of the permanent magnet synchronous motor is affected due to variable load conditions, the extended mechanical parameter sliding mode observer provided by the invention can observe the mechanical parameters of the permanent magnet synchronous motor in real time and feed back the mechanical parameters to the controller to adjust the parameters in time, so that system disturbance is inhibited, and good dynamic performance and stability of the system are ensured; compared with the prior art, the mechanical parameter identification based on the extended sliding-mode observer provides an effective way for improving the self-regulation control of the permanent magnet synchronous motor under the condition of variable load, and can be widely applied to a system of complex systems such as an electric automobile, a flywheel energy storage system, a wind energy conversion system and the like.

Claims (6)

1. A permanent magnet synchronous motor mechanical parameter identification method based on an extended sliding-mode observer is characterized by comprising the following steps:
1) establishing a mathematical model of the surface-mounted permanent magnet synchronous motor, and performing coordinate transformation by adopting vector control: firstly, establishing a mathematical model under an A-B-C three-phase static coordinate system to obtain a phase voltage equation, a flux linkage equation, a torque equation and a mechanical motion equation of each phase, and then obtaining the mathematical model under an alpha-beta two-phase static coordinate system and a d-q axis two-phase rotating coordinate system through Clark conversion and Park conversion in sequence to realize complete decoupling of exciting current and torque current so as to be convenient for control;
2) designing an extended sliding mode observer and determining observer related parameters: according to a mechanical motion equation of the PMSM, taking a damping coefficient parameter error, a rotational inertia parameter error and load torque system disturbance as an extended system state, establishing an extended state equation for designing an extended sliding mode observer, and ensuring the stability of the observer by considering a Lyapunov function;
3) observing system disturbance information by using the proposed extended sliding-mode observer, and extracting a damping coefficient, a rotational inertia and a load torque from the system disturbance information: firstly, the permanent magnet synchronous motor runs at two different steady-state speeds, the damping coefficient is estimated by utilizing the actually measured speed information and the estimated disturbance, then the rotary inertia is estimated by respectively running the permanent magnet synchronous motor at two different constant acceleration or constant deceleration, and finally, the damping coefficient and the rotary inertia are directly used for estimating the load torque.
2. The method for identifying the mechanical parameters of the permanent magnet synchronous motor based on the extended sliding-mode observer according to claim 1, wherein the specific process in the step 1) is as follows:
firstly, establishing a mathematical model under an A-B-C three-phase static coordinate system, wherein a voltage equation of the A-B-C three-phase static coordinate system is as follows:
Figure FDA0003291563980000021
wherein R issIs stator resistance, UA、UB、UCIs a stator three-phase voltage iA、iB、iCFor stator three-phase currents, ΨsA、ΨsB、ΨsCThe stator is a stator three-phase full magnetic linkage, and the vector form is as follows:
Figure FDA0003291563980000022
in the formula psisFor stator flux linkage, USIs the stator voltage; i.e. isIs the stator current;
the flux linkage equation is:
Figure FDA0003291563980000023
in the formula, ΨfA、ΨfB、ΨfCFlux linkage i for three-phase interlinkage of permanent magnet field and statorsA、isB、isCRespectively representing A, B, C three-phase stator currents, LsIn order to realize the synchronous inductance of the stator,
Figure FDA0003291563980000024
LA=LB=LC=L+Lmlwherein L isA、LB、LCFor self-induction of each phase, LmlExciting inductance for each phase; l isIs the leakage inductance of each phase;
the torque equation:
the electromagnetic torque is considered as a result of the interaction between the stator, rotor and armature, and is expressed as:
Figure FDA0003291563980000025
wherein p is the number of pole pairs, TeBeing electromagnetic torque, ΨfMagnetic linkage psi linking the permanent magnet field with the statorsA stator flux linkage;
mechanical equation of motion:
Figure FDA0003291563980000026
ωmas mechanical angular speed of the rotor, TLThe load torque is J, the moment of inertia is J, and the damping coefficient of the motor is B;
and finally obtaining a mathematical model under a d-q axis two-phase rotating coordinate system through Clark transformation and Park transformation to realize complete decoupling of exciting current and torque current so as to perform vector control:
the flux linkage equation:
Figure FDA0003291563980000031
in the formula, Ψd、ΨqIs the quadrature-direct component of the stator flux linkage, Ld、LqComponent of stator inductance in d, q axes, id、iqD and q axis stator currents;
voltage equation:
Figure FDA0003291563980000032
in the formula of Ud、UqThe quadrature-direct component, omega, of the stator voltageeRepresenting the rotor angular velocity;
the torque equation:
Figure FDA0003291563980000033
in the formula, Ld、LqThe inductance values are direct axis and quadrature axis inductance values; therefore, the decoupling of the torque is realized under the d-q axis, and the torque equation is as follows:
Figure FDA0003291563980000034
3. the method for identifying the mechanical parameters of the permanent magnet synchronous motor based on the extended sliding-mode observer according to claim 2, wherein in the step 2), the kinetic equation of the motor is expressed as follows:
Figure FDA0003291563980000035
wherein B is0、J0Is a rough estimate of the true parameter, the moment of inertia J ═ J0+ Δ J, damping coefficient B ═ B0+ Δ B, Δ J are the parameter errors between the real system and its rough estimate, ω is the rotation speed, d represents the disturbance, including parameter errors and load disturbances;
considering the disturbance d as an extended system state, in order to perform parameter estimation, the extended sliding mode observer is designed as follows:
Figure FDA0003291563980000041
wherein,
Figure FDA0003291563980000042
is an estimate of the disturbance,
Figure FDA0003291563980000043
is an estimate of the velocity, m is a sliding mode parameter, usmoRepresents the sliding mode observer signal and usmoEta, sgn (S), eta is real number, S is sliding mode surface and negative value, and is designed as
Figure FDA0003291563980000044
The error equation is now obtained as follows:
Figure FDA0003291563980000045
error of the measurement
Figure FDA0003291563980000046
Intermediate volume
Figure FDA0003291563980000047
In order to ensure the occurrence of the sliding mode, the stable condition of the variable structure of the sliding mode must be met; thus, consider the following Lyapunov function: v is 0.5s2V and s represent the speed versus time in a motion system, differentiating V from time t:
Figure FDA0003291563980000048
namely:
Figure FDA0003291563980000049
in order to ensure the stability of the extended sliding-mode observer, a stability condition must be satisfied
Figure FDA00032915639800000410
Thus η<-|e2-B0e1In practical application, the following parameter adaptation law is adopted: eta ═ l | e2-B0e1|,l>1, l is the safety factor of the sliding mode; error e1And derivatives thereof
Figure FDA00032915639800000411
Can converge to zero along a sliding mode occurring within a limited time, i.e.
Figure FDA00032915639800000412
Therefore, the error equation is simplified to
Figure FDA00032915639800000413
Namely, it is
Figure FDA00032915639800000414
Thus, e2The results of (a) are expressed as: e.g. of the type2=e-mt[C+∫r·emtdt]Where e is a mathematical constant and C is a constant, in order to ensure a disturbance estimation error e2Converging to zero, and selecting sliding mode parameters as parameters m of disturbance estimation>0; therefore, the selection of the parameters of the extended sliding-mode observer is constrained by η and m;
speed after occurrence of sliding mode
Figure FDA0003291563980000051
Observer equation simplification to
Figure FDA0003291563980000052
Thus, the equation is obtained:
Figure FDA0003291563980000053
the observer is equivalent to a low-pass filter after the sliding mode occurs; therefore, the disturbance observation effect of the extended sliding mode observer is equivalent to the output of the actual system disturbance of the filtering.
4. The method for identifying the mechanical parameters of the permanent magnet synchronous motor based on the extended sliding-mode observer according to claim 3, wherein in the step 3), the permanent magnet synchronous motor is operated at two different stable speeds, so as to estimate the damping coefficient B;
when the motor runs at a first steady-state rotating speed omega (t), based on the extended sliding-mode observer and the disturbance d, the following disturbance estimation is obtained:
Figure FDA0003291563980000054
Figure FDA0003291563980000055
respectively, the estimated values of delta B and delta J;
after a time delay τ of the motor, the disturbance of the second steady speed ω (t + τ) is estimated as follows,
Figure FDA0003291563980000056
when the permanent magnet synchronous motor control system reaches the steady-state rotating speed, the load torque TLRegarded as a constant, the above two equations are subtracted:
Figure FDA0003291563980000057
therefore, the temperature of the molten metal is controlled,
Figure FDA0003291563980000058
writing:
Figure FDA0003291563980000059
thus, an estimated parameter is obtained
Figure FDA00032915639800000510
Comprises the following steps:
Figure FDA0003291563980000061
5. the method for identifying the mechanical parameters of the permanent magnet synchronous motor based on the extended sliding-mode observer according to claim 4, wherein in the step 3), when the motor is in the first constant acceleration state ac1When in operation, based on the extended sliding-mode observer and the disturbance equation, the previous steps are performedThe resulting B estimate is substituted, and the following disturbance estimate can be obtained:
Figure FDA0003291563980000062
then, by ensuring that the two different constant accelerations correspond exactly to two different time periods, respectively, after a time delay τ, the second constant acceleration state ac2The interference estimate of (c) is given by:
Figure FDA0003291563980000063
the estimation parameters can be obtained by subtracting and sorting the two formulas
Figure FDA0003291563980000064
Figure FDA0003291563980000065
6. The method for identifying the mechanical parameters of the permanent magnet synchronous motor based on the extended sliding-mode observer according to claim 5, wherein in the step 3), when the parameters B and J are accurately estimated, the original estimation B in the state equation is adopted0And J0Estimated parameters
Figure FDA0003291563980000066
And
Figure FDA0003291563980000067
instead, the extended sliding-mode observer is then rewritten as:
Figure FDA0003291563980000068
at this time, the updated estimated disturbance
Figure FDA0003291563980000069
Expressed as:
Figure FDA00032915639800000610
it can be seen that the updated extended sliding-mode observer has the capability of estimating the load torque, and therefore, in practical applications, the proposed extended sliding-mode observer is used to estimate B and J, or to estimate the load torque T on-line in real time with B and J being knownL
CN202111165630.0A 2021-09-30 2021-09-30 Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding mode observer Active CN114006557B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111165630.0A CN114006557B (en) 2021-09-30 2021-09-30 Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding mode observer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111165630.0A CN114006557B (en) 2021-09-30 2021-09-30 Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding mode observer

Publications (2)

Publication Number Publication Date
CN114006557A true CN114006557A (en) 2022-02-01
CN114006557B CN114006557B (en) 2023-11-10

Family

ID=79922232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111165630.0A Active CN114006557B (en) 2021-09-30 2021-09-30 Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding mode observer

Country Status (1)

Country Link
CN (1) CN114006557B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114531083A (en) * 2022-03-15 2022-05-24 北京理工大学 Permanent magnet synchronous motor single current sensor prediction control and parameter identification method
CN114744938A (en) * 2022-04-11 2022-07-12 华中科技大学 Kalman filtering-based full-parameter observer and full-parameter identification method
CN114915227A (en) * 2022-06-10 2022-08-16 华中科技大学 DC speed prediction control method of permanent magnet synchronous motor
CN115664283A (en) * 2022-11-18 2023-01-31 南京信息工程大学 Sliding mode control method and system based on generalized parameter estimation observer
CN118646301A (en) * 2024-08-15 2024-09-13 浙江大学 Permanent magnet synchronous motor mechanical parameter identification method based on extended disturbance observer

Citations (4)

* 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
CN105827168A (en) * 2016-05-10 2016-08-03 华中科技大学 PMSM control method and system based on sliding mode observation
WO2018077593A1 (en) * 2016-10-27 2018-05-03 Siemens Aktiengesellschaft Method for determining a deadzone angle of a backlash in a mechanical drive-train system, method for controlling a drive motor controller as well as drive-train system
CN110557070A (en) * 2019-09-30 2019-12-10 山东深川变频科技股份有限公司 permanent magnet synchronous motor parameter identification method based on second-order sliding-mode observer

Patent Citations (4)

* 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
CN105827168A (en) * 2016-05-10 2016-08-03 华中科技大学 PMSM control method and system based on sliding mode observation
WO2018077593A1 (en) * 2016-10-27 2018-05-03 Siemens Aktiengesellschaft Method for determining a deadzone angle of a backlash in a mechanical drive-train system, method for controlling a drive motor controller as well as drive-train system
CN110557070A (en) * 2019-09-30 2019-12-10 山东深川变频科技股份有限公司 permanent magnet synchronous motor parameter identification method based on second-order sliding-mode observer

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
侯利民;徐越;何佩宇;阎馨;: "辨识转动惯量扰动观测器的PMSM滑模控制", 电机与控制学报, no. 09 *
张明晖;杨家强;陈磊;楼佳羽;: "基于扩张状态观测器的永磁电机电流预测控制", 浙江大学学报(工学版), no. 07 *
张细政;王耀南;杨民生;: "永磁同步电机无位置传感器双滑模鲁棒控制", 电机与控制学报, no. 06 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114531083A (en) * 2022-03-15 2022-05-24 北京理工大学 Permanent magnet synchronous motor single current sensor prediction control and parameter identification method
CN114531083B (en) * 2022-03-15 2024-04-12 北京理工大学 Permanent magnet synchronous motor single current sensor prediction control and parameter identification method
CN114744938A (en) * 2022-04-11 2022-07-12 华中科技大学 Kalman filtering-based full-parameter observer and full-parameter identification method
CN114744938B (en) * 2022-04-11 2024-06-04 华中科技大学 Full-parameter observer based on Kalman filtering and full-parameter identification method
CN114915227A (en) * 2022-06-10 2022-08-16 华中科技大学 DC speed prediction control method of permanent magnet synchronous motor
CN114915227B (en) * 2022-06-10 2024-09-03 华中科技大学 DC speed prediction control method of permanent magnet synchronous motor
CN115664283A (en) * 2022-11-18 2023-01-31 南京信息工程大学 Sliding mode control method and system based on generalized parameter estimation observer
CN118646301A (en) * 2024-08-15 2024-09-13 浙江大学 Permanent magnet synchronous motor mechanical parameter identification method based on extended disturbance observer

Also Published As

Publication number Publication date
CN114006557B (en) 2023-11-10

Similar Documents

Publication Publication Date Title
CN114006557B (en) Permanent magnet synchronous motor mechanical parameter identification method based on extended sliding mode observer
CN110429881B (en) Active-disturbance-rejection control method of permanent magnet synchronous motor
CN101938246A (en) Fuzzy fusion identification method of rotating speed of sensorless motor
CN110707981A (en) Permanent magnet synchronous motor speed controller based on novel extended state observer
CN107911057A (en) A kind of rotor-position robust observation procedure for flywheel energy storage system
CN108390597A (en) Permanent magnet synchronous motor nonlinear predictive controller design with disturbance observer
CN109728755A (en) A kind of PMSM inverting TSM control method
CN110995102A (en) Direct torque control method and system for permanent magnet synchronous motor
CN109639200B (en) Rotational inertia online identification method based on motor load torque detection
CN111193448A (en) Surface-mounted permanent magnet synchronous motor load torque observation method based on extended Kalman filter
CN114726278A (en) Permanent magnet synchronous motor self-adaptive control method based on mechanical parameter identification
CN106208859B (en) Permanent magnet synchronous motor control method based on interference observer and repetitive controller
Wang et al. Rotor position estimation method for permanent magnet synchronous motor based on super-twisting sliding mode observer
Xiong et al. Sensor-less complex system control of pmsm based on improved smo
Liao et al. An integrated observer framework based mechanical parameters identification for adaptive control of permanent magnet synchronous motor
CN117895851A (en) Full-speed domain control method for surface-mounted permanent magnet synchronous motor
Shi et al. Current controller based on active disturbance rejection control with parameter identification for PMSM servo systems
CN108429501B (en) Method for observing load disturbance of permanent magnet synchronous motor
CN109687792A (en) Traction electric machine rotor parameter on-line identification optimization method towards vector control system
CN115566954B (en) Embedded motor speed regulation control compensation method and system
Shao et al. Parameter identification of permanent magnet synchronous motor based on multi-innovation least square method
CN105871277A (en) Minimum variance-based nonlinear model prediction controller design method for permanent magnet servo system
Che et al. A New SMO for Speed Estimation of Sensorless Induction Motor Drives at Zero and Low Speed
Wang et al. Design of Full-Order State Observer for Motor Drive Systems Based on the Fixed Gain Filter
Ju et al. PMSM speed control method based on Kalman filter and dynamic fuzzy control in electric vehicle

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