ALEGRE BUJ Susana Tesis PDF
ALEGRE BUJ Susana Tesis PDF
ALEGRE BUJ Susana Tesis PDF
INGENIERA INDUSTRIAL
2017
TESIS DOCTORAL
MODELADO DEL VEHÍCULO ELÉCTRICO E
HÍBRIDO PARALELO POR MEDIO DE
MATLAB/SIMULINK Y PLANIFICACIÓN DE
ESTACIONES DE CARGA MEDIANTE SISTEMAS
DE INFORMACIÓN GEOGRÁFICA Y
ALGORITMOS GENÉTICOS
INGENIERA INDUSTRIAL
2017
AGRADECIMIENTOS
A mis padres, María y César, por haberme forjado como soy en la actualidad y por
su apoyo, su ayuda y su paciencia, sin la cual no hubiera conseguido mi sueño.
Quiero agradecer también a los amigos que han contribuido a la realización de este
trabajo por su constante apoyo moral. Y así como a cuantos me han escuchado y
aconsejado, pese a que no constan sus nombres, en este reducidísimo espacio.
I
ÍNDICE
Resumen........................................................................................................................XXIV
3.5.- Localización de las estaciones de carga de los vehículos eléctricos ............. 132
Aceleración angular
b Rendimiento de la batería
g Rendimiento de la transmisión
Velocidad angular
a Aceleración líneal
AG Algoritmos genéticos
CARB Junta de Recursos del Aire de California (California’s Air Resources Board)
Cd Coeficiente aerodinámico
D Difusión
G Ratio de reducción
g Aceleración de la gravedad
k Constante de Peukert
l Longitud de difusión
Lb Distancia desde el centro de la gravedad del vehículo hasta las ruedas traseras
N Número de cliente
PC Pila de combustible
PL Potencia que demanda el conductor por medio del pedal del acelerador
q Resistencia de frenado
r Radio de la rueda
T Par motor
t Tiempo [h]
U0 Tensión en vacío
v Velocidad lineal
w Frecuencia angular
Zw Impedancia de Warburg
LISTA DE FIGURAS
Figura 1. Hitos históricos del vehículo eléctrico durante los siglos XIX y XX .......... 3
Figura 21. Fuerzas que actúan en un vehículo ascendiendo una cuesta .................... 59
Figura 45. Representación ráster y vectorial de la información del mundo real ...... 99
Figura 46. Gráfico resumen de los 4 modos de carga según IEC 61851 ................ 110
Figura 50. Compatibilidad de los diferentes tipos de tomas de corriente ............... 114
Figura 53. Circuito simple equivalente del modelo de la batería .............................. 122
Figura 55. Vehículo hibrido eléctrico usando un divisor EVT ................................. 125
Figura 60. Diagrama de bloques del diseño Smart Roadster eléctrico..................... 131
Figura 61. Disposición del sistema para cargar coches eléctricos con alternador .. 132
Figura 62. Localización óptima de las estaciones de carga realizado por Sara Mehar
[46]................................................................................................................................................... 133
Figura 63. Cobertura de las particiones y la capa de las estaciones de carga [49] .. 134
Figura 72. Autonomía del vehículo eléctrico en función del ciclo de conducción 150
Figura 75. Diagrama del flujo de energía del modo sólo eléctrico ........................... 156
Figura 76. Diagrama del flujo de energía del modo híbrido...................................... 157
Figura 77. Diagrama del flujo de energía del modo recarga de baterías .................. 158
Figura 78. Diagrama del flujo de energía del modo frenada regenerativa ............... 159
Figura 79. Diagrama del flujo de energía del modo térmico puro ........................... 160
Figura 80. Organigrama de los modos de control del vehículo eléctrico híbrido en
paralelo ........................................................................................................................................... 161
Figura 81. Diagrama de bloques del modelo del vehículo eléctrico híbrido en
paralelo ........................................................................................................................................... 161
Figura 83. Análisis del motor térmico: batería de ácido – plomo............................. 165
Figura 84. Análisis del motor eléctrico: batería de ácido – plomo ........................... 165
Figura 87. Análisis del motor térmico: batería de ión – litio ..................................... 167
Figura 88. Análisis del motor eléctrico: batería de ión – litio .................................... 167
Figura 89. Análisis de la frenada regenerativa: batería de ión – litio ........................ 168
Figura 98. Esquema del programa de planificación de estaciones de carga ............ 181
Figura 99. Relación de las hojas del fichero de Excel ................................................ 182
Figura 100. Ventana gráfica del programa de algoritmos genéticos donde da a elegir
la manera de obtener la solución inicial del programa de optimización ............................... 182
Figura 101. Ventana gráfica del programa de algoritmos genéticos donde hay que
especificar la longitud y número de cromosomas .................................................................... 183
Figura 102. Ventana gráfica del programa de algoritmos genéticos donde hay que
especificar el número de iteraciones ........................................................................................... 183
Figura 105. Mapa temático donde se representa la zona que abarca cada estación de
carga y la localización de las mismas .......................................................................................... 187
Figura 107. Ciclo de conducción UNECE Elementary Urban Cycle (Part 1) ....... 200
Figura 108. Ciclo de conducción UNECE Extra-Urban Driving Cycle (Part 2) ... 201
Figura 114. Vista general del modelo del vehículo eléctrico en Simulink ............... 215
Figura 115. Detalle de la dinámica del vehículo eléctrico en Simulink .................... 216
Figura 116. Detalle de los bloques de eficacia de la transmisión en Simulink ........ 217
Figura 119. Modelo del vehículo hibrido en paralelo en Simulink ........................... 220
Figura 128. Distribución de las fuerzas de frenado en las ruedas delanteras y las
traseras ............................................................................................................................................ 225
Figura 129. Fuerzas que actúan en un vehículo durante el frenado en una superficie
plana ................................................................................................................................................ 226
Figura 133. Bloque Comparador de consumos con diferentes marchas ................. 231
Figura 136. Línea de par de máxima eficiencia del motor térmico .......................... 233
LISTA DE TABLAS
Tabla 12: Ventajas y desventajas de los SIG vectoriales y los SIG ráster ............... 102
Tabla 13: Resumen de las principales características según el modo de carga ....... 109
Tabla 15: Tipos de toma de corriente estándar según la norma IEC 62196........... 113
Tabla 16: Mejoras propuestas por Millner [40] para los PHEV ............................... 128
Tabla 22: Especificaciones técnicas de las baterías para los vehículos híbridos
eléctricos ......................................................................................................................................... 163
Tabla 25: Parámetros básicos del ciclo urbano elemental ......................................... 200
Tabla 27: Parámetros básicos del ciclo de conducción extra-urbano (Vehículos con
baja potencia) ................................................................................................................................. 202
Tabla 28: Parámetros básicos del ciclo tráfico fluido ................................................. 203
Tabla 29: Parámetros básicos del ciclo tráfico no fluido ........................................... 203
Tabla 31: Variables del módulo1 de los algoritmos genéticos .................................. 259
RESUMEN
En las últimas décadas, el consumo de petróleo del sector del transporte ha crecido
a un ritmo mayor que en cualquier otro sector. El transporte representa la cuarta parte de
las emisiones de gases de efecto invernadero y el 36% del consumo de energía en España.
A consecuencia de ello, se ha dictado una norma europea sobre emisiones, la cual regula los
límites aceptables para las emisiones de gases de combustión de los vehículos nuevos
vendidos en los Estados Miembros de la Unión Europea.
El trabajo realizado versa sobre el vehículo eléctrico e hibrido para hacer frente a la
creciente contaminación ambiental y reducir el consumo de energía. En este documento se
presenta un modelado de los vehículos eléctricos y paralelos híbridos usando el entorno
Matlab / Simulink, el cual nos permite investigar rápidamente aspectos del vehículo, como
la potencia del motor, el tipo y tamaño de la batería, el peso, etc., y ver cómo los cambios
pueden afectar al rendimiento y a la distancia recorrida. El modelo es simulado para
obtener la autonomía del vehículo eléctrico. Mediante un sistema de información
geográfica y el uso de un algoritmo matemático basado en los algoritmos genéticos se
obtiene la planificación de las estaciones de carga minimizando el coste de la inversión de la
instalación y proporcionando la mejor distribución geográfica de las mismas para
incrementar la calidad de servicio al mejorar la fiabilidad.
ABSTRACT
In the last decades, the consumption of petrol in the transport sector has increased
at a higher pace than in any other sector.
Transport represents a quarter of the total amount of greenhouse gases and 36% of
energy consumption in Spain. As a consequence, a new law regarding emissions has been
approved. Such law regulates the acceptable emission levels of combustion gases for new
vehicles sold in any of the European Union Member States. This work deals with parallel-
hybrid electric vehicles in order to face the increasing environmental pollution and reduce
energy consumption. The present paper shows a modelling of electric and parallel-hybrid
electric vehicle using Matlab/Simulink environment which allows us to access different
aspects of the vehicle such as engine power, type and size of the battery or weight and to
observe how changes can affect the performance and the distance travelled. The model is
simulated in order to obtain the electric vehicle’s autonomy. Through the use of a
Geographic Information System together with a mathematic algorithm based on genetic
algorithms the planning of charging stations is obtained, where the installation investment
cost was minimized and the geographic distribution is improved in order to increase the
quality of the service by improving reliability.
Este capítulo consta de cuatro partes. En primer lugar, se desarrollan los objetivos;
a continuación, se analiza la evolución del vehículo eléctrico a lo largo de la historia;
después se analiza el impacto de éste sobre el medio ambiente; y finalmente, se hace una
breve presentación de los contenidos de la presente tesis doctoral.
1.1.- Objetivos
Las fuentes de energía son una batería y/o un motor de combustión interna.
El coche eléctrico fue uno de los primeros automóviles que se desarrollaron, hasta
el punto que existieron vehículos eléctricos anteriores al motor de cuatro tiempos sobre el
que Diesel (motor diésel) y Otto (gasolina), basaron el automóvil actual. En la siguiente
figura se muestran los principales hitos de la historia del vehículo eléctrico a lo largo de los
siglos XIX y XX.
Robert
Davidson
consiguió
mover una
locomotora General
a 6 km/h sin Motors
usar carbón presenta
ni vapor, EV-1.
por lo Iniciativa
tanto, fue "La Jamais legislativ en
el primer Primeras Contente" EEUU: Zero
vehículo baterías superó los Casi Emission
eléctrico recargables 100 Km/h desaparecen Mandate
Los Los
1838 1839 1880 1890 1899 1912 60 70 1990
Figura 1. Hitos históricos del vehículo eléctrico durante los siglos XIX y XX
La mejora en la tecnología de las baterías está consiguiendo que sean una alternativa
a considerar, pues la pila de combustible está todavía en una fase inicial que no interesa
para el público masivo. El coche convencional comienza su decadencia ahora.
Es por ello por lo que, en un principio, la utilización supuesta del vehículo eléctrico
es como vehículo de flota o como segundo vehículo para uso particular. Esto supone una
pequeña parte del total de automóviles, debido a lo cual no puede esperarse un cambio
radical de las condiciones medioambientales y energéticas.
Se debe pensar, eso sí, en el vehículo eléctrico como una opción de futuro, a largo
plazo. Por esta razón, se encuentran en pleno desarrollo nuevas soluciones aplicables a
vehículos de uso masivo y que actúen como “puente” hasta la consolidación del vehículo
eléctrico en el parque automovilístico.
sector del transporte ha crecido a un ritmo mayor que en cualquier otro sector. Este
incremento se ha debido principalmente a las nuevas exigencias de los vehículos de uso
personal con motores convencionales de combustión interna (ICE). Debido a que algunos
problemas ambientales, como el efecto invernadero, están directamente relacionados con
las emisiones de vehículos, agencias de gobierno y organizaciones han desarrollado normas
más estrictas para el consumo de combustible y las emisiones.[2]
Los inconvenientes de este modelo, entre los que destacan la contaminación del
aire, el consumo excesivo de energía, los efectos sobre la salud o la saturación de las vías de
circulación, han generado una preocupación creciente por encontrar alternativas que
ayuden a evitar o minimizar los efectos negativos de este modelo y encontrar uno nuevo.
Figura 2. Curva de producción del petróleo, como sugirió originalmente M. King Hubbert en 1956
Gases Derivados
Energías Carbón del Carbón
renovables 1,60% 0,30%
6,20%
Electricidad
23,40%
P. Petrolíferos
50,80%
Gas
17,70%
1.3.2.1.1.- Contaminantes
del aire produciéndose una oxidación rápida, que es como químicamente se conoce a la
combustión.
Los principales contaminantes emitidos, tanto por fuentes móviles como fijas, son:
Efecto invernadero
Lluvia ácida
“Smog” fotoquímico
Efecto invernadero
Lluvia ácida
El ‘smog’ fotoquímico
A estos valores deberemos añadir otra fuente de emisiones importante, sobre todo
en los vehículos de gasolina, y es la evaporación que se produce al repostar, tanto en el
interior del depósito (ya solucionado en los automóviles modernos) como en el surtidor.
Otro factor determinante para el grado de emisiones es la calidad del combustible.
En los motores de gasolina alimentados con una mezcla rica tienden a aparecer
hidrocarburos (HC) y CO, mientras que, si la mezcla es pobre, existe más tendencia a
producir NOx. Al aumentar la temperatura de la reacción, la emisión de este último
compuesto crece.
Un motor diesel, en cambio, funciona siempre con más aire del que necesita para
quemar el gasóleo, lo que significa que hay oxígeno de sobra. Esto produce un escaso
índice de CO, pero eleva la proporción de NOx, debido a la elevada temperatura del
proceso. Estos motores emiten, además, una gran cantidad de partículas sólidas,
responsables de los humos negros.
18%
1% CO2
9% CONTAMINANTES
71% H2O
O2 Y GASES INERTES
1% N2
Los niveles de emisión dependen del ratio de aceleración y deceleración, así como
del tiempo que se pase a bajas velocidades y en paradas con el motor en ralentí. Podemos
añadir que, en condiciones de tráfico que obliguen a una marcha muy lenta, el consumo de
combustible puede llegar a incrementarse hasta un 50 por 100, lo que ocasionará un mayor
ratio de emisiones. Por todo ello, se concluye que, en entornos urbanos, precisamente
donde el nivel de contaminación empieza a resultar preocupante, las características del
tráfico contribuyen a que la aportación polutiva del VCI, independientemente de que
incorpore o no aparato catalizador, sea mayor que en un régimen de conducción menos
abrupto.
Carbón (no controladas) 955 2,92 0,03 4,34 0,14 0,02 11,82
Renovables 0 0 0 0 0 0 0
Partículas
SO2 NOx HC CO CO2
de polvo
pueden atribuirse a los vehículos eléctricos son las asociadas con la generación de la energía
necesaria para cargar las baterías del vehículo.
Precisamente las baterías son un punto negro al que se hace referencia cada vez que
se considera al vehículo eléctrico como una alternativa limpia. Una batería es un elemento
altamente contaminante si no se tiene el debido cuidado con ella. Los ácidos y metales
pesados que contiene son tóxicos (plomo, cadmio, níquel, cinc, litio, ácido sulfúrico, azufre,
etc.) y en el proceso de carga se desprenden gases que son contaminantes y peligrosos. Se
hace, pues, necesario un proceso de reciclaje que, si es bueno, puede lograr la reutilización
del 90 por 100 de una batería en desuso.
En forma gráfica, Figura 7, tomando como referencia las emisiones de los vehículos
a gasolina se obtiene:
Figura 7: Distribución relativa de las emisiones según los tipos de contaminantes y combustible utilizado
Estas estimaciones son válidas para ciclos de conducción estándares, tanto para los
vehículos de combustión interna como vehículos eléctricos. Sin embargo, cualquiera de los
dos tipos puede ser conducido bajo cualquier otro ciclo. Se ha demostrado que, a bajas
velocidades, las emisiones de los vehículos de combustión interna aumentan de manera
notable (respecto a la media de emisiones por kilómetro), mientras que las de los vehículos
eléctricos no cambian. Como consecuencia de esto, cuando los vehículos eléctricos se
conducen a velocidades bajas, como en uso urbano, las reducciones conseguidas en las
emisiones pueden ser más notables que las mencionadas aquí, y menos si se conducen a
velocidades relativamente altas.
Por lo tanto, se puede concluir que el uso de vehículos eléctricos facilitaría una
disminución de la contaminación, al menos localmente. Parece lógico pensar que lo que se
lograría, en cierto modo, sería desplazar el foco contaminante desde la ciudad, donde las
emisiones serían mínimas, si no nulas, a zonas más alejadas del ambiente urbano, donde se
encuentran ubicadas las centrales generadoras de electricidad. Que la reducción de la
contaminación se manifieste también globalmente, depende por tanto de las emisiones
producidas por esas plantas de generación, e, indiscutiblemente, será más beneficioso para
el medio ambiente, que el mix de centrales incluya más plantas respetuosas con él.
Por otra parte, el hecho de que las responsables de las emisiones sean únicamente
las plantas de generación eléctrica (en el caso de vehículos eléctricos), trae consigo la
ventaja de un mejor control de las emisiones, al tratarse de instalaciones grandes y
estacionarias cuya localización es conocida, al contrario que la gran cantidad de coches que
circulan por las carreteras, que constituye una fuente móvil difícil de vigilar. Además,
debemos tener en cuenta que, a pesar de que las emisiones de los vehículos de combustión
interna estén reguladas, con el paso del tiempo, los ingenios que las controlan pierden
precisión y no pueden garantizar el cumplir los estándares. No obstante, el aumento
progresivo del uso de vehículos eléctricos conlleva una cantidad cada vez mayor de baterías
gastadas, por lo que se hace necesario el desarrollo de tecnologías para su tratamiento y
reciclaje, con el fin de evitar problemas medioambientales.
Los vehículos eléctricos son los únicos vehículos prácticos que, por ahora, no
producen emisiones mientras funcionan, es decir, son los únicos que pueden considerarse
como vehículos de emisión cero, ZEV.
La CARB ha hecho una apuesta fuerte por ellos, ocupando una parte importante de
su normativa.
Sin embargo, el futuro de los vehículos eléctricos no parece tan cercano como se
quisiera y el motor térmico, al que todavía se le augura una larga vida, debe ser estudiado y
acondicionado para cumplir con las exigencias sociales de la actualidad.
La ley aprobada establece que el 2 por 100 de los vehículos con peso menor de
3750 lb (1700 kg) vendidos por cada fabricante en el estado de California en 1998 deben
ser de emisión cero (ZEV); y esta proporción debe incrementarse al 5 por 100 para el 2001,
y al 10 por 100 para el 2003.
El vehículo eléctrico constituye, ciertamente, una opción que podría gozar de gran
aceptación popular en un futuro. Si se logra equiparar sus prestaciones a las de un vehículo
de combustión interna, su casi total inocuidad para con el medio ambiente, acompañada de
las esperadas ayudas legales, conseguirá sin duda desbancar al automóvil térmico. Sin
embargo, no parece, al menos por el momento, que ese futuro esté muy cercano, por lo
que los fabricantes dedican la mayor parte de sus esfuerzos a la mejora de las tecnologías de
combustión existentes. El objetivo básico que se persigue es el mejor aprovechamiento de
la energía térmica, o sea, la reducción significativa del consumo en los motores, lo que
conlleva, a su vez, una reducción de las emisiones.
Los campos de actuación a los que se debe prestar atención para aprovechar mejor
la energía térmica están señalados en la Tabla 6.
- mejorar la aerodinámica
- aumentar la seguridad
- menores pérdidas
Rodadura
- menor nivel de ruidos
emisiones. Sin embargo, siguen utilizando los mismos combustibles, por lo tanto, lo único
que consiguen es retrasar, en el tiempo, y amortiguar los efectos negativos. Los
combustibles alternativos se presentan como una posible solución. Los más viables son:
Gas natural
Reduce de forma aparente ciertas emisiones (CO2, NOx y otros) pero requiere la
incorporación de depósitos especiales. Además de su alto precio, presenta el inconveniente
de ofrecer una cuarta parte del contenido energético del petróleo, lo cual puede reducir las
prestaciones. Actualmente, se utiliza casi exclusivamente en vehículos de flota, con acceso a
estaciones centrales de recarga.
Etanol (Biomasa)
Su producción es muy cara y consume casi tanto combustible como el que puede
generar. Se debe utilizar mezclado con otros combustibles. La ventaja principal es que
resulta, hasta cierto punto, abundante (se pueden utilizar tierras en desuso) y no es
petróleo.
Metanol
Se consigue, principalmente, del gas natural y ofrece un alto índice de octano, que
permite aumentar el rendimiento del motor. Sin embargo, emite ciertos ácidos perjudiciales
para la atmósfera y el motor. Produce menor energía que el petróleo, por lo que requiere
depósitos mayores y, además, algunas modificaciones en el motor.
Hidrógeno
Como ventajas ofrece una gran energía y una combustión relativamente limpia.
Además, es abundante y no tiene nada que ver con el petróleo. Sin embargo, resulta caro
(requiere mucha energía para su producción), es altamente inflamable, poco denso y su
transporte y utilización resultan complicados.
1.4.- Estructura
2.1.- Introducción
Por ejemplo, la California's Air Resources Board (CARB), define a los vehículo de
emisión cero como aquellos que no emiten contaminantes por su tubo de escape, no
emiten contaminantes por evaporación, no cuentan con sistemas de control de emisiones
que se deterioran con el tiempo y no conllevan emisiones provenientes del refinado de
gasolina.
Las baterías han tenido altos costes de fabricación, peso, tiempo de recarga, y escasa
vida útil y autonomía, lo que ha limitado la adopción masiva de vehículos eléctricos de
batería. Los adelantos tecnológicos actuales en baterías han resuelto algunos de estos
problemas.
En el vehículo eléctrico, el 46% de la energía liberada por las baterías sirve para
mover el vehículo, lo que indica una eficiencia entre 10 – 30% superior de éste, respecto al
vehículo convencional con motor de combustión.
Existen varios motivos por los que este tipo de vehículos no son viables para su
comercialización: no pueden operar por un tiempo prolongado en ausencia de luz solar y
tienen un precio muy elevado, pero pese a estas limitaciones, la tecnología solar podría
contribuir al desarrollo de coches eléctricos de batería o incluso híbridos que utilicen en
parte energía eléctrica y en parte gasolina dependiendo de las circunstancias. En estos
casos, los paneles solares ocuparían mucho menor espacio y ayudarían a alargar la
autonomía de la batería, respetando el diseño y la habitabilidad interior del vehículo.
La recarga del coche se hará, una vez esté el mercado desarrollado, en gasolineras
adaptadas para suministrar aire comprimido. En 2 ó 3 minutos el coche estará listo.
2.3.3.1.- Biocombustibles
2.3.3.1.1.- Biodiésel
2.3.3.1.2.- Bioetanol
gasolina, según la norma de calidad europea EN 228, sin necesidad de modificación alguna
del motor. Algunos fabricantes establecen el máximo de bioetanol en el 5 %, mientras otros
garantizan el buen funcionamiento del motor hasta el 10 %. Dos mezclas comunes son
E10 y E85, que contienen el etanol al 10 % y al 85 %, respectivamente. Si se emplea al
100 % en motores de encendido provocado, se deben realizar ciertas modificaciones en los
motores, además de que sus filtros deben ser cambiados con más frecuencia, ya que el
bioetanol favorece la liberación de los depósitos de partículas del tanque y de los conductos
del combustible.
Las emisiones que produce son muy similares a las del gas natural, pero si tenemos
en cuenta el ciclo de vida completo son inferiores, ya que es un combustible renovable, que
además hace que el CH4 (un potente gas de efecto invernadero) de los vertederos y plantas
depuradoras de aguas residuales no escape a la atmósfera.
Las emisiones del GtL son similares a las de la gasolina, mientras que las que
presenta el BtL son inferiores.
Estos carburantes cuentan con cuatro ventajas respecto al gas licuado de petróleo
(LPG) y el gas natural comprimido (CNG):
o Bio-hidrógeno
o Bio-dimetl-eter (bio-DME)
o Biometanol
El gas natural vehicular es el mismo gas que se utiliza en las casas, en las
termoeléctricas, en la industria en general. Lo que lo diferencia de los otros usos, son las
presiones, que en el caso del uso vehicular se comprime entre 200 y 250 bar, según la
normativa de cada país, con el objetivo de almacenar la máxima cantidad de volumen en los
cilindros dispuestos para tal fin, y de esta, forma brindar la mayor autonomía posible al
vehículo. Este gas natural es principalmente metano, que al tener un alto índice de
hidrógeno por carbono produce menos CO2 por unidad de energía entregada, en
comparación con otros hidrocarburos más pesados (con más átomos de carbono y un
menor ratio H/C).
El GNL tiene una densidad energética menor que el diésel o que la gasolina. Se
necesitan 1,7 litros (en promedio) de GNL para obtener la misma energía que se obtendría
de 1 litro de diésel, y 1,55 litros de GNL para 1 litro de gasolina.
El gas licuado del petróleo (GLP) es una mezcla de propano (C3H8) y de butano
(C4H10), de distinta proporción en función del país y del tipo de vehículo. Se extrae a partir
de procesos de refinado, y de los yacimientos de gas natural húmedo. En condiciones
normales de presión, está en estado gaseoso, pero se licua al someterlo a una presión de
unos 10 bar.
Los vehículos a GLP son muy similares a sus equivalentes de gasolina, radicando las
diferencias fundamentales en los sistemas de almacenamiento (en estado líquido) y
alimentación de combustible al motor (combustión en el motor en estado gaseoso).
Estos vehículos logran la fuerza motriz gracias a la ignición del hidrógeno dentro de
la cámara de combustión. El motor de combustión interna de hidrógeno se parece mucho a
los motores de gasolina, exceptuando que el combustible es gaseoso en lugar de líquido.
A partir del análisis de las prestaciones de los vehículos que cuentan con estas
tecnologías, es posible extraer algunas conclusiones. [9]
H2-
Eléctrico PC Híbrido Biodiésel Bioetanol Biometano GN GLP
MCI
Turismos 86,67 422,58 636,6 960 888,86 1000 337,5 546,67 255
Vehículos
75 450
industriales
Tabla 7: Resumen de la autonomía media obtenida con las diferentes tecnologías estudiadas [km]
Autonomía media
1200
1000
800
Km
600 Turismos
Autobuses
400
Vehículos industriales
200
Figura 10. Comparativa de la autonomía media proporcionada por las diferentes tecnologías (2009)
H2-
Eléctrico PC Híbrido Biodiésel Bioetanol Biometano GN GLP
MCI
Turismos 0,085 0,273 0,435 0,647 0,565 0,833 0,248 0,397 0,153
Vehículos
0,096 0,22
industriales
Tabla 8: Autonomías específicas obtenidas con las diferentes tecnologías estudiadas [km/kg]
Autonomía Específica
H2-MCI
GLP
GN
Biometano
Vehículos industriales
Bioetanol
Autobuses
Biodiesel Turismos
Híbrido
PC
Eléctrico
Figura 11. Comparativa de la autonomía específica proporcionada por las diferentes tecnologías (2009)
Durante más de cien años, la fuente de energía utilizada en los coches ha sido el
motor de combustión interna.
Los diseños actuales del vehículo eléctrico están encaminados a reducir parcial o
totalmente la participación del motor de combustión como fuente de energía en los coches.
Entre los componentes de los coches híbridos y de los coches eléctricos figuran una
batería que almacena la energía, un motor eléctrico de propulsión, un generador, una
trasmisión mecánica y un sistema de control.
lo que cada unidad debe hacer. Usando una batería de condensadores en lugar de
únicamente una batería se reduce significativamente la potencia y los ciclos de energía
requeridos, lo que aumenta la vida de la batería. Sin embargo, la eficiencia se reduce a
medida que el EMC realiza la mayor parte de la transferencia de energía entre las baterías y
los súper condensadores. El EMC se comunica con los controladores del condensador y la
batería, el motor de combustión interna, y los controles de conducción para decidir el par
requerido de los motores de la batería y los condensadores.
La clave para llevar un buen sistema de gestión de la batería, y con ello, alargar la
vida útil de la misma consiste en controlar y cargar las baterías de forma individual. Por
ello, las baterías en los vehículos híbridos deben poseer un sistema inteligente para
mantener y prolongar los ciclos de vida de las mismas ya que en estos vehículos tienen
muchos ciclos de carga y descarga.
B: Batería
E: ICE
F: Tanque de Fuel
G: Generador
M: Motor
P: convertidor de potencia
Conexión eléctrica
Conexión hidráulica
Conexión eléctrica
Debido a la disociación entre el motor y las ruedas motrices, tiene la clara ventaja
de la flexibilidad para localizar el conjunto generador del ICE. Aunque tiene una ventaja
añadida, la simplicidad de su transmisión, necesita tres dispositivos de propulsión, el ICE,
el generador y el motor eléctrico. Por lo tanto, la eficiencia del vehículo híbrido eléctrico en
serie es generalmente más baja. Otra desventaja es que todos estos dispositivos de
propulsión deben poseer un tamaño para dar la máxima potencia sostenida, por ejemplo, si
el vehículo está diseñado para subir por una pendiente larga hace que el HEV en serie sea
caro. Por otra parte, si se utiliza para viajes cortos el generador del ICE puede adoptar una
calificación más baja.
Las ventajas que posee este tipo de vehículos son las siguientes:
Los tres componentes del tren de mando deben estar dimensionados para la
potencia máxima de distancias largas, sostenido y conducción a alta velocidad.
Esto es debido a que las baterías se agotan rápidamente, dejando al motor de
combustión interna que suministre toda la potencia.
Una ventaja del híbrido paralelo frente al serie es que el híbrido paralelo sólo
necesita dos dispositivos de propulsión, la ICE y el motor eléctrico. Otra ventaja es que el
ICE es más pequeño y el motor eléctrico es también más pequeño y se utiliza para
conseguir mejor rendimiento hasta que la batería se agote. Incluso para operaciones a largo
viaje, sólo el ICE tiene que ser clasificado para la máxima potencia sostenida, mientras que
el motor eléctrico puede ser de media potencia. A continuación, se muestran los diferentes
modos posibles de operación del híbrido en paralelo:
Las ventajas del sistema híbrido mixto con respecto a los demás sistemas son las
siguientes:
Las actuales capacidades de las baterías hacen que los Vehículos Totalmente
Eléctricos (EV) se perciban como de limitada utilidad por su autonomía y tiempos de
recarga necesarios.
En un escenario de alta penetración de EVs, con alta fosilización del mix eléctrico,
se podrían ahorrar un 46 % de emisiones de CO2 en 2050.
Se considera el posible uso de los EVs como almacén energético o como una
importante sinergia con las renovables.
Ha sido bien reconocido que los vehículos eléctricos híbridos (HEV) son mucho
más eficientes y menos contaminantes que los vehículos impulsados por gasolina o diésel.
En la Figura 20 se muestra la comparativa de eficiencia energética entre el vehículo
eléctrico y el convencional. También se ha reconocido que los HEV tienen mayor
rendimiento y la aceptabilidad del usuario es mayor frente a los vehículos eléctricos puros
alimentados con baterías (EVs). Sin embargo, toda la energía utilizada en los HEV con la
estrategia de funcionamiento de la batería de carga continua (CS) proviene de la quema de
combustibles fósiles.
Desde el punto de vista del medio ambiente y de la oferta de energía, los vehículos
eléctricos tienen ciertas ventajas sobre los HEV: en su mayoría cero emisiones dispersas,
independencia de los combustibles fósiles y bajo costo de operación. Por el contrario, las
principales desventajas de los EV son la carga de la batería y el tiempo de carga de la
misma.
PHEVs son un tipo de HVE, que está diseñado y controlado para hacer pleno uso
de la energía de la batería para cierta distancia de EV puros, o la operación de carga-
agotamiento (CD). PHEV tiene las ventajas de ambos, del EV puro y del CS HEV. Cuando
la batería se encuentra en estado de alta energía después de cargar durante toda la noche, el
Las ventajas típicas de los vehículos eléctricos puros sobre los híbridos en modo CS
son las siguientes:
Yimin Gao y Merhrdad Ehsani en la referencia [15] nos muestran las estrategias de
control en toda la gama de vehículos eléctricos y los agotamientos de la carga en las
baterías. La simulación realizada se llevó a cabo en un turismo. Los resultados obtenidos
indican que una cantidad significativa de combustible puede ser desplazado por la energía
eléctrica en la conducción urbana típica.
Para modelar la dinámica del vehículo se recurre a las expresiones de la física que
rigen los principales fenómenos a los que están sometidos cualquier vehículo de ruedas.
Estos fenómenos son: el rozamiento de las ruedas, la resistencia aerodinámica, las inercias y
la resistencia a la ascensión debida al peso.
Esta parte de la fuerza es debido a la fricción del cuerpo del vehículo por moverse a
través del aire. Esencialmente, es debido a dos componentes: la forma del vehículo y el
coeficiente aerodinámico.
En el movimiento del vehículo, se generan dos zonas de presión: una que se opone
al movimiento por empuje (alta presión en el frente) y otra que tira de él hacia atrás (baja
presión en la parte trasera) como se muestra en la Figura 22.
A pesar de todo, es necesario considerar las fuerzas necesarias para acelerar las
piezas mecánicas que rotan. Esta es una aceleración rotacional debido no especialmente al
momento de inercia del motor, sino a sus altas velocidades angulares.
Y, por lo tanto,
Fta y Fwa serán negativas si el vehículo decelera y Fhc será negativa si se desciende una
cuesta. Por lo tanto, la fuerza total puede llegar a ser negativa o de sentido opuesto, si
alguna o varias de estas fuerzas prevalece en sentido contrario a todas las demás.
2.6.1.- Introducción
VOLTAJE: Por definición, una batería se compone de dos o más celdas con
cables. Una célula de tipo plomo – ácido produce aproximadamente 2,1 voltios.
Por lo tanto, una batería de plomo – ácido de tres celdas produce 6,3 V
(6,3 = 2,1 * 3). Para una batería con tapones de llenado, el número de células
puede determinarse contando el número de vasos de relleno. La tensión nominal
es la de una batería completamente cargada, y su tensión disminuirá a medida
que la batería se descargue.
Auto-
Tensión por Duración
Energía/peso Tiempo descarga por
Tipo de batería elemento (número de
(Wh/kg) de carga mes (% del
(V) recargas)
total)
Ni – Cd 48 – 80 1,25 500 10 – 14 30
Los valores de los parámetros del circuito equivalente no son constantes. El voltaje
del circuito abierto es el primer parámetro que hay que estabilizar ya que éste cambia con el
estado de carga de la batería.
La variación del voltaje en función del estado de carga de la batería en el caso de las
baterías ión-litio se muestra en la siguiente figura, Figura 27. [19]
Hay que tener en mente que ésta es la carga que se elimina de los platos de la
batería y no la carga que realmente es suministrada por la batería al vehículo eléctrico. La
carga que es suministrada es dada por la siguiente fórmula:
Otra de las situaciones a las que están sujetas es la autodescarga, que consiste en la
pérdida de características de la batería cuando se encuentra en circuito abierto (fuera de uso
o almacenada). Las causas que afectan a la autodescarga son: el diseño de la celda, el
sistema electroquímico, la temperatura (a mayor temperatura, mayor autodescarga), el
tiempo de almacenamiento y las impurezas.
Los dos fenómenos anteriores están más o menos presentes en otras tecnologías de
baterías, pero un fenómeno específico de las baterías de plomo – ácido es la carga
superficial. Esta carga superficial es debida a que las reacciones que producen o almacenan
energía tienen lugar en la interfaz entre el electrodo y el electrolito. Por ello, durante la
carga, ésta se acumula en la superficie de los electrodos y debe difundir dentro de ellos,
pero este proceso de difusión es lento.
Por ello, es importante realizar una carga adecuada de las baterías. Se pueden cargar
de diferentes formas: a tensión constante, a corriente constante o con un sistema de
regulación automática. En la Figura 28 se muestra la evolución de corriente y tensión para
cada tipo.
También hay sistemas de carga que combinan los dos métodos anteriores. Uno de
los más extendidos es cargar a corriente constante hasta alcanzar la tensión de gasificación.
A partir de ese momento, la carga pasa a ser a tensión constante [20].
Como regla general para la vida de las baterías, podemos decir que por cada 10 ºC
de aumento de la temperatura ambiente por encima de la de referencia, la vida útil se
reduce a la mitad. Por ejemplo, una batería de cinco años de duración a 25ºC, solo durará
30 meses si la temperatura en el ambiente es de 35ºC. Por ejemplo, en el caso de las
baterías de plomo – ácido, si se utilizan con inferiores a 27 °C (80 °F), su capacidad
promedio será menor. En el caso de utilizarlas a -18 °C (0 °F) la batería tendrá el 50 % de
su capacidad y a 27 °C (80 °F) tendrá el 100 % de su capacidad. Cuando las baterías se
utilicen en temperaturas superiores a 27 °C (80 °F), su capacidad promedio será superior a
la estimada pero su vida útil será más corta, siendo ésta difícil de predecir, ya que puede
variar según la aplicación, la frecuencia de uso y el nivel de mantenimiento.
100
(%) Capacidad
80
60
40
20
0
-20 0 20 40
Temperatura (ºC)
Capacidad (%)
100
80 Ácido Plomo
60 Níquel Cadmio
40 Ión Litio
20
0
-20 -10 0 10 20 30 40
Temperatura (ºC)
Desde la aparición del primer Sistema de Información Geográfica (SIG) en los años
sesenta, esta técnica ha pasado por diversas fases de desarrollo, alcanzando hoy en día áreas
tan diversas como el geomarketing o la navegación por redes de transporte asistida por
satélite. Sin embargo, actualmente se establece un objetivo fundamental para los SIG: la
solución de problemas espaciales complejos, de igual manera que en el entorno general de
los Sistemas de Información (SI) se persigue resolver problemas, por lo general, de índole
no espacial, o bien contestar rápidamente a las preguntas planteadas al sistema.
A través de las distintas etapas en la evolución de los SIG, se han alcanzado varios
logros que hacen de estos sistemas una potente herramienta y, más aún, una creciente
tecnología para todos aquellos sectores que requieren la gestión de información espacial de
manera rápida y eficaz. Si bien los actuales SIG son capaces de gestionar, almacenar,
representar e incluso desempeñar papeles fundamentales en aplicaciones avanzadas como
las antes citadas, ciertas cuestiones referentes al análisis espacial aún no se encuentran
resueltas por completo en su ámbito.
Los SIG pueden entenderse como una caja de experimentación lo que permite al
analista o al gestor territorial trabajar o plantearse diferentes escenarios virtuales de una
determinada región: por una parte, los que se producirían con la ejecución de ciertas
políticas, o los que ocurrirían siguiendo determinadas tendencias. Todo esto hace de los
SIG sean una potente herramienta de planificación cuando se dispone de una base de datos
suficiente amplia para los fines que se plantean.
Desde que la idea de los SIG fue concebida hace unos cuarenta años, han
evolucionado de igual manera que su conceptualización. Así, tenemos que definiciones
como base de datos georreferenciados o base de datos espaciales, hoy son obsoletas.
Actualmente, son frecuentes definiciones que resaltan la operatividad de los SIG, tales
como:
1. Entrada de datos
3. Manipulación y análisis
4. Salida
Sin embargo creemos que, tal y como sostienen Burrough y Bouillé, un SIG debe
verse también como un modelo del mundo real, por lo que se podría definir como: “Modelo
informatizado del mundo real, en un sistema de referencia ligado a la Tierra para satisfacer unas
necesidades de información concretas”.
Un SIG está integrado por cinco componentes principales: hardware del ordenador
(parte física del sistema), software del ordenador (programas de aplicación), datos, usuarios
y métodos.
USUARIO
SOFTWARE DATOS
SIG
HARDWARE PROCEDIMIETO
El primer elemento representa la parte física donde se asienta el SIG, y suele estar
representado por alguna plataforma de la computadora, que pueden variar desde servidores
(computador central) a computadores desktop (escritorio) o Laptop (portátil) que se
utilizan en las configuraciones de red o desconectado. Asimismo, es requerido un conjunto
de periféricos, como tabletas digitalizadoras, plotters, lectores ráster (scanner), así como
unidades de almacenamiento y procesamiento de datos, para poder desarrollar la potencia
operativa de los SIG.
desarrollar en un SIG, además de ser el aspecto que requiere un mayor esfuerzo para su
implementación en un proyecto SIG.
En este sentido, los usuarios de estos sistemas pueden optar por alguna de estas dos
soluciones ante la carencia de datos espaciales:
A pesar de estas dos posibilidades, los inconvenientes para lograr una base de datos
capaz de alimentar un proyecto de planificación territorial son múltiples, e incluso
insalvables a costes bajos o medios.
a) entrada de información
b) gestión de datos
d) salida de datos
Con estos cuatro conjuntos de funciones se puede representar todas las operaciones
posibles en un SIG.
Esta etapa es fundamental para disponer de una base de datos potente, operativa,
libre de errores y versátil, lo que permite posteriormente un adecuado funcionamiento del
SIG.
que implica que esta información debe ser homogeneizada y corregida para poder ser
introducida en el sistema.
todos los objetos geográficos con un conjunto particular de atributos, ya bien sea una
localización (espaciales) o bien el valor de una variable (temáticos).
Dado el nombre(s) de uno o varios objetos espaciales se obtiene una tabla y/o un
mapa que muestren dichos objetos espaciales y sus características temáticas y/o espaciales
(Figura 35). Este ejemplo muestra dicha búsqueda. Se quiere obtener todos los objetos
geográficos (secciones censales de una ciudad) cuyo nombre empiece por el número 1,
porque todos pertenecen al distrito 1. La búsqueda proporciona dos resultados: un mapa,
donde se puede observar la localización y una tabla en la cual da los valores de las variables
temáticas asociadas a cada sección censal.
Este tipo de búsqueda consiste en determinar qué valor temático aparece en una
localización, en el caso de un SIG vectorial, en uno o varios objetos geográficos localizados
en el mapa. También existen varias formas distintas de búsqueda espacial.
todas las secciones censales que tienen alguna extensión de terreno dentro de dicho
dominio, cumpliendo las dos condiciones aritméticas simultáneamente.
casas que, en el estrato ‘edificaciones aisladas’, están a menos de un kilometro de los ríos
del tipo seleccionado anteriormente.
Figura 38. Búsqueda espacial empleando dos estratos mediante condición geométrica
Esta función es el aspecto fundamental de los SIG, donde radica todo su potencial
operativo. Las funciones de transformación y análisis de datos son las que proveen nuevos
datos a partir de los existentes originalmente. Es aquí donde el usuario define los datos que
utilizará y cómo para resolver problemas espaciales determinados, estableciéndose así
soluciones a través del SIG con las operaciones que utilizan los datos espaciales de
diferentes maneras.
El estudio de la ‘intersección de líneas’ busca hallar los lugares donde dos redes o
sistemas de objetos lineales se cruzan y, al mismo tiempo, generar el conjunto de
nuevas aristas y nudos creados por dicha intersección.
Una vez generados los nuevos polígonos con su topología y sus identificadores, es
preciso otorgar valores en cada uno de ellos a la o a las variables temáticas que tienen
asociadas. El método es diferente de acuerdo con el tipo de las variables relacionadas con
los polígonos de origen.
1. Las variables iniciales son, las dos, de tipo nominal/ordinal: en este caso es
posible obtener los valores de la variable temática de los nuevos polígonos
de acuerdo con los siguientes procedimientos:
En este caso son posibles dos tipos diferentes de operaciones: álgebra de mapas y
enmascarado.
Enmascarado de mapas: en este caso, se elige a uno cualquiera de los dos como
mapa máscara que cubre los valores del mapa inferior, a excepción de los polígonos
dotados del valor nulo en el mapa superior.
c) Interpolación areal
exactamente iguales, se producen numerosos polígonos que no existen y que sólo son el
producto de los errores de la digitalización (ver Figura 40). Evidentemente, la operación de
superposición requiere disponer de métodos que eviten su aparición en el resultado final.
Las posibles soluciones para ello son de tres tipos.
En primer lugar, realizar la digitalización de modo que se eviten los errores citados.
Para ello, es imprescindible que cualquier nuevo mapa que se incluya en la base de datos
utilice como elementos espaciales propios los que existen en otros mapas previamente
incluidos en el SIG. Con lo cual, se debe contar con programas que busquen en la base de
datos espaciales los elementos que se van a incorporar y que tomen sus coordenadas (sí ya
han sido digitalizados) de las ya existentes. En muchas ocasiones, esta solución no es
posible, por ejemplo, si se añade a una base de datos mapas tomados de otras fuentes de
información digital no será posible evitar las diferencias en la digitalización.
Finalmente, una tercera posibilidad es eliminar los polígonos ficticios una vez
terminada la fase geométrica – topológica de la superposición, aprovechando para ello su
forma muy especifica que los identifica claramente. También ayuda a su eliminación la
circunstancia habitual de que sean de su tamaño muy reducido. Por ello, basta con buscar
polígonos menores de una extensión superficial prefijada y con la forma característica
(redondeada), se les detecta y a continuación, se les elimina, incluyéndolos en los dos
polígonos contiguos, dividiéndolos por la mitad.
Este problema es sólo uno, aunque muy importante, de los errores que pueden
afectar a las operaciones de superposición y en general a todas las tareas desarrolladas por
un SIG.
usual es el área de influencia de una línea eléctrica, es decir, la anchura que hay que dejar al
lado de dicha línea para que no se produzcan averías.
En un SIG existen diversas formas de salida de datos, las cuales dependen de los
requerimientos del usuario, las más frecuentes son: mapas analógicos, tablas de valores,
gráficos, representaciones tridimensionales, simulaciones de vuelo sobre ciertas zonas, etc.;
con estas salidas podemos representar la información contenida en la base de datos o bien
mostrar el resultado de determinadas aplicaciones.
espaciales de un mapa de acuerdo con un tema determinado; así, cada capa almacena un
tipo particular y homogéneo de objetos espaciales.
La estructuración de la información del mundo real en capas con lleva cierto nivel
de dificultad para la representación de los datos espaciales como objetos espaciales
pertenecientes a una capa, ya que cada objeto espacial está representado por información
acerca de su posición (localización absoluta X, Y), relaciones topológicas (localización
relativa en función de otros objetos espaciales y en sus relaciones) y atributos de los objetos
(características temáticas de cada uno). Siendo los aspectos topológicos y espaciales
relacionados con cada objeto espacial lo que diferencia a los SIG de otros sistemas de
información.
En este sentido, Bosque [25] indica que un SIG debe de estar en condiciones de
manejar tanto las características espaciales de los objetos geográficos (la geometría o
localización absoluta y la topología o relaciones cualitativas entre ellos) como los aspectos
temáticos asociados a los objetos o unidades de observación. Pudiendo así establecer con
los tres tipos de información que caracterizan a cada objeto espacial las posteriores
capacidades de análisis espacial que puede desarrollar un SIG.
Con relación a la representación de los hechos del mundo real en los SIG, la escala
desempeña un papel fundamental en la selección de los objetos espaciales, ya que
elementos espaciales de la realidad que en una escala pueden representarse por áreas o
líneas, en una menos detallada pueden representarse por puntos, y viceversa. De esta
manera, la representación de los elementos del mundo real en objetos espaciales se realiza
como una abstracción de cada elemento espacial del mundo real.
luego, dependiendo del modelo de datos que utilicemos, representarla a través de puntos,
líneas, áreas o celdas.
Las relaciones que se pueden establecer entre los tres objetos espaciales básicos
fueron definidas por Gatrell entre pares de objetos.
- Cruza - Interseca
Líneas
- Une - Es límite
- Superpuesto por
Áreas
- Es adyacente a
En la tabla anterior, Tabla 11, podemos ver las relaciones posibles de establecer
entre cada par de objetos espaciales del modelo vectorial; sin embargo, en caso de disponer
de un SIG ráster, es posible llevar a cabo operaciones de superposición a partir de las
celdas de cada capa temática que produzcan un resultado similar, desde un punto de vista
operativo, al del modelo vectorial aquí comentado.
En el modelo de datos vectoriales, las áreas están descritas por el conjunto de líneas
de frontera que las delimitan; dichas líneas se representan por un conjunto de coordenadas
correspondientes a sus vértices, mientras que los puntos se representan por una
coordenada (X, Y).
En este modelo de datos, los elementos del mundo real están representados por un
conjunto de celdas adyacentes, pero que en realidad no mantienen una relación mutua entre
ellas, ya que las fronteras de las regiones o áreas están definidas implícitamente, y no se
La comparación entre los dos modelos de SIG se refiere sobre todo al volumen de
almacenamiento, calidad de la representación gráfica, exhaustividad de la representación y
facilidad del proceso.
El modelo ráster tiene, por su parte, una organización muy simple de los datos, que
permite realizar con gran facilidad procesos de análisis. Por ejemplo, la superposición de
mapas es muy sencilla y fácil de programar mediante operaciones con matrices. Sus
gráficos, aunque deficientes, se pueden realizar con dispositivos baratos, como una
impresora matricial. Su principal inconveniente es el gran volumen de almacenamiento que
requiere, en especial si es necesario disponer de una representación muy precisa, lo que
exige disminuir el tamaño del ‘pixel’ y, en consecuencia, aumentar el número de filas y
columnas del mapa. Por este tipo de cuestiones resulta muy importante el desarrollo de
procedimientos con el ‘run-length’; el cual disminuye el volumen de almacenamiento
necesario. Por último, este modelo no reconoce explícitamente la existencia de objetos
geográficos, y por tanto, en las aplicaciones en que sea esencial su empleo, este modelo
tiene pocas posibilidades de ser utilizado.
Utiliza una estructura de datos muy Posee una estructura de datos muy
simple. compacta.
Tabla 12: Ventajas y desventajas de los SIG vectoriales y los SIG ráster
Las aplicaciones de un SIG se puede clasificar de acuerdo con varios criterios. Por
un lado, las tecnologías de diverso tipo que interactúan con los SIG y a las cuales éstos
aportan procedimientos e ideas, como son, por ejemplo:
Teledetección
Topografía y geodesia
Las más importantes son la clasificación por las áreas o temas donde estos
dispositivos han encontrado un campo de uso práctico y cotidiano, que serían los
siguientes, ordenados de mayor a menor desarrollo:
1. Gestión de instalaciones
2. Catastro
3. Planificación urbana
4. Transporte
5. Marketing geográfico
A. GESTIÓN DE INSTALACIONES
y cartografía a escalas grandes (Ej, 1:500) y en la que se incluyen las líneas de la topografía,
los edificios, y, sobre todo, las redes de distribución a mantener.
Son aplicaciones que permiten la cartografía del lugar donde está el automóvil y la
determinación de la ruta óptima al punto de destino. A partir de un SIG que incluya la
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 2 Conceptos básicos 105
E. GEODEMOGRAFÍA Y MARKETING
Localizar comercios
Para ello, se usan divisiones del espacio muy precisas, por ejemplo, bloques de
edificios o manzanas urbanas, o al menos secciones censales. En ellas se recoge
información (procedente habitualmente del Censo o Padrón) de las personas que allí
residen, aspectos demográficos (edades, sexo), económicos (profesiones y niveles de
ingreso), sociales (costumbres y nivel cultural, etc.)
En la actualidad existe una gran oferta de programas sobre SIG, que cubre
cualquier posibilidad de trabajo.
Por último, también existen algunos programas de uso más habitual en la enseñanza
de los SIG, aquí destaca con claridad el programa IDRISI, de la Universidad Clark (USA),
por su facilidad de uso, coste reducido, puesta al día anual, etc..
Carga semi-rápida
Carga rápida
Esta solución es la que, desde el punto de vista del cliente, se asemeja a sus
hábitos actuales de repostaje con un vehículo de combustión.
Los modos de carga tienen que ver con el nivel de comunicación entre el vehículo
eléctrico y la infraestructura de recarga (y por consiguiente, la red eléctrica), y el control que
se puede tener del proceso de carga, para programarla, ver el estado, pararla, reanudarla, o
incluso volcar electricidad a la red.
Modo 1, sin comunicación con la red. Sería el que se aplica a una toma de
corriente convencional con conector schuko.
Ambientes
Potencia
públicos o
monofasica/ Recarga
Modo Corriente (A) Carga Recarga en privados
trifasica controlada
abiertos a
(kw)
terceros
Max. 16 A por
Modo 1 3,7 /11 Lenta AC No No
fase
Max. 32 A por
Modo 2 7/ 22 Lenta AC Si No
fase
Figura 46. Gráfico resumen de los 4 modos de carga según IEC 61851
Con el fin de solucionar el problema de autonomía que tienen los VE, y permitir
que los usuarios puedan utilizar estos con la mayor libertad posible, nacen los distintos
tipos de infraestructuras de recarga del VE, para poner al alcance de todos los usuarios una
energía accesible y en cualquier lugar y momento.
Existen tres tipos de toma de corriente. En la siguiente tabla, Tabla 14, se muestran
las características [27]:
70 A (monofásica)
Intensidad máxima (A) 32 32
63 A (trifásica)
La toma de corriente Tipo 3 posee obturadores que evitan el contacto directo con
las partes en tensión, bajo un grado de protección IPXXD, denominada ‘child protection’.
Tabla 15: Tipos de toma de corriente estándar según la norma IEC 62196
Donde, L1: Línea 1; L2: Línea 2; L3: Línea 3; N: Neutro; PE: Puesta a tierra; CP:
protección de control; CS: control de seguridad; y PP: proximidad.
El estándar IEC 62196-3 adoptado por la Comisión Europea tiene dos conectores
posibles:
Los dos conectores son compatibles con la base de enchufe instalada en el vehículo.
Vía pública
En este caso, el poste de recarga se situara en la acera con el fin de dar suministro a
los vehículos eléctricos estacionados en las plazas de aparcamiento reservadas para esta
actividad. Dentro de este grupo podemos encontrar dos modalidades:
Suelen tener entre tres y seis tomas de corrientes, y debe disponer al igual
que los postes de recarga, de un sistema de protección ante riesgos
eléctricos y actos vandálicos que impida su acceso sin identificación.
Entorno controlado
- N puntos de recarga.
La autentificación de los usuarios que quieren hacer uso del sistema de recarga,
viene dada por un sistema de lectura de tarjetas RFID. Una vez realizada la autenticación,
se electrificará el punto de recarga elegido o asignado.
Interior
1. Monopunto
2. Sistemas multipunto
3.1.- Introducción
En la exposición que realiza Chan en la referencia [12], muestra que las principales
claves tecnológicas del EVs, HEVs y FCVs son:
4) La batería y ultracondensadores
Hay tres tipos principales de motores eléctricos que son adecuados para
aplicaciones de EV, HEV y FCV: motor síncrono o motores sin escobillas, los motores de
inducción; y, los motores de reluctancia. Los requisitos típicos de motor son: alto par y
potencia, amplia gama de velocidades incluyendo par constante y operaciones de potencia
constante, alta eficiencia en todo el rango de velocidades, de alta fiabilidad y robustez, todo
ello, a un costo razonable.
Prochazka [29] en su artículo nos define que la mejor opción es la del motor con
transmisión ya que si se ubica el motor directamente en las ruedas sin transmisión, éste no
posee ninguna ventaja. Por ello, se deja está ubicación para las bicicletas eléctricas.
Si se comparan motores asíncronos frente los síncronos, los primeros poseen varias
ventajas frente a los síncronos:
Las pérdidas son más bajas para accionamiento con pequeño par y alta
velocidad debido a la desexcitación del motor (al tener menor par, las pérdidas
de hierro son inferiores)
3.2.2.- Convertidores
Para la realización del modelo de la batería, se debe cumplir con los siguientes
requisitos:
El esfuerzo para parametrizar el modelo debe ser bajo tanto para las
mediciones, así como para el posterior análisis necesario
El coche eléctrico, cada día, está más en auge con lo que va a conllevar tanto un
cambio social como técnico – económico, por ello, las empresas eléctricas deben estar
preparadas para poder hacer frente a la carga masiva de vehículos eléctricos utilizando
horas – valle, implementando nuevas políticas, ampliando las instalaciones, construyendo
puntos de recarga, etc.
Figura 54. El diagrama de bloques funcional de la parte eléctrica del HEV propuesto por Habib Ullah
En la referencia [33], se muestra el trabajo realizado por Ibrahim el cual realiza una
simulación en la cual analiza el fallo de cortocircuito en el vehículo eléctrico híbrido
tipo – serie – paralelo con la finalidad de observar cómo influye en el rendimiento. Para la
realización de la simulación utiliza los siguientes softwares: Matlab/Simulink,
SimPowerSystems y SimDriveline. Demostró una posible reacción del vehículo a fallo de
cortocircuito durante la aceleración y desaceleración.
Cheng nos muestra en la referencia [34] una manera de optimizar la operación del
motor de combustión interna minimizando las emisiones y maximizando el ahorro
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 3 Estado del Arte 125
Jalalifar en la referencia [35] nos muestra el uso de una técnica de linealización por
realimentación de entrada – salida combinada con un observador backstepping de
adaptación en el marco de referencia del estator del motor de inducción utilizando en los
vehículos eléctricos híbridos en serie. Una de las mejores ventajas de este método de
control es la eliminación del sensor de flujo y disminuye el costo de control, además, el
sistema de control es robusto respecto a las variaciones de las resistencias y el par de carga
externo.
En los últimos años, se han propuesto varias estrategias para la gestión de la energía
del vehículo híbrido eléctrico. Estas estrategias de optimizar la economía de combustible se
realizan mediante la resolución de un problema en el que el consumo de combustible
aparece explícitamente en la función de coste. Sin embargo, en lugar de optimizar
explícitamente el consumo de combustible, el objetivo de Di Cairano y compañía [37] es la
mejora de la eficiencia del funcionamiento del motor. Dado que en la configuración del
vehículo híbrido eléctrico en serie, el motor está desacoplado mecánicamente de las ruedas;
para una solicitud de potencia, dado el punto de funcionamiento del motor en estado
estacionario, la velocidad del motor y el par del mismo pueden ser arbitrariamente elegidos
dentro de una curva. Por lo tanto, una curva de puntos de funcionamiento óptimo como la
función de la potencia solicitada del generador puede ser calculada, y el punto de
funcionamiento del motor puede estar controlado para localizarse en la curva de estado
estacionario. Sin embargo, durante los transitorios de la dinámica del motor, el punto de
trabajo se encuentra lejos de la curva óptima. Por ello, proponen el uso de la batería para
“suavizar” los transitorios debidos a que el punto de operación del motor se desliza en una
estrecha banda a lo largo de la curva óptima de eficiencia. Para ello, implementan un
algoritmo de control predictivo del modelo (MPC). Los resultados obtenidos al realizar el
test UDDS es que se mejora la economía de combustible.
Millner, en la referencia [40], propone tres mejoras en PHEV las cuales pueden
mejorar la economía. (Ver Tabla 16)
Mejoras Ventajas
Uso de la batería del vehículo mientras está Compensa los picos cortos en planta a escala
estacionado comercial, la demanda eléctrica para reducir los
cargos por demanda, puede proporcionar
ingresos adicionales para pagar la batería
El ciclo de vida de la batería debe ser Se presenta un modelo de batería para realizar
maximizada para evitar los costos de reemplazo las baterías de iones de litio y se utiliza para
de ésta confirmar que las estrategias anteriores son
compatibles con la larga duración de la misma
Tabla 16: Mejoras propuestas por Millner [40] para los PHEV
Knoedler nos muestra, en la referencia [42], las soluciones para lograr la máxima
eficiencia energética y minimizar las pérdidas de energía en todas las condiciones de
contorno de tráfico que varían dinámicamente. Además, muestra que las tecnologías
desarrolladas se pueden adaptar a HEV y vehículos convencionales.
En las ciudades, con una concentración cada vez mayor de vehículos, el aire y la
contaminación acústica, todos ellos representan un problema permanente. Un pequeño
coche eléctrico no produce gases y disminuye el ruido, lo que le hace ideal para los centros
de las ciudades. Además, al poseer un tamaño reducido, es más fácil encontrar una plaza de
aparcamiento y maniobrar en espacios reducidos. Dado que la mayoría de las personas que
trabajan en el centro de la ciudad, viven cerca de ésta, la carga de las baterías puede en
muchos casos estar limitada con la finalidad de que los propietarios recarguen durante la
noche. También, si el vehículo posee tres ruedas se reduce el coste de producción y
mantenimiento. El único inconveniente es la estabilidad en las curvas si se compara con los
vehículos de cuatro ruedas. A partir de todas estas hipótesis, Vasiljevic en la referencia [43]
nos muestra el modelo dinámico y la simulación de un coche eléctrico de tres ruedas. En la
Figura 58 se muestra el coche experimental.
Figura 59. Estructura típica de un vehículo eléctrico con accionamiento de las cuatro ruedas independientes
El principio de esta idea es el uso de un alternador que está conectado al eje de las
ruedas traseras. A medida que el eje de rotación del alternador genera la energía eléctrica,
éste se utiliza para cargar las baterías. La disposición de los elementos se puede ver en la
Figura 61.
Figura 61. Disposición del sistema para cargar coches eléctricos con alternador
Sara Mehar y compañía [46] proponen un algoritmo optimizado para localizar las
estaciones de carga de los vehículos eléctricos. En primer lugar, determinan las
limitaciones, y a continuación, proponen un algoritmo de optimización, al cual denominan
OLoCs (Optimized Location Scheme for electric charging stations), es decir, un algoritmo
que permite la localización óptima de la estaciones de carga.
Figura 62. Localización óptima de las estaciones de carga realizado por Sara Mehar [46]
ubicación de la estación de carga. Se demostró que el modelo y el método tienen una gran
corrección y eficacia con la ayuda de un ejemplo de planificación de una estación de carga
para una ciudad.
Figura 63. Cobertura de las particiones y la capa de las estaciones de carga [49]
Por todo esto, Oliveira en la referencia [51] investiga los efectos de los vehículos
eléctricos enchufables y los vehículos híbridos eléctricos enfuchables, la generación de
energía a partir de las mareas y del viento para recargar la batería de los vehículos híbridos y
así introducir la energía renovable en los sistemas de distribución. Además, presenta un
sistema de distribución con la generación de energía renovable. Los resultados mostraron
que la penetración de energía renovable aumenta debido a que los vehículos se consideran
como sistemas de almacenamiento de energía y ayudan en el equilibrio de la energía del
sistema. Para los fines de los servicios auxiliares de V2G (Vehicle-to-grid, traducido del
inglés como "del vehículo a la red”), describe un sistema en el cual la energía puede ser
vendida a la red eléctrica por el conductor de un vehículo todo-eléctrico o híbrido
enchufable, cuando éste es conectado a la red en los momentos que no se use para el
transporte. Alternativamente, cuando las baterías del coche necesiten ser recargadas, el flujo
se invertirá y la electricidad fluirá de la red al vehículo. Los resultados muestran que los
PHEV y PEVS encajan mejor en la regulación del mercado, por lo que es necesaria una
planificación de la capacidad precisa para nuevos estacionamientos y las fuentes
distribuidas.
Kadar en la referencia [52] nos muestra diferentes soluciones para la carga y recarga
de las baterías de los coches eléctricos. Una de ellas es la transmisión de energía inalámbrica
a alta frecuencia entre la cinta de acero de las ruedas de los coches y los electrodos
enterrados debajo de la superficie de la carretera. Otra solución es la carga por medio de un
sistema de energía fotovoltaica el cual puede generar la suficiente energía para el coche.
La nueva conciencia por el medio ambiente se hace notar por el uso de los
vehículos eléctricos, y la generación de electricidad renovable. Por supuesto, la electricidad
generada por energía fotovoltaica es una buena fuente de combustible para el tráfico y
favorable al medio ambiente ya que al generarla no se producen emisiones a la atmosfera.
Hay más opciones de cómo combinar la fuente fotovoltaica con las baterías de los
vehículos eléctricos:
A. Película delgada fotovoltaica pintada en el chasis del coche: otra solución es las
pinturas resistivas de luz en la superficie del coche. Las soluciones de la
tecnología fotovoltaica con película delgada son baratas y pueden ser ’pintadas’
en el cuerpo del coche.
Para mejorar la eficiencia del sistema eléctrico es muy importante que la demanda
eléctrica se desplace hacia las horas de menor consumo; y es ahí, donde la recarga lenta
nocturna del coche eléctrico puede jugar un papel fundamental en el aplanamiento de la
curva de demanda. Además, la recarga del vehículo eléctrico en las horas nocturnas
permitirá optimizar las infraestructuras eléctricas y aprovechar la energía renovable que en
ocasiones se deja de producir por la noche al no haber suficiente demanda para poder
integrarla en el sistema. A largo plazo, el coche eléctrico puede convertirse en un sistema
reversible de almacenamiento distribuido de energía, pudiendo verter de nuevo a la red, en
los momentos de máxima demanda, la energía que se ha almacenado durante la noche. Esta
interacción exigirá la implantación de sistemas de comunicación entre los usuarios y el
operador del sistema.[53]
4.1.- Introducción
En las secciones previas, se han definido las ecuaciones de la dinámica del vehículo
y de la batería. Para ver cuanta distancia puede recorrer un vehículo antes que la batería se
descargue, se va a simular paso a paso el proceso a través de un ciclo de conducción. La
manera de realizar este proceso está representado en el siguiente flujo grama en la Figura
68.
Figura 66. Mensajes en Matlab/Simulink para la elección del ciclo de conducción: Reglado o experimental
En el caso de elegir ciclo reglado, se pueden optar por los siguientes ciclos:
FTP75
Figura 67. Mensajes en Matlab/Simulink para la elección del ciclo de conducción: Ciclos Reglados
Inicio
NO
¿Fin de ciclo?
SI
Actualizar valores de final de
¿Bateria descargada?
ciclo y hacer uno nuevo
NO
SI
Parar
Figura 69. Mensajes en Matlab/Simulink para la elección del ciclo de conducción: Experimental
El siguiente paso es establecer los parámetros del vehículo como la masa, el tipo de
batería, potencia consumida por los accesorios, etc. El programa da la opción de poder
elegir el tipo de batería y a partir de aquí el programa actualiza los valores dependiendo del
tipo de batería elegido.
Radio 1 30 30 Aleatorio
Motor y Caja
Batería Carretera
controlador reductora
Ruedas
b m g
Potencia media
de los Energía para
accesorios=Pac mover el
vehículo
Figura 71. Flujo de energía en un sistema de tracción eléctrico con freno regenerativo
Tras cargar los datos de entrada, se obtiene el valor de la aceleración para cada
instante como diferencia de las velocidades:
La potencia mecánica requerida por motor será la potencia necesaria en las ruedas
condicionada por las pérdidas en la transmisión,
en caso de tracción:
en caso de retención:
La potencia eléctrica requerida por el motor será la mecánica condicionada por las
pérdidas calculadas en la anterior ecuación,
en caso de tracción:
en caso de retención:
En este punto, se debe considerar otros sistemas eléctricos del vehículo, las luces,
indicadores, calentadores, radio, etc. Denominaremos la potencia para los accesorios como
Pac. Se obtiene así la potencia requerida de la batería:
En caso de tracción:
En caso de retención:
En primer lugar, se va a realizar una comparativa entre los tres tipos de baterías
elegidas en la simulación, cuyas especificaciones son las siguientes:
Distancia
75,93 98,44 108,67 125,33 140,14 144,53 141,24 245,74 265,97
recorrida (Km)
Profundidad de
0,9001 0,9002 0,9000 0,9000 0,9001 0,9001 0,9003 0,9000 0,9000
descarga
140
120
FTP75
100
Distancia recorrida (Km)
20 Circulación no fluida
0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
Profundidad de descarga
Figura 72. Autonomía del vehículo eléctrico en función del ciclo de conducción
Potencia
20 760 20 760 20 760
accesorios (W)
Distancia
100,38 91,30 143,08 130,56 250,28 228,90
recorrida (Km)
Profundidad de
0,9000 0,9003 0,9001 0,9000 0,9000 0,9001
descarga
Tabla 20: Variación de la autonomía con respecto a la potencia de los accesorios en un vehículo eléctrico
La masa del vehículo eléctrico en vacio y un conductor ronda los 1605 kg. El peso
máximo admitido, PMA, para el cual está diseñado el vehículo es de 1985 kg. Como se ha
visto, el peso del vehículo tiene una gran influencia en las prestaciones del mismo. De igual
manera, sucede para la autonomía.
Distancia
98,441 81,533 140,142 119,332 245,745 209,958
recorrida (Km)
Profundidad de
0,9002 0,9001 0,9001 0,9000 0,9000 0,9000
descarga
Tabla 21: Variación de la autonomía con respecto a la carga del vehículo en un vehículo eléctrico
Los resultados obtenidos indican que el peso de la carga transportada tiene un gran
impacto en la autonomía del vehículo.
5.1.- Introducción
siempre que el motor tenga un exceso de energía necesaria para la propulsión, sin que el
estado de carga sea superior a su límite de carga completa [56].
La estrategia de control máximo se explica en la Figura 74. En esta figura, las curvas
de potencia máximas de tracción híbrida (motor térmico más el motor eléctrico), la tracción
del motor térmico por sí solo, de tracción del motor eléctrico solo y el frenado regenerativo
se representan frente a la velocidad del vehículo. Las demandas de energía también se
representan en diferentes condiciones, representadas por los puntos A, B, C, y D.
vehículo
Pe: Potencia del motor
1. Eléctrico a baja velocidad. Este modo está activado cuando la velocidad del
vehículo es inferior al mínimo que garantiza un rendimiento aceptable del
motor térmico. En este modo, el motor térmico está apagado y el vehículo
es propulsado únicamente por el motor eléctrico. Por ello, las baterías
necesitan una cierta carga inicial para propulsar el vehículo hasta que el
motor térmico entre en funcionamiento.
Energía Eléctrica
Energía Mecánica
Batería
Inversor
Reductora
Eje tractor/Ruedas
Figura 75. Diagrama del flujo de energía del modo sólo eléctrico
Energía Eléctrica
Energía Mecánica
Batería
Inversor
Reductora
Eje tractor/Ruedas
Energía Eléctrica
Energía Mecánica
Batería
Inversor
Reductora
Eje tractor/Ruedas
Figura 77. Diagrama del flujo de energía del modo recarga de baterías
4. Eléctrico puro. Este modo está activo cuando las baterías se han cargado
por completo y se mantiene activo hasta que las baterías llegan a la carga
mínima, excepto cuando el par exigido es superior al que puede dar el
motor eléctrico y en caso de frenada regenerativa. La única fuente de
energía son las baterías, por lo tanto, la única fuente de potencia es el motor
eléctrico, mientras que el motor térmico permanece apagado.
Energía Eléctrica
Energía Mecánica
Batería
Inversor
Reductora
Eje tractor/Ruedas
Figura 78. Diagrama del flujo de energía del modo frenada regenerativa
En este caso, forzar al motor térmico a dar más potencia para transformar
la parte sobrante en energía eléctrica (con las pérdidas en transmisión,
conversión, almacenamiento y posterior uso que ello conlleva) supone una
merma en la eficiencia, aparte de un desgaste innecesario de ciertos
elementos de la planta propulsora.
Para ello, este modo evalúa si circulando con la marcha más alta el
rendimiento del motor sería igual o superior al máximo que se podría
obtener haciendo trabajar tanto al motor térmico como al eléctrico en su
punto de máxima eficiencia, en cuyo caso el motor eléctrico se apaga y el
vehículo se comporta como un vehículo convencional con motor térmico.
Energía Eléctrica
Energía Mecánica
Batería
Inversor
Reductora
Eje tractor/Ruedas
Figura 79. Diagrama del flujo de energía del modo térmico puro
Si
Velocidad del Si Modo de frenada
vehículo, V mixto
Si Modo eléctrico a
¿V<Vmin?
baja velocidad
No
Potencia de la
batería, Pbat
Si Modo eléctrico
¿Ptc=Pbat?
puro
No
No
SOC de la
batería No Modo térmico puro
¿SOC<SOCmax?
Si
Modo de recarga de
la batería
Figura 80. Organigrama de los modos de control del vehículo eléctrico híbrido en paralelo
CAJA DE
CAMBIOS
CONTROL
DINÁMICA ANÁLISIS
Figura 81. Diagrama de bloques del modelo del vehículo eléctrico híbrido en paralelo
Este margen, en el nivel mínimo de carga, depende del escenario en el que se esté
trabajando. Si el sistema detecta que está circulado en ciudad, éste reserva la energía
necesaria para realizar una aceleración de 0 a 50 km/h. De lo contrario, si detecta que está
en una vía interurbana, reserva la energía necesaria para realizar una aceleración de 80 a
120 km/h.
En primer lugar, se va a realizar una comparativa entre los dos tipos de baterías
elegidas en la simulación, cuyas especificaciones son las siguientes:
Tabla 22: Especificaciones técnicas de las baterías para los vehículos híbridos eléctricos
1,8 Híbrido
1,6 Recarga
Eléctrico
Porcentaje de diferencia (%)
1,4
Frenada regenerativa
1,2
Térmico
1
0,8
0,6
0,4
0,2
0
217
271
325
379
649
703
757
811
865
109
163
433
487
541
595
919
973
1
1027
1081
1135
1189
1243
1297
1351
1405
1459
1513
1567
1621
1675
1729
1783
1837
55
Tiempo (s)
Figura 90. Distribución de modos: porcentaje de diferencia entre la batería de ácido – plomo y la batería de
ión – litio
0,30
0,25
0,20
DoD
0,15
Ácido plomo
0,10
Ión litio
0,05
0,00
100
200
300
500
700
900
400
600
800
0
1000
1100
1200
1300
1400
1500
1600
1800
1700
Tiempo (s)
Figura 91. Comparativa DoD de la batería de ácido – plomo y la batería de ión – litio
6.1.- Introducción
Los algoritmos genéticos (AG) son métodos utilizados para resolver problemas de
optimización. Están basados en los procesos genéticos de los organismos biológicos. A lo
largo de las generaciones, sobreviven los que mejor se adaptan al sistema de acuerdo con
los principios de la selección natural y supervivencia. Los algoritmos genéticos son capaces
de resolver problemas reales.
En la naturaleza, los individuos de una población compiten con otros por los
recursos como la comida, el agua y el abrigo. También, miembros de una misma especie a
menudo compiten por atraer a la pareja. Estos individuos, los cuales son más afortunados
en la supervivencia y en la atracción de la pareja, tendrán más número de descendientes.
Esto significa que los genes de los mejores adaptados, o con capacidad, estarán en los genes
de los individuos de generaciones sucesivas. La combinación de buenas características de
diferentes antecesores puede producir, a veces, descendencia supercapacitada, en los cuales
el estado físico es mejor que el de los padres. De esta manera, las especies se adaptaran
mejor a su entorno.
El poder de los algoritmos genéticos viene del hecho de que la técnica es robusta, y
puede tratar con existo un rango ancho de problemas, incluyendo aquellos que son difíciles
de resolver por otros métodos. Los algoritmos genéticos no garantizan encontrar una
solución global al problema, pero normalmente son buenos para encontrar soluciones
aceptables para problemas rápidos.
Antes que los algoritmos genéticos puedan ser ejecutados, una adecuada
codificación (o representación) para el problema debe de ser ideada. También se requiere
una función de evaluación, la cual asigna el valor que le corresponde a la solución dentro de
la población. Durante la ejecución, los padres deben ser seleccionados para la reproducción
y la recombinación, generando la descendencia.
6.2.1.1.1.- Codificación
Es asumido que una solución al problema puede ser representada como un grupo
de parámetros. Estos parámetros (conocidos como genes) están unidos para formar una
cadena de valores, a menudo referida como cromosoma.
END
TERMINAR := True
END
END
6.2.1.1.3.- Reproducción
Padres 1 0 1 0 0 0 1 1 1 0 0 0 1 1 0 1 0 0 1 0
Descendencia 1 0 1 1 0 1 0 0 1 0 0 0 1 0 0 0 1 1 1 0
Punto de mutación
Descendencia 1 0 1 0 0 1 0 0 1 0
Descendencia mutada 1 0 1 1 0 1 0 0 1 0
6.2.1.1.4.- Convergencia
Hay una gran cantidad de técnicas de optimización, algunas de las cuales solo son
aplicables a dominios limitados. Esto es un método para resolver problemas de control
multipaso, los cuales solo son aplicables donde la función de evaluación global es la suma
de la función de evaluación para cada estado del problema, y no hay interacción entre
estados.
La búsqueda iterativa y la búsqueda del gradiente pueden estar combinadas para dar
una búsqueda iterativa ‘hillclimbing’. Una vez que un pico ha sido localizado, el ‘hillclimb’
empieza otra vez, pero con otro punto de comienzo elegido aleatoriamente. Esta técnica
tiene la ventaja de la simplicidad, y puede trabajar bien si la función no tiene demasiados
máximos locales.
Sin embargo, en cada ensayo aleatorio se lleva a cabo un aislamiento del dominio y
no del dibujo global. Los procesos de búsqueda aleatorios continúan para asignarles
ensayos incluso sobre el espacio de búsqueda. Esto significa que podrán ser evaluados
tantos puntos en la región encontrada para ser de una baja capacidad como en las regiones
encontradas para ser de una alta capacidad.
una ventaja de tener más puntos de cruzamiento es que el problema del espacio puede ser
buscado más concienzudamente.
En esta vista, un punto de corte puede ser parecido a dos puntos de corte con uno
de los puntos de corte unido al principio de la cadena.
Máscara de cruzamiento 1 0 0 1 0 1 1 1 0 0
Padre 1 1 0 1 0 0 0 1 1 1 0
Descendencia 1 1 0 0 0 0 1 1 1 1
Padre 2 0 1 0 1 0 1 0 0 1 1
Otros dicen que el cruzamiento con dos puntos se realizará pobremente cuando la
población ha convergido en gran parte, reduciendo el rendimiento del cruzamiento.
Durante la ejecución del proyecto se generan varias hojas Excel, las cuales guardan
la información generada durante dicha ejecución. A continuación, se puede observar la
relación que existe entre el usuario y el programa.
ALGORITMOS GENÉTICOS
La relación entre las hojas del libro de Excel depende de si son entrada de datos o
salidas de programa se muestran en la siguiente figura, Figura 99:
ENTRADAS
•Datos_EV
•Datos_ EC
•Datos
ALGORITMOS GENÉTICOS
SALIDAS
•Costes
•Soluciones
•Solucion_optima
El programa interacciona con el usuario por medio de ventanas en las cuales les
permite realizar una serie de acciones. En primer lugar, el programa da la opción de leer los
datos de la ubicación de las estaciones de carga originales a partir de las cuales se obtendrá
el óptimo o bien que el mismo programa genere unas soluciones iniciales aleatorias.
Figura 100. Ventana gráfica del programa de algoritmos genéticos donde da a elegir la manera de obtener la
solución inicial del programa de optimización
Figura 101. Ventana gráfica del programa de algoritmos genéticos donde hay que especificar la longitud y
número de cromosomas
Figura 102. Ventana gráfica del programa de algoritmos genéticos donde hay que especificar el número de
iteraciones
Una vez introducidos estos datos, el programa realiza los cálculos hasta que obtiene
la solución óptima.
Las soluciones que se obtiene en este programa son a través de hojas ubicadas en el
fichero de Excel. Estas son:
Costes: muestra la variación de los costes. Esta hoja muestra dos columnas:
o Número de iteración.
La convergencia a este mínimo puede ser más rápida o más lenta dependiendo de
cómo se crucen o muten las soluciones y del número de cruces que se realicen en cada
iteración.
Los cromosomas, es decir, las soluciones posibles que se pueden llevar a cabo, se
cruzan de una forma aleatoria. Esto es, las soluciones que se cruzan y el punto de cruce son
aleatorios. Por todo ello, habrá ocasiones en las cuales se obtendrán mejores soluciones con
menos iteraciones que en otras ejecuciones del programa.
En las mutaciones ocurre lo mismo pero todavía más aleatorio porque se realiza un
bucle que se repite tantas veces como estaciones de carga existan en las soluciones. En el
ejemplo que se ha desarrollado, el número de veces que se repite el bucle es de catorce.
Dentro del bucle hay una condición, si el número aleatorio que se ha originado al comienzo
del bucle es mayor de cinco, entonces se realiza la mutación, sino se repite el bucle.
ejecución del programa depende del número de cruces que se realicen y por supuesto del
número de iteraciones.
Costes
180
160
140
120
100
80
60
40
20
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200
Figura 103. Representación de los costes de las pérdidas frente al número de iteraciones (dos cruces)
Costes
180
160
140
120
100
80
60
40
20
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200
Figura 104. Representación de los costes de las pérdidas frente al número de iteraciones (diez cruces)
El coste de pérdidas obtenido es menor que los iniciales. El resultado obtenido es:
Figura 105. Mapa temático donde se representa la zona que abarca cada estación de carga y la localización de
las mismas
Masa del vehículo: el peso del vehículo en vacío más el conductor y/o la
carga del maletero y viajeros
Una vez obtenidos todos los valores de entrada incluidos en el fichero de Matlab
habilitado para ello, se ha desarrollado un método para obtener de esta aplicación concreta
lo siguiente: los valores de las prestaciones, el consumo, la autonomía y la vida de las
baterías del vehículo eléctrico seleccionado y el análisis económico de la inversión.
Otro factor que puede tener una influencia significativa en el alcance del
vehículo es la cantidad de accesorios auxiliares que están consumiendo
energía. Entre el caso más favorable, que consiste en llevar el mínimo de
accesorios en funcionamiento (20 W) y el más desfavorable que consiste en
llevar el máximo de accesorios funcionando (unos 760 W entre
climatizador, iluminación de la carretera, señalización del vehículo, radio,
etc.) esto puede suponer una disminución del 9,05% en la autonomía en el
caso de la batería ácido – plomo, 8,75% para la batería de níquel – cadmio
y 8,54% para el caso de la tecnología de ión – litio.
Se pretende que este programa sea una herramienta informática con la cual se
pueda mejorar el sistema de distribución de las estaciones de carga para los vehículos
eléctricos. La labor investigadora realizada es importante para el desarrollo tecnológico y
para el diseño de los sistemas de distribución, y sobre todo, para la planificación de las
estaciones de carga de una zona.
8.1.- Introducción
En este apartado se van desarrollar y analizar las diferentes partes en las que consta
los programas de simulación de los diferentes tipos de vehículos eléctricos, es decir, el
vehículo puramente eléctrico y el híbrido en paralelo y la localización de las estaciones de
carga de los vehículos eléctricos.
En esta tesis se han tenido en cuenta seis tipos diferentes de ciclos de conducción
los cuales se han agrupado en dos grupos:
Reglado
o FTP75
Experimentales
o Tráfico fluido
8.2.1.1.- FTP75
El ciclo tiene dos fases: fase de ‘arranque en frío’ de 505 segundos en una distancia
proyectada de 5,78 kilómetros, con una velocidad media 41,2 km/h, y ‘fase transitoria’ 864
segundos. Las dos fases se separan por detener el motor durante 10 minutos. En los
EE.UU. se aplican unos factores de ponderación de 0,43 y 0,57 para la primera y segunda
fase, respectivamente.
El ciclo simula un recorrido urbano de 12,07 kilómetros (7,5 millas) con paradas
frecuentes. La velocidad máxima es de 91,2 km/h (56,7 millas/h) y la velocidad media es
del 31,5 km/h (19,6 millas/h). Un promedio de peso de la ciudad (55 %) y de la carretera
(45 %) de las economías de combustible se utiliza para determinar el impuesto por alto
consumo.
Ciclo US FTP75
100
90
80
70
Velocidad (Km/h)
60
50
40
30
20
10
0
0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1050 1100 1150 1200 1250 1300 1350 1400 1450 1500 1550 1600 1650 1700 1750 1800 1850 1900
Tiempo (segundos)
Los parámetros básicos del ciclo urbano elemental, parte 1 del Test 1, (UNECE
Elementary Urban Cycle) son los siguientes:
50
40
Velocidad (Km/h)
30
20
10
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200
Tiempo (segundos)
La prueba continúa sin detener el coche en una segunda parte que busca reproducir
las condiciones de circulación en carretera. En este caso, tras un periodo al ralentí de 20
segundos, se acelera hasta 70 km/h en cuarta, velocidad que se mantiene durante 50
segundos. Después, se baja a 50 km/h y se mantiene también durante un tiempo antes de
subir de nuevo a 70 km/h, y mantener esa velocidad otros 50 segundos. Luego se acelera
hasta 100 km/h, se mantiene esa velocidad unos segundos y se sigue acelerando hasta
alcanzar los 120 km/h. Tras 10 segundos a esa velocidad, se detiene el coche. La distancia
teórica recorrida es de casi 7 kilómetros en 6 minutos y 40 segundos a una media de
62,6 km/h.
Los parámetros básicos del ciclo extra-urbano, parte 2 del Test 1 (UNECE Extra-
Urban Driving Cycle) son los siguientes:
120
100
Velocidad (Km/h)
80
60
40
20
0
0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 380 400
Tiempo (segundos)
En el caso de vehículos con poca potencia, la prueba continúa sin detener el coche
en una segunda parte que busca reproducir las condiciones de circulación en carretera. En
este caso, tras un periodo al ralentí de 20 segundos, se acelera hasta 70 km/h en cuarta,
velocidad que se mantiene durante 50 segundos. Después, se baja a 50 km/h y se mantiene
también durante un tiempo antes de subir de nuevo a 70 km/h y mantener esa velocidad
otros 50 segundos. Luego se acelera hasta 90 km/h, se mantiene esa velocidad 84 segundos
y se detiene el coche.
Tabla 27: Parámetros básicos del ciclo de conducción extra-urbano (Vehículos con baja potencia)
90
80
70
Velocidad (Km/h)
60
50
40
30
20
10
0
0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 380 400
Tiempo (segundos)
Figura 109. Ciclo de conducción UNECE Extra-Urban Driving Cycle (Low-Powered Vehicles)
Borja Pintos [57] desarrolló una metodología que permitió generar ciclos de
conducción para bancos de rodillos, los cuales representan la conducción en tráfico real
urbano. Para ello, registró los datos mediante un equipo de captación de datos equipado a
Los parámetros básicos del ciclo con tráfico fluido son los siguientes:
50
40
Velocidad (Km/h)
30
20
10
0
0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000
Tiempo (segundos)
Los parámetros básicos del ciclo con tráfico no fluido son los siguientes:
50
40
Velocidad (Km/h)
30
20
10
0
0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1050 1100 1150 1200 1250 1300 1350 1400 1450 1500
Tiempo (segundos)
A continuación se muestran las variables utilizadas para las simulaciones tanto del
vehículo eléctrico como del vehículo eléctrico híbrido y en cada apartado se mostrarán los
datos de cada simulación.
%MOTOR ELECTRICO
%MOTOR TÉRMICO
%BATERÍA
rendimiento_max_teorico=ren_b*rend_e_m*rend_t_m;
%COMBUSTIBLE
energia_especifica_gasoil=d_gasoil*pc_gasoil;
energia_especifica_gasolina=d_gasolina*pc_gasolina;
Datos_EV
Datos_EC
Datos
8.3.2.1.- Datos_EV
Estos datos lo que representan son la demanda de los vehículos eléctricos. Para la
obtención de estos datos se ha obtenido el área de cada manzana y con la densidad de
población se ha obtenido la población que le corresponde a cada manzana. Se ha supuesto
que las familias están formadas por cuatro miembros y que poseen dos vehículos. De estos
vehículos se ha supuesto que un 20 % son eléctricos.
8.3.2.2.- Datos_EC
Estos datos son los que ha introducido el usuario como solución inicial a partir de
la cual se genera la descendencia con mejores cualidades, en este caso en concreto, se
obtienen soluciones con menor coste.
8.3.2.3.- Datos
Las variables que se utilizan durante el programa de algoritmos genéticos son las
siguientes:
Fracción EV ............................................................................. 20 %
En primer lugar, se muestra el código en Matlab a partir del cual se obtienen los
valores de las variables, como masa del vehículo, área frontal, coeficiente de rozamiento,
velocidad dependiendo del ciclo elegido, tipo de batería, etc.
if strcmp(Velocidad,'Reglado')
end
if strcmp(Velocidad,'Más')
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 8 Anexos 211
end
if strcmp(Velocidad,'Experimental')
end
if strcmp(Velocidad,'FTP75')
v=xlsread('FTP75.xlsx','D:D');
end
end
end
end
if strcmp(Velocidad,'Tráfico fluido')
v=xlsread('Experimental.xlsx','C:C');
end
if strcmp(Velocidad,'Tráfico no fluido')
v=xlsread('Experimental.xlsx','B:B');
end
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 8 Anexos 212
t=ones(N,1);
for i=2:N
t(i,1)=t(i-1,1)+1;
end
V=[t,V1];
Tipo_bateria=questdlg('Elegir tipo de
bateria:','Bateria','Acido-Plomo','NiCd', 'Ion-Litio',
'Acido-Plomo');
if strcmp(Tipo_bateria,'Acido-Plomo')
Tipo_bat=1;
end
if strcmp(Tipo_bateria,'NiCd')
Tipo_bat=2;
end
if strcmp(Tipo_bateria,'Ion-Litio')
Tipo_bat=3;
end
find_system('Name','VE_Completo');
CY=1;
Distancia_acumulativa=0;
DoD_acumulativo=0;
CR_acumulativo=0;
while DD<0.9
open_system('VE_Completo');
set_param(gcs,'simulationCommand','Start');
DoD_end(CY)=DoD(end,2);
CR_end(CY)=CR(end,2);
D_end(CY)=D(end,2)/1000;
Distancia_acumulativa=D(end,2); %D(1);
DoD_acumulativo=DoD(end,2); %DoD(1);
CR_acumulativo=CR(end,2); %CR(1);
set_param(gcs,'simulationCommand','Stop');
CY=CY+1;
end
close_system('VE_Completo');
plot(D_end,DoD_end,'k+');
ylabel('Profundidad de descarga');
xlabel('Distancia recorrida**/Km');
Figura 114. Vista general del modelo del vehículo eléctrico en Simulink
Este bloque es el encargado de calcular la fuerza motriz necesaria para cumplir los
requerimientos de la marcha del vehículo.
En primer lugar, hay que tener en cuenta las pérdidas mecánicas de la reductora y
del diferencial así como otros elementos como las juntas cardan. Un coeficiente de
recuperación es aplicado en el caso de frenada regenerativa que representa el nivel de
recuperación de energía por el motor y la batería. Los interruptores (bloque switch) sirven
para seleccionar el camino que atraviesa la energía en función de si la energía es positiva
(tracción) o negativa (frenada) y así aplicar correctamente los factores de pérdida.
La última función desarrollada en este bloque calcula el factor de intensidad que hay
que añadir a la carga eliminada de la batería para conseguir la carga eliminada en el siguiente
paso.
Los valores de las variables se obtienen de un fichero, cuyos datos de entrada están
descritos en el apartado 8.3.- Datos de entrada. Este fichero se ejecuta en Matlab, el cual te
permite elegir el tipo de ciclo de conducción y el tipo de batería como se ha explicado en el
modelo de vehículo eléctrico, a continuación, carga las variables en el espacio de trabajo de
Matlab y a partir de aquí, Simulink toma los valores.
8.5.1.1.- Desmultiplicador
Con la ayuda de este bloque, se obtiene el valor del par que el motor eléctrico
puede entregar en función de la velocidad de giro. La curva de par máximo de un motor
eléctrico tiene dos zonas claramente diferenciadas: una, en la que el par es constante, y la
otra, en la que trabaja a potencia constante. La velocidad a partir de la cual se pasa de una
zona a otra es la velocidad base, la cual en esta tesis se ha supuesto de 1500 rpm.
Figura 121. Regulación de velocidad a potencia constante y a par constante de un motor eléctrico
Las entradas a este bloque son el par y el régimen de giro obteniendo el valor del
rendimiento del motor térmico. Para ello, a partir de la gráfica de eficiencia del motor
térmico elegido, se divide en cuatro sectores. Estos sectores están representados por los
cuatro interruptores, los cuales dejan pasar el valor del rendimiento en tanto por cien
obtenido de los valores de las tablas de cada sector.
entonces se activa la salida del bloque y permanece activa durante el intervalo de tiempo
que la señal de velocidad se halle dentro del intervalo indicado anteriormente.
Figura 128. Distribución de las fuerzas de frenado en las ruedas delanteras y las traseras
Las líneas oblicuas indican la máxima deceleración (j) que se puede alcanzar
dependiendo del coeficiente de adherencia entre el neumático y el asfalto. Dicha
Figura 129. Fuerzas que actúan en un vehículo durante el frenado en una superficie plana
La mínima fuerza de frenado en las ruedas traseras, y por lo tanto, máxima fuerza
de frenado en las ruedas delanteras, cuando las ruedas delanteras están trabadas, estipulado
por la regulación ECE, está expresado como
Despejando ,
En esta ecuación se puede observar que el límite de max para cumplir con la
regulación ECE está en función de la resistencia en el frenado, q. Por lo tanto, el valor de
este índice, el cual permite conocer el valor de max, puede ser obtenido a través de la
siguiente fórmula:
Para que las señales de salida se activen es necesario que se cumplan dos señales a la
vez; y una de éstas tiene tres opciones, es decir, en el caso que la señal subir marcha se
active tiene que ser: la marcha que circule el vehículo menor al número de marchas que
posea el vehículo y a la vez que una de estas tres esté activa: el régimen de giro del motor
sea mayor al valor máximo del régimen de giro estipulado; que la salida del bloque
“Comparador de consumos con diferentes marchas” posea activa la salida “subir marcha”, y por
último, que el bloque “comparador consumos en modo recarga” tenga activa la salida en el caso
que sea mejor subir una marcha.
Para el caso de reducir marcha, las señales serán: la marcha que circule el vehículo
mayor a uno y, además, que se cumpla una de estas tres condiciones: que el régimen del
giro del motor sea menor al valor mínimo del régimen de giro estipulado; que la salida del
bloque “Comparador de consumos con diferentes marchas” posea activada la salida “reducir marcha”
y por último, que el bloque “comparador consumos en modo recarga” tenga activa la salida en el
caso que sea mejor reducir una marcha.
Las entradas a este bloque son: la marcha a la que circula el vehículo; el par del
motor térmico; y el régimen de giro. Obteniendo como salidas, dos señales booleanas:
“subir marcha” y “reducir marcha”.
El rendimiento del motor se calcula a partir del par entregado y del régimen de giro.
Los valores que se introducen a este subsistema no están situados en los mismos puntos de
la cadena cinemática, por lo que hay reubicarlos. La señal de entrada “rpma” está medida en
la salida del acoplador, por lo que hay multiplicar su valor por el índice de transmisión para
calcular el régimen de giro a la salida del motor.
Sin embargo, la señal de entrada “Nmt” se refiere al par motor entregado por el
mismo, por lo que hay que multiplicar por el inverso del índice de transmisión para
reubicarlos en la salida del acoplador. Recíprocamente, al multiplicar el valor de par
entregado por el motor por el índice de transmisión obtenemos el par aplicado en el
acoplador, y así obtenemos el valor “Nma”. Si éste, lo multiplicamos por la inversa del
índice de transmisión de las marchas superior e inferior nos proporciona el par exigido al
motor con esas marchas. A partir de las cuales se calcula el rendimiento térmico con una
marcha más o menos de la actual dependiendo en que caso se esté.
En primer lugar se calcula el par que podría dar el motor eléctrico. A continuación,
se calcula el par que entrega a determinado régimen de giro siguiendo la línea de máxima
eficiencia, bloque “par rendimiento cte”. A partir de estos dos valores, se selecciona el mínimo
para calcular el rendimiento del motor térmico.
Donde P es el par entregado y n es el régimen de giro en rpm del eje de salida del
motor.
Este bloque tiene como entrada la señal de velocidad y como salida la señal
numérica “cond”, la cual indica el modo de conducción. Si ésta posee un valor 1, indica que
se está realizando una conducción urbana y si posee un valor 2, indica que está en el caso
de conducción interurbana.
Las señales de entrada a este bloque son: la potencia de recarga efectiva de las
baterías, el par motor requerido a la salida del acoplador de par, el régimen de giro del eje
de salida del acoplador, la marcha engranada, la señal de frenado y el tipo de conducción.
Los bloques “modo on” y “modo off” no están nunca a la vez activados ya que cuando
no está activado envía al otro bloque una señal que lo desactiva. El bloque “modo on” se
desactiva cuando la señal “carga PPS” es mayor al valor máximo de la carga aceptada en la
batería, mientras que el “modo off” se desactiva cuando llega al valor mínimo. Al iniciar la
simulación y para que no exista un conflicto de activación entre ambos bloques, se ha
incluido una entrada escalón, la cual durante un segundo permanece inactiva y así no
permite que el bloque “modo off” no se active.
Este bloque está compuesto por tres subsistemas: Recarga PPS, Propulsión híbrida
y eléctrico a baja velocidad. Cada bloque tiene sus entradas específicas pero los tres poseen
las mismas señales de salida: el par motor del motor térmico, la potencia de recarga de
salida, el par del motor eléctrico, el modo de funcionamiento y el régimen de giro del
motor térmico. Las salidas de cada bloque se unen en un bus de datos, los cuales se suman
para obtener una única señal. Esta señal pasará a la salida del bloque “modo on”, “out”, si el
interruptor está activado. La señal booleana denominada “señal” es la que lleva la
información de la activación y desactivación de este bloque. Si “señal” posee el valor de
“false”, entonces la salida del bloque “modo on” denominada “off” se activa.
Las entradas de datos que posee este bloque son: el nivel de carga de las baterías, la
velocidad del vehículo, la velocidad mínima para encender el motor térmico, el par motor
requerido en el acoplador, el régimen de giro a la salida del acoplador y la marcha a la que
circula el vehículo.
Para el cálculo del par motor térmico a máxima eficiencia se multiplica el régimen
de giro del acoplador por el índice de transmisión según la marcha engranada. Esto se
realiza para situar la velocidad de giro en el motor térmico, valor de la señal “rmpt”. Una
vez obtenida esta velocidad y con la función polinómica vista en el apartado 8.5.5.2.1.-
Bloque Rend_marcha+1 se calcula el par térmico con la máxima eficiencia. Este valor se
multiplica por el índice de desmultiplicación para obtener su valor en el acoplador de par.
A este par obtenido se le resta el par motor del acoplador resultando el par motor
sobrante que se puede utilizar para recargar la batería. Este par está limitado por la
capacidad del motor térmico, por lo tanto, hay que calcular el par máximo que puede
soportar el motor eléctrico. Para ello, se multiplica la velocidad de giro del acoplador por el
índice de transmisión del motor eléctrico obteniendo la velocidad de giro de éste. A partir
de este valor y con la ayuda del bloque “motor eléctrico” se calcula el par de éste, el cual es
multiplicado por el índice de transmisión; consiguiendo el par máximo que el motor
eléctrico puede entregar al acoplador o el acoplador al motor eléctrico para trabajar como
generador. Para avalar que el par sobrante del motor térmico no sea en ningún caso
superior al par que puede soportar el motor térmico se ha colocado el bloque “Min”.
El par motor recibido por el motor eléctrico, el cual en este caso trabaja como
generador, se calcula a partir del producto de la señal que se obtiene a la salida del bloque
“Min” y la inversa del índice de éste. Multiplicando este valor por el régimen de giro del
acoplador convertido a radianes por segundo y por el índice de transmisión del motor
eléctrico se obtiene la potencia de recarga.
La señal de modo obtendrá el valor numérico de tres, el cual corresponde con este
modo de funcionamiento.
Este bloque lo que realiza, como su propio nombre indica, es la activación del
modo híbrido donde trabajan conjuntamente el motor térmico y el motor eléctrico.
Este bloque posee cinco entradas, las cuales son: la velocidad del vehículo, la
velocidad mínima para encender el motor térmico, el par motor requerido en el acoplador,
el régimen de giro a la salida del acoplador y la marcha a la que circula el vehículo.
La diferencia entre el par motor solicitado por la dinámica del sistema y el par
motor entregado por el motor térmico siguiendo la curva de máxima eficiencia;
multiplicada, a su vez, por el inverso del índice de transmisión del motor eléctrico nos
proporciona el par suministrado por este último para satisfacer los requerimientos de
potencia, correspondiente a la señal de salida “Nme”. Siendo ésta el par motor entregado
por el motor eléctrico.
La potencia de recarga será la opuesta del producto del par motor entregado por el
motor eléctrico y el régimen de giro de éste, convirtiéndolo previamente a rad/s, debido a
que en esta situación se descarga la batería.
La señal de modo obtendrá el valor numérico de dos, el cual corresponde con este
modo de funcionamiento.
Este bloque tiene cuatro entradas: la velocidad del vehículo, la velocidad mínima, el
par motor y el régimen de giro solicitados en el acoplador. Se activa cuando la velocidad del
vehículo es menor a la velocidad mínima necesaria para que se active el motor térmico.
El par motor eléctrico se obtiene multiplicando el par de salida del acoplador por el
inverso del índice de transmisión del motor eléctrico; mientras que la potencia de recarga se
calcula como el opuesto del producto de dicho par por el régimen de giro del motor
eléctrico convertido a rad/s.
En este bloque el par motor y el régimen de giro del motor térmico son nulos
debido a que está apagado el motor térmico.
La señal de modo obtendrá el valor numérico de uno, el cual corresponde con este
modo de funcionamiento.
El par motor del motor eléctrico se consigue por medio del producto del régimen
de giro del acoplador por el inverso de la relación de transmisión del motor eléctrico. El
valor de la potencia de recarga se calcula como el opuesto al producto del régimen de giro
del motor eléctrico convertido a rad/s y el par entregado por el mismo.
Este bloque se activa cuando el motor térmico trabaja con mejor eficiencia que si
trabaja en modo híbrido.
Las entradas de datos que posee este bloque son: la velocidad del vehículo, el par
motor requerido en el acoplador, el régimen de giro a la salida del acoplador, la marcha a la
que circula el vehículo y el modo de conducción.
3. El régimen de giro del motor térmico con la marcha más alta que posee el
vehículo no sea inferior al régimen de giro mínimo admitido.
Las salidas de este bloque son: la potencia de recarga y el par motor del motor
eléctrico los cuales son nulos ya que el motor eléctrico está apagado, el modo de circulación
que en este caso posee el valor numérico de seis, la señal booleana “out” la cual indica si
está el bloque activado o no, el par motor del motor térmico y el régimen de giro de éste.
Las entradas de datos que posee este bloque son: el nivel de carga de las baterías, la
velocidad del vehículo, el par motor requerido en el acoplador y el régimen de giro a la
salida del acoplador.
En este bloque hay un interruptor el cual deja pasar el valor del par motor en el
motor eléctrico si el valor de la carga de la batería es menor que la carga máximo admitida y
la velocidad del vehículo es mayor a 15 Km/h ya que debajo de esta velocidad el
El par motor en el motor eléctrico se obtiene del valor mínimo entre el par aplicado
al motor eléctrico, ya que éste trabaja como generador, y el par motor obtenido a partir del
régimen de giro del acoplador.
El par motor en motor térmico y el régimen de giro de éste son nulos, ya que éste
no trabaja en estas condiciones.
El primer cálculo que realiza este bloque es el del rendimiento eléctrico; el cual es
obtenido a partir del par motor y el régimen de giro del acoplador multiplicado por el
índice de transmisión del motor eléctrico para trasladarlo a la posición del motor eléctrico.
Este rendimiento es multiplicado por el rendimiento de la batería.
Como en este bloque el flujo de energía va en los dos sentidos, para calcular la
potencia de recarga de salida hay que tener en cuenta el sentido de la energía, para ello, se
ha introducido un interruptor el cual diferencia si la potencia de recarga de salida es
positiva o negativa. En el caso de que sea positiva, el motor eléctrico está convirtiendo la
energía mecánica en energía eléctrica por lo que la potencia efectiva de recarga, en este
caso, es el producto de la potencia de recarga de salida por el rendimiento global
anteriormente calculado. Por el contrario, si la potencia de recarga de salida es negativa, el
motor eléctrico está convirtiendo la energía eléctrica en energía mecánica, por lo tanto, la
potencia efectiva de recarga en este caso es el cociente de la potencia de recarga de salida
entre el rendimiento de la conversión.
En primer lugar, hay que calcular la potencia bruta entregada por el motor de
combustión interna el cual se calcula multiplicando el par motor por el régimen de giro del
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 8 Anexos 248
motor térmico en radianes por segundo y el resultado se divide por el rendimiento de dicho
motor, que se habrá calculado a partir del bloque rendimiento térmico, cuyas entradas son
el par motor y el régimen de giro de dicho motor.
Una vez obtenido la potencia bruta, se divide por la energía específica del
combustible expresada en J/l para obtener el consumo instantáneo de combustible. En la
siguiente tabla se muestran las características de los combustibles:
Con este bloque, se cuentan los ciclos de carga y descarga que ha realizado la
batería. El ciclo de carga comprende desde que la carga de la batería está por debajo del
52 % de la carga máxima más el margen para la energía de aceleración dependiendo del
modo en que circule, o bien urbano o bien interurbano; y esté sobre el 68 % de la carga
máxima.
Esta señal de entrada es enviada a cuatro subsistemas, los cuales analizan los
rendimientos: del motor térmico y eléctrico y en los casos de frenada regenerativa;
únicamente motor eléctrico; solo motor térmico; y por último, la distribución de modos.
Este análisis se activa si la señal de modos posee el valor numérico de cinco, el cual
corresponde a la situación de frenada regenerativa.
El rendimiento mecánico se obtiene del cociente entre la integral del producto del
par y el régimen de giro del motor eléctrico (la energía bruta mecánica que llega al motor
eléctrico) y la integral del producto del par y el régimen de giro del acoplador siendo la
energía de frenada en las ruedas delanteras.
Este bloque analiza el rendimiento del motor térmico en los diferentes modos en
los cuales puede estar funcionando.
Este bloque analiza el rendimiento del motor eléctrico en los diferentes modos en
los cuales puede estar funcionando.
tiempo de marcha en cada modo de funcionamiento. Dividiendo este valor por el tiempo
total que el circula el vehículo y multiplicando el resultado por cien, se obtiene el porcentaje
de tiempo que cada modo de funcionamiento permanece activo.
En primer lugar, se calcula el gasto del combustible que realiza el vehículo como si
fuera un vehículo convencional. Para ello, se multiplica el consumo realizado durante
100 km, la distancia que puede recorrer la batería y el gasto del combustible.
En segundo lugar, se calcula el gasto de la batería. Con el fin, de obtener este valor,
se divide el numero de ciclos de la batería entre la distancia recorrida (se integra la señal de
velocidad). Este valor, se multiplica por el número de ciclos que posee de vida dicha batería
se le suma uno para tener en cuenta el valor de la batería inicial; y el resultado se multiplica
por el precio de la batería, obteniendo así el gasto en la misma.
El índice de carga del motor térmico se obtiene dividiendo el par motor del motor
térmico y el par motor máximo del mismo. Esta señal está activada si el modo de
funcionamiento corresponde a térmico puro, valor numérico de la señal modo igual a seis.
El índice de carga del motor eléctrico, se obtiene dividiendo el par motor del motor
eléctrico y el par motor máximo del mismo. Esta señal está activada si el modo de
funcionamiento corresponde a eléctrico a baja velocidad o eléctrico puro, valor numérico
de la señal modo igual a uno o cuatro, respectivamente.
8.6.1.- Módulo1
8.6.1.1.- Variables
8.6.1.2.- Geneticos
Sub geneticos()
Dim iter As Integer
Dim num_sol As Integer
Worksheets("Costes").Activate
Range("A2:C2000").Select
Selection.ClearContents
Worksheets("Soluciones").Activate
Range("A2:C2000").Select
Selection.ClearContents
Worksheets("Solucion_optima").Activate
Range("b2:d2000").Select
Selection.ClearContents
Load Datos_entrada
Datos_entrada.Show
If opcion = 0 Then
Call leer_datos_fichero
Else
Load Datos
Datos.Show
Call genetic_aleatorio
End If
Worksheets("Costes").Activate
ActiveSheet.Cells(3, 5) = num_cromosomas
ActiveSheet.Cells(2, 5) = size_cromosomas
Call leer_datos_EV
Call leer_datos_iniciales
Call sort_by_fitness1
Dim num_fila As Integer
num_fila = 2
For iter = 0 To num_inter
Call ordenar
coste_total(iter) = es(0).fitness_values
Worksheets("Costes").Activate
For i = 0 To iter
Dim c As Integer
c=i+2
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 8 Anexos 261
ActiveSheet.Cells(c, 1).Value = i
ActiveSheet.Cells(c, 2).Value = coste_total(i)
Next i
Call do_crossovers
Call do_mutations
Call do_mutations_potencia
Call ordenar
num_sol = 0
Call get_chromosome(num_sol)
For i = 0 To size_cromosomas - 1
Worksheets("Soluciones").Activate
ActiveSheet.Cells(num_fila + i, 1).Value = cromosomas(num_sol, i, 0)
ActiveSheet.Cells(num_fila + i, 2).Value = cromosomas(num_sol, i, 1)
ActiveSheet.Cells(num_fila + i, 3).Value = cromosomas(num_sol, i, 2)
Next i
For i = 0 To size_cromosomas - 1
Worksheets("Solucion_optima").Activate
ActiveSheet.Cells(2 + i, 2).Value = cromosomas(num_sol, i, 0)
ActiveSheet.Cells(2 + i, 3).Value = cromosomas(num_sol, i, 1)
ActiveSheet.Cells(2 + i, 4).Value = cromosomas(num_sol, i, 2)
Next i
num_fila = num_fila + i
Next iter
Call asignar_EC
End Sub
8.6.1.3.- Leer_datos_fichero
Esta función, como antes se ha comentado, lo que hace es leer a través de la hoja
del fichero de Excel denominada “Datos_EC” las soluciones de las estaciones de carga. Las
coordenadas X e Y se guardan en es(i).poblacion(j,0) y es(i).poblacion(j,1), respectivamente.
La potencia de cada estación de carga se guarda en es(i).poblacion(j,2).
8.6.1.4.- Genetic_aleatorio
Esta función realiza de forma aleatoria las soluciones iniciales de las estaciones de
carga. A priori, se ha llamado a un formulario el cual le pide al usuario el número de
soluciones iniciales y la cantidad de estaciones de carga por cada solución, Figura 101.
Estos valores son cargados a las variables correspondientes. El número de soluciones
iniciales se carga a la variable “num_cromosomas” y en el caso del número de estaciones de
carga por solución a “size_cromosomas”. Una vez obtenidos estos datos se procede a la
generación de las soluciones de forma aleatoria con la función llamada Aleatoria. Para que
la función Aleatoria funcione, hay que introducir un valor por referencia a la variable para
que ésta te devuelva la solución. En este ejemplo, se han introducido el valor de 9300 en el
caso de es(i).poblacion(j,0) (coordenada X) porque la longitud de la geobase es de 9300
metros y su origen es de 679500 aproximadamente, por ello, los valores obtenidos no se
salen de la geobase, y por lo tanto, son válidos. En el segundo caso, es(i).poblacion(j,1)
(coordenada Y) ocurre lo mismo; el origen está en 4608600 y el tamaño es de 9300 metros.
Coinciden ambos datos porque la geobase es un cuadrado. Estos datos son elegidos para
que se sitúen las estaciones de carga dentro de la zona de estudio. Finalmente, se elige,
también, aleatoriamente la potencia de las estaciones de carga, en este ejemplo se han
elegido las potencias de 3700 W y 7400 W.
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 8 Anexos 263
El programa es el siguiente.
8.6.1.5.- Leer_datos_EV
8.6.1.6.- Leer_datos_iniciales
8.6.1.7.- Sort_by_fitness1
Este procedimiento evalúa las soluciones y las ordena según el valor que le devuelve
la función fitness_func. Este valor es almacenado en la variable es(i).fitness_values, la cual
almacena el coste de las pérdidas.
8.6.1.8.- Ordenar
Esta función ordena los costes de las pérdidas de menor a mayor. El programa
desarrollado es el siguiente:
8.6.1.9.- Do_crossovers
Si la variable “igual” posee el valor ‘False’ se aumenta en uno el índice. Este índice
indica el sitio donde se almacena esta solución. Si la variable “igual” posee el valor ‘True’,
entonce, no avanza y en la siguiente solución se sobrescribe perdiendo dicha solución.
8.6.1.10.- Do_mutations
Como su propio nombre indica, esta función realiza las mutaciones en la población
para que así converja antes a la solución óptima.
Dentro del bucle de las mutaciones se van a elegir dos numeros aleatorios, los
cuales están comprendidos entre cero y el número de casos posibles de localización de las
estaciones de carga introducidas al principio del programa. La primera variable, “mut1”,
indica la solución que se va a mutar y la segunda, “mut2”, el elemento a mutar.
Si la variable “igual” posee el valor ‘False’ se aumenta en uno el índice que indica el
sitio donde se almacena esta solución y si posee el valor ‘True’ no avanza y en la siguiente
solución se sobrescribe perdiendo dicha solución.
igual = True
End If
Next asd
If igual = False Then
es(num_cromosomas + m).fitness_values = cos
m=m+1
End If
End If
Next j
End Sub
8.6.1.11.- Do_mutations_potencia
8.6.1.12.- Get_cromosome
8.6.1.13.- Fitness_func
donde:
coste_perdida_pot = 0
coste_perdida_ind = 0
For i = 0 To num_EV - 1
If poblacion_EV(i, 2) = n Then
coste_perdida_pot = coste_perdida_pot + poblacion_EV(i, 0) * poblacion_EV(i, 1) * p
/g
coste_perdida_ind = coste_perdida_ind + poblacion_EV(i, 0) * poblacion_EV(i, 1) * k
/v
End If
If poblacion_EV(i, 0) > dmax Then
coste_perdida_pot = coste_perdida_pot * 2
coste_perdida_ind = coste_perdida_ind * 2
End If
Next i
Else
coste_perdida_pot = 0
For i = 0 To num_EV - 1
If poblacion_EV(i, 2) = n Then
coste_perdida_pot = coste_perdida_pot + poblacion_EV(i, 0) * poblacion_EV(i, 1) * p
/g
coste_perdida_ind = coste_perdida_ind + poblacion_EV(i, 0) * poblacion_EV(i, 1) * k
/v
End If
If poblacion_EV(i, 0) > dmax Then
coste_perdida_pot = coste_perdida_pot * 2
coste_perdida_ind = coste_perdida_ind * 2
End If
Next i
coste_perdida_pot = coste_perdida_pot * 2
coste_perdida_ind = coste_perdida_ind * 2
End If
coste_global = coste_global + coste_perdida_pot + coste_perdida_ind
Next n
fitness_func = coste_global
End Function
8.6.1.14.- Aleatorio
8.6.1.15.- Asignar_EC
8.6.2.- Módulo2
8.6.2.1.- Asignar_EC
Sub asignar_EC()
Call leer_solucion_optima
Call designacion
End Sub
8.6.2.2.- Solucion_optima
Sub leer_solucion_optima()
Worksheets("Solucion_optima").Activate
For i = 0 To num_cromosomas
poblacion_EV_op(i, 0) = Cells(i + 2, 2)
poblacion_EV_op(i, 1) = Cells(i + 2, 3)
poblacion_EV_op(i, 2) = Cells(i + 2, 4)
Next i
End Sub
8.6.2.3.- Designacion
Sub designacion()
Dim ch_size As Integer
ch_size = size_cromosomas
For bb = 0 To num_EVs - 1
poblacion_EV_ES(bb, 0) = 0
poblacion_EV_ES(bb, 1) = 0
poblacion_EV_ES(bb, 2) = 0
Next bb
For a = 0 To num_EVs - 1
dist_op(a, 0) = 0
dist_op(a, 1) = 0
Next a
For l = 0 To num_EV - 1
For j = 0 To ch_size - 1
dist_op(j, 0) = Sqr((poblacion_EV_op(j, 0) - demanda_EV(l, 0)) ^ 2 + (poblacion_EV_op(j,
[2] CHAN, C.C.; BOUSCAYROL, A.; CHEN, K., “Electric, Hybrid, and Fuel-Cell
Vehicles: Architectures and Modeling”, Vehicular Technology, IEEE Transactions on, Febrero
2010, vol.59, no.2, p.589-598
http://www.minetur.gob.es/energia/balances/Balances/LibrosEnergia/Energia_en_espan
a_2013.pdf
[4] IDAE, “Mapa tecnológico Movilidad Eléctrica”, [en línea] [Consulta: 31 de mayo de
2013] Disponible en web:
http://www.idae.es/index.php/mod.documentos/mem.descarga?file=/documentos_Movi
lidad_Electrica_ACC_c603f868.pdf
[6] RED ELECTRICA DE ESPAÑA, “Bolentín mensual Septiembre 2014”, [en línea]
[Consulta: 2 de noviembre de 2014] Disponible en web:
http://www.ree.es/sites/default/files/downloadable/ree_septiembre_2014.pdf
[12] CHAN, C.C., “The State of the Art of Electric, Hybrid, and Fuel Cell Vehicles”,
Proceedings of the IEEE, Abril 2007, vol.95, no.4, p.704-718
[13] LI-QIANG JIN; XIA-HUA ZENG; WEI WANG, “The control strategy and cost
analysis for series Plug-in hybrid electric vehicle”, Advanced Computer Control (ICACC), 2010
2nd International Conference on, marzo 2010, vol.2, p.350-354
[14] QINGSONG LI; ZHENG ZHANG; WEIJUN ZHAO, “Dynamic Control for
Four-Wheel Independent Drive Electric Vehicle”, Computer Science and Electronics Engineering
(ICCSEE), 2012 International Conference on, Marzo 2012, vol.3, p.252-256
[15] YIMIN GAO; EHSANI, M., “Design and Control Methodology of Plug-in Hybrid
Electric Vehicles”, Industrial Electronics, IEEE Transactions on, Febrero 2010, vol.57, no.2,
p.633-640
[18] LARMINIE, J.; LOWRY, J., “Electric Vehicle technology explained, 2nd Edition”,
Wiley, ISBN: 978-1-119-94273-3
Julio 2008, [en línea] [Consulta: 4 de febrero de 2014]. Disponible en web: http://e-
archivo.uc3m.es/handle/10016/7260
[21] TROJAN BATTERY COMPANY, “Guía del usuario”, [en línea] [Consulta: 10
de octubre de 2013] Disponible en web:
http://www.trojanbattery.com/pdf/TRJN0147_UsersGuideSP_SPA.pdf
[23] PANDAY, A.; BANSAL, H.O., “Temperature dependent circuit-based modeling of high
power Li-ion battery for plug-in hybrid electrical vehicles”, Advances in Technology and
Engineering (ICATE), 2013 International Conference on, 23-25 Enero 2013, p.1-6
[29] PROCHAZKA, P.; PAZDERA, I.; VOREL, P.; CERVINKA, D., “Design of
small electric car”, Power Electronics, Electrical Drives, Automation and Motion
(SPEEDAM), 2012 International Symposium on, Junio 2012, p.359-364
[31] MAS, F.B.; OLIVER, B.; ISERN, E.; PICOS, R. y ROCA, M., “Monitorización
de vehículos eléctricos para el desarrollo de modelos de consumo”, Seminario Anual de Automática,
Electrónica Industrial e Instrumentación, SAAEI, Julio 2012, p. 482-487
[32] HABIB ULLAH, M.; GUNAWAN, T.S.; SHARIF, M. R.; MUHIDA, R.,
“Design of environmental friendly hybrid electric vehicle”, Computer and Communication
Engineering (ICCCE), 2012 International Conference on, Julio 2012, p.544-548
[33] IBRAHIM, S.; ASFANI, D.A.; TAKASHI, H., “Simulation-based analysis of short
circuit fault in parallel-series type hybrid electric vehicle”, Advanced Power System Automation and
Protection (APAP), 2011 International Conference on, Octubre 2011, vol.3, p.2045-2049
[34] YUAN CHENG; CHEN, K.; CHAN, C.C.; BOUSCAYROL, A.; SHUMEI
CUI, “Global modeling and control strategy simulation”, Vehicular Technology Magazine, IEEE,
Junio 2009, vol.4, no.2, p.73-79
[36] JUNJIE LI; XIUMIN YU; HUAJIE DING; PING SUN, “Design And Control
Of The Propulsion For A Series Hybrid Electric Vehicle”, Intelligent Control and Automation
(WCICA), 2010 8th World Congress on, Julio 2010, p.2159-2162
[38] RAZAVIAN, R.; AZAD, N.L.; MCPHEE, J., “On real-time optimal control of a
series Hybrid Electric Vehicle with an ultra-capacitor”, American Control Conference (ACC),
2012, Junio 2012, p.547-552
Modelado del VE por medio de Matlab/Simulink y María Susana Alegre Buj
planificación de EC mediante SIG y algoritmos genéticos
Capítulo 9 Referencias 286
[39] EVANGELOU, S.A.; SHUKLA, A., “Advances in the modelling and control of series
hybrid electric vehicles”, American Control Conference (ACC), 2012, Junio 2012, p.527-534
[40] MILLNER, A.; JUDSON, N.; REN, B.; JOHNSON, E.; ROSS, W., “Enhanced
plug-in hybrid electric vehicles”, Innovative Technologies for an Efficient and Reliable
Electricity Supply (CITRES), 2010 IEEE Conference on, Septiembre 2010, p.333-340
[42] KNOEDLER, K.; STEINMANN, J.; JONES, S.; HUSS, A.; KURAL, E.;
BRINGMANN, O.; ZIMMERMANN, J.; LAVERSANNE, S.; SANCHEZ, D., “Optimal
energy management and recovery for FEV”, Design, Automation & Test in Europe Conference
& Exhibition (DATE), 2012, Marzo 2012, p.683-684
[43] VASILJEVIC, G.; VRHOVSKI, Z.; BOGDAN, S., “Dynamic modeling and
simulation of a three-wheeled electric car”, Electric Vehicle Conference (IEVC), 2012 IEEE
International, Marzo 2012, p.1-8
[44] CARDELLI, E.; MORETTINI, V., “Modeling and simulation of a retrofitted electric
car in urban and extra urban driving cycles”, Energy Conference and Exhibition
(ENERGYCON), 2012 IEEE International, Septiembre 2012, p.1065-1070
[46] MEHAR, S.; SENOUCI, S.M., “An optimization location scheme for electric charging
stations”, Smart Communications in Network Technologies (SaCoNeT), 2013 International
Conference on, vol.01, Junio 2013, p.1-5
[47] JUN HE; BUXIANG ZHOU; CHAO FENG; HENGXIN JIAO; JINHUA
LIU, “Electric Vehicle Charging Station Planning Based on Multiple-Population Hybrid Genetic
Algorithm”, Control Engineering and Communication Technology (ICCECT), 2012
International Conference, Diciembre 2012, p.403-406
[49] SHAOYUN GE; LIANG FENG; HONG LIU, “The planning of electric vehicle
charging station based on Grid partition method”, Electrical and Control Engineering (ICECE),
2011 International Conference, Septiembre 2011, p. 2726-2730
[52] KADAR, P.; VARGA, A., “PhotoVoltaic EV charge station”, Applied Machine
Intelligence and Informatics (SAMI), 2013 IEEE 11th International Symposium on,
Enero-Febrero 2013, p.57-60
[53] RED ELÉCTRICA DE ESPAÑA, “El vehículo eléctrico Una prometedora alternativa
de futuro”, [en línea] [Consulta: 9 de junio de 2012]. Disponible en web:
http://www.ree.es/operacion/pdf/FolletoCocheElectrico.pdf
[54] ROIG ESTRADA, MIGUEL, “Perspectiva actual de la tecnología del coche eléctrico.
Análisis de los cambios en los procesos de producción con la llegada del vehículo eléctrico”, Sancho
Rodríguez, José y Berg, Christina, Universidad Pública de Navarra. Departamento de
Ingeniería Mecánica, Energética y de Materiales, Febrero 2012, [en línea] [Consulta: 12 de
mayo de 2013] Disponible en web: http://hdl.handle.net/2454/6668
[56] EHSANI, MEHRDAD; GAO, YIMIN; EMADI, ALI; “Modern Electric, Hybrid
Electric, and Fuel Cell Vehicles Fundamentals, Theory and Design 2nd Edition”, CRC Press Taylor &
Francis Group, ISBN: 978-1-4200-5398-2