To reduce apace extraction of natural resources, to plummet the toxic emissions, and to increase the fuel economy for road transportation, hybrid vehicles are found to be promising. Hybrid vehicles use batteries and engine to propel the vehicle which minimizes dependence on liquid fuels. Battery is an important component of hybrid vehicles and is mainly characterized by its state of charge level. Here a modified state of charge estimation algorithm is applied, which includes not only coulomb counting but also open circuit voltage, weighting factor, and correction factor to track the run time state of charge efficiently. Further, presence of battery and engine together needs a prevailing power split scheme for their efficient utilization. In this paper, a fuel efficient energy management strategy for power-split hybrid electric vehicle using modified state of charge estimation method is developed. Here, the optimal values of various governing parameters are firstly computed with genetic algorithm and then fed to Pontryagin’s minimum principle to decide the threshold power at which engine is turned on. This process makes the proposed method robust and provides better chance to improve the fuel efficiency. Engine efficient operating region is identified to operate vehicle in efficient regions and reduce fuel consumption.
1. Introduction
The invention of the automobile is one of the most groundbreaking advancements in technology. Today, it is impossible to imagine the world without it anymore. The automobile industry contributes significantly to the growth of the world’s economy and affects each level of population. The present transportation structure heavily relies on internal combustion engine (ICE) based transportation, which uses fossil fuels as a source of energy. But due to toxic emissions of carbon dioxide (CO2), carbon monoxide (CO), nitrogen oxides (NOx), and unburned hydrocarbons (HCs) in large amount, they have caused environmental pollution and global warming as well. Exponential rise in population and personal transportation resulted in multifold increase in automobiles around the globe. It has been causing severe environmental problems and a threat to human life. Air pollution is a major environmental jeopardy to health due to emissions of CO2 [1]. 23 percent of total CO2 emissions in the world are caused by the transport sector [2], of which roughly 73 percent was generated by road transport [3]. In future oil production will fall, but its consumption will continue to rise so transport sector should eradicate dependence on oil by adapting the new transportation mediums like electric or hybrid vehicles which are green and sustainable. Hybrid vehicles are clean, efficient, and environment friendly transportation means. Hybrid electric vehicles (HEVs) use battery to store the electrical energy for propelling the vehicle with good fuel economy and less toxic emissions [4].
The presence of two power sources focuses on the need of designing an energy management strategy to split power between them to minimize the fuel consumption and maximize the power utilization. Complex structure of HEVs makes it challenging to design the control strategies. The preliminary objective of the control strategy is to satisfy the driver’s power demand with minimum fuel consumption. Moreover, fuel economy and emissions minimization are conflicting objectives; a smart control strategy should satisfy a trade-off between them.
In optimization based control strategies, the goal of a controller is to minimize the cost function. The cost function for an HEV may include the emission, fuel consumption, and torque, depending on the application. Global optimum solutions can be obtained by performing optimization over a fixed driving cycle. Due to causal nature of global optimization techniques, they are not suitable for real-time analysis. Therefore, global criterion is reduced to an instantaneous optimization, by introducing a cost function that depends only on the present state of the system parameters. Global optimization techniques do not consider variations of battery state of charge (SOC). Hence, a real-time optimization is performed for power split while maintaining the battery charge. The concept of real-time control strategy for efficiency and emission optimization of a parallel HEV is proposed in [5]. It considers all engine-motor torque pairs which forecast the energy consumption and emissions for every given point. An instantaneous fuel efficiency optimization strategy was developed for parallel hybrid vehicle with the charge sustaining mode in [6]. Also, to implement the global constraint, the authors developed a nonlinear penalty function in terms of battery SOC deviation from its desired value. Literature speaks that real-time optimization techniques like ECMS [6, 7], model predictive control (MPC) [8–10], Neural Network (NN) [11, 12], particle swarm optimization (PSO) [13–15], and Pontryagin’s minimum principle (PMP) [16, 17] are used extensively. Table 1 compares different real-time strategies with its pros and cons. In the presence of a priori knowledge, PMP with structural simplicity and limited computation time gives local solution to the optimization problem which in case of particular assumption can provide the global optima [17].
Comparison chart for real time optimization algorithms.
Methods
Structural complexity
Computation time
Type of solution
Requirement of a priori knowledge
Particle swarm optimization
No
More
Global
No
Energy consumption minimization strategy
Yes
Small
Local
No
Pontryagin’s minimum principle
No
Small
Local
Yes
Model predictive
No
Small
Global
No
Neural network
Yes
Small
Global
Yes
Stockar et al. used PMP to build an optimal supervisory controller by reducing a global optimization problem into local. It reduces computational requirement and gives the freedom to solve the problem in the continuous time domain [18]. Stockar et al. proposed a model based control strategy to minimize the CO2 emission. A supervisory energy management strategy is implemented as a global optimization problem and then converted into local and, using PMP, optimal energy utilization for PHEVs is obtained. A real-time optimal control can be obtained using PMP as it uses instantaneous minimization of the Hamiltonian function [19]. Kim et al. state that solution based on PMP can be global optimal under some certain assumptions [20]. Kim et al. applied PMP based control strategies to the PHEVs and found that it gives a number of alternative solutions [21]. PMP provides a near-optimal solution for optimal power management of HEVs if future driving conditions are known. It is suggested to find the proper costate, to keep SOC at a desired and predefined level.
In this paper PMP is applied to solve the performance index of the HEV in terms of fuel consumption. The requested power threshold is analyzed using PMP above which engine should be on, to work in its optimal efficient range, and generator can run to charge the battery to maintain the appropriate SOC level in the battery.
In HEVs, speed and power required by vehicle, SOC level, and engine off time play an important role to design energy management controller. PMP computes threshold power level, but initially designer needs to define the other parameters like optimal speed and torque ranges of engine, motor, and generator, highest and lowest level of SOC, target SOC, speed threshold, and engine off threshold to minimize the fuel consumption using vehicle controller. Optimal speed and torque ranges of engine, motor, and generator are defined by experiments done by the National Renewable Energy Laboratory (NREL) and provided in the package of Advanced Vehicle Simulator (ADVISOR) [22]. But speed threshold, engine off time threshold, and engine on SOC level are determined here using genetic algorithm (GA) which in turn, along with power threshold, decides engine on threshold.
GA is a heuristic search algorithm to solve optimization and search problems. This is a branch of artificial intelligence inspired by Darwin’s theory of evolution. GA is a robust and feasible approach with a wide range of search space and rapidly optimizes the parameters using simple operations. They are proven to be effective to solve complex engineering optimization problems, characterized by nonlinear, multimodal, nonconvex objective functions. GA is efficient at searching the global optima, without getting stuck in local optima. Unlike the conventional gradient based method, GA does not require any strong assumption or additional information about objective parameters. GA can also explore the solution space very efficiently. Piccolo et al. utilize GA for energy management of an on road vehicle and minimize the cost function containing fuel consumption and emission [23]. Wang and Yang implemented a robust, easy, and real-time implementable FL based energy management strategy and used GA to tune and optimize the same [24]. To optimize the fuel consumption and emissions in a series HEV, GA based control strategy has been used by [25]. It is a flexible and global optimal multiobjective control strategy which is found to be better than thermostatic and divide rectangle (DIRECT) algorithm. Wimalendra et al. applied GA in parallel HEV to find the optimal power split for improved vehicle performance and also promises to give maximum fuel economy for known driving cycle for a parallel HEV using GA [26].
This paper aims to develop an optimal controller based on PMP and GA to get the optimal power split between engine and battery to fulfill the driver’s speed and torque demand while compelling the engine to work in its efficient region. PMP is a powerful tool in optimal control theory which provides the set of necessary conditions to get the global solution of a constrained control problem.
The paper is divided into different sections. Section 2 explains the vehicle dynamics with vehicle model and planetary gear set operation. Section 3 explains power train control methodology involving engine speed control and traction torque control schemes. Section 4 elaborates problem statement, engine operating range description, and proposed control strategy. Section 5 discusses simulation and result discussion and Section 6 concludes the paper.
2. Vehicle Dynamics
The vehicle movement behavior depends upon different forces (aerodynamic drag, rolling resistance, and grading resistance) along its moving direction. Aerodynamic drag force is encountered by air in the direction of vehicle movement at a particular speed. Rolling resistance is a horizontal force which acts on the wheel center in the opposite movement direction of the wheel. Grade force acts on the vehicle either in opposite or in the same direction when a vehicle goes up or down over a slope: (1)Fr=12ρAfCDV-VW2+Pfr+Mgsinα,where Af is vehicle frontal area, CD is aerodynamic drag that characterizes the shape of the vehicle body, ρ is air density, V is vehicle speed and VW is component of wind speed with vehicle moving direction, P is force acting on the center of a standstill tire, fr is rolling resistance, and α is road angle.
Figure 1 shows the main components of the HEVs, that is, motor, generator, battery, and engine [4]. Presence of the engine and battery together in vehicle demands for coupler to add their speeds. In Toyota hybrid system (THS), planetary gear system (PGS) is used as a speed coupler. PGS contains carrier, sun, ring gear, and several pinion gears as shown in Figure 2. The ring gear is attached to the motor and final drive, engine to the carrier, and generator to the sun. Governing equations between different gear speeds and radii are given as follows:(2)ωr∗rr=-ωs∗rs+ωcrs+rr,where ωr, ωs, and ωc are ring, sun, and carrier angular speeds, respectively, and rr, rs are ring and sun radii, respectively. Neglecting energy losses in steady state operation and torques acting on sun, ring, and carrier have the relationship as follows:(3)Tc=-kysTs=-kyrTr,Tc, Ts, and Tr are the torques acting on carrier, sun, and ring gear, kyr=(1+ig)/ig and kys=(1+ig), and ig is gear ratio. While moving, engine speed ωe, motor speed ωm, and generator speed ωg are related as follows:(4)NrNs+Nr∗ωm+NsNs+Nr∗ωg=ωe,where Nr and Ns are tooth number in ring and sun gear, respectively, in Toyota Prius. As Nr=78 and Ns=30, (2) becomes(5)7.2222∗ωm+0.2778∗ωg=ωe.This equation describes that ωm is directly proportional to the linear speed of the vehicle with a quantitative change due to tire radius and final drive ratio.
Power split hybrid architecture.
Operation of a planetary gear.
2.1. Battery Modeling and SOC Estimation
In general, classical SOC estimation is performed using ampere hour counting method (in ADVISOR also), but open circuit voltage (OCV) also plays an important role in determining the SOC. Tang et al. [27] and Verbrugge and Tate [28] identified the contribution of both coulomb counting method (SOCi) and open circuit voltage method (SOCV) together to estimate the accurate SOC. References [29–31] have also identified the importance of SOCV and SOCi in calculating the run-time SOC. The SOC estimation formula proposed by the authors is given as follows:(6)SOC=wSOCV+1-wSOCi-η,where η is correction factor (CF). CF varies with the changing SOC, load L, and temperature T (i.e., CF=fSOC0,L,T) and can be formulated as in (7): (7)η298=1-SOC0100,L=0 at T=298K,ηnewSOC,T=η298+SOC/1000T-298±ε,L>0 at any T.Estimating SOC by (7) and (8) will promise a better fuel efficiency of HEV as accuracy of SOC estimation is improved.
Battery plays a vital role in HEVs. In most of the literature, energy management techniques for HEVs have used batteries with a single Rint component which consists of ohmic and polarization resistances. But, due to double-layer formation at the electrode/solution interface, capacitive effects arise [32]. This capacitance consists of purely electrical polarization capacitance and diffusion capacitance [33]. The transient response of the battery is highly influenced by double-layer and diffusion capacitance when the rates of reactions are high. This effect can be modeled using lumped capacitances in parallel with the resistances [34]. Inclusion of diffusion and double-layer resistances and capacitances (R and C components) will lead to the accurate SOC estimation. In this paper, to predict the run-time behavior of the battery, 1 RC and 2 RC models along with modified SOC estimation techniques are used to analyze the effect on fuel efficiency.
Rate of change of SOC depends on P_bat, open circuit voltage (OCV), and resistance R offered by the battery cells and capacity Qp shown in(8)SOC˙=OCV-OCV2-4∗R∗P_bat2∗R∗Qp.Required power of P_bat can be calculated as follows:(9)Pbat=ηc1kTmg1ωmg1+ηc2kTmg2ωmg2,where ηc1k and ηc2k are the efficiencies of mg1 and mg2, respectively, and are obtained from the efficiency map of mgs. Positive k represents motoring operation and negative k represents generating operations. Equations (6), (8), and (9) are applicable for different battery models proposed in the literature.
3. Powertrain Control Methodology
Power split HEVs have the potential to improve in fuel efficiency compared to series or parallel hybrids because engine speed and torque can be decoupled completely or partially from the driven wheels through speed and torque coupling. By applying suitable control strategies fuel efficiency can be improved provided it follows the control objectives, like (1) driver torque and speed demand are fulfilled, (2) engine operates in its best efficiency region, (3) target SOC level meets at the end of the trip, and (4) maximum braking energy is recuperated while braking or decelerating. While making the control strategies, different approaches can be followed as elaborated below.
3.1. Engine Speed Control Strategy
Vehicle speed ranges are divided into three regions, namely, (1) low, (2) medium, and (3) high vehicle speed, as shown in Figure 3. In low speed region motor fulfills the driver power demand and hence engine usage can be avoided which is inefficient also. Low vehicle speed, VL, threshold can be decided by the lowest engine speed allowed with zero motor/generator speed as follows:(10)VL=πkyrne_minrw30irwm/s,where ne_min is the minimum engine speed allowed, rw is the wheel radius, kyr=(1+ig)/ig, where ig is the gear ratio and is defined as rr/rs and irw is the gear ratio of the ring gear to drive train wheels. In this region motor/generator operates with a positive speed nm/g as follows:(11)nm/g=kysne_min-30irwVπkyrrw,V is the vehicle speed in m/s (V≤VL). From (3) torque produced by motor/generator, applied to the sun gear, has direction opposite to its speed. Therefore, motor/generator absorbs part of the engine power to charge the battery. Power on the motor/generator shaft Pm/g can be expressed as (12). Tm/g is torque produced by motor/generator.(12)Pm/g=2π60Tm/gnm/g=2π60Tene_min-irwkyrrwTeV.When the vehicle speed is higher than VL but lower than VH given by (13), motor/generator is deenergized and sun gear is locked to the stationary frame of the vehicle. Drive train operates in torque coupling mode. Engine speed is proportional to the vehicle speed. Consider (13)VH=πkyrne_maxrw30irwm/s,where ne_max is the maximum engine RPM allowed. In this medium speed region all the engine power is delivered to the wheels.
Various vehicle speed ranges.
When the vehicle speed is higher than the VH, for limiting the engine speed below the maximum engine allowed speed ne_max, motor/generator has to operate in the direction opposite to the engine speed. It can be expressed as follows:(14)nm/g=kysne_max-30kysirwVπkyrrw,where V≥VH. The motor generator is in motoring mode, and motoring power can be expressed as follows:(15)Pm/g=2π60Tm/gnm/g=irwkyrrwTeV-2π60irwkyrrwTene_max.
3.2. Traction Torque Control Strategy
In low vehicle speed region when sufficient SOC is available traction motor torque Tmt can be given as follows:(16)Tmt=602πPm/gntm=ne_minntm-irwkyrimwTe=-2πrw60imwne_minV-irwkyrimwTe,where imw is gear ratio from the traction motor to the driven wheels and ntm is traction motor speed. PGS, m/g, and traction motor together function as an EVT, because no energy goes into or out of the battery.
In case of medium vehicle speed range only the torque coupling mode is employed; that is, sun gear is locked to the vehicle stationary frame and engine speed is proportional to the vehicle speed. In high speed region engine speed is controlled by the engine max speed ne_max and the motor/generator works in motoring mode. If the commanded traction torque is higher than the torque that the engine can produce with its optimal throttle at the speed of ne_max and SOC of the battery is lower than SOCmin and the battery cannot be discharged any more to support motoring mode, the engine will be forced to operate at the higher speed (beyond the optimal range) to fulfill the driver power demand. In this case, engine alone mode can be used with torque coupling or engine can run at somewhat higher speed so that a motor/generator can work in generating mode to feed the traction motor to support engine by providing additional torque. For the latter case, ne can be calculated as in (17):(17)ne>30irwVπkyrrw.
If SOC is higher than the SOCmin, then the engine should be controlled at its ne_max with optimal throttle and traction motor provides additional torque to engine to support the driver torque demand.
If the commanded traction torque is smaller than the engine torque and SOC is lower than SOCmin, engine is operated according to (13) and traction motor works in generating mode. If SOC is in between range of SOCmin and SOCmax, traction motor may be de-energized and engine alone mode can be projected. If SOC is greater than the SOCmax, engine better shuts down and traction motor alone can propel the vehicle.
4. Proposed Energy Management Approach
In HEVs, presence of both motor and engine together makes it inevitable to decide engine/motor on/off condition to minimize the fuel consumption. To split the power optimally between two power sources, a cost function is derived. The cost function depends on various parameters like speed, power, SOC, and engine on/off time. The various steps involved in developing the strategy are given below.
4.1. Problem Statement
The proposed cost function involves rate of fuel consumption; that is, J=m˙ft, where m˙ft is total fuel consumption in a driving cycle. m˙f is the time rate of fuel consumption and is given by m˙f=((Pe∗ge)/(1000∗γf))(l/h) where Pe is engine power, ge is specific fuel consumption, and γf is mass density of fuel, kg/L. So total fuel consumption in a driving cycle is m˙ft=∑(Pege/1000γf)∗Δti. The cost function is minimized over ECE_EUDC driving cycle subject to the following constraints:(18)ωe,min≤ωe≤ωe,max,ωmg1,min≤ωmg1≤ωmg1,max,ωmg2,min≤ωmg2≤ωmg2,max,Te,min≤Te≤Te,max,Tmg1,min≤Tmg1≤Tmg1,max,Tmg2,min≤Tmg2≤Tmg2,max,SOCmin≤SOC≤SOCmax,where ωe,min, ωe,max, ωmg1,min, ωmg1,max, ωmg2,min, ωmg2,max, Te,min, Te,max, Tmg1,max, Tmg1,min, Tmg2,min, Tmg2,max, SOCmin, and SOCmax are the minimum and maximum values of speed and torque considered as constraints range of engine, mg1, mg2, and SOC, respectively.
Torques and speeds of mg1 and mg2 are functions of engine torque and speed, requested driving speed and torque, and gear ratios of the vehicle as follows:(19)Tmg1=-11+RTe,ωmg1=-Rζωreq+1+Rωe,Tmg2=-11+R-1+RTreqζ+RTe,ωmg2=ζωreq,where ωmg1, Tmg1, ωmg2, Tmg2, ωe, and Te are speeds and torques mg1, mg2, engine, respectively, and ωreq and Treq are the requested speed and torque. R and ζ are the gear ratio of PGS and the final drive ratio [35, 36]. As efficiency of an engine is a function of engine speed ωe and torque Te, fuel consumption will be m˙f=f(ωe,Te).
Power requested should always be delivered by either motor, engine, or generator; that is, for the successful trip completion Prequested=Pdelivered=Pengine+Pmotor+Pgenerator.
Speed, force, and torque requested by ECE_EUDC shown in Figure 4 are used to calculate power required at the wheel. Positive force/torque value shows that power is required to propel the vehicle and negative force/torque specifies that the energy will be released and regenerative braking will be applied to recuperate the released energy in the battery. Vehicle component’s and drive cycle specification are given in Table 2.
Vehicle components and drive cycle specifications.
Vehicle component specification (Toyota Prius)
Drive cycle specification (ECE_EUDC)
Components
Values
Entities
Values
Motor
31 kW
Maximum speed
74.56 mph
Engine
43 kW
Average speed
19.95 mph
Heating value of gasoline QHV
42600 J/g
Maximum acceleration
3.46 ft/s2
Generator
15 kW
Maximum deceleration
−4.56 ft/s2
Drag coefficient
0.3
No of stops
13
Battery
40 kW
Distance
6.79 miles
Final drive ratio
3.93
Time
1225 s
Frontal area
1.746 m2
Wheel radius
0.287 m
Vehicle glider mass
918 kg
ECE_EUDC driving cycle: (a) speed required, (b) force and torque required.
4.2. Determination of Efficient Operating Region of Engine
It is mandatory to identify engine’s fuel efficient regions before finding the optimal solution of the cost function. The energy management controller should keep the engine in its efficient region to minimize the liquid fuel consumption. Fuel consumption is a measure of the mass flow per unit time. Fuel flow rate per useful power output is an important parameter to determine the efficiency of the engine and is called specific fuel consumption (SFC); that is, sfc=m˙f/P. When the engine power is measured as the net power from the crankshaft, SFC is called brake specific fuel consumption (BSFC). Low values of SFC or BSFC are always desirable. The ratio of work produced to the amount of fuel energy supplied per cycle is measure of engine efficiency (fuel conversion efficiency). ηf=Wc/mfQHV=P/m˙fQHV, where Wc is work done in one cycle, mf is fuel mass consumed per cycle, and QHV is the heating value of the fuel. The efficiency can be expressed as ηf=1/(sfc∗QHV). Engine characteristics are decided by parameters like power, torque, mean effective pressure, SFC, indicated, brake power and torque, and fuel consumption characteristics.
Fuel efficient region of the engine is mainly governed by requesting power at the ring gear of PSG and maximum and minimum speeds of generator and vehicle idle speed. Based on power demand, optimal ωe∗ and Te∗ points are determined. ωe is controlled with generator torque; that is, generator torque is so adjusted that engine runs at desired speed. Engine maximum (ωe_max) and minimum (ωe_min) speed are ranged using the following equation:(20)NrNs+Nr∗ωring+NsNs+Nr∗ωgmax=ωe_max,NrNs+Nr∗ωring+NsNs+Nr∗ωgmin=ωe_min,where ωring is the speed requested at ring gear. The engine fuel efficiency map is shown in Figure 5 which infers that, below a certain speed, torque produced by the engine is less hence not efficient. ICE is rated at a specific RPM level for maximum torque and maximum power. ICE cannot produce effective torque below “some” certain speed. Maximum torque is achieved for a narrow range of speeds beyond which efficiency decreases. The characteristic of the engine is shown in Figure 6. This characteristic shows that engine’s actual horsepower is lower than the ideal lab conditions further below a certain speed, and no positive torque is achieved. For the considered engine model, maximum power of 43 kw and maximum torque of 101 N-m are provided by engine at 4000 RPM. So, it is required to operate the engine in its most efficient region for the better performance and lesser fuel consumption.
Engine efficiency map of Toyota Prius.
Generalized engine speed-torque characteristics.
4.3. Optimization Strategies
The proposed fuel efficiency optimization problem depends on various parameters of the vehicle. These parameters may have cross effects also. The proposed method uses, firstly, GA to identify optimal values of various governing parameters and then these values are fitted into PMP to produce optimum fuel efficiency.
4.3.1. Genetic Algorithm
To optimize a nonlinear problem using GA, chosen parameters will not be treated as independent variables. The combined effect of these parameters reflects on optimized output. Genetic algorithm was devised by John Holland in early 1970’s to imitate natural properties based on natural evolution. To obtain the solution of a problem, the algorithm is started with a set of solutions known as population. A new population is formed by choosing random solutions of one population and is assumed that new population is better than the old one. This course is repeated over numerous iterations or until some termination criteria is satisfied [37, 38]. The flow of the algorithm is shown in Figure 7.
Genetic algorithm process flow.
4.3.2. Pontryagin’s Minimum Principle
PMP was proposed by Russian mathematician Lev Semenovich in 1956. It gives the best possible control to take a dynamical system from one state to another in the presence of constraints for some state or input control. PMP is a special case of Euler-Lagrange equation of calculus of variations. For an optimum solution, PMP provides only necessary conditions and the sufficient conditions are satisfied by Hamilton-Jacobi-Bellman equation. In PMP, the number of nonlinear second-order differential equations linearly increases with dimension so the control based on PMP takes less computational time for getting an optimal trajectory but it could be a local optimal, not a global solution. Trajectory obtained by PMP could be considered a global optimal trajectory under certain assumptions. These are as follows: (1) trajectory obtained from PMP is unique and satisfies the necessary and boundary conditions, (2) some geometrical properties of the optimal field provide the possibility of optimality clarification, and (3) as a general statement of the second approach the absolute optimality is, mathematically, proven by clear proposition [17, 39].
To optimize any problem using PMP, the Hamiltonian is formed first and then minimized with respect to control input. Then state and costate equations are obtained by following the set procedure. The flow diagram can be corrugated as in Figure 8.
PMP process flow.
For performance measure of the form J=Sxt,ut,t+∫t0tfV(xt,ut,t) with the terminal cost Sxt,ut,t, instantaneous cost ∫t0tfV(xt,ut,t), and the state equation of the form xt=f(xt,ut,t), Hamiltonian constructions involve instantaneous cost and state equation with a time varying vector multiplier λ as follows:(21)Hxt,ut,λt,t=Vxt,ut,t+λTt∗x˙t.
According to PMP, optimal control trajectory u∗t, optimal state trajectory x∗t, and corresponding optimal costate trajectory λ∗t minimize the Hamiltonian such that(22)Hx∗t,u∗t,λ∗t,t≤Hxt,ut,λt,t.The following relations and constraints (23) must hold with the above condition:(23)x˙∗t=∂H∂λx∗t,u∗t,λ∗t,t,λ˙∗t=-∂H∂xx∗t,u∗t,λ∗t,t.Initial conditions x0 and final condition H∗+∂S/∂xtf∂tf+∂S/∂x∗-λ∗(t)tf′δxf both are assumed to be zero.
If PMP conditions are satisfied, the solution will be extrenal and if a global solution exists, it will be the global solution.
5. Strategy Analysis, Simulation, and Result Discussion
The engine in its efficient operating range and motor with sufficient SOC will lead to fuel efficient strategy. Speed, power, SOC, and engine on/off time are the deciding factors and their threshold values must be determined to run an HEV with maximum fuel efficiency.
GA first finds optimal values of engine on SOC, speed, and engine off time (cs_min_off_time, cs_eng_on_soc, cs_electric_launch_spd, and cs_eng_min_spd) thresholds while fulfilling the driver demand; that is, requested trace (road map) should meet at each instant of time over a road trip. Improper values of these parameters will reduce the fuel efficiency. After selecting threshold values of vehicular parameters using GA, they are now fed to PMP which finally reckons threshold power to turn the engine on. The effect of this hybrid control strategy is visible in terms of improved efficiency as shown in Table 3. Four different cases are analyzed here: (1) Rint battery model with conventional SOC estimation used in ADVISOR and (2) Rint, (3) 1 RC, and (4) 2 RC battery models with modified SOC estimation method [29, 31]. A considerable improvement is observed in fuel efficiency using modified SOC estimation method over conventional. Models with modified SOC estimation give 8-9 percent improvement over conventional methods. Modified SOC estimation method with Rint, 1 RC, and 2 RC models do not make much difference in efficiencies as their OCVs, resistances, and capacity variations are close to each other. To take care of the actual battery behavior one should consider R and C components instead of Rint only in HEV analysis. One RC battery model is used here further to avoid the complexity of 2 RC models. Figure 9 provides required conditions to turn the engine on/off. Here, cs_min_pwr decides minimum power commanded of the engine; below this engine should be principally shut off. cs_electric_launch_spd is a vehicle speed threshold below which engine will be off. cs_min_off time is the shortest allowed duration of the engine off period; after this time has passed, the engine may restart if high power is requested. Below cs_eng_on_soc value, the engine must be on. Below cs_eng_min_spd, fuel can be cut; that is, engine does not use fuel.
Fuel economy comparison for different battery models.
Battery model
Fuel economy (mpgge)
Trace analysis
With GA
Without GA
Percentage improvement
Rint conventional model
55.9208
44.9785
24.3278
With trace miss
Rint model∗
60.6184
51.0401
18.7662
No trace miss
1 RC model∗
60.6250
51.3732
18.00
No trace miss
2 RC model∗
60.5423
51.1247
18.428
No trace miss
∗With modified SOC estimation method.
Engine on/off decision.
To verify the correctness of proposed strategy, requested speed and delivered speed of the vehicle are compared and shown in Figure 10. The figure infers that these two match perfectly and there is no trace miss. Vehicle requested power is fulfilled by different components alone or together. Figure 4(b) signifies the time instances of negative torque; that is, kinetic energy (=1/2 MV2) stored in vehicles translating mass can be stored during these moments, if the deceleration rate is greater than 10 km/h. The traction motor operates as generator to recuperates the energy and charges battery as shown in Figure 11. Positive current flow delivers the current from the battery and negative current signifies the condition of battery getting charged.
Vehicle requested and delivered speed comparison.
Battery current over the trip.
Battery SOC variation over the trip and with engine on/off is shown in Figure 12 at 25°C with initial SOC as 80 and target as 70 percent. Figure 13 shows the motor and engine efficiency points and promise to work in most efficient range possible while acquiring the trace and maintaining SOC.
SOC status: (a) SOC variation over the trip and (b) SOC variation with engine on/off condition.
Operating points: (a) engine and (b) motor.
6. Conclusion
In this paper, a modified SOC estimation method is used to track the run-time SOC of the batteries and an optimal control based EMS is developed and implemented to control the engine on/off status. While implementing the strategy, all the important consideration like aerodynamic drag, vehicle glider mass, accessory loads, prescribed SOC level conditions, and so forth are given utmost attention. PMP along with GA and with modified SOC estimation techniques presents promising EMS. Various governing parameters of vehicle are firstly optimized using GA and then a power threshold calculation is performed using PMP. Calculation of thresholds initially using GA gives better chance to improve the fuel efficiency. Here, fuel efficiency is derived for different battery models incorporating modified and conventional SOC estimation methods. This proposed EMS yields better efficiency as compared to the default strategy available.
Conflict of Interests
The authors declare that they have no conflict of interests.
JosG. J.OlivierG. J.-M.JeroenA. H. W.2012SchipperL.FabianH.LeatherJ.Transport and carbon dioxide emissions: forecasts, options analysis, and evaluation20099Asian Development BankJapan Automobile Manufacturers Association Inc2008Japan Automobile Manufacturers Association IncEhsaniM.GaoY.EmadiA.2010chapter 2–9New York, NY, USACRC PressJohnsonV. H.WipkeK. B.RausenD. J.HEV control strategy for real-time optimization of fuel economy and emissions2000109316771690PaganelliG.ErcoleG.BrahmaA.GuezennecY.RizzoniG.General supervisory control policy for the energy optimization of charge-sustaining hybrid electric vehicles2001224511518PaganelliG.TatenoM.BrahmaA.RizzoniG.GuezennecY.Control development for a hybrid-electric sport-utility vehicle: strategy, implementation, and test resultsProceedings of the American Control ConferenceJune 2001Arlington, Va, USA50645069SciarrettaA.BackM.GuzzellaL.Optimal control of parallel hybrid electric vehicles200412335236310.1109/TCST.2004.8243122-s2.0-2442598173DebertM.ColinG.ChamaillardY.GuzzellaL.Ketfi-CherifA.BellicaudB.Predictive energy management for hybrid electric vehicles—prediction horizon and battery capacity sensitivityProceedings of the 6th IFAC Symposium Advances in Automotive Control (AAC '10)July 201027027510.3182/20100712-3-de-2013.000662-s2.0-80051907416BeckR.RichertF.BolligA.AbelD.SaengerllS.NeilK.ScholtT.NoreikatK.-E.Model predictive control of a parallel hybrid vehicle drivetrainProceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference (CDC-ECC '05)December 2005IEEE2670267510.1109/cdc.2005.15825662-s2.0-33847217283ArsieI.GraziosiM.PianeseC.RizzoG.SorrentinoM.Optimization of supervisory control strategy for parallel hybrid vehicle with provisional load estimateProceedings of the 7th International Symposium on Advanced Vehicle Control (AVEC '04)August 2004Arnhem, The Netherlands483488ProkhorovD.Toyota prius HEV neurocontrolProceedings of the International Joint Conference on Neural Networks (IJCNN '07)August 2007Orlando, Fla, USAIEEE2129213410.1109/ijcnn.2007.43712872-s2.0-51749100642HuangM.YuH.Optimal multilevel hierarchical control strategy for parallel hybrid electric vehicleProceedings of the IEEE Conference Vehicle Power and Propulsion (VPPC '06)September 2006Windsor, UK1410.1109/VPPC.2006.364369HuangM.YuH.Optimal control strategy based on PSO for powertrain of parallel hybrid electric vehicleProceedings of the IEEE International Conference on Vehicular Electronics and Safety (ICVES '06)December 2006Beijing, ChinaIEEE35235510.1109/icves.2006.3716142-s2.0-50249178125WangZ.HuangB.LiW.XuY.Particle swarm optimization for operational parameters of series hybrid electric vehicleProceedings of the IEEE International Conference Robotics and BiomimeticsDecember 2006Kunming, China682688SerraoL.RizzoniG.Optimal control of power split for a hybrid electric refuse vehicleProceedings of the American Control Conference (ACC '08)June 2008Seattle, Wash, USA4498450310.1109/acc.2008.45872042-s2.0-52449103553KimN.LeeD.Cha SW.PengH.Optimal control of a plug-in hybrid electric vehicle (PHEV) based on driving patternsProceedings of the International Battery, Hybrid and Fuel Cell Electric Vehicle SymposiumMay 2009Stavanger, Norway19StockarS.MaranoV.RizzoniG.GuzzellaL.Optimal control for plug-in hybrid electric vehicle applicationsProceedings of the American Control Conference (ACC '10)July 2010Baltimore, Md, USA502450302-s2.0-77957764911StockarS.MaranoV.CanovaM.RizzoniG.GuzzellaL.Energy-optimal control of plug-in hybrid electric vehicles for real-world driving cycles20116072949296210.1109/tvt.2011.21585652-s2.0-80052856732KimN.RousseauA.LeeD.A jump condition of PMP-based control for PHEVs201119623103801038610.1016/j.jpowsour.2011.07.0032-s2.0-80053567515KimN.ChaS. W.PengH.Optimal equivalent fuel consumption for hybrid electric vehicles201220381782510.1109/tcst.2011.21230992-s2.0-84859883656WipkeK. B.CuddyM. R.BurchS. D.ADVISOR 2.1: a user-friendly advanced powertrain simulation using a combined backward/forward approach19994861751176110.1109/25.8067672-s2.0-0033330689PiccoloA.IppolitoL.GaldiV.VaccaroA.Optimisation of energy flow management in hybrid electric vehicles via genetic algorithms1Proceedings of the IEEE/ASME International Conference on Advanced Intelligent MechatronicsJuly 2001Como, Italy43443910.1109/AIM.2001.936493WangA.YangW.Design of energy management strategy in hybrid electric vehicles by evolutionary fuzzy system Part II: tuning fuzzy controller by genetic algorithmsProceedings of the 6th World Congress on Intelligent Control and Automation (WCICA '06)2006Dalian, China8324832810.1109/WCICA.2006.1713600HuangB.ShiX.XuY.Parameter optimization of power control strategy for series hybrid electric vehicleProceedings of the IEEE Congress on Evolutionary Computation (CEC '06)July 2006Vancouver, Canada198919942-s2.0-34547252484WimalendraR. S.UdawattaL.EdirisingheE. M. C. P.KarunarathnaS.Determination of maximum possible fuel economy of HEV for known drive cycle: genetic algorithm based approachProceedings of the 4th International Conference on Information and Automation for Sustainability (ICIAFS '08)December 2008Colombo, Sri LankaIEEE28929410.1109/iciafs.2008.47839752-s2.0-64049101346TangX.MaoX.LinJ.KochB.Li-ion battery parameter estimation for state of chargeProceedings of the IEEE American Control Conference (ACC '11)July 2011IEEE9419462-s2.0-80053145193VerbruggeM.TateE.Adaptive state of charge algorithm for nickel metal hydride batteries including hysteresis phenomena20041261-223624910.1016/j.jpowsour.2003.08.0422-s2.0-0742285682PandayA.BansalH. O.Temperature dependent circuit-based modeling of high power Li-ion battery for plug-in hybrid electrical vehiclesProceedings of the International Conference on Advances in Technology and Engineering (ICATE '13)January 2013Mumbai, IndiaIEEE1610.1109/icadte.2013.65247372-s2.0-84879663122PandayA.BansalH. O.Hybrid electric vehicle performance analysis under various temperature conditions2015751962196710.1016/j.egypro.2015.07.238PandayA.BansalH. O.SrinivasanP.Thermoelectric modeling and online SOC estimation of Li-ion battery for plug-in hybrid electric vehicles2016201612235352110.1155/2016/2353521CliffsE.19912ndPrentice-HallConwayB. E.Transition from ‘Supercapacitor’ to ‘Battery’ behavior in electrochemical energy storage199113861539154810.1149/1.20858292-s2.0-0026168744ChenM.Rincón-MoraG. A.Accurate electrical battery model capable of predicting runtime and I-V performance200621250451110.1109/tec.2006.8742292-s2.0-33744765213LiuJ.PengH.FilipiZ.Modeling and analysis of the Toyota hybrid systemProceedings of the IEEE/ASME International Conference on Advanced Intelligent MechatronicsJuly 2005Monterey, Calif, USAIEEE1341392-s2.0-27644512002MiC.MasrurM. A.GaoD. W.2011London, UKJohn Wiley & SonsSumathiS.SurekhaP.2010chapter 13New York, NY, USACRC PressDebK.Practical optimization using evolutionary methods20052005008KrotovV. F.1996New York, NY, USAMarcel Dekker