Celada Padilla James Skinner
Celada Padilla James Skinner
Celada Padilla James Skinner
UNIVERSIDAD NACIONAL
“PEDRO RUIZ GALLO”
ESCUELA DE POSGRADO
TESIS
Presentada para optar el Grado Académico de Doctor en
Ciencias de la Ingeniería Mecánica y Eléctrica con
mención en Energía
AUTOR:
Celada Padilla, James Skinner
ASESOR:
Dr. Aguinaga Paz, Amado
Lambayeque- Perú
2023
ii
Autor
M.Sc. James Skinner Celada Padilla.
______________________________
Presidente
Dr. Luis Jaime Collantes Santisteban
_______________________
Secretario
Dr. Sebastian Huangal Scheineder
Vocal
Dr. Luis Antonio Pozo Suclupe
iii
AGRADECIMIENTO
DEDICATORIA
INDICE
INDICE DE TABLAS
INDICE DE FIGURAS
SIMBOLOGIA
Descripción Símbolo
RESUMEN
ABSTRACT
The method used to determine the optimal energy flows was based on the
variation of the operating parameters of the internal combustion engine and the
electric motor of a hybrid vehicle. The information on the operating parameters of
these engines was obtained from the measurements that were made to the test
vehicles with similar characteristics, which were compared with the data provided
by the manufacturers of said hybrid vehicles. It was possible to determine the
linear and non-linear relationships of the operating parameters of the internal
combustion engine and the electric motor, at different operating conditions of the
hybrid vehicle, depending on the speed of the vehicle, as well as establishing the
values in which specific fuel consumption is optimized.
The results show that, in all the operating conditions of the hybrid vehicle, the
degree of hybridization, which measures the operating times independently or
jointly of the internal combustion engine and the electric motor, decreases in all
cases as the time increases. increases the speed of the vehicle. Likewise, the
numerical modeling of the mechanical energy contribution of the internal
combustion engine and the electric motor for the drive of the hybrid vehicle was
made, this modeling was done based on the torques of said engines.
INTRODUCCIÓN
En la Unión Europea, en el trimestre 3 del año 2021, alrededor del 21% de las
ventas de vehículos nuevos fueron híbridos, es decir aquellos que no requieren
recargarse, y que solo reutilizan la energía al momento del frenado; un 9% de
las ventas correspondió a los vehículos híbridos recargables y un 10% a
vehículos totalmente eléctricos. En cuanto a la venta de vehículos nuevos que
solo utilizan motor de combustión interna, 39.5% fueron las ventas de vehículos
con motor de nafta y 17.6% de motor diésel. (Rueda. 2021).
Los países miembros de la Unión Europea, han aprobado en el año 2023, una
medida que prohíbe la venta de vehículo nuevos con motor de gasolina y/o
Diesel; con ello se busca un equilibrio perfecto entre los trabajadores, el medio
ambiente y los intereses de la Industria Automotriz; teniendo la mira del uso de
16
La china, el gigante del Asia, evalúa que, a partir del año 2035, los vehículos
nuevos que se vendan dentro del país, tienen que estar electrificados; el
Ministerio de Industria y Tecnología de la Información de China, considera que
para dicho año, el 50% de sus vehículos nuevos sean totalmente eléctricos, y el
otro 50% sean de vehículos híbridos recargables. La eliminación de la venta de
vehículo con motor térmico, será de manera paulatina, se estima que para el
2030, el 75% de vehículos será híbridos. En el año 2023, la China ha
experimentado un crecimiento explosivo en la venta de vehículos híbridos y
eléctricos, y ello aún tiene expectativas d mayor crecimiento, a medida que se
reestablezca la confianza entre consumidores y fabricantes de vehículos,
lanzando nuevos modelos al mercado. (Paul Gong, 2023).
Argentina, es otro país en América del Sur, que ha dado iniciativas en el uso de
vehículos eléctricos e híbridos, y propone que al año 2041, se deje de
comercializar vehículos híbridos; el Ministerio de Desarrollo Productivo, en
coordinación con la multinacional Toyota, propone impulsar la inserción de
sectores científicos y tecnológicos a fin de dinamizar el uso de los vehículos con
celda de combustible. (Matías Kulfas, 2021).
El fabricante Honda, planteó que para el año 2040, no fabricará vehículos con
motor de combustión interna; y a partir de dicho año los vehículos eléctricos,
serán los que comercialicen; el acuerdo es entre Honda y General Motor, en el
cual utilizarán las baterías Ultium ,con gran capacidad de almacenamiento de
energía química. Así mismo tiene el objetivo que para el año 2050, eliminar la
huella de carbono, como consecuencia de las emisiones de sus unidades, para
lo cual invertirá 3000 millones de Euros entre el 2022 y 2027. (Toshihiro Mibe,
2021).
Este trabajo fue realizado por César José Manuel Ayasta Casiano, realizando las
pruebas a fin de comparar rendimientos técnicos y económicos entre un vehículo
híbrido y un convencional de la misma marca Toyota; logró determinar que el
vehículo híbrido Toyota Prius C (con dos motores, uno de combustión interna y
el otro eléctrico), tuvo una eficiencia del 33.7% mientras que el Toyota Yaris (con
un solo motor de combustión interna) tuvo una eficiencia del 26.9%. En cuanto a
los niveles de consumo específico de combustible, el vehículo híbrido utilizado
como vehículo de prueba en esta investigación, consumió 1.12 galones de
Gasohol por 100 km de recorrido, (equivalente a 42.39 ml/km), mientras que el
vehículo convencional que se utilizó para la prueba, tuvo un consumo de 1.60
galones de Gasohol por 100 km de recorrido ((equivalente a 60.56 ml/km). El
torque y la potencia mecánica del vehículo híbrido son mayores que el
convencional, (para el vehículo híbrido el torque máximo es de 142 N.m a 2800
RPM, y 73 kW a 5200 RPM, para el vehículo convencional el torque máximo es
de 121 N.m a 4400 RPM, y 63 kW a 6000 RPM de igual manera su costo de
adquisición del vehículo híbrido es mayor que el convencional. En lo referente al
grado de contaminación por los gases de escape, para un recorrido de 100000
km el vehículo híbrido emitió 9.77 toneladas de Dióxido de Carbono (C02),
mientras que para ese mismo recorrido el vehículo convencional emitió 13.95
toneladas de Dióxido de Carbono (C02)
21
Investigación realizado por Iván Alexander Calispa Procel y Stalin Antonio García
Suárez, muestra como resultados que los costos de operación por combustible
en el vehículo híbrido está en función al modo de funcionamiento, para el Modo
Sport, en un recorrido el cual incluye 42% de superficie llana, pendientes
menores a 3° con recorrido de 17%, pendientes de 4° recorrido de 16%; el costo
fue de 3.84 dólares, y para el mismo recorrido en Modo ECO, el costo fue de
1.92 dólares. En la carga de batería para el funcionamiento en modo Eco fue
menor en 4.9 kW.h con respecto a la carga de batería en modo Sport. En lo que
respecta a la activación del motor eléctrico, este tuvo una participación de
17.44% en comparación con el funcionamiento del motor de combustión interna.
Presentada por la Autora Raquel De los Santos Salas Ibarra, la cual utilizó una
metodología basada en las mediciones de vibraciones en el motor eléctrico del
vehículo híbrido de prueba, que acciona al tren de potencia denominado Blue
Drive, siguiendo lo estipulado en la norma ISO 10816-3; consistió en la siguiente
secuencia de procesos: adquisición, ubicación del sensor, tiempo de medición.
La recolección de los datos se hizo a través del software computacional
Dewesoft, que es para el análisis de vibraciones. Para el análisis de los
resultados utilizó las herramientas de la estadística inferencial, como es el
análisis de varianza, logrando obtener el espectro de frecuencias, que es
característico de cada condición de consumo de combustible.
Las pruebas fueron estáticas, se tuvo que simular el consumo energético en
procesos de carga y descarga; se utilizó una carga adicional que fue la carga del
sistema de climatización, utilizando un tiempo de 10 segundos, con lo cual se
obtuvo mil muestras por segundo en cada tratamiento. Hizo dos ensayos por
cada tratamiento, obteniendo diferencias no significativas. Como conclusión de
la investigación es que la caracterización de cada espectro de las vibraciones,
son relativamente iguales.
23
.
Figura 4 Diagrama de Bloques Motor de Combustión Interna.
Fuente: Sagastume, 2016
Los tiempos variables del cerrado de la válvula de admisión, conocido por sus
siglas TVCVA), es la configuración que utiliza el Ciclo Atkinson, se tiene con ello
relación de compresión más pequeñas, pero hay mayor expansión; con ello el
consumo de combustible es entre 1 y 2% menor con respecto al motor con ciclo
otto. (CONUE, 2020)
28
𝛿𝑇𝑚
𝛿𝑃𝑚𝑒 4𝜋 × 𝛿𝑡
=
𝛿𝑡 𝑉𝑐
Dónde:
𝛿𝑃𝑚𝑒
: Variación de la Presión Media Efectiva.
𝛿𝑡
𝛿𝑇𝑚
: Variación del Torque Mecánico.
𝛿𝑡
𝜕𝑃𝑚
: Variación de Potencia Mecánica.
𝜕𝑡
𝑅𝑃𝑀: Número de revoluciones por minuto.
𝛿𝑇𝑚
: Variación del Torque Mecánico.
𝛿𝑡
Dónde:
eb: Fuerza electromotriz inducida
Flujo Enlazado.
N: Número de vueltas por espira
; Flujo magnético.
m: Velocidad Angular
Dónde:
F: Fuerza generada.
l: Longitud del conductor.
B: Inducción Magnética.
T: Par Motor.
Rro: Radio del rotor.
Velocidad de Sincronismo
Dónde:
n: Velocidad de giro en Revoluciones por minuto
f: Frecuencia Eléctrica.
p. Número de par de polos.
Las baterías constituyen uno de los elementos importantes dentro del flujo de
energía para el accionamiento de un vehículo híbrido; es aquí donde parte de la
energía mecánica que genera el motor de combustión interna de un vehículo
híbrido, se transforma en energía química, para su posterior uso como energía
mecánica, que finalmente permitirá el accionamiento del vehículo. La cantidad
de carga que la batería almacene es un gran reto de las tecnologías actuales, y
que los fabricantes de vehículos siguen innovando, así como también en las
prestaciones de éstos y en su dimensionamiento. El Objetivo final de una batería
en el funcionamiento de un vehículo híbrido es tener el rendimiento, en cuanto a
la autonomía, potencia y energía que esta puede adquirir y suministrar, de
acuerdo a las condiciones de funcionamiento del vehículo híbrido. (Vehículos
híbridos. Pasado, presente y futuro, 2020).
Las baterías que se utilizan para el almacenamiento de energía química son las
de Plomo ácido, que tienen baja densidad energética, son las que actualmente
predominan el mercado automotriz; las baterías Níquel – Cadmio, su principal
característica es la durabilidad, pero en el sector automotriz presentan un alto
costo; las baterías Níquel -MH, son las que más se utilizan en vehículos híbridos,
37
Las baterías de celda de Ión – Litio, utilizan al carbono en forma de grafito, como
electrodo negativo. El cátodo está fabricado de LiCoO 2 (óxido de litio-cobalto),
LiMn2O4 (también conocido como espinela) o LiFePO4 (fosfato de hierro y litio).
El electrolito de la batería, es una solución de sal de litio en un disolvente
orgánico, como es el carbonato de etileno (C3H4O3); durante el proceso de
descarga, los iones de litio hacen su desplazamiento desde el ánodo hacia el
cátodo, y en proceso de carga, lo hacen de forma inversa. (Sánchez, 2020).
𝐿𝑖 ↔ 𝐿𝑖 + + 1𝑒 −
𝐹𝑒𝑃𝑂4 + 1𝑒 − ↔ 𝐹𝑒𝑃𝑂4−
38
𝑛
𝑛
𝑠= ∑ [(𝑦𝑖𝑡 − 𝑦𝑖 ) 2]
= ∑ (𝑎𝑜 + 𝑎1 𝑎𝑖 − 𝑦𝑖 )2
𝑖=1
𝑖=1
𝑛 ∑ 𝑥𝑖 𝑦𝑖 − ∑ 𝑥𝑖 − ∑ 𝑦𝑖
𝑎1 = 2
𝑛 ∑ 𝑥𝑖 2 − (∑ 𝑥𝑖 )
∑ 𝑦𝑖 𝑎𝑖 ∑ 𝑥𝑖
𝑎0 = −
𝑛 𝑛
Dónde n, representa el número de datos.
Dónde:
𝛽𝑜 : Es la ordenada en el punto de origen del plano cartesiano; es el
valor de la variable dependiente y, cuando todos los predictores son
cero.
𝛽𝑖 : Es el valor del promedio que tiene el incremento de una unidad de
la primera variable predictora 𝑥𝑖 , sobre la variable y, manteniendo
constante el resto de variables. También se les denomina coeficientes
parciales de regresión múltiple.
𝑒𝑖 : Es el valor de residuo o de error; es la diferencia entre el valor
dado en el experimento y el estimado por el modelo.
𝑦 𝑡 (𝑥1 , 𝑥2 ) = 𝑎𝑜 + 𝑎1 𝑥1 + 𝑎2 𝑥2
2
𝑎𝑜 ∑ 𝑥1,𝑖 + 𝑎1 ∑ 𝑥1,𝑗 + 𝑎2 ∑ 𝑥1,𝑖 . 𝑥2,𝑖 = ∑ 𝑥1,𝑖 . 𝑦𝑖
2
𝑎𝑜 ∑ 𝑥2,𝑖 + 𝑎1 ∑ 𝑥1,𝑖 . 𝑥2,𝑖 + 𝑎2 ∑ 𝑥2,𝑖 + = ∑ 𝑥2,𝑖 . 𝑦𝑖
41
Regresión No Lineal
Es un método que permite encontrar un modelo de regresión no lineal, entre la
variable dependiente y un conjunto de variables independientes. Para la
generación de este modelo, se utiliza los parámetros, que son estimados por el
procedimiento de cálculo del modelo de regresión no lineal: éstos pueden ser
constantes aditivas, coeficientes multiplicativos y/o exponentes. Las
aproximaciones de los modelos no lineales pueden ser lineales, cuadrática,
exponencial, logarítmica e Inversa. (Pico, 2019)
42
∑ yi = na + b ∑ xi + c ∑ xi 2
i=1 i=1 i=1
n n n n
∑ xi yi = a ∑ xi + b ∑ xi + c ∑ xi 3 2
∑ xi yi = a ∑ xi + b ∑ xi + c ∑ xi 4
2 2 3
Es la representación de las
Variable relaciones numéricas de los Variaciones
Regresión Lineal Modelamiento
Independiente: parámetros de funcionamiento de Lineales
Regresión No Lineal computacional con
Modelación los mecanismos que transforman Variaciones No
Regresión Multivariable software libre
Energética la energía, con la finalidad de Lineales
accionar a un vehículo.
47
1.5. Hipótesis.
Hipótesis General.
La modelación energética determina los flujos óptimos de energía en un vehículo
Automotriz Liviano Híbrido.
Hipótesis Específica.
Las relaciones de la variación de la Velocidad del vehículo híbrido con los
parámetros de los motores de combustión interna y motores eléctricos son
lineales y no lineales.
48
𝑁. 𝑍𝑎 2 . 𝑝. 𝑞
𝑛= 2
𝑑 . (𝑁 − 1) + 𝑍𝑎 2 . 𝑝. 𝑞
Dónde:
n: Tamaño de la muestra.
𝑍𝑎 : Nivel de Confianza = 1.96 si la seguridad es del 95%
𝑁: Total de la Población = 42
𝑝: Proporción esperada 5%
𝑞 1−p
𝑑: Precisión 5%
Técnicas.
Observación: consistió básicamente en observar la variación de los parámetros
de funcionamiento del motor de combustión interna y motor eléctrico
Instrumentos.
En la tabla 3, se enumeran en los anexos, la ficha de recolección de datos y el
objetivo que se persigue en el registro de los valores de medición de los datos
de vuelo de los vehículos híbridos de prueba.
50
Anexo Medición de velocidad del vehículo, torque Para las 2 condiciones de funcionamiento de los
2 motor eléctrico e intensidad de corriente 3 vehículos híbridos de prueba, registre el valor
eléctrica a condiciones de funcionamiento: del torque del motor eléctrico (TE) en N.m, y de
arranque y baja velocidad intensidad de corriente eléctrica que ingresa al
motor eléctrico (IEi) en Amperios, a cada valor
de velocidad del vehículo (v) en km/h
Anexo Medición de velocidad del vehículo, torque de Para la condición de funcionamiento de los 3
3 motor de combustión interna, torque motor vehículos híbridos de prueba, registre el valor
eléctrico e intensidad de corriente eléctrica a del torque del motor eléctrico Tm en N.m, torque
condiciones de funcionamiento: velocidad de del motor eléctrico (TE) en N.m, y de intensidad
crucero de corriente eléctrica que ingresa al motor
eléctrico (IEi) en Amperios, a cada valor de
velocidad del vehículo (v) en km/h
Anexo Medición de velocidad del vehículo, torque de Para la condición de funcionamiento de los 3
4 motor de combustión interna, torque motor vehículos híbridos de prueba, registre el valor
eléctrico e intensidad de corriente eléctrica a del torque del motor eléctrico Tm en N.m, torque
condiciones de funcionamiento: aceleración y del motor eléctrico (te) en N.m, y de intensidad
pendiente de corriente eléctrica que ingresa al motor
eléctrico (IMe) en Amperios, a cada valor de
velocidad del vehículo (v) en km/h
Anexo Medición de velocidad del vehículo, intensidad Para la condición de funcionamiento de los 3
5 de corriente que genera MG2 (Ig) a vehículos híbridos de prueba, registre el valor de
condiciones de funcionamiento: deceleración o la intensidad de corriente que genera MG2 (Ig)
frenado regenerativo en Amperios a cada valor de velocidad del
vehículo (v) en km/h
ENERGÍA QUÍMICA
MOTOR DE COMBUSTION
INTERNA
ENERGÍA MECANICA
ENERGÍA MECANICA
GENERADOR DE
ELECTRICO MECANISMOS
ENERGÍA ELECTRICA
BANCO DE BATERIAS
TRANSFORMADOR
ELEVADOR
ENERGÍA ELECTRICA
ENERGÍA MECANICA
RUEDAS MOTRICES
ENERGÍA DEL
COMBUSTIBLE
ENERGÍA MECÁNICA
Fuente: Evaluación energética y exergética en un motor de combustión interna ciclo Otto de 1.6
L, 2018
ENERGÍA MECÁNICA
Generador Eléctrico
En la figura 22, se esquematiza la energía de ingreso y salida del Banco de Baterías del
vehículo, y se observa que la energía eléctrica proveniente del generador eléctrico se
almacena en forma de energía química, para luego salir como en energía eléctrica.
ENERGÍA ELÉCTRICA
BANCO DE
Carga / Descarga
BATERIAS
ENERGÍA ELÉCTRICA
ENERGÍA ELÉCTRICA
TRANSFORMADOR ELEVADOR
ENERGÍA ELÉCTRICA
ENERGÍA TERMICA
ENERGÍA ELÉCTRICA
Motor Eléctrico
Fuente: Análisis del funcionamiento de la máquina asíncrona mediante el uso del software
MATLAB/Simulink
Fuente: Análisis del funcionamiento de la máquina asíncrona mediante el uso del software
MATLAB/Simulink
Fuente: Análisis del funcionamiento de la máquina asíncrona mediante el uso del software
MATLAB/Simulink
TRANSMISIÓN
MECÁNICA
Tabla 19. Energía mecánica provenientes del motor de combustión interna que
ingresan a la transmisión mecánica
Parámetro Símbolo Unidad Tipo
Tabla 20. Energía mecánica provenientes del motor eléctrico que ingresan a la
transmisión mecánica
Parámetro Símbolo Unidad Tipo
Motor
Potencia máxima Trasera (kW) 40 40 40
Eléctrico
Voltaje (V) 650 650 650
Torque máximo (N.m) 121 121 121
Potencia máxima combinada
163/5700
(kW/RPM)
Batería Níquel Batería Níquel Batería Níquel
Tipo – Metal – Metal – Metal
Batería Hidruro (NiMh) Hidruro (NiMh) Hidruro (NiMh)
Voltaje Nominal (pack batería) (V) 244.8 244.8 244.8
Transmisión Transmisión Transmisión
automática automática automática
Transmisión Tipo eCVT, eCVT, eCVT,
continuamente continuamente continuamente
variable. variable. variable.
Fuente: Sistema Híbrido auto-recargable, HV, 2.4
65
Descripción de la
Condición de
condición de Flujo de Energía
Funcionamiento
funcionamiento
Protocolo de Pruebas.
Torque
Velocidad 𝑇𝑚 Torque 𝑇𝐸 Velocidad v Torque Torque 𝑇𝐸 Velocidad v Torque 𝑇𝑚 Torque 𝑇𝐸
v (km/h) (N.m) (N.m) (km/h) 𝑇𝑚 (N.m) (N.m) (km/h) (N.m) (N.m)
0 35 24.4 0 39.1 23.5 0 29.8 28.1
1.2 34.1 23.2 0.8 37.5 22.1 1.1 27.8 26.3
Arranque
2.1 31.1 22.1 1.8 29.7 20.7 2.1 26.1 22.1
3.1 27.8 20.7 2.9 28.1 18.7 3.1 25.1 21.2
4.2 16.4 12.2 4.1 16.5 13.8 3.8 14.5 11.2
5.1 18.5 12.4 5 17.8 14.9 4.9 14.6 11.5
6.1 20.9 12.8 5.9 20.1 15.2 6.1 14.7 11.8
7.4 21.2 13.2 7.2 21.8 15.8 6.9 14.9 12.3
8.1 22.2 13.8 8.2 22.8 16.4 7.9 15.1 12.7
8.9 23.3 14.7 9.1 23.8 16.9 9.2 15.4 13.1
9.8 24.5 15.9 10.1 25.1 17.8 10.3 15.9 13.6
10.7 25.7 17.9 10.9 26.4 18.9 11.1 17.5 14.1
12.1 27.0 18.3 11.9 27.7 19.9 11.9 17.5 14.7
12.9 28.4 19.8 12.8 29.2 20.8 12.7 19.2 15.7
Baja
13.8 29.8 21.2 14.1 31.2 22.1 13.8 24.5 16.7
Velocidad
14.7 31.3 22.7 14.9 33.5 23.2 14.7 24.5 18.1
15.7 33.8 23.9 15.6 34.1 24.7 15.4 27.2 19.8
16.4 34.5 25.7 16.4 37.1 26.1 15.9 29.2 22.2
17.8 36.2 27.7 17.2 39.1 28.4 16.8 31.2 24.8
18.5 38.0 29.7 18.1 41.8 30.2 17.9 34.5 25.4
19.1 39.4 31.7 18.8 43.1 32.5 18.7 38.5 28.4
20 41.9 34.5 19.7 45.4 34.5 19.9 41.2 31.2
21.1 44.0 35.7 20.7 47.0 37.2 21.3 47.2 34.5
22.4 46.2 37.9
23.1 48 41.2
72
31.2 66.3 60.8 52.2 29.8 61.0 59.8 51.9 27.7 59.3 60.8 52.1
31.3 66.3 60.9 52.2 30.7 61.2 59.9 51.9 28.1 60.7 60.8 52.1
32.4 66.4 60.9 52.4 31.4 61.2 59.9 51.9 28.4 61.3 60.8 52.1
32.9 66.7 61 52.4 32.1 62.9 59.9 52 28.9 62.4 60.8 52.1
32.9 66.8 61.2 52.4 32.9 63.8 59.9 52 29.1 63.4 60.9 52.1
33.2 66.8 62.2 52.4 33.2 64.0 59.9 52 30.8 63.4 60.9 52.1
33.4 66.9 61.3 52.4 33.9 64.1 59.9 52 31.5 64.3 60.9 52.2
33.5 67.2 61.3 52.5 34.5 64.3 60.0 52 32.4 64.5 60.9 52.3
33.9 67.2 61.3 52.6 34 64.3 60.0 52 33.7 64.6 60.9 52.4
34.4 67.4 61.3 52.6 34.4 64.4 60.3 52 34.1 64.9 61 52.4
34.8 67.6 61.3 52.6 34.8 64.6 60.3 52.2 34.7 65.1 61.1 52.4
35.2 67.8 61.3 52.6 35.2 64.7 60.3 52.2 35.1 65.3 61.1 52.5
35.4 67.9 61.3 52.7 35.4 64.8 60.3 52.3 35.4 65.4 61.2 52.5
73
2/2
35.9 68.2 61.4 52.7 35.9 65.1 60.3 52.3 35.9 65.2 61.2 52.5
36.1 68.3 61.5 52.7 36.1 66.5 60.3 52.3 36.1 66.8 61.3 52.5
36.7 68.6 61.5 52.8 36.7 67.4 60.3 52.3 36.7 67.7 61.3 52.6
37.8 68.7 61.5 52.9 37.8 69.1 60.3 52.3 37.8 68.2 61.3 52.6
38.1 69.1 61.6 52.9 38.1 69.3 60.5 52.3 38.1 69.1 61.3 52.6
38.5 69.3 61.7 53 39.1 68.6 60.5 52.3 39.1 68.3 61.5 52.8
40.2 70.9 61.9 53 39.9 70.1 60.5 52.5 39.9 69.9 61.5 52.8
40.9 71.8 61.9 53 40.6 73.8 60.7 52.5 40.6 73.3 61.5 53
41.8 73.7 61.9 53.2 42.1 75.8 60.7 52.6 42.1 76.5 61.6 53
42.5 75.2 62.2 53.2 42.4 77.5 60.9 52.6 42.4 77.3 61.6 53.1
42.9 77.1 62.2 53.3 42.8 80.6 60.9 52.6 42.8 79.5 61.6 53.2
43.8 79.6 62.3 53.4 43.8 81.7 61.0 52.8 43.8 81.5 61.6 53.3
44.4 81.9 62.3 53.5 44.3 82.8 61.2 52.9 44.3 83.5 61.7 53.3
44.9 84.3 62.3 53.5 44.7 86.2 61.2 52.9 44.7 87.1 61.7 53.3
45.5 86.8 62.5 53.5 45.5 87.7 61.4 52.9 45.5 86.9 61.7 53.3
45.9 89.6 62.5 53.5 45.9 91.2 61.4 53 45.9 91.2 61.7 53.3
46.5 92.0 62.7 53.8 46.5 92.9 61.6 53 46.5 94.6 61.7 53.4
47.1 95.3 62.7 53.9 46.9 96.2 61.9 53 46.9 97.7 61.7 53.5
47.8 99.6 62.7 54.1 47.7 100.8 61.9 53.1 47.4 99.8 61.7 53.6
48.4 105.2 62.8 54.4 48.1 105.7 61.9 53.3 47.9 106.2 61.9 53.7
48.9 109.3 62.8 54.4 48.4 107.5 62.1 53.5 48.3 108.4 61.9 53.7
49.5 112.8 63.0 54.7 48.6 112.5 62.1 53.7 48.6 110.8 61.9 53.7
50.2 115.9 63.1 54.9 49.1 112.8 62.3 53.8 48.9 112.5 62 53.9
50.9 121.0 63.2 55.1 49.8 116.5 62.3 53.8 49.3 116.0 62 53.9
51.4 126.0 63.2 55.3 50.1 118.8 62.3 54 49.8 119.3 62 53.1
Fuente: Mediciones Realizadas
74
2/2
2/2
Tabla 28. Registro de Medición de Consumo de Combustible para dos condiciones de funcionamiento
2/2
Tabla 31. Regresión no lineal múltiple de velocidad con torque de MCI y Torque
motor eléctrico
Vehículo de Coeficiente de
Relación Numérica
Prueba correlación múltiple
Figura 32. Modelamiento Numérico de Velocidad de vehículo y torque de MCI y Torque Motor eléctrico. Condición crucero vehículo 1
Fuente: Elaboración Propia.
88
Figura 33. Modelamiento Numérico de Velocidad de vehículo y torque de MCI y Torque Motor eléctrico. Condición crucero vehículo 2
Fuente: Elaboración Propia.
89
Figura 34. Modelamiento Numérico de Velocidad de vehículo y torque de MCI y Torque Motor eléctrico. Condición crucero vehículo 3
Fuente: Elaboración Propia.
90
Tabla 33. Regresión no lineal velocidad con torque de MCI y Torque motor
eléctrico, condición de aceleración y pendiente
Coeficiente de
Vehículo
Relación Numérica correlación
de Prueba
múltiple
Figura 36. Modelamiento Numérico de Velocidad de vehículo y torque de MCI y Torque Motor eléctrico. Condición Aceleración y
Pendiente vehículo 1.
Fuente: Elaboración Propia.
95
Figura 37. Modelamiento Numérico de Velocidad de vehículo y torque de MCI y Torque Motor eléctrico. Condición Aceleración y
Pendiente vehículo 2
Fuente: Elaboración Propia.
96
Figura 38. Modelamiento Numérico de Velocidad de vehículo y torque de MCI y Torque Motor eléctrico. Condición Aceleración y
Pendiente vehículo 3
Fuente: Elaboración Propia.
97
2/3
Velocidad v Grado de Grado de Grado de
(km/h) Hibridación HV % Velocidad (km/h) Hibridación HV % Velocidad (km/h) Hibridación HV %
44 60 44 63 44 69
45 60 45 62 45 64
46 61 46 57 46 64
47 58 47 57 47 65
48 56 48 57 48 65
49 56 49 54 49 61
50 55 50 54 50 59
51 55 51 54 51 56
52 56 52 53 52 57
53 56 53 53 53 56
54 54 54 53 54 56
55 54 55 53 55 56
56 54 56 53 56 52
57 55 57 53 57 52
58 51 58 51 58 52
59 51 59 51 59 50
60 50 60 51 60 50
61 48 61 50 61 48
62 48 62 47 62 48
63 48 63 47 63 48
64 48 64 47 64 49
65 48 65 47 65 49
66 46 66 45 66 47
67 46 67 45 67 47
68 46 68 45 68 47
69 46 69 45 69 46
70 45 70 46 70 47
71 42 71 46 71 44
72 42 72 43 72 44
73 42 73 42 73 44
74 42 74 42 74 42
75 42 75 42 75 42
76 42 76 41 76 41
77 42 77 40 77 41
78 41 78 38 78 43
79 41 79 37 79 42
80 40 80 38 80 39
81 39 81 38 81 38
82 39 82 36 82 38
83 39 83 35 83 38
84 39 84 36 84 38
85 39 85 35 85 39
86 37 86 34 86 37
87 37 87 37 87 35
88 37 88 33 88 35
89 37 89 32 89 36
103
3/3
Velocidad v Grado de Grado de Grado de
(km/h) Hibridación HV % Velocidad (km/h) Hibridación HV % Velocidad (km/h) Hibridación HV %
90 37 90 31 90 36
91 37 91 32 91 34
92 34 92 32 92 34
93 34 93 31 93 34
94 34 94 31 94 32
95 34 95 29 95 32
96 34 96 30 96 32
97 32 97 28 97 33
98 32 98 27 98 29
99 32 99 27 99 29
100 30 100 25 100 27
101 31 101 26 101 28
102 31 102 26 102 28
103 28 103 25 103 26
104 28 104 24 104 26
105 24 105 24 105 25
106 24 106 24 106 25
107 24 107 23 107 25
108 22 108 23 108 26
109 22 109 21 109 24
110 20 110 18 110 23
111 11 111 18 111 13
112 11 112 18 112 5
113 12 113 12 113 0
114 10 114 11 114 0
115 9 115 10 115 0
116 0 116 7 116 0
117 0 117 0 117 0
118 0 118 0 118 0
119 0 119 0 119 0
120 0 120 0 120 0
Fuente: Mediciones realizadas
104
Se concluye que:
v= 16.22 km/h, valor que está fuera del tango de la velocidad crucero del
vehículo 1, que es de 31.2 a 51.4 km/h.
107
v= 26 km/h, valor que está fuera del tango de la velocidad crucero que es
de 27.7 a 49.8 km/h
108
Vehículo de Coeficiente de
Relación Numérica
Prueba determinación
v= 91.52 km/h, valor que está dentro del rango de la velocidad crucero
que es de 57.2 a 111.8 km/h.
Reemplazando el valor hallado de v en la expresión de Q, se tiene:
𝑚𝑙
𝑄 = 0.0059 ∗ 91.522 − 1.081 ∗ 91.52 + 86.047 = 36.53 ,
𝑘𝑚
Los valores mínimos encontrados en cada uno de los vehículos de pruebas, tienen un
promedio de consumo específico de combustible de 38.16 ml/km, a una velocidad
promedio de 94.91 km/h.
111
+ 0.16𝑇𝑚. 𝑇𝐸 − 5.35𝑇𝐸 2
− 0.306𝑇𝑚. 𝑇𝐸 + 0.0046𝑇𝐸 2
− 0.17𝑇𝑚. 𝑇𝐸 − 5.63𝑇𝐸 2
Sustituyendo, se tiene:
Sustituyendo, se tiene:
Vehículo 1 78.52
Vehículo 2 77.45
Vehículo 3 79.97
De la tabla 41, se concluye que el aporte del torque del motor eléctrico hace
posible un menor consumo de combustible, debido a que el motor de combustión
interna del motor en esta condición de funcionamiento de aceleración y
pendiente tiene sus valores máximos de entrega de torque mecánico, llegando
a tener un valor máximo, de acuerdo a lo determinado en el item 4.3.1 que fue
de 38.16 ml/km.
118
La energía mecánica tiene dos parámetros que le dan esa característica, y son
el torque mecánico y la velocidad angular, y para la energía eléctrica es la tensión
eléctrica y la intensidad de corriente eléctrica. La energía mecánica proveniente
del motor de combustión interna, incrementa de manera no lineal el torque, sin
embargo la velocidad de rotación lo hace de manera lineal; es por ello que para
la investigación, el parámetro que tuvo mas incidencia del motor de combustión
interna fue el torque mecánico. En el caso del motor eléctrico, existe un valor
estable y constante del nivel de tensión del motor eléctrico trifásico de corriente
120
CONCLUSIONES
RECOMENDACIONES
REFERENCIAS BIBLIOGRAFICAS
..
un motor de combustión interna ciclo Otto de 1.6 L. Enfoque UTE, 2018, vol.
9, no 4, p. 221-232.
• MARTÍNEZ, Zaira Itzel Chávez. Capacidades Para la Implementación del
Sistema de Comercio de Emisiones en México: Análisis Comparado con
California. 2022. Tesis Doctoral. Centro de Investigacion y Docencia
Economicas (Mexico).
• MASCARENHAS, Tomás Bril, et al. Políticas de Desarrollo Productivo Verde
para la Argentina. 2021.
• MINELLONO, José Armando. Evaluación técnica, económica y ambiental de
la implementación de autobuses híbridos en la ciudad de Mar del Plata. 2015.
Tesis Doctoral. Universidad Nacional de Mar del Plata. Facultad de Ingeniería.
Argentina.
• MOLINA VALLEJO, Luis Carlos; VALENCIA PAZOS, Santiago
Federico. Diseño y construcción de un prototipo de vehículo híbrido ecológico
para movilidad urbana en la ciudad de Riobamba mediante el uso del software
CAD, CAM, CAE. 2020. Tesis de Licenciatura. Escuela Superior Politécnica
de Chimborazo.
• MORENO FERNÁNDEZ, Alberto. Desarrollo de un modelo del sistema de
refrigeración de un Motor de Encendido Provocado para vehículo híbrido.
2021. Tesis Doctoral. Universitat Politècnica de València.
• ORTIZ, Óscar Olarte; MARTÍNEZ, Javier Enrique Solano; QUIROGA, Óscar
Arnulfo Quiroga. GESTIÓN ENERGÉTICA BASADA EN REGLAS PARA UN
VEHÍCULO HÍBRIDO-IEEE VTS MOTOR VEHICLES CHALLENGE
2020. Encuentro Internacional de Educación en Ingeniería, 2021.
• PAGUAY QUISHPE, Alex Fernando; RAMOS SÁNCHEZ, Alexander
Sebastián. Diseño y construcción de un equipo para diagnóstico de las
bobinas del estator para MG de vehículos híbridos. 2021.
• PÉREZ, Miguel Ángel Castro, et al. El Instituto de Investigaciones Económicas
y Sociales-IIES-, 2022.
• PICO GORDÓN, José Andrés, et al. Modelado por técnicas de regresión de
los parámetros energéticos de desempeño para gasificadores tipo
downdraft. Ingeniería Energética, 2019, vol. 40, no 2, p. 138-147.
128
ANEXO 1.
Objetivo: Para las 2 condiciones de Funcionamiento de los 3 Vehículos híbridos de prueba, registre el valor del consumo de combustible (Q) en
ml/km a cada valor de velocidad del vehículo (v) en km/h
Velocidad de Aceleración y
crucero Pendiente
130
ANEXO 2
MEDICIÓN DE VELOCIDAD DEL VEHÍCULO, TORQUE MOTOR ELECTRICO E INTENSIDAD DE CORRIENTE ELECTRICA A
CONDICIONES DE FUNCIONAMIENTO: ARRANQUE Y BAJA VELOCIDAD
Objetivo: Para las 2 condiciones de Funcionamiento de los 3 Vehículos híbridos de prueba, registre el valor del Torque del Motor Eléctrico (TE)
en N. m, y de Intensidad de corriente eléctrica que ingresa al motor eléctrico (IEi) en Amperios, a cada valor de velocidad del vehículo (v) en km/h
Arranque
Baja Velocidad
131
ANEXO 3
MEDICIÓN DE VELOCIDAD DEL VEHÍCULO, TORQUE DE MOTOR DE COMBUSTIÓN INTERNA, TORQUE MOTOR ELECTRICO E
INTENSIDAD DE CORRIENTE ELECTRICA A CONDICIONES DE FUNCIONAMIENTO: VELOCIDAD DE CRUCERO
Objetivo: Para la condición de Funcionamiento de los 3 Vehículos híbridos de prueba, registre el valor del Torque del Motor Eléctrico Tm en N.M,
Torque del Motor Eléctrico (TE) en N. m, y de Intensidad de corriente eléctrica que ingresa al motor eléctrico (IEi) en Amperios, a cada valor de
velocidad del vehículo (v) en km/h
ANEXO 4
MEDICIÓN DE VELOCIDAD DEL VEHÍCULO, TORQUE DE MOTOR DE COMBUSTIÓN INTERNA, TORQUE MOTOR ELECTRICO E
INTENSIDAD DE CORRIENTE ELECTRICA A CONDICIONES DE FUNCIONAMIENTO: ACELERACIÓN Y PENDIENTE
Objetivo: Para la condición de Funcionamiento de los 3 Vehículos híbridos de prueba, registre el valor del Torque del Motor Eléctrico Tm en N.M,
Torque del Motor Eléctrico (TE) en N. m, y de Intensidad de corriente eléctrica que ingresa al motor eléctrico (IME) en Amperios, a cada valor de
velocidad del vehículo (v) en km/h
ANEXO 5
MEDICIÓN DE VELOCIDAD DEL VEHÍCULO, INTENSIDAD DE CORRIENTE QUE GENERA MG2 (IG) A CONDICIONES DE
FUNCIONAMIENTO: DECELERACIÓN O FRENADO REGENERATIVO
Objetivo: Para la condición de Funcionamiento de los 3 Vehículos híbridos de prueba, registre el valor de la intensidad de corriente que genera
MG2 (Ig) en Amperios a cada valor de Velocidad del vehículo (v) en km/h
ANEXO 6
doble="Torque";
nd=strlength(doble);
Graf_ad = strncmp(titulos,doble,nd);
[k,l]=size(Graf_ad);
e=[];
for g=1:k
nm_ad=numel(find(Graf_ad(g,:)));
if nm_ad>=2
e=[e g];
end
end
v=0;
H = numel(hojas);
for h=1:H
num_a=titulos(h,:);
eliminar=find(num_a=="0");
num_a(eliminar) = [];
repeticiones=numel(num_a)-1;
% mkdir(char(hojas(h))); %%%%%%%%
Comentar una vez creadas las carpetas
x=titulos(h,1);
nx=strlength(x);
Busqueda_Velocidad = strncmp(raw,x,nx);
135
columnaV=find(Busqueda_Velocidad==1);
numcolX=floor(columnaV/a)+1;
Cantidad_vehiculos=numel(columnaV);
Velocidad=str2double(string(raw(:,numcolX)));
RL_1=zeros(Cantidad_vehiculos,2);
RL_2=zeros(Cantidad_vehiculos,3);
RL=zeros(Cantidad_vehiculos,3);
R2_1=zeros(Cantidad_vehiculos,1);
R2_2=zeros(Cantidad_vehiculos,1);
R2=zeros(Cantidad_vehiculos,1);
Velocidad_s=rmmissing(Velocidad(:,i));
[c,d]=size(Velocidad_s);
E=[];
for r=1:repeticiones
y=titulos(h,1+r);
ny=strlength(y);
Busqueda_Y= strncmp(raw,y,ny);
columnaY=find(Busqueda_Y==1);
numcolY=floor(columnaY/a)+1;
Hibridacion=str2double(string(raw(:,numcolY)));
Hibridacion_s=rmmissing(Hibridacion(:,i));
RL_1(i,:)=polyfit(Velocidad_s,Hibridacion_s,1);
RL_2(i,:)=polyfit(Velocidad_s,Hibridacion_s,2);
Calc1 = polyval(RL_1(i,:),Velocidad_s);
Calc2 = polyval(RL_2(i,:),Velocidad_s);
SSresid1 = sum(R_1.^2);
SSresid2 = sum(R_2.^2);
R2_1(i) = 1 - SSresid1/SStotal;
R2_2(i) = 1 - SSresid2/SStotal;
if R2_1(i)<R2_2(i)
Calc = Calc2;
RL(i,:) = RL_2(i,:);
R2(i) = R2_2(i);
grf=[num2str(RL(i,1)),'v^2 + ', num2str(RL(i,2)),'v +
', num2str(RL(i,3)),' | R^2 = ', num2str(R2(i))];
136
else
Calc = Calc1;
RL(i,1:2) = RL_1(i,:);
R2(i) = R2_1(i);
grf=[num2str(RL(i,1)),'v + ', num2str(RL(i,2)),' | R^2
= ', num2str(R2(i))];
end
M=find(Hibridacion_s==max(Hibridacion_s));
m=find(Hibridacion_s==min(Hibridacion_s));
f=figure(i+Cantidad_vehiculos*(h-1));
f.WindowState = 'maximized';
scatter(Velocidad_s,Hibridacion_s,'DisplayName',y)
hold on
plot(Velocidad_s,Calc,'DisplayName',grf)
if not(or(Decre,Cres))
prim_deriv=(-RL_2(i,2))/(RL_2(i,1)*2);
res_deriv=polyval(RL_2(i,:),prim_deriv);
if RL(i,1)<0
txt = ['\downarrow Punto Máximo = (',
num2str(prim_deriv), ' ; ',num2str(res_deriv), ')'];
text(prim_deriv,res_deriv+2,txt, 'Fontsize', 10)
elseif RL(i,1)>0
txt = ['\uparrow Punto mínimo = (',
num2str(prim_deriv), ' ; ',num2str(res_deriv), ')'];
text(prim_deriv,res_deriv-0.2,txt, 'Fontsize', 10)
end
end
hold on
xlabel(x)
title({hojas(h);['Vehículo de Prueba ',num2str(i)]})
grid on
xlim([min(Velocidad_s)-2 max(Velocidad_s)+5])
legend('Location','best')
if max(h==e)
f_G=find(Graf_ad(h,:));
if max(r==(f_G-1))
E=[E Hibridacion_s];
end
end
end
if max(h==e)
137
v=v+1;
for w=1:2
if w==1
pol='poly11';
reg='RLM'; %Regresión lineal múltiple
else
pol='poly22';
reg='RNLM'; %Regresión no lineal múltiple
end
fitsurface=fit([E(:,1),E(:,2)],Velocidad_s,
pol,'Normalize','off');
Coef_sur=coeffvalues(fitsurface);
RES_EC=Coef_sur(1) + Coef_sur(2).*E(:,1) +
Coef_sur(3).*E(:,2);
if w==2
RES_EC=RES_EC + Coef_sur(4).*E(:,1).^2 +
Coef_sur(5).*E(:,1).*E(:,2) + Coef_sur(6).*E(:,2).^2;
end
R = Velocidad_s - RES_EC;
SSresid = sum(R.^2);
SStotal = (length(Velocidad_s)-1) * var(Velocidad_s);
R2m = 1 - SSresid/SStotal;
ec_rm=[num2str(Coef_sur(1)),' + ',
num2str(Coef_sur(2)),' Tm + ', num2str(Coef_sur(3)),' TE '];
if w==2
ec_rm=[ec_rm,'+ ', num2str(Coef_sur(4)),' Tm^2 +
', num2str(Coef_sur(5)),' Tm*TE + ', num2str(Coef_sur(6)),' TE^2 '];
end
ec_rm=[ec_rm, '| R^2 = ', num2str(R2m)];
f=figure(Cantidad_vehiculos*H+w+2*(v-1));
f.WindowState = 'maximized';
plot(fitsurface, [E(:,1),E(:,2)],Velocidad_s)
hold on
xlabel(titulos(h,f_G(1)))
ylabel(titulos(h,f_G(2)))
zlabel(x)
view(50,10)
title({hojas(h);reg;['Vehículo de Prueba
',num2str(i)];ec_rm})
hold off
% name=[char(hojas(h)),'\',reg,'_Vehículo de Prueba
',num2str(i)]; %%%%%%%% Comentar una vez creadas las carpetas
% saveas(gcf,name,'png')
%%%%%%%% Comentar una vez creadas las carpetas
end
end
end
pause
end
138
Siendo las 4:00 p.m. del día martes 09 de mayo de 2023, se dio inicio a la
Sustentación Virtual de Tesis soportado por el sistema Google Meet, preparado
y controlado por la Unidad de Tele Educación de la Escuela de Posgrado de la
Universidad Nacional Pedro Ruiz Gallo de Lambayeque, con la participación en
la Video Conferencia de los miembros del Jurado, nombrados con Resolución Nº
0936-2019-EPG, de fecha 24 de julio de 2019, conformado por:
Para evaluar el informe de tesis del tesista JAMES SKINNER CELADA PADILLA,
candidato a optar el grado de DOCTOR EN CIENCIAS DE LA INGENIERIA MECANICA
Y ELECTRICA CON MENCION EN ENERGIA, con la tesis titulada “MODELACION
ENERGETICA PARA DETERMINAR LOS FLUJOS OPTIMOS DE ENERGIA EN UN VEHICULO
AUTOMOTRIZ LIVIANA HIBRIDO.”.
Yo; Amado Aguinaga Paz, Docente Asesor de tesos del trabajo de investigación
del estudiante James Skinner Celada Padilla
Titulada:
MODELACION ENERGETICA PARA DETERMINAR LOS FLUJOS OPTIMOS
DE ENERGIA EN UN VEHICULO AUTOMOTRIZ LIVIANA HIBRIDO, luego de
la revisión exhaustiva del documento, constato que la misma tiene un índice de
similitud del 16%, verificable en el reporte de similitud del programa Turnitin.
El suscrito analizó dicho reporte y concluyó que cada una de las coincidencias
no constituye plagio. A mi leal saber y entender la tesis cumple con todas las
normas para el uso de citas y referencias establecidas por la Universidad
Nacional Pedro Ruiz Gallo
DNI 16527442
141
142
143
144
145
146
147
148