Drivability-Related Discrete-Time Model Predictive Control of Mode Transition in Pre-Transmission Parallel Hybrid Powertrains

During the mode transition from the pure electric propulsion mode to the hybrid propulsion mode, clutch-based pre-transmission parallel hybrid electric vehicles are subject to drivability issues. These issues originate from the fact that in the clutch-based pre-transmission parallel hybrid powertrain (CPPHP) configuration, the clutch connects the engine and the motor. Without a carefully designed mode transition control that coordinates the engine torque, clutch torque and motor torque, torque sluggishness and surges occur during the mode transition, and residual torque oscillation occurs after the mode transition. In this paper, a discrete-time model predictive control (DMPC)-based controller is proposed to address these drivability-related issues. Modeling improvements and novel drivability-related indices and constraints are all taken into consideration in the design of the discrete-time model predictive controller. Furthermore, by using discrete-time Laguerre functions and introducing the equilibrium state and the ranking of constraints, an explicit solution of the discrete-time model predictive controller is obtained. The calculation results demonstrate that the proposed controller can ensure a smooth and rapidly decaying torque difference during the mode transition, alleviating the residual torque oscillation after the mode transition and guaranteeing that the mode transition is completed within an acceptable duration.


Introduction
In comparison with conventional powertrains driven only by internal combustion engines, hybrid powertrains augment fuel efficiency and significantly reduce emissions [1].Three typical hybrid powertrain architectures are widely used in hybrid electric vehicles (HEVs) on the market: planetary gear set based series-parallel hybrids [2], clutch-based series-parallel hybrids [3] and clutch-based parallel hybrids [4].These configurations are schematically depicted in Figure 1.
The clutch-based pre-transmission parallel hybrid is the most cost-effective type of hybrid powertrain configuration among the aforementioned three [5].Numerous studies on clutch-based pre-transmission parallel hybrid powertrains (CPPHPs) have focused on energy management [6] and emissions control [7].Nevertheless, during the mode transition from the pure electric propulsion mode to the hybrid propulsion mode, drivability-related issues emerge in CPPHPs due to the clutch engagement.Within the finite duration of this process, multiple objectives need to be accomplished simultaneously: adjusting the engine and motor torques to their respective target values determined by the energy management system, sustaining a low level of torque difference between the supplied by the hybrid power source and the demanded by the driver during the mode transition, achieving a smooth clutch engagement at the ending of the mode transition and ensuring a low level of residual torque oscillation after the mode transition.These are the most significant challenges when designing a mode transition controller.
Energies 2016, 9, 740 2 of 30 by the energy management system, sustaining a low level of torque difference between the supplied by the hybrid power source and the demanded by the driver during the mode transition, achieving a smooth clutch engagement at the ending of the mode transition and ensuring a low level of residual torque oscillation after the mode transition.These are the most significant challenges when designing a mode transition controller.Conventionally, vehicle manufacturers define drivability as the driver's subjective perception of the vehicle's powertrain operation.Hence, studies have been devoted to identifying the correlation between measurable vehicle powertrain response indices and the quantification of the driver's subjective perception [8].Furthermore, commercial test tools such as AVL Drive (AVL LIST GmbH, Hans-List-Platz 1, Graz, Austria) can automatically measure a vehicle's powertrain response indices and then infer the overall drivability rating [9].In this paper, we primarily discuss the design of the mode transition controller, and therefore the term 'drivability' directly denotes a vehicle's powertrain response indices.These indices reflect a vehicle's powertrain response characteristics to the driver's maneuvers: sluggishness, surge and oscillation [10].The most commonly used vehicle powertrain response indices are the duration and the responsiveness and smoothness of the vehicle's longitudinal acceleration [11,12].Given that the external resistance can be considered to be constant during the mode transition, the responsiveness and smoothness of the vehicle's longitudinal acceleration are equivalent to those of the propulsion torque exerted on the drivetrain by the hybrid power source.Hence, in this paper, we choose the responsiveness and smoothness of the propulsion torque output by the hybrid power source in combination with the duration of the mode transition as the vehicle's powertrain response indices instead of those of the vehicle's longitudinal acceleration.The process of clutch engagement is composed of the preparation phase, the inertia phase and the post phase, as illustrated in Figure 2.
In a CPPHP, the inertia phase is the mode transition from the pure electric propulsion mode to the hybrid propulsion mode.Hence, the key to excellent mode transition control is coordinating the engine torque, clutch torque and motor torque to accomplish the aforementioned multiple objectives simultaneously during the inertia phase.Various approaches based on classical and modern control theories have been used to solve the torque coordination control problem.Koprubasi et al. proposed an output feedback control based asymptotic stable controller to smooth the clutch engagement in the CPPHP [13].Kim et al. presented a disturbance compensation based torque coordination control strategy to alleviate the torque oscillation during the inertia phase in a CPPHP [14].Chen et al. introduced an output feedback control based torque coordination controller to achieve a smooth clutch engagement at the ending of the inertia phase in the clutch-based pre-transmission series- Conventionally, vehicle manufacturers define drivability as the driver's subjective perception of the vehicle's powertrain operation.Hence, studies have been devoted to identifying the correlation between measurable vehicle powertrain response indices and the quantification of the driver's subjective perception [8].Furthermore, commercial test tools such as AVL Drive (AVL LIST GmbH, Hans-List-Platz 1, Graz, Austria) can automatically measure a vehicle's powertrain response indices and then infer the overall drivability rating [9].In this paper, we primarily discuss the design of the mode transition controller, and therefore the term 'drivability' directly denotes a vehicle's powertrain response indices.These indices reflect a vehicle's powertrain response characteristics to the driver's maneuvers: sluggishness, surge and oscillation [10].The most commonly used vehicle powertrain response indices are the duration and the responsiveness and smoothness of the vehicle's longitudinal acceleration [11,12].Given that the external resistance can be considered to be constant during the mode transition, the responsiveness and smoothness of the vehicle's longitudinal acceleration are equivalent to those of the propulsion torque exerted on the drivetrain by the hybrid power source.Hence, in this paper, we choose the responsiveness and smoothness of the propulsion torque output by the hybrid power source in combination with the duration of the mode transition as the vehicle's powertrain response indices instead of those of the vehicle's longitudinal acceleration.The process of clutch engagement is composed of the preparation phase, the inertia phase and the post phase, as illustrated in Figure 2.
In a CPPHP, the inertia phase is the mode transition from the pure electric propulsion mode to the hybrid propulsion mode.Hence, the key to excellent mode transition control is coordinating the engine torque, clutch torque and motor torque to accomplish the aforementioned multiple objectives simultaneously during the inertia phase.Various approaches based on classical and modern control theories have been used to solve the torque coordination control problem.Koprubasi et al. proposed an output feedback control based asymptotic stable controller to smooth the clutch engagement in the CPPHP [13].Kim

et al. presented a disturbance compensation based torque coordination control
Energies 2016, 9, 740 3 of 31 strategy to alleviate the torque oscillation during the inertia phase in a CPPHP [14].Chen et al. introduced an output feedback control based torque coordination controller to achieve a smooth clutch engagement at the ending of the inertia phase in the clutch-based pre-transmission series-parallel hybrid powertrain [15].Kum et al. designed a dynamic programming based engine-start controller to quickly start the engine during the inertia phase in the CPPHP [16].Torque coordination during the mode transition involves clutch control issues, and the partial objectives are identical to those of clutch control in automated manual transmissions and dual clutch transmissions.Many studies have proposed different solutions for clutch control issues.Dolcini et al. implemented a quadratic programming based lurch avoidance clutch control strategy in an automated manual transmission prototype [17].Walker and Zhang designed a proportional-integralderivative controller using the speed difference between the differential and the wheel as the input to suppress the residual torque oscillation after the inertia phase of an up shift [18].Lu et al. proposed a data-driven predictive control based gearshift control strategy to track the clutch speed difference reference trajectory by manipulating the oncoming clutch torque and engine output torque during the inertia phase of a clutch-to-clutch shift [19].
Model predictive control is adept at managing multiple objective optimization problems with constraints.Some studies have presented model predictive controllers to solve the torque coordination control problem during the mode transition in CPPHPs.Beck et al. designed a model predictive control based triple controller to coordinate the engine and motor torques during the clutch engagement [20].Pisaturo et al. utilized a composite model predictive controller to track the predefined trajectories of the engine speed and clutch speed during the inertia phase in an automated manual transmission.Numerical results validated this controller in critical operating conditions with constraints [21].
The aforementioned studies give valuable insights into the torque coordination control problem during the mode transition in CPPHPs.However, the following details can be improved: the stiffness and damping of the elastic components included in the hybrid powertrain should be modeled; the response characteristics of the engine, clutch and motor should be modeled; in the clutch inertia phase control, the responsiveness of the propulsion torque output by the hybrid power source should be regarded as one drivability-related index; and before the instant of the clutch engagement, the clutch sliding speed and speed difference between the differential and the wheel should be minimized simultaneously because the speed difference between the differential and the wheel results in residual torque oscillation.
These modeling improvements and novel drivability-related indices and constraints are all taken into consideration in the design of the discrete-time model predictive controller in this paper.These elements haven't been comprehensively and systematically taken into account in other relevant studies before, whereas these details indeed influence drivability.
Modeling the stiffness and damping of the elastic components included in the hybrid powertrain is the key to capturing the residual torque oscillation during the post phase; modeling the response characteristics of the engine, clutch and motor is the key to describing the change of the propulsion Numerical results validated this controller in critical operating conditions with constraints [21].
The aforementioned studies give valuable insights into the torque coordination control problem during the mode transition in CPPHPs.However, the following details can be improved: the stiffness and damping of the elastic components included in the hybrid powertrain should be modeled; the response characteristics of the engine, clutch and motor should be modeled; in the clutch inertia phase control, the responsiveness of the propulsion torque output by the hybrid power source should be regarded as one drivability-related index; and before the instant of the clutch engagement, the clutch sliding speed and speed difference between the differential and the wheel should be minimized simultaneously because the speed difference between the differential and the wheel results in residual torque oscillation.
These modeling improvements and novel drivability-related indices and constraints are all taken into consideration in the design of the discrete-time model predictive controller in this paper.These Energies 2016, 9, 740 4 of 31 elements haven't been comprehensively and systematically taken into account in other relevant studies before, whereas these details indeed influence drivability.
Modeling the stiffness and damping of the elastic components included in the hybrid powertrain is the key to capturing the residual torque oscillation during the post phase; modeling the response characteristics of the engine, clutch and motor is the key to describing the change of the propulsion torque output by the hybrid power source during the inertia phase; the responsiveness of the propulsion torque output by the hybrid power source is the key index to evaluate drivability during the inertia phase; the clutch sliding speed influences the torque step at the instant of the clutch engagement; and the speed difference between the differential and the wheel influences the residual torque oscillation during the post phase.These are important and necessary for drivability during the whole process of the clutch engagement, so should be contemplated comprehensively.
Furthermore, by using discrete-time Laguerre functions and introducing the equilibrium state and the ranking of constraints, an explicit solution of the discrete-time model predictive controller is obtained in this paper.The explicit solution can simplify the programming and guarantee the existence of the solution of the optimization problem.Consequently, the explicit solution benefits the attempt to implement the designed discrete-time model predictive controller on an embedded micro control unit in future studies.The key parameters of the discrete-time model predictive controller are determined by comparison with a discrete-time linear quadratic regulator (DLQR) minimizing the almost identical cost function but without constraints.
The outline of this paper is as follows: Section 2 describes the configuration of the CPPHP discussed in this paper; Section 3 describes the simulation-oriented and control-oriented models of the CPPHP; Section 4 details the design of the discrete-time model predictive controller; the parameters and performance of the designed controller are discussed in Section 5 and conclusions are drawn in Section 6.

Powertrain Configuration
The CPPHP schematically depicted in Figure 1c is primarily composed of the engine, clutch, motor, transmission, final drive, differential, axle shafts, tires and related accessories.In comparison with the conventional powertrain configuration with an automatic transmission, the CPPHP configuration substitutes a clutch for the torque converter and incorporates a motor into the primary shaft of the transmission.Other necessary components that must be added to the powertrain include the battery, motor drive and corresponding controllers.Obviously, most components and technologies of the conventional powertrain apply to the CPPHP, and structure-related issues include how to arrange the clutch and the motor in the compact space.Bosch and FEV have completely addressed this issue and produced respective transmission integration solutions.

Simulation-Oriented Model
The simulation-oriented model describes the dynamics of the main components of the CPPHP as comprehensively as possible and functions as the plant in the simulation.The simulation-oriented model will be simplified later to derive the control-oriented model.The dynamics model of the CPPHP is illustrated in Figure 3.
Energies 2016, 9, 740 4 of 30 torque output by the hybrid power source during the inertia phase; the responsiveness of the propulsion torque output by the hybrid power source is the key index to evaluate drivability during the inertia phase; the clutch sliding speed influences the torque step at the instant of the clutch engagement; and the speed difference between the differential and the wheel influences the residual torque oscillation during the post phase.These are important and necessary for drivability during the whole process of the clutch engagement, so should be contemplated comprehensively.Furthermore, by using discrete-time Laguerre functions and introducing the equilibrium state and the ranking of constraints, an explicit solution of the discrete-time model predictive controller is obtained in this paper.The explicit solution can simplify the programming and guarantee the existence of the solution of the optimization problem.Consequently, the explicit solution benefits the attempt to implement the designed discrete-time model predictive controller on an embedded micro control unit in future studies.The key parameters of the discrete-time model predictive controller are determined by comparison with a discrete-time linear quadratic regulator (DLQR) minimizing the almost identical cost function but without constraints.
The outline of this paper is as follows: Section 2 describes the configuration of the CPPHP discussed in this paper; Section 3 describes the simulation-oriented and control-oriented models of the CPPHP; Section 4 details the design of the discrete-time model predictive controller; the parameters and performance of the designed controller are discussed in Section 5 and conclusions are drawn in Section 6.

Powertrain Configuration
The CPPHP schematically depicted in Figure 1c is primarily composed of the engine, clutch, motor, transmission, final drive, differential, axle shafts, tires and related accessories.In comparison with the conventional powertrain configuration with an automatic transmission, the CPPHP configuration substitutes a clutch for the torque converter and incorporates a motor into the primary shaft of the transmission.Other necessary components that must be added to the powertrain include the battery, motor drive and corresponding controllers.Obviously, most components and technologies of the conventional powertrain apply to the CPPHP, and structure-related issues include how to arrange the clutch and the motor in the compact space.Bosch and FEV have completely addressed this issue and produced respective transmission integration solutions.

Simulation-Oriented Model
The simulation-oriented model describes the dynamics of the main components of the CPPHP as comprehensively as possible and functions as the plant in the simulation.The simulation-oriented model will be simplified later to derive the control-oriented model.The dynamics model of the CPPHP is illustrated in Figure 3.The response characteristics of the internal combustion engine can be modeled by a first order system with a time constant .Thus, the engine output torque is calculated by: 1 1 The response characteristics of the internal combustion engine can be modeled by a first order system with a time constant τ e .Thus, the engine output torque M e is calculated by: .
with the engine time constant τ e and the engine output torque command M ec .The rotational components of the engine and the associated driving parts of the clutch can be modeled as a lumped moment of inertia J epc .The rotational speed of this lumped moment of inertia ω epc is computed by: .
with the lumped moment of inertia J epc , the engine output torque M e and the torque transmitted by the clutch M c .

Clutch
The clutch functions as a torque coupler in the CPPHP.The clutch actuator can be modeled by a first order system with a time constant τ c .Then, the normalized normal force r c is calculated by: .
with the clutch actuator time constant τ c and the normalized normal force command r cc .Here r c ∈ [0, 1] and r cc ∈ [0, 1].The clutch can be modeled by a LuGre friction model.Different from the hybrid models usually adopted in the literature, the LuGre friction model is a continuous non-linear model [22].Thus, the clutch is modeled as: .
where z c denominates the average deflection of the bristles, R c the equivalent radius of the clutch, ω epc the rotational speed of the lumped moment of inertia J epc , ω cpmpp the rotational speed of J cpmpp , J cpmpp the lumped moment of inertia of the driven parts of the clutch, the rotor of the motor and the primary shaft of the transmission, σ 0 the stiffness coefficient, F Coulomb the Coulomb force, F stiction the stiction force, v Stribeck the Stribeck velocity, M c the torque transmitted by the clutch, r c the normalized normal force, σ 1 the microscopic damping coefficient and σ 2 the macroscopic damping coefficient.By using the Matlab functions 'max' and 'min', the normalized normal force command r cc is computed by: with the clutch torque command M cc .

Motor
The motor between the clutch and the transmission in the pre-transmission parallel hybrid powertrain architecture is the other power source of the hybrid powertrain.The response characteristics of the motor can be modeled by a first order system with a time constant τ m .Then, the motor output torque M m is calculated by: .
with the motor time constant τ m and the motor output torque command M mc .

Transmission
The transmission used in the CPPHP can be modeled as rigid gear pairs.The driven parts of the clutch, the rotor of the motor and the primary shaft of the transmission can be modeled as a lumped moment of inertia J cpmpp .The rotational speed of this lumped moment of inertia ω cpmpp is computed by: .
with the lumped moment of inertia J cpmpp , the torque transmitted by the clutch M c , the motor output torque M m and the torque exported from the primary shaft M s .The secondary shaft of the transmission can be modeled as a lumped moment of inertia J s .Generally, in terms of the state of the clutch, the primary modes of the CPPHP consist of the pure electric propulsion mode, the hybrid propulsion mode and the regenerative braking mode.When torque demanded by the driver is positive, the CPPHP is switched to the pure electric propulsion mode or the hybrid propulsion mode; and when torque demanded by the driver is negative, the CPPHP is switched to the regenerative braking mode.In the pure electric propulsion mode and the hybrid propulsion mode, the torque exported from the primary shaft M s is positive and for a given M s , the external resistance torque M r results in acceleration or deceleration.In the regenerative braking mode, the clutch is open and the torque exported from the primary shaft M s is negative.
In this paper, we primarily discuss the mode transition from the pure electric propulsion mode to the hybrid propulsion mode and this is the key mode transition of the clutch-based pre-transmission parallel hybrid electric vehicles.During this mode transition, ideally, the propulsion torque exported from the primary shaft M s by the hybrid power source (the engine and the motor) should sustain unchanged by applying a well-designed controller.Hence, although the motor output torque may be negative, the total propulsion torque exported from the primary shaft M s by the hybrid power source (the engine and the motor) will sustain positive and the direction of the power flow is from the hybrid power source to the drivetrain.Therefore, the rotational speed of this lumped moment of inertia ω s is calculated by: .
with the lumped moment of inertia J s , the gear ratio i g , the transmission efficiency η g , the torque exported from the primary shaft M s , and the torque exported from the secondary shaft M f d .ω cpmpp and ω s are linked by: 3.1.5.Final Drive Assuming that driving axle shafts rotate at identical speeds, the final drive and the differential can be modeled as a pair of gears.They can be modeled as a lumped moment of inertia J f d , and given Energies 2016, 9, 740 7 of 31 that the direction of the power flow is from the hybrid power source to the drivetrain during the mode transition, the rotational speed of this lumped moment of inertia ω f d is computed by: .
with the lumped moment of inertia J f d , the final drive ratio i f d , the final drive efficiency η f d , the torque exported from the secondary shaft M f d and the torque transmitted by the axle shafts M as .ω s and ω f d are linked by:

Driving Axle Shaft
Driving axle shafts can be modeled by an equivalent elastic shaft as: .
where θ denotes the torsion of the axle shaft, ω f d the rotational speed of the lumped moment of inertia J f d , ω w the wheel rotational speed, M as the torque transmitted by the axle shafts, k the equivalent stiffness coefficient and b the equivalent damping coefficient.

Longitudinal Vehicle Dynamics
Assuming that there is no slip between the tires and the road, in accordance with the power equation, the dynamics of the vehicle is modeled as [23]: .
where ω w denominates the wheel rotational speed, m the vehicle mass, R w the wheel radius, J the entire moment of inertia of the vehicle body and the tires, M as the torque transmitted by the axle shafts, M r the external resistance torque, g the gravitational acceleration, α the road slope, C rr the rolling resistance coefficient, C D the drag coefficient, A the frontal area, ρ the air density and v the vehicle longitudinal velocity.The engine, clutch actuator and motor time constants are set to be reasonable values of the correct order of magnitude by analyzing the mean value model of the engine, by referencing data provided by the manufacturer, and based on a preliminary test of a permanent magnet synchronous motor respectively.The equivalent radius of the clutch is calculated based on the specifications provided by the manufacturer.The Coulomb force, Stiction force, Stribeck velocity and macroscopic damping coefficient are determined by fitting the LuGre friction model to the characteristic curves provided by the manufacturer.The stiffness coefficient and microscopic damping coefficient are determined by analyzing the characteristics of the LuGre friction model in the literature [22].The moment of inertia of the rotational components of the engine and the parameters of the transmission, final drive, driving axle shafts and longitudinal vehicle dynamics models are determined by referencing the values of that in the literature [24,25] and by collecting the available specifications of the vehicle.The parameters of the simulation-oriented model are shown in Table 1.The simulation result of a Matlab SimDriveline (The MathWorks, Inc., Natick, MA, USA) model is used to validate the proposed model.The model depicted in Figure 4 is based on the standard component models provided by the Matlab SimDriveline.
The scenario to be simulated is: the engine output torque command M ec steps from 0 Nm to 60 Nm at the beginning of the simulation; the initial engine output torque M e0 is 0 Nm; the initial engine rotational speed ω epc0 is 285 rad/s; the normalized normal force command r cc gradually increases from 0 to 1 during the simulation; the initial normalized normal force r c0 is 0; the motor output torque command M mc steps from 50 Nm to 10 Nm at the beginning of the simulation; the initial motor output torque M m0 is 50 Nm; and the initial motor rotational speed ω cpmpp0 is 280 rad/s.
The comparison of the simulation results is shown in Figure 5.In Figure 5a, the RMSE of ω epc is 3.292 rad/s, the RMSE of ω cpmpp is 3.574 rad/s, and the deviation is primarily due to the Matlab SimDriveline model adopting the more detailed tire slip model.In Figure 5b, the RMSE of v is 0.008 m/s, and the vehicle longitudinal velocity from the proposed model is almost identical with that from the Matlab SimDriveline model.Thus, the comparison validates the proposed model in this paper.

Control-Oriented Model
The control-oriented model homes in on the most relevant dynamics of the main components in the CPPHP and is used to design the controller.Compact yet accurate models simplify the controller design.
From the simulation-oriented model to the control-oriented model, the core simplification is related to the clutch model.Although the simulation-oriented clutch model described by Equations ( 4)-( 6) is dynamic and non-linear, the stiffness coefficient σ 0 is usually enormously large.This implies that the average deflection z c reaches the equilibrium state much faster than the clutch sliding speed ω epc − ω cpmpp .Therefore, omitting the process of the average deflection z c tending toward the equilibrium state has almost no influence on the dynamics of the clutch sliding speed ω epc − ω cpmpp and thereby the control-oriented clutch model is almost consistent with the simulation-oriented clutch model.Therefore, the control-oriented clutch model can be refined as: The torque transmitted by the clutch M c is nonnegative means that the torque transmitted by the clutch propels the drivetrain.During the mode transition, guaranteeing that M c is nonnegative can sustain the direction of the torque transmitted by the clutch unchanged.Consequently, the responsiveness and smoothness of the propulsion torque output by the hybrid power source can be improved.To guarantee that the torque transmitted by the clutch M c is nonnegative during the mode transition, the clutch sliding speed ω epc − ω cpmpp should be constrained by Additionally, to guarantee that Equation ( 7) can be simplified further, the clutch torque command M cc is restricted by: The clutch actuator time constant τ c is relatively small and therefore the normalized normal force r c reaches the equilibrium state much faster than the clutch sliding speed ω epc − ω cpmpp .Taking the two constraints (20) and ( 21) into consideration and incorporating the control-oriented clutch model (19) and the clutch actuator controller (7) into the simulation-oriented clutch actuator model (3), the entire control-oriented clutch model can be formulated as: .
In practical operation, the transmission efficiency η g and the final drive efficiency η f d are state-dependent and time-variant, and the related mechanisms are sophisticated and difficult to model.However, the variation is small, so these efficiencies are treated as constants.

Optimization Objectives
In this paper, the duration of the mode transition and the responsiveness and smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source are chosen as powertrain response indices to evaluate the drivability.The level of torque difference between the supplied by the hybrid power source and the demanded by the driver reflects the responsiveness of the propulsion torque output by the hybrid power source.Thus the responsiveness of propulsion torque exerted on the drivetrain by the hybrid power source can be quantified by: where J c1 denotes the responsiveness-related cost function, k the sampling instant, N the horizon, M c (k) the torque transmitted by the clutch at sampling instant k, M m (k) the motor output torque at sampling instant k, M ceq (k) the equilibrium state of the torque transmitted by the clutch at sampling instant k and M meq (k) the equilibrium state of the motor output torque at sampling instant k.The term 'equilibrium state' denotes the hybrid powertrain state satisfying the sufficient and necessary ideal synchronization condition.The equilibrium state will be explained in detail later.
With respect to the smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source, the GV no-lurch condition claims that the time derivative of the clutch sliding speed .ω epc − .ω cpmpp equaling zero at the clutch engagement instant guarantees the smoothness.In accordance with the basic principle of clutch engagement control, before the clutch engagement instant, the clutch sliding speed ω epc − ω cpmpp should approach zero as close as possible.Thus, the GV no-lurch condition can be expressed as: .
In [17], the GV no-lurch condition is extended to a sufficient and necessary ideal synchronization condition by adding the following conditions: .
Adding the two conditions ( 27) and ( 28) is aimed at alleviating the residual torque oscillation after the clutch engagement instant.Therefore, the smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source can be quantified by: where J c2 denominates the smoothness-related cost function, k the sampling instant, N the horizon, ω epc (k) − ω cpmpp (k) the clutch sliding speed at sampling instant k, ω 1eq (k) the equilibrium state of the clutch sliding speed at sampling instant k, − ω w (k) the time derivative of torsion of the axle shaft at sampling instant k and ω 2eq (k) the equilibrium state of the time derivative of torsion of the axle shaft at sampling instant k.

Control Principle
Given the above optimization objectives and the aforementioned constraints that are necessary for assuring the validation of the control-oriented continuous-time state space model and for guaranteeing the safe operation of the related components, the discrete-time model predictive control is suitable for this situation.The basic mechanism of the discrete-time model predictive control is that at each sampling instant, based on the measurement or estimation of the state vector and the control-oriented discrete-time state space model of the controlled object, the input vector at the current sampling instant can determine the output vector at the following sampling instant.Thus, for a prediction horizon, the trajectory of the input vector can determine the trajectory of the output vector.Because the control objective can be formulated as the cost function of the set-point trajectory and predicted trajectory of the output vector, within the prediction horizon, the trajectory minimizing the cost function and satisfying all the equality or inequality constraints is the optimum trajectory of the input vector.Because only the first sampling of the optimum trajectory is implemented as the input vector for the current sampling instant, the optimization process is repeated at the next sampling instant [26].In terms of utilization in the mode transition control, the discrete-time model predictive control has two obvious advantages [27].The first advantage is that the periodically updated optimum trajectory of the input vector can readily take the change of the driver's torque demand into account.The second advantage is that the optimization process repeated at each sampling instant takes the powertrain's actual state into account and can mitigate the influence of the model uncertainty.

Equilibrium State
In this paper, the definition of the term 'equilibrium state' is the hybrid powertrain state satisfying the sufficient and necessary ideal synchronization condition.Ideally, at the clutch engagement instant, the hybrid powertrain should reach the equilibrium state and consequently the responsiveness and smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source are guaranteed.In accordance with Equation ( 26), the equilibrium state of the clutch sliding speed ω 1eq should be a positive value slightly larger than zero.This can guarantee the constraint (20) being satisfied against the influence of model uncertainty.Similarly, in accordance with Equation (28), the equilibrium state of the time derivative of torsion of the axle shaft ω 2eq should be chosen equal to zero.The equilibrium state of the engine output torque M eeq should be chosen equal to the engine output torque required by the energy management system of the hybrid powertrain M eems .Similarly, the equilibrium state of the motor output torque M meq should be chosen equal to the motor output torque required by the energy management system of the hybrid powertrain M mems .The equilibrium state of the torsion of the axle shaft θ eq and the equilibrium state of the torque transmitted by the clutch M ceq can be calculated by: Because only the first sampling of the optimum trajectory is implemented as the input vector for the current sampling instant and, in this paper, the sampling interval is ten milliseconds, the external resistance torque M r can be considered to be constant.The value of M r is computed by (17) at each sampling instant, and the control-oriented model of the hybrid powertrain is therefore updated at each sampling instant.The equilibrium state of the state vector of the control-oriented continuous-time state space model of the hybrid powertrain x1 eq can be formulated as: The equilibrium state of the engine output torque command M eceq should be chosen equal to the equilibrium state of the engine output torque M eeq .The engine time constant τ e is approximately ten times as large as the clutch actuator time constant τ c and approximately a hundred times as large as the motor time constant τ m .Thus, the engine responds to the torque command much slower than the clutch and motor.Given these engine output torque response characteristics, in this paper, the engine output torque command M ec is set equal to the equilibrium state of the engine output torque command M eceq .The equilibrium state of the clutch torque command M cceq should be chosen equal to the equilibrium state of the torque transmitted by the clutch M ceq .Similarly, the equilibrium of the motor output torque command M mceq should be chosen equal to the equilibrium state of the motor output torque M meq .By introducing the equilibrium of the state vector and input vector, the control-oriented continuous-time state space model of CPPHP ( 23) can be refined as:

Control-Oriented Discrete-Time State Space Model
The control-oriented discrete-time state space model is a prerequisite for the discrete-time model predictive controller design.Based on Equation (32), and in accordance with matrix methods for inhomogeneous systems, the control-oriented discrete-time state space model of the CPPHP can be formulated as: with A2 = e A1T and B2 = T 0 e A1t dtB1.Where x2 (k + 1) denotes the state vector at sampling instant k + 1, x2 (k) the state vector at sampling instant k, u2 (k) the input vector at sampling time k and T the sampling interval.Because shifting is inhibited during the mode transition, A1 and B1 are constant matrices during the mode transition.Thus, based on A1 and B1, A2 and B2 can be calculated by using the Matlab function 'c2d' with the specified sampling interval T.

Optimization Problem Formulation
Combining the responsiveness-related cost function J c1 with the smoothness-related cost function J c2 , the optimization problem can be formulated as: R 22 and c1 = 1 0 0 0 0 0 .Where J c denominates the cost function, N p the prediction horizon, y (k + i|k) the predicted output vector for sampling instant k + i, Q1 the weighting matrix for the output vector, u2 (k + i|k) the chosen input vector for sampling instant k + i, R the weighting matrix for the input vector, x2 (k + 1|k) the predicted state vector for sampling instant k + 1, ω 1eq (k) the equilibrium state of the clutch sliding speed chosen at sampling instant k and considered to be constant during the optimization process of sampling instant k, M cceq (k) the equilibrium state of the clutch torque command computed at sampling instant k and considered to be constant during the optimization process of sampling instant k, M ccmax (k) the maximum allowable clutch torque command calculated at sampling instant k in accordance with (21), M mcmin and M mcmax the minimum and maximum allowable motor output torque commands and M mceq (k) the equilibrium state of the motor output torque command computed at sampling instant k and considered to be constant during the optimization process of sampling instant k.In accordance with the basic mechanism of the discrete-time model predictive control, only the first sampling of the optimum trajectory of the input vector is implemented, so only the constraints on the first sampling of the trajectory of the input vector have relatively major effects [26].Thus the optimization problem formulation only includes constraints (35) and (36).The cost function (34) can be refined as: with Q = C1 T Q1C1.
The optimization problem can be simplified by introducing discrete-time Laguerre functions.The work presented in this paper is a prerequisite for an attempt to implement the designed discrete-time model predictive controller on an embedded micro control unit in future studies, so the computational load is sensitive.Theoretically, the number of parameters to be calculated directly influences the computational load of the optimization problem.The task of the discrete-time model predictive controller is solving the aforementioned optimization problem at each sampling instant.Solving the optimization problem is equivalent to finding the optimum trajectory of the input vector u2 (k + i|k) for a specified prediction horizon N p .Obviously, the number of parameters used to represent the trajectory increases rapidly with the increase of the prediction horizon N p .Given that the duration of the mode transition is approximately one second and the sampling interval in this paper is ten milliseconds, the prediction horizon in this paper is set to be 100.Because this large prediction horizon is needed for high-performance, using discrete-time Laguerre functions can enormously decrease the number of parameters and consequently reduce the computational load.The discrete-time Laguerre functions can be formulated in a vector form as described in [26]: with: . Where a ∈ [0, 1) denotes the scaling factor and N the number of discrete-time Laguerre functions.The inherent orthonormality of discrete-time Laguerre functions can be formulated as in [26]: Due to the asymptotic stability of (39), discrete-time Laguerre functions can be used to represent the trajectory of the input vector and u2 (k + i|k) can be formulated as: where l1 (i) denominates the set of discrete-time Laguerre functions calculated in accordance with (38)  and (39) with a 1 and N 1 as parameters, l2 (i) the set of discrete-time Laguerre functions computed in accordance with (38) and (39) with a 2 and N 2 as parameters and η a column vector with N 1 + N 2 elements.In accordance with the control-oriented discrete-time state space model (33), the predicted state vector for sampling instant k + i can be calculated by: Substituting ( 41) into ( 42), (42) can be refined further as: Incorporating (43) into (37), the cost function (37) can be refined further as: with The reformulation of the second term of (37) is based on the orthonormality of discrete-time Laguerre functions.Thus the aforementioned optimization problem can be reformulated as:

Optimization Solution
The control law of the discrete-time model predictive controller can be derived from the solution of the optimization problem.Essentially, this optimization problem is a quadratic programming problem.Similarly, the work presented in this paper is a prerequisite for an attempt to implement the designed discrete-time model predictive controller on an embedded micro control unit in future studies, so the convenience of programming and the existence of the solution are sensitive.Developing an explicit solution of the optimization problem can simplify the programming and guarantee the existence of the solution of the optimization problem.Due to the existence of inequality constraints, numerical methods for solving the quadratic programming problem need to identify inactive constraints and then ignore these constraints by introducing Lagrange multipliers.Obviously, m1 and rows of M1 and M2 are not linearly independent, so there may be no solution for Lagrange multipliers.Introducing the ranking of constraints can solve this problem and reduce the computational load [26].According to their importance, the constraints on the input vector have priority over that on the state vector because the former relates to the safe operation of the related components, whereas model uncertainty could interfere with the latter.Thus the inequality constraints (47) and ( 48) have a higher priority whereas the inequality constraint (46) has a lower priority.Taking the ranking of constraints into consideration, the procedure for solving the optimization problem can be formulated as the following:

•
Solve the optimization problem without inequality constraints: • Check whether η obtained in Step 1 satisfies the inequality constraint (46).If not, by introducing a Lagrange multiplier, compute η, the clutch torque command M cc (k) and the motor output torque command M mc (k): • Check whether or not η obtained in Step 2 satisfies the inequality constraints (47) and (48).If not, calculate the clutch torque command M cc (k) and the motor output torque command M mc (k): Some parameters of the designed controller are given as the following: the equilibrium state of the clutch sliding speed ω 1eq is set to be 0.5 rad/s; the maximum motor output torque command M mcmax is set to be 50 Nm; and the minimum motor output torque command M mcmin is set to be −50 Nm.
The comparison of the computational times between the proposed quadratic programming solver and the Matlab standard quadratic programming solver 'quadprog' is shown in Table 2.The comparison of the obtained optimum solutions is shown in Table 3.The proposed solver and the Matlab standard solver provide identical optimum solutions.The optimization problem to be solved is calculating the optimum input vector for the first sampling instant of the inertial phase and the values of the parameters in the solved quadratic programming problem are set equal to those chosen in Section 5.For each solver, the solving process is repeated five times for objectivity.The comparison is conducted in the Matlab programming environment in a laptop.The comparison in Tables 2 and 3 shows the obvious advantages of the proposed quadratic programming solver in this paper.

Parameter Selection
The key parameters to be determined include a 1 and N 1 for generating discrete-time Laguerre functions used to represent the trajectory of the clutch torque command M cc and a 2 and N 2 for discrete-time Laguerre functions related to the trajectory of the motor output torque command M mc .In accordance with (38) and (39), the scaling factor a determines the rate of decay of discrete-time Laguerre functions and thereby the linear combination of these discrete-time Laguerre functions that is used to represent the trajectory of the input vector u2 (k + i|k) in accordance with (41).Ideally, the clutch torque command M cc and the motor output torque command M mc should reach their respective equilibrium states simultaneously, so a 1 is chosen to be equal to a 2 .Furthermore, for the sake of brevity, we choose N 1 to be equal to N 2 .With respect to the cost function (37), choose the prediction horizon N p to be equal to infinity, and then the cost function can be formulated as: The solution minimizing the cost function ( 59) is referred to as a discrete-time linear quadratic regulator.For arbitrary a 1 ∈ [0 , 1) and a 2 ∈ [0 , 1), as N 1 and N 2 increase, the optimum trajectory of the input vector obtained by minimizing the cost function (45) will converge to that obtained by minimizing the cost function (59) [26].DLQR is not suitable to address the constraints, whereas DLQR is the ideal state feedback controller in the case of no constraints.DMPC is inherently adept at addressing the constraints.Thus, DLQR can be used to find reasonable parameters for the designed discrete-time model predictive controller and thereby the obtained DMPC can achieve high-performance similar to DLQR when the constraints become not active.The standard procedure for designing a discrete-time linear quadratic regulator is:

•
Firstly, solve the following Riccati equation for finding the Riccati matrix P ∞ : • Secondly, compute the state feedback gain matrix K: • Thirdly, calculate the optimum input vector for the current sampling instant: Practically, the state feedback gain matrix K can be computed by using the Matlab function 'dlqr' with the state matrix A2, the input matrix B2 and the weighting matrices Q and R. In the case of Q = C1 T Q1C1 and R = I 2×2 , the effects of a 1 , a 2 , N 1 and N 2 on the difference between the optimum trajectory of the input vector determined by DMPC and that determined by DLQR are shown in Figures 6-8.
Practically, the state feedback gain matrix can be computed by using the Matlab function 'dlqr' with the state matrix , the input matrix and the weighting matrices and .In the case of and , the effects of , , and on the difference between the optimum trajectory of the input vector determined by DMPC and that determined by DLQR are shown in Figures 6-8       These results demonstrate that for a given scaling factor, as N 1 and N 2 increase, the optimum trajectory calculated by DMPC converges to that computed by DLQR, and for given N 1 and N 2 , too small or too large scaling factors result in the optimum trajectory calculated by DMPC deviating from that computed by DLQR.Based on these observations, comprehensively considering the need for accurately capturing the optimum trajectory determined by DLQR and the need for reducing the computational load, in this paper we choose a 1 = a 2 = 0.9 and N 1 = N 2 = 5.
The weighting matrix for the input vector R is determined in the analysis of the performance of the designed discrete-time model predictive controller in the following section.

Performance in Nominal Conditions
The nominal conditions mean that the parameters used in the simulation-oriented model are identical with those used in the control-oriented model.The designed discrete-time model predictive controller originates from the control-oriented model.The scenario to be simulated is: firstly, the pre-transmission parallel hybrid powertrain operates in the pure electric propulsion mode; due to increase in the driver's torque demand to address the extra resistance temporarily existing before and during the preparation phase, the energy management system claims that the hybrid powertrain should switch to the hybrid propulsion mode; thus, the engine is started and operates in the speed control mode, and the preparation phase begins; when the event that the engine speed falls into the specified range that is slightly larger than the motor speed triggers the designed discrete-time model predictive controller, the preparation phase ends and the inertia phase begins; the relevant state and input variables reaching their respective equilibrium states is the symbol of the ending of the inertia phase, and then the post phase begins; and the instant the post phase begins, the engine output torque command and the motor output torque command are set equal to their respective values required by the energy management system of the hybrid powertrain and the clutch torque command gradually increases to the maximum.The initial values of the relevant variables at the beginning of the inertia phase are shown in Table 4.The calculation results shown in Figures 9-12 verify that for various combinations of weighting matrices the clutch torque command and the motor output torque command converge to their respective predefined equilibrium states simultaneously, and the clutch sliding speed and the time derivative of torsion of the axle shaft gradually tend toward their respective predefined equilibrium states.Given that the engine output torque command is set equal to the predefined equilibrium state during the mode transition and the external resistance torque can be considered to be constant during the sampling interval, the sufficient and necessary ideal synchronization condition has been satisfied, and thus the smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source is validated.Additionally, the torque difference between the supplied by the hybrid power source and the demanded by the driver, shown in Figures 9-12, demonstrates the high level of responsiveness of the propulsion torque exerted on the drivetrain by the hybrid power source.M ceq + M meq denotes the torque demanded by the driver, and M c + M m denominates the torque supplied by the hybrid power source.With respect to the duration of the mode transition, the convergence speed of relevant objects illustrated in Figures 9-12 guarantees that the mode transition can be completed within an acceptable duration.The comparison of performance with different input vector weighting matrices confirms that the decrease of the input vector weighting matrix can augment the responsiveness and smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source.However, when the input vector weighting matrix decreases to a relatively small level, the effects of any further decrease are limited because the chosen scaling factors compromise on accurately capturing the optimum trajectory determined by DLQR.Taking these observations into account, in this paper we choose R = 0.01 0 0 0.01 .
the propulsion torque exerted on the drivetrain by the hybrid power source.However, when the input vector weighting matrix decreases to a relatively small level, the effects of any further decrease are limited because the chosen scaling factors compromise on accurately capturing the optimum trajectory determined by DLQR.Taking these observations into account, in this paper we choose 0.01 0 0 0.01 .

Comparison with Conditions Without DMPC
To validate the advantages of the designed discrete-time model predictive controller, the calculation results of the scenario without DMPC are shown in Figure 13.In this scenario, when the engine speed falls into the specified range that is slightly larger than the motor speed, the engine output torque command and the motor output torque command are set equal to the respective values required by the energy management system of the hybrid powertrain, and the clutch torque command gradually increases to the maximum in the same way that it changes in the scenario with DMPC.

Comparison with Conditions Without DMPC
To validate the advantages of the designed discrete-time model predictive controller, the calculation results of the scenario without DMPC are shown in Figure 13.In this scenario, when the engine speed falls into the specified range that is slightly larger than the motor speed, the engine output torque command and the motor output torque command are set equal to the respective values required by the energy management system of the hybrid powertrain, and the clutch torque command gradually increases to the maximum in the same way that it changes in the scenario with DMPC.

Comparison with Conditions Without DMPC
To validate the advantages of the designed discrete-time model predictive controller, the calculation results of the scenario without DMPC are shown in Figure 13.In this scenario, when the engine speed falls into the specified range that is slightly larger than the motor speed, the engine output torque command and the motor output torque command are set equal to the respective values required by the energy management system of the hybrid powertrain, and the clutch torque command gradually increases to the maximum in the same way that it changes in the scenario with DMPC.Obviously, the torque difference between the supplied by the hybrid power source and the demanded by the driver in the scenario without DMPC is approximately ten times as large as that in the scenario adopting DMPC and the time derivative of torsion of the axle shaft in the scenario without DMPC is approximately six times as large as that in the scenario adopting DMPC.Obviously, the torque difference between the supplied by the hybrid power source and the demanded by the driver in the scenario without DMPC is approximately ten times as large as that in the scenario adopting DMPC and the time derivative of torsion of the axle shaft in the scenario without DMPC is approximately six times as large as that in the scenario adopting DMPC.

Robustness
The designed discrete-time model predictive controller is based on the control-oriented model of the hybrid powertrain, and the control-oriented model may deviate from the actual system due to variations in the driving conditions and parameter uncertainty.Thus, the presented discrete-time model predictive controller should be robust to variations in the driving conditions and parameter uncertainty.
The calculation results shown in Figure 14 demonstrate the robustness of the designed controller in conditions different from the nominal conditions due to variations in the driving conditions.Clearly, during the mode transition, the steady states of the relevant state and input variables deviate from their respective predefined equilibrium values due to variations in the driving conditions.Still, during the mode transition, the designed controller guarantees that the clutch sliding speed and the time derivative of torsion of the axle shaft sustain steady states adjacent to their respective predefined equilibrium states, and the torque difference between the supplied by the hybrid power source and the demanded by the driver tends toward a steady level.Thus, the responsiveness and smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source does not suffer from variations in the driving conditions.The calculation results shown in Figure 15 demonstrate the robustness of the designed controller in conditions different from the nominal conditions due to parameter uncertainty.Obviously, parameter uncertainty of the clutch actuator and motor time constants has almost no influence on the performance of the designed discrete-time model predictive controller.These calculation results indicate that the designed discrete-time model predictive controller is still stable in spite of variations in the driving conditions and parameter uncertainty, thereby validating the robustness of the designed controller.

Robustness
The designed discrete-time model predictive controller is based on the control-oriented model of the hybrid powertrain, and the control-oriented model may deviate from the actual system due to variations in the driving conditions and parameter uncertainty.Thus, the presented discrete-time model predictive controller should be robust to variations in the driving conditions and parameter uncertainty.
The calculation results shown in Figure 14 demonstrate the robustness of the designed controller in conditions different from the nominal conditions due to variations in the driving conditions.Clearly, during the mode transition, the steady states of the relevant state and input variables deviate from their respective predefined equilibrium values due to variations in the driving conditions.Still, during the mode transition, the designed controller guarantees that the clutch sliding speed and the time derivative of torsion of the axle shaft sustain steady states adjacent to their respective predefined equilibrium states, and the torque difference between the supplied by the hybrid power source and the demanded by the driver tends toward a steady level.Thus, the responsiveness and smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source does not suffer from variations in the driving conditions.The calculation results shown in Figure 15

Conclusions
In this paper, we propose a discrete-time model predictive controller applying to the mode transition from the pure electric propulsion mode to the hybrid propulsion mode in pre-transmission parallel hybrid electric vehicles.In the case of a long prediction horizon, due to the reformulation of the trajectory of the input vector by discrete-time Laguerre functions, the number of variables that needs to be computed is greatly reduced.By comparing the optimum trajectory of the input vector determined by DMPC with that determined by DLQR, which has the almost identical cost function with DMPC, the reasonable parameters are chosen.The calculation results demonstrate that the designed discrete-time model predictive controller can guarantee the responsiveness and smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source during the mode transition, ensuring a high level of drivability.Additionally, the calculation results prove that the performance of the designed discrete-time model predictive controller does not suffer from variations in the driving conditions and parameter uncertainty.In future studies, the designed discrete-time

Conclusions
In this paper, we propose a discrete-time model predictive controller applying to the mode transition from the pure electric propulsion mode to the hybrid propulsion mode in pre-transmission parallel hybrid electric vehicles.In the case of a long prediction horizon, due to the reformulation of the trajectory of the input vector by discrete-time Laguerre functions, the number of variables that needs to be computed is greatly reduced.By comparing the optimum trajectory of the input vector determined by DMPC with that determined by DLQR, which has the almost identical cost function with DMPC, the reasonable parameters are chosen.The calculation results demonstrate that the designed discrete-time model predictive controller can guarantee the responsiveness and smoothness of the propulsion torque exerted on the drivetrain by the hybrid power source during the mode transition, ensuring a high level of drivability.Additionally, the calculation results prove that the performance of the designed discrete-time model predictive controller does not suffer from variations in the driving conditions and parameter uncertainty.In future studies, the designed discrete-time model predictive controller will be refined further by deeply evaluating the importance of the introduced drivability-related indices and constraints in this paper.
hybrid powertrain[15].Kum et al. designed a dynamic programming based engine-start controller to quickly start the engine during the inertia phase in the CPPHP[16].

Figure 5 .
Figure 5.Comparison of the simulation results.(a) Comparison of the engine rotational speeds ω epc and the motor rotational speeds ω cpmpp ; (b) Comparison of the vehicle longitudinal velocities v. .

Figure 6 . 8 .
Figure 6.Difference between the optimum trajectories in the case of 0.8.The sampling interval is ten milliseconds.(a) Difference between the optimum trajectories of the clutch torque command ; (b) Difference between the optimum trajectories of the motor output torque command .

Figure 6 .Figure 7 .
Figure 6.Difference between the optimum trajectories in the case of a 1 = a 2 = 0.8.The sampling interval is ten milliseconds.(a) Difference between the optimum trajectories of the clutch torque command M cc ; (b) Difference between the optimum trajectories of the motor output torque command M mc .Energies 2016, 9, 740 19 of 30

Figure 7 .Figure 7 . 9 .
Figure 7. Difference between the optimum trajectories in the case of a 1 = a 2 = 0.9.The sampling interval is ten milliseconds.(a) Difference between the optimum trajectories of the clutch torque command M cc ; (b) Difference between the optimum trajectories of the motor output torque command M mc .

Figure 8 .
Figure 8. Difference between the optimum trajectories in the case of 0.99.The sampling interval is ten milliseconds.(a) Difference between the optimum trajectories of the clutch torque command ; (b) Difference between the optimum trajectories of the motor output torque command .

Figure 8 .
Figure 8. Difference between the optimum trajectories in the case of a 1 = a 2 = 0.99.The sampling interval is ten milliseconds.(a) Difference between the optimum trajectories of the clutch torque command M cc ; (b) Difference between the optimum trajectories of the motor output torque command M mc .

Figure 9 .
Figure 9. Performance in nominal conditions with a 1 = a 2 = 0.9, N 1 = N 2 = 5, Q1 = I 3×3 and R = I 2×2 .(a) Engine output torque command M ec , clutch torque command M cc and motor output torque command M mc ; (b) Conditions of the constraints; (c) Equilibrium state of the clutch sliding speed ω 1eq and clutch sliding speed ω 1 ; (d) Equilibrium state of the time derivative of torsion of the axle shaft ω 2eq and time derivative of torsion of the axle shaft ω 2 ; (e) Torque demanded by the driver M ceq + M meq and torque supplied by the hybrid power source M c + M m .

Figure 9 .Figure 10 .Figure 10 . 30 Figure 9 .Figure 10 .Figure 11 .
Figure 9. Performance in nominal conditions with 0.9, 5, and .(a) Engine output torque command , clutch torque command and motor output torque command ; (b) Conditions of the constraints; (c) Equilibrium state of the clutch sliding speed and clutch sliding speed ; (d) Equilibrium state of the time derivative of torsion of the axle shaft and time derivative of torsion of the axle shaft ; (e) Torque demanded by the driver and torque supplied by the hybrid power source .

.
(a) Engine output torque command M ec , clutch torque command M cc and motor output torque command M mc ; (b) Conditions of the constraints; (c) Equilibrium state of the clutch sliding speed ω 1eq and clutch sliding speed ω 1 ; (d) Equilibrium state of the time derivative of torsion of the axle shaft ω 2eq and time derivative of torsion of the axle shaft ω 2 ; (e) Torque demanded by the driver M ceq + M meq and torque supplied by the hybrid power source M c + M m .

Figure 13 .
Figure 13.Scenario without DMPC.(a) Engine output torque command , clutch torque command and motor output torque command ; (b) Equilibrium state of the clutch sliding speed and clutch sliding speed ; (c) Equilibrium state of the time derivative of torsion of the axle shaft and time derivative of torsion of the axle shaft ; (d) Torque demanded by the driver and torque supplied by the hybrid power source .

Figure 13 .
Figure 13.Scenario without DMPC.(a) Engine output torque command M ec , clutch torque command M cc and motor output torque command M mc ; (b) Equilibrium state of the clutch sliding speed ω 1eq and clutch sliding speed ω 1 ; (c) Equilibrium state of the time derivative of torsion of the axle shaft ω 2eq and time derivative of torsion of the axle shaft ω 2 ; (d) Torque demanded by the driver M ceq + M meq and torque supplied by the hybrid power source M c + M m .

Figure 15 .
Figure 15.Performance in the case of parameter uncertainty.Case 1: τ c = 0.05 s.Case 2: τ m = 0.005 s.Engine output torque command M ec steps from 0 Nm to 60 Nm at the beginning of the inertia phase.(a) Clutch torque command M cc and motor output torque command M mc ; (b) Conditions of the constraints; (c) Equilibrium state of the clutch sliding speed ω 1eq and clutch sliding speed ω 1 ; (d) Equilibrium state of the time derivative of torsion of the axle shaft ω 2eq and time derivative of torsion of the axle shaft ω 2 ; (e) Torque demanded by the driver M ceq + M meq and torque supplied by the hybrid power source M c + M m .

Table 1 .
Parameters of the simulation-oriented model.

Table 2 .
Comparison of the computational times.

Table 3 .
Comparison of the obtained optimum solutions.

Table 4 .
Initial values of the relevant variables.