Week-11 Challenge: Braking

1.) The braking energy formula is as follows : 

$benery=\frac{\mathrm{1/}}{2}\cdot m\cdot (v_i^2-v_f^2)$

 

Where

is the braking energy is Joules (J), $m$ is the mass in kg , $v_i$is the initial velocity of the vehicle in m/s and $v_f$

is the final velocity of the vehicle in m/s. 

 

The UDDS cycle is used for this particular calculation. In this code, first braking scenario has to be identified, and then the formula has to be applied. The braking scenario happens when the velocity is decreases. In the code, first those instances where the velocity decreases is found, and then the braking energy formula is applied, which gives the instantaneous braking energy. All these values are added up to obtain the total energy required for braking in that particular drive cycle. 

The braking energy formula is as follows 

Benergy=12⋅m⋅(v2i−v2f)

Benergy=12⋅m⋅(vi2-vf2)

Where, 

Benergy

Benergyis the braking energy is Joules (J), 

m is the mass in kg , 

viis the initial velocity of the vehicle in m/s and 

vfis the final velocity of the vehicle in m/s. and the other one formula is E = 1/2 * m * (Vi -Vf)^2

 

Kinetic energy is provided by the engine to accelerate the vehicle from a standstill to the desired speed. This energy is dissipated as heat during the braking.

The Drive cycle is the set of data representing the relation between vehicle velocity and time. The drive cycle is prepared in excel as shown below.

If the vehicle is decelerating from its initial velocity(Vi) to final velocity(Vf), then the energy required for braking is given as:

E = 1/2 * m * (Vi -Vf)^2

Code:

 

clc;
clear;
% Input
m = 1500; % mass in kg
d = readtable ('Drive_c.xlsx');
t = d.Time; % time(s)
v = d.Velocity / 3.6 ; % velocity(m/s)
E = zeros(1,length(t));
BE = 0;
for i = 1 : length(v)-1
if v(i)>v(i+1)
E(i) = 0.5 * m * ((v(i))^2-(v(i+1))^2) ;
BE = BE + E(i);
end

end
disp("Breaking energy is " + BE/1000 + " KWh" )
plot(t,v,linewidth=2);
title("Drive cycle");
xlabel("Time(s)");
ylabel("Velocity(m/s)");
figure
plot(t,E,linewidth=2);
title("Breaking energy");
xlabel("Time(s)");
ylabel("Breaking energy");

 

Simulation results:

The Drive cycle is as shown below.

 

The Drive cycle is as shown below.

The below graph shows the breaking energy required.

Breaking energy is 1496.1227 KWh

 

2.) The Electric motor torque-speed characteristics is as follows : 

                                                                                         

 

During an acceleration conditions, while starting, the motor starts with high constant torque, say like the 10Nm in the above motor curve until the base speed, and beyond the base speed, the torque starts decreasing, as the operating point enters the constant power region. On the other hand while braking, the motor is operating at a higher speed (more than the base speed), in such cases the peak braking torque is not available for the motor. In the above example curve, it can be observed that when the operating point is transferred from the first quadrant to the second quadrant, the braking torque is -5Nm. As in this case the peak braking torque is always less than the peak accelerating torque. 

 

This affects the amount of energy regenerated during braking. Since peak braking torque is decreased, the amount of energy regernated also decreses. Other reasons for lesser regenerative energy during braking is the lower limit of charging current that can be given to the battery. This type of regeneration of energy is possible only with the help of bi-directional power converters.

 

The electric and mechanical brakes are coordinated in two different control strategies: 

 

a.)Series control strategy: In this case, the two braking systems work one after the other. Most commonly, the electric brakes start to operate ahead of the mechanical brakes. The mechanical braking system starts operating only when the electric brakes reaches its maximum limit (max brake torque).  

                                                                                       

 b.) Parallel control strategy: In this case, both the electric and mechanical braking work simultaneously to provide the desired braking torque. 

                                                                             

The application of these strategies depends on comfort required and energy recovery required.

Application of Braking Combination :

EDIB (Electric Driven Intelligent Brake) in NISSAN LEAF

• EDIB allows and controls optimization of the regenerative brake and regular friction brake (hydraulic brake).
• The system produces an adequate braking force keeping energy regeneration increased as much as possible.
• At higher speeds, a combination of regenerative and friction braking is applied.
• At lower speeds, regenerative braking is given a higher proportion for maximised energy conversion.

 

3.) The efficiency of the motor is given by : 

η=O.PI.P

 

Where, 

ηis the efficiency 

I.P is the input power and 

O.Pis the output power. 

 

The relation between the input power and the output power is as follows :

I.P=O.P+Losses

 

The motor generally has the following types of losses : 

 

- Copper losses : These losses occurs at the windings in the stator and the rotor, and it depends on the current flowing through the windings and the internal resistance. Since the current and torque are directly proportional, the copper losses can also be represented in terms of the output torque.

 

Copper losses = I2⋅R=Kc⋅T2

 

Iron Losses : This loss is dependent on the magnetic field, and hence this loss is also called as the magnetic loss. This loss depends on the magnetic flux, and hence also depends on the speed of the motor. The constant 

Ki

 

Kidepends on the magnetic field, due to which it is constant for a PMDC motor.

 

Iron Losses = Ki⋅ω. 

 

Windage and Friction losses: These are mechanical losses, caused due to the wind resistance and friction. This is represented by the following equation. 

 

Windage Losses = Kw⋅ω3

 

All the other losses are represented using a constant 

C

Hence the total losses in a motor is given by : 

 

Losses = Kc⋅T2+Ki⋅ω+Kw⋅ω3+C

 

The input power is the electrical power given to the motor, the output power is the mechanical power. This mechanical power is the product of the torque and speed of the motor. These formulas are used to create the matlab code. 

 

The matlab code is as follows : 

 

close all

clear all

clc

 

t = linspace(1,200);

n = linspace(0,1000); 

 

[T,N] = meshgrid(t,n);

Out_P = T.*N;

 

Kc = 0.2;

Ki = 0.008;

Kw = 0.00001; 

C = 20; 

 

CuL = Kc*(T.^2);

IrL = Ki*N;

WiL = Kw*(N.^3);

 

In_P = Out_P + CuL + IrL + WiL + C ; 

eff = Out_P./In_P;

 

eff_cont = [0.7 0.75 0.8 0.85 0.9 0.91 0.92 0.93];

 

figure(1)

grid on 

plot_line = contour(N,T,eff,eff_cont);

xlabel('Speed in rad/sec');

ylabel ('Torque in Nm');

title( 'Electric Motor Efficiency Plot');

clabel(plot_line);

The electric motor efficiency plot is as follows :                            

 

                                                                                       

                                                                                                                     ELECTRIC MOTOR EFFICIENCY PLOT

 The above plot represents the efficiency contours in the torque-speed plot of the motor. The control strategies are formed such that the motor is operated in the higher efficient region. At these torque-speed operating points the losses indicated above are at their lowest.