Next Article in Journal
Reduced-Order Electro-Thermal Battery Model Ready for Software-in-the-Loop and Hardware-in-the-Loop BMS Evaluation for an Electric Vehicle
Previous Article in Journal
Fast Charging Systems for Passenger Electric Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Assessment of Low-Cost Tractor Motorization with Main Farming Implements

1
EAS Research Team, Laboratory of Research in Engineering, National High School for Electricity and mechanics (ENSEM), Casablanca 8118, Morocco
2
Foundation of Research Development and Innovation in Sciences and Engineering (FRDISI), Casablanca 8118, Morocco
*
Author to whom correspondence should be addressed.
World Electr. Veh. J. 2020, 11(4), 74; https://0-doi-org.brum.beds.ac.uk/10.3390/wevj11040074
Submission received: 28 August 2020 / Revised: 20 October 2020 / Accepted: 20 October 2020 / Published: 18 November 2020

Abstract

:
The objective of the present work is to evaluate the performance of a low-cost tractor equipped with a parallel hybrid engine, which was simulated using AMESim software. The tractor was evaluated with three different farming implements attached to the tractor, and each implement requires a different type of power. The first simulation was executed without any implements attached. The tractor was able to run for 170 s with the electric motor only, which resulted in fuel savings during this period. The first implement, a moldboard plow, was attached for the second round of evaluation, and the electric motor ran by itself for 150 s, which also led to fuel savings during operation. During the third simulation, the tractor was attached to a Bette Harvest, which has a very high-power demand. The obtained results show that both engines were engaged to provide the required energy. During the final round of evaluation, simulations were run for a straw tub grinder. In this simulation, the electric motor ran alone until the battery was fully discharged. Thereafter, the combustion engine was activated in order to facilitate operations and to charge the battery. The results show that the parallel hybrid architecture employed for the low-cost tractor significantly decreased the CO2 emissions and minimized the consumption of fuel.

1. Introduction

The continuing increase in the global population requires an adequate increase in crop production to accommodate it [1]. Moreover, as the population grows, global meat consumption and the demand for bioenergy and raw materials are projected to increase [2,3]. Farm mechanization is particularly essential towards providing increased crop production, now and probably even more so in the future [4]. In addition, agricultural mechanization is generally a necessary first step in the general mechanization of a nation so that it can shift from poverty to security and prosperity. Tractors represent the largest segment of the agricultural machinery industry and play a central role in the mechanization of agriculture [5]. Global annual tractor sales are estimated at approximately 1.8 million units, with a sales volume of USD 45 billion in 2017 [6].
The Traction farming tasks is predominantly carried out by tractors with the several power [7]. The efficiency of a tractor is measured by the amount of work performed versus the cost of this work. Hitching work is affected by the traction and travel speed of the tractor used [5,8].
Because of the energy density of diesel fuel, it is commonly assumed that diesel engines will be predominantly used in the agricultural sector in the coming decades [9,10]. However, the U.S. Environmental Protection Agency (EPA) estimates that if off-road vehicle emissions are not controlled, then these vehicles will contribute 33% of hydrocarbon (HC) emissions, 9% of carbon monoxide (CO), 9% of nitrogen oxides (NOx), and 2% of particulate matter (PM) emissions in the United States over the coming years [11].
The automotive industry has already started to move towards the construction of green, environmentally friendly vehicles. This approach focuses mainly on the three areas: emissions, noise, and fuel consumption. This strategy aims to significantly reduce vehicle emissions and increase the fuel economy of vehicles, both of which are beneficial for the environment [12].
From this point of view, tractor hybridization is a valid mechanism to ensure that the environment has a cleaner future and that the costs of the agricultural industry are minimized, as scientifically demonstrated in a previous work [13].
In the present work, the behavior of a hybrid powertrain of a low-cost tractor was evaluated using the AMESim simulation tool with a moldboard plow, Bette Harvest, or straw tub grinder attached to it.

2. State of the Art

Hybrid technology is an advanced technology for improving fuel efficiency in the automotive industry. It has improved fuel efficiency in automobiles by 25% and has also contributed to improving the fuel efficiency of construction machinery [14].
In order to improve construction machinery, Komatsu developed a model (HB205) using hybrid technology based on excavator rotation that requires high power, which has reduced fuel consumption by 25% compared with its non-hybrid model [15]. Caterpillar also developed a hybrid bulldozer (D7E Hybrid) and increased fuel efficiency by 20 percent compared with its predecessor [16]. In contrast to the automotive industry, the construction machinery industry has applied hybrid technologies to workpieces such as boom and bucket systems, and most studies have focused on using excavators [17].
In particular, various studies on hydraulic pressure control strategies for motor operating parts have been carried out to improve work efficiency. For example, Wang and Wang (2014) developed a pressure compensation system to improve the energy efficiency of a hybrid hydraulic system and performed an evaluation using a bench test. They reported that 26–33% of the energy that would have otherwise been used by this system was recovered [18]. Shen, Jiang, Su, and Karimi (2015) proposed optimal control of the variable trajectory of hybrid excavators to reduce the fuel consumption of off-road vehicles, and Lin, Wang, Hu, and Gong (2010) improved the efficiency of hybrid excavators by 17% through the development and simulation of energy regeneration systems [19,20].
Choi, Kim, Yu, and Yi (2011) developed an optimal control system for hybrid excavators.
Additionally, Yoon, Truong, and Ahn (2013) developed a parallel hybrid excavator by applying an electro-hydraulic actuator, which reduced the energy usage of the machine by 60% [21,22].
In other studies, hybrid technologies have been applied to the driving components of construction machines. Zeng, Yang, Peng, Zhang, and Wang (2014) applied several energy management strategies for wheel loaders and obtained significant increases in fuel efficiency through optimal control of the power drive [23]. Hui and Jungqing (2010) developed a parallel hybrid hydraulic system to reduce energy consumption during frequent wheel loader stop/start phases and demonstrated improved work efficiency and fuel savings through simulations [24]. Dagci, Peng, and Grizzle (2015) applied a hybrid power distribution system with two single planetary gears (SPGs) in light-duty trucks, and Keulen, Mullem, Jager, Kessels, and Steinbuch (2012) proposed an optimal control strategy that was adaptable to different truck masses and road elevations and applied hybrid technology to heavy-duty trucks [25,26].
In the agricultural machinery industry, John Deere has developed a mild hybrid tractor (Model 7030E) that operates the engine compressor cooling and air conditioning systems with one engine, which has significantly increased fuel efficiency. The efficient use of power in the drive section is important because tractors require high traction force when operating with towing implements and performing soil loading; however, there have been no studies on the use of hybrid technologies in the drive section. Thus, this study was conducted to develop a hybrid tractor. The objectives of this study were (1) to build a parallel hybrid tractor using components of a hybrid drive system, (2) to establish power management strategies for the hybrid tractor, and (3) to evaluate the performance of the hybrid tractor by comparison with a conventional tractor through field testing [27].

3. Hybrid Architecture Adopted

The parallel hybrid architecture used in this study, which is shown in Figure 1, consists of a direct mechanical connection between the hybrid power unit and the wheels. It is equipped with an electric traction motor that drives the wheels and can recover some of the braking energy to charge the batteries (regenerative braking) or to assist the internal combustion engine during acceleration.
Furthermore, the internal combustion engine and the electric motor are coupled by a mechanical device.
The electric motor can then be designed with reduced capacity, i.e., lower cost and volume.
There are several configurations of the mechanical device used for the combination, depending on the structure between the internal combustion engine and the electric motor. There can be a coupling with a single- or two-shaft configuration, speed coupling with a planetary gearbox, or a fusion of the two [28,29].
The appropriate operating mode is programmed or manually switched on [13]:
  • At high speeds, the combustion engine is used as the drive.
  • At low speeds, the electric motor is activated to optimize fuel efficiency and consumption; the electric motor does not run at a standstill to save the battery.
  • When the driver strongly accelerates, both the electric and combustion engines work simultaneously to transmit more power.
When the driver lifts his or her foot while driving or braking, the electric motor works as a generator and recovers energy. The combustion engine is then decoupled from the transmission and thus, generates few losses.
The low-cost tractor is an agricultural tractor that weighs 4.5 t and is equipped with a four-stroke V8 combustion engine with a total displacement of 12 L and a power of 120 kW. The appropriate selection of the electric part of the hybrid vehicle is a major element of the successful marriage between the thermal and electric engines. The permanent magnet synchronous motor (MSAP) appears to be an optimal solution for automotive traction because of its technical performance and, in particular, its compactness and efficiency. This type of motor is used by Toyota in the Prius for its high efficiency and good dynamic performance, which is attributed to the low stator inductances due to the large width of the apparent air gap, large magnetic field in the air gap, and lack of a DC voltage source for excitation. The low-cost tractor in this study is equipped with 20 kW synchronous motors. The outputs of the two engines are linked by a planetary gear train. The output power is adapted by a manual gearbox with four gears. In addition, it is equipped with a generator that can recover the electric energy from the combustion engine and store it in a lithium battery.

4. Mathematical Model

In the present work, we examined a model of a parallel hybrid tractor using AMESim software. The main notation used in this paper is shown in the Table 1. We produced the 1D sketch shown in Figure 2 and then determined the appropriate mission for each attached implement in order to evaluate the tractor’s behavior during the mission.
The duty cycles of agricultural tractors are very different from the cycles known to the automotive industry. The hybrid tractor operates according to the maximum speeds for each implement, as shown in Table 2 [30].

4.1. Definition of the Sub-Models Composing the Sketch

4.1.1. Tractor

The tractor sub-model is a 1D vehicle load sub-model and is used to calculate longitudinal acceleration. The speed and displacement are defined for a vehicle with two axles, taking into consideration the rolling friction and road gradient [31].
The acceleration is calculated based on the following forces:
Engine force:
F e n g i n e = T 1 r a d 1 + T 2 r a d 2
Slope strength [32]:
F s l o p e = M c a r . g . s i n ( α s l o p e )
Rolling frictional force [33]:
F f r i c t i o n = r v e h i . V c a r + F r o l l . s i g n ( V c a r )
F r o l l = R r o l l ( M b e c a u s e . g . c o s ( α s l o p e ) W z . V c )
V c   =   ( V c a r   +   V w i n d ) .   | V c a r   +   V w i n d |
Aerodynamic force [33]:
F a e r o = W x . V c
according to Newton’s second law.
The sum of all forces applied to the vehicle is divided by the body mass, including wheel inertias and trailer mass, to obtain vehicle acceleration.
A c a r = 1 M v e h i + M t r a i l + r i n 1 ( r a d 1 ) 2 + r i n 2 ( r a d 2 ) 2 . ( F e n g i n e F a e r o F s l o p e F f r i c t i o n F t r a i l e r )
Then, this acceleration is integrated twice to obtain the speed and position of the vehicle.

4.1.2. Control Unit

The control unit receives information from the driver (acceleration, braking controls, and gearbox ratio), the electric motor (maximum/minimum speed and torque), the combustion engine (maximum speed and torque), the battery (state of charge), and the vehicle speed and the gearbox (input shaft speed). The unit analyzes these values in order to minimize battery consumption. The electric motor can be used as a generator to charge the battery when the driver brakes.
In this approach, the electric motor and the combustion engine are connected to the manual gearbox. The control unit manages the power required by the engine and the electric motor. If the battery needs to be regenerated, then the motor is used to move the vehicle forward, and if the power demand is less than its optimum power (depending on its speed), then the difference is transmitted to the electric motor/generator to charge the battery.
The strategy consists of four modes:
  • Mode 1: The battery must not be charged; the combustion engine is not used.
  • Mode 2: The battery does not need to be charged; the combustion engine is used.
  • Mode 3: The battery is fully charged; the engine is not used.
  • Mode 4: The battery is not fully charged; the engine is used.
There are two transient modes when the engine starts:
  • If the engine speed is less than idle speed/2, then the clutch is engaged and the engine does not start.
  • When the engine reaches idle speed/2, the clutch remains closed and the engine starts [34].
The braking torque required by the driver is calculated following the selected brake command, as defined in Table 3.
  • Mode 1: The battery is fully charged; the combustion engine is not used.
In this mode, the electric motor is used to drive the vehicle, the combustion engine is stopped, and the clutch is opened.
The mode changes if one or both of the following conditions apply:
  • The battery charge state is below the lower limit (parameter) and the engine can start (the rotation speed of the gearbox input shaft is higher than the engine idling speed).
  • The vehicle speed is above the speed limit (parameter).
The torque required by the driver is calculated as defined in Table 4 and Table 5.
  • Mode 2: The battery is not fully charged; the combustion engine is used.
In this mode, the engine is used to drive the vehicle (only during acceleration), and the clutch is closed.
The mode changes if one or both of the following conditions apply:
  • The battery charge status is below the lower limit (parameter).
  • Vehicle speed is below the speed limit (parameter).
The torque required by the driver is calculated as defined in Table 6.
The motor control is calculated according to Table 7.
The torque control of the electric motor and the braking control of the vehicle is calculated as shown in Table 8.
  • Mode 3: The battery is fully charged; the engine is not used.
In this mode, the electric motor is used to drive the vehicle, the combustion engine is stopped, and the clutch is opened. However, we account for the fact that the engine must be started to regenerate the battery.
The mode will change if one or more of the following conditions apply:
  • The battery’s state of charge is above the upper limit (parameter).
  • Vehicle speed is above the speed limit (parameter).
  • The engine can be started (the rotation speed of the gearbox input shaft is higher than the engine idle speed).
The torque required by the driver is calculated as shown in Table 9.
The torque control of the electric motor and the braking control of the vehicle are calculated according to Table 10.
  • Mode 4: The battery is not fully charged; the engine is used.
In this mode, the electric motor is used to charge the battery, and the engine is used to drive the vehicle and charge the battery (only during acceleration).
The mode changes if one or both of the following conditions apply:
  • The battery’s state of charge is above the upper limit (parameter).
  • The engine cannot be used (the speed of rotation of the gearbox input shaft is lower than the engine idle speed).
The torque required by the driver is calculated as shown in Table 11.
The torque control of the electric motor and the braking control of the vehicle is calculated as according to Table 12.
Engine control is calculated as defined in Table 13 [35].

4.1.3. Gearbox

This model can be used to perform simple dynamic modeling of an n-speed manual gearbox (forward and reverse). The transmission losses are the same for each gear. The inertia of the input shaft and the drive axle are taken into account.
This gearbox is composed of an input shaft, flywheel, gears, clutch/synchronizer, and drive axle [36].
The gearbox is illustrated in Figure 3.
  • Torque loss
The torque loss set by the user must always be positive. A pre-treatment is carried out to ensure that the torque is always resistant to movement.
No further pre- or post-processing is carried out by the component.
The value read from the file is applied to the torque of the secondary shaft [37].
T p r a t i o = T s + T l o s s
  • Input shaft speed
The input shaft speed ωp is calculated from the motor torque at port 3, T3, and the input shaft torque, Tp (see below for the calculation of Tp).
The output rotation speed ω 3 is determined as follows [38]:
ω 3 = ω p
  • Rotational speed of the secondary shaft
The rotational speed of the secondary shaft ω s is calculated as follows [39]:
ω s = ω p a p a r a t i o
  • Torque transmitted by the secondary shaft
The relative speed ( ω r e l ) between the primary and secondary shafts must first be calculated [40]:
ω r e l = ω p r a t i o + ω s
The torque on the secondary shaft T s is the torque transmitted by the synchronizer; its calculation depends on the chosen friction model.
Hyperbolic tangent:
T s = T s l i p t a n h ( 2 ω r e l d v e l )
  • Torque transmitted by the drive shaft
The torque calculation T p of the input shaft depends on the definition of the selected speed loss [41].
T p = 1 n u T s r a t i o   T p = T s + T l o s s r a t i o
  • Torque transmitted by the drive axle
We assume the efficiency of the drive axle to be 1.0 (no loss).
The torque transmitted by the drive axle T p a   ( = T 2 ) is calculated as follows [42]:
T p a = T s p a r a t i o
  • Power Losses
There are three sources of power loss: the viscous friction of the flywheel, clutch slippage, and the efficiency of the machines.
In the event of a loss of power, the power lost from the gearbox is calculated as follows [43]:
p o w e r l o s t = v i s ω p 2 + | p o w e r P r p o w e r P a |

4.1.4. Internal Combustion Engine

The internal combustion engine (ICE) is a sub-model with different application cases.
To apply corrections and dependencies for engine starting, cold and hot temperatures are provided and can be used.
This engine component must be used whenever information on engine performance during a cycle is required. The prediction of fuel consumption or emissions during driving cycles is accurate and has a low CPU cost (compatible with a fixed time step).
Mean effective pressures (BMEP, FMEP) and couples are related by the following relationships for a four-stroke engine [13]:
B M E P = 4 π T o u t 0.01 V
F M E P = 4 π T f r i c 0.01 V
  • Torque calculation
The engine output torque is calculated from the engine control unit (ECU) load demand as follows [44]:
T o u t = T m a x c l o a d + T m i n c ( 1 l o a d )
Several corrections can be activated/applied to the torque values read from the file.
These corrections are maintained simultaneously when they are all activated.
  • Dynamic Correction
T m a x c = T m a x f i l e ρ a i r ρ r e f T a m b 273.15 + T r e f
  • Impact of temperature
If the “use temperature influence correction for FMEP” (with torque impact) option is enabled, any variation in FMEP due to a temperature that differs from the hot reference temperature will impact the engine’s torque performance.
The correction is as follows [45]:
T m a x c = T m a x f i l e Δ T f r i c   a v e c   Δ T f r i c = { T f r i c , T e m p T f r i c , T h o t   If   the   correction   is   activated 0 otherwise
  • Fuel consumption
The engine’s fuel consumption is read and corrected if necessary.
These corrections are maintained simultaneously when all of them are activated, and they can be applied to the value read from the:
Impact of temperature;
Overconsumption when starting the engine;
Deactivation of cylinders.
The actual fuel consumption is calculated as follows [46]:
c o n s = c o n s f i l e ( 1 + k s t a r t ) e x p f c o n s
The coefficient expfcons is set to 1 when the temperature influence option is not activated.
  • Engine Emissions
  • Equivalence report
The equivalence ratio is used to determine the total exhaust mass flow rate and is calculated as follows [47]:
ϕ = ϕ f i l e ( 1 + k s t a r t ) e x p e q r a t i o
Note: The expeqratio coefficient is set to 1 when the temperature influence option is not activated.
  • Exhaust gas mass flow rate
The exhaust flow rate is calculated from the fuel consumption and the equivalence ratio, as defined in Table 14.
  • Pollutant emissions
Pollutant emissions are read from related data files and corrected if necessary [48].
e m i C O = e m i C O f i l e ( α t + β t c o e f c o l d C O ) If   temperature influence   activated ( 1 + k s t a r t c o e f s t C O ) If   start   influence activated
α t = T e n g T l o w T h i g h T l o w   a d   β t = 1 α t
When the engine temperature is below the “high engine temperature threshold”, a linear interpolation between the low and high temperature thresholds is performed to correct the emissions read from the user-defined data tables.
  • Exhaust gas temperature
T e x h a u s t = T e x h a u s t f i l e e x p e x h t e m p
The coefficient expexhtemp   is set to 1 when the temperature influence option is not activated [49].
  • Heat loss from engines
Friction mean effective pressure:
F M E P = F M E P f i l e e x p f m e p c o r
The coefficient expfmepcor is set to 1 when the temperature influence option is not activated or if the FMEP file includes an engine temperature axis.
The frictional power losses [W] are obtained from the friction torque [W]. T f r i c [Nm] and engine speed w [rad/s] are determined as follows [50]:
P f r i c = T f r i c ω
  • Losses of combustion heat
The combustion heat losses (Pcomb power [W]) are calculated using the heat loss coefficient read from the table heatwall_i.data and the quantity of fuel injected [51].
P c o m b = c o n s p c i c o e f h e a t w a l l e x p h l o s s e s
  • Energy balance
In custom mode, it is possible to activate the energy balance function. The energy balance can be used to determine exhaust gas temperature or combustion heat losses.
The energy balance is established as follows [52]:
ϕ c o m b = P w o r k + P f r i c + P p u m p + P e v a p + P c o m b + P e x h
The heat of combustion φ c o m b [W] indicates released energy. This energy takes into account the production of CO and HC (instead of only CO2 and H2O for complete combustion) [53].
ϕ c o m b = ( c o n s e m i H C e m i C O 4 ) p c i
The fuel evaporation energy P e v a p [W] is deduced from the fuel consumption and the latent heat of evaporation H e v a p [J/kg] and is equal to [54]:
P e v a p = c o n s H e v a p
The motor output power and friction losses Pwork, Pfric [W] are equal to [55]:
P w o r k = T o u t w   a n d   P f r i c = T f r i c w
The motor pumping losses (when available) P p u m p [W] are deduced from the value read in the pump PMEP file [bar] as follows [56]:
P p u m p = T p u m p w
The exhaust power of the engine P e x h   [W] is equal to [57]:
P e x h = ( d m B G C p B G + d m u n b u r n e d a i r C p a i r + d m u n b u r n e d f u e l C p c a r b ) ( T e x h a u s t T a m b )
  • Deactivation of cylinders
The engine component receives the number of deactivated cylinder(s) at port 4 and the activation (or not) of specific corrections.
When the input signal to port 4 applies the data, correction is set to 0.0 (option 1), and no correction is applied to the supplied data file.
When the input signal to port 4 applies the data, correction is set to 1.0 (option 2), and corrections are applied to the value read from the data file.
When option 2 is selected, the corrections are as follows [58]:
F M E P c = F M E P n c y l a c n c y l + F M E P 0 n c y l d e a c n c y l
B M E P m a x c = B M E P m a x n c y l a c n c y l F M E P 0 n c y l d e a c n c y l
Note: The minimum BMEP is not affected by cylinder deactivation.
An equivalent BMEP (BMEPeq) is calculated and used to read the data file. This is the BMEP that the engine will produce with all cylinders activated and with the same ECU load demand.
It is calculated as follows [58]:
B M E P e q = l o a d 1 n c y l a c ( n c y l B M E P m a x + n c y l d e a c F M E P 0 ) + ( 1 l o a d ) B M E P m i n
The engine’s fuel consumption and emission variables are also corrected if the cylinder is deactivated with option 2 selected [13].
c o n s = c o n s s t d n c y l a c n c y l  
e m i C O = e m i C O s t d n c y l a c n c y l  
The equivalence ratio and exhaust gas temperature variables are not affected by cylinder deactivation.
  • Displacement variation
As with cylinder deactivation, the data file for combustion mode, n, with displacement variation must be completed using the BMEP calculated with the standard displacement. In this case, the standard displacement (engine parameter) is the largest available engine displacement.
The variation in the displacement vdis (≤1) is an input signal. The correction is made as follows:
B M E P e q = B M E P e q v d i s
The BMEPeq is used to read the data files [59].

4.1.5. Electric Motor/Generator

The following is a sub-model of an electric motor/generator with its converter. It is bi-directional (motor/generator) and independent of the technology of the motor and its converter.
  • Torque
From the input torque demand T s e t , the couple T l i m is limited as follows:
T m i n T l i m T m a x
The output torque T m   is determined from the limited torque T l i m using a first-order offset [60]:
T m = 11 + t r s T l i m
where tr is the user-defined time constant [s].
  • Power balance
The mechanical power P m e c [ W ] is calculated as follows [61]:
P m e c = T m ω
The power loss P l o s t   [ W ] is calculated from the efficiency   η :
P l o s t = ( 1 η ) | P m e c |
Two modes can then be defined:
  • P e l e c   >   P m e c Motor mode.
  • P e l e c   <   P m e c Generator mode.
The different modes of operation of the electric motor/generator are illustrated in Figure 4.
Therefore, the relationship between mechanical power and electrical power, taking the below conditions into account, is as follows [62].
P e l e c = P m e c P l o s t
with
  • P m e c < 0 in engine mode.
  • P m e c > 0 in generator mode.
  • P l o s t > 0 .
The efficiency of the motor/generator, either user-defined or derived from P l o s t , corresponds to the following.
Motor mode [63]:
η = 2 P e l e c P m e c
Generator mode:
η = P e l e c P m e c
  • Limitation
It is the responsibility of the user to ensure that the power loss is less than the input power, i.e.,
  • P l o s t < P e l e c in engine mode.
  • P l o s t < P m e c in generator mode.
  • Electric current
The output current I5 [A] is calculated as follows [63]:
I 5 = P m e c P l o s t U
where U is the input voltage [V] and given by [64]
U = V 5 V 6

4.1.6. Battery

The following describes a battery sub-model, which includes an internal resistance model in the case of a variable voltage. Experimental data are required to define the open-circuit voltage and internal resistance. Furthermore, thermal effects can be taken into account. The open-circuit voltage and internal resistance can be temperature-dependent.
This battery is a concatenation of cells in series and in parallel.
Output potential:
If the voltage is fixed, then the potential at the positive pole is calculated as follows [65]:
V 3 = f i x e d v o l t a g e + V 2
where
V 3 = f i x e d c e l l v o l t a g e S c e l l + V 2
If the voltage is variable, then the potential at the positive pole is calculated as follows [66]:
V t a r g e t = O C V b a t R b a t I 3
V t a r g e t = S c e l l ( O C V c e l l R c e l l I c e l l )
I c e l l = I 3 P c e l l
The voltage is calculated using the target V [67]:
d V d t = V t a r g e t V t a u
where
d V d t = I 3 V t a r g e t O C V b a t R b a t C f b a t
where
d V d t = I 3 V t a r g e t O C V b a t R b a t C f c e l l . Pcell Scell
and, finally,
V 3 = V + V 2
  • State of charge (SOC)
The SOC [%] of the battery is a variable state whose derivative is calculated as follows [68]:
d S O C d t = d q d t 100 C n o m
where C n o m   is the nominal capacity [As]. This derivative is limited so that the state of charge remains within a range of [0;100%].
The DOD [%] is used as an input to read the open-circuit voltage and internal resistance data files. It is derived as follows:
D O D = 100 S O C
  • Charge used by the load
The charge used by charge q [As] (or the charge removed from the battery) is calculated as follows:
d q d t = I 3
where I3 is the battery current at port 3 [A].
Note that with the TM software conditions, the battery will discharge when the battery current at port 3, I3, is negative.
  • Output power at port 1
The output power P [W] at port 1 is calculated as follows [67,68]:
P = R b a t I 3 2
where
P = R c e l l I c e l l 2 S c e l l P c e l l

5. Results and Discussion

The simulations were carried out for a low-cost tractor. The driving cycle adopted was 300 s. Additionally, the maximum working speed was specific for each tool attached to the tractor, as shown in Table 2. According to the optimal torque distribution, we found that the hybrid tractor had the following main features. The electric motor starts the tractor up to a certain power, or the combustion engine is switched on to provide traction and simultaneously recharge the battery via a generator. At a certain stabilized and low tractor power, the traction is in pure electric mode. All decelerations of the vehicle are provided by the electric motor, thus recovering the braking energy.

5.1. Tractor without Attached Tools

In this simulation, the tractor did not have any tools attached to it and had a variable speed, as shown in Figure 5 below. The maximum speed was 30 km/h.
Figure 6 shows the torque engine variation during the driving cycle described in Figure 5 when there are no implements attached to the tractor.
Figure 7 shows the variation in the torque of the electric motor during the driving cycle described in Figure 5 when there are no implements attached to the tractor.
The states of the vehicle’s engines confirm the results of the simulation. As illustrated in Figure 7, the electric motor continuously contributes throughout the driving cycle. Additionally, Figure 6 reveals a load for the internal combustion engine only after 180 s at high power demands. The generator is activated when the internal combustion engine is running to store energy in the battery. The role of the battery is to supply power to the electric motor. It should be noted that the addition of the electric motor significantly reduces the use of the internal combustion engine, and it is even possible to realize an all-electric operation.
Figure 8 shows the variation in the battery’s state of charge during the driving cycle shown in Figure 5 when there are no implements attached to the tractor.
Figure 8 shows the battery’s state of charge (%) during the mission profile. The SOC begins at 100%, and the operating range is between 100% and 75%.
Figure 9 shows the total fuel consumption during the driving cycle shown in Figure 5 when there are no implements attached to the tractor.
Thus, fuel is not consumed until after 170 s of operation, as shown in Figure 9.
The downward trend in the curve reflects the nature of the discharge during the simulation period. The fluctuating SOC is caused by the battery being powered by the recuperative generator. Because the battery’s operating limit is at a low SOC level, the tractor reaches a point at which the engine must start. Therefore, the all-electric mode time in this driving cycle is 170 s (Thermal Engine Status).
Figure 10 illustrates the variation in the CO2 released during the adopted mission profile. The curve shows an emission of 0 g for 170 s of operation; this is the result of the 100% electric operation, which minimizes the emissions rate.

5.2. Tractor with a Moldboard Plow

In this simulation, the moldboard plow tool was attached to the tractor. The speed was variable and had a maximum of 7 km/h, as shown in Figure 11 below.
Figure 12 shows the variation in the torque of the internal combustion engine during the driving cycle described in Figure 11 when a moldboard plow is attached to the tractor.
Figure 13 shows the variation in the torque of the electric motor during the driving cycle in Figure 11 when a moldboard plow is attached to the tractor.
The simulation results in Figure 13 show that the electric motor continuously contributes throughout the driving cycle, and it operates alone for 150 s. Figure 12 reveals a load for the internal combustion engine only at high power demands. With the generator functioning during operation, the internal combustion engine ensures that energy is stored in the battery, which is responsible for supplying power to the electric motor. The addition of the electric motor significantly reduces the use of the internal combustion engine.
Figure 14 shows the variation in the battery’s state of charge during the driving cycle in Figure 11 when a moldboard plow is attached to the tractor.
Figure 14 shows the variation in total fuel consumption during the driving cycle in Figure 11 when a moldboard plow is attached to the tractor.
Figure 14 shows the battery’s state of charge (%) during the mission profile. The SOC starts at 100%. The operating range is between 100% and 10% and then increases because of the contribution of the combustion engine. Thus, fuel is consumed only after 170 s of operation, as shown in Figure 15. The downward trend in the curve reflects the nature of the discharge during the simulation period. The fluctuating SOC is caused by the battery being powered by the recuperative generator. Because the battery’s operating limit is at a low SOC level, the tractor reaches a point at which the engine must start. Therefore, the all-electric mode lasts for 150 s in this driving cycle (Thermal Engine Status).
Figure 16 illustrates the variation in the CO2 released during the mission profile illustrated in Figure 11. The curve shows an emission of 0 g for 170 s of operation; this is the result of the 100% electric operation, which minimizes the emissions rate.

5.3. Tractor with the Bette Harvest

In this simulation, the Bette Harvest variable-speed implement was attached to the tractor, as shown in Figure 17. The maximum speed was 9 km/h.
Figure 18 shows the variation in the torque of the internal combustion engine during the driving cycle described in Figure 17 when a Bette Harvest is attached to the tractor.
Figure 19 shows the variation in the torque of the electric motor during the driving cycle described in Figure 17 when a Bette Harvest is attached to the tractor.
The simulation results in Figure 19 show that the electric motor contributes little power during the driving cycle due to the high power demand of the tool. Figure 18 reveals that the internal combustion engine contributes in this cycle. With the generator functioning during the operation of the combustion engine, energy is stored in the battery, which supplies power to the electric motor. These results show that the addition of the electric motor significantly reduces the use of the combustion engine.
Figure 20 shows the variation in the battery’s state of charge during the driving cycle described in Figure 17 when a Bette Harvest is attached to the tractor.
Figure 21 shows the variation in total fuel consumption during the driving cycle described in Figure 17 when a Bette harvest is attached to the tractor.
Figure 20 shows the battery’s state of charge (%) during the mission profile. The SOC starts at 100%, and the operating range is higher than the highest engine contribution. The fuel consumption of the internal combustion engine is very high because of the high-power requirement, as shown in Figure 21.
Figure 22 illustrates the variation in CO2 emissions during the adopted mission profile. The curve illustrates continuous and high CO2 emissions, which is explained by the considerable intervention of the combustion engine.

5.4. Tractor with a Straw Tub grinder

In this simulation, the straw tub grinder tool was attached to the tractor. The tool does not require any movement, as shown in Figure 23.
Figure 24 shows the variation in the torque of the internal combustion engine during the driving cycle described in Figure 23 when a straw tub grinder is attached to the tractor.
Figure 25 shows the variation in the torque of the electric motor during the driving cycle described in Figure 23 when a straw tub grinder is attached to the tractor.
The simulation results in Figure 25 show that the electric motor contributes substantially during the driving cycle because of the low power demand of the tool. However, the combustion engine contributes after 570 s, which is after the battery is discharged, as shown in Figure 24.
Figure 26 shows the battery’s state of charge (%) during the mission profile. The SOC starts at 100%. The operating range decreases because of the large contribution of the electric motor. Fuel is not consumed until 570 s, when the battery is discharged and the internal combustion engine starts to operate, as shown in Figure 27.
Figure 28 illustrates the variation in CO2 emissions during the adopted mission profile. The curve shows that CO2 is not emitted until 570 s, which is when the combustion engine starts to operate and the battery is discharged.

6. Conclusions

Agriculture is one of the most important industries in the world, and the tractor is regarded as the main axial machine used in this sector.
The present work reports an evaluation of a low-cost tractor that is equipped with a parallel hybrid engine. Simulations were carried out with AMESim software with three different implements attached to the tractor. Each implement requires a different type of power.
First, without any attached implements, the tractor was able to run for 170 s with the electric motor only. This saved fuel during this period.
The second simulation was carried out with a moldboard plow attached to the tractor, which requires draft power. The simulation showed that the tractor worked with the electric motor for 150 s, which reflects notable fuel savings.
In the third simulation, the tractor was attached to a Bette Harvest, which requires draft power and PTO at the same time. Thus, very high power was required, and it could only be provided by placing a demand on both engines.
In the final simulation, the straw tub grinder was used, which only needs PTO. During the operation, the electric motor ran alone until the battery was discharged. Thereafter, the combustion engine started to operate and charged the battery.
Finally, the parallel hybrid architecture can minimize fuel consumption through the participation of the electric motor. H. Fathollahzadeh et al. found that when working the soil with a conventional tractor attached to a moldboard plow, the average fuel consumption was 30 L/ha, whereas the average fuel consumption by our low-cost tractor was 17 L/ha while performing the same work. This represents a significant decrease in consumption and emissions [69]. This improvement is the result of the power obtained by recovering energy from the combustion engine. The possibility of integrating renewable energy into our architecture will be the subject of our future work.

Author Contributions

Conceptualization, Y.Z., M.E.M., S.B. and H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paillard, S. Agrimonde—Scenarios and Challenges for Feeding the World in 2050; Springer Science and Business Media: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  2. Alexandratos, N.; Bruinsma, J. World Agriculture towards 2030/2050; Food and Agriculture Organization of the United Nations: Roma, Italy, 2012. [Google Scholar]
  3. Food and Agriculture Organization of the United Nations. The Future of Food and Agriculture—Trends and Challenges; FAO: Rome, Italy, 2017. [Google Scholar]
  4. Kormawa, P. La Mécanisation Agricole Durable Cadre Stratéguique pour L’Afrique; FAO: Rome, Italy, 2019. [Google Scholar]
  5. Renius, K.T. Fundamentals of Tractors Design; Springer International Publishing: Cham, Switzerland, 2020. [Google Scholar]
  6. Tordo, S.; Lorenzato, G.; Zhao, J.; McEneaney, K.; Sarmiento-Saher, S.P. Options for Increased Private Sector in Resilience Investment Focus on Agriculture; World Bank: Washington, DC, USA, 2017. [Google Scholar]
  7. Salvat, B. Tracteurs Agricoles Renault; France Agricole: Paris, France, 2012. [Google Scholar]
  8. Sakai, J. Tractors: Two-wheel tractors for wet land farming. In Handbook of Agriculture Engineering; International Commission of Agricultural and Biosystems Engineering (CIGR): Liege, Belgium, 1999; Volume III, pp. 54–95. [Google Scholar]
  9. Karner, J.; Baldinger, M.; Schober, P.; Reichl, P.H. Hybrid systems for agricultural engineering. Landtechnik 2013, 68, 22–25. [Google Scholar]
  10. Karner, M.J. Prospects of hybrid systems on agricultural machinery. J. Agric. Eng. 2014, 1, 33–37. [Google Scholar] [CrossRef]
  11. Office of Mobile Sources. Non-Road Engines and Air Pollution; Environmental Protection Agency (EPA): Washington, DC, USA, 1669.
  12. Kannusamy, M.Y.L.; Ravindran, V.; Rao, N. Analysis of Multiple Hybrid Electric Concept in Agricultural Tractor through Simulation Technique; SAE International: Warrendale, PA, USA, 2019. [Google Scholar]
  13. Zahidi, Y.; Benhadou, S.; Medromi, H. Determining the optimal motorization for low-cost tractor. Int. J. Simul. Syst. Sci. Technol. 2020, 21, 1–13. [Google Scholar]
  14. Wang, L.; Zhang, Y.; Yin, C.; Zhang, H.; Wang, C. Hardware-in-the-loop simulation for the design and verification of the control system of a series–parallel hybrid electric city-bus. Simul. Model. Pract. Theory 2012, 25, 148–162. [Google Scholar] [CrossRef]
  15. Inoue, H. Development of hybrid hydraulic excavators. In Proceedings of the 9th JFPS International Symposium, Matsue, Japan, 28–31 October 2014. [Google Scholar]
  16. Zhang, B.; Guo, S.; Zhang, X.; Xue, Q.; Teng, L. Adaptive smoothing power following control strategy based on an optimal efficiency map for a hybrid electric tracked vehicle. Energies 2020, 13, 1893. [Google Scholar] [CrossRef]
  17. Ceresoli, F. Excavator with Two or Three Arms: Dynamical Behavior and Structural Implications; Department of Mechanical and Industrial Engineering, University of Brescia: Brescia, Italy, 2020. [Google Scholar]
  18. Wang, D.; Lu, X.; Wang, H.; Chen, M. A new design of foaming agent mixing device for a pneumatic foaming system used for mine dust suppression. Int. J. Min. Sci. Technol. 2016, 26, 187–192. [Google Scholar] [CrossRef]
  19. Shen, W.; Jiang, J.; Su, X.; Karimi, H.R. Control strategy analysis of the hydraulic hybrid excavator. J. Frankl. Inst. 2014, 352, 541–561. [Google Scholar] [CrossRef]
  20. Lin, T.; Wang, Q.; Hu, B.; Gong, W. Development of hybrid powered hydraulic construction machinery. Autom. Constr. 2020, 19, 11–19. [Google Scholar] [CrossRef]
  21. Choi, J.; Kim, H.; Yu, S.; Yi, K.C. Development of integrated controller for a compound hybrid excavator. J. Mech. Sci. Technol. 2011, 25, 1557. [Google Scholar] [CrossRef]
  22. Yoon, J.I.; Truong, D.Q.; Ahn, K.K. A generation step for an electric excavator with a control strategy and verifications of energy consumption. Int. J. Precis. Eng. Manuf. 2013, 14, 755–766. [Google Scholar] [CrossRef]
  23. Zeng, X.; Yang, N.; Peng, Y.; Zhang, Y.; Wang, J. Research on energy saving control strategy of parallel hybrid loader. Autom. Constr. 2014, 38, 100–108. [Google Scholar] [CrossRef]
  24. Hui, J.J.S. Research on the system configuration and energy control strategy for parallel hydraulic hybrid loader. Autom. Constr. 2010, 19, 100–108. [Google Scholar] [CrossRef]
  25. Dagci, O.H.; Peng, H.; Grizzle, J.W. Power-Split Hybrid Electric Powertrain Design with Two Planetary Gearsets for Light-Duty Truck Applications; International Federation of Automatic Control (IFAC): New York, NY, USA, 2015; pp. 8–15. [Google Scholar]
  26. van Keulen, T.; van Mullem, D.; de Jager, B.; Kessels, J.T.; Steinbuch, M. Design, implementation, and experimental validation of optimal power split control for hybrid electric trucks. Control Eng. Pract. 2012, 20, 547–558. [Google Scholar] [CrossRef]
  27. Bodria, L. Evolution and Prospects of Agricultural Mechanization in the World; Italian Agriculture Machinery Manufacture Federation: Boulogne, Italy, 2016. [Google Scholar]
  28. Ehsani, M.; Gao, Y.; Miller, J.M. Hybrid electric vehicles: Architecture and motor drives. Proc. IEEE 2007, 95, 719–728. [Google Scholar] [CrossRef]
  29. Modern, C.K. Electric Vehicle Technology; Springer: New York, NY, USA, 2001. [Google Scholar]
  30. American Society of Agricultural and Biological Engineers. Agricultural Machinery Management Data; American Society of Agricultural Engineers: St. Joseph, MI, USA, 2000. [Google Scholar]
  31. Baek, S.W.; Lee, S.W. Aerodynamic drag reduction on a realistic vehicle using continuous. Microsyst. Technol. 2019, 26, 11–23. [Google Scholar] [CrossRef]
  32. Khanipour, A.; Ebrahimi, K.M.; Seale, W.J. Conventional design and simulation of an urban hybrid bus. Int. J. Mech. Mechatron. Eng. 2007, 1, 147. [Google Scholar]
  33. Szadkowski, B.; Chrzan, P.J.; Roye, D. A study of energy requirements for electric and hybrid vehicles in cities. In Proceedings of the 2003 International Conference on Clean, Efficient and Safe Urban Transport, Gdansk, Poland, 4–6 June 2003. [Google Scholar]
  34. Bissam, A. Étude et Commande des Véhicules Hybrides Parallèles; Université de Bejaïa: Bejaïa, Algery, 2013. [Google Scholar]
  35. Ehsani, M.; Gao, Y.; Longo, S.; Ebrahimi, K. Modern Electric, Hybride Electric and Fuel Cell Vehicles; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  36. Bouzidi, I. Etude des Vibrations Libres D’une Boite de Vitesse Par la Méthode des Éléments Finis. Ph.D. Thesis, The Universite Abou Bekr Belkaid-Tlemcenfaculte De Technologie, Chetouane, Algery, 2013. [Google Scholar]
  37. Dion, J.L.; Le Moyne, S.; Chevallier, G.; Sebbah, H. Gear impacts and idle gear noise: Experimental study and non linear dynamic model. Mech. Syst. Signal Process. 2009, 23, 2608–2628. [Google Scholar] [CrossRef] [Green Version]
  38. Mèmeteau, H.B. Technologie Fonctionnelle de L’automobile: Transmission, Train Roulant et Equipement Electrique; Dunod: Paris, France, 2014. [Google Scholar]
  39. Martin, M.R.; Hoeijmakers, J. The electrical variable transmission in a city bus. In Proceedings of the 35th Annual IEEE Power Electronics Specialists Conference, Aachen, Germany, 20–25 June 2004. [Google Scholar]
  40. Kowalak, P.; Borkowski, T.; Bonisławski, M.; Hołub, M.; Myśków, J. A statistical approach to zero adjustment in torque measurement of shippropulsion shafts. Measurement 2020, 164, 108088. [Google Scholar] [CrossRef]
  41. Tamada, S.; Bhattacharjee, D.; Dan, P.K. Review on automatic transmission control in electric and non-electric automotive powertrain. Int. J. Veh. Perform. 2020, 6, 98–128. [Google Scholar] [CrossRef]
  42. Adeline, A. Modelisation Dynamique Globale des Boites de Vitesses Automobile; L’Institut National des Sciences Appliquees de Lyon: Lyon, France, 1997. [Google Scholar]
  43. Barathiraja, K.; Devaradjane, G.; Paul, J.; Rakesh, S.; Jamadade, G. Analysis of automotive transmission gearbox synchronizer wear due to torsional vibration and the parameters influencing wear reduction. Eng. Fail. Anal. 2019, 105, 427–443. [Google Scholar]
  44. Hridoy, S.A.A. Lightweight Authenticated Encryption for Vehicle Controller Area Network. Ph.D. Thesis, The Queen’s University, Kingston, ON, Canada, 2020. [Google Scholar]
  45. Kaveh, M.; Abbaspour-Gilandeh, Y. Impacts of hybrid (convective-infrared-rotary drum) drying on the quality attributes of green pea. J. Food Process. Eng. 2020, 43, e13424. [Google Scholar] [CrossRef]
  46. Baytorun, A.N.; Zaimoğlu, Z.; Akyüz, A.; Üstün, S.; Çaylı, A. Comparison of greenhouse fuel consumption calculated using different methods with actual fuel consumption. Turk. J. Agric. Food Sci. Technol. 2018, 7, 850–857. [Google Scholar] [CrossRef]
  47. Hansson, P.A.; Lindgren, M.; Norén, O. A comparison between different methods of calculating average engine emissions for agricultural tractors. J. Agric. Eng. 2001, 1, 80. [Google Scholar]
  48. Coelho, M.C.; Farias, T.L.; Rouphail, N.M. Effect of roundabout operations on pollutant emissions. Transp. Res. D Transp. Environ. 2006, 11, 333–343. [Google Scholar] [CrossRef]
  49. Galindo, J.; Lujan, J.M.; Serrano, J.R.; Dolz, V.; Guilain, S. Description of a heat transfer model suitable to calculate transient processes of turbocharged diesel engines with one-dimensional gas-dynamic codes. Appl. Therm. Eng. 2006, 26, 66–76. [Google Scholar] [CrossRef]
  50. Açıkgöz, B.; Çelik, C.; Soyhan, H.S.; Gökalp, B.; Karabağ, B. Emission characteristics of a hydrogen–CH4 fuelled spark ignition engine. Fuel 2015, 159, 298–307. [Google Scholar] [CrossRef]
  51. Chewning, S.R. Micro Adiabatic Combustion Engine: Concept Development, Simulation and Combustion Experiments. Ph.D. Thesis, Oregon State University, Corvallis, OR, USA, 2009. [Google Scholar]
  52. Abedin, M.J.; Masjuki, H.H.; Kalam, M.A.; Sanjid, A.; Rahman, S.A.; Masum, B.M. Energy balance of internal combustion engines using alternative fuels. Int. J. Simul. Syst. Sci. Technol. 2013, 26, 20–33. [Google Scholar] [CrossRef]
  53. Bertana, M.; Truchot, B.; Marlair, G. Comparison of the fire consequences of an electric vehicle and an internal combustion engine vehicle. In Proceedings of the International Conference on Fires in Vehicles, Chicago, IL, USA, 27–28 September 2012. [Google Scholar]
  54. Janhunen, T.T. Ultra-Low Nox Hcci-Combustion in the Z Engine; SAE Technical Paper; SAE: Warrendale, PA, USA, 2012. [Google Scholar]
  55. Sim, K.; Koo, B.; Kim, C.H.; Kim, T.H. Development and performance measurement of micro-power packusing micro-gas turbine driven automotive alternators. Appl. Energy 2013, 103, 309–319. [Google Scholar] [CrossRef]
  56. Galindo, J.; Serrano, J.R.; Climent, H.; Varnier, O. Impact of two-stage turbocharging architectures on pumping losses of automotive engines based on an analytical model. Energy Convers. Manag. 2010, 51, 1958–1969. [Google Scholar] [CrossRef]
  57. Fu, J.; Liu, J.; Yang, Y.; Ren, C.; Zhu, G. A new approach for exhaust energy recovery of internal combustion engine: Steam turbocharging. Appl. Therm. Eng. 2013, 52, 150–159. [Google Scholar] [CrossRef]
  58. Boretti, A.; Scalco, J. Piston and Valve Deactivation for Improved Part Load Performances of Internal Combustion Engines; SAE International: Warrendale, PA, USA, 2011. [Google Scholar]
  59. Ibrahim, H.; Ilinca, A.; Perron, J. Moteur Diesel Suralimenté Bases et Calculs Cycles Réel, Théorique et Thermodynamique; Université du Québec à Rimouski: Rimouski, QC, Canada, 2006. [Google Scholar]
  60. Hanselman, D. Minimum torque ripple, maximum efficiency excitation of brushless permanent magnet motors. IEEE Trans. Ind. Electron. 1994, 41, 292–300. [Google Scholar] [CrossRef]
  61. Rigacci, M.; Sato, R.; Shirase, K. Experimental evaluation of mechanical and electrical power consumption of feed drive systems driven by a ball-screw. Precis. Eng. 2020, 64, 280–287. [Google Scholar] [CrossRef]
  62. Mohammed, K.G. Mechanical and electrical design calculations of hybrid vehicles. In Applied Electromechanical Devices and Machines for Electric Mobility Solutions; IntechOpen: London, UK, 2020. [Google Scholar]
  63. Chung, C.T.; Wu, C.H.; Hung, Y.H. Evaluation of driving performance and energy efficiency for a novel full hybrid system with dual-motor electric drive and integrated input- and output-split e-CVT. Energy 2020, 191, 116508. [Google Scholar] [CrossRef]
  64. Daguse, B. Modélisation Analytique Pour le Dimensionnement par Optimisation D’une Machine Dédiée à une Chaîne de Traction Hybride à Dominante Électrique. Ph.D. Thesis, Ecole Doctorale Sciences et Technologies de l’Information des Télécommunications et des Systèmes, Paris, France, 2013. [Google Scholar]
  65. Chouchane, M.; Primo, E.N.; Franco, A.A. Mesoscale effects in the extraction of the solid-state lithium diffusion coefficient values of battery active materials: Physical insights from 3D modeling. J. Phys. Chem. Lett. 2020, 11, 2775–2780. [Google Scholar] [CrossRef] [PubMed]
  66. Liu, Y.; Tang, S.; Li, L.; Liu, F.; Jiang, L.; Jia, M. Simulation and parameter identification based on electrochemical- thermal coupling model of power lithium ion-battery. J. Alloy. Compd. 2020, 844, 156003. [Google Scholar] [CrossRef]
  67. Li, A. Analyse Expérimentale et Modélisation d’Eléments de Batterie et de Leurs Assemblages: Application Aux Véhicules Electriques Et Hybrides. Ph.D. Thesis, L’Universite Claude Bernard, Lyon, France, 2006. [Google Scholar]
  68. Urbain, M. Modélisation Electrique et Energétique des Accumulateurs Li-Ion. Estimation en Ligne de la SOC et de la SOHL. Ph.D. Thesis, Institut National Polytechnique de Lorraine, Lorraine, France, 2009. [Google Scholar]
  69. Fathollahzadeh, H.; Mobli, H.; Rajabipour, A.; Minaee, S.; Jafari, A.; Tabatabaie, S.M.H. Average and instantaneous fuel consumption of Iranian conventional tractor with moldboard plow in tillage. ARPN J. Eng. Appl. Sci. 2010, 5, 30–35. [Google Scholar]
Figure 1. Principle of parallel architecture.
Figure 1. Principle of parallel architecture.
Wevj 11 00074 g001
Figure 2. Sketch of the parallel hybrid tractor.
Figure 2. Sketch of the parallel hybrid tractor.
Wevj 11 00074 g002
Figure 3. Operation modes of the gear of a gearbox.
Figure 3. Operation modes of the gear of a gearbox.
Wevj 11 00074 g003
Figure 4. Electric motor operation.
Figure 4. Electric motor operation.
Wevj 11 00074 g004
Figure 5. Tractor driving cycle without any attached implements.
Figure 5. Tractor driving cycle without any attached implements.
Wevj 11 00074 g005
Figure 6. Thermal engine torque.
Figure 6. Thermal engine torque.
Wevj 11 00074 g006
Figure 7. Electric motor torque.
Figure 7. Electric motor torque.
Wevj 11 00074 g007
Figure 8. Battery’s state of charge.
Figure 8. Battery’s state of charge.
Wevj 11 00074 g008
Figure 9. Fuel consumption.
Figure 9. Fuel consumption.
Wevj 11 00074 g009
Figure 10. The variation in the CO2 release.
Figure 10. The variation in the CO2 release.
Wevj 11 00074 g010
Figure 11. Tractor driving cycle with a moldboard plow attached.
Figure 11. Tractor driving cycle with a moldboard plow attached.
Wevj 11 00074 g011
Figure 12. Thermal engine torque.
Figure 12. Thermal engine torque.
Wevj 11 00074 g012
Figure 13. Electric motor torque.
Figure 13. Electric motor torque.
Wevj 11 00074 g013
Figure 14. Battery’s state of charge.
Figure 14. Battery’s state of charge.
Wevj 11 00074 g014
Figure 15. Fuel consumption.
Figure 15. Fuel consumption.
Wevj 11 00074 g015
Figure 16. The variation in the CO2 release.
Figure 16. The variation in the CO2 release.
Wevj 11 00074 g016
Figure 17. Tractor driving cycle with a Bette Harvest attached.
Figure 17. Tractor driving cycle with a Bette Harvest attached.
Wevj 11 00074 g017
Figure 18. Thermal motor torque.
Figure 18. Thermal motor torque.
Wevj 11 00074 g018
Figure 19. Electric motor torque.
Figure 19. Electric motor torque.
Wevj 11 00074 g019
Figure 20. Battery’s state of charge.
Figure 20. Battery’s state of charge.
Wevj 11 00074 g020
Figure 21. Fuel consumption.
Figure 21. Fuel consumption.
Wevj 11 00074 g021
Figure 22. The variation in the CO2 release.
Figure 22. The variation in the CO2 release.
Wevj 11 00074 g022
Figure 23. Tractor driving cycle with a straw tub grinder attached.
Figure 23. Tractor driving cycle with a straw tub grinder attached.
Wevj 11 00074 g023
Figure 24. Thermal engine torque.
Figure 24. Thermal engine torque.
Wevj 11 00074 g024
Figure 25. Electric motor torque.
Figure 25. Electric motor torque.
Wevj 11 00074 g025
Figure 26. Battery’s state of charge.
Figure 26. Battery’s state of charge.
Wevj 11 00074 g026
Figure 27. Fuel consumption.
Figure 27. Fuel consumption.
Wevj 11 00074 g027
Figure 28. The variation in the CO2 release.
Figure 28. The variation in the CO2 release.
Wevj 11 00074 g028
Table 1. Main notation used in this paper.
Table 1. Main notation used in this paper.
SymbolDescription
ρ a i r Air density
SxActive area of the vehicle for air resistance
CxCoefficient of drag for longitudinal axis translation
CzCoefficient of drag for translation on the vertical axis
WxStirring coefficient for the translation on the longitudinal axis
WzStirring coefficient for translation on the vertical axis
T 1 Front-wheel torque
T 2 Rear-wheel torque
r a d 1 Front-wheel radius
r a d 2 Rear-wheel radius
M c a r Mass of the vehicle
g Gravity
α s l o p e Slope angle in rad
α Angle of slope in %
V c a r Longitudinal speed of the vehicle
r v e h i Viscous coefficient of friction
F r o l l Rolling resistance force
R r o l l Coefficients of rolling friction
W z Vertical wind coefficient
V w i n d Wind speed
F t r a i l e r Resistive force from the connected tool
M t r a i l Trail mass capture
r i n 1 Rear-wheel inertia
r i n 2 Front-wheel inertia
F e n g i n e Driving force
BraketorqBraking torque required by the driver [Nm]
TorbrakdrvDriver brake control [Nm]
MaxtorqvehMaximum braking torque of the vehicle [Nm]
SommetorqTorque required by the driver [Nm]
TorvehdrvDriver acceleration control [W or Nm]
tmaxMelMaximum engine torque [Nm]
BraketorqBraking torque required by the driver [Nm]
sensorMErevSpeed of the electric motor [rpm]
TpullTorque of the electric motor for lifting [Nm]
sensorMESpeed of the electric motor [rad/s]
GearratioTransmission ratio between the electric motor and the vehicle
torMETorque control of the electric motor [Nm]
GainaccelGain on acceleration control
tmaxMelMaximum engine torque [Nm]
GearatioTransmission ratio between the electric motor and the vehicle
TorbrakBrake control
GainbrakesGain on brake control
MaxtorqvehMaximum braking torque of the vehicle [Nm]
tminMelMinimum engine torque
tmaxMthMaximum engine torque [Nm]
throttleMTThermal engine load
GainacceGain on acceleration control
ClutchsigClutch control
startengFiring: engine on/off
torMETorque control of the electric motor [Nm]
torbrakBrake control
maxtorqvehMaximum braking torque of the vehicle [Nm]
sensorMERotational speed of the electric motor [rad/s]
gainaccelGain on acceleration control
ω p a Rotational speed of the driven axle
p a r a t i o Drive axle transmission ratio
r a t i o Transmission ratio of the engaged gear
ω p Speed of rotation of the input shaft
ω s Speed of rotation of the secondary shaft
T s l i p Maximum Coulomb friction torque
ω r e l Relative rotational speed
d v e l Speed threshold of the rotating stick
T s Secondary shaft torque
n u Efficiency of the engaged gear
T l o s s Loss of torque
r a t i o The transmission ratio of the engaged gear
p a r a t i o Drive axle transmission ratio
T s Torque transmitted by the secondary shaft
v i s Coefficient of viscous friction on the flywheel
ω p Flywheel rotation speed
p o w e r P r Input shaft power ( Tp     ω p )
p o w e r P a Power to the drive axle ( Tpa     ω pa )
B M E P Average effective brake pressure [ bar ]
P M F Mean effective friction pressure [ bar ]
T o u t Engine output torque [ Nm ]
T f r i c Engine friction torque [ Nm ]
V Swept volume of the engine [ L ]
T m a x c   e t   T m i n c Maximum and minimum corrected torques
T m a x Maximum engine torque read from the file [ Nm ]
ρ r e f Reference ambient air density [ kg / m 3 ]
T a m b Ambient air temperature [ K ]
TrefReference ambient temperature [ degC ]
T m a x Maximum engine torque [Nm]
T f r i c ,   T h o t Frictional torque at high engine temperature [Nm]
T f r i c ,   T e m p Frictional torque at the current temperature [Nm]
T h o t “Hot engine temperature” parameter [degC]
T e m p Current engine or oil temperature [ degC ]
c o n t r e Actual fuel consumption
c o n s f i l e Value read from the fuel file
k s t a r t Overconsumption when starting the engine
e x p f c o n s Effect of cold temperature (user-defined expression)
φ Current equivalency report
φ f i l e Value read from the data file
e x p e q r a t i o Effect of cold temperature (user-defined expression)
d m e x h Exhaust gas mass flow rate [g/s]
c o n s Fuel consumption [g/s]
A F s t o Stoichiometric air/fuel ratio [null]
w Motor speed [rpm]
i d l e s p e e d Engine idling speed [rpm]
f l o w 0 Fresh air mass flow rate at idle [g/s]
T e n g Motor temperature [degC]
T h i g h High threshold for motor temperature parameter [degC]
T l o w Low threshold for motor temperature parameter [degC]
c o e f c o l d C O Correction coefficient for cold CO emission
c o e f s t C O Coefficient of CO when enriching the fuel when starting the engine
k s t a r t Overconsumption when starting the engine
e m i C O f i l e Value of CO
F i c h i e r   F M E P Average effective friction pressure read in the [bar] file
e x p f m e p c o r Temperature effect correction (user-defined expression)
T e x h a u s t f i l e Exhaust gas temperature read in the [degC] file
e x p e x h t e m p Temperature effect correction (user-defined expression)
c o e f h e a t w a l l Combustion heat ratio for wall heating
c o n s Fuel consumption [kg/s]
p c i Specific calorific value of the fuel [J/kg]
e x p h l o s s e s Correction of the effect of temperature
c o n s Fuel consumption [kg/s]
e m i H C HC emissions [kg/s]
e m i C O CO emissions [kg/s]
p c i Specific calorific value of the fuel [J/kg]
d m B G Burnt exhaust gas
d m u n b u r n e d a i r Exhaust unburned air mass flow rate (when φ ≤ 1)
Unburned fuelExhaust unburned fuel mass flow rate (when φ ≥ 1)
C p B G Thermal capacity of the flue gas set by the user [J/K/kg] for φ = 1
C p c a r b Heating capacity of the fuel set by the user [J/K/kg]
P a i r e Thermal capacity of fresh air and gas set at 1040.0 [J/K/kg]
n c y l Total number of engine cylinders
n c y l d e a c Number of deactivated cylinders
n c y l a c Number of activated cylinders
B M E P m a x Maximum BMEP read from the file and corrected for temperature effects
F M E P FMEP read from the file at current load and corrected for temperature effects
F M E P 0 FMEP at zero load
T m i n and T m a x Negative and positive torque corresponding to defined parameters
T m Couple [ Nm ]
ω Speed of rotation [ rad / s ]
P e l e c Electrical power [ W ]
SOCState of charge
DODDepth of discharge
CPVOOpen-circuit voltage
V2Potential at negative pole [V]
fixedvoltageFixed battery voltage [V]
fixedcellvoltageFixed cell voltage [V].
SealNumber of cells in series
CPVObatOpen-circuit voltage read from the file/expression for the whole battery [V]
OCVcellOpen-circuit voltage read from the file/expression for a cell [V]
RbatInternal resistance read from the file/expression for the whole drum set [Ohm]
RcellInternal resistance read from the file/expression for a cell [Ohm]
I3Battery current at the positive pole [A]
IcellCellular current [A]
CellNumber of cells in series
CellNumber of parallel cells
Table 2. Parameters of implements.
Table 2. Parameters of implements.
ImplementName of Implement Maximal VelocityThe Needed Force for Implements
Wevj 11 00074 i001Moldboard plow7 km/hDraft power
Wevj 11 00074 i002Bette Harvest8 km/hPTO + Draft power
Wevj 11 00074 i003Straw tub grinder0 km/hPTO power
Table 3. Brake command.
Table 3. Brake command.
Brake Command: Signal [0.1]Brake Command: Wheel Torque (Nm)
braketorq t o r b r a k d r v m a x t o r q v e h t o r b r a k d r v
Table 4. Torque command.
Table 4. Torque command.
Sommetorq
Acceleration signal command [0.1] t o r v e h d r v t m a x M e l b r a k e t o r q
Acceleration control: wheel power [W]if |sensorMErev| < 1.0 in torvehdrv > 0
t p u l l
if |sensorMErev| < 1.0 in torvehdrv ≤ 0
b r a k e t o r q
if |sensorMErev| ≥ 1.0
t o r v e h d r v / | s e n s o r M E | / g e a r a t i o b r a k e t o r q
Acceleration control: wheel torque [Nm] t o r v e h d r v g e a r a t i o b r a k e t o r q
Table 5. Torque required.
Table 5. Torque required.
TorMETorbrak
Torque requested by the driver is positive (acceleration requested) M I N ( s o m m e t o r q .
g a i n a c c e l ,
t m a x M e l )
0
Torque requested by the driver is negative or null (brake requested) and no gear is selected0 ( t o r M E . g e a r a t i s o m m e t o r q ) .
g a i n b r a k e s m a x t o r q u e h
Torque requested by the driver is negative or null (brake requested) and no gear is selected M A X ( s o m m e t o r q g e a r a t i o n ,
t m i n M e l )
t o r M e . g e a r a t i o s o m m e t o r q
g a i n b r a k e s m a x t o r q v e h .
Table 6. Motor control.
Table 6. Motor control.
Sommetorq
Acceleration command: signal [0.1] T o r v e h d r v . t m a x M t h b r a k e t o r q
Acceleration command: wheel power [W]Si |sensorMErev| < 5.0 and torvehdrv > 0
t p u l l . g e a r a t i o
Si |sensorMErev| < 5.0 en torvehdrv ≤ 0
b r a k e t o r q
Si |sensorMErev| ≥ 5.0
t o r v e h d r v / | s e n s o r M E | / g e a r a t i o b r a k e t o r q
Acceleration command: wheel power [Nm] t o r v e h d r v g e a r a t i o * b r a k e t o r q
Table 7. Motor control.
Table 7. Motor control.
ThrottleMTClutchsigStarteng
Torque requested by the driver is positive (acceleration requested). M I N ( s o m m e t o r q . g a i n a c c e l , t m a x M t h ) t m a x M t h 11
Torque requested by the driver is negative or null (acceleration requested).000 or 1 function of engdecel parameter
Table 8. Torque required.
Table 8. Torque required.
TorMETorbrak
Torque requested by the driver is positive (acceleration requested).00
Torque requested by the driver is negative or null (brake requested) and no gear is selected.0 ( t o r M E . g e a r a t i o s o m m e t o r q ) g a i n b r a k e s m a x t o r q u e h
Torque requested by the driver is negative or null (brake requested) and the gear is selected. M A X ( s o m m e t o r q g e a r r a t i o , t m i n M e l ) ( t o r M E . g e a r r a t i o s o m m e t o r q ) . g a i n b r a k e s m a x t o r q u v e h
Table 9. Torque required by the driver.
Table 9. Torque required by the driver.
Sommetorq
Acceleration command: signal [0.1] T o r v e h . t m a x M e l b r a k e t o r q
Acceleration command: wheel power [W]if |sensorMErev| < 1.0 and torvehdrv > 0
t p u l l
if |sensorMErev| < 1.0 and torvehdrv ≤ 0
b r a k e t o r q
if |sensorMErev| ≥ 1.0
t o r v e h d r v / | s e n s o r M E | / g e a r a t i o b r a k e t o r q
Acceleration command: wheel torque [Nm] t o r v e h d r v g e a r a t i o b r a k e t o r q
Table 10. Torque required.
Table 10. Torque required.
TorMETorbrak
Torque requested by the driver is positive (acceleration requested) M I N ( s o m m e t o r q .
g a i n a c c e l , t m a x M e l )
0
Torque requested by the driver is negative or null (brake requested) and no gear is selected0 ( t o r q M E . g e a r a t i o s o m m e t a r q ) .
g a i n b r a k e s m a x t o r q v e h
Torque requested by the driver is negative or null (brake requested) and no gear is selected M A X ( s o m m e t o r q g e a r a t i o ,
t m i n M e l )
( t o r q M E . g e a r a t i o s o m m e t a r q )
. g a i n b r a k e s m a x t o r q v e h
Table 11. Torque required.
Table 11. Torque required.
Sommetorq
Acceleration command: signal [0.1]. T o r v e h d r v .
t m a x M t h b r a k e t o r q
Acceleration command: wheel power [W]if |sensorMErev| < 1.0 and torvehdrv > 0
t p u l l
if |sensorMErev| < 1.0 and torvehdrv ≤ 0
b r a k e t o r q
if |sensorMErev| ≥ 1.0
Acceleration command: wheel torque [Nm] t o r v e h d r v g e a r a t i o b r a k e t o r q
Table 12. Torque control.
Table 12. Torque control.
TorMETorbrak
Torque requested by the driver is positive (acceleration requested) and the torque requested for the engine is lower than the optimum torque. M A X ( s o m m e t o r q .
g a i n a c c e l ,
t m i M e l )
0
Torque requested by the driver is positive (acceleration requested) and torque requested for the engine is higher than the optimum torque.00
Torque requested by the driver is negative or null (brake requested) and the gear is selected.0 ( t o r q M E . g e a r a t i o s o m m e t a r q ) . g a i n b r a k e s m a x t o r q v e h
Table 13. Engine control.
Table 13. Engine control.
ThrottlemtClutchsigStarteng
Torque requested by the driver is positive (acceleration requested) and the torque requested for the engine is lower than the optimum torque. M I N ( s o m m e t o r q . g a i n a c c e l t o r M E , t m a x M t h ) t m a x M t h 1
Torque requested by the driver is positive (acceleration requested) and the torque requested for the engine is higher than the optimum torque. M I N ( s o m m e t o r q . g a i n a c c e l t o r M E , t m a x M t h ) t m a x M t h 1
Torque requested by the driver is negative or null (brake requested).001 or 1 function of engdecel parameter
Table 14. Exhaust flow.
Table 14. Exhaust flow.
Positive ChargeEngine Brake
d m e x h = c o n s ( A F s t o ϕ + 1 ) d m e x h = f l o w 0 w i d l e s p e e d + c o n s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zahidi, Y.; El Moufid, M.; Benhadou, S.; Medromi, H. An Assessment of Low-Cost Tractor Motorization with Main Farming Implements. World Electr. Veh. J. 2020, 11, 74. https://0-doi-org.brum.beds.ac.uk/10.3390/wevj11040074

AMA Style

Zahidi Y, El Moufid M, Benhadou S, Medromi H. An Assessment of Low-Cost Tractor Motorization with Main Farming Implements. World Electric Vehicle Journal. 2020; 11(4):74. https://0-doi-org.brum.beds.ac.uk/10.3390/wevj11040074

Chicago/Turabian Style

Zahidi, Yassine, Mohamed El Moufid, Siham Benhadou, and Hicham Medromi. 2020. "An Assessment of Low-Cost Tractor Motorization with Main Farming Implements" World Electric Vehicle Journal 11, no. 4: 74. https://0-doi-org.brum.beds.ac.uk/10.3390/wevj11040074

Article Metrics

Back to TopTop