Next Article in Journal
Study of Reversible Platelet Aggregation Model by Nonlinear Dynamics
Next Article in Special Issue
Adaptive Fault Estimation for Hyperbolic PDEs
Previous Article in Journal
High-Speed Interval Type-2 Fuzzy Systems for Dynamic Parameter Adaptation in Harmony Search for Optimal Design of Fuzzy Controllers
Previous Article in Special Issue
A Valid Dynamical Control on the Reverse Osmosis System Using the CESTAC Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Control Method for IPMSM Based on Active Disturbance Rejection Control and Model Predictive Control

1
School of Automation, Central South University, Changsha 410083, China
2
Hunan Xiangjiang Artificial Intelligence Academy, Changsha 410083, China
3
Power Construction Corporation of China, Jiangxi Electric Power Design Institute CO., LTD, Nanchang 330096, China
4
Department of Electrical and Electronic Engineering, Chiba University, Chiba 2638522, Japan
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(7), 760; https://doi.org/10.3390/math9070760
Submission received: 8 February 2021 / Revised: 24 March 2021 / Accepted: 29 March 2021 / Published: 1 April 2021
(This article belongs to the Special Issue Model Predictive Control and Optimization for Cyber-Physical Systems)

Abstract

:
In this paper, the speed tracking problem of the interior permanent magnet synchronous motor (IPMSM) of an electric vehicle is studied. A cascade speed control strategy based on active disturbance rejection control (ADRC) and a current control strategy based on improved duty cycle finite control set model predictive control (FCSMPC) are proposed, both of which can reduce torque ripple and current ripple as well as the computational burden. First of all, in the linearization process, some nonlinear terms are added into the control signal for voltage compensation, which can reduce the order of the prediction model. Then, the dq-axis currents are selected by maximum torque per ampere (MTPA). Six virtual vectors are employed to FCSMPC, and a novel way to calculate the duty cycle is adopted. Finally, the simulation results show the validity and superiority of the proposed method.

1. Introduction

With the gradual exhaustion of resources and the aggravation of global pollution, Electric Vehicles (EVs) have become the future development trend of the automobile industry. The era of Industry 4.0, marked by Cyber-Physical Systems (CPS), promotes the rapid development of EVs. The cyber subsystem is located in the feedback loop of the electric vehicle drive motor. The cyber subsystem drives actuators and the controlled drive motor based on the current states and inputs of the sensors [1,2]. The permanent magnet synchronous motor, especially interior permanent magnet synchronous machines (IPMSM), has become the ideal driving motor for EVs due to its advantages such as high efficiency, high power factor, wide operating speed range, high torque to the current ratio, and robustness.
The nonlinear dynamic of IPMSM, its motor parameter variations, and external disturbances can deteriorate control performances, and proportiona-integral (PI) and linear-quadratic (LQ) regulators cannot achieve good performances. Accordingly, some intelligent control methods have been proposed to improve the control performances, such as robust and adaptive nonlinear approaches, sliding mode control, neural network, fuzzy controllers and genetic algorithms (GAs), backstepping control, and active disturbance rejection control (ADRC).
Robust and adaptive nonlinear approaches have been widely used for control of IPMSM. In [3], an adaptive flux observer is proposed to improve its robustness. In [4], an adaptive interconnected high gain observer has been employed for sensorless control of IPMSM. Moreover, a model reference adaptive system (MRAS) has been employed for speed estimation in the sensorless control of IPMSM [5]. However, adaptive control requires a large number of iterative updates, and real-time control will limit its practical applications.
Sliding mode control has also been utilized in IPMSM control. Higher-order sliding mode (HOSM) techniques have been used to reduce the chattering phenomenon. A hybrid terminal HOSM observer has been devised for position and speed estimation [6]. A super-twisting algorithm (STA) based on HOSM has been used to estimate the rotor position [7]. The designs mentioned above would decrease the control performances due to convergence and robustness problems. Hence, it is necessary to reduce chattering while maintaining robustness. A fractional order sliding mode controller has been considered in previous works. In [8], an adaptive fractional order sliding mode controller for IPMSM has been presented. An intelligent robust fractional surface sliding mode control has been employed for a nonlinear system [9]. The works mentioned above improve transient response and robustness. However, the chattering phenomenon still exists.
Neural-network-based approaches have been presented to estimate the rotor position and speed estimations. A linear neural network named TLS EXIN neuron has been employed for estimation of the speed of the PMSM [10]. In [11], an adaptive linear neural (ADALINE)-network-based filter has been proposed to reduce the harmonic ripple in the rotor position estimation error. However, the approaches are complex to calculate and difficult to implement.
In addition, fuzzy controllers and GAs have been employed in IPMSM. A fuzzy logic speed controller has been designed to reduce torque ripples and flux ripples in [12]. An optimization method, based on GA, for optimal trajectories calculation has been proposed [13]. In [9], GAs have been used to solve the problem of determining and optimizing the parameters of the fuzzy sliding mode controller. However, the fuzzy control has low precision and a steady state error.
Due to its systematic recursive design, backstepping control is a popular control method for IPMSM. In [4], an improved backstepping method has been developed in order to ensure the rotor speed and the rotor current to efficiently track the Benchmark trajectories references. A discrete-time fuzzy position tracking controller using a backstepping approach has been proposed for IPMSM [14]. In the designs mentioned above, computation in real-time would limit control performances.
ADRC is an attractive control technique for IPMSM due to its extended state observer, which could observe and compensate interior model disturbance and exterior load disturbance [15]. It is usually used as the speed controller in the IPMSM drive systems. In [16], an optimal control strategy based on active ADRC has been proposed to achieve high-precision control of PMSM. The controllers mentioned above also improve transient response and robustness.
Moreover, model predictive control (MPC) has received widespread attention because it can respond to the assumed future system behavior in advance and meet system constraints more easily. Combined with the control methods mentioned above, MPC becomes a more powerful and effective control method for IPMSM drive systems.
Finite control set model predictive control (FCSMPC) and continuous control set MPC (CCSMPC) are the two main categories of MPC methods for IPMSM. CCSMPC obtains a set of continuous voltages through an inverter with pulse width modulation (PWM) [17,18]. As opposed to CCSMPC, FCSMPC limits the control space to a finite set consisting of eight voltage vectors and obtains a finite voltage set through the inverter without PWM [19,20]. It has a simple structure, fast dynamic responses, and it is easier to include system constraints. However, FCSMPC usually enumerates and calculates all the candidate voltage vectors and selects an optimal output based on the cost function, which leads to a large amount of calculation and difficulty when applied in practical usage [21,22], reducing the amount of calculation by reducing the candidate switching states for prediction and evaluation. A lookup table has been proposed to reduce the number of candidate voltage vectors from eight to three [23]. In [24], the sector of the desired voltage vector is directly determined, rather than traversing all switch states. Although the approaches mentioned above can reduce the calculation time, they still have large torque ripples and current ripples. Hence, the emphasis of this paper is to decrease the torque ripple and current ripple of FCSMPC.
In this paper, a cascade speed controller based on ADRC and a current controller based on improved FCSMPC for IPMSM are proposed. A different approach to dealing with nonlinear terms in IPMSM is proposed. In the linearization process, some nonlinear terms are added into the control signal for voltage compensation, which can reduce the order of the prediction model. The ADRC controller is used to estimate the nonlinear terms in the mechanical motion equation and generate the appropriate torque reference. Then, the torque reference is converted into the dq-axis current reference of the current controller according to MTPA. Six virtual vectors are employed to FCSMPC, and a novel way to calculate the duty cycle is adopted. The proposed FCSMPC controller applies both current and voltage constraints and reduces torque ripple and current ripple as well as computational burden. Compared with traditional PI control and traditional FCSMPC, the effectiveness and superiority of the designed ADRC controller and improved duty cycle FCSMPC controller are verified.
This paper is organized as follows. Section 2 introduces the mathematical model of IPMSM. Section 3 presents the design details of the ADRC speed controller, and Section 4 presents the design details of the improved duty cycle FCSMPC method in IPMSM. The simulation result is presented in Section 5. The conclusion is given in Section 6.

2. Mathematical of IPMSM

In the three-phase static reference frame, the mathematical model of IPMSM is a nonlinear equation related to the instantaneous position of the rotor, which makes the voltage and flux linkage equations too complicated to be analyzed. The IPMSM established in the three-phase static reference frame (abc) can be converted into a synchronous reference frame (dq) by the Clark transformation and the Park transformation.
The mathematic model of the IPMSM can be described as follows:
d i d ( t ) d t   =       R s L d i d ( t )   +   L q L d ω ( t ) i q ( t )   +   1 L d u d ( t ) d i q ( t ) d t =     L d L q ω ( t ) i d ( t )     R L q i q ( t )     ψ f L q ω ( t )   +   1 L q u q ( t ) d ω ( t ) d t   =       B J ω ( t )   +   p n J ( T e ( t )     T L )
where variables u d , u q , i d ,   i q , L d ,   L q are the d-axis voltage, q-axis voltage, d-axis current, q-axis current, d-axis inductance and q-axis inductance, respectively; variables ω , R , and ψ f are the electrical angular velocity, stator resistance, and rotor flux, respectively; J is the moment of inertia, B is the viscous coefficient of the load, and T L is the disturbance torque, which is given by
T e ( t ) = 3 2 p n i q ( t ) [ ψ f   +   ( L d     L q ) i d ( t ) ]
Since ω e ( t ) , i q ( t ) and i d ( t ) can be measured, ω e ( t ) i q ( t ) and ω e ( t ) i d ( t ) can be considered as compensation for control signals and removed from the original model and then added after solving the control signal:
d i d ( t ) d t =     R s L d i d ( t )   +   1 L d u ¯ d ( t ) d i q ( t ) d t =     R s L q i q ( t )     ψ f L q ω e ( t )   +   1 L q u ¯ q ( t ) d ω ( t ) d t = 3 2 p n 2 J ψ f i q ( t )     B J ω e ( t )     p n J T L
where u ¯ d ( t ) = L q ω e ( t ) i q ( t ) + u d ( t ) , u ¯ q ( t ) = L d ω e ( t ) i d ( t ) + u q ( t ) ;
The state vector of the system can be expressed as:
x = [ i d i q ω e ] T  
Define the control vector as u   =   [ u d , u q ] T , the disturbance as d   =   T L . Then, the state equation is transformed into the following standard form:
x ˙ ( t ) = A c x ( t )   +   B c u ( t )   +   B d c d ( t )  
where,
A c   =   [     R s L d 0 0 0     R s L q     ψ f L q 0 3 p n 2 2 J ψ f     B J ] , B c   =   [ 1 L d 0 0 0 1 L q 0 ] T , B d c   =   [ 0 0     p n J ] T

3. Design of Speed Controller Based on Improved Duty Cycle FCSMPC

3.1. Structure of Control System

The system block diagram of the overall motor-control scheme is shown in Figure 1, which consists of an IPMSM, a three-phase voltage inverter, a coordinate transformation, speed controller based on ADRC and a current controller based on improved duty cycle FCSMPC.
Assuming that ω r e f is the electrical angular velocity reference, T e * is the torque reference, i d * is the d-axis current reference, i q * is the q-axis current reference, i d q is the dq-axis current, i a b c is the abc-axis current and θ is the electrical angle.
The ADRC controller is used to estimate the nonlinear terms in the mechanical motion equation and generate the appropriate torque reference. Then, the torque reference is converted to the dq-axis current reference of the current loop according to MTPA. According to the current reference, the switching signal is generated by the improved duty cycle FCSMPC. Finally, the output voltage signal is used to control IPMSM. The ADRC controller is introduced in Section 3, the MTPA and improved duty cycle FCSMPC are introduced in Section 4.

3.2. Speed Controller Based on ADRC

ADRC controller is usually composed of a tracking differentiator, extended state observer (ESO) and nonlinear error feedback control law. The tracking differentiator has no use in the control of IPMSM. Therefore, only the extended state observer and nonlinear error feedback control law need to be designed.
Part of the speed differential equation can be regarded as the model disturbance of the system, which can be written as:
d ω ( t ) d t = a ( t )   +   b u  
where,
a ( t ) =     B J ω e ( t )     p n J T L , b = p n J , u = T e
Equation (6) contains the uncertainty term a ( t ) and the control term u . If a ( t ) can be correctly observed by ADRC and can be fed back into a known quantity, then the dynamic equation can be simplified into a first-order equation.
The extended state observer corresponding to Equation (6) can be expressed as:
{ e = z 1     ω z ˙ 1 = z 2     β 1 fal ( e , α 1 , δ 1 )   +   b u z ˙ 2 =     β 2 fal ( e , α 2 , δ 1 )
where z 1 is used to observe the rotational speed state, which is used to track the actual speed. z 2 is the expanded state, which represents the estimation of the interference term a ( t ) in Equation (6). u is the control variable, e is the input deviation, δ 1 is the adjustable coefficient greater than 0, and ω is the actual speed. α 1 , α 2 , β 1 and β 2 are the coefficients that need to be adjusted.
fal ( ) in Equation (7) can be expressed as:
fal ( e , α , δ ) = { | e | α sign ( e ) , | e | > δ e / δ ( 1     α ) , | e | δ
where,
sign ( e ) = { 1 , i f e > 0 0 , i f e = 0     1 , i f e < 0
Equation (7) can estimate the speed ω and the interference term a ( t ) , and the first-order nonlinear control law is given as:
u 0 = k 1 fal ( z 1     ω , α 3 , δ 2 )
Substitute Equation (9) to Equation (6) and the final control output can be given as:
u = u 0     z 2 b

4. Design of Current Controller Based on Improved Duty Cycle FCSMPC

4.1. The Maximum Torque Per Ampere (MTPA)

The Maximum Torque Per Ampere control is used to maximize the torque when the amplitude of the current vector is the minimum. Solving the current value under the Maximum Torque Per Ampere control can be equivalently expressed as:
min i d 2   +   i q 2 s . t . T e ( t ) = 3 p n i q ( t ) [ ψ f   +   ( L d     L q ) i d ( t ) ] 2
The problem shown in Equation (11) can be solved by the Lagrange multiplier method. The derivation process of MTPA is not discussed in this paper. By derivation, the relationship between i d and i q can be described as:
i d =     ψ f 2 ( L d     L q )     ψ f 2 4 ( L d     L q ) 2   +   i q 2
The relationship between T e and i q can be described as:
  T e ( t ) = 3 2 p n i q ( ψ f 2     ( L d     L q ) ψ f 2 4 ( L d     L q ) 2   +   i q 2 )
Equation (13) can be solved numerically by the fsolve() function of MATLAB. In order to meet the time limit of online calculation, the i q corresponding to each torque can be obtained in advance according to the interval and then made into a table. Hence, the corresponding i q can be obtained only through the given torque- i q table during online calculation.
Due to the complexity of the solution, calculating i d directly through Equation (12) cannot meet the requirements of real-time applications. Therefore, using Taylor’s formula to expand near a point close to 0 (select 0.001) and ignoring the minimal term and higher-order terms in the expansion:
i d =     1 ψ f 2 ( L d     L q ) 2 ( i q     0.001 ) 2
It is also necessary to consider the issue of constraints. The dq-axis current must satisfy the constraint i d 2   +   i q 2 I N 2 . Therefore, the actual steps are as follows:
set torque T e ;
calculate i q and i d according to Equations (13) and (14);
if i d 2   +   i q 2 I N 2 , loop to Step 5; otherwise, loop to Step 4;
update i d =     sqrt ( I N 2     i q 2 ) ;
output i q and i d as current reference values for the current regulator.

4.2. Current Controller Based on FCSMPC

FCSMPC is based on the finite switching states. This method does not require a modulator such as Space Vector Pulse Width Modulation (SVPWM), and the nonlinear characteristics of the system and various constraints can also be easily integrated into the system.
For the design of FCSMPC, the current prediction model discretized by the first-order Euler discretization method can be expressed as:
[ i d ( k   +   1 ) i q ( k   +   1 ) ] = F ( k ) [ i d ( k ) i q ( k ) ]   +   G [ v d ( k ) v q ( k ) ]   +   H ( k )  
where,
F ( k ) = [ 1     R T L d T ω ( k )     T ω ( k ) 1     R T L q ] , G = [ T L d 0 0 T L q ] , H ( k ) = [ 0     T ω e ( k ) ψ f L q ]
The prediction equation of FCSMPC can be expressed as:
X i ( k   +   1 ) = F ( k ) X _ ( k )   +   G M ( k ) D u _ i ( k )   +   H ( k )
where
M ( k )   =   [ cos θ ( k ) sin θ ( k )     sin θ ( k ) cos θ ( k ) ] , u _ ( k )   =   [ u a ( k ) u b ( k ) u c ( k ) ] ,
D   =   [ v d ( k ) v q ( k ) ]   =   2 3 E M ( k ) [ 1     1 2     1 2 0 3 2     3 2 ] [ u a ( k ) u b ( k ) u c ( k ) ] =   M ( k ) D ( k ) u _ ( k ) ,
where   u _ i ( k ) is the selection of the i -th switch state, and X i ( k   +   1 ) is the current state at the next sampling time after the i -th switch state is applied.
The traditional FCSMPC method traverses all the switch states according to Equation (16), obtains all the predicted states at the next moment, and selects the switch state corresponding to the vector with the smallest difference in the expected voltage vector amplitude for control. Considering the current constraints, the corresponding cost function can be expressed as:
min i = 0 , 1 , 7 J = | X i ( k   +   1 )     X # |   +   flim , flim = { 0 , i f | X i ( k   +   1 ) | I N , i f | X i ( k   +   1 ) | I N

4.3. Improved Duty Cycle FCSMPC

The duty cycle FCSMPC is to add a zero vector in one period, so the optimal vector is only used for a part of the sampling period. Assuming that the ratio of the optimal vector action time to the sampling period final state can be expressed as:
( 1 γ ) X _ 0 ( k + 1 ) + γ X _ s e l e c t ( k + 1 )
define ε X i _ ( k + 1 ) = X _ # X _ i ( k + 1 ) , which represents the deviation of the given state vector X _ # from the i -th vector X _ i ( k   +   1 ) , ε X _ 0 ( k   +   1 )   =   ( x 0 , y 0 ) , which represents the deviation from the given value under the action of the zero vector.
Assuming that the voltage vector V _ i = ( x 1 , y 1 ) , the angle between the two voltage vectors can be expressed as:
co s _ θ = x 0 x 1   +   y 0 y 1 x 0 2   +   y 0 2   +   x 1 2   +   y 1 2  
The optimal vector can be obtained by comparing the size of each vector’s cos θ _ , assuming that the optimal vector is j ; the ideal origin is O ; the future state under the action of zero vector is A, X _ 0 ( k   +   1 ) ; the current state of the system is B, X _ ( k ) ; the given state is C, X _ # ; and the state after its action is D, X _ j ( k   +   1 ) .
The geometric relationship is shown in Figure 2:
If k   +   1 is omitted, then the Equation (18) can be expressed as:
( 1 γ ) X _ 0 + γ X _ j = ( 1 γ ) X _ 0 + γ ( X _ 0 + A D ) = X _ 0 + γ A D
Assuming that the vector of this state is O E , which satisfies A E   =   γ A D , D C = ( x 0 , y 0 ) and A C = ( x 1 , y 1 ) , then
A E   =   γ A D   =   γ ( x 1     x 0 , y 1     y 0 ) E C = A C     A E = ( ( 1     γ ) x 1   +   γ x 0 , ( 1     γ ) y 1   +   γ y 0 )  
Due to C E A E , substitute it into Equation (21) and obtain
γ = x 1 2   +   y 1 2     x 0 x 1     y 0 y 1 ( x 1     x 0 ) 2   +   ( y 1     y 0 ) 2  
However, this method still has its limitations. With only six vectors and two zero vectors to choose from, it cannot approach the possible vectors sufficiently. In order to make the optimal vector more accurate, consider adding 6 virtual voltage vectors. All virtual voltage vectors can be regarded as the result of two adjacent basic vectors acting for half of the time.
First calculate all the basic vectors to select the optimal vector and the sub-optimal vector and then compare with the virtual vector between these two vectors, so that the optimal vector can be obtained in only eight computations.
In addition, due to current constraints, the voltage vector selected according to Equation (19) may cause the stator current to exceed the limit value at the next moment. Therefore, after selecting the optimal vector and calculating the duty cycle. It can be judged according to the following process.
First, the average value of the control voltage in a sampling period can be expressed as:
v ¯ d q ( k ) = 1 T k T ( k   +   1 ) T v d q ( t ) d t  
Define ρ ( k ) _   =   [ ρ a ( k ) ρ b ( k ) ρ c ( k ) ] T , which represents the ratio of the time when each switch state is 1 to the sampling period. Substituting it into Equation (23), the control quantity of a sampling period can be expressed as:
v ¯ d q ( k ) = M ( k ) D ρ ( k )
Substituting Equation (24) into Equation (15), the predicted current value at the next moment can be obtained. If the predicted current exceeds the current limit value, then the sub-optimal vector is selected until the stator current is less than the limit.

5. Simulation Results

The parameters of IPMSM are selected from Toyota Prius Hybrid vehicles, as shown in Table 1.
The relevant parameters of ADRC speed controller are α 1   =   0.8 , δ 1   =   0.001 , β 1   =   2000 , β 2   =   8 × 10 5 , k 1   =   3800 , α 3   =   0.9 .
The simulation condition simulates the driving condition of the car. The car starts with a load of 10 N·m until it reaches the speed of 1000 r/min. The load is increased to 30 N·m at 0.4 s to simulate uphill conditions. After smooth operation, the load is reduced to 10 N·m at 0.6 s to simulate flat road driving. The entering speed is reduced to 500 r/min at 0.8 s, and it enters the deceleration driving mode.
Under the same conditions, simulation experiments were carried out for the control system of the following controllers, as shown in the control block diagram in Figure 1, assuming that C 1 is the speed controller and C 2 is the current controller.
  • Case 1: C 1 uses PI control, and C 2 also uses PI control;
  • Case 2: C 1 uses PI control, and C 2 uses improved duty cycle FCSMPC;
  • Case 3: C 1 uses ADRC control, and C 2 uses traditional FCSMPC;
  • Case 4: C 1 uses ADRC control, and C 2 uses improved duty cycle FCSMPC.
Figure 3 shows the comparison of the speed response in Cases 1, 2 and 4. It indicates that all three controllers can rapidly increase the speed until reaching the given value of 1000 r/min during the starting process with load. The first magnified portion shows the speed response of three controllers as the speed approaches steady state. The peak time in Case 1 is the same as that in Case 2, but the overshoot in Case 2 is smaller than that in Case 1, and the steady-state time is significantly less than that in Case 1. Compared with Case 1 and Case 2, the adjustment time in Case 4 is significantly lower than these two cases, and there is no overshoot. Therefore, it has significant advantages. The second magnified portion shows the speed response when entering the uphill condition at 0.4 s. Due to the increase of load torque, there is a short decrease of rotating speed in all three control cases. Compared with Case 2, the rotational speed in Case 2 decreases less than that in Case 1, and the time for the rotational speed to recover the given value is shorter. Compared with Case 1 and Case 2, the speed drop in Case 4 is the smallest, and it returns to the given value in a short time. After a period of stable operation, the system enters the flat driving condition. The third magnified portion shows the speed response in the three control situations when the load drops. There are larger fluctuations and longer adjustment times in Case 1, the adjustment times in Case 2 and Case 4 are shorter, and the overshoot in Case 4 is smaller. Afterwards, it enters the deceleration condition at 0.8 s. The fourth magnified portion shows the speed response when the target speed is about to reach 500 r/min during deceleration. It indicates that there is no overshoot in Case 4, and the time to reach the steady state is the shortest.
From the comparison of the speeds in the two cases in Figure 4, in acceleration driving condition, the peak time in Case 3 is much longer than in Case 4. Additionally, there are still fluctuations in Case 3 when reaching the steady state. Similarly, the same situation also occurs when the given speed drops. The adjustment time is significantly shorter in Case 4.
Table 2 shows the speed dynamic response in all cases. There is almost no overshoot and transient time in the acceleration and deceleration driving condition using the ADRC controller in Case 3 and Case 4. This proves that the control system using the ADRC controller has a better control effect. Compared with Case 3, Case 4 has less peak time, overshoot and transient time in all conditions. Obviously, compared with the control system using the traditional FCSMPC, the control system using the improved duty cycle FCSMPC has a better control effect.
Figure 5 and Figure 6 show the d-axis current response and the q-axis current response in Case 1, Case 2 and Case 4. When the speed drops, the currents are shown in the first magnified portion of the two figures. In the three cases, the currents are fluctuant in the steady state. The dq-axis currents have very large overshoots in Case 1 and Case 2, while there is almost no overshoot in Case 4. When entering the uphill condition at 0.4 s, the absolute value of the dq-axis currents is increased to offset the increased load torque. At this time, the time to reach the steady state is the shortest in Case 4. The time to reach the steady state in Case 4 is still the shortest when it reaches the flat driving condition at 0.6 s. When entering the deceleration driving condition at 0.8 s, the q-axis current reverses to generate the maximum torque in the negative direction to decelerate rapidly, and then it returns to constant speed driving. The currents decrease so that the output torque is equal to the load torque, and the dq-axis current responses are shown in the fourth magnified portion of the two figures. The dq-axis current has no overshoot in Case 4, and the adjustment time is the shortest. There are large overshoots in Case 1 and Case 2, and the adjustment time is relatively long.
Figure 7 shows that the dq-axis current under the control in Case 3 fluctuates greatly, which will cause the electromagnetic torque value generated to be smaller than that in Case 4. Due to the excessive current fluctuation, it is difficult for it to stabilize when it reaches the steady state. Hence, the speed still fluctuates in the steady state.
Figure 8 shows the torque response in Case 3 and Case 4. In the motor starting phase and the speed drop phase, the maximum torque in Case 3 is less than that in Case 4, and the torque in Case 3 has large torque fluctuations throughout the simulation process.
Figure 9 shows the stator current waveforms and Fourier analysis results in four cases. (a), (b), (c), and (d) respectively represent the results of Case 4, Case 3, Case 2, and Case 1. Total Harmonic Distortion (THD) indicates the degree of waveform distortion caused by all harmonics. The smaller the THD, the lower the degree of distortion between the waveform and the sine wave of a certain frequency. The proposed controller has a THD of 2.25% in Case 4, 17.00% in Case 3, 2.26% in Case 2, and 2.62% in Case 1. The sines of the waveforms in Case 2 and Case 4 are comparable, while the sine is slightly worse in Case 1 and the worst in Case 3. The smoothness of the waveform generated by the traditional FCSMPC is low. While the improved duty cycle FCSMPC is used as the current controller in Case 2 and Case 4, the smoothness of the waveforms is the highest. Therefore, the effectiveness of the proposed current controller can be further proved.

6. Conclusions

This paper presents the application of a cascade speed control strategy based on ADRC and a current control strategy based on improved duty cycle FCSMPC in an IPMSM speed tracking problem. More than six virtual vectors are used to reduce torque ripple and current ripple, and this novel way to select the optimal vector leads to a reduction in the computational burden. The proposed controllers are investigated at different driving conditions with current and voltage constraints. Compared with traditional PI control and traditional FCSMPC, the simulation results indicate that the proposed controllers exhibit better dynamic performances keeping all system variables within rated. Additionally, the THD and current ripple is lower than that of PI control and traditional FCSMPC. Our future work will be extending the results to other potential applications including energy storage management.

Author Contributions

Conceptualization, F.L. and H.L.; methodology, F.L., L.L. and R.Z.; software, F.L.; validation, F.L., R.Z. and K.L.; formal analysis, F.L. and H.L.; investigation, F.L., H.L. and L.L.; resources, F.L., H.L. and R.Z.; data curation, F.L. and H.L.; writing—original draft preparation, F.L. and H.L.; writing—review and editing, F.L. and H.L.; visualization, L.L. and K.L.; supervision, H.L.; project administration, H.L.; funding acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 61973318, and in part by Distinguished Young Foundation of Hunan Province of China under Grant 2020JJ2045.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, S.; Koren, I.; Krishna, C.M. Adaptive Workload Adjustment for Cyber-Physical Systems using Deep Reinforcement Learning. Sustain. Comput. Inform. Syst. 2021, 100525, ISSN 2210-5379. [Google Scholar] [CrossRef]
  2. Tehrani, K. A smart cyber physical multi-source energy system for an electric vehicle prototype. J. Syst. Archit. 2020, 111, 101804. [Google Scholar] [CrossRef]
  3. Sinetova, M.; Pyrkin, A.; Bobtsov, A.; Ortega, R.; Vedyakov, A. An Adaptive Flux and Position Observer for Interior Permanent Magnet Synchronous Motors. IFAC-PapersOnLine 2019, 52, 43–48. [Google Scholar] [CrossRef]
  4. Hamida, M.A.; de Leon, J.; Glumineau, A. Experimental sensorless control for IPMSM by using integral backstepping strategy and adaptive high gain observer. Control Eng. Pract. 2017, 59, 64–76. [Google Scholar] [CrossRef]
  5. Khlaief, A.; Boussak, M.; Châari, A. A MRAS-based stator resistance and speed estimation for sensorless vector controlled IPMSM drive. Electr. Power Syst. Res. 2014, 108, 1–15. [Google Scholar] [CrossRef]
  6. Zheng, J.; Feng, Y.; Yu, X. Hybrid terminal sliding mode observer design method for permanent magnet synchronous motor control system. In Proceedings of the 2008 International Workshop on Variable Structure Systems, Antalya, Turkey, 8–10 June 2008; pp. 106–111. [Google Scholar] [CrossRef]
  7. Moreno, J.A.; Osorio, M. Strict Lyapunov Functions for the Super-Twisting Algorithm. IEEE Trans. Autom. Control 2012, 57, 1035–1040. [Google Scholar] [CrossRef]
  8. Hosseini, S.H.; Tabatabaei, M. IPMSM velocity and current control using MTPA based adaptive fractional order sliding mode controller. Eng. Sci. Technol. Int. J. 2017, 20, 896–908. [Google Scholar] [CrossRef]
  9. Delavari, H.; Ghaderi, R.; Ranjbar, A.; Momani, S. Fuzzy fractional order sliding mode controller for nonlinear systems. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 963–978. [Google Scholar] [CrossRef]
  10. Accetta, A.; Cirrincione, M.; Pucci, M. TLS EXIN based neural sensorless control of a high dynamic PMSM. Control Eng. Pract. 2012, 20, 725–732. [Google Scholar] [CrossRef]
  11. Zhang, G.; Wang, G.; Xu, D.; Zhao, N. ADALINE-Network-Based PLL for Position Sensorless Interior Permanent Magnet Synchronous Motor Drives. IEEE Trans. Power Electron. 2016, 31, 1450–1460. [Google Scholar] [CrossRef]
  12. Mehedi, F.; Djilali, A.B.; Yahdou, A.; Bouyakoub, I. Fuzzy Logic Based DTC-SVM for Speed Control of Five-Phase IPMSM. In Proceedings of the 2019 7th International Renewable and Sustainable Energy Conference (IRSEC), Agadir, Morocco, 27–30 November 2019; pp. 1–6. [Google Scholar] [CrossRef]
  13. Kumagai, S.; Inoue, K.; Kato, T. GA Based Optimized Trajectories of Rotating Speed and d-q Axis Currents for an IPMSM. In Proceedings of the 2018 International Power Electronics Conference (IPEC-Niigata 2018 -ECCE Asia), Niigata, Japan, 20–24 May 2018; pp. 1264–1269. [Google Scholar] [CrossRef]
  14. Yu, J.; Shi, P.; Yu, H.; Chen, B.; Lin, C. Approximation-Based Discrete-Time Adaptive Position Tracking Control for Interior Permanent Magnet Synchronous Motors. IEEE Trans. Cybern. 2015, 45, 1363–1371. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, L.; Zhang, H.; Liu, X. Robust sensorless of ADRC controlled PMSM based on MRAS with stator resistance identification. In Proceedings of the 30th Chinese Control Conference, Yantai, China, 22–24 July 2011; pp. 3575–3579. [Google Scholar]
  16. Meng, Y.; Liu, B.; Wang, L. Speed Control of PMSM Based on an Optimized ADRC Controller. Math. Probl. Eng. 2019, 2019. [Google Scholar] [CrossRef]
  17. Errouissi, R.; Al-Durra, A.; Muyeen, S.; Leng, S. Continuous-time model predictive control of a permanent magnet synchronous motor drive with disturbance decoupling. IET Electr. Power Appl. 2017, 11, 697–706. [Google Scholar] [CrossRef]
  18. Wang, W.; Fan, Y.; Chen, S.; Zhang, Q. Finite control set model predictive current control of a five-phase PMSM with virtual voltage vectors and adaptive control set. CES Trans. Electr. Mach. Syst. 2018, 2, 136–141. [Google Scholar] [CrossRef]
  19. Lang, X.; Yang, M.; Long, J.; Li, Y.; Xu, D. A novel direct predictive speed and current controller for PMSM drive. In Proceedings of the 2016 IEEE 8th International Power Electronics and Motion Control Conference (IPEMC-ECCE Asia), Hefei, China, 22–26 May 2016; pp. 2551–2556. [Google Scholar] [CrossRef]
  20. Xie, W.; Wang, X.; Wang, F.; Xu, W.; Kennel, R.M.; Gerling, D.; Lorenz, R.D. Finite-Control-Set Model Predictive Torque Control with a Deadbeat Solution for PMSM Drives. IEEE Trans. Ind. Electron. 2015, 62, 5402–5410. [Google Scholar] [CrossRef]
  21. Li, Y.; Zhao, C.; Zhou, Y.; Qin, Y. Model predictive torque control of PMSM based on data drive. Energy Rep. 2020, 6, 1370–1376. [Google Scholar] [CrossRef]
  22. Siami, M.; Khaburi, D.A.; Rodriguez, J. Simplified Finite Control Set-Model Predictive Control for Matrix Converter-Fed PMSM Drives. IEEE Trans. Power Electron. 2018, 33, 2438–2446. [Google Scholar] [CrossRef]
  23. Habibullah, M.; Lu, D.D.; Xiao, D.; Rahman, M.F. A Simplified Finite-State Predictive Direct Torque Control for Induction Motor Drive. IEEE Trans. Ind. Electron. 2016, 63, 3964–3975. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, W.; Yang, W.; Yan, W.; Xu, D. Improved Finite Control Set Model Predictive Control for Permanent Magnet Synchronous Motor Drives. In Proceedings of the 2019 Chinese Control Conference (CCC), Guangzhou, China, 27–30 July 2019; pp. 2977–2982. [Google Scholar] [CrossRef]
Figure 1. The block diagram of interior permanent magnet synchronous machines (IPMSM) cascade control strategy based on active disturbance rejection control (ADRC) and improved finite control set model predictive control (FCSMPC).
Figure 1. The block diagram of interior permanent magnet synchronous machines (IPMSM) cascade control strategy based on active disturbance rejection control (ADRC) and improved finite control set model predictive control (FCSMPC).
Mathematics 09 00760 g001
Figure 2. The solving diagram of duty cycle FCSMPC.
Figure 2. The solving diagram of duty cycle FCSMPC.
Mathematics 09 00760 g002
Figure 3. The comparison of the speed response in Case 1, Case 2 and Case 4.
Figure 3. The comparison of the speed response in Case 1, Case 2 and Case 4.
Mathematics 09 00760 g003
Figure 4. The speed response in Case 3 and Case 4.
Figure 4. The speed response in Case 3 and Case 4.
Mathematics 09 00760 g004
Figure 5. The q-axis current response in Case 1, Case 2 and Case 4.
Figure 5. The q-axis current response in Case 1, Case 2 and Case 4.
Mathematics 09 00760 g005
Figure 6. The d-axis current response in Case 1, Case 2 and Case 4.
Figure 6. The d-axis current response in Case 1, Case 2 and Case 4.
Mathematics 09 00760 g006
Figure 7. The dq-axis current response in Case 3 and Case 4.
Figure 7. The dq-axis current response in Case 3 and Case 4.
Mathematics 09 00760 g007
Figure 8. The torque response in Case 3 and Case 4.
Figure 8. The torque response in Case 3 and Case 4.
Mathematics 09 00760 g008
Figure 9. (ad) The stator current waveforms and Fourier analysis results in four cases.
Figure 9. (ad) The stator current waveforms and Fourier analysis results in four cases.
Mathematics 09 00760 g009aMathematics 09 00760 g009b
Table 1. The Toyota Prius Hybrid vehicles’ IPMSM parameters.
Table 1. The Toyota Prius Hybrid vehicles’ IPMSM parameters.
DescriptionValueSymbol
ψ f (Rotor flux)0.035Wb
L d (d-axis inductance)0.169mH
L q (q-axis inductance)0.331mH
R s (Resistance of stator)0.07 Ω
J (Moment of inertia)0.1312 kg m 2
p n (Poles)4
U N (Rated voltage)500V
I N (Rated current)250A
f N (Rated frequency)200Hz
Table 2. Comparison of speed dynamic response in all cases.
Table 2. Comparison of speed dynamic response in all cases.
Driving ConditionCasePeak Time (s)Overshoot (r/min)Transient Time (s)
AccelerationCase 10.22394.09260.2335
Case 20.22492.28560.2286
Case 30.30630.1175 (almost 0)-(null)
Case 40.22470.0077 (almost 0)-(null)
Uphill roadCase 10.00221.62170.0051
Case 20.00150.90970.0034
Case 30.00090.62100.0017
Case 40.00090.58110.0014
Flat roadCase 10.00231.56720.0066
Case 20.00160.89140.0036
Case 30.00110.70190.0018
Case40.00090.57530.0014
DecelerationCase 10.08655.75430.0957
Case 20.08613.26250.0901
Case 30.10550.1126 (almost 0)-(null)
Case 40.08550.0114 (almost 0)-(null)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, F.; Li, H.; Liu, L.; Zou, R.; Liu, K. A Control Method for IPMSM Based on Active Disturbance Rejection Control and Model Predictive Control. Mathematics 2021, 9, 760. https://doi.org/10.3390/math9070760

AMA Style

Liu F, Li H, Liu L, Zou R, Liu K. A Control Method for IPMSM Based on Active Disturbance Rejection Control and Model Predictive Control. Mathematics. 2021; 9(7):760. https://doi.org/10.3390/math9070760

Chicago/Turabian Style

Liu, Fang, Haotian Li, Ling Liu, Runmin Zou, and Kangzhi Liu. 2021. "A Control Method for IPMSM Based on Active Disturbance Rejection Control and Model Predictive Control" Mathematics 9, no. 7: 760. https://doi.org/10.3390/math9070760

APA Style

Liu, F., Li, H., Liu, L., Zou, R., & Liu, K. (2021). A Control Method for IPMSM Based on Active Disturbance Rejection Control and Model Predictive Control. Mathematics, 9(7), 760. https://doi.org/10.3390/math9070760

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop