0

1
Electrical Vehicle Design and Modeling
Erik Schaltz
Aalborg University
Denmark
1. Introduction
Electric vehicles are by many seen as the cars of the future as they are high efficient, produces
no local pollution, are silent, and can be used for power regulation by the grid operator.
However, electric vehicles still have critical issues which need to be solved. The three main
challenges are limited driving range, long charging time, and high cost. The three main
challenges are all related to the battery package of the car. The battery package should both
contain enough energy in order to have a certain driving range and it should also have a
sufficient power capability for the accelerations and decelerations. In order to be able to
estimate the energy consumption of an electric vehicles it is very important to have a proper
model of the vehicle (Gao et al., 2007; Mapelli et al., 2010; Schaltz, 2010). The model of an
electric vehicle is very complex as it contains many different components, e.g., transmission,
electric machine, power electronics, and battery. Each component needs to be modeled
properly in order prevent wrong conclusions. The design or rating of each component is a
difficult task as the parameters of one component affect the power level of another one. There
is therefore a risk that one component is rated inappropriate which might make the vehicle
unnecessary expensive or inefficient. In this chapter a method for designing the power system
of an electric vehicle is presented. The method insures that the requirements due to driving
distance and acceleration is fulfilled.
The focus in this chapter will be on the modeling and design of the power system of a battery
electric vehicle. Less attention will therefore be put on the selection of each component
(electric machines, power electronics, batteries, etc.) of the power system as this is a very big
task in it self. This chapter will therefore concentrate on the methodology of the modeling and
design process. However, the method presented here is also suitable for other architectures
and choice of components.
The chapter is organized as follows: After the introduction Section 2 describes the modeling
of the electric vehicle, Section 3 presents the proposed design method, Section 4 provides a
case study in order to demonstrate the proposed method, and Section 5 gives the conclusion
remarks.

2. Vehicle modeling
2.1 Architecture

Many different architectures of an electric vehicle exist (Chan et al., 2010) as there are many
possibilities, e.g., 1 to 4 electric machines, DC or AC machines, gearbox/no gearbox, high or

www.intechopen.com

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

low battery voltage, one or three phase charging, etc. However, in this chapter the architecture
in Fig. 1 is chosen.
The purpose of the different components in Fig. 1 will here shortly be explained: The traction
power of the wheels is delivered by the three phase electric machine. The torque of the left
and right wheels are provided by a differential with also has a gear ratio in order to fit the high
speed of the electric machine shaft to the lower speed of the wheels. The torque and speed
of the machine are controlled by the inverter which inverts the battery DC voltage to a three
phase AC voltage suitable for the electric machine. When analyzing the energy consumption
of an electric vehicle it is important also to include the losses due to the components which
not are a part of the power chain from the grid to the wheels. These losses are denoted as
auxiliary loss and includes the lighting system, comfort system, safety systems, etc. During
the regenerative braking it is important that the maximum voltage of the battery is not
exceeded. For this reason a braking resistor is introduced. The rectifier rectifies the three
phase voltages and currents of the grid to DC levels and the boost converter makes it possible
to transfer power from the low voltage side of the rectifier to the high voltage side of the
battery.

kDE

/xX- k
X
/xY- k
Y

kTH
-

Tgevkhkgt xTH
/

Dqquv"
eqp/
xgtvgt

kDcv -

kCwz
Cwzk/
nkct{
nqcfu

kDT
Dtcmkpi"
tgukuvqt

Dcvvgt{

kKpx
Kpxgtvgt

xDcv
/

kC
kD
kE

Gngevtke"
ocejkpg

Vtcpuokuukqp

Itkf
/xW- k

u
u

Fig. 1. Architecture of the battery electric vehicle. In the figure the main components of the
vehicles which have an inuence on the energy consumption of the vehicle is shown.
2.2 Force Model

The forces which the electric machine of the vehicle must overcome are the forces due to
gravity, wind, rolling resistance, and inertial effect. These forces can also be seen in Fig. 2
where the forces acting on the vehicle are shown.

fn
fwind

vcar
fI
ft
fg

frr

Fig. 2. Free body diagram of the forces (thick arrows) acting on the car.

www.intechopen.com

33

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

The traction force of a vehicle can be described by the following two equations (Ehsani et al.,
2005):
fn




f t = Mcar v car + Mcar g sin() + sign(vcar ) Mcar g cos() crr
     



fI

fg

f rr

1
+ sign(vcar + vwind ) air Cdrag Afront (vcar + vwind )2


2

(1)

f wind


3.6
crr = 0.01 1 +
vcar ,
100


where f t
fI
f rr
fg
fn
f wind

Mcar
vcar
v car
g = 9.81
air = 1.2041
crr
Cdrag
Afront
vwind

[N ]
[N ]
[N ]
[N ]
[N ]
[N ]
[rad]
[kg]
m/s
[
]

2
m/s


2
 m/s 3
kg/m
[]
[]
 2
m
[m/s]

(2)

Traction force of the vehicle

Inertial force of the vehicle
Rolling resistance force of the wheels
Gravitational force of the vehicle
Normal force of the vehicle
Force due to wind resistance
Angle of the driving surface
Mass of the vehicle
Velocity of the vehicle
Acceleration of the vehicle
Free fall acceleration
Air density of dry air at 20 C
Tire rolling resistance coefficient
Aerodynamic drag coefficient
Front area
Headwind speed

2.3 Auxiliary loads

The main purpose of the battery is to provide power for the wheels. However, a modern car
have also other loads which the battery should supply. These loads are either due to safety,
e.g., light, wipers, horn, etc. and/or comfort, e.g., radio, heating, air conditioning, etc. These
loads are not constant, e.g., the power consumption of the climate system strongly depend on
the surrounding temperature. Even though some average values are suggested which can be
seen in Table 1. From the table it may be understood that the total average power consumption
is pAux = 857 W.
Radio
Heating Ventilation Air Condition (HVAC)
Lights
Total pAux

52 W
489 W
316 W
857 W

Table 1. Average power level of the auxiliary loads of the vehicle. The values are inspired
from (Ehsani et al., 2005; Emadi, 2005; Lukic & Emadi, 2002).

www.intechopen.com

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

2.4 Transmission

From Fig. 1 it can be understood that the torque, angular velocity, and power of the
transmission system are given by the following equations:
t = f t rw
t
2
vcar
=
rw

w =
w

pt = f t vcar ,

(3)
(4)
(5)
(6)

where t
[Nm] Traction torque
w [Nm] Torque of each driving wheel
rw
[m] Wheel radius
w [rad/s] Angular velocity of the wheels
pt
[W] Traction power
It is assumed that the power from the shaft of the electric machine to the two driving wheels
has a constant efficiency of TS = 0.95 (Ehsani et al., 2005). The shaft torque, angular velocity,
and power of the electric machine are therefore

TS Gt , pt < 0
s =
(7)
t
TS G , pt 0

where s [Nm]
s [rad/s]
ps
[W ]
G
[]

s = Gw

(8)

ps = s s ,

(9)

Shaft torque of electric machine

Shaft angular velocity of electric machine
Shaft power of electric machine
Gear ratio of differential

2.5 Electric machine

For propulsion usually the induction machine (IM), permanent magnet synchronous machine
(PMSM), and switched reluctance machine (SRM) are considered. The "best" choice is
like many other components a trade off between, cost, mass, volume, efficiency, reliability,
maintenance, etc. However, due to its high power density and high efficiency the PMSM is
selected. The electric machine is divided into an electric part and mechanic part. The electric
part of the PMSM is modeled in the DQ-frame, i.e.,
vd = R s i d + L d

diq
+ e Ld id + e pm
dt

3
v i + vq i q ,
=
2 dd

vq = R s i q + L q
pEM

www.intechopen.com

did
e Lq iq
dt

(10)
(11)
(12)

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

55

where vd
[V] D-axis voltage
vq
[V] Q-axis voltage
id
[A] D-axis current
iq
[A] Q-axis current
Rs
[ ] Stator phase resistance
Ld
[H] D-axis inductance
Lq
[H] Q-axis inductance
pm [Wb] Permanent magnet ux linkage
e [rad/s] Angular frequency of the stator
pm [Wb] Permanent magnet ux linkage
pEM [W] Electric input power
The mechanical part of the PMSM can be modeled as follows:
ds
+ Bv s + c + s
dt
ps = s s ,
e = Js

(13)
(14)

kgm
Shaft moment of inertia
where Js
e
Electromechanical torque
[Nm]
Coulomb torque
c
[Nm]
Bv [Nms/rad] Viscous friction coefficient
The coupling between the electric and mechanic part is given by



3P
pm iq + Ld Lq id iq
22
P
e = s ,
2
e =

(15)
(16)

where P [] Number of poles

2.6 Inverter

A circuit diagram of the inverter can be seen in Fig. 3. The inverter transmits power between
the electric machine (with phase voltages vA , vB , and vC ) and the battery by turning on and
off the switches QA+ , QA- , QB+ , QB- , QC+ , and QC- . The switches has an on-resistance RQ,Inv .
The diodes in parallel of each switch are creating a path for the motor currents during the
deadtime, i.e., the time where both switches in one branch are non-conducting in order to
avoid a shoot-through.
The average power losses of one switch pQ,Inv and diode pD,Inv in Fig. 3 during one
fundamental period are (Casanellas, 1994):




m
m
1
1
(17)
+ i RQ,Inv Ip2 +
+ i cos(EM ) VQ,th,Inv Ip
pQ,Inv =
8 3
2
8




1
1
m
m
pD,Inv =
i RD,Inv Ip2 +
i cos(EM ) VD,th,Inv Ip
(18)
8 3
2
8
2Vp
mi =
,
(19)
VBat

www.intechopen.com

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

kKpx

SC/
XS.vj.Kpx-/

TF.Kpx

FC/

TS.Kpx

SEXS.vj.Kpx-/
TS.Kpx

XF.vj.Kpx
TF.Kpx

SD/
XS.vj.Kpx-/
TS.Kpx

xC
/

FEXF.vj.Kpx
TF.Kpx

kD

kC

EKpx

TS.Kpx

xDcv

XF.vj.Kpx
TF.Kpx

FDXF.vj.Kpx

kE

FD/
XF.vj.Kpx
TF.Kpx

SE/
XS.vj.Kpx-/
TS.Kpx

xD
/

TS.Kpx

SDXS.vj.Kpx-/

FC-

SCXS.vj.Kpx-/

FE/
XF.vj.Kpx
TF.Kpx
-

xE
/

Fig. 3. Circuit diagram of inverter.

where pQ,Inv
[W] Power loss of one switch
pD,Inv
[W] Power loss of one diode
EM
[rad] Power factor angle
Ip
[A] Peak phase current
Vp
[V] Peak phase voltage
mi
[] Modulation index
VBat
[V] Battery voltage
RQ,Inv
[ ] Inverter switch resistance
RD,Inv
[ ] Inverter diode resistance
VQ,th,Inv [V] Inverter switch threshold voltage
VD,th,Inv [V] Inverter diode threshold voltage
If it is assumed that the threshold voltage drop of the switches and diodes are equal, i.e.,
Vth,Inv = VQ,th,Inv = VD,th,Inv , and that the resistances of the switches and diodes also are
equal, i.e., RInv = RQ,Inv = RD,Inv , the total power loss of the inverter is given by


3
6
PInv,loss = 6 PQ,Inv + PD,Inv = RInv Ip2 + Vth,Inv Ip .
2

(20)

The output power of the inverter is the motor input power pEM . The inverter input power
and efficiency are therefore
pInv = vBat iInv = pEM + pInv,loss

pEM
pInv , pEM 0
Inv =
pInv
pEM , pEM < 0,
where iInv [A] Inverter input current
pInv [W] Inverter input power
Inv [] Inverter efficiency

www.intechopen.com

(21)
(22)

77

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

2.7 Battery

The battery pack is the heart of an electric vehicle. Many different battery types exist, e.g.,
lead-acid, nickel-metal hydride, lithium ion, etc. However, today the lithium ion is the
preferred choice due to its relatively high specific energy and power. In this chapter the
battery model will be based on a Saft VL 37570 lithium ion cell. Its specifications can be
seen in Table 2.
Maximum voltage
Nominal voltage
Minimum voltage
1 h capacity
Nominal 1 h discharge current
Maximum pulse discharge current

VBat,max,cell
VBat,nom,cell
VBat,min,cell
Q1,cell
IBat,1,cell
IBat,max,cell

4.2 V
3.7 V
2.5 V
7 Ah
7A
28 A

Table 2. Data sheet specifications of Saft VL 37570 LiIon battery (Saft, 2010).
2.7.1 Electric model

The battery will only be modeled in steady-state, i.e., the dynamic behavior is not considered.
The electric equivalent circuit diagram can be seen in Fig.4. The battery model consist of an
internal voltage source and two inner resistances used for charging and discharging. The
two diodes are ideal and have only symbolics meaning, i.e., to be able to shift between the
charging and discharging resistances. Discharging currents are treated as positive currents,
i.e., charging currents are then negative.

TDcv.fku.egnn
xDcv.kpv.egnn

kDcv.egnn

TDcv.ejc.egnn

xDcv.egnn
/

Fig. 4. Electric equivalent circuit diagram of a battery cell.

From Fig. 4 the cell voltage is therefore given by

vBat,int,cell RBat,cell,disiBat,cell , iBat,cell 0
vBat,cell =
vBat,int,cell RBat,cell,chaiBat,cell , iBat,cell < 0,

(23)

where vBat,cell
[V] Battery cell voltage
vBat,int,cell [V] Internal battery cell voltage
iBat,cell
[A] Battery cell current
RBat,cell,dis [] Inner battery cell resistance during discharge mode
RBat,cell,cha [] Inner battery cell resistance during charge mode
The inner voltage source and the two resistances in Fig. 4 depend on the depth-of-discharge
of the battery. The battery cell have been modeled by the curves given in the data sheet of the
battery. It turns out that the voltage source and the resistances can be described as 10th order
polynomials, i.e.,
10
9
8
7
6
RBat,cell,dis = a10 DoDBat
+ a9 DoDBat
+ a8 DoDBat
+ a7 DoDBat
+ a6 DoDBat
5
4
3
2
+ a5 DoDBat
+ a4 DoDBat
+ a3 DoDBat
+ a2 DoDBat
+ a1 DoDBat + a0

www.intechopen.com

(24)

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

10
9
8
7
6
vBat,int,cell = b10 DoDBat
+ b9 DoDBat
+ b8 DoDBat
+ b7 DoDBat
+ b6 DoDBat
5
4
3
2
+ b5 DoDBat
+ b4 DoDBat
+ b3 DoDBat
+ b2 DoDBat
+ b1 DoDBat + b0
10
9
8
7
6
+ c9 DoDBat
+ c8 DoDBat
+ c7 DoDBat
+ c6 DoDBat
RBat,cell,cha = c10 DoDBat
5
4
3
2
+ c5 DoDBat
+ c4 DoDBat
+ c3 DoDBat
+ c2 DoDBat
+ c1 DoDBat

where a10
a4
b10
b4
c10
c4

=
=
=
=
=
=

-634.0,
-416.4,
-8848,
784,
2056,
578,

a9
a3
b9
b3
c9
c3

=
=
=
=
=
=

2942.1,
60.5,
40727,
-25,
-9176,
25,

a8
a2
b8
b2
c8
c2

=
=
=
=
=
=

-5790.6,
-4.8,
-79586,
55,
17147,
3,

a7
a1
b7
b1
c7
c1

=
=
=
=
=
=

6297.4,
0.2,
86018,
0,
-17330,
0,

a6
a0
b6
b0
c6
c0

=
=
=
=
=
=

+ c0

(25)

(26)

-4132.1, a5 = 1677.7
0.0
-56135, b5 = -5565
4
10168, c5 = -3415
0

2.7.2 Capacity model

The inner voltage source, charging resistance, and discharge resistance all depend on the
depth-of-discharge. The state-of-charge and depth-of-discharge depend on the integral of the
current drawn or delivered to the battery, i.e.,
DoDBat = DoDBat,ini +
SoCBat = 1 DoDBat

 i
Bat,eq,cell

QBat,1,cell

dt

(27)
(28)

where DoDBat
[] Depth-of-discharge
DoDBat,ini [] Initial depth-of-discharge
SoCBat
[] Battery state-of-charge
iBat,eq,cell [A] Equivalent battery cell current
The equivalent battery cell current depend on the sign and amplitude of the current (Schaltz,
2010). Therefore



iBat,cell k
IBat,1,cell IBat,1,cell
, iBat,cell 0
iBat,eq,cell =
(29)
, iBat,cell < 0
Bat,cha iBat,cell

1
, iBat,cell IBat,1,cell
(30)
k=
1.125 , iBat,cell > IBat,1,cell,
where k
[] Peukert number
Bat,cha = 0.95 [] Charging efficiency
It is seen that the peukert number has two different values depending on the amplitude of the
discharge current. For currents higher than the nominal 1 hour discharge current IBat,1,cell the
capacity is therefore reduced significant.
2.7.3 Simulation results

In order to verify the methods used to calculate the state-of-charge, internal voltage source,
and charging resistance calculations are compared to the data sheet values. The results can be
seen in Fig. 5 where the battery cell voltage is shown for different C-values (1 C is the nominal
discharge current of IBat,1,cell = 7 A, which means that C/2 is equal to 3.5 A). It is seen that
the calculated voltages almost are identical to the data sheet values. It is also noticed that the
voltage is strongly depending on the current level and the delivered Ah, and that the voltage
drops significant when the battery is almost completely discharged.

www.intechopen.com

99

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

Voltage vBat,cell [V]

3.5

C/5 Data sheet

C/5 Calculated
C/2 Data sheet
C/2 Calculated
1C Data sheet
1C Calculated
2C Data sheet
2C Calculated

2.5

Capacity QBat,cell [Ah]

Fig. 5. Data sheet values (Saft, 2010) and calculations of the battery voltage during constant
discharge currents.
2.8 Boost converter

The circuit diagram of the boost converter can be seen in Fig. 6. The losses of the boost

XF.vj.DE

TF.DE

ETH

XS.vj.DE -/

FDE

SDE

xTH

kDE

kN
/

kTH

EDcv xDcv
/

TS.DE

Fig. 6. Electric circuit diagram of the boost converter.

converter are due to the switch resistance RQ,BC and threshold voltage VQ,th,BC and the
diodes resistance RD,BC and threshold voltage VD,th,BC . In order to simplify it is assumed
that the resistances and threshold voltages of the switch QBC and diode DBC are equal, i.e.,
RBC = RQ,BC = RD,RF and Vth,BC = VQ,th,BC = VD,th,BC . The power equations of the boost
converter are therefore given by
PRF = VRF iRF = PBC + PLoss,BC

(31)

PBC = VBat iBC

(32)

PLoss,BC =

www.intechopen.com

2
RBC iRF

+ Vth,BC iRF ,

(33)

10

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

10

where PRF
PBC
PLoss,BC
VRF
Vth,BC
RBC
iRF
iBC

[W]
[W]
[W]
[V ]
[V ]
[]
[ A]
[ A]

Input power of boost converter

Output power of boost converter
Power loss of boost converter
Input voltage of boost converter
Threshold voltage of switch and diode
Resistance of switch and diode
Input current of boost converter
Output current of boost converter

2.9 Rectifier

In order to utilize the three phase voltages of the grid vU , vV , and vW they are rectified by a
rectifier as seen in Fig. 7. In the rectifier the loss is due to the resistance RD,RF and threshold
voltage VD,th,RF .
kTH

Xgjkeng

kX

/xY-

kY

TTH

FW-

FX-

FY-

Xvj.TH

Xvj.TH

kN

Xvj.TH

TTH

/xX-

Xvj.TH

TTH

kW

/xW-

Xvj.TH

Itkf

Xvj.TH

TTH

TTH

TTH

FW/

FX/

FY/

xTH

ETH

Dqquv"
eqpxgtvgt

Fig. 7. Electric circuit diagram of the rectifier.

The average rectified current, voltage, and power are given by (Mohan et al., 2003)

3
iRF = IGrid
2

3 2
VRF =
V 2RRF iRF 2Vth,RF
LL
PRF = VRF iRF = PGrid PRF,loss

3 2
V I
LL RF
2
= 2RRF iRF
+ 2Vth,RF iRF ,

PGrid =
PRF,loss
where IGrid
PGrid
Ploss,RF
RRF
Vth,RF

[ A]
[W ]
[W ]
[]
[ V]

www.intechopen.com

Grid RMS-current
Power of three phase grid
Total loss of the rectifier
Resistance of switch and diode
Threshold voltage of switch and diode

(34)
(35)
(36)

(37)
(38)

11
11

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

2.10 Simulation model

The models of each component of the power system in the electric vehicle have now been
explained. When combining all the sub models a model of the battery electric vehicle is
obtained. In Fig. 8 the implementation in a Matlab/Simulink environment can be seen. The
overall vehicle model includes the model of the forces acting on the vehicle (wind, gravity,
rolling resistance, etc.), and the individual components of the power train, i.e., transmission,
electric machine, inverter, battery, boost converter, rectifier. The wind speed vwind and road
angle have been set to zero for simplicity. The input to the simulation model is a driving
cycle (will be explained in Section 4) and the output of the model is all the currents, voltages,
powers, torques, etc, inside the vehicle.

3. Design method
3.1 Parameter determination

The parameter determination of the components in the vehicle is an iterative process. The
parameters are calculated by using the models given in Section 2 and the outputs of the
Matlab/Simulink model shown in Fig. 8.
3.1.1 Battery

The maximum rectified voltage can be calculated from Equation 35 in no-load mode, i.e.,

3 2
3 2
(39)
VRF,max =
V =
400 V = 540 V.
LL

In order to insure boost operation during charging the rectified voltage of the rectifier should
always be greater than this value. The required number of series connected cells is therefore
NBat,s =

VRF,max
540 V
=
216 cells.
VBat,cell,min
2.5 V

(40)

The number of series connected cells NBat,s is due to the voltage requirement of the battery
pack. However, in order to insure that the battery pack contains sufficient power and energy
it is probably not enough with only one string of series connected cells. The battery pack
will therefore consist of NBat,s series connected cells and NBat,p parallel strings. The number
of parallel strings NBat,p are calculated in an iterative process. The ow chart of the sizing
procedure of the battery electric vehicle can be seen in Fig. 9. In the Initialization-process the
base parameters are defined, e.g., wheel radius and nominal bus voltage, initial power ratings
of each component of the vehicle are given, and the base driving cycle is loaded into the
workspace of Matlab. In the Is the minimum number of parallel strings obtained?-decision
block it is verified if the minimum number of parallel strings that fulfills both the energy and
power requirements of the battery have been reached. If not a Simulation routine-process is
executed. This process are executed several times during the sizing procedure and its ow
chart is therefore shown separately in Fig. 9. This process consist of three sub-processes.
The first sub-process is Design components. In this process the parameters of each
component of the battery electric vehicle are determined, e.g., motor and power electronic
parameters. The next sub-process is the Vehicle simulation-process. In this process the
Simulink-model of the vehicle is executed due to the parameters specified in the previous
sub-process. In the third and last sub-process, i.e., the Calculate the power and energy of each
component-process, the energy and power of each component of the vehicle are calculated.

www.intechopen.com

12

12

w_s [rad/s]

w_s [rad/s]

I_p_peak [A]

f_t [N]

tau_s [Nm]

tau_s [Nm]

p_EM [W]

I_p_peak [A]

i_Inv [A]

p_EM [W]

[t v_car]
From
Workspace

Electric machine

Transmission system

v_Bat [V]

i_Aux [A]

v_Bat [V]

Inverter model

Electric auxiliary loads

v_vehicle [km/h]

Braking resistor
Boost converter model

Battery model
i_Aux [A]
v_Bat [V]

[t d_v_car_dt]

i_BC [A]

i_Inv [A]

d_v_vehicle_dt [m/s^2]

i_Bat [A]

From
Workspace1

f_t [N]

SoC_Bat []

SoC_Bat []
i_Bat_cha_max [A]

i_Bat_cha_max [A]

i_RF [A]
v_RF [V]

0
Constant

i_Bat [A]

i_BC [A]
i_Bat_cha_max [A]

v_wind [m/s]

v_Bat [V]

alpha [rad]

Vehicle model

i_RF [A]

v_RF [V]

Rectifier model

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

Fig. 8. Matlab/Simulink implementation of the battery electric vehicle.

www.intechopen.com

v_vehicle [km/h]

13
13

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

The three sub-processes in the Simulation routine-process are executed three times in order
to make sure that parameters converges to the same values for the same input. After the
Simulation routine-process is finish the Calculate number of parallel strings-process is
applied. In this process the number of parallel strings NBat,p is either increased or decreased.
When the minimum possible number of parallel strings that fulfills both the energy and power
requirements of the battery has been found the Simulation routine-process is executed in
order to calculate the grid energy due to the final number of parallel strings.
Dgikp

Ukowncvkqp"tqwvkpg

Kpkvkcnk|cvkqp

hqt"k?3<5

Ku"vjg"okpkowo"pwodgt"qh"
rctcnngn"uvtkpiu"qdvckpgfA

[gu

Fgukip"eqorqpgpvu

Pq
Ukowncvkqp"tqwvkpg

Xgjkeng"ukowncvkqp
Ukowncvkqp"tqwvkpg
Ecnewncvg"vjg"rqygt"cpf"
gpgti{"qh"gcej"eqorqpgpv

Ecnewncvg"pwodgt"qh"
rctcnngn"uvtkpiu

Gpf

Fig. 9. Sizing procedure of the battery electric vehicle.

In principle all the energy of a battery could be used for the traction. However, in order to
prolong the lifetime of the battery it is usually recommended not to charge it to more than
90 % of its rated capacity and not to discharge it below SoCBat,min = 20 %, i.e., only 70 %
of the available energy is therefore utilized. In Fig. 10 it can be seen how the Calculate
number of parallel strings-process finds the minimum number of parallel strings NBat,p
that fulfills both the energy and power requirements. This process is a part of the sizing
procedure shown in Fig. 9. In Fig. 10(a) the minimum state-of-charge min(SoCBat) is shown
and in Fig. 10(b) the maximum battery cell discharge current max(iBat,cell ) is shown. From
the figure it is understood that the first iteration is for NBat,p = 10. However, both the
minimum state-of-charge and maximum discharge current are satisfying their limits, i.e.,
SoCBat,min = 0.2 and IBat,max,cell = 28 A, respectively. Therefore the number of parallel strings
is reduced to NBat,p = 3 for iteration number two. However, now the state-of-charge limit is
exceeded and therefore the number of parallel strings is increased to NBat,p = 8 for iteration
three. This process continuous until iteration number six where the number of parallel strings
settles to NBat,p = 6, as this is the minimum number of parallel strings which ensures that
both the state-of-charge and maximum current requirements are fulfilled.

www.intechopen.com

14

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

14

Minimum state-of-charge min (SoCBat ) []

0.6

1
0.5

0.4
0.3

(a)

4
6

0.2

SoCBat,min

0.1
5
0

2
3

30

5
6
7
8
9
Number of parallel strings NBat,p []


Maximum cell discharge current max iBat,cell [A]

IBat,max,cell

25
20

10

2
(b)

15
5
10

4
6
3

5
3

5
6
7
8
Number of parallel strings NBat,p []

1
9

10

Fig. 10. Number of parallel strings NBat,p due to the Calculate number of parallel
strings-process in Fig. 9. The numbers in the green and yellow boxes indicate the iteration
number of the design procedure. The yellow boxes are the first and last iteration number. (a)
Minimum state-of-charge SoCBat,min. The red dashed horizontal lines indicates the minimum
allowed state-of-charge. (b) Maximum cell discharge current max(iBat,cell ). The dashed red
horizontal line indicates the maximum allowed discharge current.
3.1.2 Electric machine

In order to design the machine design constraints from UQM Technologies (UQM, 2010) are
applied. The machine from UQM Technologies is a brushless permanent magnet synchronous
machine with the specifications in Table 3.
The phase angle between the voltage and current is not specified, but is assumed to be
EM,nom = 0.55 rad, which corresponds to a power factor of cos(EM,nom ) = 0.85. The
shaft angular velocity at (maximum power, maximum torque), (maximum power, continuous

www.intechopen.com

15
15

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

Continuous shaft power

45 kW
Peak shaft power
Ps,max,UQM
75 kW
Maximum speed
ns,max 8000 rpm
Shaft speed at peak torque and power
ns,corner 3000 rpm
Number of poles (Jensen et al., 2009)
P
18
Continuous shaft torque
s,cont,UQM 150 Nm
Peak shaft torque
s,max,UQM 240 Nm
Efficiency at maximum shaft power and continuous torque
EM,a
94 %
Efficiency at maximum shaft power and maximum speed
EM,b
90 %
Table 3. Specifications of UQM PowerPhase 75 (UQM, 2010) from UQM Technologies.
torque), and maximum speed are therefore
ns,corner
2 = 314.2 rad/s
s,corner =
60 s/min
Ps,max,UQM
s,nom =
= 500 rad/s
s,cont,UQM
ns,max
s,max =
2 = 837.8 rad/s.
60 s/min

(41)
(42)
(43)

The relationship between the maximum and continuous shaft torque can be defined as
aEM =

s,cont,UQM
= 1.6.
s,max,UQM

(44)

It is only possible to have the maximum shaft torque of a PMSM as long as the product of
the shaft torque and angular velocity is below the maximum shaft power due to the voltage
induced by the permanent magnet. For a PMSM with a given shaft peak torque s,max the
speed-torque contour can be written as
Ps,max = s,max s,corner

s,max , s s,corner
s,limit =
Ps,max
, s > s,corner.
s

(45)
(46)

The peak shaft torque s,max is selected in such a way that the (s , s )-output from
the Matlab/Simulink simulation is below the shaft torque contour s,limit calculated by
Equation 46.
By trial-and-error-method it turns out that if the coulomb torque and viscous friction are
responsible for 2 % and 6 %, respectively, of the power loss at maximum speed and power,
the maximum efficiency is located around the nominal point of operation. Therefore
0.02 Ps,max
s EM,b
0.06 Ps,max
Bv = 2
.
s EM,b
c =

(47)
(48)

The nominal electro mechanical torque is:

www.intechopen.com

s,max
aEM
= c + Bv s,nom + s,cont.

s,cont =

(49)

e,cont

(50)

16

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

16

The machine will be designed at nominal speed s,nom, maximum power Ps,max , and
minimum bus voltage VBus,min . The speed is approximately proportional to the terminal
voltage. At the minimum bus voltage the machine should be able run at maximum speed
with a modulation index mi = 1. Because the machine is designed at nominal speed, but at
the minimum battery voltage, the modulation index is
mi,nom =

s,nom
= 0.3581.
s,max

(51)

The voltages of the machine are therefore:

VBat,min
Vp,nom = mi,nom
2

Vd,nom = Vp,nom sin(EM,nom )

Vq,nom = Vp,nom cos(EM,nom ).

(52)
(53)
(54)

Several control properties of the PMSM can be applied. Due to its simple implementation the
Id = 0 property is selected even though the reluctance then cannot be utilized. Therefore,
when using Id = 0 control the machine parameters can be calculated as follows:
Ps,max
EM,a

(55)

Iq,cont =

2 PEM,max
3 Vq,nom

(56)

pm =

2 2 e,cont
3 P Iq,cont

(57)

PEM,max =

P
= 4500 rad/s
2
Vd,nom
Lq =
e,nom Iq,cont

e,nom = s,nom

Rs =

Vq,nom e,nom pm
.
Iq,cont

(58)
(59)
(60)

The efficiency of the machine for different torque-speed characteristics can be seen in Fig. 11.
It is seen that the efficiency is highest at continuous torque s,cont and nominal speed ns,nom . A
common mistake in electric vehicle modeling is to assume a fixed efficiency of the components
and it can be understood from Fig. 11 that wrong conclusions therefore can be made if the
electric machines not is operating in a sufficient point of operation. The corner speed ns,corner,
nominal speed ns,nom, maximum speed ns,max , continuous torque s,cont, peak torque s,max ,
and the torque contour s,limit are also shown in the figure.
3.1.3 Transmission

The maximum speed of the electric machine is ns,max = 8000 rpm. The required gear ratio of
the differential is therefore
G=

www.intechopen.com

ns,max 2 rw
.
Vcar,max 60 1.1

(61)

17
17

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

ns,corner

ns,nom

ns,max
95 % EM
93 % EM
90 % EM
85 % EM
80 % EM
75 % EM
70 % EM
s,limit

s,max

90
80
70

Shaft torque s [Nm]

60

> 93 %
> 90 %
> 85 %
> 80 %
> 75 %
> 70 %

s,cont

50
40
30
20
10
0

1000

2000

3000

4000

5000

6000

7000

8000

Shaft speed ns [rpm]

Fig. 11. Efficiency map of the electric machine.
From Equation 61 it is seen that the differential is designed with a 10 % buffer of the maximum
speed of the car.
3.1.4 Rectifier

It is expected that most of the charging of the vehicle will take place at private homes, where
the maximum RMS-current is IGrid,max = 16 A. The maximum grid power and rectifier current
are therefore

3 3
= 10.6 kW
(62)
PGrid,max =
V I
LL Grid,max

3
= 19.6 A.
(63)
IRF,max =
I
2 Grid,max

www.intechopen.com

18

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

18

It is assumed that the rectifier has an efficiency of RF,nom = 0.98 at maximum grid power.
The switch-on resistance can therefore be calculated from Equation 36:


2
PRF,max = RF,nom PGrid,max = PGrid,max 2RRF IRF,max
+ 2Vth,RF IRF,max
(64)

RRF =

PGrid,max (1 RF,nom ) 2Vth,RF IRF,max

= 199 m.
2
2IRF,max

(65)

3.1.5 Boost converter

It is assumed that the boost converter has efficiency BC,nom = 0.98 at maximum power. The
maximum power of the boost converter is therefore
PBC,max = PRF,max BC,nom .

(66)

The threshold voltage is Vth,BC = 1.5 V. From Equation 31 the resistance of the boost converter
is therefore
2
PRF,max = BC,nom PRF,max + RBC IRF,max
+ Vth,BC IRF,max


RBC =

PRF,max (1 BC,rat ) Vth,BC IRF,max

= 475 .
2
IRF,max

(67)

(68)

It might be noticed that the values of the resistances of the rectifier and boost converter are
higher than one would expect. However, this is because it is assumed that all the loss is due
to the threshold voltage and resistance of the switches and diodes.
3.1.6 Inverter

The inverter is also designed at the maximum power of the electric machine and at minimum
battery voltage. However, the inverter is designed at the peak shaft torque instead of at the
maximum speed. The inverter is assumed to have an efficiency of Inv,nom = 0.98 at this point
of operation. The loss PInv,loss,max and resistance RInv of the inverter are therefore
PInv,loss,max =

1 Inv,nom
PEM,max
Inv,nom

e,max = c + Bv s,corner + s,max

Iq,max =

RInv =

2 2 e,max
3 P pm

(70)
(71)

6
Vth,Inv Iq,max
.
3 2
2 Iq,max

PInv,loss,max

The threshold voltage is assumed to be Vth,Inv = 1 V.

www.intechopen.com

(69)

(72)

19
19

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

3.2 Battery charging control

During the charging of the battery, i.e., both due to the regenerative braking and the grid, it is
very important that the maximum battery charging current and voltage not are exceeded. The
maximum allowed cell charging current can be calculated from the inner and outer voltage of
the battery cell, i.e.,
V
vBat,int,cell VBat,max,cell vBat,int,cell
Bat,max,cell
,
IBat,1,cell
RBat,cell,cha
RBat,cell,cha
iBat,cell,cha,max =
(73)
V
vBat,int,cell
I
>
I
.
, Bat,max,cell
Bat,1,cell
Bat,1,cell
RBat,cell,cha

In Equation 73 it is insured that neither the maximum allowed voltage or current are exceeded.
The battery pack consist of NBat,s series connected cells and NBat,p parallel connected strings.
The total voltage and current of the battery pack can therefore be calculated as
vBat = NBat,s vBat,cell

(74)

iBat = NBat,p iBat,cell

(75)

iBat,cha,max = NBat,p iBat,cell,cha,max.

(76)

During the charging of the battery the battery cell voltage vBat,cell should not exceed
VBat,max,cell = 4.2 V and the maximum cell charging current should not be higher than
IBat,1,cell = 7 A (Saft, 2010). In order to charge the battery as fast as possible either the
maximum voltage or maximum current should be applied to the battery. The requested
battery charging current, i.e., the output current of the boost converter iBC , is therefore

iBC
= iBat,cha,max ,

(77)

which means that the requested output power of the boost converter is

pBC
.
= vBat iBC

(78)

The requested charging current insures that neither the maximum allowed voltage or current
are exceeded. However, for a big battery pack the required charging power might be so high
that a special charging station is necessary.
The requested input current of the boost converter, i.e., the rectifier current iRF , can be
calculated by Equation 31 and 78:
2

 

Vth,BC vRF 4RBC pBC

V

RF
th,BC

.
(79)
iRF =
2RBC
The grid RMS-current can therefore from Equation 34 be calculated as



< I
23 iRF
, 23 iRF
Grid,max
IGrid =


IGrid,max , 2 i IGrid,max .
3 RF

www.intechopen.com

(80)

20

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

20

Glider mass
Wheel radius
Front area
Aerodynamic drag coefficient

Mglider 670 kg
rw
0.2785 m
Afront
1.68 m2
Cdrag
0.3

Table 4. Parameters of the vehicle used for the case study.

Thereby it is ensured that the maximum RMS grid current is not exceeded. The actual values
can therefore be obtained by calculating backwards, i.e.,

3
iRF =
(81)
I
2 Grid
pRF = vRF iRF
(82)
2
pBC = pRF RBC iRF
Vth,BC iRF
pBC
iBC =
.
vBat

(83)
(84)

4. Case study
4.1 Driving cycle

When different cars are compared in terms of energy consumption a standard driving cycle is
used. An often used driving cycle is the New European Driving Cycle (NEDC) as this driving
cycle contains both city driving with several start-and-stops and motorway driving, i.e., it is a
good representation of a realistic driving environment. The NEDC has a maximum speed of
120 km/h, an average speed of 33.2 km/h, a duration of 1184 s, and a length of 10.9 km. The
NEDC profile can be seen in Fig. 12. The input to the simulation will be the NEDC repeated 14
times as this should provide a driving distance of 153 km which is assumed to be an acceptable
driving distance.
4.2 Vehicle parameters

The energy consumption of a given vehicle depend on the physical dimensions and total mass
of the vehicle. For this case study the parameters in Table 4 are used. The glider mass is the
mass of the vehicle without motor, battery, power electronics, etc. It might be understood
from the parameters in Table 4 that it is a rather small vehicle, i.e., similar to a Citron C1.
4.3 Results

In Fig. 13 the battery state-of-charge, current, voltage, and the power of the grid and battery
can be seen. It is understood from Fig. 13(a) that the battery is designed due to its energy
requirement rather than the power requirement as the state-of-charge reaches the minimum
allowed value of SoCBat,min = 0.2. In Fig. 13(b) and (c) the battery current and voltage are
shown, respectively. It is seen that when the current becomes higher the voltage becomes
lower as the power should be the same. In Fig. 13(d) the battery and grid power are shown.
It is seen that the charging of the battery is limited by the maximum allowed grid power
PGrid,max . After approximately two hours the battery reaches the maximum voltage, and it is
therefore seen that the battery then is charged under constant-voltage approach, which means
that the battery current and power and grid power slowly are decreased until the battery
reaches its initial state-of-charge value.

www.intechopen.com

21
21

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

120

Speed [km/h]

100

80

60

40

20

200

400

600

800

1000

Time [s]
Fig. 12. New European Driving Cycle (NEDC). This driving cycle will be repeated 14 times
and thereby serving as the input profile of the Matlab/Simulink simulation model.
Due to the minimum battery pack voltage requirement NBat,s = 216 series connected battery
cells are required. The chosen vehicle is designed to be able to handle 14 repetitions of the
NEDC. From Fig. 10 it is understood that NBat,p = 5 parallel strings are demanded in order to
fulfill this requirement. This means that the battery pack has a capacity of
EBat =

VBat,nom,cell NBat,s QBat,1,cell NBat,s

= 28.0 kWh.
1000 Wh/kWh

(85)

The energy distribution of the vehicle can be seen in Fig. 14. During the 14 NEDC repetitions
Et = 11.2 kWh is delivered to the surface between the driving wheels and the road, but
EGrid = 22.7 kWh charging energy is taken from the grid. This means that only 49 % of the
charging energy from the grid is used for the traction and that the grid energy consumption is
148.3 Wh/km. The rest of the energy is lost in the path between the wheels and the grid. The
auxiliary loads are responsible for the biggest energy loss at 17 %. However, it is believed that
this can be reduced significant by using diodes for the light instead of bulbs, and to use heat
pumps for the heating instead of pure resistive heating.
The battery is responsible for the second largest energy waist as 14 % of the grid energy is
lost in the battery. The battery was only designed to be able to handle the energy and power
requirements. However, in order to reduce the loss of the battery it might be beneficial to
oversize the battery as the battery peak currents then will become closer to its nominal current

www.intechopen.com

22

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

22

Battery state-of-charge SoCBat []

0.8
0.6

(a)

0.4
0.2
0

Battery current iBat [A]

40
20

(b)

0
0

Battery voltage vBat [V]

VBat,max

900
800

(c)

700
0

Power [kW]
30
20

Battery pBat
Grid pGrid

PGrid,max

10

(d)

0
10
0

4
Time [h]

Fig. 13. Simulation results of the vehicle with 14 repeated NEDC cycles as input. (a) Battery
state-of-charge. (b) Battery current. (c) Battery voltage. (d) Power of the battery and grid.

www.intechopen.com

23
23

Electrical
Design and Modeling
Electrical VehicleVehicle
Design and Modeling

which will reduce the negative inuence of the peukert phenomena. However, a heavier
battery will also increase the traction power, so the gained reduction in battery loss should be
higher than the increased traction power. A bigger battery will of course also make the vehicle
more expensive, but these issues are left for future work.
ELoss,RF : 2 %
ELoss,Bat : 14 %

Et : 49 %

EAux : 17 %

ELoss,BC : 2 %
ELoss,Inv : 2 %
ELoss,EM : 10 %
ELoss,TS : 4 %
Fig. 14. Energy distribution in the vehicle relative to the grid energy.

5. Conclusion
In this chapter a battery electric vehicle have been modeled and designed. The battery of
the electric vehicle is designed in such a way that both the power and energy requirements
are fulfilled for a given driving cycle. The design procedure is an iterative process as the
power ow inside the vehicle depends on the parameters of each component of the power
system between the grid and driving wheels. The loss of each component in the vehicle
depend on the internal states of the vehicle, i.e., the voltages, currents, speed, torques, and
state-of-charge. These states have been included in the modeling in order to obtain a realistic
energy calculation of the vehicle. A case study with a small vehicle undergoing 14 driving
cycles of type NEDC resulted in a grid energy consumption of 148.3 Wh/km with an efficiency
of 49 % from the grid to the driving wheels. However, a relatively big part of the energy loss
is due to the auxiliary loads, e.g., light, safety systems, comfort systems, etc., and the battery.
For this work the only design constraint of the battery was the voltage limit, and the energy
and power requirements. For future work it is recommended also to include the cost and
overall efficiency as design parameters. It is also suggested to investigate how the loss due to
the auxiliary loads can be reduced.

www.intechopen.com

24

24

Electric Vehicles Modelling and

Simulations
Will-be-set-by-IN-TECH

6. References
Casanellas, F. (1994). Losses in pwm inverters using igbts, IEE Proceedings - Electric Power
Applications 141(5): 235 239.
Chan, C. C., Bouscayrol, A. & Chen, K. (2010). Electric, hybrid, and fuel-cell vehicles:
Architectures and modeling, IEEE Transactions on Vehicular Technology 59(2): 589
598.
Ehsani, M., Gao, Y., Gay, S. E. & Emadi, A. (2005). Modern Electric, Hybrid Electric, and Fuel Cell
Vehicles - Fundamentals, Theory, and Design, first edn, CRC Press LLC.
Emadi, A. (2005). Handbook of Automotive Power Electronics and Motor Drives, first edn, Taylor
& Francis.
Gao, D. W., Mi, C. & Emadi, A. (2007). Modeling and simulation of electric and hybrid
vehicles, Proceedings of the IEEE 95(4): 729 745.
Jensen, K. K., Mortensen, K. A., Jessen, K., Frandsen, T., Runlfsson, G. & Thorsdttir, T.
(2009). Design of spmsm drive system for renault kangoo, Aalborg University .
Lukic, S. & Emadi, A. (2002). Performance analysis of automotive power systems: effects of
power electronic intensive loads and electrically-assisted propulsion systems, Proc.
of IEEE Vehicular Technology Conference (VTC) 3: 1835 1839.
Mapelli, F. L., Tarsitano, D. & Mauri, M. (2010). Plug-in hybrid electric vehicle: Modeling,
prototype realization, and inverter losses reduction analysis, IEEE Transactions on
Industrial Electronics 57(2): 598 607.
Mohan, N., Underland, T. M. & Robbins, W. P. (2003). Power electronics, third edn, John Wiley.
Saft (2010). Saftbatteries. URL: http://www.saftbatteries.com
Schaltz, E. (2010). Design of a Fuel Cell Hybrid Electric Vehicle Drive System, Department of
Energy Technology, Aalborg University.
UQM (2010). Uqm technologies. URL: http://www.uqm.com

www.intechopen.com

Electric Vehicles - Modelling and Simulations

Edited by Dr. Seref Soylu

ISBN 978-953-307-477-1
Hard cover, 466 pages
Publisher InTech

Published online 12, September, 2011

Published in print edition September, 2011

In this book, modeling and simulation of electric vehicles and their components have been emphasized
chapter by chapter with valuable contribution of many researchers who work on both technical and regulatory
sides of the field. Mathematical models for electrical vehicles and their components were introduced and
merged together to make this book a guide for industry, academia and policy makers.

How to reference

In order to correctly reference this scholarly work, feel free to copy and paste the following:
Erik Schaltz (2011). Electrical Vehicle Design and Modeling, Electric Vehicles - Modelling and Simulations, Dr.
Seref Soylu (Ed.), ISBN: 978-953-307-477-1, InTech, Available from:
http://www.intechopen.com/books/electric-vehicles-modelling-and-simulations/electrical-vehicle-design-andmodeling

InTech Europe

University Campus STeP Ri

Slavka Krautzeka 83/A
51000 Rijeka, Croatia
Phone: +385 (51) 770 447
Fax: +385 (51) 686 166
www.intechopen.com

InTech China

Unit 405, Office Block, Hotel Equatorial Shanghai

No.65, Yan An Road (West), Shanghai, 200040, China
Phone: +86-21-62489820
Fax: +86-21-62489821