Mecatronica-Jorge Pablo Diaz Velilla
Mecatronica-Jorge Pablo Diaz Velilla
Mecatronica-Jorge Pablo Diaz Velilla
DESARROLLO DE UN PROTOTIPO
MECATRÓNICO DE SEGUIDOR SOLAR
PARABÓLICO PARA LA EDIFICACIÓN
TESIS DOCTORAL
(Ingeniero Industrial)
2018
Departamento de Tecnología de la Edificación
DESARROLLO DE UN PROTOTIPO
MECATRÓNICO DE SEGUIDOR SOLAR
PARABÓLICO PARA LA EDIFICACIÓN
(Ingeniero Industrial)
(Doctor en CC Físicas)
2018
Tribunal nombrado por el Magfco. y Excmo. Sr. Rector de la Universidad
Politécnica de Madrid, el día de de 2018
EL TRIBUNAL CALIFICADOR
Presidente:
Vocal:
Vocal:
Vocal:
Secretario:
Calificación:
El secretario
DESARROLLO DE UN PROTOTIPO
MECATRÓNICO DE SEGUIDOR SOLAR
PARABÓLICO PARA LA EDIFICACIÓN
Junio de 2018
A todos aquellos que creen en las energías
renovables y el desarrollo sostenible como
solución al problema energético y
medioambiental de nuestro planeta
RESUMEN
Son muchísimas las actividades que el ser humano desarrolla, que generan esta
insostenible contaminación medioambiental: amén de la producción energética “clásica”,
son también responsables el transporte, fabricación de bienes de consumo, minería,
telecomunicaciones, industrias alimenticias… y edificación. En efecto, se estima que en
España es responsable del 33% de consumo de la energía anual generada (imputando
la necesaria para fabricar materiales de construcción precisos para el sector y energía
térmica y eléctrica empleada por hogares y edificio). Y dicho sector es responsable, con
el actual “mix” energético nacional, del 30% de emisiones nocivas. A resultas de lo
anterior parece imprescindible que las sociedades avanzadas inicien una transición hacia
un modelo económico de crecimiento basado en el desarrollo sostenible, donde los
dos pilares sean la eficiencia energética y las energías renovables. En añadidura, están
apareciendo nuevos actores que cambiarán en las próximas décadas el paradigma
energético a escala mundial: generación distribuida, edificios inteligentes, viviendas de
consumo casi nulo y procesado masivo de datos. Un apasionante desafío para el futuro.
Energy production since the Industrial Revolution has focused on the use of fossil
fuels (coal, oil and natural gas). On the other hand, the emerging powers (China, India
and Brazil), eager to increase their economies and replicating the growth model of
industrialized countries, are equally intensive in the consumption of previous raw
materials. This inevitably creates two serious problems:
• The progressive depletion of oil and natural gas fields (estimated life of 50 years)
and coal (150 years).
• The unbridled emissions of greenhouse gases (CO2 and CH4) that go hand in
hand with the burning of fuels (particles present in the atmosphere in
concentrations of 413 ppm, already responsible for an increase of 1ºC in
atmospheric temperature). It is sadly known as climate change.
There are a great many activities that people carry out that generate this
unsustainable environmental pollution: in addition to "classic" energy production,
transport, the manufacture of consumer goods, mining, telecommunications, the food
industry... and building are also responsible. In fact, it is estimated that in Spain it is
responsible for 33% of the annual energy consumption generated (imputed to the energy
needed to manufacture precise building materials for the sector and thermal and electrical
energy used by households and buildings). And this sector is responsible, with the current
national energy mix, for 30% of harmful emissions. As a result, it seems imperative that
advanced societies begin a transition to an economic model of growth based on
sustainable development, where the two pillars are energy efficiency and renewable
energy. In addition, new players are emerging that will change the global energy paradigm
in the coming decades: distributed generation, intelligent buildings, near-zero
consumption housing and massive data processing. An exciting challenge for the future.
As with the existing research projects in environmental robotics and the confluence
of the above, this Doctoral Thesis aims to design, simulate, manufacture, validate and
experiment with a mechatronic prototype of a paraboloidal and parabolic trough solar
thermal tracker based on Arduino, capable of capturing and storing operational and
meteorological variables. At the same time, a meteorological station will be built
specifically for the purpose of establishing measurement standards for the magnitudes
measured on the first measurement. It will be necessary to take into account the current
state of the art in concentrating solar energy technology, the use of conical and quadric
optical properties, spherical trigonometry, the fundamentals of solar concentration and
the tracking of single and dual-axis trajectories. All this with the specific objective of
analysing the validity of the prototype as a tool for acquiring useful data for technicians
and engineers, both to improve the meteorological value tables and to be used as training
material. As a general objective, an ambitious one can be set: to serve as a modest
contribution to the scientific-technical community in order to speed up the transfer of
knowledge of concentrated solar thermal energy from current industrial plants to the
Building sector.
1. INTRODUCCIÓN ..............................................................................................................................2
1.1 Dedicatoria y agradecimientos .....................................................................................................2
1.2 Prolegómenos ..............................................................................................................................3
1.3 Objetivos y alcance de la Tesis ....................................................................................................5
1.4 Estructura de la Tesis ...................................................................................................................6
LISTA DE FIGURAS CAPÍTULO 1 …………………………………………………………………………. ..8
ÍNDICE CAPÍTULO 1
1. INTRODUCCIÓN ....................................................................................................................... 2
1.1 Dedicatoria y agradecimientos ........................................................................................... 2
1.2 Prolegómenos .................................................................................................................... 3
1.3 Objetivos y alcance de la Tesis .......................................................................................... 5
1.4 Estructura de la Tesis......................................................................................................... 6
LISTA DE FIGURAS CAPÍTULO 1 ................................................................................................ 8
1
Jorge Pablo Díaz Velilla Introducción
1. INTRODUCCIÓN
2
Jorge Pablo Díaz Velilla Introducción
1.2 Prolegómenos
En el mes de octubre del año 2013, el alumno autor de esta Tesis tuvo la
oportunidad de defender el ya extinto DEA (Diploma de Estudios Avanzados) en la
Escuela Técnica Superior de Ingenieros Industriales de la Universidad Politécnica
de Madrid. En aquel entonces, los objetivos y alcance del proyecto DEA eran:
3
Jorge Pablo Díaz Velilla Introducción
4. Especial motivación personal por aunar en una sola, las dos líneas
tecnológicas que siempre me fueron más atractivas como ingeniero: la
Robótica y la Energía Solar. Ambas han confluido en este trabajo y se han
fundido en el dispositivo mecatrónico ya citado. Es necesario definir la
Mecatrónica como aquella línea de conocimiento multidisciplinar que
engloba ciencias como la Mecánica, Electrónica, Control y Computación.
4
Jorge Pablo Díaz Velilla Introducción
5
Jorge Pablo Díaz Velilla Introducción
6
Jorge Pablo Díaz Velilla Introducción
7
Jorge Pablo Díaz Velilla Introducción
8
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
ÍNDICE CAPÍTULO 2
9
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
10
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
11
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
La ratio medio de producción oscila entre 6,8 y 7,2 barriles de petróleo por
tonelada, dependiendo de la densidad del petróleo. Por tanto, las reservas de crudo
se calculan entre 0,97 y 1,003 billones de barriles. Sin embargo, el límite real de las
reservas podría estar más cercano aún si se consideran modelos de previsión de
la demanda con un consumo creciente, como de hecho ha venido ocurriendo a lo
largo de todo el siglo pasado. Adicionalmente, los descubrimientos de nuevos
yacimientos se han minorado abruptamente en las últimas décadas haciendo
prácticamente insostenible por largo tiempo, los altísimos niveles de extracción
actuales, máxime sin incluir la ya mentada demanda de los consumidores asiáticos.
12
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
13
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
El año exacto del pico no podrá determinarse hasta que ya haya sucedido y
se tenga un horizonte temporal suficientemente amplio para validar la hipótesis, es
decir, se sabrá a posteriori. De hecho, basándose en datos actuales de producción,
la Asociación para el Estudio del Pico del Petróleo y el Gas (ASPO en sus siglas
inglesas), estimó que el pico del petróleo ocurrió ya en el año 2010, siendo el del
gas natural algunos años posterior. Este hecho implicaría importantísimas
consecuencias para los países desarrollados, que dependen en gran medida (en
exceso en nuestra modesta opinión) de petróleo barato y abundante, especialmente
para el sector del transporte, agricultura, industria química y la calefacción
doméstica. No en vano, un cada vez mayor número de expertos no solo estima sino
que afirma sin ambages que el pico de producción, de hecho, ya llegó. Después de
los devastadores efectos del huracán Katrina, Arabia Saudí reconoció su
imposibilidad para incrementar su producción y así atenuar la crisis por las pérdidas
en la producción y el refino sufridas en la zona del Golfo de Méjico. Por eso, son
muchos los que piensan que estamos inmersos en el inicio de la crisis definitiva del
petróleo. Y se puede afirmar que es irreversible porque será la última y la que
obligará a efectuar mayores ajustes y recortes en su consumo como nunca antes
se ha hecho.
14
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
Figura 2.3. Países con mayores emisiones totales de CO2 (fuente: WRI/CAIT, CDIAC)
15
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
16
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
para evitar que a través del transporte que suponen las corrientes atmosféricas, se
perjudique a países limítrofes (contaminación transfronteriza) que nada han tenido
que ver en la producción de tales contaminantes. La tendencia en este sentido, será
la de quemar únicamente combustibles fósiles de bajo nivel de acidez, aunque lo
más sensato es abogar por un abandono progresivo de cualquier combustible fósil
en beneficio de energías renovables respetuosas con el medio ambiente.
17
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
18
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
Figura 2.5. Huella ecológica provocada por los distintos países (fuente: Wikipedia)
20
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
Así por ejemplo, bajo el prisma de un tipo de vida básicamente agraria, bien
organizada y sin monocultivos extensivos, se estima que entre 1 y 2 hectáreas es
aproximadamente el terreno preciso para atender a las necesidades de una familia
autosuficiente. Por contra, teniendo en cuenta el desarrollo industrial físicamente
fundamentado en la disponibilidad de combustibles fósiles y partiendo de
necesidades típicas de un ciudadano urbano medio se llega a la siguiente
conclusión: serían necesarios otros 2,5 planetas como el nuestro para que los 7.000
millones de seres humanos actuales pudieran vivir todos bajo ese tipo de vida.
Estas contundentes conclusiones, hacen necesario distinguir dos elementos clave:
• La huella ecológica poco tiene que ver con el espacio físico ocupado por un
grupo humano.
Figura 2.6. Huella ecológica promedio por países (fuente: State of Ontario's Biodiversity)
21
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
22
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
23
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
24
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
26
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
Por tanto, parece que únicamente puede ser una solución adicional para
mitigar la emisión de CO2 pero en modo alguno es la panacea o respuesta a todos
los problemas medioambientales. En este sentido, una de las técnicas más
desarrolladas hasta ahora en materia de CCS es el llamado almacenamiento
subterráneo en capas acuíferas salinas profundas; en yacimientos de petróleo y/o
gas natural ya agotados; en minas de carbón y en otras formaciones geológicas.
Algunas asociaciones conservacionistas advierten del peligro y error que supondría
apostar únicamente por la captura y secuestro de carbono en detrimento de las
energías renovables, ya que es una tecnología que no conduce a una reducción
real de emisiones de CO2 ni a un cambio del modelo energético basado en
combustibles fósiles. Además, utilizar la CCS significa trasladar inexorablemente a
las generaciones futuras el problema de qué hacer con el CO2 que hoy
almacenaríamos en cuevas subterráneas, de forma similar a como ocurre con los
residuos radiactivos. Desgraciadamente, nunca se podrá eliminar la posibilidad de
fugas a la atmósfera, ya sean lentas filtraciones a largo plazo como escapes
catastróficos a corto plazo, que podrían poner en peligro vidas humanas, echar a
perder depósitos de agua dulce, restringir usos agrícolas de los suelos, contaminar
capas acuíferas y superficiales y actuar brusca y catastróficamente sobre el clima.
27
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
28
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
29
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
30
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
se estima que el uso de dicha energía per capita en EE.UU., Europa o Australia sea
dos a tres veces mayor que el de China o India, por lo que los países de la
Organización para la Cooperación y el Desarrollo Económico (OCDE) deberán
reducir sus emisiones de dióxido de carbono antes que las economías en vías de
industrialización. Ya es sabido qué reducciones se esperan para los países
desarrollados. Para los llamados países en vías de desarrollo, como China e India,
se deben estabilizar sus emisiones de CO2 para 2020 y comenzar a reducirlas a
partir de 2030, creando a la vez un suministro energético asequible y seguro. Y algo
muy importante, manteniendo un desarrollo económico constante y renovable.
Parece por tanto, que se impone un uso y producción inteligente de la energía.
• La reforma del sector eléctrico, teniendo en cuenta las propuestas del Libro
Blanco para orientarlo hacia la sostenibilidad y la renovabilidad.
32
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
• Misma demanda eléctrica per capita en todas las provincias de una misma
comunidad autónoma.
33
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
• Red Natura 2000: Zonas de Especial Protección para las Aves (ZEPA) +
Lugares de Interés para la Conservación (LIC).
Pues bien, los principales resultados del estudio para cada una de las
tecnologías consideradas en comparación con la demanda eléctrica proyectada
para 2050, adoptando un enfoque conservador (es decir, considerando la mejor
tecnología actualmente existente para la transformación de cada uno de los
recursos renovables en electricidad, incluyendo mejoras tecnológicas sólo cuando
parezca obvio que estarán disponibles para el año 2050) arrojan resultados muy
esperanzadores. En global, se puede afirmar que existen multitud opciones de
configurar un “mix” de generación de electricidad 100% renovable incluso con
exceso de producción renovable.
34
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
35
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
Figura 2.12. Defensa de Siracusa por Arquímedes con energía solar (fuente: Wikipedia)
36
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
A partir del año 1950 con el desarrollo de la llamada “carrera espacial” entre
la antigua Unión Soviética y los Estados Unidos, los ingenieros de la época se
percatan de que tienen que dotar a los satélites geoestacionarios de una fuente de
energía autónoma y fiable. Retoman entonces el antedicho efecto fotovoltaico,
dando vida a la primera aplicación industrial de la energía fotovoltaica: la
incorporación de paneles solares en tales satélites para dotarles de electricidad y
así habilitar el funcionamiento de sus componentes electrónicos.
37
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
Figura 2.14. Gran instalación de colectores de energía solar térmica (fuente: autor)
38
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
39
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
𝑇𝑓
𝜂 = 1− (1)
𝑇𝑐
donde:
40
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
41
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
42
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
Dicho motor es propulsado por el proceso de expansión del gas cuando éste
se calienta, seguido de una compresión del mismo cuando se enfría. El motor
contiene una cantidad de gas fija que trasiega constantemente entre la parte
fría y la caliente, gracias a una fuente térmica exterior. Desde el punto de
vista termodinámico, el ciclo Stirling es un ciclo cerrado formado por dos
procesos isotermos (a temperatura constante) y dos isócoros (a volumen
constante).
43
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
• La eficiencia teórica del ciclo es la misma que la del ciclo de Carnot, máxima
eficiencia teórica termodinámicamente hablando, entre dos fuentes térmicas.
• Describe un ciclo cerrado, con lo cual utiliza la misma sustancia tras cada
ciclo, pudiendo utilizarse un regenerador para intercambio térmico.
44
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
45
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
46
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
1
𝛼 𝐺 ∙ 𝑐𝑜𝑠𝜃 4 (2)
𝑇=( ∙ )
𝜀 𝜎
donde:
47
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
• Superficies cilindro-parabólicas.
• Superficies paraboloidales.
𝐴𝑎
𝐶= (3)
𝐴𝑟
donde:
Figura 2.20. Parámetros relevantes del factor de concentración (fuente: ETSII UPM)
48
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
1
𝛼𝑒𝑓 𝐺 ∙ 𝐶 4 (4)
𝑇=( ∙ )
𝜀𝑒𝑓 𝜎
donde:
𝐴𝑎 2𝐷 1
2𝐷
𝐶𝑚𝑎𝑥 =( ) = ≈ 212 (5)
𝐴𝑟 𝑠𝑒𝑛𝜃𝑠
𝐴𝑎 3𝐷 1
3𝐷
𝐶𝑚𝑎𝑥 =( ) = ≈ 45.000 (6)
𝐴𝑟 𝑠𝑒𝑛𝜃𝑠
49
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
50
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
51
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
52
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
53
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
54
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
55
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
56
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
57
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
58
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
BIBILIOGRAFÍA CAPÍTULO 2
[1]. Bayrak F.; Abu-Hamdeh N.; A. Alnefaie K.; F. Öztop H. A review on exergy
analysis of solar electricity production. Renewable and Sustainable Energy
Reviews, Volume 74, July 2017, pp. 755-770.
[3]. Shahbaz M.; Sarwar S.; Chen W.; Nasir M. M. Dynamics of electricity
consumption, oil price and economic growth: Global perspective. Energy
Policy, Volume 108, September 2017, pp. 256-270.
[4]. Luz-Sant’Ana I.; Román-Román P.; Torres-Ruiz F.; Modeling oil production and
its peak by means of a stochastic diffusion process based on the Hubbert curve.
Energy, Volume 133, 15 August 2017, pp. 455-470.
[6]. Bakhshipour Z.; Asadi A.; B.K. Huat B.B.K.; Sridharan A.; Kawasaki S. Effect
of acid rain on geotechnical properties of residual soils. Soils and Foundations,
Volume 56, Issue 6, December 2016, pp. 1008-1020.
[7]. Akitt J.W. Some observations on the greenhouse effect at the Earth's surface:
Reflection of energy near the earth's surface: A bounded random walk
approach. Spectrochimica Acta Part A: Molecular and Biomolecular
Spectroscopy, In press, accepted manuscript, Available online 5 July 2017.
[8]. Zhang L.; Dzakpasu M.; Chen R.; Wang X.C. Validity and utility of ecological
footprint accounting: A state-of-the-art review. Sustainable Cities and Society,
Volume 32, July 2017, pp. 411-416.
59
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
[13]. Landrum A.R.; Lull R.B.; Akin H.; Hasell A.; Jamieson K.H. Processing the
papal encyclical through perceptual filters: Pope Francis, identity-protective
cognition, and climate change concern. Cognition, Volume 166, September
2017, pp. 1-12.
[14]. Zou C.; Zhao Q.; Zhang G.; Xiong B. Energy revolution: From a fossil energy
era to a new energy era. Natural Gas Industry B, Volume 3, Issue 1, January
2016, pp. 1-11.
[15]. García X.; Domínguez J.; Linares P.; López O. Renovables 2050. Un informe
sobre el potencial de las energías renovables en la España peninsular. Instituto
de Investigaciones Tecnológicas, Universidad Pontificia Comillas. Madrid.
España. Julio 2005.
[16]. Bukarica V.; Tomšić Ž. Energy efficiency policy evaluation by moving from
techno-economic towards whole society perspective on energy efficiency
market. Renewable and Sustainable Energy Reviews, Volume 70, April 2017,
pp. 968-975.
[19]. Anderson W.W.; Chai Y.G. Becquerel effect solar cell. Energy Conversion,
Volume 15, Issues 3–4, 1976, pp. 85-94.
[21]. S.E. Wright. A generalized and explicit conceptual statement of the principle
of the second law of thermodynamics. International Journal of Engineering
Science, Volume 111, February 2017, pp. 12-18.
[22]. Canturk B.; Oikonomou T.; Bagci B.G. Group theory, entropy and the third
law of thermodynamics. Annals of Physics, Volume 377, February 2017, pp.
62-70.
[23]. Helvaci H.U.; Khan Z.A. Thermodynamic modelling and analysis of a solar
organic Rankine cycle employing thermofluids. Energy Conversion and
Management, Volume 138, 15 April 2017, pp. 493-510.
[24]. Hafez A.Z.; Soliman A.; El-Metwally K.A.; Ismail I.M. Solar parabolic dish
Stirling engine system design, simulation, and thermal analysis. Energy
Conversion and Management, Volume 126, 15 October 2016, pp. 60-75.
60
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
[27]. Li Q.; Shirazi A.; Zheng C.; Rosengarten G.; Taylor R.A. Energy
concentration limits in solar thermal heating applications. Energy, Volume 96,
1 February 2016, pp. 253-267.
[28]. Criado C.; Alamo N. Optical properties of conics: a method for obtaining
reflecting and focusing profiles. Optics Communications, Volume 167, Issues
1–6, 15 August 1999, pp. 83-88.
[29]. Kim D-S.; Kim Y.H. Some characterizations of spheres and elliptic
paraboloids. Linear Algebra and its Applications, Volume 437, Issue 1, 1 July
2012, pp. 113-120.
[30]. Rostamzadeh H.; Ebadollahi M.; Ghaebi H.; Amidpour M.; Kheiri R. Energy
and exergy analysis of novel combined cooling and power (CCP) cycles.
Applied Thermal Engineering, Volume 124, September 2017, pp. 152-169.
[31]. Mahoney A.R.; Price H. Solar Field Performance of New UVAC Receivers at
SEGS VI. KJC Operation Co. 2002. ASME Journal of Solar Energy
Engineering.
[32]. Beltagy H.; Semmar D.; Lehaut C.; Said N. Theoretical and experimental
performance analysis of a Fresnel type solar concentrator. Renewable Energy,
Volume 101, February 2017, pp. 782-793.
[33]. Mendoza L.S.; Carrillo G.E.; Melian V.R.; Silva E.E.; Martinez A.M.
Mathematical modeling of the geometrical sizing and thermal performance of a
Dish/Stirling system for power generation. Renewable Energy, Volume 107,
July 2017, pp. 23-35.
[34]. Thalange V.C.; Dalvi V.H.; Mahajani S.M.; Panse S.V.; Patil R.N. Design,
optimization and optical performance study of tripod heliostat for solar power
tower plant. Energy, Volume 135, 15 September 2017, pp. 610-624.
[35]. Prieto C.; Miró L.; Peiró G.; Oró E.; Cabeza L.F. Temperature distribution and
heat losses in molten salts tanks for CSP plants. Solar Energy, Volume 135,
October 2016, pp. 518-526.
[36]. Relloso S.; García E. Tower Technology Cost Reduction Approach after
Gemasolar Experience. Energy Procedia, Volume 69, May 2015, pp. 1660-
1666.
61
Jorge Pablo Díaz Velilla Fundamentos de energía termosolar
62
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
ÍNDICE CAPÍTULO 3
63
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
360 ∙ 𝑛
𝐺 = 𝐺𝑐𝑡𝑒 ∙ [1 + 0,033 ∙ cos ( )] (1)
365
donde:
64
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Entre los datos más relevantes del Sol, se estiman que los más relevantes
son los que se indican seguidamente: edad (4,5∙109 años); masa (1030 kg); radio
(696∙103 km); densidad (1,4∙103 kg/m3); distancia promedio Sol-Tierra (150∙106 km);
pérdida de masa (09 kg/s); radiación emitida (3,86∙1026 W). Ahora bien, de esa
energía incidente sobre la Tierra antes de atravesar la atmósfera, realmente a la
superficie de esta llega una cuantía menor, aproximadamente 1.000 W/m2
dependiendo de las latitudes. Dicha merma energética es debida a varios factores:
donde:
65
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
66
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Esto es, ocurre como si el Sol fuese una gigantesca fuente de radiación
electromagnética cuyo poder emisivo fuese el de un cuerpo negro a una
temperatura aproximada de 5.777 K (concepto denominado en la Ciencia como
cuerpo negro). Una ínfima porción de esa energía es la recibida por la Tierra. En
este sentido, se conoce como espectro electromagnético (Norgard J., 2007) de
la radiación como la distribución energética de ondas electromagnéticas en su
conjunto. Si se refiere a un cuerpo concreto se denomina simplemente espectro a
la radiación electromagnética que emite (espectro de emisión) o absorbe (espectro
de absorción) la sustancia de la cual está constituida dicho cuerpo.
67
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
68
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
1
𝑚= (3)
cos 𝜃𝑧
donde:
Los dispositivos más utilizados para medir la radiación solar en todos sus formatos
así como otras variables asociadas, son los que a continuación se detallan:
69
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Pues bien, todos los proyectos llevados a cabo hasta la fecha en nuestro
país para caracterizar la radiación solar, se centran en cuantificar y representar
gráficamente, con suficiente resolución espacial, el promedio tanto mensual como
anual de la radiación que llega a la superficie terrestre española. Así, el primigenio
Atlas de Radiación Solar del ya extinto Instituto Nacional de Meteorología, databa
del año 1984 y fue construido mediante una red de radiómetros en superficie, que
aunque intentaba “mapear” todo el territorio nacional, no eran suficientes para
hacerlo con suficiente precisión amén de que se elaboró con tecnologías muy
precarias comparadas con las disponibles hoy en día.
Desde entonces hasta nuestros días, han surgido en España varios atlas
monográficos sobre radiación solar, especialmente desde el advenimiento político
de las comunidades autónomas: como referencia, cabe destacar el atlas basado en
imágenes de satélite Meteosat 6 para Galicia, confeccionado con datos del período
2002-2004 junto con el método de estimación de irradiancias Heliosat-2.
70
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
71
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Figura 3.9. Mapa de radiación: España y Portugal (fuente: Institute for European Studies)
72
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
𝑅 − 𝑅𝑠𝑓𝑐
𝐶𝐴𝐿 = (4)
𝑅𝑚𝑎𝑥 − 𝑅𝑠𝑓𝑐
donde:
73
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
𝑆𝐼𝑆
𝑘= (5)
𝑆𝐼𝑆𝑐𝑙𝑠
donde:
donde:
donde:
SIS: irradiancia total que incide sobre la superficie horizontal dentro de la banda
espectral 0,2 - 4 μm o Surface Incoming Solar radiation, en W/m2
E0: irradiancia solar en la cima de la atmósfera (constante solar), en W/m2
θ: ángulo zenital, en rad
T: transmisividad atmosférica adimensional
74
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
75
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
284 + 𝑛
𝛿 = 23,45 ∙ sin (360 ∙ ) (9)
365
donde:
• Azimut solar β o γs: ángulo formado por la recta que une al observador con
el sur geográfico (S), con la recta que le une a la proyección ortogonal del
Sol sobre la horizontal. Está comprendido entre [-180º, +180º].
76
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
• Elevación solar α: ángulo formado por la recta que une al observador con
la posición del Sol, con la recta que le une a la proyección ortogonal del Sol
sobre la horizontal.
77
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
donde:
78
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Figura 3.16. Gráfico 3D del posicionamiento solar para latitud 40,4º (fuente: Sunpath3D)
Figura 3.17. Gráfico 2D del posicionamiento solar para latitud 40,4º (fuente: Solarpraxis)
cos 𝜃 = cos 𝛿
tan 𝜑 ( 14 )
tan 𝛽 =
cos 𝛾
cos 𝜃 = 1
𝜋 ( 15 )
𝛽 = − 𝜃𝑧
2
En cuanto a las tecnologías de seguimiento algorítmico solar, básicamente
las dos tendencias actuales son la que se indican seguidamente:
80
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
81
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Figura 3.20. Discos Stirling con seguidor a dos ejes (fuente: Solar-energy)
82
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
83
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
84
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
85
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
86
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Figura 3.23. Diseño 3D de Rover Curiosity para el análisis de Marte (fuente: NASA)
87
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Entre las tareas más inmediatas en las que la robótica implementa mejoras
en la construcción, se pueden enumerar un sinfín de ellas: automatización del
proceso constructivo; movimiento de tierras; colocación de ladrillos y paneles
aislantes; proyección y pulverización de aislamiento térmico y mortero; pintado de
grandes techos y superficies; instalación de vidrios; soldadura de estructuras
metálicas; fabricación de encofrados; control de calidad de producto terminado;
demolición de construcciones avejentadas y así un largo etcétera.
88
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Figura 3.25. Monitorización robotizada: agua, tuberías y paneles (fuente: ETSII UPM)
89
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
Figura 3.26. Principio topológico de generación distribuida (fuente: Smart grids info)
91
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
92
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
93
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
94
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
95
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
96
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
97
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
El informe del BPIE señala en añadidura, una importante diferencia entre las
estrategias de los Estados Miembros en cuanto a los edificios de nueva
construcción y los ya existentes que precisan una rehabilitación. Así, solo ocho
naciones tienen establecidos requisitos EECN para edificios que sean objeto de
reforma integral. Y si bien es cierto es que aunque todavía no existe una definición
completa de EECN en España, no es correcto afirmar que no se están implantando
estrategias hacia este objetivo ya que la Administración está realizando una
progresiva aproximación al concepto, como demuestra una primera modificación
que tuvo lugar en 2013 del Documento Básico de Ahorro de Energía del Código
Técnico de la Edificación (CTE) que modificó completamente la manera de plantear
criterios energéticos respecto al año 2006 en que fue promulgado dicho CTE. Así,
se ha fijado indicadores válidos para proseguir hacia la anhelada definición de
EECN, incluyendo un nuevo indicador de energía primaria no renovable, así como
un enfoque más prestacional en el cual los técnicos tienen mayor flexibilidad para
implantar soluciones tecnológicas que cumplan los requisitos precisos.
98
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
99
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
BIBILIOGRAFÍA CAPÍTULO 3
[1]. Li H.; Lian Y.; Wang X.; Weibin Ma W.; Zhao L. Solar constant values for
estimating solar radiation. Energy, Volume 36, Issue 3, March 2011, pp. 1785-
1789.
[2]. Vindel J.M.; Polo J.; Zarzalejo L.F. Modeling monthly mean variation of the
solar global irradiation. Journal of Atmospheric and Solar-Terrestrial Physics,
Volume 122, January 2015, pp. 108-118.
[4]. Polo J. Solar global horizontal and direct normal irradiation maps in Spain
derived from geostationary satellites. Journal of Atmospheric and Solar-
Terrestrial Physics, Volumes 130–131, August 2015, pp. 81-88.
[7]. Kittler R.; Darula S. Determination of time and sun position system. Solar
Energy, Volume 93, July 2013, pp. 72-79.
[9]. Wang Y.; Zheng W.; Sun S.; Li L. X-ray pulsar-based navigation system with
the errors in the planetary ephemerides for Earth-orbiting satellite. Advances in
Space Research, Volume 51, Issue 12, 15 June 2013, pp. 2394-2404.
[10]. Kwan T.H.; Wu X. High performance P&O based lock-on mechanism MPPT
algorithm with smooth tracking. Solar Energy, Volume 155, October 2017, pp.
816-828.
[11]. Barbieri G.; Goldoni G.; Borsari R.; Fantuzzi C. Modelling and Simulation for
the Integrated Design of Mechatronic Systems. IFAC-PapersOnLine, Volume
48, Issue 10, 2015, pp. 75-80.
[12]. Cardoso da Silva N.; Maurício Rosário J.; Lima C.R.E. Computer Aided
Choice for DC Motors of Antropomorphics Robots. IFAC Proceedings Volumes,
Volume 33, Issue 26, September 2000, pp. 493-497.
[13]. Hayati S.; Rankin A.; Kim W.; Leger P.; Ali K. Advanced robotics technology
infusion to the NASA Mars exploration Rover (MER) project. IFAC Proceedings
Volumes, Volume 40, Issue 15, 2007, pp. 180-185.
100
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
[15]. Kurkin A.A.; Tyugin D.Y.; Kuzin V.D.; Chernov A.G.; Zeziulin D.V.
Autonomous Mobile Robotic System for Environment Monitoring in a Coastal
Zone. Procedia Computer Science, Volume 103, 2017, pp. 459-465.
[16]. Soheyli S.; Mehrjoo M.; Mayam M.H.S. Modeling and optimal resources
allocation of a novel tri-distributed generation system based on sustainable
energy resources. Energy Conversion and Management, Volume 143, 1 July
2017, pp. 1-22.
[18]. Tummala A.; Velamati R.K.; Sinha D.K.; Indraja V.; Krishna V.H. A review on
small scale wind turbines. Renewable and Sustainable Energy Reviews,
Volume 56, April 2016, pp. 1351-1371.
[19]. Židonis A.; Aggidis G.A. State of the art in numerical modelling of Pelton
turbines. Renewable and Sustainable Energy Reviews, Volume 45, May 2015,
pp. 135-144.
[20]. An L.; Chen R. Mathematical modeling of direct formate fuel cells. Applied
Thermal Engineering, Volume 124, September 2017, pp. 232-240.
[21]. Bayindir R.; Colak I.; Fulli G.; Demirtas K. Smart grid technologies and
applications. Renewable and Sustainable Energy Reviews, Volume 66,
December 2016, pp. 499-516.
101
Jorge Pablo Díaz Velilla Seguimiento solar, robótica de monitorización y edificación sostenible
102
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
ÍNDICE CAPÍTULO 4
103
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Dicho lo anterior, es fundamental definir las fases (Schuh G. et al., 2016) del
diseño mecatrónico del prototipo de seguidor termosolar híbrido paraboloidal y
cilindro-parabólico de esta Tesis, etapas que de forma pormenorizada se describen
a continuación y sobre las cuales gravitan el resto de capítulos del escrito:
104
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
105
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
106
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
107
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
108
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
109
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
En primer lugar y al objeto de construir una idea general del prototipo global
ya diseñado, se incluye una imagen 3D completa de todos los componentes
ensamblados (a excepción del cableado para facilitar la percepción visual del
conjunto). Y progresivamente se aborda el trazado de componentes mecánicos,
ópticos y electrónicos, respectivamente.
Figura 4.7. Diseño 3D completo del prototipo de seguidor termosolar (fuente: autor)
110
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
111
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Information Modeling), en el cual un modelo de este tipo contiene, por así decirlo,
el ciclo de vida completo de cualquier construcción arquitectónica, desde el
concepto inicial hasta la edificación en sí misma. FreeCAD, como otros muchos
modernos modeladores CAD tridimensionales, tiene también un banco de trabajo
para el dibujo bidimensional, muy útil tanto para el dibujo 2D per se, como para su
posterior extrusión hasta lograr sólidos en 3D (Beale R. y Bowers P., 2018).
112
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
113
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
114
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Se adjuntan por tanto una tabla resumen que “traduce” los anteriores
símbolos con los dispositivos reales que el prototipo incorpora, la mayoría de ellos
del fabricante Seeedstudio; las entradas y salidas del Arduino Mega ocupadas; la
variable asignada posteriormente en programación; y el sumatorio de entradas y
salidas analógicas y digitales. De otra parte, se incorpora una imagen con el
esquema electrónico 2D completo de los componentes que conforman el prototipo
(sensores y actuadores) con la nomenclatura gráfica propia de la base de datos de
Fritzing. Merece la pena recalcar que en función de la complejidad en el trazado de
todas las conexiones y al objeto de simplificar las mismas, se ha recurrido a adquirir
dispositivos muy similares a los diseñados pero del fabricante ya mentado
Seeedstudio, familia de dispositivos Grove, que mediante unos ingeniosos
conectores rápidos, permite agilizar notablemente el montaje de los mismos. En el
momento de redacción de estas líneas, no existe en Fritzing una base de datos de
estos dispositivos lo suficientemente completa como para haber recurrido a ella. Y
de otra parte, sigue pareciendo más interesante representar de forma más “purista”
los componentes que el prototipo debe implementar, de cara a una trazabilidad más
clara para un futuro lector de la misma que desee reproducir el experimento.
También se incluye en tercer lugar un mapa electrónico (Alghazzawi T.F., 2016)
de diseño con la simbología estándar para esta línea de conocimiento, para una
rápida interpretación por especialistas en la materia. En último lugar, se incluye un
ágil conteo de las entradas y salidas analógicas y digitales junto con los puertos
que el prototipo precisará para su correcto gobierno.
115
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
116
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
117
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
• Generación del diagrama de flujo del (los) algoritmos así como alternancia entre
flujogramas clásicos o de Nassi-Shneiderman.
118
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
119
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
120
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Figura 4.18. Texto auxiliar de ayuda para interpretación de variables (fuente: autor)
121
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
122
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
123
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
124
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
125
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
126
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Figura 4.34. Movimiento manual en azimut y elevación con joystick (fuente: autor)
127
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
128
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
129
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
130
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
131
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Servo servos[13];
void setup()
{
E=1;
S=0;
132
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Jx=0;
Jy=0;
SN=0;
SS=0;
SE=0;
SO=0;
Sm=(SN + SS + SE + SO) / 4;
Smax=1200;
Smin=200;
alfa=0;
delta_alfa=1;
alfa_repseg=120;
beta=0;
delta_beta=1;
beta_repseg=0;
servo_pausa=10;
Tp=0;
Tcp=0;
Tmax=300;
V=0;
Vmax=100;
C=0;
I=0;
Ta=0;
Hr=0;
A=0;
Pa=0;
pinMode(23,INPUT);
pin_E=digitalRead(23);
pinMode(31,INPUT);
133
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
pin_C=digitalRead(31);
pinMode(47,INPUT);
pin_I=digitalRead(47);
pinMode(11,INPUT);
pin_servo_ALFA=digitalRead(11);
pinMode(7,INPUT);
pin_servo_BETA=digitalRead(7);
Serial.begin(9600);
void loop()
{
Serial.println("Posicionamiento en AZIMUT con joystick");
int *J=(int*)malloc(3*sizeof(int));
J[0]=0;
J[1]=0;
J[2]=0;
pinMode(,INPUT_PULLUP);
J=readJoystick_pin_pin();
if (J[0] >= 512) {
Serial.println("Movimiento del seguidor hacia el E");
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta + delta_beta);
delay(servo_pausa);
}else if (J[0] < 512) {
Serial.println("Movimiento del seguidor hacia el O");
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta - delta_beta);
delay(servo_pausa);
}
Serial.println("Posicionamiento en ELEVACIÓN con joystick");
if (J[1] >= 512) {
Serial.println("Movimiento del seguidor hacia el N");
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa + delta_alfa);
delay(servo_pausa);
}else if (J[1] < 512) {
Serial.println("Movimiento del seguidor hacia el S");
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa - delta_alfa);
delay(servo_pausa);
}
Serial.println("Autómata ENCENDIDO: nueva iteración cada 10 minutos");
Serial.println("Seguimiento en ELEVACIÓN");
Serial.println("Captura radiaciones");
SN=analogRead(pin_SN);
SS=analogRead(pin_SS);
if (SN > SS) {
Serial.println("Movimiento del seguidor hacia el N");
134
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa + delta_alfa);
delay(servo_pausa);
}else {
if (SN < SS) {
Serial.println("Movimiento del seguidor hacia el S");
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa - delta_alfa);
delay(servo_pausa);
}else {
Serial.println("El seguidor permanece con misma elevación");
}
}
while (pin_E == HIGH) {
}
Serial.println("Seguimiento en AZIMUT");
Serial.println("Captura radiaciones");
SE=analogRead(pin_SE);
SO=analogRead(pin_SO);
if (SE > SO) {
Serial.println("Movimiento del seguidor hacia el E");
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta + delta_beta);
delay(servo_pausa);
}else {
if (SE < SO) {
Serial.println("Movimiento del seguidor hacia el O");
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta - delta_beta);
delay(servo_pausa);
}else {
Serial.println("El seguidor permanece con mismo azimut");
}
}
Serial.println("Verificación condiciones REPOSO-SEGURIDAD");
Serial.println("1. Comprobación VELOCIDAD VIENTO");
pinMode(pin_V,INPUT);
V=analogRead(pin_V);
if (V >= Vmax) {
Serial.println("Posición reposo-seguridad ACTIVADA por alta velocidad VIENTO");
S=1;
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa_repseg);
delay(servo_pausa);
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta_repseg);
delay(servo_pausa);
}else {
Serial.println("No hay peligro por velocidad VIENTO");
S=0;
Serial.println("2. Comprobación caída LLUVIA");
if (pin_C == HIGH) {
Serial.println("Posición reposo-seguridad ACTIVADA por caída LLUVIA");
S=1;
alfa=alfa_repseg;
beta=beta_repseg;
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa_repseg);
delay(servo_pausa);
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta_repseg);
delay(servo_pausa);
135
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
}else {
Serial.println("No hay peligro por caída LLUVIA");
S=0;
Serial.println("3. Comprobación impacto GRANIZO");
if (pin_I == HIGH) {
Serial.println("Posición reposo-seguridad ACTIVADA por impacto GRANIZO");
S=1;
alfa=alfa_repseg;
beta=beta_repseg;
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa_repseg);
delay(servo_pausa);
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta_repseg);
delay(servo_pausa);
}else {
Serial.println("No hay peligro por impacto GRANIZO");
S=0;
Serial.println("4. Comprobación niveles MÁXIMOS RADIACIÓN");
if (Sm >= Smax) {
Serial.println("Posición reposo-seguridad ACTIVADA por alta RADIACIÓN");
S=1;
alfa=alfa_repseg;
beta=beta_repseg;
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa_repseg);
delay(servo_pausa);
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta_repseg);
delay(servo_pausa);
}else {
Serial.println("No hay peligro por alta RADIACIÓN");
S=0;
Serial.println("5. Comprobación niveles MINIMOS RADIACIÓN");
if (Sm <= Smin) {
Serial.println("Posición reposo-seguridad ACTIVADA por baja RADIACIÓN");
S=1;
alfa=alfa_repseg;
beta=beta_repseg;
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa_repseg);
delay(servo_pausa);
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta_repseg);
delay(servo_pausa);
}else {
Serial.println("La radiación AÚN es aprovechable");
S=0;
Serial.println("6. Comprobación ALTA TEMPERATURA en colector PARABOLOIDAL");
if (pin_Tp >= Tmax) {
Serial.println("Posición reposo-seguridad ACTIVADA por alta TEMPERATURA colector
PARABOLOIDAL");
S=1;
alfa=alfa_repseg;
beta=beta_repseg;
servos[pin_servo_ALFA].attach(pin_servo_ALFA);
servos[pin_servo_ALFA].write(alfa_repseg);
delay(servo_pausa);
servos[pin_servo_BETA].attach(pin_servo_BETA);
servos[pin_servo_BETA].write(beta_repseg);
delay(servo_pausa);
}else {
136
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
137
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
138
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
139
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
140
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
donde:
qe: presión estática, en kN/m2
qb: presión dinámica del viento, con 0,42 kN/m2 para zona A; 0,45 kN/m2 para
zona B y 0,52 kN/m2 para zona C según mapa de velocidades.
ce: coeficiente de exposición, adimensional y función de la altura y la
aspereza del punto de exposición. Hasta 8 alturas, toma un valor de 2, pero
un cálculo riguroso se realiza según tabla y expresión contiguas.
cp: coeficiente eólico o de presión, adimensional y en tablas según forma y
orientación respecto al viento, pudiendo ser presión o succión.
141
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
𝑐𝑒 = 𝐹 ∙ (𝐹 + 7𝑘) (2)
𝑚á𝑥 𝑍 (3)
𝐹 = 𝑘 ∙ ln ( )
𝐿
donde:
F, k, Z, L: parámetros característicos del entorno según tabla anexa.
142
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
143
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
144
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
145
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Las cargas recibidas por cada perfil en L son las que le transmite el
casquillo a él unido. Ambos perfiles se comportan exactamente igual debido
a la simetría de las cargas, que son exactamente las mismas que recibía el
casquillo correspondiente, más un momento flector debido a la distancia
vertical de 20 mm correspondiente a la longitud del casquillo. Esto produce
un momento en el perfil de valor 2.981,6 N·mm.
Figura 4.52. Esfuerzos sobre casquillos de cruceta (fuente: Julián García de León)
146
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
147
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Pre-procesado
Figura 4.54. Etapas de pre-procesado de diseño CAD (fuente: Julián García de León)
148
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
149
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
150
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Abrazadera de acero: se
observan valores contenidos de la
tensión máxima y la distribución de
tensiones es bastante coherente con
la solicitación y la geometría. La zona
de mayor concentración se produce
justo en el borde de la sujeción de la
arandela y el máximo se encuentra
en la cara orientada al sentido de la
fuerza de tracción (sobre la
superficie curva interna de la
abrazadera). Considerando un acero Figura 4.60. Tensiones von Mises de abrazadera
(fuente: Julián Garcia de León)
comercial genérico fácilmente
localizable en cualquier ferretería (el
material real seleccionado seguirá este criterio) se van a obtener los valores
máximos de las tensiones, como punto de partida de los requisitos que
habría que verificar que cumple la pieza. Se concluye que la pieza va a
soportar adecuadamente las solicitaciones calculadas, pues la tensión
máxima es solo de 148,42 MPa < 240 MPa que puede tomarse como valor
promedio de la tensión máxima admisible para un acero convencional.
Soporte de polietileno:
aunque el programa ofrece un
máximo en un nodo concreto, no hay
diferencia real entre las dos zonas
simétricamente situadas respecto de
la carga (la disposición de los ejes es
arbitraria, el sistema coordenado no
coincidirá con el resto de modelos de
esta pieza debido a que se ha
manipulado la posición de la misma
en el preproceso para facilitar la
aplicación de las cargas, pero esto Figura 4.61. Tensiones von Mises de soporte (fuente:
no tiene la mayor incidencia en los Julián Garcia de León)
cálculos). La distribución de
tensiones que presentan los resultados es coherente con la geometría del
problema. Dadas las restricciones de movimiento consideradas y las cargas
aplicadas, es perfectamente plausible que se produzcan las máximas
tensiones en las zonas con cambio de geometría más abrupta. Con estos
resultados, y teniendo en cuenta las propiedades del material que se escogió
como posibilidad para la pieza real (polietileno), se llega a la conclusión de
que, en las condiciones extremas de carga supuestas para llevar a cabo la
simulación, la pieza fallaría: se supera la tensión de rotura, 57,17 MPa > 33
MPa (resistencia a la tracción del antedicho plástico).
151
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Casquillos de latón: se
aprecia que la mayor tensión
corresponde a la superficie inferior,
donde está unido el casquillo al perfil,
en la dirección de aplicación de la
carga (eje Y). El tornillo de acero
interno únicamente le confiere rigidez
al sistema y reduce las tensiones que
tendría el casquillo si no estuviera
roscado a aquel. Atendiendo a los
valores habituales para el límite
elástico del latón (material estimado Figura 4.62. Tensiones von Mises de casquillo (fuente:
Julián García de León)
para esta pieza) se está
razonablemente lejos de
comprometer la estructura por el mismo pues 238,18 MPa < 360 MPa, siendo
además los desplazamientos pequeños, por lo que puede admitirse un
correcto comportamiento de la pieza.
Perfil de aluminio: se
aprecian zonas de elevadísima
concentración de tensiones,
correspondientes a las zonas de la
ligadura con uno de los tornillos. A la
vista de los resultados y las
propiedades de distintos tipos de
aluminios (75 MPa) se llega a la
conclusión de que la pieza va a fallar
pues 500 MPa > 75 MPa. Pero es
factible pensar que esas tensiones
tan elevadas son debidas a las
Figura 4.63. Tensiones von Mises de perfil de aluminio
singularidades introducidas en el (fuente: Julián García de León)
modelo por la aplicación de las
cargas puntuales en los nodos para simular el sistema equivalente al
momento flector, lo que no deja de ser una burda aproximación. El modelo
parece, a todas luces, insuficientemente preciso. Sin embargo, un dato
importante, las reacciones sí se pueden considerar por buenas, ya que para
este caso las singularidades no influyen. Queda, pues, la opción de mejorar
el modelo en un futuro, como se comentará en las conclusiones. Una de las
opciones más adecuadas a priori sería el modelado mediante la teoría de
láminas.
152
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
153
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Desde finales de la década de los años 80, la industria relacionada con las
tecnologías de concentración ha desarrollado códigos de diseño óptico para el
diseñar y modelar sistemas de esta índole, entre ellos obviamente los de
concentración solar (Jorgensen G., 1994). Estos primigenios códigos fueron
escritos para modelar geometrías ópticas muy particulares de forma que
posteriores adaptaciones fueron construidas a partir de las anteriores de forma
evolutiva. A modo de ejemplos históricos de tales códigos se pueden citar
OPTDSH (concentradores de disco parabólico de apertura circular); ODMF
(concentradores de disco parabólico multifacético); SOLFUR (horno solar).
Conforme la energía termosolar de concentración comenzó a expandirse, el sector
se convenció de que se necesitaba una herramienta de cálculo y diseño más
generalista. Y ese fue justamente el punto de partida de la metodología de trazado
de rayos de Montecarlo como método de simulación en la que están basados los
actuales códigos de programación. No en vano, otros códigos de modelado de
aquel entonces estaban limitados tanto en la definición geométrica como en la
representación óptica de las superficies de reflexión (Wendelin T., 2003). Y en
añadidura, el tiempo y esfuerzo precisos para aprender estos complejos paquetes
informáticos (de entornos de desarrollo y uso poco “amigables”) provocó que
quedaran en desuso. A resultas de todo ello, surgió una nueva generación de
aplicaciones gratuitas y/o de pago tales como Tonatiuh, Comsol Multiphysics (Wave
Optics Module) y SolTrace, entre otros.
Figura 4.65. Haz de rayos por método Montecarlo (fuente: Comsol, Soltrace, Tonatiuh)
154
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
155
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Pues bien, una vez fijados conocidos los orígenes de esta aplicación,
se procede sin más dilación a iniciar el proceso de parametrización de
SolTrace. Dos características definen al Sol como fuente de irradiación: la
distribución angular de la intensidad de radiación a través del disco solar
(denominada forma del Sol) y la posición solar. El primer paso a dar es
determinar la dirección del Sol: la mejor opción que el programa ofrece es
definir un punto en el sistema de coordenadas global de forma que un vector
director desde el origen hasta un punto final, defina la dirección del Sol. De
esta forma SolTrace calcula la posición del Sol en azimut y elevación y
determina un vector unitario correspondiente basado en las ecuaciones de
posicionamiento solar basadas en el tiempo solar que provienen de la
relación geométrica esférica de la Tierra y el Sol, que se recuerdan:
156
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
157
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
2 2
𝜎ó𝑝𝑡𝑖𝑐𝑜 = √4𝜎𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 + 𝜎𝑒𝑠𝑝𝑒𝑐𝑢𝑙𝑎𝑟𝑖𝑑𝑎𝑑 (10)
donde σ se sustituirá por la mitad del ancho del pulso rectangular en el caso
de emplear dicha distribución. Y donde conforme las leyes de Snell, un error
de pendiente de θ se traduce en un error del vector reflejado de 2θ.
158
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
159
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
160
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
161
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
162
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
163
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Figura 4.77. Parametrización trazado rayos sobre colector paraboloidal (fuente: autor)
Una vez realizado el trazado de los haces de rayos para ambos tipos
de colectores, las conclusiones que se pueden obtener en cuanto a
irradiancia concentrada en W/m2 se pueden representan en términos de
trazado de contorno (a modo de visualización termográfica) o bien en
términos de trazado de superficie (a modo de representación tridimensional).
Seguidamente se muestran las salidas por pantalla que el programa
SolTrace genera tras la simulación para el colector cilindro-parabólico y el
paraboloidal, respectivamente.
164
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Figura 4.79. Irradiancia del receptor cilindro-parabólico en trazado contorno (fuente: autor)
Figura 4.80. Irradiancia del receptor cilindro-parabólico en trazado superficie (fuente: autor)
165
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Figura 4.81. Irradiancia del receptor paraboloidal trazado contorno (fuente: autor)
Figura 4.82. Irradiancia del receptor paraboloidal trazado superficie (fuente: autor)
166
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
𝑃 4 𝑃 1
= 𝜎 ∙ 𝑇 4 ⟹ 𝑇 = √( ) ∙ ( ) (11)
𝐴 𝐴 𝜎
donde:
P/A: irradiancia recibida por cada receptor como cuerpo negro, en W/m 2
σ=5,67∙10-8 W/m2K4: constante de Stefan-Boltzmann
T: temperatura de cada receptor, en K
167
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
BIBILIOGRAFÍA CAPÍTULO 4
[1]. Schuh G.; Rudolf S.; Breunig S. Modular Platform Design for Mechatronic
Systems using Axiomatic Design and Mechatronic Function Modules. Procedia
CIRP, Volume 50, 2016, pp. 701-706.
[2]. González-Lluch C.; Company P.; Contero M.; Camba J.D.; Plumed R. A survey
on 3D CAD model quality assurance and testing tools. Computer-Aided Design,
Volume 83, February 2017, pp. 64-79.
[4]. Burleyson C.D.; Nathalie Voisin, Z.; Taylor T.; Yulong X. I.K. Simulated building
energy demand biases resulting from the use of representative weather
stations. Research article. Applied Energy, Volume 209, 1 January 2018, pp.
516-528.
[5]. Camba J.D.; Contero M.; Company P. Parametric CAD modeling: An analysis
of strategies for design reusability. Computer-Aided Design, Volume 74, May
2016, pp. 18-31.
[6]. Beale R.; Bowers P. Chapter Two: CAD and Design Automation in Piping
Design. Book chapter. The Planning Guide to Piping Design (Second Edition),
2018, pp. 39-59.
[7]. Bieker U.; Neumann A. Using logic programming and coroutining for electronic
CAD. Open archive, Research article. The Journal of Logic Programming,
Volume 26, Issue 2, February 1996, pp.199-215.
[8]. Sun X.; Robaldo L. On the complexity of input/output logic. Research article.
Journal of Applied Logic, Volume 25, December 2017, pp. 69-88.
168
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
[13]. Woodward M.R. The use of Nassi-Shneiderman charts and supporting tools
in software engineering education. Research article. Computers & Education,
Volume 11, Issue 4, 1987, pp. 267-279.
[14]. Barišić A; Amaral, V.; Goulão M. Usability driven DSL development with USE-
ME. Research article. Computer Languages, Systems & Structures, Volume
51, January 2018, pp. 118-157.
[15]. Ishikawa M.; Maruta I. Rapid Prototyping for Control Education using Arduino
and Open-Source Technologies. Research article. IFAC Proceedings
Volumes, Volume 42, Issue 24, 2010, pp. 317-321.
[16]. Townsend L.W.; PourArsalan M.; Hall M.I.; Anderson J.A.; Adamczyk A.M.
Estimates of Carrington-class solar particle event radiation exposures on Mars.
Research article. Acta Astronautica, Volume 69, Issues 7–8, September–
October 2011, pp. 397-405.
[18]. Kater A.; Wriedt J.P.; Meurer T. Finite Element Modeling and Motion Planning
of an Adaptive Elastic Wingsail. Research article. IFAC-PapersOnLine, Volume
50, Issue 1, July 2017, pp. 6768-6773.
[19]. Singiresu S.R. Chapter 22: Finite Element Analysis Using ANSYS. Book
chapter. The Finite Element Method in Engineering (Sixth Edition), 2018, pp.
703-721.
[21]. Liu D.; Wang W.; Liu J. Sensitivity Analysis of Meteorological Parameters on
Building Energy Consumption. Open access, Research article. Energy
Procedia, Volume 132, October 2017, pp. 634-639.
[23]. Habte F.; Chowdhury A.G.; Zisis I. Effect of wind-induced internal pressure
on local frame forces of low-rise buildings. Research article. Engineering
Structures, Volume 143, 15 July 2017, pp. 455-468.
[24]. Dassi F.; Kamenski L.; Farrell P.; Si H. Tetrahedral mesh improvement using
moving mesh smoothing, lazy searching flips, and RBF surface reconstruction.
Research article. Computer-Aided Design, In press, corrected proof, Available
online 8 December 2017.
169
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
[25]. Khanna K.; V.K. Gupta; Nigam S.P. Creep Analysis in Functionally Graded
Rotating Disc Using Tresca Criterion and Comparison with Von-Mises
Criterion. Research article. Materials Today: Proceedings, Volume 4, Issue 2,
Part A, 2017, pp. 2431-2438.
[29]. Neumann A.; Witzke, A.; Jones, S.A.; Schmitt, G. Representative Terrestrial
Solar Brightness Profiles, Journal of Solar Energy Engineering, Vol. 124, May
2002, pp. 198-204.
[30]. Steele C.R.; Balch, C.D.; Jorgensen, G.J.; Wendelin, T.; Lewandowski, A.
Membrane Dish Analysis: A Summary of Structural and Optical Analysis
Capabilities, National Renewable Energy Laboratory, Golden, CO. 1991.
NREL/TP-253-3432.
[31]. Ratzel A.C.; Boughton, B.D. CIRCE.001: A Computer Code for Analysis of
Point-Focus Concentrators with Flat Targets, SAND86-1866, 1987, Sandia
National Laboratories, Albuquerque, NM.
170
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Figura 4.1. Fases del diseño mecatrónico (fuente: National Instruments) .................................... 105
Figura 4.2. Planta del prototipo mecatrónico (fuente: autor).......................................................... 107
Figura 4.3. Alzado del prototipo mecatrónico (fuente: autor) ......................................................... 107
Figura 4.4. Vista latera del prototipo mecatrónico (fuente: autor) .................................................. 108
Figura 4.5. Modelización de tareas: concentración y seguimiento solar (fuente: autor) ............... 109
Figura 4.6. Variables meteorológicas capturables por el prototipo (fuente: autor) ........................ 109
Figura 4.7. Diseño 3D completo del prototipo de seguidor termosolar (fuente: autor) .................. 110
Figura 4.8. Interfaz de trabajo de FreeCAD 0.16 (fuente: autor) ................................................... 111
Figura 4.9. Diseño 3D de varios dispositivos mecánicos (fuente: autor) ....................................... 112
Figura 4.10. Diseño 3D de dispositivos ópticos (fuente: autor) ..................................................... 113
Figura 4.11. Diseño 3D de varios dispositivos electrónicos (fuente: autor) ................................... 113
Figura 4.12. Interfaz de trabajo de Fritzing 0.9.3 (fuente: autor) ................................................... 114
Figura 4.13. Esquema Fritzing conceptual del prototipo (fuente: autor) ........................................ 115
Figura 4.14. Esquema Fritzing electrónico del prototipo (fuente: autor) ........................................ 116
Figura 4.15. Contabilización de entradas y salidas en el autómata (fuente: autor) ....................... 116
Figura 4.16. Interfaz de trabajo de PSeInt (fuente: autor) ............................................................ 119
Figura 4.17. Flujograma completo (fuente: autor) .......................................................................... 120
Figura 4.18. Texto auxiliar de ayuda para interpretación de variables (fuente: autor) .................. 121
Figura 4.19. Definición e inicialización de variables (fuente: autor) ............................................... 121
Figura 4.20. Posicionamiento manual con joystick (fuente: autor) ................................................ 121
Figura 4.21. Seguimiento en elevación (fuente: autor) .................................................................. 122
Figura 4.22. Seguimiento en azimut (fuente: autor) ....................................................................... 122
Figura 4.23. Comprobación velocidad de viento (fuente: autor) .................................................... 122
Figura 4.24. Comprobación caída de lluvia (fuente: autor) ............................................................ 123
Figura 4.25. Comprobación caída de granizo (fuente: autor) ........................................................ 123
Figura 4.26. Comprobación máximos de radiación (fuente: autor) ................................................ 123
Figura 4.27. Comprobación mínimos de radiación (fuente: autor) ................................................ 123
Figura 4.28. Comprobación alta temperatura colector paraboloidal (fuente: autor) ...................... 124
Figura 4.29. Comprobación alta temperatura colector cilindro-parabólico (fuente: autor) ............ 124
Figura 4.30. Captura de datos meteorológicos (fuente: autor) ...................................................... 124
Figura 4.31. Muestra de datos y condición de salida (fuente: autor) ............................................. 124
Figura 4.32. Interfaz de programación de Visualino (fuente: autor) .............................................. 125
Figura 4.33. Definición e inicialización de variables (fuente: autor) ............................................... 127
Figura 4.34. Movimiento manual en azimut y elevación con joystick (fuente: autor) .................... 127
Figura 4.35. Seguimiento en elevación (fuente: autor) .................................................................. 128
Figura 4.36. Seguimiento en azimut (fuente: autor) ....................................................................... 128
Figura 4.37. Comprobación velocidad de viento (fuente: autor) .................................................... 129
Figura 4.38. Comprobación caída de lluvia (fuente: autor) ............................................................ 129
Figura 4.39. Comprobación caída de granizo (fuente: autor) ........................................................ 129
Figura 4.40. Comprobación máximos de radiación (fuente: autor) ................................................ 130
Figura 4.41. Comprobación mínimos de radiación (fuente: autor) ................................................ 130
Figura 4.42. Comprobación alta temperatura colector paraboloidal (fuente: autor) ...................... 130
Figura 4.43. Comprobación alta temperatura colector cilindro-parabólico (fuente: autor) ............ 130
Figura 4.44. Captura de datos meteorológicos (fuente: autor) ...................................................... 131
Figura 4.45. Muestra global de datos operativos y meteorológicos (fuente: autor) ....................... 131
Figura 4.46. Metodología de elementos finitos (fuente: Julián García de León) .......................... 140
Figura 4.47. Valores básicos de velocidad de viento en España (fuente: CTE) ........................... 141
Figura 4.48. Dimensiones y fuerzas en colectores (fuente: Julián García de León) ..................... 144
Figura 4.49. Reacciones en tubo cilindro-parabólico (fuente: Julián García de León) .................. 144
Figura 4.50. Reacciones en soporte (fuente: Julián García de León) ........................................... 145
Figura 4.51. Esfuerzos sobre cruceta (fuente: Julián García de León) ......................................... 145
Figura 4.52. Esfuerzos sobre casquillos de cruceta (fuente: Julián García de León) ................... 146
Figura 4.53. Esfuerzos sobre polea (fuente: Julián García de León) ............................................ 147
Figura 4.54. Etapas de pre-procesado de diseño CAD (fuente: Julián García de León) .............. 148
Figura 4.55. Mallado y fuerzas en abrazadera (fuente: Julián García de León)............................ 149
Figura 4.56. Mallado y fuerzas en soporte (fuente: Julián García de León) .................................. 149
171
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Figura 4.57. Mallado y fuerzas en casquillo (fuente: Julián García de León) ................................ 150
Figura 4.58. Mallado y fuerzas en perfil en L (fuente: Julián García de León) .............................. 150
Figura 4.59. Mallado y fuerzas en polea (fuente: Julián García de León) ..................................... 150
Figura 4.60. Tensiones von Mises de abrazadera (fuente: Julián Garcia de León) ...................... 151
Figura 4.61. Tensiones von Mises de soporte (fuente: Julián Garcia de León) ............................ 151
Figura 4.62. Tensiones von Mises de casquillo (fuente: Julián García de León) .......................... 152
Figura 4.63. Tensiones von Mises de perfil de aluminio (fuente: Julián García de León) ............. 152
Figura 4.64. Tensiones von Mises en polea (fuente: autor)........................................................... 153
Figura 4.65. Haz de rayos por método Montecarlo (fuente: Comsol, Soltrace, Tonatiuh) ............ 154
Figura 4.66. Parametrización del factor de forma solar (fuente: autor) ......................................... 157
Figura 4.67. Parametrización óptica del concentrador cilindro-parabólico (fuente: autor) ............ 159
Figura 4.68. Parametrización óptica del absorbedor cilindro-parabólico (fuente: autor) ............... 160
Figura 4.69. Parametrización parte interior tubo cilindro-parabólico (fuente: autor) ..................... 160
Figura 4.70. Parametrización parte exterior tubo cilindro-parabólico (fuente: autor) .................... 161
Figura 4.71. Parametrización dimensional colector cilindro-parabólico (fuente: autor) ................. 161
Figura 4.72. Parametrización óptica concentrador paraboloidal (fuente: autor) ............................ 161
Figura 4.73. Parametrización dimensional concentrador paraboloidal (fuente: autor) .................. 162
Figura 4.74. Parametrización dimensional absorbedor paraboloidal (fuente: autor) ..................... 162
Figura 4.75. Parametrización trazado rayos colector cilindro-parabólico (fuente: autor) .............. 163
Figura 4.76. Trazado de rayos sobre colector cilindro-parabólico (fuente: autor) ......................... 163
Figura 4.77. Parametrización trazado rayos sobre colector paraboloidal (fuente: autor) .............. 164
Figura 4.78. Trazado de rayos sobre colector paraboloidal (fuente: autor) ................................... 164
Figura 4.79. Irradiancia del receptor cilindro-parabólico en trazado contorno (fuente: autor) ....... 165
Figura 4.80. Irradiancia del receptor cilindro-parabólico en trazado superficie (fuente: autor) ..... 165
Figura 4.81. Irradiancia del receptor paraboloidal trazado contorno (fuente: autor) ..................... 166
Figura 4.82. Irradiancia del receptor paraboloidal trazado superficie (fuente: autor) .................... 166
172
Jorge Pablo Díaz Velilla Diseño, modelado y simulación del prototipo
Tabla 1. Mapeo de señales del prototipo termosolar (fuente: autor) ............................................. 117
Tabla 2. Coeficiente de exposición ce según aspereza y altura (fuente: CTE) .............................. 142
Tabla 3. Coeficientes F, k, Z y L según tipo de entorno (fuente: CTE) .......................................... 142
Tabla 4. Coeficiente de presión cp o succión cs (fuente: CTE) ...................................................... 142
Tabla 5. Cálculos de factor F y coeficientes de exposición (fuente: Julián García de León) ........ 143
173
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
ÍNDICE CAPÍTULO 5
174
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
175
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
Figura 5.3. Engranajes, eje, polea y correa de transmisión (fuente: Model Reyna)
176
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
177
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
178
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
179
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
180
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
181
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
• Botón de encendido/apagado
182
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
Figura 5.11. Cotas dimensionales del Arduino Mega 2560 (fuente: www.arduino.cc)
• Entradas analógicas: 16
• Memoria SRAM: 8 KB
• Memoria EEPROM: 4 KB
183
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
184
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
185
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
186
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
Figura 5.19. Velocidad lineal y amperaje consumido por actuador (fuente: Firgelli)
187
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
188
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
5.3.3 Sensórica
189
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
190
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
191
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
192
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
193
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
donde:
a: altitud, en m
p: presión medida por el sensor, en hPa
p0: presión a nivel del mar, igual a 1.013,25 hPa
194
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
Figura 5.33. Esquema electrónico del sensor presión barométrica (fuente: Seeedstudio)
195
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
Figura 5.36. Esquema electrónico del sensor de calidad del aire (fuente: Seeedstudio)
197
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
198
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
Figura 5.40. Respuesta del sensor frente longitud de onda (fuente: Seeedstudio)
199
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
200
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
1.000
𝑖= (2)
𝑑
𝑖 ∙ 10
𝑣 = (3)
4
Sustituyendo la ecuación ( 2 ) en la ecuación ( 3 ) se tendrá:
1.000
∙ 10 2.500
𝑣= 𝑑 = (4)
4 𝑑
201
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
donde:
i: impulsos, adimensional
d: duración, en s
v: velocidad, en km/h
202
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
203
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
204
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
205
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
206
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
207
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
208
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
209
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
210
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
211
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
212
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
213
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
x10
x10
x10
x10
x10
Figura 5.63. Comparativa costes dispositivos prototipo vs industrial (fuente: autor)
214
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
215
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
216
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
217
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
A priori podría parecer que acometer las tareas de fabricación del prototipo
de seguidor termosolar no debería representar especial esfuerzo: nada más lejos
de la realidad, tal y como la experiencia empírica ha demostrado en este trabajo.
Son múltiples y variadas las herramientas, equipos y dispositivos que han sido
precisos para construirlo, toda vez que es un equipo construido a medida para este
proyecto, sin contar por tanto con manuales o videotutoriales de montaje (Klein T.P.
y Reinhart G., 2016). No en vano, en algunos casos ha sido preciso acudir al taller
de Fabricación Mecánica del centro de formación Salesianos Carabanchel, centro
en el cual el autor desarrolla su tarea profesional como profesor del Ciclo formativo
de grado superior denominado Eficiencia energética y energía solar térmica.
218
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
219
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
220
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
221
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
222
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
224
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
225
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
226
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
227
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
BIBILIOGRAFÍA CAPÍTULO 5
[2]. Liu D.; Wang G.; Ren T. Transmission principle and geometrical model of
eccentric face gear. Mechanism and Machine Theory, Volume 109, March
2017, pp 51-64.
[3]. Sánchez M.; Martínez I.G.; Rincón E.A.; Duran M.D.; Design and Thermal-optic
Analysis of an Ultra-solar Concentrator. Energy Procedia, Volume 57, 2014, pp
311-320.
[4]. Jin J.; Ling Y.; Hao Y. Similarity analysis of parabolic-trough solar collectors.
Applied Energy, Volume 204, 15 October 2017, pp 958-965.
[5]. Candelas F.A.; García G.J.; Puente S.; Pomares J.; Torres F. Experiences on
using Arduino for laboratory experiments of Automatic Control and Robotics.
IFAC-PapersOnLine, Volume 48, Issue 29, 2015, pp 105-110.
[6]. Frenzel L.E. Chapter 13: Software, Programming and Electronics: Today, Most
Electronics Is a Combination of Hardware and Software. Electronics Explained
(Second Edition), 2017, pp 325-346.
[11]. Andreev Y.S.; Demkovich N.A.; Isaev R.M.; Tselischev A.A.; Vasilkov S.D.
Effect of the Method for Treatment of the Functional Surface of a Piezoelectric
Vibration Sensor on its Sensitivity. Procedia Engineering, Volume 176, 2017,
pp 96-106.
228
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
[12]. Bord I.; Tardy P.; Menil F. Temperature influence on a differential capacitive
rain sensor performances. Sensors and Actuators B: Chemical, Volume 125,
Issue 1, 16 July 2007, pp 262-267.
[15]. Bognár G.; Szabó P.G.; Takács G. Generalization of the thermal model of
infrared radiation sensors. Microelectronics Journal, Volume 46, Issue 6, June
2015, pp 543-550.
[16]. Yina L.; Liang Z.; Chunliu Z.; Dongning W. High-Temperature sensor based
on peanut flat-end reflection structure. Optik - International Journal for Light
and Electron Optics, Volume 148, November 2017, pp 293-299.
[17]. Ross M.; Pister K.S.J. Micro-windmill for optical scanning and flow
measurement. Sensors and Actuators A: Physical, Volume 47, Issues 1–3,
March–April 1995, pp 576-579.
[18]. Dongsheng Y.; Hertle R.W. Digital head posture analysis using an electronic
compass. Journal of American Association for Pediatric Ophthalmology and
Strabismus, Volume 11, Issue 1, February 2007, pp 105.
[19]. Yinsheng W.; Shudong W.; Hongcai Z.; Hairong G.; Xueyong W. A high
resolution tilt measurement system based on multi-accelerometers.
Measurement, Volume 109, October 2017, pp 215-222.
[20]. Dicianno B.E.; Mahajan H.; Cooper R.A. Advanced Joystick Algorithms for
Computer Access Tasks. PM&R, Volume 7, Issue 6, June 2015, pp 555-561.
[21]. Leinonen M.; Juuti J.; Jantunen H. Interface circuit for resistive sensors
utilizing digital potentiometers. Sensors and Actuators A: Physical, Volume
138, Issue 1, 20 July 2007, pp 97-104.
[22]. Qureshi F.Z.; Terzopoulos D. Chapter 21: Virtual Vision for Camera Networks
Research. Academic Press Library in Signal Processing, Volume 4, 2014, pp
609-625.
[23]. Liang Ch.; Qile Z.; Zhigang H.; Xinyuan J.; Chuang S. GNSS Global Real-
time Augmentation Positioning: Real-time Precise Satellite Clock Estimation,
Prototype System Construction and Performance Analysis. Advances in Space
Research, In press, accepted manuscript, Available online 9 September 2017.
229
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
[26]. Spicer R.A.; Valdes P.J.; Spicer T.E.V.; Craggs H.J.; Yang J. New
developments in CLAMP: Calibration using global gridded meteorological data.
Palaeogeography, Palaeoclimatology, Palaeoecology, Volume 283, Issues 1–
2, 1 December 2009, pp 91-98.
[27]. Yonghui Z.; Hongxu L.; Qian S.; Kai W.; Guoliang C. Real time remote
monitoring and pre-warning system for Highway landslide in mountain area.
Journal of Hydrology, Volume 245, Issues 1–4, 1 May 2001, pp 32-47.
[29]. Olivieri E.; Rotter M.; De Combarieu M.; Forget P.; Pari P. Full range resistive
thermometers. Cryogenics, Volume 72, Part 2, December 2015, pp 148-152.
[30]. Polo J.; Fernandez-Neira W.G.; Alonso-García M.C. On the use of reference
modules as irradiance sensor for monitoring and modelling rooftop PV systems.
Renewable Energy, Volume 106, June 2017, pp 186-191.
[31]. Margaret R. W. Weather Vanes. Mayo Clinic Proceedings, Volume 90, Issue
1, January 2015, pp e7.
[32]. Devaraju J.T.; Suhas K.R.; Mohana H.K.; Vijaykumar A.P. Wireless Portable
Microcontroller based Weather Monitoring Station. Measurement, Volume 76,
December 2015, pp 189-200.
230
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
231
Jorge Pablo Díaz Velilla Construcción de prototipo y estación de validación
232
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
ÍNDICE CAPÍTULO 6
233
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
234
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
235
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
236
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.6. Trayectoria solar en polares del 21jun17 y 22sep17 (fuente: SunEarthTools)
Figura 6.7. Trayectoria solar en polares del 21dic17 y 20mar18 (fuente: SunEarthTools)
Figura 6.8. Trayectoria solar en cartesianas del 21jun17 y 22sep17 (fuente: SunEarthTools)
Figura 6.9. Trayectoria solar en cartesianas del 21dic17 y 20mar18 (fuente: SunEarthTools)
237
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
238
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
239
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.10. Fecha e instantes temporales para capturas de datos (fuente: autor)
240
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.12. Movimiento del seguidor hacia direcciones del espacio (fuente: autor)
241
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.16. Temperatura ambiente y humedad relativa del aire (fuente: autor)
• Presión barométrica, en hPa (se recuerda 1 hPa = 10-3 bar ≈ 10-3 atm).
242
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Las únicas variables que no se han representado de forma gráfica, son las
denominadas C (sensor de lluvia) e I (sensor de granizo), por motivos
constructivos obvios: ni los sensores, ni el prototipo al completo, están realmente
preparados para recibir agua de lluvia y/o el impacto del granizo. Cualquiera de los
antedichos efectos meteorológicos, pondría en peligro la integridad física de la
electrónica implicada, malogrando todo el conjunto de forma inexorable. Es por ello
que se tuvo especial cuidado en la toma de datos de los días señalados, para evitar
al máximo la posible caída de agua y/o granizo sobre el seguidor: afortunadamente,
la meteorología acompañó y no se presentaron ninguna de las antedichas
circunstancias. Los valores digitales que ambas variables adquieren en todas las
tomas de datos, son siempre 0 (no hay presencia ni de agua ni de granizo), con lo
cual carece de sentido su representación gráfica por no aportar nada al respecto.
243
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
244
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
245
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
246
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
247
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
248
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
249
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
250
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
251
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
252
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.32. Temperaturas focos (Tp vs Tcp) del 21jun17 (fuente: autor)
Figura 6.33. Temperaturas focos (Tp vs Tcp) del 22sep17 (fuente: autor)
253
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.34. Temperaturas focos (Tp vs Tcp) del 21dic17 (fuente: autor)
Figura 6.35. Temperaturas focos (Tp vs Tcp) del 20mar18 (fuente: autor)
254
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
255
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
En este primer epígrafe se trazan las irradiancias (en W/m2) capturadas por
los sensores SO y SE, para cada día (Hinkelman L.M. y Schaeffer N., 2017).
.
Figura 6.36. Irradiancias sensores (SO vs SE) del 21jun17 (fuente: autor)
Figura 6.37. Irradiancias sensores (SO vs SE) del 22sep17 (fuente: autor)
256
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.38. Irradiancias sensores (SO vs SE) del 21dic17 (fuente: autor)
Figura 6.39. Irradiancias sensores (SO vs SE) del 20mar18 (fuente: autor)
257
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
258
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.40. Irradiancias sensores (SN vs SS) del 21jun17 (fuente: autor)
Figura 6.41. Irradiancias sensores (SN vs SS) del 22sep17 (fuente: autor)
259
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.42. Irradiancias sensores (SN vs SS) del 21dic17 (fuente: autor)
Figura 6.43. Irradiancias sensores (SN vs SS) del 20mar18 (fuente: autor)
260
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
261
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
262
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
263
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
264
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
265
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Sea como fuere, en todos los casos analizados, dicha velocidad de viento
oscila entre valores comprendidos entre 0 y 6 m/s (21,6 km/h), cifras muy
contenidas incluso para encontrarnos a una altura de unos 10 m sobre el nivel del
suelo en el taller de energías renovables del centro de formación Salesianos
Carabanchel. Especialmente bajas parecen ser las velocidades de viento el 21dic17
en contraposición al 20mar18 que es el día de mayores velocidades elevadas.
266
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
267
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
268
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
269
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Ahora bien, el incremento de temperatura (ΔT) para cada uno de los días
señalados, parece muy similar en todos los casos: oscilando entre 7-10ºC
aproximadamente, lo cual no deja de ser un dato “extraño” frente a lo que cabría
esperarse (mayor incremento en verano que en invierno, por ejemplo).
Volviendo una vez más a cotejar los datos capturados por el prototipo frente
a datos meteorológicos contrastados, las mediciones del prototipo parecen más
“suaves” (al menos visualmente) que las propias de dichos datos estadísticos.
Cierto es que es preciso extraer de dichos datos el intervalo coincidente con las
horas de trabajo del prototipo (se recuerda, de 7:30-19:30h). Y que además es
necesario escalar el factor de forma de ambos gráficos para que su aspecto sea
coincidente. Además, se sigue manteniendo el ΔTmáx ≈ 10ºC. Luego haciendo una
reflexión sobre ambas circunstancias se vuelve a dar validez a los resultados
obtenidos con el prototipo mecatrónico.
270
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
271
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
La humedad relativa de las fechas indicadas, en tanto por ciento (%), son las
indicadas en los gráficos adjuntos.
272
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
273
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
274
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
275
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
276
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
277
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
La calidad del aire en grandes ciudades es uno de los temas que cobra cada
vez mayor relevancia. La creciente conciencia medioambiental de la población, la
proliferación de enfermedades respiratorias e incluso las muertes prematuras de
ciudadanos por causas de difícil identificación (30.000 personas/año) pero
imputables a la mala calidad del aire, justifican con creces la inclusión de sensores
capacitados para identificar dicha calidad (Van Wart J. et al., 2015). Los valores
mostrados indican una adecuada calidad general del aire (buena, nivel 4) que solo
se ve minorada (aceptable, nivel 3), en determinadas franjas horarias cuya pauta
parece repetirse. Y dicha pauta no es otra que las horas de entrada y salida de
alumnos al centro, especialmente de la sección Enseñanza Secundaria
Obligatoria-Bachiller-Formación Profesional (8:00h) y Primaria (9:00h) o (14:15h
para el primer bloque de los antedichos alumnos y 16:30h para el segundo).
Especialmente importante para el caso de los alumnos que en función de su
temprana edad, son traídos por sus progenitores en transporte privado con motores
de combustión interna alternativos (ver imágenes siguientes). Se puede deducir que
esas bajadas de la calidad del aire se pueden deber a los tubos de escape de los
motores de combustión de dichos automóviles.
Figura 6.76. Hora punta de entrada de automóviles en Salesianos Carabanchel (fuente: autor)
278
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
279
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
280
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
281
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Sí que parece cierto, que los días 21dic17 y 20mar18 son días de mayor
presión barométrica (alcanzando máximos cercanos a los 948 hPa) frente a las
fechas de 21jun17 y 22sep17. Y sea como fuere, las diferencias de presión intradía
tampoco son extremadamente altas, de escasamente 5 hPa el 20mar18.
282
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
283
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
284
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.90. Comparativa imagen real vs termográfica de colector paraboloidal (fuente: autor)
285
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
286
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
287
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
288
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
289
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
290
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
291
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
292
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
293
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
294
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
295
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
296
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
297
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.114. Posibles fases del despliegue de colector paraboloidal (fuente: Wolfram)
298
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
• Díaz J.P.; Morón, C.; Ferrández D.; Saiz P.; Vega G.; (2017). “Mechatronic
Prototype of Parabolic Solar Tracker”. Sensors, 16(6), 882;
doi:10.3390/s16060882.
• Morón, C.; Ferrández D.; Saiz P.; Vega G.; Díaz J.P. (2017). “New Prototype
of Photovoltaic Solar Tracker Based on Arduino”. Energies, 10(9), 1298;
doi:10.3390/en10091298.
• Díaz, J.P.; Ferrández, D.; Morón C. (2017). “Frío solar como alternativa a
climatización con refrigerantes: análisis de aplicabilidad a caso real mediante
simulación por ordenador”, CITE 2017, Madrid, España. ISBN: 978-84-
16397-56-3.
• Morón C.; Díaz, J.P.; Ferrández, D.; Morón C. (2017); Saiz P. “Electrical
energy generation in building by means of Stirling dish technology“, Anales
de Edificación. Vol. 3, Nº2, 9-16 (2017). ISSN: 2444-1309. Doi:
10.20868/ade.2017.3567.
• Morón C.; Díaz, J.P.; Ferrández, D.; Morón C. (2017); Saiz P. “Protocol for
technical evaluation of a Solar Thermal Installation “, Anales de Edificación.
Vol. 2, Nº2, 21-27 (2016). ISSN: 2444-1309. Doi: 10.20868/ade.2016.3307.
299
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
• Díaz, J.P.; Ferrández, D.; Morón C. (2018); Vega G. “El software libre en la
enseñanza: aplicación en los ciclos formativos de grado superior”, CINIE
2017, Madrid, España. ISBN: 978-84-16397-55-6.
• Ferrández, D.; Díaz, J.P.; Morón C.; Ramos Mª Paz. “Learning and
knowledge technologies in professional training”. INTED 2016.
• Ramos Mª Paz. Ferrández, D.; Saiz P.; Díaz, J.P.; Morón C.; “Mindfulness in
vocational education training”. Edulearn 2016.
300
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Herramientas hardware:
Herramientas software:
• Microsoft Windows 10
• FreeCAD 0.16
• Visualino 1.0
• Evernote 7.17.1
• PSeInt 2017
• Fritzing 9.3b
• SolTrace 2012.7.9
• SmartView 4.3
• AR Compass 3D
301
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
BIBLIOGRAFÍA CAPÍTULO 6
[1]. Farah S.; Saman W.; Boland J. Development of robust meteorological year
weather data. Research article. Renewable Energy, Volume 118, April 2018,
pp 343-350.
[2]. Qing X.; Niu Y. Hourly day-ahead solar irradiance prediction using weather
forecasts by LSTM. Research article. Energy, Volume 148, 1 April 2018, pp.
461-468.
[3]. Cheng W.Y.Y.; Liu Y.; Bourgeois A.; Wu Y.; Haupt S.E. Short-term wind
forecast of a data assimilation/weather forecasting system with wind turbine
anemometer measurement assimilation. Research article. Renewable Energy,
Volume 107, July 2017, pp. 340-351.
[4]. Hong T.; Chang W-K; Lin H-W. A fresh look at weather impact on peak
electricity demand and energy use of buildings using 30-year actual weather
data. Research article. Applied Energy, Volume 111, November 2013, pp. 333-
350.
[5]. Hinkelman L.M.; Schaeffer N. Relating solar resource and its variability to
weather and climate across the northwestern United States. Research article.
Solar Energy, Volume 157, 15 November 2017, pp. 966-978.
[6]. Soares E.; Costa P.; Costa B.; Leite D. Ensemble of evolving data clouds and
fuzzy models for weather time series prediction. Research article. Applied Soft
Computing, Volume 64, March 2018, pp. 445-453.
[7]. Tsoka S.; Tolika K.; Theodosiou T.;Tsikaloudaki K.; Bikas D. A method to
account for the urban microclimate on the creation of ‘typical weather year’
datasets for building energy simulation, using stochastically generated data.
Research article. Energy and Buildings, Volume 165, 15 April 2018, pp. 270-
283.
[8]. Burleyson C.D.; Voisin N.; Taylor Z.T.; Xie Y.; Kraucunas I. Simulated building
energy demand biases resulting from the use of representative weather
stations. Research article. Applied Energy, Volume 209, 1 January 2018, pp.
516-528.
[9]. Hosseini M.; Lee B.; Vakilinia S. Energy performance of cool roofs under the
impact of actual weather data. Research article. Energy and Buildings, Volume
145, 15 June 2017, pp. 284-292.
[11]. Anaya K.L.; Pollitt M.G. Using stochastic frontier analysis to measure the
impact of weather on the efficiency of electricity distribution businesses in
302
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
[12]. Hosseini M.; Tardy F.; Lee B. Cooling and heating energy performance of a
building with a variety of roof designs; the effects of future weather data in a
cold climate. Research article. Journal of Building Engineering, Volume 17,
May 2018, pp. 107-114.
[13]. Van Wart J.; Grassini P.; Yang H.; Claessens L.; Cassman K.G. Creating
long-term weather data from thin air for crop simulation modeling. Open
access, Research article. Agricultural and Forest Meteorology, Volumes 209–
210, 15 September 2015, pp. 49-58.
[14]. Tao S.; Corcoran J.; Hickman M.; Stimson R. The influence of weather on
local geographical patterns of bus usage. Research article. Journal of
Transport Geography, Volume 54, June 2016, pp. 66-80.
303
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
304
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
Figura 6.58. Temperatura ambiente del 21dic17 (fuente: autor) ................................................... 269
Figura 6.59. Temperatura ambiente del 20mar18 (fuente: autor) .................................................. 269
Figura 6.60. Temperatura ambiente estadística del 21jun17 (fuente: AEMET)............................. 270
Figura 6.61. Temperatura ambiente estadística del 22sep17 (fuente: AEMET)............................ 271
Figura 6.62. Temperatura ambiente estadística del 21dic17 (fuente: AEMET) ............................. 271
Figura 6.63. Temperatura ambiente estadística del 20mar18 (fuente: AEMET) ........................... 271
Figura 6.64. Humedad relativa del 21jun17 (fuente: autor) ........................................................... 272
Figura 6.65. Humedad relativa del 22sep17 (fuente: autor) .......................................................... 272
Figura 6.66. Humedad relativa del 21dic17 (fuente: autor) ........................................................... 273
Figura 6.67. Humedad relativa del 20mar18 (fuente: autor) .......................................................... 273
Figura 6.68. Humedad relativa estadística del 21jun17 (fuente: AEMET) ..................................... 274
Figura 6.69. Humedad relativa estadística del 22sep17 (fuente: AEMET) .................................... 275
Figura 6.70. Humedad relativa estadística del 21dic17 (fuente: AEMET) ..................................... 275
Figura 6.71. Humedad relativa estadística del 20mar18 (fuente: AEMET) ................................... 275
Figura 6.72. Calidad del aire 21jun17 (fuente: autor) .................................................................... 276
Figura 6.73. Calidad del aire 22sep17 (fuente: autor) ................................................................... 276
Figura 6.74. Calidad del aire 21dic17 (fuente: autor)..................................................................... 277
Figura 6.75. Calidad del aire 20mar18 (fuente: autor) ................................................................... 277
Figura 6.76. Hora punta de entrada de automóviles en Salesianos Carabanchel (fuente: autor). 278
Figura 6.77. Calidad aire estadística del 21jun17 (fuente: Munimadrid) ....................................... 278
Figura 6.78. Calidad aire estadística del 22sep17 (fuente: Munimadrid) ...................................... 279
Figura 6.79. Calidad aire estadística del 21dic17 (fuente: Munimadrid) ........................................ 279
Figura 6.80. Calidad aire estadística del 20mar18 (fuente: Munimadrid) ...................................... 279
Figura 6.81. Presión barométrica 21jun17 (fuente: autor) ............................................................. 280
Figura 6.82. Presión barométrica 22sep17 (fuente: autor) ............................................................ 280
Figura 6.83. Presión barométrica 21dic17 (fuente: autor) ............................................................. 281
Figura 6.84. Presión barométrica 20mar18 (fuente: autor) ............................................................ 281
Figura 6.85. Presión barométrica estadística del 21jun17 (fuente: AEMET) ................................. 282
Figura 6.86. Presión barométrica estadística del 22sep17 (fuente: AEMET) ................................ 283
Figura 6.87. Presión barométrica estadística del 21dic17 (fuente: AEMET) ................................. 283
Figura 6.88. Presión barométrica estadística del 20mar18 (fuente: AEMET) ............................... 283
Figura 6.89. Interfaz de trabajo de SmartView 4.3 (fuente: Autor) ................................................ 284
Figura 6.90. Comparativa imagen real vs termográfica de colector paraboloidal (fuente: autor) .. 285
Figura 6.91. Representación 3D de termografía de colector paraboloidal (fuente: autor) ............ 285
Figura 6.92. Histograma térmico de colector paraboloidal (fuente: autor) ..................................... 285
Figura 6.93. Imagen real vs termográfica de colector cilindro-parabólico (fuente: autor) ............. 286
Figura 6.94. Representación 3D de termografía de colector cilindro-parabólico (fuente: autor) ... 286
Figura 6.95. Histograma térmico de colector cilindro-parabólico (fuente: autor) ........................... 286
Figura 6.96. Alumnos de Salesianos Carabanchel con el prototipo (fuente: autor) ...................... 288
Figura 6.97. Adición de motor Stirling a prototipo termosolar (fuente: autor) ................................ 288
Figura 6.98. Alumnos de Robótica educativa en Salesianos Carabanchel (fuente: autor) ........... 289
Figura 6.99. Correa de tracción alternativa (fuente: Edison) ......................................................... 290
Figura 6.100. Identificación de sección real de cilindro-parábola (fuente: autor) .......................... 290
Figura 6.101. Arduino MEGA industrial Controllino (fuente: autor) ............................................... 292
Figura 6.102. Soluciones termográficas portátiles (fuente: Seekthermal, Flir y Thermapp).......... 292
Figura 6.103. Simulación de maniobras de programación (fuente: autor) ..................................... 293
Figura 6.104. Modelización de seguidor mediante funciones de transferencia (fuente: autor) ..... 293
Figura 6.105. Soluciones SCADA de gestión remota (fuente: MyOpenlab y Labview) ................. 294
Figura 6.106. Pantalla TFT de visualización de datos (fuente: autor) ........................................... 294
Figura 6.107. Programa gratuito de control difuso (fuente: Xfuzzy) .............................................. 295
Figura 6.108. Panel fotovoltaico con batería integrada (fuente: autor) .......................................... 295
Figura 6.109. Brújula en RA desde perspectiva del seguidor (fuente: autor) ................................ 296
Figura 6.110. Trayectoria en RA desde perspectiva del seguidor (fuente: autor) ......................... 296
Figura 6.111. Cálculo de engranajes con Relatran 3.5 (fuente: Relatran) .................................... 297
Figura 6.112. Seguidor fotovoltaico retráctil (fuente: Smartflower) ................................................ 297
Figura 6.113. Fases del despliegue de paneles fotovoltaicos (fuente: Smartflower) .................... 298
Figura 6.114. Posibles fases del despliegue de colector paraboloidal (fuente: Wolfram) ............. 298
305
Jorge Pablo Díaz Velilla Pruebas de funcionamiento y conclusiones
306