Design of an Aftermarket Hybridization Kit: Reducing Costs and Emissions Considering a Local Driving Cycle

: For decades, drivers and ﬂeet managers have been impacted by the instability of fuel prices, the need to save resources and the duty to meet and attain environmental regulations and certiﬁcations. Aiming to increase performance and efﬁciency and reduce emissions and mileage costs, plug-in electric vehicles (PHEVs) have been pointed out as a viable option, but there are gaps related to tools that could improve the numerous existing conventional vehicles. This study presents the design of an aftermarket hybridization kit that converts a vehicle originally driven by a combustion engine into a PHEV. To achieve this goal, an optimization was conducted with the objective of decreasing the cost (regarding fuel consumption and battery charging) to perform a local driving cycle, while attenuating the tailpipe emissions and reducing the battery mass. The torque curves of the electric motors, the battery capacity, the parameters for a gear shifting strategy and the parameters for a power split control were the design variables in the optimization process. This study used the Campinas driving cycle, which was experimentally obtained in a real-world driving scenario. The use of a local driving cycle to tune the design variables of an aftermarket optimization kit is important to achieve a customized product according to the selling location. Among the optimum solutions, the best trade-off conﬁguration was able to decrease the mileage cost in 22.55%, and reduce the tailpipe emissions by 28.4% CO, 33.55% NOx and 19.11% HC, with the addition of a 137 kg battery.


Introduction
Local and global environmental issues have resulted in the establishment of regulations to reduce vehicle emissions worldwide [1,2]. On the other hand, the instability of fuel prices leads customers to react by reducing consumption [3] or by changing energy source [4] when it is possible. Plugin hybrid electric vehicles are driven by electric motors (EMs) and internal combustion engines (ICE) and can be connected to the grid [5]. They are one of the solutions for mitigating local emissions [6] and absorbing fuel and electric energy price instabilities.
The relevance of the plug-in electric vehicles (PHEVs) in the current engineering is reflected by the recent scientific activities. Du et al. [7] optimized a plug-in hybrid bus to avoid excessive changes of operation modes, reducing the number of engine startups, the clutch abrasion and the jerk. Golpîra Therefore, the aim of the present study is to conduct a multi-objective optimization of the aftermarket hybridization kit. The optimization problem is formulated to decrease the mileage cost (in terms of fuel consumption and battery charging), mitigating the tailpipe emissions and decreasing the battery mass. The selected design variables are the EMs' torque curves, the battery capacity, the parameters of the gear shifting strategy and the parameters of the power split control. Moreover, this study is based on the Campinas driving cycle [38], which was experimentally obtained in a real-world driving scenario. Local driving cycles allow a more realistic analysis of the possible gains of the aftermarket hybridization vehicle, focusing on the local costumer.

Simulation Model
In the literature, models of vehicle longitudinal dynamics and vehicle components are often implemented in MATLAB TM [37,[39][40][41]. All the equations presented here have been implemented by the authors by means of MATLAB/Simulink TM block diagrams. In this study, the vehicle longitudinal dynamics are modeled based on the equations presented by Gillespie [42]. However, some modifications were required to include the electric drive train in the model.
The vehicle required traction torque T req (Nm) is defined by Equation (1) as a function of the movement resistance forces caused by the gravitational acceleration g (m/s 2 ), vehicle mass M (kg) and required acceleration a req (m/s 2 ), the aerodynamic drag D A (N), the tires rolling resistance R x (N) and the climbing resistance that varies according to the road angle α (rad).
The tire dynamic radius r (m) can be defined as proposed by Genta and Morello [43] as 98% of the geometric radius r g (m), that varies according to the tire standard geometry. The vehicle speed effects are also considered by means of the kv factor [44].
The rolling resistance R x (Equation (3)) and the aerodynamic drag D A (Equation (4)) increase with the vehicle speed V (m/s). Moreover, D A varies as a function of the air density ρ (kg/m 3 ), the projection of the vehicle frontal area A (m 2 ) and the C D coefficient that represents the vehicle shape.
The required acceleration a req represents the driver behavior. In this paper, a local driving cycle of the city of Campinas, Brazil [38], is applied ( Figure 2). The target speed V c (m/s) is defined as the cycle speed at one time step ∆ t (s) ahead of the current simulation time. Therefore, the requested acceleration can be calculated by Equation (5).  As aforementioned, the main objective of the aftermarket hybridization kit is to add an extra propelling system in the conventional vehicle. Therefore, the power management control (PMC) defines the traction torque required by the engine T eng (Nm) and the electric motors T el (Nm). In the same way, as presented by Eckert et al. [34], the extra electric motors act only as an auxiliary propelling system. Thus, the throttle pedal signal is acquired and the PMC splits the requested traction torque T eng according to the P E factor (defined in the optimization problem). The engine torque T eng and torque of the electric motors T el are defined by Equations (6) and (7) respectively, according to the engine inertia I e (kgm 2 ); gearbox and differential gear ratios N t and N d ; inertia I t (kgm 2 ) and I d (kgm 2 ); the frontal and real wheels inertia I w f (kgm 2 ) and I wr (kgm 2 ); and the overall powertrain efficiency η td .
T el = P E T req + I wr a req (7)

ICE System Restrictions
Once the required traction torque parcel exerted by the engine is defined, it has to be compared to the available traction torque T av (Nm) that corresponds to the 100% throttle curve (Figure 3a). Moreover, the clutch transmissible torque limit T cl (Nm) (that acts during the gear shifting) also is considered in the definition of the available traction torque. The clutch transmissible torque is calculated according to its friction coefficient µ cl , the number of clutch faces N (-) and internal and external disks radii R i (m) and R o (m) (see Equation (8)). Moreover, the spring force F n (N) (Figure 3b) controls the torque transmission according to the clutch pedal displacement (CPD).
It can be seen in Figure 3b that the clutch maximum F n value is applied when it is closed (CPD = 0). In this condition, the total engine output torque is transmitted to the gearbox. On the other hand, during the gear shifting process, the clutch disconnects the ICE from the gearbox (CPD = 0 to CPD = 100 mm) in 0.3 s. After the clutch fully opened (CPD = 100 mm), the gear ratio is changed in 0.2 s; then, the ICE is reconnected to the gearbox (CPD = 100 to CPD = 0 mm) in 0.5 s [32]. Therefore, the gearbox input torque T gb (Nm) is defined according to the rules presented in Equation (9) and the available traction torque at the vehicle frontal wheels T a f (Nm) is then defined by Equation (10).

Electric Drivetrain Restrictions
As the frontal propelling system, the electric drivetrain is also limited by the maximum torque of the EMs curves T ev (Nm) according to their current speed. Therefore, if the required torque T el is higher than the sum of available torques T ev of both in-wheel EMs, there will be performance losses. The available traction torque of the rear-drive system T ar (Nm) is defined by the combination of Equations (11) and (12).
T ar = T rl − I wr a req (12)

Traction Restrictions
Once the available traction torques T a f and T ar are defined, it is possible to verify whether the available torque can be transmitted through the tire/ground contact. The maximum transmissible torque for the frontal T F(max) (Nm) and rear T R(max) (Nm) wheels are defined by Equations (13) and (14) respectively [46]. µ represents the tire-ground peak friction coefficient, and L (m), h (m), b (m) and c (m) correspond respectively to the vehicle's wheelbase; the height of its center of gravity; and the longitudinal distance between its front and rear axles.
Therefore, the effective traction torque of the frontal T F (Nm) and rear T R (Nm) propelling system are defined by Equations (15) and (16).

Acceleration Iterative Process
With both effective traction torques, it is possible to define the vehicle current acceleration a x (m/s 2 ) by means of Equation (17), which is integrated by ODE5 of the Simulink TM , to find out the current vehicle speed V and also the displacement.
Due to the many constraints presented, the vehicle acceleration a x may become lower than the requested a req . Therefore, it is necessary to perform an iterative process among Equations (1), (6), (7) and (9)-(17), considering a req = a x until the convergence of the a x value. Moreover, the tire slipping has to be corrected in this iterative process. The coefficient e (Figure 4) estimates the difference between the tire tangential speed and the vehicle displacement speed as a function of the traction torque T F/R . This parameter changes the engine and EMs speeds, and therefore the maximum available torques according to the drivetrain respective curves.

Fuel Consumption and Emissions
After the convergence, the engine effective torque T e (Nm) and speed ω e are defined by the Equations (18) and (19).
These values correspond to the engine operation point, and are input to the ADVISOR TM fuel converter block, which defines the ICE fuel consumption F C (l) and the tailpipe emissions HC (g/km) NOx (g/km) and CO (g/km) according to the ICE maps shown in Figure 5. ADVISOR TM [47] is a free vehicular simulation tool that considers the ICE transient regimes, simulating the combustion heat and its influence in the catalyst efficiency that significantly changes the generated tailpipe emissions [32].

Battery and Equivalent Fuel Consumption
As initially proposed by Eckert et al. [34], the hybridization kit contains a lead-acid battery pack with a nominal voltage of 48 V. This voltage value was selected because the hybridization kit must be adapted to conventional vehicles. Electrical systems above 60 V can represent a potential danger to the passenger [48], and therefore require a shock protection system [49] that would increase the cost of the kit. Furthermore, the lead-acid batteries were selected due to their maintenance-free, good recycling capacity [50] associated with easy manufacture and high-volume production [51] that decrease the final cost of the kit.
In the simulations, the lead-acid battery block from Simulink TM database [34] is used. The developed control stops the use of the auxiliary electric drive train system when the battery state of charge SoC is under 40% to avoid fatigue caused by deep discharge [34]. Aiming to define the best battery pack for the hybridization kit, its capacity B c (Ah) has to be defined and included as input to the Simulink TM battery block. Moreover, the battery mass M bat (kg) is also calculated by Equation (20) as a function of the lead-acid specific energy S E = 40 (Wh/kg) and added to the vehicle mass M. Thus, the resulting M bat is added to the vehicle overall mass M.
Once the driving cycle is simulated, the final battery SoC is used as an input to the charge simulation, which infers the cost of recharging the battery from the electric grid. As in Eckert et al. [34], the battery is recharged at constant current [52]-the maximum I (A) value that does not overcome the allowed battery voltage V bat (V). The power consumption P c (W) to return the battery SoC to 100% is defined by Equation (21), considering a charge efficiency η c that concerns the conversion from AC to DC and the heat losses.
Differences in the cost of the energy sources may impact the result of the optimization of vehicle solutions [53]. To allow a fair comparison between the conventional vehicles and the hybridized one, the cost of battery charge is converted in a equivalent fuel consumption F eq (l) by Equation (22) according to the charging time C t (s) and the ratio between the fuel F cost ($/l) and electric energy E cost ($/Ws) costs. Finally, the F eq is added to the engine fuel consumption F C , resulting in the total fuel consumption F T (l) used to compare the results.

Electric Motors
One of the goals of this study was to find out optimum in-wheel EMs to be applied in the aftermarket hybridization kit; a generic torque curve and efficiency map are presented in Figure 6 [34,54,55]. This approach allows for the definition of a theoretical torque curve, based on three main points. The first one is the EM maximum torque T max (Nm) and its respective speed ω Tc (rad/s) at the constant torque phase. The second point is defined at the constant power phase, according to Tong [56], which defines the EM best operating region between 0.1T max and 0.3T max . In this paper, we adopt the upper limit of 0.3T max , as shown in Equation (24), and the EM speed of this point is defined by Equation (25). Finally, the last torque curve point represents the speed at which the EM reaches null torque, and it is defined by linear progression of the previously defined points (T max , ω Tc ) and (T Pc , ω Pc ).  [34,54,55].
With the EM torque curve fully defined, the map that defines the EM efficiency η EM is interpolated based on the data presented in Figure 6. Moreover, the EM inertia is interpolated based on some values presented by Corrêa et al. [57]. Finally, the inverter efficiency η inv is obtained from Table 1 and the electric current I (A) results from the Equation (26). The EMs are able to regenerate some part of the kinetic energy during the braking. In this paper, the regenerative braking is limited to 10% of the EM maximum torque [54,59]. If the required braking torque exceeds the bounded, the remaining parcel is dissipated by the conventional brake of the vehicle [59].

Vehicle Parameters
The simulated vehicle is based on the 1.0L engine frontal traction Brazilian Chevrolet Celta TM .
The parameters used to model this conventional vehicle and the aftermarket hybridization kit in the simulations are presented in Table 2.

Gear Shifting Strategy
Finally, the simulation parameters related to the applied gear shifting strategy have to be defined; they significantly influence the energy management [60], the vehicle performance, the fuel consumption and the emissions [32,33]. The use of an adequate gear shifting strategy improves the gains reached by the vehicle hybridization, once the electric motors decrease the ICE torque demand, allowing for the anticipation of the upshifts, moving the engine to a lower fuel consumption/better efficiency operation point [34]. On the other hand, anticipating the upshifts increases the ICE warm-up period, leading to poor catalyst efficiency and higher tailpipe emissions.
Therefore, the applied shifting strategy has to be carefully defined. The hybridization kit was developed to be adapted to a conventional vehicle a manual gearbox; the most adequate gear will be shown to the driver by means of a user interface adapted to the vehicle. As in previous works [32,34,45], the gear shifting strategy is based on speed limits; in other words, when the vehicle reaches one of these limits, the gear will be changed to the next (upshift) or previous gear (downshift). In the simulation the shifting process is controlled by the vector V shift (Equation (27)) that contains the upshift speeds V u(i) (km/h) and the speed interval V D (km/h) used to define the downshift speeds V d(i) (km/h) (Equation (28)). This difference between V u(i) and V d(i) is necessary to prevent gear shifting instabilities, as highlighted by Xi et al. [61].
The gearbox transmission ratio N t is then defined according to the shifting rule present in Equation (29).

Optimization
As mentioned before, the hybridization kit is composed of two in-wheel electric motors assembled directly in the vehicle rear wheels and powered by a lead-acid battery. The main objective of this study is, focusing on a local drive cycle (Campinas, Brazil), to define the best configuration for the battery and electric motors and also the power split control and gear shifting strategy associated with each proposed configuration. To do that, a multi-objective optimization problem is formulated and solved by means of the interactive adaptive-weight genetic algorithm (i−AWGA) proposed by Gen, Cheng and Lin [62]. The i−AWGA technique performs a wide search for the best solution, and it is not limited to false minimums, as occurs in some other optimization techniques. Moreover, this method was applied in several previous works, regarding vehicle multi-objective optimization [32][33][34]59], reaching satisfactory results.

Problem Formulation
The first optimization criterion f 1 (X) is the minimization of the total fuel consumption F T (Equation (23)) that represents the ICE fuel consumption and the conversion of the battery charge cost in hypothetical fuel (Equation (22)).
The second optimization criterion f 2 (X) is the minimization of the engine carbon monoxide (CO (g/km)), nitrogen oxide (NOx (g/km)) and hydrocarbon (HC (g/km)) emissions. These emissions are combined in a single parameter named emission factor E F [32,33]. The E F value is defined by the adaptive-weight technique [62] according to Equation (30) as a function of the maximum and minimum values of each gas emission presented in the solution database.
Finally, the third optimization criterion f 3 (X) is the minimization of the battery mass M bat that represents the majority of the extra weight added by the hybridization kit to the vehicle and reduces the available space in the vehicle trunk where it is meant be assembled.
Therefore, the optimization criteria are: subjected to the constraints C presented in Equation (34).
Besides the EMs and battery parameters used to narrow the optimization process, some extra constraints were included to ensure the correct behavior of the vehicle. The first one is the ICE operation range, which cannot be below the idle speed (ω min ≈ 84 rad/s) and cannot overcome the maximum allowed speed (ω max ≈ 680 rad/s).
The final constraints are included to avoid any configuration that presents poor acceleration performance caused by the fuel consumption minimization, and to prevent configurations that increase the tailpipe emissions due to the ICE warm-up delay. To be considered viable, the PHEV configuration cannot generate more CO, NOx and HC emissions than the conventional vehicle. Moreover, the vehicle performance is evaluated as the correlation coefficient R calculated by the Equation (35) [30,63] in which the simulated speed profile V s is compared to the standard driving cycle V cs in discrete time steps of 0.1 s, and the analyzed PHEV must reach at least the same performance of the conventional vehicle R = 0.9997, using a standard gear shifting strategy proposed by its manufacturer [64].

Genetic Algorithm
The optimization starts with the definition of chromosomes X (Equation (36)) that constrain the design variables of a specific potential solution.
The initial solution database is composed of randomly defined chromosomes (considering the constrains C) that are simulated and included in the population if they reach the minimum required performance.

Selection, Crossover and Mutation
Once an initial solutions population is developed, the results are classified by the i−AWGA technique (Equation (37)), which compares the current result of each optimization criterion f k (X) with the maximum f max k and minimum f min k values of the analyzed criterion presented in the population. Moreover, a penalty value P p = 0 is added to the dominated solutions, aiming to decrease their probability of reproduction. On the other hand, the non-dominated solutions (Pareto frontier) receive P p = 1 value, which increases their fitness value Ft(X), improving their selection probability as a function of the sum of all Ft(X) values according to the population size PM as shown by Equation (38).
S P (X) = Ft(X) The crossover process starts with the selection of two members of the population, which are named Member 1 (M 1 ) and Member 2 (M 2 ). The values of the design variables of these selected chromosomes are then randomly combined, generating a new chromosome denominated X cr which is simulated. It is included in the population if it respects the constraints C and if the minimum performance criteria are reached.
To provide diversity among the values of the design variables present in the population, the mutation process randomly changes some of the chromosomes values of the selected members and of the crossover chromosome X cr , generating three new chromosomes X mt1 , X mt2 and X mtc that are evaluated by the constrains C. Table 3 shows the mutation rules according to the operator 0 ≤ Mut ≤ 1 defined by the Matlab TM function rand.

Population Control and Convergence Criterion
The initial population is composed of 100 members (randomly generated chromosomes), which are combined by the crossover and mutation operators. The addition of new solutions increases the population size until it reaches a limit value P lim = 200. Once the population limit is reached, the worst solutions (higher Pareto ranking) are eliminated from the population. If a population composed of exclusively non-dominated solutions (Pareto frontier) reaches P lim , the limit is increased P lim = P lim + 100 to avoid the elimination of the whole population.
The convergence is defined by the stagnation of the evolution process, which is characterized by the repetition of the Pareto frontier [33,59,65] for over 20 generations of 80 new simulated solutions (crossover and mutation). Figure 7 shows the optimization flowchart according to the applied method.

Results and Discussion
After the convergence of the algorithm, the Pareto frontiers of non-dominated solutions are defined as optimum configurations with different compromise among the optimization criteria. Among these solutions, some of them are selected to be analyzed: the minimum energy cost solution min f 1 and the minimum emissions min f 2 , the minimum battery min f 3 . Moreover, the min f 1 solution presents the best trade-off (max Ft value) among the analyzed criteria.
However, the min f 1 (best trade-off) and min f 2 present large batteries (up to ≈144 kg), representing the advantages of the electric powertrain system against the combustion one. On the other hand, the minimum battery min f 3 solution was able to save fuel, but the tailpipe emissions were close to the ones of the conventional vehicle. In addition to these solutions, the Pareto frontier presents several other optimum solutions, and two of them were selected based on the decrease of the battery mass (criterion f 3 ). These extra solutions were named To f f 1.5 (making Ft( f 3 ) 50% more valuable: 1.5 × Ft( f 3 )) and To f f 2 (making Ft( f 3 ) 100% more valuable: 2 × Ft( f 3 )), increasing the importance of the Ft value related to the battery mass minimization, selecting the lighter hybridization kit configurations. Table 4 shows the results of the selected configurations and their respective chromosomes. Moreover, Figure 8 shows the Pareto frontier and highlights the selected configurations.

Standard Vehicle
The standard solution to be compared to the optimized PHEV configurations is the conventional vehicle configuration propelled only by the ICE and using a standard gear shifting strategy proposed by the vehicle manufacturer [64]. Figure 9 shows the ICE running time for each operation point (torque/speed), where it is possible to observe that the ICE presents several points close to the idle/low-speed regime, which is a characteristic of the Campinas driving cycle (Figure 2). The standard gear shifting strategy is shown in Figure 9; it represents the higher speed section of the cycle. As it can be seen in Figure 8c and Table 4, the conventional vehicle fuel consumption and resulting emissions are fully dominated by the optimum PHEV configurations.

Minimum Battery Size Solution min f 3
The minimum battery size solution focuses on minimizing the EMs and battery, resulting in a more compact hybridization kit, with a 24.14 kg battery and low power EMs responsible only for 8.22% of the vehicle power demand. This configuration allows for the upshifting anticipation (Figure 10c) compared to the standard shifting procedure, which saves fuel by moving the ICE operation point to a low-speed/high-torque [32,33], as can be observed in Figure 10a. The EMs operation is presented in Figure 10b where it is possible to observe that in the majority of the time, the electric driving system is used only as a supplementary propelling system, which ensures the vehicle performance even with the upshifting anticipation, and decreases the ICE torque, saving fuel. This configuration was able to reduce the equivalent fuel consumption (F T ) by 3.8%, and present tailpipe emissions close to the standard vehicle results (−3.03% CO, −2.3% NOx and −1.57% HC), as is shown in Table 4.

Minimum Emissions Solution min f 2
The minimum emission PHEV configuration(min f 2 ) presents a high level of electrification of the vehicle propelling system, in which the EMs are responsible for 55.4% of the power demand. Therefore, the EMs present high torque, as shown in Figure 11b, in which it is possible to observe that these EMs operate at maximum torque for a large amount of time. Moreover, these high torque EMs also expressively decrease the ICE required torque (Figure 11a), associated with the upshifting anticipation (Figure 11c), keeping the ICE running under ≈ 2700 rpm. On the other hand, this configuration needs a large battery of ≈144 kg to ensure the availability of the electric propelling system during the 42 km of the analyzed driving cycle. Due to the mentioned minimization of the ICE using, it was possible to save 1.22 L of fuel, compared to the conventional vehicle, which represents 21.74% economy regarding the overall cost to perform the desired path. Due to the heavy use of the electric drivetrain system, this PHEV configuration results in −33.88% CO, −35.36% NOx and −18.19% HC.
Regarding the cost relation between fuel and electricity, the presented results have shown a favorable scenario for vehicle electrification. However, the use of large batteries may lead to other issues, such as a higher cost of the aftermarket hybridization kit, difficulty of assembling the battery in the vehicle and possible interference in vehicle handling [34].

Minimum Cost Solution min f 1 (Best Trade-Off)
Regarding the optimized scenario, the minimum cost solution min f 1 , which represents the minimum equivalent fuel consumption (F T ) among the Pareto frontier is rated as the best-compromised solution (higher fitness value Ft). This PHEV configuration presents EMs that provide up to 105 Nm, which act similarly to the min f 2 EMs, keeping the ICE running at a low speed and torque regime, which together, lead to fuel economy. The electric propelling system is responsible to fulfill 66.20% of the power demand, which results in a high concentration of EM operation points close to the maximum torque curve, as shown in Figure 12b. This configuration was able to decrease the cost to perform the driving cycle in 22  Besides the satisfactory reached results regarding cost end emissions of the min f 1 configuration, this aftermarket hybridization kit needs a 137 kg battery, which implies in the same issues listed before regarding large batteries. Therefore, the configurations To f f 1.5 and To f f 2 analyzed below, are included to represent some alternative hybridization kits composed of lighter batteries and EMs, which can be more suitable in cases in which a lower powertrain electrification degree is desired. The first of the selected trade-off solutions was defined by assigning 50% more weight to the battery mass minimization at Equation (37) that determines the fitness Ft value, used to classify the population. This configuration is powered by a ≈ 90 kg battery and EMs able to perform up to 83 Nm, which operates the majority of the time at maximum torque regime, as shown in Figure 13b. Moreover, the power split control (P E ) was defined as 37.17% due to the lower EMs torque capacity compared to the min f 1 and min f 2 PHEV configurations. Similarly to the previously mentioned configurations, the To f f 1.5 keeps the ICE running at low torque and speed, which enables 950 mL fuel-saving and an equivalent fuel consumption 18.75% lower compared to the conventional vehicle. Furthermore, the To f f 1.5 solution also enhanced the ICE emissions, by decreasing 23.78% of the CO, 17.76% NOx and 11.12% HC.  The To f f 2 solution, in the same way as the minimum battery size PHEV configuration (min( f 3 )), represents a lower powertrain electrification degree. Even with the lower gains in terms of energy consumption, these configurations may become good alternatives to provide cheaper hybridization kits due to their smaller batteries and electric motors, in addition to their lower interference in the vehicle weight distribution and drivability.

Conclusions
In this paper, an aftermarket hybridization kit that converts a vehicle originally driven by a combustion engine into a PHEV was proposed. An optimization problem was formulated, aiming to decrease the overall cost (regarding fuel consumption and battery charging) to perform a local driving cycle, attenuating the tailpipe emissions and reducing the battery mass, which leads to a low price hybridization kit. Some design variables were taken into account; they could change the capacity and the weight of the batteries, the power and the torque of the EMs, the gear shifting strategy and the PMC strategy. Considering the needs of the customer, the optimization of the kit was suitable to a local driving cycle, reflecting the real service loads that it would face.
After the driving cycle was set, an optimization process was conducted by means of the i−AWGA, which was able to find out a Pareto frontier of optimum configurations of the kit. Some of these solutions could be selected, such as the ones that further reduced the travel cost or emissions, resulting in vehicles with a high degree of electrification. These solutions were suitable for customers who were prone to adding more weight to their vehicles and buying a kit with a higher purchase price (due to the cost of large battery packs).
On the other hand, an optimum configuration with the minor battery weight could be found, which would result in a minor purchase price and would be less effective at reducing travel costs and emissions. If battery size was used as a rule, it would be possible to find other optimum solutions in the Pareto frontier to significantly reduce the wanted index but with a non-extravagant additional weight of the kit.
Furthermore, the optimization algorithm led to another aftermarket hybridization kit configuration which presented good trade-off among the optimization criteria. Each of these optimum configurations presented improvements in the overall energy cost (fuel and electricity) and engine emissions according to the applied electrification degree resulting from the assembled electric motors and battery.
Finally, the presented results have shown that it is possible to obtain expressive enhancements by converting a conventional vehicle propelled only by a combustion engine into a PREV. However, this study is only theoretical, and some other steps are necessary to ensure the feasibility of the proposed concept.
Our future work will guide the decisions among all the Pareto solutions. A study regarding the aftermarket hybridization kit cost, the purchase price and the consumer acceptation will contribute to define the kit setup for a certain local. Regarding the emissions, the electric energy mix and the kit material, manufacture and disposal must be taken into account in future analyses to ensure environmental gains. Moreover, the effectiveness of the hybridization kits should be evaluated under a set of different driving scenarios, to ensure the robustness of the designed devices. Finally, the real-time control for vehicle power management control also needs to be developed, aiming to make the hybridization kit feasible.