# Drivability Improvement Control for Vehicle Start-Up Applied to an Automated Manual Transmission

Authors: Danna Jiang; Ying Huang; Zhe Zuo; Huan Li
Journal: Shock and Vibration
 (2017)

---

## Abstract

Drivability is the key factor for the automated manual transmission. It includes fast response to the driver’s demand and the driving comfort. This paper deals with a control methodology applied to an automated manual transmission vehicle for drivability enhancement during vehicle start-up phase. Based on a piecewise model of powertrain, a multiple-model predictive controller (mMPC) is designed with the engine speed, clutch disc speed, and wheel speed as the measurable input variables and the engine torque reference and clutch friction torque reference as the controller’s output variables. The model not only includes the clutch dynamic, the flexible shaft dynamic, but also includes the actuators’ delay character. Considering the driver’s intention, a slipping speed trajectory is generated based on the acceleration pedal dynamically. The designed control strategy is verified on a complete powertrain and longitudinal vehicle dynamic model with different driver’s torque demands.


---

## Body

## 1. Introduction

Automated manual transmission (AMT) is widely used in modern vehicles especially in trucks due to the advantage of its low weight, low cost, and high efficiency just like a manual transmission (MT) and driving comfort without directly operating the clutch and gearshift just like an automated transmission (AT) [1]. Another transmission using a dual clutch system also has the advantage of low weight, high efficiency, and driving comport, but it is not as simple or as cost-effective as the AMT. Therefore, AMT is considered as an inexpensive add-on solution for classical MT [2].In AMT, the start-up or gearshift phase is managed by an actuator driven clutch which is controlled by the transmission control unit (TCU) with communication to engine management system (EMS) for coordination control. It can be acceptable to the customers provided that the control strategies are able to limit the variation of the vehicle acceleration during vehicle start-up or gearshifts. That is one of the most important aspects of vehicle drivability [3]. Therefore, the drivability control in the clutch control is the key factor to guarantee the AMT performance.Many solutions for dry clutch control in AMT have been presented in recent years. There are many different approaches proposed for dry clutch control in the literatures, for example, the classical control [4], fuzzy control [5], linear-quadratic based optimal control [6–8], decoupling control [9], robust control [10], predictive control [2, 11], and the hybrid control [6, 12].However, the clutch control is expected to satisfy the different and sometimes conflicting objectives, such as fast and smooth start-up. Hence, optimization based algorithm becomes a potential solution for this problem. Because the clutch dynamics in slipping and engagement phases are typical nonlinear processes, nonlinear based control methods are also needed for this problem. In this paper, we choose multiple-model predictive control strategy to design the clutch controller algorithm considering that it can handle multivariable control problem naturally and deal with nonlinear system with switching the subcontrollers. And the calculation speed is higher than hybrid predictive control strategy when the system has high order.Unlike other researches of clutch control, this paper aims at drivability enhancement in the start-up phase of a heavy-duty truck with an automated manual transmission. Not only the slipping time and jerk considered in the controller design but also the torsion vibration and actuator delays which have bad effects on the drivability of the vehicle and driver’s intentions which directly affect the driver’s feeling are considered as well.The approach proposed in this paper deals with the clutch control of multiphases defined by the clutch position. In different phases, different control objectives are used. Torsion vibration which is associated with drivability and comfort is selected as one of the control objectives in every phase. In clutch going to engagement slipping phase, the other control objective is the driver’s intention to decide slipping speed trajectory tracking. A hierarchical approach is used. The high level is torque reference control for engine and clutch actuator, considering the modern engine controller widely adopting torque-based method. Three angular speeds are measured as the controller inputs. They are the engine speed, clutch disk speed, and the wheel speed. The controller outputs are the optimized engine torque requirement and the clutch friction torque requirement. In this paper, we just focus on the high level of the clutch control. And the low level controller is assumed as an ideal one.This paper is organized as follows. In Section2, the structure of an automotive driveline is discussed, and a simplified driveline is introduced, with the dry clutch, flexible drive shaft and actuators, engine, and the clutch pneumatic actuator dynamic considered in the model. In Section 3, the different control objectives for different phase of the AMT are stated. And driver’s intention is also considered in the control objectives. In this section, the control scheme based on constrained MPC strategy is presented. The switching logic is designed for different MPC controllers. In the scheme, a Kalman estimator is designed to estimate the vehicle load torque. With the estimated load torque, the prediction can be much more accurate. In Section 4, the proposed strategy is tested and discussed. Finally, the concluding remarks are given in Section 5.


## 2. System Modeling

Figure1 shows a rear-driven heavy-duty truck. The driveline mainly consists of the dry clutch, gearbox, propeller shaft, differential, half driven shafts, and wheels. The clutch electropneumatic circuit is the actuator to drive the clutch close or open. The engine is the power source which is also an actuator working together with clutch actuators.Figure 1
Overall structure of the AMT truck.Figure2 shows a simplified driveline model, whereas a three-inertia model is used in this study. One inertia indicates the engine flywheel Je with viscous friction de. Another inertia indicates the equivalent mass moment from the clutch disc to the differential gear, which is summarized to J1 with viscous friction d1 between the clutch disc and gearbox, and the last inertia J2 indicates the equivalent mass moment of remaining driveline and the vehicle with viscous friction d2.Figure 2
Simplified model of the driveline.The clutch friction torque is represented byTc and the engine torque is represented by Te. The gearbox ratio together with the final gear is described by the ratio r. It is assumed that the flexibility of drivetrain is equivalent to a spring and damper represented by ks and ds. Wheel slip is neglected for simplicity. The driving resistance torque is described by Tl, which consists of the rolling resistance torque, grade resistance torque, and the aerodynamics drag torque. Usually, the resistance torques is unknown.The equivalent inertia ofJ1 can be calculated using (1) as follows:(1)J1=Jc+Jg1+Jg2rg2+Jdrgrd2,where Jc is the inertia of clutch; Jg1 is the inertia of the input gear of the transmission; Jg2 is the inertia of the output gear of the transmission; Jd is inertia of differential input gear; rg is the gear transmission ratio; and rd is the differential ratio.The equivalent inertia ofJ2 can be calculated using (2) as follows:(2)J2=Jw+mRw2,where Jw is the wheel inertia, m is the mass of vehicle body, and Rw is the tire radius.The equivalent stiffness ofks can be calculated using (3) as follows. (The damping coefficient of ds is calculated by the same equation with replacing the stiffness coefficient by the damping coefficient.) The model acts as if there are several damped springs in series and parallel:(3)1ks=1kcr2+1kprd2+1kh1+kh2,where kc is the stiffness coefficient of clutch torsion damper, kp is the stiffness coefficient of propeller shaft, kh1 is the stiffness coefficient of left half shaft, and kh2 is the stiffness coefficient of right half shaftThe dynamic of the drive train is decided by the clutch status. When the clutch is slipping, it is shown in(4)Jeω˙e=Te-Tc-deω˙eJ1ω˙c=Tc-d1ω˙c-ksθs+dsθ˙srJ2ω˙w=-Tl-d2ω˙w+ksθs+dsθ˙θ˙s=θ˙cr-θ˙w=ωcr-ωw,where θs is the drive shaft torsion angle; ωe is the engine speed; ωc is the clutch speed; ωw is the wheel speed.At this time, the clutch friction torque is generated by Coulomb friction and decided mainly by the actuator forceFn, as shown in(5)Tc=FnμRasign⁡ωe-ωc,where Fn is the press force actuating on the clutch plate, μ is the dynamic friction coefficient of the clutch surface material, and Ra is the effective radius of the clutch plate.When the clutch is engaged, the engine is rigidly coupled to the driveline. At this time, the clutch slipping speedωsl=ωe-ωc=0. Consequently, the two equations of motion of the engine and the clutch are merged into a single equation:(6)Je+J1ω˙e/c=Te-de+d1ω˙e/c-ksθs+dsθ˙sr.Moreover, during engaged phase, the torque through the clutch cannot be altered by the actuator force anymore. Instead of a controlled input, it becomes a constrained variable.(7)Tc≤Tc,max=Fnμ0Ra,where Tc,max is the maximum friction torque of the clutch disc and μ0 is the static friction coefficient.So the switch logic model is shown in Figure3. We define two statuses of the clutch as Slipping and Engaged.Figure 3
Switching logic of the clutch statuses.Considering the effect of the engine and clutch actuator dynamics, the engine torque and the clutch friction torque are modeled as first-order dynamics with time delaysτe and τc, respectively. The equations are as follows:(8)τeT˙e+Te=Te,sp+Te,reqτcT˙c+Tc=Tc,sp,where Te,sp and Tc,sp are the set points of the clutch controller for the engine controller and clutch actuator controller. And Te,req is the driver demand torque. The torque-based controls are now commonly used in modern engine controllers. The torque requirement from the clutch controller will be added to the driver demand torque in the engine controller. Therefore, the real engine torque is the sum of these requirements.The engine torque delay is related to the engine speed and engine cylinder number, which is often presented by (9), where the cylinder number of the engine is nclinder and the engine speed is ne in revolutions per second.(9)τe=2nclinderne.And for the clutch actuator time delay parameter, we can use simulation result of the detailed clutch actuator model to identify it. Literature [9] has shown that the first-order linear function with the identified parameters has a similar simulation result to the detailed model.In the continuous state-space representation, the driveline model can be written as follows (10). In the equation, the controller outputs and disturbance torque are distinguished.(10)x˙=Acx+Bcuu+BcrTe,req+BclTl=Ac1x+Bcu1u+BcrTe,req+BclTlclutchstatus=SlippingAc2x+Bcu2u+BcrTe,req+BclTlclutchstatus=Engagedy=Ccx,where the state variables, input variables, and output variables are defined in as follows:(11)x=ωeθsωcωwTeTcTu=Te,spTc,spy=ωeωslωsTeTc.State-space equation parameters are defined as follows:(12a)Ac1=-deJe0001Je-1Je001r-1000-ks/rJ1-d1+ds/r2J1ds/rJ101J10ksJ2ds/rJ2-d1+dsJ20000001τe0000001τcAc2=0-ks/rJe+J1-de+d1+ds/r2Je+J1ds/rJe+J11Je+J10001r-1000-ks/rJe+J1-de+d1+ds/r2Je+J1ds/rJe+J11Je+J100ksJ2ds/rJ2-d1+dsJ20000001τe0000000(12b)Bcu1=00001τe0000001τcTBcu2=00001τe0000000TBcr=00001τe0TBcl=000-1J200T


## 3. Controller Design

### 3.1. Control Objectives

For the AMT clutch controller design, different phases need different control objectives. Usually, when vehicle is starting up, the clutch status is Slipping and Engaged in turn. Because the slipping time and the driving comfort are conflicted, a driver intention will be used to change the slipping time dynamically. When the driver wants a faster dynamic process, the slipping time will be decreased. But when the driver does not need so fast dynamic process, then the slipping time could be increased, so as to increase the driving comfort.Because the torsion vibration influences the driving comfort badly, the objective of torsion speedθ˙s should be controlled in every phase. Good drivability also means that vehicle responds to the driver’s demand exactly with no delay. In this sense, the controller output should not have a big effect on the driver’s intended torque. Therefore, the controller is expected to minimize the engine torque requirement from the clutch controller in clutch Engaged phase.Moreover, the controller should satisfy the constraints such as that the engine torque and clutch torque set points should be limited in the permitted minimum and maximum values. And in order to keep the engine in the so-called no-kill condition [9], the engine speed should be bigger higher than the idle speed.For every phase, the control objective is defined as follows.


#### 3.1.1. Slipping Phase

In Slipping phase, the conversion condition isωsl=0. We design a reference trajectory for clutch slipping speed as follows (13). Typical trajectory is showed in Figure 4.(13)ωslref=2ωsl,0t-t0,slippingtf,slipping3-3ωsl,0t-t0,slippingtf,slipping2+ωsl,0.Figure 4
Reference trajectory of clutch slipping speed.It obeys the constraint as follows:(14)ωslreft0,slipping=ωsl,0,ωslreftf,slipping=0,ω˙slreft0,slipping=ω˙slreftf,slipping=0,where ωslref is the reference slipping speed, ωsl,0 is the initial slipping speed, t0,slipping is the time of slipping start, and tf,slipping is the duration of slipping; it is a variable which is a function of driver acceleration pedal α, shown in (15), where tf,slippingmax and tf,slippingmin are the maximum and minimum permitted opening time, respectively:(15)tf,slipping=tf,slippingmax-tf,slippingmax-tf,slippingminα.Meanwhile the torsion vibration should be kept small. And the engine torque should be kept close to the driver torque intended. Therefore, the control objective can be written as (16). Jslipping is the cost function of this phase.(16)min⁡Jslipping=min⁡λ1,slippingy2-ωslref2+λ2,slippingy32+λ3,slippingu12,where λ1,slipping, λ2,slipping, and λ3,slipping are weighting factors of slipping time, torsion vibration, and driver torque intended.


#### 3.1.2. Engaged Phase

In Engaged phase, the engine is synchronized with the transmission. The clutch torque does not influence the dynamics of the driveline. But the clutch maximum static friction torque should be kept not smaller than the engine torque, so that slipping is avoided. Meanwhile, the torsion vibration should be reduced, and the engine torque should be close to the driver torque intended. So the control objective can be written as (17). Jengaged is the cost function of this phase.(17)min⁡Jengaged=min⁡λ2,engagedy32+λ3,engagedu12(18)Tc,sp>Te,sp,where λ2,engaged and λ3,engaged are weighting factors of torsion vibration and driver torque intended.


### 3.2. Control Scheme

Multiple-MPC (mMPC) strategy is used for clutch control, due to the ability to handle multivariable system, to take time-domain constraints into account explicitly and deal with multiple objectives in a somehow optimal sense. What is more, there are multiphases with different control targets, so we design different controllers for different phases and switch the controllers by the clutch status.The principle of MPC is to calculate a sequence of control actions over a finite receding horizon by optimizing the certain optimization function, while only implementing the first element on the control plant each time. In this application, the clutch control scheme is shown in Figure5.Figure 5
mMPC based clutch controller scheme.An estimator uses the measured engine speed, clutch speed, and wheel speed to estimate the unmeasured disturbanceTl and other state variables of the state-space equation in (10). Therefore, all of the state variables and disturbance input variables are measured. We can predict the future output variables of the equation. And an optimization function which is defined in Section 3.1 could be solved to calculate the target engine torque and clutch friction torque.


### 3.3. Load and State Variables Estimation

The vehicle is subjected to various load disturbances. Since the load mainly depends on unmeasured entities, such as slope, it is difficult to include it in the model. The literature [13] presented a method to estimate the load torque Tl. The method is based on the assumption that the load disturbance is constant or slowly varying, so it can be described or modeled as follows:(19)T˙l=0.TreatingTl as a state variable, the system can be described as follows:(20)x˙T˙l=AcBcl00xTl+BcuBcr00uTe,req=Alxl+Bluly=Cc0xTl=Clxl.The state estimation is given by(21)x˙^l=Alx^l+Blul+Kly-Clx^l.ThenKl is a Kalman estimator coefficient and can be calculated using the algebraic Riccati equation to minimize the estimation error.By now, the disturbance torquesTe,req and Tl can be treated the same. So we define two new symbols as follows:(22)Td=Te,reqTlBcd=BcrBclT.So the state-space equation (10) can be rewritten as follows:(23)x˙=Acx+Bcuu+BcdTd.


### 3.4. Multiple-Model Predictive Controller

Because the system state-space equation considers the engine time delay which is not a constant, the coefficients of the state equation are changed step by step. If we use method of nonlinear MPC, it will be very complex and not suitable for real-time application. Here, we introduce a new method for solving this problem. Because the delay time of the engine torque generation does not change much in a prediction horizon, we assume that the delay time is a constant in a prediction horizon. Then we can deal with the equation as a linear system in every prediction process and keep the delay time updating step by step. So in every step, we need to firstly calculate the delay timeτe(k) using the measured state ωe(k) and select the right model based on the clutch status in (10). Then we update the coefficients of the state-space equation Ac(k), Bcu(k), and Bcr(k), which have relationship with τe(k). We convert the updated continuous state-space equation to discrete state-space equation using forward Euler method. Using the discrete state-space equation, we estimate the load torque T^l(k) and state variables which are not measured directly based on Kalman filter strategy. After that, we can predict the output vector of the system Y(k+1∣k) in predictive horizon. By solving the right optimization function based on the clutch status defined in Section 3.1 with constraints, the control vector in control horizon U(k) can be calculated. The first element of U(k) is selected to be the controller output. The whole process of the mMPC is described in Figure 6.Figure 6
mMPC clutch controller process.According to the main principle of predictive control [14], with the measured state vector x^(k) as the initial condition at time k, the future torsion speed is predicted on the discrete system model. In this paper, the prediction horizon is illustrated with p. The control horizon is illustrated with m, satisfying m≤p. In order to reduce the system prediction equations, the assumptions are set as follows.Out of control horizon, the control variable keeps unchanged.Disturbance torque keeps unchanged in prediction horizon.We define that the predicted output vector in horizonp is Y(k+1∣k) and the predicted input vector in horizon m is U(k) as shown in(24)Yk+1∣k=Yk+1∣kYk+2∣k⋮yk+p∣kp×1,Uk=ukuk+1⋮uk+m-1m×1,where k+1∣k represents the prediction of time k+1 at time k.According to the basis of MPC and exploiting above contents, we can infer the sequences of outputs to be predicted and present them in the form of(25)Yk+1∣k=Sxx^k+SuUk+SdT^dk,where Sx, Sd, and Su are calculated by (26) and T^d(k) is calculated using the estimation variable T^l(k), which is obtained from (27).(26)Sx=CACA2⋮CApp×1Sd=CBd∑i=12CAdi-1Bd⋮∑i=1pCAdi-1Bdp×1Su=CBu0⋯0CABuCBu⋯0⋮⋮⋮⋮CAp-1BuCAp-2Bu⋯CAp-mBup×m(27)T^dk=TrkT^lkT.Because the main control requirement is to minimize the cost function, the optimization function can be defined as(28)minUk⁡JUk,Yk+1∣k=WyYk+1∣k-Rek+1∣k2+WuUk2,where Wy and Wu are the weighting factor matrixes of the two objectives, respectively. Re⁡(k+1∣k) is the reference trajectory matrix. In different phases, these two weighting factors are composed with λ defined in Section 3.1. For the reference trajectory, because it includes time information, in every prediction step, the value can be calculated based on the time step.The optimization problem (28) can be formulated as a quadratic programming (QP) problem in(29)minUk⁡UkTHUk-Gk+1∣kTUk,where(30a)H=SuTWyTWySu+WuTWu(30b)Gk+1∣k=2SuTWyTWyEpk+1∣k(30c)Epk+1∣k=Rek+1-Sxxk-SdT^dk.We assume that the constraints are also kept unchanged in prediction horizon. Solving the optimization function with the constraints can get the control sequenceU(k). Applying the first element u(k) to the plant and this process is repeated at every time step.


## 3.1. Control Objectives

For the AMT clutch controller design, different phases need different control objectives. Usually, when vehicle is starting up, the clutch status is Slipping and Engaged in turn. Because the slipping time and the driving comfort are conflicted, a driver intention will be used to change the slipping time dynamically. When the driver wants a faster dynamic process, the slipping time will be decreased. But when the driver does not need so fast dynamic process, then the slipping time could be increased, so as to increase the driving comfort.Because the torsion vibration influences the driving comfort badly, the objective of torsion speedθ˙s should be controlled in every phase. Good drivability also means that vehicle responds to the driver’s demand exactly with no delay. In this sense, the controller output should not have a big effect on the driver’s intended torque. Therefore, the controller is expected to minimize the engine torque requirement from the clutch controller in clutch Engaged phase.Moreover, the controller should satisfy the constraints such as that the engine torque and clutch torque set points should be limited in the permitted minimum and maximum values. And in order to keep the engine in the so-called no-kill condition [9], the engine speed should be bigger higher than the idle speed.For every phase, the control objective is defined as follows.


### 3.1.1. Slipping Phase

In Slipping phase, the conversion condition isωsl=0. We design a reference trajectory for clutch slipping speed as follows (13). Typical trajectory is showed in Figure 4.(13)ωslref=2ωsl,0t-t0,slippingtf,slipping3-3ωsl,0t-t0,slippingtf,slipping2+ωsl,0.Figure 4
Reference trajectory of clutch slipping speed.It obeys the constraint as follows:(14)ωslreft0,slipping=ωsl,0,ωslreftf,slipping=0,ω˙slreft0,slipping=ω˙slreftf,slipping=0,where ωslref is the reference slipping speed, ωsl,0 is the initial slipping speed, t0,slipping is the time of slipping start, and tf,slipping is the duration of slipping; it is a variable which is a function of driver acceleration pedal α, shown in (15), where tf,slippingmax and tf,slippingmin are the maximum and minimum permitted opening time, respectively:(15)tf,slipping=tf,slippingmax-tf,slippingmax-tf,slippingminα.Meanwhile the torsion vibration should be kept small. And the engine torque should be kept close to the driver torque intended. Therefore, the control objective can be written as (16). Jslipping is the cost function of this phase.(16)min⁡Jslipping=min⁡λ1,slippingy2-ωslref2+λ2,slippingy32+λ3,slippingu12,where λ1,slipping, λ2,slipping, and λ3,slipping are weighting factors of slipping time, torsion vibration, and driver torque intended.


### 3.1.2. Engaged Phase

In Engaged phase, the engine is synchronized with the transmission. The clutch torque does not influence the dynamics of the driveline. But the clutch maximum static friction torque should be kept not smaller than the engine torque, so that slipping is avoided. Meanwhile, the torsion vibration should be reduced, and the engine torque should be close to the driver torque intended. So the control objective can be written as (17). Jengaged is the cost function of this phase.(17)min⁡Jengaged=min⁡λ2,engagedy32+λ3,engagedu12(18)Tc,sp>Te,sp,where λ2,engaged and λ3,engaged are weighting factors of torsion vibration and driver torque intended.


## 3.1.1. Slipping Phase

In Slipping phase, the conversion condition isωsl=0. We design a reference trajectory for clutch slipping speed as follows (13). Typical trajectory is showed in Figure 4.(13)ωslref=2ωsl,0t-t0,slippingtf,slipping3-3ωsl,0t-t0,slippingtf,slipping2+ωsl,0.Figure 4
Reference trajectory of clutch slipping speed.It obeys the constraint as follows:(14)ωslreft0,slipping=ωsl,0,ωslreftf,slipping=0,ω˙slreft0,slipping=ω˙slreftf,slipping=0,where ωslref is the reference slipping speed, ωsl,0 is the initial slipping speed, t0,slipping is the time of slipping start, and tf,slipping is the duration of slipping; it is a variable which is a function of driver acceleration pedal α, shown in (15), where tf,slippingmax and tf,slippingmin are the maximum and minimum permitted opening time, respectively:(15)tf,slipping=tf,slippingmax-tf,slippingmax-tf,slippingminα.Meanwhile the torsion vibration should be kept small. And the engine torque should be kept close to the driver torque intended. Therefore, the control objective can be written as (16). Jslipping is the cost function of this phase.(16)min⁡Jslipping=min⁡λ1,slippingy2-ωslref2+λ2,slippingy32+λ3,slippingu12,where λ1,slipping, λ2,slipping, and λ3,slipping are weighting factors of slipping time, torsion vibration, and driver torque intended.


## 3.1.2. Engaged Phase

In Engaged phase, the engine is synchronized with the transmission. The clutch torque does not influence the dynamics of the driveline. But the clutch maximum static friction torque should be kept not smaller than the engine torque, so that slipping is avoided. Meanwhile, the torsion vibration should be reduced, and the engine torque should be close to the driver torque intended. So the control objective can be written as (17). Jengaged is the cost function of this phase.(17)min⁡Jengaged=min⁡λ2,engagedy32+λ3,engagedu12(18)Tc,sp>Te,sp,where λ2,engaged and λ3,engaged are weighting factors of torsion vibration and driver torque intended.


## 3.2. Control Scheme

Multiple-MPC (mMPC) strategy is used for clutch control, due to the ability to handle multivariable system, to take time-domain constraints into account explicitly and deal with multiple objectives in a somehow optimal sense. What is more, there are multiphases with different control targets, so we design different controllers for different phases and switch the controllers by the clutch status.The principle of MPC is to calculate a sequence of control actions over a finite receding horizon by optimizing the certain optimization function, while only implementing the first element on the control plant each time. In this application, the clutch control scheme is shown in Figure5.Figure 5
mMPC based clutch controller scheme.An estimator uses the measured engine speed, clutch speed, and wheel speed to estimate the unmeasured disturbanceTl and other state variables of the state-space equation in (10). Therefore, all of the state variables and disturbance input variables are measured. We can predict the future output variables of the equation. And an optimization function which is defined in Section 3.1 could be solved to calculate the target engine torque and clutch friction torque.


## 3.3. Load and State Variables Estimation

The vehicle is subjected to various load disturbances. Since the load mainly depends on unmeasured entities, such as slope, it is difficult to include it in the model. The literature [13] presented a method to estimate the load torque Tl. The method is based on the assumption that the load disturbance is constant or slowly varying, so it can be described or modeled as follows:(19)T˙l=0.TreatingTl as a state variable, the system can be described as follows:(20)x˙T˙l=AcBcl00xTl+BcuBcr00uTe,req=Alxl+Bluly=Cc0xTl=Clxl.The state estimation is given by(21)x˙^l=Alx^l+Blul+Kly-Clx^l.ThenKl is a Kalman estimator coefficient and can be calculated using the algebraic Riccati equation to minimize the estimation error.By now, the disturbance torquesTe,req and Tl can be treated the same. So we define two new symbols as follows:(22)Td=Te,reqTlBcd=BcrBclT.So the state-space equation (10) can be rewritten as follows:(23)x˙=Acx+Bcuu+BcdTd.


## 3.4. Multiple-Model Predictive Controller

Because the system state-space equation considers the engine time delay which is not a constant, the coefficients of the state equation are changed step by step. If we use method of nonlinear MPC, it will be very complex and not suitable for real-time application. Here, we introduce a new method for solving this problem. Because the delay time of the engine torque generation does not change much in a prediction horizon, we assume that the delay time is a constant in a prediction horizon. Then we can deal with the equation as a linear system in every prediction process and keep the delay time updating step by step. So in every step, we need to firstly calculate the delay timeτe(k) using the measured state ωe(k) and select the right model based on the clutch status in (10). Then we update the coefficients of the state-space equation Ac(k), Bcu(k), and Bcr(k), which have relationship with τe(k). We convert the updated continuous state-space equation to discrete state-space equation using forward Euler method. Using the discrete state-space equation, we estimate the load torque T^l(k) and state variables which are not measured directly based on Kalman filter strategy. After that, we can predict the output vector of the system Y(k+1∣k) in predictive horizon. By solving the right optimization function based on the clutch status defined in Section 3.1 with constraints, the control vector in control horizon U(k) can be calculated. The first element of U(k) is selected to be the controller output. The whole process of the mMPC is described in Figure 6.Figure 6
mMPC clutch controller process.According to the main principle of predictive control [14], with the measured state vector x^(k) as the initial condition at time k, the future torsion speed is predicted on the discrete system model. In this paper, the prediction horizon is illustrated with p. The control horizon is illustrated with m, satisfying m≤p. In order to reduce the system prediction equations, the assumptions are set as follows.Out of control horizon, the control variable keeps unchanged.Disturbance torque keeps unchanged in prediction horizon.We define that the predicted output vector in horizonp is Y(k+1∣k) and the predicted input vector in horizon m is U(k) as shown in(24)Yk+1∣k=Yk+1∣kYk+2∣k⋮yk+p∣kp×1,Uk=ukuk+1⋮uk+m-1m×1,where k+1∣k represents the prediction of time k+1 at time k.According to the basis of MPC and exploiting above contents, we can infer the sequences of outputs to be predicted and present them in the form of(25)Yk+1∣k=Sxx^k+SuUk+SdT^dk,where Sx, Sd, and Su are calculated by (26) and T^d(k) is calculated using the estimation variable T^l(k), which is obtained from (27).(26)Sx=CACA2⋮CApp×1Sd=CBd∑i=12CAdi-1Bd⋮∑i=1pCAdi-1Bdp×1Su=CBu0⋯0CABuCBu⋯0⋮⋮⋮⋮CAp-1BuCAp-2Bu⋯CAp-mBup×m(27)T^dk=TrkT^lkT.Because the main control requirement is to minimize the cost function, the optimization function can be defined as(28)minUk⁡JUk,Yk+1∣k=WyYk+1∣k-Rek+1∣k2+WuUk2,where Wy and Wu are the weighting factor matrixes of the two objectives, respectively. Re⁡(k+1∣k) is the reference trajectory matrix. In different phases, these two weighting factors are composed with λ defined in Section 3.1. For the reference trajectory, because it includes time information, in every prediction step, the value can be calculated based on the time step.The optimization problem (28) can be formulated as a quadratic programming (QP) problem in(29)minUk⁡UkTHUk-Gk+1∣kTUk,where(30a)H=SuTWyTWySu+WuTWu(30b)Gk+1∣k=2SuTWyTWyEpk+1∣k(30c)Epk+1∣k=Rek+1-Sxxk-SdT^dk.We assume that the constraints are also kept unchanged in prediction horizon. Solving the optimization function with the constraints can get the control sequenceU(k). Applying the first element u(k) to the plant and this process is repeated at every time step.


## 4. Evaluation the Controllers by Simulation

### 4.1. Simulation Model

In this section, some simulation results are given and analyzed to evaluate the controllers. In order to simulate, a simulation model shown in Figure7 is built in Matlab/Simulink/Simscape. It includes the longitudinal vehicle dynamics, lumped mass driveline, mean value engine model, pneumatic clutch actuator, and the engine controller and clutch actuator controller model.. The constructed simulation model can capture the important transient dynamics of the driveline, such as the delay of the engine torque generation, the delay of the clutch actuator, clutch damper vibration, flexible drive shaft, and half drive shaft oscillation. As well as tire slip, even road slope and road surface which may have an effect on the clutch controller are all included in the model. Some main vehicle parameters are listed in Table 1.Table 1
Normal values of parameters of simulation model.
Symbol
Parameters
Value & unit
J
e
Engine flywheel inertia
2.1 kg/m2
J
c
Clutch disk inertia
0.1 kg/m2
J
g
1
Equivalent inertia of gearbox primary shaft
0.1 kg/m2
J
g
2
Equivalent inertia of gearbox secondary shaft
0.1 kg/m2
J
d
Equivalent inertia of propeller shaft and differential gear inertia
0.1 kg/m2
J
2
Equivalent inertia of half drive shaft and wheels and vehicle
4020 kg/m2
d
e
Damping of engine shaft to bearings
0.1 Nm/rad/s
d
c
Damping of clutch shaft to bearings
0.1 Nm/rad/s
d
2
Damping of wheels to bearings
2 Nm/rad/s
d
s
Damping of the drive shaft
3 kNm/rad/s
k
s
Stiffness of the drive shaft
175 kNm/rad
r
Transmission ratio of drive train
35.04
n
Engine cylinder number
6
τ
c
Clutch actuator delay time
0.3 s
m
Vehicle mass
16 ton
R
w
Radius of tire
0.5 mFigure 7
Simulation model of a heavy-duty truck.For the engine model, the input is fuel mass and, for the clutch actuator, the input is solenoid On/Off valves; an engine controller and clutch actuator controller are also modeled with simplified control strategy. And we assume that these two controllers have high control accuracy. In order to simulate the gearshift, a simple gear controller is modeled. It just includes gearshift timing. And the clutch control strategy described in this paper is also built as a model in Simulink. In this simulation, the sampling period for MPC control updating is set to 10 ms, and the simulation updating period is 1 ms.


### 4.2. Simulation Results

In order to fully evaluate the controllers, there are several maneuvers defined as follows:(1)
Open loop control for the start-up with low torque demand (M1)(2)
Open loop control for the start-up with high torque demand (M2)(3)
Close loop control for the start-up with low torque demand (M3)(4)
Close loop control for the start-up with high torque demand (M4)Figures8 and 9 are the simulation results of maneuvers M1 and M2. M1 and M2 are maneuvers without clutch optimization controller. M1 starts with a low driver demand engine torque which is 100 Nm. And clutch friction torque is also set to 100 Nm. M2 starts with a high driver demand engine torque which is 800 Nm. And clutch friction torque is also set to 800 Nm. From the simulation results, we can see that, for both M1 and M2, there are torsion vibrations in both Slipping and Engaged phase. And the vibration frequency in Slipping phase is higher than that in Engaged phase. Higher torque can evoke serious oscillation but can shorten the slipping time. From Figures 8(c) and 9(c), we can see that the acceleration curves also show many shakes and jerks in the whole process. The shake and jerk evoke a bad drivability.Figure 8
M1 maneuver simulation results.
(a)
(b)
(c)Figure 9
M2 maneuver simulation results.
(a)
(b)
(c)Figure10 is the simulation results of maneuver M3. M3 is the start-up maneuver with optimization controller with low driver torque demand which is 100 Nm. It is a low driver pedal position input. Therefore, the controller generates a reference slipping speed trajectory as shown in Figure 10(d) lasted about 5 s. From the simulation results in Figures 10(a)–10(d), we can see that the torsion vibration is reduced and the engine speed tends to increase at a nearly fixed rate. The acceleration of the whole process does not have many shakes and jerks, and the slipping speed has a good tracking to the reference trajectory. From Figures 10(e)-10(f), we can see that the controller outputs are in the area of constraints. We set the constraint of engine torque requirement from clutch control as [−100 Nm, 200 Nm]. The minimum value is based on that the engine real torque cannot be lower than zero. The maximum value is based on the low driver pedal maximum engine torque limit. We set the constraint of the clutch friction torque as [0, 300 Nm], with considering that the clutch torque should not be larger or lower than engine torque.Figure 10
M3 maneuver simulation results.
(a)
(b)
(c)
(d)
(e)
(f)Figure11 shows the simulation results of maneuver M4. M4 is the start-up maneuver with optimization controller with high driver torque demand which is 800 Nm. It is a high driver pedal position input. So the controller generates a reference slipping speed trajectory as shown in Figure 11(d) lasted about 1.5 s. From the simulation results in Figures 11(a)–11(d), we can see that the torsion vibration is reduced and the engine speed tends to increase at nearly fixed rate. The acceleration of the whole process does not have many shakes and jerks, and the slipping speed has a good tracking to the reference trajectory. From Figures 11(e)-11(f), we can see the controller outputs are in the area of constraints. We set the constraint of engine torque requirement from clutch control to [−800 Nm, 300 Nm]. The minimum value is based on that the engine real torque cannot be lower than zero. The maximum value is based on the high driver pedal maximum engine torque limit. We set the constraint of the clutch friction torque to [0,  1100 Nm], with considering that the clutch torque should not be larger or lower than the engine torque.Figure 11
M4 maneuver simulation results.
(a)
(b)
(c)
(d)
(e)
(f)


## 4.1. Simulation Model

In this section, some simulation results are given and analyzed to evaluate the controllers. In order to simulate, a simulation model shown in Figure7 is built in Matlab/Simulink/Simscape. It includes the longitudinal vehicle dynamics, lumped mass driveline, mean value engine model, pneumatic clutch actuator, and the engine controller and clutch actuator controller model.. The constructed simulation model can capture the important transient dynamics of the driveline, such as the delay of the engine torque generation, the delay of the clutch actuator, clutch damper vibration, flexible drive shaft, and half drive shaft oscillation. As well as tire slip, even road slope and road surface which may have an effect on the clutch controller are all included in the model. Some main vehicle parameters are listed in Table 1.Table 1
Normal values of parameters of simulation model.
Symbol
Parameters
Value & unit
J
e
Engine flywheel inertia
2.1 kg/m2
J
c
Clutch disk inertia
0.1 kg/m2
J
g
1
Equivalent inertia of gearbox primary shaft
0.1 kg/m2
J
g
2
Equivalent inertia of gearbox secondary shaft
0.1 kg/m2
J
d
Equivalent inertia of propeller shaft and differential gear inertia
0.1 kg/m2
J
2
Equivalent inertia of half drive shaft and wheels and vehicle
4020 kg/m2
d
e
Damping of engine shaft to bearings
0.1 Nm/rad/s
d
c
Damping of clutch shaft to bearings
0.1 Nm/rad/s
d
2
Damping of wheels to bearings
2 Nm/rad/s
d
s
Damping of the drive shaft
3 kNm/rad/s
k
s
Stiffness of the drive shaft
175 kNm/rad
r
Transmission ratio of drive train
35.04
n
Engine cylinder number
6
τ
c
Clutch actuator delay time
0.3 s
m
Vehicle mass
16 ton
R
w
Radius of tire
0.5 mFigure 7
Simulation model of a heavy-duty truck.For the engine model, the input is fuel mass and, for the clutch actuator, the input is solenoid On/Off valves; an engine controller and clutch actuator controller are also modeled with simplified control strategy. And we assume that these two controllers have high control accuracy. In order to simulate the gearshift, a simple gear controller is modeled. It just includes gearshift timing. And the clutch control strategy described in this paper is also built as a model in Simulink. In this simulation, the sampling period for MPC control updating is set to 10 ms, and the simulation updating period is 1 ms.


## 4.2. Simulation Results

In order to fully evaluate the controllers, there are several maneuvers defined as follows:(1)
Open loop control for the start-up with low torque demand (M1)(2)
Open loop control for the start-up with high torque demand (M2)(3)
Close loop control for the start-up with low torque demand (M3)(4)
Close loop control for the start-up with high torque demand (M4)Figures8 and 9 are the simulation results of maneuvers M1 and M2. M1 and M2 are maneuvers without clutch optimization controller. M1 starts with a low driver demand engine torque which is 100 Nm. And clutch friction torque is also set to 100 Nm. M2 starts with a high driver demand engine torque which is 800 Nm. And clutch friction torque is also set to 800 Nm. From the simulation results, we can see that, for both M1 and M2, there are torsion vibrations in both Slipping and Engaged phase. And the vibration frequency in Slipping phase is higher than that in Engaged phase. Higher torque can evoke serious oscillation but can shorten the slipping time. From Figures 8(c) and 9(c), we can see that the acceleration curves also show many shakes and jerks in the whole process. The shake and jerk evoke a bad drivability.Figure 8
M1 maneuver simulation results.
(a)
(b)
(c)Figure 9
M2 maneuver simulation results.
(a)
(b)
(c)Figure10 is the simulation results of maneuver M3. M3 is the start-up maneuver with optimization controller with low driver torque demand which is 100 Nm. It is a low driver pedal position input. Therefore, the controller generates a reference slipping speed trajectory as shown in Figure 10(d) lasted about 5 s. From the simulation results in Figures 10(a)–10(d), we can see that the torsion vibration is reduced and the engine speed tends to increase at a nearly fixed rate. The acceleration of the whole process does not have many shakes and jerks, and the slipping speed has a good tracking to the reference trajectory. From Figures 10(e)-10(f), we can see that the controller outputs are in the area of constraints. We set the constraint of engine torque requirement from clutch control as [−100 Nm, 200 Nm]. The minimum value is based on that the engine real torque cannot be lower than zero. The maximum value is based on the low driver pedal maximum engine torque limit. We set the constraint of the clutch friction torque as [0, 300 Nm], with considering that the clutch torque should not be larger or lower than engine torque.Figure 10
M3 maneuver simulation results.
(a)
(b)
(c)
(d)
(e)
(f)Figure11 shows the simulation results of maneuver M4. M4 is the start-up maneuver with optimization controller with high driver torque demand which is 800 Nm. It is a high driver pedal position input. So the controller generates a reference slipping speed trajectory as shown in Figure 11(d) lasted about 1.5 s. From the simulation results in Figures 11(a)–11(d), we can see that the torsion vibration is reduced and the engine speed tends to increase at nearly fixed rate. The acceleration of the whole process does not have many shakes and jerks, and the slipping speed has a good tracking to the reference trajectory. From Figures 11(e)-11(f), we can see the controller outputs are in the area of constraints. We set the constraint of engine torque requirement from clutch control to [−800 Nm, 300 Nm]. The minimum value is based on that the engine real torque cannot be lower than zero. The maximum value is based on the high driver pedal maximum engine torque limit. We set the constraint of the clutch friction torque to [0,  1100 Nm], with considering that the clutch torque should not be larger or lower than the engine torque.Figure 11
M4 maneuver simulation results.
(a)
(b)
(c)
(d)
(e)
(f)


## 5. Conclusions

A multiple-MPC for dry clutch control for AMT truck has been proposed. Two driveline models are built based on the clutch position. Two different control objectives are designed for different phases including Slipping and Engaged phases. Therefore, two MPC controllers with two prediction models are designed for the AMT clutch control for the start-up process. In these controllers, the clutch slipping time is decided by the driver’s intention. The torsion vibration is controlled not only in the Engaged phase but also in the Slipping phase. Simulation results have shown the mMPC controllers have good performance in each phase during start-up process in AMT truck. The drivability of the vehicle is enhanced with the optimization controller.The proposed scheme is only for the start-up process control so far, it will be enhanced in the future for the gearshift process control to form a complete high level controller for the AMT clutch control.



---

*Source: 5783527-2017-05-11.xml*