#### 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.

#### 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 Figure 1. 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.

##### 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

After equivalent transformation, the electromechanical characteristics of BLDC can be described by the following equations [12]:

The electromagnetic torque of BLDC, , can be obtained as follows:

Depending on (3) and (4), the torque constant, , is equal to back EMF constant, , after the dimension change.

When BLDC is working on traction, the electromagnetic torque of the motor can also be described by

Based on the Laplace transform of (1)–(5), the electric-mechanical performance of BLDC can be expressed by the following transfer function:

Consequently, the steady state characteristic can be established:

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, , 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]:

According to (4), the value of the current is directly proportional to the value of electromagnetic torque. Thus, we can use a coefficient, , to describe the relationship between the copper loss and the electromagnetic torque. The copper loss can be rewritten as follows:

The copper loss coefficient, , can be calculated by

Icon loss power, , 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

In the BLDC, can be approximately seen as a constant and can be measured by an experimental method.

The mechanical loss , output from the rotor, is largely wasted on air resistance. It can be obtained by

The mechanical power output and total power input of the electric motor can be described by the following equations:

According to (7)–(13), the electric conversion energy efficiency, , may be established by

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:

When the tractor plows, the force balance of the drivetrain can be described by the following equations:

Driving speed and tractive force can be described as follows:

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 particles, can be built as . The th particle is a vector with dimensionalities and can be built as to express the position of this particle, namely, a potential solution of the objective function. The fitness values of each may be calculated by the objective function. The speed of can be built as , the individual best solution is , and the swarm optimal solution is . During the iterative process, every particle of the swarm uses its and to update the and , which can be described as follows:where ; ; and and are random numbers in .

##### 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.

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

When the tractor works under the light load gear or transport gear, the relationship can be described by the following equation:

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, and and 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 Figure 3. 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.

When the position sensor measures the depth of AP, load torque can be calculated by

During every sampling time, the internal task actions of LTCS are built in Figure 4, 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.

The best working point () that corresponds to the best electric energy efficiency of the drive motor is selected as the searching target of the PSO searching process. Because should be greater than , the enclosure space is defined as the constraint. 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.

The membership functions that belong to the footplate depth, the footplate deepening rate, and the collection mode are described by Figures 5–7.

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 , the quantization factor of footplate depth, , is calculated by

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, , is calculated by

Depending on the test of McGehee et al. [17], 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 Figure 8.

#### 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, , 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.

The working condition in the simulation can be described by Figure 9 (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.

According to our previous research [3], Table 3 shows the simulative parameters of the 1804 hybrid drivetrain, which acts as the controlled object.

Based 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:

In (24), and 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 . The slip efficiency used in the simulation is established in Figure 10.

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

Accordingly, the transmission efficiency, , is equal to 94.12%, and the rolling efficiency, , is dynamically equal to 11.2/.

Based on (22) and (23), and 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 velocity , , and is mainly based on the formula presented by Zhang et al. as follows [21]:

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 . The adjustment torque 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, . Thus, is described by

In every control step, the value of 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

Figure 11 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 Nm, and the electric energy conversion efficiency is 92.38%.

Figure 12 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 Nm, 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.

The control single of the load torque controller in the simulation is expressed by Figure 13. 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 14 depicts the simulative results of the voltage-current performance controlled by LTCS. According to the discussion of Figure 4, the current, , is defined as , which follows the relationship of (4). And the rotate speed of EM, , is adjusted by changing the voltage.

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.

Figure 15 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.

By comparing the heavy load control performance of LTCS and LTFC in Figure 15, 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.

Figure 16 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 17 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/kNh, which can enhance the seeding uniformity.

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.

Figure 18 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.

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 Figure 19. 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 20 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.

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.

#### Nomenclature

: | Maximum opening angle of AP (°) |

: | Distance from barycenter to front axle (m) |

: | Depth of AP (°) |

: | Width of single ploughshare (cm) |

: | Constant loss of EM (kW) |

, : | Acceleration factors (-) |

: | Per phase back EMF (V) |

: | Back electromotive force (V) |

: | Air resistance (kN) |

: | Rolling resistance (kN) |

: | Plowing resistance (kN) |

: | Acceleration resistance (kN) |

: | Slope resistance (kN) |

: | Tractive force (kN) |

: | Driving force (kN) |

: | Rolling resistance coefficient (-) |

: | Tractor gravity (kN) |

: | Plowing depth (cm) |

: | Hitch point height (m) |

: | Per phase current (A) |

: | Winding current (A) |

: | Gear ratio of AI (-) |

: | Gear ratio of GT (-) |

: | Gear ratio between EM and PTO (-) |

: | Gear ratio of MD (-) |

: | Rotational inertia of drivetrain (kgm^{2}) |

: | Soil specific resistance (kN/cm^{2}) |

: | Back EMF constant (Vmin/r) |

: | Torque constant (Nm/A) |

: | Copper loss coefficient (kW/N^{2}m^{2}) |

: | Iron loss coefficient (kWmin/r) |

: | Air resistance coefficient (kWmin^{3}/r^{3}) |

: | Wheelbase (m) |

: | Per phase self-inductance (H) |

: | Mutual inductance (H) |

: | Self-inductance (H) |

: | Leakage inductance (H) |

: | Free travel ratio of AP (-) |

: | Copper loss (kW) |

: | Icon loss (kW) |

: | Electric power input of EM (kW) |

: | Mechanical loss (kW) |

: | Mechanical power output of EM (kW) |

: | Rating power of EM (kW) |

: | Searching range of rotate speed (r/min) |

: | Stator winding resistance (Ω) |

: | Searching range of torque (Nm) |

: | Rolling radius of drive wheel (m) |

: | Electromagnetic torque (Nm) |

: | Load torque of EM (Nm) |

: | Maximum torque of EM (Nm) |

: | Torque output of driving wheel (Nm) |

: | Torque output of AI (Nm) |

: | AP release time (s) |

: | Per phase voltage (V) |

: | Input voltage (V) |

: | Viscous resistance coefficient (Nmmin/r) |

: | Tractor speed (km/h) |

: | Number of ploughshares (-) |

: | Rotate speed of EM (r/min) |

: | Rotate speed of driving wheel (r/min) |

: | Rotate speed of AI (r/min) |

: | Mechanical efficiency of AI (%) |

: | Transmission efficiency (%) |

: | Rolling efficiency (%) |

: | Mechanical efficiency of GT (%) |

: | Mechanical efficiency between EM and PTO (%) |

: | Electric conversion energy efficiency of EM (%) |

: | Tractive efficiency (%) |

: | Slip efficiency (%) |

: | Mechanical efficiency of MD (%) |

: | Inertia weight (-) |

: | Quantization factor of footplate depth (-) |

: | Footplate deepening rate quantization factor (-) |

: | Number of iterations (-) |

: | Slip rate (-). |

#### Abbreviations

AI: | Agricultural implement |

AP: | Accelerator pedal |

BLDC: | Brushless direct current motor |

EM: | Electric motor |

EMF: | Electromotive force |

ES: | Energy system |

GT: | Gearbox transmission |

ICE: | Internal combustion engine |

LTCS: | Load torque control strategy |

LTFC: | Load torque followed control method |

MC: | Motor controller |

MD: | Main drive |

PSO: | Particle swarm optimization algorithm |

PTO: | Power take-off shaft. |

#### Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank the YTO Group Corporation for the study permission and technical support. The project is supported by National Natural Science Foundation of China (no. 51375145), Science and Technique Foundation of Henan Province (Grant no. 142102210424), and Research Program of Application Foundation and Advanced Technology of Henan Province (no. 152300410080). The authors thank Accdon for its linguistic assistance during the preparation of this paper.

#### Supplementary Materials

Description of Supplementary Material A:

Corresponding to Figure 9 in the text, data of plowing parameter numbering, time, equivalent resistance and soil specific resistance in the tractive simulation can be seen in the table of Supplementary Material A. In this table, different colors are used to discriminate each set of plowing parameters in Table 2. And the main independent variable, soil specific resistance, is changed from 0.007 to 0.01, which contains a wide soil humidity scope of common soil conditions in the Northeast China region.

Description of Supplementary Material B

The figure in Supplementary Material B describes the initialization signals’ superposition process during the simulation. As stated in the text, there are two initialization signal source, time trigger and load torque controller, in load torque control strategy (LTCS). Time trigger outputs the initialization signal periodically, and the load controller produces the initialization signal following the fuzzy rules of Figure 8. The initialization signal of simulation is obtained by the superposition of the two initialization signals.