# Design of a Load Torque Based Control Strategy for Improving Electric Tractor Motor Energy Conversion Efficiency

Authors: Mengnan Liu; Liyou Xu; Zhili Zhou
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

In order to improve the electrical conversion efficiency of an electric tractor motor, a load torque based control strategy (LTCS) is designed in this paper by using a particle swarm optimization algorithm (PSO). By mathematically modeling electric-mechanical performance and theoretical energy waste of the electric motor, as well as the transmission characteristics of the drivetrain, the objective function, control relationship, and analytical platform are established. Torque and rotation speed of the motor’s output shaft are defined as manipulated variables. LTCS searches the working points corresponding to the best energy conversion efficiency via PSO to control the running status of the electric motor and uses logic and fuzzy rules to fit the search initialization for load torque fluctuation. After using different plowing forces to imitate all the common tillage forces, the simulation of traction experiment is conducted, which proves that LTCS can make the tractor use electrical power efficiently and maintain agricultural applicability on farmland conditions. It provides a novel method of fabricating a more efficient electric motor used in the traction of an off-road vehicle.


---

## Body

## 1. Introduction

A potential development trend of farming power is that sustainable energies, such as electricity and biodiesel, will gradually and eventually replace the nonrenewable fossil fuels. An electric tractor is one of the transformation directions. Carlini et al. presented a series hybrid electric powertrain for a small crawler tractor used for logging in forests. It reached a 30% increase in fuel economy and better ecological features in comparison with the traditional tractor [1]. Florentsev et al. developed an electromechanical drive train for a 300 hp class tractor. The plowing production was 2% higher compared to the traditional tractor with power-shift drive train, while fuel consumption was lower by 18% [2]. Based on our research in the static analysis of a 180 hp class electric tractor, the tractive efficiency average increased by 2.98% [3]. These achievements show that the electric tractor has the structural advantage in energy conservation. However, in its drivetrain, an electric motor acts as the engine whose energy efficiency tightly links to the running status. In fact, owing to the diversity of farming works, complexity of soil characteristics, and variation of driving resistance, the load torque of the drive wheel changes frequently and suddenly when tractors work on the farm [4]. Consequently, the motor cannot always maintain high-level energy conversion efficiency. Focused on this problem, control strategies designed for promoting the energy efficiency of electric motors meet the developing requirement of resource-saving and environment friendly agricultural machinery technology.Related studies have been performed for improving the energy conversion economy of electric vehicles. By analyzing five types of control strategies used in high-power vehicles, Garcia et al. proved that the equivalent fuel consumption minimization strategy had better applicability [5]. Via comparison with a batch Genetic Algorithm-based optimization, Sorrentino et al. designed a rule-based control strategy for on-board energy management of series hybrid vehicles, which could improve the fuel economy of hybrid solar vehicles [6]. However, researches like these only refer to the management optimization of energy systems that included the battery, the electric generator powered by an internal combustion engine (ICE), and other assistance energy sources. They also do not take electric energy-using device into consideration, which means such control strategies are not suitable for the energy conversion economy improvement of the electric vehicle with a simple energy system, such as a purely electric vehicle. To account for this, Shabbir and Evangelou investigated a real-time control strategy, which uses a control map computed off-line to maximize the drivetrain efficiency of hybrid electric vehicles [7]. Another control strategy was proposed by Osornio-Correa et al. based on the heuristic map created to analyze the restrictions and benefits of using either of the on-board power plants under different driving conditions for maximum energy economy of a parallel hybrid electric vehicle [8]. However, those researches lack discussion of how to improve the EM energy conversion efficiency. In addition, the establishment of a control map requires extensive experimental data as a practical foundation (which means it would be difficult to use in product development phases of the vehicle). Wu et al. optimized principal parameters of both the powertrain and the control strategy using a PSO and achieved an acceptable result. However, the optimization of electric motor control strategy was also ignored [9]. Elwer et al. optimized the scaling factors of a fuzzy controller used for enhancing the rapidity, accuracy, and robustness of the electric vehicle motor [10]. However, the controller is not helpful in the improvement of motor energy conversion efficiency. By using an adaptive network in a fuzzy inference system optimization algorithm, the strategy put forth by Qian et al. made the ICE of series hybrid electric vehicle work with high energy conversion efficiency [11]. Unfortunately, in a series hybrid electric vehicle, ICE is mechanically decoupled to the drive wheel, and its mechanical power output is not constrained by load characteristic. Therefore, the control strategy cannot be imitatively used in traction motor. In an ordinary electric drivetrain or pure electric vehicle, an electric motor works as the sole energy conversion device, so the energy conversion efficiency of the drivetrain (or the vehicle) is only promoted by means of controlling the working state of the electric motor. Nevertheless, the energy conversion efficiency of the electric motor tightly couples with its rotate speed and torque, which are also dependent upon the actual driving demand of the vehicle. For the road vehicle, complex and frequently changing road condition means it is impractical that control of the mechanical characteristics of the electric motor is dependent upon its energy conversion efficiency, while neglecting the speed requirements of safety, practicability, and traffic laws. However, in the agricultural conditions of the tractor, the frequent change of speed is not required, and the load force is normally stable. Therefore, we can try to control the electric motor’s running state for improved energy conversion efficiency. In summary, research on an efficient control strategy for an electric tractor motor is both needed and innovative.This paper begins with establishing the mathematical model that describes the electric-mechanical performance and theoretical energy consumption of the electric motor, as well as the transmission characteristics of the electric tractor. Then, it introduces the design of LTCS and how it can control the electric motor to maintain high energy conversion efficiency. Finally, simulation of the traction experiment is conducted to verify the effect and agricultural applicability of LTCS.


## 2. Mathematical Model

### 2.1. Introduction of the Electric Tractor Drivetrain

Although the electric tractor shares many common structural features with the common electric vehicle, there are some significant differences between the two drivetrains. Particularly, the electric vehicle layout is organized as shown in Figure1. The electric tractor drivetrain is powered by the energy system (ES). The ES can contain the storage battery in a pure electric tractor or the ICE, generator, battery, super-capacitor, and so forth in a hybrid electric tractor. The electric motor (EM) provides the mechanical power for the drive wheel and power take-off shaft (PTO). PTO joints with the EM via the clutch. When the tractor only drives the agricultural implement (AI), the clutch decouples the mechanical connection between the PTO and the drive wheel, which makes the mechanical power flows of them be independent of each other. The gearbox transmission (GT) can transform the dynamic power performance field of the drive wheel via the gear shifting operation of driver. And the GT directly transmits the mechanical power to PTO, which can make the rotate speed of PTO match with either the rotate speed of the EM or the forward velocity of the electric tractor. The main drive (MD) acts as the differential mechanism and final reducer of the transmission system. By stepping on the accelerator pedal (AP), the driver can control the EM through electrical adjustment of the motor controller (MC). Generally, there are two types of load force exerted upon the electric tractor: one is driving force and the other is tilling force.Figure 1
Scheme of an electric tractor’s drivetrain.


### 2.2. Electric-Mechanical Performance of the Electric Motor (EM)

According to the previous research, the EM is a brushless direct current motor (BLDC) [3]. The voltage per phase can be expressed by(1)
U
a
=
R
s
I
a
+
L
a
-
L
m
d
I
a
d
t
+
e
a
,
U
b
=
R
s
I
b
+
L
b
-
L
m
d
I
b
d
t
+
e
b
,
U
c
=
R
s
I
c
+
L
c
-
L
m
d
I
c
d
t
+
e
c
.After equivalent transformation, the electromechanical characteristics of BLDC can be described by the following equations [12]:(2)
U
t
=
R
s
I
s
+
L
n
-
L
m
d
I
s
d
t
+
e
s
,
L
s
=
L
n
-
L
m
,
e
s
=
k
E
ω
r
.The electromagnetic torque of BLDC,T
e, can be obtained as follows:(3)
T
e
=
9.55
e
a
I
a
+
e
b
I
b
+
e
c
I
c
ω
r
=
9.55
e
s
I
s
ω
r
,
(4)
T
e
=
k
T
I
s
.Depending on (3) and (4), the torque constant, k
T, is equal to back EMF constant, k
E, after the dimension change.When BLDC is working on traction, the electromagnetic torque of the motor can also be described by(5)
T
e
=
T
L
+
J
d
ω
r
d
t
+
V
ω
r
.Based on the Laplace transform of (1)–(5), the electric-mechanical performance of BLDC can be expressed by the following transfer function:(6)
ω
r
s
=
k
T
R
s
+
s
L
s
s
J
+
V
+
k
E
k
T
U
t
s
-
R
s
+
s
L
s
R
s
+
s
L
s
s
J
+
V
+
k
E
k
T
T
L
s
.Consequently, the steady state characteristic can be established:(7)
ω
r
=
k
T
R
s
V
+
k
E
k
T
U
t
-
R
s
R
s
V
+
k
E
k
T
T
L
.According to (4), (6), and (7), the mechanical performance of the EM can be linked to the controllable electric parameter.


### 2.3. Energy Consumption Model of the EM

When the EM is running, electric power loss contains copper loss, icon loss, mechanical loss, and constant power used for driving the MC.Copper loss,P
c, is the prime waste of the BLDC. This is due to the fact that the stator winding impedance converts the electric power to the thermal dissipated to the surroundings. It can be obtained as follows [13]: (8)
P
c
=
I
s
2
R
s
.According to (4), the value of the current is directly proportional to the value of electromagnetic torque. Thus, we can use a coefficient, k
c, to describe the relationship between the copper loss and the electromagnetic torque. The copper loss can be rewritten as follows:(9)
P
c
=
k
c
T
e
2
.The copper loss coefficient,k
c, can be calculated by(10)
k
c
=
R
s
k
T
2
.Icon loss power,P
i, mainly contains magnetic hysteresis loss and current loss, and its mechanisms are complex and changing. It mainly depends on the speed of rotor and is described by(11)
P
i
=
k
i
ω
r
.In the BLDC,k
i can be approximately seen as a constant and can be measured by an experimental method.The mechanical lossP
M, output from the rotor, is largely wasted on air resistance. It can be obtained by(12)
P
M
=
k
w
ω
r
3
.The mechanical power output and total power input of the electric motor can be described by the following equations:(13)
P
o
u
t
=
T
L
ω
r
9550
P
i
n
=
P
o
u
t
+
P
c
+
P
i
+
P
M
+
C
.According to (7)–(13), the electric conversion energy efficiency, η
m, may be established by(14)
η
m
=
P
o
u
t
P
i
n
=
T
L
ω
r
T
L
ω
r
+
9550
k
c
T
e
2
+
k
i
+
k
w
ω
r
2
ω
r
+
C
.According to (14), which is the objective function of PSO searching, the theoretical relation of EM energy conversion efficiency to the mechanical performance can be established.


### 2.4. Transmission Characteristics of the Drivetrain

When the electric tractor works in the farm, the torque output of the EM distributes to the driving wheel and PTO and overcomes the driving and working resistance, respectively. The torque coupling relationship between the electric motor shaft, driving wheel, and PTO can be obtained by applying the following equations:(15)
T
L
ω
r
=
T
T
N
ω
T
N
η
g
η
0
+
T
w
ω
w
η
g
′
η
a
,
T
L
=
T
T
N
i
g
i
0
η
g
η
0
+
T
w
i
g
′
i
a
η
g
′
η
a
,
ω
r
=
ω
T
N
i
g
i
0
=
ω
w
i
g
′
i
a
.When the tractor plows, the force balance of the drivetrain can be described by the following equations:(16)
F
T
N
=
F
g
+
F
f
+
F
p
+
F
A
f
+
F
i
,
F
T
N
=
T
T
N
1000
r
T
N
,
F
g
=
Z
b
1
h
k
k
.Driving speed and tractive force can be described as follows:(17)
v
=
0.377
r
ω
r
η
δ
i
g
i
0
,
F
T
=
F
T
N
η
T
,
η
T
=
η
c
η
f
η
δ
.According to (15)–(17), the relationship between load and mechanical output of the EM can be expressed.


## 2.1. Introduction of the Electric Tractor Drivetrain

Although the electric tractor shares many common structural features with the common electric vehicle, there are some significant differences between the two drivetrains. Particularly, the electric vehicle layout is organized as shown in Figure1. The electric tractor drivetrain is powered by the energy system (ES). The ES can contain the storage battery in a pure electric tractor or the ICE, generator, battery, super-capacitor, and so forth in a hybrid electric tractor. The electric motor (EM) provides the mechanical power for the drive wheel and power take-off shaft (PTO). PTO joints with the EM via the clutch. When the tractor only drives the agricultural implement (AI), the clutch decouples the mechanical connection between the PTO and the drive wheel, which makes the mechanical power flows of them be independent of each other. The gearbox transmission (GT) can transform the dynamic power performance field of the drive wheel via the gear shifting operation of driver. And the GT directly transmits the mechanical power to PTO, which can make the rotate speed of PTO match with either the rotate speed of the EM or the forward velocity of the electric tractor. The main drive (MD) acts as the differential mechanism and final reducer of the transmission system. By stepping on the accelerator pedal (AP), the driver can control the EM through electrical adjustment of the motor controller (MC). Generally, there are two types of load force exerted upon the electric tractor: one is driving force and the other is tilling force.Figure 1
Scheme of an electric tractor’s drivetrain.


## 2.2. Electric-Mechanical Performance of the Electric Motor (EM)

According to the previous research, the EM is a brushless direct current motor (BLDC) [3]. The voltage per phase can be expressed by(1)
U
a
=
R
s
I
a
+
L
a
-
L
m
d
I
a
d
t
+
e
a
,
U
b
=
R
s
I
b
+
L
b
-
L
m
d
I
b
d
t
+
e
b
,
U
c
=
R
s
I
c
+
L
c
-
L
m
d
I
c
d
t
+
e
c
.After equivalent transformation, the electromechanical characteristics of BLDC can be described by the following equations [12]:(2)
U
t
=
R
s
I
s
+
L
n
-
L
m
d
I
s
d
t
+
e
s
,
L
s
=
L
n
-
L
m
,
e
s
=
k
E
ω
r
.The electromagnetic torque of BLDC,T
e, can be obtained as follows:(3)
T
e
=
9.55
e
a
I
a
+
e
b
I
b
+
e
c
I
c
ω
r
=
9.55
e
s
I
s
ω
r
,
(4)
T
e
=
k
T
I
s
.Depending on (3) and (4), the torque constant, k
T, is equal to back EMF constant, k
E, after the dimension change.When BLDC is working on traction, the electromagnetic torque of the motor can also be described by(5)
T
e
=
T
L
+
J
d
ω
r
d
t
+
V
ω
r
.Based on the Laplace transform of (1)–(5), the electric-mechanical performance of BLDC can be expressed by the following transfer function:(6)
ω
r
s
=
k
T
R
s
+
s
L
s
s
J
+
V
+
k
E
k
T
U
t
s
-
R
s
+
s
L
s
R
s
+
s
L
s
s
J
+
V
+
k
E
k
T
T
L
s
.Consequently, the steady state characteristic can be established:(7)
ω
r
=
k
T
R
s
V
+
k
E
k
T
U
t
-
R
s
R
s
V
+
k
E
k
T
T
L
.According to (4), (6), and (7), the mechanical performance of the EM can be linked to the controllable electric parameter.


## 2.3. Energy Consumption Model of the EM

When the EM is running, electric power loss contains copper loss, icon loss, mechanical loss, and constant power used for driving the MC.Copper loss,P
c, is the prime waste of the BLDC. This is due to the fact that the stator winding impedance converts the electric power to the thermal dissipated to the surroundings. It can be obtained as follows [13]: (8)
P
c
=
I
s
2
R
s
.According to (4), the value of the current is directly proportional to the value of electromagnetic torque. Thus, we can use a coefficient, k
c, to describe the relationship between the copper loss and the electromagnetic torque. The copper loss can be rewritten as follows:(9)
P
c
=
k
c
T
e
2
.The copper loss coefficient,k
c, can be calculated by(10)
k
c
=
R
s
k
T
2
.Icon loss power,P
i, mainly contains magnetic hysteresis loss and current loss, and its mechanisms are complex and changing. It mainly depends on the speed of rotor and is described by(11)
P
i
=
k
i
ω
r
.In the BLDC,k
i can be approximately seen as a constant and can be measured by an experimental method.The mechanical lossP
M, output from the rotor, is largely wasted on air resistance. It can be obtained by(12)
P
M
=
k
w
ω
r
3
.The mechanical power output and total power input of the electric motor can be described by the following equations:(13)
P
o
u
t
=
T
L
ω
r
9550
P
i
n
=
P
o
u
t
+
P
c
+
P
i
+
P
M
+
C
.According to (7)–(13), the electric conversion energy efficiency, η
m, may be established by(14)
η
m
=
P
o
u
t
P
i
n
=
T
L
ω
r
T
L
ω
r
+
9550
k
c
T
e
2
+
k
i
+
k
w
ω
r
2
ω
r
+
C
.According to (14), which is the objective function of PSO searching, the theoretical relation of EM energy conversion efficiency to the mechanical performance can be established.


## 2.4. Transmission Characteristics of the Drivetrain

When the electric tractor works in the farm, the torque output of the EM distributes to the driving wheel and PTO and overcomes the driving and working resistance, respectively. The torque coupling relationship between the electric motor shaft, driving wheel, and PTO can be obtained by applying the following equations:(15)
T
L
ω
r
=
T
T
N
ω
T
N
η
g
η
0
+
T
w
ω
w
η
g
′
η
a
,
T
L
=
T
T
N
i
g
i
0
η
g
η
0
+
T
w
i
g
′
i
a
η
g
′
η
a
,
ω
r
=
ω
T
N
i
g
i
0
=
ω
w
i
g
′
i
a
.When the tractor plows, the force balance of the drivetrain can be described by the following equations:(16)
F
T
N
=
F
g
+
F
f
+
F
p
+
F
A
f
+
F
i
,
F
T
N
=
T
T
N
1000
r
T
N
,
F
g
=
Z
b
1
h
k
k
.Driving speed and tractive force can be described as follows:(17)
v
=
0.377
r
ω
r
η
δ
i
g
i
0
,
F
T
=
F
T
N
η
T
,
η
T
=
η
c
η
f
η
δ
.According to (15)–(17), the relationship between load and mechanical output of the EM can be expressed.


## 3. Design of Load Torque Based Control Strategy (LTCS)

### 3.1. Particle Swarm Optimization Algorithm (PSO) Introduction

In comparison with classical optimization methods, PSO is an efficient range-searching algorithm. By selecting the search threshold, the population behavior can be restricted to only running under working conditions and practical requirements. It is a helpful advantage to ensure that the control parameters produced by PSO conform to the speed and torque requirement of farming. Moreover, PSO can deal with nonsmooth, noncontinuous, and nondifferentiable functions simply with objective function information. Compared with other intelligence algorithms such as the Genetic Algorithm, there is one simple operator: velocity calculation in the PSO. Using the simpler algorithm means shorter computing time, less memory, and improved robustness [14]. Therefore, PSO is chosen here as the primary algorithm of LTCS.By means of simulating bird flock preying behavior, particle swarm optimization used widely in the engineering field has the advantage of high robustness and low complexity [15]. In a searching space, the particles population, which contains n particles, can be built as X
=
(
X
1
,
X
2
,
…
,
X
n
). The ith particle is a vector with m dimensionalities and can be built as X
i
=
X
i
1
,
X
i
2
,
…
,
X
i
m
T to express the position of this particle, namely, a potential solution of the objective function. The fitness values of each X
i may be calculated by the objective function. The speed of X
i can be built as V
i
=
V
i
1
,
V
i
2
,
…
,
V
i
n
T, the individual best solution is Y
i
=
Y
i
1
,
Y
i
2
,
…
,
Y
i
n
T, and the swarm optimal solution is Y
g
=
Y
g
1
,
Y
g
2
,
…
,
Y
g
n
T. During the iterative process, every particle of the swarm uses its Y
i and Y
g to update the V
i and X
i, which can be described as follows:(18)
V
i
d
γ
+
1
=
ε
V
i
d
γ
+
c
1
R
1
Y
i
d
γ
-
X
i
d
γ
+
c
2
R
2
Y
g
d
γ
-
X
i
d
γ
,
X
i
d
γ
+
1
=
X
i
d
γ
+
V
i
d
γ
+
1
,where d
=
1
,
2
,
…
,
n; i
=
1
,
2
,
…
,
m; and R
1 and R
2 are random numbers in [
0
,
1
].


### 3.2. Working Characteristic

The 1804 series hybrid electric tractor was selected as the study object of this paper in a project designed by taking the YTO-1804 as a prototype [3]. Its drivetrain conforms to the description in Figure 1, and the speed-traction force performance curve is built in Figure 2.Figure 2
Speed-traction force performance curve of the 1804 series hybrid electric tractor.When an 1804 hybrid works under heavy load gear or moderate load gear, the relationship between the driving speed and the traction force could be described by(19)
d
F
T
d
v
>
d
v
d
F
T
.When the tractor works under the light load gear or transport gear, the relationship can be described by the following equation:(20)
d
v
d
F
T
≈
d
F
T
d
v
,
d
F
T
d
v
<
d
v
d
F
T
.Since the 1804 hybrid always tills with a stable speed and main farm works are done by using heavy, moderate, and light load gear, meeting the requirement of the traction force is more important than meeting the requirement of the speed. Additionally, the control result can match the change of traction force well, andT
L and ω
r are, respectively, defined as the independent variable and dependent variable of the PSO.


### 3.3. Control Strategy

Control function of the LTCS in the vehicle is shown in Figure3. Compared with the electric tractor in Figure 1, the LTCS is added between the driver and the MC. Two input signals, which are, respectively, acquired from the position sensor on the accelerator pedal (AP) and the shift lever, are deduced and calculated by LTCS. Afterward, based on internal rules, LTCS provides control signals to the MC and energy system (ES) to control the running status of the EM.Figure 3
Schematic representation of the control function for the electric drivetrain.When the position sensor measures the depth of AP, load torque can be calculated by(21)
T
L
=
B
-
A
·
N
A
·
1
-
N
·
T
m
a
x
.During every sampling time, the internal task actions of LTCS are built in Figure4, which includes two parts. One is the search process of the best energy-using working point, and the other is the timely control of the search process initialization.Figure 4
Schematic description of internal task actions within the LTCS.The best working point (ω
r
′
,
T
L
′) that corresponds to the best electric energy efficiency of the drive motor is selected as the searching target of the PSO searching process. Because T
L
′ should be greater than T
L, the enclosure space [
T
L
,
T
L
+
Q
] is defined as the constraint. Q is the adjustment torque.After every initialization signal, the search process runs only once. Via (4) and (6), the searched best working point was used as the control parameter of the torque based current feedback control of the MC and the DC-AC converter of the energy system (ES) to, respectively, control the expected phase current and input voltage until the next initialization occurs, which follows these rules:(1)
In every sampling step, the initialization signals sent by the time trigger and load torque controller are synchronous and mutually independent. When the initialization is triggered, the search process and following control process occurred immediately.
(2)
Using the sampling step as the period, the time trigger periodically sends the initialization signal. It is used to deal with the accumulation of smaller AP position changes, which cannot trigger the initialization signal of the load torque controller during every period.
(3)
In every period, the initialization can also be controlled by the torque controller, which judges the load torque changing status. If the load change that occurred during this period is large enough, the PSO searching is initialized immediately.


### 3.4. The Design of Load Torque Controller

The load change degree is judged by the gearshift and footplate positions, which can be measured by the gear sensor and the AP position sensor.Shifting gears is an operation greatly affected by the subjective intention of the driver. In other words, when the tractor meets visible topographic changes or the working type changes, where the load torque will likely experience a significant change, the driver usually shifts the gear. Logical algorithms are used to process the signal of the gearshift position. If the driver shifts the gear (which means the load torque condition change is large enough), the controller sends the initialization signal.Under any gear, the driver can estimate the load force situation by viewing the tachometer reading. When the AP position is stable, changing the tachometer reading means that the force balance on the EM shaft has been broken. The driver should adjust the AP depth based on the change trend in the tachometer reading until the tractor speed is similar to what it was previously. A fuzzy controller is designed to deal with the AP signal, inferring the change degree of working condition and providing the initialization signal.Wang et al. discussed the fuzzy rules that reflect the relationship between diver intention and AP operation [16]. Based on that, two of the fuzzy rules used for deducing the initialization of LTCS are shown here in Table 1.Table 1
Fuzzy rules to infer initialization.
Deeping rate
Footplate depth
S
M
B
−B
INIT
INIT
INIT
−M
HOLD
INIT
INIT
S
HOLD
HOLD
HOLD
+M
HOLD
INIT
INIT
+B
INIT
INIT
INITThe membership functions that belong to the footplate depth, the footplate deepening rate, and the collection mode are described by Figures5–7.Figure 5
Membership functions of the footplate depth.Figure 6
Membership functions of the footplate deepening rate.Figure 7
Membership functions of the collection mode.In the case of the footplate depth, it has three chosen membership functions: S, small depth; M, moderate depth; and B, big depth. They cover the full range of footplate depth. When the AP controls BLDC, its control function is equivalent to the voltage control process of the potentiometer. For the torque control of BLDC, the torque controlled by AP is normally proportional to the expected electric current, so the footplate depth membership functions are equally distributed. To transform the real domain into fuzzy domain scope[
0
,
1
], the quantization factor of footplate depth, ε
d, is calculated by(22)
ε
d
=
A
1
-
N
.In the case of the footplate deepening rate, five membership functions are considered: S, small change rate; +M, moderate deepening rate; −M, moderate returning rate; +B, big deepening rate; and −B, big returning rate. When the tractor speed increases with an unchanged AP position, which means the load torque reduces, the driver may return the pedal until the speed changes back and vice versa. Thus, the formulation of footplate deepening rate membership functions should contain this case. The quantization factor of footplate deepening rate,ε
d
r, is calculated by(23)
ε
d
r
=
A
1
-
N
t
r
.Depending on the test of McGehee et al. [17], t
r can be set as 1.28 s.Regarding the output variables, the collection mode contains two membership functions: INIT, the initialization of LTCS, and HOLD, nothing else changing in the load torque control process. Because the opposition between initialization and hold operations is definite and complete, their fuzzy domains are complementary.Surface of the fuzzy rules can be viewed in Figure8.Figure 8
Fuzzy surface.


## 3.1. Particle Swarm Optimization Algorithm (PSO) Introduction

In comparison with classical optimization methods, PSO is an efficient range-searching algorithm. By selecting the search threshold, the population behavior can be restricted to only running under working conditions and practical requirements. It is a helpful advantage to ensure that the control parameters produced by PSO conform to the speed and torque requirement of farming. Moreover, PSO can deal with nonsmooth, noncontinuous, and nondifferentiable functions simply with objective function information. Compared with other intelligence algorithms such as the Genetic Algorithm, there is one simple operator: velocity calculation in the PSO. Using the simpler algorithm means shorter computing time, less memory, and improved robustness [14]. Therefore, PSO is chosen here as the primary algorithm of LTCS.By means of simulating bird flock preying behavior, particle swarm optimization used widely in the engineering field has the advantage of high robustness and low complexity [15]. In a searching space, the particles population, which contains n particles, can be built as X
=
(
X
1
,
X
2
,
…
,
X
n
). The ith particle is a vector with m dimensionalities and can be built as X
i
=
X
i
1
,
X
i
2
,
…
,
X
i
m
T to express the position of this particle, namely, a potential solution of the objective function. The fitness values of each X
i may be calculated by the objective function. The speed of X
i can be built as V
i
=
V
i
1
,
V
i
2
,
…
,
V
i
n
T, the individual best solution is Y
i
=
Y
i
1
,
Y
i
2
,
…
,
Y
i
n
T, and the swarm optimal solution is Y
g
=
Y
g
1
,
Y
g
2
,
…
,
Y
g
n
T. During the iterative process, every particle of the swarm uses its Y
i and Y
g to update the V
i and X
i, which can be described as follows:(18)
V
i
d
γ
+
1
=
ε
V
i
d
γ
+
c
1
R
1
Y
i
d
γ
-
X
i
d
γ
+
c
2
R
2
Y
g
d
γ
-
X
i
d
γ
,
X
i
d
γ
+
1
=
X
i
d
γ
+
V
i
d
γ
+
1
,where d
=
1
,
2
,
…
,
n; i
=
1
,
2
,
…
,
m; and R
1 and R
2 are random numbers in [
0
,
1
].


## 3.2. Working Characteristic

The 1804 series hybrid electric tractor was selected as the study object of this paper in a project designed by taking the YTO-1804 as a prototype [3]. Its drivetrain conforms to the description in Figure 1, and the speed-traction force performance curve is built in Figure 2.Figure 2
Speed-traction force performance curve of the 1804 series hybrid electric tractor.When an 1804 hybrid works under heavy load gear or moderate load gear, the relationship between the driving speed and the traction force could be described by(19)
d
F
T
d
v
>
d
v
d
F
T
.When the tractor works under the light load gear or transport gear, the relationship can be described by the following equation:(20)
d
v
d
F
T
≈
d
F
T
d
v
,
d
F
T
d
v
<
d
v
d
F
T
.Since the 1804 hybrid always tills with a stable speed and main farm works are done by using heavy, moderate, and light load gear, meeting the requirement of the traction force is more important than meeting the requirement of the speed. Additionally, the control result can match the change of traction force well, andT
L and ω
r are, respectively, defined as the independent variable and dependent variable of the PSO.


## 3.3. Control Strategy

Control function of the LTCS in the vehicle is shown in Figure3. Compared with the electric tractor in Figure 1, the LTCS is added between the driver and the MC. Two input signals, which are, respectively, acquired from the position sensor on the accelerator pedal (AP) and the shift lever, are deduced and calculated by LTCS. Afterward, based on internal rules, LTCS provides control signals to the MC and energy system (ES) to control the running status of the EM.Figure 3
Schematic representation of the control function for the electric drivetrain.When the position sensor measures the depth of AP, load torque can be calculated by(21)
T
L
=
B
-
A
·
N
A
·
1
-
N
·
T
m
a
x
.During every sampling time, the internal task actions of LTCS are built in Figure4, which includes two parts. One is the search process of the best energy-using working point, and the other is the timely control of the search process initialization.Figure 4
Schematic description of internal task actions within the LTCS.The best working point (ω
r
′
,
T
L
′) that corresponds to the best electric energy efficiency of the drive motor is selected as the searching target of the PSO searching process. Because T
L
′ should be greater than T
L, the enclosure space [
T
L
,
T
L
+
Q
] is defined as the constraint. Q is the adjustment torque.After every initialization signal, the search process runs only once. Via (4) and (6), the searched best working point was used as the control parameter of the torque based current feedback control of the MC and the DC-AC converter of the energy system (ES) to, respectively, control the expected phase current and input voltage until the next initialization occurs, which follows these rules:(1)
In every sampling step, the initialization signals sent by the time trigger and load torque controller are synchronous and mutually independent. When the initialization is triggered, the search process and following control process occurred immediately.
(2)
Using the sampling step as the period, the time trigger periodically sends the initialization signal. It is used to deal with the accumulation of smaller AP position changes, which cannot trigger the initialization signal of the load torque controller during every period.
(3)
In every period, the initialization can also be controlled by the torque controller, which judges the load torque changing status. If the load change that occurred during this period is large enough, the PSO searching is initialized immediately.


## 3.4. The Design of Load Torque Controller

The load change degree is judged by the gearshift and footplate positions, which can be measured by the gear sensor and the AP position sensor.Shifting gears is an operation greatly affected by the subjective intention of the driver. In other words, when the tractor meets visible topographic changes or the working type changes, where the load torque will likely experience a significant change, the driver usually shifts the gear. Logical algorithms are used to process the signal of the gearshift position. If the driver shifts the gear (which means the load torque condition change is large enough), the controller sends the initialization signal.Under any gear, the driver can estimate the load force situation by viewing the tachometer reading. When the AP position is stable, changing the tachometer reading means that the force balance on the EM shaft has been broken. The driver should adjust the AP depth based on the change trend in the tachometer reading until the tractor speed is similar to what it was previously. A fuzzy controller is designed to deal with the AP signal, inferring the change degree of working condition and providing the initialization signal.Wang et al. discussed the fuzzy rules that reflect the relationship between diver intention and AP operation [16]. Based on that, two of the fuzzy rules used for deducing the initialization of LTCS are shown here in Table 1.Table 1
Fuzzy rules to infer initialization.
Deeping rate
Footplate depth
S
M
B
−B
INIT
INIT
INIT
−M
HOLD
INIT
INIT
S
HOLD
HOLD
HOLD
+M
HOLD
INIT
INIT
+B
INIT
INIT
INITThe membership functions that belong to the footplate depth, the footplate deepening rate, and the collection mode are described by Figures5–7.Figure 5
Membership functions of the footplate depth.Figure 6
Membership functions of the footplate deepening rate.Figure 7
Membership functions of the collection mode.In the case of the footplate depth, it has three chosen membership functions: S, small depth; M, moderate depth; and B, big depth. They cover the full range of footplate depth. When the AP controls BLDC, its control function is equivalent to the voltage control process of the potentiometer. For the torque control of BLDC, the torque controlled by AP is normally proportional to the expected electric current, so the footplate depth membership functions are equally distributed. To transform the real domain into fuzzy domain scope[
0
,
1
], the quantization factor of footplate depth, ε
d, is calculated by(22)
ε
d
=
A
1
-
N
.In the case of the footplate deepening rate, five membership functions are considered: S, small change rate; +M, moderate deepening rate; −M, moderate returning rate; +B, big deepening rate; and −B, big returning rate. When the tractor speed increases with an unchanged AP position, which means the load torque reduces, the driver may return the pedal until the speed changes back and vice versa. Thus, the formulation of footplate deepening rate membership functions should contain this case. The quantization factor of footplate deepening rate,ε
d
r, is calculated by(23)
ε
d
r
=
A
1
-
N
t
r
.Depending on the test of McGehee et al. [17], t
r can be set as 1.28 s.Regarding the output variables, the collection mode contains two membership functions: INIT, the initialization of LTCS, and HOLD, nothing else changing in the load torque control process. Because the opposition between initialization and hold operations is definite and complete, their fuzzy domains are complementary.Surface of the fuzzy rules can be viewed in Figure8.Figure 8
Fuzzy surface.


## 4. Simulation of the Traction Experiment

### 4.1. Parameters and Conditions of the Simulation

The simulation of the traction experiment is designed as a backward type. Working conditions, which act as the input of the simulation model, are required to contain all the common working resistances of the 1804 hybrid tractor. According to the common soil conditions of the Northeast China region, the rolling resistance coefficient,f, is 0.1, the adhesion coefficient is 0.7, and the clayed ground, whose soil specific resistances randomly changed from 0.007 to 0.01, is tilled. Parameters of the plowing resistances, which imitate resistances of common working conditions, are expressed in Table 2.Table 2
Plowing parameters.
Number
Z
b
l
h
k
0
7
35
30
1
6
35
30
2
5
32
25
3
4
35
25
4
4
40
31
5
3
30
20
6
3
35
22
7
1
50
20The working condition in the simulation can be described by Figure9 (see data in Supplementary Material A in Supplementary Material available online at http://dx.doi.org/10.1155/2016/2548967). The equivalent resistance ranged from 5 kN to 62.5 kN, which agrees with the tractive force scope in Figure 2.Figure 9
Tractive force of the simulation. Note: numbering 0–7 corresponds to the numbers of plowing parameters in Table2.According to our previous research [3], Table 3 shows the simulative parameters of the 1804 hybrid drivetrain, which acts as the controlled object.Table 3
Simulative parameters of the electric drivetrain.
Part
Parameter name
Value and unit
Vehicle
Overall length
5.2 m
Overall width
2.69 m
Overall height
2.97 m
L
2.85 m
h
T
1.6 m
G
112.64 kN
r
T
N
0.9 m
a
1.97 m
GT
Heavy load gear ratio
5.38
Moderate load gear ratio
3.88
Light load gear ratio
2.63
Transport gear ratio
1.65
η
g
96.04%
EM
Voltage rating
380 V
Peak voltage
540 V
Torque rating
4138 N·m
Rotate speed rating
300 r/min
Reduction ratio
23
P
r
130 kW
R
s
0.04 Ω
J
0.0001 kg·m2
V
0.0368 N·m·min/r
k
c
0.15 kW/N2·m2
k
i
0.1 kW·min/r
k
w
1.2e
-
7 kW·min3/r3
k
E
0.0543 V·min/r
k
T
0.5186 N·m/A
C
0.21 kW
MD
i
0
3.2
η
0
98%
AP
A
54°
N
0.15Based on mass experimental data, Zhou et al. presented the theoretical models used to predict theoretical tractive performance of a wheeled tractor [18]. According to this model, the slip rate, δ, can be calculated by the following equation:(24)
δ
=
δ
∗
L
n
ϕ
m
a
x
ϕ
m
a
x
-
L
F
T
N
/
a
G
+
h
T
F
T
N
.In (24), δ
∗ and ϕ
m
a
x are fitting coefficients whose values are recommended as 0.0757 and 0.624 under the common soil conditions of the Northeast China region. The relationship between η
δ and δ is η
δ
=
1
-
δ. The slip efficiency used in the simulation is established in Figure 10.Figure 10
Slip efficiency in the simulation.In the simulation of plowing, the clutch between GT and PTO is separated. Therefore, the mechanical loss of the agricultural implement (AI) is irrelevant. The rest of the parameters of efficiency are calculated as(25)
η
c
=
η
0
η
g
,
η
f
=
G
f
F
T
N
.Accordingly, the transmission efficiency,η
c, is equal to 94.12%, and the rolling efficiency, η
f, is dynamically equal to 11.2/F
T
N.Based on (22) and (23), ε
d and ε
d
r of the load torque controller are, respectively, assigned as 44.8 and 35.86.In the simulation, the tractor working system is simplified as a longitudinal model, which tills in a straight line, and all the side forces can be disregarded. A driver response time to the peak force is set as 0.1 s. It is assumed that the driver can accurately adjust the AP depth in accordance with the load force change. The sampling step length is 5 s. The load force fluctuation owed to the AI and gear switching is ignored.


### 4.2. PSO Setting in the Simulation

In the simulation, the PSO searching process runs on the platform developed by Birge [19]. Shi and Eberhart researched the set of inertia weights by the simulation on the benchmark problems and found that inertia weight starting with a value from 0.9 to 0.4 through the course of the run greatly improved the performance of PSO [20]. Thus, the start weight value selected here is 0.9 and the end weight value is 0.4.The determination of maximum velocityV
i
m
a
x, c
1, and c
2 is mainly based on the formula presented by Zhang et al. as follows [21]:(26)
V
i
m
a
x
=
λ
X
i
m
a
x
,
c
1
=
c
2
=
φ
2
.In this case, the searching space of PSO is 2-dimensional and the objective function, (14), is unimodal. Therefore, it is recommended that λ is equal to 0.05 and φ is set as 4.1.As stated earlier, the PSO is searching for the best working point in constraint[
T
L
,
T
L
+
Q
]. The adjustment torque Q is set as 0.5 kN for the purpose of promoting of the search accuracy. In addition, in order to prevent the mechanical power of the searched working point from exceeding the rating power of the EM, the rotate speed of the EM should be bounded by the range, X
2
m
a
x. Thus, X
i
m
a
x is described by(27)
X
1
m
a
x
=
R
T
=
T
L
,
T
L
+
0.5
,
X
2
m
a
x
=
R
M
A
V
=
0
,
9550
P
r
T
L
.In every control step, the value ofT
L is dependent upon (21).Shi and Eberhart advised that probably the most common population sizes of PSO are from 20 to 50 [20]. Experimental verification has found that there is no visible difference in searching results when the population size is, respectively, set as 20, 30, 40, and 50. To simplify, we set the population size at 20.


### 4.3. Analysis of the Simulation Result

Figure11 depicts the searching result when the tractive force is 55 kN (which belongs to the heavy load scope). After 38 evolutions, the population fitness degree becomes stable. The best working point is searched as 256.339 r/min, 3676.97 N·m, and the electric energy conversion efficiency is 92.38%.Figure 11
PSO searching on the heavy load situation.Figure12 depicts the searching result when the load torque is 14 kN (which belongs to the light load scope). After 20 evolutions, the population fitness degree becomes stable. The best working point was searched as 218.8 r/min, 2864 N·m, and the electric energy conversion efficiency was 92.83%. The searching results on two situations show that the LTCS makes the result of PSO searching self-adaptive for the tractive force. In addition, the convergent evolution times testify that the set of PSO parameters make the searching process sufficiently rapid.Figure 12
PSO searching on the light load situation.The control single of the load torque controller in the simulation is expressed by Figure13. In comparison with Figure 9, when the hybrid works with heavy, moderate, and light load gear, initialization signals can frequently be activated by strong fluctuations in the load force. After 680 s, when the load force begins to meet the transport condition with a relatively stable load force, there is no signal output. In conclusion, the rules of the load torque controller can judge the load force with acceptable accuracy and adaptivity. The initialization signals’ superposition process during the simulation is described in Supplementary Material B.Figure 13
Output of load torque controller in the simulation.Figure14 depicts the simulative results of the voltage-current performance controlled by LTCS. According to the discussion of Figure 4, the current, I
s, is defined as T
L
′, which follows the relationship of (4). And the rotate speed of EM, ω
r, is adjusted by changing the voltage.Figure 14
The electrical characteristic controlled in the simulation.The load torque followed control method (LTFC) is used popularly in the BLDC control of the electric tractor. According to the excepted torque, LTFC controls the mechanical output of the EM by adjusting the phase current [12]. Compared with another electric vehicle electric motor control strategy, the speed followed control method (which controls the BLDC with expected speed); according to the discussion in Figure 2, LTFC is more suitable for the farming working conditions, which have high-level load force. Thus, we selected the LTFC as the comparison for analyzing the control performance of LTCS.Figure15 depicts the drivetrain working points of LTCS when the tractor works with heavy loads. The figure shows that some of working points are distributed on the edge of the speed-traction force performance curve, and others are in the central area near 92.5%. Both of the two distributions are near the best efficiency area under an equivalent traction force. The speed range of working points is 4.1 km/h~5.8 km/h, which can meet the speed scope of plowing.Figure 15
Working points distribution of the high load gear.By comparing the heavy load control performance of LTCS and LTFC in Figure15, all of the working points of both control strategies can meet the speed requirement of plowing, but the working points of LTCS present more efficient energy conversion.Figure16 depicts the drivetrain working points when the tractor works with moderate load. The energy conversion efficiency of the working points controlled by LTCS exceeds 92.1%. When the traction force is less than 33.3 kN, the efficiency can be greater than 92.5%. In addition, the speed range is 5.1 km/h~7.4 km/h, which meets the speed requirement of cultivating and hoeing. Compared with the case of LTFC, LTCS presents more efficient speed and control performance.Figure 16
Working points distribution of moderate load gear.Figure17 depicts the working points’ distribution when the tractor works with a light load. As shown, the energy conversion efficiency of EM controlled by LTCS is almost always superior. The speed range is 8.5 km/h~10.7 km/h, which can meet the speed scope of sowing, harrowing, and stubble chopping. In addition, the maximum speed change rate is 0.2 km/kN·h, which can enhance the seeding uniformity.Figure 17
Working points distribution of light load gear.By contrast, the speeds of all of the working points controlled by LTFC are higher than the speed range of sowing, harrowing, and stubble chopping. This is because when the tractor is working with AI, the rotate speed of PTO is commonly changeless. In the 1804 series hybrid electric tractor, the rotate speed is 1000 r/min or 540 r/min. When the tractor speed exceeds the speed scope of its farming type, the farming quality may be decreased (e.g., for the sowing, a faster tractor speed means longer seed spacing, which will reduce the seeding rate and consequently decrease the crop production). Therefore, LTFC is not suitable for the light working conditions of electric tractor.Figure18 shows the distribution of working points when LTCS is used for the transport condition. Although the efficiencies of all of the working points can be approximated to the maximum under an equivalent traction force, the speed is relatively slow for transporting, which may reduce the conveying efficiency. Furthermore, when the 1804 hybrid transports out of the farm, the control result cannot attain adequate traffic security and drive flexibility, which means LTCS is not appropriate for the road transport condition.Figure 18
Working points distribution of transport gear.In the design of LTCS, load torque is used as the independent variable of target function, and the enclosure space based on load torque is the only constraint of the PSO search. Consequently, the speed changes dependent on the status of the load force, and that is why the speed cannot be adjusted independently and intentionally.Simulation results of energy conversion efficiency are depicted in Figure19. In each simulation, the electric tractor controlled by LTCS is more efficient than when controlled by LTFC in electric use. In addition, its mean energy conversion efficiency is 92.46%, and there are small improvements in regions A and B. The reason for this is that further improvement of energy conversion efficiency is limited by the rating power of electric motor, which proves that the set of PSO searching ranges are appropriate.Figure 19
The comparison of energy conversion efficiency under the control of LTCS and LTFC.Figure20 depicts the electric power consumption of the ES. Due to the improvement of the EM energy conversion efficiency and the speed-limited function of LTCS, the mean electric power consumption can reduce by 23 kW, and the energy-saving effect of farming work can increase by 21.3% in comparison with the LTFC.Figure 20
The comparison of ES power consumption under the control of LTCS and LTFC.In summary, when LTCS is used for working the farm, the electric tractor possesses superior agricultural adaptation and an obvious energy-saving effect. In off-road conditions, the electric tractor can also be engaged in some transport operations where speed performance is not required. In the future, when the LTCS is used in a realistic environment, the following is recommended:(1)
The objective function should be further amended by considering the change of temperature or be built by experimental method.
(2)
The mechanical dynamic characteristics of the motor output shaft should be controlled by the PSO method presented by Elwer et al. [10].
(3)
The control parameters of PSO should be optimized by considering the hardware performance of the electric control unit selected; the principle is based on the premise of appropriate precision, so the convergence rate should be increased to the maximum speed.


## 4.1. Parameters and Conditions of the Simulation

The simulation of the traction experiment is designed as a backward type. Working conditions, which act as the input of the simulation model, are required to contain all the common working resistances of the 1804 hybrid tractor. According to the common soil conditions of the Northeast China region, the rolling resistance coefficient,f, is 0.1, the adhesion coefficient is 0.7, and the clayed ground, whose soil specific resistances randomly changed from 0.007 to 0.01, is tilled. Parameters of the plowing resistances, which imitate resistances of common working conditions, are expressed in Table 2.Table 2
Plowing parameters.
Number
Z
b
l
h
k
0
7
35
30
1
6
35
30
2
5
32
25
3
4
35
25
4
4
40
31
5
3
30
20
6
3
35
22
7
1
50
20The working condition in the simulation can be described by Figure9 (see data in Supplementary Material A in Supplementary Material available online at http://dx.doi.org/10.1155/2016/2548967). The equivalent resistance ranged from 5 kN to 62.5 kN, which agrees with the tractive force scope in Figure 2.Figure 9
Tractive force of the simulation. Note: numbering 0–7 corresponds to the numbers of plowing parameters in Table2.According to our previous research [3], Table 3 shows the simulative parameters of the 1804 hybrid drivetrain, which acts as the controlled object.Table 3
Simulative parameters of the electric drivetrain.
Part
Parameter name
Value and unit
Vehicle
Overall length
5.2 m
Overall width
2.69 m
Overall height
2.97 m
L
2.85 m
h
T
1.6 m
G
112.64 kN
r
T
N
0.9 m
a
1.97 m
GT
Heavy load gear ratio
5.38
Moderate load gear ratio
3.88
Light load gear ratio
2.63
Transport gear ratio
1.65
η
g
96.04%
EM
Voltage rating
380 V
Peak voltage
540 V
Torque rating
4138 N·m
Rotate speed rating
300 r/min
Reduction ratio
23
P
r
130 kW
R
s
0.04 Ω
J
0.0001 kg·m2
V
0.0368 N·m·min/r
k
c
0.15 kW/N2·m2
k
i
0.1 kW·min/r
k
w
1.2e
-
7 kW·min3/r3
k
E
0.0543 V·min/r
k
T
0.5186 N·m/A
C
0.21 kW
MD
i
0
3.2
η
0
98%
AP
A
54°
N
0.15Based on mass experimental data, Zhou et al. presented the theoretical models used to predict theoretical tractive performance of a wheeled tractor [18]. According to this model, the slip rate, δ, can be calculated by the following equation:(24)
δ
=
δ
∗
L
n
ϕ
m
a
x
ϕ
m
a
x
-
L
F
T
N
/
a
G
+
h
T
F
T
N
.In (24), δ
∗ and ϕ
m
a
x are fitting coefficients whose values are recommended as 0.0757 and 0.624 under the common soil conditions of the Northeast China region. The relationship between η
δ and δ is η
δ
=
1
-
δ. The slip efficiency used in the simulation is established in Figure 10.Figure 10
Slip efficiency in the simulation.In the simulation of plowing, the clutch between GT and PTO is separated. Therefore, the mechanical loss of the agricultural implement (AI) is irrelevant. The rest of the parameters of efficiency are calculated as(25)
η
c
=
η
0
η
g
,
η
f
=
G
f
F
T
N
.Accordingly, the transmission efficiency,η
c, is equal to 94.12%, and the rolling efficiency, η
f, is dynamically equal to 11.2/F
T
N.Based on (22) and (23), ε
d and ε
d
r of the load torque controller are, respectively, assigned as 44.8 and 35.86.In the simulation, the tractor working system is simplified as a longitudinal model, which tills in a straight line, and all the side forces can be disregarded. A driver response time to the peak force is set as 0.1 s. It is assumed that the driver can accurately adjust the AP depth in accordance with the load force change. The sampling step length is 5 s. The load force fluctuation owed to the AI and gear switching is ignored.


## 4.2. PSO Setting in the Simulation

In the simulation, the PSO searching process runs on the platform developed by Birge [19]. Shi and Eberhart researched the set of inertia weights by the simulation on the benchmark problems and found that inertia weight starting with a value from 0.9 to 0.4 through the course of the run greatly improved the performance of PSO [20]. Thus, the start weight value selected here is 0.9 and the end weight value is 0.4.The determination of maximum velocityV
i
m
a
x, c
1, and c
2 is mainly based on the formula presented by Zhang et al. as follows [21]:(26)
V
i
m
a
x
=
λ
X
i
m
a
x
,
c
1
=
c
2
=
φ
2
.In this case, the searching space of PSO is 2-dimensional and the objective function, (14), is unimodal. Therefore, it is recommended that λ is equal to 0.05 and φ is set as 4.1.As stated earlier, the PSO is searching for the best working point in constraint[
T
L
,
T
L
+
Q
]. The adjustment torque Q is set as 0.5 kN for the purpose of promoting of the search accuracy. In addition, in order to prevent the mechanical power of the searched working point from exceeding the rating power of the EM, the rotate speed of the EM should be bounded by the range, X
2
m
a
x. Thus, X
i
m
a
x is described by(27)
X
1
m
a
x
=
R
T
=
T
L
,
T
L
+
0.5
,
X
2
m
a
x
=
R
M
A
V
=
0
,
9550
P
r
T
L
.In every control step, the value ofT
L is dependent upon (21).Shi and Eberhart advised that probably the most common population sizes of PSO are from 20 to 50 [20]. Experimental verification has found that there is no visible difference in searching results when the population size is, respectively, set as 20, 30, 40, and 50. To simplify, we set the population size at 20.


## 4.3. Analysis of the Simulation Result

Figure11 depicts the searching result when the tractive force is 55 kN (which belongs to the heavy load scope). After 38 evolutions, the population fitness degree becomes stable. The best working point is searched as 256.339 r/min, 3676.97 N·m, and the electric energy conversion efficiency is 92.38%.Figure 11
PSO searching on the heavy load situation.Figure12 depicts the searching result when the load torque is 14 kN (which belongs to the light load scope). After 20 evolutions, the population fitness degree becomes stable. The best working point was searched as 218.8 r/min, 2864 N·m, and the electric energy conversion efficiency was 92.83%. The searching results on two situations show that the LTCS makes the result of PSO searching self-adaptive for the tractive force. In addition, the convergent evolution times testify that the set of PSO parameters make the searching process sufficiently rapid.Figure 12
PSO searching on the light load situation.The control single of the load torque controller in the simulation is expressed by Figure13. In comparison with Figure 9, when the hybrid works with heavy, moderate, and light load gear, initialization signals can frequently be activated by strong fluctuations in the load force. After 680 s, when the load force begins to meet the transport condition with a relatively stable load force, there is no signal output. In conclusion, the rules of the load torque controller can judge the load force with acceptable accuracy and adaptivity. The initialization signals’ superposition process during the simulation is described in Supplementary Material B.Figure 13
Output of load torque controller in the simulation.Figure14 depicts the simulative results of the voltage-current performance controlled by LTCS. According to the discussion of Figure 4, the current, I
s, is defined as T
L
′, which follows the relationship of (4). And the rotate speed of EM, ω
r, is adjusted by changing the voltage.Figure 14
The electrical characteristic controlled in the simulation.The load torque followed control method (LTFC) is used popularly in the BLDC control of the electric tractor. According to the excepted torque, LTFC controls the mechanical output of the EM by adjusting the phase current [12]. Compared with another electric vehicle electric motor control strategy, the speed followed control method (which controls the BLDC with expected speed); according to the discussion in Figure 2, LTFC is more suitable for the farming working conditions, which have high-level load force. Thus, we selected the LTFC as the comparison for analyzing the control performance of LTCS.Figure15 depicts the drivetrain working points of LTCS when the tractor works with heavy loads. The figure shows that some of working points are distributed on the edge of the speed-traction force performance curve, and others are in the central area near 92.5%. Both of the two distributions are near the best efficiency area under an equivalent traction force. The speed range of working points is 4.1 km/h~5.8 km/h, which can meet the speed scope of plowing.Figure 15
Working points distribution of the high load gear.By comparing the heavy load control performance of LTCS and LTFC in Figure15, all of the working points of both control strategies can meet the speed requirement of plowing, but the working points of LTCS present more efficient energy conversion.Figure16 depicts the drivetrain working points when the tractor works with moderate load. The energy conversion efficiency of the working points controlled by LTCS exceeds 92.1%. When the traction force is less than 33.3 kN, the efficiency can be greater than 92.5%. In addition, the speed range is 5.1 km/h~7.4 km/h, which meets the speed requirement of cultivating and hoeing. Compared with the case of LTFC, LTCS presents more efficient speed and control performance.Figure 16
Working points distribution of moderate load gear.Figure17 depicts the working points’ distribution when the tractor works with a light load. As shown, the energy conversion efficiency of EM controlled by LTCS is almost always superior. The speed range is 8.5 km/h~10.7 km/h, which can meet the speed scope of sowing, harrowing, and stubble chopping. In addition, the maximum speed change rate is 0.2 km/kN·h, which can enhance the seeding uniformity.Figure 17
Working points distribution of light load gear.By contrast, the speeds of all of the working points controlled by LTFC are higher than the speed range of sowing, harrowing, and stubble chopping. This is because when the tractor is working with AI, the rotate speed of PTO is commonly changeless. In the 1804 series hybrid electric tractor, the rotate speed is 1000 r/min or 540 r/min. When the tractor speed exceeds the speed scope of its farming type, the farming quality may be decreased (e.g., for the sowing, a faster tractor speed means longer seed spacing, which will reduce the seeding rate and consequently decrease the crop production). Therefore, LTFC is not suitable for the light working conditions of electric tractor.Figure18 shows the distribution of working points when LTCS is used for the transport condition. Although the efficiencies of all of the working points can be approximated to the maximum under an equivalent traction force, the speed is relatively slow for transporting, which may reduce the conveying efficiency. Furthermore, when the 1804 hybrid transports out of the farm, the control result cannot attain adequate traffic security and drive flexibility, which means LTCS is not appropriate for the road transport condition.Figure 18
Working points distribution of transport gear.In the design of LTCS, load torque is used as the independent variable of target function, and the enclosure space based on load torque is the only constraint of the PSO search. Consequently, the speed changes dependent on the status of the load force, and that is why the speed cannot be adjusted independently and intentionally.Simulation results of energy conversion efficiency are depicted in Figure19. In each simulation, the electric tractor controlled by LTCS is more efficient than when controlled by LTFC in electric use. In addition, its mean energy conversion efficiency is 92.46%, and there are small improvements in regions A and B. The reason for this is that further improvement of energy conversion efficiency is limited by the rating power of electric motor, which proves that the set of PSO searching ranges are appropriate.Figure 19
The comparison of energy conversion efficiency under the control of LTCS and LTFC.Figure20 depicts the electric power consumption of the ES. Due to the improvement of the EM energy conversion efficiency and the speed-limited function of LTCS, the mean electric power consumption can reduce by 23 kW, and the energy-saving effect of farming work can increase by 21.3% in comparison with the LTFC.Figure 20
The comparison of ES power consumption under the control of LTCS and LTFC.In summary, when LTCS is used for working the farm, the electric tractor possesses superior agricultural adaptation and an obvious energy-saving effect. In off-road conditions, the electric tractor can also be engaged in some transport operations where speed performance is not required. In the future, when the LTCS is used in a realistic environment, the following is recommended:(1)
The objective function should be further amended by considering the change of temperature or be built by experimental method.
(2)
The mechanical dynamic characteristics of the motor output shaft should be controlled by the PSO method presented by Elwer et al. [10].
(3)
The control parameters of PSO should be optimized by considering the hardware performance of the electric control unit selected; the principle is based on the premise of appropriate precision, so the convergence rate should be increased to the maximum speed.


## 5. Conclusion

In this paper, the load torque based control strategy used for maintaining high energy conversion efficiency of the traction motor in an electric tractor is presented and studied. The control process contains a PSO search that seeks the best working point with maximum energy conversion efficiency and initiation control and whose target is to adjust the searching results dynamically based on peak torque. First, by the mathematical modeling, the relation between the mechanical characteristic and energy conversion is established, and then the target function is performed. Second, by analyzing the speed-traction force performance of the 1804 hybrid, the load torque is selected as the constraint of the PSO search. By describing the schematic representation and the internal task action, the control method and flow are obtained. By making logic and fuzzy rules, the torque controller based on driver intention was designed. Finally, after the working condition collection, which covers most of the tillage force, simulation of traction experiment is conducted. The distribution of initialization signal output from the load torque controller fits the load force fluctuation, which means the initialization control can be self-adaptive to the working condition, and the search results can be time-sensitive. Compared with LTFC, LTCS maintains the energy conversion efficiency of the EM near the maximum during the simulation time. Meanwhile, the working speed is sufficiently stable for the requirement of agriculture uniformity and suitable for the speed scope of tillage.



---

*Source: 2548967-2016-04-28.xml*