Fuzzy-Energy-Management-Based Intelligent Direct Torque Control for a Battery—Supercapacitor Electric Vehicle

: This paper presents a proposed fuzzy energy management strategy developed for a battery − super capacitor electric vehicle. In addition to providing different driving modes, the proposed strategy delivers the suitable type and amount of power to the vehicle. Furthermore, the proposed strategy takes into account possible failures in vehicle power sources. The speed and torque of the HEV traction machine are simultaneously controlled using a genetic algorithm that provides simultaneous tuning via the use of newly proposed cost functions that give the designer the ability to tradeoff and prioritize between the design variables to be minimized. The simulation results show that the intelligent speed and torque control and the fuzzy power management strategy improved the vehicle’s performance in terms of ripple minimization. The real-time simulation is conducted using the RT LAB simulator, and the results obtained correspond to those obtained in the numerical simulation using MATLAB/Simulink.


Introduction
The road transportation industry depends on oil, an exhausting source of polluting energy. To reduce or eliminate this dependence, researchers have used electric motors to replace the conventional EV traction chain, which relies on an internal combustion engine. Several research axes were born out of the increasing importance given to electric vehicles, such as control of electric propulsion systems [1,2], energy management [3], vehicular communication systems [4], autonomous vehicles [5], vehicle chassis systems [6], traction motor control [7], power electronics [8,9], power train architecture [10], vehicle safety [11], fault tolerance [12], and grid connected electric vehicles [13], which are among the biggest trends nowadays in automotive engineering.
During the last decade, several bioinspired algorithms, such as the genetic algorithm [14], particle swarm optimization [15], ant colony optimization [16], and adaptive quantuminspired evolutionary algorithm [17], have been used in several engineering disciplines. These algorithms are employed because of their ability to automatically find optimal and near-optimal solutions in contrast with classical methods, which are tedious and timeconsuming. While an HEV is undergoing a driving cycle, its torque and speed must be precisely controlled. In this paper, the DTC technique was chosen because of its fast torque response, which is highly recommended in HEV applications. PMSM ensures HEV traction because of its high efficiency, high torque to mass ratio, and fast recovery time during transitions. PMSM flux and torque are measured using an estimator that optimizes hardware as it requires fewer sensors than the classical estimator. The classical DTC technique requires the designer to fix the bounds of two hysteresis controllers and to tune one PI controller. DTC-SVM needs three PI controllers to be adjusted. Both of the previously mentioned control techniques force the user to use classical tuning methods or trial and error. When a fuzzy logic controller is used with DTC, the designer must select appropriate membership functions and a good system. It is worth noting that classical methods are too approximate. The authors in [18] used a classical method for finding controller gains. However, in their design approach, they assumed only a linearized plant in series with a PI controller, neglecting many nonlinear components, such as inverter and hysteresis controllers. Usually, the gains obtained from this method are not consistent.
A battery−supercapacitor is a good power source combination for use in HEV. First, both sources are rechargeable, which is an important feature that allows for energy recovery. Secondly, these two energy storage sources have a complementary characteristic, which makes their combination reduce the drawbacks of each source. For example, the discharge time of supercapacitors is very small compared with the battery discharge time. Supercapacitors have a low power density, whereas batteries have a high energy density. Supercapacitors have higher charge and discharge cycles in contrast with batteries [19]. In [20,21], an electric vehicle powered using battery and fuel cell is studied. In these works, the HEV was simulated for only a few seconds. This short simulation period does not allow for rating the effectiveness of the power management strategy, because the power sources' state of charge will not be affected in a short period. In [22], the performance and lifetime of a hybrid electric vehicle is studied. However, the vehicle in that work is controlled using field-oriented control, which is very sensitive to parameter variation. In [23], improvement of the energy management strategy of hybrid electric vehicle is studied. Even though advanced artificial intelligence techniques are used, important power peaks in the order of kilowatts are still noted. The authors in [24] discussed energy management and faulttolerant strategies for a fuel cell−supercapacitor HEV. The energy management proposed in that work resulted in a relatively deep discharge of the supercapacitor. In [25], the authors resumed the HEV system to design a vehicle speed controller into a first order transfer function, neglecting the presence of nonlinear devices such as AC machines, inverters, and hysteresis controllers. In that work, the PMSM stator flux and torque were estimated to be very sensitive to DC components, resulting in very poor estimations at low speeds.
In this paper, the HEV powered using a battery and supercapacitor is controlled using DTC-based space vector modulation. The PMSM stator flux and torque are measured using a robust estimator that requires a minimum number of sensors. An approach is proposed that makes it possible to adapt the genetic algorithm to HEV. In addition to self-tuning HEV regulators, this approach allows the designer to tradeoff or prioritize between the variables via an adequate proposed cost function. The use of vehicle power sources is managed by a proposed fuzzy power management strategy, which is intended to reduce the power source switching and hence reduce the DC bus ripples associated with each switching. A fuzzy set of rules is designed to provide the HEV with the quantity and quality of power required for traction. The effectiveness of the proposed power management strategy is verified through an extended 660 s simulation period.
This paper is structured as follows: Section 2 deals with electric vehicle modeling and presents a detailed analysis of the different forces exerted on the vehicle as it undergoes a driving cycle. The direct torque control technique used to control the torque of the PMSM, ensuring vehicle traction, is detailed in Section 3. Section 4 points out how a genetic algorithm is adapted for use with the electric vehicle system. Furthermore, the different steps that precede the formulation of the cost function are explained. Section 5 highlights the proposed fuzzy power management strategy used in this paper to control the battery and supercapacitor power flow. Section 6 concludes this paper and confirms the effectiveness of the proposed intelligent direct torque control technique and fuzzy power management using numerical simulation and a RT LAB simulator.

Hybrid Electric Vehicle Dynamics
The motion of an EV is governed by Newton's second law, given below: where F ext is the sum of all of the external forces acting on the EV, m represents vehicle mass, and "a" is the EV acceleration.
Equation (1) can be rewritten as the difference between the EV traction force F traction and the net resistive force Fr, as it is shown in the equation below: The net resistive force acting on the EV while it is undergoing a given driving cycle is given by the following equation: where F aero is the aerodynamic force, F g is the road slope force, F tire is the friction force between the vehicle tires and ground surface, and F acc represents the acceleration force required to achieve maximum vehicle speed from rest. The aerodynamic force to which the hybrid electric vehicle is subjected is given by Equation (4), as shown below [26]: where ρ represents air density, A stands for EV front area, C d represents the aerodynamic coefficient, V wheel is the EV longitudinal speed, and V wind is the wind speed. The slope force due to road inclination is given by Equation (5), where g represents gravity acceleration and α represents the road inclination angle.
The friction force between the vehicle tires and the ground surface is given by Equation (6), which predicts the values of the rolling resistance force, with acceptable accuracy for speeds up to 128 km/h [27]: Acceleration force is modeled by Equation (7), where ta is the EV time constant characteristic. The data of the HEV used in this work are shown in Table 1. The parameters of the PMSM used for HEV traction are presented in Table 1. This machine is chosen because of its robustness and its high torque-to-mass ratio. Its stator flux components and developed electromagnetic torque are given by Equations (8) and (9), respectively.

DTC Algorithm
When the driver pushes the accelerator pedal, they give the motor a torque reference to follow. As DTC provides a fast and precise torque response, it was chosen from among other control techniques in order to control the EV in this work. DTC is a vector control technique that was developed by Takahashi and M. Depenbrock in the 1980s. In classic DTC, the stator flux in the αβ frame, stator flux magnitude, electromagnetic torque, and stator flux angle are estimated using Equations (10)-(13), respectively.
V sα and V sβ are calculated using Equation (14), as shown below: From the above equations, it can be seen that the estimator used in classical DTC requires three sensors for flux and torque measurements (i a and i b current sensors and V DC voltage sensor). In addition, the position sensor is needed to sense and control the EV speed. So, the total number of required sensors is equal to four. In this work, an estimator that requires less sensors and that does not involve any integral terms was used, and it is shown in the figure below ( Figure 1).
To implement DTC with a speed loop, the torque and flux estimator shown above requires only three sensors, i a , i b , and θ (IC is deduced from i a , i b ), whereas the classical estimator needs another extra sensor for the DC bus voltage measurement, as already stated in Equation (14).
GA is a stochastic universal adaptive search optimization technique based on inherent selection [28][29][30][31]. In this algorithm, the chromosomes will mate through a process called crossover, producing offspring that combine their parent's genes. In generation, few chromosomes will undergo mutation in their genes. The number of chromosomes that will undergo crossover and mutation is controlled via crossover and mutation rates. The basic steps of the GA algorithm are highlighted in Figure 2, and the GA parameters used are shown in Table 2. To implement DTC with a speed loop, the torque and flux estimator shown above requires only three sensors, ia, ib, and θ (IC is deduced from ia, ib), whereas the classical estimator needs another extra sensor for the DC bus voltage measurement, as already stated in Equation (14).
GA is a stochastic universal adaptive search optimization technique based on inherent selection [28][29][30][31]. In this algorithm, the chromosomes will mate through a process called crossover, producing offspring that combine their parent's genes. In generation, few chromosomes will undergo mutation in their genes. The number of chromosomes that will undergo crossover and mutation is controlled via crossover and mutation rates. The basic steps of the GA algorithm are highlighted in Figure 2, and the GA parameters used are shown in Table 2. Table 2. GA parameters used.

Parameter Population Iterations Selection Selection Rate Mutation Rate
Value 100 50 Roulette wheel 0.5 0.5 Figure 2. Convergence of test function over iterations.  To implement DTC with a speed loop, the torque and flux estimator shown above requires only three sensors, ia, ib, and θ (IC is deduced from ia, ib), whereas the classical estimator needs another extra sensor for the DC bus voltage measurement, as already stated in Equation (14).
GA is a stochastic universal adaptive search optimization technique based on inherent selection [28][29][30][31]. In this algorithm, the chromosomes will mate through a process called crossover, producing offspring that combine their parent's genes. In generation, few chromosomes will undergo mutation in their genes. The number of chromosomes that will undergo crossover and mutation is controlled via crossover and mutation rates. The basic steps of the GA algorithm are highlighted in Figure 2, and the GA parameters used are shown in Table 2. Table 2. GA parameters used.

Parameter Population Iterations Selection Selection Rate Mutation Rate
Value 100 50 Roulette wheel 0.5 0.5 Figure 2. Convergence of test function over iterations. Figure 2. Convergence of test function over iterations. It is worth noting that the convergence of the GA algorithm used was tested using the benchmark test function [32][33][34][35] shown in Figure 3, whose expression is given by Equation (15). This benchmark function is full of local minima and has a unique global optimum, which may trap the algorithm if it is not well conceived. Figure 2 shows the result of minimizing Equation (15) using the GA algorithm. It can be seen that the GA algorithm succeeded in avoiding all of the local optimums and converged to the same global optimal solution already found antically. result of minimizing Equation (15) using the GA algorithm. It can be seen that the GA algorithm succeeded in avoiding all of the local optimums and converged to the same global optimal solution already found antically.

Cost Function and GA Adaptation
The global HEV scheme is shown in Figure 4. The outer controller PI1 is used for speed regulation and the two internal regulators, PI2 and PI3, are dedicated to torque and flux regulation. Two power sources-the battery and supercapacitor-feed the DC bus inverter through a bidirectional DC−DC converter. Information about the power sources used are highlighted in Table 3.

Cost Function and GA Adaptation
The global HEV scheme is shown in Figure 4. The outer controller PI1 is used for speed regulation and the two internal regulators, PI2 and PI3, are dedicated to torque and flux regulation. Two power sources-the battery and supercapacitor-feed the DC bus inverter through a bidirectional DC−DC converter. Information about the power sources used are highlighted in Table 3.  The objective is to minimize the torque, speed, and flux ripples simultaneously. So, the best solution is the one that significantly minimizes all of the previously mentioned quantities. In any optimization algorithm, the cost function is responsible for selecting the population's fittest particles. In this work, the cost function given by Equation (16) is used along with GA to simultaneously minimize the speed torque and flux errors.
α, β, and γ are weighting factors and their sum is equal to 1. These weighting factors represent the importance attributed to the different design variables. For example, torque ripples will be minimized more than speed and flux ripples if their corresponding weighting factor is greater than the other factors. The proposed cost function also allows the designer to tradeoff between design variables by putting α = β = γ. The average absolute speed error, average absolute torque error, and average absolute flux error are calculated using Equations (17)- (19), respectively: ^ n  The objective is to minimize the torque, speed, and flux ripples simultaneously. So, the best solution is the one that significantly minimizes all of the previously mentioned quantities. In any optimization algorithm, the cost function is responsible for selecting the population's fittest particles. In this work, the cost function given by Equation (16) is used along with GA to simultaneously minimize the speed torque and flux errors.
α, β, and γ are weighting factors and their sum is equal to 1. These weighting factors represent the importance attributed to the different design variables. For example, torque ripples will be minimized more than speed and flux ripples if their corresponding weighting factor is greater than the other factors. The proposed cost function also allows the designer to tradeoff between design variables by putting α = β = γ. The average absolute speed error, average absolute torque error, and average absolute flux error are calculated using Equations (17)- (19), respectively: where Ω i and T i represent the HEV reference speed and torque, respectively. T i and Ω i represent the measured HEV speed and torque, respectively. ϕ ref is the reference flux and it is equal to the permanent magnet flux, and ϕ i is the measured PMSM flux at a given instant, n represents the length of speed, and torque absolute error vectors. Figure 5 depicts the genetic algorithm pseudo code used in this paper. It is worth noting that while the HEV is running, the absolute error between its speed and the reference speed is calculated after each sampling time using Equation (17), as shown in Figure 6 (same applies for torque and flux errors). instant, n represents the length of speed, and torque absolute error vectors. Figure 5 depicts the genetic algorithm pseudo code used in this paper. It is worth noting that while the HEV is running, the absolute error between its speed and the reference speed is calculated after each sampling time using Equation (17), as shown in Figure 6 (same applies for torque and flux errors).  It is worth noting that while the HEV is running, the absolute error between its speed and the reference speed is calculated after each sampling time using Equation (17), as shown in Figure 6 (same applies for torque and flux errors).     The internal structure of the fuzzy power management block shown in the figure above is detailed in Figure 8. This block fuzzifies its inputs using adequate membership functions, and then outputs the supercapacitor reference power based on a set of welldefined rules. It is worth noting that establishing rules governing the fuzzy logic controller's operation is critical for good power management. The internal structure of the fuzzy power management block shown in the figure above is detailed in Figure 8. This block fuzzifies its inputs using adequate membership functions, and then outputs the supercapacitor reference power based on a set of well-defined rules. It is worth noting that establishing rules governing the fuzzy logic controller's operation is critical for good power management. The internal structure of the fuzzy power management block shown in the figure above is detailed in Figure 8. This block fuzzifies its inputs using adequate membership functions, and then outputs the supercapacitor reference power based on a set of welldefined rules. It is worth noting that establishing rules governing the fuzzy logic controller's operation is critical for good power management.     The fuzzy logic controller represents the SC reference power. The membership shown in Figure 11 illustrates the output fuzzy sets used for deffuzification.   The fuzzy logic controller represents the SC reference power. The membership shown in Figure 11 illustrates the output fuzzy sets used for deffuzification. The fuzzy logic controller represents the SC reference power. The membership shown in Figure 11 illustrates the output fuzzy sets used for deffuzification. The fuzzy logic controller represents the SC reference power. The membership shown in Figure 11 illustrates the output fuzzy sets used for deffuzification.

Simulation and Results
The global HEV simulation scheme built in the MATLAB/Simulink environment is shown in Figure 12. The outer controller is used for speed regulation and the two internal regulators are dedicated to torque and flux regulation. It can be seen that electric vehicle wheels are driven using PMSM, which is subjected to the different load torque already presented in Section 2.

Simulation and Results
The global HEV simulation scheme built in the MATLAB/Simulink environment is shown in Figure 12. The outer controller is used for speed regulation and the two internal regulators are dedicated to torque and flux regulation. It can be seen that electric vehicle wheels are driven using PMSM, which is subjected to the different load torque already presented in Section 2. When GA is used for regulator tuning, a population of chromosomes will explore the search space in order to find the adequate PI combination that minimizes the cost function. Figure 13a,c,e shows the distribution of the first generation of potential solutions in the search space. It can be seen that the search space dimension of each controller is equal to two (two axes), and the coordinates of each potential solution in the search space represents a possible PI combination, which may result in the minimization of the cost function, given by Equation (16). Figure 13b,d,f shows that the last generation of chromosomes generated by GA for each PI controller has converged to an almost unique location, which represents the optimum PI gains. From these three last-mentioned figures, When GA is used for regulator tuning, a population of chromosomes will explore the search space in order to find the adequate PI combination that minimizes the cost function. Figure 13a,c,e shows the distribution of the first generation of potential solutions in the search space. It can be seen that the search space dimension of each controller is equal to two (two axes), and the coordinates of each potential solution in the search space represents a possible PI combination, which may result in the minimization of the cost function, given by Equation (16). Figure 13b,d,f shows that the last generation of chromosomes generated by GA for each PI controller has converged to an almost unique location, which represents the optimum PI gains. From these three last-mentioned figures, one can state that the selected optimum (Kp, Ki) combinations for PI1, PI2, and PI3 are (119.2, 20.3), (25.03, 24.97), and (7.64 × 10 5 , 15.77), respectively. The speed of the HEV is presented in Figure 14. One can notice that HEV speed precisely follows its reference, which contains slow, medium, and aggressive acceleration and decelerations.     The electromagnetic torque developed by the HEV traction machine is shown in Figure 15. It can be seen that the PMSM torque response is fast and approximately matches the reference torque. From this figure, it can be seen that the traction machine operates in both the motor and generator mode, as the torque has positive and negative values.   A magnification of the HEV power, speed, and torque during the period [0; 70 s] is provided in Figure 16. It can be seen that the DTC-SVM-based GA resulted in a narrow PMSM torque hysteresis band of 1 N·m, which represents almost one-tenth of the hysteresis band obtained with a two-level inverter in [22]. It can be seen that the battery and the supercapacitor-generated power is equal to the power developed by the PMSM for traction. Note that the power sources and motor ripples are very narrow compared with [22,23].  Figure 16. It can be seen that the DTC-SVM-based GA resulted in a narrow PMSM torque hysteresis band of 1 N·m, which represents almost one-tenth of the hysteresis band obtained with a two-level inverter in [22]. It can be seen that the battery and the supercapacitor-generated power is equal to the power developed by the PMSM for traction. Note that the power sources and motor ripples are very narrow compared with [22,23].  Figure 17 shows the PMSM stator flux on the alpha-beta frame. The circular shape means that the flux is constant in the time domain and its amplitude is equal to the circle radius, which is equal to the reference flux. The magnification in Figure 17 confirms that the PMSM flux follows its reference, which is equal to 0.192 Wb.  Figure 17 shows the PMSM stator flux on the alpha-beta frame. The circular shape means that the flux is constant in the time domain and its amplitude is equal to the circle radius, which is equal to the reference flux. The magnification in Figure 17 confirms that the PMSM flux follows its reference, which is equal to 0.192 Wb. Figure 18 shows that the inverter DC bus voltage follows its reference value, which is equal to 560 V. The magnified area on the same figure shows that the ripples band is less than or equal to 10 V, and this is twice as small as what was obtained in [20]. This is because the established fuzzy rules suppressed many of the power source switchings.  Figure 17 shows the PMSM stator flux on the alpha-beta frame. The circular shape means that the flux is constant in the time domain and its amplitude is equal to the circle radius, which is equal to the reference flux. The magnification in Figure 17 confirms that the PMSM flux follows its reference, which is equal to 0.192 Wb.   Figure 18 shows that the inverter DC bus voltage follows its reference value, which is equal to 560 V. The magnified area on the same figure shows that the ripples band is less than or equal to 10 V, and this is twice as small as what was obtained in [20]. This is because the established fuzzy rules suppressed many of the power source switchings.      Figure 19 confirms the efficiency of the power management strategy. One can see that the power provided by the battery and/or supercapacitor is equal to the power developed by the HEV at every simulation instant. Furthermore, the magnified area on the same figure shows that the power has a very tolerable hysteresis band of almost 200 W.  Figure 18 shows that the inverter DC bus voltage follows its reference value, which is equal to 560 V. The magnified area on the same figure shows that the ripples band is less than or equal to 10 V, and this is twice as small as what was obtained in [20]. This is because the established fuzzy rules suppressed many of the power source switchings.     Figure 20 shows the battery and the supercapacitor power. At each instant, one can notice that their sum is equal to the reference HEV required power shown in the figure below. In addition, this figure clearly shows the different HEV operation modes, namely: the traction mode using the battery only, the traction mode using the supercapacitor only, the traction using both sources, supercapacitor charging, battery charging, and the charging of both power sources.  The restored power will be used to charge the battery, because its state of charge (SOC) is lower than that of the supercapacitor, as shown in Figure 22. From the previous figure, it can be seen that even though the SOC of the two power sources is initially different, their final SOC is nearly the same, which reflects good power management. From the figure, it can also be seen that the supercapacitor and/or battery are charging during the regenerative braking periods.  The restored power will be used to charge the battery, because its state of charge (SOC) is lower than that of the supercapacitor, as shown in Figure 22. From the previous figure, it can be seen that even though the SOC of the two power sources is initially different, their final SOC is nearly the same, which reflects good power management. From the figure, it can also be seen that the supercapacitor and/or battery are charging during the regenerative braking periods.

Real Time Simulation
The proposed soft transition strategy was executed in real-time using the OPAL-RT LAB simulator. As shown in Figure 23, the HEV simulation system was first split into a master computation bloc (SM_computation) and user interface console (SC_GUI) bloc, which show the effect control signal variations. The computation subsystem was executed in real-time on one of the CPU cores of the real-time target. Data between the computation subsystem and GUI subsystem were exchanged asynchronously through the TCP/IP link. Figure 24 shows the real time simulation bench established in our research laboratory, in which element 1 represents the host PC, element 2 shows the FPGA-based real time simulator (OP 5700), element 3 highlights the unit measurement of the data acquisition interface (OP 8660), and element 4 is the oscilloscope.

Real Time Simulation
The proposed soft transition strategy was executed in real-time using the OPAL-RT LAB simulator. As shown in Figure 23, the HEV simulation system was first split into a master computation bloc (SM_computation) and user interface console (SC_GUI) bloc, which show the effect control signal variations. The computation subsystem was executed in real-time on one of the CPU cores of the real-time target. Data between the computation subsystem and GUI subsystem were exchanged asynchronously through the TCP/IP link. Figure 24 shows the real time simulation bench established in our research laboratory, in which element 1 represents the host PC, element 2 shows the FPGA-based real time simulator (OP 5700), element 3 highlights the unit measurement of the data acquisition interface (OP 8660), and element 4 is the oscilloscope.   Figure 25 shows the load power required by the HEV, as well as the power developed by the battery and supercapacitor. It can be seen that the power developed by the HEV power sources follows its reference precisely. Furthermore, it is worth noting that a simple load power profile was chosen to clearly show each power source's power contribution.   Figure 25 shows the load power required by the HEV, as well as the power developed by the battery and supercapacitor. It can be seen that the power developed by the HEV power sources follows its reference precisely. Furthermore, it is worth noting that a simple load power profile was chosen to clearly show each power source's power contribution.  Figure 26 shows the power developed by each power source along with the FC control signal that represents the amount of total power that FC delivers at each time instant. From Figure 24, at P1, when = 0, FC is off and delivers no power. At P2 = 0.3 and this will make the FC delivers 30% of the required power for traction and the SC will ensure the remaining power. At P2 = 1 meaning that FC will deliver all the required power for traction.  Figure 26 shows the power developed by each power source along with the FC control signal k FC f uzzy that represents the amount of total power that FC delivers at each time instant. From Figure 24, at P1, when k FC f uzzy = 0, FC is off and delivers no power. At P2 k FC f uzzy = 0.3 and this will make the FC delivers 30% of the required power for traction and the SC will ensure the remaining power. At P2 k FC f uzzy = 1 meaning that FC will deliver all the required power for traction. Figure 26 shows the power developed by each power source along with the FC control signal that represents the amount of total power that FC delivers at each time instant. From Figure 24, at P1, when = 0, FC is off and delivers no power. At P2 = 0.3 and this will make the FC delivers 30% of the required power for traction and the SC will ensure the remaining power. At P2 = 1 meaning that FC will deliver all the required power for traction.

Conclusions
In this work, hybrid electric vehicle controllers were simultaneously tuned using a genetic algorithm with minimization of PMSM torque and speed ripples as an optimization goal. The simulation results show a significant reduction in PMSM torque and speed ripples. This reflects the effectiveness of the cost function formulated in this paper, mainly to tradeoff between torque and speed ripple minimization. The developed fuzzy energy management strategy proposed in this paper operates the battery and the supercapacitor at safe operating points and enables different driving modes. The fuzzy power management used here equilibrated the battery and the supercapacitor's use and prevented exhausting one power source over another. Both numerical results obtained using MATLAB/Simulink and the real time simulation results obtained using the RT LAB simulator are satisfactory. As further work, we intend to take advantage of the superiority

Conclusions
In this work, hybrid electric vehicle controllers were simultaneously tuned using a genetic algorithm with minimization of PMSM torque and speed ripples as an optimization goal. The simulation results show a significant reduction in PMSM torque and speed ripples. This reflects the effectiveness of the cost function formulated in this paper, mainly to tradeoff between torque and speed ripple minimization. The developed fuzzy energy management strategy proposed in this paper operates the battery and the supercapacitor at safe operating points and enables different driving modes. The fuzzy power management used here equilibrated the battery and the supercapacitor's use and prevented exhausting one power source over another. Both numerical results obtained using MATLAB/Simulink and the real time simulation results obtained using the RT LAB simulator are satisfactory. As further work, we intend to take advantage of the superiority of metaheuristic algorithms in order to ensure optimal use of vehicle power sources, which will result in increasing their lifespan and vehicle autonomy.

Conflicts of Interest:
The authors declare no conflict of interest.