Real-Time Co-optimization of Gear Shifting and Engine Torque for Predictive Cruise Control of Heavy-Duty Trucks

Fuel consumption is one of the main concerns for heavy-duty trucks. Predictive cruise control (PCC) provides an intriguing opportunity to reduce fuel consumption by using the upcoming road information. In this study, a real-time implementable PCC, which simultaneously optimizes engine torque and gear shifting, is proposed for heavy-duty trucks. To minimize fuel consumption, the problem of the PCC is formulated as a nonlinear model predictive control (MPC), in which the upcoming road elevation information is used. Finding the solution of the nonlinear MPC is time consuming; thus, a real-time implementable solver is developed based on Pontryagin’s maximum principle and indirect shooting method. Dynamic programming (DP) algorithm, as a global optimization algorithm, is used as a performance benchmark for the proposed solver. Simulation, hardware-in-the-loop and real-truck experiments are conducted to verify the performance of the proposed controller. The results demonstrate that the MPC-based solution performs nearly as well as the DP-based solution, with less than 1% deviation for testing roads. Moreover, the proposed co-optimization controller is implementable in a real-truck, and the proposed MPC-based PCC algo-rithm achieves a fuel-saving rate of 7.9% without compromising the truck’s travel time.


Introduction
Heavy-duty trucks have an irreplaceable role in freight transportation in the daily life of human beings.As a commercial transport, one of the main concerns for the heavy-duty truck is energy consumption [1].One way of reducing energy consumption is by the development of a smarter cruise control (CC) system [2].At present, CC system which releases the driver from operating either on the accelerator or brake pedal is widely equipped to reduce driver workload [3] for some heavy-duty trucks.However, the traditional CC system usually brakes under slope conditions; this mechanism is not conducive to energy saving.Predictive cruise control (PCC) technologies, which can avoid the unnecessary braking by allowing the vehicle speed to vary around the set speed, has attracted a tremendous attention for heavy-duty trucks, considering the increasing demand for fuel economy [4][5][6].As an optimal control problem, PCC uses the upcoming road traffic information to optimize the longitudinal driving strategies [7,8].
Given that dynamic programming (DP) can provide a global optimal control solution of the optimization problem of PCC, it is used as a benchmark to analyze the optimal trajectories [9].However, DP cannot be directly applied to real-time control because of the heavy computation burden (known as the curse of dimensionality) [10].Similar to other advanced driving assistance system, PCC should be run on vehicle control unit in real time.Biswas et al. divided the real-time controllers in energy management systems into three categories, namely, (i) precomputed control laws, (ii) online computed control laws, and (iii) blended pre-and online computed control laws [11].The precomputed control laws, such as rule-based and fuzzy logic controllers [12,13], need tremendous calibration of the control parameters.Traditionally, these approaches are less robust to model inaccuracies and the change of the driving conditions.An example of the blended pre-and online computed control laws is explicit model predictive control (MPC) [14].However, the explicit MPC involves the exponential growth of the number of control regions for the long predictive horizon and many constraints [15].
MPC is expected to generate online computed control laws of the real-time controllers [16,17].However, the considerable computational burden from finding an optimal solution restricts the application of MPC in real-time manner [18,19].If the optimization problem can be formulated as a linear convex quadratic programming (QP) problem, the fast computation time can be achieved [20].Unfortunately, the PCC problem is a nonlinear and constrained optimization problem.Finding an analytical solution for such problem is often impossible or extremely difficult.To avoid solving complex optimization problems directly, some hierarchical optimization methods are utilized, for example, double-layer scheme for energy-optimal braking control [21], bi-level MPC for eco-driving applications [22].
At present, many studies are focused on numerical solution methods instead of the analytical solution for MPC problems [23,24].These numerical approaches fall into two broad categories [25], namely, (i) direct methods, and (ii) indirect methods.In the direct method, the original optimization problem is transcribed to a nonlinear programming problem by parameterizing the state and control input, while in the indirect method, the original optimization problem is converted into a boundary-value problem by using optimality conditions [26].Pontryagin's maximum principle (PMP) can be used as an optimal control theory to obtain the optimality conditions [27].Then, the shooting method is used to solve a boundary value problem by guessing an initial value of the costate variable.In Ref. [28], the authors used the indirect shooting method to derive an online optimal control law of the gearshift command for multispeed electric vehicles.Many heavy-duty trucks are equipped with automated manual transmission with discrete gear ratios.Thus, co-optimization of gear shifting and engine torque for PCC problem involves a mixed-integer programming problem [29].Considering that finding a fast optimal solution for the mixed-integer programming problem is extremely challenging, two sub-step optimization method is proposed by converting an original optimization problem into two sub-optimization problems in Ref. [30].
Based on the above previous works, a real-time implementable PCC is proposed to simultaneously optimizes engine torque and gear shifting for heavy-duty trucks.
The main contributions of this study are devising a realtime implementable co-optimization method and deriving an MPC-based solution.The MPC-based solution performs nearly the same as a global optimal control solution for the PCC optimization problem.The procedure to achieve this goal is summarized as follows.First, the PCC problem for heavy-duty trucks is formulated as a nonlinear MPC with discrete control input of the gear shifting command.Second, after simplifying that gear shifting only occurs in the first step of the predictive horizon, a fast solution is derived using PMP and an indirect shooting method, in which the state constraints on the truck speed and gear position are treated in practical implementation approaches.Third, hardware-in-theloop (HIL) experiments on a embedded computer and a vehicle control unit are carried out to show the real-time capability of the proposed co-optimization controller.Fourth, real-truck experiments are conducted to evaluate the fuel-saving performance of the proposed MPC-based controller, and a fuel-saving rate of 7.9% is achieved.
The remainder of this paper is organized as follows.The PCC system and system model of heavy-duty trucks are presented in Section 2. The optimization formulation of the PCC problem is discussed in Section 3. The solution methods are derived based on MPC and DP in Section 4. The proposed PCC solutions are evaluated through simulation, HIL and real-truck experiments, and some main results are represented in Section 5 and Section 6.Finally, the conclusions are provided in Section 7.

PCC System and System Model
In this section, the fuel-saving reasons of PCC system are firstly demonstrated by comparison to the conventional CC.Then, the longitudinal dynamic model of a heavyduty truck is constructed.Thereafter, the fuel consumption rate is fitted by a 2D polynomial function.

Fuel-Saving Reasons of PCC System
The conventional CC is designed to maintain the driver's set speed.It can be regarded as a speed tracking problem.PCC uses the information on the road profile to optimize the velocity profiles to reduce fuel consumption.It can be regarded as a speed optimization problem with upcoming road information.The difference between PCC and CC is shown in Figure 1.Compared with CC, the main fuel-saving reason for typical uphill and downhill conditions can be summarized as follows.
I. When approaching an ascent, the truck with PCC accelerates in advance, but the truck with CC does not have any uphill preparation.When climbing the ascent, the engine for CC works in a low-speed and high-torque region, and the truck often gains more torque by downshifting.The truck with PCC can maintain a high gear due to the early acceleration.Because working in a high gear is conducive to fuel saving, PCC saves more fuel than CC.II.When approaching a descent, the truck with PCC decelerates in advance, but the truck with CC does not have any downhill preparation.When on the descent, CC often brakes the truck to avoid excessive speed, while PCC can reduce the brake operation due to the early deceleration.Brake dissipates the energy in the form of heat; thus, PCC saves more fuel by avoiding unnecessary braking.

Longitudinal Dynamic Model
The vehicle under discussion is a heavy-duty truck.Its power comes from an internal combustion engine (ICE) mated to an automated manual transmission (AMT).Engine torque, brake torque, and gear number are the control inputs to determine the longitudinal motion of the truck.Traditionally, the engine and the brake system do not work simultaneously.The engine works in the traction mode, and the brake system works in brake mode.When the truck is in the brake mode, engine drag occurs, and the engine drag resistance force is equal to the minimum (negative).

Traction Mode
The longitudinal dynamic motion of the truck can be described using the truck distance s, the truck speed v and the longitudinal acceleration a, as follows: The force diagram of the heavy-duty truck moving on a slope is shown in Figure 2. According to the force balance, the engine traction force is the sum of various external forces, as follows: (1) where F t is the traction force, F f is the rolling resistance, F w is the aerodynamic friction, F i is the gravitational resistance, and F j is the acceleration resistance.
The discrete equation of the engine traction force is as follows: where η t is the combined transmission efficiency, which accounts for both the efficiency of the transmission and the main gearbox, i f is the final drive ratio, G(k) is the gear number at step k, i g (G(k)) is the transmission ratio in gear G(k), T e (k) is the engine torque at step k, and r w is the wheel radius.The transmission ratio i g is determined by the selected gear G(k).
The discrete equations of the resistance force are given as follows: where m is the mass of the truck, g is the gravity constant, f is the rolling resistance coefficient, α(k) is the slope angle at step k that varies along with the position of the truck, C D is the aerodynamic drag coefficient, A is the frontal area, ρ is the air density, and δ is the lumped rota- tional inertia coefficient.
The discrete equation of the longitudinal dynamic motion in the traction mode can be expressed as follows: where t is the time interval from step k to step k + 1. (2)

Brake Mode
According to the force balance, the longitudinal dynamic motion in the brake mode is given as follows: where F drag is the engine drag resistance force, and F b is the braking force.
The discrete equation of the longitudinal dynamic motion in the brake mode can be expressed as follows:

Gear Shifting
The transmission system for the investigated truck is 12-speed AMT.A gear shift input u g (k) is introduced to describe the gearshift operation, as follows: Gear skipping in the gearshift operation is not allowed.The gearshift operation only includes downshift, hold, and upshift.Thus, the gear shift input u g (k) has three entries, u g (k) ∈ {−1, 0, 1} , which respectively represent downshift, hold, and unshift.

Fuel Consumption Rate
Fuel consumption rate (g/s) is used to model the fuel consumption of ICE.Traditionally, the engine manufacturer provides a map of brake specific fuel consumption (BSFC) (g/kW•h).The fuel consumption rate can be obtained from the BSFC map by unit conversion.For a convenient mathematical formulation, the fuel consumption rate Q t is fitted by a 2D polynomial function where ι i,j represent the fitting coefficients, and ω e is the engine speed.The fitting performance of the fuel consumption rate is shown in Figure 3. Black dots in the upper subfigure of Figure 3 are experimental data points, which are obtained by a large number of steady state experiments.The bottom subfigure of Figure 3 shows the model residuals of the 2D polynomial function (9).The root mean square error (RMSE) is 0.1493.
The relationship between the engine speed ω e and the vehicle speed v(k) can be described as follows: (6)

Problem Formulation
This section formulates the optimization problem of PCC.PCC allows the speed of the truck to vary around the driver's set speed.The speed variation range provides opportunity to reduce fuel consumption by using the upcoming road slope information.PCC can be regarded as a speed optimization problem.Compared with our previous work in Refs.[16,31], the optimal control input of the gear selection is considered.

Objective Function
Evidently, the goals of PCC are to track the driver's set speed and to reduce the fuel consumption.Thus, the problem of the PCC can be described as a nonlinear multiobjective optimization.The performance index to be minimized is formulated as follows: with where N is the predictive horizon, v ref is the driver's set speed, and κ 1 and κ 2 are weighting coefficients.The item Q t (k)�t is penalized to minimize the fuel consumption.
The item v(k) − v ref 2 �t is penalized to minimize the tracking error of the driver's set speed.The terminal constraint (v(N ) − v ref ) 2 is penalized to guarantee that the vehicle reaches the reference speed at step N. The terminal constraint ϕ[x(N ), N ] is a soft constraint with respect to the terminal state v(N).

System Dynamics
The state equation of the discrete system is given as follows:

Capacity Constraints
When the engine is working, the engine torque T e (k) has a capacity constraint, expressed as follows: where T e,l and T e,u are the minimal and maximal torques, respectively.PCC allows the speed of the truck to vary within a specified range.Thus, the speed of the truck has a band constraint, where v min and v max are the minimal and maximal speeds, respectively.
The truck under investigation has a 12-speed transmission.The transmission ratio i g (G(k)) in different gears is designed to meet not only the maximum climbing slope, but also the maximum speed.The lowest gear is related to the maximum climbing slope, and the highest gear is related to the maximum speed.The applicable gear ratio for a specified truck speed is constrained by the engine speed range because the engine has a speed range.According to Eq. ( 10), the constraints of the applicable gears G(k) can be expressed as follows: where ω e,min and ω e,max are the minimal and maximal engine speeds, respectively.To determine the suitable ( 13) (14) T e,l (ω e (k)) ≤ T e (k) ≤ T e,u (ω e (k)), (15)  gears for the current vehicle speed, we can use the given velocity v(k) and Eq. ( 16), along with the maximum speed ω e,max , to calculate the transmission ratio i g (G(k)) of the large value, which corresponds to the low gear.Similarly, by combining v(k) with the minimum speed ω e,min , we can calculate the transmission ratio i g (G(k)) of the small value, corresponding to the high gear.These calculations help us determine the range of the transmission ratio, and subsequently, we can check if the transmission ratio of gears 1 to 12 falls within this range.On this basis, we can determine which gears are appropriate for the present vehicle speed.As a result, only some of the twelve gears can satisfy the constraint ( 16) for a given vehicle speed.

Formulation Summary
In summary, the PCC problem is formulated as a constrained nonlinear optimization problem, written as follows: In the formulation of Eq. ( 17), the optimal control inputs are the engine torque and the gear shifting command, i.e., u(k) = T e (k) u g (k) T .Note that the engine torque is a continuous variable and the gear shifting is a discrete variable.The actuator system of the gear shifting u g (k) is a transmission control unit, and the actuator system of the engine torque T e (k) is an engine control unit (ECU).
Traditionally, the ECU receives the position signal of the accelerator pedal, then a two-dimensional map of the pedal position and the engine speed is used to determine an engine torque command.The engine responds to the torque demand via many control loops, such as the air/fuel ratio feedback loop and the ignition angle feedforward loop.The ECU under investigation can directly receive the torque demand by a controller area network.Thus, the engine torque demand, instead of the pedal position, is selected as the optimal control input.
The system states are the truck distance, the truck speed, and the gear position, i.e., x(k The optimization problem is subject to state and input constraints.

Solution
Bearing the real-vehicle implementation in mind, this section proposes a real-time implementable solution for PCC optimization problem with consideration of gear shifting.First, a real-time solution is derived based on PMP under the framework of MPC.Then, a nonreal-time solution is given based on DP as a benchmark.

Solution Method Based on MPC
In this part, an indirect method is discussed firstly on the basis of PMP for optimization problem (17).Then, the implementation of the fast solver is summarized.Thereafter, the treatment on state constraint and gear shifting implementation are discussed from an engineering implementation perspective.

Pontryagin's Minimum Principle
The PCC optimization problem considering gear selection is a challenging task due to the discrete nature of gears.To obtain rapid solution method, a simplification that the gear position in the predictive horizon remains unchanged is made.In other word, the gear shifting only occurs in the first step of the predictive horizon, and the gear position does not change for the remaining steps of the predictive horizon.Thus, the gear position G(k) in the predictive horizon has three states, namely, downshift, hold, and upshift.
Benefit from this simplification, only two of the system states need to be well treated for finding the optimal solution of the optimization problem (17), i.e., states of the truck distance and the truck speed.Following the idea of Ref. [32], the Hamiltonian function is defined by the following: where 1 (k + 1) and 2 (k + 1) are the costate variables with respect to states s(k) and v(k), respectively.
According to the statement of Pontryagin's Minimum Principle, the optimal control input u * (k) , the optimal state trajectory x * (k) , and the corresponding costate T must minimize the hamiltonian function so that ( 18) for all steps, k = 1, 2, • • • , N and for all permissible con- trol inputs u ∈ .In addition, the update equation for costate variables must be satisfied.The boundary conditions of the costate variables are as follows: Given that the slope angle α(k) varies slowly along the truck distance s(k), i.e., ∂α(k) ∂s(k) ≈ 0 , the update equation for the costate 1 (k) becomes as follows: In addition, the boundary condition of the costate 1 (k) is as follows: Thus, the following relationship can be obtained: for all steps k = 1, 2, • • • , N .The Hamiltonian func- tion (18) then becomes as follows: For clarity of expression, the Hamiltonian is reformulated by a 2D polynomial using control variable T e (k) and state variable v(k) as follows: (19) where By substituting Eq. ( 26) into Eq.(20b), the update equation for the costate 2 (k) becomes as follows: Using the expansion formula of ρ 4 (k) and ρ 5 (k) , the transformation relationship between the costate variables 2 (k) and 2 (k + 1) can be expressed as follows: From Eq. (21b), the boundary condition of the costate 2 (k) is as follows: We now turn our attention to treat the necessary condition (19) of the optimal control.The gear shifting command u g (k) is a discrete variable, while the engine torque T e (k) is a continuous variable.In the predictive horizon, the gear shifting command u g (k) is used only in the first step, and the gear position G(k) has three possible states.Thus, the three different gear shift operation (downshift, hold, and upshift) have three independent Hamiltonian functions, defined as where subscripts −1, 0, 1 represent downshift, hold, and upshift, respectively.The optimal control T * e (k), u * g (k) must minimize the Hamiltonian function H at every step instant, written as follows: (29) Given that each Hamiltonian function in Eq. ( 30) is the quadratic function of the optimal input T e (k) (see Eq. ( 26)), the explicit optimal solution can be solved to minimize the Hamiltonian function, as follows:

Implementation of the Fast Solver
It can be seen from Eq. ( 32) that if the optimal state variable v * (k) and the optimal costate variable * 2 (k) are known, the optimal control inputs T * e (k), u * g (k) can be (32) obtained.When the initial co-variate 2 (0) is given, the corresponding terminal costate 2 (N ) can be obtained step by step by using the update Eq. ( 28).If the obtained the terminal costate 2 (N ) satisfies the boundary condi- tion (29), then, the optimal control input for the initial co-variate 2 (0) is the optimal solution of the optimiza- tion problem (17).If the obtained terminal costate 2 (N ) does not satisfy the boundary condition (29), then, we need to adjust the initial co-variate 2 (0) .This descrip- tion is the so-called shooting method.The original optimization problem ( 17) is formulated as a problem of finding an initial value.Here, the bisection method is used to solve the problem of finding an initial value.Detailed information on using the bisection method can be found in our previous work [16].

Discussion on the Braking Torque Input and the Band Constraint of Vehicle Speed
In the solution derivation, the state constraint v min ≤ v(k) ≤ v max is inefficiently considered mainly because considering this constraint directly complicates the solution method, and the computational burden becomes huge.Instead of directly considering ths constraint, an engineering method is used to deal with the constraint.If the state constraint v min ≤ v(k) ≤ v max is violated, the control input is determined by a prescribed lookup table, i.e., By defining F req = F w + F i + F j , e 1 = v − v min and e 2 = v − (v max + ς) , the prescribed lookup tables are written as where ς is the speed tolerance, k p0 (e 1 ) and k i0 (e 1 ) are the controller gains of the traction mode, and k p1 (e 2 ) and k i1 (e 2 ) are the controller gains of the brake mode.When the state constraint v min ≤ v(k) ≤ v max is not violated, the fast solver is used to determine the optimal control input.

Discussion on Gear Shifting Implementation
Eq. ( 16) shows that the gear G(k) has twelve candidates which are constrained by the velocity v(k) and the engine speed ω e (k) .Therefore, the optimal control laws for downshift, hold, and upshift operations from the three ( 33) (34) Hamiltonian functions Eq. ( 31) should be met the gear constraint Eq. ( 16).If the upshift and downshift operations violate the constraint, then the gear position will hold the last step's value.In addition, a gear hold time, such as 2 s, is introduced to prevent frequent shifting.
Once the gear is shifted, the current gear position needs to be held for 2 s.

Solution Method Based on DP (Benchmark)
The DP algorithm is selected as a benchmark because it can provide a global optimal control.The illustration of the DP calculation is shown in Figure 4 for the optimization problem (17).The driving distance is divided into N equal segments.Every segment of the distance has a truck speed v(k).The truck speed v(k) is discretized as follows: (35) where σ is the discretization interval.
According to Bellman's optimality principle, the recursive equation of DP is given for the optimization problem (17), as follows: (36) where . Starting from the terminal state, the optimal solution and the corresponding value function can be obtained by a backward iterative procedure.Setting k = N − 1 , Eq. ( 36) becomes (37) The corresponding optimal input at step N − 1 can be obtained by the following argument, Then, setting k = N − 2, • • • , 1, 0 , the value functions and optimal control sequence {u * (N − 2), • • • , u * (1), u * (0)} can also be obtained by the backward iterative procedure.
Given an initial state x(0), the optimal control input u * (0) at step k = 0 is generated.Then, applying the opti- mal control input u * (0) on the system model ( 13) yields a new state x (1).For the new state x(1), the optimal control input u * (1) at step k = 1 is generated.Similarly, the optimal state trajectories and the optimal control policy can be generated by a forward iterative procedure.

Simulation Results
The simulations are conducted in a joint simulation environment to demonstrate the effectiveness of the proposed PCC scheme.The high-fidelity model of the heavy-duty truck is constructed using the commercial software Trucksim, and the proposed PCC is built in Matlab/Simulink.The important parameters of the truck are given in Table 1.More detailed information on parameters can be found in Ref. [33].The complete control architecture of the proposed predictive cruise control is shown in Figure 5.The reference speed is set by the driver.The vehicle speed is measured by speed (38) (39) sensor.The position and road-slope information are provided by GPS and HD map, respectively.The proposed controller receives this information and provides the co-optimization control law.Then, the actuators of engine control unit, electronic stability control system and AMT respond to the co-optimization control law by controlling the heavy-duty truck.
In the simulation validation, three methods, namely, PID controller, DP algorithm, and MPC algorithm, are used.PID controller represents the traditional constantspeed cruise.DP algorithm as the benchmark of optimal control problem provides a global optimal solution.The artificial and authentic roads are used to evaluate the control performance of the proposed predictive cruise controller.The cruise speed is set as 90 km/h to simulate a high-speed cruise condition.
The performance index of the fuel saving rate will be used to evaluate the performance of each control algorithm, which is defined as where f PID and S PID are the fuel consumption and travel distance of the PID controller, f DP,MPC and S DP,MPC are the fuel consumption and travel distance of the MPC controller or DP controller.

Artificial Road
A single ascent road is first used to conduct the simulation to show the fuel-saving mechanism.The initial speed of the truck is set to the driver's set speed, i.e., v ref = 90 km/h.The allowed range of the speed is set as [80 100] km/h.Simulation results for the single ascent road are shown in Figure 6. Figure 6 shows that the truck with MPC algorithm and DP algorithm accelerates in advance when approaching an ascent.Acceleration in advance can avoid excessive reduction in the truck speed for the uphill process.The times to downshift gear for MPC algorithm and DP algorithm are later than that of the PID controller.The proposed controller saves more fuel than PID controller because working in a high gear is conducive to fuel saving.
Similarly, the simulation is conducted with a descent road.The corresponding results are shown in Figure 7.It is can be seen from Figure 7 that the truck with MPC algorithm and DP algorithm decelerates in advance when approaching a descent, while the truck with PID controller does not have any downhill preparation.When on the descent, the truck with PID controller has brake operation, whereas the proposed PCC has no Table 2 summarizes the performance comparison for the three different controllers.Compared with the PID controller, the fuel saving rate for the PCC with MPC algorithm and DP algorithm is approximately 1% for the ascent condition, approximately 17% for the descent condition.The deviation of fuel saving rates between MPC algorithm and DP algorithm is small.Note that the fuel saving-performance of Table 2 depends on how long the flat road before and after the slopes.The greater the proportion of flat roads is, the smaller fuel savings the PCC controller obtains.

Authentic Road
Two expressway routes are used to persuasively compare the control algorithms.The selected test roads are authentic roads, as shown in Figure 8.One route is from Yicheng service area to Chenan toll station (Route I), and the other is from Caopuyuan to Xijin (Route II).Their mileages are approximately 33 km and 40 km, respectively.Route I is an expressway section between two cities.Route II is an expressway across the city of Quzhou.The slope information on these two routes is provided by the map vendor of NavInfo.Simulation results are shown in Figure 9 and Figure 10.From Figure 9, a discernible discrepancy exists in vehicle speed between the DP algorithm and the MPC algorithm within the range of 25000 m to 30000 m.This disparity arises because the DP algorithm, serving as a global optimal controller, employs a more extensive set of global road slope information.Consequently, it initiates deceleration earlier when approaching a descent.In contrast, the MPC algorithm relies solely on the road slope information within the predictive horizon, resulting in less advanced deceleration.The findings within the 25000 m to 30000 m range align with those illustrated in typical descent roads in Figure 7.To further compare the fuel efficiencies for the three different controllers on authentic roads, the distribution of engine operating points for Figure 9 and Given that the engine does not inject fuel under zero or negative torque conditions, this portion of the data is not utilized as the statistical data.Judging from Figure 11 and Figure 12, it can be seen that the number of operating points in high efficiency area are near for DP algorithm and MPC algorithm, greater than that of the PID controller.The difference in the distribution of working points leads to the difference in fuel saving.
Table 3 summarizes the performance comparison for these two expressway routes.Compared with the PID controller, the fuel saving rate for the PCC with MPC algorithm and DP algorithm is approximately 4%−6% for the authentic roads.As revealed in Table 3, the control performance of the proposed MPC-based solution is nearly the same as that of the DP-based solution with less than 1% deviation and considerably better than PID controller.

Experimental Results
HIL and real-truck experiments are conducted to verify the real-time capability of the proposed co-optimization controller.The experiments are first implemented in an embedded computer with high computation capability, then in an in-vehicle control unit, and finally on a real truck driving on real roads.

Embedded Computer
The HIL test bench of the embedded computer is shown in Figure 13.The series of the embedded computer is  14.It can be seen that the simulation results and HIL experimental results are consistent.This agreement is primarily attributed to the utilization of identical algorithms, parameters, and initial conditions in both environments.Additionally, both the HIL platform and Matlab/Simulink software employ floating-point calculations, ensuring the absence of computational deviation.

Vehicle Control Unit
The HIL test bench of the vehicle control unit is shown in Figure 15.The vehicle control unit is an automotive-grade processing chip, SAK-TC233L, which is developed by Infineon.The vehicle control unit communicates with the upper computer via LAUTER-BACH Debugger and ValueCAN.The C source code of the proposed controller was generated via MATLAB toolbox C/C++ Coder, after which it was compiled by the HighTec's compiler.An upper computer was used to monitor and record the inputs and outputs of the vehicle control unit.
The proposed control method is integrated into the system as a functional component, constituting only a fraction of the entire framework.Consequently, the choice of the sampling time necessitates consideration of other functional components and the real-time characteristics of the entire system.A small sampling time may impede real-time execution, leading to system overrun prompts.Inheriting the original engineering requirements of the developed control unit, the sampling time of the vehicle control unit is set at 0.02 s.In this case, the proposed control method can be implemented in a real-time vehicle control unit, and the HIL experimental results of the vehicle control unit are shown in Figure 16.The test scenario is from an authentic road.The road slope information is provided by the map vendor.The set reference speed is 78 km/h.The initial speed differs from that in Figure 14.The divergence in initial speeds allows us to assess the adaptive performance of the proposed algorithms under various initial conditions, simulating real-world scenarios where vehicles may initiate tasks at different vehicle speeds.The results in both Figure 14 and Figure 16 show that the proposed controller works well for different initial speeds.As demonstrated in Figure 16 that the simulation results are nearly consistent with the results of the HIL test.The zoomed subfigure reveals slight deviations around 410 s.This is mainly due to the effect of data type conversion.Typically, the vehicle control unit does not support floating-point operations due to their large memory occupation and computational burden.In HIL test, SAK-TC233L employs fixed-point operations, whereas Matlab/Simulink software utilizes double-type floating-point operations.

Real Truck Implementation
To validate the capability of implementing the proposed algorithm in real trucks and its fuel-saving performance, a real-truck experiment is conducted in Hubei Province using the test truck shown in Figure 17.The experiment covers a test route from HeDianzhen to Gucheng toll station, as illustrated in Figure 17.The test distance is 100 km, and the speed limit for trucks is 100 km/h.During the experiments, a digital map is equipped to the test truck to provide up-coming road information.Moreover, the driver is required to intervene in emergency circumstances by applying brakes.The rule-based PCC developed by a component company is employed as the benchmark controller.
The experimental results, as depicted in Figure 18, show that the PCC algorithm outperforms the

Figure 1
Figure 1 Demonstration on the difference between PCC and CC

Figure 2
Figure 2 Force diagram of the truck moving on a road

Figure 3
Figure 3 Fitting performance: fuel consumption rate (upper sub-figure) and model residuals (bottom sub-figure)

Figure 4
Figure 4 Illustration of the DP calculation

Figure 5
Figure 5 Control architecture of the proposed predictive cruise control

Figure 6
Figure 6 Simulation results for the ascent condition

Figure 7
Figure 7 Simulation results for the descent condition

Figure 9
Figure 9Simulation results for the authentic road: route I

Figure 10
Figure 10Simulation results for the authentic road: route II

Figure 11 Figure 12
Figure 11 Distribution of working points in BSFC map: route I

Figure 13 Figure 14 Figure 14 Figure 15
Figure 13 HIL test bench of the embedded computer

Figure 16
Figure 16 HIL experimental results on the vehicle control unit for the authentic road

Figure 17
Figure 17 Experimental truck and test route

Figure 18
Figure 18 Real-truck experimental results

Table 1
Parameter values for the vehicle model

Table 2
Performance comparison for ascent and descent conditions

Table 3
Performance comparison for authentic roads Qingdao Automotive Research Institute of Jilin University, China.His research interests include intelligent connected vehicle, new energy vehicle and vehicle dynamics.Zhuo Yang received the B.Sc. degree in automation of electrical engineering machines from Jiangsu University, China, in 2008.He is currently working at Dongfeng Commercial Vehicle Technical Center, China, and is the section chief of the vehicle electronic control development office of the electronic control development department.His current work direction is vehicle electronic control system matching and application software development, including body domain (BDCU), vehicle domain (VDCU), and intelligent cabin (IDCU).Lin Zhang received the PhD degree in automotive engineering from Jilin University, China, in 2019.He is currently an Associate Professor at School of Automotive Studies, Tongji University, China.His research mainly includes the vehicle control in terms of safety, energy-saving, and intelligence, including vehicle dynamics and control, HEV control, and trajectory planning.Hong Chen received the B.S. and M.S. degrees in process control from Zhejiang University, China, in 1983 and 1986, respectively, and the PhD degree in system dynamics and control engineering from the University of Stuttgart, Germany, in 1997.In 1986, she joined Jilin University of Technology, China.From 1993 to 1997, she was a Wissenschaftlicher Mitarbeiter with the Institut fuer Systemdynamik und Regelungstechnik, University of Stuttgart, Germany.Since 1999, she has been a Professor with Jilin University, China, and hereafter a Tang Aoqing Professor.From 2015 to 2019, she served as the Director of State Key Lab-oratory Automotive Simulation and Control, Jilin University, China.She joins Tongji University, China, as a Distinguished Professor.Her current research interests include model predictive control, nonlinear control, and applications in mechatronic systems focusing on automotive systems.