Aplicação de Modelos Elétricos de Bateria na Predição do ... - Unijuí
Aplicação de Modelos Elétricos de Bateria na Predição do ... - Unijuí
Aplicação de Modelos Elétricos de Bateria na Predição do ... - Unijuí
Transforme seus PDFs em revista digital e aumente sua receita!
Otimize suas revistas digitais para SEO, use backlinks fortes e conteúdo multimídia para aumentar sua visibilidade e receita.
<strong>Aplicação</strong> <strong>de</strong> Mo<strong>de</strong>los <strong>Elétricos</strong> <strong>de</strong> <strong>Bateria</strong> <strong>na</strong><br />
<strong>Predição</strong> <strong>do</strong> Tempo <strong>de</strong> Vida <strong>de</strong> Dispositivos Móveis<br />
Cleber Mateus Duarte Porciuncula<br />
Dissertação <strong>de</strong> Mestra<strong>do</strong> submetida ao Programa <strong>de</strong> Pós-Graduação<br />
em Mo<strong>de</strong>lagem Matemática da Universida<strong>de</strong> Regio<strong>na</strong>l <strong>do</strong> Noroeste<br />
<strong>do</strong> Esta<strong>do</strong> <strong>do</strong> Rio Gran<strong>de</strong> <strong>do</strong> Sul - <strong>Unijuí</strong> - como parte <strong>do</strong>s requisi-<br />
tos necessários para a obtenção <strong>do</strong> Grau <strong>de</strong> Mestre em Mo<strong>de</strong>lagem<br />
Matemática.<br />
Paulo Sérgio Sausen, Dsc.<br />
Orienta<strong>do</strong>r<br />
Airam Teresa Zago Romcy Sausen, Dsc.<br />
Co-Orienta<strong>do</strong>ra<br />
Ijuí, RS, Brasil<br />
c⃝Cleber Mateus Duarte Porciuncula, Abril <strong>de</strong> 2012
<strong>Aplicação</strong> <strong>de</strong> Mo<strong>de</strong>los <strong>Elétricos</strong> <strong>de</strong> <strong>Bateria</strong> <strong>na</strong><br />
<strong>Predição</strong> <strong>do</strong> Tempo <strong>de</strong> Vida <strong>de</strong> Dispositivos Móveis<br />
Cleber Mateus Duarte Porciuncula<br />
Dissertação <strong>de</strong> Mestra<strong>do</strong> apresentada em Abril <strong>de</strong> 2012<br />
Paulo Sérgio Sausen, Dsc.<br />
Orienta<strong>do</strong>r<br />
Airam Teresa Zago Romcy Sausen, Dsc.<br />
Co-Orienta<strong>do</strong>ra<br />
Robinson Figueire<strong>do</strong> De Camargo, D.sc.<br />
Componente da Banca<br />
Sandro Sawicki, D.sc.<br />
Componente da Banca<br />
Ijuí, RS, Brasil, Abril <strong>de</strong> 2012<br />
ii
À minha mãe a qual <strong>de</strong>vo a vida.<br />
Agra<strong>de</strong>cimentos<br />
Aos meus professores orienta<strong>do</strong>res, Paulo Sérgio Sausen e Airam Sausen, pelas orien-<br />
tações, atenção e disponibilida<strong>de</strong>.<br />
Aos <strong>de</strong>mais professores <strong>do</strong> Mestra<strong>do</strong> pelo auxílio e pelos ensi<strong>na</strong>mentos.<br />
Aos colegas <strong>de</strong> turma (especialmente ao André Be<strong>de</strong>n<strong>do</strong>, Alan Oliveira e Marília Boés-<br />
sio) pela amiza<strong>de</strong>, ajuda, reflexões e troca <strong>de</strong> conhecimento.<br />
À colega Geni, pela ajuda em momentos <strong>de</strong> dúvidas.<br />
À UNUJUÍ pela infra-estrutura oferecida.<br />
À CAPES pelo apoio fi<strong>na</strong>nceiro recebi<strong>do</strong>.<br />
i
Resumo<br />
Nos últimos anos, houve um aumento significativo <strong>na</strong> utilização <strong>de</strong> dispositivos eletrôni-<br />
cos móveis tais como: notebooks, câmeras digitais, celulares inteligentes (i.e., smart-<br />
phones), entre outros. A utilização <strong>de</strong>stes dispositivos, <strong>na</strong> maior parte <strong>do</strong> tempo, está<br />
condicio<strong>na</strong>da ao esta<strong>do</strong> <strong>de</strong> carga <strong>de</strong> suas baterias, então surge a necessida<strong>de</strong> <strong>de</strong> predi-<br />
zer por quanto tempo este dispositivo po<strong>de</strong>rá ficar operacio<strong>na</strong>l. Uma forma <strong>de</strong> realizar<br />
esta predição é a partir <strong>do</strong> uso <strong>de</strong> mo<strong>de</strong>los matemáticos que simulam a <strong>de</strong>scarga <strong>de</strong> uma<br />
bateria, e consequentemente, estimam o seu tempo <strong>de</strong> vida.<br />
Neste trabalho é apresentada a avaliação <strong>de</strong> um mo<strong>de</strong>lo elétrico, presente <strong>na</strong> ferramenta<br />
computacio<strong>na</strong>l MatLab/Simulink, <strong>de</strong>nomi<strong>na</strong><strong>do</strong> Battery, que foi escolhi<strong>do</strong> em virtu<strong>de</strong> da<br />
praticida<strong>de</strong> no que se refere a extração <strong>de</strong> seus parâmetros <strong>de</strong> configuração e por consi<strong>de</strong>rar<br />
parte <strong>do</strong>s efeitos não-lineares que ocorrem em um processo real <strong>de</strong> <strong>de</strong>scarga, que influ-<br />
enciam diretamente no tempo <strong>de</strong> vida da bateria. Todas as simulações computacio<strong>na</strong>is<br />
com o mo<strong>de</strong>lo Battery foram realizadas <strong>na</strong> ferramenta computacio<strong>na</strong>l MatLab/Simulink,<br />
a partir <strong>de</strong> um circuito implementa<strong>do</strong> em diagrama <strong>de</strong> blocos, consi<strong>de</strong>ran<strong>do</strong> os parâme-<br />
tros da bateria <strong>de</strong> íon lítio Nokia mo<strong>de</strong>lo BL-5F, e os parâmetros da bateria <strong>de</strong> íon lítio<br />
polímero PL-383562.<br />
A avaliação <strong>do</strong> mo<strong>de</strong>lo Battery foi realizada a partir da comparação <strong>de</strong> seus resulta-<br />
<strong>do</strong>s simula<strong>do</strong>s, com da<strong>do</strong>s experimentais obti<strong>do</strong>s a partir <strong>de</strong> uma plataforma <strong>de</strong> testes<br />
(i.e., testbed), e com os resulta<strong>do</strong>s simula<strong>do</strong>s <strong>de</strong> um mo<strong>de</strong>lo elétrico <strong>de</strong> alta acurácia en-<br />
contra<strong>do</strong> <strong>na</strong> literatura, <strong>de</strong>nomi<strong>na</strong><strong>do</strong> Mo<strong>de</strong>lo para Predizer Runtime e Características V-I<br />
(i.e., tensão e corrente) <strong>de</strong> uma bateria.<br />
A partir da análise <strong>do</strong>s resulta<strong>do</strong>s obti<strong>do</strong>s em simulações compara<strong>do</strong>s com os da<strong>do</strong>s<br />
experimentais da bateria <strong>de</strong> íon lítio Nokia mo<strong>de</strong>lo BL-5F, verificou-se que o menor erro<br />
<strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida apresenta<strong>do</strong> pelo mo<strong>de</strong>lo Battery foi <strong>de</strong> 2,12% em <strong>de</strong>scarga<br />
contínua, e <strong>de</strong> 1,79% em <strong>de</strong>scarga variável, indican<strong>do</strong> resulta<strong>do</strong>s satisfatórios. Na com-<br />
paração realizada entre o mo<strong>de</strong>lo Battery e o mo<strong>de</strong>lo elétrico <strong>de</strong> alta acurácia, o erro <strong>na</strong><br />
predição <strong>do</strong> tempo <strong>de</strong> vida apresenta<strong>do</strong> pelo mo<strong>de</strong>lo Battery foi <strong>de</strong> 0,534%, em <strong>de</strong>scarga<br />
contínua, e o erro <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida <strong>do</strong> mo<strong>de</strong>lo <strong>de</strong> alta acurácia foi <strong>de</strong> 0,395%.<br />
Em <strong>de</strong>scarga variável, o erro <strong>do</strong> mo<strong>de</strong>lo Battery foi <strong>de</strong> 1,621%, e o erro <strong>do</strong> mo<strong>de</strong>lo <strong>de</strong> alta<br />
acurácia foi <strong>de</strong> 0,338%. Portanto, as diferenças <strong>de</strong> erros entre os mo<strong>de</strong>los, <strong>na</strong> predição <strong>do</strong><br />
tempo <strong>de</strong> vida, foram <strong>de</strong> ape<strong>na</strong>s 0,139%, em <strong>de</strong>scarga contínua, e <strong>de</strong> 1,283%, em <strong>de</strong>scarga<br />
variável, indican<strong>do</strong> que os resulta<strong>do</strong>s apresenta<strong>do</strong>s pelo mo<strong>de</strong>lo Battery são satisfatórios,<br />
consi<strong>de</strong>ran<strong>do</strong> que no seu processo <strong>de</strong> extração <strong>de</strong> parâmetros, não foi realiza<strong>do</strong> nenhum<br />
teste experimental com a bateria utilizada, já que os parâmetros foram extraí<strong>do</strong>s direta-<br />
mente <strong>do</strong> datasheet da bateria.<br />
ii
Abstract<br />
In recent years there has been a significant increase in the use of mobile electronics<br />
<strong>de</strong>vices such as laptops, digital cameras, smartphones, among others. The use of these<br />
<strong>de</strong>vices, most of the time, <strong>de</strong>pends of the state of charge of batteries, then, It arises the<br />
need for predicting how long this <strong>de</strong>vice can be operatio<strong>na</strong>l. One way to perform this<br />
prediction is based on the use of mathematical mo<strong>de</strong>ls that simulate the discharge of a<br />
battery, and therefore estimate its lifetime.<br />
In this work is presented the evaluation of an electrical mo<strong>de</strong>l called Battery, avai-<br />
lable in the Matlab/Simulink computatio<strong>na</strong>l tool, this mo<strong>de</strong>l was chosen because of the<br />
practicality at the extraction of configuration parameters and for consi<strong>de</strong>ring the effects<br />
nonlinear occurred in a real process of discharge, which directly influence the lifetime<br />
of the battery. All computer simulations with the Battery mo<strong>de</strong>l were performed in<br />
Matlab/Simulink tool, from a circuit implemented in block diagram, consi<strong>de</strong>ring the pa-<br />
rameters of the Nokia BL-5F lithium ion battery, and the parameters of the PL-383562<br />
polymer lithium ion battery.<br />
The evaluation of the Battery mo<strong>de</strong>l was performed comparing the simulated results<br />
with the experimental data obtained from a testbed, and with simulated results of a<br />
high accuracy electric mo<strong>de</strong>l, called Mo<strong>de</strong>l for Predicting Runtime and V-I (voltage and<br />
current) characteristics of a battery.<br />
From the a<strong>na</strong>lysis of the simulation results compared with the experimental data of<br />
the Nokia BL-5F lithium ion battery, it observed the lowest error in the prediction of life-<br />
time presented by Battery mo<strong>de</strong>l was 2.12% in continuous discharge and 1.79% in variable<br />
discharge, indicating satisfactory results. In the comparison between the Battery mo<strong>de</strong>l<br />
and the high accuracy electric mo<strong>de</strong>l, the error in the prediction of lifetime presented<br />
by Battery mo<strong>de</strong>l was 0.534%, in continuous discharge, and the error in the prediction<br />
of lifetime of the highly accurate mo<strong>de</strong>l was 0.395%. In variable discharge, the error of<br />
the Battery mo<strong>de</strong>l was 1.621%, and the error of the high accuracy mo<strong>de</strong>l was 0.338%.<br />
Therefore, the differences of errors between the mo<strong>de</strong>ls in the prediction of lifetime, were<br />
just 0.139% in continuous discharge, and 1.283% in variable discharge, indicating that<br />
the results presented by Battery mo<strong>de</strong>l are satisfactory, consi<strong>de</strong>ring that in its parame-<br />
ters extraction process, it were not performed experimental tests with the used battery,<br />
because the parameters were extracted directly from the battery datasheet.<br />
iii
Lista <strong>de</strong> Símbolos<br />
i(t) - corrente <strong>de</strong> <strong>de</strong>scarga;<br />
K - Constante <strong>de</strong> polarização (Ah) ou resistência <strong>de</strong> polarização (Ohms);<br />
α - parâmetro que representa a capacida<strong>de</strong> da bateria no mo<strong>de</strong>lo a<strong>na</strong>lítico <strong>de</strong> Rakhmatov-<br />
Vrudhula;<br />
β - parâmetro que representa a não-linearida<strong>de</strong> da bateria no mo<strong>de</strong>lo a<strong>na</strong>lítico <strong>de</strong><br />
Rakhmatov-Vrudhula;<br />
C(x, t) - função concentração <strong>de</strong> espécies eletroativas <strong>do</strong> mo<strong>de</strong>lo a<strong>na</strong>lítico <strong>de</strong> Rakhmatov-<br />
Vrudhula;<br />
L - tempo <strong>de</strong> vida da bateria;<br />
w - comprimento <strong>do</strong> eletrólito da bateria;<br />
C - capacida<strong>de</strong> da bateria;<br />
C ′ - capacida<strong>de</strong> da bateria no início da operação <strong>do</strong> mo<strong>de</strong>lo a<strong>na</strong>lítico Linear;<br />
I - corrente constante <strong>de</strong> <strong>de</strong>scarga;<br />
td - tempo <strong>de</strong> duração da corrente <strong>de</strong> <strong>de</strong>scarga <strong>do</strong> mo<strong>de</strong>lo a<strong>na</strong>lítico Linear;<br />
a - parâmetro relacio<strong>na</strong><strong>do</strong> ao tipo <strong>de</strong> bateria da Lei <strong>de</strong> Peukert;<br />
b - parâmetro relacio<strong>na</strong><strong>do</strong> ao tipo <strong>de</strong> bateria da Lei <strong>de</strong> Peukert;<br />
iv
J(x, t) - fluxo <strong>de</strong> espécies eletroativas;<br />
D - constante <strong>de</strong> difusão;<br />
v - número <strong>de</strong> elétrons envolvi<strong>do</strong>s <strong>na</strong> reação eletroquímica;<br />
F - constante <strong>de</strong> Faraday;<br />
A - área da superfície <strong>do</strong> eletro<strong>do</strong>;<br />
Capacity - capacida<strong>de</strong> nomi<strong>na</strong>l em Ah;<br />
Ik - corrente, on<strong>de</strong> k = 0, ..., n, e k ∈ N;<br />
Ik−1 - corrente, on<strong>de</strong> k = 0, ..., n, e k ∈ N;<br />
tk - tempo, on<strong>de</strong> k = 0, ..., n, e k ∈ N;<br />
tk−1 - tempo, on<strong>de</strong> k = 0, ..., n, e k ∈ N;<br />
DC - corrente contínua;<br />
AC - corrente alter<strong>na</strong>da;<br />
Rseries - resistência em série;<br />
Rtransient - resistência transiente;<br />
Ctransient - capacitância transiente;<br />
RC - re<strong>de</strong> resistiva capacitiva;<br />
SOC - esta<strong>do</strong> <strong>de</strong> carga;<br />
VOC - tensão <strong>de</strong> circuito aberto;<br />
Zac - impedância que mo<strong>de</strong>la o equivalente eletroquímico da bateria;<br />
v
Lseries - indutância em série;<br />
Vlost - queda <strong>de</strong> tensão para perdas inter<strong>na</strong>s;<br />
Ccapacity - capacida<strong>de</strong> total da bateria;<br />
Rself−discharge - resistência <strong>de</strong> auto-<strong>de</strong>scarga;<br />
Rtransient_S - resistência transiente <strong>de</strong> curta duração;<br />
Rtransient_L - resistência transiente <strong>de</strong> longa duração;<br />
Ctransient_S - capacitância transiente <strong>de</strong> curta duração;<br />
Ctransient_L - capacitância transiente <strong>de</strong> longa duração;<br />
f1 - número <strong>de</strong> ciclos;<br />
f2 - temperatura;<br />
Ri - resistência inter<strong>na</strong>;<br />
Ibatt - corrente elétrica da bateria;<br />
Vbatt - tensão elétrica da bateria;<br />
Q - capacida<strong>de</strong> máxima da bateria;<br />
E0 - tensão constante;<br />
i ∗ - corrente dinâmica em baixa frequência;<br />
it - capacida<strong>de</strong> extraída da bateria;<br />
A - tensão exponencial;<br />
vi
B - capacida<strong>de</strong> exponencial;<br />
vii
Lista <strong>de</strong> Tabelas<br />
2.1 Comparação entre as categorias [1]. . . . . . . . . . . . . . . . . . . . . . . . 23<br />
3.1 Parâmetros da bateria <strong>de</strong> níquel metal hidreto Pa<strong>na</strong>sonic HHR650D [22]. . . . . 39<br />
3.2 Parâmetros da bateria <strong>de</strong> íon lítio Nokia BL-5F. . . . . . . . . . . . . . . . . . 42<br />
4.1 Perfis <strong>de</strong> <strong>de</strong>scargas variáveis. . . . . . . . . . . . . . . . . . . . . . . . . . . . 48<br />
4.2 Perfis <strong>de</strong> <strong>de</strong>scargas variáveis, seus tempos, média e <strong>de</strong>svio padrão. . . . . . . . . 49<br />
4.3 Perfis <strong>de</strong> <strong>de</strong>scargas contínuas, seus tempos, média e <strong>de</strong>svio padrão [2]. . . . . . . 50<br />
4.4 Parâmetros utiliza<strong>do</strong>s <strong>na</strong>s simulações <strong>do</strong> mo<strong>de</strong>lo elétrico Battery. . . . . . . . . 51<br />
4.5 Tempos <strong>de</strong> vida experimentais (Tve) e tempos <strong>de</strong> vida simula<strong>do</strong>s (Tvs) pelo<br />
mo<strong>de</strong>lo elétrico Battery com <strong>de</strong>scargas contínuas. . . . . . . . . . . . . . . . . 54<br />
4.6 Tempos <strong>de</strong> vida experimentais (Tve) e tempos <strong>de</strong> vida simula<strong>do</strong>s (Tvs) pelo<br />
mo<strong>de</strong>lo elétrico Battery com <strong>de</strong>scargas variáveis. . . . . . . . . . . . . . . . . . 55<br />
5.1 Parâmetros da bateria <strong>de</strong> íon lítio polímero PL-383562 utiliza<strong>do</strong>s <strong>na</strong>s simulações<br />
<strong>do</strong> mo<strong>de</strong>lo Battery. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63<br />
5.2 Comparação entre o mo<strong>de</strong>lo para Predizer Runtime e Característica V-I e o<br />
mo<strong>de</strong>lo Battery. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65<br />
1
Lista <strong>de</strong> Figuras<br />
2.1 Esquema <strong>de</strong> uma célula eletroquímica [5]. . . . . . . . . . . . . . . . . . . . . 11<br />
2.2 Tensão elétrica em função da <strong>de</strong>scarga [5]. . . . . . . . . . . . . . . . . . . . . 12<br />
2.3 Esquema <strong>do</strong>s diferentes esta<strong>do</strong>s <strong>de</strong> operação da bateria [4]. . . . . . . . . . . . 13<br />
2.4 Capacida<strong>de</strong> da bateria em função <strong>de</strong> sua taxa <strong>de</strong> <strong>de</strong>scarga [5]. . . . . . . . . . . 14<br />
2.5 Densida<strong>de</strong> <strong>de</strong> energia e ano <strong>de</strong> implantação comercial das tecnologias <strong>de</strong> baterias<br />
[7]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16<br />
2.6 Seção <strong>de</strong> perfil <strong>de</strong> <strong>de</strong>scarga constante [5]. . . . . . . . . . . . . . . . . . . . . . 20<br />
2.7 <strong>Bateria</strong> totalmente carregada [4]. . . . . . . . . . . . . . . . . . . . . . . . . . 21<br />
2.8 Mo<strong>de</strong>lo basea<strong>do</strong> em Thevenin [1]. . . . . . . . . . . . . . . . . . . . . . . . . 23<br />
2.9 Mo<strong>de</strong>lo basea<strong>do</strong> em Impedância [21]. . . . . . . . . . . . . . . . . . . . . . . . 24<br />
2.10 Mo<strong>de</strong>lo basea<strong>do</strong> em Runtime [6]. . . . . . . . . . . . . . . . . . . . . . . . . . 24<br />
2.11 Mo<strong>de</strong>lo para Predizer Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong> [1]. . . . . 25<br />
2.12 Curvas características <strong>de</strong> capacida<strong>de</strong> utilizável <strong>de</strong> baterias [1]. . . . . . . . . . . 26<br />
2.13 Tensão em circuito aberto em relação ao SOC [1]. . . . . . . . . . . . . . . . . 27<br />
2.14 Resposta transiente para um evento <strong>de</strong> <strong>de</strong>grau <strong>de</strong> corrente [1]. . . . . . . . . . . 28<br />
2.15 Parâmetros extraí<strong>do</strong>s da bateria <strong>de</strong> íon lítio polímero [1]. . . . . . . . . . . . . 30<br />
2.16 Diagrama esquemático <strong>do</strong> mo<strong>de</strong>lo Battery <strong>do</strong> MatLab/Simulink [22]. . . . . . . 31<br />
3.1 Diagrama em blocos simplifica<strong>do</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico Battery. . . . . . . . . . . 35<br />
3.2 Subsistema Mo<strong>de</strong>l Contiuous. . . . . . . . . . . . . . . . . . . . . . . . . . . . 35<br />
3.3 Curva característica <strong>de</strong> <strong>de</strong>scarga [22]. . . . . . . . . . . . . . . . . . . . . . . 36<br />
3.4 Tempo <strong>de</strong> resposta da bateria <strong>de</strong> 30 segun<strong>do</strong>s [22]. . . . . . . . . . . . . . . . . 37<br />
3.5 Curva característica <strong>de</strong> <strong>de</strong>scarga da bateria <strong>de</strong> Ni-MH Pa<strong>na</strong>sonic HHR650D [22]. 38<br />
3.6 Interface <strong>do</strong> mo<strong>de</strong>lo elétrico Battery <strong>na</strong> forma genérica <strong>de</strong> registro [22]. . . . . . 39<br />
3.7 Interface <strong>do</strong> mo<strong>de</strong>lo elétrico Battery <strong>na</strong> forma não genérica <strong>de</strong> registro [22]. . . . 40<br />
3.8 Curva característica <strong>de</strong> <strong>de</strong>scarga para uma <strong>de</strong>scarga contínua <strong>de</strong> 150 mA <strong>na</strong><br />
bateria <strong>de</strong> íon lítio Nokia BL-5F. . . . . . . . . . . . . . . . . . . . . . . . . . 41<br />
3.9 Matriz N × 3 gerada após simulação com o mo<strong>de</strong>lo elétrico Battery. . . . . . . . 43<br />
4.1 Diagrama <strong>do</strong> testbed [24]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45<br />
1
Lista <strong>de</strong> Figuras 2<br />
4.2 Foto <strong>do</strong> hardware. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45<br />
4.3 Interface <strong>do</strong> software <strong>de</strong> controle <strong>do</strong> testbed. . . . . . . . . . . . . . . . . . . . 46<br />
4.4 Janela <strong>de</strong> configuração <strong>do</strong> teste a ser aplica<strong>do</strong> <strong>na</strong> bateria. . . . . . . . . . . . . 47<br />
4.5 <strong>Bateria</strong>s conectadas à fonte exter<strong>na</strong>. . . . . . . . . . . . . . . . . . . . . . . . 48<br />
4.6 <strong>Bateria</strong>s conectadas ao testbed. . . . . . . . . . . . . . . . . . . . . . . . . . . 48<br />
4.7 Curva real <strong>de</strong> <strong>de</strong>scarga variável para o perfil P1 da bateria <strong>de</strong> íon lítio Nokia BL-5F. 50<br />
4.8 Curva real <strong>de</strong> <strong>de</strong>scarga para o perfil <strong>de</strong> <strong>de</strong>scarga contínua <strong>de</strong> 150 mA <strong>na</strong> bateria<br />
<strong>de</strong> íon lítio Nokia BL-5F. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50<br />
4.9 Programa implementa<strong>do</strong> em diagramas <strong>de</strong> blocos no Simulink para as simulações<br />
com o mo<strong>de</strong>lo elétrico Battery. . . . . . . . . . . . . . . . . . . . . . . . . . . 52<br />
4.10 Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais <strong>na</strong>s <strong>de</strong>scargas con-<br />
tínuas <strong>de</strong> (a) 100 mA, (b) 350 mA, (c) 650 mA e (d) 950 mA, <strong>na</strong> bateria <strong>de</strong> íon<br />
lítio BL-5F. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58<br />
4.11 Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais <strong>do</strong>s perfis <strong>de</strong><br />
<strong>de</strong>scargas variáveis (a) P3, (b) P5 e (c) P7, <strong>na</strong> bateria <strong>de</strong> íon lítio BL-5F. . . . . 59<br />
5.1 Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais para a <strong>de</strong>scarga<br />
contínua <strong>de</strong> 80 mA, <strong>na</strong> bateria <strong>de</strong> íon lítio polímero PL-383562 [1]. . . . . . . . 61<br />
5.2 Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais para a <strong>de</strong>scarga<br />
<strong>de</strong> quatro <strong>de</strong>graus <strong>de</strong> corrente, <strong>na</strong> bateria <strong>de</strong> íon lítio polímero PL-383562 [1]. . . 61<br />
5.3 Curvas reais <strong>de</strong> <strong>de</strong>scarga da bateria <strong>de</strong> íon lítio polímero PL-383562 para dife-<br />
rentes correntes contínuas <strong>de</strong> <strong>de</strong>scarga [25]. . . . . . . . . . . . . . . . . . . . . 62<br />
5.4 Curva simulada para uma <strong>de</strong>scarga contínua <strong>de</strong> 80 mA. . . . . . . . . . . . . . 64<br />
5.5 Curva simulada para uma <strong>de</strong>scarga <strong>de</strong> quatro <strong>de</strong>graus <strong>de</strong> corrente. . . . . . . . . 64
Sumário<br />
1 Apresentação da Dissertação 5<br />
1.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5<br />
1.2 Motivação . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6<br />
1.3 Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7<br />
1.3.1 Objetivo Geral . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7<br />
1.3.2 Objetivos Específicos . . . . . . . . . . . . . . . . . . . . . . . . . . 8<br />
1.4 Contribuições . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8<br />
1.5 Estrutura <strong>do</strong> Documento . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9<br />
2 Revisão Bibliográfica 10<br />
2.1 <strong>Bateria</strong>s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10<br />
2.1.1 Nível <strong>de</strong> Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11<br />
2.1.2 Tempo <strong>de</strong> Vida . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11<br />
2.1.3 Características Não-lineares . . . . . . . . . . . . . . . . . . . . . . 11<br />
2.1.4 Esta<strong>do</strong> <strong>de</strong> Carga da <strong>Bateria</strong> . . . . . . . . . . . . . . . . . . . . . . 14<br />
2.1.5 Capacida<strong>de</strong> . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15<br />
2.2 Tipos <strong>de</strong> <strong>Bateria</strong>s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15<br />
2.3 Mo<strong>de</strong>los <strong>de</strong> <strong>Bateria</strong>s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17<br />
2.3.1 Mo<strong>de</strong>los Eletroquímicos . . . . . . . . . . . . . . . . . . . . . . . . 17<br />
2.3.2 Mo<strong>de</strong>los Estocásticos . . . . . . . . . . . . . . . . . . . . . . . . . . 18<br />
2.3.3 Mo<strong>de</strong>los A<strong>na</strong>líticos . . . . . . . . . . . . . . . . . . . . . . . . . . . 19<br />
2.3.4 Mo<strong>de</strong>los <strong>Elétricos</strong> . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22<br />
3 Mo<strong>de</strong>lo Elétrico Battery 33<br />
3.1 Hipóteses <strong>do</strong> Mo<strong>de</strong>lo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33<br />
3.2 Equação Matemática <strong>do</strong> Mo<strong>de</strong>lo . . . . . . . . . . . . . . . . . . . . . . . . 34<br />
3.3 Mo<strong>de</strong>lo em Diagrama <strong>de</strong> Blocos . . . . . . . . . . . . . . . . . . . . . . . . 34<br />
3.4 Parâmetros <strong>do</strong> Mo<strong>de</strong>lo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35<br />
3.4.1 Extração <strong>de</strong> Parâmetros . . . . . . . . . . . . . . . . . . . . . . . . 37<br />
3
Sumário 4<br />
3.5 Saída <strong>do</strong> Mo<strong>de</strong>lo via Simulink . . . . . . . . . . . . . . . . . . . . . . . . . 41<br />
3.6 Avaliação <strong>do</strong> Mo<strong>de</strong>lo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42<br />
4 Ambiente <strong>de</strong> Simulação 44<br />
4.1 Testbed . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44<br />
4.1.1 Hardware . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45<br />
4.1.2 Software . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45<br />
4.2 Meto<strong>do</strong>logia A<strong>do</strong>tada nos Testes Experimentais . . . . . . . . . . . . . . . 47<br />
4.3 Da<strong>do</strong>s Experimentais Obti<strong>do</strong>s . . . . . . . . . . . . . . . . . . . . . . . . . 49<br />
4.4 Simulações Computacio<strong>na</strong>is . . . . . . . . . . . . . . . . . . . . . . . . . . 51<br />
4.5 Análise <strong>do</strong>s Resulta<strong>do</strong>s <strong>de</strong> Simulações . . . . . . . . . . . . . . . . . . . . . 54<br />
5 Análise Comparativa <strong>do</strong>s Mo<strong>de</strong>los 60<br />
5.1 <strong>Aplicação</strong> <strong>do</strong> Mo<strong>de</strong>lo Elétrico para Predizer Runtime e Características V-I<br />
<strong>de</strong> uma bateria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60<br />
5.2 <strong>Aplicação</strong> <strong>do</strong> Mo<strong>de</strong>lo Elétrico Battery . . . . . . . . . . . . . . . . . . . . . 62<br />
5.3 Comparação entre os Mo<strong>de</strong>los <strong>Elétricos</strong> . . . . . . . . . . . . . . . . . . . . 63<br />
6 Conclusões e Trabalhos Futuros 66<br />
Referências Bibliograficas 69<br />
A Lista das Publicações Relacio<strong>na</strong>das com a Dissertação 72<br />
A.1 Artigo Aceito em Periódico Nacio<strong>na</strong>l . . . . . . . . . . . . . . . . . . . . . 72<br />
A.2 Artigos Publica<strong>do</strong>s em A<strong>na</strong>is . . . . . . . . . . . . . . . . . . . . . . . . . . 72<br />
A.3 Resumos Expandi<strong>do</strong>s Publica<strong>do</strong>s em A<strong>na</strong>is . . . . . . . . . . . . . . . . . . 72<br />
A.4 Resumos Publica<strong>do</strong>s em A<strong>na</strong>is . . . . . . . . . . . . . . . . . . . . . . . . . 73<br />
A.5 Artigo Submeti<strong>do</strong> . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Capítulo 1<br />
Apresentação da Dissertação<br />
1.1 Introdução<br />
As baterias eletroquímicas têm gran<strong>de</strong> importância nos mais diversos tipos <strong>de</strong> sistemas<br />
elétricos, pois a sua energia química armaze<strong>na</strong>da é convertida em energia elétrica, po<strong>de</strong>n<strong>do</strong><br />
ser fornecida para estes sistemas on<strong>de</strong> e quan<strong>do</strong> a energia é necessária. Atualmente,<br />
são inúmeras as áreas <strong>de</strong> aplicação <strong>de</strong>stas baterias, por exemplo, <strong>na</strong> área <strong>de</strong> saú<strong>de</strong>, em<br />
marca passos; <strong>na</strong> área automotiva, em veículos híbri<strong>do</strong>s-elétricos; <strong>na</strong> área <strong>de</strong> dispositivos<br />
eletrônicos móveis, entre outras. Consi<strong>de</strong>ran<strong>do</strong>, nos últimos anos, o crescimento <strong>de</strong>sta<br />
última área, tem ocorri<strong>do</strong> um aumento <strong>na</strong> necessida<strong>de</strong> <strong>de</strong> baterias cada vez menores,<br />
mais leves, e <strong>de</strong> melhor <strong>de</strong>sempenho. Por outro la<strong>do</strong>, observa-se que a popularida<strong>de</strong><br />
<strong>do</strong>s dispositivos móveis como: celulares, câmeras digitais, notebooks, entre outros, tem<br />
impulsio<strong>na</strong><strong>do</strong> a evolução das baterias, mas esta evolução não é ainda capaz <strong>de</strong> aten<strong>de</strong>r<br />
a progressiva <strong>de</strong>manda <strong>de</strong> energia e a limitação <strong>de</strong> tamanho, exigidas pelos dispositivos<br />
móveis atuais [1].<br />
O funcio<strong>na</strong>mento <strong>de</strong>stes dispositivos móveis está condicio<strong>na</strong><strong>do</strong> ao tempo <strong>de</strong> vida <strong>de</strong><br />
suas baterias. Este tempo, é por <strong>de</strong>finição o tempo que a bateria <strong>de</strong>mora para atingir um<br />
<strong>de</strong>termi<strong>na</strong><strong>do</strong> nível inferior <strong>de</strong> capacida<strong>de</strong> <strong>de</strong> carga, <strong>de</strong>nomi<strong>na</strong><strong>do</strong> nível <strong>de</strong> cutoff. Quan<strong>do</strong> o<br />
nível <strong>de</strong> cutoff é atingi<strong>do</strong> ocorre a impossibilida<strong>de</strong> da bateria <strong>de</strong> fornecer energia elétrica<br />
para o sistema, sen<strong>do</strong> consi<strong>de</strong>rada <strong>de</strong>scarregada. Neste contexto, surge a importância <strong>de</strong><br />
<strong>de</strong>finir um méto<strong>do</strong> que possa predizer o tempo <strong>de</strong> vida da bateria, e consequentemente,<br />
<strong>do</strong> sistema que é alimenta<strong>do</strong> por ela. Uma forma <strong>de</strong> realizar esta predição é a partir<br />
da aplicação <strong>de</strong> mo<strong>de</strong>los matemáticos <strong>de</strong> baterias que simulam a <strong>de</strong>scarga <strong>de</strong> energia <strong>do</strong><br />
sistema.<br />
Nos últimos anos, diferentes mo<strong>de</strong>los matemáticos <strong>de</strong> baterias foram <strong>de</strong>senvolvi<strong>do</strong>s,<br />
<strong>de</strong>ntre eles po<strong>de</strong>m ser cita<strong>do</strong>s: os a<strong>na</strong>líticos, os estocásticos, os elétricos e os eletroquími-<br />
cos, cada um com suas características e níveis <strong>de</strong> complexida<strong>de</strong>, que serão <strong>de</strong>scritos mais<br />
5
Capítulo 1. Apresentação da Dissertação 6<br />
<strong>de</strong>talhadamente no <strong>de</strong>correr <strong>de</strong>ste trabalho. Mas é importante <strong>de</strong>stacar, que os mo<strong>de</strong>los<br />
que apresentam os melhores resulta<strong>do</strong>s <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida da bateria, são<br />
os mo<strong>de</strong>los que consi<strong>de</strong>ram aspectos físicos das operações <strong>de</strong> <strong>de</strong>scarga da bateria. Estas<br />
operações <strong>de</strong> <strong>de</strong>scarga são compostas por perfis <strong>de</strong>fini<strong>do</strong>s como um conjunto <strong>de</strong> operações<br />
que um <strong>de</strong>termi<strong>na</strong><strong>do</strong> dispositivo po<strong>de</strong> realizar em um <strong>de</strong>termi<strong>na</strong><strong>do</strong> intervalo <strong>de</strong> tempo.<br />
Por exemplo, po<strong>de</strong>-se escutar música MP3, ou conversar a partir <strong>de</strong> um celular, porém,<br />
para cada uma <strong>de</strong>stas operações tem-se uma taxa <strong>de</strong> <strong>de</strong>scarga diferente <strong>na</strong> bateria [2].<br />
Estu<strong>do</strong>s em baterias [3,4] revelam que as taxas <strong>de</strong> <strong>de</strong>scarga são não-lineares no tempo<br />
e <strong>de</strong>pen<strong>de</strong>m da capacida<strong>de</strong> residual da bateria, ou seja, para diferentes perfis <strong>de</strong> <strong>de</strong>scarga<br />
são obti<strong>do</strong>s diferentes tempos <strong>de</strong> vida. Isto significa que a capacida<strong>de</strong> efetiva da bateria<br />
não é a mesma para diferentes perfis <strong>de</strong> <strong>de</strong>scarga. Além disto, em perío<strong>do</strong>s em que a<br />
corrente <strong>de</strong> <strong>de</strong>scarga é nula ou reduzida significativamente, ocorre um efeito <strong>de</strong> recuperação<br />
<strong>na</strong> capacida<strong>de</strong> da bateria. Então, para realizar a predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> uma<br />
bateria, é muito importante que o mo<strong>de</strong>lo a ser utiliza<strong>do</strong> consi<strong>de</strong>re os efeitos não-lineares<br />
<strong>de</strong>nomi<strong>na</strong><strong>do</strong>s efeito <strong>de</strong> recuperação e efeito da taxa <strong>de</strong> capacida<strong>de</strong>, ou parte <strong>de</strong>les, que<br />
ocorrem em um processo real <strong>de</strong> <strong>de</strong>scarga, e seus efeitos <strong>na</strong> capacida<strong>de</strong> da mesma.<br />
Portanto, neste trabalho é apresenta<strong>do</strong> o estu<strong>do</strong> e aplicação <strong>de</strong> um mo<strong>de</strong>lo elétrico<br />
<strong>de</strong> bateria que consi<strong>de</strong>ra parte <strong>do</strong>s efeitos não-lineares (i.e., efeito <strong>de</strong> recuperação), sen<strong>do</strong><br />
<strong>de</strong> fácil implementação, especialmente no que se refere a extração e <strong>de</strong>finição <strong>de</strong> seus<br />
parâmetros <strong>de</strong> configuração, <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> dispositivos móveis. Este<br />
mo<strong>de</strong>lo elétrico é <strong>de</strong>nomi<strong>na</strong><strong>do</strong> Battery e está presente <strong>na</strong> ferramenta computacio<strong>na</strong>l Mat-<br />
Lab/Simulink, sua aplicação é realizada em <strong>do</strong>is momentos: inicialmente é realizada a<br />
comparação <strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s a partir <strong>do</strong> mo<strong>de</strong>lo, com os resulta<strong>do</strong>s experimentais<br />
obti<strong>do</strong>s <strong>de</strong> uma plataforma <strong>de</strong> testes; em seguida é realizada a comparação <strong>do</strong>s resulta<strong>do</strong>s<br />
simula<strong>do</strong>s a partir <strong>do</strong> mo<strong>de</strong>lo, com resulta<strong>do</strong>s simula<strong>do</strong>s consi<strong>de</strong>ran<strong>do</strong> o mo<strong>de</strong>lo elétrico<br />
para Predizer Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong>, que é um mo<strong>de</strong>lo elétrico<br />
<strong>de</strong> alta acurácia encontra<strong>do</strong> <strong>na</strong> literatura.<br />
O restante <strong>de</strong>ste capítulo está organiza<strong>do</strong> da seguinte forma. Na Seção 1.2 é apre-<br />
sentada a motivação para a realização <strong>de</strong>ste trabalho. Na Seção 1.3 são apresenta<strong>do</strong>s os<br />
objetivos, geral e específicos. As contribuições são apresentadas <strong>na</strong> Seção 1.4. E fi<strong>na</strong>l-<br />
mente, <strong>na</strong> Seção 1.5, é apresentada a estrutura <strong>de</strong>ste trabalho.<br />
1.2 Motivação<br />
A cada dia é maior o número <strong>de</strong> dispositivos móveis utiliza<strong>do</strong>s <strong>na</strong>s mais diversas áreas,<br />
tais como: telefones celulares, re<strong>de</strong>s wireless, bluetooth, smartphones, câmeras digitais,<br />
notebooks, entre outros. Estes dispositivos tor<strong>na</strong>ram-se parte importante <strong>do</strong> cotidiano das
Capítulo 1. Apresentação da Dissertação 7<br />
pessoas e muitas vezes passam a ser indispensáveis <strong>na</strong> vida diária das mesmas, tanto no<br />
lazer, quanto no trabalho, pois são capazes <strong>de</strong> modificar suas roti<strong>na</strong>s e as suas formas <strong>de</strong><br />
tomar <strong>de</strong>cisões.<br />
Desta maneira, a mobilida<strong>de</strong> <strong>de</strong>ixa <strong>de</strong> ser uma facilida<strong>de</strong>, tor<strong>na</strong>n<strong>do</strong>-se uma necessida<strong>de</strong>,<br />
permitin<strong>do</strong> o acesso a da<strong>do</strong>s e informações em qualquer lugar a qualquer momento. Assim,<br />
tor<strong>na</strong>-se não somente uma opção para facilitar tarefas particulares, mas também uma<br />
oportunida<strong>de</strong> <strong>de</strong> melhora <strong>na</strong> gestão <strong>de</strong> negócios, possibilitan<strong>do</strong> a integração <strong>de</strong> dispositivos<br />
móveis com sistemas <strong>de</strong> gestão e e-business.<br />
Em virtu<strong>de</strong> das suas características, muitos <strong>de</strong>stes dispositivos não possuem qualquer<br />
tipo <strong>de</strong> conexão com uma re<strong>de</strong> elétrica que lhes permita um constante atendimento <strong>de</strong> suas<br />
necessida<strong>de</strong>s <strong>de</strong> consumo <strong>de</strong> energia, e acabam sen<strong>do</strong> supri<strong>do</strong>s por fontes <strong>de</strong> alimentação<br />
próprias e individuais, sen<strong>do</strong> a principal <strong>de</strong>las a utilização <strong>de</strong> baterias. Deste mo<strong>do</strong>, a<br />
utilização <strong>de</strong>stes dispositivos está diretamente limitada ao tempo <strong>de</strong> vida das baterias que<br />
os alimentam. Sen<strong>do</strong> assim, surge um novo <strong>de</strong>safio: <strong>de</strong>finir um méto<strong>do</strong> eficaz que permita<br />
predizer o tempo <strong>de</strong> vida das baterias, e consequentemente, <strong>do</strong> sistema como um to<strong>do</strong>.<br />
Conforme já mencio<strong>na</strong><strong>do</strong>, uma forma <strong>de</strong> realizar esta predição é a partir da utilização<br />
<strong>de</strong> mo<strong>de</strong>los matemáticos <strong>de</strong> baterias, cada um com suas características e grau <strong>de</strong> com-<br />
plexida<strong>de</strong>. Logicamente, <strong>de</strong>staca-se a importância <strong>de</strong> utilizar um mo<strong>de</strong>lo que possua bom<br />
nível <strong>de</strong> acurácia, que consi<strong>de</strong>re pelo menos algumas das características não-lineares da<br />
bateria e ao mesmo tempo seja prático e <strong>de</strong> fácil implementação. Sen<strong>do</strong> assim, através <strong>do</strong><br />
presente trabalho, será investiga<strong>do</strong>, a partir <strong>de</strong> estu<strong>do</strong>s comparativos, se o mo<strong>de</strong>lo elétrico<br />
Battery apresenta estas características.<br />
1.3 Objetivos<br />
Nesta seção são apresenta<strong>do</strong>s os objetivos <strong>de</strong>ste trabalho. Para facilitar a compreensão,<br />
optou-se em dividir os Objetivos em Objetivo Geral e Específicos, os quais são <strong>de</strong>talha<strong>do</strong>s<br />
a seguir.<br />
1.3.1 Objetivo Geral<br />
O objetivo geral <strong>de</strong>ste trabalho é estudar/aplicar a utilização <strong>de</strong> mo<strong>de</strong>los elétricos, que<br />
sejam capazes <strong>de</strong> simular a <strong>de</strong>scarga <strong>de</strong> energia <strong>de</strong> baterias que alimentam dispositivos<br />
móveis, possibilitan<strong>do</strong> a predição <strong>de</strong> seus tempos <strong>de</strong> vida.
Capítulo 1. Apresentação da Dissertação 8<br />
1.3.2 Objetivos Específicos<br />
Buscan<strong>do</strong> alcançar o objetivo geral <strong>de</strong>ste trabalho, os seguintes objetivos específicos foram<br />
traça<strong>do</strong>s:<br />
- Realizar uma revisão bibliográfica <strong>do</strong> esta<strong>do</strong> da arte <strong>do</strong>s diferentes tipos <strong>de</strong> mo<strong>de</strong>los<br />
matemáticos que <strong>de</strong>screvem a <strong>de</strong>scarga <strong>de</strong> energia em baterias, enfatizan<strong>do</strong> os mo<strong>de</strong>los<br />
elétricos;<br />
- Escolher <strong>de</strong>ntre os diferentes tipos <strong>de</strong> mo<strong>de</strong>los elétricos estuda<strong>do</strong>s, um mo<strong>de</strong>lo a<strong>de</strong>-<br />
qua<strong>do</strong> e prático, especialmente no que se refere a extração e <strong>de</strong>finição <strong>de</strong> seus parâmetros,<br />
para a predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> baterias;<br />
- Realizar o estu<strong>do</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico escolhi<strong>do</strong> a partir da comparação <strong>de</strong> seus resul-<br />
ta<strong>do</strong>s simula<strong>do</strong>s, com os resulta<strong>do</strong>s obti<strong>do</strong>s a partir <strong>de</strong> uma plataforma <strong>de</strong> testes;<br />
- Realizar o estu<strong>do</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico escolhi<strong>do</strong> a partir da comparação <strong>de</strong> seus re-<br />
sulta<strong>do</strong>s simula<strong>do</strong>s, com os resulta<strong>do</strong>s simula<strong>do</strong>s <strong>de</strong> um mo<strong>de</strong>lo elétrico <strong>de</strong> alta acurácia<br />
encontra<strong>do</strong> <strong>na</strong> literatura.<br />
1.4 Contribuições<br />
Nesta dissertação são introduzidas contribuições para a realização da predição <strong>do</strong> tempo<br />
<strong>de</strong> vida <strong>de</strong> dispositivos móveis que são alimenta<strong>do</strong>s por baterias. Estas contribuições são<br />
apresentada a seguir:<br />
1. Estu<strong>do</strong> e aplicação <strong>do</strong> mo<strong>de</strong>lo elétrico Battery para predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong><br />
baterias.<br />
2. <strong>Aplicação</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico Battery basea<strong>do</strong> em resulta<strong>do</strong>s reais obti<strong>do</strong>s a partir<br />
<strong>de</strong> uma plataforma <strong>de</strong> testes e não em resulta<strong>do</strong>s obti<strong>do</strong>s a partir <strong>do</strong> emprego <strong>de</strong><br />
simula<strong>do</strong>res <strong>de</strong> baterias (e.g., simula<strong>do</strong>r DualFoil) [2].<br />
3. Comparação <strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s a partir <strong>do</strong> mo<strong>de</strong>lo elétrico Battery com os<br />
resulta<strong>do</strong>s coleta<strong>do</strong>s a partir <strong>de</strong> uma plataforma <strong>de</strong> testes.<br />
4. Comparação <strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s a partir <strong>do</strong> mo<strong>de</strong>lo elétrico Battery com os<br />
resulta<strong>do</strong>s simula<strong>do</strong>s a partir <strong>de</strong> um mo<strong>de</strong>lo elétrico <strong>de</strong> alta acurácia encontra<strong>do</strong> <strong>na</strong><br />
literatura <strong>de</strong>nomi<strong>na</strong><strong>do</strong> mo<strong>de</strong>lo elétrico para Predizer Runtime e Características V-I<br />
<strong>de</strong> uma bateria.<br />
Durante o <strong>de</strong>senvolvimento <strong>de</strong>ste trabalho, seus resulta<strong>do</strong>s parciais já foram publica<strong>do</strong>s<br />
em eventos regio<strong>na</strong>is, <strong>na</strong>cio<strong>na</strong>is e periódico <strong>na</strong>cio<strong>na</strong>l. A relação completa <strong>do</strong>s trabalhos<br />
publica<strong>do</strong>s e <strong>do</strong>s que estão em processo <strong>de</strong> avaliação, po<strong>de</strong> ser verificada no Apêndice A.
Capítulo 1. Apresentação da Dissertação 9<br />
1.5 Estrutura <strong>do</strong> Documento<br />
Este trabalho está organiza<strong>do</strong> da seguinte forma:<br />
- No Capítulo 2 são abordadas as principais proprieda<strong>de</strong>s e características <strong>de</strong> uma<br />
bateria, os tipos <strong>de</strong> baterias mais utiliza<strong>do</strong>s atualmente, assim como alguns mo<strong>de</strong>los <strong>de</strong><br />
<strong>de</strong>scarga <strong>de</strong> baterias encontra<strong>do</strong>s <strong>na</strong> literatura, visan<strong>do</strong> uma melhor compreensão <strong>do</strong> con-<br />
texto <strong>de</strong> baterias no qual este trabalho está inseri<strong>do</strong>;<br />
- No Capítulo 3 é apresenta<strong>do</strong> o mo<strong>de</strong>lo elétrico Battery, escolhi<strong>do</strong> para ser utiliza<strong>do</strong><br />
neste trabalho. Este mo<strong>de</strong>lo utiliza equações matemáticas que <strong>de</strong>screvem o <strong>de</strong>caimento<br />
da tensão elétrica <strong>de</strong> baterias recarregáveis (e.g., íon lítio), sen<strong>do</strong> necessária a informação<br />
<strong>de</strong> parâmetros, como por exemplo: tipo <strong>de</strong> bateria, capacida<strong>de</strong> típica, resistência inter<strong>na</strong>,<br />
entre outros, on<strong>de</strong> os mesmos são extraí<strong>do</strong>s diretamente <strong>de</strong> uma curva característica <strong>de</strong><br />
<strong>de</strong>scarga da bateria simulada, em conjunto com da<strong>do</strong>s <strong>de</strong> seu datasheet;<br />
- No Capítulo 4 é apresentada a plataforma <strong>de</strong> testes, bem como a meto<strong>do</strong>logia uti-<br />
lizada <strong>na</strong> realização <strong>do</strong>s testes experimentais e também <strong>na</strong>s simulações com o mo<strong>de</strong>lo<br />
elétrico Battery. Também são apresenta<strong>do</strong>s os resulta<strong>do</strong>s das simulações <strong>de</strong>ste mo<strong>de</strong>lo e<br />
<strong>do</strong>s testes experimentais, juntamente com sua análise;<br />
- No Capítulo 5 é apresentada uma avaliação comparativa entre os resulta<strong>do</strong>s simula<strong>do</strong>s<br />
<strong>do</strong> mo<strong>de</strong>lo elétrico Battery e os resulta<strong>do</strong>s simula<strong>do</strong>s <strong>do</strong> mo<strong>de</strong>lo elétrico para Predizer<br />
Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong>;<br />
- E por fim, no Capítulo 6 são apresentadas as conclusões e os trabalhos futuros<br />
sugeri<strong>do</strong>s para este trabalho.
Capítulo 2<br />
Revisão Bibliográfica<br />
Neste capítulo são aborda<strong>do</strong>s alguns conceitos básicos relacio<strong>na</strong><strong>do</strong>s à bateria, suas princi-<br />
pais características e proprieda<strong>de</strong>s, com o objetivo <strong>de</strong> facilitar a leitura e a compreensão<br />
<strong>de</strong>ste trabalho. Também é realizada uma revisão bibliográfica <strong>do</strong> esta<strong>do</strong> da arte <strong>do</strong>s di-<br />
ferentes tipos <strong>de</strong> baterias utilizadas atualmente, e <strong>do</strong>s mo<strong>de</strong>los matemáticos encontra<strong>do</strong>s<br />
<strong>na</strong> literatura, usa<strong>do</strong>s para predição <strong>do</strong> seu tempo <strong>de</strong> vida.<br />
O restante <strong>de</strong>ste capítulo está organiza<strong>do</strong> como segue. Na Seção 2.1 são apresentadas<br />
as principais características das baterias. Na Seção 2.2 são <strong>de</strong>scritos os principais tipos<br />
<strong>de</strong> baterias. Na Seção 2.3 são <strong>de</strong>scritos alguns mo<strong>de</strong>los <strong>de</strong> baterias da literatura.<br />
2.1 <strong>Bateria</strong>s<br />
Uma bateria é constituída por uma ou mais células eletroquímicas ligadas em série ou em<br />
paralelo, ou ainda através <strong>de</strong> uma combi<strong>na</strong>ção mista <strong>de</strong> ambas. Nestas células a energia<br />
química armaze<strong>na</strong>da é convertida em energia elétrica, a partir <strong>de</strong> reações eletroquímicas,<br />
po<strong>de</strong>n<strong>do</strong> fornecer corrente elétrica para um circuito externo através <strong>de</strong> seus eletro<strong>do</strong>s 1<br />
<strong>de</strong>nomi<strong>na</strong><strong>do</strong>s <strong>de</strong> âno<strong>do</strong> (polarida<strong>de</strong> negativa) e <strong>de</strong> cáto<strong>do</strong> (polarida<strong>de</strong> positiva), separa<strong>do</strong>s<br />
por um eletrólito 2 .<br />
Durante o fornecimento <strong>de</strong> corrente elétrica da bateria para um circuito ou sistema,<br />
o âno<strong>do</strong> libera elétrons para o circuito e o cáto<strong>do</strong> recebe elétrons <strong>do</strong> circuito, sen<strong>do</strong> que<br />
estes elétrons são origi<strong>na</strong><strong>do</strong>s por reações eletroquímicas <strong>de</strong>ntro da bateria, e são chama<strong>do</strong>s<br />
<strong>de</strong> espécies eletroativas [4, 5]. O esquema <strong>de</strong> uma célula eletroquímica é apresenta<strong>do</strong> <strong>na</strong><br />
Figura 2.1.<br />
Em uma bateria eletroquímica, o produto <strong>de</strong> duas gran<strong>de</strong>zas importantes fornece a<br />
quantida<strong>de</strong> <strong>de</strong> energia armaze<strong>na</strong>da <strong>na</strong> mesma [3]. Estas gran<strong>de</strong>zas são: Tensão Elétrica<br />
1 condutores metálicos por on<strong>de</strong> a corrente elétrica se movimenta<br />
2 condutor <strong>de</strong> eletricida<strong>de</strong> sóli<strong>do</strong> ou líqui<strong>do</strong><br />
10
Capítulo 2. Revisão Bibliográfica 11<br />
Figura 2.1: Esquema <strong>de</strong> uma célula eletroquímica [5].<br />
ou Voltagem, medida em Volts (V) e Capacida<strong>de</strong>, medida em Ampère-hora (Ah). Teori-<br />
camente, uma bateria <strong>de</strong> 50 Ah po<strong>de</strong> fornecer a um sistema 5 A durante 10 horas ou<br />
50 A durante 1 hora. Na prática observa-se que as operações reais <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> uma<br />
bateria são influenciadas por efeitos não-lineares ocorri<strong>do</strong>s <strong>na</strong> mesma, que influenciam<br />
diretamente no seu tempo <strong>de</strong> vida [5].<br />
rias.<br />
A seguir serão <strong>de</strong>scritas algumas proprieda<strong>de</strong>s e características encontradas em bate-<br />
2.1.1 Nível <strong>de</strong> Cutoff<br />
O nível <strong>de</strong> cutoff é <strong>de</strong>fini<strong>do</strong> como o valor limite inferior <strong>de</strong> carga em que a bateria consegue<br />
fornecer uma tensão necessária para um sistema. Quan<strong>do</strong> este valor é atingi<strong>do</strong>, a bateria<br />
não é mais capaz <strong>de</strong> fornecer energia ao sistema <strong>de</strong>vi<strong>do</strong> à impossibilida<strong>de</strong> <strong>de</strong> ocorrerem<br />
reações eletroquímicas [4], isto não significa que ela está completamente <strong>de</strong>scarregada,<br />
mas sim, indisponível para alimentar o sistema.<br />
2.1.2 Tempo <strong>de</strong> Vida<br />
O tempo <strong>de</strong> vida da bateria é o tempo que ela <strong>de</strong>mora para atingir um <strong>de</strong>termi<strong>na</strong><strong>do</strong> nível<br />
inferior <strong>de</strong> carga (i.e., nível <strong>de</strong> cutoff ), on<strong>de</strong> a bateria não é mais capaz <strong>de</strong> fornecer energia<br />
elétrica para o sistema.<br />
2.1.3 Características Não-lineares<br />
Para a mo<strong>de</strong>lagem matemática <strong>do</strong> comportamento <strong>de</strong> baterias e posteriormente a predição<br />
<strong>do</strong> seu tempo <strong>de</strong> vida, é importante o conhecimento <strong>do</strong>s efeitos não-lineares que ocorrem
Capítulo 2. Revisão Bibliográfica 12<br />
durante um perío<strong>do</strong> <strong>de</strong> <strong>de</strong>scarga. O comportamento i<strong>de</strong>al seria que a capacida<strong>de</strong> da bateria<br />
fosse constante para qualquer corrente <strong>de</strong> <strong>de</strong>scarga, e que toda a sua energia armaze<strong>na</strong>da<br />
pu<strong>de</strong>sse ser utilizada. No entanto, em um processo real <strong>de</strong> <strong>de</strong>scarga, <strong>de</strong>vi<strong>do</strong> aos efeitos<br />
não-lineares, a tensão da bateria é reduzida durante a <strong>de</strong>scarga e a sua capacida<strong>de</strong> efetiva é<br />
reduzida para altas correntes. Destaca-se que, <strong>de</strong>pen<strong>de</strong>n<strong>do</strong> <strong>do</strong> tipo <strong>de</strong> bateria, estes efeitos<br />
têm maior ou menor influência em sua capacida<strong>de</strong> [3, 5]. A seguir são <strong>de</strong>scritas algumas<br />
características não-lineares, que se fazem presentes em um processo real <strong>de</strong> <strong>de</strong>scarga.<br />
Efeito <strong>de</strong> Recuperação<br />
Define-se por efeito <strong>de</strong> recuperação <strong>de</strong> uma bateria a reorganização <strong>do</strong>s elétrons no<br />
eletrólito em um intervalo <strong>de</strong> tempo em que a corrente <strong>de</strong> <strong>de</strong>scarga é nula ou reduzida<br />
significativamente. Assim, a capacida<strong>de</strong> efetiva da bateria é aumentada, pois uma maior<br />
quantida<strong>de</strong> <strong>de</strong> carga tor<strong>na</strong>-se disponível antes da bateria alcançar o nível <strong>de</strong> cutoff. Na<br />
Figura 2.2 são apresenta<strong>do</strong>s <strong>do</strong>is tipos <strong>de</strong> <strong>de</strong>scarga, para uma <strong>de</strong>scarga intermitente ou<br />
variável, o efeito <strong>de</strong> recuperação é observa<strong>do</strong>, resultan<strong>do</strong> em um aumento da capacida<strong>de</strong><br />
da bateria e consequentemente <strong>do</strong> seu tempo <strong>de</strong> vida, diferente <strong>do</strong> que ocorre no caso <strong>de</strong><br />
uma <strong>de</strong>scarga contínua, on<strong>de</strong> a bateria não possui tempo para recuperar-se [5].<br />
Figura 2.2: Tensão elétrica em função da <strong>de</strong>scarga [5].<br />
Na Figura 2.3 são apresentadas, <strong>de</strong> forma simplificada, as operações <strong>de</strong> uma bateria.<br />
Na Figura 2.3 (A), percebe-se que a bateria está totalmente carregada e a concentração<br />
<strong>de</strong> espécies eletroativas é constante por to<strong>do</strong> comprimento w <strong>do</strong> eletrólito. Quan<strong>do</strong> uma<br />
<strong>de</strong>scarga é iniciada (Figura 2.3 (B)), ocorrem reações eletroquímicas que resultam <strong>na</strong><br />
redução <strong>de</strong> espécies eletroativas próximas ao eletro<strong>do</strong>, geran<strong>do</strong> um gradiente <strong>de</strong> concen-<br />
tração no eletrólito, ocasio<strong>na</strong>n<strong>do</strong> a difusão <strong>de</strong> espécies eletroativas. No momento em que<br />
a corrente <strong>de</strong> <strong>de</strong>scarga é reduzida significativamente ou é nula, a bateria encontra-se em<br />
um perío<strong>do</strong> <strong>de</strong> relaxação (Figura 2.3 (C)), possibilitan<strong>do</strong> a reorganização <strong>do</strong>s elétrons,<br />
reequilibran<strong>do</strong> o sistema, e geran<strong>do</strong> um gradiente <strong>de</strong> concentração nulo no eletrólito, tor-
Capítulo 2. Revisão Bibliográfica 13<br />
<strong>na</strong>n<strong>do</strong> disponível uma maior quantida<strong>de</strong> <strong>de</strong> carga <strong>na</strong> superfície <strong>do</strong> eletro<strong>do</strong>, aumentan<strong>do</strong><br />
assim a capacida<strong>de</strong> efetiva da bateria (Figura 2.3 (D)). Destaca-se que a concentração <strong>de</strong><br />
espécies eletroativas <strong>na</strong> superfície <strong>do</strong> eletro<strong>do</strong> será menor <strong>do</strong> que a concentração inicial.<br />
À medida que a concentração <strong>na</strong> superfície <strong>do</strong> eletro<strong>do</strong> diminui, a tensão da bateria é<br />
reduzida. Quan<strong>do</strong> a tensão da bateria atinge um limite inferior ao valor <strong>de</strong> corte (i.e.,<br />
cutoff ), as reações eletroquímicas não po<strong>de</strong>m mais ocorrer e a bateria não é mais capaz <strong>de</strong><br />
fornecer energia ao sistema. Neste momento a bateria po<strong>de</strong> ser consi<strong>de</strong>rada <strong>de</strong>scarregada,<br />
apesar <strong>de</strong> ainda existirem espécies eletroativas no eletrólito (Figura 2.3 (E)) [2, 4].<br />
Figura 2.3: Esquema <strong>do</strong>s diferentes esta<strong>do</strong>s <strong>de</strong> operação da bateria [4].<br />
Efeito da Taxa <strong>de</strong> Capacida<strong>de</strong><br />
O efeito da taxa <strong>de</strong> capacida<strong>de</strong> [4,5] é <strong>de</strong>pen<strong>de</strong>nte da capacida<strong>de</strong> atual da bateria e da<br />
intensida<strong>de</strong> da corrente <strong>de</strong> <strong>de</strong>scarga, ou seja, para altas correntes <strong>de</strong> <strong>de</strong>scarga a capacida<strong>de</strong><br />
efetiva da bateria é baixa, já que não há tempo suficiente para a reorganização <strong>do</strong>s elétrons<br />
no eletrólito, (i.e., efeito <strong>de</strong> recuperação) reduzin<strong>do</strong> <strong>de</strong>sta forma, a capacida<strong>de</strong> da bateria.<br />
Já para correntes alter<strong>na</strong>das, a capacida<strong>de</strong> efetiva da bateria é aumentada, pois <strong>na</strong> troca<br />
<strong>de</strong> uma corrente alta para uma corrente baixa, ou em perío<strong>do</strong>s sem corrente, ocorre o<br />
efeito <strong>de</strong> recuperação, aumentan<strong>do</strong> a capacida<strong>de</strong> efetiva da bateria. Esta capacida<strong>de</strong><br />
nunca será maior ou igual a capacida<strong>de</strong> inicial [4, 5]. Então, observa-se que para altas<br />
taxas <strong>de</strong> <strong>de</strong>scarga, a capacida<strong>de</strong> efetiva da bateria diminui, assim como para baixas taxas,<br />
a sua capacida<strong>de</strong> efetiva aumenta [5].<br />
Na Figura 2.4 é apresentada uma relação entre a capacida<strong>de</strong> efetiva da bateria e a taxa<br />
<strong>de</strong> <strong>de</strong>scarga, observa-se que para altas taxas <strong>de</strong> <strong>de</strong>scarga a capacida<strong>de</strong> efetiva da bateria
Capítulo 2. Revisão Bibliográfica 14<br />
é reduzida.<br />
Figura 2.4: Capacida<strong>de</strong> da bateria em função <strong>de</strong> sua taxa <strong>de</strong> <strong>de</strong>scarga [5].<br />
2.1.4 Esta<strong>do</strong> <strong>de</strong> Carga da <strong>Bateria</strong><br />
O emprego <strong>de</strong> mo<strong>de</strong>los matemáticos em baterias permite a <strong>de</strong>termi<strong>na</strong>ção <strong>de</strong> informações<br />
importantes da mesma, tais como: tempo <strong>de</strong> vida, <strong>de</strong>sempenho, entre outras. Um<br />
parâmetro importante é o chama<strong>do</strong> esta<strong>do</strong> <strong>de</strong> carga da bateria (i.e., SOC ), que <strong>de</strong>fine a<br />
capacida<strong>de</strong> e a energia útil que po<strong>de</strong>m ser utilizadas em um da<strong>do</strong> momento. Os principais<br />
fatores que influenciam o SOC são: (1) resistência inter<strong>na</strong> (i.e., oposição à passagem <strong>de</strong><br />
corrente elétrica), (2) tipo <strong>de</strong> <strong>de</strong>scarga, (3) mo<strong>do</strong> <strong>de</strong> <strong>de</strong>scarga e (4) taxa <strong>de</strong> carga/<strong>de</strong>scarga<br />
[6], que serão <strong>de</strong>scritos a seguir.<br />
Resistência Inter<strong>na</strong><br />
Resistência <strong>de</strong> sobre-carga e sobre-<strong>de</strong>scarga: quan<strong>do</strong> a bateria é sobre-carregada ou<br />
sobre-<strong>de</strong>scarregada, sua resistência inter<strong>na</strong> aumenta significativamente <strong>de</strong>vi<strong>do</strong> a difusão<br />
no eletrólito.<br />
Resistência <strong>de</strong> auto-<strong>de</strong>scarga: consi<strong>de</strong>ra uma eletrólise (i.e., uma reação eletroquímica)<br />
em alta tensão e uma perda lenta <strong>de</strong> cargas elétricas através <strong>do</strong>s termi<strong>na</strong>is da bateria em<br />
baixa tensão. Este fenômeno é sensível à temperatura e varia inversamente com ela.<br />
Resistência para carga e <strong>de</strong>scarga: esta resistência <strong>de</strong>fine a resistência <strong>do</strong> eletrólito, <strong>do</strong><br />
fluí<strong>do</strong> e das placas da bateria.<br />
Tipo <strong>de</strong> Descarga<br />
Descarga contínua: ocorre quan<strong>do</strong> a bateria fornece energia continuamente para o<br />
sistema, sem que haja por parte da bateria, perío<strong>do</strong>s <strong>de</strong> relaxação, fazen<strong>do</strong> com que a<br />
capacida<strong>de</strong> da mesma diminua continuamente.<br />
Descarga intermitente: ocorre quan<strong>do</strong> a bateria fornece energia para o sistema por<br />
algum intervalo <strong>de</strong> tempo. Em outro intervalo o sistema é <strong>de</strong>sconecta<strong>do</strong> <strong>de</strong> mo<strong>do</strong> que a<br />
capacida<strong>de</strong> da bateria é recuperada (i.e., efeito <strong>de</strong> recuperação), e após algum tempo o
Capítulo 2. Revisão Bibliográfica 15<br />
sistema é conecta<strong>do</strong> e a <strong>de</strong>scarga recomeça. Neste tipo <strong>de</strong> <strong>de</strong>scarga, a capacida<strong>de</strong> efetiva<br />
da bateria é aumentada.<br />
Mo<strong>do</strong> <strong>de</strong> Descarga<br />
Sistema constante: ocorre quan<strong>do</strong> a bateria alimenta um sistema que possui a resistên-<br />
cia constante. Assim, durante o processo <strong>de</strong> <strong>de</strong>scarga da bateria, sua corrente diminui<br />
conforme diminui a sua tensão.<br />
Corrente constante: ocorre quan<strong>do</strong> a corrente <strong>de</strong> <strong>de</strong>scarga se mantém constante, ou<br />
seja, o sistema diminui continuamente sua resistência para manter constante a corrente<br />
<strong>de</strong> <strong>de</strong>scarga. A tensão da bateria <strong>de</strong>cai rapidamente <strong>de</strong>vi<strong>do</strong> a média alta da corrente <strong>de</strong><br />
<strong>de</strong>scarga.<br />
Potência constante: ocorre quan<strong>do</strong> uma potência elétrica constante é estabelecida<br />
pelo sistema alimenta<strong>do</strong> pela bateria, <strong>de</strong> mo<strong>do</strong> que a corrente <strong>de</strong> <strong>de</strong>scarga aumenta para<br />
compensar a diminuição da tensão. Neste mo<strong>do</strong> <strong>de</strong> <strong>de</strong>scarga, as baterias apresentam um<br />
menor tempo <strong>de</strong> vida.<br />
Taxa <strong>de</strong> Carga/Descarga<br />
A taxa <strong>de</strong> carga/<strong>de</strong>scarga indica a corrente <strong>de</strong> carga ou <strong>de</strong>scarga pela qual a bateria<br />
foi submetida. Para prolongar o tempo <strong>de</strong> vida útil <strong>de</strong> uma bateria, suas taxas <strong>de</strong> carga<br />
e <strong>de</strong>scarga não po<strong>de</strong>m ser muito altas, pois sen<strong>do</strong> altas, o seu tempo <strong>de</strong> vida po<strong>de</strong> ser<br />
reduzi<strong>do</strong> significativamente.<br />
2.1.5 Capacida<strong>de</strong><br />
A capacida<strong>de</strong> <strong>de</strong> uma bateria po<strong>de</strong> ser expressa <strong>de</strong> três formas diferentes [4, 7]:<br />
Capacida<strong>de</strong> teórica: baseia-se <strong>na</strong> quantida<strong>de</strong> <strong>de</strong> energia armaze<strong>na</strong>da, sen<strong>do</strong> o limite<br />
máximo <strong>de</strong> energia que po<strong>de</strong> ser extraí<strong>do</strong> da bateria <strong>na</strong> prática;<br />
Capacida<strong>de</strong> padrão: é a energia que po<strong>de</strong> ser extraída sob condições especificadas pelo<br />
fabricante;<br />
Capacida<strong>de</strong> atual: é aquela que po<strong>de</strong> exce<strong>de</strong>r a capacida<strong>de</strong> padrão, mas não po<strong>de</strong><br />
exce<strong>de</strong>r a capacida<strong>de</strong> teórica <strong>de</strong> uma bateria.<br />
Consi<strong>de</strong>ran<strong>do</strong> as características e proprieda<strong>de</strong>s <strong>de</strong> uma bateria, a sua capacida<strong>de</strong> em<br />
relação ao perfil <strong>de</strong> corrente <strong>de</strong> <strong>de</strong>scarga <strong>de</strong>pen<strong>de</strong> <strong>de</strong> <strong>do</strong>is efeitos, o efeito da taxa <strong>de</strong><br />
capacida<strong>de</strong> e o efeito <strong>de</strong> recuperação [4].<br />
2.2 Tipos <strong>de</strong> <strong>Bateria</strong>s<br />
Nesta seção são apresentadas as tecnologias utilizadas <strong>na</strong>s últimas décadas para aten<strong>de</strong>r<br />
a crescente <strong>de</strong>manda <strong>de</strong> baterias recarregáveis para dispositivos eletrônicos portáteis. Ao
Capítulo 2. Revisão Bibliográfica 16<br />
comparar tecnologias <strong>de</strong> baterias, os seguintes aspectos <strong>de</strong>vem ser consi<strong>de</strong>ra<strong>do</strong>s: <strong>de</strong>nsi-<br />
da<strong>de</strong> <strong>de</strong> energia (carga armaze<strong>na</strong>da por unida<strong>de</strong> <strong>de</strong> peso da bateria), ciclo <strong>de</strong> vida (o<br />
número <strong>de</strong> ciclos <strong>de</strong> carga/<strong>de</strong>scarga antes <strong>do</strong> <strong>de</strong>scarte da bateria), impacto ambiental,<br />
segurança, custo, tensão <strong>de</strong> alimentação disponível e características <strong>de</strong> carga/<strong>de</strong>scarga<br />
[7]. Na Figura 2.5 é apresentada a evolução das tecnologias <strong>de</strong> baterias recarregáveis em<br />
relação a <strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia.<br />
Figura 2.5: Densida<strong>de</strong> <strong>de</strong> energia e ano <strong>de</strong> implantação comercial das tecnologias <strong>de</strong> baterias<br />
[7].<br />
Dentre as tecnologias mais populares <strong>de</strong> baterias recarregáveis para dispositivos eletrôni-<br />
cos portáteis existem:<br />
Níquel Cádmio (Ni-Cd)<br />
Esta é uma tecnologia consolidada, utilizada por várias décadas para <strong>de</strong>senvolver ba-<br />
terias recarregáveis para aparelhos eletrônicos portáteis. Suas vantagens incluem o baixo<br />
custo e a possibilida<strong>de</strong> <strong>de</strong> serem utilizadas em operações <strong>de</strong> altas taxas <strong>de</strong> <strong>de</strong>scarga. Em-<br />
bora a tecnologia Ni-Cd venha per<strong>de</strong>n<strong>do</strong> espaço nos últimos anos <strong>de</strong>vi<strong>do</strong> a sua baixa<br />
<strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia e toxicida<strong>de</strong>, esta tecnologia é ainda bastante empregada em apli-<br />
cações <strong>de</strong> baixo custo, como por exemplo, rádios portáteis [7, 8].<br />
Níquel Metal Hidreto (Ni-MH)<br />
Estas baterias foram bastante utilizadas nos últimos anos para alimentar computa<strong>do</strong>res<br />
portáteis. Elas têm aproximadamente duas vezes a <strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia das baterias Ni-<br />
Cd. No entanto, elas possuem um ciclo <strong>de</strong> vida mais curto, são mais caras e ineficientes<br />
em altas taxas <strong>de</strong> <strong>de</strong>scarga [7, 8].
Capítulo 2. Revisão Bibliográfica 17<br />
Íon Lítio (Li-íon)<br />
Esta é a tecnologia <strong>de</strong> bateria que mais cresce atualmente, possui <strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia<br />
significativamente superior e ciclo <strong>de</strong> vida aproximadamente duas vezes maior <strong>do</strong> que os<br />
ciclos <strong>de</strong> vida das baterias Ni-MH. As baterias <strong>de</strong> Li-íon são mais sensíveis às caracte-<br />
rísticas da corrente <strong>de</strong> <strong>de</strong>scarga e mais caras que as baterias Ni-MH. Por outro la<strong>do</strong>,<br />
elas apresentam um longo tempo <strong>de</strong> vida, por isso são mais populares sen<strong>do</strong> usadas em<br />
notebooks, Perso<strong>na</strong>l Digital Assistants (PDAs) e celulares [7, 8].<br />
Alcali<strong>na</strong> Recarregável<br />
Enquanto as baterias alcali<strong>na</strong>s <strong>de</strong>scartáveis têm si<strong>do</strong> utilizadas por muitos anos, a<br />
tecnologia reutilizável alcali<strong>na</strong> manganês foi <strong>de</strong>senvolvida como uma alter<strong>na</strong>tiva <strong>de</strong> baixo<br />
custo em que a <strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia e o ciclo <strong>de</strong> vida são comprometi<strong>do</strong>s [7]. Embora a<br />
<strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia inicial das baterias alcali<strong>na</strong>s reutilizáveis seja superior a <strong>de</strong>nsida<strong>de</strong><br />
das baterias <strong>de</strong> Ni-Cd, verificou-se uma rápida diminuição <strong>de</strong>sta <strong>de</strong>nsida<strong>de</strong> com relação<br />
aos ciclos <strong>de</strong> vida. Por exemplo, após 10 ciclos, é normalmente observada uma redução<br />
<strong>de</strong> 50% <strong>na</strong> <strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia <strong>de</strong>ste tipo <strong>de</strong> bateria, e em 50 ciclos, observa-se uma<br />
redução <strong>de</strong> 75% [7].<br />
Lítio Polímero<br />
Esta tecnologia emergente permite a fabricação <strong>de</strong> baterias ultra-fi<strong>na</strong>s (i.e., menos <strong>de</strong><br />
1 mm <strong>de</strong> espessura) e espera-se através <strong>de</strong>la, aten<strong>de</strong>r as necessida<strong>de</strong>s da próxima geração<br />
<strong>de</strong> computação e comunicação <strong>de</strong> dispositivos portáteis em relação ao tamanho e peso.<br />
Além disto, elas têm melhorias em relação à tecnologia <strong>de</strong> íon lítio no que se refere a<br />
<strong>de</strong>nsida<strong>de</strong> <strong>de</strong> energia e segurança. Como pontos negativos, estas baterias possuem um<br />
alto custo <strong>de</strong> fabricação e enfrentam <strong>de</strong>safios no gerenciamento térmico interno [7, 8].<br />
2.3 Mo<strong>de</strong>los <strong>de</strong> <strong>Bateria</strong>s<br />
Nesta seção são apresenta<strong>do</strong>s alguns mo<strong>de</strong>los <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> baterias encontra<strong>do</strong>s <strong>na</strong> lit-<br />
eratura. Conforme os trabalhos correlatos, existe uma varieda<strong>de</strong> <strong>de</strong> mo<strong>de</strong>los matemáticos<br />
<strong>de</strong> baterias com diferentes características. Estes mo<strong>de</strong>los capturam as características reais<br />
<strong>de</strong> operação da bateria e po<strong>de</strong>m ser utiliza<strong>do</strong>s para prever o comportamento da mesma<br />
sob várias condições <strong>de</strong> carga e <strong>de</strong>scarga. Também são úteis para projetos <strong>de</strong> sistemas<br />
alimenta<strong>do</strong>s, porque permitem a análise <strong>do</strong> comportamento <strong>de</strong> <strong>de</strong>scarga da bateria sob<br />
diferentes especificações <strong>do</strong> projeto [4]. A seguir são apresenta<strong>do</strong>s alguns <strong>de</strong>stes mo<strong>de</strong>los.<br />
2.3.1 Mo<strong>de</strong>los Eletroquímicos<br />
Os mo<strong>de</strong>los eletroquímicos baseiam-se nos processos químicos que ocorrem <strong>na</strong> bateria, e<br />
são <strong>de</strong>scritos com gran<strong>de</strong> <strong>de</strong>talhamento, necessitan<strong>do</strong> a informação <strong>de</strong> uma gran<strong>de</strong> quanti-
Capítulo 2. Revisão Bibliográfica 18<br />
da<strong>de</strong> <strong>de</strong> parâmetros da bateria. Por isto, estes mo<strong>de</strong>los são altamente complexos e difíceis<br />
<strong>de</strong> implementar, são consi<strong>de</strong>ra<strong>do</strong>s <strong>na</strong> literatura os mo<strong>de</strong>los <strong>de</strong> maior acurácia [5, 7].<br />
Doyle, Fuller e Newman [5, 9] <strong>de</strong>senvolveram um mo<strong>de</strong>lo eletroquímico para células<br />
<strong>de</strong> lítio e íon lítio. A resolução das equações que compõem o mo<strong>de</strong>lo fornecem a tensão<br />
e a corrente em função <strong>do</strong> tempo; as fases <strong>de</strong> potencial no eletrólito e no eletro<strong>do</strong>; a<br />
concentração sali<strong>na</strong>; a taxa <strong>de</strong> reação e a <strong>de</strong>nsida<strong>de</strong> <strong>de</strong> corrente no eletrólito em função<br />
<strong>do</strong> tempo e da posição <strong>na</strong> célula [5, 9, 10].<br />
Este mo<strong>de</strong>lo po<strong>de</strong> ser encontra<strong>do</strong> no programa Fortran Dualfoil, disponível gratuita-<br />
mente <strong>na</strong> Internet [11], para simular baterias <strong>de</strong> íon lítio. O programa computa a mudança<br />
<strong>de</strong> todas as proprieda<strong>de</strong>s da bateria ao longo <strong>do</strong> tempo para um perfil <strong>de</strong> carga <strong>de</strong>fini<strong>do</strong><br />
pelo usuário. A partir <strong>do</strong>s da<strong>do</strong>s <strong>de</strong> saída <strong>do</strong> programa, é possível obter o tempo <strong>de</strong> vida<br />
da bateria. Além <strong>do</strong> perfil <strong>de</strong> carga, o usuário tem que <strong>de</strong>finir mais <strong>de</strong> cinquenta parâme-<br />
tros relacio<strong>na</strong><strong>do</strong>s à bateria, como por exemplo: a espessura <strong>do</strong>s eletro<strong>do</strong>s, a concentração<br />
inicial <strong>de</strong> sal no eletrólito e a capacida<strong>de</strong> global <strong>de</strong> calor. Para a <strong>de</strong>finição <strong>de</strong> to<strong>do</strong>s estes<br />
parâmetros, <strong>de</strong>ve-se ter um conhecimento <strong>de</strong>talha<strong>do</strong> sobre a bateria que está sen<strong>do</strong> mo<strong>de</strong>-<br />
lada. Sen<strong>do</strong> assim, <strong>de</strong>vi<strong>do</strong> à alta precisão <strong>do</strong>s resulta<strong>do</strong>s gera<strong>do</strong>s pelo programa Dualfoil,<br />
ele acaba sen<strong>do</strong> frequentemente utiliza<strong>do</strong> como uma referência <strong>de</strong> comparação com outros<br />
mo<strong>de</strong>los da literatura, em substituição a utilização <strong>de</strong> resulta<strong>do</strong>s experimentais [5].<br />
2.3.2 Mo<strong>de</strong>los Estocásticos<br />
Os mo<strong>de</strong>los estocásticos <strong>de</strong>screvem a bateria <strong>de</strong> uma forma abstrata, on<strong>de</strong> o processo <strong>de</strong><br />
<strong>de</strong>scarga e o efeito <strong>de</strong> recuperação são <strong>de</strong>scritos como processos estocásticos [5, 12].<br />
Um <strong>do</strong>s primeiros mo<strong>de</strong>los estocásticos <strong>de</strong> bateria foi <strong>de</strong>senvolvi<strong>do</strong> por Chiasserini e<br />
Rao [5, 13, 14] utilizan<strong>do</strong> ca<strong>de</strong>ias <strong>de</strong> Markov para mo<strong>de</strong>lar a <strong>de</strong>scarga <strong>de</strong> uma bateria. O<br />
estu<strong>do</strong> <strong>de</strong> predição é elabora<strong>do</strong> com base em esta<strong>do</strong>s que correspon<strong>de</strong>m ao número <strong>de</strong><br />
unida<strong>de</strong>s <strong>de</strong> carga que a bateria é capaz <strong>de</strong> efetuar, on<strong>de</strong> estes esta<strong>do</strong>s po<strong>de</strong>m assumir<br />
valores, referentes à transmissão <strong>de</strong> um pacote <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> energia, ou pela transmissão<br />
contínua <strong>de</strong> corrente [13]. Neste mo<strong>de</strong>lo, em cada etapa <strong>de</strong> consumo <strong>de</strong> energia é efetua<strong>do</strong><br />
um cálculo probabilístico, cujo objetivo é verificar a recuperação da bateria em função da<br />
mudança <strong>de</strong> corrente. Este cálculo também é feito <strong>na</strong> transmissão <strong>do</strong> pacote <strong>de</strong> <strong>de</strong>scarga,<br />
em que é verificada a <strong>de</strong>scarga da bateria.<br />
Os resulta<strong>do</strong>s apresenta<strong>do</strong>s pelos mo<strong>de</strong>los estocásticos <strong>de</strong> Chiasserini e Rao [5], têm<br />
um <strong>de</strong>svio máximo <strong>de</strong> 4%, com um <strong>de</strong>svio médio <strong>de</strong> 1% em relação ao mo<strong>de</strong>lo eletro-<br />
químico. Estes resulta<strong>do</strong>s mostram que o mo<strong>de</strong>lo estocástico fornece uma boa <strong>de</strong>scrição<br />
<strong>do</strong> comportamento da bateria com relação a <strong>de</strong>scargas variáveis. Porém, estes mo<strong>de</strong>los<br />
possuem limitações, pois ape<strong>na</strong>s o efeito <strong>de</strong> recuperação é consi<strong>de</strong>ra<strong>do</strong>.
Capítulo 2. Revisão Bibliográfica 19<br />
2.3.3 Mo<strong>de</strong>los A<strong>na</strong>líticos<br />
Nos mo<strong>de</strong>los a<strong>na</strong>líticos, as principais proprieda<strong>de</strong>s da bateria são mo<strong>de</strong>ladas utilizan<strong>do</strong>-se<br />
um conjunto menor <strong>de</strong> equações, o que tor<strong>na</strong> a implementação <strong>de</strong>ste tipo <strong>de</strong> mo<strong>de</strong>lo mais<br />
simples quan<strong>do</strong> comparada a outros mo<strong>de</strong>los [5]. Os mo<strong>de</strong>los a<strong>na</strong>líticos po<strong>de</strong>m incluir<br />
mo<strong>de</strong>los <strong>de</strong> carga constante ou carga variável, e alguns conseguem capturar o efeito da<br />
taxa <strong>de</strong> capacida<strong>de</strong> e o efeito <strong>de</strong> recuperação. Eles são computacio<strong>na</strong>lmente eficientes e<br />
flexíveis, requeren<strong>do</strong> avaliação <strong>de</strong> simples expressões a<strong>na</strong>líticas e po<strong>de</strong>m ser facilmente<br />
estendi<strong>do</strong>s para diferentes tipos <strong>de</strong> baterias [7, 15].<br />
A seguir são <strong>de</strong>scritos alguns mo<strong>de</strong>los a<strong>na</strong>líticos encontra<strong>do</strong>s <strong>na</strong> literatura.<br />
Mo<strong>de</strong>lo Linear<br />
O mo<strong>de</strong>lo a<strong>na</strong>lítico mais simples é o Mo<strong>de</strong>lo Linear [3–5], on<strong>de</strong> a bateria é consi<strong>de</strong>rada<br />
como um recipiente linear <strong>de</strong> corrente. A capacida<strong>de</strong> restante C <strong>de</strong> uma bateria é dada<br />
por:<br />
C = C ′ − I·td, (2.1)<br />
on<strong>de</strong>: C ′ é a capacida<strong>de</strong> no início da operação, I é a corrente constante <strong>de</strong> <strong>de</strong>scarga<br />
durante a operação, e td é o tempo <strong>de</strong> duração da operação. A capacida<strong>de</strong> remanescente<br />
é calculada sempre que a corrente <strong>de</strong> <strong>de</strong>scarga for alterada. Porém, este mo<strong>de</strong>lo não é<br />
capaz <strong>de</strong> capturar os efeitos não-lineares da bateria, que ocorrem <strong>na</strong>s operações físicas <strong>de</strong><br />
<strong>de</strong>scarga, e que influenciam diretamente em sua capacida<strong>de</strong>, e consequentemente no seu<br />
tempo <strong>de</strong> vida.<br />
Lei <strong>de</strong> Peukert<br />
O mo<strong>de</strong>lo a<strong>na</strong>lítico mais simples, que consi<strong>de</strong>ra parte <strong>de</strong> suas proprieda<strong>de</strong>s não-lineares,<br />
para a estimação <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> baterias é a Lei <strong>de</strong> Peukert [5,15]. Esta lei captura<br />
a relação não-linear entre o tempo <strong>de</strong> vida da bateria e sua taxa <strong>de</strong> <strong>de</strong>scarga, mas não<br />
consi<strong>de</strong>ra o efeito <strong>de</strong> recuperação. Conforme a lei <strong>de</strong> Peukert, o tempo <strong>de</strong> vida L da<br />
bateria é da<strong>do</strong> por:<br />
L = a<br />
,<br />
Ib (2.2)<br />
on<strong>de</strong>: I é a corrente <strong>de</strong> <strong>de</strong>scarga, e a e b são parâmetros que <strong>de</strong>pen<strong>de</strong>m da bateria e são<br />
obti<strong>do</strong>s a partir <strong>de</strong> experimentos.<br />
Os resulta<strong>do</strong>s da aplicação da Lei <strong>de</strong> Peukert, para a predição <strong>do</strong> tempo <strong>de</strong> vida da<br />
bateria, são razoáveis para cargas constantes, mas não para cargas variáveis [2]. Conforme<br />
[16], Rakhmatov e Vrudhula propõem um mo<strong>de</strong>lo que é uma extensão da Lei <strong>de</strong> Peukert<br />
para cargas variáveis, on<strong>de</strong> a corrente I da equação (2.2) é substituída pela média da<br />
corrente até um tempo t = L. Para uma seção <strong>de</strong> perfil <strong>de</strong> <strong>de</strong>scarga constante, com
Capítulo 2. Revisão Bibliográfica 20<br />
pontos no tempo tk e mudanças <strong>na</strong> corrente Ik, conforme apresenta<strong>do</strong> <strong>na</strong> Figura 2.6,<br />
tem-se que o tempo <strong>de</strong> vida L da bateria é da<strong>do</strong> por:<br />
L =<br />
a<br />
( ∑ n<br />
k=1 Ik(tk−tk−1)<br />
L<br />
Figura 2.6: Seção <strong>de</strong> perfil <strong>de</strong> <strong>de</strong>scarga constante [5].<br />
) b . (2.3)<br />
Para n = 1, a equação (2.3) reduz-se a equação (2.2). Apesar da extensão da Lei <strong>de</strong><br />
Peukert tratar <strong>de</strong> perfis <strong>de</strong> <strong>de</strong>scargas variáveis, o mo<strong>de</strong>lo é consi<strong>de</strong>ra<strong>do</strong> simples porque<br />
ape<strong>na</strong>s a média das correntes <strong>de</strong> <strong>de</strong>scarga é utilizada e o efeito <strong>de</strong> recuperação não é<br />
mo<strong>de</strong>la<strong>do</strong>. Na prática, embora perfis <strong>de</strong> <strong>de</strong>scarga apresentem a mesma média, o efeito<br />
<strong>de</strong> recuperação <strong>na</strong> bateria po<strong>de</strong>rá ser diferente, <strong>de</strong>pen<strong>de</strong>n<strong>do</strong> da or<strong>de</strong>m <strong>do</strong>s <strong>de</strong>graus das<br />
correntes <strong>de</strong> <strong>de</strong>scarga, o que vai influenciar no tempo <strong>de</strong> vida da bateria [5].<br />
Mo<strong>de</strong>lo <strong>de</strong> Difusão <strong>de</strong> Rakhmatov e Vrudhula<br />
O mo<strong>de</strong>lo <strong>de</strong> difusão <strong>de</strong> Rakhmatov e Vrudhula [16] <strong>de</strong>screve a evolução da concentração<br />
das espécies eletroativas no eletrólito, sen<strong>do</strong> utiliza<strong>do</strong> para predizer o tempo <strong>de</strong> vida <strong>de</strong><br />
uma bateria a partir <strong>de</strong> uma corrente <strong>de</strong> <strong>de</strong>scarga constante ou variável [2]. Observa-se que<br />
os processos químicos ocorri<strong>do</strong>s em ambos os eletro<strong>do</strong>s são consi<strong>de</strong>ra<strong>do</strong>s idênticos, mas<br />
por simplificação <strong>de</strong> mo<strong>de</strong>lagem, ape<strong>na</strong>s o fluxo em um <strong>de</strong>les é consi<strong>de</strong>ra<strong>do</strong>. Este mo<strong>de</strong>lo<br />
captura os efeitos não-lineares (i.e., efeito da taxa <strong>de</strong> capacida<strong>de</strong> e efeito <strong>de</strong> recuperação)<br />
que ocorrem durante um processo <strong>de</strong> <strong>de</strong>scarga da bateria, e que afetam diretamente sua<br />
capacida<strong>de</strong> e seu tempo <strong>de</strong> vida. É também consi<strong>de</strong>ra<strong>do</strong> <strong>de</strong> fácil implementação quan<strong>do</strong><br />
compara<strong>do</strong> com os <strong>de</strong>mais mo<strong>de</strong>los [4], em especial aos mo<strong>de</strong>los eletroquímicos.<br />
Conforme a Figura 2.7, a difusão é consi<strong>de</strong>rada unidimensio<strong>na</strong>l em uma região <strong>de</strong><br />
comprimento w <strong>do</strong> eletrólito. Seja C(x, t) a concentração <strong>de</strong> espécies eletroativas no
Capítulo 2. Revisão Bibliográfica 21<br />
Figura 2.7: <strong>Bateria</strong> totalmente carregada [4].<br />
tempo t ϵ [0, L] e <strong>na</strong> distância x ϵ [0,w] <strong>do</strong> eletro<strong>do</strong>, quan<strong>do</strong> uma bateria está completa-<br />
mente carregada, a concentração é constante através <strong>do</strong> comprimento <strong>do</strong> eletrólito, isto é,<br />
C(x, 0) = C ∗ , on<strong>de</strong> C ∗ representa a capacida<strong>de</strong> inicial da bateria. A bateria é consi<strong>de</strong>rada<br />
<strong>de</strong>scarregada quan<strong>do</strong> C(0, t) é inferior ao nível <strong>de</strong> cutoff. A evolução da concentração é<br />
<strong>de</strong>scrita pelas Leis <strong>de</strong> Fick [15,17], dadas pelo sistema <strong>de</strong> Equações Diferenciais Parciais<br />
(EDPs), apresenta<strong>do</strong> a seguir:<br />
{ −J(x, t) = D ∂C(x,t)<br />
∂x<br />
∂C(x,t)<br />
∂t<br />
= D ∂2 C(x,t)<br />
∂ 2 x<br />
(2.4)<br />
on<strong>de</strong>: J(x, t) é o fluxo <strong>de</strong> espécies eletroativas em função <strong>do</strong> tempo t e em função <strong>de</strong> uma<br />
distância x <strong>do</strong> eletro<strong>do</strong>, D é a constante <strong>de</strong> difusão, e C(x, t) é a função concentração<br />
<strong>de</strong> espécies eletroativas no tempo t ∈ [0, L] e <strong>na</strong> distância x ∈ [0, w]. Para uma bateria<br />
completamente carregada (i.e., t = 0), a concentração <strong>de</strong> espécies eletroativas é constante<br />
no comprimento <strong>do</strong> eletrólito, proporcio<strong>na</strong>n<strong>do</strong> a seguinte condição inicial:<br />
C(x, 0) = C ∗ . (2.5)<br />
A bateria é consi<strong>de</strong>rada <strong>de</strong>scarregada quan<strong>do</strong> C(0, t) é inferior ao nível <strong>de</strong> cutoff. De<br />
acor<strong>do</strong> com a Lei <strong>de</strong> Faraday, o fluxo <strong>de</strong> espécies eletroativas <strong>na</strong> superfície <strong>do</strong> eletro<strong>do</strong><br />
(x = 0) é proporcio<strong>na</strong>l à corrente i(t) (i.e., carga exter<strong>na</strong> aplicada). Por outro la<strong>do</strong>, por<br />
simplificação <strong>de</strong> mo<strong>de</strong>lagem, o fluxo <strong>na</strong> extremida<strong>de</strong> x = w é consi<strong>de</strong>ra<strong>do</strong> zero. Estas<br />
suposições fornecem as seguintes condições <strong>de</strong> fronteira:<br />
i(t)<br />
vF A<br />
∂C(x, t)<br />
= D |x=0, (2.6)<br />
∂x<br />
∂C(x, t)<br />
0 = D |x=w, (2.7)<br />
∂x<br />
on<strong>de</strong>: A é a área da superfície <strong>do</strong> eletro<strong>do</strong>, F é a constante <strong>de</strong> Faraday, e v é o número<br />
<strong>de</strong> elétrons envolvi<strong>do</strong>s <strong>na</strong> reação eletroquímica <strong>na</strong> superfície <strong>do</strong> eletro<strong>do</strong>.<br />
A partir <strong>de</strong> manipulações matemáticas, <strong>de</strong>scritas em [17], obtém-se uma solução<br />
a<strong>na</strong>lítica para o sistema <strong>de</strong> EDPs, apresenta<strong>do</strong> <strong>na</strong> equação (2.4), que relacio<strong>na</strong> o tempo
Capítulo 2. Revisão Bibliográfica 22<br />
<strong>de</strong> vida L para um perfil <strong>de</strong> <strong>de</strong>scarga i(t) com os parâmetros α e β que necessitam ser<br />
estima<strong>do</strong>s, dada pela expressão geral:<br />
α =<br />
∫ L<br />
0<br />
i(τ)<br />
√ L − τ dτ + 2<br />
∞∑<br />
n=1<br />
∫ L<br />
0<br />
i(τ)<br />
√ L − τ e − β2 n 2<br />
(L−τ) dτ, (2.8)<br />
on<strong>de</strong>: α relacio<strong>na</strong>-se com a capacida<strong>de</strong> da bateria, β relacio<strong>na</strong>-se com o comportamento<br />
não-linear da bateria, L é o tempo <strong>de</strong> vida da bateria e i(τ) é o perfil <strong>de</strong> <strong>de</strong>scarga.<br />
Rakhmatov-Vrudhula [15] comparam o seu mo<strong>de</strong>lo com o programa <strong>de</strong> simulação Dual-<br />
foil, e com a versão estendida da Lei <strong>de</strong> Peukert, <strong>na</strong> qual é possível utilizar cargas variáveis.<br />
Os resulta<strong>do</strong>s <strong>de</strong> simulação <strong>do</strong> Dualfoil são usa<strong>do</strong>s como valores <strong>de</strong> referência, uma vez<br />
que possuem boa acurácia. Para cargas contínuas, o mo<strong>de</strong>lo <strong>de</strong> Rakhmatov-Vrudhula<br />
prediz o tempo <strong>de</strong> vida com um erro médio <strong>de</strong> 3%, e um erro máximo <strong>de</strong> 6% em compara-<br />
ção com os resulta<strong>do</strong>s obti<strong>do</strong>s utilizan<strong>do</strong>-se <strong>do</strong> programa Dualfoil [5]. Por outro la<strong>do</strong>, a<br />
Lei <strong>de</strong> Peukert apresentou um erro médio <strong>de</strong> 14% e um erro máximo <strong>de</strong> 43%. A Lei <strong>de</strong><br />
Peukert tem si<strong>do</strong> utilizada <strong>de</strong> forma satisfatória para cargas baixas, mas os erros aumen-<br />
tam significativamente para cargas altas. Para cargas variáveis e interrompidas, a análise<br />
<strong>do</strong> mo<strong>de</strong>lo a<strong>na</strong>lítico <strong>de</strong> Rakhmatov-Vrudhula apresenta melhores resulta<strong>do</strong>s, ou seja, um<br />
erro máximo <strong>de</strong> 2,7% e um erro médio abaixo <strong>de</strong> 1%. Neste cenário, a Lei <strong>de</strong> Peukert<br />
não apresenta bons resulta<strong>do</strong>s, principalmente por não consi<strong>de</strong>rar um efeito não-linear<br />
importante <strong>na</strong> bateria, que é o efeito <strong>de</strong> recuperação [2].<br />
2.3.4 Mo<strong>de</strong>los <strong>Elétricos</strong><br />
Os mo<strong>de</strong>los elétricos [1], também <strong>de</strong>nomi<strong>na</strong><strong>do</strong>s mo<strong>de</strong>los <strong>de</strong> circuitos elétricos, empregam<br />
uma combi<strong>na</strong>ção <strong>de</strong> fontes <strong>de</strong> tensão, resistores e capacitores para a simulação <strong>de</strong> <strong>de</strong>scarga<br />
da bateria. A acurácia <strong>de</strong>stes mo<strong>de</strong>los, em relação a predição <strong>do</strong> tempo <strong>de</strong> vida da bateria,<br />
situa-se entre a acurácia <strong>do</strong>s mo<strong>de</strong>los a<strong>na</strong>líticos e eletroquímicos (i.e., entre 1% e 5% [1]).<br />
Na área <strong>de</strong> engenharia elétrica, estes mo<strong>de</strong>los são mais intuitivos e fáceis <strong>de</strong> manusear,<br />
principalmente quan<strong>do</strong> utiliza<strong>do</strong>s a partir <strong>de</strong> simula<strong>do</strong>res <strong>de</strong> circuito.<br />
A forma básica <strong>do</strong>s mo<strong>de</strong>los elétricos para diferentes tipos <strong>de</strong> baterias é a mesma,<br />
um capacitor representa a capacida<strong>de</strong> da bateria, uma taxa <strong>de</strong> <strong>de</strong>scarga normaliza<strong>do</strong>ra<br />
<strong>de</strong>termi<strong>na</strong> a perda <strong>de</strong> capacida<strong>de</strong> em altas correntes <strong>de</strong> <strong>de</strong>scarga, um circuito representa<br />
o consumo da capacida<strong>de</strong> da bateria, uma tabela <strong>de</strong> pesquisa representa a tensão versus<br />
o esta<strong>do</strong> da carga, e um resistor representa a resistência inter<strong>na</strong> da bateria.<br />
Nesta seção são <strong>de</strong>scritas características <strong>de</strong> alguns mo<strong>de</strong>los elétricos encontra<strong>do</strong>s <strong>na</strong><br />
literatura, como por exemplo, mo<strong>de</strong>los basea<strong>do</strong>s em Thevenin, mo<strong>de</strong>los basea<strong>do</strong>s em<br />
Impedância, mo<strong>de</strong>los basea<strong>do</strong>s em Runtime, mo<strong>de</strong>lo para Predizer Runtime e Caracterís-<br />
ticas V-I <strong>de</strong> uma <strong>Bateria</strong>, e o mo<strong>de</strong>lo Battery.
Capítulo 2. Revisão Bibliográfica 23<br />
Conforme [1], existem muitos mo<strong>de</strong>los elétricos <strong>de</strong> baterias, <strong>de</strong>s<strong>de</strong> chumbo-áci<strong>do</strong> até<br />
íon lítio polímero. Porém, a maioria <strong>de</strong>stes mo<strong>de</strong>los elétricos po<strong>de</strong> ser classificada em três<br />
categorias básicas: mo<strong>de</strong>los basea<strong>do</strong>s em Thevenin, em Impedância, e em Runtime. Na<br />
Tabela 2.1, é apresentada uma comparação entre estas categorias e a seguir é realizada a<br />
sua <strong>de</strong>scrição.<br />
Tabela 2.1: Comparação entre as categorias [1].<br />
Capacida<strong>de</strong> <strong>de</strong> previsão Mo<strong>de</strong>lo basea<strong>do</strong> Mo<strong>de</strong>lo basea<strong>do</strong> Mo<strong>de</strong>lo basea<strong>do</strong><br />
em Thevenin em Impedância em Runtime<br />
corrente contínua (DC) Não Não Sim<br />
corrente variável (AC) Limita<strong>do</strong> Sim Não<br />
Transiente Sim Limita<strong>do</strong> Limita<strong>do</strong><br />
Runtime Não Não Sim<br />
Categoria 1: Mo<strong>de</strong>lo Basea<strong>do</strong> em Thevenin<br />
Em sua forma básica [1], o mo<strong>de</strong>lo basea<strong>do</strong> em Thevenin, apresenta<strong>do</strong> <strong>na</strong> Figura 2.8,<br />
utiliza um resistor Rseries e uma re<strong>de</strong> resistiva capacitiva paralela formada pelo resistor<br />
Rtransient e pelo capacitor Ctransient objetivan<strong>do</strong> prever a resposta da bateria para cargas<br />
transientes em um esta<strong>do</strong> particular <strong>de</strong> carga, consi<strong>de</strong>ran<strong>do</strong> constante a tensão <strong>de</strong> circuito<br />
aberto VOC(SOC). No entanto, este mo<strong>de</strong>lo não captura as variações <strong>de</strong> tensão da bateria<br />
no esta<strong>do</strong> estacionário (i.e., resposta DC), assim como a informação <strong>de</strong> tempo <strong>de</strong> vida [1].<br />
No mo<strong>de</strong>lo, o resistor Rself−discharge representa o fenômeno <strong>de</strong> auto-<strong>de</strong>scarga da bateria,<br />
o resistor Rseries representa a resistência inter<strong>na</strong> da bateria, e a re<strong>de</strong> RC representa o<br />
comportamento transiente da bateria e a constante <strong>de</strong> tempo para condições transientes.<br />
Figura 2.8: Mo<strong>de</strong>lo basea<strong>do</strong> em Thevenin [1].<br />
Conforme [1], os mo<strong>de</strong>los <strong>de</strong>riva<strong>do</strong>s <strong>do</strong> mo<strong>de</strong>lo basea<strong>do</strong> em Thevenin, possuem me-<br />
lhorias <strong>de</strong>vi<strong>do</strong> a adição <strong>de</strong> componentes para a predição <strong>do</strong> tempo <strong>de</strong> vida da bateria,<br />
mas ainda continuam com algumas <strong>de</strong>svantagens. Por exemplo, em [18] é utiliza<strong>do</strong> um<br />
capacitor variável em vez <strong>do</strong> VOC(SOC) para representar a tensão não-linear <strong>de</strong> circuito<br />
aberto, necessitan<strong>do</strong> <strong>de</strong> uma integral sobre a tensão da bateria para a obtenção <strong>do</strong> SOC.<br />
Em [19], é mo<strong>de</strong>lada a relação não-linear entre a tensão <strong>de</strong> circuito aberto e SOC, mas é<br />
<strong>de</strong>sconsi<strong>de</strong>ra<strong>do</strong> o comportamento transiente da bateria. Em [20] são necessárias equações
Capítulo 2. Revisão Bibliográfica 24<br />
matemáticas adicio<strong>na</strong>is para obtenção <strong>do</strong> SOC e predição <strong>do</strong> tempo <strong>de</strong> vida da bateria.<br />
Portanto, existem várias <strong>de</strong>rivações <strong>do</strong> mo<strong>de</strong>lo basea<strong>do</strong> em Thevenin, mas nenhuma <strong>de</strong>las<br />
po<strong>de</strong> predizer o tempo <strong>de</strong> vida da bateria <strong>de</strong> forma simples e prática [1].<br />
Categoria 2: Mo<strong>de</strong>lo Basea<strong>do</strong> em Impedância<br />
Mo<strong>de</strong>los basea<strong>do</strong>s em impedância, conforme apresenta<strong>do</strong> <strong>na</strong> Figura 2.9, empregam o<br />
méto<strong>do</strong> <strong>de</strong> espectroscopia <strong>de</strong> impedância eletroquímica para obter um mo<strong>de</strong>lo equivalente<br />
<strong>de</strong> impedância AC no <strong>do</strong>mínio <strong>de</strong> frequência e usam uma re<strong>de</strong> equivalente complicada<br />
(Zac) para ajustar o espectro <strong>de</strong> impedância, sen<strong>do</strong> este processo <strong>de</strong> ajuste difícil e não<br />
intuitivo [1]. Além disto, estes mo<strong>de</strong>los funcio<strong>na</strong>m somente para um SOC constante e<br />
uma temperatura <strong>de</strong>finida, por isto não po<strong>de</strong>m prever a resposta DC, ou o tempo <strong>de</strong><br />
vida da bateria [1]. O mo<strong>de</strong>lo apresenta<strong>do</strong> <strong>na</strong> Figura 2.9 é composto pela combi<strong>na</strong>ção <strong>do</strong><br />
resistor Rseries e <strong>do</strong> indutor Lseries para representar a resistência inter<strong>na</strong> da bateria, e a<br />
impedância Zac para mo<strong>de</strong>lar o equivalente eletroquímico da bateria. A tensão VOC(SOC)<br />
representa o SOC da bateria que é mo<strong>de</strong>lada como uma fonte <strong>de</strong> tensão [21].<br />
Figura 2.9: Mo<strong>de</strong>lo basea<strong>do</strong> em Impedância [21].<br />
Categoria 3: Mo<strong>de</strong>lo Basea<strong>do</strong> em Runtime<br />
Segun<strong>do</strong> [1], mo<strong>de</strong>los basea<strong>do</strong>s em Runtime, conforme apresenta<strong>do</strong> <strong>na</strong> Figura 2.10,<br />
utilizam uma re<strong>de</strong> <strong>de</strong> circuito complexa para simular o tempo <strong>de</strong> vida da bateria e a<br />
resposta DC para o caso <strong>de</strong> <strong>de</strong>scargas contínuas. No entanto, estes mo<strong>de</strong>los não po<strong>de</strong>m<br />
simular com acurácia, nem o tempo <strong>de</strong> vida, nem a resposta DC, para o caso <strong>de</strong> <strong>de</strong>scargas<br />
variáveis [1].<br />
Figura 2.10: Mo<strong>de</strong>lo basea<strong>do</strong> em Runtime [6].
Capítulo 2. Revisão Bibliográfica 25<br />
Na Figura 2.10, observa-se que o mo<strong>de</strong>lo basea<strong>do</strong> em Runtime é composto por três<br />
partes. Na primeira parte (Figura 2.10 (a)) é apresenta<strong>do</strong> o comportamento transiente<br />
da bateria, há um resistor Rtransient e um capacitor Ctransient. Na segunda parte (Figura<br />
2.10 (b)) é apresentada a resistência <strong>de</strong> auto-<strong>de</strong>scarga da bateria, a capacida<strong>de</strong> total<br />
da bateria e a queda <strong>de</strong> tensão para perdas inter<strong>na</strong>s, há um resistor Rself−discharge, um<br />
capacitor Ccapacity e uma tensão Vlost, a corrente Ibatt representa a corrente <strong>de</strong> <strong>de</strong>scarga<br />
da bateria. Na última parte (Figura 2.10 (c)) é apresentada a tensão nos termi<strong>na</strong>is e o<br />
SOC da bateria, há um resistor Rseries que mo<strong>de</strong>la a resistência inter<strong>na</strong> da bateria e uma<br />
tensão VOC(SOC) que representa o seu esta<strong>do</strong> <strong>de</strong> carga.<br />
Mo<strong>de</strong>lo para Predizer Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong><br />
O mo<strong>de</strong>lo para Predizer Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong>, é um mo<strong>de</strong>lo<br />
<strong>de</strong> bateria abrangente, intuitivo e <strong>de</strong> alta acurácia, que combi<strong>na</strong> as capacida<strong>de</strong>s tran-<br />
sientes <strong>do</strong>s mo<strong>de</strong>los basea<strong>do</strong>s em Thevenin, as características AC <strong>do</strong>s mo<strong>de</strong>los basea<strong>do</strong>s<br />
em Impedância, e a informação <strong>de</strong> tempo <strong>de</strong> vida <strong>do</strong>s mo<strong>de</strong>los basea<strong>do</strong>s em Runtime [1].<br />
No mo<strong>de</strong>lo elétrico apresenta<strong>do</strong> da Figura 2.11, o capacitor Ccapacity e a fonte <strong>de</strong><br />
corrente controlada mo<strong>de</strong>lam a capacida<strong>de</strong>, o SOC, e o tempo <strong>de</strong> vida da bateria. A<br />
re<strong>de</strong> RC simula a resposta transiente. A tensão gerada pela fonte controlada é usada<br />
para relacio<strong>na</strong>r o SOC com a tensão <strong>de</strong> circuito aberto VOC. Este mo<strong>de</strong>lo prevê o tempo<br />
<strong>de</strong> vida da bateria, o esta<strong>do</strong> estacionário e a resposta transiente <strong>de</strong> forma acurada. Ele<br />
também captura todas as características elétricas e dinâmicas da bateria como: capacida<strong>de</strong><br />
utilizável, tensão em circuito aberto e resposta transiente, que serão <strong>de</strong>scritas a seguir.<br />
Figura 2.11: Mo<strong>de</strong>lo para Predizer Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong> [1].<br />
Capacida<strong>de</strong> Utilizável<br />
A capacida<strong>de</strong> utilizável é a energia extraída quan<strong>do</strong> uma bateria é <strong>de</strong>scarregada a partir<br />
<strong>de</strong> um esta<strong>do</strong> carrega<strong>do</strong> para uma tensão fi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga, sen<strong>do</strong> que esta capacida<strong>de</strong><br />
diminui conforme aumenta o número <strong>de</strong> ciclos, a corrente <strong>de</strong> <strong>de</strong>scarga e o tempo <strong>de</strong><br />
armaze<strong>na</strong>mento (Figura 2.12(a), (c) e (d)) [1]. Esta capacida<strong>de</strong> utilizável aumenta com o
Capítulo 2. Revisão Bibliográfica 26<br />
aumento <strong>de</strong> temperatura (Figura 2.12(b)). O fenômeno da capacida<strong>de</strong> utilizável po<strong>de</strong> ser<br />
mo<strong>de</strong>la<strong>do</strong> por um capacitor carrega<strong>do</strong>, um resistor <strong>de</strong> auto-<strong>de</strong>scarga Rself−discharge e um<br />
resistor equivalente (i.e., a soma <strong>de</strong> Rseries, Rtransient_S e Rtransient_L).<br />
Figura 2.12: Curvas características <strong>de</strong> capacida<strong>de</strong> utilizável <strong>de</strong> baterias [1].<br />
Um capacitor carrega<strong>do</strong> representa a carga total armaze<strong>na</strong>da <strong>na</strong> bateria, por converter<br />
a capacida<strong>de</strong> nomi<strong>na</strong>l da bateria <strong>de</strong> Ah, para carga em Coulomb. O valor <strong>do</strong> capacitor<br />
Ccapacity é da<strong>do</strong> por:<br />
Ccapacity = 3600·Capacity·f1(ciclo)·f2(temp), (2.9)<br />
on<strong>de</strong>: Capacity é a capacida<strong>de</strong> nomi<strong>na</strong>l em Ah e f1(ciclo) e f2(temperatura) são fatores<br />
<strong>de</strong> correção <strong>de</strong>pen<strong>de</strong>ntes <strong>do</strong> número <strong>de</strong> ciclos e da temperatura. Ao <strong>de</strong>finir a tensão<br />
inicial VSOC em Ccapacity igual a 1 V ou 0 V , a bateria é inicializada em seu esta<strong>do</strong><br />
totalmente carregada (SOC <strong>de</strong> 100%), ou totalmente <strong>de</strong>scarregada (SOC <strong>de</strong> 0%). Em<br />
outras palavras, a tensão VSOC representa o SOC da bateria quantitativamente.<br />
A variação da capacida<strong>de</strong> utilizável <strong>de</strong>pen<strong>de</strong>nte da corrente, apresentada <strong>na</strong> Figura<br />
2.12(c), ocorre a partir <strong>de</strong> diferentes valores <strong>de</strong> SOC no fi<strong>na</strong>l da <strong>de</strong>scarga para diferentes<br />
correntes, <strong>de</strong>vi<strong>do</strong> a diferentes quedas <strong>de</strong> tensão sobre o resistor interno (i.e., a soma <strong>de</strong><br />
Rseries, Rtransient_S e Rtransient_L), e a mesma tensão fi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga. Quan<strong>do</strong> a bateria<br />
está sen<strong>do</strong> carregada ou <strong>de</strong>scarregada, a fonte <strong>de</strong> corrente controlada Ibatt é utilizada para<br />
carregar ou <strong>de</strong>scarregar o capacitor Ccapacity, <strong>de</strong> mo<strong>do</strong> que o SOC, representa<strong>do</strong> pela tensão<br />
VSOC, muda di<strong>na</strong>micamente. Portanto, o tempo <strong>de</strong> vida da bateria é obti<strong>do</strong> quan<strong>do</strong> a<br />
tensão da bateria atinge a tensão fi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga.<br />
O resistor <strong>de</strong> auto-<strong>de</strong>scarga Rself−discharge é usa<strong>do</strong> para caracterizar a perda <strong>de</strong> energia<br />
<strong>de</strong> auto-<strong>de</strong>scarga quan<strong>do</strong> as baterias são armaze<strong>na</strong>das por um longo tempo [1]. Teori-<br />
camente, este resistor é uma função <strong>de</strong> SOC, temperatura e frequentemente número <strong>de</strong>
Capítulo 2. Revisão Bibliográfica 27<br />
ciclos. Ele po<strong>de</strong> ser simplifica<strong>do</strong> como um resistor <strong>de</strong> alto valor ou ignora<strong>do</strong>, <strong>de</strong> acor<strong>do</strong><br />
com a curva <strong>de</strong> retenção <strong>de</strong> capacida<strong>de</strong> apresentada <strong>na</strong> Figura 2.12(d). Conforme po<strong>de</strong><br />
ser observa<strong>do</strong> <strong>na</strong> Figura 2.12(d), a capacida<strong>de</strong> utilizável diminui lentamente com o tempo,<br />
quan<strong>do</strong> nenhum circuito está conecta<strong>do</strong> à bateria.<br />
Tensão em Circuito Aberto<br />
A tensão <strong>de</strong> circuito aberto VOC é alterada para diferentes níveis <strong>de</strong> capacida<strong>de</strong>, i.e.,<br />
SOC, conforme apresenta<strong>do</strong> <strong>na</strong> Figura 2.13. A relação não-linear entre a tensão <strong>de</strong> circuito<br />
aberto e o SOC é importante para ser incluída no mo<strong>de</strong>lo. Assim, a fonte <strong>de</strong> tensão con-<br />
trolada VOC(VSOC) é utilizada para representar esta relação. A tensão <strong>de</strong> circuito aberto<br />
é normalmente medida como a tensão termi<strong>na</strong>l <strong>de</strong> circuito aberto no esta<strong>do</strong> estacionário<br />
em vários pontos SOC [1].<br />
Figura 2.13: Tensão em circuito aberto em relação ao SOC [1].<br />
Resposta Transiente<br />
Em um evento <strong>de</strong> <strong>de</strong>grau <strong>de</strong> corrente (i.e., mudança brusca <strong>de</strong> um esta<strong>do</strong> <strong>de</strong> corrente<br />
em um curto espaço <strong>de</strong> tempo), a tensão da bateria Vbatt respon<strong>de</strong> lentamente, conforme<br />
apresenta<strong>do</strong> <strong>na</strong> Figura 2.14. Por isto, esta resposta transiente é caracterizada pela re<strong>de</strong><br />
RC <strong>de</strong>stacada <strong>na</strong> Figura 2.11. A re<strong>de</strong> RC consiste <strong>de</strong> um resistor Rseries e duas re<strong>de</strong>s<br />
paralelas RC compostas pelo resistor Rtransient_S e o capacitor Ctransient_S, e pelo resistor<br />
Rtransient_L e o capacitor Ctransient_L. O resistor Rseries é responsável pela queda <strong>de</strong> tensão<br />
instantânea <strong>de</strong> resposta ao <strong>de</strong>grau <strong>de</strong> corrente. Os resistores Rtransient_S e Rtransient_L, e<br />
os capacitores Ctransient_S e Ctransient_L, são responsáveis pelas constantes <strong>de</strong> tempo, longa<br />
e curta, <strong>de</strong> resposta ao <strong>de</strong>grau <strong>de</strong> corrente, mostradas pelos <strong>do</strong>is círculos pontilha<strong>do</strong>s <strong>na</strong><br />
Figura 2.14. Com base em numerosas curvas experimentais, usar duas constantes <strong>de</strong><br />
tempo RC é a melhor opção entre precisão e complexida<strong>de</strong>, pois elas mantém erros <strong>de</strong>ntro<br />
<strong>de</strong> 1 mV para to<strong>do</strong>s os ajustes <strong>de</strong> curva [1].<br />
Teoricamente, to<strong>do</strong>s os parâmetros <strong>de</strong>ste mo<strong>de</strong>lo são funções multivariáveis <strong>de</strong> SOC,<br />
corrente, temperatura e número <strong>de</strong> ciclos. No entanto, <strong>de</strong>ntro <strong>de</strong> um erro <strong>de</strong> tolerân-<br />
cia, alguns parâmetros po<strong>de</strong>m ser simplifica<strong>do</strong>s [1]. Por exemplo, uma bateria <strong>de</strong> baixa<br />
capacida<strong>de</strong> em uma aplicação <strong>de</strong> temperatura constante, po<strong>de</strong> ter seus efeitos <strong>de</strong> tempe-
Capítulo 2. Revisão Bibliográfica 28<br />
Figura 2.14: Resposta transiente para um evento <strong>de</strong> <strong>de</strong>grau <strong>de</strong> corrente [1].<br />
ratura <strong>de</strong>sconsi<strong>de</strong>ra<strong>do</strong>s, assim como uma bateria frequentemente utilizada po<strong>de</strong> também<br />
ter a sua taxa <strong>de</strong> auto-<strong>de</strong>scarga (5% por mês [1]) <strong>de</strong>sconsi<strong>de</strong>rada.<br />
A partir da análise da Figura 2.11, po<strong>de</strong>-se concluir que em condições <strong>de</strong> esta<strong>do</strong> esta-<br />
cionário, os capacitores Ctransient_S e Ctransient_L funcio<strong>na</strong>m como um circuito aberto para<br />
DC, pois eles estão totalmente carrega<strong>do</strong>s, e neste caso, oferecem uma alta resistência<br />
à DC. Além disso, a resistência <strong>de</strong> auto-<strong>de</strong>scarga Rself−discharge po<strong>de</strong> ser <strong>de</strong>sconsi<strong>de</strong>rada<br />
[1]. Sen<strong>do</strong> assim, a tensão da bateria Vbatt para a condição <strong>de</strong> esta<strong>do</strong> estacionário é dada<br />
por:<br />
Vbatt = VOC − Ibatt(Rseries + Rtransient_S + Rtransient_L), (2.10)<br />
on<strong>de</strong>: VOC representa a tensão em circuito aberto da bateria, Ibatt representa a corrente da<br />
bateria, Rseries representa a sua resistência inter<strong>na</strong>, Rtransient_S representa a resistência<br />
transiente <strong>de</strong> curta duração e Rtransient_L representa a resistência transiente <strong>de</strong> longa<br />
duração.<br />
Quan<strong>do</strong> o mesmo mo<strong>de</strong>lo é consi<strong>de</strong>ra<strong>do</strong> em condições transientes (i.e., primeiros 30<br />
s [6]), on<strong>de</strong> uma mudança rápida <strong>na</strong> carga acontece, o efeito da capacitância ocorre nos<br />
capacitores Ctransient_S e Ctransient_L. Neste caso, os capacitores vão se comportar como<br />
um curto-circuito para o momento transitório, até que eles fiquem totalmente carrega<strong>do</strong>s.<br />
Sen<strong>do</strong> assim, a tensão da bateria Vbatt para condições transientes é dada por:<br />
Vbatt = VOC − IbattRseries −<br />
1<br />
Ctransient_S + Ctransient_L<br />
∫ 30<br />
0<br />
Ibattdt, (2.11)<br />
on<strong>de</strong>: VOC representa a tensão em circuito aberto da bateria, Ibatt representa a corrente da<br />
bateria, Rseries representa a sua resistência inter<strong>na</strong>, Ctransient_S representa a capacitância<br />
transiente <strong>de</strong> curta duração e Ctransient_L representa a capacitância transiente <strong>de</strong> longa<br />
duração. De acor<strong>do</strong> com [1], é importante <strong>de</strong>stacar que to<strong>do</strong>s os parâmetros presentes <strong>na</strong>s<br />
equações (2.10) e (2.11) são expressos em função <strong>do</strong> esta<strong>do</strong> <strong>de</strong> carga (SOC ) da bateria.<br />
Apesar <strong>do</strong> mo<strong>de</strong>lo para Predizer Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong> apre-<br />
sentar uma boa acurácia, pois ele engloba as características <strong>do</strong>s mo<strong>de</strong>los elétricos basea<strong>do</strong>s
Capítulo 2. Revisão Bibliográfica 29<br />
em Thevenin, Impedância e Runtime, já <strong>de</strong>scritos neste capítulo, o seu processo <strong>de</strong> ex-<br />
tração <strong>de</strong> parâmetros não é prático. Em [1], para aplicar este mo<strong>de</strong>lo <strong>na</strong> predição <strong>do</strong> tempo<br />
<strong>de</strong> vida <strong>de</strong> uma bateria <strong>de</strong> íon lítio polímero mo<strong>de</strong>lo PL-383562, foi necessária a realização<br />
<strong>de</strong> quarenta testes experimentais, sen<strong>do</strong> utilizadas quatro correntes <strong>de</strong> <strong>de</strong>scarga pulsantes<br />
(i.e., 80 mA, 160 mA, 320 mA e 640mA). Para cada uma <strong>de</strong>stas correntes <strong>de</strong> <strong>de</strong>scarga<br />
pulsantes foram realiza<strong>do</strong>s <strong>de</strong>z testes, obten<strong>do</strong>-se então, <strong>de</strong>z curvas reais <strong>de</strong> <strong>de</strong>scarga para<br />
cada corrente. Logo em seguida, para as <strong>de</strong>z curvas reais <strong>de</strong> <strong>de</strong>scarga obtidas para cada<br />
uma das correntes consi<strong>de</strong>radas, foi escolhida a curva média <strong>de</strong> <strong>de</strong>scarga, a partir da qual<br />
foram medi<strong>do</strong>s os valores <strong>do</strong>s parâmetros (VOC, Rseries, Rtransient_S, Rtransiente_L, Ctransient_S<br />
e Ctransiente_L) em diferentes pontos <strong>de</strong> SOC. Como resulta<strong>do</strong> <strong>de</strong>sta extração <strong>de</strong> parâ-<br />
metros, foram obti<strong>do</strong>s seis gráficos conforme apresenta<strong>do</strong>s <strong>na</strong> Figura 2.15, on<strong>de</strong> cada um<br />
representa quatro comportamentos <strong>de</strong> um <strong>do</strong>s 6 (seis) parâmetros em função <strong>do</strong> SOC, já<br />
que os parâmetros foram medi<strong>do</strong>s, consi<strong>de</strong>ran<strong>do</strong>-se quatro correntes <strong>de</strong> <strong>de</strong>scarga.<br />
A partir <strong>de</strong> interpolações, foram encontradas sete funções para representar as curvas<br />
<strong>do</strong>s parâmetros, conforme apresentadas a seguir <strong>na</strong>s equações (2.12)-(2.17). Cada função<br />
representa uma média das quatro funções obtidas a partir das correntes <strong>de</strong> <strong>de</strong>scarga<br />
pulsantes <strong>de</strong> 80 mA, 160 mA, 320 mA e 640mA, em cada gráfico. Na verda<strong>de</strong>, estas<br />
sete funções representam os sete parâmetros <strong>do</strong> mo<strong>de</strong>lo. Desta forma, é possível observar<br />
que os parâmetros <strong>do</strong> mo<strong>de</strong>lo são não-lineares, pois as funções que os representam são<br />
não-lineares. Os parâmetros neste caso são consi<strong>de</strong>ra<strong>do</strong>s <strong>de</strong>pen<strong>de</strong>ntes somente <strong>do</strong> SOC<br />
da bateria, já que não nota-se nos gráficos uma diferença consi<strong>de</strong>rável entre as curvas <strong>de</strong><br />
parâmetros para diferentes correntes <strong>de</strong> <strong>de</strong>scarga. Isto indica que os parâmetros po<strong>de</strong>m<br />
ser consi<strong>de</strong>ra<strong>do</strong>s in<strong>de</strong>pen<strong>de</strong>ntes das correntes <strong>de</strong> <strong>de</strong>scarga [1].<br />
VOC(SOC) = −1, 031 exp −35SOC +3, 685 + 0, 2156SOC − 0, 1178SOC 2 + 0, 3201SOC 3 .<br />
(2.12)<br />
RSeries(SOC) = 0, 1562 exp −24,37SOC +0, 07446. (2.13)<br />
RT ransient_S(SOC) = 0, 3208 exp −29,14SOC +0, 04669. (2.14)<br />
CT ransient_S(SOC) = −752, 9 exp −13,51SOC +703, 6. (2.15)<br />
RT ransient_L(SOC) = 6, 603 exp −155,2SOC +0, 04984. (2.16)<br />
CT ransient_L(SOC) = −6056 exp −27,12SOC +4475. (2.17)<br />
O processo <strong>de</strong> extração <strong>de</strong> parâmetros <strong>do</strong> mo<strong>de</strong>lo para Predizer Runtime e Caracte-<br />
rísticas V-I <strong>de</strong> uma <strong>Bateria</strong> foi <strong>de</strong>scrito <strong>de</strong> forma sucinta, maiores <strong>de</strong>talhes po<strong>de</strong>m ser<br />
encontra<strong>do</strong>s em [1]. Conforme já mencio<strong>na</strong><strong>do</strong> para a utilização <strong>do</strong> mo<strong>de</strong>lo para Predizer<br />
Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong>, <strong>de</strong>vem ser realiza<strong>do</strong>s, obrigatoriamente,
Capítulo 2. Revisão Bibliográfica 30<br />
Figura 2.15: Parâmetros extraí<strong>do</strong>s da bateria <strong>de</strong> íon lítio polímero [1].<br />
testes experimentais com a bateria a ser simulada, no processo <strong>de</strong> extração <strong>de</strong> parâmetros<br />
<strong>do</strong> mo<strong>de</strong>lo. Neste processo, é necessária a utilização <strong>de</strong> quatro curvas experimentais <strong>de</strong>
Capítulo 2. Revisão Bibliográfica 31<br />
<strong>de</strong>scarga, o que totaliza a necessida<strong>de</strong> <strong>de</strong> quarenta testes experimentais, e ainda, a uti-<br />
lização <strong>de</strong> uma ferramenta <strong>de</strong> ajuste <strong>de</strong> curva para concluir o processo <strong>de</strong> <strong>de</strong>finição <strong>do</strong>s<br />
parâmetros <strong>do</strong> mo<strong>de</strong>lo.<br />
No entanto, o principal objetivo <strong>de</strong>ste trabalho é utilizar um mo<strong>de</strong>lo elétrico <strong>de</strong><br />
<strong>de</strong>scarga que seja prático e <strong>de</strong> fácil implementação, com relação a extração <strong>de</strong> seus parâme-<br />
tros <strong>de</strong> configuração, <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> baterias que alimentam dispositivos<br />
móveis. Portanto, é <strong>de</strong>seja<strong>do</strong> um mo<strong>de</strong>lo elétrico que não necessite, obrigatoriamente, <strong>de</strong><br />
testes experimentais para a obtenção <strong>de</strong> seus parâmetros <strong>de</strong> configuração, ou ainda, que<br />
no caso da necessida<strong>de</strong> <strong>de</strong>stes testes, seja realiza<strong>do</strong> o menor número possível <strong>do</strong>s mesmos.<br />
Logicamente, são <strong>de</strong>sejadas todas estas características mencio<strong>na</strong>das em conjunto com<br />
uma boa acurácia <strong>do</strong> mo<strong>de</strong>lo, que po<strong>de</strong>m gerar uma otimização <strong>de</strong> tempo e <strong>de</strong> resulta<strong>do</strong>s,<br />
especialmente nos casos <strong>de</strong> <strong>de</strong>senvolvimento <strong>de</strong> novos projetos.<br />
Mo<strong>de</strong>lo Battery<br />
O mo<strong>de</strong>lo elétrico Battery é um mo<strong>de</strong>lo <strong>de</strong> bateria parametriza<strong>do</strong>, genérico e dinâmico, que<br />
representa os mais populares tipos <strong>de</strong> baterias recarregáveis (i.e., íon lítio, níquel cádmio,<br />
níquel metal hidreto e chumbo áci<strong>do</strong>). Com este mo<strong>de</strong>lo elétrico é possível simular o<br />
comportamento <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> um tipo específico <strong>de</strong> bateria, mediante a introdução <strong>de</strong><br />
parâmetros no mo<strong>de</strong>lo, tais como: tipo <strong>de</strong> bateria, resistência inter<strong>na</strong>, capacida<strong>de</strong> máxima,<br />
entre outros [22]. Seu diagrama esquemático é apresenta<strong>do</strong> <strong>na</strong> Figura 2.16.<br />
Figura 2.16: Diagrama esquemático <strong>do</strong> mo<strong>de</strong>lo Battery <strong>do</strong> MatLab/Simulink [22].<br />
Nas simulações realizadas a partir <strong>do</strong> emprego <strong>do</strong> mo<strong>de</strong>lo elétrico Battery, é possível<br />
observar uma boa resposta transiente e uma acurácia satisfatória nos tempos <strong>de</strong> vida<br />
simula<strong>do</strong>s, o que indica que o mo<strong>de</strong>lo possui resposta transiente como os mo<strong>de</strong>los elétricos
Capítulo 2. Revisão Bibliográfica 32<br />
basea<strong>do</strong>s em Thevenin e capacida<strong>de</strong> <strong>de</strong> previsão <strong>do</strong> tempo <strong>de</strong> vida, como os mo<strong>de</strong>los<br />
basea<strong>do</strong>s em Runtime.<br />
O mo<strong>de</strong>lo elétrico Battery, está presente <strong>na</strong> ferramenta computacio<strong>na</strong>l MatLab/Simulink,<br />
sen<strong>do</strong> que o MatLab (Matrix Laboratory) é um software, on<strong>de</strong> os elementos básicos <strong>de</strong> tra-<br />
balho são matrizes, no qual problemas po<strong>de</strong>m ser expressos em notação próxima a notação<br />
matemática e solucio<strong>na</strong><strong>do</strong>s por meio <strong>de</strong> cálculos computacio<strong>na</strong>is eficientes e confiáveis [23].<br />
As tarefas típicas realizadas pelo MatLab incluem:<br />
- Cálculo matemático e computacio<strong>na</strong>l;<br />
- Aquisição <strong>de</strong> da<strong>do</strong>s;<br />
- Análise <strong>de</strong> da<strong>do</strong>s;<br />
- Visualização <strong>de</strong> resulta<strong>do</strong>s;<br />
- Mo<strong>de</strong>lagem, simulação e prototipação;<br />
- Desenvolvimento <strong>de</strong> aplicativos.<br />
O Simulink é uma extensão <strong>do</strong> MatLab, on<strong>de</strong> po<strong>de</strong>m ser cria<strong>do</strong>s mo<strong>de</strong>los para a<strong>na</strong>lisar<br />
e controlar o comportamento <strong>de</strong> sistemas dinâmicos, a partir da implementação em diagra-<br />
mas <strong>de</strong> blocos. No aplicativo Simulink, existe uma biblioteca chamada SimPowerSystems<br />
que possui ferramentas para mo<strong>de</strong>lar e simular a geração, transmissão, distribuição e con-<br />
sumo <strong>de</strong> energia elétrica, e é nesta biblioteca que encontra-se o mo<strong>de</strong>lo elétrico Battery.<br />
No próximo capítulo, o mo<strong>de</strong>lo elétrico Battery é apresenta<strong>do</strong> e <strong>de</strong>scrito mais <strong>de</strong>talha-<br />
damente, em virtu<strong>de</strong> <strong>de</strong>le ser o mo<strong>de</strong>lo elétrico utiliza<strong>do</strong> neste trabalho, objetivan<strong>do</strong>-se<br />
verificar se é um mo<strong>de</strong>lo a<strong>de</strong>qua<strong>do</strong> para a predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> baterias que<br />
alimentam dispositivos móveis.
Capítulo 3<br />
Mo<strong>de</strong>lo Elétrico Battery<br />
Neste capítulo é apresenta<strong>do</strong> o mo<strong>de</strong>lo elétrico Battery, utiliza<strong>do</strong> neste trabalho para a<br />
predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> baterias utilizadas em dispositivos móveis. Este mo<strong>de</strong>lo foi<br />
escolhi<strong>do</strong> pela praticida<strong>de</strong> no processo <strong>de</strong> extração <strong>de</strong> parâmetros, que po<strong>de</strong>m ser obti<strong>do</strong>s<br />
diretamente a partir <strong>de</strong> uma única curva real <strong>de</strong> <strong>de</strong>scarga da bateria em conjunto com<br />
os da<strong>do</strong>s <strong>de</strong> seu datasheet, caso já exista, no datasheet, uma curva real <strong>de</strong> <strong>de</strong>scarga da<br />
bateria, não é necessária a realização <strong>de</strong> testes experimentais. A não realização <strong>de</strong> testes<br />
experimentais, vai proporcio<strong>na</strong>r uma otimização <strong>de</strong> tempo importante, principalmente<br />
quan<strong>do</strong> esta vem juntamente com uma acurácia satisfatória <strong>do</strong> mo<strong>de</strong>lo.<br />
Este capítulo está organiza<strong>do</strong> como segue. Na Seção 3.1 são apresenta<strong>do</strong>s as hipóte-<br />
ses <strong>do</strong> mo<strong>de</strong>lo. Na Seção 3.2 é apresentada a equação matemática <strong>do</strong> mo<strong>de</strong>lo elétrico<br />
Battery para baterias <strong>de</strong> íon lítio. Na Seção 3.3 é apresenta<strong>do</strong> o mo<strong>de</strong>lo elétrico Battery<br />
implementa<strong>do</strong> <strong>na</strong> forma <strong>de</strong> diagrama <strong>de</strong> blocos no MatLab/Simulink. Na Seção 3.4 são<br />
apresenta<strong>do</strong>s os parâmetros <strong>de</strong> configuração <strong>do</strong> mo<strong>de</strong>lo. E por fim, <strong>na</strong> Seção 3.5 é <strong>de</strong>scrito<br />
como o mo<strong>de</strong>lo será avalia<strong>do</strong>.<br />
3.1 Hipóteses <strong>do</strong> Mo<strong>de</strong>lo<br />
O mo<strong>de</strong>lo elétrico Battery possui as seguintes hipóteses <strong>de</strong> mo<strong>de</strong>lagem [22]:<br />
- A resistência inter<strong>na</strong> da bateria é consi<strong>de</strong>rada constante durante sua <strong>de</strong>scarga e não<br />
varia com a amplitu<strong>de</strong> da corrente;<br />
- A capacida<strong>de</strong> efetiva da bateria não se altera com as variações <strong>de</strong> amplitu<strong>de</strong> das<br />
correntes <strong>de</strong> <strong>de</strong>scarga;<br />
- A temperatura é negligenciada;<br />
- A auto-<strong>de</strong>scarga da bateria não é representada;<br />
- O mo<strong>de</strong>lo não consi<strong>de</strong>ra efeito memória (i.e., perda <strong>de</strong> capacida<strong>de</strong> da bateria, <strong>de</strong>vi<strong>do</strong><br />
ao fato <strong>de</strong>la não ser submetida a ciclos completos <strong>de</strong> carga/<strong>de</strong>scarga).<br />
33
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 34<br />
3.2 Equação Matemática <strong>do</strong> Mo<strong>de</strong>lo<br />
O mo<strong>de</strong>lo elétrico Battery possui a característica <strong>de</strong> simular os mais populares tipos <strong>de</strong><br />
baterias recarregáveis. Logo, para cada tipo <strong>de</strong> bateria há uma equação matemática que<br />
<strong>de</strong>screve o <strong>de</strong>caimento <strong>de</strong> tensão <strong>do</strong> sistema. Neste trabalho são utilizadas baterias <strong>de</strong> íon<br />
lítio, então, a equação matemática que representa o mo<strong>de</strong>lo elétrico Battery, para este<br />
tipo <strong>de</strong> bateria é dada a seguir [22]:<br />
f(it, i ∗ , i) = E0 − K Q<br />
Q − it i∗ − K Q<br />
Q − it it + Aexp(−Bit) , (3.1)<br />
on<strong>de</strong>: E0 é a tensão constante, K é a constante <strong>de</strong> polarização ou resistência <strong>de</strong> polari-<br />
zação, i* é a corrente dinâmica em baixa frequência, it é a capacida<strong>de</strong> extraída, Q é a<br />
capacida<strong>de</strong> máxima da bateria, A é a tensão exponencial e B é a capacida<strong>de</strong> exponencial.<br />
3.3 Mo<strong>de</strong>lo em Diagrama <strong>de</strong> Blocos<br />
Nesta seção é apresenta<strong>do</strong> o diagrama <strong>de</strong> blocos simplifica<strong>do</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico Battery,<br />
conforme a Figura 3.1. Observa-se que o diagrama possui um resistor e uma fonte <strong>de</strong><br />
tensão controlada, no bloco Mo<strong>de</strong>l Continuous, são utiliza<strong>do</strong>s os valores <strong>do</strong>s parâmetros da<br />
bateria para resolver as equações matemáticas, e a partir <strong>de</strong>stas resoluções são simuladas<br />
as curvas <strong>de</strong> <strong>de</strong>scarga da bateria. Na Figura 3.1, os pontos 1 e 2, <strong>do</strong> circuito em diagrama<br />
<strong>de</strong> blocos, <strong>de</strong>finem o local on<strong>de</strong> <strong>de</strong>ve ser aplica<strong>do</strong> o perfil <strong>de</strong> <strong>de</strong>scarga no mo<strong>de</strong>lo. Este<br />
perfil é aplica<strong>do</strong> a partir <strong>do</strong> componente Controlled Current Source, que será apresenta<strong>do</strong><br />
<strong>na</strong> Seção 5.4 <strong>do</strong> Capítulo 5.<br />
Na Figura 3.2, é apresenta<strong>do</strong> o subsistema Mo<strong>de</strong>l Continuous <strong>do</strong> diagrama <strong>de</strong> blocos<br />
simplifica<strong>do</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico Battery. Neste subsistema é realizada a simulação <strong>do</strong> fenô-<br />
meno <strong>de</strong> <strong>de</strong>scarga da bateria, consi<strong>de</strong>ran<strong>do</strong> os seus parâmetros informa<strong>do</strong>s, a partir da<br />
combi<strong>na</strong>ção <strong>de</strong> vários compontes como amplifica<strong>do</strong>res, satura<strong>do</strong>res, blocos <strong>de</strong> equações,<br />
entre outros. O resulta<strong>do</strong> fi<strong>na</strong>l da interação entre estes componentes, origi<strong>na</strong> quatro parce-<br />
las <strong>de</strong> tensão, que compõem o mo<strong>de</strong>lo matemático fi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga, equação (3.2), <strong>do</strong><br />
mo<strong>de</strong>lo elétrico Battery, para um <strong>de</strong>termi<strong>na</strong><strong>do</strong> tipo <strong>de</strong> bateria recarregável, que simula o<br />
<strong>de</strong>caimento <strong>de</strong> tensão <strong>do</strong> sistema. Observa-se que o subsistema Mo<strong>de</strong>l Continuous simula<br />
três si<strong>na</strong>is, que constituem a saída <strong>do</strong> mo<strong>de</strong>lo elétrico Battery, que são: tensão, SOC e<br />
perfil <strong>de</strong> <strong>de</strong>scarga da bateria.
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 35<br />
Figura 3.1: Diagrama em blocos simplifica<strong>do</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico Battery.<br />
Figura 3.2: Subsistema Mo<strong>de</strong>l Contiuous.<br />
3.4 Parâmetros <strong>do</strong> Mo<strong>de</strong>lo<br />
Os parâmetros <strong>do</strong> mo<strong>de</strong>lo elétrico Battery são informa<strong>do</strong>s para representar um tipo par-<br />
ticular <strong>de</strong> bateria, basea<strong>do</strong> em suas características <strong>de</strong> <strong>de</strong>scarga. Uma curva característica<br />
<strong>de</strong> <strong>de</strong>scarga é composta por três seções, conforme apresenta<strong>do</strong> <strong>na</strong> Figura 3.3.<br />
A primeira seção representa a queda <strong>de</strong> tensão exponencial quan<strong>do</strong> a bateria está<br />
totalmente carregada. Depen<strong>de</strong>n<strong>do</strong> da bateria, esta área é maior ou menor. A segunda<br />
seção representa a carga que po<strong>de</strong> ser extraída da bateria até a tensão <strong>de</strong>cair abaixo da<br />
tensão <strong>de</strong> cutoff. E a terceira seção representa a <strong>de</strong>scarga total da bateria, quan<strong>do</strong> a<br />
tensão diminui rapidamente.<br />
Para que o mo<strong>de</strong>lo elétrico Battery execute sua função <strong>de</strong> simulação <strong>de</strong> <strong>de</strong>scarga da<br />
bateria, e consequentemente possa predizer o seu tempo <strong>de</strong> vida, os seguintes parâmetros
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 36<br />
<strong>de</strong>vem ser informa<strong>do</strong>s [22]:<br />
Figura 3.3: Curva característica <strong>de</strong> <strong>de</strong>scarga [22].<br />
- Tipo <strong>de</strong> bateria: informa a tecnologia da bateria.<br />
- Tensão <strong>de</strong> cutoff (V ): representa o fi<strong>na</strong>l da zo<strong>na</strong> linear das características <strong>de</strong><br />
<strong>de</strong>scarga.<br />
- Capacida<strong>de</strong> típica (Ah): é a capacida<strong>de</strong> mínima efetiva da bateria.<br />
- SOC inicial (%): representa o esta<strong>do</strong> inicial <strong>de</strong> carga da bateria. O esta<strong>do</strong> <strong>de</strong><br />
100% indica que a bateria está totalmente carregada e 0% indica que ela está sem carga.<br />
Este parâmetro é usa<strong>do</strong> como condição inicial para a simulação.<br />
- Capacida<strong>de</strong> máxima (Ah): é a capacida<strong>de</strong> Q, quan<strong>do</strong> uma <strong>de</strong>scontinuida<strong>de</strong> ocorre<br />
<strong>na</strong> tensão da bateria. Segun<strong>do</strong> [22], este valor po<strong>de</strong> ser consi<strong>de</strong>ra<strong>do</strong> como aproximada-<br />
mente 105% da capacida<strong>de</strong> típica da bateria.<br />
- Tensão com bateria carregada (V ): é representada para uma <strong>de</strong>termi<strong>na</strong>da cor-<br />
rente <strong>de</strong> <strong>de</strong>scarga. Esta tensão da bateria totalmente carregada não é a tensão quan<strong>do</strong> a<br />
bateria não possui nenhum sistema conecta<strong>do</strong> a ela.<br />
- Corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga (A): é a corrente <strong>de</strong> <strong>de</strong>scarga a qual a bateria<br />
é submetida, e a partir <strong>de</strong>ste processo, é obtida uma curva característica <strong>de</strong> <strong>de</strong>scarga da<br />
bateria. Por exemplo, uma corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga típica para uma bateria Ni-MH<br />
<strong>de</strong> 1,5 Ah é 20% da capacida<strong>de</strong> típica, ou seja, 0,3 A [22].<br />
- Resistência inter<strong>na</strong> (Ohms): representa a resistência inter<strong>na</strong> da bateria. Esta<br />
resistência permanece constante durante a simulação, ou seja, não varia com a amplitu<strong>de</strong><br />
da corrente elétrica. Uma opção para este parâmetro é um valor genérico correspon<strong>de</strong>n<strong>do</strong><br />
a 1% <strong>do</strong> produto da tensão <strong>de</strong> cutoff pela capacida<strong>de</strong> típica da bateria [22].<br />
- Capacida<strong>de</strong> <strong>na</strong> tensão <strong>de</strong> cutoff (Ah): é a capacida<strong>de</strong> extraída da bateria até a<br />
sua tensão atingir o valor da tensão <strong>de</strong> cutoff.<br />
- Zo<strong>na</strong> exponencial (tensão, capacida<strong>de</strong>): a tensão exponencial e a capacida<strong>de</strong> ex-<br />
ponencial correspon<strong>de</strong>m ao fi<strong>na</strong>l da zo<strong>na</strong> exponencial da curva característica <strong>de</strong> <strong>de</strong>scarga.<br />
A tensão <strong>de</strong>ve estar entre a tensão <strong>de</strong> cutoff e a tensão com a bateria carregada, e a
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 37<br />
capacida<strong>de</strong> entre zero e a capacida<strong>de</strong> nomi<strong>na</strong>l.<br />
- Unida<strong>de</strong>: durante a simulação, representa a <strong>de</strong>finição da unida<strong>de</strong> no eixo x, que<br />
po<strong>de</strong> ser em segun<strong>do</strong>s (s) ou em Ampère-hora (Ah).<br />
- Tempo <strong>de</strong> resposta da bateria (s): este valor representa a dinâmica da tensão<br />
da bateria e po<strong>de</strong> ser observada no caso <strong>de</strong> uma corrente <strong>de</strong> <strong>de</strong>scarga <strong>do</strong> tipo <strong>de</strong>grau,<br />
conforme apresenta<strong>do</strong> <strong>na</strong> Figura 3.4.<br />
Figura 3.4: Tempo <strong>de</strong> resposta da bateria <strong>de</strong> 30 segun<strong>do</strong>s [22].<br />
No exemplo <strong>de</strong> simulação apresenta<strong>do</strong> <strong>na</strong> Figura 3.4, é possível perceber que o mo<strong>de</strong>lo<br />
consi<strong>de</strong>ra aspectos físicos importantes ocorri<strong>do</strong>s em um processo real <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> uma<br />
bateria. Quan<strong>do</strong> a bateria é submetida a um <strong>de</strong>grau <strong>de</strong> corrente, a sua tensão começa<br />
a <strong>de</strong>cair <strong>de</strong> forma exponencial e não instantaneamente, conforme po<strong>de</strong> ser observa<strong>do</strong><br />
no tempo <strong>de</strong> 4000 s, este comportamento exponencial da tensão ocorre por 30 s. No<br />
instante <strong>de</strong> 4100 s, ocorre uma aumento exponencial <strong>na</strong> tensão da bateria, que representa<br />
um efeito não-linear importante a ser consi<strong>de</strong>ra<strong>do</strong>, que é o efeito <strong>de</strong> recuperação. Este<br />
aumento exponencial <strong>de</strong> tensão dura por 30 s, e <strong>de</strong>pois a tensão estabiliza.<br />
3.4.1 Extração <strong>de</strong> Parâmetros<br />
Para que o mo<strong>de</strong>lo elétrico Battery execute sua função <strong>de</strong> simulação <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> uma<br />
bateria e consequentemente, <strong>de</strong> seu tempo <strong>de</strong> vida, os parâmetros apresenta<strong>do</strong>s <strong>na</strong> seção<br />
anterior <strong>de</strong>vem ser informa<strong>do</strong>s. Conforme já mencio<strong>na</strong><strong>do</strong> no início <strong>de</strong>ste capítulo, estes<br />
parâmetros po<strong>de</strong>m ser obti<strong>do</strong>s a partir <strong>de</strong> uma única curva característica <strong>de</strong> <strong>de</strong>scarga da<br />
bateria em conjunto com especificações encontradas em seu datasheet. Em alguns casos,
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 38<br />
esta curva característica <strong>de</strong> <strong>de</strong>scarga já está presente no datasheet, em outros é necessário<br />
obtê-la <strong>de</strong> forma experimental.<br />
Na Figura 3.5 é apresenta<strong>do</strong> um exemplo <strong>de</strong> uma curva característica <strong>de</strong> <strong>de</strong>scarga<br />
presente no datasheet <strong>de</strong> uma bateria <strong>de</strong> Ni-MH [22], com os seus respectivos parâme-<br />
tros, on<strong>de</strong>, o parâmetro c representa a tensão elétrica quan<strong>do</strong> a bateria está totalmente<br />
carregada; o parâmetro e representa a tensão e a capacida<strong>de</strong> da bateria no fi<strong>na</strong>l da zo<strong>na</strong><br />
exponencial da curva; o parâmetro a representa o fi<strong>na</strong>l da zo<strong>na</strong> linear das característi-<br />
cas <strong>de</strong> <strong>de</strong>scarga da bateria; o parâmetro b representa a capacida<strong>de</strong> máxima da bateria,<br />
quan<strong>do</strong> ocorre uma <strong>de</strong>scontinuida<strong>de</strong> <strong>na</strong> sua tensão; o parâmetro d representa a corrente<br />
nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga da bateria, a partir da qual esta curva característica <strong>de</strong> <strong>de</strong>scarga<br />
foi obtida. O valor <strong>de</strong>sta corrente nomi<strong>na</strong>l para este tipo <strong>de</strong> bateria é normalmente 20%<br />
da capacida<strong>de</strong> típica da bateria [22]. Estes parâmetros são aproxima<strong>do</strong>s e <strong>de</strong>pen<strong>de</strong>m da<br />
precisão <strong>do</strong>s pontos obti<strong>do</strong>s da curva real <strong>de</strong> <strong>de</strong>scarga [22]. Os parâmetros resistência<br />
inter<strong>na</strong> e capacida<strong>de</strong> típica, não foram extraí<strong>do</strong>s da curva característica <strong>de</strong> <strong>de</strong>scarga da<br />
bateria, sen<strong>do</strong> obti<strong>do</strong>s a partir <strong>do</strong> datasheet da mesma.<br />
Figura 3.5: Curva característica <strong>de</strong> <strong>de</strong>scarga da bateria <strong>de</strong> Ni-MH Pa<strong>na</strong>sonic HHR650D [22].<br />
Como resulta<strong>do</strong> da extração <strong>de</strong> parâmetros da curva característica <strong>de</strong> <strong>de</strong>scarga apre-<br />
sentada <strong>na</strong> Figura 3.5 e <strong>do</strong>s da<strong>do</strong>s <strong>do</strong> datasheet, os parâmetros obti<strong>do</strong>s são apresenta<strong>do</strong>s<br />
<strong>na</strong> Tabela 3.1.<br />
Após a obtenção <strong>de</strong> to<strong>do</strong>s os parâmetros necessários para o mo<strong>de</strong>lo elétrico Battery,<br />
é condição informá-los em uma interface gráfica, para que possam ser realizadas as simu-<br />
lações <strong>de</strong> <strong>de</strong>scarga da bateria. Esta extração <strong>de</strong> parâmetros da curva real <strong>de</strong> <strong>de</strong>scarga da<br />
bateria é necessária somente uma vez, para a calibração <strong>do</strong> mo<strong>de</strong>lo. Depois <strong>de</strong> feita esta<br />
calibração, basta submeter o mo<strong>de</strong>lo a diferentes tipos <strong>de</strong> perfis <strong>de</strong> <strong>de</strong>scarga, para que as<br />
simulações sejam realizadas.<br />
O registro <strong>do</strong>s parâmetros <strong>na</strong> interface <strong>do</strong> mo<strong>de</strong>lo po<strong>de</strong> ser realiza<strong>do</strong> a partir <strong>de</strong> duas
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 39<br />
Tabela 3.1: Parâmetros da bateria <strong>de</strong> níquel metal hidreto Pa<strong>na</strong>sonic HHR650D [22].<br />
Parâmetro Valor<br />
Resistência inter<strong>na</strong> 0,002 Ω<br />
Tensão cutoff 1,18 V<br />
Capacida<strong>de</strong> típica 6,5 Ah<br />
Capacida<strong>de</strong> máxima 7 Ah (5,38 h * 1,3 A)<br />
Tensão com bateria carregada 1,39 V<br />
Corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga 1,3 A<br />
Capacida<strong>de</strong> <strong>na</strong> tensão cutoff 6,25 Ah<br />
Tensão exponencial 1,28 V<br />
Capacida<strong>de</strong> exponencial 1,3 Ah<br />
formas diferentes. Uma forma é a genérica, on<strong>de</strong> são informa<strong>do</strong>s somente três parâmetros<br />
da bateria, a tensão <strong>de</strong> cutoff, a capacida<strong>de</strong> típica e o SOC inicial da bateria, os <strong>de</strong>mais,<br />
são preenchi<strong>do</strong>s <strong>de</strong> forma genérica pelo programa. Nesta forma <strong>de</strong> registro, a opção Use<br />
parameters based on battery type and nomi<strong>na</strong>l values <strong>de</strong>ve estar marcada, a interface <strong>de</strong>ste<br />
mo<strong>do</strong> <strong>de</strong> registro é apresentada <strong>na</strong> Figura 3.6. A outra forma <strong>de</strong> registro <strong>de</strong> parâmetros é<br />
a não genérica, que é a<strong>do</strong>tada neste trabalho, on<strong>de</strong> é informada uma quantida<strong>de</strong> maior <strong>de</strong><br />
parâmetros. Além <strong>do</strong>s parâmetros informa<strong>do</strong>s <strong>na</strong> forma genérica, é necessário o registro<br />
<strong>de</strong> parâmetros como capacida<strong>de</strong> máxima, resistência inter<strong>na</strong>, entre outros. Nesta forma<br />
<strong>de</strong> registro, a opção Use parameters based on Battery type and nomi<strong>na</strong>l values <strong>de</strong>ve estar<br />
<strong>de</strong>smarcada. Na Figura 3.7 é apresentada a sua interface.<br />
Tanto <strong>na</strong> Figura 3.6, quanto <strong>na</strong> Figura 3.7, observa-se que a interface possui três<br />
guias. A guia Parameters, on<strong>de</strong> são informa<strong>do</strong>s os parâmetros relevantes para a simula-<br />
ção <strong>de</strong> <strong>de</strong>scarga da bateria. A guia View Discharge Characteristics que permite a geração<br />
automática <strong>de</strong> curvas <strong>de</strong> <strong>de</strong>scarga da bateria simulada, mas sem a possibilida<strong>de</strong> <strong>de</strong> ma-<br />
nipulação <strong>do</strong>s da<strong>do</strong>s gera<strong>do</strong>s. E fi<strong>na</strong>lmente, a guia Battery Dy<strong>na</strong>mics, on<strong>de</strong> é <strong>de</strong>fini<strong>do</strong> o<br />
Figura 3.6: Interface <strong>do</strong> mo<strong>de</strong>lo elétrico Battery <strong>na</strong> forma genérica <strong>de</strong> registro [22].
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 40<br />
Figura 3.7: Interface <strong>do</strong> mo<strong>de</strong>lo elétrico Battery <strong>na</strong> forma não genérica <strong>de</strong> registro [22].<br />
tempo <strong>de</strong> resposta da bateria, ou seja, a dinâmica da tensão da bateria simulada.<br />
É importante <strong>de</strong>stacar que estas interfaces servem para a entrada <strong>de</strong> parâmetros <strong>do</strong><br />
mo<strong>de</strong>lo, e também para a simulação <strong>de</strong> uma ou mais curvas <strong>de</strong> <strong>de</strong>scarga, geradas a partir<br />
<strong>de</strong> <strong>de</strong>scargas contínuas, informadas diretamente <strong>na</strong> interface. Estas <strong>de</strong>scargas contínuas<br />
precisam ser informadas <strong>na</strong> guia View Discharge Characteristics. Mas para a simulação<br />
<strong>de</strong> curvas no caso <strong>de</strong> <strong>de</strong>scargas variáveis, e também para a possibilida<strong>de</strong> <strong>de</strong> manipulação<br />
<strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s, tanto <strong>na</strong>s simulações em <strong>de</strong>scargas contínuas quanto variáveis, é<br />
necessária a implementação <strong>de</strong> um circuito em diagrama <strong>de</strong> blocos no MatLab/Simulink,<br />
além da entrada <strong>de</strong> parâmetros <strong>na</strong> interface. Este circuito é apresenta<strong>do</strong> e <strong>de</strong>scrito <strong>de</strong><br />
forma mais <strong>de</strong>talhada no próximo capítulo.
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 41<br />
3.5 Saída <strong>do</strong> Mo<strong>de</strong>lo via Simulink<br />
A saída <strong>do</strong> mo<strong>de</strong>lo elétrico Battery via Simulink é uma matriz <strong>de</strong> or<strong>de</strong>m N ×3 que contém<br />
três si<strong>na</strong>is, o SOC, a corrente e a tensão. O SOC, que representa o esta<strong>do</strong> <strong>de</strong> carga da<br />
bateria, é da<strong>do</strong> por:<br />
(<br />
SOC = 100 1 − 1<br />
Q<br />
∫ t<br />
0<br />
)<br />
i(t)dt , (3.2)<br />
on<strong>de</strong>: Q é a capacida<strong>de</strong> máxima da bateria e i(t) é a corrente <strong>de</strong> <strong>de</strong>scarga. A corrente<br />
é a própria corrente <strong>de</strong> <strong>de</strong>scarga aplicada ao mo<strong>de</strong>lo, e a tensão, representa a tensão nos<br />
termi<strong>na</strong>is da bateria, que é obtida pela equação (3.1).<br />
A seguir, é apresenta<strong>do</strong> um exemplo <strong>do</strong> processo <strong>de</strong> extração <strong>de</strong> parâmetros, uti-<br />
lizan<strong>do</strong> a bateria <strong>de</strong> íon lítio Nokia BL-5F que possui a capacida<strong>de</strong> nomi<strong>na</strong>l <strong>de</strong> 950 mAh,<br />
e também, a matriz resultante da simulação <strong>de</strong> <strong>de</strong>scarga <strong>de</strong>sta bateria, com o emprego<br />
<strong>do</strong> mo<strong>de</strong>lo elétrico Battery. Esta bateria <strong>de</strong> íon lítio não possui datasheet, por isto foi<br />
necessária a realização <strong>de</strong> testes experimentais, para a obtenção <strong>de</strong> uma curva caracterís-<br />
tica <strong>de</strong> <strong>de</strong>scarga da bateria, on<strong>de</strong> foi consi<strong>de</strong>rada uma corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>de</strong><br />
150 mA, conforme apresenta<strong>do</strong> <strong>na</strong> Figura 3.8. Esta foi a curva utilizada para a calibração<br />
<strong>do</strong> mo<strong>de</strong>lo, já que o mesmo po<strong>de</strong> ser calibra<strong>do</strong> a partir <strong>de</strong> uma única curva <strong>de</strong> <strong>de</strong>scarga<br />
da bateria [22].<br />
Figura 3.8: Curva característica <strong>de</strong> <strong>de</strong>scarga para uma <strong>de</strong>scarga contínua <strong>de</strong> 150 mA <strong>na</strong> bateria<br />
<strong>de</strong> íon lítio Nokia BL-5F.<br />
A curva característica <strong>de</strong> <strong>de</strong>scarga da bateria é necessária para que seja possível a<br />
<strong>de</strong>finição <strong>de</strong> alguns parâmetros <strong>de</strong> configuração <strong>do</strong> mo<strong>de</strong>lo. Os parâmetros que não po<strong>de</strong>m<br />
ser obti<strong>do</strong>s a partir da curva, como por exemplo a resistência inter<strong>na</strong> da bateria, po<strong>de</strong>m<br />
ser <strong>de</strong>termi<strong>na</strong><strong>do</strong>s conforme sugeri<strong>do</strong> por [22]. Como resulta<strong>do</strong> <strong>de</strong>ste processo, os seguintes<br />
parâmetros foram <strong>de</strong>fini<strong>do</strong>s para a bateria <strong>de</strong> íon lítio Nokia BL-5F, conforme apresenta<strong>do</strong><br />
<strong>na</strong> Tabela 3.2.<br />
Quan<strong>do</strong> estes parâmetros são informa<strong>do</strong>s ao mo<strong>de</strong>lo elétrico Battery, o mesmo fica<br />
calibra<strong>do</strong> a partir da corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 150 mA e <strong>do</strong>s parâmetros extraí<strong>do</strong>s
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 42<br />
Tabela 3.2: Parâmetros da bateria <strong>de</strong> íon lítio Nokia BL-5F.<br />
Parâmetros Valores<br />
Resistência inter<strong>na</strong> 0,02945 Ω<br />
Tensão <strong>de</strong> cutoff 3,10 V<br />
Capacida<strong>de</strong> típica 0,950 Ah<br />
Capacida<strong>de</strong> máxima 0,9975 Ah<br />
Tensão com bateria carregada 4,21 V<br />
Corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga 0,15 A<br />
Capacida<strong>de</strong> <strong>na</strong> tensão cutoff 0,8129 Ah<br />
Tensão exponencial 3,7 V<br />
Capacida<strong>de</strong> exponencial 0,5032 Ah<br />
da curva característica <strong>de</strong> <strong>de</strong>scarga, origi<strong>na</strong>da por esta corrente nomi<strong>na</strong>l. Nesta situação<br />
<strong>de</strong> calibração <strong>do</strong> mo<strong>de</strong>lo, foi então realizada uma simulação com o mesmo, consi<strong>de</strong>ran<strong>do</strong><br />
a aplicação <strong>de</strong> uma corrente <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 950 mA. Como resulta<strong>do</strong> da simulação, a<br />
saída <strong>do</strong> mo<strong>de</strong>lo é dada pela matriz apresentada <strong>na</strong> Figura 3.9, on<strong>de</strong> a primeira colu<strong>na</strong><br />
representa o SOC da bateria, em percentual; a segunda colu<strong>na</strong> representa a corrente <strong>de</strong><br />
<strong>de</strong>scarga da bateria, em Ampére; e a terceira colu<strong>na</strong>, representa a tensão da bateria em<br />
Volts.<br />
Durante a simulação <strong>de</strong> <strong>de</strong>scarga pelo mo<strong>de</strong>lo, os valores <strong>de</strong> SOC e da tensão da bateria<br />
são atualiza<strong>do</strong>s a cada segun<strong>do</strong>, já que neste exemplo, as simulações foram realizadas a<br />
partir da utilização <strong>do</strong> méto<strong>do</strong> numérico <strong>de</strong> Bogacki Shampine, presente no Simulink, com<br />
um passo fixa<strong>do</strong> em 1 segun<strong>do</strong>. Logicamente, por tratar-se <strong>de</strong> uma simulação <strong>de</strong> <strong>de</strong>scarga,<br />
os valores <strong>de</strong> SOC e tensão da bateria vão diminuir com o passar <strong>do</strong> tempo, conforme<br />
po<strong>de</strong> ser observa<strong>do</strong> <strong>na</strong>s colu<strong>na</strong>s 1 e 3 da matriz. É importante <strong>de</strong>stacar, que o mo<strong>de</strong>lo<br />
recalcula os valores <strong>de</strong> SOC e tensão da bateria a cada segun<strong>do</strong>, durante to<strong>do</strong> o tempo <strong>de</strong><br />
simulação <strong>de</strong>fini<strong>do</strong>. Então, quanto maior este tempo <strong>de</strong> simulação, mais linhas a matriz<br />
vai possuir.<br />
3.6 Avaliação <strong>do</strong> Mo<strong>de</strong>lo<br />
Em um primeiro momento, o mo<strong>de</strong>lo elétrico Battery será avalia<strong>do</strong> a partir da comparação<br />
<strong>de</strong> seus resulta<strong>do</strong>s simula<strong>do</strong>s com os resulta<strong>do</strong>s experimentais obti<strong>do</strong>s a partir <strong>de</strong> uma<br />
plataforma <strong>de</strong> testes (i.e., testbed). Posteriormente, o mo<strong>de</strong>lo será avalia<strong>do</strong> a partir da<br />
comparação <strong>de</strong> seus resulta<strong>do</strong>s simula<strong>do</strong>s com os resulta<strong>do</strong>s simula<strong>do</strong>s <strong>do</strong> mo<strong>de</strong>lo elétrico<br />
para Predizer Runtime e características V-I <strong>de</strong> uma bateria.
Capítulo 3. Mo<strong>de</strong>lo Elétrico Battery 43<br />
Figura 3.9: Matriz N × 3 gerada após simulação com o mo<strong>de</strong>lo elétrico Battery.
Capítulo 4<br />
Ambiente <strong>de</strong> Simulação<br />
Neste capítulo são apresentadas a plataforma <strong>de</strong> testes (i.e., testbed), e a meto<strong>do</strong>logia<br />
a<strong>do</strong>tada para os processos <strong>de</strong> carga e <strong>de</strong>scarga das baterias utilizadas nos experimentos<br />
realiza<strong>do</strong>s. Observa-se que os resulta<strong>do</strong>s obti<strong>do</strong>s a partir <strong>do</strong>s experimentos, consi<strong>de</strong>ran<strong>do</strong><br />
baterias <strong>de</strong> íon lítio, são usa<strong>do</strong>s para comparação com os resulta<strong>do</strong>s das simulações obti<strong>do</strong>s<br />
a partir <strong>do</strong> mo<strong>de</strong>lo elétrico Battery.<br />
O restante <strong>de</strong>ste capítulo está organiza<strong>do</strong> como segue. Na Seção 5.1 é realizada uma<br />
<strong>de</strong>scrição da plataforma <strong>de</strong> testes utilizada neste trabalho. Na Seção 5.2 é apresentada<br />
a meto<strong>do</strong>logia utilizada para a realização <strong>do</strong>s testes experimentais. Na Seção 5.3 são<br />
apresenta<strong>do</strong>s os resulta<strong>do</strong>s experimentais obti<strong>do</strong>s a partir da plataforma. Na Seção 4.4 é<br />
realizada uma <strong>de</strong>scrição <strong>de</strong> como são efetuadas as simulações computacio<strong>na</strong>is. E por fim,<br />
<strong>na</strong> Seção 4.5 é realizada a análise <strong>do</strong>s resulta<strong>do</strong>s.<br />
4.1 Testbed<br />
Nesta seção é apresentada a plataforma <strong>de</strong> testes, também <strong>de</strong>nomi<strong>na</strong>da testbed, utilizada<br />
neste trabalho para a realização <strong>do</strong>s experimentos [24]. Este testbed, a partir da qual é<br />
possível capturar as curvas reais <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> uma bateria, é constituí<strong>do</strong> <strong>de</strong> três partes:<br />
(i) sistema <strong>de</strong> controle (software), (ii) hardware e (iii) bateria, conforme apresenta<strong>do</strong> <strong>na</strong><br />
Figura 4.1.<br />
O hardware realiza a comunicação com o computa<strong>do</strong>r e a administração <strong>do</strong>s módulos<br />
<strong>de</strong> sensoriamento e controle <strong>de</strong> <strong>de</strong>scarga. O software, <strong>de</strong>senvolvi<strong>do</strong> em C++, apresenta<br />
uma interface intuitiva para a informação <strong>do</strong>s parâmetros da bateria, o mesmo envia ao<br />
hardware as configurações <strong>do</strong> tipo <strong>de</strong> <strong>de</strong>scarga o qual a bateria <strong>de</strong>ve ser submetida. E<br />
por fim, a bateria, que neste trabalho é consi<strong>de</strong>rada uma bateria <strong>de</strong> íon lítio amplamente<br />
utilizadas em telefones celulares. A seguir o hardware e o software são <strong>de</strong>scritos mais<br />
<strong>de</strong>talhadamente.<br />
44
Capítulo 4. Ambiente <strong>de</strong> Simulação 45<br />
4.1.1 Hardware<br />
Figura 4.1: Diagrama <strong>do</strong> testbed [24].<br />
O hardware é composto por três placas eletrônicas, sen<strong>do</strong> que o principal componente<br />
presente em uma <strong>de</strong>las é um microcontrola<strong>do</strong>r, que possui a função <strong>de</strong> efetuar aquisições<br />
<strong>de</strong> temperatura e <strong>do</strong>s módulos <strong>de</strong> tensão e <strong>de</strong> corrente. O microcontrola<strong>do</strong>r, através <strong>de</strong> um<br />
algoritmo <strong>de</strong> controle Proporcio<strong>na</strong>l Integral (PI) faz o controle da corrente <strong>de</strong> <strong>de</strong>scarga<br />
a que é submetida a bateria, com referência aos coman<strong>do</strong>s <strong>do</strong> software. Mais <strong>de</strong>talhes<br />
po<strong>de</strong>m ser obti<strong>do</strong>s em [24]. Na Figura 4.2 é apresentada esta parte <strong>do</strong> testbed.<br />
4.1.2 Software<br />
Figura 4.2: Foto <strong>do</strong> hardware.<br />
O software possibilita a obtenção <strong>do</strong>s resulta<strong>do</strong>s gera<strong>do</strong>s pelos testes. Após o preenchi-<br />
mento <strong>de</strong> da<strong>do</strong>s (i.e., parâmetros e tipo <strong>de</strong> <strong>de</strong>scarga) <strong>na</strong> interface, o software administra<br />
o controle da <strong>de</strong>scarga a que a bateria é submetida. Além <strong>do</strong> software possuir recursos<br />
<strong>de</strong> proteção da bateria, no caso <strong>de</strong> problemas no sistema, ele possibilita também salvar
Capítulo 4. Ambiente <strong>de</strong> Simulação 46<br />
as imagens <strong>do</strong>s gráficos em Bitmap e os relatórios em formato texto. Na Figura 4.3 é<br />
apresentada a interface <strong>do</strong> software que é dividida em <strong>do</strong>is ca<strong>na</strong>is, o primeiro se encontra<br />
<strong>na</strong> meta<strong>de</strong> superior da interface, e o segun<strong>do</strong> está <strong>na</strong> meta<strong>de</strong> inferior. Devi<strong>do</strong> ao fato <strong>de</strong><br />
existirem <strong>do</strong>is ca<strong>na</strong>is, é possível realizar testes com duas baterias simultaneamente.<br />
Figura 4.3: Interface <strong>do</strong> software <strong>de</strong> controle <strong>do</strong> testbed.<br />
A seguir é realizada a <strong>de</strong>scrição <strong>do</strong>s blocos 1, 2, 3 e 4 <strong>de</strong>staca<strong>do</strong>s <strong>na</strong> interface apresen-<br />
tada <strong>na</strong> Figura 4.3 consi<strong>de</strong>ran<strong>do</strong> o ca<strong>na</strong>l I da mesma:<br />
1. Apresentação da duração <strong>do</strong> teste aplica<strong>do</strong> à bateria;<br />
2. Apresentação <strong>do</strong>s valores instantâneos obti<strong>do</strong>s da bateria;<br />
3. Configuração, controle e armaze<strong>na</strong>mento;<br />
4. Representação gráfica <strong>do</strong> teste aplica<strong>do</strong> à bateria.<br />
Esta <strong>de</strong>scrição <strong>do</strong>s blocos <strong>do</strong> ca<strong>na</strong>l I é válida também para o ca<strong>na</strong>l II da interface.<br />
No bloco 3, o botão Configura<strong>do</strong>r abre a janela apresentada <strong>na</strong> Figura 4.4. Nesta janela<br />
faz-se o registro <strong>de</strong> alguns parâmetros da bateria, como tensão nomi<strong>na</strong>l, capacida<strong>de</strong>, entre<br />
outros. E <strong>na</strong> peque<strong>na</strong> tabela da janela, é <strong>de</strong>fini<strong>do</strong> o perfil <strong>de</strong> <strong>de</strong>scarga a ser aplica<strong>do</strong> <strong>na</strong><br />
bateria, ou seja, seu tempo em s e seu valor em mA.<br />
Ao fi<strong>na</strong>l <strong>do</strong> experimento, é possível gerar um relatório em formato texto com os parâ-<br />
metros já registra<strong>do</strong>s no software, e também, as informações <strong>de</strong> tensão elétrica, corrente<br />
elétrica e temperatura da bateria, obtidas a partir <strong>do</strong>s sensores <strong>de</strong> corrente, tensão e<br />
temperatura, presentes no hardware <strong>do</strong> testbed.
Capítulo 4. Ambiente <strong>de</strong> Simulação 47<br />
Figura 4.4: Janela <strong>de</strong> configuração <strong>do</strong> teste a ser aplica<strong>do</strong> <strong>na</strong> bateria.<br />
É importante mencio<strong>na</strong>r que durante um experimento no testbed, no caso <strong>de</strong> falhas<br />
<strong>na</strong> comunicação serial entre o hardware e o software, ou quan<strong>do</strong> a curva <strong>de</strong> tensão da<br />
bateria simulada pelo testbed atinge o nível <strong>de</strong> tensão <strong>de</strong> cutoff configura<strong>do</strong> no software,<br />
a plataforma pára <strong>de</strong> operar imediatamente.<br />
4.2 Meto<strong>do</strong>logia A<strong>do</strong>tada nos Testes Experimentais<br />
A meto<strong>do</strong>logia utilizada para os ensaios experimentais foi a <strong>de</strong> a<strong>do</strong>tar um padrão único<br />
em to<strong>do</strong>s os experimentos, objetivan<strong>do</strong> minimizar qualquer alteração no resulta<strong>do</strong> fi<strong>na</strong>l<br />
<strong>do</strong>s testes. Observa-se que foram adquiridas baterias novas, que inicialmente foram sub-<br />
metidas ao mesmo processo <strong>de</strong> carregamento utiliza<strong>do</strong> em [2], ou seja, cada uma <strong>de</strong>las<br />
foi conectada a uma fonte exter<strong>na</strong> <strong>de</strong> carregamento (ver Figura 4.5) para ser submetida<br />
a uma corrente <strong>de</strong> carga lenta e constante correspon<strong>de</strong>n<strong>do</strong> a 20% <strong>de</strong> sua capacida<strong>de</strong> no-<br />
mi<strong>na</strong>l, até atingir um valor máximo <strong>de</strong> carga equivalente a 4,2 V . Portanto, as baterias<br />
foram carregadas com uma corrente inicial <strong>de</strong> carga <strong>de</strong> 190 mA, já que suas capacida<strong>de</strong>s<br />
nomi<strong>na</strong>is são <strong>de</strong> 950 mAh. Após a bateria estar completamente carregada (i.e., tensão<br />
<strong>de</strong> 4,2 V e corrente <strong>de</strong> carga igual a 0), ela foi <strong>de</strong>sconectada da fonte <strong>de</strong> carga e então,<br />
conectada à plataforma para o início <strong>do</strong> processo <strong>de</strong> <strong>de</strong>scarga. Na Figura 4.6 po<strong>de</strong>m ser<br />
observadas duas baterias conectadas à plataforma.<br />
Para a realização <strong>do</strong>s testes foram <strong>de</strong>fini<strong>do</strong>s <strong>de</strong>z perfis <strong>de</strong> <strong>de</strong>scargas variáveis. Con-<br />
si<strong>de</strong>ran<strong>do</strong> o caso <strong>do</strong>s perfis <strong>de</strong> <strong>de</strong>scargas contínuas, não foram realiza<strong>do</strong>s testes, já que<br />
foram utiliza<strong>do</strong>s os da<strong>do</strong>s experimentais coleta<strong>do</strong>s em [2]. Em ambos os casos as correntes<br />
foram medidas em mA. Para os testes experimentais optou-se por utilizar baterias <strong>de</strong> íon<br />
lítio que alimentam dispositivos móveis, como telefones celulares mo<strong>de</strong>lo BL-5F fabrica<strong>do</strong><br />
pela Nokia, que equipam o telefone celular Nokia mo<strong>de</strong>lo N95. Para a configuração <strong>do</strong><br />
software <strong>do</strong> testbed, foram <strong>de</strong>fini<strong>do</strong>s os seguintes parâmetros para o <strong>de</strong>scarregamento da<br />
bateria: (i) tipo <strong>de</strong> bateria: íon lítio; (ii) tensão nomi<strong>na</strong>l: 3,7 V ; (iii) capacida<strong>de</strong> nomi<strong>na</strong>l:
Capítulo 4. Ambiente <strong>de</strong> Simulação 48<br />
Figura 4.5: <strong>Bateria</strong>s conectadas à fonte exter<strong>na</strong>.<br />
Figura 4.6: <strong>Bateria</strong>s conectadas ao testbed.<br />
950 mAh; (iv) corrente e tempo <strong>de</strong> <strong>de</strong>scarga: conforme Tabela 4.1; (v) tensão <strong>de</strong> cutoff :<br />
3,10 V .<br />
Tabela 4.1: Perfis <strong>de</strong> <strong>de</strong>scargas variáveis.<br />
Perfil Degraus <strong>de</strong> corrente <strong>do</strong> perfil (mA) Tempo por <strong>de</strong>grau (s)<br />
P1 300-900-600 600<br />
P2 300-900-600 60<br />
P3 500-950-750 60<br />
P4 100-500-950 60<br />
P5 250-450-650 60<br />
P6 450-0-800 60<br />
P7 800-200-500 60<br />
P8 850-750-500 60<br />
P9 500-200-100 60<br />
P10 470-920 60<br />
Na Tabela 4.1 são apresenta<strong>do</strong>s os perfis <strong>de</strong> <strong>de</strong>scargas variáveis. Verifica-se que <strong>na</strong>
Capítulo 4. Ambiente <strong>de</strong> Simulação 49<br />
primeira colu<strong>na</strong> está a i<strong>de</strong>ntificação <strong>do</strong> perfil; <strong>na</strong> segunda colu<strong>na</strong> há nove perfis <strong>de</strong> <strong>de</strong>scarga<br />
com três <strong>de</strong>graus <strong>de</strong> corrente e ape<strong>na</strong>s um perfil <strong>de</strong> <strong>de</strong>scarga com <strong>do</strong>is <strong>de</strong>graus. Na terceira<br />
colu<strong>na</strong> é indica<strong>do</strong> o tempo que dura cada um <strong>do</strong>s <strong>de</strong>graus <strong>do</strong> perfil. Por exemplo, no caso<br />
da primeira linha da Tabela 4.1, on<strong>de</strong> há um perfil <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 300-900-600 mA para<br />
um tempo por <strong>de</strong>grau <strong>de</strong> 600 s, significa que o <strong>de</strong>grau <strong>de</strong> 300 mA dura 600 s, o <strong>de</strong>grau<br />
<strong>de</strong> 900 mA dura 600 s e o <strong>de</strong>grau <strong>de</strong> 600 mA também dura 600 s.<br />
O tempo <strong>de</strong> vida ou <strong>de</strong> duração <strong>do</strong> processo <strong>de</strong> <strong>de</strong>scarga, é o tempo que a bateria<br />
<strong>de</strong>mora até atingir a tensão <strong>de</strong> cutoff, <strong>de</strong>finida no software em 3,10 V , que é a tensão limite<br />
inferior que possibilita ao dispositivo móvel Nokia N95 permanecer ainda operacio<strong>na</strong>l.<br />
Assim, após cada <strong>de</strong>scarga <strong>de</strong> uma bateria, a mesma era submetida novamente a um<br />
carregamento, sen<strong>do</strong> a<strong>do</strong>tada a meto<strong>do</strong>logia já <strong>de</strong>scrita anteriormente, como forma <strong>de</strong><br />
garantir que a bateria possuísse sua carga completa no início <strong>de</strong> um novo processo <strong>de</strong><br />
<strong>de</strong>scarga.<br />
Deste mo<strong>do</strong>, foram realiza<strong>do</strong>s os testes para os <strong>de</strong>z perfis <strong>de</strong> <strong>de</strong>scargas variáveis<br />
<strong>de</strong>fini<strong>do</strong>s sempre com duas baterias simultaneamente, sen<strong>do</strong> que para cada perfil foram<br />
realiza<strong>do</strong>s <strong>de</strong>z testes, objetivan<strong>do</strong> obter um número significativo <strong>de</strong> amostras, para então<br />
ser calculada a média e <strong>de</strong>svio padrão para cada perfil <strong>de</strong> <strong>de</strong>scarga. Os resulta<strong>do</strong>s <strong>de</strong>stes<br />
testes são apresenta<strong>do</strong>s <strong>na</strong> próxima seção.<br />
4.3 Da<strong>do</strong>s Experimentais Obti<strong>do</strong>s<br />
Nesta seção, são apresenta<strong>do</strong>s os da<strong>do</strong>s experimentais obti<strong>do</strong>s a partir <strong>do</strong> testbed. Na<br />
Tabela 4.2 são apresenta<strong>do</strong>s os da<strong>do</strong>s experimentais obti<strong>do</strong>s consi<strong>de</strong>ran<strong>do</strong> os <strong>de</strong>z perfis <strong>de</strong><br />
<strong>de</strong>scargas variáveis apresenta<strong>do</strong>s <strong>na</strong> Tabela 4.1. Na Tabela 4.3 são apresenta<strong>do</strong>s os da<strong>do</strong>s<br />
experimentais obti<strong>do</strong>s por [2], consi<strong>de</strong>ran<strong>do</strong> <strong>de</strong>z perfis <strong>de</strong> <strong>de</strong>scargas contínuas.<br />
Tabela 4.2: Perfis <strong>de</strong> <strong>de</strong>scargas variáveis, seus tempos, média e <strong>de</strong>svio padrão.<br />
Perfil Tempos <strong>de</strong> vida (s) Média Desvio<br />
padrão<br />
P1 4269 4225,2 4210,8 4498,2 4201,8 4396,8 4221 4488 4435,8 4713 4366,2 2,83<br />
P2 4423,2 4744,8 4048,2 4429,2 4582,2 3883,8 4432,2 4594,8 4431 4765,2 4433,4 4,67<br />
P3 3693 3142,8 3856,2 3684 3676,2 3334,8 3145,2 3343,8 3324 3337,2 3453,6 4,20<br />
P4 5014,2 4672,8 4858,8 5374,8 4843,2 5001 5185,2 5355 4843,2 4993,2 5014,2 3,83<br />
P5 5889 6067,8 5526 6097,2 5883 6247,2 5347,8 6265,8 5710,2 6112,8 5914,8 5,11<br />
P6 6997,8 6834 6285 5938,2 6097,2 6636 6445,2 6837 5896,8 6795 6476,4 6,72<br />
P7 4888,2 5235 5460 5098,2 4738,2 5236,2 5406 5086,8 5067 5059,2 5127,6 3,65<br />
P8 3835,2 3784,8 3697,8 3799,2 3639 3966 3258 3639 3706,2 2905,2 3622,8 5,22<br />
P9 10309,8 10306,2 10480,2 10663,2 10489,8 10669,8 10129,2 10482 11185,2 10483,8 10519,8 4,76<br />
P10 3819 3939 3555 3547,8 2953,2 3466,2 3718,2 3550,2 3808,8 3436,8 3579,6 4,60
Capítulo 4. Ambiente <strong>de</strong> Simulação 50<br />
Tabela 4.3: Perfis <strong>de</strong> <strong>de</strong>scargas contínuas, seus tempos, média e <strong>de</strong>svio padrão [2].<br />
Perfis (mA) Tempos <strong>de</strong> vida (s) Média Desvio<br />
padrão<br />
100 30055,2 28908 29209,2 30381 29779,8 27715,8 28540,8 29158,2 28324,8 28155 29022,6 14,41<br />
150 19662 19743 19510,8 19113 19780,8 19870,2 20263,2 19588,2 18820,8 19047 19540,2 7,22<br />
250 10963,2 11631 11712 11536,2 11952 11587,2 11794,2 11451 11431,8 11332,2 11539,2 4,55<br />
350 8350,8 8263,8 8635,8 8343 7936,8 8446,2 7350 8328 8104,8 7990,2 8175 5,97<br />
450 6424,8 6555 6235,8 6394,8 6160,8 6583,2 6117 6298,8 6241,2 6406,2 6342 2,63<br />
550 5277 5419,2 5056,8 5275,2 5182,2 4906,8 5068,8 5083,8 5331 4881 5148 2,99<br />
650 4309,8 4342,8 4299 3850,2 4419 4533 4264,2 3903 4159,2 4354,8 4243,8 3,61<br />
750 3907,2 3645 3630 3958,2 3007,8 3804 3727,2 3787,2 3688,8 3286,2 3644,4 4,83<br />
850 3079,2 3375 3183,6 3228 2748 3195 3370,2 2754 3016,8 3067,8 3102 3,66<br />
950 2698,8 2794,2 3001,2 2959,2 2509,2 2725,2 2815,8 2908,2 2691 2914,8 2802 2,50<br />
Como exemplo, é apresentada <strong>na</strong> Figura 4.7 uma curva real <strong>de</strong> <strong>de</strong>scarga da bateria <strong>de</strong><br />
íon lítio Nokia BL-5F, obtida a partir <strong>do</strong> testbed, referente ao perfil <strong>de</strong> <strong>de</strong>scarga variável<br />
P1 da Tabela 4.2. Na Figura 4.8, é apresentada uma curva real <strong>de</strong> <strong>de</strong>scarga da mesma<br />
bateria, também obtida <strong>do</strong> testbed, referente ao perfil <strong>de</strong> <strong>de</strong>scarga contínua <strong>de</strong> 150 mA da<br />
Tabela 4.3.<br />
Figura 4.7: Curva real <strong>de</strong> <strong>de</strong>scarga variável para o perfil P1 da bateria <strong>de</strong> íon lítio Nokia BL-5F.<br />
Figura 4.8: Curva real <strong>de</strong> <strong>de</strong>scarga para o perfil <strong>de</strong> <strong>de</strong>scarga contínua <strong>de</strong> 150 mA <strong>na</strong> bateria <strong>de</strong><br />
íon lítio Nokia BL-5F.<br />
A partir das curvas reais <strong>de</strong> <strong>de</strong>scarga obtidas <strong>do</strong>s resulta<strong>do</strong>s experimentais, foi possível
Capítulo 4. Ambiente <strong>de</strong> Simulação 51<br />
realizar simulações computacio<strong>na</strong>is <strong>do</strong> mo<strong>de</strong>lo elétrico Battery, as quais serão <strong>de</strong>scritas e<br />
apresentadas <strong>na</strong> próxima seção.<br />
4.4 Simulações Computacio<strong>na</strong>is<br />
Para a realização das simulações <strong>do</strong> tempo <strong>de</strong> vida das baterias <strong>de</strong> íon lítio foi escolhi<strong>do</strong><br />
o mo<strong>de</strong>lo elétrico Battery, <strong>de</strong>vi<strong>do</strong> a vantagens já mencio<strong>na</strong>das no Capítulo 3, on<strong>de</strong> foi<br />
<strong>de</strong>staca<strong>do</strong> que este mo<strong>de</strong>lo precisa ser calibra<strong>do</strong> uma única vez a partir <strong>de</strong> parâmetros<br />
<strong>de</strong> uma curva real <strong>de</strong> <strong>de</strong>scarga da bateria a ser simulada, bastan<strong>do</strong> posteriormente, a<br />
aplicação <strong>de</strong> diferentes perfis <strong>de</strong> <strong>de</strong>scarga no mo<strong>de</strong>lo, para que ele realize as simulações<br />
<strong>de</strong> <strong>de</strong>scarga da bateria.<br />
Objetivan<strong>do</strong> verificar a partir <strong>de</strong> qual curva real <strong>de</strong> <strong>de</strong>scarga o mo<strong>de</strong>lo elétrico Battery<br />
apresenta os melhores resulta<strong>do</strong>s <strong>de</strong> simulações, o mo<strong>de</strong>lo foi calibra<strong>do</strong> <strong>de</strong>z vezes, sen<strong>do</strong><br />
utiliza<strong>do</strong>s parâmetros <strong>de</strong> <strong>de</strong>z curvas reais <strong>de</strong> <strong>de</strong>scarga diferentes da bateria <strong>de</strong> íon lítio<br />
Nokia BL-5F, a partir das respectivas correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga: 100 mA, 150 mA,<br />
250 mA, 350 mA, 450 mA, 550 mA, 650 mA, 750 mA, 850 mA e 950 mA. O mo<strong>de</strong>lo ainda<br />
foi calibra<strong>do</strong> com os parâmetros obti<strong>do</strong>s a partir da média <strong>do</strong>s parâmetros das <strong>de</strong>z curvas<br />
reais <strong>de</strong> <strong>de</strong>scarga, logo sen<strong>do</strong> utilizada a média das <strong>de</strong>z correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga,<br />
o que resultou em uma corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga média <strong>de</strong> 505 mA. Os parâmetros<br />
obti<strong>do</strong>s a partir das <strong>de</strong>z curvas reais <strong>de</strong> <strong>de</strong>scarga, em conjunto com informações da bateria<br />
e da<strong>do</strong>s <strong>de</strong>termi<strong>na</strong><strong>do</strong>s, usa<strong>do</strong>s para a calibração <strong>do</strong> mo<strong>de</strong>lo em cada corrente nomi<strong>na</strong>l <strong>de</strong><br />
<strong>de</strong>scarga, são apresenta<strong>do</strong>s <strong>na</strong> Tabela 4.4.<br />
Tabela 4.4: Parâmetros utiliza<strong>do</strong>s <strong>na</strong>s simulações <strong>do</strong> mo<strong>de</strong>lo elétrico Battery.<br />
Parâmetros Calibração em uma das correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga<br />
Tipo <strong>de</strong> bateria íon lítio<br />
Tensão cutoff 3,10 V<br />
Capacida<strong>de</strong> típica 0,950 Ah<br />
Esta<strong>do</strong> <strong>de</strong> carga 100%<br />
Capacida<strong>de</strong> máxima 0,9975 Ah<br />
Tensão com carga completa 4,20 V<br />
Corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga Uma das correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga (Ind) A<br />
Resistência inter<strong>na</strong> 0,02945 Ω<br />
Capacida<strong>de</strong> <strong>na</strong> tensão cutoff (Ind·tempo <strong>na</strong> tensão cutoff ) Ah<br />
Tensão exponencial (Tensão observada no fi<strong>na</strong>l <strong>do</strong> comportamento exponencial da curva <strong>de</strong> <strong>de</strong>scarga) V<br />
Capacida<strong>de</strong> exponencial (Ind·tempo <strong>na</strong> tensão exponencial) Ah<br />
Tempo <strong>de</strong> resposta 30 s<br />
Unida<strong>de</strong> no eixo x s<br />
Para as simulações <strong>do</strong> mo<strong>de</strong>lo elétrico Battery, <strong>na</strong> ferramenta computacio<strong>na</strong>l Mat-
Capítulo 4. Ambiente <strong>de</strong> Simulação 52<br />
lab/Simulink, foi utiliza<strong>do</strong> o méto<strong>do</strong> numérico <strong>de</strong> Bogacki Shampine com o tamanho <strong>do</strong><br />
passo fixa<strong>do</strong> em 1 segun<strong>do</strong>. É importante <strong>de</strong>stacar alguns aspectos importantes sobre os<br />
parâmetros da bateria <strong>de</strong> íon lítio Nokia BL-5F, consi<strong>de</strong>ra<strong>do</strong>s durante as simulações. Em<br />
virtu<strong>de</strong> <strong>de</strong> ser uma empresa privada, a Nokia não disponibiliza datasheets técnicos especí-<br />
ficos <strong>de</strong> seus dispositivos ou componentes. Por isto, foram necessários obrigatoriamente,<br />
testes experimentais com esta bateria, para a obtenção <strong>de</strong> curvas reais <strong>de</strong> <strong>de</strong>scarga. Para<br />
a <strong>de</strong>finição <strong>do</strong> parâmetro <strong>de</strong> resistência inter<strong>na</strong>, foi utiliza<strong>do</strong> um valor que correspon<strong>de</strong><br />
a 1% <strong>do</strong> produto da tensão <strong>de</strong> cutoff da bateria pela sua capacida<strong>de</strong> típica, conforme<br />
proposto em [22]. O valor <strong>do</strong> parâmetro capacida<strong>de</strong> máxima foi consi<strong>de</strong>ra<strong>do</strong> como 105%<br />
da capacida<strong>de</strong> típica da bateria, conforme também proposto em [22]. A capacida<strong>de</strong> típica<br />
da bateria é a sua capacida<strong>de</strong> nomi<strong>na</strong>l, que no caso da bateria utilizada é <strong>de</strong> 950 mAh.<br />
O esta<strong>do</strong> <strong>de</strong> carga foi configura<strong>do</strong> como 100%, pois no início <strong>do</strong> processo <strong>de</strong> <strong>de</strong>scarga,<br />
a bateria está totalmente carregada. A tensão <strong>de</strong> cutoff foi <strong>de</strong>finida como 3,10 V , pois<br />
este é o valor <strong>de</strong> limite inferior <strong>de</strong> tensão elétrica que ainda mantém o celular Nokia N95<br />
operacio<strong>na</strong>l. E fi<strong>na</strong>lmente, o parâmetro tempo <strong>de</strong> resposta foi <strong>de</strong>fini<strong>do</strong> como 30 s, por<br />
representar <strong>de</strong> forma satisfatória a dinâmica <strong>de</strong> tensão da bateria [22].<br />
Tanto para as simulações em <strong>de</strong>scargas contínuas, quanto para as simulações em<br />
<strong>de</strong>scargas variáveis, foi necessário utilizar-se <strong>de</strong> um programa <strong>na</strong> forma <strong>de</strong> diagrama <strong>de</strong><br />
blocos, implementa<strong>do</strong> no Simulink, para a geração e posterior tratamento das curvas si-<br />
muladas pelo mo<strong>de</strong>lo. Na Figura 4.9 é apresenta<strong>do</strong> este programa em diagrama <strong>de</strong> blocos.<br />
Figura 4.9: Programa implementa<strong>do</strong> em diagramas <strong>de</strong> blocos no Simulink para as simulações<br />
com o mo<strong>de</strong>lo elétrico Battery.
Capítulo 4. Ambiente <strong>de</strong> Simulação 53<br />
A parte grifada no diagrama, representa os <strong>de</strong>z perfis <strong>de</strong> <strong>de</strong>scargas variáveis, conforme<br />
apresenta<strong>do</strong>s <strong>na</strong> Tabela 4.1. Estes perfis foram cria<strong>do</strong>s a partir <strong>de</strong> um somatório <strong>de</strong> si<strong>na</strong>is<br />
gera<strong>do</strong>s por componentes chama<strong>do</strong>s Pulse generator, presentes no Simulink. Para o caso<br />
das simulações em <strong>de</strong>scargas contínuas, foi utiliza<strong>do</strong> o componente chama<strong>do</strong> Constant,<br />
bastan<strong>do</strong> alterar os seus valores <strong>de</strong> acor<strong>do</strong> com os valores das <strong>de</strong>scargas contínuas. Em<br />
momentos diferentes, os perfis <strong>de</strong> <strong>de</strong>scarga são aplica<strong>do</strong>s ao mo<strong>de</strong>lo elétrico Battery, re-<br />
presenta<strong>do</strong> pelo componente Battery, a partir <strong>do</strong> componente Controlled Current Source.<br />
Conforme já <strong>de</strong>scrito no Capítulo 3, a saída <strong>do</strong> mo<strong>de</strong>lo elétrico Battery via Simulink é<br />
uma matriz <strong>de</strong> or<strong>de</strong>m N × 3 constituída por três si<strong>na</strong>is. Estes si<strong>na</strong>is são: SOC, corrente e<br />
tensão. Durante as simulações com o mo<strong>de</strong>lo, a partir <strong>do</strong>s componentes <strong>de</strong>nomi<strong>na</strong><strong>do</strong>s To<br />
Workspace presentes no diagrama, estes três si<strong>na</strong>is simula<strong>do</strong>s são envia<strong>do</strong>s para a área <strong>de</strong><br />
trabalho <strong>do</strong> Matlab, juntamente com um si<strong>na</strong>l <strong>de</strong> tempo gera<strong>do</strong> pelo compontente Digital<br />
Clock, on<strong>de</strong> to<strong>do</strong>s estes si<strong>na</strong>is po<strong>de</strong>m ser fi<strong>na</strong>lmente manipula<strong>do</strong>s, para em seguida ocorrer<br />
a geração <strong>de</strong> gráficos conten<strong>do</strong> curvas simuladas <strong>de</strong> <strong>de</strong>scarga.<br />
No primeiro momento, as simulações com o mo<strong>de</strong>lo elétrico Battery foram realizadas<br />
consi<strong>de</strong>ran<strong>do</strong> a calibração <strong>do</strong> mo<strong>de</strong>lo a partir <strong>de</strong> uma corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 100<br />
mA, ou seja, consi<strong>de</strong>ran<strong>do</strong> os parâmetros extraí<strong>do</strong>s da curva real <strong>de</strong> <strong>de</strong>scarga origi<strong>na</strong>da<br />
<strong>de</strong>sta corrente nomi<strong>na</strong>l. Após a calibração <strong>do</strong> mo<strong>de</strong>lo, ele foi submeti<strong>do</strong>, em momentos<br />
diferentes, a nove perfis <strong>de</strong> <strong>de</strong>scargas contínuas (i.e., 150 mA, 250 mA, 350 mA, 450 mA,<br />
550 mA, 650 mA, 750 mA, 850 mA e 950 mA). Destaca-se que cada um <strong>de</strong>stes nove<br />
perfis <strong>de</strong> <strong>de</strong>scarga utiliza<strong>do</strong>s <strong>na</strong>s simulações <strong>do</strong> mo<strong>de</strong>lo, calibra<strong>do</strong> em 100 mA, vai ser em<br />
outro momento, a corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>do</strong> mo<strong>de</strong>lo utilizada para calibração <strong>do</strong><br />
mesmo.<br />
Portanto, este procedimento foi repeti<strong>do</strong> para o mo<strong>de</strong>lo calibra<strong>do</strong> <strong>na</strong>s correntes no-<br />
mi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 150 mA, 250 mA, 350 mA, 450 mA, 550 mA, 650 mA, 750 mA,<br />
850 mA e 950 mA. No caso <strong>de</strong> simulações em <strong>de</strong>scargas contínuas, para cada calibração<br />
<strong>do</strong> mo<strong>de</strong>lo, foram realizadas nove simulações, pois observa-se que não é realizada a si-<br />
mulação com a corrente nomi<strong>na</strong>l a partir da qual o mo<strong>de</strong>lo é calibra<strong>do</strong>. Observa-se que<br />
simulações realizadas a partir da aplicação ao mo<strong>de</strong>lo, da corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga<br />
<strong>do</strong> mo<strong>de</strong>lo, e consequentemente a partir <strong>do</strong>s parâmetros obti<strong>do</strong>s da curva real <strong>de</strong> <strong>de</strong>scarga<br />
gerada por esta corrente, apresentam tempos <strong>de</strong> vida simula<strong>do</strong>s com um erro bastante<br />
próximo <strong>de</strong> 0%. Similarmente, para o caso da corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga média, ou<br />
seja, calibração <strong>do</strong> mo<strong>de</strong>lo em 505 mA, não foram realizadas as simulações com os perfis<br />
<strong>de</strong> <strong>de</strong>scargas contínuas em questão, justamente porque estes perfis foram utiliza<strong>do</strong>s para<br />
a <strong>de</strong>termi<strong>na</strong>ção <strong>do</strong>s parâmetros <strong>de</strong> calibração <strong>do</strong> mo<strong>de</strong>lo em 505 mA.<br />
No caso das simulações com <strong>de</strong>scargas variáveis, para cada calibração <strong>do</strong> mo<strong>de</strong>lo,<br />
incluin<strong>do</strong> a calibração média (i.e., calibração em 505 mA), foram aplica<strong>do</strong>s <strong>de</strong>z perfis <strong>de</strong>
Capítulo 4. Ambiente <strong>de</strong> Simulação 54<br />
<strong>de</strong>scarga variáveis ao mo<strong>de</strong>lo, totalizan<strong>do</strong> <strong>de</strong>z simulações para cada calibração <strong>do</strong> mo<strong>de</strong>lo.<br />
Portanto, para o caso <strong>de</strong> <strong>de</strong>scargas contínuas houve um total <strong>de</strong> <strong>de</strong>z calibrações, enquanto<br />
que para o caso <strong>de</strong> <strong>de</strong>scargas variáves, houve um total <strong>de</strong> onze calibrações <strong>do</strong> mo<strong>de</strong>lo.<br />
4.5 Análise <strong>do</strong>s Resulta<strong>do</strong>s <strong>de</strong> Simulações<br />
Nesta seção são apresenta<strong>do</strong>s os resulta<strong>do</strong>s das simulações <strong>do</strong> mo<strong>de</strong>lo elétrico Battery<br />
consi<strong>de</strong>ran<strong>do</strong> uma bateria <strong>de</strong> íon lítio Nokia mo<strong>de</strong>lo BL-5F. Em seguida é realizada uma<br />
análise comparativa, objetivan<strong>do</strong> verificar em qual calibração o mo<strong>de</strong>lo elétrico apresenta<br />
melhores resulta<strong>do</strong>s para a predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> uma bateria <strong>de</strong> íon lítio.<br />
Na Tabela 4.5 são apresenta<strong>do</strong>s os tempos <strong>de</strong> vida médios da bateria <strong>de</strong> íon lítio obti<strong>do</strong>s<br />
<strong>de</strong> forma experimental (Tve) por [2], e os tempos <strong>de</strong> vida simula<strong>do</strong>s (Tvs), obti<strong>do</strong>s a partir<br />
<strong>do</strong> mo<strong>de</strong>lo elétrico Battery, consi<strong>de</strong>ran<strong>do</strong> o mo<strong>de</strong>lo calibra<strong>do</strong> <strong>na</strong>s correntes nomi<strong>na</strong>is <strong>de</strong><br />
<strong>de</strong>scarga <strong>de</strong> 100 mA, 150 mA, 250 mA, 350 mA, 450 mA, 550 mA, 650 mA, 750 mA,<br />
850 mA e 950 mA.<br />
Tabela 4.5: Tempos <strong>de</strong> vida experimentais (Tve) e tempos <strong>de</strong> vida simula<strong>do</strong>s (Tvs) pelo mo<strong>de</strong>lo<br />
elétrico Battery com <strong>de</strong>scargas contínuas.<br />
Cal. Perfis 100 150 250 350 450 550 650 750 850 950 Erro médio<br />
(mA) (mA) (%)<br />
- Tve (s) 29022,6 19540,2 11539,2 8175 6342 5148 4243,8 3644,4 3102 2802 -<br />
100 Tvs (s) - 19055 11170 7790 5911 4715 3887 3279 2814 2447 -<br />
- Erro (%) - 2,48 3,20 4,71 6,80 8,41 8,41 10,03 9,28 12,67 7,33<br />
150 Tvs (s) 29565 - 11467 8019 6103 4882 4037 3417 2942 2567 -<br />
- Erro (%) 1,87 - 0,63 1,91 3,77 5,17 4,87 6,24 5,16 8,39 4,22<br />
250 Tvs (s) 29715 19615 - 8068 6143 4916 4067 3443 2966 2589 -<br />
- Erro (%) 2,39 0,38 - 1,31 3,14 4,51 4,17 5,53 4,38 7,60 3,71<br />
350 Tvs (s) 29810 19681 11577 - 6172 4942 4090 3465 2986 2608 -<br />
- Erro (%) 2,71 0,72 0,33 - 2,68 4,00 3,62 4,92 3,74 6,92 3,29<br />
450 Tvs (s) 30237 19980 11772 8253 - 5052 4190 3556 3072 2689 -<br />
- Erro (%) 4,18 2,25 2,02 0,95 - 1,86 1,27 2,43 0,97 4,03 2,22<br />
550 Tvs (s) 30732 20324 11997 8428 6444 - 4305 3663 3171 2783 -<br />
- Erro (%) 5,89 4,01 3,97 3,09 1,61 - 1,44 0,51 2,22 0,68 2,60<br />
650 Tvs (s) 30566 20207 11919 8366 6391 5133 - 3622 3132 2745 -<br />
- Erro (%) 5,32 3,41 3,29 2,34 0,77 0,29 - 0,61 0,97 2,03 2,12<br />
750 Tvs (s) 30697 20297 11976 8408 6425 5161 4286 - 3151 2762 -<br />
- Erro (%) 5,77 3,87 3,79 2,85 1,31 0,25 0,99 - 1,58 1,43 2,43<br />
850 Tvs (s) 30447 20116 11848 8302 6329 5072 4200 3560 - 2679 -<br />
- Erro (%) 4,91 2,95 2,68 1,55 0,20 1,48 1,03 2,32 - 4,39 2,39<br />
950 Tvs (s) 30850 20404 12046 8463 6471 5202 4323 3677 3183 - -<br />
- Erro (%) 6,30 4,42 4,39 3,52 2,03 1,05 1,87 0,89 2,61 - 3,01<br />
Cal. = Calibração, Tve = Tempo <strong>de</strong> vida experimental, Tvs = Tempo <strong>de</strong> vida simula<strong>do</strong>
Capítulo 4. Ambiente <strong>de</strong> Simulação 55<br />
Na Tabela 4.6 são apresenta<strong>do</strong>s os tempos <strong>de</strong> vida médios da bateria <strong>de</strong> íon lítio<br />
obti<strong>do</strong>s <strong>de</strong> forma experimental para perfis <strong>de</strong> <strong>de</strong>scargas variáveis e os tempos <strong>de</strong> vida<br />
simula<strong>do</strong>s, obti<strong>do</strong>s a partir <strong>do</strong> mo<strong>de</strong>lo elétrico Battery, consi<strong>de</strong>ran<strong>do</strong> o mo<strong>de</strong>lo calibra<strong>do</strong><br />
<strong>na</strong>s correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 100 mA, 150 mA, 250 mA, 350 mA, 450 mA, 550<br />
mA, 650 mA, 750 mA, 850 mA, 950 mA e 505 mA.<br />
Tabela 4.6: Tempos <strong>de</strong> vida experimentais (Tve) e tempos <strong>de</strong> vida simula<strong>do</strong>s (Tvs) pelo mo<strong>de</strong>lo<br />
elétrico Battery com <strong>de</strong>scargas variáveis.<br />
Cal. Perfis P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 Erro médio<br />
(mA) (mA) (%)<br />
- Tve (s) 4366,2 4433,4 3453,6 5014,2 5914,8 6476,4 5127,48 3622,8 10520 3579,6 -<br />
100 Tvs (s) 4234 4040 3174 4499 5713 5910 4882 3444 9924 3441 -<br />
- Erro (%) 3,03 8,87 8,10 10,27 3,41 8,75 4,79 4,94 5,67 3,87 6,17<br />
150 Tvs (s) 4344 4220 3340 4821 5897 6104 5068 3618 10137 3569 -<br />
- Erro (%) 0,51 4,81 3,29 3,85 0,30 5,75 1,16 0,13 3,64 0,30 2,37<br />
250 Tvs (s) 4366 4224 3356 4826 5907 6264 5083 3624 10287 3579 -<br />
- Erro (%) 0,00 4,72 2,83 3,75 0,13 3,28 0,87 0,03 2,21 0,02 1,79<br />
350 Tvs (s) 4386 4229 3493 4834 5920 6267 5232 3632 10296 3594 -<br />
- Erro (%) 0,45 4,61 1,14 3,59 0,09 3,23 2,04 0,25 2,13 0,40 1,79<br />
450 Tvs (s) 4469 4400 3510 5006 6085 6446 5260 3794 10488 3707 -<br />
- Erro (%) 2,35 0,75 1,63 0,16 2,88 0,47 2,58 4,73 0,30 3,56 1,94<br />
550 Tvs (s) 4566 4576 3677 5184 6261 6627 5433 3839 10824 3831 -<br />
- Erro (%) 4,58 3,22 6,47 3,39 5,85 2,33 5,96 5,97 2,89 7,02 4,77<br />
650 Tvs (s) 4527 4417 3672 5176 6251 6470 5418 3810 10662 3808 -<br />
- Erro (%) 3,68 0,37 6,32 3,23 5,68 0,10 5,67 5,17 1,35 6,39 3,79<br />
750 Tvs (s) 4545 4428 3674 5179 6256 6623 5424 3822 10821 3816 -<br />
- Erro (%) 4,10 0,12 6,38 3,29 5,77 2,26 5,78 5,50 2,86 6,60 4,27<br />
850 Tvs (s) 4462 4399 3506 5003 6091 6446 5260 3794 10642 3701 -<br />
- Erro (%) 2,19 0,78 1,52 0,22 2,98 0,47 2,58 4,73 1,16 3,39 2,00<br />
950 Tvs (s) 4576 4577 3679 5188 6269 6630 5443 3869 10831 3816 -<br />
- Erro (%) 4,81 3,24 6,53 3,47 5,99 2,37 6,15 6,80 2,96 6,60 4,89<br />
505 Tvs (s) 4459 4398 3506 5003 6082 6444 5254 3793 10481 3700 -<br />
- Erro (%) 2,13 0,80 1,52 0,22 2,83 0,50 2,47 4,70 0,37 3,36 1,89<br />
Cal. = Calibração, Tve = Tempo <strong>de</strong> vida experimental, Tvs = Tempo <strong>de</strong> vida simula<strong>do</strong><br />
A aproximação entre os tempos <strong>de</strong> vida experimentais e simula<strong>do</strong>s apresenta<strong>do</strong>s <strong>na</strong>s<br />
Tabelas 4.5 e 4.6 indica que o mo<strong>de</strong>lo elétrico Battery prediz o tempo <strong>de</strong> vida da bateria<br />
com uma acurácia satisfatória. Observa-se pelos resulta<strong>do</strong>s, que para todas as calibrações<br />
realizadas no caso <strong>de</strong> <strong>de</strong>scargas contínuas, somente <strong>na</strong> calibração <strong>de</strong> 100 mA, o mo<strong>de</strong>lo<br />
apresentou um erro médio que ultrapassou 5%, atingin<strong>do</strong> um erro <strong>de</strong> 7,33%, da mesma<br />
forma, no caso das <strong>de</strong>scargas variáveis, atingin<strong>do</strong> um erro <strong>de</strong> 6,17%, o que indica que<br />
tanto para as <strong>de</strong>scargas contínuas, quanto para as <strong>de</strong>scargas variáveis, o maior erro médio<br />
apresenta<strong>do</strong> pelo mo<strong>de</strong>lo ocorreu nesta calibração.
Capítulo 4. Ambiente <strong>de</strong> Simulação 56<br />
Conforme po<strong>de</strong> ser verifica<strong>do</strong> <strong>na</strong> Tabela 4.3, foi justamente o perfil <strong>de</strong> <strong>de</strong>scarga con-<br />
tínua <strong>de</strong> 100 mA que apresentou o maior <strong>de</strong>svio padrão. Isto ocorre <strong>de</strong>vi<strong>do</strong> ao fato <strong>de</strong><br />
que a corrente <strong>de</strong> <strong>de</strong>scarga é baixa, pois nesta situação, os efeitos não-lineares, em espe-<br />
cial o efeito <strong>de</strong> recuperação, se fazem mais presentes, origi<strong>na</strong>n<strong>do</strong> maiores variações nos<br />
tempos <strong>de</strong> <strong>de</strong>scarga e consequentemente, <strong>na</strong>s curvas reais <strong>de</strong> <strong>de</strong>scarga. Desta forma, o<br />
processo <strong>de</strong> extração <strong>de</strong> parâmetros <strong>na</strong> curva <strong>de</strong> <strong>de</strong>scarga que é necessário para o mo<strong>de</strong>lo<br />
elétrico Battery, tor<strong>na</strong>-se menos preciso, o que vai influenciar diretamente <strong>na</strong> acurácia<br />
<strong>do</strong> mo<strong>de</strong>lo. No entanto, nota-se que para as <strong>de</strong>mais calibrações, consi<strong>de</strong>ran<strong>do</strong> tanto as<br />
<strong>de</strong>scargas contínuas quanto as <strong>de</strong>scargas variáveis, o erro médio apresenta<strong>do</strong> pelo mo<strong>de</strong>lo<br />
não ultrapassou os 5%, valor este espera<strong>do</strong> para o caso <strong>de</strong> mo<strong>de</strong>los elétricos [1].<br />
Consi<strong>de</strong>ran<strong>do</strong> que a bateria <strong>de</strong> íon lítio possui uma capacida<strong>de</strong> típica <strong>de</strong> 950 mAh, para<br />
as simulações, foram feitas as calibrações <strong>do</strong> mo<strong>de</strong>lo com parâmetros basea<strong>do</strong>s em uma<br />
corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga baixa (100 mA), até uma corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga alta<br />
(950 mA). Para <strong>de</strong>scargas contínuas, o mo<strong>de</strong>lo apresentou o menor erro médio quan<strong>do</strong><br />
sua calibração foi realizada com parâmetros da curva real <strong>de</strong> <strong>de</strong>scarga baseada <strong>na</strong> corrente<br />
nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 650 mA, igual a 2,12%. Já para <strong>de</strong>scargas variáveis, o mo<strong>de</strong>lo<br />
apresentou o mesmo e o menor erro médio em duas calibrações, ou seja, <strong>na</strong>s calibrações <strong>de</strong><br />
250 mA e 350 mA. Este erro foi <strong>de</strong> 1,79%. Porém, <strong>de</strong>staca-se que <strong>na</strong> corrente nomi<strong>na</strong>l <strong>de</strong><br />
<strong>de</strong>scarga intermediária <strong>de</strong> 450 mA, tanto em <strong>de</strong>scargas contínuas, quanto em <strong>de</strong>scargas<br />
variáveis, o mo<strong>de</strong>lo apresentou resulta<strong>do</strong>s satisfatórios, próximos <strong>do</strong> melhor resulta<strong>do</strong><br />
apresenta<strong>do</strong> pelo mo<strong>de</strong>lo.<br />
A<strong>na</strong>lisan<strong>do</strong> os resulta<strong>do</strong>s da Tabela 4.5, percebe-se que os erros no tempo <strong>de</strong> vida<br />
apresenta<strong>do</strong>s pelo mo<strong>de</strong>lo, variam consi<strong>de</strong>ravelmente à medida em que os valores das<br />
<strong>de</strong>scargas contínuas se afastam <strong>do</strong> valor da corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga utilizada para<br />
calibrar o mo<strong>de</strong>lo. Nota-se também, que os erros apresenta<strong>do</strong>s em cada uma das duas<br />
tabelas possuem uma diferença consi<strong>de</strong>rável entre o mínimo e o máximo <strong>de</strong> erro no tempo<br />
<strong>de</strong> vida, consi<strong>de</strong>ran<strong>do</strong> a mesma calibração. Isto ocorre <strong>de</strong>vi<strong>do</strong> ao mo<strong>de</strong>lo elétrico Battery<br />
não capturar um efeito não-linear importante, que ocorre <strong>na</strong> bateria durante seu processo<br />
<strong>de</strong> <strong>de</strong>scarga, chama<strong>do</strong> <strong>de</strong> efeito da taxa <strong>de</strong> capacida<strong>de</strong>.<br />
É importante <strong>de</strong>stacar que todas as simulações executadas pelo mo<strong>de</strong>lo elétrico Bat-<br />
tery po<strong>de</strong>m ser baseadas em uma única curva real <strong>de</strong> <strong>de</strong>scarga da bateria simulada. Os<br />
parâmetros obti<strong>do</strong>s a partir <strong>de</strong>sta curva são utiliza<strong>do</strong>s para inicializar o mo<strong>de</strong>lo elétrico, e<br />
após esta inicialização, basta aplicar diferentes perfis <strong>de</strong> <strong>de</strong>scarga ao mo<strong>de</strong>lo, para que ele<br />
realize as simulações <strong>de</strong> <strong>de</strong>scarga da bateria. Porém, neste trabalho, as simulações foram<br />
baseadas em <strong>de</strong>z curvas reais <strong>de</strong> <strong>de</strong>scarga, objetivan<strong>do</strong>-se verificar a partir <strong>de</strong> qual <strong>de</strong>las,<br />
o mo<strong>de</strong>lo apresentou os melhores resulta<strong>do</strong>s simula<strong>do</strong>s.<br />
Como exemplo, <strong>na</strong> Figura 4.10, é apresentada a comparação entre os resulta<strong>do</strong>s simu-
Capítulo 4. Ambiente <strong>de</strong> Simulação 57<br />
la<strong>do</strong>s pelo mo<strong>de</strong>lo elétrico Battery e os da<strong>do</strong>s experimentais, consi<strong>de</strong>ran<strong>do</strong> as <strong>de</strong>scargas<br />
contínuas <strong>de</strong> 100 mA, 350 mA, 650 mA e 950 mA aplicadas <strong>na</strong> bateria <strong>de</strong> íon lítio Nokia<br />
BL-5F. Neste exemplo, as simulações foram realizadas com o mo<strong>de</strong>lo calibra<strong>do</strong> <strong>na</strong> cor-<br />
rente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 150 mA. Nas quatro comparações po<strong>de</strong>m ser observa<strong>do</strong>s<br />
os pontos A e B. O ponto A está localiza<strong>do</strong> no início <strong>do</strong> processo <strong>de</strong> <strong>de</strong>scarga e o ponto<br />
B, no fi<strong>na</strong>l, quan<strong>do</strong> é atingida a tensão <strong>de</strong> cutoff <strong>de</strong> 3,10 V . No ponto A, <strong>na</strong>s quatro<br />
comparações, nota-se uma diferença <strong>de</strong> tensão entre as curvas (i.e., experimental e simu-<br />
lada) nos instantes iniciais, que ocorre em <strong>de</strong>corrência <strong>do</strong> comportamento real da tensão<br />
elétrica nos termi<strong>na</strong>is da bateria nos ensaios com a plataforma, pois em um processo real<br />
<strong>de</strong> <strong>de</strong>scarga, a tensão elétrica diminui bruscamente nos instantes iniciais e <strong>de</strong>pois ten<strong>de</strong> a<br />
estabilizar. Ao mesmo tempo, verifica-se também nos instantes iniciais das comparações,<br />
que esta diferença entre as curvas vai aumentan<strong>do</strong>, conforme vai aumentan<strong>do</strong> a corrente<br />
<strong>de</strong> <strong>de</strong>scarga, o que evi<strong>de</strong>ncia a presença <strong>do</strong> efeito da taxa <strong>de</strong> capacida<strong>de</strong> no processo real<br />
<strong>de</strong> <strong>de</strong>scarga, ou seja, quanto maior a corrente <strong>de</strong> <strong>de</strong>scarga em uma bateria, menor será sua<br />
capacida<strong>de</strong> efetiva. Nas quatro comparações, o ponto B representa a diferença <strong>de</strong> tempo<br />
para atingir a tensão <strong>de</strong> cutoff, ou seja, a diferença no tempo <strong>de</strong> vida entre as curvas reais<br />
e simuladas. Nota-se que quanto mais o valor da corrente <strong>de</strong> <strong>de</strong>scarga aplicada ao mo<strong>de</strong>lo<br />
vai se afastan<strong>do</strong> <strong>do</strong> valor da corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>do</strong> mo<strong>de</strong>lo (i.e., 150 mA neste<br />
caso), ou seja, corrente <strong>na</strong> qual o mo<strong>de</strong>lo está calibra<strong>do</strong>, maior vai fican<strong>do</strong> o erro no tempo<br />
<strong>de</strong> vida simula<strong>do</strong>, conforme já evi<strong>de</strong>ncia<strong>do</strong> <strong>na</strong> análise da Tabela 4.5.<br />
No exemplo da Figura 4.11, é apresentada uma comparação entre os resulta<strong>do</strong>s simu-<br />
la<strong>do</strong>s pelo mo<strong>de</strong>lo elétrico Battery e os da<strong>do</strong>s experimentais, consi<strong>de</strong>ran<strong>do</strong> os perfis <strong>de</strong><br />
<strong>de</strong>scargas variáveis P3, P5 e P7 aplica<strong>do</strong>s à bateria <strong>de</strong> íon lítio Nokia BL-5F. Neste exem-<br />
plo, as simulações também foram realizadas com o mo<strong>de</strong>lo elétrico Battery calibra<strong>do</strong> em<br />
150 mA. Neste caso <strong>de</strong> <strong>de</strong>scargas variáveis, também ocorre uma diferença nos instantes<br />
iniciais, entre a curva real e a curva simulada, e esta diferença vai <strong>de</strong>pen<strong>de</strong>r <strong>do</strong> valor <strong>do</strong><br />
<strong>de</strong>grau <strong>de</strong> corrente inicial, ou seja, quanto maior for o valor <strong>de</strong>ste <strong>de</strong>grau, maior será a<br />
diferença. Na questão <strong>de</strong> diferença no tempo <strong>de</strong> vida, entre a curva real e a curva simu-<br />
lada, no fi<strong>na</strong>l <strong>do</strong> processo <strong>de</strong> <strong>de</strong>scarga, não foi observada nenhuma relação com a corrente<br />
nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>do</strong> mo<strong>de</strong>lo.
Capítulo 4. Ambiente <strong>de</strong> Simulação 58<br />
Figura 4.10: Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais <strong>na</strong>s <strong>de</strong>scargas con-<br />
tínuas <strong>de</strong> (a) 100 mA, (b) 350 mA, (c) 650 mA e (d) 950 mA, <strong>na</strong> bateria <strong>de</strong> íon lítio BL-5F.
Capítulo 4. Ambiente <strong>de</strong> Simulação 59<br />
Figura 4.11: Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais <strong>do</strong>s perfis <strong>de</strong> <strong>de</strong>scar-<br />
gas variáveis (a) P3, (b) P5 e (c) P7, <strong>na</strong> bateria <strong>de</strong> íon lítio BL-5F.
Capítulo 5<br />
Análise Comparativa <strong>do</strong>s Mo<strong>de</strong>los<br />
Neste capítulo, é apresenta<strong>do</strong> um estu<strong>do</strong> comparativo, a partir <strong>de</strong> da<strong>do</strong>s experimentais e<br />
resulta<strong>do</strong>s simula<strong>do</strong>s, entre <strong>do</strong>is mo<strong>de</strong>los elétricos <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> baterias: o mo<strong>de</strong>lo Bat-<br />
tery, e o mo<strong>de</strong>lo para Predizer Runtime e Características V-I <strong>de</strong> uma <strong>Bateria</strong> consi<strong>de</strong>ra<strong>do</strong><br />
<strong>na</strong> literatura o mo<strong>de</strong>lo elétrico mais acura<strong>do</strong> [1].<br />
Conforme já aborda<strong>do</strong> no Capítulo 3 o mo<strong>de</strong>lo elétrico Battery possui um processo<br />
prático <strong>de</strong> extração <strong>de</strong> parâmetros. Observa-se que uma vantagem <strong>de</strong>ste mo<strong>de</strong>lo resi<strong>de</strong> <strong>na</strong><br />
presença, no datasheet da bateria a ser simulada, <strong>de</strong> uma curva característica <strong>de</strong> <strong>de</strong>scarga,<br />
quan<strong>do</strong> este fato ocorre uma fase muito importante <strong>do</strong> processo <strong>de</strong> extração <strong>de</strong> parâme-<br />
tros po<strong>de</strong> ser ignorada, que é a fase <strong>do</strong>s testes experimentais. Por outro la<strong>do</strong>, conforme já<br />
<strong>de</strong>scrito no Capítulo 2, verifica-se que o mo<strong>de</strong>lo para Predizer Runtime e Características<br />
V-I, possui um processo <strong>de</strong> extração <strong>de</strong> parâmetros mais <strong>de</strong>mora<strong>do</strong>, <strong>de</strong>mandan<strong>do</strong> obriga-<br />
toriamente, a fase <strong>de</strong> testes experimentais com a bateria a ser simulada. Neste contexto,<br />
é apresentada, neste capítulo, uma comparação entre estes <strong>do</strong>is mo<strong>de</strong>los elétricos encon-<br />
tra<strong>do</strong>s <strong>na</strong> literatura, consi<strong>de</strong>ran<strong>do</strong> uma bateria <strong>de</strong> íon lítio polímero e perfis <strong>de</strong> <strong>de</strong>scarga<br />
utiliza<strong>do</strong>s em [1].<br />
Este capítulo está organiza<strong>do</strong> como segue. Na Seção 5.1 é apresentada a aplicação<br />
<strong>do</strong> mo<strong>de</strong>lo elétrico para Predizer Runtime e Características V-I consi<strong>de</strong>ran<strong>do</strong> <strong>do</strong>is casos<br />
<strong>de</strong> simulação, o primeiro com corrente contínua e o segun<strong>do</strong> com corrente variável. Na<br />
Seção 5.2 é apresentada a aplicação <strong>do</strong> mo<strong>de</strong>lo Battery também consi<strong>de</strong>ran<strong>do</strong> <strong>do</strong>is casos<br />
<strong>de</strong> simulação. Na Seção 5.3 é realizada uma comparação entre os mo<strong>de</strong>los.<br />
5.1 <strong>Aplicação</strong> <strong>do</strong> Mo<strong>de</strong>lo Elétrico para Predizer Run-<br />
time e Características V-I <strong>de</strong> uma bateria<br />
O mo<strong>de</strong>lo elétrico para Predizer Runtime e Características V-I é avalia<strong>do</strong>, após a <strong>de</strong>termi-<br />
<strong>na</strong>ção <strong>de</strong> seus parâmetros, a partir <strong>de</strong> uma comparação entre seus resulta<strong>do</strong>s simula<strong>do</strong>s,<br />
60
Capítulo 5. Análise Comparativa <strong>do</strong>s Mo<strong>de</strong>los 61<br />
e os da<strong>do</strong>s experimentais <strong>de</strong> uma bateria <strong>de</strong> íon lítio polímero PL-383562 com capacida<strong>de</strong><br />
<strong>de</strong> 850 mAh, obti<strong>do</strong>s a partir <strong>de</strong> um testbed. Para esta avaliação, foram aplica<strong>do</strong>s ao<br />
mo<strong>de</strong>lo <strong>do</strong>is perfis <strong>de</strong> <strong>de</strong>scarga, o primeiro foi um perfil <strong>de</strong> <strong>de</strong>scarga contínua <strong>de</strong> 80 mA,<br />
o segun<strong>do</strong> foi um perfil <strong>de</strong> <strong>de</strong>scarga variável com quatro <strong>de</strong>graus <strong>de</strong> corrente <strong>de</strong> 0 mA,<br />
400 mA, 160 mA e 640 mA para os intervalos <strong>de</strong> tempo <strong>de</strong> 615 s, 738 s, 1785 s e 431 s,<br />
respectivamente.<br />
Na Figura 5.1 são apresentadas as curvas <strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s e <strong>do</strong>s da<strong>do</strong>s expe-<br />
rimentais consi<strong>de</strong>ran<strong>do</strong> a <strong>de</strong>scarga contínua. O erro obti<strong>do</strong> <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida<br />
pelo mo<strong>de</strong>lo, para uma bateria <strong>de</strong> íon lítio polímero, foi <strong>de</strong> 0,395% [1]. Na Figura 5.2 são<br />
apresentadas as curvas <strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s e <strong>do</strong>s da<strong>do</strong>s experimentais consi<strong>de</strong>ran<strong>do</strong><br />
a <strong>de</strong>scarga variável, bem como os <strong>de</strong>graus <strong>de</strong> correntes consi<strong>de</strong>ra<strong>do</strong>s para simulação. O<br />
erro obti<strong>do</strong> <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida, pelo mo<strong>de</strong>lo, para uma bateria <strong>de</strong> íon lítio<br />
polímero foi <strong>de</strong> 0,338% [1].<br />
Figura 5.1: Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais para a <strong>de</strong>scarga<br />
contínua <strong>de</strong> 80 mA, <strong>na</strong> bateria <strong>de</strong> íon lítio polímero PL-383562 [1].<br />
Figura 5.2: Comparação entre resulta<strong>do</strong>s simula<strong>do</strong>s e da<strong>do</strong>s experimentais para a <strong>de</strong>scarga <strong>de</strong><br />
quatro <strong>de</strong>graus <strong>de</strong> corrente, <strong>na</strong> bateria <strong>de</strong> íon lítio polímero PL-383562 [1].
Capítulo 5. Análise Comparativa <strong>do</strong>s Mo<strong>de</strong>los 62<br />
A aproximação observada <strong>na</strong>s Figuras 5.1 e 5.2, entre os resulta<strong>do</strong>s das simulações e<br />
os da<strong>do</strong>s experimentais, indica que este mo<strong>de</strong>lo elétrico <strong>de</strong> bateria proposto em [1] prediz<br />
o tempo <strong>de</strong> vida da bateria com uma acurácia satisfatória.<br />
5.2 <strong>Aplicação</strong> <strong>do</strong> Mo<strong>de</strong>lo Elétrico Battery<br />
Consi<strong>de</strong>ran<strong>do</strong> que o objetivo <strong>de</strong>ste capítulo é apresentar um estu<strong>do</strong> comparativo entre<br />
o mo<strong>de</strong>lo elétrico Battery e o mo<strong>de</strong>lo elétrico para Predizer Runtime e Características<br />
V-I <strong>de</strong> uma bateria. Neste contexto, o mo<strong>de</strong>lo Battery foi emprega<strong>do</strong> consi<strong>de</strong>ran<strong>do</strong> os<br />
mesmos cenários utiliza<strong>do</strong>s para as simulações <strong>do</strong> mo<strong>de</strong>lo elétrico para Predizer Runtime<br />
e Características V-I, ou seja, uma <strong>de</strong>scarga contínua <strong>de</strong> 80 mA e uma <strong>de</strong>scarga variável<br />
com os <strong>de</strong>graus <strong>de</strong> corrente <strong>de</strong> 0 mA, 400 mA, 160 mA e 640 mA, uma tensão com carga<br />
completa <strong>de</strong> 4,1 V e uma tensão <strong>de</strong> cutoff <strong>de</strong> 3 V .<br />
Para a realização das simulações, a partir <strong>do</strong> emprego <strong>do</strong> mo<strong>de</strong>lo Battery, este precisa<br />
ser calibra<strong>do</strong> com parâmetros extraí<strong>do</strong>s <strong>de</strong> uma curva real <strong>de</strong> <strong>de</strong>scarga da bateria simulada.<br />
Na Figura 5.3 são apresentadas algumas curvas reais <strong>de</strong> <strong>de</strong>scarga da bateria <strong>de</strong> íon lítio<br />
polímero PL-383562. Para a calibração foram extraí<strong>do</strong>s os parâmetros da curva real <strong>de</strong><br />
<strong>de</strong>scarga <strong>de</strong>sta bateria origi<strong>na</strong>da a partir da corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> 360 mA, em<br />
conjunto com da<strong>do</strong>s <strong>do</strong> seu datasheet e com parâmetros <strong>de</strong>termi<strong>na</strong><strong>do</strong>s conforme proposto<br />
em [22]. Esta corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga foi escolhida, porque representa <strong>de</strong> forma<br />
satisfatória o parâmetro, tensão com carga completa, que neste caso <strong>de</strong>ve ser <strong>de</strong> 4,1 V ,<br />
Figura 5.3: Curvas reais <strong>de</strong> <strong>de</strong>scarga da bateria <strong>de</strong> íon lítio polímero PL-383562 para diferentes<br />
correntes contínuas <strong>de</strong> <strong>de</strong>scarga [25].
Capítulo 5. Análise Comparativa <strong>do</strong>s Mo<strong>de</strong>los 63<br />
e também por ser uma corrente intermediária. Na comparação realizada no Capítulo<br />
4, entre os resulta<strong>do</strong>s simula<strong>do</strong>s <strong>do</strong> mo<strong>de</strong>lo elétrico Battery e resulta<strong>do</strong>s experimentais,<br />
verificou-se que o mo<strong>de</strong>lo calibra<strong>do</strong> em uma corrente nomi<strong>na</strong>l intermidiária, apresenta<br />
resulta<strong>do</strong>s satisfatórios. Então, os seguintes parâmetros ficaram <strong>de</strong>fini<strong>do</strong>s para a simulação<br />
<strong>de</strong> <strong>de</strong>scarga da bateria consi<strong>de</strong>rada, conforme apresenta<strong>do</strong>s <strong>na</strong> Tabela 5.1.<br />
Tabela 5.1: Parâmetros da bateria <strong>de</strong> íon lítio polímero PL-383562 utiliza<strong>do</strong>s <strong>na</strong>s simulações <strong>do</strong><br />
mo<strong>de</strong>lo Battery.<br />
Calibração<br />
Parâmetros <strong>de</strong> 360 mA<br />
Tipo <strong>de</strong> bateria íon lítio<br />
Tensão <strong>de</strong> cutoff 3,00 V<br />
Capacida<strong>de</strong> típica 0,850 Ah<br />
Esta<strong>do</strong> <strong>de</strong> carga 100%<br />
Capacida<strong>de</strong> máxima 0,8925 Ah<br />
Tensão com carga completa 4,1 V<br />
Corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga 0,36 A<br />
Resistência inter<strong>na</strong> 0,0255 Ω<br />
Capacida<strong>de</strong> <strong>na</strong> tensão <strong>de</strong> cutoff 0,7551 Ah<br />
Tensão exponencial 3,82 V<br />
Capacida<strong>de</strong> exponencial 0,3 Ah<br />
Tempo <strong>de</strong> resposta 30 s<br />
Unida<strong>de</strong> no eixo x s<br />
Após serem informa<strong>do</strong>s os parâmetros da Tabela 5.1, ao mo<strong>de</strong>lo Battery, foram realiza-<br />
das as simulações computacio<strong>na</strong>is. Tanto para a simulação em <strong>de</strong>scarga contínua, quanto<br />
para a simulação em <strong>de</strong>scarga variável, foi utiliza<strong>do</strong> um programa <strong>na</strong> forma <strong>de</strong> diagrama<br />
<strong>de</strong> blocos, implementa<strong>do</strong> no MatLab/Simulink, para a geração e posterior tratamento das<br />
curvas simuladas pelo mo<strong>de</strong>lo. Este programa foi apresenta<strong>do</strong> e <strong>de</strong>scrito <strong>de</strong>talhadamente<br />
no capítulo anterior.<br />
Na Figura 5.4 são apresentadas as curvas <strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s consi<strong>de</strong>ran<strong>do</strong> a<br />
<strong>de</strong>scarga contínua, bem como a curva da corrente contínua utilizada. Na Figura 5.5 são<br />
apresentadas as curvas <strong>do</strong>s resulta<strong>do</strong>s simula<strong>do</strong>s consi<strong>de</strong>ran<strong>do</strong> a <strong>de</strong>scarga variável, bem<br />
como os <strong>de</strong>graus <strong>de</strong> correntes consi<strong>de</strong>ra<strong>do</strong>s para simulação. Observa-se a partir das curvas<br />
encontradas que os resulta<strong>do</strong>s simula<strong>do</strong>s pelo mo<strong>de</strong>lo Battery estão muito próximos <strong>do</strong>s<br />
resulta<strong>do</strong>s obti<strong>do</strong>s pelo mo<strong>de</strong>lo para Predizer Runtime e Características V-I.<br />
5.3 Comparação entre os Mo<strong>de</strong>los <strong>Elétricos</strong><br />
Uma comparação entre os resulta<strong>do</strong>s <strong>de</strong> simulações <strong>de</strong>stes mo<strong>de</strong>los elétricos está apre-<br />
sentada <strong>na</strong> Tabela 5.2, consi<strong>de</strong>ran<strong>do</strong> uma <strong>de</strong>scarga contínua <strong>de</strong> 80 mA e uma <strong>de</strong>scarga
Capítulo 5. Análise Comparativa <strong>do</strong>s Mo<strong>de</strong>los 64<br />
Figura 5.4: Curva simulada para uma <strong>de</strong>scarga contínua <strong>de</strong> 80 mA.<br />
Figura 5.5: Curva simulada para uma <strong>de</strong>scarga <strong>de</strong> quatro <strong>de</strong>graus <strong>de</strong> corrente.<br />
variável com quatro <strong>de</strong>graus <strong>de</strong> corrente. A partir da análise <strong>do</strong>s resulta<strong>do</strong>s observa-se<br />
que para o perfil <strong>de</strong> <strong>de</strong>scarga contínua, o erro no tempo <strong>de</strong> vida apresenta<strong>do</strong> pelo mo<strong>de</strong>lo
Capítulo 5. Análise Comparativa <strong>do</strong>s Mo<strong>de</strong>los 65<br />
Battery quan<strong>do</strong> compara<strong>do</strong> com os da<strong>do</strong>s experimentais <strong>de</strong> [1] foi <strong>de</strong> 0,534%, enquanto<br />
que para o mo<strong>de</strong>lo para Predizer Runtime e Característica V-I o erro foi <strong>de</strong> 0,395%. Para<br />
o perfil <strong>de</strong> <strong>de</strong>scarga variável, o erro no tempo <strong>de</strong> vida apresenta<strong>do</strong> pelo mo<strong>de</strong>lo Battery<br />
quan<strong>do</strong> compara<strong>do</strong> com os da<strong>do</strong>s experimentais <strong>de</strong> [1] foi <strong>de</strong> 1,621%, enquanto que para<br />
o mo<strong>de</strong>lo para Predizer Runtime e Característica V-I o erro foi <strong>de</strong> 0,338%. Verifica-se<br />
então, uma diferença nos erros <strong>do</strong> tempo <strong>de</strong> vida <strong>do</strong>s mo<strong>de</strong>los <strong>de</strong> 0,139% para o caso da<br />
<strong>de</strong>scarga contínua <strong>de</strong> 80 mA, e <strong>de</strong> 1,283% para o caso da <strong>de</strong>scarga em quatro <strong>de</strong>graus <strong>de</strong><br />
corrente.<br />
Tabela 5.2: Comparação entre o mo<strong>de</strong>lo para Predizer Runtime e Característica V-I e o mo<strong>de</strong>lo<br />
Battery.<br />
Mo<strong>de</strong>lo para Predizer Runtime e Característica V-I Mo<strong>de</strong>lo Battery<br />
Perfis (mA) Tve (s) Tvs (s) Erro (%) Tvs (s) Erro (%)<br />
Descarga contínua 35586 35446 0,395 35396 0,534<br />
Descarga variável 11846 11886 0,338 11654 1,621<br />
Tve = Tempo <strong>de</strong> vida experimental, Tvs = Tempo <strong>de</strong> vida simula<strong>do</strong><br />
Conforme o espera<strong>do</strong>, o mo<strong>de</strong>lo para Predizer Runtime e Características V-I <strong>de</strong> uma<br />
<strong>Bateria</strong> apresentou melhores resulta<strong>do</strong>s, pois trata-se <strong>de</strong> um mo<strong>de</strong>lo <strong>de</strong> alta acurácia, pois<br />
engloba as características <strong>de</strong> três categorias <strong>de</strong> mo<strong>de</strong>los elétricos, <strong>de</strong>mandan<strong>do</strong> um pro-<br />
cesso <strong>de</strong> extração <strong>de</strong> parâmetros mais <strong>de</strong>talha<strong>do</strong>, (i.e., obrigatorieda<strong>de</strong> <strong>de</strong> experimentos<br />
físicos e interpolações), e portanto, mais refi<strong>na</strong><strong>do</strong>. No entanto, os resulta<strong>do</strong>s apresenta-<br />
<strong>do</strong>s pelo mo<strong>de</strong>lo Battery são satisfatórios pois ficaram bastante próximos <strong>do</strong>s resulta<strong>do</strong>s<br />
apresenta<strong>do</strong>s pelo mo<strong>de</strong>lo para Predizer Runtime e Características V-I, consi<strong>de</strong>ran<strong>do</strong> que<br />
para as simulações <strong>do</strong> mo<strong>de</strong>lo Battery não foi necessária a realização <strong>de</strong> nenhum teste ex-<br />
perimental com a bateria simulada, sen<strong>do</strong> feita a extração <strong>de</strong> parâmetros para o mo<strong>de</strong>lo<br />
diretamente <strong>de</strong> uma curva característica <strong>de</strong> <strong>de</strong>scarga presente no datasheet da bateria, o<br />
que po<strong>de</strong> ocasio<strong>na</strong>r uma maior possibilida<strong>de</strong> <strong>de</strong> erros durante a <strong>de</strong>finição <strong>do</strong>s parâmetros<br />
<strong>de</strong> configuração <strong>do</strong> mo<strong>de</strong>lo.<br />
Um fator importante a ser mencio<strong>na</strong><strong>do</strong>, que teve influência nos bons resulta<strong>do</strong>s apre-<br />
senta<strong>do</strong>s pelo mo<strong>de</strong>lo Battery, são as características consistentes da bateria <strong>de</strong> íon lítio<br />
polímero PL-383562, conforme cita<strong>do</strong> em [1]. Devi<strong>do</strong> a esta consistência, a <strong>de</strong>finição <strong>de</strong><br />
parâmetros para o mo<strong>de</strong>lo elétrico Battery, a partir <strong>de</strong> uma curva <strong>de</strong> <strong>de</strong>scarga presente em<br />
datasheet, tor<strong>na</strong>-se mais precisa, contribuin<strong>do</strong> <strong>de</strong>sta forma, para uma acurácia satisfatória<br />
<strong>do</strong> mo<strong>de</strong>lo.
Capítulo 6<br />
Conclusões e Trabalhos Futuros<br />
Neste trabalho é apresenta<strong>do</strong> um estu<strong>do</strong> <strong>de</strong> um mo<strong>de</strong>lo elétrico <strong>de</strong>nomi<strong>na</strong><strong>do</strong> Battery,<br />
presente <strong>na</strong> ferramenta computacio<strong>na</strong>l MatLab/Simulink, que simula a <strong>de</strong>scarga <strong>do</strong>s mais<br />
populares tipos <strong>de</strong> baterias recarregáveis (i.e., íon lítio, níquel cádmio, níquel metal hidreto<br />
e chumbo áci<strong>do</strong>). O principal motivo da escolha <strong>de</strong>ste mo<strong>de</strong>lo, se <strong>de</strong>ve ao fato <strong>do</strong> mesmo<br />
ser <strong>de</strong> fácil implementação no que se refere a extração <strong>de</strong> seus parâmetros <strong>de</strong> configuração,<br />
e principalmente, porque em alguns casos, a fase <strong>de</strong> testes experimentais para a extração<br />
<strong>de</strong> parâmetros <strong>do</strong> mo<strong>de</strong>lo po<strong>de</strong> ser ignorada. Além disto, o mo<strong>de</strong>lo consi<strong>de</strong>ra parte <strong>do</strong>s<br />
efeitos não-lineares ocorri<strong>do</strong>s em um processo real <strong>de</strong> <strong>de</strong>scarga <strong>de</strong> uma bateria.<br />
Todas as simulações computacio<strong>na</strong>is com o emprego <strong>do</strong> mo<strong>de</strong>lo Battery, foram realiza-<br />
das <strong>na</strong> ferramenta computacio<strong>na</strong>l MatLab/Simulink, a partir <strong>de</strong> um circuito implementa<strong>do</strong><br />
<strong>na</strong> forma <strong>de</strong> diagrama <strong>de</strong> blocos. Primeiramente para a obtenção <strong>do</strong>s resulta<strong>do</strong>s das si-<br />
mulações foram consi<strong>de</strong>ra<strong>do</strong>s os parâmetros <strong>de</strong> uma bateria <strong>de</strong> íon lítio Nokia BL-5F, e<br />
posteriormente, os parâmetros da bateria <strong>de</strong> íon lítio polímero PL-383562.<br />
O estu<strong>do</strong> <strong>do</strong> mo<strong>de</strong>lo elétrico Battery foi dividi<strong>do</strong> em <strong>do</strong>is momentos: primeiramente foi<br />
realizada uma comparação <strong>do</strong>s seus resulta<strong>do</strong>s simula<strong>do</strong>s, consi<strong>de</strong>ran<strong>do</strong> os parâmetros <strong>de</strong><br />
uma bateria <strong>de</strong> íon lítio Nokia BL-5F, com os resulta<strong>do</strong>s experimentais obti<strong>do</strong>s a partir <strong>de</strong><br />
uma plataforma <strong>de</strong> testes; em um segun<strong>do</strong> momento foi realizada uma comparação <strong>de</strong> seus<br />
resulta<strong>do</strong>s simula<strong>do</strong>s, consi<strong>de</strong>ran<strong>do</strong> os parâmetros <strong>de</strong> uma bateria <strong>de</strong> íon lítio polímero<br />
PL-383562, com resulta<strong>do</strong>s simula<strong>do</strong>s a partir <strong>do</strong> mo<strong>de</strong>lo elétrico para Predizer Runtime<br />
e Característica V-I, que é consi<strong>de</strong>ra<strong>do</strong> <strong>na</strong> literatura um mo<strong>de</strong>lo elétrico <strong>de</strong> alta acurácia.<br />
Nos resulta<strong>do</strong>s das simulações <strong>do</strong> mo<strong>de</strong>lo elétrico Battery a partir da comparação <strong>de</strong><br />
seus resulta<strong>do</strong>s simula<strong>do</strong>s com os da<strong>do</strong>s experimentais, consi<strong>de</strong>ran<strong>do</strong> tanto <strong>de</strong>scargas con-<br />
tínuas quanto <strong>de</strong>scargas variáveis para uma bateria <strong>de</strong> íon lítio Nokia BL-5F, verificou-se,<br />
para <strong>de</strong>scargas contínuas, que o mo<strong>de</strong>lo apresentou resulta<strong>do</strong>s satisfatórios <strong>na</strong> maioria das<br />
correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga ou calibração <strong>do</strong> mo<strong>de</strong>lo, resulta<strong>do</strong>s estes com um erro no<br />
tempo <strong>de</strong> vida abaixo <strong>de</strong> 5%, o que é espera<strong>do</strong> para mo<strong>de</strong>los elétricos. O menor erro no<br />
66
Capítulo 6. Conclusões e Trabalhos Futuros 67<br />
tempo <strong>de</strong> vida, para o caso <strong>de</strong> <strong>de</strong>scargas contínuas, foi <strong>de</strong> 2,12% e ocorreu <strong>na</strong> corrente<br />
nomi<strong>na</strong>l <strong>de</strong> calibração <strong>de</strong> 650 mA, indican<strong>do</strong> que o mo<strong>de</strong>lo apresentou o melhor resulta<strong>do</strong><br />
a partir da corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga equivalente a aproximadamente 68% da capaci-<br />
da<strong>de</strong> nomi<strong>na</strong>l da bateria <strong>de</strong> íon lítio Nokia BL-5F, que é <strong>de</strong> 950 mAh. Já para o caso <strong>de</strong><br />
<strong>de</strong>scargas variáveis, o menor erro no tempo <strong>de</strong> vida foi <strong>de</strong> 1,79% e ocorreu <strong>na</strong>s correntes<br />
nomi<strong>na</strong>is <strong>de</strong> calibração <strong>de</strong> 250 mA e 350 mA, indican<strong>do</strong> que os melhores resulta<strong>do</strong>s <strong>do</strong><br />
mo<strong>de</strong>lo foram obti<strong>do</strong>s a partir das correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga equivalentes a aproxi-<br />
madamente 26% e 37% da capacida<strong>de</strong> nomi<strong>na</strong>l da bateria utilizada. Os piores resulta<strong>do</strong>s<br />
apresenta<strong>do</strong>s pelo mo<strong>de</strong>lo, tanto <strong>na</strong>s <strong>de</strong>scargas contínuas quanto <strong>na</strong>s <strong>de</strong>scargas variáveis,<br />
ocorreram <strong>na</strong> menor corrente nomi<strong>na</strong>l <strong>de</strong> <strong>de</strong>scarga <strong>do</strong> mo<strong>de</strong>lo (i.e., 100 mA), pois é justa-<br />
mente em correntes baixas <strong>de</strong> <strong>de</strong>scarga que os efeitos não-lineares se fazem mais presentes<br />
em um processo real <strong>de</strong> <strong>de</strong>scarga, fazen<strong>do</strong> com que o processo <strong>de</strong> extração <strong>de</strong> parâmetros,<br />
para o mo<strong>de</strong>lo Battery, que <strong>de</strong>pen<strong>de</strong> <strong>de</strong> uma curva característica <strong>de</strong> <strong>de</strong>scarga, tor<strong>na</strong>-se<br />
menos preciso, diminuin<strong>do</strong> <strong>de</strong>sta forma, a acurácia <strong>do</strong> mo<strong>de</strong>lo. No entanto, verificou-<br />
se que em uma corrente nomi<strong>na</strong>l <strong>de</strong> calibração intermediária (i.e., 450 mA), tanto <strong>na</strong>s<br />
<strong>de</strong>scargas contínuas, quanto <strong>na</strong>s <strong>de</strong>scargas variáveis, o mo<strong>de</strong>lo apresentou um resulta<strong>do</strong><br />
satisfatório, próximo <strong>do</strong> melhor resulta<strong>do</strong> apresenta<strong>do</strong> pelo mo<strong>de</strong>lo. Portanto, não é su-<br />
geri<strong>do</strong> que o mo<strong>de</strong>lo elétrico Battery seja calibra<strong>do</strong> com correntes nomi<strong>na</strong>is <strong>de</strong> <strong>de</strong>scarga<br />
altas e baixas, mas sim, com uma corrente nomi<strong>na</strong>l <strong>de</strong> calibração intermediária em um<br />
valor próximo <strong>de</strong> 50% da capacida<strong>de</strong> nomi<strong>na</strong>l da bateria simulada.<br />
A partir da comparação realizada entre os resulta<strong>do</strong>s das simulações <strong>do</strong> mo<strong>de</strong>lo Bat-<br />
tery e os resulta<strong>do</strong>s das simulações <strong>do</strong> mo<strong>de</strong>lo para predizer Runtime e Características V-I<br />
<strong>de</strong> uma <strong>Bateria</strong>, consi<strong>de</strong>ran<strong>do</strong> os mesmos cenários simula<strong>do</strong>s no caso anterior, verificou-se<br />
que apesar <strong>do</strong>s parâmetros <strong>do</strong> mo<strong>de</strong>lo Battery terem si<strong>do</strong> obti<strong>do</strong>s diretamente <strong>do</strong> datasheet<br />
da bateria a ser simulada, ou seja, não foi necessária a realização <strong>de</strong> nenhum teste ex-<br />
periemental com a bateria <strong>de</strong> íon lítio polímero PL-383562, o mo<strong>de</strong>lo apresentou resulta<strong>do</strong>s<br />
satisfatórios. As diferenças, entre os mo<strong>de</strong>los, nos erros <strong>de</strong> tempo <strong>de</strong> vida foram <strong>de</strong> ape-<br />
<strong>na</strong>s 0,139% para <strong>de</strong>scarga contínua e <strong>de</strong> 1,283% para <strong>de</strong>scarga variável. Estes resulta<strong>do</strong>s<br />
satisfatórios apresenta<strong>do</strong>s pelo mo<strong>de</strong>lo Battery, sem a necessida<strong>de</strong> da realização <strong>de</strong> testes<br />
experimentais com a bateria simulada, representam uma vantagem <strong>do</strong> mo<strong>de</strong>lo <strong>na</strong> questão<br />
<strong>de</strong> otimização <strong>de</strong> tempo, que é tão <strong>de</strong>sejada nos dias atuais <strong>na</strong>s mais diversas áreas <strong>de</strong><br />
engenharia.<br />
Como trabalhos futuros, sugere-se esten<strong>de</strong>r este trabalho realizan<strong>do</strong> simulações con-<br />
si<strong>de</strong>ran<strong>do</strong> outras baterias <strong>de</strong> íon lítio, objetivan<strong>do</strong> verificar se os resulta<strong>do</strong>s apresenta<strong>do</strong>s<br />
pelo mo<strong>de</strong>lo Battery vão continuar ocorren<strong>do</strong> <strong>na</strong>s mesmas correntes nomi<strong>na</strong>is <strong>de</strong> cali-<br />
bração. E posteriormente, esten<strong>de</strong>r este trabalho realizan<strong>do</strong> simulações com outras tec-<br />
nologias <strong>de</strong> baterias recarregáveis, também objetivan<strong>do</strong> verificar em qual calibração são
Capítulo 6. Conclusões e Trabalhos Futuros 68<br />
apresenta<strong>do</strong>s os melhores resulta<strong>do</strong>s pelo mo<strong>de</strong>lo. E por fim, sugere-se esten<strong>de</strong>r este<br />
trabalho realizan<strong>do</strong> uma análise estatística <strong>do</strong>s resulta<strong>do</strong>s encontra<strong>do</strong>s.
Referências Bibliograficas<br />
[1] M. Chen and G. Rincón-Mora, “Accurate electrical battery mo<strong>de</strong>l capable of predict-<br />
ing runtime and i-v performance,” IEEE Transactions on Energy Conversion, vol. 21,<br />
no. 2, june 2006.<br />
[2] K. K. Schnei<strong>de</strong>r, “Avaliação <strong>de</strong> mo<strong>de</strong>los a<strong>na</strong>líticos <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong><br />
baterias utilizadas em dispositivos móveis,” Dissertação <strong>de</strong> Mestra<strong>do</strong>, Universida<strong>de</strong><br />
Regio<strong>na</strong>l <strong>do</strong> Noroeste <strong>do</strong> Esta<strong>do</strong> <strong>do</strong> Rio Gran<strong>de</strong> <strong>do</strong> Sul, Ijuí, Rio Gran<strong>de</strong> <strong>do</strong> Sul,<br />
Brasil, 2011.<br />
[3] M. Jonger<strong>de</strong>n and B. Haverkort, “Which battery mo<strong>de</strong>l to use?” Imperial College<br />
Lon<strong>do</strong>n, pp. 76–88, July 2008.<br />
[4] P. S. Sausen, “Gerenciamento integra<strong>do</strong> <strong>de</strong> energia e controle <strong>de</strong> topologia em re<strong>de</strong>s<br />
<strong>de</strong> sensores sem fio,” Master’s Thesis, Universida<strong>de</strong> Fe<strong>de</strong>ral <strong>de</strong> Campi<strong>na</strong> Gran<strong>de</strong>,<br />
Campi<strong>na</strong> Gran<strong>de</strong>, Paraíba, Brasil, july 2008.<br />
[5] M. Jonger<strong>de</strong>n and B. Haverkort, “Battery mo<strong>de</strong>ling,” Thecnical Report in Faculty<br />
Electrical Engineering, Mathematics and Computer Science, 2008.<br />
[6] N. Parthasarathy, “Adaptive battery monitoring using parameter estimation,” Mas-<br />
ter’s Thesis, The Ohio State University, Columbus, Ohio, EUA, 2009.<br />
[7] S. D. K. Lahiri, A. Raghu<strong>na</strong>than and D. Panigrahi, “Battery-driven system <strong>de</strong>sign:<br />
A new frontier in low power <strong>de</strong>sign,” Chemical Engineering Science.<br />
[8] N. Kularat<strong>na</strong>, Power Electronics Design Handbook. EDN Series for Design Engineers,<br />
1998.<br />
[9] T. F. F. M. Doyle and J. Newman, “Mo<strong>de</strong>ling of galvanostatic charge and discharge<br />
of the lithium, polymer, insertion cell,” Jour<strong>na</strong>l of the Electrochemical Society, vol.<br />
140, no. 6, pp. 1526–1533, 1993.<br />
69
REFERÊNCIAS BIBLIOGRAFICAS 70<br />
[10] M. D. T. F. Fuller and J. Newman, “Simulation and optimization of the dual lithium<br />
ion insertion cell,” Jour<strong>na</strong>l of the Electrochemical Society, vol. 141, no. 1, pp. 1–10,<br />
1994.<br />
[11] Simula<strong>do</strong>r Fortran Dualfoil, Disponível em http://www.cchem.berkeley.edu/jsngrp/<br />
fortran.html.<br />
[12] S. D. R. R. A. R. D. Panigrahi, C. Chiasserini and K. Lahiri, “Battery life estimation<br />
of mobile embed<strong>de</strong>d systems,” NEC USA Inc., and the California Micro Program.<br />
[13] C. Chiasserini and R. Rao, “Pulsed battery discharge in communication <strong>de</strong>vices,” Pro-<br />
ceedings of the 5th Inter<strong>na</strong>tio<strong>na</strong>l Conference on Mobile Computing and Networking,<br />
pp. 88–95, 1999.<br />
[14] “A mo<strong>de</strong>l for battery pulsed discharge with recovery effect,” Wireless Communica-<br />
tions and Networking Conference, pp. 636–639, 1999.<br />
[15] D. Rakhmatov and S. Vrudhula, “An a<strong>na</strong>lytical high-level battery mo<strong>de</strong>l for use in<br />
energy ma<strong>na</strong>gement of portable electronic systems,” Natio<strong>na</strong>l Science Foundation’s<br />
State/Industry/University Cooperative Research Centers’ (NSFS/IUCRC) Center for<br />
Low Power Electronics (CLPE), pp. 1–6, 2001.<br />
[16] P. S. S. e. A. S. K. K. Schnei<strong>de</strong>r, “Análise comparativa <strong>de</strong> mo<strong>de</strong>los a<strong>na</strong>líticos para<br />
predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> baterias em dispositivos móveis,” 9th Conference on<br />
Dy<strong>na</strong>mics, Control and their Applications, june 2010.<br />
[17] D. Rakhmatov and S. Vrudhula, “Energy ma<strong>na</strong>gement for battery-powered embed<strong>de</strong>d<br />
systems,” ACM Transactions on Embed<strong>de</strong>d Computing Systems, vol. 2, no. 3, pp.<br />
277–324, August 2003.<br />
[18] W. A. L. Z. M. Salameh, M. A. Casacca, “A mathematical mo<strong>de</strong>l for lead-acid bat-<br />
teries,” IEEE Transactions on Energy Conversion, vol. 7, no. 1, pp. 93–98, march<br />
1992.<br />
[19] D. R. S. R. M. Valvo, F. E. Wicks, “Development and application of an improved<br />
equivalent circuit mo<strong>de</strong>l of a lead acid battery,” Proceedings Energy Conversion Eng.<br />
Conference, vol. 2, pp. 1159–1163, august 1996.<br />
[20] M. C. S. Barsali, “Dy<strong>na</strong>mical mo<strong>de</strong>ls of lead-acid batteries: Implementation issues,”<br />
IEEE Transactions on Energy Conversion, vol. 17, no. 1, pp. 16–23, march 2002.
REFERÊNCIAS BIBLIOGRAFICAS 71<br />
[21] R. D. S. Buller, M. Thele and E. Kar<strong>de</strong>n, “Impedance based simulation mo<strong>de</strong>ls of<br />
super-capacitors and li-ion batteries for power electronic applications,” Industrial<br />
Application Conference, vol. 3, 2003.<br />
[22] MathWorks, "Implement generic battery mo<strong>de</strong>l", Disponível em<br />
http://www.mathworks.nl/access/help<strong>de</strong>sk/help/toolbox/physmod/powersys/ref/<br />
battery.html, 2010.<br />
[23] Y. E. Matsumoto, Simulink 7.2 - Guia Prático, 1st ed. Belo Horizonte: Editora<br />
Érica Ltda, 2008.<br />
[24] L. M. H. B. Nonemacher and P. S. Sausen, “Desenvolvimento <strong>de</strong> um testbed para<br />
avaliação <strong>de</strong> mo<strong>de</strong>los matemáticos utiliza<strong>do</strong>s <strong>na</strong> predição <strong>do</strong> tempo <strong>de</strong> vida <strong>de</strong> bate-<br />
rias,” Congresso Regio<strong>na</strong>l <strong>de</strong> Iniciação Científica e Tecnológica em Engenharia, pp.<br />
2–4, 2010.<br />
[25] Lithium Ion Rechargeable Batteries - Technical Handbook, Disponível em<br />
http://www.sony.com.cn/products/ed/battery/<strong>do</strong>wnload.pdf, 2011.
Apêndice A<br />
Lista das Publicações Relacio<strong>na</strong>das com<br />
a Dissertação<br />
A.1 Artigo Aceito em Periódico Nacio<strong>na</strong>l<br />
1. C. M. D. Porciuncula, A. V. Oliveira, P. S. Sausen, A. Sausen, Avaliação Comparativa<br />
entre o Mo<strong>de</strong>lo Elétrico Battery e os Mo<strong>de</strong>los A<strong>na</strong>líticos Linear e Lei <strong>de</strong> Peukert. Revista<br />
Brasileira <strong>de</strong> Computação Aplicada (ISSN 2176-6649), Passo Fun<strong>do</strong>, v. 3, n. 2, p. xx-xx,<br />
mar. 2012.<br />
A.2 Artigos Publica<strong>do</strong>s em A<strong>na</strong>is<br />
1. C. M. D. Porciuncula, A. V. Oliveira, P. S. Sausen, A. Sausen, Mo<strong>de</strong>lo Elétrico Battery<br />
Versus Mo<strong>de</strong>los A<strong>na</strong>líticos. A<strong>na</strong>is <strong>do</strong> EREMATSUL 2011 - XVII Encontro Regio<strong>na</strong>l <strong>de</strong><br />
Estudantes <strong>de</strong> Matemática da Região Sul. Curitiba, PR, Brasil, Setembro <strong>de</strong> 2011.<br />
2. C. M. D. Porciuncula, A. V. Oliveira, P. S. Sausen, A. Sausen, Avaliação Compar-<br />
ativa entre o Mo<strong>de</strong>lo Elétrico Battery e os Mo<strong>de</strong>los A<strong>na</strong>líticos Linear e Lei <strong>de</strong> Peukert.<br />
A<strong>na</strong>is <strong>do</strong> SCA 2011 - III Simpósio <strong>de</strong> Computação Aplicada. Passo Fun<strong>do</strong>, RS, Brasil,<br />
Outubro <strong>de</strong> 2011.<br />
A.3 Resumos Expandi<strong>do</strong>s Publica<strong>do</strong>s em A<strong>na</strong>is<br />
1. C. M. D. Porciuncula, A. V. Oliveira, P. S. Sausen, A. Sausen, <strong>Aplicação</strong> <strong>do</strong> Mo<strong>de</strong>lo<br />
Elétrico Battery <strong>na</strong> Estimação <strong>do</strong> Tempo <strong>de</strong> Vida <strong>de</strong> <strong>Bateria</strong>s Utilizadas em Dispositivos<br />
Móveis. A<strong>na</strong>is <strong>do</strong> DINCON 2011 - 1 a Conferência Brasileira <strong>de</strong> Dinâmica, Controle e<br />
72
Apêndice A. Lista das Publicações Relacio<strong>na</strong>das com a Dissertação 73<br />
Aplicações. Águas <strong>de</strong> Lindóia, SP, Brasil, 28 <strong>de</strong> agosto a 01 <strong>de</strong> setembro <strong>de</strong> 2011.<br />
2. A. V. Oliveira, C. M. D. Porciuncula, P. S. Sausen, A. Sausen, Estimação <strong>do</strong>s Parâ-<br />
metros <strong>do</strong> Mo<strong>de</strong>lo A<strong>na</strong>lítico <strong>de</strong> Difusão <strong>de</strong> Rakhmatov-Vrudhula para a <strong>Predição</strong> <strong>do</strong> Tempo<br />
<strong>de</strong> Vida <strong>de</strong> <strong>Bateria</strong>s. A<strong>na</strong>is <strong>do</strong> DINCON 2011 - 1 a Conferência Brasileira <strong>de</strong> Dinâmica,<br />
Controle e Aplicações. Águas <strong>de</strong> Lindóia, SP, Brasil, 28 <strong>de</strong> agosto a 01 <strong>de</strong> setembro <strong>de</strong><br />
2011.<br />
A.4 Resumos Publica<strong>do</strong>s em A<strong>na</strong>is<br />
1. C. M. D. Porciuncula, P. S. Sausen, A. Sausen, A. V. Oliveira, <strong>Aplicação</strong> <strong>de</strong> Mo<strong>de</strong>los<br />
<strong>Elétricos</strong> <strong>na</strong> Estimação <strong>do</strong> Tempo <strong>de</strong> Vida <strong>de</strong> <strong>Bateria</strong>s Utilizadas em Dispositivos Móveis.<br />
A<strong>na</strong>is <strong>do</strong> Seminário Inter<strong>na</strong>cio<strong>na</strong>l <strong>de</strong> Educação <strong>do</strong> Mercosul 2011 - XIII Seminário Inter-<br />
<strong>na</strong>cio<strong>na</strong>l <strong>de</strong> Educação <strong>do</strong> Mercosul, X Seminário Interinstitucio<strong>na</strong>l, I Curso <strong>de</strong> Práticas<br />
Sócioculturais Interdiscipli<strong>na</strong>res. Cruz Alta, RS, Brasil, Junho <strong>de</strong> 2011.<br />
2. C. M. D. Porciuncula, P. S. Sausen, A. Sausen, Estu<strong>do</strong> e <strong>Aplicação</strong> <strong>de</strong> Mo<strong>de</strong>los<br />
<strong>Elétricos</strong> <strong>de</strong> Descarga <strong>na</strong> Estimação <strong>do</strong> Tempo <strong>de</strong> Vida <strong>de</strong> <strong>Bateria</strong>s Utilizadas em Dispo-<br />
sitivos Eletrônicos Móveis. A<strong>na</strong>is <strong>do</strong> CNEM 2011 - II Congresso Nacio<strong>na</strong>l <strong>de</strong> Educação<br />
Matemática e IX Encontro Regio<strong>na</strong>l <strong>de</strong> Educação Matemática. Ijuí, RS, Brasil, Junho <strong>de</strong><br />
2011.<br />
3. A. V. Oliveira, P. S. Sausen, A. Sausen, C. M. D. Porciuncula, Estimação <strong>do</strong>s Pa-<br />
râmetros Não-Lineares <strong>do</strong> Mo<strong>de</strong>lo Matemático <strong>de</strong> Difusão <strong>de</strong> Rakhmatov-Vrudhula Uti-<br />
lizan<strong>do</strong> Técnicas <strong>de</strong> I<strong>de</strong>ntificação <strong>de</strong> Sistemas. A<strong>na</strong>is CT&I e SOCIEDADE - XVIII<br />
Seminário <strong>de</strong> Iniciação Científica, XV Jor<strong>na</strong>da <strong>de</strong> Pesquisa, XI Jor<strong>na</strong>da <strong>de</strong> Extensão.<br />
Ijuí, RS, Brasil, Outubro <strong>de</strong> 2010.<br />
4. C. M. D. Porciuncula, P. S. Sausen, A. Sausen, Avaliação Comparativa entre o<br />
Mo<strong>de</strong>lo Elétrico Battery e Da<strong>do</strong>s Experimentais <strong>de</strong> uma <strong>Bateria</strong> <strong>de</strong> Íon Lítio. A<strong>na</strong>is <strong>do</strong><br />
Salão <strong>do</strong> Conhecimento 2011 - XVI Jor<strong>na</strong>da <strong>de</strong> Pesquisa <strong>Unijuí</strong>. Ijuí, RS, Brasil, Outubro<br />
<strong>de</strong> 2011.
Apêndice A. Lista das Publicações Relacio<strong>na</strong>das com a Dissertação 74<br />
A.5 Artigo Submeti<strong>do</strong><br />
1. Foi submeti<strong>do</strong> um artigo no XIX Congresso Brasileiro <strong>de</strong> Automática (CBA 2012),<br />
com os resulta<strong>do</strong>s fi<strong>na</strong>is <strong>de</strong>sta dissertação