Finite Control Set Model Predictive Speed Control of a DC Motor

Thepaper describes the design procedure for a finite control setmodel predictive control (FCS-MPC) of brushed permanentmagnet DC (PMDC)machine supplied fromDC-DC converter. Full order linear Kalman filter is used for estimation of an unmeasured load torque and reduction of speed measurement noise. A new cost function has been introduced with a feedforward dynamic current component and a feedforward static load current component.The performance of the proposed control strategy is compared to the conventional PI-PWM cascade speed control through the experimental verification on the 250W laboratory prototype. Obtained results show excellent dynamic behaviour and indicate possible energy savings of the proposed speed control.


Introduction
Permanent magnet DC (PMDC) machines with brushes are being of less importance today and they are continuously replaced by permanent magnet synchronous machines (PMSM) or brushless DC (BLDC) machines.Despite this fact, PMDC machines still have their field of use in various types of applications, from the machine tools and small home appliances to photovoltaic, robotic, and medical applications.They are widely used as auxiliary drives in automotive solutions due to presence of a DC grid in modern vehicles, and they are especially suitable for a battery operated equipment.The biggest drawback of these machines is a commutator, which needs regular maintenance and causes sparkles, and, therefore, the machines cannot be used in explosive environment.On the other hand, linear speed-torque curve, high torque at low speed, and the simplicity of torque/speed control are well-known advantages of the conventional PMDC machines.Coreless PMDC machines have almost no cogging torque, less vibrations, and smoother motor running at low speed against conventional PMDC machines.
In order to work in all four quadrants, these machines are often supplied by H-bridge DC-DC converter, which changes mean value of the voltage on the machine terminals.This converter uses pulse-width modulation (PWM) with 2-level or 3-level type of control.
The PMDC machines have been extensively studied and many advanced control strategies have been presented in the past, but they are studied by many researchers even today [1][2][3][4][5][6][7].Among the new approaches, a new and attractive alternative for the control of H-bridge supplied PMDC machine is finite control set model predictive control (FCS-MPC) [8].A framework for the analytical derivation of FCS-MPC can be found in [9].FCS-MPC has been applied to the control of several types of electrical drives.Techniques for reducing a steady-state error of FCS-MPC controlled H-bridge have been presented in [10].Detailed description of the algorithm design for the commutation torque ripple minimization applied to a BLDC machine can be found in [11].Direct speed control of a PMSM for a single-mass system has been developed in [12].Predictive strategy for the speed control of a two-mass system driven by PMSM has been proposed in [13], introducing reduced order extended Kalman filter for mitigation of oscillations and compensation of load torque disturbances.In authors' opinion and supported by the literature review, a FCS-MPC for the PMDC machine has not been reported yet.
Novelty of this paper lies in the following: the predictive control strategy is, for the first time, proposed for the combination of DC-DC converter and PMDC machine.The proposed predictive control is compared to the conventional approach based on PI controllers.In addition, we have proposed a new cost function, which is different compared to the commonly used ones [9].In the novel cost function, two current components are introduced: a current component covering a dynamic feedforward torque,   , and a current component covering a static load torque,   .We will also show how these components affect overall control performance of FCS-MPC controller.

Drive Model
where   is the armature resistance,   is the armature inductance,   and V  are the armature current and voltage, respectively, and V emf is the back-emf voltage.The electric torque generated by the motor is given by where   is the magnitude of motor torque and   is the motor torque constant.From (2) it is obvious that nonlinearity of magnetization curve is neglected in the machine model.The mechanical part of the machine is given by where   and  are the motor inertia and viscous friction, respectively;   is the shaft angular speed; and   is the load torque acting as an unknown external disturbance.For the sake of simplicity, it is assumed that characteristic of the external torque is nearly constant and the sampling time is short enough.Then, the disturbance torque is independent to the motor state variables and can be included in the system model as an augmentative state variable [14].The system equations are summarized in the state space representation as where The truncation should be performed until the input explicitly appears in the expansion [15].For the current equation (1), the order used is   = 1, but, for the angular speed equation (3),   = 2.It results in the following sampled-data model: where

State Observer
In this section, we introduce a full order Kalman filter to estimate the slowly varying load torque and to reduce impact of a quantization noise on the measurement of angular speed.KF proposed here is based on the predictor-corrector form [16].It is assumed that state variables are excited by random disturbances (e.g., inaccurately modeled dynamics), and so we consider the sampled-data model ( 7) and (8) in the presence of noise as follows: where k() is the vector of zero-mean process noise and w() is the vector of zero-mean measurement noise.These noises are assumed uncorrelated as follows: Process noise auto-covariance matrix, Q, is assumed to be where  11 is the variance of V 1 () and so forth.Measurement noise auto-covariance matrix, R, is assumed to be where  11 is the variance of  1 () and so forth.Then, the state estimate of the machine, which is further used in the control system (e.g., posteriori estimate), is given by K is the Kalman gain given by Predicted state estimate (e.g., a priori estimate) is given by The model ( 7) and ( 8) is assumed to be linear and time invariant.Thus, Kalman filter gain, K, will converge towards a steady-state gain.We precalculated value of the gain matrix, K, offline and then we use it in the online predictive control strategy.It is obvious that measurement noise is introduced by the incremental encoder and the current sensor.The incremental encoder with  pulses per revolution (ppr) is evaluated within the sampling time,  ENC .With quadrature evaluation in the control system, it yields to the quantization error that Regarding a standard deviation,  22 , of and for  ENC = 500 s and  = 1024 ppr, the variance of the speed measurement noise is then obtained by

Power Converter
PMDC motor is driven by conventional H-bridge (Figure 1) fed by DC link with the voltage value  DC .The H-bridge has four admissible switching states, generating three output voltage values, + DC , − DC , and zero voltage (Table 1).
In fact, there are two possibilities of PWM control of such DC-DC converter: two-level or three-level voltage control.In the first case, both legs of converter are complementary controlled, and thus only positive or negative DC bus voltage is applied.In the second case, one leg holds defined switching state and the other leg is switching according to voltage demand.FCS-MPC does not use exactly one or another approach.Instead, the switching state is directly chosen by a FCS-MPC controller, as explained in Section 6.

Speed Control Structure with PI Controller and Feedforward Compensation
In this section, we provide only a short note of the conventional speed control based on PI controller, shown in Figure 2.
It is presented here for two reasons: we will compare it to the novel predictive controller and, furthermore, it creates a basic theory to introduction of novel current components.Later, these components will be implemented into a cost function.
It is expected that torque response of PMDC machine is several times faster than the closed loop mechanical dynamics.Thus, the electrical dynamics can be neglected (  =  ref ), and torque setpoint,  ref , can be described by the following: where   is the output of the PI controller,   is the dynamic feedforward torque, and T is the feedforward compensation of estimated load torque.With the use of the standard PI controller having transfer function, and simple mechanical dynamics (neglecting often a very small value of ), the closed loop transfer function of the system in Figure 2 without feedforward (i.e., at   () = 0) will be In servo applications, when the system is driven from one position to another, trapezoidal speed references are used instead of stepwise speed change.Speed reference during acceleration and deceleration must be tracked without any steady-state error [17].This can be achieved with two-degreeof-freedom (2DOF) control structure.In other words, the feedforward   () must be properly set up.It can be easily shown that if the feedforward term is used in Figure 2, closed loop transfer function (22) changes to an ideal unity gain.Because we neglected the dynamics of the electrical circuit and parameters of the machine may vary with time, the unity gain is only an ideal case.The control structure in Figure 2 with the feedforward is capable of assuring a tracking of trapezoidal speed reference with the high precision.
It is obvious from (23) that if the feedforward is to be applied, a value of reference angular acceleration,  ref , must be available.In servo applications, this value is generated from a trajectory generator or a ramp generator [18].Based on ( 2) and (23), we can write an equation for a dynamic feedforward component of current as follows: In a similar way, a static load torque component of current can be written as We formulated these current components for later use in the novel cost function and we will also show their influence on the control performance.

Novel FCS-MPC Based Control Structure
where î and ω are the current and the speed filtered by KF, T is the estimated load torque, considered as the second system input introduced in order to predict correctly the future states,   is the vector of predicted states for th admissible voltage, and is the admissible voltage value.The predicted values are compared using cost function and the voltage, for which the cost function is evaluated as minimal, is then applied to the converter.Figure 4 presents the control algorithm for FCS-MPC of H-bridge driven PMDC machine.
It consists of steps that can be summarized as follows: (1) Measurement of the machine current and speed.
(2) Filtration of the current and speed by KF and estimation of a load torque.(3) Calculation of required current components   and   .(4) Prediction of the future values of speed and current for each admissible voltage.(5) Evaluation of the cost function, , for three voltages, using variable state as an argument for which the cost function  is evaluated as minimal (compare to [17]).( 6) Application of the optimal actuation to the H-bridge.

Controller Design. The following criteria have been proposed for the predictive controller:
(i) Precise reference tracking for the step and trapezoidal speed reference.(ii) Limitation of the maximal current.(iii) Smooth behaviour of the electrical torque.(iv) Fast disturbance rejection.
The cost function implementation itself can be considered as a controller, which can fulfill abovementioned demands at once.The following cost function has been introduced: where   and   are predicted values of the machine angular speed and current, respectively.The cost function consists of the sum of three terms: (a) This term favors the voltage vectors that achieve reference speed tracking and it is weighted with weighting factor,  1 .
(b) This term provides the current control of the machine and is weighted with weighting factor,  2 .The summation of the aforementioned currents   and   is the current reference.
(c) This term in the cost function presents current limitation, where  lim is nonlinear function defined as in [17]: The current limitation works on the simple principle that when predicted value of the current for any voltage is above the limit,  MAX , the cost function gets very high value (10 12 ), forcing the controller to choose another voltage.

Experimental Results
Proposed For 0:  p2 , i ap2 implemented into simulator.The second PMDC machine was used as a load machine with the same parameters.Experimental results were taken for both control approaches under the same conditions, such as the current value limitation  MAX = 10 A and the sampling time   = 50 s.
A standard cascade speed control structure with the inner current control loop and the superimposed speed control loop with antiwindup and dynamic feedforward is denoted in experimental results as PI.It uses 2-level voltage control with 10 kHz PWM frequency.FCS-MPC speed control, as presented in Section 6, is denoted in experimental results as MPC.For a good comparability of both methods, the controllers were tuned to reach the same settling time and overshoot.It resulted in the values   = 0.5 Nms and   = 11 Nm for PI controller and  1 = 150 and  2 = 1 for MPC controller.
The experimental results comparing these control approaches are shown in Figures 5-8. Figure 5 shows the results of a reference step from 0 to 80 rad/s.Looking at Figure 5, it can be concluded that MPC controller leads to a faster disturbance rejection than PI controller.This is also supported by Figure 6, where reference tracking of a ramp signal from 0 to 80 rad/s with angular acceleration  ref = 1066.6rad/s 2 is shown.Both figures also show that MPC controller has higher oscillations of machine current than PI controller.These oscillations can be suppressed by increasing of  2 in the cost function (28), but this would lead to the deterioration of the speed performance.This is well-known disadvantage of FCS-MPC approaches [12].
The performance of the Kalman filter for load torque estimation is shown in Figure 7. Measured value (blue line) is obtained by the torque transducer.Estimated value (red line) during the acceleration (t = 0.025-0.1 s) does not track the measured value, because the torque sensor also measures a  dynamic torque, whereas the torque estimator was designed only for estimation of a static torque and load.Figure 8 shows the influence of the introduced current components on the speed control quality.If the   current component in the cost function ( 28) is zero, we can observe considerable dynamic tracking error in the acceleration and deceleration, but the steady-state value of the reference speed is reached.On the other hand, if the   current component is zero, dynamic tracking of the trapezoidal reference during acceleration is adhered, but, after the loading by external torque, the steady-state accuracy is not preserved.Therefore, a significance of proposed current components   and   in the cost function has been proved.We consider this as the main contribution of this paper since these current components can be easily applied to the other types of FCS-MPC controllers, even for AC drives.

Efficiency Issues and Parameter Sensitivity
8.1.Efficiency Issues.Now we compare number of switching states for both approaches.One transistor turn-on and turnoff yields the   = 2 number of switching states.A standard cascade PI controller in Figure 2 uses PWM with the fixed frequency, where the number of switching states (for Hbridge) can be defined as follows: where   is PWM carrier frequency and  RUN is a time period during which the number of switching states is calculated.For example, for H-bridge running with 10 kHz frequency PWM, the overall number of switching states during 1 second is   = 80 000.On the other hand, MPC controller has a variable switching frequency; thus, exact explicit formula for the number of switching states cannot be performed.It has been found that number of switching states of MPC differs when the drive runs at steady state and transients.Therefore, we proposed a dynamic factor,   , as the following: where   is a period, where the drive runs at transients (accelerating or decelerating), and   is a time, where drive runs at steady state; that is, if   = 100%, drive is all the time accelerating or decelerating.The number of switching states for the different values of angular acceleration,  REF , obtained by Matlab simulation is shown in Figure 9.It is obvious that, for MPC controlled drive, the number of switching states decreases with increasing dynamic factor.For example, for   = 100%, the drive needs 28 000 switching states, but, for   = 5%, the drive needs 44 000 switching states (the red line).Note that number of switching states was recalculated to the period of 1 second in order to obtain the fairest possible comparison with PI-PWM controller.
The number of switching states is directly connected to the switching losses of converter transistors and these losses could present a significant share of the total converter losses and, therefore, the total efficiency of the converter.
For DC-DC converters, the total losses are typically obtained by simulation or by calculation [19] and, in practice, efficiencies of 70% to 95% are achieved [20].In this contribution, an experiment has been provided instead of the calculation.The results are shown in Figure 10.Total efficiency has been measured during the interval of 1 s with various values of dynamic factor.The input and output voltage and current of the converter have been measured and then the average value of input and output power has been calculated.Note that, for the valid measurement of the output PWM voltage, there is a need for very short sampling time, which is beyond the computational limits of the used simulator; for example, for 1% precision measurement of 10 kHz output PWM voltage, a voltage sample for calculation has to be stored in memory every 1 s.For that reason, PWM voltage frequency in Figure 10 was lowered to 1 kHz.
The theoretical maximum of switching frequency using FCS-MPC is 100 kHz (using 10 s as sampling time) in steady state.However, during transients, it can be lowered to almost single switching state for a certain time period.Therefore, it can be assumed that, with the higher dynamic factor, the switching losses using FCS-MPC will be lower.
The resulting, experimentally measured efficiency using 1 kHz PI-PWM and FCS-MPC control is depicted in Figure 10.It is obvious that, with lower values of dynamic factor, the efficiency of PI-PWM is better than the efficiency of FCS-MPC.On the other hand, with dynamic factor rising to 100%, the performance of both control approaches is comparable.However, it is expected that, with the higher switching frequency of PWM, the switching losses will rise, and so the efficiency of FCS-MPC will overcome the efficiency of PI-PWM, even for the lower dynamic factor values.

Parameters' Variation.
In following experiments, the parameters of the FCS-MPC controller remained the same as in Section 7 and real values of inertia and armature resistance have been varied.The variation of the armature resistance,   , for 3 different values of real resistance is shown in Figure 11.It can be observed that, for 180% change in resistance, speed tracking is not preserved, because the maximum speed error is almost 5% (green line).In the same case, the current oscillations are mitigated, probably due to  the parasitic inductance of the used resistor.Note that the ripple observed in average speed error response is mainly caused by using a low cost encoder with only 1000 ppr, with quantization error, as was shown in Section 3. Variation of total inertia value is shown in Figure 12.Note that speed reference has been modified in order to keep the current into unsaturated values (i.e., below 10 A).The load torque has been applied at 0.6 s.It can be observed that, for the same controller parameters, the speed tracking during acceleration from t = 0.1 s to 0.4 is preserved for both original and raised inertia value.Even more, the disturbance rejection is excellent and minimal influence of the applied load torque can be seen in the speed error response.

Conclusion
Finite control set model predictive control has already been used in the control of various electrical drives, but its application for the speed control is still an open problem.An experimental validation of the finite control set model predictive speed control of PMDC machine has been presented.
A new cost function, with the consideration of additional feedforward current components, has been developed.It was shown that the new cost function significantly reduces steady-state error and improves dynamic speed tracking.In addition, proposed MPC controller fluently combines 2level and 3-level voltage control.This controller has been compared with the PI-PWM controller.For comparability, both controllers were tuned to the same settling time and   overshoot and the measurements were taken under the same boundary conditions.Presented type of PMDC machine speed control is especially suitable for the drives with very high dynamic demands.It can be concluded that FCS-MPC controller has a fast disturbance rejection and it needs less number of switching states than PI-PWM controller.But its complexity, computational burden, and torque oscillations may be a problem.Finally, it can be confirmed that FCS-MPC controllers are worthy of attention against conventional controllers.

6. 1 .
Control Algorithm.Block diagram of the FCS-MPC speed control is shown in Figure3.The proposed algorithm uses predictions of the future machine state in ( + 1) step for each value of admissible voltage.It means that three predictions are needed in each cycle.The prediction itself is based on the modified sampled-data model described in Section 2:

Figure 5 :Figure 6 :Figure 7 :
Figure 5: Experimental results for the step speed reference.

Figure 8 :
Figure 8: Experimental results for the trapezoidal speed referenceinfluence of different current components of the cost function.

Figure 9 :
Figure 9: Simulation results-number of switching states as a function of the dynamic factor.

Figure 10 :
Figure 10: Experimental results-measured values of converter efficiency as a function of the dynamic factor.
2.1.Continuous Model.In this section, a continuous time state space model and discrete model for the PMDC machine are obtained.The stator dynamics of the PMDC machine is described by