Next Article in Journal
A Systematic Literature Review on Fake News in the COVID-19 Pandemic: Can AI Propose a Solution?
Next Article in Special Issue
A Study on Longitudinal Motion Scenario Design for Verification of Advanced Driver Assistance Systems and Autonomous Driving Systems
Previous Article in Journal
High-Precision Numerical Research on Flow and Structure Noise of Underwater Vehicle
Previous Article in Special Issue
Cooperative Roundabout Control Strategy for Connected and Autonomous Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New HEV Power Distribution Algorithm Using Nonlinear Programming

1
Department of Electrical Engineering, Hanyang University, Seoul 04763, Republic of Korea
2
Department of Electrical and Biomedical Engineering, Hanyang University, Seoul 04763, Republic of Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(24), 12724; https://doi.org/10.3390/app122412724
Submission received: 15 November 2022 / Revised: 3 December 2022 / Accepted: 8 December 2022 / Published: 12 December 2022

Abstract

:
An equivalent consumption minimization strategy (ECMS) is one of the most powerful and practical ways to improve the fuel efficiency of hybrid electric vehicles (HEVs). In an ECMS, it is important to determine the optimal equivalent factor to reach a global optimal solution. The optimal equivalent factor is determined by driving conditions. Previous studies have used an adaptive ECMS (A-ECMS) to determine the appropriate equivalent factor according to changing driving conditions. An A-ECMS adjusts the equivalent factor by controlling the battery’s state of charge (SOC) to follow a reference SOC trajectory. It is therefore critical to identify a reference SOC trajectory that reflects real-world driving conditions. These conditions, which are composed of the HEV’s nonlinear dynamics and complex constraints, can be formulated into a nonlinear optimal control problem (NOCP). Here, we propose applying nonlinear programming (NLP) to an A-ECMS. The NLP-based ECMS algorithm can be divided into two parts: the use of an NLP to solve an NOCP to obtain the reference SOC trajectory and the application of an NLP solution (the result of the first part) to an A-ECMS. Simulation results demonstrate that the proposed NLP-based ECMS closely resembles a global optimal solution for dynamic programming in a relatively brief calculation time.

1. Introduction

As global warming accelerates, environmental regulations are being introduced and enforced around the world. In the United States, the National Highway Traffic Safety Administration has developed corporate average fuel economy standards that require an average fuel efficiency of 23.2 km/L by 2025. In Europe, the Euro 7 regulations for vehicle emissions also require nitrogen oxide emission to be lowered to 40–120 mg/kWh by 2025 [1]. In response to strong environmental regulations, the automobile industry is making great efforts to electrify the vehicle powertrain.
Hybrid vehicles that provide propulsion power from two or more sources are typical examples of electrified powertrains. Hybrid electric vehicles (HEVs) are generally powered by both an internal combustion engine (ICE) and an electric motor (EM). HEV can reduce fuel consumption by operating ICE and EM in high-efficiency areas. When combined with an ICE-based powertrain, an EM can support the ICE in low-speed sections that are associated with low efficiencies and can regenerate energy during braking. As a result, an effective energy management strategy (EMS) that reflects driving and vehicle conditions is important to improve the efficiency of a HEV.
The study of EMSs of HEVs can be classified into two main categories. The first comprises rule-based EMSs, in which rules based on the experience of engineers determine the power distribution ratio between the ICE and EM. The advantage of a rule-based EMS is its intuitiveness. However, due to its innate rigidity, it is impossible to design a rule that anticipates every possible driving situation, and fuel efficiency may deteriorate in situations that are not reflected in the rule [2,3]. In addition, a rule-based EMS requires considerable time and resources as intensive parameter-tuning is required to anticipate diverse driving situations [4]. The second category comprises optimization-based EMSs, which use optimization theories to determine the power distribution ratio between the ICE and EM. An optimization-based EMS can determine the power distribution ratio closest to a global optimal solution. However, such a strategy can be complicated and require lengthy computational times due to complex optimization problems.
Various studies have been conducted on optimization-based EMSs because the hardware performance of the hybrid control unit (HCU) has been increased and fuel consumption can be minimized by using a global optimal solution in various driving situations. In an optimization-based EMS, an optimal control problem (OCP) can be defined and calculated by using optimization theory. Depending on how it is defined, or which optimization theory is applied, an OCP can be divided into offline and online optimization algorithms.
An offline optimization algorithm produces a global optimal solution that can minimize fuel consumption in all expected driving conditions, assuming that all conditions are known in advance. The most representative algorithm is dynamic programming (DP), which is based on the principle of optimality proposed by Richard Bellman. After splitting the complex problem into several simpler subproblems, the solution of the subproblems is sequentially obtained in a recursive manner, so that a complete algorithm obtains the global optimal solution. Although DP can identify a global optimal solution, it is difficult to apply such solutions to a real-world HCU because DP requires lots of computational resources. A global optimal solution obtained by DP is used primarily in offline analysis in the development of other EMS options, such as rule-based designs or as machine-learning training data [5].
Online optimization produces a solution that can minimize fuel consumption for a given driving condition in a limited time. In general, it can be applied to a real-world HCU by reducing the computational load by using a simplification method (e.g., linearization) or local optimization. Pontryagin’s minimum principle (PMP), model predictive control (MPC), and equivalent consumption minimization strategy (ECMS) are representative algorithms. A PMP performs local optimization by applying a numerical method to a Hamiltonian function, which considers the battery’s state of charge (SOC) and the fuel consumption characteristics of the ICE [6,7,8]. A PMP needs to efficiently design the Hamiltonian function to reduce computational loads. An MPC obtains an optimal solution for a given driving condition on a finite horizon. For an infinite horizon, an OCP is defined by designing cost functions and constraints that represent the characteristics of the system. A numerical solution can be obtained by applying a solver, such as quadratic programming and nonlinear programming (NLP) to a defined OCP [9,10]. An MPC needs to design system models and constraints carefully to reduce computational loads and create a coping mechanism for infeasible regions. Recently, research on the use of nonlinear-MPC (NMPC) by NLP solver has also been actively conducted. NMPC directly determines the torque of the engine and motor. NLP solvers used in NMPC include the sequential quadratic programming (SQP) method [11,12], mixed-integer nonlinear programming (MINLP) method [13], and the interior point method [14]. NMPC is one of the most efficient and powerful online optimization methods. When designing an NMPC, two analyses are most important. The first is an analysis for model simplification to reduce long computational time. The second is an analysis of the infeasible region to design a strategy for the infeasible solution. An ECMS derives a solution that minimizes equivalent energy consumption, which is the equivalent energy consumption in the form of a cost function that reflects fuel energy and electrical energy. This cost function defines the equivalent factor, which confirms the relative importance of electrical energy and fuel energy, according to driving conditions. The equivalent factor is a therefore an important tuning factor when designing an ECMS. A traditional ECMS uses an equivalent factor with a fixed value [15,16]. However, the equivalent factor can be sensitive to driving conditions [17,18,19,20]. An adaptive-ECMS (A-ECMS) has been proposed to regulate the equivalent factor according to changes in driving conditions. An A-ECMS adjusts the equivalent factor so that the SOC can follow the reference SOC trajectory while considering driving conditions. To obtain the reference SOC trajectory, researchers have applied particle swarm optimization (PSO) [21,22,23], approximate DP [24,25,26], and neural networks (NN) [27,28,29,30,31].
An A-ECMS is one of the most suitable EMS options to use in real-world HCUs [32,33,34,35,36]. However, A-ECMS designs have some drawbacks. First, computational time must be minimized to operate in a real-world HCU. An A-ECMS compares the control input candidates presented in a complex solution space to determine the power distribution ratio. Because the calculation of many control input candidates should reflect the characteristics of the HEV, system models and constraints need to be efficiently designed for computation. If the system model and constraints are too complex or poorly designed, the computational time increases significantly. Determining the reference SOC trajectory accurately and rapidly is also challenging. An A-ECMS adjusts the equivalent factor according to the reference SOC trajectory. If the accuracy of the reference SOC trajectory is insufficient, fuel consumption cannot be minimized. Even an accurate reference SOC trajectory cannot be used for an A-ECMS if it requires a long computational time. Only accurate and rapid calculation of a reference SOC trajectory can maximize the performance of an A-ECMS.
In this paper, an NLP-based ECMS is proposed. NLP-based ECMS does not directly determine the torque of the engine and motor. The equivalent factor is determined by NLP, and the torque of the engine and motor is determined by the A-ECMS scheme using the equivalent factor. Therefore, the OCP size of the NLP can be reduced, and stable control of the vehicle is possible, even if the NLP has derived an infeasible solution. As mentioned above, it is important to obtain the reference SOC trajectory accurately and rapidly, and NLP can provide rapid and accurate calculations. To reduce NLP computation time, an optimal system model and constraints for an HEV with a P2 structure (TMED) was designed. Results of simulations showed that the proposed NLP-based ECMS closely resembled the global optimal solution produced by DP in a relatively short calculation time.
The main contributions of this paper are as follows. First, a suitable required torque preprocessing method for a P2 structure is proposed. In this process, the required torque demanded by the ICE and EM according to the change in gear ratio is calculated. Secondly, a method for calculating the reference SOC trajectory using NLP is proposed. System models, constraints, and cost functions were designed to reflect the characteristics of the HEV with a P2 structure and driving conditions. To reduce the computational time, the NLP model was simplified. Finally, an A-ECMS that can efficiently follow the reference SOC trajectory is proposed as an optimization-based EMS resembling the global optimal solution of DP.
This paper consists of the following structure. Section 2 presents a system model of a HEV with a P2 structure, Section 3 presents analysis result of the NLP-based ECMS, Section 4 discusses the significance of our findings, and Section 5 is discussion of the result.

2. System Model

In this paper, the NLP-based ECMS was applied to HEVs with a P2 structure, which includes a four-cylinder ICE, a permanent-magnet synchronous motor (PMSM), a lithium-ion battery, and a five-speed automatic transmission. The configuration of the P2 hybrid powertrain is shown in Figure 1. An HEV with a P2 powertrain was simulated by using Autonomie [37], a forward-looking simulation tool developed by Argonne. The NLP-based ECMS algorithm was designed to use Autonomie’s data. In addition, an Autonomie simulation environment was used to analyze the performance of the NLP-based ECMS. Table 1 lists the specifications of the major components of the Autonomie model.

2.1. Internal Combustion Engine

An EMS is designed to minimize fuel consumption of an HEV. The ICE of an HEV is modeled as quasistationary, represented by a brake-specific fuel consumption (BSFC) map configured by using the ICE’s steady-state test data, as shown in Figure 2. A BSFC map is formulated as a function of torque ( T e n g ) and speed ( w e n g ), as shown in Equation (1). The relationship between BSFC and the power ( P e n g ) is shown in Equation (2). Q L H V is the lower heating value. We have
b s f c = f T e n g ,   w e n g
P e n g = b s f c × Q L H V .
The quasistationary model cannot consider the transient response of the ICE. The transient response of an ICE adversely affects fuel consumption. The transient response of an ICE has complex dynamics, which are difficult to reflect in an EMS accurately. The transient response of the ICE is simplified as a first-order delay function to reflect it in the EMS, as shown in Equation (3). As a result, when the input torque ( T e , c m d ) is determined, the actual torque ( T e n g ) passes through the low-efficiency area of the BSFC map:
T ˙ e n g = 1 τ e n g T e , c m d T e n g .

2.2. Electrical Motor

An HEV is being driven while maintaining the SOC at a certain level. The SOC is determined by the current of the EM, when excluding the influence of accessories. The current of the EM is determined by power. The power of the EM was represented in a quasistationary model in the form of a map. The power map was expressed as a function of torque ( T m o t ) and speed ( w m o t ), as shown in Equation (4). The power of an EM also has a nonlinear characteristic, as shown in Figure 3. We have
P m o t = f T m o t ,   w m o t .
Because the power map of an EM is also quasistationary, it cannot consider a transient response. A transient response is simplified as a first-order delay function, as shown in Equation (5). The time delay ( τ e n g , τ m o t ) should be decided by reflecting how quickly the EM is compared to the ICE. As a result, the input torque ( T m , c m d ) passes through the low-efficiency area of the power map when the actual torque is determined. We have
T ˙ m o t = 1 τ m o t T m , c m d T m o t .

2.3. Battery

The battery is a complex nonlinear system that is affected by temperature, voltage, and SOC. It is difficult to build a model that accurately represents the battery’s dynamic characteristics in an EMS because the formula is complicated and takes a long time to compute. Therefore, the internal resistance model, which is one of the control-oriented models, was applied to the EMS. This model calculates the battery SOC by defining the equivalent circuit expressed by internal resistance ( R i n t ), open-circuit voltage ( V o c ), and EM power ( P m o t ). The formula is shown in Equation (6). In this formula, R i n t and V o c are functions of SOC:
I m o t = V o c V o c 2 4 · R i n t · P m o t 2 · R i n t .
The SOC of the battery can be calculated integrating the current used by the EM, as shown in Equation (7). This formula is calculated by the initial SOC ( S O C i n i t ) and battery capacity ( Q b a t ):
S O C = S O C i n i t 1 Q b a t · I m o t t d t .

3. NLP-Based ECMS Algorithm

In this paper, the NLP-based ECMS algorithm, a combination of NLP and an A-ECMS, is proposed. A conceptual diagram of the NLP-based ECMS algorithm is shown in Figure 4. Based on the reference SOC trajectory obtained by NLP, the A-ECMS operates to minimize fuel consumption of the HEV. The process of performing the NLP-based ECMS algorithm is as follows. First, the required torque demanded by ICE and EM is calculated according to the driving condition. The gear ratio must be properly considered when calculating the required torque, because the gear ratio of the P2 powertrain is a discrete variable. Secondly, the reference SOC trajectory corresponding to the required torque through NLP is calculated. NLP can cause long calculation times and infeasible solutions if the system model and constraints are poorly designed. Thirdly, the A-ECMS determines the control input so that the measured SOC can follow the reference SOC trajectory.

3.1. Required Torque Preprocessing Algorithm

The required torque preprocessing algorithm determined the required torque used to calculate the reference SOC trajectory in NLP. The required torque preprocessing algorithm consisted of three steps. In first step, the speed was calculated along with the acceleration of the wheel by using the vehicle-speed and road-grade profiles. Inertial force, grade resistance force, air resistance force, rolling resistance torque, and final drive efficiency were reflected. In the second step, the gear ratio trajectory was determined by using the gear shift rule, which is a function of the speed and acceleration determined in the first step. A map for gear shift rules is provided in Figure 5. In third step, the required torque was determined by the gear ratio and transmission efficiency. The algorithm calculated the gear ratio and the torque which will be distributed to the EM and ICE.

3.2. Nonlinear Programming

The NLP-based ECMS algorithm calculated the reference SOC trajectory. NLP then calculated the solution for a nonlinear optimal control problem (NOCP) by using a solver. The NOCP was derived from the vehicle model and driving conditions. The IPOPT solver calculated the solution for the NOCP. The solution for the NOCP, the result of IPOPT, was the reference SOC trajectory. The more accurately and efficiently the NOCP was defined, the higher the quality of the solution provided by NLP derived from the solver [38,39,40,41,42].
The NOCP can be explained in three parts. The first part involves a nonlinear system dynamics model. The state vector, input vector, and known-disturbance vector of the HEV model is shown in Equation (8). The system dynamics equations of the components are represented by state space by using vectors, as shown in Equation (9). As part of the simplification work to accelerate the NLP computation, the speed of the EM was treated as a known disturbance. This is because the solution of the NLP, which was used indirectly in the form of the SOC, was determined to satisfy the required torque derived from the EM speed. We have
x = T e n g T m o t S O C ,   u = T c m d ,   e n g T c m d ,   m o t ,   and   d = w m o t = w r e f
x ˙ = T ˙ e n g T ˙ m o t S O C ˙ = 1 τ e n g · T e n g + T c m d ,   e n g 1 τ m o t · T m o t + T c m d , m o t 1 Q b a t · V o c V o c 2 4 · R i n t · T m o t · w m o t 2 · R i n t = f x , u , d .
The second part involves nonlinear constraints. The system model was used in NLP to predict how the HEV would behave within the time horizon. That is, the system model operated as a nonlinear constraint when the NLP calculated the solution. The NLP-based ECMS includes constraints of the time horizon by applying the system model to a fourth order Runge–Kutta method, which is a discretization method. The fourth-order Runge–Kutta method is shown in Equation (10). We have
x k + 1 = x k + T s t e p 6 · k 1 + 2 · k 2 + 2 · k 3 + k 4 k 1 = f x k ,   u k ,   d k k 2 = f x k + 0.5 · T s t e p · k 1 ,   u k ,   d k k 3 = f x k + 0.5 · T s t e p · k 2 ,   u k ,   d k k 4 = f x k + T s t e p · k 3 ,   u k ,   d k .
The Runge–Kutta method, which requires a relatively large amount of computation, was applied rather than the Euler method, which has a small amount of computation, because the step time ( T s t e p ) can be relatively large. A large step time can predict longer intervals using the same number of receiving horizons. Equation (10) can be reformulated to equality constraints, as shown in Equation (11). The first nonlinear constraint was derived from a system dynamics model. We have
g 1 x = x k + 1 x k + T s t e p 6 · k 1 + 2 · k 2 + 2 · k 3 + k 4 = 0 .
If the described Runge–Kutta method is used, the initial value of the state vector is important because it is accumulated from the initial value. That means the initial values of the state vector are also equality constraints. The initial values of engine torque, motor torque, and SOC were set as shown in Equation (12). Because the SOC must be situated within a certain range at the end of driving, the final value was also important for the SOC. The final value of the SOC was also set to an equality constraint as in Equation (13). We have
g 2 x = x 1 0 = T e n g 0 = 0 g 3 x = x 2 0 = T m o t 0 = 0   g 4 x = x 3 0 = S O C 0 = 0.6
g 5 x = x 3 T f i n = S O C f i n .
The maximum/minimum torque is determined by the speed ( w r e f ). Because the speed ( w r e f ) is a known disturbance, the maximum/minimum torque can be known in advance and is an inequality constraint, as shown in Equation (14):
g 6 x f e n g , m i n w r e f T e n g f e n g , m a x w r e f g 7 x f m o t , m i n w r e f T m o t f m o t , m a x w r e f .
The third part involves a cost function. The algorithm includes a term for fuel consumption and a term for tracking the required torque by using the sum of ICE torque ( T e n g )   and EM torque ( T m o t ) . The cost function of an NLP-based ECMS is shown in Equation (15). Q f u e l is the weight cost of fuel consumption ( c o s t f u e l ), and Q t r q is the weight cost of the reference torque ( c o s t t r q ). We have
c o s t = Q f u e l · c o s t f u e l + Q t r q · c o s t t r q .
The cost of fuel consumption can be calculated by using Equation (16). Surface-fitting for a fuel-consumption rate map was performed to include fuel consumption in the cost function of NLP. To derive a polynomial with an appropriate degree, various models were analyzed as shown in Figure 6. We have
c o s t f u e l = s u m p 20 · T e n g 2 + p 10 · T e n g + p 11 · T e n g · w r e f + p 01 · w r e f .
The cost for tracking the required torque using the sum of ICE torque ( T e n g ) and EM torque ( T m o t ) can be calculated by using Equation (17). The required torque is the reference torque given to the ICE and EM. The cost function was designed to use reference torque, not speed, to decrease the computational time by reducing the size of the system model. We have
c o s t t r q = s u m T e n g + T m o t T r e f 2 .
The constraint and cost functions were defined in the form of an OCP for receding horizons. The defined problem was solved by an NLP solver. The IPOPT solver, which is a rapid and powerful NLP solver, was used for the NLP-based ECMS algorithm.

3.3. Adaptive Equivalent Consumption Minimization Strategy

The NLP-based ECMS algorithm was controlled by an A-ECMS based on the reference SOC trajectory, which is the result of NLP. An A-ECMS computes the cost function for input grids and obtains the optimal power distribution ratio, which minimize the cost function. A schematic diagram of A-ECMS consisting of five steps is shown in Figure 7.
The first step determined the number of grids for the control input candidates ( u g r i d ) for the power distribution ratio, as described by Equation (18). The A-ECMS in this paper used 201 grids (N = 201). We have
u g r i d = u N × 1 = 1 ,   1 .
The second step calculated the torque candidate of the ICE and EM. The required torque ( T r e q ), calculated by a preprocessing algorithm, was distributed to the EM torque candidates ( T m o t ,   g r i d ) by the input grid. The engine torque candidates ( T e n g ,   g r i d ) were based on the difference between the required torque and EM torque candidates. The formulas are found in Equations (19) and (20). In this process, the maximum/minimum torque of the EM and the ICE was considered as shown in Equation (21):
T m o t ,   g r i d = T r e q · u g r i d
T e n g ,   g r i d = T r e q T m o t , g r i d
s u b j .   t o .   T m o t , m i n T m o t , g r i d T m o t , m a x T e n g , m i n T e n g , g r i d T e n g , m a x .
The third step calculated the power candidate of the ICE and EM. The ICE power candidates ( P e n g ,   g r i d ) were obtained from the BSFC ( m ˙ b s f c ) and lower heating value ( Q L H V ). The EM power candidates ( P m o t ,   g r i d ) were obtained from the internal resistance ( R i n t ), open-circuit voltage ( V o c ), and EM current ( I m o t ). Internal resistance and open-circuit voltage are functions of the battery’s SOC, and the EM current can be calculated by using Equation (22). SOC candidates were calculated from the EM power candidates because the EM power candidates affected the SOC. The formulas are as follows Equations (23)–(25):
I m o t , g r i d = V o c V o c 2 4 · R i n t · T m o t , g r i d · w m o t 2 · R i n t
P m o t ,   g r i d = V o c · I m o t , g r i d R i n t · I m o t , g r i d 2
P e n g   g r i d = Q L H V · m ˙ b s f c T e n g , g r i d ,   w e n g
S O C   g r i d = ( 1 / Q b a t ) · I m o t , g r i d .
The fourth step calculates the equivalent factor ( S e q ), which was used to equalize the power candidates of the ICE and EM. This process determines whether to use power from the EM or the ICE. An A-ECMS regulates the equivalent factor, so that the battery’s SOC follows the reference SOC trajectory, which can be calculated by NLP. If the battery’s SOC exceeded the reference trajectory, the equivalent factor was decreased. If the battery’s SOC was less than reference trajectory, the equivalent factor was increased. This is a form of a proportional–integral–derivative controller ( S p i d ), as shown in Equations (26) and (27). Equation (28) was used to prevent the equivalent factor from changing rapidly [33,38]. We have
S O C e r r = S O C r e f S O C g r i d
S p i d = S 0 + K p · S O C e r r + K i · S O C e r r + K d · d S O C e r r d t
S e q = S p i d l 1 · t a n π 2 · l s · S O C e r r Δ S O C .
The final step calculates the cost function and derives the optimal input, which is the power-split ratio between the ICE and EM that minimizes the cost function from the candidates. The cost function ( J e c m s ) of an A-ECMS reflecting the equivalent factor is described by Equation (29). The power candidates of the EM were equalized by using the equivalent factor as the power candidates of the ICE. We have
J e c m s = S e q · P m o t , g r i d + P e n g , g r i d .
An A-ECMS is a quasistationary method that uses an input grid. Because the transient response of the ICE was not considered, the ICE torque changed rapidly or was turned on and off frequently. To prevent this problem, Equations (30) and (31) were used to reflect the cost function. Equation (30) is a cost function for engine torque change in previous ( T e n g ) and current steps ( T e n g , g r i d ). Equation (31) is the cost of starting an engine that was in EV mode in the previous step ( u c m d ):
J t r q = T e n g , g r i d · w e n g T e n g · w e n g 2
J s t r = u g r i d 1 · u c m d = 1 .
In addition, the torque candidates of the ICE and EM are constrained by maximum/minimum torque conditions. If the engine torque and motor torque do not satisfy the required torque because of the constraints, the vehicle cannot follow the speed profile. To prevent this problem, Equation (32) was applied to reflect the cost function by using the difference between the power components (EM and ICE) and the required power:
J r e f = a b s T r e q T m o t T e n g · w g b .
The cost function of A-ECMS ( J ) is shown in Equation (33), and the optimal solution of A-ECMS ( u * ) can be obtained from Equation (34). Each term that constitutes the cost function must be tuned by using a weight ( W 1 ,   W 2 ,   W 3 ) . When setting the weight, it is important to prevent contamination of the original cost function ( J e c m s ) with the ICE power candidate, EM power candidate, and equivalent factor. We have
J = J e c m s + W 1 · J t r q + W 2 · J s t r + W 3 · J r e f
u * = m i n i m i z e u g r i d   J .

4. Results

In this paper, the performance of an NLP-based ECMS was validated in a simulation environment by using Autonomie, a forward-looking simulation tool developed by Argonne. An HEV model (P2 powertrain) and an NLP-based EMS were connected in a Simulink environment, as shown in Figure 8. The FTP-75 driving cycle was used to carry out the simulation and verify the performance. Simulations were run on a desktop computer equipped with an AMD Ryzen 5 5600X 6-Core Processor (3.70 GHz) and 16 GB of RAM.
Dynamic programming was used to compare the performance of the NLP-based ECMS. As mentioned above, DP was used to calculate a global optimal solution. That is, an optimal solution that can minimize fuel consumption for a given speed profile. If NLP-based ECMS ensures the HEV follows the optimal trajectory, the result of the NLP-based ECMS will get close to that of DP.
In this paper, the reductions in fuel consumption in the NLP-based ECMS were compared with the results of a rule-based EMS, DP, and DP-based ECMS. The computational times of each EMS were also compared.
Section 4.1 describes the DP model, which is a comparative group, and analyzes the results of the DP model. Section 4.2 analyzes the performance of NLP’s reference SOC trajectory calculation. Section 4.3 analyzes the performance of the NLP-based ECMS by using the HEV vehicle model of Autonomie.

4.1. Comparison Group (Global Optimal Solution/Dynamic Programming)

First, we analyzed the results of deriving the reference SOC trajectory of the NLP-based ECMS. Because an NLP-based ECMS derives the reference SOC trajectory by using NLP, the results of NLP and DP could be compared. The DP setting as a comparison target is shown in Table 2. The input grid involves the power-split ratio, which ranges from −1 to 1. The state grid is for the SOC, which ranges from 0.3 to 0.9. As the number of input grids and state grids increased, the calculation time increased. As the step time decreased, the calculation time increased significantly.
As shown in Figure 9 and Figure 10, the DP result changed significantly according to the setting for the input grid, the state grid, and the step time. Figure 9 displays the results of DP calculations for the SOC trajectory, which differs greatly depending on the case. Figure 10 depicts the result of DP calculations for ICE torque. If the grid and step time were set incorrectly, the calculated engine torque would include ripple.
Case 0 is the DP result with the largest grid and the smallest step time set. Case 0 had the longest calculation time because of the setting, but the solution was the closest to the global optimal. Case 0 was therefore used as the reference result, which is an ideal result that minimizes fuel consumption during the speed profile. Case 0 cannot be applied to the ECMS because it took too much time to calculate. Accordingly, Case 2, Case 5, and Case 6, which required relatively short computation times, were further compared. The results of Case 2, Case 5, and Case 6 were applied to the ECMS to implement a DP-based ECMS and compared it with the NLP-based ECMS proposed in this paper. Case 2, Case 5, and Case 6, which have relatively short calculation times, were also included in the comparison group. The results of Case 2, Case 5, and Case 6 were applied to the ECMS to design a DP-based ECMS and compared them with the NLP-based ECMS proposed in this paper.

4.2. Reference SOC Trajectory

The NLP-based ECMS calculated the reference SOC trajectory by using NLP. Accuracy and computational time were important factors. Accuracy indicates how close the results of NLP were to the global optimal solution.
The results of NLP are compared with DP models in the comparison group. Table 3 compares the NLP and DP results. The calculation time is the time spent calculating the reference SOC trajectory. The sum of error is the difference between the global optimal solution and the reference SOC trajectory obtained by DP and NLP. As a result, NLP was the quickest and the most accurate to calculate the reference SOC trajectory.
Case 0, Case 2, Case 5, and Case 6 are results for the reference framework obtained by the DP model. Case 0 takes the longest time because it has the largest number of input/state grids and the smallest step size. Because the most accurate DP model is used, Case 0 is used as a reference because it is the closest result to a global optimal solution.
In order to be applied to DP based ECMS, the DP model of Case 2, Case 5, and Case 6 has been simplified. Fewer input/state girds and larger step sizes were used, which can reduce the computation time. Because the model is simplified, however, the results were moved away from the global optimal solution. The cost function of the DP is the fuel consumption of the HEV as shown in Equation (1). The cost function used in the DP of each case and the system model (the equations, parameters, maps, etc. shown in Section 2) are also identical.
The NLP case is a result obtained by NLP. The cost function and system models of NLP and DP are different. The cost function and system model of NLP are shown in Section 3.2. However, because the target vehicle is same, same vehicle parameters are used as shown in Section 2.
The tendency of the reference SOC trajectory is important because it determines the charging or discharging of the EM. As shown in Figure 11, the reference SOC trajectory obtained by NLP was similar to the global optimal solution. As shown in Figure 12, the states of ICE ( x 1 ) and EM torque ( x 2 ) were found by the NLP input to be first-order delay functions. Additionally, the states of ICE and EM torque were limited by constraints on the maximum/minimum ICE and EM torque.

4.3. NLP-Based ECMS

Figure 13 depicts the results of the simulation by applying the NLP-based ECMS to Autonomie’s HEV model. The NLP-based ECMS properly determined the ICE torque and EM torque so that the HEV follows the speed profile. In addition, the NLP-based ECMS appropriately determined the charging and discharging mode of the EM, so that the battery SOC follows the reference SOC trajectory calculated by NLP. An error occurred between the engine torque and motor torque and command because of the influence of clutch and transmission.
Table 4 reveals that the NLP-based ECMS reduced fuel consumption in a manner similar to the Case-0 DP-based ECMS with the shortest simulation time. The Case-0 DP-based ECMS is the closest to the global optimal solution because the A-ECMS used the most accurate DP result. An NLP-based ECMS is an EMS that optimally controls HEV.

5. Discussion

An NLP-based ECMS that can determine an appropriate equivalent factor for an A-ECMS was proposed. The NLP-based ECMS used a reference SOC trajectory to control the equivalent factor. The NOCP was designed to calculate the reference SOC trajectory. The NOCP includes ICE torque, EM torque, and battery SOC as state variables, and the system models for these variables were applied in a simplified form to increase calculation speed. The NOCP also includes a cost function, which was designed through surface-fitting to the fuel-consumption-rate map. The solution closest to the global optimal with a short computation time through NLP-based ECMS was preferred. In addition, the proposed method was able to respond to infeasible solutions that occurred in NLP by using a strategy involving the equivalent factor. The performance of the proposed method, which was verified in an Autonomie simulation environment, was superior to those of a DP-based ECMS and rule-based EMS in terms of calculation time and fuel consumption.
The proposed NLP-based ECMS completed its tasks rapidly in a desktop-computer simulation. However, the hybrid or electric control unit of an actual vehicle likely has far fewer computational resources. To apply an NLP-based ECMS to a real-world environment, a solver that can be applied to embedded system is therefore required. In this paper, the possibility that NLP-based ECMS can reduce fuel consumption close to the global optimal solution (DP result) is confirmed. However, empirical research on NLP-based ECMS is also important. For empirical verification, it is necessary to compare NLP-based ECMS with other online optimization approaches, such as PSO [21,22,23] or NN-based method [27,28,29,30,31]. We are conducting a follow-up study for empirical verification of NLP-based ECMS. In conclusion, follow-up studies are required to apply NLP-based ECMS to embedded systems. NLP-based ECMS includes NLP, which is a challenge when applied to embedded systems. To apply NLP-based ECMS to embedded systems, two representative follow-up studies are needed.
The first is to improve the computational speed and computational stability of NLP. To improve the computational performance of NLP, it is necessary to simplify OCP. Simplifying OCP means softening constraints on road information (speed profile, grade profile, etc.) or vehicle models (degree of surface fitting polynomial for engine/motor map). It is important to determine the appropriate prediction/control horizon through this process. Improving computational stability means considering the infeasible region of the NLP. NLP can derive infeasible solutions because of the exceptional driving situation, like a cut-in.
The second is to develop a NLP solver suitable for embedded systems. Various studies have been actively conducted on the research subject [43,44,45], such as SQP or interior-point method. NLP requires a very large amount of computation because it uses numerical and recursive methods. Therefore, it is necessary to optimize the amount of computation through solution space analysis. The solution space is determined by an operation characteristic or a driving situation of the HEV.
The proposed NLP-based ECMS does not consider clutch and transmission effects. Because the simulation results show that clutch and transmission effects degrade the performance of the proposed method, the influence of clutch and transmission must be considered when defining the OCP in NLP. A MINLP method may be used because the clutch and transmission involve mathematically discrete variables.
In this paper, previously supplied driving conditions were assumed. However, it is impossible to accurately anticipate actual driving conditions. An EMS must therefore reflect changing driving conditions, which depend on traffic conditions.

Author Contributions

Formal analysis, J.L.; project administration, J.L. and H.L.; software, J.L.; validation, J.L.; writing—original draft preparation, J.L.; writing—review and editing, J.L. and H.L.; supervision, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by the Autonomous Driving Development Innovation Program (20014476, Development of Mixed Reality-based Autonomous Driving Parts and System Evaluation Technology) funded by the Ministry of Trade, Industry & Energy (MOTIE, Korea).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Y.; Huang, Z. Optimization-based energy management strategy for a 48-V mild parallel hybrid electric power system. J. Energy Resour. Technol. 2020, 142, 052002. [Google Scholar]
  2. Larsson, V.; Johannesson, L.; Egardt, B.; Lassson, A. Benefit of route recognition in energy management of plug-in hybrid electric vehicles. In Proceedings of the 2012 American Control Conference (ACC), Montréal, QC, Canada, 27–29 June 2012; pp. 1314–1320. [Google Scholar]
  3. Overington, S.; Rajakaruna, S. Review of PHEV and HEV operation and control research for future direction. In Proceedings of the 2012 3rd IEEE International Symposium on Power Electronics for Distributed Generation Systems (PEDG), Aalborg, Denmark, 25–28 June 2012; pp. 385–392. [Google Scholar]
  4. Lee, J.; Park, B. Development and evaluation of a cooperative vehicle intersection control algorithm under the connected vehicles environment. IEEE Trans. Intell. Transp. Syst. 2012, 13, 81–90. [Google Scholar] [CrossRef]
  5. Zeng, Y.; Cai, Y.; Kou, G.; Gao, W.; Qin, D. Energy management for plug-in hybrid electric vehicle based on adaptive simplified-ECMS. Sustainability 2018, 10, 2060. [Google Scholar] [CrossRef] [Green Version]
  6. Hou, C.; Ouyang, M.; Xu, L.; Wang, H. Approximate Pontryagin’s minimum principle applied to the energy management of plug-in hybrid electric vehicles. Appl. Energy 2014, 115, 174–189. [Google Scholar] [CrossRef]
  7. Mechichi, O.; Trigui, R.; El Amraoui, L. Adaptive λ-Control Strategy for Plug-In HEV Energy Management Using Fast Initial Multiplier Estimate. Appl. Sci. 2022, 12, 10543. [Google Scholar] [CrossRef]
  8. Tormos, B.; Pla, B.; Bares, P.; Pinto, D. Energy Management of Hybrid Electric Urban Bus by Off-Line Dynamic Programming Optimization and One-Step Look-Ahead Rollout. Appl. Sci. 2022, 12, 4474. [Google Scholar] [CrossRef]
  9. Borhan, H.; Vahidi, A.; Phillips, A.M.; Kuang, M.L.; Kolmanovsky, I.V.; Di Cairano, S. MPC-based energy management of a power-split hybrid electric vehicle. IEEE Trans. Control Syst. Technol. 2011, 20, 593–603. [Google Scholar] [CrossRef]
  10. Vu, T.M.; Moezzi, R.; Cyrus, J.; Hlava, J.; Petru, M. Parallel Hybrid Electric Vehicle Modelling and Model Predictive Control. Appl. Sci. 2021, 11, 10668. [Google Scholar] [CrossRef]
  11. Zhang, B.; Xu, F.; Shen, T. A Real-Time Energy Management Strategy for Parallel HEVs with MPC. In Proceedings of the 2019 IEEE Vehicle Power and Propulsion Conference (VPPC), Hanoi, Vietnam, 14–17 October 2019; pp. 1–5. [Google Scholar]
  12. Shivappriya, S.N.; Karthikeyan, S.; Prabu, S.; Pérez de Prado, R.; Parameshachari, B.D. A Modified ABC-SQP-Based Combined Approach for the Optimization of a Parallel Hybrid Electric Vehicle. Energies 2020, 13, 4529. [Google Scholar] [CrossRef]
  13. Yang, L.; Wang, W.; Yang, C.; Du, X.; Zhang, W. Online mixed-integer optimal energy management strategy for connected hybrid electric vehicles. J. Clean. Prod. 2022, 374, 133908. [Google Scholar] [CrossRef]
  14. Liu, H.; Liu, R.; Xu, R.; Han, L.; Ruan, S. Hierarchical energy management strategy considering switching schedule for a dual-mode hybrid electric vehicle. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022, 236, 938–949. [Google Scholar] [CrossRef]
  15. Onori, S.; Serrao, L.; Rizzoni, G. Hybrid Electric Vehicles: Energy Management Strategies; Springer: Cham, Switzerland, 2016. [Google Scholar]
  16. Simmons, K.; Guezennec, Y.; Onori, S. Modeling and energy management control design for a fuel cell hybrid passenger bus. J. Power Sources 2014, 246, 736–746. [Google Scholar] [CrossRef]
  17. Zhang, P.; Wu, X.; Du, C.; Xu, H.; Wang, H. Adaptive equivalent consumption minimization strategy for hybrid heavy-duty truck based on driving condition recognition and parameter optimization. Energies 2020, 13, 5407. [Google Scholar] [CrossRef]
  18. Lin, W.; Zhao, H.; Zhang, B.; Wang, Y.; Xiao, Y.; Xu, K.; Zhao, R. Predictive Energy Management Strategy for Range-Extended Electric Vehicles Based on ITS Information and Start–Stop Optimization with Vehicle Velocity Forecast. Energies 2022, 15, 7774. [Google Scholar] [CrossRef]
  19. Zanelli, A.; Servetto, E.; De Araujo, P.; Vankayala, S.N.; Vondrak, A. Numerical Assessment of Auto-Adaptive Energy Management Strategies Based on SOC Feedback, Driving Pattern Recognition and Prediction Techniques. Energies 2022, 15, 3896. [Google Scholar] [CrossRef]
  20. Kommuri, N.K.; McGordon, A.; Allen, A.; Truong, D.Q. A Novel Adaptive Equivalence Fuel Consumption Minimisation Strategy for a Hybrid Electric Two-Wheeler. Energies 2022, 15, 3192. [Google Scholar] [CrossRef]
  21. Li, T.; Cui, W.; Cui, N. Soft Actor-Critic Algorithm-Based Energy Management Strategy for Plug-In Hybrid Electric Vehicle. World Electr. Veh. J. 2022, 13, 193. [Google Scholar] [CrossRef]
  22. Wang, Z.; Jiao, X. Optimization of the powertrain and energy management control parameters of a hybrid hydraulic vehicle based on improved multi-objective particle swarm optimization. Eng. Optim. 2021, 53, 1835–1854. [Google Scholar] [CrossRef]
  23. Oubelaid, A.; Taib, N.; Nikolovski, S.; Alharbi, T.E.; Rekioua, T.; Flah, A.; Ghoneim, S.S. Intelligent speed control and performance investigation of a vector controlled electric vehicle considering driving cycles. Electronics 2022, 11, 1925. [Google Scholar] [CrossRef]
  24. Pérez, W.; Tulpule, P.; Midlam-Mohler, S.; Rizzoni, G. Data-Driven Adaptive Equivalent Consumption Minimization Strategy for Hybrid Electric and Connected Vehicles. Appl. Sci. 2022, 12, 2705. [Google Scholar] [CrossRef]
  25. Wang, W.; Cai, Z.; Liu, S. Design of Real-Time Control based on DP and ECMS for PHEVs. Math. Probl. Eng. 2021, 2021, 6667614. [Google Scholar] [CrossRef]
  26. Han, J.; Park, Y.; Kum, D. Optimal adaptation of equivalent factor of equivalent consumption minimization strategy for fuel cell hybrid electric vehicles under active state inequality constraints. J. Power Sources 2014, 267, 491–502. [Google Scholar] [CrossRef]
  27. Sun, X.; Xue, M.; Cai, Y.; Tian, X.; Jin, Z.; Chen, L. Adaptive ECMS Based on EF Optimization by Model Predictive Control for Plug-in Hybrid Electric Buses. IEEE Trans. Transp. Electrif. 2022. early access. [Google Scholar] [CrossRef]
  28. Chen, Z.; Liu, Y.; Zhang, Y.; Lei, Z.; Chen, Z.; Li, G. A neural network-based ECMS for optimized energy management of plug-in hybrid electric vehicles. Energy 2022, 243, 122727. [Google Scholar] [CrossRef]
  29. Han, L.; Jiao, X.; Zhang, Z. Recurrent neural network-based adaptive energy management control strategy of plug-in hybrid electric vehicles considering battery aging. Energies 2020, 13, 202. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, B.; Meng, Q.; Wu, J.; Ni, Y. An Adaptive Adjustment Method of Equivalent Factor Considering Speed Predict Information. World Electr. Veh. J. 2021, 12, 211. [Google Scholar] [CrossRef]
  31. Shanmuganathan, J.; Victoire, A.A.; Balraj, G.; Victoire, A. Deep Learning LSTM Recurrent Neural Network Model for Prediction of Electric Vehicle Charging Demand. Sustainability 2022, 14, 10207. [Google Scholar] [CrossRef]
  32. Musardo, C.; Rizzoni, G.; Guezennec, Y.; Staccia, B. A-ECMS: An adaptive algorithm for hybrid electric vehicle energy management. Eur. J. Control 2005, 11, 509–524. [Google Scholar] [CrossRef]
  33. Sivertsson, M. Adaptive control using map-based ECMS for a PHEV. IFAC Proc. Vol. 2012, 45, 357–362. [Google Scholar] [CrossRef] [Green Version]
  34. Onori, S.; Serrao, L.; Rizzoni, G. Adaptive equivalent consumption minimization strategy for hybrid electric vehicles. In Proceedings of the Dynamic Systems and Control Conference, Cambridge, MA, USA, 12–15 September 2010; pp. 499–505. [Google Scholar]
  35. Gu, B.; Rizzoni, G. An adaptive algorithm for hybrid electric vehicle energy management based on driving pattern recognition. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, Chicago, IL, USA, 5–10 November 2006; pp. 249–258. [Google Scholar]
  36. Kessels, J.T.; Koot, M.W.; Van Den Bosch, P.P.; Kok, D.B. Online energy management for hybrid electric vehicles. IEEE Trans. Veh. Technol. 2008, 57, 3428–3440. [Google Scholar] [CrossRef]
  37. AUTONOMIE. Available online: https://vms.taps.anl.gov/ (accessed on 26 October 2022).
  38. Balazs, A.; Morra, E.; Pischinger, S. Optimization of electrified powertrains for city cars. SAE Int. J. Altern. Powertrains 2012, 1, 381–394. [Google Scholar] [CrossRef]
  39. CasADi. Available online: https://web.casadi.org/ (accessed on 26 October 2022).
  40. IPOPT. Available online: https://coin-or.github.io/Ipopt/ (accessed on 26 October 2022).
  41. Biegler, L.T.; Zavala, V.M. Large-scale nonlinear programming using IPOPT: An integrating framework for enterprise-wide dynamic optimization. Comput. Chem. Eng. 2009, 33, 575–582. [Google Scholar] [CrossRef]
  42. Pirnay, H.; López-Negrete, R.; Biegler, L.T. Optimal sensitivity based on IPOPT. Math. Program. Comput. 2012, 4, 307–331. [Google Scholar] [CrossRef]
  43. Jabr, R.A.; Coonick, A.H.; Cory, B.J. A primal-dual interior point method for optimal power flow dispatching. IEEE Trans. Power Syst. 2002, 17, 654–662. [Google Scholar] [CrossRef]
  44. Debert, M.; Yhamaillard, G.; Ketfi-herifellicaud, G.A. Predictive energy management for hybrid electric vehicles-Prediction horizon and battery capacity sensitivity. IFAC Proc. Vol. 2010, 43, 270–275. [Google Scholar] [CrossRef] [Green Version]
  45. Dueri, D.; Açıkmeşe, B.; Scharf, D.P.; Harris, M.W. Customized real-time interior-point methods for onboard powered-descent guidance. J. Guid. Control Dyn. 2017, 40, 197–212. [Google Scholar] [CrossRef]
Figure 1. Structure of an HEV with a P2 powertrain.
Figure 1. Structure of an HEV with a P2 powertrain.
Applsci 12 12724 g001
Figure 2. Brake-specific fuel consumption (BSFC) map.
Figure 2. Brake-specific fuel consumption (BSFC) map.
Applsci 12 12724 g002
Figure 3. Electric motor power map.
Figure 3. Electric motor power map.
Applsci 12 12724 g003
Figure 4. Conceptual diagram of the NLP-based ECMS algorithm.
Figure 4. Conceptual diagram of the NLP-based ECMS algorithm.
Applsci 12 12724 g004
Figure 5. Gear shift rule. (a) Upshift map. (b) Downshift map.
Figure 5. Gear shift rule. (a) Upshift map. (b) Downshift map.
Applsci 12 12724 g005
Figure 6. Surface fitting results for a fuel-consumption rate map. (a) surface fitting result with various degree. (b) Fitting error according to the degree of polynomial.
Figure 6. Surface fitting results for a fuel-consumption rate map. (a) surface fitting result with various degree. (b) Fitting error according to the degree of polynomial.
Applsci 12 12724 g006
Figure 7. Conceptual diagram of A-ECMS in NLP-ECMS.
Figure 7. Conceptual diagram of A-ECMS in NLP-ECMS.
Applsci 12 12724 g007
Figure 8. Autonomie simulation environment.
Figure 8. Autonomie simulation environment.
Applsci 12 12724 g008
Figure 9. Comparison of DP results (1): state of charge.
Figure 9. Comparison of DP results (1): state of charge.
Applsci 12 12724 g009
Figure 10. Comparison of DP results (2): engine torque.
Figure 10. Comparison of DP results (2): engine torque.
Applsci 12 12724 g010
Figure 11. Comparison of NLP results to DP results (2): reference SOC trajectory.
Figure 11. Comparison of NLP results to DP results (2): reference SOC trajectory.
Applsci 12 12724 g011
Figure 12. Result of NLP: (a) ICE torque; (b) EM torque. (blue = state/red = input/black = maximum/minimum torque constraints.
Figure 12. Result of NLP: (a) ICE torque; (b) EM torque. (blue = state/red = input/black = maximum/minimum torque constraints.
Applsci 12 12724 g012
Figure 13. Autonomie (P2 HEV) simulation result with an EMS (NLP-based ECMS): (a) vehicle speed; (b) state of charge; (c) ICE torque; (d) EM torque.
Figure 13. Autonomie (P2 HEV) simulation result with an EMS (NLP-based ECMS): (a) vehicle speed; (b) state of charge; (c) ICE torque; (d) EM torque.
Applsci 12 12724 g013
Table 1. Specification of HEV.
Table 1. Specification of HEV.
ComponentParameterValue
Engine (4-cylinder)Displacement
Maximum power
2.2 [L]
85.0 [kW]
Motor (PMSM)Maximum power29.2 [kW]
Battery (Lithium-ion)Capacity
Nominal voltage
7.03 [Ah]
324 [V]
Transmission (5-speed automatic)Gear ratio2.563/1.552/1.022/0.727/0.52
Final driveGear ratio4.438
WheelWheel radius0.30115 [m]
ChassisMass1680 [kg]
Table 2. Settings for comparison groups (DP).
Table 2. Settings for comparison groups (DP).
CaseInput GridState GridStep TimeCalculation Time
0 (reference)2016010.1 [s]2481.9 [s]
1201610.1 [s]210.7 [s]
2216010.1 [s]210.1 [s]
321610.1 [s]89.2 [s]
4201611.0 [s]20.4 [s]
5216011.0 [s]20.8 [s]
621611.0 [s]8.8 [s]
Table 3. Comparison of NLP results and DP results (1): calculation time and sum of error.
Table 3. Comparison of NLP results and DP results (1): calculation time and sum of error.
CaseInput Grid State GridStep TimeCalculation TimeSum of Error
0 (reference)2016010.1 [s]2481.9 [s]0
2216010.1 [s]210.1 [s]45.187
5216011.0 [s]20.8 [s]508.71
621611.0 [s]8.8 [s]717.21
NLP---2.6 [s]156.77
Table 4. Analysis of Autonomie simulation results: Fuel Improvement and Simulation Time.
Table 4. Analysis of Autonomie simulation results: Fuel Improvement and Simulation Time.
CaseFuel
Consumption
Fuel
Improvement
Final
SOC
Simulation
Time
DP0.512 [kg]
Ruel-based0.595 [kg]0.0 [%]0.61 [-]73.79 [s]
00.527 [kg]11.4 [%]0.61 [-]2553.18 [s]
20.527 [kg]11.4 [%]0.61 [-]281.49 [s]
50.533 [kg]10.4 [%]0.61 [-]90.39 [s]
60.536 [kg]9.9 [%]0.61 [-]81.00 [s]
NLP0.531 [kg]10.8 [%]0.61 [-]74.18 [s]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, J.; Lee, H. A New HEV Power Distribution Algorithm Using Nonlinear Programming. Appl. Sci. 2022, 12, 12724. https://doi.org/10.3390/app122412724

AMA Style

Lee J, Lee H. A New HEV Power Distribution Algorithm Using Nonlinear Programming. Applied Sciences. 2022; 12(24):12724. https://doi.org/10.3390/app122412724

Chicago/Turabian Style

Lee, Jooin, and Hyeongcheol Lee. 2022. "A New HEV Power Distribution Algorithm Using Nonlinear Programming" Applied Sciences 12, no. 24: 12724. https://doi.org/10.3390/app122412724

APA Style

Lee, J., & Lee, H. (2022). A New HEV Power Distribution Algorithm Using Nonlinear Programming. Applied Sciences, 12(24), 12724. https://doi.org/10.3390/app122412724

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

Article Metrics

Back to TopTop