A Modified ABC-SQP-Based Combined Approach for the Optimization of a Parallel Hybrid Electric Vehicle

In this paper, an improved fuel consumption and emissions control strategy based on a mathematical and heuristic approach is presented to optimize Parallel Hybrid Electric Vehicles (HEVs). The well-known Sequential Quadratic Programming mathematical method (SQP-Hessian approach) presents some limitations to achieve fuel consumption and emissions control optimization, as it is not able to find the global minimum, and it generally shows efficient results in local exploitation searches. The usage of a combined Modified Artificial Bee Colony algorithm (MABC) with the SQP approach is proposed in this work to obtain better optimal solutions and overcome these limitations. The optimization is performed with boundary conditions, considering that the optimized vehicle performance has to satisfy Partnership for a New Generation of Vehicles (PNGV) constraints. The weighting factor of the vehicle’s performance parameters in the objective function is varied, and optimization is carried out for two different driving cycles, namely Federal Test Procedure (FTP) and Economic commission Europe—Extra Urban Driving Cycle (ECE-EUDC), using the MABC and MABC with SQP approaches. The MABC with SQP approach shows better performance in terms of fuel consumption and emissions than the pure heuristic approach for the considered vehicle with similar boundary conditions. Moreover, it does not present significant penalties for final battery charging and it offers an optimized size of the key vehicle’s components for different driving cycles.


Introduction
Nowadays, the modern ways are initialized to minimize the impact on natural resources and to increase the dependency over the non-renewable resources, due to the increasing fuel prices and environmental issues [1]. The improvement in the automobile science helps to create the modern ways for minimizing the usage of gasoline [2]. Many years of devoted research has culminated in the improvement of Hybrid Electric Vehicle (HEV) technology. The increase in vehicle numbers subsequently increases the pollution, which causes global warming [3]. In order to save the earth balance, it is mandatory that automotive industries have to produce low-emission and better energy-efficient vehicles. HEV performance with reduced emissions and improved fuel efficiency seems to be promising as one kind of an alternative technology [4]. In hybrid cars, the battery in the vehicle acts as an additional energy source to supplement the main engine, so that the vehicle operates at its maximum possible efficiency [5]. The battery in the vehicle acts as an additional electrical form of energy source to supplement the main engine, such that the vehicle operating efficiency is shifted towards its maximum possible efficiency and the vehicle operation range is extended. Also, the control strategy [6] for the vehicle has to be properly adapted. The emission reduction and increment in fuel economy are the most difficult goals in the automotive industry. The cooperative agreements, legislative efforts and society causes the major constraints in the automotive industry. The aforementioned requirements are solved by using the capacity of the HEV as well as these HEV are considered as alternative for the conventional vehicles. Hence, the minimum fuel optimal control issue is solved by using the novel control strategy in parallel HEV. The major contributions of this research are indicated as follows. This paper studies the optimization of parameters of an electric assist control strategy for a parallel HEV. From this study, a Sequential Quadratic Programming (SQP) mathematical method for a parallel HEV has been developed with a Modified Artificial Bee Colony algorithm (MABC) to minimize the instantaneous fuel consumption of the vehicle. The developed control method mainly concentrates on the fuel consumption minimization along with the improvement in power demand and drivability. Furthermore, the hybrid combination of MABC with SQP is used to optimize the size of the components. Results show that the MABC with SQP approach shows better performance in terms of fuel consumption and emission than the pure heuristic approach for the considered vehicle with similar boundary conditions. This research paper can be summarized as follows. It presents a literature review of recent techniques for optimization of HEV in Section 2, the proposed method with control parameters configurations is described in Section 3, simulation results are provided in Section 4, and the main conclusions of the research are drawn in Section 5.

Literature Review
In the literature, several associated works are available which are based on optimization of HEVs. Certain critical works are assessed here to justify the contribution of this work.
Enang et al. [7] presents comprehensive, quasi-static longitudinal modeling and parallel HEV validation. Further research was carried out over the New European Driving Cycle (NEDC) using the validated model for investigating the impact of an early gear upshift on cumulative fuel consumption. In this work, without any change in the final battery charging state, better real-time fuel saving is achieved. In spite of the fuel saving potential of this control strategy, the technique employed is unable to guarantee cost function (fuel consumption) optimality. It also cannot satisfy itself in real time, in the end.
Zhang et al. [1] have thoroughly analyzed the various HEV energy management strategies based on bibliometric and propose that driving cycle recognition would help with better energy management of HEVs. For most common types of energy management strategies, which are used in depth in HEVs, the statistical analysis is carried out, with essential characteristics containing pros and cons. Interconnections and scope for change between the different energy management strategies are also exposed from the control point of view. But, this study has not given any basic design as an illustration and relevance for an optimized vehicle design and output for different energy management strategies. Similarly, Shankar has analyzed the design of parallel hybrid electric vehicles for Indian road conditions [8].
Long, and Nhan [9] used the Pheromone-Based Bees Algorithm (PBA) to optimize the size of key component and controlling strategy of parallel HEV. The developed PBA modifies the parameters of control strategy and size of key component for minimizing the emission and Fuel Consumption (FC), while vehicle efficiency meets PNGV constraints. In this analysis, MATLAB software has been employed as a simulation tool, and driving cycles Federal Test Procedure (FTP) and European Driving Cycle (ECE-EUDC) are utilized to evaluate FC, emissions, and dynamic performances [10]. However, in this paper, the effect and change in weightage of objective function parameters and its effect on vehicle design and change in its performance are not discussed [11]. The key components of the vehicle size, which would be considered for optimization, would be internal combustion engine, electric motor, and battery. The key objective function of this research is to minimize the Fuel Consumption ( ), Carbon Monoxide ( ), Nitrous Oxide ( ), and Hydrocarbons ( ), so as to obtain maximum vehicle performance, subject to Partnership for a New Generation of Vehicles (PNGV) constraint limits. The objective function is as follows in Equation (1): Distance is >140 ft in 5 s Delta State of Charge (∆SOC) is ≤0.5% Also, subjected to the following operating conditions: a. Engine power is Peng > 0 (10) The key components of the vehicle size, which would be considered for optimization, would be internal combustion engine, electric motor, and battery. The key objective function of this research is to minimize the Fuel Consumption (FC), Carbon Monoxide (CO), Nitrous Oxide (NOX), and Hydrocarbons (HC), so as to obtain maximum vehicle performance, subject to Partnership for a New Generation of Vehicles (PNGV) constraint limits. The objective function is as follows in Equation (1): where 'Y' is the feasible solution search space. 'y' is a variable vector that contains a vehicle component size and parameters of control strategy, and 'd 1 ' to 'd 4 ' are termed as weighting factors of different parameters considered in the objective function, which can vary between 0 and 1, indicating degree of weightage of the parameter optimized and to investigate the effect of different parameters on the optimization results. Subjected to v(y) > 0 for y = 1 to 8, the PNGV constraint limits as defined by the US Consortium are as stated below in Equations (2)-(9): Acceleration time is ≤12 s for 0-60 mph (2) Acceleration time is ≤5.3 s for 40-60 mph Acceleration time is ≤23.4 s for 0-85 mph (4) Gradeability is ≥6.5% at 55 mph for 1200 s Maximum speed is ≥85.1 mph Maximum acceleration is >16.4 ft/s 2 Distance is >140 ft in 5 s Delta State of Charge (∆SOC) is ≤0.5% Also, subjected to the following operating conditions: a. Engine power is Peng > 0 (10) b. The power output of the engine is restricted to the max power rating of the engine Peng < Peng_max_power (11) c. For the charge-sustaining control system response, the State of Charge (SOC) of the battery should be always within: the defined minimum and maximum limits. In this work, an analysis has been performed on HEV with two different driving cycles, by varying the values of weighting factors in Equation (1), to analyze their effect on vehicle performances. The control strategy is considered as an important role on the performance and fuel economy of a vehicle. This research paper uses the Electric Assist Control Strategy (EACS). This is a type of charge-sustaining strategy, which means that the battery charge is always maintained between low and high State of Charge (SOC) levels [4]. The primary energy source is the ICE in this technique and the secondary source is electric motor. EACS uses the electric motor at low speeds, where either the ICE does not operate efficiently or where the power required goes beyond the maximum deliverable torque of the engine. If the State of Charge (SOC) of battery is low, the engine supports the battery charging through the electric motor.

Driving Cycles
The FTP and ECE-EUDC transient cycles are considered as the base driving cycles for the analysis. These vehicle driving cycles are considered due to the fact that the vehicle under study is a light-duty passenger car and they cover the majority part of the different driving conditions that a vehicle would face during its usage [21]. The corresponding driving cycle details are as shown in Figures 2 and 3.
Energies 2020, 13, x FOR PEER REVIEW 5 of 21 c. For the charge-sustaining control system response, the State of Charge (SOC) of the battery should be always within: the defined minimum and maximum limits. In this work, an analysis has been performed on HEV with two different driving cycles, by varying the values of weighting factors in Equation (1), to analyze their effect on vehicle performances. The control strategy is considered as an important role on the performance and fuel economy of a vehicle. This research paper uses the Electric Assist Control Strategy (EACS). This is a type of charge-sustaining strategy, which means that the battery charge is always maintained between low and high State of Charge (SOC) levels [4]. The primary energy source is the ICE in this technique and the secondary source is electric motor. EACS uses the electric motor at low speeds, where either the ICE does not operate efficiently or where the power required goes beyond the maximum deliverable torque of the engine. If the State of Charge (SOC) of battery is low, the engine supports the battery charging through the electric motor.

Driving Cycles
The FTP and ECE-EUDC transient cycles are considered as the base driving cycles for the analysis. These vehicle driving cycles are considered due to the fact that the vehicle under study is a light-duty passenger car and they cover the majority part of the different driving conditions that a vehicle would face during its usage [21]. The corresponding driving cycle details are as shown in    c. For the charge-sustaining control system response, the State of Charge (SOC) of the battery should be always within: the defined minimum and maximum limits. In this work, an analysis has been performed on HEV with two different driving cycles, by varying the values of weighting factors in Equation (1), to analyze their effect on vehicle performances. The control strategy is considered as an important role on the performance and fuel economy of a vehicle. This research paper uses the Electric Assist Control Strategy (EACS). This is a type of charge-sustaining strategy, which means that the battery charge is always maintained between low and high State of Charge (SOC) levels [4]. The primary energy source is the ICE in this technique and the secondary source is electric motor. EACS uses the electric motor at low speeds, where either the ICE does not operate efficiently or where the power required goes beyond the maximum deliverable torque of the engine. If the State of Charge (SOC) of battery is low, the engine supports the battery charging through the electric motor.

Driving Cycles
The FTP and ECE-EUDC transient cycles are considered as the base driving cycles for the analysis. These vehicle driving cycles are considered due to the fact that the vehicle under study is a light-duty passenger car and they cover the majority part of the different driving conditions that a vehicle would face during its usage [21]. The corresponding driving cycle details are as shown in     Schedule, and represents city driving conditions of a passenger car vehicle. It consists of a city route with 17 stops for a distance of 7.5 miles. The average and the maximum speeds are 19.6 miles per hour and 56.7 miles per hour. The idle time and the total cycle time are 259 and 1369 s, respectively.

Control Strategy
To obtain the minimum fuel consumption and emission targets during vehicle optimization [6], the control strategy is used to identify the ideal operating point of the vehicle's engine and motor. The flexibility of this strategy lies in allowing a vehicle to adjust its controls by considering the driving location or limits for local control of vehicle, real-time adjustment to driving cycles, and incorporating the temperature effects on fuel use, engine-out emissions, and catalyst behavior. In this analysis, the parallel EACS has been utilized [22]: the EACS uses the motor for additional power, when needed by the vehicle, and maintains the charge in the batteries. This strategy uses the electric motor in different ways: 1.
The motor is utilized in all driving torque, when the vehicle is operated in certain minimum speed.

2.
The motor is utilized for torque assist at the operating speed of engine, when the required torque is greater than max torque which is delivered by the engine.

3.
Through regenerative braking, the motor charges the batteries.

4.
The motor will be used to produce the required torque when the engine runs inefficiently for a required engine torque at the given speed. Also, the engine shuts off at this condition.

5.
When there is low-battery SOC condition, the motor will charge the battery through the excess torque provided by the engine.

Vehicle Configuration for Optimization
The vehicle considered for the optimization is a parallel hybrid electric vehicle, with a gasoline engine and battery as sources of power for vehicle driving [2]. The key parameters of the vehicle considered are as depicted in Table 1. The optimization is performed using the MATLAB simulation tool [23]. In this analysis, the key input parameters considered are fuel converter scaling factor for torque range, motor/controller torque scaling range, number of battery modules, and a set of control strategy parameters, as shown in Table 2 below.
The mentioned parameters are given as input to the MABC and combined MABC-SQP algorithm and the resultant optimized outputs in terms of fuel consumption and emissions are obtained. The approaches are discussed in detail below.

Details on SQP-Hessian Approach
A bisection-based Sequential Quadratic Programming (SQP) method is used to determine the required component sizes, fuel usage, and emission performance based on the performance criteria specified as a part of the analytical approach. For the given minimization problem, an approximation of the Lagrangian function Hessian using a quasi-Newton updating method is made at each big iteration. Given for a general minimization problem description in a general problem, the main idea is to formulate a QP sub-problem based on a quadratic approximation of the Lagrangian function.
where f (x) is the main objective function, g i (x) are the inequality constraints, λ i is Lagrange multipliers under the non-negativity constraint, and m is the total number of restrictions. Implementation of SQP consists of four main stages, as outlined below.
Updating the Hessian Matrix. The definite positive quasi-Newton approximation of the Lagrangian Hessian equation, H, is calculated in each major iteration, where λ i , i = 1, . . . , m, is an estimate of the Lagrange multipliers. where Quadratic Programming (QP) Solution. The solution for the QP is computed using 2 steps: 1.
Generation of feasible points.

2.
Detection of the feasible point in which the convergence takes place with number of iterative sequences.
A search path generated by a linear summation of any combination of the Zk columns is guaranteed to remain at the limits of the active constraints. The matrix Zk is generated from the last m−l columns where If Zk is found, a direction of new searchd k is sought that minimizes q(d), whered k is in the null space of the active constraints.
Initialization. The algorithm calls for a feasible starting point. If the current point of the SQP approach is not feasible, the Linear Programming (LP) problem can be solved as is shown below: min γ, such that: The notation A i indicates the ith row of the matrix 'A', and b is the coefficient vector. Line Search and Merit Function. The QP sub-problem solution generates a vector, dk, that is used to create a new iterate: xk The phase length parameter αk is set to cause an appropriate decrease in a merit function. Powell (1978) recommends setting the penalty parameter as: This allows for a positive involvement from constraints which are inactive in the QP solution, but which have been active recently.
where || || represents the Euclidean norm. It means greater contributions from constraints with lower gradients to the penalty parameter, which would be the case for active constraints at the point of solution.

Heuristic Approach
Numerous sources describe different evolutionary approaches to the problem. These are used due to the fact that in most cases, for the problem, it is very hard to get a solution, and therefore, in many cases, it is unsolvable. Some of the heuristic approaches available in the literature include Particle Swarm Optimization (PSO), Genetic Algorithm (GA) [5,24] and Genetic Programming (GP), and Artificial Bee Colony (ABC) [25].
PSO has a general disadvantage that it will fall into local optimal solution in high-dimension space and very low convergence rate during iteration or premature convergence [26]. With GA, choosing the right number of the population size and number of generations is also tough. Also, GA is a type of non-deterministic algorithm and in most of the cases, it provides a sub-optimal solution and may not be best suitable for multi-constrained optimization problems. With GP, the generation of mathematical formulations may be too difficult to interpret, as its combination of multiple elementary functions often results in extremely complex models. Another major issue is the selection of suitable parameter settings to control the algorithm run. Also, the time required for GP to converge to the global optimum is dependent on the parameters that govern the evolution process. One of the most popular heuristic algorithms used for multi-constrained optimization with a multi-input and multi-output system is Artificial Bee Colony optimization (ABC). Some of the advantages of ABC over other algorithms Energies 2020, 13, 4529 9 of 22 include simplicity, flexibility, robustness, ability to explore local solutions, ease of implementation, and ability to handle the objective cost. Hence, this algorithm is used for the problem optimization and the details of the ABC algorithm are given below.

Modified Artificial Bee Colony (MABC)-Based Optimisation
The heuristic approach used for the optimization is the Modified Artificial Bee Colony (MABC) algorithm. In the MABC algorithm, a new search equation described by the global best guided ABC inspired by PSO [27] has been utilized.
MABC overcomes the drawbacks of ABC and considers the global best information from the below equation. vi where vi,j is the new neighboring food source, yj is the jth element of the global best solution, Ψi,j is a uniform random number in the range of 0 to 1.5, Φi,j is a random number in the range [−1, 1], and j ∈ {1, 2, . . . , n} is a randomly chosen index. The major steps of the MABC algorithm are outlined as follows: Initialize all parameters, Repeat while Termination criterions are not met.
Step 1. New food sources are computed in the Employee bee phase.
Step 2. Location of the food sources is updated based on their amount of nectar in the Onlooker Bee phase.
Step 3. Search the new food sources in place of rejected food sources in the Scout bee phase.
Step 4. Remember the best food source identified so far.
The output of the algorithm would be the best solution identified so far. The Modified ABC [25] is good at exploration but has a disadvantage of exploitation problems, which causes the algorithm to get trapped to local optima and it may not achieve the global optimal solution [23].

Simulation and Analysis
Wipke et al. [17] developed a National Renewable Energy Laboratory (NREL) MATLAB-based vehicle simulator called the ADVISOR tool for setting vehicle parameters for acceleration and grade tests with the required conditions. This tool provides the results for those tests and in addition to this, vehicle Fuel Consumption (FC) and emission parameters are also derived: HC, CO, and NOX. Gowrishankar and Nirmal [20] show how to fix the zero delta SOC correction technique with the tolerance limit, which is a very important step in PNGV constraints. This work analyzes the ECE-EUDC and FTP driving cycles.
The analysis is done with four different cases for two driving cycles, ECE-EUDC and FTP. Case 1 denotes results for the MABC algorithm with the FTP driving cycle, while Case 2 denotes the results for the same FTP driving cycle for the MABC-SQP algorithm. Whereas, Case 3 evaluates the MABC algorithm for the ECE-EUDC driving cycle and Case 4 explains the process of the ECE-EUDC driving cycle for the MABC-SQP algorithm. More weightage (50%) is given to minimization of fuel utilization in Cases 1 and 2, and the remaining 50% weightage is distributed among the other three emission vehicle parameters. In Cases 3 and 4, maximum weightage (70%) is given for emission and the remaining 30% weightage is kept for the fuel usage. The performance of the algorithm for various requirements of the vehicle is analyzed under these various objective functions with different weightage. Tony and Wipke [3] show the influences of the driving cycle on the peak speed and distance, with relation to the fuel consumption and emission. The optimal control problem setup is solved over the ECE-EUDC and FTP driving cycles, and the results are compared below to benchmark the proposed MABC + SQP Energies 2020, 13, 4529 10 of 22 algorithm over the MABC algorithm in finding the optimal solution. The analysis is done in two steps, namely MABC only and the MABC with SQP-Hessian approach. Table 3 shows the corresponding initial input and optimized values obtained for the considered driving cycles for various cases for the FTP driving cycle for the MABC algorithm.  Figure 4 shows the battery SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 59%.

Modified Artificial Bee Colony Algorithm (MABC) FTP Driving Cycle
Energies 2020, 13, x FOR PEER REVIEW 10 of 21     Also, Figures 7 and 8 show energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 28,000 kJ, and energy utilized by the motor controller is higher, 1200 kJ, during regeneration mode, in order to  Figure 5 shows fuel converter efficiency and as it is seen from the figure, there is no big impact and the efficiency is maintained at around 27%. Because of the large variation of speeds ranging from idle to nearing 60 mph top speed in this driving cycle, the loss at lower speeds is compensated. Figure 6 shows the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation.
Energies 2020, 13, x FOR PEER REVIEW 11 of 21 recharge the battery and keep it within the intended SOC level. Additionally, it can be seen that although the efficiency of the fuel converter is around 27% only, which is similar to the previous case, the energy loss has been reduced after optimization. Also, energy usage by braking components is also near to 1400 kJ in this mode, to support higher energy regeneration.   Also, Figures 7 and 8 show energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 28,000 kJ, and energy utilized by the motor controller is higher, 1200 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. Additionally, it can be seen that although the efficiency of the fuel converter is around 27% only, which is similar to the previous case, the energy loss has been reduced after optimization. Also, energy usage by braking components is also near to 1400 kJ in this mode, to support higher energy regeneration.

Modified Artificial Bee Colony + Sequential Programming Algorithm (MABC + SQP) FTP Driving Cycle
In order to obtain a balance between the exploration and exploitation issues and to overcome the same, a combined approach of using MABC as a global optimizer for exploration and SQP as a deterministic local search for exploitation has been used. The results show that with reduced emissions and still with an improved fuel economy compared to the initial design, as shown in the results for most of the cases, the proposed MABC with SQP-Hessian approach shows its superiority in finding an optimal solution compared to the other individual heuristic approach. Table 4 shows the corresponding initial input and optimized values obtained for the considered driving cycles for various cases of the MABC + SQP algorithm.

Modified Artificial Bee Colony + Sequential Programming Algorithm (MABC + SQP) FTP Driving Cycle
In order to obtain a balance between the exploration and exploitation issues and to overcome the same, a combined approach of using MABC as a global optimizer for exploration and SQP as a deterministic local search for exploitation has been used. The results show that with reduced emissions and still with an improved fuel economy compared to the initial design, as shown in the results for most of the cases, the proposed MABC with SQP-Hessian approach shows its superiority in finding an optimal solution compared to the other individual heuristic approach. Table 4 shows the corresponding initial input and optimized values obtained for the considered driving cycles for various cases of the MABC + SQP algorithm.

Modified Artificial Bee Colony + Sequential Programming Algorithm (MABC + SQP) FTP Driving Cycle
In order to obtain a balance between the exploration and exploitation issues and to overcome the same, a combined approach of using MABC as a global optimizer for exploration and SQP as a deterministic local search for exploitation has been used. The results show that with reduced emissions and still with an improved fuel economy compared to the initial design, as shown in the results for most of the cases, the proposed MABC with SQP-Hessian approach shows its superiority in finding an optimal solution compared to the other individual heuristic approach. Table 4 shows the corresponding initial input and optimized values obtained for the considered driving cycles for various cases of the MABC + SQP algorithm.
It is evident from the above Tables 3 and 4 that the solutions obtained for different cases are according to the weightage of FC for different cases. As for Cases 1 and 2, they have been set with a higher weightage value for FC, and the improvement in FC in mpg is also high when compared to Cases 3 and 4. Figures 9 and 10 show the vehicle performance change after optimization of its design with the MABC approach for the FTP driving cycle, Case 1. Figure 9 shows the battery SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 57%. Figure 10 shows fuel converter efficiency and as it is seen from the figure, there is no big impact and the efficiency is maintained around 26%. Because of the large variation of speeds ranging from idle to nearing 60 mph top speed in this driving cycle, the loss at lower speeds is compensated. Figure 11 depicts the electric motor torque-speed characteristics and also, the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation. It is evident from the above Tables 3 and 4 that the solutions obtained for different cases are according to the weightage of FC for different cases. As for Cases 1 and 2, they have been set with a higher weightage value for FC, and the improvement in FC in mpg is also high when compared to Cases 3 and 4. Figures 9 and 10 show the vehicle performance change after optimization of its design with the MABC approach for the FTP driving cycle, Case 1. Figure 9 shows the battery SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 57%.   Figure 10 shows fuel converter efficiency and as it is seen from the figure, there is no big impact and the efficiency is maintained around 26%. Because of the large variation of speeds ranging from idle to nearing 60 mph top speed in this driving cycle, the loss at lower speeds is compensated. Figure 11 depicts the electric motor torque-speed characteristics and also, the limits of maximum motoring about 27,000 kJ, and energy utilized by the motor controller is higher, 1300 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. Moreover, it can be seen that although the efficiency of the fuel converter is around 26% only, which is similar to the previous case, the energy loss has been reduced after optimization. Also, energy usage by braking components is also near to 1400 kJ in this mode, to support higher energy regeneration.   Additionally, Figures 12 and 13 show the energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 27,000 kJ, and energy utilized by the motor controller is higher, 1300 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. Moreover, it can be seen that although the efficiency of the fuel converter is around 26% only, which is similar to the previous case, the energy loss has been reduced after optimization. Also, energy usage by braking components is also near to 1400 kJ in this mode, to support higher energy regeneration. Figure 10 shows fuel converter efficiency and as it is seen from the figure, there is no big impact and the efficiency is maintained around 26%. Because of the large variation of speeds ranging from idle to nearing 60 mph top speed in this driving cycle, the loss at lower speeds is compensated. Figure 11 depicts the electric motor torque-speed characteristics and also, the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation.
Additionally, Figures 12 and 13 show the energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 27,000 kJ, and energy utilized by the motor controller is higher, 1300 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. Moreover, it can be seen that although the efficiency of the fuel converter is around 26% only, which is similar to the previous case, the energy loss has been reduced after optimization. Also, energy usage by braking components is also near to 1400 kJ in this mode, to support higher energy regeneration.    Table 5 shows the ECE-EUDC driving cycle for the MBAC algorithm.   Table 5 shows the ECE-EUDC driving cycle for the MBAC algorithm. Figures 14 and 15 show the vehicle performance change after optimization of its design with the MABC + SQP approach for the ECE-EUDC driving cycle, Case 1. Figure 14 shows the battery Energies 2020, 13, 4529 15 of 22 SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 39%.  Figures 14 and 15 show the vehicle performance change after optimization of its design with the MABC + SQP approach for the ECE-EUDC driving cycle, Case 1. Figure 14 shows the battery SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 39%.    Figure 15 shows fuel converter efficiency and as it is seen from the figure, there is no big impact and the efficiency is maintained at around 24%. Because of the large variation of speeds ranging from idle to nearing 75 mph top speed in this driving cycle, the loss at lower speed is compensated. Figure 16 shows the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation. Also, Figures 17 and 18 show energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 20,000 kJ, and energy utilized by the motor controller is higher, 1000 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. It can be seen that although the efficiency of the fuel converter is around 24% only, which is similar to Energies 2020, 13, 4529 16 of 22 the previous case, the energy loss has been reduced after optimization [28]. Also, energy usage by braking components is also near to 500 kJ in this mode, to support higher energy regeneration.

Modified Artificial Bee Colony Algorithm (MABC) ECE EUDC Driving Cycle
Energies 2020, 13, x FOR PEER REVIEW 15 of 21 Figure 15 shows fuel converter efficiency and as it is seen from the figure, there is no big impact and the efficiency is maintained at around 24%. Because of the large variation of speeds ranging from idle to nearing 75 mph top speed in this driving cycle, the loss at lower speed is compensated. Figure 15. Fuel converter efficiency. Figure 16 shows the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation. Also, Figures 17 and 18 show energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 20,000 kJ, and energy utilized by the motor controller is higher, 1000 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. It can be seen that although the efficiency of the fuel converter is around 24% only, which is similar to the previous case, the energy loss has been reduced after optimization [28]. Also, energy usage by braking components is also near to 500 kJ in this mode, to support higher energy regeneration.     Figure 16 shows the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation. Also, Figures 17 and 18 show energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 20,000 kJ, and energy utilized by the motor controller is higher, 1000 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. It can be seen that although the efficiency of the fuel converter is around 24% only, which is similar to the previous case, the energy loss has been reduced after optimization [28]. Also, energy usage by braking components is also near to 500 kJ in this mode, to support higher energy regeneration.     Figure 16 shows the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation. Also, Figures 17 and 18 show energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 20,000 kJ, and energy utilized by the motor controller is higher, 1000 kJ, during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. It can be seen that although the efficiency of the fuel converter is around 24% only, which is similar to the previous case, the energy loss has been reduced after optimization [28]. Also, energy usage by braking components is also near to 500 kJ in this mode, to support higher energy regeneration.    Figures 19 and 20 show the vehicle performance change after optimization of its design with the MABC approach for the ECE-EUDC driving cycle, Case 1. Figure 19 shows the battery SOC variation  Figures 19 and 20 show the vehicle performance change after optimization of its design with the MABC approach for the ECE-EUDC driving cycle, Case 1. Figure 19 shows the battery SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 64%.  Figures 19 and 20 show the vehicle performance change after optimization of its design with the MABC approach for the ECE-EUDC driving cycle, Case 1. Figure 19 shows the battery SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 64%.   Figure 21 illustrates the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation.

Modified Artificial Bee Colony Algorithm (MABC + SQP) ECE EUDC Driving Cycle
Moreover, Figures 22 and 23 represent the energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 19,000 kJ, and energy utilized by the motor controller is close to 800 kJ during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. Also, it can be seen   Figures 19 and 20 show the vehicle performance change after optimization of its design with the MABC approach for the ECE-EUDC driving cycle, Case 1. Figure 19 shows the battery SOC variation with respect to the driving cycle, and it can be observed from the figure that variation between the initial and final value of battery SOC is less than 0.5% and SOC value lies close to 64%.   Figure 21 illustrates the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation.

Modified Artificial Bee Colony Algorithm (MABC + SQP) ECE EUDC Driving Cycle
Moreover, Figures 22 and 23 represent the energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 19,000 kJ, and energy utilized by the motor controller is close to 800 kJ during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. Also, it can be seen  Figure 21 illustrates the electric motor torque-speed characteristics and also the limits of maximum motoring and regeneration torque. At lower rpm, the torque points lie more on the positive region, supplementing the engine operation, and at rpm greater than 2000, torque points also shift to the negative region due to regeneration operation.
Energies 2020, 13, x FOR PEER REVIEW 17 of 21 that although the efficiency of the fuel converter is around 24% only, which is similar to the previous case, the energy loss has been reduced after optimization. Furthermore, energy usage by braking components is also near to 550 kJ in this mode, to support higher energy regeneration.  Moreover, Figures 22 and 23 represent the energy utilization of different key components of the vehicle during power mode and regeneration mode of the vehicle. It can be observed from these figures that the energy utilized by the fuel converter is higher during the power mode, which is about 19,000 kJ, and energy utilized by the motor controller is close to 800 kJ during regeneration mode, in order to recharge the battery and keep it within the intended SOC level. Also, it can be seen that although the efficiency of the fuel converter is around 24% only, which is similar to the previous case, the energy loss has been reduced after optimization. Furthermore, energy usage by braking components is also near to 550 kJ in this mode, to support higher energy regeneration.

Inference
The result for BABC has been referred from this work [20] as it was also using similar vehicle configuration. To illustrate the improvement in the vehicle performance after optimization, Case 1 results of all the three methods discussed above are tabulated below for comparison. Tables 6-8 show the comparative analysis for the two driving cycles and its units for the constraints are already mentioned in the respective tables. As seen from Tables 6, 7 and 8, there is a significant increase in the fuel consumption in terms of miles per gallon. This shows that the combined MABC + SQP method can obtain superior results in terms of vehicle performance compared to other methods. Similarly, for the other cases, based on the weightage of the variables, either the fuel efficiency would be increased, or the emissions would be reduced, as discussed in the previous section.

Inference
The result for BABC has been referred from this work [20] as it was also using similar vehicle configuration. To illustrate the improvement in the vehicle performance after optimization, Case 1 results of all the three methods discussed above are tabulated below for comparison. Tables 6-8 show the comparative analysis for the two driving cycles and its units for the constraints are already mentioned in the respective tables. As seen from Tables 6, 7 and 8, there is a significant increase in the fuel consumption in terms of miles per gallon. This shows that the combined MABC + SQP method can obtain superior results in terms of vehicle performance compared to other methods. Similarly, for the other cases, based on the weightage of the variables, either the fuel efficiency would be increased, or the emissions would be reduced, as discussed in the previous section.

Inference
The result for BABC has been referred from this work [20] as it was also using similar vehicle configuration. To illustrate the improvement in the vehicle performance after optimization, Case 1 results of all the three methods discussed above are tabulated below for comparison. Tables 6-8 show the comparative analysis for the two driving cycles and its units for the constraints are already mentioned in the respective tables. As seen from Tables 6-8, there is a significant increase in the fuel consumption in terms of miles per gallon. This shows that the combined MABC + SQP method can obtain superior results in terms of vehicle performance compared to other methods. Similarly, for the other cases, based on the weightage of the variables, either the fuel efficiency would be increased, or the emissions would be reduced, as discussed in the previous section.

Conclusions
The design of a parallel HEV involves a number of variables that must be optimized for better fuel economy and vehicle performance. Driving cycle pattern also plays a key role in assessing the fuel economy and emissions. The objective of this study was to minimize the overall fuel consumption of a parallel HEV on FTP and ECE-EUDC driving cycles to improve the vehicle performance. The weighting factors also have the influence on determining the parameters for the vehicle performance. In this research, a hybrid MABC with SQP approach was proposed to obtain better optimal solutions. This method was analyzed with various boundary conditions and PNGV constraints, which occur in every driving cycle and test case. From the test cases, the performance of HEV and its PNGV constraints were analyzed. From the results, it was clearly indicated that the proposed approach achieved an increased fuel economy of 12% for FTP and 20% for ECE-EUDC driving cycles. In most of the cases, the size of the Energy Storage System (ESS) components were minimized, which will decrease the overall cost of the HEV in parallel. Furthermore, the proposed MABC-SQP approach decreased one of the engine emission parameters, HC, up to 19% for the ECE-EUDC driving cycle. While considering the FTP driving cycle, the engine emission parameter CO was decreased up to 19%. While considering the SOC portion, the change was within 0.5% tolerance on every solution with reduced battery modules, which satisfies the required vehicle performance. While comparing the results of the individual heuristic approach [20], the proposed MABC + SQP approach achieved better results, as indicated above. Also, in comparison to a MABC algorithm, the hybrid MABC + SQP solution was found to be outstanding in that it achieved remarkable real-time fuel savings and reduced emissions, with no significant penalties for final battery charging and reduced size of key vehicle components for different driving cycles.
Author Contributions: Methodology, formal analysis, and investigation, S.N.S.; software, resources, and data curation, S.K.; data collection and validation, S.P.; writing-original draft preparation and review and editing, R.P.d.P.; supervision and project administration, B.D.P. All authors have read and agreed to the published version of the manuscript.