Adaptive Backstepping Based MTPA Sensorless Control of PM-Assisted SynRM with Fully Uncertain Parameters

1Hebei Key Laboratory of Distributed Energy Storage and Micro-Grid, North China Electric Power University, Baoding, Hebei 071003, China 2State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Baoding 071003, China 3Department of Technology and Information, State Grid Hebei Electric Power Company, Shijiazhuang, Hebei 050022, China


Introduction
Drive system is vital for various industrial applications.One of the important trends for the drive field is that direct current (DC) motor drive systems are being replaced by alternate current (AC) motor drive system due to low cost, convenient control, and superior performance of AC motor drive systems.A permanent magnet-assisted synchronous reluctance motor (PMa-SynRM) based drive system is a unique state-ofthe-art technology, where permanent magnets are equipped in magnetic barrier of rotor.Permanent magnet synchronous motor (PMSM) incorporated with synchronous reluctance motor makes PMa-SynRM appeal to the public in high efficiency, good utilization of the inverter, and wide range of flux-weak regulation.The uniqueness in structure of PMa-SynRM attracts the investigator's attention of optimizing the design of motor body and putting forward more diverse and fascinating structures [1,2].
As for its control, it is a bottleneck problem researched insufficiently.The synchronous change of speed and frequency for PMa-SynRM makes speed be regulated by adjusting voltage or frequency, which leads to three commonly used approaches to control: voltage-to-frequency control (VFC), direct torque control (DTC), and vector control (VC).VFC is a scalar and open-loop control, and its voltage regulated by space vector modulation always follows the change of given value.Low cost and ease of use are the prominent superiorities of VFC [3].However, the weaknesses of slow response, poor performance of control, and insufficient utilization of torque hinder its successful development.Takahashi and Noguchi [4] and Depenbrock [5] studied DTC about 30 years ago, which was applied to control asynchronous motor first and was gradually introduced in PMSM until it had been developed and matured [6].In terms of space vector and stator flux orientation, a hysteresis comparator is designed to control flux and torque.Through comparison of feedback value and given value in flux and torque, respectively, the states of converter switches are produced and torque output with high performance and rapid response is gained ultimately. Definite physical meaning of control signals and quick torque response are the remarkable assets for DTC.But torque ripple, current impulse, and massive online computations restrict its application [7].VC is originated from the principle of flux orientation and enables controlling flux and torque to be decoupled like DC motor [8].Therefore, it is famous for flexibility, high accuracy, and stable performance in low speed.For PMa-SynRM, maximum torque per ampere (MTPA) vector control has great advantages of large torque in low speed, small capacity in inverter unit, and facility to utilize reluctance torque compared to   = 0 vector control [9,10].The essence of MTPA control is to allocate --axis currents appropriately.The relationship of constrained -axis currents in MTPA principle is multiple, high-order, and nonlinear.The --axis reference currents predetermined through offline method are a good solution.One of the representative offline methods is look-up table method [11], which is easy to execute but adds burden to store and requires extra offline effort [12].Additionally, conventional MTPA control relies heavily on precise modeling parameters.Variation of parameters directly worsens its control performance.Especially for PMa-SynRM, previous investigations clearly demonstrate that its resistance, inductance, and flux linkage are directly related to operating temperature [13].Hence, research on MTPA control with uncertain parameters is a hot and thorny issue.Niazi and Toliyat [14] estimated --axis inductances and flux linkage of PMa-SynRM under multiple reference frame and [15] proposed a robust MTPA control scheme of PMa-SynRM for variation of --axis inductances.The actuality is that only partial system uncertainties are considered, which desires the other system parameters to be known.
Backstepping control (BC) is a new style recursive technology for uncertain nonlinear plant [16,17].With the aid of virtual control variables and plenty of recursive steps, the ultimate controller can be completed systematically and original high-order plant is reduced to a lower dimension.Furthermore, integration of backstepping control with adaptive control and uncertain plant parameters can be estimated through selecting a suitable Lyapunov function.A robust nonlinear controller based on BC and MTPA has been proposed for speed control of interior PMSM [18,19].Nevertheless, to the best of our knowledge, previous researches on employment of BC into PMa-SynRM have not been found.A majority of control schemes utilize conventional PID control [20].Moreover, three intrinsic deficiencies in traditional BC should be mentioned as follows: (A1) Linear treatment of uncertain plant parameters results in solving highly complicated regression matrices in the design of control; utilization of various linearization theories may ignore some beneficial nonlinearities of the plant.
(A2) Partial parameters uncertainties considered cause to the designed control scenario vulnerable to variation of extra values.
(A3) Determination of control parameters in the last time derivative of ultimate Lyapunov function requires the positive definite terms which only include partial model information.Once the derivative is induced to be positive definite by these positive terms, the asymptotic stability of the controller will be damaged and unfortunately it is nearly impossible to discover.
The innovation of the study is firstly to propose a new nonlinear and fully adaptive BC approach with no problems of overparameterization and singularity for an uncertain PMa-SynRM and hence (A1) is excluded.Previous researches such as that in [21] cannot resolve full parameter uncertainties, and the problems of overparameterization and singularity cannot be complete to be coped with.Furthermore, the paper enables overcoming the shortcomings of conventional BC, and totally seven parameter uncertainties of six motor parameters and one load parameter are taken into account, which construct a full adaptive structure and eliminate (A2) and (A3).Merely the number of pole pairs is considered to be known, due to the fact that it does not change in operation and can be acquired from nameplate.The other novelty in the context is that adaptive BC is introduced into MTPA control through implicit and symbolic computation methodology.Additionally, an active flux based improved reduced-order Luenberger speed observer is developed to estimate the speed of PMa-SynRM, resisting its strong and high nonlinearity of dynamic model.The remainder of the paper is arranged as follows: the dynamics of PMa-SynRM with fully uncertain parameters are introduced in Section 2. The relation constrained in --axis currents under MTPA control is explicitly derived in Section 3. Section 4 discloses the design scenario of the proposed controller with relation to Lyapunov stability analysis.An improved Luenberger speed observer on the basis of active flux is developed in Section 5 to estimate the speed.The effectiveness and correctness of the proposed algorithm are validated in Section 6 and the simulation results are discussed.The final section sums up the paper.

Dynamics of PMa-SynRM with Fully Uncertain Parameters
The model of a typical PMa-SynRM is expressed in rotating --axis coordinate system as follows [22]: Flux Linkage Equations Electromagnetic Torque Thus, the dynamic model of a PMa-SynRM can be described as where Due to the existence of product terms between electrical stator currents and mechanical rotor speed in ( 5) and ( 6), PMa-SynRM as an electromechanical coupling system belongs to a highly nonlinear plant.With regard to the various electrical parameters   ,   ,   , and   , it should be noted that they could be directly measured or calculated, but their values may vary with operating and experimental conditions such as temperature, humidity, and flux saturation.The remaining mechanical parameters   and   are even more difficult, which are improbable to measure or calculate.The final load torque   is frequently uncertain.It follows that PMa-SynRM suffers from inevitable parameter uncertainties and immeasurable disturbances.For this purpose, nonlinear adaptive control should be introduced to eliminate uncertainties and reject disturbances.

Parameter Constraint Relation in MTPA Control
The key of MTPA is to pursue electromagnetic torque as large as possible through reasonable distribution of --axis stator currents.This problem can be transformed into an extreme problem in the following: where   is the stator current.Through introduction of a Lagrangian multiplier , the extreme problem given in ( 8) is converted into an auxiliary function solution problem in the following: Using partial derivatives of   ,   , and  in ( 9), respectively, we can obtain Let ( 10) be equal to zero; then the relation constrained in --axis currents in MTPA control can be derived as Equation (11) indicates that if   is acquired,  ref will be determined.

Nonlinear Adaptive Backstepping Controller
The control target of the study is to guarantee all the signals to be bounded and ensures the speed and currents to track their respective references precisely and rapidly in spite of full parameter uncertainties in PMa-SynRM and load disturbance.For achievement of the goal, the designed controller enables tracking the variations of parameters uncertainties.Hence, online parameter estimation laws for fully uncertain parameters with no problems of overparameterization and singularity should be conducted and vary with actual parameters adaptively.

Designed Nonlinear and Fully Adaptive Backstepping
Controller.The tracking errors   ,   , and   of   ,   , and   are defined as follows: where  ref ,  ref , and  ref are the references of   ,   , and   , respectively.
Then, (21) can be rewritten as By derivative of   in ( 14), the dynamics of -axis current tracking error can be obtained as where In terms of   =  5 ,   =  6 , and   =  8 , ( 24) can be arranged as For reaching the control goal and stabilizing the tracking error terms, the control inputs can be designed as follows: where   and   are two positive feedback gains.Substituting ( 27) and ( 28) into ( 23) and ( 26), respectively, we can obtain Step 3. The positive definite Lyapunov function candidate is defined as follows: Derivative of the Lyapunov function yields By substituting (32) into (31), the derivative of the Lyapunov function can be simplified as follows after some mathematical computations: where   ,   , and   are the positive control gains.Suppose that where The final Lyapunov function is selected as where  is an eighth-order positive definite diagonal matrix; the diagonal elements are described as  1 ,  2 ,  3 ,  4 ,  5 ,  6 ,  7 , and  8 .
Then, the derivative of the Lyapunov function yields Letting we can obtain In terms of (36) and (40), the adaptive law depicting the uncertainties for the fully uncertain parameters can be written as By substituting (40) into (38), the derivative of the final Lyapunov function can be rewritten as Equation ( 42) represents that there exists a positive definite Lyapunov function  and its derivative V with respect to time being nonpositive for the closed-loop system.
Additionally, in terms of (41), the fully uncertain parameters of PMa-SynRM and load torque disturbance can be estimated and represented as (43)

Stability Analysis
Theorem 1.For the dynamics of PMa-SynRM described in ( 5)- (7) with parameter uncertainties and external load disturbance, the proposed controller shown in ( 27) with the adaptive law of fully uncertain parameters given in (41) enables ensuring the tracking error signals in ( 12)-( 14) to converge to (0, 0, 0) asymptotically through Barbalat's Lemma.
To prove this theorem, by integration of (42) in its both sides from 0 to +∞, we can give Thus, It follows that Furthermore, () is uniformly continuous and V() is finite.Through Barbalat's Lemma [23], it gives that which proves that the tracking errors (  ,   ,   ) can converge to (0, 0, 0) asymptotically under full parameters' uncertainties and uncertain load torque disturbance for any initial condition.It is worth mentioning that the choices of control parameter and adaption gain are no restriction for use, except for positive definite terms.Previous studies have shown that control gains are frequently necessary to satisfy inequality relation with uncertain parameters [21] and selection of adaptation gains in multiparameter estimation situation is complicated [24] as (A3) encountered.The relaxation of control and adaptation gains is beneficial for implementation of the proposed control scheme apparently.

Proposed Luenberger Speed Observer
Luenberger speed observer is an excellent linear estimation method.The product term of     in (4) represents the model nonlinearity of PMa-SynRM, which results in the fact that Luenberger speed observer cannot be used directly.Hence, an active flux   is firstly introduced to linearize (4) as follows: Then, the torque equation ( 4) can be rewritten as By substituting (48) into (3), we can obtain In terms of ( 2) and ( 50), the components of   in rotating --axis coordinate system are Equation ( 51) indicates that active flux is oriented on -axis, which is shown in Figure 1.Furthermore, for the mathematical model (see ( 5) and ( 6)), new input variables are introduced as By substituting (52) and ( 51) into ( 5)-( 7), we can get Equations ( 53)-( 55) manifest that the nonlinear model of the original system has been linearized.Additionally,   is only related to   and the linearized reduced-order observer can be constructed to estimate the speed.Through utilization of the simplified rotor motion equation ( 55) and -axis current equation ( 53), the reduced-order matrix equation can be given as ) .
(56) For the above matrix equation, linear Luenberger observer is designed as where î and ω are the estimations of -axis current and speed, respectively, and  is the feedback gain of designed Luenberger observer.By subtracting (57) from (56), we can get where u q =   6 i qＬ？＠ +   5 i q + n p   8   r i d + k q e q + 3n p  Equation (59) shows that the observed error of speed enables converging to zero through selecting a suitable parameter  to configurate the eigenvalue of matrix.Moreover, due to the fact that the designed speed observer is of first order, the computational burden is decreased distinctly.

Digital Simulation Analysis and Discussion
The overall configuration diagram of the proposed control scenario and concrete nonlinear backstepping controller with full adaptation law are shown in Figures 2 and 3, respectively.
To testify the superiority of the controller, adaption laws, and Luenberger speed observer adequately, the simulations consist of three groups.
Case 1. Due to the fact that very low speed regulation of PMa-SynRM is a complicated issue in control field, to validate the performance of the proposed controller and Luenberger    speed observer, the speed reference is set to a very low and changing value, which is 1 rad/s between 0 and 1 s, 2 rad/s between 1 s and 4 s, and 1 rad/s between 4 s and 6 s.That is to say that the simulation time in the first simulation is 6 s.Simultaneously, the sudden variation of load torque is taken into account and is abruptly applied to the motor as 4 N⋅m from 0 to 2 s, 2 N⋅m from 2 to 5 s, and 4 N⋅m from 5 s to 6 s.
The simulation results responding to the speed variation and torque disturbance are indicated in Figures 4(a Case 2. To illustrate the capability of the adaption law and testify the robustness of the control scheme retaining the full parameter perturbations furthermore, the parameter perturbations are described as follows: per phase stator resistance increases from 2.875 Ω to 4 Ω at 2 s and changes to 2.875 Ω at 5 s, -axis and -axis inductances decrease from 0.133 H and 0.058 H to 0.1 H and 0.05 H, respectively, at 2 s and return to their respective initial value at 5 s, viscous damping coefficient varies from 0.005 N/rad/s to 0.01 N/rad/s at 2 s and returns to 0.005 N/rad/s at 5 s, permanent magnet flux changes from 0.38 Wb to 0.5 Wb at 2 s and changes to 0.38 Wb at 5 s, and moment of inertia increases from 0.03 kg m 2 to 0.05 kg m 2 at 2 s and changes to 0.03 kg m 2 at 5 s.Additionally, in terms of the persistency of excitation condition, a good result with parameter estimation requires abundant signals, where the speed reference is chosen as  ref = 2 sin(), and simulation time is also 6 s.The load torque is the same as in Case 1. Figures 5(a)-5(g) plot the estimated errors of fully uncertain parameters, which clarify that all the uncertain and varied parameters are estimated by the adaption law rightly.As a consequence, the simulation results shown in Figure 5 exhibit the accuracy of the adaption law and reveal the good robustness of the presented control scenario.
Simultaneously, for demonstrating the influence of each parameter variation on the presented control scheme, the dynamic responses of speed tracking are displayed in Figure 6 when changing each parameter of the following parameters:   ,   ,   ,   ,   ,   , and   .The range of variation for each parameter is set to ±50% rated value.Figure 6 indicates that the speed of the motor tracks the reference precisely with the aid of the controller and the tracking errors approach zero.Hence, the proposal controller goes against all the parameter perturbations.Merely, when the parameter deviates from its rated value, the waveform responses of speed at startup will become larger, especially for   ,   , and   .In other words, the variations of these parameters   ,   , and   have a more significant influence on the speed response.However, all the speed responses in Figure 6 converge to their expected values within 0.3 s.Case 3.For various applications in usual, the motor operates at medium speed or high speed is indispensable.Therefore, this case indicates the simulation results in high speed condition.The speed reference is selected as 300 rad/s from zero to 1 s and 600 rad/s from 1 s to 4 s.After 4 s, the speed reference returns to 300 rad/s.Meanwhile, the load torque changes from 4 N⋅m to 2 N⋅m at 2 s and returns to 4 N⋅m at 5 s.
The results in this case are demonstrated in Figure 7.It is implied that, under the high speed condition, the estimated  Consequently, the presented control algorithm in the context achieves the closed-loop control of rotor speed and stator current and performs the estimations of fully unknown parameters for PMa-SynRM.

Conclusions
In this paper, on the basis of nonlinear dynamic model and parameter relation constrained in MTPA control, a nonlinear and robust adaptive backstepping speed sensorless control scheme with mismatched fully uncertain parameters is proposed for a PMa-SynRM.All the parameters of the motor and load torque are thought as uncertain except for the number of pole pairs of rotor.The design of the control scheme avoids solving complex regression matrices and the only prerequisite for all controller and adaption gains is positive definite, which reduces the workload of choice of feedback gains tremendously and ensures the asymptotic stability of the controller readily.Additionally, an active flux based Luenberger speed observer is developed.The speed observer resolves the nonlinearity of dynamic model for PMa-SynRM effectively.The proposed controller fulfills the global asymptotic tracking of a varying speed reference and --axis current references constrained by MTPA control.Simultaneously, the designed controller is of great robustness against uncertainties in both PMa-SynRM and load torque perturbation.Digital simulations have been implemented to verify the correctness and feasibility of the proposed control scheme well.

Figure 1 :
Figure 1: Distribution of active flux in 0 coordinated system.

Figure 2 :
Figure 2: The overall configuration diagram of the proposed control scenario.

Figure 3 :
Figure 3: The concrete nonlinear backstepping controller with full adaptation law.

Figure 4 :
Figure 4: The simulation results in Case 1.(a) The estimated speed and its reference; (b) -axis current and its reference; (c) -axis current and its reference; (d) --axis control voltages; (e) torque and its estimation.

Figure 5 :
Figure 5: The simulation results in Case 2. (a)   and its estimation; (b)   and its estimation; (c)   and its estimation; (d)   and its estimation; (e)   and its estimation; (f)   and its estimation; (g)   and its estimation.
)-4(e).The estimated speed and its reference are shown in Figure 4(a).The dynamic tracking responses of --axis currents and their references are in Figures 4(b) and 4(c) individually.The control inputs of --axis voltages are displayed in Figure 4(d).The estimated load torque and its reference are demonstrated in Figure 4(e).Figures 4(a)-4(c) indicate that the proposed controller can guarantee the output signals to track their respective references correctly and rapidly.Furthermore, Figure 4(a) verifies the effectiveness of the designed Luenberger speed observer.From Figure 4(e), the load torque is estimated precisely.

Figure 6 :
Figure 6: Speed tracking responses with the parameters variation.(a) Variation of   .(b) Variation of   .(c) Variation of   .(d) Variation of   .(e) Variation of   .(f) Variation of   .(g) Variation of   .

Figure 7 :
Figure 7: The simulation results in Case 3. (a) The estimated speed and its reference; (b) -axis current and its reference; (c) -axis current and its reference; (d) torque and its estimation.
and   are the --axis currents,   and   are the --axis voltages,   and   are the --axis flux linkages,   and   are the --axis inductors,   is the stator resistance,   is the number of pole pairs,   is the mechanical speed of rotor,   is the flux linkage of permanent magnet,   is the moment of inertia of rotor,   is the viscous damping coefficient, and   is the load torque denoting external torque disturbance.