Universidad Nacional de San Agustı́n

Escuela Profesional de Ingenierı́a Electrónica

Estrategia de gestión del lado de la

demanda en un hogar inteligente
utilizando un vehı́culo eléctrico y un
sistema hı́brido de energı́a renovable

Proyecto de Fin de Carrera

Integrantes:
Taipe Sacasqui, Leonardo – 20210285
Huertas Valverde, José Jesús

Docente:
Dr. Germán Alberto Echaiz Espinoza

Arequipa – 2024
 Resumen
En esta tesina se presenta una estrategia de gestión de energı́a orientada
al sector residencial, combinando fuentes de energı́a renovable con un vehı́culo
eléctrico, dentro del contexto de un hogar inteligente. El objetivo principal es
optimizar el uso de la energı́a disponible —ya sea solar, eólica o proveniente
de la red eléctrica— mediante un sistema hı́brido que permita reducir la de-
pendencia del suministro convencional y mejorar la eficiencia energética del
hogar.
El sistema propuesto considera una arquitectura compuesta por paneles
solares, turbina eólica, baterı́as de almacenamiento y un vehı́culo eléctrico
conectado a la red doméstica. Este vehı́culo no solo consume energı́a, sino
que también puede actuar como fuente de respaldo mediante el modo vehicle-
to-home (V2H). Para lograr una gestión eficiente, se utiliza un modelo de
programación lineal con restricciones no lineales, implementado en MATLAB,
que permite tomar decisiones inteligentes sobre cuándo cargar, descargar o
intercambiar energı́a entre los distintos componentes del sistema.
A lo largo del desarrollo se simulan diferentes escenarios de consumo y
producción, teniendo en cuenta factores como la radiación solar, el viento, y
el comportamiento tı́pico de los usuarios. Los resultados permiten visualizar
cómo la cooperación entre el vehı́culo eléctrico y las fuentes renovables puede
minimizar el consumo desde la red eléctrica y asegurar una mayor autonomı́a
energética.
Esta propuesta busca aportar una alternativa realista para enfrentar los
desafı́os actuales del consumo energético residencial, haciendo uso de tecno-
logı́as limpias, automatización y gestión inteligente de recursos.

Palabras clave: Sistemas hı́bridos renovables, Vehı́culo eléctrico, Red a vehı́cu-

lo.

1. Objetivos
a. Objetivos del tema elegido
Los objetivos derivados del estudio base sobre la gestión energética en hogares
inteligentes con vehı́culos eléctricos y sistemas hı́bridos de energı́a renovable, adap-
tados a las posibilidades de este proyecto, son:

Analizar el comportamiento energético de un hogar inteligente con integración

de un sistema hı́brido compuesto por paneles fotovoltaicos, turbina eólica,
baterı́a de almacenamiento y vehı́culo eléctrico.

Comprender los flujos de energı́a entre los distintos componentes del sistema
(energı́a solar, eólica, baterı́a, red eléctrica y vehı́culo eléctrico) en los modos
home-to-vehicle, vehicle-to-home y grid-to-vehicle.

Estudiar la estrategia general de gestión de energı́a basada en programación

lineal y restricciones no lineales.

Identificar los momentos de mayor consumo energético residencial y cómo el

sistema puede responder ante ellos mediante la optimización del uso de energı́as
renovables.

1
b. Objetivos que se aumentarán al tema elegido
Como parte del enfoque particular de esta tesina, se añaden los siguientes obje-
tivos complementarios, pensados para ser alcanzables dentro del tiempo y recursos
disponibles:
Sustentar teóricamente las principales ecuaciones utilizadas en el modelo energéti-
co, detallando su origen, variables involucradas y fundamentos fı́sicos o ma-
temáticos.
Organizar y presentar una interpretación accesible de la arquitectura propues-
ta en el artı́culo base, explicando sus componentes y conexiones desde una
perspectiva académica clara.
Diseñar un esquema conceptual en Simulink que represente gráficamente el
flujo energético entre los elementos principales del sistema, permitiendo vi-
sualizar las interacciones entre las fuentes renovables, el vehı́culo eléctrico, el
sistema de almacenamiento y la red eléctrica.

2. Resumen de avance
Durante la primera etapa del desarrollo de esta tesina se han llevado a cabo
una serie de actividades fundamentales que sientan las bases teóricas y prácticas del
trabajo. A continuación, se detallan los avances realizados:
Se realizó un análisis profundo de las ecuaciones presentadas en el artı́culo
base, enfocándose en comprender su origen y aplicación dentro del sistema
hı́brido de gestión energética.
Se investigaron y contrastaron diversas fuentes académicas confiables para
validar cada una de las fórmulas utilizadas en el modelo, lo que permitió ase-
gurar su veracidad y consistencia con fundamentos de electrónica de potencia
y energı́as renovables.
A partir de fórmulas más simples de generación y conversión de energı́a, se
reconstruyó el camino lógico-matemático hasta llegar a las ecuaciones emplea-
das en el artı́culo original. Este desarrollo incluyó el estudio de modelos de
paneles solares (como el modelo de dos diodos), turbinas eólicas y baterı́as.
Se explicó el propósito de cada ecuación dentro del sistema propuesto, permi-
tiendo entender cómo cada componente (paneles solares, turbina eólica, baterı́a
y vehı́culo eléctrico) aporta al flujo energético del hogar inteligente.
Se ha elaborado una pequeña simulación preliminar en Simulink, en la cual
se representó el comportamiento de una fuente fotovoltaica conectada a una
carga doméstica simple, como primer acercamiento al modelado del sistema
completo. Esta simulación servirá como base para futuras ampliaciones del
esquema hı́brido completo.
Estos avances permiten contar con una comprensión sólida de la base teórica del
sistema y sentar las bases para una implementación progresiva del modelo energético
propuesto.

2
 Ecuación 1: Modelo equivalente de un panel
fotovoltaico (dos diodos)
Esta ecuación representa el comportamiento fı́sico de una celda fotovoltaica realista.
Se basa en un modelo de dos diodos que describe con gran fidelidad la corriente
generada en función del voltaje aplicado, considerando recombinaciones internas y
la fotocorriente.

! qV
" ! qV
"
I = Iph → I01 e n1 Ns kB T → 1 → I02 e n2 Ns kB T → 1 (1)
Donde:

Iph : Corriente generada por la luz (fotocorriente),

I01 , I02 : Corrientes de saturación de los diodos 1 y 2,

n1 , n2 : Factores de idealidad de cada diodo,

Ns : Número de celdas en serie,

T : Temperatura absoluta (K),

kB : Constante de Boltzmann,

q: Carga del electrón.

Para facilitar el análisis, se simplifica a un solo diodo equivalente:

! qV "
I = Iph → I0 e nNs kB T → 1 (2)
La fotocorriente Iph se modela en función de la irradiancia y la temperatura:
G
Iph = [Iscn + µi (T → Tref )] (3)
Gref
La corriente en el punto de máxima potencia (Vpv→mpp ) se calcula como:
# qV $
pv→mpp
Impp = Iph → I0 e nNs Nss kB T
→1 (4)

Sustituyendo Iph en Impp :

# qV $
G pv→mpp
Impp = [Iscn + µi (T → Tref )] → I0 e nNs Nss kB T
→1 (5)
Gref
La potencia generada se obtiene multiplicando por Vpv→mpp , número de ramas en
paralelo y eficiencia del convertidor:

Ppv = Np · Impp · Vpv→mpp · ωDC/DC (6)

Sustituyendo Impp :

3
 #
G
Ppv = Np · Vpv→mpp · ωDC/DC · [Iscn + µi (T → Tref )]
Gref
# qV $$
pv→mpp
→I0 e nNs Nss kB T
→1 (7)

Considerando I0 como función del voltaje de circuito abierto Vocn :

I + µi (T → Tref )
I0 = !scn " (8)
q(Vocn +µv (T →Tref ))
exp nNs kB T
→1

Sustituyendo esta expresión de I0 en la fórmula de potencia:

% # $
G
Ppv = Np [Iscn + µi (T → Tref )]
Gref

# # $ $
Np [Iscn + µi (T → Tref )] qVpv→mpp
→ ! " · exp →1
exp
q(Vocn +µv (T →Tref ))
→ 1 nNss Ns kB T 
nNs kB T
nNs kB T Nss
· · Vpv→mpp · ωDC/DC (9)
q

Resultado final:

 
 # $ # # $ $
G Np [Iscn + µi (T → Tref )] qVpv→mpp nN
Ppv = Np [Iscn + µi (T → Tref )] → ! " · exp →1 ·
 Gref exp
q(Vocn +µv (T →Tref ))
→ 1 nNss Ns kB T 
nNs kB T
(10)

Importancia: Esta expresión permite determinar con precisión la potencia eléctrica

entregada por el campo fotovoltaico. Es esencial para el control de flujos energéticos
en el hogar inteligente de la tesina, ya que permite decidir cuándo cargar baterı́as,
alimentar cargas o inyectar energı́a a la red, asi entonces, la potencia generada por
un panel solar se define como el producto de la corriente y el voltaje generados,
como se muestra en la Ecuación (1) [1, 2].
—

Ecuación 2: Cálculo del voltaje en el punto de

máxima potencia Vpv→mpp mediante la función de
Lambert W
Esta ecuación permite obtener una expresión cerrada del voltaje Vpv→mpp , punto
crucial para el cálculo de la potencia máxima que puede entregar un sistema foto-
voltaico. Su derivación se basa en el modelo de un solo diodo y en la aplicación de
la función de Lambert W , utilizada para resolver ecuaciones del tipo xex = A.

4
 Vpv→mpp : Voltaje en el punto de máxima potencia (V).

Impp : Corriente en el punto de máxima potencia (A).

Iph : Fotocorriente generada por irradiancia solar (A).

I0 : Corriente de saturación inversa del diodo (A).

Vocn : Voltaje de circuito abierto nominal (V).

Iscn : Corriente de cortocircuito nominal (A).

G: Irradiancia solar actual (W/m²).

Gref : Irradiancia de referencia (W/m²), tı́picamente 1000 W/m2 .

T : Temperatura del módulo (°C o K).

Tref : Temperatura de referencia (25°C o 298.15 K).

q: Carga del electrón, 1,602 ↑ 10→19 C.

kB : Constante de Boltzmann, 1,38 ↑ 10→23 J/K.

n: Factor de idealidad del diodo (sin unidad).

Ns : Número de celdas en serie.

Np : Número de celdas en paralelo.

Nss : Número de submódulos (strings) en serie.

ki : Coeficiente de temperatura de la corriente (A/°C).

kv : Coeficiente de temperatura del voltaje (V/°C).

vt : Tensión térmica equivalente, dada por:

nNs kB T
vt =
q
! qV "
I = Iph → I0 e nNs kB T → 1 (11)
En el punto de máxima potencia:

I = Impp , V = Vpv→mpp

La fotocorriente se define como:

G
Iph = [Iscn + ki (T → Tref )] (12)
Gref
Iscn : corriente de cortocircuito en condiciones estándar,

ki : coeficiente de temperatura de la corriente (A/°C).

5
Sustituyendo Iph en la ecuación de corriente:
# qV $
G pv→mpp
Impp = [Iscn + ki (T → Tref )] → I0 e nNs kB T
→1 (13)
Gref
Despejamos para obtener la forma:
qVpv→mpp
G
I0 · e nNs kB T
= [Iscn + ki (T → Tref )] → Impp + I0
Gref
Definimos la tensión térmica:
nNs kB T
vt =
q
Reorganizamos la expresión:
# $
Vpv→mpp
1 G
e vt = [Iscn + ki (T → Tref )] → Impp + I0
I0 Gref
Luego:
# $
Vpv→mpp Vpv→mpp vt G
· e vt = [Iscn + ki (T → Tref )] → Impp + I0
vt I0 Gref
Aplicando la función inversa:
# # $$
Vpv→mpp vt G
=W [Iscn + ki (T → Tref )] → Impp + I0
vt I0 Gref
Finalmente:

# # $ $
1 G
Vpv→mpp = vt · W [Iscn + ki (T → Tref )] → Impp + I0 · vt (14)
I0 Gref

Si en lugar de I0 usamos su aproximación experimental basada en Vocn :

I + ki (T → Tref )
I0 = !scn "
q(Vocn +kv (T →Tref ))
exp nNs kB T
→1
y reagrupando los factores, se obtiene una expresión computacional más directa:

# # # $ ,
G (Vocn + kv (T → Tref ))Nss
Vpv→mpp = LambertW 0, [Iscn + ki (T → Tref )] Np · exp · exp(1) [(Iscn
Gref vt
(15)

Importancia: Esta fórmula permite calcular Vpv→mpp , valor indispensable para esti-
mar la máxima potencia eléctrica entregada por los paneles. Es esencial en la gestión
energética ya que permite determinar con precisión el punto óptimo de operación
del sistema fotovoltaico bajo condiciones variables de temperatura e irradiancia, por
lo tanto, la estimación del punto de máxima potencia requiere conocer tanto la co-
rriente como el voltaje en ese punto, lo cual se puede calcular considerando factores
como temperatura e irradiancia, según el modelo de dos diodos [2, 3].

6
Ecuaciones 3, 4 y 5: Modelado de potencia eólica
en el sistema hı́brido

Estas ecuaciones permiten estimar la energı́a que se puede obtener del viento a
través de una turbina eólica. Se consideran tanto principios fı́sicos como correcciones
empı́ricas y pérdidas eléctricas, fundamentales para una estimación realista en el
sistema hı́brido propuesto en la tesina.

Sı́mbolo Descripción
Ecin Energı́a cinética del viento
m Masa del aire que atraviesa el rotor
ε Densidad del aire [kg/m3 ]
S Área barrida por las palas [m2 ]
R Radio de las palas [m]
Vwind Velocidad del viento [m/s]
Pwind,raw Potencia cinética total del viento
Cp (ϑ, ϖ) Coeficiente de potencia
ϖ Ángulo de pitch de las palas [°]
ϑ Tip Speed Ratio (TSR)
ϑi TSR corregido
! Velocidad angular del rotor [rad/s]
ωAC/DC Eficiencia del convertidor AC/DC
ωDC/DC Eficiencia del convertidor DC/DC

Paso 1. Energı́a cinética del viento

1 2
Ecin = mVwind
2

m = ε · V = ε · S · Vwind · t

1 3
↓ Ecin = · ε · S · Vwind ·t
2

Paso 2. Potencia cinética instantánea

dEcin 1 3
Pwind,raw = = · ε · S · Vwind
dt 2

S = ϱR2

1
↓ Pwind,raw = · ε · ϱR2 · Vwind
3
2

7
Paso 3. Aplicación del coeficiente de potencia
1
Protor = Cp (ϑ, ϖ) · Pwind,raw = · ε · ϱR2 · Cp (ϑ, ϖ) · Vwind
3
2

Paso 4. Consideración de pérdidas eléctricas

PWT = Protor · ωAC/DC · ωDC/DC

Ecuación (3): Potencia útil del aerogenerador

1
PWT = · ε · ϱR2 · Cp (ϑ, ϖ) · Vwind
3
· ωAC/DC · ωDC/DC
2

Paso 5. Relación de velocidad de punta (TSR)

vtip = R · !

vtip R·!
ϑ= =
Vwind Vwind
Ecuación (5): Relación de velocidad de punta

R·!
ϑ=
Vwind

Paso 6. Corrección empı́rica de ϑ

Ecuación (4): TSR corregido

1 1 0,035
= → 3
ϑi ϑ + 0,08ϖ ϖ + 1

Paso 7. Modelo del coeficiente de potencia

# $
116
Cp (ϑ, ϖ) = 0,517 + 0,08ϖ → 4,06 · e→21/ωi + 0,0068ϑ
ϑi

Paso 8. Sustitución total en PWT

- # $ .
1 2 116 →21/ωi 3
PWT = ·ε·ϱR · 0,517 + 0,08ϖ → 4,06 e + 0,0068ϑ ·Vwind ·ωAC/DC ·ωDC/DC
2 ϑi

Importancia: Estas ecuaciones permiten determinar la producción eléctrica del

sistema eólico en función de las condiciones del viento. En el contexto de la tesina,

8
son fundamentales para estimar el aporte renovable en la arquitectura hı́brida del
hogar inteligente, especialmente en momentos de baja irradiancia solar.
(Referencia: Serway y Jewett, Cap. 10: Rotation of a Rigid Object About a Fixed Axis))

Ecuaciones 6, 7 y 8: Modelado de la energı́a en la

baterı́a del sistema hı́brido

Estas ecuaciones permiten modelar la dinámica de la energı́a almacenada en la

baterı́a, considerando el flujo de potencia, la eficiencia del convertidor y las pérdidas
durante carga y descarga. Son clave para simular y controlar la interacción energética
entre el sistema fotovoltaico, el hogar y el vehı́culo eléctrico.

Ecuación (6)

dE(t)
P (t) =
dt
dEbatt (t)
Pbatt (t) =
dt
dEbatt
P =
dt
Resultado final - Ecuación (6)

dEbatt
=P
dt

Ecuación (7)

Ebatt (t) → Ebatt (t → 1)

P ↔
”t

”Ebatt = Ebatt (t) → Ebatt (t → 1)

Resultado final - Ecuación (7)

”Ebatt = Ebatt (t) → Ebatt (t → 1)

Ecuación (8)

9
Parámetros involucrados:

P: Potencia en el enlace baterı́a-bus (W)

ωDC/DC : Eficiencia del convertidor DC/DC bidireccional
ωch : Eficiencia durante la carga de la baterı́a
ωdis : Eficiencia durante la descarga de la baterı́a
”Ebatt : Variación de energı́a en la baterı́a (Wh)

Cálculo base con paso de tiempo unitario:

”Ebatt = P · ”t

”t = 1 h

”Ebatt = P

”Ebatt = P · ωch · ωDC/DC , si P > 0

P
”Ebatt = , si P < 0
ωDC/DC · ωdis

”Ebatt = 0, si P = 0

Considerando eficiencia de los procesos:



 P · ωch · ωDC/DC si P > 0 (carga)

 P
”Ebatt = si P < 0 (descarga)

 ωDC/DC · ωdis

0 si P = 0 (sin flujo)

Resultado final - Ecuación (8)



 P · ωch · ωDC/DC si P > 0

 P
”Ebatt = si P < 0

 ωDC/DC · ωdis

0 si P = 0

10
 Situación Energı́a útil en baterı́a Razón
P >0 P · ωDC/DC · ωch Entra energı́a, pero con pérdidas
P
P <0 Se necesita más energı́a para compensar pérdidas
ωDC/DC · ωdis
P =0 0 No hay flujo de energı́a

Importancia: Estas ecuaciones permiten representar el comportamiento dinámico

y realista del almacenamiento en baterı́a dentro del sistema hı́brido. Son claves para
la gestión eficiente de la energı́a, y para simular cómo la baterı́a se carga desde
fuentes renovables o descarga hacia el hogar o vehı́culo eléctrico.
(Referencia: Serway y Jewett, Cap. 8, Secc. 8.5)

Referencias
[1] G. M. Masters, Renewable and E!cient Electric Power Systems, 2nd ed. John
Wiley & Sons, 2013.

[2] M. G. Villalva, J. R. Gazoli, and E. R. Filho, “Comprehensive approach to

modeling and simulation of photovoltaic arrays,” IEEE Transactions on Power
Electronics, vol. 24, no. 5, pp. 1198–1208, 2009.

[3] G. Walker, “Evaluating mppt converter topologies using a matlab pv model,”

Journal of Electrical & Electronics Engineering, Australia, vol. 21, no. 1, pp.
49–56, 2001.

11