Batteries, Flywheels and Supercapacitors
73The first step is to find
an equation for the current I from a battery when it is operating at a power
P watts. In general we know that
P =
V ×
IIf we then combine this with the basic equation for the terminal voltage of a battery,
which we have written as Equation (3.1), we get
P =
V ×
I =
(E −
IR) ×
I =
EI −
RI2
This
is a quadratic equation for I . The normal useful solution
6
to this equation is
I =
E −
√
E2
− 4
RP2
R(3.26)
This equation allows us easily to use MATLAB® or similar mathematical software to simulate the constant power discharge of a battery. The MATLAB® script file below shows this done fora lead acid battery. The graph of voltage against time is shown in
Figure 3.21.
0 500 1000 1500 2000 2500 3000 3500 4000
Discharged
Time/Seconds
Constant power discharge of a lead acid battery 105 110 115 120 125 130 135 Battery voltage/Volts
Figure 3.21Graph of voltage against time fora constant power discharge of a lead acid battery at 5000 W. The nominal ratings of the battery are 120 V, 50 Ah
6
As with all quadratics, there are two solutions. The other corresponds to a lunatic way of operating
the battery at a huge current, so large that the internal resistance causes the voltage to drop to a low value, so that the power is achieved with a low voltage and very high current. This is immensely inefficient.
74Electric Vehicle Technology Explained, Second Edition A constant P discharge experiment fora lead acid battery The system has 10 batteries,
each V lead acid, Ah We uses steps, as these are sufficiently small fora constant power discharge. We setup arrays to store the data.
T=(0:10:10000);
% Time goes up to 10,000 ins steps This is 1001 values.
CR=zeros(1,1001); % Charge rem. from bat. Peukert corrected.
DoD=zeros(1,1001); % Depth of discharge, start fully charged.
V=zeros(1,1001);
% Battery voltage, initially set to zero.
NoCells=60;
% 10 of 6 cell (12 V) batteries.
Capacity=50;
% 50
Ah batteries, 10 h rate capacity.
k=1.12;
% Peukert coefficient deltaT = 10;
% Takes steps, OK for constant power.
P = 5000;
% We will drain the battery at 5 kW Calculated values
Rin= (0.022/Capacity)*NoCells;
% Internal res, Equation (3.2)
PeuCap= ((Capacity/10)^k)*10;
% See Equation (3.18)
% Starting
voltage set outside loopE=open_circuit_voltage_LA(0,NoCells);
I = (E - (E*E - (4*Rin*P))^0.5)/(2*Rin);
%Equation (V E - I*Rin;
%Equation (for n=2:1001
E=open_circuit_voltage_LA(DoD(n-1),NoCells); Equation IEEE- (4*Rin*P))^0.5)/(2*Rin);
%Equation (3.26)
CR(n) = CR(n-1) +((deltaT * I^k)/3600);
%Equation (3.23)
DoD(n)= CR(n)/PeuCap;
%Equation (if DoD(n)>1
DoD(n)=1;
end
% We will say that the battery is "dead" if the depth of discharge exceeds 99%.
V(n)=open_circuit_voltage_LA(DoD(n),NoCells) - I*Rin; Equation (if DoD(n)>0.99
V(n)=0;
end end plot(T,V,'b.');
YLABEL('Battery voltage/Volts');
XLABEL('Time/Seconds');
TITLE('Constant power discharge of a lead acid battery');
axis([0 4000 100 When we come to simulate the battery being used in a vehicle, the issue of regenerative braking will arise. Here a certain power is dissipated
into the battery. If we look again at
Figure 3.1, and consider
the situation if the current I is flowing
into the battery, then the equation becomes
Share with your friends: