Next Article in Journal
Feasibility of Acoustic Print Head Monitoring for Binder Jetting Processes with Artificial Neural Networks
Next Article in Special Issue
Influence of the Final Ratio on the Consumption of an Electric Vehicle under Conditions of Standardized Driving Cycles
Previous Article in Journal
Dairy Waste and Potential of Small-Scale Biogas Digester for Rural Energy in India
Previous Article in Special Issue
Processor-in-the-Loop Architecture Design and Experimental Validation for an Autonomous Racing Vehicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parallel Hybrid Electric Vehicle Modelling and Model Predictive Control

1
Institute for Nanomaterials, Advanced Technologies and Innovation, Technical University of Liberec, 460 01 Liberec, Czech Republic
2
Faculty of Mechatronics, Informatics and Interdisciplinary Studies, Technical University of Liberec, 460 01 Liberec, Czech Republic
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(22), 10668; https://doi.org/10.3390/app112210668
Submission received: 22 September 2021 / Revised: 3 November 2021 / Accepted: 9 November 2021 / Published: 12 November 2021
(This article belongs to the Special Issue Frontiers in Mechatronics Systems for Automotive)

Abstract

:
This paper presents the modelling and calculations for a hybrid electric vehicle (HEV) in parallel configuration, including a main electrical driving motor (EM), an internal combustion engine (ICE), and a starter/generator motor. The modelling equations of the HEV include vehicle acceleration and jerk, so that simulations can investigate the vehicle drivability and comfortability with different control parameters. A model predictive control (MPC) scheme with softened constraints for this HEV is developed. The new MPC with softened constraints shows its superiority over the MPC with hard constraints as it provides a faster setpoint tracking and smoother clutch engagement. The conversion of some hard constraints into softened constraints can improve the MPC stability and robustness. The MPC with softened constraints can maintain the system stability, while the MPC with hard constraints becomes unstable if some input constraints lead to the violation of output constraints.

1. Introduction

Controllers for HEVs powertrains and speeds can be included model-free or model-based. Model-free controllers are mostly used with heuristic, fuzzy, neuro, AI, or human virtual and augmented reality. The use of model-free methods will be presented in the next part of this study. Model-based controllers can be used with a conventional adaptive PID, H2, H, or sliding mode. However, all conventional control methods cannot include the real-time dynamic constraints of the vehicle physical limits, the surrounding obstacles, and the environment (road and weather) conditions. Therefore, a MPC with horizon state and open loop control prediction subject to dynamic constraints are mainly used to control as real-time the HEV speeds and torques. Due to the limit size of this paper, we have reviewed some of the most recent research of MPC applications for HEVs. In this paper, vehicle dynamic formulas and calculations are referred to the reference [1].
A recent modelling and control of the dual clutch transmission for HEVs are presented in [2], in which a new controller was designed for synchronizing the dual clutch transmission (DCT) with higher performances and lower fuel consumption. Another MPC for an autonomous driving vehicle is developed in [3], in which the MPC was used to drive the HEV to exactly track the given feasible trajectories. Moreover, a controller for a hybrid dual clutch transmission powertrain for a HEV is introduced in [4], in which the ICE and EM were driven by a DCT powertrain. A MPC for a HEV with linear parameters and a varying model is presented in [5], in which the MPC controller was designed to improve the fuel economy of the power-split HEV.
A MPC for a HEV is necessary not only to control the torque and speed, but also to control the gas emission and improve the fuel economy. The authors of [6] developed a MPC with a multi-objective function for HEVs for fuel economy and exhaust emission and collision detection as well as to optimize the vehicle speed and engine torque. A new MPC design for HEVs with adaptive cruise of autonomous electric vehicles was presented in [7]. A hybrid MPC to optimize the HEV mode selection was introduced in [8], in which the vehicle thermal management was controlled by this MPC subject to decision-making algorithms. Fuel economy and lower emission were also controlled by a MPC with outer approximation and semi-convex cut generation, as presented in [9].
Due to the recent world commitment to limit the increase in global warming and stop using fossil fuels, plug-in and the pure electric vehicles are expanding. MPC algorithms are also being developed to control the plug-in hybrid vehicle (PHEV), as shown in [10]. In this study, a non-linear MPC is designed to control the torque-split and to optimize fuel management. The authors of [11] also present a data-based scenario MPC framework to optimize power consumption.
Non-linear model predictive control (NMPC) has been widely used due to the rapid increase in computer capacity and speeds. Computers now can calculate a real-time solution from complex non-linear functions. Therefore, the authors of [12] provide a MPC for the non-linear energy management of the power split HEV. Energy efficiency management for HEVs is now also extended in communication among vehicles, as shown in [13], in which a MPC framework is proposed to generate the optimal torque and velocity by connecting the information from vehicle to vehicle.
The authors of [14] reviewed the latest model-based controllers in the market to assess the improvement of energy management for HEVs, and, in this study, the MPC was used to calculate optimal energy, torque, and speed. Since a MPC is a model-based algorithm, difficulties will arise when there are existing mismatches between the model and the plant or the plant uncertainties. These mismatches and uncertainties may lead to the instability of the controller. Robust model predictive control (RMPC) algorithms are, therefore, developed to deal to these uncertainties. The authors of [15] present a new method using matrix inequalities based on RMPC for HEVs considering the external disturbances, the time varying delays, and the model uncertainties. The authors of [16] introduce a real-time NMPC for the energy management of HEVs using sequential quadratic programing.
The use of a MPC for pure electric vehicles is also mentioned in [17] regarding full battery consumption and road slope condition. The authors of [18] presented a decentralized MPC of plug-in electric vehicles charged based on the alternative direction method of multipliers. A real-time control for HEVs, longitudinal tracking, jaw movement, dual-mode power split, and minimizing energy were presented by the authors of [19,20,21,22,23,24,25,26,27,28]. However, none of the recent MPC methods mentioned is concerned with the MPC with softened constraints. The MPC is always subject to many strict constraints on states, outputs, and inputs; therefore, the MPC may not find a feasible solution and it may become unstable. Since the MPC is a real-time optimizer, any failure solution cannot be tolerated. We propose to convert some physical strict constraints into softened constraints, while adding some large penalty values into the MPC objective function. In this way, we can increase the stability and the robustness of the system dealing with uncertainties and some initial conditions, which may lead the input to violate output constraints.
The layout of this paper is as following: Section 2 presents the modelling of the parallel HEV; Section 3 introduces the design of the MPC; Section 4 develops the MPC algorithms with softened constraints; Section 5 illustrates the simulations of the MPC for the HEV; and Section 6 is the conclusion.

2. Modelling of the Parallel HEV

In parallel hybrid electric vehicles, both combustion engine and main electric motor are installed in parallel and work in independent configurations. The vehicle can run and switch in four driving modes: pure main electric motor (EV1) at a low speed and/or low load; pure combustion engine (EV2) at a high speed and/or high load; both main electric motor and combustion engine (EV3) at a very high and/or very high load; and a combination of the main electric motor, generator motor, and combustion engine (EV4) at extreme high speed and/or load. In 2021, Hyundai introduced a new version of Sonata Hybrid, the middle size family passenger HEV that combines updated technologies for a typical parallel electric vehicle, as shown in Figure 1.
This parallel HEV consists of one internal combustion engine (ICE) with four cylinders with multiple point injection, a volume of 2.4 litters, a maximum power of 156 kW at 6000 rpm, and a peak torque of 265 Nm; one electric motor starter (EM2) with a maximum power of 8 kW and a maximum torque of 43 Nm; the main electric motor (EM1) with a maximum power of 35 kW and a maximum torque of 205 Nm; the battery HEV Li-ion with a capacity of 6.1 Ah; a transmission gearbox with fully automated 6 speeds; and a friction clutch engagement. The vehicle curb weight is 1569 kg. This Sonata Hybrid vehicle is used to simulate our system modelling and test the new MPC scheme with softened constraints.
The schematic architecture of the 2021 Hyundai Sonata Hybrid in Figure 1 can be modelled with a simple drivetrain and is shown in Figure 2. The first part of this mechanical structure consists of an internal combustion engine (ICE) and the electric starter/generator motor (EM2) can be grouped into one inertia J 1 , including the left clutch disk, the sharp 1, EM2, and ICE. J 1 is modelled as one rigid inertia. M I C E and M E V 1 are the torques on the ICE and EM2, respectively. θ 1 and ω 1 are the angular position and the velocity of the sharp 1, respectively. Similarly, J 2 is modelled as the lumped rigid inertia of the main electric motor EM1 and the right clutch disk. θ 2 and ω 2 are the angular position and velocity of the sharp 2, respectively. The third powertrain part connecting the gearbox and the vehicle’s driven wheels can be modelled by a gear ratio i via a torsional spring and damper with k θ , k β , and k α as the stiffness, damping, and acceleration coefficient, respectively, of which the acceleration has not been studied before. The third part, the lumped inertia J 3 , consists of the rest of the vehicle, including the gearbox, differential gear, shaft 3, and the driven wheels. θ 3 and ω 3 are the angular position and velocity of shaft 3, respectively. r r is the rolling radius of the vehicle’s wheels.
In this paper, all vehicle dynamic formulas and constraints were taken from the technical book in [1]. The vehicle resistance torque is the approximation of the air density ρ , air drag coefficient c w , the vehicle crossing area A , the wheel rolling radius r r , vehicle friction resistant coefficient f r , natural gravity g , vehicle mass m , and the polynomial coefficients of a 0 , a 1 and a 2 The vehicle rolling resistance torque M v can be calculated as:
M v = ( ρ 2 c w A ( r ω 3 ) 2 + f r m g ) r + a 0 + a 1 ω 3 + a 2 ω 3 2
In Equation (1), the additional road conditions, such as the road dynamics, the road increase, and other environment conditions, can be added as disturbances that lead to some reduction of or increase in the vehicle rolling resistance torque. Changes of vehicle velocity, depending on the road conditions as well as the vehicle dynamic constraints between the vehicle speed and vehicle steering wheel, are referred to in [1].
At a low speed of less than 40 km/h, the clutch is open, and only the main electric motor EM1 propels the HEV. The contribution of some other exponential coefficients is small and can be ignored. The vehicle rolling resistance torque at a low speed can be simplified as:
M v = M v 0 + k v ω 3
where M v 0 is the initial resistance constant of air drag and rolling friction. k v is a linear coefficient that depends on the gear ratio.
On the first part, the torque applied is:
M 1 o = J 1 ω ˙ 1
This torque can be calculated as:
M 1 o = M I C E + M M 2 M C
where M I C E is the torque from ICE; M M 2 is the torque from motor ME2; and M C is the torque from the clutch.
When the clutch is locked, the clutch torque M C is the maximum static clutch friction:
M C = 2 3 r C F N C μ S   when   ( M C = M f m a x S t a t i c )
where r C is the clutch radius; F N C is the normal force; and µS is the clutch friction coefficient.
When the clutch is in the transitional period, M C < M f m a x S t a t i c , the clutch torque is:
M C = r C F N C s i g n ( ω 1 ω 2 ) μ K   when   ( M C < M f m a x S t a t i c )
where μ K is the clutch slipping coefficient.
On the second part, the torque applied on the main motor ME1 is:
M 2 o = k θ θ 2 + k θ i θ 3 + k v ω 3
The sum of inertias is calculated as:
M 2 o = J 2 ω ˙ 2 i + J 3 ω ˙ 3 + k v ω 3
The torque changing is calculated as:
M ˙ 2 o = J 2 ω ¨ 2 i + k α ( ω ˙ 2 i ω ˙ 3 ) + k β ( ω 2 i ω 3 )
The balance of torque M 2 o is calculated as:
M 2 o = ( M E M 2 + M c ) η i M v 0
where η is the transmission efficiency of the gearbox and the differential gear.
The above torque equations can be transformed to the following dynamic equations:
θ ˙ 1 = ω 1
The angular acceleration of the shaft 1 is calculated as:
ω ˙ 1 = k β 1 ω 1 J 1 + M I C E J 1 + M M 1 J 1 + M C J 1
where k β 1 is the shaft 1 friction coefficient.
The angular acceleration of the shaft 2 is calculated as:
ω ˙ 2 = k β 2 ω 3 J 2 i J 3 ω ˙ 3 J 2 i η M M 2 J 2 + η M C J 2 M v 0 J 2 i
where k β 2 is the shaft 2 friction coefficient.
Finally, the angular acceleration of the shaft 3 is calculated as:
ω ˙ 3 = k β 3 ω 3 J 3 + M v 0
where k β 3 is the shaft 3 friction coefficient.
The jerk on the drivetrain is calculated as:
ω ¨ 3 = k β 2 ω 2 J 3 i ( k β 2 J 2 i 2 + k α k v ) ω 3 J 2 J 3 i 2 ( k v + k α J 3 + k α J 2 i 2 ) ω ˙ 3 + k α η ( M M 2 + M C ) J 2 J 3 i k α M v 0 J 2 J 3 i 2
The torque generated on the main motor is calculated as:
M D C _ M O T O R = k T R I V k E k T R I ω M D C _ M O T O R = k T R I V k E k T R I ω
where M D C _ M O T O R is the main motor torque; k T is the motor constant, k T = M T o r q u e I C u r r e n t (Nm/A); k E is the electromotive force (EMF) constant, k E = k T ; R I is the resistance; V is the voltage supply; and ω is the angular velocity.
Now we proceed and transform all the above equations into a first order linear system as:
θ ˙ 1 = [ 0 + ω 1 + 0 + 0 + 0 + 0 + 0 ] + [ 0 + 0 + 0 + 0 + 0 ]
ω ˙ 1 = [ 0 + ( k β 1 + k E 1 k T 1 R I 1 ) ω 1 J 1 + 0 + 0 + 0 + 0 + 0 ] + [ M I C E J 1 + k T 1 V 1 R I 1 J 1 + 0 + M C J 1 + 0 ]
θ ˙ 2 = [ 0 + 0 + 0 + ω 2 + 0 + 0 + 0 ] + [ 0 + 0 + 0 + 0 + 0 + ]
ω ˙ 2 = [ 0 + 0 + 0 + 0 + 0 + ( k β 2 + k E 2 k T 2 R I 2 ) ω 3 J 2 i + J 3 ω ˙ 3 J 2 i ] + [ 0 + 0 + η k T 2 V 2 R I 2 J 2 + η M C J 2 + M v 0 J 2 i ]
θ ˙ 3 = [ 0 + 0 + 0 + 0 + 0 + ω 3 + 0 ] + [ 0 + 0 + 0 + 0 + 0 ]
ω ˙ 3 = [ 0 + 0 + 0 + 0 + 0 + k β 3 ω 3 J 3 + 0 ] + [ 0 + 0 + 0 + 0 + M v 0 ]
ω ¨ 3 = [ 0 + 0 + 0 + ( k β 2 + k E 2 k T 2 R I 2 ) ω 2 J 3 i + 0 + ( k β 2 J 2 i 2 + k α k v ) ω 3 J 2 J 3 i 2 + ( k v + k α J 3 + k α J 2 i 2 ) ω ˙ 3 ] + [ 0 + 0 + k α η k T 2 V 2 R I 2 J 2 J 3 i + k α η M C J 2 J 3 i + k α M v 0 J 2 J 3 i 2 ]
If we put the space vector as x 0 = [ θ 1 ω 1 θ 2 ω 2 θ 2 ω 3 ω ˙ 3 ] , the input variables as u 0 = [ M I C E V 1 V 2 M C M v 0 ] for the torque on the combustion engine (ICE), the input voltage for motor EM1 and EM2, torque on clutch, and the initial air-drag load, a linear space state of the vehicle dynamics system can be expressed as:
x ˙ 0 = [ 0 1 0 0 0 0 0 0 ( k β 1 + k E 1 k T 1 R 1 ) J 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 ( k β 2 + k E 2 k T 2 R 2 ) J 2 i J 3 ω ˙ 3 J 2 i 0 0 0 0 0 1 0 0 0 0 0 0 k β 3 ω 3 J 3 0 0 0 0 ( k β 2 + k E 2 k T 2 R 2 ) J 3 i 0 ( k β 2 J 2 i 2 + k α k v ) J 2 J 3 i 2 ( k v + k α J 3 + k α J 2 i 2 ) ] x 0 + [ 0 0 0 0 0 1 J 1 k T 1 R 1 J 1 0 1 J 1 0 0 0 0 0 0 0 0 η k T 2 R 2 J 2 η J 2 1 J 2 i 0 0 0 0 0 0 0 0 0 1 0 0 k α η k T 2 R 2 J 2 J 3 i k α η J 2 J 3 i k α J 2 J 3 i 2 ] u 0
The linear first order state space model in Equation (24) can be used to form the MPC algorithms in the next part. The system in Equation (24) includes acceleration ω ˙ 3 and jerk ω ¨ 3 , which can be used to simulate and regulate the HEV driving comfortability.
When the HEV runs in a low speed of less than 40km/h, only the main motor EM1 is activated. Considering the inputs of M I C E = 0 , V 1 = 0 , and M C = 0 , and the state variables of θ 1 = 0 and ω 1 = 0 , then the above linear system can be simplified as:
{ x ˙ p = A p x p + B p u p y p = C p x p + D p u p A p = [ 0 1 0 0 0 k β 2 + k E 2 k T 2 R I 2 J 2 0 0 0 0 0 1 0 0 0 ( k β 3 k v ) J 3 ] ,   B p = [ 0 0 k T 2 R I 2 J 2 0 0 0 0 1 ] C p = [ 0 0 0 1 k θ 0 k θ i 0 ] ;   D p = [ 0 0 0 1 ]
where the states x p = [ θ 2   ω 2   θ 3   ω 3 ] , inputs u p = [ V 2 M v 0 ] , outputs y p = [ ω 3 T T o r q u e 3 ] . The output T T o r q u e 3 is the unmeasured torque at shaft 3. In this equation, k θ is the torsional rigidity, k θ = M T o r q u e φ = G J l , and φ is the twist angle, φ = θ 2 θ 3 i . G is the rigidity modulus. l is the shaft length. J is the lumped inertia moment, J = J 2 + J 3 .
When the HEV runs in a high speed greater than 40 km/h, the starter motor EM2 activates the combustion engine ICE while the friction clutch is still open, the state equations of the first part can be written as:
θ ˙ 1 = [ 0 + ω 1 ] + [ 0 + 0 ]
ω ˙ 1 = [ 0 ( k β 1 + k E 1 k T 1 R I 1 ) J 1 ω 1 ] + [ ς k T 1 R I 1 J 1 V 1 + 1 J 1 M I C E ]
where ς is the additional coefficient for starting motor EM2 as a compensation load for the starting period. The linear state space system in the first part is:
{ x ˙ e = A e x e + B e u e y e = C e x e + D e u e ,   with A s = [ 0 1 0 ( k β 1 + k T 1 2 R I 1 ) J 1 ] ;   B s = [ 0 0 ς k T 1 R I 1 J 1 1 J 1 ] ;   C s = [ 0 1 0 ( k β 1 + k E 1 k T 1 R 1 ) J 1 ] ;   D s = [ 0 0 0 0 ]
where x e = [ θ 1 ω 1 ] , u e = [ V 1 M I C E ] , and y e = [ ω 1 T T o r q u e 1 ] . The output T T o r q u e 1 is the unmeasured torque at shaft 1.

3. Model Predictive Control for the HEV

A MPC is an open loop, infinite horizon prediction, and optimization subject to dynamic constraints. The continuous first order linear space state equation in Equation (24) can be discretized into time intervals with a discrete k and k + 1 = k + Δt ( Δ t is the computer scanning speed or the time sampling interval). Now, the continuous time form in Equation (24) can be discretized into:
{ x t + 1 = A x t + B u t y t = C x t + D u t
which is subject to the states, inputs, outputs, and the input increase constraints
u t U ,   Δ u t = u t u t 1 Δ U ,   and   y t Y
A MPC calculates the open loop input and output prediction horizon. For the calculation simplicity, we assumed the input prediction length is always equal to the output prediction length, or N u = N y . The objective function of the MPC for the HEV is:
m i n U { Δ u 1 , , Δ u 1 + N u 1 } { J ( U , x ( t ) ) = k = 0 N y 1 [ ( y t + k | t r ) Q ( y t + k | t r ) + Δ u t + k R Δ u t + k ] }
By substituting x k + j | k = A k x ( k ) + j = 0 k 1 A j B u k + j 1 j , Equation (31) can be transformed as
V ( x ( k ) ) = 1 2 x ( k ) Y x ( k ) + m i n U { 1 2 U H U + x ( k ) F U }
subject to the linear matrices inequality (LMI), G U W + E x ( t ) , where the column vector U [ u k , , u k + N 1 ] s ,   s m N u is the optimization vector, H = H > 0 , and H, F, Y, G, W, E are obtained from Q, R and in (9) as only the optimizer vector U is needed, the term involving Y is usually removed from (10). The optimization problem (10) is a quadratic program (QP). The MPC optimizer will calculate the optimal input vector U { Δ u k , , Δ u k + N u 1 } subject to the hard constraints of the inputs, u k U , and u k + i [ u m a x m i n ] ; of the outputs y k Y , and y k + i | k [ y m a x m i n ] ; and of the input increments Δ u k + i [ Δ u m a x m i n ] . But only the first input increment, Δ u k , is taken into the implementation. Then, the optimizer will update the outputs and states variables with the new update input and repeat the calculation for the next time interval. Therefore, the MPC is also called as the receding time horizon control. A diagram control system for this NMPC is shown in Figure 3.
The MPC scheme for the HEV in Figure 3 calculates the real-time optimal control action, Δ u k , and feeds into the vehicle dynamic equations and updates the current states, inputs, and outputs. The updated states, inputs, and outputs will feedback and compare to the reference desired trajectory data for generating the next optimal control action, Δ u k , in the next interval.
When the system is non-linear and has a general derivative nonlinear form, it is calculated as:
X ˙ = f ( x , u )
where x is the state variables and u is the inputs. The non-linear equation in (33) can be approximated in a Taylor series at referenced positions of ( x r , u r ) for X ˙ r = f ( x r , u r ) , so that:
X ˙ f ( x r , u r ) + f x , r ( x x r ) + f u , r ( u u r )
where f x . r and f r . x are the Jacobian function calculating approximation of x and u , respectively, moving around the referenced positions ( x r , u r ) .
Substituting Equation (34) for X ˙ r = f ( x r , u r ) , we can obtain an approximation linear form in continuous time ( t ) :
X ˜ ˙ ( t ) = A ( t ) X ˜ ( t ) + B ( t ) u ˜ ( t )
The linearized system in Equation (35) can be used as the linear system in Equation (24) for the MPC calculation. However, the MPC real-time optimal control action Δ u k + i | k must be fed into the original non-linear system in Equation (33) for the updated states, outputs, and inputs.

4. The MPC with Softened Constraints for the HEV

The conventional MPC objective function in Equation (31) subject to the constraints in Equation (30) regarding states, outputs, inputs, and input increase may deal with so many hard constraints. The MPC optimizer may not find out a solution that satisfies all constraints. Thus, we considered to widen the MPC feasibility by converting some possible hard constraints from Equation (30) into softened constraints to increase the possibility of finding a solution. The new MPC scheme subject to the softened constraints has the following form:
m i n U { Δ u k , , Δ u k + N u 1 } { J ( U , x ( k ) = i = 0 N y 1 [ ( y k + i | k r k + i | k ) Q ( y k + i | k r k + i | k ) + Δ u k + i | k R Δ u k + i | k + ε i ( k ) Λ ε i ( k ) + 2 μ ε k + i | k ] }
subject to
[ 1 z i z i X + μ ε i I ] 0 { m i n j X j j x m a x 2 z i v e r t { χ u * ( . | k ) k + i | k ( x ( k ) ) } ,   i { 1 , , N }
where μ is assigned as big values as a weighting factor ( μ > 0 ), and ε i is the constraints penalty terms ( ε i 0 ) added into the MPC objective function. X and z i are the corresponding matrix of the hard constraints.
The new items in Equation (37) are softened constraints selected from hard constraints in u k U , and u k U , Δ u k + i [ Δ u m a x m i n ] , for i = 0 , 1 , , N u 1 ,
y k Y , and y k + i | k [ y m a x m i n ] , for i = 0 , 1 , , N y 1 , Δ u k = u k u k 1 Δ U , and Δ u k + i = 0 , for i N u , x k | k = x ( k ) , x k + i + 1 | k = A ( k ) x k + i | k + B ( k ) u k + i , u k + i | k = u k + i 1 | k + Δ u k + i | k , y k + i | k = C ( k ) x k + i | k , where, ε i ( k ) = [ ε y ; ε u ] , y k + i | k y m a x m i n , and u k + i | k u m a x m i n ; and Λ = Λ 0 is the additional penalty matrix (generally Λ > 0 and assign to small values). In this new MPC scheme, the penalty term of the softened constraints i = 0 N p [ ε k + i | k Λ ε k + i | k + 2 μ ε k + i | k ] is added into the objective function with the positive definite and symmetric matrix Λ . This term penalizes the violations of constraints and, where possible, the free constrained solution is returned.
This MPC calculates the new optimization vector U S = [ U ε ] and the new MPC computational algorithms are:
Ψ S ( x ( t ) ) = m i n U S { 1 2 U S H S U S + x ( t ) F S U S } ,
subject to G S U S W S + E S x ( k ) , where U S [ u k , u k + 1 , , u k + N p 1 , ε k , ε k + 1 , , ε k + N p ] is the new optimization input vector; H S = [ H 0 0 M ] ; F S = [ F μ ] ; and the matrices for inequality constraints H, F, G, W, and E are obtained from Equation (38),
G S = [ G 0 g S I 0 I ] with g S = [ 0 0 0 0 Z B 0 0 0 Z A B Z B 0 0 Z A N p 1 B Z A N p 2 B Z B ] ,
W S = [ W w S 0 ] with w S = [ z z ] , and E S = [ E e S 0 ] with e S = [ Z Z A Z A 2 Z A N p ] .
To illustrate the ability of this controller, we test the two MPC schemes in Equations (31) and (36) by the following simple example as considering the non-linear system shown below:
x ˙ 1 = 2 x 2 + u ( 1 + x 1 ) x ˙ 2 = 2 x 1 + u ( 1 3 x 2 )
It is assumed that the system in Equation (39) is subjected to the hard state and input constraints x m i n = [ 1 1 ] and 2 u 2 . The linearized approximation of this system in (35) is: x ˙ = A x + B u , in which A = [ 0 2 2 0 ] and B = [ 1 1 ] . The weighting matrices are chosen as Q = [ 1 0 0 1 ]   and   R = 1 . The weighting matrices for softened constraints are chosen as A = [ 0 2 2 0 ] and B = [ 1 1 ] . It is assumed that the system is starting form an initial state position, x 0 = [ 0.72 0.35 ] . Figure 4 shows the performance of two NMPC schemes: this initial state position, x0, does not lead to any violation of states and input ( x m i n = [ 1 1 ] and 2 u 2 ). In this x 0 , the solutions of the two control schemes are always available. We can see that the NMPC with a softened state approaches the asymptotic point faster than the hard constraints. It means that, if we somehow loosen some of the constraints, the optimizer can generate easier optimal inputs and the system will be more stable.
It is interesting to see in Figure 4 that both schemes have x 1   m i n H a r d = 0.8475 and x 1   m i n S o f t e n e d = 0.8483 , which almost reach the hard constraint of x m i n = [ 1 1 ] . These states still have not violated the state constraints but selecting other initial positions x 0 may lead to state and input violations.
If we select x 0 = [ 0.9 0.8 ] , this initial condition will lead to the violations of the state and the input constraints as x 1   m i n = 1.0441 and u m a x = 2.2303 . These violations will make the RMPC with hard constraints infeasible. Meanwhile, the RMPC scheme with softened constraints is still running well and it is still easy to find optimal input solutions, as shown in Figure 5. After a very short transitional period, the fully constrained solution is returned, or there is no more constrained violation.
The new MPC scheme with softened constraints for the HEV will be further analyzed and simulated in the next section.

5. The MPC with Softened Constraints for the HEV

5.1. The MPC for the HEV in Pure Electrical Drive

The main motor ME1 was used to run the HEV at a low speed. In this mode, the clutch is open. ICE and ME2 are off. We run the MPC in this mode with the discrete time interval of 0.05 s. ME1 has a maximum power of 35 kW, a maximum torque of 205 Nm, rigidity torque k θ = 1158 , inertia J 2 = 1 , constants k E 2 = k T 2 = 10 , inertia J 3 = 2 , gear ratio i = 2.34 , damping k β 2 = 0.5 and k β 3 = 12 , and resistance R I 2 = 5 .
Some softened constraints were converted, as input constraints for the DC voltage applied for the vehicle is | V 2 | 300   V , u ( t ) = + / 5   V . The output softened constraints are also set on the shaft with the shear strength for carbon steel of τ = 25 MPa or N/mm2. The output torque on the shaft 2 is constrained as | T | = τ π d 3 16 , with the diameter d = 0.05   m . Then, the torque softened constraint on shaft 2 is | T 2 | = 455   Nm .
The MPC parameters were set up with the predictive horizon of N u = N y = N p = 5 ; the weighting matrices are set at Q = [ 1 0 0 1 ] and R = [ 1 ] . The MPC performance with softened constraints is shown in Figure 6.
It is worth noting that the weighting matrix for output Q and input R can be varied according to the desired variation on outputs or inputs. If we want to limit the errors or keep the output variation in a small value, we have to pay for more input energy or increase the input variation. With this aim in mind, we increased Q and reduced R. It means that any small variation in output will lead to a big penalty amount adding to the MPC objective function. Figure 7 shows the MPC for the HEV performance with Q = 100 and R = 1.
As shown in Figure 6 and Figure 7, we set up softened constraints on the input voltage of | V 2 | 300   V . The MPC allows some input voltage violation at the starting time to ensure the controller stability and feasibility. Then, after a very short transitional period, the solution is returned without constraint violation. In these cases, the MPC with hard constraints becomes infeasible and unstable.

5.2. The MPC for the HEV in High Speed with ICE

When the HEV runs at a high speed, the starter/generator ME2 starts the ICE. Depending on the required output torque, the ICE alone, or the ICE and ME1, or all ICE, ME1 and ME2 will be running and together providing torque.
In this mode, we assumed that the vehicle runs at ω 3 = 2000 rpm, and the torque of the air drag resistance at this speed of M v 0 = 30 Nm. The parameters of the starter motor EM2 are as constants k E 2 = k T 2 = 15 , inertia J 1 = 1 , damping coefficient k β 1 = 0.5 , resistance R I 1 = 7 , compensation ς = 0.5 , and a discrete time of 0.05 s.
The softened constraints were imposed on the input voltage constraints for the starter of | V 1 | 48   V , u ( t ) = + / 5   V / i n t e r v a l , and the output constrained torque on shaft 1 of | T 2 | 455   Nm .
For the MPC parameters, we selected the predictive horizon length of N u = N y = N p = 5 and the weighting matrices Q = [ 10 0 0 10 ] and R = [ 1 0 0 1 ] . The MPC performance with starting EM2 is shown in Figure 8.
Figure 8 shows that the EM2 starts in 1 s and the ICE is fully ignited and runs in 2.3 s; the ICE speed reaches the setpoints of 2000 rpm and steadily runs at 6.2 kW, providing the output torque of 31 Nm.
In the next simulation, we run the EM2 and the ICE to track the speed desired setpoints and ignite the clutch engagement. It was assumed that the main motor EM1 runs at 1500 rpm and the starter EM2 starts the ICE and is engaged into the system. The clutch engagement must take place at ω 1 ω 2 or ω 1 = 1.05 ω 2 for driving comfortability and a low jerk. The ICE and ME2 must track on the EM1 speed at +5% offset. The full engagement is completed in 2.3 s, as shown in Figure 9.
In Figure 9, we see the ICE and ME2 tracking ME1 on the desired setpoints in 1.9 s. In the normal mode, at speed higher than 40 km/h, the starter ME2 ignites the ICE and is turned as a generator charging to battery. The main motor EM1 now is also turned off. Only the ICE propels the HEV.
In Figure 10, it can be seen that the ME2 is turned off and becomes the generator after igniting the ICE. The main motor EM1 is also turned off, and the ICE alone propels the HEV. The HEV reaches and tracks the desired speed setpoints after 3.5 s.
Finally, we compared the performances of the MPC with hard constraints and he MPC with softened constraints. We run the MPC with hard constraints in Equation (31) and the MPC with softened constraints in Equation (36) to track the desired speed setpoints, as shown in Figure 11.
Figure 11 shows that the MPC with hard constraints generates smaller inputs and, hence, it needs a longer time to track into the speed setpoint. The MPC with hard constraints reaches the speed setpoint after 4.5 s, while the MPC with softened constraints needs only 3.5 s to fully track into the speed setpoint.

6. Conclusions

In this study, we have presented the modelling of the HEV and the MPC algorithms for controlling the HEV. In the HEV modelling, we have included the system acceleration and jerk into the equations to investigate and compare the vehicle driving comfortability with different control parameters. The MPC scheme with softened constraints has proved its superiority over the MPC with only hard constraints. The control system now becomes more flexible, stable, and robust against model uncertainties, time variant, and constraint violations. The new MPC scheme can control the HEV with faster clutch engagement and lower jerk reduction. The MPC with softened constraint still stable and robust, while the MPC with only hard constraints becomes unstable and infeasible because of the constraint violations. In a following study, we will investigate the control of the HEV friction clutch for a smooth and fast engagement with high comfortability and low jerk, and apply these algorithms in the real HEV.

Author Contributions

Conceptualization, T.M.V.; methodology, T.M.V.; formal analysis, T.M.V. and R.M.; writing—original draft preparation, T.M.V. and R.M.; writing—review and editing, R.M. and J.C.; supervision, J.H. and M.P. All authors have read and agreed to the published version of the manuscript.

Funding

The results were obtained with the financial support of the Ministry of Education, Youth, and Sports in the Czech Republic and the European Union, in the framework of the project “Modular platform for autonomous chassis of specialized electric vehicles for freight and equipment transportation”, Reg. No. CZ.02.1.01/0.0/0.0/16_025/0007293.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available from the corresponding author upon request.

Acknowledgments

The authors would like to thank Technical University of Liberec for Student Grant Scheme, project number SGS-2021-3059.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Minh, V.T. Advanced Vehicle Dynamics; Universiti of Malaya Press: Kuala Lumpur, Malaysia, 2012; p. 265. [Google Scholar]
  2. Moustafa, A.; Minh, V.T. Modeling and simulation of dual clutch transmission and hybrid electric vehicles. In Proceedings of the 11th International Conference of DAAAM Baltic Industrial Engineering, Tallinn, Estonia, 20–22 April 2016; pp. 168–174. [Google Scholar]
  3. Minh, V.T.; Afzulpurkar, N.; Muhamad, W.W. Fault detection model-based controller for process systems. Asian J. Control 2011, 13, 382–397. [Google Scholar] [CrossRef]
  4. Mrochen, M.; Sawodny, O. Modeling and simulation of a hybrid dual-clutch transmission powertrain. IFAC-PapersOnLine 2018, 51, 886–891. [Google Scholar] [CrossRef]
  5. Takahashi, Y.; Hidaka, K. Model predictive control for hybrid electric vehicles with linear parameter-varying model. In Proceedings of the 2018, 18th International Conference on Control, Automation and Systems (ICCAS), PyeongChang, Korea, 17–20 October 2018; pp. 1501–1506. [Google Scholar]
  6. Hua, X.; Zhanga, X.; Tanga, X.; Lin, X. Model predictive control of hybrid electric vehicles for fuel economy, emission reductions, and inter-vehicle safety in car-following scenarios. Energy 2020, 196, 117101. [Google Scholar] [CrossRef]
  7. Chen, Y.; Feng, G.; Wu, S.; Tan, X. A new hybrid model predictive controller design for adaptive cruise of autonomous electric vehicles. J. Adv. Transp. 2021, 2021, 1–25. [Google Scholar] [CrossRef]
  8. Glos, J.; Solc, F.; Otava, L.; Vaclavek, P. Hybrid model predictive control for fully electric vehicle thermal management system optimal mode selection. In Proceedings of the IECON 2020 the 46th Annual Conference of the IEEE Industrial Electronics Society, Singapore, 18–21 October 2020; pp. 2036–2043. [Google Scholar]
  9. Massimo, D.M.; Joris, G.; Jan, S.; Goele, P. Real-time model predictive control for a parallel hybrid electric vehicle using outer approximation and semi-convex cut generation. In Proceedings of the IEEE 16th International Workshop on Advanced Motion Control, Kristiansand, Norway, 22 February–20 April 2020; pp. 198–203. [Google Scholar]
  10. Oncken, J.; Chen, B. Real-time model predictive powertrain control for a connected plug-in hybrid electric vehicle. IEEE Trans. Veh. Technol. 2020, 69, 8420–8432. [Google Scholar] [CrossRef]
  11. East, S.; Cannon, M. Scenario model predictive control for data-based energy management in plug-in hybrid electric vehicles. arXiv 2021, arXiv:2104.07343. [Google Scholar]
  12. Shi, D.; Wang, S.; Cai, Y.; Chen, L.; Yuan, C.; Yin, C. Model predictive control for nonlinear energy management of a power split hybrid electric vehicle. Intell. Autom. Soft Comput. 2020, 26, 27–39. [Google Scholar] [CrossRef]
  13. Zhang, F.; Hu, X.S.; Liu, T.L.; Xu, K.; Duan, Z.D.; Pang, H. Computationally efficient energy management for hybrid electric vehicles using model predictive control and vehicle-to-vehicle communication. IEEE Trans. Veh. Technol. 2020, 70, 237–250. [Google Scholar] [CrossRef]
  14. Torreglosa, J.P.; Garcia-Triviño, P.; Vera, D.; López-García, D.A. Analyzing the improvements of energy management systems for hybrid electric vehicles using a systematic literature review: How far are these controls from rule-based controls used in commercial vehicles? Appl. Sci. 2020, 10, 8744. [Google Scholar] [CrossRef]
  15. Liu, W.; Chen, G.; Knoll, A. Matrix inequalities based robust model predictive control for vehicle considering model uncertainties, external disturbances, and time-varying delay. Front. Neurorobot. 2021, 14, 617293. [Google Scholar] [CrossRef]
  16. Schmitt, L.; Keller, M.; Albin, T.; Abel, D. Real-time nonlinear model predictive control for the energy management of hybrid electric vehicles in a hierarchical framework. In Proceedings of the 2020 American Control Conference, Denver, CO, USA, 1–3 July 2020; pp. 1961–1967. [Google Scholar]
  17. Ibrahim, M.H.; Maghfiroh, H.; Nizam, M.; Apribowo, C.H.B. Plug-in hybrid electric vehicle mode selection strategy for full battery consumption and known road slope condition. AIP Conf. Proc. 2020, 2217, 030149. [Google Scholar]
  18. Germanà, R.; Liberati, F.; di Giorgio, A. Decentralized model predictive control of plug-in electric vehicles charging based on the alternating direction method of multipliers. In Proceedings of the 28th Mediterranean Conference on Control and Automation (MED), Saint-Raphaël, France, 15–18 September 2020; pp. 739–745. [Google Scholar]
  19. Qi, Y.; Xiang, C.; Wang, W.; Wen, B.; Ding, F. Model predictive coordinated control for dual mode power split hybrid electric vehicle. Int. J. Automot. Technol. 2018, 19, 345–358. [Google Scholar] [CrossRef]
  20. Wang, Y.; Cao, S.; Yang, H.; Zuo, Z.; Wang, L.; Luo, X. Model predictive longitudinal control for autonomous electric vehicles with tracking differentiator. Int. J. Syst. Sci. 2021, 52, 2564–2579. [Google Scholar] [CrossRef]
  21. Liu, H.; Yan, S.; Shen, Y.; Li, C.; Zhang, Y.; Hussain, F. Model predictive control system based on direct yaw moment control for 4WID self-steering agriculture vehicle. Int. J. Agric. Biol. Eng. 2021, 14, 175–181. [Google Scholar] [CrossRef]
  22. Minh, V.T.; Pumwa, J. Feasible path planning for autonomous vehicles. Math. Probl. Eng. 2014, 2014, 317494. [Google Scholar] [CrossRef]
  23. Minh, V.T.; Hashim, F.B.M.; Awang, M. Development of a real-time clutch transition strategy for a parallel hybrid electric vehicle. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2011, 226, 188–203. [Google Scholar] [CrossRef]
  24. Minh, V.T.; Hashim, F.B.M. Tracking setpoint robust model predictive control for input saturated and softened state constraints. Int. J. Control Autom. Syst. 2011, 9, 958–965. [Google Scholar] [CrossRef]
  25. Minh, V.T.; Afzulpurkar, N.; Muhamad, W.M.W. Fault detection and control of process systems. Math Probl. Eng. 2007, 2007, 080321. [Google Scholar] [CrossRef]
  26. Minh, V.T.; Moezzi, R.; Dhoska, K.; Pumwa, J. Model predictive control for autonomous vehicle tracking. Int. J. Innov. Technol. Interdiscip. Sci. 2021, 4, 560–603. [Google Scholar]
  27. Ovchinnikov, I.; Kovalenko, P. Predictive control model to simulate humanoid gait: Predictive control model to simulate humanoid gait. Int. J. Innov. Technol. Interdiscip. Sci. 2019, 1, 9–17. [Google Scholar]
  28. Pumwa, J. Time variant predictive control of autonomous vehicles: Time variant pedictive control of autonomous vehicles. Int. J. Innov. Technol. Interdiscip. Sci. 2019, 2, 62–77. [Google Scholar]
Figure 1. The 2021 Hyundai Parallel HEV.
Figure 1. The 2021 Hyundai Parallel HEV.
Applsci 11 10668 g001
Figure 2. Simplified structure of the parallel HEV.
Figure 2. Simplified structure of the parallel HEV.
Applsci 11 10668 g002
Figure 3. Diagram of the MPC system.
Figure 3. Diagram of the MPC system.
Applsci 11 10668 g003
Figure 4. Performance of the two NMPC schemes.
Figure 4. Performance of the two NMPC schemes.
Applsci 11 10668 g004
Figure 5. Softened constraint NMPC.
Figure 5. Softened constraint NMPC.
Applsci 11 10668 g005
Figure 6. The MPC for the HEV with N p = 5 and Q = R .
Figure 6. The MPC for the HEV with N p = 5 and Q = R .
Applsci 11 10668 g006
Figure 7. The MPC for the HEV with N p = 5 , Q = 100, and R = 1.
Figure 7. The MPC for the HEV with N p = 5 , Q = 100, and R = 1.
Applsci 11 10668 g007
Figure 8. The MPC for the HEV with ICE and ME2.
Figure 8. The MPC for the HEV with ICE and ME2.
Applsci 11 10668 g008
Figure 9. The MPC for the HEV with ICE and ME2.
Figure 9. The MPC for the HEV with ICE and ME2.
Applsci 11 10668 g009
Figure 10. The MPC for the HEV propelled by only ICE.
Figure 10. The MPC for the HEV propelled by only ICE.
Applsci 11 10668 g010
Figure 11. Comparison of two MPC performances.
Figure 11. Comparison of two MPC performances.
Applsci 11 10668 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vu, T.M.; Moezzi, R.; Cyrus, J.; Hlava, J.; Petru, M. Parallel Hybrid Electric Vehicle Modelling and Model Predictive Control. Appl. Sci. 2021, 11, 10668. https://doi.org/10.3390/app112210668

AMA Style

Vu TM, Moezzi R, Cyrus J, Hlava J, Petru M. Parallel Hybrid Electric Vehicle Modelling and Model Predictive Control. Applied Sciences. 2021; 11(22):10668. https://doi.org/10.3390/app112210668

Chicago/Turabian Style

Vu, Trieu Minh, Reza Moezzi, Jindrich Cyrus, Jaroslav Hlava, and Michal Petru. 2021. "Parallel Hybrid Electric Vehicle Modelling and Model Predictive Control" Applied Sciences 11, no. 22: 10668. https://doi.org/10.3390/app112210668

APA Style

Vu, T. M., Moezzi, R., Cyrus, J., Hlava, J., & Petru, M. (2021). Parallel Hybrid Electric Vehicle Modelling and Model Predictive Control. Applied Sciences, 11(22), 10668. https://doi.org/10.3390/app112210668

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop