See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/336731902

MODELO SIMPLIFICADO PARA UN MÓDULO FOTOVOLTAICO USANDO MATLAB

Conference Paper · October 2019

CITATIONS READS
0 1,906

2 authors:

Jose Alexis Toledo Hector Alfredo Aguilar-Burgoin

Autonomous University of Baja California Sur Autonomous University of Queretaro
6 PUBLICATIONS 0 CITATIONS 5 PUBLICATIONS 0 CITATIONS

SEE PROFILE SEE PROFILE

All content following this page was uploaded by Jose Alexis Toledo on 23 October 2019.

The user has requested enhancement of the downloaded file.

 MODELO SIMPLIFICADO PARA UN MÓDULO FOTOVOLTAICO USANDO MATLAB

José A. Toledo Cornejo | Alfredo Flores Irigollen1 | Héctor A. Aguilar Burgoin | Jesús A. Martínez Guerra
Departamento Académico de Ingeniería en Pesquerías, Universidad Autónoma de Baja California Sur, Carretera al Sur km 5.5, La Paz,
Baja California Sur, 23080, México,
+52 (612) 1238800, [email protected]

RESUMEN
En este artículo se presenta un modelo que describe el comportamiento dinámico de un panel fotovoltaico cuya característica esencial es
que no requiere de procesos iterativos para la estimación de parametros como las resistencias en serie y paralelo y tampoco de los
coeficientes de idealidad de los diodos. El modelo permite caracterizar las gráficas “corriente – voltaje” y “potencia – voltaje” del panel
fotovoltaico. Dicho modelo es flexible y además permite realizar simulaciones para irradiación y temperatura variables.
El presente modelo fue validado mediante mediciones realizadas en un panel fotovoltaico de la marca ASTROPOWER AP-1106, con
resultados que se aproximan favorablemente a los datos experimentales con un coeficiente de determinación superior al 90%.
ABSTRACT
This article presents a model that describes the dynamic behavior of a photovoltaic panel whose essential characteristic is that it does not
require iterative processes for the specification of parameters such as series and parallel resistors and neither of the ideality coefficients of
the diodes. The model allows characterizing the “current - voltage” and “power - voltage” graphs of the photovoltaic panel. This model is
flexible and also allows simulations for variable irradiation and temperature.
The present model was validated by measurements made on a photovoltaic panel of the ASTROPOWER AP-1106 brand, with results
that favorably approximate the experimental data with a determination coefficient greater than 90%.
Palabras clave: Simulación, simulink, arduino, labVIEW, panel fotovoltaico.

INTRODUCCIÓN
Las energías renovables a finales del año 2018 contaron con una potencia mundial instalada de 2378 GW. Con lo cual, representaron
cerca del 27% de la generación de electricidad en todo el mundo. De dicha potencia instalada, 505 GW fueron pertenecientes a la energía
solar fotovoltaica [1]. Esta fuente de energía renovable es catalogada como no contaminante e inagotable [2], ya que convierte los fotones
provenientes de los rayos del sol en energía eléctrica. A nivel nacional para este mismo año, la producción total de energía eléctrica en
México fue de 317,278 GWh. Donde, el 79% fue a partir de combustibles fósiles, 4.3% por medio de energía nuclear y tan solo un solo
0.7 % por medio de la energía solar, el porcentaje restante fue distribuido entre las demás fuentes de energía renovable [3].

A pesar del crecimiento exponencial que la energía solar ha ido teniendo a lo largo de los años, aún posee puntos débiles e importantes
para su investigación y desarrollo [4]. Por ejemplo, en temas de eficiencia, la cual suele ser baja al no aprovechar al máximo la energía
incidente[5]. Además, la eficiencia no solo esta determinada por la cantidad de radiación solar incidente, sino que también por la temperatura
de trabajo del módulo, así pues, mientras la temperatura en el panel sea mayor, menor será su rendimiento [4]. Por tanto, es necesario el uso
de herramientas que permitan analizar y simular el comportamiento eléctrico de los paneles solares bajo distintitas condiciones de
temperatura y radiación [6].

Dicho lo anterior este articulo aborda la simulación de un panel solar fotovoltaico utilizando las ecuaciones propuestas por Gabriel Villa y
Carlos Robles que permite simular de una manera muy acertada las curvas características de corriente-voltaje.[6] Para ello se utilizó el
modelo matemático propuesto por ellos dos utilizando las herramientas pertinentes para realizar la captura de los datos y el análisis de ellos.
Para la obtención de los datos fue necesario utilizar hardware especializado para la captura de las variables necesarias, dichas variables
alimentan el modelo anteriormente mencionado y se comparan las curvas obtenidas por medio del modelo con las curvas que los fabricantes
de los paneles brindan en sus fichas de especificaciones. El panel utilizado fue un ASTROPOWER AP-1106, un panel monocristalino de
110 W, el cual en su ficha técnica muestra sus curvas características de corriente-voltaje para diferentes valores de temperatura y radiación,
tratando de aproximarse a la realidad, fuera de las condiciones ideales de pruebas de laboratorio.

Para facilitar los numerosos cálculos que se requieren se utilizó el software especializado “MATLAB” el cual permite realizar una gran
cantidad de cálculos en muy poco tiempo, además permite visualizar de manera muy amigable, por medio de gráficas, los datos obtenidos
por medio de las mediciones. En dicho software se escribieron las ecuaciones necesarias para realizar los cálculos del modelo matemático
y la comparación de los datos reales obtenidos por medio de mediciones.
ANTECEDENTES
La literatura en torno al análisis, modelación y simulación de sistemas fotovoltaicos es muy amplia, a saber: Gradella et al. [7] presentan
un método para modelar y simular arreglos fotovoltaicos. Su método consiste en encontrar los parámetros de la ecuación no lineal I-V
mediante ajuste de la curva a tres puntos: voltaje a circuito abierto, potencia máxima, corriente a circuito cerrado. Para realizar esto, se
basan en el uso de un modelo fotovoltaico de simple diodo que incluye los efectos de las resistencias en serie y paralelo, y garantiza que la
potencia máxima del modelo coincide con la potencia máxima del arreglo. Validan su modelo con datos experimentales de arreglos
comerciales.
Por otro lado, diversos autores han empleado Simulink, un ambiente de programación en bloque de Matlab, para modelar y simular el
comportamiento dinámico de módulos fotovoltaicos. Así, por ejemplo, Selmi y Belghouthi [8] presentaron un modelo para una celda solar
de unión triple InGaN, para lo cual utilizaron ecuaciones representadas en bloque en el ambiente de Simulink. Así mismo, los autores Resk
y Hasaneen [9] analizaron el modelo de una celda solar de unión triple pero lo extendieron a un panel. Su modelo incluye un algoritmo
MPPT que utiliza redes neuronales artificiales. Es también muy destacable el trabajo de Acevedo-Luna y Morales-Acevedo [10] , quienes
crearon un algoritmo que permite obtener los parámetros de los modelos de diodo simple y doble diodo de una celda solar con base en los
datos experimentales de las curvas I-V. Los turcos Setiuawan y Siregar [11] validaron el modelo de bloques Simulink de un sistema solar
fotovoltaico y evaluaron diferentes configuraciones del inversor para clima tropical. Finalmente, Kapoor et al. [12] crearon un control
adaptativo para estudiar el comportamiento dinámico de cualquier carga en un panel solar.

METODOLOGÍA
De manera general, el primer paso fue conocer el voltaje máximo del panel fotovoltaico, con el fin de calcular las resistencias a utilizar y
así hacer un divisor resistivo, el cual baja el voltaje al valor que se necesite. Es decir, nuestro panel tiene un valor máximo de 16.7 voltios,
mientras que Arduino Uno, nuestro sistema para hacer el registro automático de los datos, sólo puede leer señales de hasta 5 V, por tanto,
es necesario reducir el voltaje máximo del panel y así no generar daños en el Arduino. La Figura 1 presenta el diagrama esquemático de la
conexión entre el panel (representado como fuente de poder), el divisor resistivo y el Arduino Uno.

Figura 1. Diagrama esquemático entre el panel solar, el divisor resistivo y la tarjeta Arduino. Elaboración propia.

Por otro lado, también se revisó en la tabla de especificaciones del panel el valor de la corriente máxima del mismo, y así poder elegir el
sensor de corriente del rango adecuado. El valor máximo de la corriente para este trabajo fue de 6.6 A, por lo tanto, el sensor que se utilizó
fue el ACS712, capaz de soportar hasta 10 A en corriente directa. La Figura 2 muestra el diagrama esquemático de la conexión serie del
panel con una resistencia, la cual es necesaria para representar la carga conectada al panel, subsecuentemente este subsistema se conecta al
sensor de amperaje, el cual a su vez se acopló a la tarjeta Arduino.
 Figura 2. Diagrama de conexión entre el panel solar, la carga, el sensor de corriente y la tarjeta Arduino. Elaboración propia.

Así mismo, para realizar las mediciones de temperatura ambiente, se montó el sensor LM35 con el Arduino Uno tal como se aprecia en
la Figura 3. Este sensor es capaz de soportar temperaturas entre el rango que va desde los -55ºC a los 150ºC.

Figura 3. Diagrama de conexión entre el sensor de temperatura ambiente y la tarjeta Arduino. Elaboración propia.

Finalmente, el circuito completo que permitió la lectura de los datos se muestra en la Figura 4, es importante destacar, que la tarjeta
Arduino se enlazó con LabVIEW para realizar la adquisición y registro de datos.

Figura 4. Diagrama completo de las conexiones a la tarjeta Arduino.

MODELO MATEMATICO DEL PANEL FOTOVOLTAICO
Lo esencial de este modelo matemático es que es mucho más sencillo de manipular que otros modelos, debido a que solo es necesario
encontrar un parámetro de ajuste de curva que puede ser obtenido directamente de la ecuación I-V [6] .De modo que los demás parámetros
necesarios para el modelo son sencillamente obtenidos de la hoja de especificaciones del panel fotovoltaico deseado. Como se verá en el
modelo, estas ecuaciones podrán representar al panel FV para valores dinámicos de irradiación (sobre plano inclinado u horizontal) y
temperatura, sin necesidad de hacer correcciones por el área de la celda, el número de ellas por panel o correcciones por conexiones serie-
paralelo de las células FV internas del panel. La ecuación propuesta para representar la relación I-V del panel FV es la siguiente [6]:
V 1
Ix ( − )
I(V) = (
−1
)
[1 − e b Vx b ] (1)
1−e b

Donde:
I(V); es la corriente en función del voltaje, 𝑏 es el parámetro de ajuste, que como se mencionó con anterioridad es el único parámetro
que deberá ser calculado, Ix es la fotocorriente generada para los valores dinámicos de irradiancia y temperatura, dichos valores también
afectan a Vx que corresponden al voltaje de circuito abierto.

La ecuación que nos describe el comportamiento del voltaje es la siguiente [6]:

𝐺 𝑉 −𝑉
𝐺 [ 𝑙𝑛( 𝑚𝑎𝑥 𝑜𝑐 )]
𝐺𝑟𝑒𝑓 𝑉𝑚𝑎𝑥 −𝑉𝑚𝑖𝑛
𝑉𝑥 = 𝑁𝑠 𝐾𝑣(𝑇 − 𝑇𝑛 ) + 𝑁𝑠 𝑉𝑚𝑎𝑥 − 𝑁𝑠 (𝑉𝑚𝑎𝑥 − 𝑉𝑚𝑖𝑛 )𝑒 (2)
𝐺𝑟𝑒𝑓
Donde:
𝑁𝑠 es el número deseado de paneles en serie, 𝐺𝑟𝑒𝑓 es una constante de irradiancia solar para cualquier panel, toma el valor de 1000
W/m2, 𝐺 es nuestro vector de irradiancia, 𝐾𝑣 es el coeficiente de voltaje de circuito abierto, este nos indica el decaimiento de voltaje por
cada aumento de ° C de temperatura en el panel (deberá ser un valor negativo), 𝑇 es nuestro de vector de temperaturas, 𝑇𝑛 es la
temperatura nominal de pruebas, toma el valor de 25 °C, 𝑉𝑜𝑐 es el voltaje de circuito abierto, específico para cada panel, 𝑉𝑚𝑎𝑥 es el voltaje
de circuito abierto máximo, el cual corresponde al 103% del Voc, asimismo, 𝑉𝑚𝑖𝑛 , es el voltaje de circuito abierto mínimo, el cual
corresponde al 85% del Voc.

Al mismo tiempo, la ecuación que nos describe el comportamiento de la fotocorriente es la siguiente [6]:

G
Ix = Np [Isc + Ki(T − Tn)] (3)
Gref

Donde Np es el número deseado de paneles en paralelo, Isc es la corriente de cortocircuito, Ki es el coeficiente de temperatura de
cortocircuito, este nos indica el aumento de amperaje por cada aumento de ° C de temperatura en el panel (las unidades deberán ser en
A/°C).

Como se ha mencionado, la bondad de este modelo matemático es que solo es necesario calcular el parámetro “b” adimensional, el cual
fluctúa entre el rango de 0.01 a 0.18. Así, la aproximación de la ecuación (4) calcula este parámetro necesario [6].

1
(− )
1−𝑒 𝑏 ≈ 1. (4)
RESULTADOS Y DISCUSIÓN

La Figura 5 muestra el comportamiento I-V del panel para tres condiciones de radiación y temperaturas distintas. Las líneas continuas
representan a los resultados teóricos, mientras que los puntos corresponden a las mediciones experimentales. Se aprecian tres grupos de
colores diferentes, cada grupo es la comparación entre el valor teórico y el valor experimental.

Figura 5. Comparación entre los datos resultados y experimentales.

Para la simulación con una radiación de 1000 W/m2 y 25°C de temperatura ambiente, se obtuvo una correlación de 0.99. Para las
condiciones de 1000 W/m2 con 60°C ambiente, la correlación resultó ser de 0.90. Finalmente, para las condiciones de 800 W/m2 y 45°C de
temperatura ambiente la correlación fue de 0.99.

La Figura 6 muestra las curvas teóricas P-V correspondientes a las tres condiciones de radiación. Para la simulación con una radiación
de 1000 W/m2 y 25°C de temperatura ambiente, se obtuvo una potencia máxima de 110 W. Para las condiciones de 1000 W/m2 con 60°C
ambiente, la potencia máxima resultó ser de 96 W. Finalmente, para las condiciones de 800 W/m2 y 45°C de temperatura ambiente la
potencia pico fue de 82 W.

Figura 6. Potencia teórica bajo tres condiciones de radiación y temperatura.

 CONCLUSIONES

• La implementación en Simulink del modelo matemático del panel fotovoltaico permitió estimar su comportamiento con una muy
buena precisión.
• La principal ventaja de este modelo radica en que no es necesario conocer los factores de idealidad de los diodos del panel, lo que
en muchas ocasiones es difícil de conocer ya que estos no se encuentran en las tablas de especificaciones de los paneles, y para
calcularlos es necesario el uso de ciclos iterativos.
• Se puede apreciar que la temperatura es un factor que afecta drásticamente la potencia pico del panel fotovoltaico. El modelo
manifiesta claramente este efecto.
• El valor del parámetro de ajuste “b” resultó ser de 0.091, el cual está dentro del rango especificado por la literatura.
• La bondad del modelo quedo demostrada al reproducir las mediciones experimentales con un coeficiente de determinación
superior al 90 %.

REFERENCIAS

[1] REN 21, «Renewables 2019 Global Status Report». Paris, p. 336, 2019.
[2] ACCIONA, «Energía solar fotovoltaica», 2015. .
[3] CENACE, «Programa de ampliación y modernización de la red nacional de transmisión y redes generales de
distribución del mercado eléctrico mayorista 2019-2033», México, 2019.
[4] SEBASTIÁN ANDRÉS MAGGI SILVA, «MODELAMIENTO TÉRMICO DE UN PANEL FOTOVOLTAICO CON
DISIPADOR DE CALOR OPERANDO EN EL NORTE DE CHILE.», UNIVERSIDAD DE CHILE, Chile, 2013.
[5] Anderson Guillermo Vera-Dávila, Jhan Carlos Delgado-Ariza, y Sergio Basilio Sepúlveda-Mora, «Validación del
modelo matemático de un panel solar empleando la herramienta Simulink de Matlab», 2, vol. 8, 2018.
[6] Carlos Robles y Gabriel Villa, «CONTROL DEL PUNTO DE MÁXIMA POTENCIA DE UN PANEL SOLAR
FOTOVOLTAICO, UTILIZANDO LÓGICA DIFUSA», Revista electrónica de estudios telemáticos, vol. 10, n.o 2,
2011.
[7] Gradella Villalva M., Jonas Rafael Gazoli, y Ernesto Ruppert Filho, «Comprehensive Approach to Modeling and
Simulation of Photovoltaic Arrays», vol. 24, n. o 5, pp. 1198-1208, 2009.
[8] T. Selmi y R. Belghouthi, «A novel widespread Matlab/Simulink based modeling of InGaN double hetero-junction p-i-
n solar cell», International Journal of Energy and Environmental Engineering, vol. 4, n.o 8, pp. 273-281, 2017.
[9] H. Rezk y E.S. Hasaneen, «A new MATLAB/Simulink model of triple-junction solar cell and MPPT based on artificial
neural networks for photovoltaic energy systems», A new MATLAB/Simulink model of triple-junction solar cell and
MPPT based on artificial neural networks for photovoltaic energy systems, vol. 3, n.o 6, pp. 873-881, 2015.
[10] A. Acevedo-Luna y A. Morales-Acevedo, «Study of validity of the single-diode model for solar cells by I-V curves
parameters extraction using a simple numerical method.», Journal of Materials Science: Materials in Electronics, pp.
1-7.
[11] E. A. Setiawan, A. Setiawan, y D. Siregar, «Analysis on solar panel performance and PV-inverter configuration for
tropical region», Journal of Thermal Engineering, vol. 3, n.o 3, pp. 1259-1270, 2017.
[12] D. Kapoor, P. Sodhi, y D. Deb, «Solar panel simulation using adaptive control», In International Conference on
Control Applications, pp. 1124-1130, 2012.

View publication stats