# Real-Time Control Strategy of Elman Neural Network for the Parallel Hybrid Electric Vehicle

Authors: Ruijun Liu; Dapai Shi; Chao Ma
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

Through researching the instantaneous control strategy and Elman neural network, the paper established equivalent fuel consumption functions under the charging and discharging conditions of power batteries, deduced the optimal control objective function of instantaneous equivalent consumption, established the instantaneous optimal control model, and designs the Elman neural network controller. Based on the ADVISOR 2002 platform, the instantaneous optimal control strategy and the Elman neural network control strategy were simulated on a parallel HEV. The simulation results were analyzed in the end. The contribution of the paper is that the trained Elman neural network control strategy can reduce the simulation time by 96% and improve the real-time performance of energy control, which also ensures the good performance of power and fuel economy.


---

## Body

## 1. Introduction

Under the dual pressure of environmental pollution and energy crisis, hybrid vehicles have advantages of both conventional vehicles and electric vehicles, which have characteristics of energy conservation, environmental protection, diverse shapes, and strong implementation. Hybrid vehicles have become an effective way to solve the problem of energy crisis and environmental protection and also have been one of the most perspective vehicle models.According to different connective ways of power system, hybrid electric vehicle (HEV) can be mainly divided into four styles: series, parallel, series-parallel, and complex. The dynamic structure diagram is shown in Figure1.Figure 1
Classification of hybrid electric vehicles.Engine output energy of series HEV is transformed two times and the efficiency of the motor and generator is relatively low, so series HEV loses more energy and leads to lower efficiency than vehicles of internal combustion engine. Parallel HEV (PHEV) is equipped with series and parallel power systems, and their structures and control systems are more complex and have higher cost. Complex HEV structures and control systems are most complex and have highest cost. However, parallel hybrid power system can adapt to various road conditions and is widely used by enterprises [1, 2].As the core of multiple energy control system, the energy control strategy determines performances of PHEV. Based on vehicle’s torque, energy control strategies of PHEV are mainly divided into four types [3, 4]: static logic threshold energy control strategy, instantaneous optimal energy control strategy, global optimal energy control strategy, and neural network energy control strategy. The static logic threshold control strategy cannot guarantee the optimal fuel economy of PHEV, does not adapt to dynamic conditions, and cannot make the whole system to achieve maximum efficiency. Besides, its threshold parameters are set by engineering experience [5, 6]. The global optimal control strategy can achieve global optimal fuel consumption of HEV. Its defects are complex algorithm, large amount of calculation, and knowing the whole condition in advance [7, 8]. Neural network energy control strategy can adapt to diverse conditions with good robustness and obtain global fuel consumption optimum by engineering experience [9, 10]. Instantaneous optimal energy control strategy can realize minimum equivalent fuel consumption of PHEV in each control cycle which is widely used to distribute PHEV energy [11].The basic principle of instantaneous optimal control strategy is based on the model of the optimal curve of engines; the object function of the whole power system was optimized on the specific operating points of parallel HEV. On the basis of the instantaneous optimal operating points, it makes power of variable states redistributed and make the loss of energy minimized in the energy flow process at any time (see Figure2). Instantaneous optimal control strategy has good fuel economy at any time and bad real-time performance. Its real-time performance is influenced by these factors which are the accuracy of various components battery ages and engine and motor characteristics [12–14]. So it is difficult to improve the real-time performance of instantaneous optimal control strategy by changing these factors.Figure 2
Energy flow diagram of parallel hybrid electric vehicles.The hybrid vehicles possess good power performance and fuel economy and obtain rapid allocation energy by finding a new energy control strategy. Elman neural network is a feedback neural network and has a very strong computing ability and stability [15]. The instantaneous rules of the instantaneous optimal control strategy are used to train Elman neural network, establish Elman neural network controller, and improve the real-time performance of energy control [14–16].Based on the research of the instantaneous optimal control strategy, the strategy possesses good fuel economy and makes energy distributed reasonably. However, its real-time performance is poor. In order to solve bad real-time defects of instantaneous control strategy, instantaneous optimal control rules are used to train the Elman neural network control strategy and improve the real-time performance of the trained Elman energy control strategy on the premise that it can guarantee advantages of the instantaneous optimal control strategy. The results show that the trained Elman neural network control strategy can replace the instantaneous optimal control strategy, optimize power distribution, and make the simulation time reduced by 60%.


## 2. Research Energy Control Strategy

### 2.1. Instantaneous Optimal Control Strategy

Instantaneous optimal control strategy is defined as follows. In order to achieve the minimum fuel consumption of HEV, the optimal output power of the engine and electric motor is calculated in each control cycle of hybrid power system. Working conditions of HEV and calculation expressions of the equivalent fuel consumption are different in every time. So an optimal objective function should be established [17]. Working condition of hybrid system is divided into two cases of power battery charging and discharging. Then, objective functions of the instantaneous equivalent minimum fuel consumption were established on two working conditions.Here, the full line represents the circulation and transformation of fuel chemical energy in the hybrid power system. The dotted line represents electric current circulation and transformation in the hybrid power system.


#### 2.1.1. Calculate the Equivalent Fuel Consumption of Battery Discharging

When the power battery takes part in driving hybrid cars, its SOC value will reduce and deviate from the target of SOC value. In order to compensate for the used electricity and restore SOC value of power batteries, the engine drives the motor to charge power batteries in the future time [18]. The charging time is divided into n control cycles, and the motor power of each cycle is P
mc
_
chg
_
i  (
i
=
1,2
,
3
…
,
n
) [5].The relationship between the motor power (P
mc) with driving vehicle and the motor power (P
mc
_
chg
_
i) with power batteries charging is
(1)
P
mc
=
η
-
mc
2
η
-
dischg
η
-
chg
∑
i
=
1
n
P
mc
_
chg
_
i
,
where η
-
mc is the average efficiency of motor; η
-
chg is the average efficiency of power battery charging; η
-
dischg is the average efficiency of power battery discharging.When the motor drives the vehicle, the energy consumption of power batteries can be converted into the engine fuel consumption. The equivalent fuel consumption rate of the motor is(2)
b
mc
_
eq
=
∑
i
=
1
n
(
P
fc
Y
_
i
b
Y
_
i
-
P
fc
N
_
i
b
N
_
i
)
η
¯
mc
2
η
¯
dischg
η
¯
chg
∑
i
=
1
n
P
mc
_
chg
_
i
,
where P
fc
N
_
i, when the engine does not charge power batteries, is the engine power of the i control period; b
N
_
i, when the engine does not charge power batteries, is the fuel consumption rate of the engine; P
fc
Y
_
i, when the engine charges power batteries, is the engine power of the i control period; b
Y
_
i, when the engine charges power batteries, is the fuel consumption rate of the engine.Letb
-
chg be the average fuel consumption rate when the engine charges power batteries:
(3)
b
¯
chg
=
∑
i
=
1
n
P
fc
Y
_
i
b
Y
_
i
-
P
fc
N
_
i
b
N
_
i
P
mc
_
chg
_
i
.Merge (3) and (2):
(4)
b
mc
_
eq
=
b
¯
chg
η
¯
mc
2
η
¯
dischg
η
¯
chg
.When the motor drives the vehicle after a period ofΔ
t, equivalent instantaneous fuel consumption of the motor is
(5)
m
mc
_
eq
=
b
¯
chg
P
mc
η
¯
mc
2
η
¯
dischg
η
¯
chg
.


#### 2.1.2. Calculate the Fuel Consumption of Batteries Charging

When the power battery is charged by the engine, its SOC value will rise and even exceed the target of the SOC value. In order to maintain SOC values, power battery energy will be consumed in future [18]. Discharging time is divided into n control cycles, and the motor power of the each control cycle is P
mc
2
_
i  (
i
=
1,2
,
3
…
,
n
) [17].In a certain period of discharging time, the relationship between motor power (P
mc
2
_
i) with driving vehicle and motor power (P
mc
_
chg
_
2) with charging power batteries is
(6)
∑
i
=
1
n
P
mc
2
_
i
=
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
chg
2
η
¯
dischg
2
,
where P
mc
_
chg
2 is the motor power when the power battery is charged; η
-
mc
2 is the average efficiency of motor; η
-
chg
2 is the average efficiency of the power battery charging; η
-
dischg
2 is the average efficiency of the power battery discharging.When the motor drives the vehicle, the relationship between the motor power battery energy consumption and the fuel consumption rate is(7)
b
mc
_
eq
2
=
P
fc
Y
2
b
Y
2
-
P
fc
N
2
b
N
2
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
chg
2
η
¯
dischg
2
,
where P
fc
Y
2, when the engine charges the power battery, is the power of the engine; b
Y
2, when the engine charges the power battery, is the fuel consumption rate of the engine; P
fc
N
2, when the engine does not charge the power battery, is the power of the engine; b
N
2, when the engine does not charge the power battery, is the fuel consumption rate of the engine.Letb
chg
2 be the fuel consumption rate when the engine charges the power battery:
(8)
b
chg
2
=
P
fc
Y
2
b
Y
2
-
P
fc
N
2
b
N
2
P
mc
_
chg
2
.Simplify the (8) formula:
(9)
b
mc
_
eq
2
=
b
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
.When the motor charges power batteries after a period ofΔ
t, the objective function of equivalent instantaneous fuel consumption of the motor is
(10)
m
mc
_
eq
2
=
b
¯
chg
2
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
.


#### 2.1.3. Deduce the Objective Function of the Instantaneous Optimal Control Strategy

Set two new variables:(11)
f
eq
_
dischg
=
b
¯
chg
η
¯
mc
2
η
¯
dischg
η
¯
chg
,
f
eq
_
chg
=
b
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
,
where b
chg
2 is the fuel consumption rate when the engine charges power batteries at the present moment; b
-
chg is the average fuel consumption rate when the engine charges power batteries in the future time.The instantaneous control objective function of the lowest fuel consumption is(12)
M
=
∑
Min
⁡
{
[
T
mc
(
t
)
,
m
mc
_
eq
ω
mc
(
t
)
]
Δ
t
m
fc
[
T
fc
(
t
)
,
ω
(
t
)
]
Δ
t
f
i
f
f
f
f
f
f
f
f
+
m
mc
_
eq
[
T
mc
(
t
)
,
ω
mc
(
t
)
]
Δ
t
}
,
where ω
fc
(
t
) is engine speed; ω
mc
(
t
) is motor speed; T
fc is output torque of the engine; T
mc is output torque of the motor. Consider
(13)
m
mc
_
eq
[
P
mc
(
t
)
]
=
{
f
eq
_
dischg
P
mc
T
mc
>
0
(
power
batteries
discharge
)
f
eq
_
chg
P
mc
_
chg
2
T
mc
<
0
(
power
batteries
charge
)
.P
mc and P
mc
_
chg
2 are both motor power, so they can be unified as P
mc. By calculating, the improved instantaneous control objective function of the minimum fuel consumption is
(14)
M
=
∑
Min
⁡
{
[
λ
f
eq
_
dischg
+
(
1
-
λ
)
f
eq
_
chg
]
m
fc
[
P
fc
(
t
)
]
Δ
t
f
f
f
f
f
f
f
i
i
f
+
[
λ
f
eq
_
dischg
+
(
1
-
λ
)
f
eq
_
chg
]
P
mc
Δ
t
}
,
where λ
=
(
1
+
sign
⁡
(
T
mc
)
)
/
2.


#### 2.1.4. Improve the Objective Function of the Instantaneous Optimal Control Strategy

SOC value change of batteries and braking energy recovery both have a certain effect on energy control. The optimal function of the working point needs to be improved.(
1) Revise the SOC Value Function of Power Batteries. When power batteries work, their SOC value is maintained at the high efficient range by the reset function, in order to reduce the loss energy in the process of charging and discharging power batteries and make hybrid system keep better performances. The working principle of SOC reset function is as follows. When the SOC value is more than the target region, the hybrid power system will give priority to consuming power battery energy. It does not stop until SOC value decreases to the target region under the effect of the reset function. When the battery SOC value is lower than the target region, the hybrid system will give priority to consuming fuel energy to drive the vehicle and recover the value of SOC. It does not stop until SOC value returns to the target region under the effect of the reset function.SetK
SOC be the variable in the reset function, and the value table between K
SOC and SOC value is shown in Table 1.Table 1
The value table betweenK
soc and SOC.
SOC
0.11
0.14
0.2
0.25
0.3
0.35
K
soc
7
2
1.05
1.05
1.03
1.04
SOC
0.4
0.45
0.5
0.55
0.59
0.64
K
soc
1.02
1.02
1
0.99
0.98
0.96
SOC
0.7
0.74
0.81
0.85
0.9
K
soc
0.96
0.94
0. 94
0.8
0.2Based on the Matlab platform and Table1 data, the fitting curve between K
SOC and SOC value is constructed by using fitting curve toolbox, as shown in Figure 3.Figure 3
SOC andK
SOC relationship curve.Polynomial function of the fitting curve is(15)
K
soc
=
-
8866
x
7
+
33460
x
6
-
52310
x
5
+
43730
x
4
-
21010
x
3
+
5766
x
2
-
833.5
x
+
49.77
,
where x: is SOC value.Considering the influence of power battery SOC, the formula of instantaneous equivalent fuel consumption of HEV can be expressed asK
soc
m
mc
_
eq.(
2
) Revise the Objective Function of Braking Recovery Energy. The adopted method which revises the equivalent fuel consumption function of power battery energy is as follows: the average braking power is calculated at a period of time before the current moment. The power is used as the standard of fuel consumption correction of the braking recovery power in the next moment [19–22].The statistical time range is divided inton  (
n
>
0
) cycles. In each cycle, let braking power be a fixed value. Therefore, the average braking power of the whole time can be expressed as
(16)
P
-
braking
=
∑
i
=
1
n
P
braking
_
i
n
.When power batteries discharge, the objective function of instantaneous equivalent fuel consumption is(17)
m
mc
_
eq
=
b
-
chg
(
P
mc
+
P
-
braking
)
η
-
mc
2
η
-
dischg
η
-
chg
.In summary, taking the influence of power battery SOC and brake energy recovery into consideration, the final objective function of instantaneous optimal control strategy is(18)
M
=
∑
Min
⁡
{
[
P
mc
(
t
)
+
P
-
braking
(
t
)
]
m
fc
[
P
fc
(
t
)
]
Δ
t
fffffffffffffff
+
K
soc
m
mc
_
eq
[
P
mc
(
t
)
+
P
-
braking
(
t
)
]
Δ
t
}
.


### 2.2. Elman Neural Network Control Strategy

#### 2.2.1. The Structure of Elman Neural Network

Elman neural network is put forward by Jeffrey L. Elman in 1990 and is a typical local recessionary grid, as shown in Figure4 [23].Figure 4
Elman neural network model.P is the input of the neural network and its size is R
×
1; b
1 is the neuronal threshold vector of the feedback layer and its size is S
1
×
1; I
W
1 is the connective weight vector of the neurons and the input vector in the input layer and its size is S
1
×
R; n
1 is the middle operational result of the neurons in the feedback layer, namely, weighted sum of the connective weight vector and the threshold vector, and its size is S
1
×
1; a
1 is the output vector of the feedback layer’s neuron in the K iteration and its size is S
1
×
1; D is feedback node. It is the same way that  b
2, L
W
2, n
2, and a
2 are related to parameters of the output layer.The inputs of structure diagram of Elman neural network is the required torque, speed, and SOC value of power battery and its output is motor torque. The structure diagram of Elman neural network is shown in Figure5 [24].Figure 5
Structure diagram of Elman neural network.The structure diagram of Elman neural network contains input layer, hidden layer, undertaken layer, and output layer. Let the input vector of the input layer be three-dimensional vectoru; the output vector of the output layer is one-dimensional vector y; the output vector of the hidden layer is n-dimensional vector x; the output vector of the undertaken layer is n-dimensional vector r; w
1, w
2, and w
3 are respective connective weights of the hidden layer to the output layer, the input layer to the hidden layer, and the undertaken layer to the hidden layer; g
(
·
) is the driving function of the output neurons; f
(
·
) is the driving function of hidden layer; h
(
·
) is the driving function of undertaken layer; net(
·
) is the net input driving function of a certain layer; A shows the input layer; B shows the undertaken layer; K shows the iterative sequence.Define two functions:(19)
v
i
(
k
)
=
{
u
n
(
k
)
,
if
i
∈
A
,
r
n
(
k
)
,
if
i
∈
B
,
w
i
(
k
)
=
{
w
2
,
if
i
∈
A
,
w
3
,
if
i
∈
B
.The input and output functions of theN neuron of the hidden layer are
(20)
net
n
(
k
+
1
)
=
∑
i
∈
A
∪
B
w
i
(
k
)
v
i
(
k
)
,
x
n
(
n
+
1
)
=
f
(
net
n
(
k
+
1
)
)
.The input and output functions of theN neuron of the hidden layer are
(21)
net
n
(
k
)
=
∑
i
∈
A
∪
B
w
i
(
k
-
1
)
v
i
(
k
-
1
)
,
r
(
k
)
=
h
(
net
n
(
k
)
)
.The input and output functions of the output layer’s neuron are(22)
net
n
(
k
+
1
)
=
∑
i
∈
A
∪
B
w
i
(
k
+
1
)
x
n
(
k
+
1
)
,
y
(
k
+
1
)
=
g
(
net
(
k
+
1
)
)
.


#### 2.2.2. Select the Parameters of Elman Neural Network

The neuron number is determined by following formula [25]:
(23)
k
=
m
+
n
+
β
,
where m is the number of the input vector; n is the neuron number of the output vector; β is a constant, (1~10).The excitation function of Elman neural network of the feedback layer selects the Tansig function [26]:
(24)
tansig
(
x
)
=
2
1
+
e
-
2
x
-
1
.


#### 2.2.3. Learning and Training Mechanism of Elman Neural Network

Elman neural network is trained by Levenberg-Marquardt algorithm. The error index function of Levenberg-Marquardt arithmetic is(25)
E
(
w
)
=
1
2
∑
i
=
1
p
∥
y
i
-
y
i
′
∥
2
=
1
2
∑
i
=
1
p
e
i
2
(
w
)
,
where p is the sample number; e
i is the systemic error; y
i
′ is the actual output of the network.The formula of the adjusting weight is(26)
w
k
+
1
=
w
k
+
Δ
w
.The computing formula of the increment weight is(27)
Δ
w
=
[
J
T
(
w
)
J
(
w
)
+
u
I
]
-
1
J
T
(
w
)
e
(
w
)
,
where u is learning rate; I is the unit matrix; J
(
w
) is the Jacobian matrix. Consider
(28)
J
(
w
)
=
[
∂
e
1
(
w
)
∂
w
1
∂
e
1
(
w
)
∂
w
2
⋯
∂
e
1
(
w
)
∂
w
n
∂
e
2
(
w
)
∂
w
1
∂
e
2
(
w
)
∂
w
1
⋯
∂
e
2
(
w
)
∂
w
n
⋮
⋮
⋱
⋮
∂
e
n
(
w
)
∂
w
1
∂
e
n
(
w
)
∂
w
2
⋯
∂
e
n
(
w
)
∂
w
n
]
.Neural network weight is adjusted by the Levenberg-Marquardt algorithm, and the adjusting process is shown in Figure6.Figure 6
Process map of the adjusting network weight.


## 2.1. Instantaneous Optimal Control Strategy

Instantaneous optimal control strategy is defined as follows. In order to achieve the minimum fuel consumption of HEV, the optimal output power of the engine and electric motor is calculated in each control cycle of hybrid power system. Working conditions of HEV and calculation expressions of the equivalent fuel consumption are different in every time. So an optimal objective function should be established [17]. Working condition of hybrid system is divided into two cases of power battery charging and discharging. Then, objective functions of the instantaneous equivalent minimum fuel consumption were established on two working conditions.Here, the full line represents the circulation and transformation of fuel chemical energy in the hybrid power system. The dotted line represents electric current circulation and transformation in the hybrid power system.


### 2.1.1. Calculate the Equivalent Fuel Consumption of Battery Discharging

When the power battery takes part in driving hybrid cars, its SOC value will reduce and deviate from the target of SOC value. In order to compensate for the used electricity and restore SOC value of power batteries, the engine drives the motor to charge power batteries in the future time [18]. The charging time is divided into n control cycles, and the motor power of each cycle is P
mc
_
chg
_
i  (
i
=
1,2
,
3
…
,
n
) [5].The relationship between the motor power (P
mc) with driving vehicle and the motor power (P
mc
_
chg
_
i) with power batteries charging is
(1)
P
mc
=
η
-
mc
2
η
-
dischg
η
-
chg
∑
i
=
1
n
P
mc
_
chg
_
i
,
where η
-
mc is the average efficiency of motor; η
-
chg is the average efficiency of power battery charging; η
-
dischg is the average efficiency of power battery discharging.When the motor drives the vehicle, the energy consumption of power batteries can be converted into the engine fuel consumption. The equivalent fuel consumption rate of the motor is(2)
b
mc
_
eq
=
∑
i
=
1
n
(
P
fc
Y
_
i
b
Y
_
i
-
P
fc
N
_
i
b
N
_
i
)
η
¯
mc
2
η
¯
dischg
η
¯
chg
∑
i
=
1
n
P
mc
_
chg
_
i
,
where P
fc
N
_
i, when the engine does not charge power batteries, is the engine power of the i control period; b
N
_
i, when the engine does not charge power batteries, is the fuel consumption rate of the engine; P
fc
Y
_
i, when the engine charges power batteries, is the engine power of the i control period; b
Y
_
i, when the engine charges power batteries, is the fuel consumption rate of the engine.Letb
-
chg be the average fuel consumption rate when the engine charges power batteries:
(3)
b
¯
chg
=
∑
i
=
1
n
P
fc
Y
_
i
b
Y
_
i
-
P
fc
N
_
i
b
N
_
i
P
mc
_
chg
_
i
.Merge (3) and (2):
(4)
b
mc
_
eq
=
b
¯
chg
η
¯
mc
2
η
¯
dischg
η
¯
chg
.When the motor drives the vehicle after a period ofΔ
t, equivalent instantaneous fuel consumption of the motor is
(5)
m
mc
_
eq
=
b
¯
chg
P
mc
η
¯
mc
2
η
¯
dischg
η
¯
chg
.


### 2.1.2. Calculate the Fuel Consumption of Batteries Charging

When the power battery is charged by the engine, its SOC value will rise and even exceed the target of the SOC value. In order to maintain SOC values, power battery energy will be consumed in future [18]. Discharging time is divided into n control cycles, and the motor power of the each control cycle is P
mc
2
_
i  (
i
=
1,2
,
3
…
,
n
) [17].In a certain period of discharging time, the relationship between motor power (P
mc
2
_
i) with driving vehicle and motor power (P
mc
_
chg
_
2) with charging power batteries is
(6)
∑
i
=
1
n
P
mc
2
_
i
=
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
chg
2
η
¯
dischg
2
,
where P
mc
_
chg
2 is the motor power when the power battery is charged; η
-
mc
2 is the average efficiency of motor; η
-
chg
2 is the average efficiency of the power battery charging; η
-
dischg
2 is the average efficiency of the power battery discharging.When the motor drives the vehicle, the relationship between the motor power battery energy consumption and the fuel consumption rate is(7)
b
mc
_
eq
2
=
P
fc
Y
2
b
Y
2
-
P
fc
N
2
b
N
2
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
chg
2
η
¯
dischg
2
,
where P
fc
Y
2, when the engine charges the power battery, is the power of the engine; b
Y
2, when the engine charges the power battery, is the fuel consumption rate of the engine; P
fc
N
2, when the engine does not charge the power battery, is the power of the engine; b
N
2, when the engine does not charge the power battery, is the fuel consumption rate of the engine.Letb
chg
2 be the fuel consumption rate when the engine charges the power battery:
(8)
b
chg
2
=
P
fc
Y
2
b
Y
2
-
P
fc
N
2
b
N
2
P
mc
_
chg
2
.Simplify the (8) formula:
(9)
b
mc
_
eq
2
=
b
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
.When the motor charges power batteries after a period ofΔ
t, the objective function of equivalent instantaneous fuel consumption of the motor is
(10)
m
mc
_
eq
2
=
b
¯
chg
2
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
.


### 2.1.3. Deduce the Objective Function of the Instantaneous Optimal Control Strategy

Set two new variables:(11)
f
eq
_
dischg
=
b
¯
chg
η
¯
mc
2
η
¯
dischg
η
¯
chg
,
f
eq
_
chg
=
b
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
,
where b
chg
2 is the fuel consumption rate when the engine charges power batteries at the present moment; b
-
chg is the average fuel consumption rate when the engine charges power batteries in the future time.The instantaneous control objective function of the lowest fuel consumption is(12)
M
=
∑
Min
⁡
{
[
T
mc
(
t
)
,
m
mc
_
eq
ω
mc
(
t
)
]
Δ
t
m
fc
[
T
fc
(
t
)
,
ω
(
t
)
]
Δ
t
f
i
f
f
f
f
f
f
f
f
+
m
mc
_
eq
[
T
mc
(
t
)
,
ω
mc
(
t
)
]
Δ
t
}
,
where ω
fc
(
t
) is engine speed; ω
mc
(
t
) is motor speed; T
fc is output torque of the engine; T
mc is output torque of the motor. Consider
(13)
m
mc
_
eq
[
P
mc
(
t
)
]
=
{
f
eq
_
dischg
P
mc
T
mc
>
0
(
power
batteries
discharge
)
f
eq
_
chg
P
mc
_
chg
2
T
mc
<
0
(
power
batteries
charge
)
.P
mc and P
mc
_
chg
2 are both motor power, so they can be unified as P
mc. By calculating, the improved instantaneous control objective function of the minimum fuel consumption is
(14)
M
=
∑
Min
⁡
{
[
λ
f
eq
_
dischg
+
(
1
-
λ
)
f
eq
_
chg
]
m
fc
[
P
fc
(
t
)
]
Δ
t
f
f
f
f
f
f
f
i
i
f
+
[
λ
f
eq
_
dischg
+
(
1
-
λ
)
f
eq
_
chg
]
P
mc
Δ
t
}
,
where λ
=
(
1
+
sign
⁡
(
T
mc
)
)
/
2.


### 2.1.4. Improve the Objective Function of the Instantaneous Optimal Control Strategy

SOC value change of batteries and braking energy recovery both have a certain effect on energy control. The optimal function of the working point needs to be improved.(
1) Revise the SOC Value Function of Power Batteries. When power batteries work, their SOC value is maintained at the high efficient range by the reset function, in order to reduce the loss energy in the process of charging and discharging power batteries and make hybrid system keep better performances. The working principle of SOC reset function is as follows. When the SOC value is more than the target region, the hybrid power system will give priority to consuming power battery energy. It does not stop until SOC value decreases to the target region under the effect of the reset function. When the battery SOC value is lower than the target region, the hybrid system will give priority to consuming fuel energy to drive the vehicle and recover the value of SOC. It does not stop until SOC value returns to the target region under the effect of the reset function.SetK
SOC be the variable in the reset function, and the value table between K
SOC and SOC value is shown in Table 1.Table 1
The value table betweenK
soc and SOC.
SOC
0.11
0.14
0.2
0.25
0.3
0.35
K
soc
7
2
1.05
1.05
1.03
1.04
SOC
0.4
0.45
0.5
0.55
0.59
0.64
K
soc
1.02
1.02
1
0.99
0.98
0.96
SOC
0.7
0.74
0.81
0.85
0.9
K
soc
0.96
0.94
0. 94
0.8
0.2Based on the Matlab platform and Table1 data, the fitting curve between K
SOC and SOC value is constructed by using fitting curve toolbox, as shown in Figure 3.Figure 3
SOC andK
SOC relationship curve.Polynomial function of the fitting curve is(15)
K
soc
=
-
8866
x
7
+
33460
x
6
-
52310
x
5
+
43730
x
4
-
21010
x
3
+
5766
x
2
-
833.5
x
+
49.77
,
where x: is SOC value.Considering the influence of power battery SOC, the formula of instantaneous equivalent fuel consumption of HEV can be expressed asK
soc
m
mc
_
eq.(
2
) Revise the Objective Function of Braking Recovery Energy. The adopted method which revises the equivalent fuel consumption function of power battery energy is as follows: the average braking power is calculated at a period of time before the current moment. The power is used as the standard of fuel consumption correction of the braking recovery power in the next moment [19–22].The statistical time range is divided inton  (
n
>
0
) cycles. In each cycle, let braking power be a fixed value. Therefore, the average braking power of the whole time can be expressed as
(16)
P
-
braking
=
∑
i
=
1
n
P
braking
_
i
n
.When power batteries discharge, the objective function of instantaneous equivalent fuel consumption is(17)
m
mc
_
eq
=
b
-
chg
(
P
mc
+
P
-
braking
)
η
-
mc
2
η
-
dischg
η
-
chg
.In summary, taking the influence of power battery SOC and brake energy recovery into consideration, the final objective function of instantaneous optimal control strategy is(18)
M
=
∑
Min
⁡
{
[
P
mc
(
t
)
+
P
-
braking
(
t
)
]
m
fc
[
P
fc
(
t
)
]
Δ
t
fffffffffffffff
+
K
soc
m
mc
_
eq
[
P
mc
(
t
)
+
P
-
braking
(
t
)
]
Δ
t
}
.


## 2.1.1. Calculate the Equivalent Fuel Consumption of Battery Discharging

When the power battery takes part in driving hybrid cars, its SOC value will reduce and deviate from the target of SOC value. In order to compensate for the used electricity and restore SOC value of power batteries, the engine drives the motor to charge power batteries in the future time [18]. The charging time is divided into n control cycles, and the motor power of each cycle is P
mc
_
chg
_
i  (
i
=
1,2
,
3
…
,
n
) [5].The relationship between the motor power (P
mc) with driving vehicle and the motor power (P
mc
_
chg
_
i) with power batteries charging is
(1)
P
mc
=
η
-
mc
2
η
-
dischg
η
-
chg
∑
i
=
1
n
P
mc
_
chg
_
i
,
where η
-
mc is the average efficiency of motor; η
-
chg is the average efficiency of power battery charging; η
-
dischg is the average efficiency of power battery discharging.When the motor drives the vehicle, the energy consumption of power batteries can be converted into the engine fuel consumption. The equivalent fuel consumption rate of the motor is(2)
b
mc
_
eq
=
∑
i
=
1
n
(
P
fc
Y
_
i
b
Y
_
i
-
P
fc
N
_
i
b
N
_
i
)
η
¯
mc
2
η
¯
dischg
η
¯
chg
∑
i
=
1
n
P
mc
_
chg
_
i
,
where P
fc
N
_
i, when the engine does not charge power batteries, is the engine power of the i control period; b
N
_
i, when the engine does not charge power batteries, is the fuel consumption rate of the engine; P
fc
Y
_
i, when the engine charges power batteries, is the engine power of the i control period; b
Y
_
i, when the engine charges power batteries, is the fuel consumption rate of the engine.Letb
-
chg be the average fuel consumption rate when the engine charges power batteries:
(3)
b
¯
chg
=
∑
i
=
1
n
P
fc
Y
_
i
b
Y
_
i
-
P
fc
N
_
i
b
N
_
i
P
mc
_
chg
_
i
.Merge (3) and (2):
(4)
b
mc
_
eq
=
b
¯
chg
η
¯
mc
2
η
¯
dischg
η
¯
chg
.When the motor drives the vehicle after a period ofΔ
t, equivalent instantaneous fuel consumption of the motor is
(5)
m
mc
_
eq
=
b
¯
chg
P
mc
η
¯
mc
2
η
¯
dischg
η
¯
chg
.


## 2.1.2. Calculate the Fuel Consumption of Batteries Charging

When the power battery is charged by the engine, its SOC value will rise and even exceed the target of the SOC value. In order to maintain SOC values, power battery energy will be consumed in future [18]. Discharging time is divided into n control cycles, and the motor power of the each control cycle is P
mc
2
_
i  (
i
=
1,2
,
3
…
,
n
) [17].In a certain period of discharging time, the relationship between motor power (P
mc
2
_
i) with driving vehicle and motor power (P
mc
_
chg
_
2) with charging power batteries is
(6)
∑
i
=
1
n
P
mc
2
_
i
=
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
chg
2
η
¯
dischg
2
,
where P
mc
_
chg
2 is the motor power when the power battery is charged; η
-
mc
2 is the average efficiency of motor; η
-
chg
2 is the average efficiency of the power battery charging; η
-
dischg
2 is the average efficiency of the power battery discharging.When the motor drives the vehicle, the relationship between the motor power battery energy consumption and the fuel consumption rate is(7)
b
mc
_
eq
2
=
P
fc
Y
2
b
Y
2
-
P
fc
N
2
b
N
2
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
chg
2
η
¯
dischg
2
,
where P
fc
Y
2, when the engine charges the power battery, is the power of the engine; b
Y
2, when the engine charges the power battery, is the fuel consumption rate of the engine; P
fc
N
2, when the engine does not charge the power battery, is the power of the engine; b
N
2, when the engine does not charge the power battery, is the fuel consumption rate of the engine.Letb
chg
2 be the fuel consumption rate when the engine charges the power battery:
(8)
b
chg
2
=
P
fc
Y
2
b
Y
2
-
P
fc
N
2
b
N
2
P
mc
_
chg
2
.Simplify the (8) formula:
(9)
b
mc
_
eq
2
=
b
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
.When the motor charges power batteries after a period ofΔ
t, the objective function of equivalent instantaneous fuel consumption of the motor is
(10)
m
mc
_
eq
2
=
b
¯
chg
2
P
mc
_
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
.


## 2.1.3. Deduce the Objective Function of the Instantaneous Optimal Control Strategy

Set two new variables:(11)
f
eq
_
dischg
=
b
¯
chg
η
¯
mc
2
η
¯
dischg
η
¯
chg
,
f
eq
_
chg
=
b
chg
2
η
¯
mc
2
2
η
¯
dischg
2
η
¯
chg
2
,
where b
chg
2 is the fuel consumption rate when the engine charges power batteries at the present moment; b
-
chg is the average fuel consumption rate when the engine charges power batteries in the future time.The instantaneous control objective function of the lowest fuel consumption is(12)
M
=
∑
Min
⁡
{
[
T
mc
(
t
)
,
m
mc
_
eq
ω
mc
(
t
)
]
Δ
t
m
fc
[
T
fc
(
t
)
,
ω
(
t
)
]
Δ
t
f
i
f
f
f
f
f
f
f
f
+
m
mc
_
eq
[
T
mc
(
t
)
,
ω
mc
(
t
)
]
Δ
t
}
,
where ω
fc
(
t
) is engine speed; ω
mc
(
t
) is motor speed; T
fc is output torque of the engine; T
mc is output torque of the motor. Consider
(13)
m
mc
_
eq
[
P
mc
(
t
)
]
=
{
f
eq
_
dischg
P
mc
T
mc
>
0
(
power
batteries
discharge
)
f
eq
_
chg
P
mc
_
chg
2
T
mc
<
0
(
power
batteries
charge
)
.P
mc and P
mc
_
chg
2 are both motor power, so they can be unified as P
mc. By calculating, the improved instantaneous control objective function of the minimum fuel consumption is
(14)
M
=
∑
Min
⁡
{
[
λ
f
eq
_
dischg
+
(
1
-
λ
)
f
eq
_
chg
]
m
fc
[
P
fc
(
t
)
]
Δ
t
f
f
f
f
f
f
f
i
i
f
+
[
λ
f
eq
_
dischg
+
(
1
-
λ
)
f
eq
_
chg
]
P
mc
Δ
t
}
,
where λ
=
(
1
+
sign
⁡
(
T
mc
)
)
/
2.


## 2.1.4. Improve the Objective Function of the Instantaneous Optimal Control Strategy

SOC value change of batteries and braking energy recovery both have a certain effect on energy control. The optimal function of the working point needs to be improved.(
1) Revise the SOC Value Function of Power Batteries. When power batteries work, their SOC value is maintained at the high efficient range by the reset function, in order to reduce the loss energy in the process of charging and discharging power batteries and make hybrid system keep better performances. The working principle of SOC reset function is as follows. When the SOC value is more than the target region, the hybrid power system will give priority to consuming power battery energy. It does not stop until SOC value decreases to the target region under the effect of the reset function. When the battery SOC value is lower than the target region, the hybrid system will give priority to consuming fuel energy to drive the vehicle and recover the value of SOC. It does not stop until SOC value returns to the target region under the effect of the reset function.SetK
SOC be the variable in the reset function, and the value table between K
SOC and SOC value is shown in Table 1.Table 1
The value table betweenK
soc and SOC.
SOC
0.11
0.14
0.2
0.25
0.3
0.35
K
soc
7
2
1.05
1.05
1.03
1.04
SOC
0.4
0.45
0.5
0.55
0.59
0.64
K
soc
1.02
1.02
1
0.99
0.98
0.96
SOC
0.7
0.74
0.81
0.85
0.9
K
soc
0.96
0.94
0. 94
0.8
0.2Based on the Matlab platform and Table1 data, the fitting curve between K
SOC and SOC value is constructed by using fitting curve toolbox, as shown in Figure 3.Figure 3
SOC andK
SOC relationship curve.Polynomial function of the fitting curve is(15)
K
soc
=
-
8866
x
7
+
33460
x
6
-
52310
x
5
+
43730
x
4
-
21010
x
3
+
5766
x
2
-
833.5
x
+
49.77
,
where x: is SOC value.Considering the influence of power battery SOC, the formula of instantaneous equivalent fuel consumption of HEV can be expressed asK
soc
m
mc
_
eq.(
2
) Revise the Objective Function of Braking Recovery Energy. The adopted method which revises the equivalent fuel consumption function of power battery energy is as follows: the average braking power is calculated at a period of time before the current moment. The power is used as the standard of fuel consumption correction of the braking recovery power in the next moment [19–22].The statistical time range is divided inton  (
n
>
0
) cycles. In each cycle, let braking power be a fixed value. Therefore, the average braking power of the whole time can be expressed as
(16)
P
-
braking
=
∑
i
=
1
n
P
braking
_
i
n
.When power batteries discharge, the objective function of instantaneous equivalent fuel consumption is(17)
m
mc
_
eq
=
b
-
chg
(
P
mc
+
P
-
braking
)
η
-
mc
2
η
-
dischg
η
-
chg
.In summary, taking the influence of power battery SOC and brake energy recovery into consideration, the final objective function of instantaneous optimal control strategy is(18)
M
=
∑
Min
⁡
{
[
P
mc
(
t
)
+
P
-
braking
(
t
)
]
m
fc
[
P
fc
(
t
)
]
Δ
t
fffffffffffffff
+
K
soc
m
mc
_
eq
[
P
mc
(
t
)
+
P
-
braking
(
t
)
]
Δ
t
}
.


## 2.2. Elman Neural Network Control Strategy

### 2.2.1. The Structure of Elman Neural Network

Elman neural network is put forward by Jeffrey L. Elman in 1990 and is a typical local recessionary grid, as shown in Figure4 [23].Figure 4
Elman neural network model.P is the input of the neural network and its size is R
×
1; b
1 is the neuronal threshold vector of the feedback layer and its size is S
1
×
1; I
W
1 is the connective weight vector of the neurons and the input vector in the input layer and its size is S
1
×
R; n
1 is the middle operational result of the neurons in the feedback layer, namely, weighted sum of the connective weight vector and the threshold vector, and its size is S
1
×
1; a
1 is the output vector of the feedback layer’s neuron in the K iteration and its size is S
1
×
1; D is feedback node. It is the same way that  b
2, L
W
2, n
2, and a
2 are related to parameters of the output layer.The inputs of structure diagram of Elman neural network is the required torque, speed, and SOC value of power battery and its output is motor torque. The structure diagram of Elman neural network is shown in Figure5 [24].Figure 5
Structure diagram of Elman neural network.The structure diagram of Elman neural network contains input layer, hidden layer, undertaken layer, and output layer. Let the input vector of the input layer be three-dimensional vectoru; the output vector of the output layer is one-dimensional vector y; the output vector of the hidden layer is n-dimensional vector x; the output vector of the undertaken layer is n-dimensional vector r; w
1, w
2, and w
3 are respective connective weights of the hidden layer to the output layer, the input layer to the hidden layer, and the undertaken layer to the hidden layer; g
(
·
) is the driving function of the output neurons; f
(
·
) is the driving function of hidden layer; h
(
·
) is the driving function of undertaken layer; net(
·
) is the net input driving function of a certain layer; A shows the input layer; B shows the undertaken layer; K shows the iterative sequence.Define two functions:(19)
v
i
(
k
)
=
{
u
n
(
k
)
,
if
i
∈
A
,
r
n
(
k
)
,
if
i
∈
B
,
w
i
(
k
)
=
{
w
2
,
if
i
∈
A
,
w
3
,
if
i
∈
B
.The input and output functions of theN neuron of the hidden layer are
(20)
net
n
(
k
+
1
)
=
∑
i
∈
A
∪
B
w
i
(
k
)
v
i
(
k
)
,
x
n
(
n
+
1
)
=
f
(
net
n
(
k
+
1
)
)
.The input and output functions of theN neuron of the hidden layer are
(21)
net
n
(
k
)
=
∑
i
∈
A
∪
B
w
i
(
k
-
1
)
v
i
(
k
-
1
)
,
r
(
k
)
=
h
(
net
n
(
k
)
)
.The input and output functions of the output layer’s neuron are(22)
net
n
(
k
+
1
)
=
∑
i
∈
A
∪
B
w
i
(
k
+
1
)
x
n
(
k
+
1
)
,
y
(
k
+
1
)
=
g
(
net
(
k
+
1
)
)
.


### 2.2.2. Select the Parameters of Elman Neural Network

The neuron number is determined by following formula [25]:
(23)
k
=
m
+
n
+
β
,
where m is the number of the input vector; n is the neuron number of the output vector; β is a constant, (1~10).The excitation function of Elman neural network of the feedback layer selects the Tansig function [26]:
(24)
tansig
(
x
)
=
2
1
+
e
-
2
x
-
1
.


### 2.2.3. Learning and Training Mechanism of Elman Neural Network

Elman neural network is trained by Levenberg-Marquardt algorithm. The error index function of Levenberg-Marquardt arithmetic is(25)
E
(
w
)
=
1
2
∑
i
=
1
p
∥
y
i
-
y
i
′
∥
2
=
1
2
∑
i
=
1
p
e
i
2
(
w
)
,
where p is the sample number; e
i is the systemic error; y
i
′ is the actual output of the network.The formula of the adjusting weight is(26)
w
k
+
1
=
w
k
+
Δ
w
.The computing formula of the increment weight is(27)
Δ
w
=
[
J
T
(
w
)
J
(
w
)
+
u
I
]
-
1
J
T
(
w
)
e
(
w
)
,
where u is learning rate; I is the unit matrix; J
(
w
) is the Jacobian matrix. Consider
(28)
J
(
w
)
=
[
∂
e
1
(
w
)
∂
w
1
∂
e
1
(
w
)
∂
w
2
⋯
∂
e
1
(
w
)
∂
w
n
∂
e
2
(
w
)
∂
w
1
∂
e
2
(
w
)
∂
w
1
⋯
∂
e
2
(
w
)
∂
w
n
⋮
⋮
⋱
⋮
∂
e
n
(
w
)
∂
w
1
∂
e
n
(
w
)
∂
w
2
⋯
∂
e
n
(
w
)
∂
w
n
]
.Neural network weight is adjusted by the Levenberg-Marquardt algorithm, and the adjusting process is shown in Figure6.Figure 6
Process map of the adjusting network weight.


## 2.2.1. The Structure of Elman Neural Network

Elman neural network is put forward by Jeffrey L. Elman in 1990 and is a typical local recessionary grid, as shown in Figure4 [23].Figure 4
Elman neural network model.P is the input of the neural network and its size is R
×
1; b
1 is the neuronal threshold vector of the feedback layer and its size is S
1
×
1; I
W
1 is the connective weight vector of the neurons and the input vector in the input layer and its size is S
1
×
R; n
1 is the middle operational result of the neurons in the feedback layer, namely, weighted sum of the connective weight vector and the threshold vector, and its size is S
1
×
1; a
1 is the output vector of the feedback layer’s neuron in the K iteration and its size is S
1
×
1; D is feedback node. It is the same way that  b
2, L
W
2, n
2, and a
2 are related to parameters of the output layer.The inputs of structure diagram of Elman neural network is the required torque, speed, and SOC value of power battery and its output is motor torque. The structure diagram of Elman neural network is shown in Figure5 [24].Figure 5
Structure diagram of Elman neural network.The structure diagram of Elman neural network contains input layer, hidden layer, undertaken layer, and output layer. Let the input vector of the input layer be three-dimensional vectoru; the output vector of the output layer is one-dimensional vector y; the output vector of the hidden layer is n-dimensional vector x; the output vector of the undertaken layer is n-dimensional vector r; w
1, w
2, and w
3 are respective connective weights of the hidden layer to the output layer, the input layer to the hidden layer, and the undertaken layer to the hidden layer; g
(
·
) is the driving function of the output neurons; f
(
·
) is the driving function of hidden layer; h
(
·
) is the driving function of undertaken layer; net(
·
) is the net input driving function of a certain layer; A shows the input layer; B shows the undertaken layer; K shows the iterative sequence.Define two functions:(19)
v
i
(
k
)
=
{
u
n
(
k
)
,
if
i
∈
A
,
r
n
(
k
)
,
if
i
∈
B
,
w
i
(
k
)
=
{
w
2
,
if
i
∈
A
,
w
3
,
if
i
∈
B
.The input and output functions of theN neuron of the hidden layer are
(20)
net
n
(
k
+
1
)
=
∑
i
∈
A
∪
B
w
i
(
k
)
v
i
(
k
)
,
x
n
(
n
+
1
)
=
f
(
net
n
(
k
+
1
)
)
.The input and output functions of theN neuron of the hidden layer are
(21)
net
n
(
k
)
=
∑
i
∈
A
∪
B
w
i
(
k
-
1
)
v
i
(
k
-
1
)
,
r
(
k
)
=
h
(
net
n
(
k
)
)
.The input and output functions of the output layer’s neuron are(22)
net
n
(
k
+
1
)
=
∑
i
∈
A
∪
B
w
i
(
k
+
1
)
x
n
(
k
+
1
)
,
y
(
k
+
1
)
=
g
(
net
(
k
+
1
)
)
.


## 2.2.2. Select the Parameters of Elman Neural Network

The neuron number is determined by following formula [25]:
(23)
k
=
m
+
n
+
β
,
where m is the number of the input vector; n is the neuron number of the output vector; β is a constant, (1~10).The excitation function of Elman neural network of the feedback layer selects the Tansig function [26]:
(24)
tansig
(
x
)
=
2
1
+
e
-
2
x
-
1
.


## 2.2.3. Learning and Training Mechanism of Elman Neural Network

Elman neural network is trained by Levenberg-Marquardt algorithm. The error index function of Levenberg-Marquardt arithmetic is(25)
E
(
w
)
=
1
2
∑
i
=
1
p
∥
y
i
-
y
i
′
∥
2
=
1
2
∑
i
=
1
p
e
i
2
(
w
)
,
where p is the sample number; e
i is the systemic error; y
i
′ is the actual output of the network.The formula of the adjusting weight is(26)
w
k
+
1
=
w
k
+
Δ
w
.The computing formula of the increment weight is(27)
Δ
w
=
[
J
T
(
w
)
J
(
w
)
+
u
I
]
-
1
J
T
(
w
)
e
(
w
)
,
where u is learning rate; I is the unit matrix; J
(
w
) is the Jacobian matrix. Consider
(28)
J
(
w
)
=
[
∂
e
1
(
w
)
∂
w
1
∂
e
1
(
w
)
∂
w
2
⋯
∂
e
1
(
w
)
∂
w
n
∂
e
2
(
w
)
∂
w
1
∂
e
2
(
w
)
∂
w
1
⋯
∂
e
2
(
w
)
∂
w
n
⋮
⋮
⋱
⋮
∂
e
n
(
w
)
∂
w
1
∂
e
n
(
w
)
∂
w
2
⋯
∂
e
n
(
w
)
∂
w
n
]
.Neural network weight is adjusted by the Levenberg-Marquardt algorithm, and the adjusting process is shown in Figure6.Figure 6
Process map of the adjusting network weight.


## 3. Modeling of Energy Control Strategy

### 3.1. Modeling of Instantaneous Optimal Control Strategy

On the platform of Matlab/simulink, the instantaneous control model is established and mainly includes two parts: the calculation model of average braking energy of recovered power and the calculative model of optimal working point [25, 26].First, the calculation model of average braking energy of the recovery power is established as shown in Figure7. Inputs of the model are the required torque and rotational speed and its output is the calculated average recovery power of braking energy in the transmission system.Figure 7
Flow diagram of average recovery power of braking energy.Second, the calculation model of the revised SOC value of power battery is established, as shown in Figure8. The model can maintain the SOC value of power batteries at a high efficient range.Figure 8
Flow diagram of the revised SOC value of power batteries.Third, the calculation model of the optimal operating point is established, as shown in Figure9. Inputs of the model are the required torque, rotational speed, modified S
target-motor value of power batteries, and the average braking recovery power and the output torque of the motor is the optimal allocation power between the engine and motor in the transmission system.Figure 9
Flow diagram of the optimal operating point.Finally, the whole simulation diagram of the instantaneous control strategy is established, as shown in Figure10.Figure 10
Total simulation diagram of the instantaneous control strategy.


### 3.2. The Elman Neural Network Model

Elman neural network is gradually learning to do something by imitating human brain. Its essence is to improve the learned knowledge and the neurons weight [27–30]. The flow diagram of the trained Elman neural network is shown in Figure 11.Figure 11
The training flow diagram of the trained Elman neural network.


## 3.1. Modeling of Instantaneous Optimal Control Strategy

On the platform of Matlab/simulink, the instantaneous control model is established and mainly includes two parts: the calculation model of average braking energy of recovered power and the calculative model of optimal working point [25, 26].First, the calculation model of average braking energy of the recovery power is established as shown in Figure7. Inputs of the model are the required torque and rotational speed and its output is the calculated average recovery power of braking energy in the transmission system.Figure 7
Flow diagram of average recovery power of braking energy.Second, the calculation model of the revised SOC value of power battery is established, as shown in Figure8. The model can maintain the SOC value of power batteries at a high efficient range.Figure 8
Flow diagram of the revised SOC value of power batteries.Third, the calculation model of the optimal operating point is established, as shown in Figure9. Inputs of the model are the required torque, rotational speed, modified S
target-motor value of power batteries, and the average braking recovery power and the output torque of the motor is the optimal allocation power between the engine and motor in the transmission system.Figure 9
Flow diagram of the optimal operating point.Finally, the whole simulation diagram of the instantaneous control strategy is established, as shown in Figure10.Figure 10
Total simulation diagram of the instantaneous control strategy.


## 3.2. The Elman Neural Network Model

Elman neural network is gradually learning to do something by imitating human brain. Its essence is to improve the learned knowledge and the neurons weight [27–30]. The flow diagram of the trained Elman neural network is shown in Figure 11.Figure 11
The training flow diagram of the trained Elman neural network.


## 4. Simulation Results and Analysis

The basic vehicle parameters are shown in Table2.Table 2
Parallel hybrid electric vehicle parameters.
Vehicle
Curb weight
1605 kg
Face area
2.65 m²
Wheel base
2.775 m
Height of the center of mass
0.5 m
Front axle load distribution ratio
0.51
Coefficient drag
0.32
Engine
Peak power
118 kW
Displacement
2.5 L
Power battery pack
Voltage
244.8/650 V
Style
NI-MH
Volume
6.5 Ah
Mold number
34
Motor
Peak power
105 kW
Style
PMSMTraffic parameters of simulation experiments are described in Table3.Table 3
Traffic parameters of simulation experiments.
Parameter
NEDC
HWFET
Idle time (s)
298
6
Top speed (km/h)
10.93
16.51
Cycle time (s)
1184
765
Average speed (km/h)
33.21
77.58
Maximum acceleration (m/s2)
1.06
1.43
Maximum deceleration (m/s2)
−1.39
−1.48
Park time (time)
13
1
Traveling distance (km)
120
96.4The original control model is replaced by the instantaneous optimal control model in ADVISOR 2002. Then the trained Elman neural network controller is imported to the software [31, 32].Simulation results are shown in Figure13. HWFET_mdf and NEDC_mdf show simulation results of the instantaneous optimal control strategy of the PHEV model (Figure 12); HWFET_net and NEDC_net show simulation results of the Elman neural network controller of the PHEV model.Figure 12
Simulation model of the parallel hybrid electric vehicle.Contrast power battery SOC.
(a)
NEDC working condition
(b)
HWFET working conditionCompared with the instantaneous optimal control strategy, the Elman neural network strategy can slow down the decline of SOC value and maintain SOC value at high efficient range on the NEDC working condition in Figure13(a). Meanwhile, the Elman neural network strategy has approximate SOC value change compared with the instantaneous optimal control strategy on the HWFET working in Figure 13(b).As shown in Figure13(a), Elman neural network strategy can make the engine produce more torque than the instantaneous optimal control strategy at the beginning of 600 s on the NEDC working condition. It can make the vehicle start, accelerate, and climb better. After the vehicle starts, the two control strategies play the same role on the vehicle energy control. As seen in Figure 13(b), Elman neural network strategy can make the engine produce slightly more torque than instantaneous optimal control strategy on the HWFET working condition, while these two strategies have similar effect on the engine torque control.As seen in Figure15, Elman neural network strategy can make motor produce slightly more torque than the instantaneous optimal control strategy at some moments on the NEDC and HWFET working condition, while the two strategies have the similar effect on the motor torque control.As shown in Table4, compared with the instantaneous optimal neural strategy, fuel consumption of Elman neural network strategy only increases about 0.5 (L/100 km) on the NEDC and HWFEF working condition, which implies that Elman neural network controller can also have the advantage of low fuel consumption. The slight increase in fuel consumption can be accepted since it has a little effect on the fuel economy of the whole vehicle. As seen in Table 5, Elman neural network strategy makes the simulation time decreased greatly compared with the instantaneous optimal neural strategy and improves the response time of the vehicle greatly.Table 4
Fuel consumption of 100 km (L/100 km).
Strategy
Road
NEDC
HWFEF
Instantaneous optimal control
9.4
6.5
Elman neural network
9.8
7Table 5
Simulation time (s).
Strategy
Road
NEDC
HWFEF
Instantaneous optimal control
471.3
315.8
Elman neural network
15.6
10.2In conclusion, as seen in Figures13, 14, and 15, Elman neural network strategy can replace the instantaneous optimal control strategy to maintain SOC value at the high efficient range and achieve a reasonable distribution of the torque between the engine and the motor. The significance of the paper is that the simulation time of energy control is reduced by 96%.Contrast engine torque.
(a)
NEDC working condition
(b)
HWFET working conditionContrast the motor torque.
(a)
NEDC working condition
(b)
HWFET working condition


## 5. Conclusion

Through the research on the instantaneous optimal strategy and Elman neural network control strategy, we deduce the objective functions of instantaneous optimal control and establish the instantaneous control model and design the Elman controller. Based on the ADVISOR 2002 platform, two control strategies were simulated on a hybrid electric vehicle.It is seen from the simulation results that the trained Elman neural network strategy shows similar control ability on the vehicle energy distribution compared with the instantaneous optimal control strategy, which ensures good performances of power and fuel economy of HEV, reduces the control reaction time greatly, and overcomes the disadvantage of poor real-time performance of the instantaneous optimal control strategy. The research significance of the paper is that the simulation time of energy control is reduced by 96%.Future works are listed as below.(1)
Simulation and experiment should be improved by adding more design parameters, such as vehicle emission.
(2)
It is necessary to do lots of experiments to enrich simulation results.
(3)
Actual road condition is more complex than the simulation road condition, so control strategies need to be tested in the actual road conditions.



---

*Source: 596326-2014-08-26.xml*