Integrated Optimization of Structure and Control Parameters for the Height Control System of a Vertical Spindle Cotton Picker

Vertical picking method is a predominate method used to harvest cotton crop. However, a vertical picking method may cause spindle bending of the cotton picker if spindles collide with stones on the cotton field. Thus, how to realize a precise height control of the cotton picker is a crucial issue to be solved. The objective of this study is to design a height control system to avoid the collision. To design it, the mathematical models are established first. Then a multi-objective optimization model represented by structure parameters and control parameters is proposed to take the pressure of chamber without piston, response time and displacement error of the height control system as the optimization objectives. An integrated optimization approach that combines optimization via simulation, particle swarm optimization and simulated annealing is proposed to solve the model. Simulation and experimental test results show that the proposed integrated optimization approach can not only reduce the pressure of chamber without piston, but also decrease the response time and displacement error of the height control system.


Introduction
Cotton is one of the most important crops [1]. According to the data from International Cotton Advisory Committee, the world cotton production was more than 24 million tons in 2020. Increasing cotton production creates an urgent need for high efficiency harvesting machines. Due to the high efficiency in cotton harvesting, in recent years, the mechanical cotton pickers have gained more and more application. Among these mechanical cotton pickers, the most representative one is the horizontal spindle cotton pickers produced by John Deere and CaseIH, which can work not only reliably and stably, but also harvest with low cotton loss and less cotton leaf [2].
However, during the cotton harvesting process, the traditional horizontal spindle cotton picker is parallel to the cotton crop, which may lead to a low fiber quality of the cotton. To solve this problem, vertical spindle cotton pickers are developed in recent years [3]. It adopts a telescopic cotton removal method to avoid the reverse pulling by the cotton removal tray, which can efficiently improve the cotton fiber quality. Nevertheless, when using the vertical spindle cotton pickers, the spindles of cotton pickers may touch the soil. This may cause a deformation of spindles and decrease the cotton picking efficiency. Therefore, it is important to develop a height control system for the vertical spindle cotton picker.
Usually, there are two kinds of height control system, i.e., mechanical height control system and electrohydraulic height control system. In early agricultural practice, the mechanical height control system is widely used in the no-till seeding [4] and transplanting machines due to its simple structure, easy control logic and low cost [5][6][7][8]. Unfortunately, as the control accuracy and response time of the mechanical height control system are poor, it cannot be used in cases that need high control accuracy and rapid response. To solve this problem, electro-hydraulic height control systems are developed by researchers [9,10].
However, although the electro-hydraulic height control system has higher control accuracy than the mechanical height control system, the control of its hydraulic cylinder is very complex. In the past, many researchers have been engaged in this area and a variety of control methods have been proposed, such as proportion integration differentiation (PID) control algorithm [11], sliding-mode fuzzy control method [12], etc. Among these methods, PID control algorithm is widely used because of its features of simplicity, robustness and reliability [13]. Many improved PID algorithms are proposed based on the basic PID to enhance the probability of the suitable selection of the controller gains [14]. For example, in the work presented by Peng et al., they combined traditional PID with fuzzy control algorithm for planning trajectory of an electro-hydraulic robot [15].
A difficulty in a fuzzy PID control method is the design of fuzzy rules because it highly depends on the experience of experts or large amounts of experimental data [16]. Hence, the combinations with various natureinspired optimization algorithms (NIOAs) have a significant impact on the performance of PID control. Other researchers used the NIOAs to optimize PID parameters and fuzzy controller, which can improve the control efficiency of the system. Examples of the approaches can be found in the works reported in References [17][18][19][20], such as Grey Wolf Optimizer algorithm [17], Particle Swarm Optimization [18,19] and Simulated Annealing algorithm [20].
Summarizing the findings of the above discussion, it can be found that PID parameters have been optimized for hydraulic actuators to obtain accurate and precise control performance. However, the height control system of the cotton picker is a mechanical-electric-hydraulic coupled system. Its design should consider interactions between mechanical and electrical parts of the complete system. Several works about integrated optimization approach have been carried out in the past decades, such as parallel robots [21,22], mesh reflector deployable space antennas [23], pipeline leak detection [24] and multistage synchronous induction coilgun [25]. It is confirmed that there is a lack of work in integrated optimization approach for the height control system.
Generally, for a height control system, its performance is dependent on its mechanical structure design and controller design. In designing such a system, the mechanical structure cannot be optimized without considering its influence on control. Conversely, the control performance, which can be optimized by using control design methods, can further be improved by modifying the mechanical structure. Hence, optimization of mechanical structure parameters and control parameters in an integrated manner can further improve the performance of the height control system [26,27]. Despite the success achieved by integrated optimization approach, little work has been devoted to height control system. Motivated by those findings, this paper aims to bridge the research gap and make the following contributions. Firstly, an integrated optimization model of mechanical structure parameters and control parameters is proposed to take the pressure of the chamber without piston, response time and displacement error of height control system as the optimization objectives. Secondly, this model is solved by Optimization via Simulation method (OvS), Particle Swarm Optimization (PSO) and Simulated Annealing (SA). To our best knowledge, it is the first time that both mechanical structure parameters and control parameters are simultaneously considered in the proposed integration problem for the height control system.
The rest of the paper is organized as follows. Section 2 establishes mathematical models of this system. Section 3 presents an integrated optimization method. Section 4 gives the case study, followed by the conclusion in Section 5. Figure 1 shows the schematic of the proposed height control system, which consists of signal acquisition device, control system, hydraulic system and lifting mechanism. The signal of height is obtained via the angle sensor installed on the signal acquisition device. Then the signal of height is used to govern the electro-hydraulic servo valve. The hydraulic cylinder controlled by electrohydraulic servo valve is used as the power sector of lifting mechanism. With cylinder controlled by valve to lift the spindles of the cotton picker, it can avoid spindles colliding with stones on the cotton field.

Modeling of the Lifting Mechanism
The diagram of lifting mechanism of the height control system is shown in Figure 2. Line AB represents the hydraulic cylinder. Point B and point C are articulated on the frame of the spindle cotton picker corresponding to the installation position. L 1 represents the total length of the hydraulic cylinder, L 2 is the length of articulated rod AC, L 3 is the installation position on the frame, and L 4 is the distance of picking head centroid from articulated point A. L 1 can be further expressed as follows: where L 0 and L represent cylinder length and piston displacement of the hydraulic cylinder, respectively. For convenience of analyzing the dynamics of the lifting mechanism, it is simplified to an equivalent diagram as shown in Figure 3. Then, the dynamic model of the equivalent lifting mechanism is formulated as follows: where A 1 is the area of the end face of piston, A 2 is the effective area of the piston in cylinder chamber with piston, P 1 is the input oil pressure of the chamber without piston, P 2 is the pressure of chamber with piston, M is the equivalent mass of the lifting mechanism, t is the elapsed (1) time of the lifting mechanism from starting to move to the end of lifting, B is viscous damping coefficient of the piston, F L is the external load of the lifting mechanism. Equivalent mass M is calculated based on the theory that kinetic energy of the equivalent lifting mechanism is equal to the original lifting mechanism. Mass of hydraulic cylinder, piston and articulated rod AC is negligible compared to the mass of picking unit. Besides, velocity of hydraulic cylinder, piston and articulated rod AC is same order of magnitude with velocity of picking unit. So kinetic energy of hydraulic cylinder, piston and articulated rod AC is ignored. Equivalent mass M is calculated based on the following Eq. where m is the mass of picking unit, m = G g , G is gravitational force of picking unit, g is gravitational acceleration, ẋ and ẏ are two components of barycenter velocity of picking unit. x and y are calculated as follows: External load F L is calculated as shown in Eq. (6): Based on Eqs. (1)-(6), the dynamic model of the equivalent lifting mechanism is expressed as follows: The lifting mechanism of the height control system realizes the lifting function of the picking unit by extending the piston of the hydraulic cylinder. Thus, it can prevent the spindle from being bended in the collision with stones and metals. In Figure 4, point A is the articulation of the hydraulic cylinder and the picking unit. When the articulation is in the state of point A, the picking unit works in the normal position. The piston of the hydraulic cylinder does not extend and its stroke is zero. When the piston of the hydraulic cylinder moves to point A 1 , the picking unit is lifted. In this process, the profile height h can be calculated as: When the surface of cotton field changes, the angle variation of profiling slider is detected by the angle sensor. The angle sensor transmits the angle variation into the controller. The controller will output signal to make the piston rod of the hydraulic cylinder to extend out. The profile height h of the lifting mechanism is determined by the angle change of the angle sensor, which can be modelled as shown in Eq. (9): where l is the distance between angle sensor and touchpoint of the profiling slide, Δθ is the angle change of angle sensor.

Modeling of the Hydraulic Control System
When the hydraulic fluid is pumped into the cylinder of the height control system, the hydraulic pressure will force the piston to move back and forth. The displacement and speed of the hydraulic cylinder is determined by the opening size of the electro-hydraulic servo valve. Therefore, the height and rising speed of the picking unit are adjusted by the electro-hydraulic servo valve. In order to achieve precise position tracking control, a PID controller, which is widely used in various industries due to (8) its high reliability and robustness [28][29][30][31], is adopted to control the electro-hydraulic servo valve.
As shown in Figure 5, the input of the PID controller is ΔU, whose value is the difference between the actual displacement signal U 2 and the input displacement signal U 1 . After a proportional, integral and differential processing of the PID controller, the new displacement signal U is the output. Then the displacement signal U is converted into control current I of the electro-hydraulic servo valve by the servo amplifier. The transfer function of the servo amplifier is expressed in Eq. (10): where K A is the amplification factor of servo amplifier, ω A is the natural frequency of the amplifier.
The control current I of the electro-hydraulic servo valve controls the displacement of the valve core, thus resulting in the change of valve opening, and finally realizing the control of the flow q v of the electro-hydraulic servo valve. The transfer function of the electro-hydraulic servo valve is [32]: where K sv is the gain of the electro-hydraulic servo valve, ω sv is the natural frequency of the amplifier, ξ sv is the damping ratio of the electro-hydraulic servo valve.
When the output fluid of the electro-hydraulic servo valve flows into the hydraulic cylinder, the piston of the hydraulic cylinder will extend out of the cylinder block. So that the picking unit will be lifted to the position of the height required by the system. The transfer function of the hydraulic cylinder is: where K h is the gain of hydraulic cylinder, ω h is the natural frequency of hydraulic cylinder, ξ h is the damping ratio of hydraulic cylinder.
When the displacement L of the piston rod is detected by the sensor installed on the hydraulic cylinder, it will be converted into actual displacement signal. The difference between the actual displacement (10) signal and the input displacement signal are sent into the PID controller to adjust the displacement of the piston of the hydraulic cylinder. The transfer function of the displacement sensor is: where K se is the gain of displacement sensor, ω se is the natural frequency of displacement sensor.

Multi-objective Optimization of Height Control System
The height control system is a mechanics-electronicshydraulics coupled system that consisting of signal acquisition device, control system, hydraulic system and lifting mechanism. For such a system, pressure in the chamber without piston, response time and displacement error are the most important performance indicators. A reduction of the pressure in the chamber without piston will decrease the pressure required by the hydraulic system and hence improve the reliability of the hydraulic system. Moreover, minimization the response time and displacement error can improve response speed and precision of the height control system. As a result, the probability of the spindle cotton picker, whose off-ground height can be adjusted by the height control system, colliding to the ground will be reduced. Hence, in order to reduce pressure in the chamber without piston and minimize response time and displacement error of the hydraulic system, a multi-objective optimization is adopted in this paper. Moreover, for the height control system, the pressure of the chamber without piston is highly dependent on the mechanical structure parameters, i.e., inner diameter of hydraulic cylinder D, length of hydraulic cylinder barrel L 0 and the length of rod AC L 2 . The response time and displacement error is related to the combination of the mechanical structure parameters (i.e., inner diameter of hydraulic cylinder D, length of hydraulic cylinder barrel L 0 ) and control parameters. (i.e., proportion coefficient K P , integration coefficient K I and differential coefficient K D ). Hence, during the optimization process, the mechanical structure parameters and control parameters are both taken as the optimization variables. With the obtained optimal mechanical structure parameters and control parameters, a better height control system can be designed to reduce the probability of the spindle cotton picker colliding to the ground. In this section, a multi-objective optimization model is proposed firstly. Then an algorithm combined OvS, PSO and SA is proposed to solve the model.

Pressure of the Chamber without Piston
In a hydraulic system, the lower the pressure in the chamber without piston is, the lower the pressure required by the hydraulic system is. As a result, the reliability of the hydraulic system will be improved. Hence, in this paper, the pressure of the chamber without piston P 1 is set as an optimization objective. Oil return circuit is directly connected to the oil tank from the hydraulic schematic diagram, so P 2 = 0. According to Eq. (7), P 1 is expressed as follows:

Response Time and Displacement Error
For the height control system of the spindle cotton picker, response time t 0 and displacement error e are two vital performance indicators. With the aim to improve the profiling precision, reduce response time of hydraulic cylinder and decrease the probability of spindle colliding to the ground, the response time t 0 and displacement error e are chosen to be optimization objectives. For the sake of algorithmic convergence, the response time t 0 and displacement error e are merged into one normalized objective H, which can be seen in Eq. (15): where t opt and e opt are the response time and displacement error required by the height control system, ω 1 and ω 2 are the corresponding weights.
The hydraulic system of the height control system is complex with nonlinear characteristics. Its response time t 0 and displacement error e are influenced by multiple factors, such as the mechanical structure parameters and control parameters. During the optimization process, it is difficult to establish a precise mathematical model of the hydraulic system to acquire the response time and displacement error. Hence, a simulation model based on AMESim is proposed to acquire the response time and displacement error of the hydraulic system.
The AMESim is a kind of software serving for advanced modeling and simulation of engineering systems based on intuitive graphical interface. It is specialized in multi-disciplinary collaborative simulation especially for mechanics, electronics and hydraulics. As the height control system is a typical mechanical-electro-hydraulic (14) system, several libraries are involved to build the whole model, like hydraulic component design library for the hydraulic pump, planner mechanical library for mechanical attachments and control library for control operations. Meanwhile, the required parameters such as the inner diameter of hydraulic cylinder D and length of hydraulic cylinder barrel L 0 are input in the model. The simulation model of the height control system can be seen in Figure 6.

Variables
The decision variables considered for the integrated optimization problem can be classified into two categories. One is related to the structure and the variables needed to be decided are the inner diameter of hydraulic cylinder D, length of hydraulic cylinder barrel L 0 and the length of rod AC L 2 . Another is the PID control parameters, i.e., proportion coefficient K P , integration coefficient K I and differential coefficient K D .

Constraints
(1) The inner diameter of hydraulic cylinder D should be selected from the Standard GB/T 2348, which can be seen in Eq. (2) The overall length of hydraulic cylinder L 1 , length of rod AC L 2 and installation dimension L 3 are restricted to the geometry constraints of a triangle: where the overall length of hydraulic cylinder L 1 is the sum of the length of hydraulic cylinder barrel L 0 and travel range of piston L, installation dimension L 3 is constant. The length of hydraulic cylinder barrel L 0 is: where B is width of rod, B = (0.6-0.8)D, L A is the guiding distance of hydraulic rod, L A = (0.6-1.5)D, L M is sealing distance of hydraulic rod determined by sealing method, L C is the rest length of hydraulic cylinder barrel.
(3) The length of hydraulic cylinder barrel L 0 is within twenty times the length of its inner diameter: (4) The lifting height should be more than 100 mm when the hydraulic cylinder reaches its maximum travel range: where L max is the maximum travel range of hydraulic cylinder. (5) The PID parameters are restricted by the reasonable ones: where K Pmax , K Imax and K Dmax are the maximum values of K P , K I and K D . The widely used method for acquiring these values is simulation tests. During each test, one of the PID parameters is changed while another two parameters remain unchanged. If the response time and displacement error of the height control system are within a reasonable range, then the value of the changed parameter will be recorded. The tests will last for several times until all the maximum values of the three PID parameters are settled. Figure 6 AMESim model of the height control system

Optimization Model
Based on the analysis above, a multi-objective optimization model with the aim to minimize the pressure of chamber without piston, response time and displacement error of the height control system, is then formulated as shown in Eq. (22):

Optimization Solution
In this paper, the optimization model is solved via a simulation based optimization method. The response time and displacement error of the hydraulic cylinder are firstly acquired by the AMESim platform. Then (22) min [P 1 (L 0 , D, L 2 ), H(L 0 , D, K P , K I , K D )], 63, 80, 90, 100, 110, 125, 140, 160, 180], the simulation results are chosen to be fitness functions of the multi-objective optimization model. For the discrete structure parameters, i.e., inner diameter of hydraulic cylinder D, length of hydraulic cylinder barrel L 0 and length of rod AC L 2 , they are firstly searched using SA and the optimal control parameters K P , K I and K D for each set of structure parameters are then obtained through PSO algorithm. The reason for choosing SA and PSO is that the SA is capable of high searching efficiency and fast convergence speed for solving discrete optimization problem [33] and so does the PSO for continuous optimization problems [34,35]. The flowchart of the algorithm is shown in Figure 7 and the detailed procedure is described as follows.
Step 1: Set the parameters of SA and PSO. Generate initial solution x_current(D, L 0 , L 2 , K P , K I , K D ) and calculate objective function object (P 1 , H). Set x SA _current as x_current(1:3).
Step 2: Enter the outer loop of SA. Judge whether the temporary temperature temp is greater than the minimum temperature temp min set in Step 1 or not. If so, go to Step 3. Otherwise, go to Step 10.

Input SA and PSO parameters
Generate initial solutions x_current(D, L 0 , L 2 , K P , K I , K D )  Step 3: Judge whether the current iteration count g is less than the maximum iteration count g M or not. If so, go to Step4. Otherwise, go to Step 9.
Step 4: Randomly generate the adjacent solution x SA _new of x SA _current in the outer loop of SA and calculate the pressure object SA _new(P 1 ) of the chamber without piston in the hydraulic cylinder.
Step 5: Enter the inner loop of PSO and calculate the optimal PID parameters corresponding to the structural parameters x SA _new generated in Step 4. The details are as follows.
Step 5.1: Initialize the PSO, generate the initial solution x pos _current(K Pi , K Ii , K Di ). Invoke the model in AMESim and get the output results of response time and displacement error of hydraulic cylinder, namely the value of objective function object pso _cur rent(H).
Step 5.2: Update the velocity v id and position x id of every particle according to Eqs. (23) and (24) and generate the adjacent solution x pos _new(K Pi , K Ii , K Di ) of PSO. Invoke the model in AMESim and obtain the output results of response time and displacement error of profiling hydraulic cylinder, namely the value of objective function object pso _new(H).
where v T id and x T id are velocity and position of the dth (d = 1, 2 or 3) dimension of the ith particle in the Tth iteration, p id is the best position of each particle during the iterative process, p gd is the best position of all particles in the temporary iterative procedure, c 1 and c 2 are learning factors, r 1 and r 2 are two random numbers, ω is inertia weight.
Step 5.3: Update the iteration count of PSO, T = T + 1.
Step 5.4: Judge whether the current iteration count T is greater than the maximum iteration count T max or not. If so, go back to Step 5.2, otherwise go to Step 5.5.
Step 5.5: Generate the optimal PID parameters x pos _new that are related to the structure parameters from the outer loop of SA.
Step 7: Judge the dominating relationship of the solutions in SA. The details are as follows.
Step 7.1: Judge whether x_new is dominated by x_current or not. If not, set x_current as x_new and add x_new to Pareto Archive, otherwise, go to Step 7.2.
Step 7.2: Calculate the probability of acceptance Prob according to Eq. (25): , E(f(x), λ) = λ 1 P 1 + λ 2 H. λ 1 and λ 2 are set as 0.5. Step 7.3: Generate a random number Rand and judge whether Rand is less than Prob or not. If so, set x_current as x_new and add x_new to the Pareto Archive. Otherwise, retain x_current.
Step 8: Update the iteration count of SA, g = g + 1. Go back to Step 3.
Step 9: A commonly-used exponent cooling strategy is adopted in the outer loop of SA: where a is the cooling factor of SA and faster cooling speed is achieved with less value of a. Go back to Step 2.
Step 10: Output the Pareto Archive.

Case Study
To validate the proposed multi-objective optimization model and algorithm and gain a better understanding of the influence of structure parameters and control parameters on pressure of the chamber without piston, response time and displacement error, a series of experiments are carried out in this section.

Optimization Results
To verify the effectiveness of the proposed approach, three optimization cases are studied. Hydraulic parameters and structure parameters of the height control system for the cotton picker are shown in Table 1. Parameters of SA and PSO are given in Table 2. Optimization results are shown in Table 3.
Case 1: Height control system with original structure parameters and PID control parameters.
Case 2: Height control system with optimal PID control parameters but original structure parameters. Case 3: Height control system with optimal structure parameters and PID control parameters obtained by the proposed integration optimization approach.

Simulation Results under an Identical Height
According to the structural parameters of the height control system and displacement of the hydraulic cylinder, the pressure of the chamber without piston P 1 for the three cases is calculated and plotted in Figure 8.
From Figure 8, it can be seen that the maximum pressure of Case 1, Case 2 and Case 3 is 29.07 MPa, 29.07 MPa and 24.81 MPa, respectively. Compared with Case 1 and Case 2, the maximum pressure of Case 3 is reduced by 14.65%. Furthermore, it also can be found that the average pressure of Case 1, Case 2 and Case 3 is 11.82 MPa, 11.82 MPa and 12.31 MPa, respectively. The average pressure is slightly increased by 4.15% when compared with Case 1 and Case 2. Note that the maximum pressure and average pressure of Case 1 and Case 2 are the same. This is because the load of the hydraulic cylinder will not be changed if only control parameters are optimized but structural parameters remains unchanged.
To study the response characteristics of the three cases under the same height, an AMESim simulation is performed. During the simulation process, the output load of the hydraulic cylinder is set to be the maximum load.    Figure 8 Pressure of hydraulic cylinder for different cases As the step signal is the most common and useful reference signal for judging the response characteristics of a controlled system, during the simulation, the step signal is used to study the response characteristics of the height control system. The lifting height of spindle cotton picker is set to be 100 mm. The displacement of the hydraulic cylinder for the three cases is calculated by Eq. (8) and the piston displacement of Case 1, Case 2 and Case 3 is 80 mm, 80 mm and 72 mm. Simulation is carried out with step signal as shown in Figure 9. In Figure 10, the response curves of Case 1, Case 2 and Case 3 is given. The response time of the three cases are 0.98s, 0.75s and 0.48s respectively. Moreover, the displacement error of the three cases is 4.0%, 2.6% and 1.1%. The response time and displacement error of Case 3 are less than that of Case 1 or Case 2. This is because the displacement of the piston rod L for Case 3 is smaller than that of Case 1 and Case 2 under the same required height. Moreover, from Figure 10, it can be found that the control parameters has a great influence on the response time and displacement error of the height control system.

Simulation Results under an Identical Displacement
In order to verify the performance of the height control system, comparative simulations are made on the three cases with the same displacement of 80 mm. From Figure 11, it can be seen that the maximum pressure of the chamber without piston for Case 1 and Case 2 is 29.07 MPa, while the maximum pressure of Case 3 is only 24.81 MPa. Compared to Case 1 and Case 2, the maximum pressure of Case 3 is reduced by 14.65%. From Figure 12, it can be found that the response time of Case 3 is 0.62 s, which is reduced by 18.42% compared with Case 1  and Case 2. Besides, the displacement error of Case 3 is reduced by 1.4%. It can be found that significant improvement is achieved with the proposed approach. From the analysis above, it can be found that the optimal control parameters will vary with different structure parameters, optimization of structure parameters and control parameters in an integrated manner can further reduce the pressure, shorten the response time, and improve the positioning accuracy of the height control system.

Experimental Validation
According to the optimal structure and PID control parameters, as shown in Figure 13, a height control system for a cotton picker is constructed. In order to verify its performance, a practical test is carried out. Obstacles with the height of 100 mm, 80 mm and 50 mm are put on a smooth road to simulate the fluctuation of cotton field ground. The cotton picker travels through obstacles at a working speed of 5.8 km/h. The response time of the height control system is recorded as shown in Table 4. From Table 4, it can be found that the response time is 0.55 s, 0.43 s and 0.28 s under the height of 100 mm, 80 mm and 50 mm. All the response time is less than 0.93 s, which is the maximum response time to lift the picking unit to avoid colliding the obstacles at a speed of 5.8 km/h. This proves that the optimal structure and control parameters of the height control system can meet the requirements.

Conclusions
In this paper, a vertical spindle cotton picker has been presented to solve existing problems in the traditional spindle cotton picker. A mechanical-electric-hydraulic coupled height control system is designed for the vertical cotton picker to avoid collision between spindles and bumps on the cotton field, and improve the efficiency of cotton picking. The conclusions are as follows.
(1) Both structure and control parameters are optimized to improve terrain-tracking performance of the height control system. A multi-objective integrated optimization model is established, which takes the pressure, response time and displacement error as the optimization objectives. (2) The multi-objective integrated optimization model is solved by using the optimization method based on AMESim simulation. Experimental results show that the maximum pressure of the hydraulic cylinder, the response time and the displacement error can be effectively reduced by 14.65%, 18.42% and 1.4%, respectively. (3) Future work will be focused on the following aspects. Effects of the SA and PSO algorithm parameters on the optimization results of the height control system are not analyzed in the current study. It is expected that these works will be the future work. More improved optimization algorithms such as Artificial Bee Colony and Cuckoo Search algorithm [36][37][38] will be proposed to solve the proposed multi-objective optimization model.