1. Introduction
The tractor is one of the most important agricultural machinery [
1,
2] and faces complicated working conditions. Different operations, such as plowing, deep scarification, winter tillage, rotary tillage, seeding, transplanting, harvest, transportation, ditching and so on, require tractors to meet different operational requirements [
3,
4,
5,
6]. The performance of speed regulating transmission system of tractors, such as transmission ratio, the variation range of transmission ratio, transmission efficiency and so on, has a great influence on the fuel economy, dynamic property and other key performance indexes of the full tractor. The HMCVT (Hydro-Mechanical Continuously Variable Transmission) can realize continuously variable speeds in a large range and also inherits the advantages of mechanical transmission’s high efficiency and hydraulic transmission’s high power [
7,
8,
9,
10] while overcoming the disadvantages, and thus, can meet the requirements of agricultural machinery like tractors better. Therefore, the HMCVT is applied widely to agricultural and engineering machinery, such as tractors, and the research on HMCVT has great significance [
11,
12,
13,
14].
A reasonable design is the primary premise of the application of HMCVT transmission parameters to agricultural machinery like tractors. A good design for HMCVT transmission parameters of the tractor needs to take the following aspects into consideration: (1) the requirement of full tractor’s driving and working conditions on HMCVT speed variation range; (2) the continuity of transmission ratio variation intervals of variable speed stages of multi-stage HMCVT, to avoid any discontinuity point in the variation process of HMCVT transmission ratio from the minimum value to the maximum value, reduce the impact of stage shifting and improve the ride comfort of the tractor in operation and the reliability of parts; (3) the maximization of work efficiency of HMCVT for energy conservation and emission reduction.
Currently, there is relatively little research on transmission parameter optimization design based on HMCVT speed regulation characteristics and efficiency characteristics (that is, the transmission efficiency changes with the system operating conditions, efficiency can be calculated from the transmission output power and input power ratio). Ince and Guler [
15] made detailed theoretical analysis and derivation of characteristic models for a new continuously variable speed system, including kinematic analysis, power flows and efficiencies at different stages. This study provides a valuable reference for the optimization design of a power-split infinitely variable power transmission system. Cheng and Lu [
16] measured the speed regulation characteristics of the hydraulic system of HMCVT. Four models of speed regulation characteristics were established and compared. In this study, the improved particle swarm optimization algorithm and the modified HMCVT speed-regulating characteristic model are applied to the optimal design of transmission parameters. This is helpful in improving the practical application effect of the optimized design results. However, this study ignored the effect of the efficiency characteristics of HMCVT. Cheng et al. [
17] used an improved genetic algorithm to optimize the transmission parameters of the HMCVT, which has a three-planetary alignment cooperative motion. This study takes into account the efficiency characteristics of HMCVT. However, the efficiency of the hydraulic system (mainly composed of pumps and motors) is assumed to be a constant value. Zhu et al. [
18] designed a multi-stage HMCVT by taking the ratio of each gear’s maximum and minimum output speed as the fixed value (namely, equal-ratio design method). The approximate method was also adopted by He et al. [
19] and Zhang et al. [
20] (similar to the equal-ratio design method, the difference between each gear’s maximum and minimum output speed is taken as the fixed value in this study) and other studies. These methods have certain limitations in terms of design flexibility. Additionally, it is often difficult to take into account the efficiency characteristics of the system. This is not friendly in terms of energy conservation.
To sum up, the current research on HMCVT design has some shortcomings, including the following aspects: (1) The design only considers the speed regulation characteristics of HMCVT. This results in limited transmission efficiency during HMCVT operation. (2) The design method is not flexible. In the classical design method, the ratio of the maximum and minimum output speed of each gear is fixed by solving equations, or the difference between the maximum and minimum output speed of each gear is fixed. (3) The changing characteristics of hydraulic system efficiency are ignored. The hydraulic system efficiency is often set to a fixed value in the study. This leads to a certain error between the analysis results and the actual system. (4) The working characteristics of agricultural machinery in the whole life cycle affect the performance of HMCVT. Agricultural machinery often operates in certain speed ranges when traveling or operating. The design of HMCVT can further improve the working efficiency of the whole agricultural machinery by matching the service characteristics of agricultural machinery. No reports have been found on the introduction of full life cycle working characteristics of agricultural machinery in design. In addition, the heuristic intelligent optimization algorithm is used to optimize the HMCVT transmission parameters. The decision variables are many, as well as the objective functions are many and nonlinear. Therefore, it is difficult to optimize the design based on the characteristics of HMCVT transmission ratio change, the working characteristics of agricultural machinery in the whole life cycle, and the maximization of HMCVT transmission efficiency.
To solve the problems above, this paper proposes a transmission parameter optimization design method based on the verified speed regulation characteristic model and the verified hydraulic system efficiency characteristic model aiming at improving the performance (mainly referring to speed regulation characteristics and efficiency characteristics) of tractor HMCVT. The research in the paper considers a type of 5-stage tractor HMCVT as the object. According to the pump-motor efficiency model, the paper builds an HMCVT efficiency model by using the engaging power method. The paper uses the I-GA improved in four aspects (i.e., avoiding any super individual in the initial population, the adaptive scaling of population size in the iteration process, the adaptive changes of crossover probability and the adaptive changes of mutation probability) to optimize and improve HMCVT performance in speed regulating and efficiency characteristics. In the optimization, the paper considers the following three points as the comprehensive objective: the HMCVT speed regulating characteristics and efficiency characteristics matching the driving and working requirements of the tractor, the continuous speed regulating characteristic of HMCVT and the maximization of mean transmission efficiency of HMCVT. According to the HMCVT transmission parameter optimization design method for tractors proposed in the paper, the optimal matching design of HMCVT transmission parameters and performance improvement can be realized by obtaining the efficiency data samples of the “pump-motor” system selected in the HMCVT design and combining with the requirements of tractor operation and driving conditions.
4. The New HMCVT Optimization Design Method
The paper considers the HMCVT of the tractor as the design object. Suppose the maximum speed of a tractor running in the fields is 32 km/h, the radius of the wheel is 0.9 m, engine speed changes in the range of 750~2300 r/min, and the other train transmission ratio is 28. The calculation of tractor speed refers to the calculation formula of vehicle speed [
25,
26,
27]. As the conceptual design phase for HMCVT is studied in this paper, tire slip is ignored. Then,
, the minimum transmission ratio of HMCVT to be designed, can be calculated and obtained, which is 0.87 (the minimum transmission ratio in stage HM
4). Next, suppose the minimum steady speed in stage HM
1 required by a tractor in field operations is 2 km/h (stage H
0, i.e., the hydraulic transmission system, offers the speed for the tractor to start up and the speed less than 2 km/h). Additionally, to meet the tractor’s field driving requirement to large torque, the engine’s working speed is about 1500 r/min. Therefore, we can calculate and obtain that
, the maximum transmission ratio of stage HM
1 of HMCVT, is 9.09.
Different working conditions have different requirements for the tractor’s speed. For instance, some rotary tillage and seeding tasks of the tractor are low-speed operations, while plowing, intertillage, and field transportation are medium-high-speed operations. The conventional transmission ratio design of the speed control system of tractors only considers the requirements of the full-vehicle dynamic property. However, the conventional HMCVT design generally adopts equal-difference or equal-ratio changing transmission design schemes for the convenience of design calculation. Therefore, the conventional design principle neglects not only the main working speed intervals in the life cycle of the tractor but also the influence of the transmission system’s efficiency characteristics on the full-vehicle energy saving of the tractor.
The paper proposes an optimization design method of transmission parameters considering the HMCVT efficiency maximization in the life cycle of the tractor. The steps are as follows:
Step 1. Analyze the statistical data of operating frequency of speed in the whole life cycle of the tractor.
Our design takes scholar Wang’s analysis of the research results of Resch and Renius for reference [
28]. In the whole life cycle, the tractor spends 61~68% of working time in the speed interval of 4~12 km/h and 15~25% of working time in the speed interval of 12~20 km/h. Therefore, the paper designs the percentage of time of each time interval of the tractor in the whole life cycle as follows:
Step 2. Calculate the variation range of transmission ratio corresponding to the high-frequency used speed interval of the tractor.
With the minimum working speed and the maximum working speed of the engine, calculate the minimum transmission ratios and the maximum transmission ratios of HMCVT corresponding to different speed intervals in calculation Equation (13). Suppose the minimum transmission ratio and the maximum transmission ratio of HMCVT are and , respectively, when . Suppose the minimum transmission ratio and the maximum transmission ratio of HMCVT are and , respectively, when .
Step 3. Calculate the weight of service time of each hydro-mechanical power dividing stage (i.e., HM stage) of HMCVT.
According to the design results of a group of transmission parameters, calculate the proportion of the transmission ratio of each HM stage of HMCVT to the transmission ratio interval in Step 2 in the current design results. Next, according to the proportion and combined with the proportion ratio of each section’s service time in Step 1, estimate the proportion of service time of each HM stage in the whole life cycle of the tractor. Then, consider the proportion as the service time weight
required by the optimization of the objective function in the subsequent optimization design of HMCVT.
in which
is the service time weight of the
kth HM stage;
is the proportion of transmission ratio of the
kth HM stage in interval
;
is the proportion of transmission ratio of the
kth HM stage in interval
;
is the proportion of transmission ratio of the
kth HM stage in other transmission ratio intervals;
is the rest of transmission ratio range of all HM stages in the intervals except
and
.
Step 4. Set objective function
with the comprehensive optimization objective of transmission ratio continuity, transmission ratio variation range matching required speed of tractor and service efficiency maximization in the whole life cycle of the tractor.
in which
is the variable to be optimized and
(in which
act as a whole to transmit power in the whole HMCVT transmission system, so they are combined with being one variable to be optimized in the optimization design;
is a reverse gear and the HMCVT has a consistent working principle for moving forward and backward, so the paper sets
in the optimization design;
is the transmission ratio value required by the purely hydraulic system stage in the case of starting up, and because the purely hydraulic system transmission ratio can realize continuously variable regulation by itself, the determination of
’s value only needs to match the maximum transmission ratio of stage HM
1, so the optimization design in the paper doesn’t take the parameter into account.);
and
are the weight coefficients of the objective function;
is the objective function realizing the continuity of transmission ratio and the variation range of transmission ratio matching the required vehicle speed of the tractor, as shown in Equation (16);
is the objective function realizing the service efficiency maximization of the tractor in the whole life cycle, as shown in Equation (17).
in which
is the displacement ratio of the variable pump;
is the minimum transmission ratio of stage HM
4 set in the optimization design of HMCVT. According to the design requirement,
should take a value of 0.87. However, to ensure that HMCVT’s speed regulating characteristics meet the maximum vehicle speed requirement after optimization design, and also for the expectation of getting a vehicle speed as high as possible, the optimization design stage in the paper takes the value of
less than 0.87.
The paper adopts the mean value of transmission efficiencies in working conditions of different engine speeds and different motor output end loads as the transmission efficiency of a fixed displacement ratio and writes it as
the mean transmission efficiency in the displacement ratio. The research in the paper calculates the ratio of
’s integral value of variation interval of displacement ratio to the variation length of displacement ratio and considers it as the man value of transmission efficiency of each HM stage. Additionally, the research considers the sum value of the mean value of transmission efficiencies of stages and considers the sum value of the mean efficiency of HM stages and the product of efficiency weight as the service efficiency of HMCVT with currently designed parameters. Besides, the optimization process adopted in the paper is a process searching the minimum value of the objective function, so the paper considers the loss value of the system deducting service efficiency as the objective function of
.
in which
is the mean value of efficiency of the
kth HM stage with displacement ratio of
(
Section 2 and
Section 3 give the efficiency models of HM stages of the HMCVT); Constant 2 presents the variable length of displacement ratio of each HM stage.
The purpose of using the objective function
’s weight coefficients
and
is to balance the orders of magnitude of
and
, and, on the other hand, reflect the degree of importance tendency between
and
in the optimization design. In the optimization design of the paper, it’s recommended to substitute the original design parameters (
Table 1) into Equations (16) and (17) to obtain the values of
and
through calculation and make
and
.
Step 5. Use a heuristic intelligent optimization algorithm for the optimization design of transmission parameters.
The heuristic optimization algorithm is widely used to solve problems in complex engineering projects [
29,
30,
31,
32,
33,
34], which can effectively solve nonlinear, multi-objective or multi-decision variable optimization problems. The case of optimization design in the paper has many variables to be optimized (a total of 11 variables to be optimized), and the calculation model of the objective function is complicated (i.e., it’s required to calculate the speed regulating characteristics in different working conditions of HMCVT for each individual in each iteration process of algorithm and consider the service time weight and efficiency of each HM stage in the whole life cycle of the tractor), so the research in the paper adopts the I-GA for the optimization design of transmission parameters. The GA has improved in the following four aspects: the super individual screening, the population size adaptive change, the adaptive changes of crossover probability and mutation probability. See reference [
17] for the specific improvement process and optimization results after improvements.
5. Results and Discussion
In the optimization design, the paper sets the I-GA’s initial population size to be 500, and, in the process of iteration, the maximum value of adaptive change of population size to be 1000 and the maximum number of iterations to be 100.
Figure 2 shows the iteration evolution curve.
According to
Figure 2, the HMCVT’s optimization process has a fast rate of descent in the initial period (in the first 20 iterations, the objective function decreases from 4.5593 to 1.1491, decreasing by about 74.80%), indicating the I-GA can avoid the phenomenon of prematurity well in the early period of optimization. Meanwhile, because of the introduction of population size adaptive change and crossover & mutation probability adaptive change, the optimization process still has an evolution ability in the middle and later periods to get out of the optimal local solution.
Table 1 shows the optimization design results of HMCVT transmission parameters.
Figure 3 shows the transmission characteristics of HMCVT after the optimization design.
According to
Figure 3,
after optimization has an error of about 0.62% with the expected value in design;
and
have an error of about 0.62%;
and have an error of about 0.98%;
and
have an error of about 1.94%;
, less than 0.87 the expected value in design (indicating the maximum vehicle speed meets design requirement).
After optimization, the transmission ratio characteristics of each HMCVT working stage (a total of 4 HM stages) change continuously with the variable pump displacement ratio, without discontinuity point or transmission ratio mutation phenomenon. The total transmission ratio of the transmission system can be continuously changed from 9.15 to 0.86, meeting the design requirements of the tractor in
Section 4 (the maximum speed required for the tractor to drive in the fields is 32 km/h, and the minimum steady speed of stage HM
1 is 2 km/h). Specifically, the transmission ratio varies from 5.44 to 9.15 (stage HM
1), 2.86 to 5.47 (stage HM
2), 2.83 to 1.68 (stage HM
3), and 1.65 to 0.86 (stage HM
4) for each HMCVT working stage. The transmission ratio of HMCVT presents a nonlinear change characteristic, which corresponds to Equations (1)–(4), and the degree of nonlinearity decreases with the increase of the HM stage.
Figure 4 shows the comparison of mean transmission efficiencies of HM stages of HMCVT (i.e., the mean value of efficiencies with different engine speeds and different motor output end loads in each displacement ratio condition. The calculation in the paper sets the engine speed to be 750~2300 r/min and the motor output end load to be 34.71~138.84 Nm. Meanwhile, the paper chooses 50 evenly spaced working conditions in the engine speed interval and the output end load interval. Therefore, we calculate the mean of transmission efficiencies in a total of 2500 working conditions in the case of one displacement ratio) before and after the optimization.
According to
Figure 4, after the optimization, the mean transmission efficiencies of HM stages of HMCVT in the displacement ratio working conditions (i.e., calculating the mean value of transmission efficiencies in a total of 2500 working conditions with a certain displacement ratio) all improve to a certain extent. For stage HM
1 and stage HM
3, the mean transmission efficiency is improved by 4.08% to the most and 1.72% on average. For stage HM
2 and stage HM
4, the mean transmission efficiency is improved by 1.64% to the most and 0.71% on average. As for the improvement result considering HMCVT service efficiency in the whole life cycle of the tractor (i.e., objective function
), the value after optimization is improved by 19.93% compared with the value before optimization.
Figure 5 (corresponding to stage HM
1 and stage HM
3) and
Figure 6 (corresponding to stage HM
2 and stage HM
4) show the gross transmission efficiencies of the system of HMCVT in different working conditions (for example, in the case that engine speed is 750~2300 r/min, motor output end load is 34.71~138.84 Nm, and the displacement ratio is in −1~1)
In order to further demonstrate the effectiveness of the proposed design method, it is compared with the traditional design method (i.e., the design method that ignores the efficiency characteristics of HMCVT). Taking reference [
16] as an example, the traditional design method adopts the improved particle swarm optimization algorithm to optimize the design of HMCVT. The design results are shown in
Figure 7.
According to
Figure 7, the design results of the traditional design method can well meet the requirements of HMCVT speed regulation characteristics. However, the HMCVT efficiency characteristics are ignored in the design process, which leads to relatively poor efficiency characteristics after optimization design and high energy consumption of the system. For stage HM
1 and stage HM
3, the mean transmission efficiency is reduced by 24.30% to the most and 13.05% on average. For stage HM
2 and stage HM
4, the mean transmission efficiency is reduced by 5.35% to the most and 2.21% on average. Considering the whole life cycle of the tractor, the design result of HMCVT service efficiency decreases by 20.42%.