A Causal and Real-Time Capable Power Management Algorithm for Off-Highway Hybrid Propulsion Systems

Hybrid propulsion systems allow for a reduction of fuel consumption and pollutant emissions of future off-highway applications. A challenging aspect of a hybridization is the larger number of system components that further increases both the complexity and the diversification of such systems. Hence, beside a standardization on the hardware side for off-highway systems, a high flexibility and modularity of the control schemes is required to employ them in as many different applications as possible. In this paper, a causal optimization-based power management algorithm is introduced to control the power split between engine and electric machine in a hybrid powertrain. The algorithm optimizes the power split to achieve the maximum power supply efficiency and, thereby, considers the energy cost for maintaining the battery charge. Furthermore, the power management provides an optional function to control the battery state of charge in such a way that a target value is attained. In a simulation case study, the potential and the benefits of the proposed power management for the hybrid powertrain—aiming at a reduction of the fuel consumption of a DMU (diesel multiple unit train) operated on a representative track—will be shown.


Motivation
In the automotive industry, hybrid electric and full electric vehicles play a significant role in efforts to meet future legislated emission targets.Beside stricter limits for pollutant emissions of nitrogen oxides and particulates, European regulations also aim at reducing the average CO 2 emissions of passenger cars to 95 g/km [1].For off-highway applications-such as diesel multiple unit trains (DMUs) or marine applications-so far no legislative regulation to reduce CO 2 emissions has been implemented.Nevertheless, lowering fuel consumption and thus CO 2 emissions is the major goal in developing future off-highway propulsion systems.This is comprehensible considering that the fuel share of total life cycle costs for off-highway applications can be up to over 90% and, additionally, fuel prices are very likely to increase in the mid-and long-term future [2].Hybrid powertrain technology can help to significantly lower fuel consumption of rail and marine applications [3][4][5][6].Compared to conventional non-hybrid powertrains, additional operational functionalities of hybrids like recuperation of kinetic energy, engine operation point shifting or electric boosting improve the system performance and reduce pollutant emissions [7]. Figure 1 shows an overview of the functionalities of hybrid powertrains and their impact on emissions, fuel consumption and system performance.an overview of the functionalities of hybrid powertrains and their impact on emissions, fuel consumption and system performance.On the other hand, electrification of propulsion systems further increases the number of system components and, therefore, the complexity of already highly sophisticated off-highway architectures consisting of multiple engines, generators, gearboxes, etc.This leads to an even higher diversification of system configurations.To cope with this variety in series production, it is compulsory to standardize components of hybrid propulsion systems and to use them for various applications.Besides the hardware standardization, it is also necessary to harmonize software functions and control algorithms in order to apply them to as many different system configurations and architectures as possible.Therefore, a power management algorithm (PMA) is required to control the power split between combustion engine and electric motor in a hybrid propulsion system.The PMA has to offer the flexibility and modularity to be applicable to a vast variety of systems, including multi-engine system topologies.

Classification of Power Management Algorithms
In the past, a wide range of different academic approaches of PMAs to optimize the power split in hybrid propulsion systems for on-and non-road applications have been published in the literature.In order to classify multiple PMAs by the basis of their data, one can differentiate between causal and non-causal strategies [8,9].In causal strategies, the decisions made by the PMAs are based upon present and past state variables.On the contrary, in a non-causal approach, knowledge of the whole drive cycle is processed.Considering data of the complete drive cycle enables an optimization of the fuel consumption for a global optimum [9,10].For DMUs, where a detailed knowledge of the track and timetable is available, a non-causal PMA in combination with a predictive strategy is feasible.From a functional perspective, heuristic, often called rule-based or sub-optimal, and optimizationbased, or also known as optimal, algorithms can be distinguished [9,[11][12][13].Heuristic strategies include deterministic rule-based and fuzzy logic control algorithms [14][15][16][17].Those approaches are state of the art in most prototype and production hybrid vehicles [9].They are real-time capable and are processed online.The optimization-based or optimal PMAs are mainly processed offline.Those strategies are normally non-causal and, therefore, process information about the whole drive cycle.The optimal PMAs are divided into numerical methods such as dynamic programming (DP) and analytical approaches, like Pontryagins minimum principle [11,18,19].DP enables to find a global optimum but is not real-time capable and, therefore, frequently used as benchmark PMA.Other optimization-based PMAs are the so called equivalent consumption minimization strategies (ECMS) [9,12].The approach associates the stored electric energy to a future increase or decrease of fuel consumption.Therefore, usually an equivalence factor is imposed to convert battery into fuel power, and based on that into an equivalent fuel consumption.In the meantime, many different ECMS On the other hand, electrification of propulsion systems further increases the number of system components and, therefore, the complexity of already highly sophisticated off-highway architectures consisting of multiple engines, generators, gearboxes, etc.This leads to an even higher diversification of system configurations.To cope with this variety in series production, it is compulsory to standardize components of hybrid propulsion systems and to use them for various applications.Besides the hardware standardization, it is also necessary to harmonize software functions and control algorithms in order to apply them to as many different system configurations and architectures as possible.Therefore, a power management algorithm (PMA) is required to control the power split between combustion engine and electric motor in a hybrid propulsion system.The PMA has to offer the flexibility and modularity to be applicable to a vast variety of systems, including multi-engine system topologies.

Classification of Power Management Algorithms
In the past, a wide range of different academic approaches of PMAs to optimize the power split in hybrid propulsion systems for on-and non-road applications have been published in the literature.In order to classify multiple PMAs by the basis of their data, one can differentiate between causal and non-causal strategies [8,9].In causal strategies, the decisions made by the PMAs are based upon present and past state variables.On the contrary, in a non-causal approach, knowledge of the whole drive cycle is processed.Considering data of the complete drive cycle enables an optimization of the fuel consumption for a global optimum [9,10].For DMUs, where a detailed knowledge of the track and timetable is available, a non-causal PMA in combination with a predictive strategy is feasible.From a functional perspective, heuristic, often called rule-based or sub-optimal, and optimization-based, or also known as optimal, algorithms can be distinguished [9,[11][12][13].Heuristic strategies include deterministic rule-based and fuzzy logic control algorithms [14][15][16][17].Those approaches are state of the art in most prototype and production hybrid vehicles [9].They are real-time capable and are processed online.The optimization-based or optimal PMAs are mainly processed offline.Those strategies are normally non-causal and, therefore, process information about the whole drive cycle.The optimal PMAs are divided into numerical methods such as dynamic programming (DP) and analytical approaches, like Pontryagins minimum principle [11,18,19].DP enables to find a global optimum but is not real-time capable and, therefore, frequently used as benchmark PMA.Other optimization-based PMAs are the so called equivalent consumption minimization strategies (ECMS) [9,12].The approach associates the stored electric energy to a future increase or decrease of fuel consumption.Therefore, usually an equivalence factor is imposed to convert battery into fuel power, and based on that into an equivalent fuel consumption.In the meantime, many different ECMS approaches, primarily for automotive applications, were introduced.They use rule-or map-based routines as well as online adaption to determine the equivalent factor and optimize power split [20][21][22][23][24].

Requirements for the Off-Highway Application Power Management
In this paper a new power management algorithm (PMA) which controls the power split between combustion engine and electric motor in a hybrid propulsion system is presented.It is an online optimization-based causal strategy.Only present and past state variables are used, and no knowledge of the whole drive cycle is necessary.The PMA is developed to be applied to various applications and system configurations, including systems with multiple engines and a variety of electrical components.It is suited for use in hardware-in-the-loop (HiL) experiments as well as field testing of systems for off-highway applications.Furthermore it has to fulfill, amongst others, the following requirements:

•
Modular structure   It is important to point out that some of these requirements significantly differ from demands for PMAs used in automotive applications.Especially the necessity of a high modularity and the flexibility of the algorithm for various applications, topologies and system components are challenging concerning algorithm development and code implementation.
The paper comprises a simulation case study where the PMA is exemplarily applied to control a hybrid propulsion system for a DMU application.In Section 2, boundary conditions for the simulation case study are presented, including the specification of the system hardware in Section 2.1, as well as a short description of the drive strategies in Section 2.2, which are used for the simulations.Afterwards, in Section 3, a brief overview over the simulation model approach is given along with examples for a model validation.For the model validation, measured data from a hardware-in-the-loop (HiL) test rig is used.It represents a realistic prototype of a hybrid propulsion system driving a virtual vehicle model of the DMU.In Section 4, the PMA is described comprehensively.At first, the general problem for the control of hybrid systems is defined.Section 4.2 points out the possible operation modes of parallel hybrid powertrains which are relevant for the PMA.Thereafter, Section 4.3 presents the different functions and their interaction with the PMA.The results of the simulation case study are finally discussed in Section 5.

Hybrid Diesel Multiple Unit Train (DMU) System Specification
The PMA introduced in this paper has been developed for a Siemens Desiro VT642 DMU dedicated for regional passenger transportation.In a simulation case study, the vehicle is equipped with two identical parallel hybrid power units (PU), which form the propulsion system.Each PU consists of a diesel internal combustion engine (ICE) and an electric motor/generator unit (MG), which is placed on the engine output shaft center, see Figure 2. The angular velocities of ICE and MG can be decoupled by a clutch (CL).The gearbox (GR) is a six-speed automatic transmission with a hydrodynamic torque-converter on its input shaft.As energy storage unit, an electrochemical lithium-ion accumulator (BAT) is used in each PU.Furthermore, Figure 2 shows the arrangement of the auxiliary loads, comprising electric auxiliaries (AUX) and mechanical power take-offs (PTO) at the ICE, like the hydraulic fan.Table 1 shows a more detailed specification of the vehicle and system components.In the simulation case study, a system configuration is analyzed, where two identical PUs operate synchronously.It means, for example, that the power demand for the ICE calculated by the PMA is identical in both PUs.It is possible and, moreover, from an efficiency perspective rather beneficial to consider two PUs that are logically linked but operated individually.Such an operation strategy, however, is not state-of-the-art yet because a high voltage link between the electric circuits of both PUs would be required.Figure 3 outlines a system configuration with two PUs, which are linked on their electric circuit and share a joint battery.The second combustion engine is masked, which implies that such systems enable to restore one ICE.This is possible because-compared to a conventional non-hybrid system-the hybridization increases the system power considerably.Thus, considering track profiles and mission scenarios with a relatively low traction power demand, it is feasible to remove one of the ICEs and compensate the loss of traction power with the MGs.It becomes clear that even quite conventional propulsion system configurations of DMUs with two PUs can differ significantly by their design and operation, not to mention systems layouts with more than two PUs.This shows the necessity of a highly flexible PMA approach for off-highway hybrid systems.In the simulation case study, a system configuration is analyzed, where two identical PUs operate synchronously.It means, for example, that the power demand for the ICE calculated by the PMA is identical in both PUs.It is possible and, moreover, from an efficiency perspective rather beneficial to consider two PUs that are logically linked but operated individually.Such an operation strategy, however, is not state-of-the-art yet because a high voltage link between the electric circuits of both PUs would be required.Figure 3 outlines a system configuration with two PUs, which are linked on their electric circuit and share a joint battery.The second combustion engine is masked, which implies that such systems enable to restore one ICE.This is possible because-compared to a conventional non-hybrid system-the hybridization increases the system power considerably.Thus, considering track profiles and mission scenarios with a relatively low traction power demand, it is feasible to remove one of the ICEs and compensate the loss of traction power with the MGs.It becomes clear that even quite conventional propulsion system configurations of DMUs with two PUs can differ significantly by their design and operation, not to mention systems layouts with more than two PUs.This shows the necessity of a highly flexible PMA approach for off-highway hybrid systems.In the simulation case study, a system configuration is analyzed, where two identical PUs operate synchronously.It means, for example, that the power demand for the ICE calculated by the PMA is identical in both PUs.It is possible and, moreover, from an efficiency perspective rather beneficial to consider two PUs that are logically linked but operated individually.Such an operation strategy, however, is not state-of-the-art yet because a high voltage link between the electric circuits of both PUs would be required.Figure 3 outlines a system configuration with two PUs, which are linked on their electric circuit and share a joint battery.The second combustion engine is masked, which implies that such systems enable to restore one ICE.This is possible because-compared to a conventional non-hybrid system-the hybridization increases the system power considerably.Thus, considering track profiles and mission scenarios with a relatively low traction power demand, it is feasible to remove one of the ICEs and compensate the loss of traction power with the MGs.It becomes clear that even quite conventional propulsion system configurations of DMUs with two PUs can differ significantly by their design and operation, not to mention systems layouts with more than two PUs.This shows the necessity of a highly flexible PMA approach for off-highway hybrid systems.

Drive Strategy
Optimizations of the drive strategy to operate a DMU on a given track are the object of research in many publications.Often dynamic programming (DP) algorithms are used to optimize the velocity Energies 2017, 10, 10 5 of 17 trajectory of hybrid vehicles [19,25].Due to the fact that DP requires a high computational power, they are not real-time capable.Recently, approaches have been proposed where the DP algorithm is outsourced and computed on a cloud computing platform.The optimized velocity profile is sent back to the vehicle afterwards [26].Such approaches have significant potential to reduce fuel consumption.
However, for the simulation case study in this paper, a more simple drive strategy is used.Initially, a mission specification for the track, on which the DMU is operated, is defined.It specifies the driving distance, drive and stop durations as well as velocity limits of each track sequence.Based on this information, a driver simulation model, implemented in Simulink and Stateflow, computes the requested velocity profile.The objective of the driving strategy is the reduction of the DMU's energy consumption traveling on the tracks considering all the specified constraints (driving duration, velocity limits etc.).For this purpose, the algorithm aims at minimizing the cruising speed, considering the quadratic influence of velocity on the required traction force, see Equation (1).To operate the vehicle at a cruising speed as low as possible, the driving strategy asks for maximum acceleration and deceleration in the particular parts of each track sequence.Acceleration and deceleration demands are restricted by the performance capabilities of the PU, vehicle limitations or even operational boundaries, like a deceleration limit of −1 m/s 2 due to comfort considerations.In DMUs with conventional PUs, such a drive strategy, where the energy consumption is minimized, is considerably simple and leads to good results.Due to its simplicity it is fairly realistic that a human DMU operator in real world driving can easily track such a velocity profile.
For a hybrid system with energy storage unit, the possibility of a recovery of kinetic energy arises.For the hybrid drive strategy, hence, the algorithm of the driver model retards the deceleration phase.The strategy demands the vehicle to decelerate with a constant brake power, equivalent to the maximum mechanic brake power of the MG.This leads to slightly higher cruising speeds and, thus, energy consumption, but also significantly increases the amount of recuperated energy.Therefore, this drive strategy is used for hybrid operation mode.

System Simulation Model and Validation
The virtual vehicle and propulsion system are built up as a dynamic simulation model in the simulation environment GT-SUITE.GT-SUITE library blocks are used to model the physical components as well as to implement the basic control functions belonging to the components.Physical components and functions are combined to so-called compounds and archived in a compound library.The compounds possess standardized interfaces to enable a high degree of modularity, which allows to build simulation models of various system configurations for diverse applications.Each component provides the actual values of its physical state variables as well as their maximum and minimum values for each of them, as model outputs.More complex control algorithms, like the PMA, are realized in MATLAB/Simulink and Stateflow.The Simulink models can be compiled to dll-code and implemented in the GT-SUITE environment.The simulation models are suited for performance simulations, controls development, system parameter and configuration optimizations as well as hardware-in-the-loop testing.

Modelling System Components
This section gives a brief overview over the physical modeling theory and the equations, which the component models are based on.Furthermore, the predefined objects of the GT-SUITE library, which are used for the particular component models, are presented.Additionally, the most important model parameters are specified.

Vehicle Model
The dynamic vehicle simulation model uses the "GT VehicleBody", which is based on the equation of motion.The traction force F Tr is equal to the sum of all drive resistance forces: with track inclination ∝, air density ρ, vehicle velocity v and the reduced mass moment of inertia θ red .The required input traction torque M Drive to the vehicle differential is calculated by: Table 2 gives an overview over the most important vehicle model parameters.

Engine Model
The engine simulation model used in this case study is based on the "GT EngineState" object.It computes fuel consumption using a static brake specific fuel consumption (bs f c) map, see Figure 4a.The fuel consumption FC ICE is calculated by: The angular velocity ω ICE of the ICE in rad sec (or n ICE in rpm) follows from .
with the engine torque M ICE , the mass moment of inertia θ ICE of the ICE and the mechanical auxiliary torque M PTO .The mechanical auxiliaries or PTOs (power-take-offs) are directly linked to the engine output shaft.In the real vehicle, M PTO is a function of both angular velocity and PTO load, which depends on coolant or boost temperatures etc.In the simulation model, M PTO or P PTO is a function of the engine angular velocity, whereas the load is set constant, see Figure 4b.
The dynamic vehicle simulation model uses the "GT VehicleBody", which is based on the equation of motion.The traction force is equal to the sum of all drive resistance forces: with track inclination ∝, air density , vehicle velocity and the reduced mass moment of inertia .The required input traction torque to the vehicle differential is calculated by: Table 2 gives an overview over the most important vehicle model parameters.The engine simulation model used in this case study is based on the "GT EngineState" object.It computes fuel consumption using a static brake specific fuel consumption ( ) map, see Figure 4a.The fuel consumption is calculated by: The angular velocity of the ICE in (or in ) follows from with the engine torque , the mass moment of inertia of the ICE and the mechanical auxiliary torque .The mechanical auxiliaries or PTOs (power-take-offs) are directly linked to the engine output shaft.In the real vehicle, is a function of both angular velocity and PTO load, which depends on coolant or boost temperatures etc.In the simulation model, or is a function of the engine angular velocity, whereas the load is set constant, see Figure 4b

Gearbox Model
The gearbox is a six-speed transmission with a hydrodynamic torque converter (HTC) and transmits the input torque M GR to the vehicle differential input shaft.The model consists of a continuously-variable gear set, a HTC including a model of the torque converter clutch, and the controls module.It manages the clutch actuation and gear shifts.The controls are implemented in Simulink and integrated into the GT-SUITE model as dll-file.The required gear is computed externally by one module of the power management.The gearbox output torque is calculated by: (valid at closed torque converter clutch) with the MG effective torque M MG and the gear ratio i GR .
For simplicity, the gearbox efficiency η GR is set constant.The HTC enables the railway vehicle to accelerate with high traction forces.In gears with a gear ratio of z > 1 it is "deactivated" and locked with a torque converter clutch.To calculate the turbine torque of the HTC, M Turb the impeller torque M Imp is multiplied by the torque ratio i HTC , which is taken from the torque ratio table in Figure 5a.The impeller torque is calculated using the capacity factor K HTC .

Motor/Generator Model
To determine the electro-mechanic losses in the MG model, the efficiency η MG is interpolated linearly from a map, see Figure 5b.Beside the losses of the MG itself, it also includes the electric efficiency of the electric converter.From the mechanic power P MG and the efficiency, the electric power in the generator mode becomes Pe Gen = η MG • P MG , whereas the motor mode is characterized by Pe Mot = P MG η MG .Generally, electric machines are capable of running on overload for a restricted period of time.The overload duration is thermally restricted.Therefore, the MG involves a simple heat model to simulate heating and cooling sequences: Energies 2017, 10, 10 7 of 17

Gearbox Model
The gearbox is a six-speed transmission with a hydrodynamic torque converter (HTC) and transmits the input torque to the vehicle differential input shaft.The model consists of a continuously-variable gear set, a HTC including a model of the torque converter clutch, and the controls module.It manages the clutch actuation and gear shifts.The controls are implemented in Simulink and integrated into the GT-SUITE model as dll-file.The required gear is computed externally by one module of the power management.The gearbox output torque is calculated by: (valid at closed torque converter clutch) with the MG effective torque and the gear ratio .For simplicity, the gearbox efficiency is set constant.The HTC enables the railway vehicle to accelerate with high traction forces.In gears with a gear ratio of z > 1 it is "deactivated" and locked with a torque converter clutch.To calculate the turbine torque of the HTC, the impeller torque is multiplied by the torque ratio , which is taken from the torque ratio table in Figure 5a.The impeller torque is calculated using the capacity factor .

Motor/Generator Model
To determine the electro-mechanic losses in the MG model, the efficiency is interpolated linearly from a map, see Figure 5b.Beside the losses of the MG itself, it also includes the electric efficiency of the electric converter.From the mechanic power and the efficiency, the electric power in the generator mode becomes = • , whereas the motor mode is characterized by = .Generally, electric machines are capable of running on overload for a restricted period of time.The overload duration is thermally restricted.Therefore, the MG involves a simple heat model to simulate heating and cooling sequences:  The product of mass m and specific heat c (here for iron) represents the time constant of the differential Equation ( 9) and depends on the MG size.The heat flows follow from: . .
The parameters m , α w , T Coolant and A for this simplified approach are set to match pre-defined heating-up and cool-down times as well as to reach thermal steady state at the rated continuous power of the MG.The internal control routine of the MG model deactivates the overload mode and reduces the maximum power as soon as a temperature threshold is exceeded.MG power in overload mode is increased by a factor of 1.5 compared to the continuous operation mode, see Figure 5b.

Battery Model
The battery model is based on the equivalent electric circuit shown in Figure 6a.The voltage of the battery U BAT follows from Kirchhoff's voltage law for the illustrated equivalent circuit: Energies 2017, 10, 10 8 of 17 The product of mass and specific heat (here for iron) represents the time constant of the differential Equation ( 9) and depends on the MG size.The heat flows follow from: The parameters , , and for this simplified approach are set to match predefined heating-up and cool-down times as well as to reach thermal steady state at the rated continuous power of the MG.The internal control routine of the MG model deactivates the overload mode and reduces the maximum power as soon as a temperature threshold is exceeded.MG power in overload mode is increased by a factor of 1.5 compared to the continuous operation mode, see Figure 5b.
The open-voltage circuit represents the equilibrium potential of the battery, which corresponds to the number of battery cells connected in series (13).The battery stack capacity follows from the number of parallel linked battery cells , with: The internal resistance of the battery in ( 12) is taken from 2-D lookup tables parametrized for each battery cell type.For charging and discharging events, different tables of ( ) as a function of battery state of charge are employed.From the power sum at the electric circuit, the battery current can be calculated.Beside the MG, a further power sink of the electric circuit is the electric auxiliary power : The integration of according to (15) leads to the battery .It is one of the most important input variables of the PMA of a hybrid propulsion system, because it represents the energetic status of the energy storage.It is often the reference to trigger mode switches between various system operation modes.
is the initial state of charge value of the battery: The open-voltage circuit U OC represents the equilibrium potential of the battery, which corresponds to the number n ser of battery cells connected in series (13).The battery stack capacity C BAT follows from the number of parallel linked battery cells n par , with: The internal resistance R i of the battery in ( 12) is taken from 2-D lookup tables parametrized for each battery cell type.For charging and discharging events, different tables of R i (SOC) as a function of battery state of charge SOC are employed.From the power sum at the electric circuit, the battery current I BAT can be calculated.Beside the MG, a further power sink of the electric circuit is the electric auxiliary power Pe Aux : The integration of I BAT according to (15) leads to the battery SOC.It is one of the most important input variables of the PMA of a hybrid propulsion system, because it represents the energetic status of the energy storage.It is often the reference to trigger mode switches between various system operation modes.SOC 0 is the initial state of charge value of the battery:

Model Validation
This section gives a brief overview of the model validation.The component models outlined in the previous section are used for open-and closed-loop component tests.The measured data of those tests contributes to a further improvement of the model accuracy.

Battery Model Validation
For model parameterization and validation purposes, measured data from open-loop battery tests was used.In the test procedure, a battery prototype on a battery test bench with a DC/DC power source was stressed with a battery current profile simulated beforehand.The current profile results from a simulation of a DMU with two hybrid PUs operating in regional passenger transportation.Battery state variables as voltage, current and SOC are measured.Afterwards, simulations of the battery were executed to optimize the internal resistance tables in the battery model for a minimum deviation of measured and simulated values for voltage and SOC.The optimization is based on Brent's algorithm [27].Figure 6b shows a good matching between the calculated and observed SOC.

Powertrain Model Validation
In more extensive test campaigns of the whole propulsion system, validation data for the PU model was recorded.For this purpose, a prototype of a MTU hybrid PU was installed on a system test bench.In hardware-in-the-loop tests (HiL tests), the powertrain on the test bed is coupled with the simulation model of a DMU running on a virtual track.The test bed automation system is linked via a real-time interface to the simulation model of the vehicle and driver.In such a HiL test environment, performance tests can be executed or software functions are validated etc. Figure 7a gives a schematic overview of the experimental set-up.In Figure 7b, the simulated and measured results for the accumulated fuel consumptions over the whole drive cycle-for conventional and hybrid operation-are compared.The stated results are normalized on the fuel consumption for conventional mode.For the measurements campaign and the model validation, a previous version of the PMA, comprising a rule-based approach, was used.Obviously, the results underline and validate the achieved high accuracy of the simulation models.

Model Validation
This section gives a brief overview of the model validation.The component models outlined in the previous section are used for open-and closed-loop component tests.The measured data of those tests contributes to a further improvement of the model accuracy.

Battery Model Validation
For model parameterization and validation purposes, measured data from open-loop battery tests was used.In the test procedure, a battery prototype on a battery test bench with a DC/DC power source was stressed with a battery current profile simulated beforehand.The current profile results from a simulation of a DMU with two hybrid PUs operating in regional passenger transportation.Battery state variables as voltage, current and are measured.Afterwards, simulations of the battery were executed to optimize the internal resistance tables in the battery model for a minimum deviation of measured and simulated values for voltage and .The optimization is based on Brent's algorithm [27].Figure 6b shows a good matching between the calculated and observed .

Powertrain Model Validation
In more extensive test campaigns of the whole propulsion system, validation data for the PU model was recorded.For this purpose, a prototype of a MTU hybrid PU was installed on a system test bench.In hardware-in-the-loop tests (HiL tests), the powertrain on the test bed is coupled with the simulation model of a DMU running on a virtual track.The test bed automation system is linked via a real-time interface to the simulation model of the vehicle and driver.In such a HiL test environment, performance tests can be executed or software functions are validated etc. Figure 7a gives a schematic overview of the experimental set-up.In Figure 7b, the simulated and measured results for the accumulated fuel consumptions over the whole drive cycle-for conventional and hybrid operation-are compared.The stated results are normalized on the fuel consumption for conventional mode.For the measurements campaign and the model validation, a previous version of the PMA, comprising a rule-based approach, was used.Obviously, the results underline and validate the achieved high accuracy of the simulation models.

Problem Definition
The PMA controls the power flow in the propulsion system and optimizes the power split between the ICE and the MG.Based on the optimal control problem [26], a suitable cost functional or performance index J is stated in (17).The cost functional J substitutes a corrected amount of fuel consumption FC Cor of the system in a time interval t ∈ [t 0 , t end ].The objective of the PMA is to minimize J and FC Cor over the given drive cycle.FC Cor represents the quantity to evaluate the performance of the hybrid system according to its energy consumption and efficiency: J is calculated by the integration of fuel mass flow .
m Fuel in kg sec , as a function of the control inputs u(t) and state variables x(t).The first term in (17) stands for a correction term to penalize the deviation of energy stored in the battery between t 0 and t end : Here, λ represents a fuel equivalent factor, which converts the deviation of SOC into a value of fuel consumption in kg % .In contrast to some ECMS control algorithms, where an online adaption of λ is part of the PMA, λ is set constant in this case.It is determined by the amount of fuel, which is needed to balance battery charge deviation, assuming that charging is performed while the vehicle is at standstill (halt) at the end of each drive cycle: E 100 BAT stands for the maximum amount of energy, which can be stored in the battery at SOC = 100%.H u is the lower heating value of diesel fuel, and η op chrg denotes the charging efficiency at a predefined ICE operation point.This operating point for charging at standstill (halt) is optimized offline for maximum charging efficiency and set as a parameter.

Hybrid Operation Modes
The configuration of the parallel hybrid system, see Figure 2, enables different operation modes.Table 3 outlines the various operation modes and assigns them to drive events.The drive events result from the requested traction power P Trac .This table shows the signs of ICE and MG power demands, P ICE and P MG , as well as the time derivative .SOC of the state of charge for the different drive modes.The signs take the auxiliary powers into account: P PTO < 0 at the ICE output shaft and Pe AUX > 0 at the electric circuit.In Table 3 the plus sign, +, stands for a power request higher than zero.The minus sign, -, represents a negative power request, whereas, 0, refers to a constant value.The double plus and minus signs indicate very high or low power requests, which for example arise at acceleration or braking events.
During acceleration, two modes are possible.At Pure ICE mode, the traction power is exclusively provided by the ICE.The MG runs as generator and supplies electric auxiliaries.In the Combined Mode, the traction power is split between ICE and MG.During cruising, Pure Electric driving is possible.Furthermore, the battery can be charged by increasing the ICE load.In the case of a negative traction power demand, the deceleration event, the MG serves as generator to recuperate as much energy as possible.The clutch between ICE and MG is opened, to operate the ICE at idle speed.At halt, the vehicle velocity is zero.In Normal Halt mode, the ICE powers the PTO and supplies power for the AUX.If necessary, the ICE load can be increased to charges the battery.Additionally, the ICE may be switched off in Start/Stop mode, and the AUX power is supplied by the BAT.At deceleration and halt events, the different operation modes are triggered according to pre-defined thresholds of the SOC and other state variables.In those events, the PMA is a heuristic rule-based approach, where no optimization is executed.During acceleration and cruising, the optimization routine for P ICE and P MG is active.

PMA Optimization Routine
Figure 8 shows a signal flow scheme of the PMA optimization routine including its most important input and output signals.Each block of the scheme represents one function of the routine, which is described in this section.In case of a DMU with more than one PU, P ICE and P MG represent the power demand from the particular power path.This means, for example, that if more than one engine is present, P ICE substitutes the sum of engine power demands.A further function block splits P ICE and distributes the particular share of power to the different engines.In case of PUs running synchronously, as assumed in this paper, power is equally split to the PUs' machines.In this section, for simplicity, the PMA refers to the ICE and MG, neglecting the fact that more than one of each component is involved.
Energies 2017, 10, 10 11 of 17 Mode, the traction power is split between ICE and MG.During cruising, Pure Electric driving is possible.Furthermore, the battery can be charged by increasing the ICE load.In the case of a negative traction power demand, the deceleration event, the MG serves as generator to recuperate as much energy as possible.The clutch between ICE and MG is opened, to operate the ICE at idle speed.At halt, the vehicle velocity is zero.In Normal Halt mode, the ICE powers the PTO and supplies power for the AUX.If necessary, the ICE load can be increased to charges the battery.Additionally, the ICE may be switched off in Start/Stop mode, and the AUX power is supplied by the BAT.At deceleration and halt events, the different operation modes are triggered according to pre-defined thresholds of the and other state variables.In those events, the PMA is a heuristic rule-based approach, where no optimization is executed.During acceleration and cruising, the optimization routine for and is active.

PMA Optimization Routine
Figure 8 shows a signal flow scheme of the PMA optimization routine including its most important input and output signals.Each block of the scheme represents one function of the routine, which is described in this section.In case of a DMU with more than one PU, and represent the power demand from the particular power path.This means, for example, that if more than one engine is present, substitutes the sum of engine power demands.A further function block splits and distributes the particular share of power to the different engines.In case of PUs running synchronously, as assumed in this paper, power is equally split to the PUs' machines.In this section, for simplicity, the PMA refers to the ICE and MG, neglecting the fact that more than one of each component is involved.

Power Vector Definition
From the traction power required and the power request of the mechanical auxiliaries , this function block calculates the propelling power demand .It always equals the sum of the power demands, plus .Moreover, the gearbox input power limitations , as well as further vehicle constraints, have to be considered.In a next functional step, the inappropriate operation modes for the current drive event are rejected.This is necessary because the battery always has to be in a valid state of charge range, ≤ ≤ .If, for instance, is above its maximum value, Charge Mode is prohibited.Based on the feasible modes and constraints regarding the ICE, MG and BAT the maximum and minimum values for the power split ratio are defined.From those extremes, the operation point vectors , and for the split ratios of the length k are determined.They follow from:

Power Vector Definition
From the traction power required and the power request of the mechanical auxiliaries P PTO , this function block calculates the propelling power demand P D .It always equals the sum of the power demands, P ICE plus P MG .Moreover, the gearbox input power limitations [P min GR , P max GR ] as well as further vehicle constraints, have to be considered.In a next functional step, the inappropriate operation modes for the current drive event are rejected.This is necessary because the battery SOC always has to be in a valid state of charge range, SOC min ≤ SOC ≤ SOC max .If, for instance, SOC is above its maximum value, Charge Mode is prohibited.Based on the feasible modes and constraints regarding the ICE, MG and BAT the maximum and minimum values for the power split ratio ψ are defined.From those extremes, the operation point vectors ψ, P ICE and P MG for the split ratios of the length k are determined.They follow from: subject to

Power Split Optimization
The optimization function block determines the control output signals for P ICE and P MG .Therefore, the algorithm calculates a power supply efficiency vector of length k, Ω = (Ω 1 . . .Ω k ), where each element of the vector assigns a power supply efficiency Ω i to the particular elements of the operation point vectors.The power supply efficiency correlates to the sum of requested powers, P D , P PTO and Pe AUX , divided by the sum of power consumption, P Fuel plus Pe BAT− (electric discharge power): with: It is differentiated between electric battery power for charging Pe BAT+ and discharging Pe BAT− .P k BAT for each vector element k is calculated by: hence: The values for the conversion efficiencies η ICE , η MG and η BAT are taken from efficiency maps (e.g., Figure 5b) or, in case of η ICE , are calculated based on lookup tables of bs f c(P ICE (k), ω ICE ).The variable (t) in ( 23) stands for an energy cost equivalent charge factor and is described more detailed in the next subsection.Finally, the control output signals for P ICE and P MG are determined by the maximum value Ω max of the power supply efficiency vector Ω: arg max{Ω(P ICE , P MG , P D , P PTO , Pe AUX )} → P ICE , P MG . (27)

Energy Cost Equivalent Charge Factor
The energy cost equivalent charge factor (t) accounts for the amount of fuel energy that was used in the time interval t ∈ [t 0 , t] to charge the battery.Its valid range is η min chrg ≤ (t) ≤ 1.If the battery is charged only by recuperation during braking, without any consumption of fuel, (t) has the value 1.But if the energy stored in the battery is provided by the ICE, (t) equals the medium efficiency of the charging procedure.Hence, for .SOC > 0, (t) is computed by: with: and: In the case of discharging, .
SOC ≤ 0, the PMA modifies (28) and holds (t) constant.It follows from ( 28) that (t) substitutes the ratio of energy E * BAT , stored in the battery, and the fuel energy needed to charge the battery E * ChrgFuel .Please note that P ChrgFuel represents the share of fuel consumed for charging.P ChrgFuel can be calculated from the electric charging power at the battery terminal and the product of conversion efficiencies for the ICE as well as MG, leading to: The index * implies that E * BAT represents the amount of energy in between the admissible boundaries of the battery state of charge, SOC min ≤ SOC ≤ SOC max .This makes sense because only this partition of the battery is used during operation.
From a mathematical point of few, (t) set in ( 23) is handled like an efficiency value, which increases the energy costs, or in other words, reduces the efficiency of the electrical path.In case of (t) being close to one, it means that most of the energy in the battery was provided by "free" sources like recuperation or even external sources, for instance, external charging in a depot.As a consequence, a high portion of P D is provided by the electric path.This results in ψ close to zero due to its higher conversion efficiency as compared to the ICE path.In applications, system configurations or even drive cycles, where not much of the consumed energy can be recuperated during the deceleration and only a relatively small BAT is used, Charge Mode is employed quite regularly to maintain SOC within its admissible interval.In this case, (t) falls below one and, thus, increases the electric energy cost.Due to its implementation, (t) adapts itself dependent on system parameters, the drive cycle and other operational boundaries.

SOC Control Function
There are system configurations or constraints for specific applications, where SOC does not only have to be within its bounds but have to match a specific value at the end of the drive cycle or even among single track sequences.This is the case if the additional installed electric power is intended to improve system performance or, for example, to replace one combustion engine, see Figure 3. Therefore, the PMA features a function block that determines an expected amount of energy that is recuperated during braking down from the current vehicle velocity.The computation addresses the driving resistance forces according to (1), estimates a brake trajectory in compliance with the drive strategy, see Section 2.2, and determines the recuperation duration t Rec = f (v).From t Rec and an assumed mean electric recuperation power Pe mean Rec , a recuperated SOC calc Rec is calculated: The factor F i Rec in (32) represents a correction factor and is adapted online for a better match of the predicted SOC calc Rec and the real SOC Rec .It is calculated by a comparably simple approach according to (33) but leads to relatively good results (see Section 5): subject to − 0.15 ≤ dF Rec ≤ 0.15. (34) The computation of F i Rec is processed at discrete states i, after each recuperation event, where the velocity is reduced from cruising speed to a halt.Based on SOC calc Rec and a target value for SOC Trg , a finite-state machine sets a trigger to initiate the Charge Mode.The PMA strives for equality of SOC and SOC Trg at the start of the next cycle sequence.Hence, the Charge Mode is activated until SOC ≥ (SOC Trg − SOC calc Rec ) is achieved.The SOC control function is a straightforward approach to control the state of charge of the battery BAT in order to attain a predefined target value.By a proper choice of the target value, it allows for implementing different operating strategies for hybrid systems.If SOC Trg is close to the upper SOC limit, a performance oriented strategy is applied.Here, it is ensured that at any time a high battery SOC is available for electrically boosted acceleration events.For a low SOC Trg , the PMA allows a wider range of possible operation modes and, moreover, to optimize a power split only from an efficiency point of view.Furthermore, this charge triggering function is intended to be the interface for further optional predictive algorithms, which process the knowledge of the whole drive cycle.

Case Study Results
In the simulation case study, the DMU specified in Section 2 with two identical PUs is employed.The vehicle parameters are shown in Table 2.The electric auxiliary power is set to Pe Aux = 30 kW for the vehicle.In the reference case, where the DMU runs on Pure ICE mode, the auxiliaries are supplied by a generator driven by the ICE.For evaluation purposes, the corrected fuel consumption FC Cor , calculated with ( 15)-( 17), is used.The DMU drives on a track which was generated from an evaluation of various DMUs' field data.This track is representative for the real world operation of such a vehicle.
The fuel consumption of the reference diesel-driven vehicle adds up to 32.8 L. In hybrid operation, the DMU consumes 25.8 L, which denotes a reduction of FC Cor by 21.3%.Figure 9 shows the vehicle velocity profile over time for the diesel and hybrid system as well as the altitude profile.The velocity profiles differ due to the modified drive strategies and their particular brake trajectories (Section 2.2).
Energies 2017, 10, 10 14 of 17 The computation of is processed at discrete states , after each recuperation event, where the velocity is reduced from cruising speed to a halt.Based on and a target value for , a finite-state machine sets a trigger to initiate the Charge Mode.The PMA strives for equality of and at the start of the next cycle sequence.Hence, the Charge Mode is activated until ≥ ( − ) is achieved.The control function is a straightforward approach to control the state of charge of the battery BAT in order to attain a predefined target value.By a proper choice of the target value, it allows for implementing different operating strategies for hybrid systems.If is close to the upper limit, a performance oriented strategy is applied.Here, it is ensured that at any time a high battery is available for electrically boosted acceleration events.For a low , the PMA allows a wider range of possible operation modes and, moreover, to optimize a power split only from an efficiency point of view.Furthermore, this charge triggering function is intended to be the interface for further optional predictive algorithms, which process the knowledge of the whole drive cycle.

Case Study Results
In the simulation case study, the DMU specified in Section 2 with two identical PUs is employed.The vehicle parameters are shown in Table 2.The electric auxiliary power is set to = 30 kW for the vehicle.In the reference case, where the DMU runs on Pure ICE mode, the auxiliaries are supplied by a generator driven by the ICE.For evaluation purposes, the corrected fuel consumption , calculated with ( 15)-( 17), is used.The DMU drives on a track which was generated from an evaluation of various DMUs' field data.This track is representative for the real world operation of such a vehicle.
The fuel consumption of the reference diesel-driven vehicle adds up to 32.8 L. In hybrid operation, the DMU consumes 25.8 L, which denotes a reduction of by 21.3%.Figure 9 shows the vehicle velocity profile over time for the diesel and hybrid system as well as the altitude profile.The velocity profiles differ due to the modified drive strategies and their particular brake trajectories (Section 2.2).In Figure 10a, the battery is shown.Due to the activated control function, the at the end of the drive cycle and also at the end of each cycle sequence are fairly close to the target value of 55%.The maximum deviation of the and its target value at halt events occurs after the seventh drive sequence and adds up to 0.59%.Considering the comparatively simple approach from Section 4.3.4,Figure 10b shows the energy cost-equivalent charge factor ( ) over cycle time.It becomes obvious that ( ) falls during cruising and raises throughout recuperation.This is the typical behavior and points out that Charge Mode is active quite often.The fact that ( ) is not below one during the whole cycle shows that a relatively high amount of electric charging power for the battery is gained by recuperation.As a result, Pure Electric operation is executed rather often during cruising, because electric energy is quite "cheap".Figure 11a shows the percentages of the In Figure 10a, the battery SOC is shown.Due to the activated SOC control function, the SOCs at the end of the drive cycle and also at the end of each cycle sequence are fairly close to the target value SOC Trg of 55%.The maximum deviation of the SOC and its target value at halt events occurs after the seventh drive sequence and adds up to 0.59%.Considering the comparatively simple approach from Section 4.3.4,Figure 10b shows the energy cost-equivalent charge factor (t) over cycle time.It becomes obvious that (t) falls during cruising and raises throughout recuperation.This is the typical behavior and points out that Charge Mode is active quite often.The fact that (t) is not below one during the whole cycle shows that a relatively high amount of electric charging power for the battery is gained by recuperation.As a result, Pure Electric operation is executed rather often during cruising, because electric energy is quite "cheap".Figure 11a shows the percentages of the individual operation modes during acceleration and cruising.The Combined Mode extensively occurs during vehicle acceleration.The reason is that a combined use of ICE and MG increases the vehicle performance and, therefore, results in lower cruising speeds.About 56% of the time, the Pure Electric mode is activated, whereas the ICE is running on idle speed decoupled from the MG by the clutch.This leads to the question, whether it is feasible to completely switch off the ICE during driving.The employment of Start/Stop mode at halt further reduces the fuel consumption by 2.5%, which leads to a FC Cor of 25.0 L. If Start/Stop mode is expanded to cruising events as well, the FC Cor is further lowered by 1.2%, see Figure 11b.

Conclusions and Outlook
In this contribution, an optimization-based causal PMA for hybrid propulsion systems is presented.It optimizes the power split between ICE and MG according to vectors of possible operation points.For each vector element, which represents one feasible operational power distribution, the corresponding element of a power supply efficiency vector is calculated.The power split that is related to the highest efficiency determines the power demand for ICE and MG as control output.The computation of the power supply efficiency uses an energy-cost-equivalent charge factor, which accounts for the energy consumption to charge the battery.It differentiates between regenerative charging during braking and "energy expensive" charging via the ICE.Moreover, as an optional feature, the power split can be controlled in such a way that the battery matches a predefined target value at the end of the drive cycle.In a case study, the presented PMA is employed for the control a hybrid propulsion system of a DMU consisting of two identical PUs.Thereby, the DMU was operated on a track which was generated from logging field data of comparable vehicles.In

Conclusions and Outlook
In this contribution, an optimization-based causal PMA for hybrid propulsion systems is presented.It optimizes the power split between ICE and MG according to vectors of possible operation points.For each vector element, which represents one feasible operational power distribution, the corresponding element of a power supply efficiency vector is calculated.The power split that is related to the highest efficiency determines the power demand for ICE and MG as control output.The computation of the power supply efficiency uses an energy-cost-equivalent charge factor, which accounts for the energy consumption to charge the battery.It differentiates between regenerative charging during braking and "energy expensive" charging via the ICE.Moreover, as an optional feature, the power split can be controlled in such a way that the battery matches a predefined target value at the end of the drive cycle.In a case study, the presented PMA is employed for the control a hybrid propulsion system of a DMU consisting of two identical PUs.Thereby, the DMU was operated on a track which was generated from logging field data of comparable vehicles.In

Conclusions and Outlook
In this contribution, an optimization-based causal PMA for hybrid propulsion systems is presented.It optimizes the power split between ICE and MG according to vectors of possible operation points.For each vector element, which represents one feasible operational power distribution, the corresponding element of a power supply efficiency vector is calculated.The power split that is related to the highest efficiency determines the power demand for ICE and MG as control output.The computation of the power supply efficiency uses an energy-cost-equivalent charge factor, which accounts for the energy consumption to charge the battery.It differentiates between regenerative charging during braking and "energy expensive" charging via the ICE.Moreover, as an optional feature, the power split can be controlled in such a way that the battery SOC matches a pre-defined target value at the end of the drive cycle.In a case study, the presented PMA is employed for the control a hybrid propulsion system of a DMU consisting of two identical PUs.Thereby, the DMU was operated on a track which was generated from logging field data of comparable vehicles.In comparison to the conventional diesel-driven system, a fuel consumption reduction of 21.3% was achieved.Employing Start/Stop at halt and, additionally, during Pure Electric driving, the fuel consumption was further reduced by 3.7%.It was stated that the PMA leads to high proportion of Pure Electric driving during cruising events.This has a significant influence on the engine load profile because especially areas of low engine loads, where engine efficiency is weak, disappear.This is even more true as soon as the engine is switched off during Pure Electric driving.Beside the shown fuel consumption benefits, the Pure Electric mode leads to shorter engine operation time and extends the time-before overhaul (TBO), which further reduces operational costs for the system operator.Another aspect which has to be considered is the influence of hybridization on exhaust emissions.DMUs' modern diesel engines are compliant to EPA Tier 4i emission standards.Those are commonly equipped with SCR (selective catalytic reduction) exhaust after-treatment systems.Investigations in the HiL test environment show that the changed engine operation in hybrid systems, compared to conventional ones, considerably influence exhaust after-treatment performance, hence real driving emissions [7].The implementation of Start/Stop during Pure Electric driving are one possibility to influence exhaust after-treatment temperatures and, for example, prevent the SCR-catalyst from cooling down too quickly during idling.
In future investigations, the introduced PMA will be applied to system configurations with multiple PUs, where the ICEs and MGs are operated individually.Furthermore, the PMA's suitability for other applications will be proven.To demonstrate the "real-world" capability of the PMA, it will be applied to a PU prototype of a DMU propulsion system in hardware-in-the-loop testing campaigns on a system test bench.In those experiments, the impact and the benefits of the proposed PMA and its different operation modes on fuel consumption and emissions will be investigated for different DMU propulsion system configurations.

Figure 2 .
Figure 2. Scheme of a parallel hybrid power unit.

Figure 3 .
Figure 3. Scheme of a parallel hybrid propulsion system for a diesel multiple unit train (DMU) consisting of two PUs and a combined electric circuit with the potential to replace one ICE.

Figure 2 .
Figure 2. Scheme of a parallel hybrid power unit.

Figure 2 .
Figure 2. Scheme of a parallel hybrid power unit.

Figure 3 .
Figure 3. Scheme of a parallel hybrid propulsion system for a diesel multiple unit train (DMU) consisting of two PUs and a combined electric circuit with the potential to replace one ICE.

Figure 3 .
Figure 3. Scheme of a parallel hybrid propulsion system for a diesel multiple unit train (DMU) consisting of two PUs and a combined electric circuit with the potential to replace one ICE. .

Figure 4 .
Figure 4. (a) Generic bsfc map for a diesel engine; (b) Mechanical auxiliary power over engine speed .

Figure 4 .
Figure 4. (a) Generic bsfc map for a diesel engine; (b) Mechanical auxiliary power P PTO over engine speed n ICE .

Figure 5 .
Figure 5. (a) Torque converter torque ratio and capacity factor over the speed ratio; (b) MG efficiency map.Figure 5. (a) Torque converter torque ratio i HTC and capacity factor K HTC over the speed ratio; (b) MG efficiency map.

Figure 5 .
Figure 5. (a) Torque converter torque ratio and capacity factor over the speed ratio; (b) MG efficiency map.Figure 5. (a) Torque converter torque ratio i HTC and capacity factor K HTC over the speed ratio; (b) MG efficiency map.

Figure 6 .
Figure 6.(a) Equivalent electric circuit for the battery; (b) Battery validation cycle: simulated and measured over time.

Figure 6 .
Figure 6.(a) Equivalent electric circuit for the battery; (b) Battery validation cycle: simulated SOC Sim and measured SOC Mes over time.

Figure 7 .
Figure 7. (a) Schematic overview over experimental HiL layout; (b) Measured and simulated normalized fuel consumption for conventional and hybrid mode.

Figure 7 .
Figure 7. (a) Schematic overview over experimental HiL layout; (b) Measured and simulated normalized fuel consumption for conventional and hybrid mode.

Figure 8 .
Figure 8. Signal flow scheme of the PMA optimization routine.

Figure 8 .
Figure 8. Signal flow scheme of the PMA optimization routine.

Figure 9 .
Figure 9. Altitude and velocity profiles for conventional and hybrid mode.

Figure 9 .
Figure 9. Altitude and velocity profiles for conventional and hybrid mode.

Energies 2017 ,
10, 10 15 of 17 individual operation modes during acceleration and cruising.The Combined Mode extensively occurs during vehicle acceleration.The reason is that a combined use of ICE and MG increases the vehicle performance and, therefore, results in lower cruising speeds.About 56% of the time, the Pure Electric mode is activated, whereas the ICE is running on idle speed decoupled from the MG by the clutch.This leads to the question, whether it is feasible to completely switch off the ICE during driving.The employment of Start/Stop mode at halt further reduces the fuel consumption by 2.5%, which leads to a of 25.0 L. If Start/Stop mode is expanded to cruising events as well, the is further lowered by 1.2%, see Figure 11b .

Figure 10 .Figure 11 .
Figure 10.(a) and its target value over cycle time; (b) Energy cost equivalent charge factor over cycle time.

Figure 10 .
Figure 10.(a) SOC and its target value over cycle time; (b) Energy cost equivalent charge factor over cycle time. .

Figure 10 .Figure 11 .
Figure 10.(a) and its target value over cycle time; (b) Energy cost equivalent charge factor over cycle time.

Figure 11 .
Figure 11.(a) Percental proportion of operation modes during acceleration and cruising (hybrid without Start/Stop); (b) Overview over percental reductions of corrected fuel consumption.

•
Suitability for different mission scenarios and velocity/load control modes (including driver-controlled traction torque demand and automatically controlled drive strategies)•Consideration of variable auxiliary or external load requests in the PMA.

Table 1 .
System specification of vehicle and components.

Table 1 .
System specification of vehicle and components.

Table 1 .
System specification of vehicle and components.

Table 3 .
Drive Events and Hybrid Operation Modes.