ANÁLISE DO SOFTWARE PSIM POR MEIO DA IMPLEMENTAÇÃO

COMPUTACIONAL DE MÉTODOS CLÁSSICOS DE MPPT PARA PAINEIS

FOTOVOLTAICOS
Leonardo Rosenthal Caetano Silva1, Marcelo Pansani Freitas2, Marcio Augusto Tamashiro3, Anderson
Rodrigo Piccini4, Geraldo Caixeta Guimarães5, Jaqueline Oliveira Rezende6, Larissa Marques Peres7
Universidade Federal de Uberlândia, Uberlândia – MG, 2 Instituto Federal de Educação Ciências e Tecnologia do
1,5,6,7

Triângulo Mineiro, Paracatu – MG, 3,4 Instituto Federal de Educação Ciências e Tecnologia do Tocantins, Palmas – TO
1
leonardo.rosenthal@hotmail.com, 2marcelopansani@iftm.edu.br, 3tamashiro@ifto.edu.br, 4anderson@ifto.edu.br,
5
gcaixeta@ufu.br, 6jaqueline.oliveirarezende@gmail.com, 7larissaa_mp@yahoo.com.br

Resumo - Este artigo apresenta a implementação de aqui é mínima. Conforme observado nos informes gerais da
dois métodos clássicos de MPPT bastante conhecidos, por Agência Nacional de Energia Elétrica (ANEEL) de junho de
meio do software comercial PSIM. Pretende-se assim 2013, tal geração não chega a 0,1% da demanda geral de
testar as principais funcionalidades desse programa, além energia do Brasil [3]. Outros países do mundo, mesmo em
de servir como um guia prático de uso inicial destinado a condições naturais mais desfavoráveis estão mais avançados
alunos e professores de Engenharia Elétrica. na exploração dessa fonte renovável, como exemplo:
Alemanha, Itália e Espanha, onde em 2011 apresentaram um
Palavras-Chave – Método da condutância incremental, consumo anual de energia solar de 37,5 Terawatt-hora
método perturba e observa, MPPT, PSIM. (TWh), representando 67,4% de todo o consumo mundial
[4].
ANALYSIS OF PSIM SOFTWARE VIA Uma resolução publicada pela ANEEL, em abril de 2012,
COMPUTER IMPLEMENTATION OF deverá estimular a produção de energia elétrica no Brasil por
meio da geração fotovoltaica. Trata-se da Resolução
CLASSICAL METHODS OF MPPT FOR Normativa (REN) 482/2012 da ANEEL [5], que define o
PHOTOVOLTAIC PANELS Sistema de Compensação como um arranjo no qual a energia
ativa injetada por alguma unidade consumidora, com
Abstract - This paper presents the implementation of microgeração ou minigeração distribuída, é cedida à
two well-known classical methods of MPPT, through of distribuidora local e posteriormente compensada com o
the commercial software PSIM. The aim is to test the consumo de energia elétrica ativa dessa mesma unidade
main features of this program, and also to serve as a consumidora ou de outra unidade de mesma titularidade.
practical guide for initial use for students and faculty of Apesar das vantagens da geração fotovoltaica através do
Electrical Engineering. uso de painéis fotovoltaicos, a eficiência desse tipo de
1 conversão de energia ainda é baixa, e desta forma, torna-se
Keywords – Incremental conductance method, perturb necessária a utilização de técnicas para a obtenção da
and observe method, MPPT, PSIM. máxima potência possível (MPPT – Maximum Power Point
Tracking) desses painéis [6][7].
Vários métodos para a extração da máxima potência são
I. INTRODUÇÃO encontrados na literatura [6][7][8][9]. Como exemplo, tem-se
os métodos da Tensão Constante, da Tensão de Circuito
A quantidade de energia produzida por um sistema
Aberto, de Curto-Circuito por Pulsos, Perturba e Observa, da
fotovoltaico depende da taxa de irradiação solar e da
Condutância Incremental, Hill Climbing, Beta, Correlação de
temperatura do local onde é instalado. No Brasil, as regiões
Ripple, métodos baseados em inteligência artificial, entre
Nordeste, Centro-Oeste e Norte são as que possuem o maior
outros. Dentre esses, destacam-se dois métodos bastante
potencial de aproveitamento da energia solar. A região Sul
eficazes, quando comparados com os demais, [10] além de
dentre todas as regiões é a menos privilegiada, entretanto
serem de fácil implementação:
ainda possui taxas de irradiação melhores do que aquelas
 Método Perturba e Observa (P&O), e
encontradas na Alemanha [1] ou em outros países que
empregam largamente a energia solar fotovoltaica, conforme  Método da Condutância Incremental.
verificado no Atlas Solarimétrico do Brasil do ano 2000 [2].
Apesar da alta incidência de raios solares em todo o Este artigo irá apresentar os dois métodos clássicos acima
território nacional ao longo do ano, a energia solar gerada mencionados, bem como as suas implementações no
software PSIM. Pretende-se assim demonstrar as vantagens
desse software para a realização de estudos sobre sistemas
fotovoltaicos, além de servir, para estudantes e professores,
como guia inicial de utilização do PSIM.

1
 II. MÉTODOS CLÁSSICOS DE MPPT

A. Método Perturba e Observa

O método Perturba e Observa (P&O) opera de forma
periódica, incrementando ou decrementando a tensão de
saída terminal do painel fotovoltaico (PV), e comparando a
potência obtida no ciclo atual com a potência do ciclo
anterior [11]. Caso a potência aumente, o sistema de controle
muda o ponto de operação naquela direção; caso contrário,
muda o ponto de operação na direção oposta [11]. Além
disso, a tensão terminal do PV é sempre alterada e quando o
painel opera no ponto de máxima potência (MPP) ocorrem Fig. 2 Curva P-V com ponto de máxima potência.
oscilações na potência de saída em torno do seu valor
máximo, o que resulta em perdas de potência [11]. Os O fluxograma do algoritmo desse método de MPPT é
métodos de P&O podem ser divididos em P&O clássico, mostrado na Figura 3 [10].
otimizado e de três pontos [8]. O P&O clássico opera com
perturbações fixas; em relação ao método otimizado, uma
média de diversas amostras de potência é usada para ajustar
dinamicamente a magnitude da perturbação; já no P&O de
três pontos, são usados três pontos distintos da curva
potência versus tensão, para determinar a direção e
magnitude da próxima perturbação [8].
A Figura 1 [12] apresenta o fluxograma do algoritmo do
P&O clássico.

Fig. 3 Fluxograma do algoritmo do método de condutância

incremental clássico.

III. SOFTWARE PSIM

PSIM, segundo o desenvolvedor, é um dos simuladores

mais rápidos para o campo da eletrônica de potência. Esse
software utiliza um forte algoritmo dedicado a circuitos
Fig. 1 Fluxograma do algoritmo do método P&O clássico.
elétricos [13]. É capaz de simular circuitos de controle nas
B. Método da Condutância Incremental várias formas: circuito analógico, diagrama de blocos da
Este método é baseado no fato de que a inclinação da função de transferência no domínio s ou z, código C
curva de potência do painel solar é nula no MPP, positiva à personalizado, ou no Matlab/Simulink. A sua biblioteca de
esquerda e negativa à direita (vide Figura 2) [7]. Uma vez controle fornece uma lista abrangente de componentes e
alcançado o MPP, as perturbações são encerradas até que se blocos de funções, o que torna possível construir
observem alterações na corrente do painel. Desta forma, não praticamente qualquer esquema de controle de forma rápida e
ocorrem oscilações na potência do PV e o passo de conveniente.
incremento determina a velocidade do método [7]. Também Possui diversos módulos adicionais opcionais para atender
pode ser usado passo variável para melhorar a eficácia deste às necessidades específicas em diversas aplicações, tais como
método. Este método também é muito utilizado pela sua acionamentos de motores, controle digital, energia renovável,
eficácia e apresenta bom desempenho diante de variações ou suporte a DSP (Digital Signal Processing) e FPGA
rápidas das condições meteorológicas [8]. (Field-programmable Gate Array). Esses módulos fornecem
aos usuários a flexibilidade para adaptar o PSIM para suas

2
necessidades, e aumentar significativamente a sua
capacidade. Desta função de transferência, pelo método de Cohen-
Entre outras características do PSIM destacam-se [13]: Coon, foi calculado os valores da constante proporcional e da
• Interface gráfica de usuário intuitiva e de fácil constante de tempo integral de um controlador PI na forma
utilização; dependente, também utilizando o programa Calculadora de
• Interatividade que permite a alteração de parâmetros e PID. Os parâmetros do controlador PI encontrados foram,
visualização de tensões/correntes durante a simulação; respectivamente, Kp = 1722,494 e Tr = 0.00001. O gráfico
• Simulação de sistemas conversores de potência e da função controlada é visto na Figura 6.
controle grandes e complexos em um curto espaço de
tempo;
• Compilador C embutido que permite a digitação de
código C diretamente, sem compilar.
• Tem como requisitos de sistema um mínimo de 128MB
de memória RAM e cerca de 135MB de espaço em
disco disponível para a sua instalação.
• Roda no Microsoft Windows XP, Vista, 7 ou 8.

IV. IMPLEMENTAÇÃO COMPUTACIONAL DOS

MÉTODOS DE MPPT

A construção de circuitos elétricos, eletrônicos e de Fig. 6 Sistema fotovoltaico sob degrau sintonizado.
controle é semelhante, e tão fácil quanto a maneira utilizada
no Matlab/Simulink. Para a implementação dos dois métodos Tabela I – Blocos utilizados na implementação de cada MPPT
de MPPT foi utilizado o esquema indicado na Figura 4. Bloco Nome Função
Rotula e vincula pontos
comuns de sinais elétricos
Label
dos esquemas
implementados
Utilizado para definir a
frequência de atualização
Zero-order
do código do método P&O,
hold
cujo valor simulado é de
1000 Hz
Fig. 4 Esquema comum da implementação dos métodos de MPPT.
Disponível para escrita
simplificada do algoritmo
Os principais blocos utilizados estão listados na Tabela I. para solução de um
A implementação de cada um dos algoritmos de MPPT Simplified C determinado problema em
Block código C. No caso do
(ver Figuras 1 e 3) em código C encontram-se listados na
esquema implementado são
Tabela II do Anexo. utilizados duas entradas
Para a sintonização do bloco PI, nos modelos de P&O e de (x1 e x2) e uma saída (y1)
condutância incremental, foi aplicado um degrau de zero a Somador de sinais. No
mil unidades na irradiação solar do sistema. Pela medida da caso do esquema
Summer implementado, este bloco
potência máxima gerada nas fotocélulas e com a resposta de está funcionando como um
potência na saída da linha, pode se descobrir a função de subtrador
transferência do sistema, visto na Figura 5 e na Equação 1. Proportional-
Implementação
Foi utilizado o método de Hägglund, utilizando o software computacional do
Integral
controlador proporcional-
Calculadora de PID [14]. controller
integral (PI)

Limiter Limitador de sinal

Comparator Comparador de sinais

Fonte de tensão de onda

Triangular- triangular. No caso do
wave voltage esquema foi utilizado
source como referência de sinal
para o comparador
Terminal de referência de
Ground
Fig. 5 Sistema fotovoltaico sob degrau sem sintonia. tensão

(1)

3
 V. SIMULAÇÃO E RESULTADOS Tabela III – Blocos utilizados na implementação do sistema
fotovoltaico
Para a realização das simulações foi utilizado o mesmo Bloco Nome Função
sistema fotovoltaico com conversor CC-CC do tipo Buck Disponível para definir
configurações básicas da
(abaixador de tensão), apresentado na Figura 7. Esse circuito, simulação, como passo de
já com todos os parâmetros configurados, é um dos exemplos Simulation tempo e tempo total. Para a
disponíveis no módulo Renewable Energy do PSIM. Houve Control realização das simulações
apenas a substituição do bloco de controle original pelo esses tempos foram
ajustados, respectivamente,
apresentado na Seção IV, Figura 4. para 1us e 0.4s
Utilizado para a captura e a
impressão posterior do
Voltage probe
sinal elétrico, em um ponto
determinado do circuito
Utilizado para a captura e a
impressão posterior da
Current probe corrente elétrica, em um
ponto determinado do
circuito
Fonte de tensão de onda
quadrada. No sistema
simulado foi utilizado para
Square-wave
gerar valores de irradiação
voltage source
alternada e na forma
quadrada, com frequência
de 5Hz
Fonte de tensão DC. No
Fig. 7 Sistema fotovoltaico simulado. sistema simulado foi
utilizado para gerar o valor
DC voltage
fixo da temperatura e
A Tabela III apresenta os principais blocos utilizados na source
também a tensão de saída
montagem desse sistema fotovoltaico. de referência do conversor
Os blocos do tipo label nomeados de Vin e In foram Buck
inseridos para a obtenção, respectivamente, da tensão e Medição da corrente em
Current sensor um ponto determinado do
corrente de saída do módulo fotovoltaico, e aplicação de seus
circuito
valores na entrada do circuito de cada MPPT.
Diversos blocos do tipo Voltage probe foram colocados Medição da diferença de
para a obtenção dos sinais de irradiação (S), potência Voltage sensor potencial elétrico em dois
pontos distintos do circuito
máxima gerada (Pmax), tensão de saída (Vcell) e potência
real gerada (Po) pelo painel fotovoltaico. Este último foi
obtido por meio da utilização do bloco Multiplier que realiza Multiplier Multiplicador de sinais
a multiplicação dos valores de tensão e corrente de saída do
painel. A Figura 8 apresenta o sinal de irradiação simulado. Módulo fotovoltaico
As Figuras 9 e 10 apresentam, respectivamente, as potências (modelo físico), com
reais geradas (em azul) em cada um dos métodos de MPPT, entradas para valores fixos
em comparação a potência máxima gerada (em vermelho) ou variáveis de irradiação
Solar module
(S) e de temperatura (T).
pelo painel fotovoltaico. (physical mode)
Dispões ainda de terminais
de saída DC e para
verificação de sua potência
máxima gerada

Capacitor Capacitor

Inductor Indutor
Fig. 8 Sinal da irradiação simulado.
Diode Diodo

On-off switch
Chave controlada on-off
controller

Chave MOSFET com

MOSFET
acionamento pelo terminal
switch
de gate
Fig. 9 Gráfico das potências geradas pelo método P&O.

4
 [4] OS 10 países que mais usam energia solar no mundo.
Exame.com, 2012. Disponivel em:
<http://exame.abril.com.br/meio-ambiente-e-energia/
noticias/os-10-paises-que-mais-usam-energia-solar-no-
mundo#1>. Acesso em: 17 jul. 2013.
[5] ANEEL. Procedimentos de Distribuição de Energia
Elétrica no Sistema Elétrico Nacional – PRODIST.
Fig. 10 Gráfico das potências geradas pelo método da condutância [S.l.]: [s.n.], 2012.
incremental [6] DESAI, H. P.; PATEL, H. K. Maximum point algorithm
in PV generation: An overview. In: INTERNATIONAL
VI. CONCLUSÕES CONFERENCE ON POWER ELECTRONICS AND
DRIVE SYSTEMS - PEDS, 7., 2007, Bangkok.
A utilização do software PSIM, pelos menos para esse Proceedings… Bangkok: IEEE, 2007. p. 624-630.
tipo de simulação, apresenta-se como uma ótima opção. A [7] ESRAM, T,; CHAPMAN, P. L. Comparison of
forma disponível para a criação de circuitos elétricos, photovoltaic array maximum power point tracking
eletrônicos e de controle, semelhante ao Matlab/Simulink, techniques. IEEE Transactions on Energy Conversion,
realmente permite a realização dessa tarefa de forma rápida e New York, v.24, n.2, p. 439-449, 2007.
precisa. A possibilidade da escrita e inserção de código C na [8] FARANDA, R.; LEVA, S.; MAUGERI, V. MPPT
simulação ampliam em muito as funcionalidades desse techniques for PV systems: Energetic and cost
software. comparison. In: POWER AND ENERGY SOCIETY
Outra constatação a partir dos gráficos das Figuras 9 e 10, GENERAL MEETING - PESGM, 9., 2008, Pittsburgh.
obtidos com o PSIM, é que os dois métodos são satisfatórios Proceedings… Pittsburgh: IEEE, 2008. p. 1-6.
para a se encontrar a melhor transferência de potência para a [9] JAEN, C.; MOYANO, C.; SANTACRUZ, X.; POU, J.;
unidade consumidora. Apresentam resultados semelhantes, ARIAS, A. Overview of maximum power point tracking
rápidos e com bom aproveitamento de potência, desempenho control techniques used in photovoltaic systems. In:
e fácil aplicação. O rastreamento de potência é extremamente INTERNATIONAL CONFERENCE ON
eficiente, com perda de aproximadamente 3,5% [10] nos dois ELECTRONICS CIRCUITS AND SYSTEMS - ICECS,
métodos. 15., 2008, Malta. Proceedings… Malta: IEEE, 2008. p.
1099-1102.
AGRADECIMENTOS [10] BRITO, Moacyr A. G.; Junior, Luigi G.; SAMPAIO,
Leonardo P.;CANESIN, Carlos A, Avaliação das
Os autores agradecem o suporte financeiro da CAPES e Principais Técnicas para Obtenção de MPPT de Painéis
também o apoio fornecido pelas instituições de ensino Fotovoltaicos. UNESP, 2010.
envolvidas neste trabalho (IFTO, IFMG e UFU). [11] HUSSEIN, K. H.; HOSHINO, T. OSAKADA, M.
Maximum photovoltaic power tracking: an algorithm
REFERÊNCIAS for rapidly changing atmospheric conditions. IEE
Generation, Transmission and Distribution, United
[1] SALAMONI, I. T.; RÜTHER, R. O potencial brasileiro Kingdom, v. 142, n. 1, p. 59-64, 1995.
da geração solar fotovoltaica conectada à rede elétrica: [12] SEGUEL, Julio Igor López. Projeto de um sistema
análise de paridade de rede. IX Encontro Nacional e V fotovoltaico autônomo de suprimento de energia usando
Latino Americano de Conforto no Ambiente técnica MPPT e controle digital. UFMG, 2009.
Construtivo - ENCAC. Ouro Preto, MG: [s.n.]. 2007. p. [13] PSIM Simulation Software. Disponivel em:
1658 a 1667. <http://www.psim-europe.com/psim.php>. Acesso em:
[2] COORDENADOR CHIGUERU TIBA. Atlas 11 jun. 2014.
Solarimétrico do Brasil: Banco de Dados Terrestre. [14] ROSENTHAL, Leonardo C. S.; RODRIGUES, R. B.
Recife: Universitária da UFPE, 2000. Calculadora de PID. UNIUBE, 2012.
[3] ANEEL. Informações Gerais 2° Trimestre de 2013.
Agencia Nacional de Energia Elétrica - ANEEL.
Disponivel em: <http://www.aneel.gov.br/arquivos/
PDF/Z_IG_Jun_2013.pdf>. Acesso em: 17 out. 2013.

5
 ANEXO

Tabela II – Implementação em código C de cada método de MPPT

P&O Condutância incremental
static float V = 0;
static float Vant = 0;
static float P = 0;
static float deltaV = 0;
static float Pant = 0;
static float VRef = 16;
static float deltaP = 0;
static float I = 0;
static float V = 0;
static float Iant = 0;
static float Vant = 0;
static float deltaI = 0;
static float deltaV = 0;
static float delta = 0.1;
static float VRef = 16;
static float I = 0;
V = x1;
static float delta = 0.1;
I = x2;
V = x1;
deltaV = V - Vant;
I = x2;
deltaI = I - Iant;
P = V*I;
if (deltaV != 0)
{
deltaP = P - Pant;
if ((deltaI/deltaV) != (-I/V))
deltaV = V - Vant;
{
if((deltaI/deltaV) > (-I/V))
if (deltaP > 0)
{VRef = VRef + delta;}
{
else
if(deltaV > 0)
{VRef = VRef - delta;}
{VRef = VRef + delta;}
}
else
}
{VRef = VRef - delta;}
else
}
{
else
if (deltaI != 0)
{
{
if(deltaV > 0)
if(deltaI > 0)
{VRef = VRef - delta;}
{VRef = VRef + delta;}
else
else
{VRef = VRef + delta;}
{VRef = VRef - delta;}
}
}
}
Vant = V;
Pant = P;
Vant = V;
y1 = VRef;
Iant = I;
y1 = VRef;