Next Article in Journal
Power and Frequency Control in the National Power System of the 370 MW Coal Fired Unit Superstructured with a Gas Turbine
Previous Article in Journal
Development of a Model for the Estimation of the Energy Consumption Associated with the Transportation of CO2 in Pipelines
Previous Article in Special Issue
Boundary Identification for Traction Energy Conservation Capability of Urban Rail Timetables: A Case Study of the Beijing Batong Line
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electric Powertrain Topology Analysis and Design for Heavy-Duty Trucks

Department of Mechanical Engineering, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
*
Author to whom correspondence should be addressed.
Submission received: 30 March 2020 / Revised: 1 May 2020 / Accepted: 5 May 2020 / Published: 12 May 2020

Abstract

:
Powertrain system design optimization is an unexplored territory for battery electric trucks, which only recently have been seen as a feasible solution for sustainable road transport. To investigate the potential of these vehicles, in this paper, a variety of new battery electric powertrain topologies for heavy-duty trucks is studied. Thereby, topological design considerations are analyzed related to having: (a) a central or distributed drive system (individually-driven wheels); (b) a single or a multi-speed gearbox; and finally, (c) a single or multiple electric machines. For reasons of comparison, each concurrent powertrain topology is optimized using a bilevel optimization framework, incorporating both powertrain components and control design. The results show that the combined choice of powertrain topology and number of gears in the gearbox can result in a 5.6% total-cost-of-ownership variation of the vehicle and can, significantly, influence the optimal sizing of the electric machine(s). The lowest total-cost-of-ownership is achieved by a distributed topology with two electric machines and two two-speed gearboxes. Furthermore, results show that the largest average reduction in total-cost-of-ownership is achieved by choosing a distributed drive over a central drive topology (−1.0%); followed by using a two-speed gearbox over a single speed (−0.6%); and lastly, by using two electric machines over using one for the central drive topologies (−0.3%).

Graphical Abstract

1. Introduction

Driven by stringent demands to further reduce CO2 emissions of road transport, an increasing number of truck manufacturers are announcing full-electric versions of their vehicles [1]. This has unavoidably led to a trend towards the development of battery electric trucks, for both city and highway applications. Among all the vehicles being announced and being sold, there is a large variation in the powertrain design: the electric machine used, the transmission technology, and the drivetrain topology. This large variation in the powertrain design (also referred to as architecture) raises the question of what is the best powertrain design. To answer this question, one has to understand and try to find the powertrain topology and component specifications that minimize the objective at hand, be that total-cost-of-ownership (TCO), energy efficiency, construction costs, or any other. This is equally relevant for customers, suppliers, and manufacturers.

1.1. System-Level Design

To find this optimal powertrain design, one must define and solve the problem by looking both at the powertrain’s hardware design and its control design due to the coupling between both [2]. Previous research has shown that investigating these problems separately leads to sub-optimal solutions [3]. For integrated plant and control design, common practice is to split the problem into four different design layers: (i) topology, (ii) component technology, (iii) component sizing, and (iv) powertrain controls [2]. This is sometimes referred to as a system-level design (SLD), which is graphically represented in Figure 1.
Considering all four design levels leads to a multi-level optimization problem, which requires a large number of function evaluations, which quickly increases when more design parameters are added. Due to this high computational burden, most powertrain design research has focused on one (e.g., control) or two design levels at maximum (e.g., sizing and control). In the context of battery electric trucks, most literature focuses on the feasibility [4,5,6,7,8] or performance evaluation [9] analyses. In addition, some recent research has focused on the type of technology used for the components, in particular focusing on: (i) the electric machine technology [10], (ii) battery technology [11,12], or (iii) the transmission technology [10,13,14,15]. A few publications considered optimizing the components’ sizing for a single or a few topologies [10,15,16], often combined with the control optimization. Only limited research seems to focus specifically on the variation of the powertrain topology for all-electric vehicles; with limited examples that can be found either focused on hybrid electric vehicles [17,18,19,20] (e.g., where the focus is limited to the type of hybrid powertrain: series, parallel or power split) or the research is limited to the light-duty vehicle domain [21,22,23,24,25]. Therefore, there is little knowledge on the influence of topological design choices for full battery electric powertrains for heavy-duty trucks, which is crucial in order to evaluate the full potential of these type of vehicles. Next to that, this knowledge is also highly valuable for vehicle manufacturers to support them in the choice of the most suitable powertrain topology for their vehicle(s).

1.2. Contribution and Outline of the Work

In this paper, the potential of novel electric powertrain topologies is investigated for heavy-duty trucks by analyzing the potential of various design options. In particular, the focus is on understanding the effects of the following three topological design choices: (a) a central or distributed drive, via the use of a final drive transmission; (b) incorporating a single or a multi-speed (discrete) gearbox, via the use of one or more gears; and finally, (c) the number of electric machines.
Although huge amounts of possible topologies can be constructed even with a small set of component elements, in this work only a relatively small set of relevant topologies is studied. For each of the feasible topologies, an integrated powertrain design optimization is performed. The performance metrics, i.e., total-cost-of-ownership and energy consumption, with constraints on drivability (e.g, acceleration, top speed, uphill driving), are compared in order to understand the influence of certain powertrain design choices on the topological level (a)–(c) explained above.
In the remainder of the paper, the powertrain design problem and the design method used are discussed in Section 2. The simulation model, including the sub-models, used in the design optimization is discussed in Section 3. The results of the design optimization are discussed in detail in Section 4, followed by the conclusions in Section 5, respectively.

2. Problem Description

For any given powertrain topology, as shown in Figure 1, finding an optimal design is typically a bilevel optimization problem, where both the powertrain components and the control algorithm are optimized. Throughout the rest of the paper, for most parameters, the domain to which the parameter belongs is mentioned. In the case of equal types of parameters and parameters for which this is self-evident, this is omitted.

2.1. System-Level Objective Function

The main design objective for commercial vehicles, which is often used, is the total-cost-of-ownership, C TCO R + (euro),
C TCO = C e + C v .
where C e R + (euro) are the electricity running costs of the vehicle and C v R + (euro) is the vehicle cost. Throughout the rest of this paper, for all cost parameters, indicated by the symbol C, C R + holds. The TCO represents the cost to operate the vehicle, which preferably should be as low as possible. The electricity costs are calculated as:
C e = 100 c e E v D y y .
Here, c e R + (euro/kWh) is the electricity cost per unit energy, E v R + (kWh/100 km) is the average energy consumption, D y (km/year) the yearly mileage of the vehicle, and y (year) the years of vehicle ownership. In the running costs, C e , only the costs related to the energy consumption of the vehicle are included. The costs related to maintenance, road taxes, driver wages, and loss of interest are left out of consideration. Therefore, also the running costs only include the cost of the electricity; hence, C e will also be referred to as the electricity cost. Especially, road taxes and driver wages are not dependent on the vehicles’ powertrain design. The average energy consumption, E v , is defined as energy consumed over a drive cycle starting at time, t 0 (s), and ending at final time, t end (s),
E v = 100 D c t o t end P s ( t ) d t .
Here, D c (km) is the distance traveled over the drive cycle and P s R (W) the internal battery power. The electricity cost, C e , depends on the size of the battery, electric machine(s), and gear ratio(s), since these values determine the limitations on the controlled inputs and states.
Consequently, the vehicle cost, C v , in the form of the depreciation cost of the vehicle over the lifetime, is modeled as a function of the costs of the powertrain components. The depreciation is the difference between the initial acquisition cost of the vehicle and the price for which the vehicle is sold to the second-hand market at the end of its lifetime. Linear depreciation is assumed with an economic life time, y el , at which the vehicle no longer has a resale value. This results in:
C v = y y el C pt + C 0 ,
with C pt (euro) the costs of the powertrain, which depend on the powertrain design, and C 0 (euro) a system base cost for all other remaining components and parts of the vehicle. Note that minimizing the TCO does not necessarily imply that both C e and C v are minimized. An example is when an extra electric machine is added to the powertrain. In this case, the vehicle cost C v will increase; however, the energy consumption and therewith the cost C e will decrease, resulting in a decrease in C TCO .

2.2. Topology Design Space

A topology refers to the layout in which the components in the powertrain are connected to each other. Depending on the set of components, degree of connection per component, and number of instances per component, a set of possible topologies, T p , can be generated. For the problem considered here, with the component set listed in Table 1, the total size of this design set is in the order of O ( 10 20 ) .
In this study, topologies are considered that are constructed with up to four electric machines EM k , two transmission units TR l , and the choice of whether to use a final drive transmission FD or not. A single battery pack BA is assumed to be connected to the electric machine(s) via integrated power electronics (not modeled separately). The powertrain built by these components is connected to one set of truck wheels that drive the vehicle. The subscript k { 1 , , n m } N refers to the index of the electric machine, with n m N (-) the total number of electric machines present in the powertrain build. The subscript l { 1 , , n t } N refers to the index of the transmission with n t N 0 (-) the total number of transmission units in the powertrain.
Using functional and application-specific constraints, this large set of possible topologies can be reduced to a smaller set of feasible topologies T f T p , which can be done using constraint programming [26]. In this work, a sub-set of the feasible topologies T f is chosen, based on three topological design choices combined with an overview from the most recent announcements from vehicle manufacturers. This set, in total eight separate topologies, is graphically shown in Figure 2.
Each of these topologies is analyzed with different types of transmission units. Only discrete multi-speed gearboxes will be considered here, with different numbers of gears in the gearbox. Herein, the combination of a chosen powertrain topology and chosen number of gears for the gearbox(es) is referred to as a powertrain configuration. The choice of eight topologies and four different gearbox options results in a unique set of powertrain configurations of (distributed + central) = 12 + 44 = 56 options. Due to identical options for two topologies and the fact that for one topology, only up to three speeds are considered, a total set of 44 powertrain configurations is studied. Each configuration in the set is indicated by T i . j l T f . The first sub-script, i, refers to the name of the topology,
i { D 1 , , D 3 , C 1 , , C 5 } ,
where D refers to distributed and C to central. Typically, different classifications of topologies can be made, with often the largest division being between central and distributed drive powertrains. Central indicates that the power to drive the vehicle is generated centrally and is divided over the wheels by using a final drive transmission, also called a differential. This is the lower set of five topologies in Figure 2. In the case of a distributed powertrain, the power is generated separately for each wheel. In this type of topology, the electric machines are most of the time located in or close to the driven wheels. Next, the subscript j l , defined by:
j l { 1 , 2 , 3 , 4 } ,
refers to the total number of gears of the discrete gearbox l, where the subscript l refers to the index of the gearbox. In the case of topologies with multiple gearboxes, each subscript j l refers to the total number of gears for each gearbox, e.g., T C 4.14 refers to powertrain configuration C 4 where TR1 is a one-speed and TR2 is a four-speed gearbox. Whereas the technology of the transmission is varied, the component technologies of the electric machine(s) with power electronics and the battery pack are considered to be fixed.

2.3. System Design Problem

For each topology, the system design problem requires finding the optimal sizing of the battery, electric machine(s), and gear ratios. The multi-level optimization problem aims to find the optimal powertrain design parameters, x p , and control inputs, x c ( t ) , respectively, by solving:
min x p , x c ( t ) C TCO ( x p , x c ( t ) ) = min x p , x c ( t ) C e ( x p , x c ( t ) ) + C v ( x p ) s . t . h p ( x p ) = 0 , g p ( x p ) 0 , h c ( x p , x c ( t ) ) = 0 , g c ( x p , x c ( t ) ) 0 .
Here, g and h are the inequality and equality constraints, respectively. Due to the unidirectional coupling between the electricity and the acquisition cost of the vehicle, the problem is solved using a nested optimization approach [3], where the control problem can be seen as an additional optimization problem in the set of constraints. By doing that, (7) becomes:
Input : topology T i , j l ,   drive   cycle ,   vehicle parameters Outer loop   ( component sizing problem ) : min x p C TCO ( x p | x c ( t ) ) = C e ( x p | x c ( t ) ) + C v ( x p ) , s . t . h p ( x p ) = 0 , g p ( x p ) 0 , Inner loop   ( control problem ) : min x c ( t ) C e ( x c ( t ) | x p ) , h c ( x c ( t ) | x p ) = 0 , g c ( x c ( t ) | x p ) 0 .
Note that bold notation indicates variables that are either a vector or a matrix. The superscript ∗ refers to the optimal solution of the control problem solved. Through (8), it is implied that in the inner loop, the optimal control signals are found only for the powertrain component sizing parameters x p that satisfy the (in)equality constraints in the outer loop.

2.4. Outer Loop-Component Sizing Problem

The set of powertrain component sizing variables x p R n p depends on the specific vehicle topology that is optimized, but it generally consists of:
x p = { E b , P m , R } ,
where E b R + (kWh) is the total battery energy content and P m R + n m (W) is a vector containing the mechanical peak power of each electric machine, P m k P m (W). R R + n r contains the gear ratios of all gearboxes. R is subdivided into multiple vectors r l R (-), where each vector r l contains the gear ratios for one of the gearbox units present in the powertrain, which is described as:
r l = { r l ( 1 ) , , r l ( j l ) } .
Here, each entry is the gear ratio of one of the gears of gearbox unit l, with j l the total number of gears for this gearbox. To reduce the number of design variables, the value of the gear ratios in the case of a central drive topology includes the ratio of the final drive. Thus, the final drive ratio is r f = 1 . Furthermore, for these topologies, the total ratio between the wheels and the electric machine(s) is optimized. The total number of component sizing parameters, n p , which depends on the powertrain configuration being optimized, is:
n p = 1 + n m + n r ,
where n m is the total number of electric machines and n r N 0 is calculated by:
n r = l = 1 n t j l ,
and denotes the total number of gear ratios present in the powertrain, where n t is the total number of gearboxes.

Component Constraints

The set of constraints defined under the component sizing problem of (8) consists of:
g p g p 1 : top speed 85 km / h , g p 2 : kinematic top speed 85 km / h , g p 3 : acceleration time for 0 80 km / h 80 s , g p 4 : gradeability at 85 km / h 1 % , g p 5 : gradeability from standstill 10 % , g p 8 : 10 P m k 550 k , g p 9 : 100 E b 450 , g p 10 : 0.1 r l ( x l ) 40 l , h p h p 1 : r f = 1 .
The first five inequality constraints, g p 1 g p 5 , are target performance and drivability constraints. The first constraint, g p 1 , is a constraint on the minimal peak power of the electric machines to able to drive at the required top speed. The second constraint g p 2 is a kinematic constraint on the top speed, where in at least one of the gear positions, the maximum speed of all motors should not be exceeded when driving at the required top speed. The last three inequality constraints are constraints on the value of the design variables.
The equality constraint, h p 1 , is related to the assumption that the combined ratio of the gearbox and final drive is optimized, in the case of a central drive topology. The set of equality constraints can also include symmetry constraints, depending on the powertrain topology; for example for T D 1 in Figure 2, where equally sized electric machines (and gearboxes) must be connected to each wheel of the same driven axle. However, this only holds for distributed topologies. Keep in mind that the formulation of the constraints in (13) can vary per topology, due to the location of the gearbox(es) and number of electric machines.

2.5. Inner Loop Control Optimization

Consequently, for each choice of values for the component sizing parameters, x p , in the inner loop of (8), one needs to solve the following control problem:
min x c ( t ) C e ( x c ( t ) | x p , Λ ( t ) ) , subject to g c , ξ ( t ) [ ξ ̲ , ξ ¯ ] , h c 1 : ξ ˙ ( t ) f ( ξ ( t ) , x p , x c ( t ) , Λ ( t ) ) = 0 , h c 2 : ξ ( t 0 ) ξ 0 = 0 .
Here, ξ ( t ) R + (-) refers to the state of the system, in this case the battery state-of-charge (SOC). In this description, the bar above ( . ) ¯ (or below ( . ) ̲ ) a symbol indicates the maximum (or minimum) value, respectively. The first equality constraint refers to the system dynamics that have to be satisfied, and the second constraint refers to the initial SOC of the battery. The control problem is solved over a drive cycle, mentioned by Λ ( t ) R 2 , for which the vehicle speed v R + (km/h) and the road slope α R (rad) are given by:
Λ ( t ) = v ( t ) α ( t ) , t [ t 0 , t end ] ,
where t 0 (s) is the start time of the drive cycle and t end (s) the end time.
As is the case for the component sizing variables, which control variables are present depends on the specific powertrain topology T i . j l . When all control variables are present, the vector of control variables consists of:
x c ( t ) = { u ts ( t ) , x g ( t ) } ,
where u ts ( t ) [ 0 , 1 ] R + n m 1 (-) is the torque split over two electric machines, defined as:
u ts ( t ) = τ m 1 ( t ) τ m 1 ( t ) + τ m 2 ( t ) ,
where τ m 1 R (Nm) is the torque generated by EM1 and τ m 2 R (Nm) by the second electric machine EM2. In case u ts = 0 , only EM2 is used, whereas the case u ts = 1 means only using EM1. If more than two electric machines are present in the powertrain, for each extra electric machine, an extra torque split variable is added.
The gear position as the controlled state becomes x g ( t ) N n t for all gearbox units,
x g ( t ) = { x 1 ( t ) , , x n t ( t ) } ,
where x l ( t ) N (-) is the chosen gear position for each gearbox unit. The controlled input for the gearbox unit relates to the gear shift command u l ( t ) Z (-),
x l ( t + ) = x l ( t ) + u l ( t ) ,
whereby there are no limitations on the gear shift command. t indicates the time before the gear shift event and t + after. x l can take the following values,
x l ( t ) { 1 , , j l } ,
where j l is the total number of gears per gearbox.

Constraints

The set of constraints, g c , defined in the control design problem, includes the physical limitations of the components such as the maximum torque and speed for each electric machine ( g c 1 , g c 2 , g c 3 ) and constraints on the battery currents ( g c 4 , g c 5 ). These constraints are take into account inside the simulation model.
g c g c 1 : max . electric machine torque k , g c 2 : min . electric machine torque k , g c 3 : max . electric machine speed k , g c 4 : max . battery charge current , g c 5 : max . battery discharge current .

2.6. Optimization Framework

To solve the (nested) design optimization problem defined in (8) for each feasible powertrain configuration, the optimal component sizing problem is solved using particle swarm optimization (PSO) [27], and the control problem is solved using a local minimization of the control objective with respect to its input(s). The same algorithm settings were used for all powertrain configurations. To increase the chance of finding the global optimal solution for each design optimization, the optimization was run multiple times, and the results giving the lowest objective value were used.

3. System Modeling

To optimize the control variables for a given drive cycle, a simulation model was required. For this, a backward facing quasi-static simulation model of the vehicle was used, since the predictions on the overall energy consumption of these type of models were sufficiently accurate [28]. The model was implemented in MATLAB©, as was the optimization process and the analysis of the results. In the rest of this section, the powertrain component models used in the simulation model are discussed separately.

3.1. Vehicle Road-Load Model

The tractive force required to drive the vehicle, F r R (N), taking into account the rolling friction, aerodynamic friction, and gradient resistance, can be calculated as [28],
F r ( t ) = m v c r g cos ( α ( t ) ) + g sin ( α ( t ) ) + d v d t + 1 2 ρ A c d v ( t ) 2 ,
where the variables and their values can be found in Table 2. Furthermore, the angular speed ω w R + (rad/s) and torque demand τ w R (Nm) at the wheels are calculated by:
ω w ( t ) = v ( t ) r w ,
τ w ( t ) = F r ( t ) r w ,
where r w (m) is the wheel radius as defined in Table 2.

3.2. Final Drive Model

In case a final drive transmission is present in the powertrain topology, the final drive is modeled using a fixed efficiency, η f = 0.97 R + (-), for which the following holds,
if FD present ω f ( t ) = ω w ( t ) r f , τ f ( t ) = τ w ( t ) r f η f sign ( τ w ( t ) ) .
Here, ω f (rad/s) and τ f (Nm) are the angular speed and torque before the final drive unit, which depending on the topology are either equal to the torque and speed of the electric machine connected to the final drive or of the gearbox connected to the final drive. When no final drive is present, (25) becomes:
if FD not ̲ present ω f ( t ) = ω w ( t ) , τ f ( t ) = τ w ( t ) 2 .

3.3. Transmission Model

A transmission unit (TR), i.e., a multi-speed gearbox, is modeled using a fixed efficiency for each gear, which varies with the ratio of the gear. It is assumed that for larger gear ratios, multiple gear-stages are required to realize the ratio. Herein, it is assumed that one gear-stage, i.e., one pair of gear wheels, can realize maximally a ratio of four [29]. Furthermore, it is assumed that each gear-stage has a constant efficiency of η g p = 0.985 (-) [29], including losses from bearings and gear teeth contact. In case a final drive is present, it is taken into account that with a single stage final drive, a maximum ratio of seven can be achieved [30], which implies lower gear ratios for the gearbox. The efficiency per gear, η l ( x l ) , is:
η l ( x l ) = η g p if r l ( x l ) r f 4 , η g p 2 if 4 < r l ( x l ) r f 16 , η g p 3 if r l ( x l ) r f > 16 ,
where r l ( x l ) is the gear for which the efficiency is calculated. The value of r f equals one in the case of a distributed topology and seven in the case of a central drive topology. The angular velocity, ω ti l , and torque, τ ti l , at the gearbox inlet for gearbox unit l is given by:
ω ti l ( t ) = ω to l ( t ) r l ( x l ( t ) ) ,
τ ti l ( t ) = τ to l ( t ) r l ( x l ( t ) ) η l ( x l ( t ) ) sign ( τ to l ( t ) ) ,
where ω to l (rad/s) and τ to l (Nm) are the angular speed and torque at the gearbox output. Note that the gear ratio and efficiency of gearbox l depend on the chosen gear position x l at any given time t. Depending on the specific topology of the vehicle, the torque and speed at the gearbox outlet are either equal to the torque and angular speed of another gearbox or equal to ω f (rad/s) and τ f (Nm).

3.4. Electric Machine Model

Each electric machine (EM) was modeled using a lookup map of the efficiency, η m k (-), which depended on the electric machines’ torque, ω m k (rad/s), and angular speed, τ m k (Nm). An efficiency map of a passenger car was used as a base for the model (see Figure 3), as reliable open-source efficiency maps for higher power commercial vehicle electric machines were unavailable. Here, a “medium” type of electric machine was selected, as opposed to, typically, a low-speed-high-torque or a high-speed-low-torque machine, to investigate in a more balanced way the effect of a transmission and topology choice. This would otherwise lead to (largely) compromising the vehicle performance or reversibly would lead to possible excessive machine or transmission scaling. The losses of the power electronics (Inverter) are included in this efficiency map.
The electric power consumption of the electric machine P m k , e (W) is calculated as,
P m k , e ( t ) = ω m k ( t ) τ m k ( t ) η m k sign ( τ m k ( t ) ) ,
η m k = f m k ( ω m k , τ m k s m k ) .
Depending on the specific powertrain topology, the electric machine torque is determined from the model of the gearbox connected to the electric machine or from the final drive model.
The base efficiency map was scaled along the torque axis to scale the performance of the electric machine for different sizes [28]. This method of scaling is common practice in powertrain design studies. The scaling of the efficiency map was performed using the scaling variable s m k R + (-),
s m k = P m k P m k , base ,
where P m k is the peak mechanical power of the resized electric machine and P m k , base (W) is the peak mechanical power of the base electric machine. The maximum torque of the scaled machine τ ¯ m k is calculated as:
τ ¯ m k = s m k τ ¯ m .
The same should be applied to the minimum machine torque τ ̲ m k . In the case of topologies without a final drive and where the wheels were individually driven by the electric machines, i.e., using hub electric machines, the same efficiency map of the electric machine was used.

3.5. Battery Model

The battery pack of the vehicle was modeled using an equivalent circuit model [28]. The open circuit voltage U oc R + (V) and resistance R R + ( Ω ) of the complete battery pack depended on the configuration of the battery, i.e., the number of cells in series, n se N (-), and in parallel, n pa N (-),
U oc = n se U oc , cl ,
R = n se n pa R cl ,
where U oc , cl (V) and R cl ( Ω ) are the open-circuit voltage and the internal resistance of a single cell, respectively. It was assumed that these values were equal for all cells in the battery pack. Both U oc , cl and R cl were modeled as a function of the battery state-of-charge ξ . The influence of the battery temperature and battery aging were not considered in this model. The state-of-charge of the battery was determined using Coulomb counting, as:
ξ ( t ) = t = t 0 t I b ( t ) 3600 Q 0 d t + ξ ( t 0 ) ,
where Q 0 R + (Ah) is the nominal battery capacity, which depends on the nominal cell capacity Q 0 , cl (Ah) as:
Q 0 = n pa Q 0 , cl .
The battery current I b R (A) is calculated using:
I b = η c U oc U oc 2 4 R P b 2 R ,
where η c (-) is the Coulombic efficiency that differs for charging and discharging of the battery [28],
η c = 1 if I b 0 , 0.98 if I b < 0 .
The power at the battery terminal, P b , is:
P b = k = 1 n m P m k , e + P a .
where P a (W) denotes the power consumption of the auxiliaries. Now, the internal battery power P s , which is needed to calculate the average energy consumption, E v , in (3), is calculated using:
P s = U oc I b .
Scaling of the battery size is done by scaling total energy content, which is scaled with the total number of parallel branches in the battery n pa according to:
n pa = E b E b , 0 .
where E b , 0 (kWh) is the energy content of one parallel battery branch:
E b , 0 = Q 0 , c l U o c .
The number of cells in series n se was kept fixed, meaning the battery voltage was kept fixed during scaling. The battery parameters are given in Table 3, and a plot of the battery cell resistance and open-circuit voltage versus the state-of-charge is shown in Figure 4.

3.6. Auxiliary Units

The auxiliaries considered here were the electric air compressor (EAC), electro-hydraulic power steering (EHPS), electronic control unit (ECU), air conditioning unit (HVAC), and a cooling circuit. All were modeled by a fixed average power consumption and added up to a total auxiliary power of P a = 4.86 kW.

3.7. Vehicle Mass Model

The change in the total mass of the vehicle, m v R + (kg), by scaling the powertrain components, is taken into account using,
m v = m ta + m pt + m tr + m ca ,
where m ta is the mass of the tractor without powertrain, m tr the trailer mass, and m ca the transported cargo mass. Moreover, the mass of the powertrain m pt (kg) scales with the size of the powertrain components according to:
m pt = m b + k = 1 n m ( m m k + m i k ) + l = 1 n t m g l .
Here, m b is the mass of the battery, m m k the mass for electric machine k, m i k the inverter mass belonging to this electric machine, and m g l the mass per gearbox unit as detailed in Table 4. The gearbox mass model depends on the maximum torque capability of the electric machine connected to that specific gearbox [30].

3.8. Cost Model Parameters

The cost of the powertrain of the vehicle C pt (euro), from (4), is given as:
C pt = C b + k = 1 n m ( C m k + C i k ) + l = 1 n t C g l ,
with C b the cost of the battery, C m k the cost per electric machine k, and C i k the cost per inverter. The cost of the transmission unit, or the gearbox, is denoted by C g l . The linear cost model parameters of the electric machine, inverter, and battery models were estimated based on multiple references [11,33,34,35,36,37] and are listed in Table 5. Hereto, the values of all references are plotted versus the year the cost models applied. Based on this, estimates were made of the value of these parameters in 2019. As no reliable estimates could be made on the constant part of the cost models, they were omitted. However, they can easily be added in the future. The price of electricity, c e , was based on the European average of non-household energy prices of 2018 [38].

3.9. Drive Cycle

Typically, for large-scale powertrain design studies, to constrain the simulation times, a shorter, yet representative, long-haul drive cycle is used. Similarly, in this work, the drive cycle used (Figure 5) was obtained from the Vehicle Energy consumption Calculation Tool (VECTO) [39], which is a simulation tool in development for CO2 declaration of commercial vehicles above 7.5 tons, and it contains no on-route charging.

4. Design Results and Analysis

For each of the 44 powertrain configurations, the results of the optimization problem (8) are discussed next. It was important to investigate each design choice separately and to identify the potential trade-offs between them.

4.1. Total-Cost-of-Ownership

One of the most important parameters when designing a new vehicle concept is the total-cost-of-ownership, as defined in (1), which should be be as low as possible. The optimal value of the TCO, denoted by C TCO , for each of the 44 powertrain configurations optimized is shown in Figure 6. In this figure, markers are used to indicate the different types of gearboxes used.
The topology with the lowest TCO was T D 1 , which had one electric machine and gearbox for each of the two driven wheels. This topology was equipped with a two-speed; hence, T D 1.2 , gave overall the lowest TCO with a value of C TCO = 187.7 k euro. The topology with the highest TCO was T C 5 , which was a central drive topology with in total two electric machines and two gearboxes that were placed in series. This topology also showed the largest variation in C TCO with the number of gears in the gearboxes. This topology was equipped with two three-speeds; hence, T C 5.33 , resulted in being the highest TCO overall with a value of C TCO = 199 k euro. Thus, between the powertrain configurations with the highest and lowest C TCO , the difference was −11,200 euro, meaning that for the set of topologies studied, choosing a different topology, and gearbox type, could result in a variation of −5.6% in the optimal TCO value. Other performance metrics and the values of the optimal sizing parameters x p for each topology can be found in Appendix B. The group of powertrain configurations with topology T C 5 , which are the three most right vertical lines in Figure 6, was clearly the least favorable option regarding TCO. Therefore, this topology will be left out for the rest of the analyses discussed in this section.

4.2. Optimal Component Sizing

Among the investigated set of configurations, the optimal battery size E b , shown in Figure 7, varied from 221 kWh to 210 kWh, which was a relative variation of −5.1%. The results favored the smallest battery size possible. This was a trade-off between the battery size influencing the (i) battery efficiency, (ii) the vehicle mass, and thereby vehicle energy consumption; and, (iii) the vehicle cost. The optimal battery size reached a boundary optimum, where the battery size was constrained by the lower state-of-charge limit at the end of the drive cycle. The distributed drive configurations having smaller optimal battery sizes are located more towards the left on the horizontal axis in Figure 7, due to the lower energy consumption of this group of topologies.
Another important parameter was the total electric machine power, P ^ m (Figure 7, vertical axis), which was the sum of all electric machines’ mechanical peak powers. The value of P ^ m varied between 497 kW ( T C 2.4 ) and 268 kW ( T D 1.3 ), which was a reduction of −46%. The topologies T D 1 and T C 1 equipped with a multi-speed gearbox seemed to require smaller total electric machine powers.

4.3. Central Drive Versus Distributed Drive Topology

The first of three topological design choices to discuss is the choice of a central or distributed drive powertrain. Therefore, in Figure 8, the C TCO of the configurations with a distributed (Figure 8, left) and central drive (Figure 8, right) powertrain are grouped into two separate plots. The value of C TCO is plotted against the number of gears in the gearbox.
The group of distributed configurations (Figure 8, left) had an overall lower C TCO (−1.0%) than the group of central configurations (Figure 8, right). In order to find the cause of this difference, the two parts of which C TCO consisted, the electricity cost C e and vehicle cost C v , are plotted against one another in Figure 9. As can be seen in Figure 9, both groups are located at the same height on the vertical axis, meaning the vehicle cost C v was similar for both groups. Thus, the difference in C TCO between both groups was mainly related to the difference in electricity cost, C e (Figure 9, horizontal axis), for which the distributed configurations are located more towards lower values. This difference in electricity cost was largely coupled to the fact that the distributed powertrains did not include a final drive. Therefore, they had a higher powertrain efficiency and inherently a lower energy consumption.

4.4. Gearbox Type and Location

The second topological design choice investigated was the choice of the total number of gears in the gearbox. In Figure 8, it can be observed that the influence of the chosen number of gears, visible as the variation over the horizontal axis of Figure 8, varied per topology. Using a two-speed had a reduced C TCO over using a one-speed for all topologies, except for T C 4 (which was unexpected and where the multi-layer algorithm had most likely reached a local minimum). For all configurations that showed a reduction in C TCO from a one-speed to a two-speed, the average reduction in C TCO was −0.6%. The reduction in C TCO was slightly larger for the central drive topologies (−0.9%) compared to the distributed topologies (−0.4%).

4.4.1. Two-Speed Gearboxes

The configuration with a two-speed had the lowest C TCO for five out of eight topologies studied. Adding more gears to the gearbox, i.e., using a three-speed or four-speed, did not significantly reduce C TCO (for most topologies). Thus, the largest impact on C TCO was achieved by using a two-speed instead of a single-speed. This was related to the fact that using a second gear increased the freedom in placing the working points of the electric machine. In the case of a single-speed gearbox, the choice in the gear ratio was limited as both the top speed and maximum torque constraints had to be satisfied by the choice of this gear ratio. However, in the case of a two-speed gearbox, only one gear had to satisfy these constraint or each gear only one constraint. Therefore, there was more freedom in the choice of the gear ratio value, which facilitated that the working points could be placed closer to the region of high efficiency of the electric machine.

4.4.2. Three-Speed Gearboxes

Adding a third gear to the gearbox reduced the electricity cost compared to a two-speed for some of the topologies. However, for most topologies, the reduced electricity cost did not outweigh the increased vehicle cost due to increased transmission cost because of the use of an extra gear. This was only the case for T D 2 and T C 4 , for which a three-speed resulted in the lowest C TCO . The average change in C TCO when going from a two-speed to a three-speed, for the group of configurations studied, was 0.0%. The only benefit of adding a third gear was improved gradeability. Topologies T D 2 , T C 2 , and T C 3 had a significantly higher gradeability performance when using a three-speed than when using a two-speed (Table A2 in Appendix B). In the case of a three-speed, the third gear was used to place the working points of the electric machine at low vehicle speeds into the region of the efficiency map with high efficiencies. As this required one of the gears to have a higher gear ratio, the gradeability performance was improved as well.

4.4.3. Four-Speed Gearboxes

When adding a fourth gear to the gearbox, for most topologies, there was an increase in C TCO , an average increase of +0.8% among the configurations that showed an increase in TCO. The steep increase for T D 2 and T D 3 might be due to the convergence to local minima. Therefore, no clear conclusion could be drawn on the effect of adding a fourth gear on the TCO, and the difference in this between distributed and central topologies. This would require more in-depth analysis; however, the general trend seemed to indicate an increase in TCO when adding a fourth gear.

4.4.4. Gearbox Location

In the case that a topology incorporated multiple electric machines, the question arose where to locate the gearbox unit in the powertrain and how many gearboxes should be added. For example, in the case of a central drive powertrain with two EMs, should one gearbox be incorporated that is used by one electric machine ( T C 3 ) or used by both electric machines ( T C 2 ), or should two gearboxes be used ( T C 4 )? Results showed that the lowest C TCO was achieved for T C 2 , when both EMs shared the same gearbox. This resulted in a lower C TCO , on average −1.0%, as when only one EM was equipped with a gearbox, T C 3 . The reason a lower C TCO was achieved compared to supplying each electric machine with its own gearbox, T C 4 , was largely related to the vehicle cost, C v . The vehicle cost was lower for T C 2 due to lower gearbox and electric machine costs. The latter was due to a lower total EM machine peak power P ^ m for T C 2 .

4.5. Single Versus Multiple Electric Machines

The last design choice involved the number of electric machines in the powertrain. Hereto, topology T D 1 (two EMs) was compared to T D 2 (four EMs), both distributed powertrains. Next to that, T C 1 (one EM) and T C 2 (two EMs) were compared for the central drive powertrains. Whereas topology T D 1 contained two EMs, both were equally sized due to the symmetry constraint between the left and right wheel of the vehicle. Therefore, there was only one design variable for the EM size for this topology.
In the case of a distributed topology, the option with the least number of electric machines, T D 1 , had the lowest C TCO . This could be concluded from the left plot in Figure 8 where the line of T D 1 was located below the line of T D 2 . T D 1 (two EMs) had a higher electricity cost, C e , compared to T D 2 (four EMs). However, due to lower vehicle cost, T D 1 resulted in a lower C TCO . Hence, the reduced electricity cost due to the use of extra electric machines did not outweigh the required increase in vehicle cost. Hence, for the distributed topologies, adding extra EMs increased the TCO, with on average +0.7%. A different observation was made for the central topologies. Among the group of central drive topologies, the option with the most electric machines, T C 2 (two EMs, different sizing) had the lowest C TCO . Thus, for the central drive topologies, the reduced electricity cost outweighed the required investment cost of adding an extra electric machine to the power. Hence, for the central topologies, adding an extra EM reduced the TCO, by on average −0.3%.
The total machine power P ^ m , for both comparing T D 1 to T D 2 and T C 1 to T C 2 , was higher for the options with more EMs in the powertrain, T D 2 and T C 2 . Adding extra EMs to the powertrain increased the optimal total machine power. When looking at the individual sizes of the electric machines, plotted in Figure 10, it can be observed that the electric machines were not equally sized. For two powertrain configurations, the sizes were close to equal, a difference of around 10 kW; however, for most configurations, the difference was larger than 10 kW. The difference in sizing could be observed in Figure 10 by the deviation of the points from the red dotted line, which indicated where both EM sizing parameters would be equal. A possible explanation of the difference in EM sizing was related to the existence of multiple local minima with (almost) equal for the topologies with multiple EMs, e.g., topologies T D 2 , T C 2 , and T C 4 . Either using equally sized or a large difference in size of electric machines resulted in similar objective function values. Therefore, from the objective point of view, there was no preference for either solution, and some optimizations converged to one solution, e.g., T D 2.4 and T C 4.11 , while some to the other solution, e.g., T D 2.3 and T C 4.22 .
Note that P m 2 had values close to 10 kW for T D 3 and T C 3 . This was the value of the lower bound of the EM sizing parameter. In both topologies, the machine(s) with sizing parameter P m 2 was directly connected to the final drive, and therefore had a transmission ratio of one to the driven wheels. Due to this, the work points for these machines were outside the feasible working range of the machine; this machine was minimally used on the drive cycle, and its size was minimized in order to reduce the vehicle cost as much as possible.

4.6. Comparison of the Influence of Design Choices

The question now arose about which of the three topological design choices had the largest impact on C TCO . This comparison was performed from the perspective of the most widely used powertrain configuration by truck manufacturers, T C 1.1 . In Figure 11, the influence of each of the three topological design choices is depicted, showing the effect of: adding one gear to the gearbox ( T C 1.2 ); adding an extra electric machine to the topology ( T C 2.1 ); or moving to a distributed type of topology ( T D 1.1 ). For each of these design choices, the change in C TCO , C e and C v is indicated by Δ C TCO (%), e.g., this value for the design choice from T C 1.1 to T D 1.1 is calculated as:
Δ C TCO = C TCO ( T D 1.1 ) C TCO ( T C 1.1 ) C TCO ( T C 1.1 ) · 100 ,
The values of Δ C e and Δ C v are calculated in the exact same way.
Seen from the perspective of T C 1.1 , the largest reduction in C TCO was realized by going to a distributed topology (−1.1%), followed by the addition of an extra gear (−0.8%), and last, by adding an extra EM (−0.1%). This was in line with the results from Section 4.3, Section 4.4 and Section 4.5, where also the largest reduction on average on the TCO was shown for going to a distributed topology (on average −1.0%), followed by adding a gear (on average −0.6%), and the lowest impact for adding an EM, on average −0.3% for central topologies. The latter even showed an increase in TCO for distributed topologies (on average +0.7%). Therefore, it seemed that for the application studied, the largest impact on TCO was realized when going from a central to a distributed topology, followed by the use of a multi-speed gearbox and last by the addition of extra electric machines.

4.7. Results Discussion

Note that each optimal set of parameters x p for each powertrain configuration was subject to the grid defined for the components and the parameters chosen for the optimization algorithm. The same optimization algorithm parameters were used for all powertrain configurations, and this (e.g., optimization stopping criteria) could also influence x p . For example, the population size was slightly increased for the optimization of T C 4 and T C 5 , due to the higher number of design parameters. The use of PSO did not guarantee finding the exact global optima; however, PSO converged to a solution in the proximity of the optimum, as was shown in an example in [27].
In order to check whether the algorithm converged to a minimum, a sensitivity study of the objective to the sizing parameters was performed. Hereto, the change in the objective value was calculated for a change of 5%, in both directions, in each of the sizing parameters x p . Each time, only one parameter was varied, it being a one-at-a-time sensitivity study. The maximum variation in the objective value due to variation in the design parameters was on average −0.09%, with a maximum of −0.31% for T C 3.4 , for the distributed configurations and most central configurations. Mainly further down-scaling the battery or electric machine(s) had the largest influence on the objective value. The group of configurations of T C 4 and T C 5 showed larger variations in the objective, ranging from −0.11% to −0.8%. These topologies more often converged to local minima. This might be related to the high number of design parameters and the possible existence of multiple local minima when multiple electric machines were present in the powertrain, as discussed in Section 4.5. Apparently, the tuned PSO algorithm was not robust enough to be properly applied to all design problems. The large variation in the shape of the objective functions among different powertrain configurations would require re-tuning of the algorithm; something that should be investigated more in future research. For this paper, this meant that the exact values of the results of especially topologies T C 4 and T C 5 should be observed with more care.
Next to that, the value of the optimal battery size in kWh, discussed in Section 4.2 and shown in Figure 7, should be interpreted with some care. The focus should be more on the variation of the battery size among different powertrain configurations, rather than the absolute value. The optimization of the battery in this paper focused on minimizing system efficiency and did not include effects related to battery aging and charging of the battery. These effects will influence the exact value of the battery size and will most likely result in larger batteries when being implemented in a real vehicle. To determine the exact influence of these effects, more extensive research is required.

4.7.1. Optimal Gear Use

Next to the design parameters, also the optimal control trajectories were studied. In Figure 12, the use of the gears and electric machines is shown for two powertrain configurations. In this figure, the road-load force at the wheels is plotted versus the speed of the vehicle, showing the load at the wheels of the vehicle.
Regarding the use of gears, shown in the lower two plots of Figure 12, the most used gear was the last gear having the lowest gear ratio. This gear was engaged for around 90% of the drive cycle. This also held for the other topologies studied. As can be observed in Figure 12, the most used gear was mainly used at higher vehicle speeds, above 65 km/h. This coincided with the vehicle being driven at speeds above 65 km/h for 88% of the time on the drive cycle. Using the second gear, the work points, corresponding with higher vehicle speeds, of both EMs were located as close as possible to the high efficiency region of the map (see the location of the dark crosses in Figure 13). In Figure 13, the corresponding work points for each of the electric machines in the powertrain are shown, with different colors indicating the gear engaged for each work point.

4.7.2. Torque Split

Regarding the use of the electric machines during motoring mode, when the vehicle was driven, plots were made to show whether only the first EM was used (EM1 only), the second EM was used (EM2 only), or both EMs were used simultaneously (dual motor). Keep in mind that for distributed topologies in the “EM1 only” mode, both EM1 machines were used, as shown in Figure 2. The use of these three modes is indicated with different markers in the upper two plots of Figure 12 and in percentages in Figure 11. The most used mode during driving was the dual motor mode for both configurations studied in this section, followed by the use of the largest EM. During braking, mainly the smallest EM was used to recuperate energy.
To conclude on the results, when only looking at total-cost-of-ownership for choosing the design of an electric powertrain for a long-haul truck, T D 1.2 was the most optimal. Furthermore, from a central drive vs. distributed drive point of view, the later was preferred. Regarding the number of gears in the gearbox, the option with a two-speed seemed to result in the lowest overall TCO. Adding a third gear would improve the gradeability of the vehicle if required; however, for most topologies, this came at the cost of an increased TCO. Regarding choosing the number of electric machines, for both central and distributed drive topologies, the option with two electric machines resulted in the lowest TCO, where for the distributed topologies, they were equal in size. As can be observed throughout this section, in most cases, the conclusion and effects of a design choice were different for the central and distributed drive powertrains. Hence, there existed an interaction between the design choice of a distributed or central drive topology and the other two topological design choices.

4.8. Alternative Optimization Objective: Energy Consumption

As C TCO highly depended on the cost model assumptions, the optimization described in Section 2 was re-performed with only the energy consumption as objective; thus, C TCO = C e in (7) and (8). The optimal energy consumption E v , shown in Figure 14, varied from 181.7 kWh/100 km ( T C 5.33 ) to 169.5 kWh/100 km ( T D 2.2 ), which was a variation of −6.7%. In the case of energy consumption as the optimization objective, the topology with four electric machines was a better choice ( T D 2.2 ) instead of the two EM topology ( T D 1.2 ) in the case of C TCO as the objective. This was due to the fact that in the case of energy consumption as the objective, there was no cost penalty in the use of more powertrain components.
Similar to the case of C TCO being the objective to minimize, the largest impact on energy consumption regarding the choice of number of gears in the gearbox was when going from a one-speed to a two-speed. Regarding the best choice of the number of gears, here as well, the two-speed had the lowest TCO for the distributed topologies. However, for T D 1 , the three-speed had a marginally lower energy consumption. However, for the central topologies, no a clear conclusion could be drawn on the best gearbox type, but the preference seemed to be more towards four speeds. Hence, more gears as when C TCO was the objective.
Another observation was that in the case of energy consumption as the objective, higher total EM powers were found. Compared to when TCO was the objective function to minimize, there was no direct penalty on the objective value when increasing the size of the electric machine, as electric machine cost was not part of the objective function. Larger electric machine sizes increased the average electric machine efficiency over the drive cycle and, therewith, resulted in a lower energy consumption.

5. Conclusions

The main goal was to investigate the effect of discrete topological design choices for electric powertrains for heavy-duty trucks. Therefore, a powertrain design optimization study, centered on the costs of operation and components, was performed for a group of various powertrain configurations. A fixed set of 44 feasible configurations was considered, constructed from eight topologies for which each a different number of discrete gear ratio values in the gearbox was considered. In particular, this set was considered in order to analyze the topological design choices related to: (a) a central or distributed drive (individually-driven wheels) powertrain, (b) a single or a multi-speed gearbox, and (c) using single or multiple electric machines. Results showed that the choice of topology and discrete gear ratio in the gearbox (for the set studied) could result in a variation in the optimal total-cost-of-ownership of at maximum 5.6% and a variation in the sum of the electric machines’ peak power specification of 46%. The topology resulting in the lowest total-cost-of-ownership was a distributed topology with in total two electric machines and equipped with a two-speed gearbox per machine. Additionally, studying the set of powertrain configurations showed that the largest average reduction in total-cost-of-ownership was achieved by: (i) choosing a distributed drive powertrain topology over a central drive topology (−1.0%); followed by (ii) using a multi-speed gearbox, preferably a two-speed, over a single-speed gearbox (−0.6%); and the smallest impact was shown for (iii) using multiple electric machines over using one single machine (−0.3%) for central drive topologies, which even showed an increase in total-cost-of-ownership for distributed topologies (+0.7%). The framework, used in this paper, could easily be applied to different (vehicle) applications by changing the model parameters and the drive cycle. Therefore, in future research, the effect of the chosen model parameters on the results will be further examined in order to understand what the influence is of the vehicle, cost, mass, and drive cycle parameters on the three topological design choices and on the choice of the optimal topology ultimately.

List of Sub-Scripts

Author Contributions

Conceptualization, F.J.R.V. and T.H.; formal analysis, F.J.R.V.; software, F.J.R.V.; supervision, E.S. and T.H.; visualization, F.J.R.V.; writing, original draft, F.J.R.V.; writing, review and editing, E.S. and T.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This research was conducted at Eindhoven University of Technology.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BABattery
EACElectric air compressor
ECUElectronic control unit
EHPSElectric hydraulic power steering
EMElectric machine
FDFinal drive
HVACair conditioning unit
PSOParticle swarm optimization
TCOTotal-cost-of-ownership
TRTransmission
VECTOVehicle Energy Consumption Calculation Tool
WDriven wheel
List of Symbols
Afrontal aream2Q0battery capacityAh
ceelectricity costeuro/kWhrffinal drive ratio-
Ccost parameterseurorlvector with gear ratios of gearbox l-
cdair drag coefficient-rwwheel radiusm
crrolling resistance coefficient-Rinternal battery resistanceΩ
DdistancekmRMatrix with all gear ratios in the powertrain-
Ebbattery energy contentkWhSmkelectric machine scaling parameter-
Evenergy consumption, energykWh/100 kmttimes
Frroad-load forceNtaccacceleration times
ffunction indication-Ttopology-
ginequality constraints-Tset of topologies-
ggravitational accelerationm2/sutstorque split control variable-
hequality constraints-Uocopen circuit voltageV
itopology name-vvehicle speedm/s
IcurrentAxvector with design variables-
jlnumber of gears in gearbox uni l-xlgear position of gearbox l-
kindex of the electric machine-xgvector with gear positions of all gearboxes-
lindex of the gearbox unit-yyearyear
mmasskgαroad sloperad
ntotal number of instances-ηefficiency-
PpowerWΛdrive cycle-
Pmkmechanical peak power for electric machine kWξbattery state-of-charge-
Pmk,eelectric power for electric machine kWρair densitykg/m3
Pmvector with all electric machine peak powersWτtorqueNm
P ˙ m sum of electric machine peak powersWωangular speedrad/s
List of Sub-Scripts
0begin time, nominal, base vehicle costlgearbox index
85related to 85 km/h constraintmelectric machine
accaccelerationocopen-circuit
bbatterypplant
cdrive cyclepaparallel
ccontrolptpowertrain
Ccentralrtotal number of gear ratios
cacargosinternal battery power
chchargingseseries
clsingle battery cellstrelated to standstill constraint
ddrag coefficienttnumber of transmission units
Ddistributedtatractor
didischargeTCOtotal-cost-of-ownership
eelectricity cost, electric powertitransmission inlet
eleconomical life timetotransmission output
endend timetoptop speed
ffinal drivetrtrailer
ggearboxtstorque split
gpgear pairvvehicle
iinverterwwheel
itopology nameyyearly
kelectric machine index

Appendix A. Design Parameters Per Topology

Table A1 shows the specific plant design and control parameters per topology.
Table A1. Overview of component sizing and control design parameters per topology.
Table A1. Overview of component sizing and control design parameters per topology.
TopologyComponent Sizing ParametersControl Parameters
T i . j l
D1. j 1 x p = { E b , P m 1 , r 1 } x c ( t ) = { x 1 ( t ) }
D2. j 1 x p = { E b , P m 1 , P m 2 , r 1 } x c ( t ) = { u ts ( t ) , x 1 ( t ) }
D3. j 1 x p = { E b , P m 1 , P m 2 , r 1 } x c ( t ) = { u ts ( t ) , x 1 ( t ) }
C1. j 1 x p = { E b , P m 1 , r 1 } x c ( t ) = { x 1 ( t ) }
C2. j 1 x p = { E b , P m 1 , P m 2 , r 1 } x c ( t ) = { u ts ( t ) , x 1 ( t ) }
C3. j 1 x p = { E b , P m 1 , P m 2 , r 1 } x c ( t ) = { u ts ( t ) , x 1 ( t ) }
C4. j 1 j 2 x p = { E b , P m 1 , P m 2 , r 1 , r 2 } x c ( t ) = { u ts ( t ) , x 1 ( t ) , x 2 ( t ) }
C5. j 1 j 2 x p = { E b , P m 1 , P m 2 , r 1 , r 2 } x c ( t ) = { u ts ( t ) , x 1 ( t ) , x 2 ( t ) }

Appendix B. Optimization Results Data

In this Appendix, all the optimization results are shown. Table A2 shows the calculated objective and performance for each optimized powertrain configuration, and Table A3 shows the optimal values of the design parameters.
Table A2. Additional performance parameters for the powertrain configurations optimized for total-cost-of-ownership, as discussed in Section 4.
Table A2. Additional performance parameters for the powertrain configurations optimized for total-cost-of-ownership, as discussed in Section 4.
TopologyPerformance Parameters
T i . j l E v Top SpeedGradeabilityAcceleration
0–80 km/h
Range
(kWh/100 km)(km/h)(%)(s)(km)
D1.1173.2143113798
D1.2171.7131114698
D1.3172.0126125198
D1.4171.4128114798
D2.1171.6148113598
D2.2170.9150113498
D2.3170.6135324099
D2.4170.7156382598
D3.1173.2111113698
D3.2173.2101254398
D3.3172.0104373798
D3.4172.71205032100
C1.1175.7145113598
C1.2175.7127124898
C1.3175.6125184998
C1.4174.7132154298
C2.1174.0155123099
C2.2172.61411237100
C2.3173.7129304498
C2.4174.4131454198
C3.1176.0147113698
C3.2174.3149173598
C3.3173.5150382998
C3.4175.61302144101
C4.11173.9153113298
C4.21175.694114399
C4.22174.6147403098
C4.24173.4136173998
C4.31173.4142173699
C4.32174.4136403899
C4.33173.8142433398
C4.34176.0129334499
C4.41173.6148373198
C4.44175.2135393899
C5.11177.2126116298
C5.12175.1151374198
C5.13173.81644829100
C5.14173.7146903398
C5.21175.5134115098
C5.22175.81355757106
C5.23176.8184902098
C5.31177.6141203698
C5.32177.81329041101
C5.33175.8148902998
C5.41175.217843498
Table A3. Optimal component sizing parameters for the powertrain configurations optimized for total-cost-of-ownership, as discussed in Section 4.
Table A3. Optimal component sizing parameters for the powertrain configurations optimized for total-cost-of-ownership, as discussed in Section 4.
TopologySizing
T i . j l E b P m 1 P m 2 r 1 ( 1 ) r 1 ( 2 ) r 1 ( 3 ) r 1 ( 4 ) r f r 2 ( 1 ) r 2 ( 2 ) r 2 ( 3 ) r 2 ( 4 )
(kWh)(kW)(kW)(-)(-)(-)(-)(-)(-)(-)(-)(-)
D1.1213182 7.6
D1.2211149 9.35.3
D1.3211134 11.35.10.9
D1.4210145 10.19.88.35.1
D2.1211111877.0
D2.2210122856.76.0
D2.32111114724.26.85.1
D2.421012211219.815.96.86.2
D3.1213188117.6
D3.22121471421.36.5
D3.32111641027.018.35.4
D3.42151757133.729.420.45.4
C1.1216391 7.5 1
C1.2216285 10.56.7 1
C1.3215272 16.812.86.1 1
C1.4214313 12.99.49.36.51
C2.12152781886.7 1
C2.22152491148.75.1 1
C2.321416413025.514.85.1 1
C2.42142971335.629.17.05.21
C3.1216368627.6 1
C3.22143487312.56.2 1
C3.32134151022.78.55.2 1
C3.42212951018.814.714.56.71
C4.112142462037.3 15.4
C4.212186526822.911.8 15.6
C4.222153931624.96.1 118.43.9
C4.242132736110.65.2 124.222.212.67.6
C4.312163561712.510.85.1 10.3
C4.3221518514926.719.75.1 132.27.8
C4.3321321016134.222.96.8 120.617.55.1
C4.3421814415123.614.39.4 133.310.68.55.1
C4.4121322518436.115.17.44.515.2
C4.442172735929.416.915.77.2130.229.919.67.0
C5.112181742100.7 18.8
C5.122151233140.1 127.85.2
C5.132171264070.6 123.512.45.1
C5.142132231791 120.118.910.65.3
C5.21215652573.41.4 16.0
C5.222331182131.10.9 138.86.3
C5.232172935500.50.4 139.36.86.1
C5.312183373128.717.310.3 10.5
C5.322242892523.610.51.2 14.40.6
C5.3321620821617.65.51.3 118.810.95.6
C5.4121510429928.816.412.91.315.2

References

  1. Verbruggen, F.J.R.; Hoekstra, A.E.; Hofman, T. Evaluation of the state-of-the-art of full-electric medium and heavy-duty trucks. In Proceedings of the 2018 31th International Electric Vehicle Symposium & Exhibition (EVS), Kobe, Japan, 30 September–3 October 2018. [Google Scholar]
  2. Silvas, E. Integrated Optimal Design for Hybrid Electric Vehicles. Ph.D. Thesis, Eindhoven University of Technology, Eindhoven, The Netherlands, 2015. [Google Scholar]
  3. Fathy, H.K.; Reyer, J.A.; Papalambros, P.Y.; Ulsov, A. On the coupling between the plant and controller optimization problems. In Proceedings of the 2001 American Control Conference (ACC), Arlington, VA, USA, 25–27 June 2001; Volume 3, pp. 1864–1869. [Google Scholar]
  4. Sen, B.; Ercan, T.; Tatari, O. Does a battery-electric truck make a difference?—Life cycle emissions, costs, and externality analysis of alternative fuel-powered Class 8 heavy-duty trucks in the United States. J. Clean. Prod. 2017, 141, 110–121. [Google Scholar] [CrossRef]
  5. Sripad, S.; Viswanathan, V. Performance metrics required of next-generation batteries to make a practical electric semi truck. ACS Energy Lett. 2017, 2, 1669–1673. [Google Scholar] [CrossRef] [Green Version]
  6. Çabukoglu, E.; Georges, G.; Küng, L.; Pareschi, G.; Boulouchos, K. Battery electric propulsion: An option for heavy-duty vehicles? Results from a Swiss case-study. Transp. Res. Part C Emerg. Technol. 2018, 88, 107–123. [Google Scholar] [CrossRef]
  7. Huismans, H. Electric Trucks: Wishful Thinking or the Real Deal? Master’s Thesis, Delft University of Technology, Delft, The Netherlands, 2018. [Google Scholar]
  8. Zhou, T.; Roorda, M.J.; MacLean, H.L.; Luk, J. Life cycle GHG emissions and lifetime costs of medium-duty diesel and battery electric trucks in Toronto, Canada. Transp. Res. Part D Transp. Environ. 2017, 55, 91–98. [Google Scholar] [CrossRef]
  9. Karlsson, A.; Domingues-Olavarría, G.; Alaküla, M. Alternative EV powertrain topologies designed for operation in a conductive electric road system. In Proceedings of the 2018 IEEE International Conference on Electrical Systems for Aircraft, Railway, Ship Propulsion and Road Vehicles & International Transportation Electrification Conference (ESARS-ITEC), Nottingham, UK, 7–9 November 2018; pp. 1–6. [Google Scholar]
  10. Morozov, A.; Humphries, K.; Zou, T.; Rahman, T.; Angeles, J. Design, Analysis, and Optimization of a Multi-Speed Powertrain for Class-7 Electric Trucks. SAE Int. J. Altern. Powertrains 2018, 7, 27–42. [Google Scholar] [CrossRef]
  11. Fries, M.; Wolff, S.; Horlbeck, L.; Kerler, M.; Lienkamp, M.; Burke, A.; Fulton, L. Optimization of hybrid electric drive system components in long-haul vehicles for the evaluation of customer requirements. In Proceedings of the 2017 IEEE 12th International Conference on Power Electronics and Drive Systems, Honolulu, HI, USA, 12–15 December 2017; pp. 1–141. [Google Scholar]
  12. Mareev, I.; Becker, J.; Sauer, D. Battery dimensioning and life cycle costs analysis for a heavy-duty truck considering the requirements of long-haul transportation. Energies 2017, 11, 55. [Google Scholar] [CrossRef] [Green Version]
  13. Kollmeyer, P.; McFarland, J.; Jahns, T. Comparison of class 2a truck electric vehicle drivetrain losses for single-and two-speed gearbox systems with IPM traction machines. In Proceedings of the 2015 IEEE International Electric Machines & Drives Conference (IEMDC), San Diego, CA, USA, 11–15 May 2015; pp. 1501–1507. [Google Scholar]
  14. Morozov, A.; Humphries, K.; Zou, T.; Martins, S.; Angeles, J. Design and optimization of a drivetrain with two-speed transmission for electric delivery step van. In Proceedings of the 2014 IEEE International Electric Vehicle Conference (IEVC), Florence, Italy, 17–19 December 2014; pp. 1–8. [Google Scholar]
  15. Tan, S.; Yang, J.; Zhao, X.; Hai, T.; Zhang, W. Gear Ratio Optimization of a Multi-Speed Transmission for Electric Dump Truck Operating on the Structure Route. Energies 2018, 11, 1324. [Google Scholar] [CrossRef] [Green Version]
  16. Verbruggen, F.J.R.; Rangarajan, V.; Hofman, T. Powertrain design optimization for a battery electric heavy-duty truck. In Proceedings of the 2019 American Control Conference (ACC), Philadelphia, PA, USA, 10–12 July 2019; pp. 1488–1493. [Google Scholar]
  17. Katrašnik, T. Analytical framework for analyzing the energy conversion efficiency of different hybrid electric vehicle topologies. Energy Convers. Manag. 2009, 50, 1924–1938. [Google Scholar] [CrossRef]
  18. Soriano, F.; Moreno-Eguilaz, M.; Alvarez, J.; Riera, J. Topological analysis of powertrains for refuse-collecting vehicles based on real routes–Part II: Hybrid electric powertrain. Int. J. Automot. Technol. 2016, 17, 883–894. [Google Scholar] [CrossRef] [Green Version]
  19. Xu, C.; Guo, K.; Yang, F. A Comparative study of Different Hybrid Electric Powertrain Architectures for Heavy-Duty Truck. IFAC-PapersOnLine 2018, 51, 746–753. [Google Scholar] [CrossRef]
  20. Dagci, O.H.; Peng, H.; Grizzle, J.W. Hybrid electric powertrain design methodology with planetary gear sets for performance and fuel economy. IEEE Access 2018, 6, 9585–9602. [Google Scholar] [CrossRef]
  21. Song, Z.; Hofmann, H.; Li, J.; Han, X.; Zhang, X.; Ouyang, M. A comparison study of different semi-active hybrid energy storage system topologies for electric vehicles. J. Power Sources 2015, 274, 400–411. [Google Scholar] [CrossRef]
  22. Bianchi, E.L.; Polinder, H.; Bandyopadhyay, S. Energy consumption of electric powertrain architectures: A comparative study. In Proceedings of the 19th European Conference on Power Electronics and Applications (EPE’17 ECCE Europe), Warsaw, Poland, 11–14 September 2017; p. 1. [Google Scholar]
  23. Li, K.; Bouscayrol, A.; Han, S.; Cui, S. Comparisons of electric vehicles using modular cascade machines system and classical single drive electric machine. IEEE Trans. Veh. Technol. 2018, 67, 354–361. [Google Scholar] [CrossRef]
  24. Wang, B.; Hung, D.L.S.; Zhong, J.; Teh, K.Y. Energy Consumption Analysis of Different Bev Powertrain Topologies by Design Optimization. Int. J. Automot. Technol. 2018, 19, 907–914. [Google Scholar] [CrossRef]
  25. Domingues-Olavarría, G.; Márquez-Fernández, F.J.; Fyhr, P.; Reinap, A.; Andersson, M.; Alaküla, M. Optimization of Electric Powertrains Based on Scalable Cost and Performance Models. IEEE Trans. Ind. Appl. 2019, 55, 751–764. [Google Scholar] [CrossRef]
  26. Silvas, E.; Hofman, T.; Serebrenik, A.; Steinbuch, M. Functional and cost-based automatic generator for hybrid vehicles topologies. IEEE/ASME Trans. Mechatron. 2015, 20, 1561–1572. [Google Scholar] [CrossRef] [Green Version]
  27. Ebbesen, S.; Kiwitz, P.; Guzzella, L. A generic particle swarm optimization Matlab function. In Proceedings of the 2012 American Control Conference (ACC), Montreal, QC, Canada, 27–29 June 2012; pp. 1519–1524. [Google Scholar]
  28. Guzzella, L.; Sciaretta, A. Vehicle Propulsion Systems, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  29. Steinhilper, W.; Sauer, B. Structural Elements of Mechanical Engineering 2: Fundamentals of Machine Elements for Powertrains, 7th ed.; Springer: Berlin/Heidelberg, Germany, 2012. (In German) [Google Scholar]
  30. Naunheimer, H.; Bertsche, B.; Ryborz, J.; Novak, W. Automotive Transmissions: Fundamentals, Selection, Design and Application, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  31. Staunton, R.H.; Ayers, C.W.; Marlino, L.; Chiasson, J.; Burress, B. Evaluation of 2004 Toyota Prius Hybrid Electric Drive System; Technical Report ONRL/TM-2006/423; Oak Ridge National Lab: Oak Ridge, TN, USA, 2006. [Google Scholar]
  32. Argonne National Laboratory. Autonomie Vehicle System Simulation Tool. Available online: https://www.anl.gov/es/autonomie-vehicle-system-simulation-tool (accessed on 20 July 2018).
  33. Thielmann, A.; Sauer, A.; Wietschel. Gesamt-Roadmap fur die Elektromobilitat 2030; Technical Report; Fraunhofer ISI: Karlsruhe, Germany, 2015. [Google Scholar]
  34. Den Boer, E.; Aarnink, S.; Kleiner, F.; Pagenkopf, J. Zero Emissions Trucks: An Overview of State-of-the-Art Technologies and Their Potential; Technical Report 13.4841.21; CE Delft: Delft, The Netherlands, 2013. [Google Scholar]
  35. Summerton, P.; Harrison, P. An Economic Assessment of Low Carbon Vehicles; Technical Report; Cambridge Econometrics and Ricardo-AEA: Cambridge, UK, 2013. [Google Scholar]
  36. Pourabdollah, M.; Egardt, B.; Murgovski, N.; Grauers, A. Effect of driving, charging, and pricing scenarios on optimal component sizing of a PHEV. Control Eng. Pract. 2017, 61, 217–228. [Google Scholar] [CrossRef]
  37. Silvas, E.; Bergshoeff, E.; Hofman, T.; Steinbuch, M. Comparison of bi-level optimization frameworks for sizing and control of a hybrid electric vehicle. In Proceedings of the 2014 IEEE Vehicle Power and Propulsion Conference (VPPC), Coimbra, Portugal, 27–30 October 2014; pp. 1–6. [Google Scholar]
  38. Eurostat. Electricity Price Statistics. Available online: https://ec.europa.eu/eurostat/statistics-explained/index.php?title=Electricity_price_statistics (accessed on 13 August 2018).
  39. Fontaras, G.; Rexeis, M.; Dilara, P.; Hausberger, S.; Anagnostopoulos, K. The Development of a Simulation Tool for Monitoring Heavy-Duty Vehicle CO2 Emissions and Fuel Consumption in Europe; Technical Report No. 2013-24-0150; SAE: Warrendale, PA, USA, 2013. [Google Scholar]
Figure 1. Graphical layout of the powertrain system-level design problem and its main design sub-levels.
Figure 1. Graphical layout of the powertrain system-level design problem and its main design sub-levels.
Energies 13 02434 g001
Figure 2. Set of eight different feasible powertrain topologies studied, classified into distributed and central drive systems. The first subscript (i) indicates the topology name; the second subscript after the dot ( j l ) indicates the total number of gears in the gearbox. For each topology, different gearbox types are considered, resulting in a set of in total 44 powertrain configurations studied. W = driven wheel, EMk = electric machine w.power electronics, TR1 = transmission unit, BA = battery pack, FD = final drive transmission.
Figure 2. Set of eight different feasible powertrain topologies studied, classified into distributed and central drive systems. The first subscript (i) indicates the topology name; the second subscript after the dot ( j l ) indicates the total number of gears in the gearbox. For each topology, different gearbox types are considered, resulting in a set of in total 44 powertrain configurations studied. W = driven wheel, EMk = electric machine w.power electronics, TR1 = transmission unit, BA = battery pack, FD = final drive transmission.
Energies 13 02434 g002
Figure 3. Efficiency map of PMSMused for scaling the electric machine efficiency with the peak power [31]. The map includes both machine and inverter efficiencies. Peak efficiency η ¯ m = 0.937, peak power P m k , base = 50 kW, maximum torque τ ¯ m = 375 Nm, maximum speed ω ¯ m = 629 rad/s, and base speed ω m , b = 116 rad/s.
Figure 3. Efficiency map of PMSMused for scaling the electric machine efficiency with the peak power [31]. The map includes both machine and inverter efficiencies. Peak efficiency η ¯ m = 0.937, peak power P m k , base = 50 kW, maximum torque τ ¯ m = 375 Nm, maximum speed ω ¯ m = 629 rad/s, and base speed ω m , b = 116 rad/s.
Energies 13 02434 g003
Figure 4. Data on open-circuit voltage ( U oc , cl ) and internal resistance ( R cl ) for the base battery cell used [32].
Figure 4. Data on open-circuit voltage ( U oc , cl ) and internal resistance ( R cl ) for the base battery cell used [32].
Energies 13 02434 g004
Figure 5. Drive cycle velocity and height profile for the VECTO regional long-haul drive cycle [39].
Figure 5. Drive cycle velocity and height profile for the VECTO regional long-haul drive cycle [39].
Energies 13 02434 g005
Figure 6. Optimal total-cost-of-ownership, C TCO , for each configuration shown in Figure 2. The different markers indicate the gearbox type present in the optimized topology. The dark gray arrow indicates the improvement from the powertrain configuration with the highest to the one with the lowest total-cost-of-ownership.
Figure 6. Optimal total-cost-of-ownership, C TCO , for each configuration shown in Figure 2. The different markers indicate the gearbox type present in the optimized topology. The dark gray arrow indicates the improvement from the powertrain configuration with the highest to the one with the lowest total-cost-of-ownership.
Energies 13 02434 g006
Figure 7. The sum of all electric machines mechanical peak powers P ^ m = k = 1 n m P m k , plotted versus the optimal battery size E b for each powertrain configuration optimized. The markers indicate the gearbox type present in the topology for that specific data point. The data points belonging to the same topology are connected with lines. Gray colored lines indicate distributed topologies, and black lines indicate central drive topologies.
Figure 7. The sum of all electric machines mechanical peak powers P ^ m = k = 1 n m P m k , plotted versus the optimal battery size E b for each powertrain configuration optimized. The markers indicate the gearbox type present in the topology for that specific data point. The data points belonging to the same topology are connected with lines. Gray colored lines indicate distributed topologies, and black lines indicate central drive topologies.
Energies 13 02434 g007
Figure 8. The optimal total-cost-of-ownership ( C TCO ) versus the number of gears in the gearbox. The results are grouped with on the left the distributed configurations and on the right the central drive configurations. The results belonging to the same topology are connected with lines. The red dotted lines and the red number indicate the average total-cost-of-ownership for each of the two groups.
Figure 8. The optimal total-cost-of-ownership ( C TCO ) versus the number of gears in the gearbox. The results are grouped with on the left the distributed configurations and on the right the central drive configurations. The results belonging to the same topology are connected with lines. The red dotted lines and the red number indicate the average total-cost-of-ownership for each of the two groups.
Energies 13 02434 g008
Figure 9. The vehicle cost C v plotted versus the electricity running cost C e for each optimized powertrain configuration. The data points belonging to the same topology are connected with lines, and the different markers indicate the gearbox type present. Gray colored lines are used for distributed topologies and black lines for central drive topologies. “Combinations” refers to the configurations of T C 4 equipped with two different types of gearboxes, for example a one-speed in combination with a two-speed.
Figure 9. The vehicle cost C v plotted versus the electricity running cost C e for each optimized powertrain configuration. The data points belonging to the same topology are connected with lines, and the different markers indicate the gearbox type present. Gray colored lines are used for distributed topologies and black lines for central drive topologies. “Combinations” refers to the configurations of T C 4 equipped with two different types of gearboxes, for example a one-speed in combination with a two-speed.
Energies 13 02434 g009
Figure 10. Electric machine sizing parameters plotted against one another. The markers indicate the gearbox type present in the topology for that specific data point. The data points belonging to the same topology are connected with lines. Gray colored lines indicated distributed topologies, and black lines indicate central drive topologies. “Combinations” refers to the configurations of T C 4 equipped with two different types of gearboxes, for example a one-speed in combination with a two-speed.
Figure 10. Electric machine sizing parameters plotted against one another. The markers indicate the gearbox type present in the topology for that specific data point. The data points belonging to the same topology are connected with lines. Gray colored lines indicated distributed topologies, and black lines indicate central drive topologies. “Combinations” refers to the configurations of T C 4 equipped with two different types of gearboxes, for example a one-speed in combination with a two-speed.
Energies 13 02434 g010
Figure 11. Design choice tree from the perspective of T C 1.1 , showing the impact on the optimal total-cost-of-ownership and optimal sizing parameters for the three topological design choices studied: (top) moving to a distributed drive; (middle) adding an extra gear to the gearbox; (bottom) adding an extra electric machine to the powertrain. The change in sizing parameters due to each choice is shown in the boxes above and below the configurations. In case control variables are present, the use of the gears and electric machine modes for driving the vehicle are mentioned for each configuration.
Figure 11. Design choice tree from the perspective of T C 1.1 , showing the impact on the optimal total-cost-of-ownership and optimal sizing parameters for the three topological design choices studied: (top) moving to a distributed drive; (middle) adding an extra gear to the gearbox; (bottom) adding an extra electric machine to the powertrain. The change in sizing parameters due to each choice is shown in the boxes above and below the configurations. In case control variables are present, the use of the gears and electric machine modes for driving the vehicle are mentioned for each configuration.
Energies 13 02434 g011
Figure 12. Visualization of the work points of the vehicle over the drive cycle for both T D 2.2 and T C 2.2 . The upper plots indicate the use of the electric machines, where EM1 only means that only EM1 supplies power to the wheels and EM2 only that this is only done by EM2. Dual motor indicates that both machines together drive the wheels of the vehicle. The lower two plots indicate the engaged gear for each work point.
Figure 12. Visualization of the work points of the vehicle over the drive cycle for both T D 2.2 and T C 2.2 . The upper plots indicate the use of the electric machines, where EM1 only means that only EM1 supplies power to the wheels and EM2 only that this is only done by EM2. Dual motor indicates that both machines together drive the wheels of the vehicle. The lower two plots indicate the engaged gear for each work point.
Energies 13 02434 g012
Figure 13. Visualization of the workpoints of the electric machines over the drive cycle for both T D 2.2 and T C 2.2 . The workpoints correspond with the workpoints shown in Figure 12, however, here plotted for each electric machine in the powertrain. The upper plots indicate the use of the electric machines, where (EM1 only) means that only EM1 supplies power to the wheels, and (EM2 only) that this is only done by EM2. (Dual motor) indicates that both machines together drive the wheels of the vehicle. The lower two plots indicate the engaged gear for each workpoint.
Figure 13. Visualization of the workpoints of the electric machines over the drive cycle for both T D 2.2 and T C 2.2 . The workpoints correspond with the workpoints shown in Figure 12, however, here plotted for each electric machine in the powertrain. The upper plots indicate the use of the electric machines, where (EM1 only) means that only EM1 supplies power to the wheels, and (EM2 only) that this is only done by EM2. (Dual motor) indicates that both machines together drive the wheels of the vehicle. The lower two plots indicate the engaged gear for each workpoint.
Energies 13 02434 g013
Figure 14. Plot of the optimal energy consumption E v versus the number of gears in the gearbox present for the optimization results with energy consumption as the optimization objective. The results are divided into two figures with, on the left, the distributed and, on the right, the central drive. The red line and number indicate the average energy consumption for each of the two groups.
Figure 14. Plot of the optimal energy consumption E v versus the number of gears in the gearbox present for the optimization results with energy consumption as the optimization objective. The results are divided into two figures with, on the left, the distributed and, on the right, the central drive. The red line and number indicate the average energy consumption for each of the two groups.
Energies 13 02434 g014
Table 1. Powertrain system platform.
Table 1. Powertrain system platform.
Component TypeNo. InstancesDegree
1:EMk41
2:TR122
3:FD13
4:BA11
Table 2. Vehicle parameters.
Table 2. Vehicle parameters.
NameSymbolValueUnit
Wheel radius r w 0.492m
Gravitational acceleration g 9.81m/s2
Rolling friction coefficient c r 0.006-
Aerodynamic drag coefficient c d 0.73-
Frontal areaA9.75m2
Air density ρ 1.225kg/m3
Table 3. Battery pack specs. Battery type: Li-ion.
Table 3. Battery pack specs. Battery type: Li-ion.
NameSymbolValueUnit
Cell capacity Q 0 , cl 41Ah
Nominal cell voltage U oc , cl ( ξ = 0.5 ) 3.62V
Max. cell charge current I ¯ ch 150A
Max. cell discharge current I ¯ di 150A
Battery nominal voltage U oc ( ξ = 0.5 ) 520V
Max. battery state of charge ξ ¯ 1.0-
Min. battery state of charge ξ ̲ 0.2-
Initial battery state of charge ξ ( t 0 ) 1.0-
Table 4. Mass models of powertrain components.
Table 4. Mass models of powertrain components.
ComponentMass Equation (kg)
Tractor mass m ta = 5400
Trailer mass m tr = 7500
Cargo mass m ca = 25,000
Battery m b = 6.7 E b
Inverter m i k = 0.1 ( P m k 1000 )
Electric machine m m k = 0.8 ( P m k 1000 )
Transmission [30] m gl = 1.723 ( r ¯ l τ ¯ m k ) 0.439 j l 0.219
Table 5. Cost models of powertrain components and assumed simulation input parameters.
Table 5. Cost models of powertrain components and assumed simulation input parameters.
ComponentCost ModelUnit
Electric machine cost model C m k = 16 ( P m k 1000 ) euro
Inverter cost model C i k = 15 ( P m k 1000 ) euro
Battery cost model C b = 250 E b euro
Transmission cost model [30] C g l = 16.73 m g l euro
Base vehicle cost C 0 = 75 , 000 euro
Years of ownership y = 4year
Years of economical lifetime y el = 8year
Yearly mileage D y = 150 , 000 km
Electricity cost [38] c e = 0.114 euro/kWh

Share and Cite

MDPI and ACS Style

Verbruggen, F.J.R.; Silvas, E.; Hofman, T. Electric Powertrain Topology Analysis and Design for Heavy-Duty Trucks. Energies 2020, 13, 2434. https://0-doi-org.brum.beds.ac.uk/10.3390/en13102434

AMA Style

Verbruggen FJR, Silvas E, Hofman T. Electric Powertrain Topology Analysis and Design for Heavy-Duty Trucks. Energies. 2020; 13(10):2434. https://0-doi-org.brum.beds.ac.uk/10.3390/en13102434

Chicago/Turabian Style

Verbruggen, Frans J. R., Emilia Silvas, and Theo Hofman. 2020. "Electric Powertrain Topology Analysis and Design for Heavy-Duty Trucks" Energies 13, no. 10: 2434. https://0-doi-org.brum.beds.ac.uk/10.3390/en13102434

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