DISEÑO Y CONSTRUCCIÓN DE UN MOTOR STIRLING SOLAR TIPO ALFA DE

PEQUEÑO PORTE

BRAYAN LEANDRO ZAPATA GUZMÁN

CÓDIGO: 2131213

MARÍA FERNANDA BERNAL CASTAÑEDA

CÓDIGO: 2130243

UNIVERSIDAD AUTÓNOMA DE OCCIDENTE

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE ENERGÉTICA Y MECÁNICA
PROGRAMA INGENIERÍA MECÁNICA
SANTIAGO DE CALI
2018
DISEÑO Y CONSTRUCCIÓN DE UN MOTOR STIRLING SOLAR TIPO ALFA DE
PEQUEÑO PORTE

BRAYAN LEANDRO ZAPATA GUZMÁN

CÓDIGO: 2131213
MARÍA FERNANDA BERNAL CASTAÑEDA
CÓDIGO: 2130243

PROYECTO DE GRADO para optar al título de

Ingeniero Mecánico

Director
CARLOS EDUARDO CASTANG MONTIEL
Ingeniero Mecánico, Mg

UNIVERSIDAD AUTÓNOMA DE OCCIDENTE

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE ENERGÉTICA Y MECÁNICA
PROGRAMA INGENIERÍA MECÁNICA
SANTIAGO DE CALI
2018
 Nota de aceptación:

Aprobado por el Comité de Grado en

cumplimiento de los requisitos exigidos
por la Universidad Autónoma de
Occidente para optar al título de
Ingeniero Mecánico

___________________________
FELIX GONZALEZ PEREZ

___________________________
JUAN RICARDO VIDAL

Santiago de Cali, 05 de octubre de 2018

3
 AGRADECIMIENTOS

Primero que todo para poder haber llegado a estas instancias de nuestras vidas
debemos hacer un repaso de todos aquellos que hicieron parte de ellas, ya que a
partir de cada amigo y cada profesor pudimos aprender algo nuevo.

Darle las gracias a Dios por darnos la vida y permitirnos a pesar de las dificultades
culminar nuestro proceso.

Queremos agradecerle a nuestros padres y familiares que hicieron posible nuestro

desarrollo como profesionales con su apoyo, amor, dedicación y paciencia.

Les agradecemos a nuestros profesores por compartirnos sus conocimientos y

enseñanzas, las cuales nos permitieron tener unas buenas bases para culminar
nuestra carrera.

Le agradecemos a nuestro director de trabajo de grado, Carlos Castang, por confiar

en nosotros y ponernos este gran reto. Sus enseñanza, apoyo y contribución en la
finalización de este proyecto fueron importantes.

A nuestros amigos (los mecánicos, multimedia, ambientales, biomédicos,

mecatrónicos) queremos expresarles nuestro agradecimiento por su gran amistad,
pues fueron parte de ese gran momento en nuestras vidas que fue el pregrado,
esperamos seguir siempre incondicionales para sus vidas, éxitos.

4
 CONTENIDO
pág.

RESUMEN 16

INTRODUCCIÓN 18

1. PLANTEAMIENTO DEL PROBLEMA 20

2. JUSTIFICACIÓN 22

3. OBJETIVOS 23

3.1 OBJETIVO GENERAL 23

3.2 OBJETIVOS ESPECÍFICOS 23

4. ESTADO DEL ARTE 24

5. MARCO TEÓRICO 29

5.1 MOTOR STIRLING 29

5.1.1 Ventajas y desventajas del motor Stirling 29

5.1.2 Clasificación de los motores Stirling: 30

5.1.2.1 Motor Stirling tipo Alfa 30

5.1.3 Mecanismo del motor Stirling tipo alfa 31

5.1.4 Ciclo Stirling ideal: 35

5.1.5 Ciclo Stirling real: 37

5.1.6 Modelos termodinámicos 39

5.1.6.1 Modelos de Primer Orden 39

5.1.6.2 Modelos de Segundo y Tercer Orden 46

5
5.2 COLECTOR SOLAR 46

5.2.1 Sistema disco parabólico 46

6. METODOLOGÍA 54

7. DESARROLLO DE PARÁMETROS INICIALES DEL DISEÑO

DEL MOTOR STIRLING Y EL COLECTOR SOLAR 55

7.1 DISEÑO DEL MOTOR STIRLING TIPO ALFA 55

7.1.1 Cálculo para el volumen barrido: 55

7.1.2 Cálculo de diámetro de los pistones 56

7.1.3 Diseño de las camisas de los pistones 56

7.2 DISEÑO DEL COLECTOR SOLAR 60

8. DISEÑO DEL MOTOR STILING SOLAR 63

8.1 COMPONENTES DEL MOTOR STIRLING 63

8.2 COMPONENTES DEL COLECTOR SOLAR 79

9. CONSTRUCCIÓN DEL MOTOR STIRLING SOLAR 82

9.1 FABRICACIÓN DE COMPONENTES DEL MOTOR STIRLING 82

9.2 FABRICACIÓN DE COMPONENTES DEL DISCO

PARABÓLICO 87

9.3 ENSAMBLE DEL MOTOR STIRLING SOLAR 88

9.4 COSTOS DE ELABORACIÓN DEL MOTOR STIRLING

SOLAR: 89

10. RESULTADOS Y VALIDACIÓN CON EL MODELO

MATEMÁTICO 90

10.1 RESULTADOS 90

6
10.2 VALIDACIÒN CON EL MODELO MATEMÀTICO 95

11. CONCLUSIONES 107

12. RECOMENDACIONES 109

BIBLIOGRAFIA 110

ANEXOS 114

7
 LISTA DE CUADROS

pág.

Cuadro 1. Ciclo Stirling del motor tipo Alfa 33

8
 LISTA DE FIGURAS
pág.

Figura 1. Motores Stirling tipo Alpha-Beta-Gamma 30

Figura 2. Configuración alfa de motor Stirling 31

Figura 3. Etapas del ciclo Stirling y su variación de volumen 35

Figura 4. Ciclo ideal Stirling 37

Figura 5. Ciclo ideal vs Ciclo real 38

Figura 6. Esquema de la teoría de Schmidt 41

Figura 7. Sección transversal de un concentrador parabólico 47

Figura 8. Dimensión de la imagen generada en el punto focal de un

colector de concentración. 52

Figura 9. Parámetros de las dimensiones del motor Stirling solar tipo

Alfa 58

Figura 10. Dimensiones del cilindro caliente (cotas en mm) 59

Figura 11. Dimensiones del cilindro frio (cotas en mm) 59

Figura 12. Relación de concentración dependiendo el ángulo de

borde (ϕr). 61

Figura 13. Diseño cigüeñal en CAD 63

Figura 14. Dimensiones del cigüeñal (cotas en mm) 64

Figura 15. Diseño de las bielas en CAD 65

Figura 16. Secciones de las bielas 65

Figura 17. Dimensiones de las bielas (cotas en mm) 66

Figura 18. Diseño de pistones en CAD 67

Figura 19. Dimensiones de los pistones (cotas en mm) 67

Figura 20. Diseño del cilindro de la zona caliente en CAD 69

9
Figura 21. Dimensiones del cilindro de la zona caliente (cotas en
mm) 69

Figura 22. Diseño del cilindro de la zona fría en CAD 71

Figura 23. Dimensiones del cilindro de la zona fría 72

Figura 24. Perfil de una aleta (cotas en metros) 73

Figura 25. Diseño del regenerador en CAD 76

Figura 26. Diseño de la estructura del Motor Stirling Solar en CAD 77

Figura 27. Dimensiones de la estructura del Motor Stirling Solar 77

Figura 28. Diseño del volante de inercia en CAD 78

Figura 29. Diámetro del volante de inercia 79

Figura 30. Croquis del diseño del colector solar con parámetros
calculados 80

Figura 31. Diseño del colector solar en CAD 80

Figura 32. Cigüeñal del mecanismo 82

Figura 33. Rodamientos seleccionados para el mecanismo 83

Figura 34. Bielas del mecanismo 83

Figura 35. Pistones del mecanismo 84

Figura 36. Cilindros de la zona fría y caliente del mecanismo 84

Figura 37. Buje de latón para el cilindro de la zona fría 85

Figura 38. Regenerador del mecanismo 85

Figura 39. Caja del mecanismo 86

Figura 40. Plato solar del motor Stirling 87

Figura 41. Motor Stirling Solar 88

Figura 42. Ensamble del motor Stirling con el plato solar 88

10
Figura 43. Prueba en seco del motor Stirling con motor trifásico 92

Figura 44. Prueba de temperatura de los cilindros con el plato solar 93

Figura 45. Prueba del motor Stirling con soplete 94

Figura 46. Presión vs Volumen del motor Stirling solar tipo alfa 102

Figura 47. Variación del volumen de expansión vs el ángulo del

cigüeñal 103

Figura 48. Variación del volumen de expansión vs el ángulo del

cigüeñal 103

Figura 49. Variación de la presión vs el ángulo del cigüeñal 104

Figura 50. Variación de potencia con respecto al diámetro de

diseño 104

11
 LISTA DE TABLAS

pág.

Tabla 1. Estado del arte 24

Tabla 2. Ventajas y desventajas del motor Stirling 29

Tabla 3. Errores típicos para concentradores parabólicos. 50

Tabla 4. Simbología de ecuación para la eficiencia óptica 53

Tabla 5. Simbología de ecuación para la energía absorbida. 53

Tabla 6. Simbología de las dimensiones de los cilindros de las zonas

de expansión y compresión 57

Tabla 7. Parámetros de diseño del colector según los diámetros de

apertura (Da) 62

Tabla 8. Valores para determinar el calor del receptor (Qrecep). 70

Tabla 9. Costos de fabricación 89

Tabla 10. Valores de temperaturas registrados en cada prueba 90

Tabla 11. Datos obtenidos de la prueba en seco 92

Tabla 12. Valores de temperaturas en las zonas de expansión y

compresión con diferentes fuentes de calor 94

Tabla 13. Parámetros de operación iniciales para el desarrollo del

modelo termodinámico del motor Stirling solar 96

Tabla 14. Magnitud de volúmenes con sus respectivas unidades 97

Tabla 15. Valores de Presión y Volumen en función del ángulo del

cigüeñal 101

Tabla 16. Delta de temperatura vs eficiencia del ciclo Stirling 105

Tabla 17. Calor captado por el receptor teórica vs experimental 106

Tabla 18. Potencia alcanzada teóricamente por medio de Beale y

Schmidt 106
12
Tabla 19. Eficiencia teórica con las temperaturas alcanzadas 106

13
 LISTA DE DIAGRAMAS

Diagrama 1. Metodología 54

Diagrama 2. Diseño del motor 55

14
 LISTA DE ANEXOS
pág.

Anexo A. Tabla de emisividad normal de los metales 114

Anexo B. Tabla de absortividad en superficies de los metales 115

Anexo C. Tabla de propiedades térmicas de algunos elementos

metálicos 116

Anexo D. Eficiencia de aletas circulares de espesor constante 117

Anexo E. Planos 118

15
 RESUMEN

En este estudio se realiza el diseño y construcción de un motor Stirling solar tipo

alfa empleando como fuente de alimentación un disco parabólico. Este prototipo
utiliza aire como fluido de trabajo a 1 atm. El desarrollo del proyecto se dividió en
tres fases. Inicialmente se realizó el modelamiento termodinámico del motor a través
de métodos de primer orden utilizando Matlab con el fin de analizar teóricamente la
eficiencia y potencia entregada. Seguido a ello se realizó el estudio de movimiento
en Solidworks para detallar en CAD el diseño del motor y posteriormente seleccionar
los materiales y procesos de manufactura adecuados. Finalmente, se construyó el
motor Stirling y el colector solar para la validación experimental de la eficiencia y
potencia alcanzada. El proceso de diseño permitió establecer que el volumen óptimo
de desplazamiento es de 78 cm3 con un diámetro de pistón de 5 cm, alcanzando un
gradiente de temperatura entre la zona de compresión y expansión de 250°C. Se
realizaron pruebas con tres tipos de fuente de calor, con el fin de medir las
temperaturas alcanzadas en las zonas de expansión y compresión. Se utilizaron
como fuentes de calor, el plato solar, un mechero y un soplete. El proyecto forma
parte del semillero de investigación de sistemas térmicos de la U.A.O. el cual
pretende aportar en la solución de la problemática de las ZNI (Zonas No
Interconectadas) del país haciendo uso de fuentes de energía renovables por medio
de sistemas alternativos de generación de energía eléctrica.

Palabras clave: energía, eficiencia, disco parabólico, motor Stirling, métodos de

primer orden, potencia neta, tipo alfa, ZNI.

16
 ABSTRACT

This study is about the design and construction of a solar Stirling engine type alpha
using as a power source a parabolic disc. This prototype uses air as a working fluid
at 1 atm. The development of the project has been carried out in three phases.
Initially, the thermodynamic modeling of the motor was carried out by means of first
order methods using Matlab in order to analyze theoretically the efficiency and
delivered power. Then it is done the movement study in Solidworks to detail in CAD
the design of the engine and then select the appropriate materials and manufacturing
processes. Finally, the Stirling engine and the solar collector were built for the
experimental validation of the efficiency and power reached. The design process
was allow to establish the optimum displacement volume that is 78 cm3 with a piston
diameter of 5 cm reaching a temperature gradient between the area of compression
and expansion of 250°C. Three types of heat source were tested in order to measure
the temperatures reached in the expansion and compression zones. They were used
as heat sources, the solar dish, a lighter and a blowtorch. The project is part of a
group of investigation of thermal systems of the U.A.O which tries to contribute in
the solution of the problem of the NIZ (Non-interconnected zones) of the country
making use of renewable energy sources by means of alternative systems of
generation of electric power.

Keywords: energy, efficiency, parabolic disc, Stirling engine, first order methods,
net power, type alpha, NIZ.

17
 INTRODUCCIÓN

En la actualidad Colombia cuenta con una matriz energética estable donde

predomina la producción de energía por medio de centrales hidroeléctricas y
centrales térmicas 1. Sin embargo, cerca de del 58% del territorio nacional pertenece
a las Zonas No Interconectas (ZNI), debido al costo que conlleva conectar estos
lugares a la red nacional ya que por lo general son apartados y de difícil acceso.
Según el Plan Indicativo de Expansión de Cobertura de Energía Eléctrica, realizado
por la Unidad de Planeación Minero-Energética (UPME) del Ministerio de Minas y
Energía 2, se considera que se requieren cerca de $4,3 billones de pesos
colombianos para poder suplir la necesidad de toda la población.

En el año 2017 el Instituto de Planificación y Promoción de Soluciones Energéticas

para las zonas no interconectadas (IPSE) informó que durante ese año se
dispusieron cerca de 100.000 millones de pesos en infraestructura 3 . De tal modo,
el escenario tanto para la población que habita estas zonas como para el país no
es muy optimo, ya que en un mundo donde el acceso a la tecnología es
determinante para el desarrollo socioeconómico de la población, el uso de la energía
es de vital importancia.

Las zonas no interconectadas se encuentran repartidas en 17 departamentos, cinco

capitales departamentales, 37 cabeceras municipales, 112 municipios y 1448
localidades 4. De la cuales cerca de 800 localidades (48% de las localidades) según
el Centro Nacional de Monitoreo (CNM) registran valores de radiación solar

1
UPME. Energía eléctrica SIN (Sistema Interconectado Nacional) [En línea]. Unidad de Planeación
Minero Energética. [Consultado: 20 de junio del 2018]. Disponible en Internet:
http://www1.upme.gov.co/InformacionCifras/Paginas/PETROLEO.aspx
2
SÁENZ, Jorge. En penumbras: 2,5 millones de colombianos no tienen energía eléctrica [En línea].
En: El Espectador. Agosto. 2017. [Consultado: 19 de abril del 2018]. Disponible en Internet:
https://www.elespectador.com/economia/en-penumbras-25-millones-de-colombianos-no-tienen-
energia-electrica-articulo-706892
3
Instituto de Planificación y Promoción de Soluciones Energéticas para las Zonas No
Interconectadas (IPSE). Informe de gestión[En línea] ipse 2017. p. 5. [Consultado: 20 de junio del
2018]. Disponible en Internet: http://www.ipse.gov.co/ipse/informes-de-gestion
4
UNIVERSIDAD NACIONAL DE COLOMBIA. Microrred generaría energía a bajo costo en zonas no
interconectadas [En línea]. Agenciadenoticias.unal.edu.co. 18 de abril de 2017. [Consultado: 20 de
junio del 2018]. Disponible en Internet: http://agenciadenoticias.unal.edu.co/detalle/article/microrred-
generaria-energia-a-bajo-costo-en-zonas-no-interconectadas.html

18
promedio entre 4 a 5,5 kWh/m2 5. Ante estos datos es relevante considerar la
generación de energía a través de sistemas que usen la energía solar, como lo son
las celdas fotovoltaicas y el motor Stirling. Siendo este último capaz de alcanzar una
eficiencia cercana a la teórica del ciclo de Carnot, el cual funciona por medio de un
fluido confinado que está sometido a un proceso de compresión y expansión a
través de una fuente de calor externa. Por consiguiente, observando las
características geográficas de las zonas no interconectadas y el modo de operación
del equipo este trabajo tiene como objetivo la construcción de un motor Stirling tipo
alfa que utiliza como fuente de alimentación un disco parabólico.

5
Instituto de Planificación y Promoción de Soluciones Energéticas para las Zonas No
Interconectadas, Op. cit, p. 54.

19
 1. PLANTEAMIENTO DEL PROBLEMA

En la actualidad Colombia cuenta con una sólida matriz de generación energética

supliendo en gran medida la energía requerida por la población. Según las cifras de
XM 6 administrador del sistema eléctrico nacional, Colombia durante el 2016 tuvo
una capacidad de generación de 16.594 MW mientras que la demanda máxima fue
de 9.904 MW, estando la capacidad de generar un 40% por encima de lo requerido.
No obstante, a pesar de que el país cuenta con un buen panorama de generación,
no toda la población tiene el acceso a la red eléctrica nacional debido a la dificultad
que existe en la transmisión de energía a zonas apartadas y con difícil acceso. De
tal modo según Jorge Valencia 7, exdirector de la UPME, existen alrededor de
470.000 viviendas sin servicio a lo largo del territorio nacional.

En Colombia, las principales fuentes de energía son la hidroeléctrica y la

termoeléctrica, abarcando aproximadamente un 90% de la generación. Estas
fuentes de energía eléctrica no habían tenido ningún problema hasta los últimos
años, ya que el país cuenta con abundantes fuentes hídricas. Sin embargo, producto
de las altas emisiones de CO2 se han generado cambios climáticos los cuales han
provocado que los fenómenos como el del Niño, se hayan acrecentado generando
temporadas de sequías más largas, y por ende racionamiento de energía 8. La
gerente general de XM, María Nohemí Arboleda Arango, señalo que la dependencia
a la hidroelectricidad exige al máximo la operación de la generación térmica y por
ello es importante diversificar la matriz del sistema energético, para disminuir la
vulnerabilidad a los impactos que sufre la principal fuente de energía 9. Es por esto

6
SAAVEDRA CRESPO, Mónica Andrea. Generación eléctrica en Colombia está holgada ante la
demanda [en línea]. EL MUNDO.COM. Agosto, 2017. [Consultado el 15 de mayo de 2018] Disponible
en Internet: http://www.elmundo.com/noticia/Generacion-electrica-en-Colombia-esta-holgada-ante-
la-demanda/357381
7
SÁENZ, Op. cit. Disponible en Internet: https://www.elespectador.com/economia/en-penumbras-
25-millones-de-colombianos-no-tienen-energia-electrica-articulo-706892
8
AHUMADA Omar. Recomiendan al Gobierno racionamiento programado de energía del 5 %. 2016
[en línea]. En: El Tiempo, 6 de marzo de 2016. [Consultado: 5 de mayo del 2018]. Disponible en
http://www.eltiempo.com/archivo/documento/CMS-16531353
9
SAAVEDRA CRESPO Mónica A. Op. cit. Disponible en Internet:
http://www.elmundo.com/noticia/Generacion-electrica-en-Colombia-esta-holgada-ante-la-
demanda/357381

20
que se tienen proyectos que buscan la generación de energía a partir de fuentes
como la eólica y solar.

Una de las metas del gobierno es poder llevar el servicio de energía eléctrica a una
de cada tres viviendas que no cuenta con el servicio de red eléctrica. De tal modo
la Unidad de Planeación Minero Energética (UPME), con el apoyo del Ministerio de
Minas y Energía (MME) y el Instituto de Planificación y Promoción de Soluciones
Energéticas para las ZNI (IPSE), lideraron en el 2016 un proceso de consultoría con
el objetivo de estructurar esquemas empresariales que permitan ampliar la
cobertura 10. El resultado obtenido de los modelos financieros de los esquemas
empresariales propuestos por la consultoría no logro un cierre financiero del
proyecto, por ende, se requiere de soluciones técnicas menos robustas
económicamente.

Con base al escenario que se vive en el país se concluye que se requiere un sistema
que pueda satisfacer la necesidad que tienen las personas que habitan en dichas
zonas no interconectadas. De tal manera, ante esta situación este trabajo tiene
como objetivo diseñar y construir un Motor Stirling que utiliza como fuente de calor
energía solar, con el fin de analizar la eficiencia que tiene el sistema para ser una
solución alternativa al déficit energético. Sugiriendo así la pregunta problema de
este proyecto de grado: ¿Qué capacidad energética puede generar un prototipo de
Motor Stirling Solar tipo alfa?

10
UNIDAD DE PLANEACIÓN MINERO ENERGÉTICA. Informe de rendición de cuentas UPME 2016
2017 [en línea]. Unidad de Planeación Minero Energética: Octubre, 2017, p.29. [Consultado: 15 de
mayo de 2018]. Disponible en Internet:
http://www1.upme.gov.co/InformesGestion/Memorias_Rendicion_2017.pdf. ISSN 2463-123X

21
 2. JUSTIFICACIÓN

Colombia cuenta con los recursos naturales necesarios para generar electricidad a
partir de las fuentes no convencionales de energía renovable. Según la Unidad de
Planeación Minero Energética 11, el país cuenta con una irradiación solar promedio
de 194 W/m2 por día de luz, con vientos localizados de velocidades medias en el
orden de 9 metros por segundo (m/s) (a 80 metros de altura para el caso particular
de La Guajira), y posibilidades energéticas del orden de 450.000 Terajules (TJ) por
año en residuos de biomasa.

Teniendo en cuenta lo anterior, es importante pensar en el aprovechamiento de esta

energía limpia, renovable y sin costo; e intentar darle un uso para la generación de
energía mecánica la cual puede ser posteriormente transformada en energía
eléctrica.

Es por esto que surge la necesidad de plantear en este trabajo de grado el diseño
y la construcción de un motor Stirling solar, el cual es una máquina capaz de
convertir el calor de la energía solar en energía mecánica. El funcionamiento se
debe al confinamiento de un gas sometido a un delta de temperatura que provoca
la expansión y compresión constante del gas y, por consiguiente, el desplazamiento
de un pistón de potencia el cual junto con otro pistón pueden obtener un movimiento
rotacional que entrega potencia. Este motor cuenta con las principales
características de utilizar energías renovables, no generar emisiones a la atmósfera,
y tener un diseño versátil.

Por lo tanto, la construcción de este prototipo permitirá conocer el funcionamiento

del motor Stirling desde una perspectiva real (física), la cual servirá como referencia
para futuras investigaciones, ya que los resultados que se obtengan en este
proyecto serán los primeros pasos para seguir analizando este sistema de
generación de energía que tiene gran potencial, puesto que representa una opción
interesante en esta época de búsqueda de fuentes alternas a las tradicionales que
sean eficientes, económicas y sobre todo limpias.

11
El Tiempo. Energías renovables, con gran potencial para crecer en el país [en línea]. En: El
Tiempo: 22 de abril de 2018. [Consultado: 20 de mayo del 2018]. Disponible en Internet:
https://www.eltiempo.com/vida/medio-ambiente/las-energias-renovables-tienen-un-gran-potencial-
para-crecer-en-colombia-208116

22
 3. OBJETIVOS

3.1 OBJETIVO GENERAL

Diseñar y construir un Motor Stirling tipo Alfa alimentado por energía solar con
materiales de bajo costo como alternativa de solución para la generación de energía
eléctrica renovable.

3.2 OBJETIVOS ESPECÍFICOS

• Realizar los cálculos pertinentes para establecer las dimensiones que adoptaría
el Motor Stirling Solar.

• Diseñar el sistema mecánico y el componente de concentración de la energía

solar usando las herramientas CAD para realizar un estudio de movimiento del
Motor Stirling Solar.

• Construir el sistema compuesto por el motor y el colector solar con materiales de

bajo costo para obtener la información de la cantidad de energía generada.

23
 4. ESTADO DEL ARTE

Tabla 1. Estado del arte

Parámetros
Título Objetivos Resultados
de diseño

General: Motor:
A partir de la ejecución del
Diseñar y fabricar a • Volumen
proyecto, los autores
escala un motor Stirling de
obtuvieron una eficiencia
accionado por energía desplazado
de Carnot de 53%, al igual
solar. : 49cm3
que una salida de potencia
• Fluido: Aire del modelo isotérmico de
Específicos: • Temperatur 65W. Por otro lado, del
a máxima: análisis simple se obtuvo
Diseño y - Definir los parámetros 638K una eficiencia de 7% y
fabricación a de funcionamiento del • Temperatur salida de potencia de 20W,
escala de un motor. a mínima: ya que los
motor Stirling - Seleccionar el 300K intercambiadores de calor
accionado concentrador de • Tipo: Alfa no son ideales, sumado a
por energía energía solar.
las pérdidas por bombeo.
solar 12 - Diseñar el motor
Fuente de Con respecto a la
Stirling. energía: funcionalidad consideran
• Lente que es necesario procesos
- Fabricar y ensamblar Fresnel. de mecanizados altamente
el prototipo.
precisos.
- Especificar el costo
del motor.
- Evaluar
comportamiento y
funcionalidad.

12
AGREDA Víctor, MARASCA Tony. Diseño y fabricación a escala de un motor Stirling accionado
por energía solar [en línea]. Trabajo de investigación Ingeniero Mecánico. Bárbula, Venezuela
Universidad de Carabobo. Facultad de Ingeniería. Escuela de Ingeniería Mecánica. Mayo, 2013.
[Consultado en 15 de abril de 2018]. Disponible en internet:
https://es.scribd.com/document/207121829/DISENO-Y-FABRICACION-A-ESCALA-DE-UN-
MOTOR-STIRLING-ACCIONADO-POR-ENERGIA-SOLAR

24
 Tabla 1. (Continuación)
General Motor:
• Fluido:
Construir un sistema Aire
prototipo conformado • Temperatur
básicamente por un a máxima:
concentrador de energía 150°C El motor Stirling construido
solar y un motor Stirling • Temperatur entregó una potencia
alimentado por dicho a mínima: máxima (estimada) de 7,74
concentrador, que 27°C Wa
permita obtener energía • Tipo: 859 RPM y con una
mecánica a partir de Beta temperatura máxima de
energía solar. 300°C en las pruebas con
Fuente de gas natural como fuente de
Específicos: energía: energía. En las pruebas
con el concentrador solar,
- Conocer las Concentrador la
Diseño de un
características que parabólico de temperatura máxima fue de
posee la energía solar satélite TV. 233°C con lo que se obtuvo
sistema una potencia de 7,35 W a
que llega a la ciudad
prototipo de
conversión
de Medellín para • Recubierto 844 RPM.
establecer su en cromo
de energía
potencial de uso. pulido.
solar en
potencia • Reflectivida
- Diseñar un sistema d: 0,65.
mecánica
que permita convertir
compuesto
la energía solar en
por un
energía térmica y
concentrador
disponer de una
solar y motor
fuente de alta
Stirling 13.
temperatura.

- Diseñar el prototipo de
una máquina térmica
tipo Stirling que
permita convertir la
energía térmica

13
GARCÍA José Luis. Diseño de un sistema prototipo de conversión de energía solar en potencia
mecánica compuesto por un concentrador solar y motor Stirling [en línea]. Proyecto de grado
Ingeniero Mecánico. Medellín. Universidad EAFIT. Facultad de Ingeniería. Departamento de
Ingeniería, 2010. [Consultado el 10 de abril de 2018]. Disponible en internet:
https://repository.eafit.edu.co/handle/10784/4273

25
 Tabla 1. (Continuación)
obtenida en trabajo
mecánico.

- Construir un prototipo
funcional conformado
por un concentrador
de energía solar y un
motor Stirling.

- Realizar pruebas de
los componentes por
separado y del
sistema en conjunto
para verificar su
operación y
rendimiento.

- Evaluar y concluir.

General: Motor: Según lo establecido por el

Diseño y autor, el prototipo no dio
construcción Diseñar y construir un • Tipo: ningún resultado que se
de un motor motor Stirling tipo beta Beta pudiera medir, a pesar de
Stirling tipo acoplado a un • Fluido: que el sistema de
beta concentrador solar Aire alimentación alcanzó a
acoplado a parabólico. brindar una temperatura de
un Fuente de 300°C.
concentrador Específicos: energía.
solar
parabólico - Generar energía • Concentra
cóncavo 14. mecánica con el uso dor
del motor Stirling parabólico
mediante .
concentración solar.

- Convertir la radiación
solar en energía de

14
RODRÍGUEZ José Arturo. Diseño y construcción de un motor Stirling tipo beta acoplado a un
concentrador solar parabólico cóncavo. [En línea]. Trabajo de investigación Ingeniero Mecánico.
Universidad de San Carlos de Guatemala. Facultad de Ingeniería. Escuela de Ingeniería mecánica,
2015. [Consultado el 5 de abril de 2018]. Disponible en internet:
http://www.repositorio.usac.edu.gt/3173

26
 Tabla 1. (Continuación)
alta temperatura como
fuente de calor
externa del motor
Stirling.

- Contribuir como
proyecto innovador en
el campo de la energía
termo solar en
Guatemala y ser
referencia para futuros
estudios y desarrollos
tecnológicos aplicados
al campo de la
termodinámica.
General: Motor:

Utilizar las herramientas • Potencia

adecuadas y la estimada:
información necesarias 25W
para crear el prototipo de
un motor Stirling que • Temperatur
estará diseñado para la a máx.: Se concluye por parte de
generación alternativa de 500°C. los autores que no se logró
energía eléctrica, con el el funcionamiento del
Diseño y propósito de proporcionar • Temperatur prototipo, sin embargo, se
prototipo de un medio para obtener a min: obtuvo un mayor
un Motor energía eléctrica con 50°C. conocimiento de las
Stirling 15 menos emisiones fuentes alternativas.
contaminantes y un motor • Velocidad:
altamente eficiente. 650 RPM

Específicos: Fuente de
energía:
- Elaborar una lista de
ventajas y desventajas • Gas
que muestre los butano

15
HERNANDEZ Miguel, SANCHEZ José Salvador. Diseño y prototipo de un motor Stirling. [En línea].
Trabajo de investigación Ingeniero Mecánico. Universidad Nacional Autónoma de México. Facultad
de Ingeniería. 2013 [Consultado el 20 de marzo de 2018]. Disponible en internet:
http://132.248.52.100:8080/xmlui/handle/132.248.52.100/4961

27
 Tabla 1. (Continuación)
beneficios de un motor
Stirling.

- Diseñar y crear un
prototipo de motor
Stirling a partir de un
análisis
termodinámico.

- Mostrar aplicaciones
próximas para
motores Stirling como
generación alternativa
de energía eléctrica.

Para la
realización del
proyecto: Por medio de las
conclusiones el autor
• Se establece que el prototipo
Diseño y Objetivo General: establece es funcional y cumple con
construcción partes del el objetivo planteado.
de un Diseñar y construir un colector Durante las pruebas
prototipo de prototipo de un solar experimentales se lograron
concentrador concentrador solar PDR. temperaturas alrededor de
solar parabólico de disco
• Relación 400°C, a pesar de que se
parabólico de reflector con dos ejes de
de realizaron durante días
disco seguimiento y de un
concentra nublados. De tal modo se
reflector para diámetro de 1.5m, que
ción en un estima temperaturas
generación proporcione energía
colector mayores a 600°C en
de energía térmica de alta calidad
PDR. mejores condiciones.
térmica 16. (elevada temperatura).
• Concentra
dor.
• Sistema
de control.
• Sistema
mecánico.

16
DURÁN MEDINA, Pino; BARBOSA SALDAÑA, Juan Gabriel. Diseño y construcción de un
prototipo de concentrador solar parabólico de disco reflector para generación de energía térmica [En
línea]. Instituto Politécnico Nacional. Escuela superior de Ingeniería Mecánica y Eléctrica, 2012. 54
p. [Consultado: 15 de mayo del 2018]. Disponible en
http://tesis.ipn.mx/jspui/handle/123456789/17648

28
 5. MARCO TEÓRICO

5.1 MOTOR STIRLING

5.1.1 Ventajas y desventajas del motor Stirling

Tabla 2. Ventajas y desventajas del motor Stirling

VENTAJAS DESVENTAJAS

El aporte de calor es externo por lo que

Baja densidad de potencia debido a la
sus condiciones de combustión son
flexibles combustión externa

Funciona con cualquier fuente de calor Lo condiciona su tamaño

La construcción puede resultar difícil

Reduce la mayor parte de las emisiones
para que tenga una larga vida útil

Largo tiempo de encendido y apagado

Requieren de poca lubricación
del motor

Lento tiempo de respuesta, es decir

no tiene un arranque inmediato, sino
Cuenta con mecanismos sencillos que pasado un tiempo de
calentamiento comienza a ganar
velocidad.

Usa un fluido de trabajo de una única

fase, manteniendo las presiones internas Se requieren grandes superficies de
cercanas a la presión de diseño y reduce intercambios de calor
los riesgos de explosión

29
5.1.2 Clasificación de los motores Stirling:

Existen 3 tipos de motores o configuraciones de motores Stirling: Alfa, Beta y

Gamma. La diferencia que manejan estas configuraciones se basa en la
complejidad del mecanismo con el cual se diseñe, pero tienen la misma finalidad. A
continuación, se muestra en la Figura 1 los diseños de los mecanismos de las
diferentes configuraciones.

Figura 1. Motores Stirling tipo Alpha-Beta-Gamma

alfa beta gamma

Fuente: PASCUAL TOLEDANO, Juan Ignacio. Configuraciones Alfa, Beta y Gamma

[Imagen]. Simulación de un motor Stirling mediante modelos de primer orden [Figura].
España: Universidad Carlos III de Madrid, 2012. p. 24. [Consultado: 10 de febrero del
2018]. Disponible en Internet: https://e-archivo.uc3m.es/handle/10016/18099#preview

5.1.2.1 Motor Stirling tipo Alfa

Esta configuración cuenta con dos pistones que se encuentran en cilindros

independientes que están unidos por medio de un ducto en el cual se desplaza el
fluido de la zona caliente a la zona fría y viceversa. No cuenta con un desplazador,
pero sus pistones están desfasados a 90°. Normalmente este tipo de motor es el
más simple a comparación de las otras configuraciones, sin embargo, su dificultad
se basa en que los pistones deben de tener sellos para contener el gas y evitar
fugas. Adicional en la zona fría cuenta con aletas para mejor su transferencia de
calor. 17

17
PASCUAL TOLEDANO, Juan Ignacio. Configuraciones Alfa, Beta y Gamma [Imagen]. Simulación
de un motor Stirling mediante modelos de primer orden [en línea]. Trabajo de investigación para
Ingeniero Técnico Mecánico. España. Universidad Carlos III de Madrid. Departamento de ingeniería
30
Figura 2. Configuración alfa de motor Stirling

5.1.3 Mecanismo del motor Stirling tipo alfa

Teniendo en cuenta la investigación que se describe en el título “marco teórico”

sobre el funcionamiento del motor Stirling, los tipos de configuraciones y los
parámetros que se deben considerar para su diseño, el motor que se elige para el
desarrollo del proyecto de grado es el tipo Alfa, por las siguientes características:

• El motor Stirling tipo Alfa se compone de dos cilindros-pistones y un mecanismo

de biela manivela, lo cual lo hace más sencillo para su fabricación. Para el desarrollo
del proyecto se utiliza como ejemplo el mecanismo del sistema del cabezal de un
compresor de aire ya que se compone de los mismos elementos y tiene un desfase
de 90°.

• Teniendo en cuenta que el motor tipo Alfa maneja en las dos zonas (fría y
caliente) el mismo mecanismo (cilindro-pistón), el volumen desplazado es igual en
ambas zonas, es decir que se tiene una relación de volumen mayor en comparación
de los otros tipos de motores Stirling.

térmica y fluidos, 2012. p. 24. [Consultado: 10 de febrero del 2018]. Disponible en Internet: https://e-
archivo.uc3m.es/handle/10016/18099#preview

31
• Para estimar teóricamente la potencia que obtendrá el motor Stirling a diseñar
se requiere del modelo termodinámico de primer orden de Schmidt, el cual para la
configuración Alfa es más sencilla y precisa.

A continuación, se ilustra por medio del Cuadro 1 el comportamiento del fluido de

trabajo con el mecanismo del motor Stirling el cual genera el ciclo para una
configuración tipo Alfa:

32
Cuadro 1. Ciclo Stirling del motor tipo Alfa

Por medio de la energía solar que

es captada por medio de un disco
parabólico, se calienta el gas de
trabajo hasta llegar a su
temperatura máxima generando
una expansión de este. El pistón
se desplaza hacia abajo
transmitiendo movimiento al
cigüeñal. El gas de trabajo se
desplaza por el regenerador hacia
la zona fría

Cuando el gas de trabajo va

ocupando el espacio en el cilindro
frio, el mecanismo de biela
manivela genera un empuje del
cigüeñal el cual permite desplazar
el pistón caliente de su punto
inferior a su punto superior. En
este punto del proceso la presión
cae ya que la temperatura
disminuye por la transferencia de
calor.

33
 Cuadro 1. (Continuación)

El pistón de la zona caliente se

desplaza hacia la parte superior
(compresión) cuando el gas de
trabajo ocupa todo el volumen del
cilindro frio, gracias al proceso de
transferencia de calor ya que su
temperatura disminuye y la
presión cae. En este proceso el
volante de inercia ayuda a la
compresión del gas.

El volumen del gas de trabajo en

la zona fría ya es mínimo y el
pistón caliente se encuentra en su
punto superior. El gas se caliente
de nuevo para iniciar el empuje del
pistón caliente hacia el punto
inferior. Nuevamente inicia el ciclo.

Fuente: SCHLAICH BERGERMANN, &. Partner GbR. EuroDish - Stirling System

Description [Tabla]. Stuttgart, junio 2001. p. 6. [Consultado: 5 marzo del 2018]. Disponible
en Internet: http://pointfocus.com/images/pdfs/eurodish.pdf.

34
5.1.4 Ciclo Stirling ideal:

Los motores térmicos se encargan de transformar el trabajo producido por un flujo

de calor por medio de una seria de procesos termodinámicos. Estos conjuntos de
procesos generan lo denominado “ciclo”. El rendimiento máximo al cual puede llegar
un ciclo termodinámico es el ciclo de Carnot. Este ciclo es una configuración ideal
en donde no existen perdidas por ningún fenómeno por lo tanto el ciclo Stirling es
con base al ciclo de Carnot. Los procesos de evolución del ciclo están compuestos
por 18:

Figura 3. Etapas del ciclo Stirling y su variación de volumen

Fuente: GARCÍA J. Luis. Etapas del ciclo Stirling y su respectiva variación en el

volumen dentro del cilindro [Imagen]. Diseño de un sistema prototipo de conversión
de energía solar en potencia mecánica compuesto por un concentrador solar y
motor Stirling [Figura]. Trabajo de investigación Ingeniero Mecánico. Medellín,
Colombia. Universidad EAFTI. Facultad de Ingeniería. Departamento de Ingeniería
Mecánica, 2010. [Consultado el 10 de abril de 2018]. Disponible en internet:
https://repository.eafit.edu.co/handle/10784/4273

18
AGÜERO ZAMORA, Víctor Raúl. Diseño y construcción de un motor Stirling para la generación de
energía eléctrica. [en liena] Proyecto de grado Ingeniero Mecánico. Lima Perú.Universidad Nacional
de Ingeniería. Facultad de Ingeniería Mecánica. 2011. P. 42. [Consultado: 5 marzo del 2018].
Disponible en Internet: http://cybertesis.uni.edu.pe/handle/uni/756

35
1->2: Expansión isotérmica con absorción de calor (Isoterma)
𝑉𝑉2
𝑄𝑄1−2 = 𝑅𝑅 ∗ 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 ∗ 𝐿𝐿𝐿𝐿 � �
𝑉𝑉1 (1)

𝑉𝑉2
𝑊𝑊1−2 = 𝑃𝑃1 ∗ 𝑉𝑉1 ∗ 𝐿𝐿𝐿𝐿 � �
𝑉𝑉1 (2)

𝑉𝑉2
𝑠𝑠2 − 𝑠𝑠1 = 𝑅𝑅 ∗ 𝐿𝐿𝐿𝐿 � �
𝑉𝑉1 (3)

2->3: Extracción de calor de regeneración a volumen constante (Isocora)

𝑄𝑄𝑟𝑟𝑟𝑟𝑟𝑟= 𝑄𝑄1−2 = 𝑐𝑐𝑐𝑐 ∗ (𝑇𝑇3 − 𝑇𝑇2)
(4)

𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇
𝑠𝑠3 − 𝑠𝑠2 = 𝑐𝑐𝑐𝑐 ∗ 𝐿𝐿𝐿𝐿 � �
𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 (5)

3->4: Compresión isotérmica con extracción de calor (Isoterma)

𝑃𝑃3
𝑄𝑄3−4 = 𝑅𝑅 ∗ 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 ∗ 𝐿𝐿𝐿𝐿 � �
𝑃𝑃4 (6)

𝑃𝑃3
𝑊𝑊3−4 = 𝑃𝑃3 ∗ 𝑉𝑉3 ∗ 𝐿𝐿𝐿𝐿 � �
𝑃𝑃4 (7)

𝑉𝑉4
𝑠𝑠4 − 𝑠𝑠3 = 𝑅𝑅 ∗ 𝐿𝐿𝐿𝐿 � �
𝑉𝑉3
(8)

4->1: Aportación de calor de regeneración a volumen constante (Isocora)

𝑄𝑄4−1 = 𝑐𝑐𝑐𝑐 ∗ (𝑇𝑇1 − 𝑇𝑇4)
(9)

36
 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇
𝑠𝑠3 − 𝑠𝑠4 = 𝑐𝑐𝑐𝑐 ∗ 𝐿𝐿𝐿𝐿 � �
𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 (10)

Figura 4. Ciclo ideal Stirling

Fuente: PASCUAL TOLEDANO, Juan Ignacio. Simulación de un motor Stirling

mediante modelos de primer orden. Proyecto de grado Ingeniería Técnica
Mecánica. España. Universidad Carlos III de Madrid. Facultad de Ingeniería. 2012.
Pág. 10.

5.1.5 Ciclo Stirling real:

En los sistemas desarrollados térmicamente la eficiencia más alta a la que se puede

llegar es a la eficiencia por el ciclo de Carnot. El ciclo real se diferencia del ideal en
algunos aspectos:

• El movimiento del pistón es usualmente sinusoidal en vez de discontinuo,

llevando a la manipulación no óptima del fluido del trabajo.

37
• Los procesos de expansión y compresión en el motor Stirling real son más bien
politrópicos en vez de isotérmicos. Esto causa fluctuaciones en la presión y
temperatura en el gas de trabajo lo que resulta en procesos adiabáticos y perdidas
transitorias de calor.

• Se presentan perdidas debido a la fricción del gas durante el desplazamiento

particularmente durante el flujo del regenerador.

• Al regenerador y a los intercambiadores de calor de la máquina Stirling real sí

tienen volumen a diferencia de los de la máquina ideal.

• Debido a que el mecanismo de los motores Stirling es complejo en su proceso

de manufactura, el movimiento no es totalmente ideal entre el pistón y el
desplazador y que además el proceso de transferencia de calor no llega a realizar
ciclos netamente isotérmicos, la gráfica del ciclo real es un ciclo “redondeado” en
forma de elipse.

Figura 5. Ciclo ideal vs Ciclo real

Fuente: RODRIGUEZ, J. Arturo. Ciclo ideal vs ciclo ideal [imagen]. Diseño y

construcción de un motor Stirling tipo beta acoplado a un concentrador solar
parabólico cóncavo [en línea]. Trabajo de investigación Ingeniero Mecánico.
Guatemala. Universidad de San Carlos de Guatemala. Facultad de Ingeniería. 2015.
30 p. [Consultado: 12 marzo de 2018]. Disponible en Internet:
http://www.repositorio.usac.edu.gt/3173/

38
5.1.6 Modelos termodinámicos

Los modelos termodinámicos permiten la estimación de la cantidad de trabajo que

va a realizar el motor Stirling proporcionando las dimensiones de cada parte. Para
ello se cuenta con tres modelos: primer, segundo y tercer orden, para los cuales
corresponden a un análisis de conversación de masa, energía y movimiento lineal.
A continuación, se describen los modelos que permiten el desarrollo del diseño del
motor Stirling:

5.1.6.1 Modelos de Primer Orden

Estos modelos están basados en condiciones ideales ya que utilizan factores de

ajuste para poder estimar la potencia y la eficiencia. Estos resultados no son
suficientemente confiables para determinar el rendimiento neto ni dimensiones
finales puesto que sería conveniente hacer los otros modelos termodinámicos. Sin
embargo, para desarrollar prototipos estos modelos de primer orden permiten dar
una estimación de las dimensiones iniciales y potencia generada como primera
instancia para identificar las variables críticas y poder establecer mejoras.

• Fórmula de Beale:

Un criterio tradicional consiste en estimar la potencia máxima mediante el número

de Beale. William Beale observó que el número de Beale indica que
aproximadamente todos los modelos de motores Stirling tienen un número
adimensional que está relacionado con la temperatura del foco caliente, como se
muestra en el gráfico 1. La expresión matemática que involucra el número de Beale
es:
𝑃𝑃 = 𝐵𝐵𝐵𝐵. 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃. 𝑉𝑉𝑉𝑉. 𝑓𝑓 (11)

P es la potencia (Watts),
B es el número de Beale
Vo volumen en el cilindro de expansión (cm3)
f es la frecuencia del motor (Hertz)
P media es la presión media del ciclo (bar)

En este gráfico el valor de Beale puede verse afectado también por el tamaño,
diseño y eficiencia que presente el motor. Por lo tanto, la línea punteada inferior es
la guía para determinar este número en motores pequeños y con eficiencia baja. La

39
línea superior es para motores grande con eficiencia alta. Y la línea central es la
estándar para la mayoría de motores.

Gráfico 1. Curva para determinar el número de Beale en función de la temperatura

Fuente: WALKER Graham. Número de Beale en función de la temperatura

[Imagen]. Stirling engine [Grafica]. Lugar de publicación: New York 1980, Oxford
University Press. p. 74. [Consultado: 10 de marzo de 2018]. Disponible en Internet:
https://es.scribd.com/doc/45062130/Stirling-Engines-G-Walker-Oxford-1980-WW

• Teoría de Schmidt:

El análisis del ciclo Stirling anteriormente descrito, es un análisis térmico simple el

cual solo se tiene dos fases de compresión isotérmicas y dos procesos
regenerativos con volúmenes constantes. Es por esto que Gustav Schmidt (1871)
desarrolló un análisis clásico que suministra una aproximación más razonable del
rendimiento del motor. Sin embargo, el método sigue conservando parámetros
idealizados ya que la teoría estudia el comportamiento del movimiento del
mecanismo con base a la suposición de una compresión y expansión isotérmica y
una regeneración perfecta.

40
No obstante, el análisis del ciclo de Schmidt es una herramienta útil y sencilla para
poder determinar la variación de los volúmenes de los espacios de trabajo en
función del ángulo del eje del cigüeñal. A continuación, se muestra la teoría para el
desarrollo de los cálculos por el método de Schmidt:

Figura 6. Esquema de la teoría de Schmidt

Fuente: URIELI Israel; Berchowitz David M. Stirling Cycle Machine Analysis

[Figura]. Bristol: A. Hilger, 1984 [Consultado: 10 de abril de 2018]. Disponible en
Internet: https://www.ohio.edu/mechanical/stirling/isothermal/Schm_summary.html

El ciclo de Schmidt es un ciclo el cual el pistón de potencia y el desplazador

(dependiendo de la configuración) tienen un movimiento sinusoidal, como se puede

41
observar en la figura 6. Para el desarrollo de este se tienen las siguientes
suposiciones 19:

• El proceso de regeneración es perfecto

• La presión en cada instante es igual en todo el sistema

• El fluido de trabajo obedece a la ley de gases ideales, PV=mRT

• No hay fugas y la masa se mantiene constante

• La variación de volumen se mueve sinusoidalmente

• No existen gradientes de temperatura en los intercambiadores de calor

• La temperatura del pistón y el cilindro son constantes

• La temperatura del fluido de trabajo en el espacio de compresión y expansión son

constantes

Para desarrollar el análisis de Schmidt lo primero que se debe determinar son los
volúmenes en la zona fría y en la caliente para un ángulo de cigüeñal. El volumen
actual, que es igual al volumen en cada instante del ciclo, está determinado por un
ángulo X. (el cual es X=0 cuando el pistón de expansión se encuentra en el límite
superior).

El volumen actual de expansión Ve y el volumen actual de compresión están

determinados por las siguientes ecuaciones:

𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉
𝑉𝑉𝑉𝑉 = (1 + cos(𝑥𝑥)) + 𝑉𝑉𝑉𝑉𝑙𝑙𝑙𝑙 (12)
2

19
WALKER Graham. Stirling engine [en línea]. New York, Oxford University Press. 1980.
[Consultado: 10 de marzo de 2018]. Disponible en Internet:
https://es.scribd.com/doc/45062130/Stirling-Engines-G-Walker-Oxford-1980-WW

42
 𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉
𝑉𝑉𝑉𝑉 = (1 − cos(𝑥𝑥 − 𝑑𝑑𝑑𝑑)) + 𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 (13)
2

El volumen actual total está definido por la siguiente ecuación:

𝑉𝑉 = 𝑉𝑉𝑉𝑉 + 𝑉𝑉𝑉𝑉 + 𝑉𝑉𝑉𝑉 (14)

La masa total al interior del motor m es calculada usando la presión del motor P y
las temperaturas en las zonas, con la siguiente ecuación:

𝑃𝑃𝑃𝑃𝑃𝑃 𝑃𝑃𝑃𝑃𝑃𝑃 𝑃𝑃𝑃𝑃𝑃𝑃

𝑚𝑚 = + + (15)
𝑅𝑅𝑅𝑅𝑅𝑅 𝑅𝑅𝑅𝑅𝑅𝑅 𝑅𝑅𝑅𝑅𝑅𝑅

Los principales parámetros de diseño de motor Stirling son la relación de

temperaturas t, la relación de volúmenes de barrido v y las relaciones de volúmenes
muertos ya que por medio de estos parámetros se establece la presión media y la
potencia de salida.

𝑇𝑇𝑇𝑇
𝑡𝑡 = (16)
𝑇𝑇ℎ

𝑉𝑉𝑉𝑉𝑉𝑉
𝑣𝑣 = (17)
𝑉𝑉𝑉𝑉𝑉𝑉

𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉
𝑋𝑋𝑋𝑋𝑋𝑋 = (18)
𝑉𝑉𝑉𝑉𝑉𝑉

𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉
𝑋𝑋𝑋𝑋𝑋𝑋 = (19)
𝑉𝑉𝑉𝑉𝑉𝑉

𝑉𝑉𝑉𝑉
𝑋𝑋𝑋𝑋 = (20)
𝑉𝑉𝑉𝑉𝑉𝑉

43
La temperatura del regenerador se considera como el promedio entre la temperatura
de la zona fría y la zona caliente.

(𝑇𝑇ℎ − 𝑇𝑇𝑇𝑇)
𝑇𝑇𝑇𝑇 =
𝑇𝑇ℎ (21)
ln � �
𝑇𝑇𝑇𝑇

Para determinar la masa total se tiene en cuenta los volúmenes con sus respectivas
temperaturas y se sustituyen en la ecuación de gases ideales como se puede
evidenciar de la siguiente forma:

𝑀𝑀 = 𝑚𝑚𝑚𝑚 + 𝑚𝑚𝑚𝑚 + 𝑚𝑚𝑚𝑚 + 𝑚𝑚ℎ (22)

𝑃𝑃𝑃𝑃
𝑀𝑀 = (23)
𝑅𝑅𝑅𝑅

Donde:

𝑉𝑉𝑉𝑉 𝑉𝑉𝑉𝑉 𝑉𝑉ℎ

𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 ∗ � �+� �+� �
𝑀𝑀 = 𝑇𝑇𝑇𝑇 𝑇𝑇𝑇𝑇 𝑇𝑇ℎ (24)
𝑅𝑅

A continuación, se presentan las ecuaciones que determinan los coeficientes a, S,

B y c:

𝑣𝑣. sin(𝑑𝑑𝑑𝑑)
𝑎𝑎 = 𝑡𝑡𝑡𝑡𝑡𝑡−1 � �
𝑡𝑡 + cos(𝑑𝑑𝑑𝑑) (25)

4. 𝑡𝑡. 𝑋𝑋𝑋𝑋
𝑆𝑆 = 𝑡𝑡 + 2. 𝑡𝑡. 𝑋𝑋𝑋𝑋𝑋𝑋 + + 𝑣𝑣 + 2. 𝑋𝑋𝑋𝑋𝑋𝑋 (26)
1 + 𝑡𝑡

44
 𝐵𝐵 = �𝑡𝑡 2 + 2. 𝑡𝑡. 𝑣𝑣. cos(𝑑𝑑𝑑𝑑) + 𝑣𝑣 2 (27)

𝐵𝐵
𝑐𝑐 = (28)
𝑆𝑆

Teniendo en cuenta el cálculo de la presión P, se determina por medio de un arreglo

matemático la presión promedio que es equivalente a la presión media que
posteriormente se utilizara para los cálculos de potencia con respecto a un análisis
de primer orden. Esta presión se calcula con el fin de determinar un valor entero en
todo el ciclo considerando que la presión es constante.

𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃√1 − 𝑐𝑐 2
𝑃𝑃 = (29)
1 − 𝑐𝑐. cos(𝑥𝑥 − 𝑎𝑎)

1 2. 𝑚𝑚. 𝑅𝑅. 𝑇𝑇𝑇𝑇

𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 = ∮ 𝑃𝑃𝑃𝑃𝑃𝑃 =
2𝜋𝜋 𝑉𝑉𝑉𝑉𝑉𝑉√𝑆𝑆 2 − 𝐵𝐵2 (30)

La potencia y la eficiencia del motor Stirling está determinada por las siguientes
ecuaciones:

𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃. 𝑉𝑉𝑉𝑉𝑉𝑉. 𝜋𝜋. 𝑐𝑐. sin(𝑎𝑎)

𝑊𝑊𝑊𝑊 = (31)
1 + √1 − 𝑐𝑐 2

𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃. 𝑉𝑉𝑉𝑉𝑉𝑉. 𝜋𝜋. 𝑐𝑐. 𝑡𝑡. sin(𝑎𝑎)

𝑊𝑊𝑊𝑊 = (32)
1 + √1 − 𝑐𝑐 2

𝑊𝑊𝑊𝑊 = 𝑊𝑊𝑊𝑊 + 𝑊𝑊𝑊𝑊 (33)

𝑊𝑊𝑊𝑊 ∗ 𝑛𝑛
𝐿𝐿𝐿𝐿 = (34)
60

45
 𝑊𝑊𝑊𝑊
𝑒𝑒 = (35)
𝑊𝑊𝑊𝑊

5.1.6.2 Modelos de Segundo y Tercer Orden

Estos métodos son más complejos, desde el punto de vista físico y matemático ya
que analizan los procesos termodinámicos que suceden dentro del motor utilizando
la ecuación de estado, balance de energía y continuidad, con el fin de poder
optimizar el rendimiento. Para ello se requiere de un software de cálculo que
solucione las ecuaciones modelándolo para los diferentes instantes de tiempo del
ciclo del motor.

La ventaja que tienen estos modelos de Segundo y Tercer Orden es que el análisis
es un poco más detallado para determinar la magnitud de las pérdidas de
temperatura tanto en él foco frio como en él caliente y la caída de presión.

5.2 COLECTOR SOLAR

5.2.1 Sistema disco parabólico

El sistema de disco parabólico tiene los siguientes componentes:

• Concentrador: Parte del colector que dirige la radiación hacia el receptor. La

temperatura proveniente de la energía acumulada en el colector se puede
incrementar por medio de la reducción del área por la cual ocurren perdidas de calor.

• Receptor: Elemento del sistema donde la radiación es absorbida y convertida en

algún otro tipo de energía.

• Sistema de generación, (motor de combustión externa Stirling).

• Estructura soporte y mecanismos.

46
Teóricamente la radiación incidente en la dirección normal de la superficie del disco
es concentrada en un punto focal que está a una distancia f del concentrador, el
receptor. Entre el colector y el receptor se maneja un margen de tolerancia del
ángulo comprendido entre ambos, ya que se supone que los rayos de Sol llegan
linealmente a la Tierra.

La geometría del colector solar parabólico tipo disco y los aspectos geométricos
requeridos para el diseño se muestran en la figura 7. Las dimensiones del disco
parabólico se obtienen a partir de la parábola que forma la sección transversal del
concentrador como se puede observar en la misma figura. La geometría de una
parábola con eje focal sobre el eje de las ordenadas esta dado por la siguiente
ecuación:

𝑥𝑥 2 = 4𝑓𝑓𝑓𝑓 (36)
La distancia f que se muestra en la figura 7 es denominada longitud focal, es la
distancia del punto focal al vértice, por otro lado, Da se denomina como el diámetro
de apertura de la parábola. Entre los puntos F y A existe una distancia r que es igual
a la distancia focal, esta distancia incrementa hasta alcanzar un valor máximo rr,
denominado radio de apertura que corresponde a la distancia que existe entre en el
punto F y B.

Figura 7. Sección transversal de un concentrador parabólico

El cálculo de los factures necesarios para el diseño y la construcción del disco se

hacen teniendo en cuenta el diámetro de apertura (Da), foco (f) y el diámetro del

47
receptor. De tal modo el diámetro de apertura y el ángulo de apertura (ϕr) vienen
relacionados por la ecuación (37).

⎡ 8 � 𝑓𝑓 � ⎤
𝐷𝐷𝐷𝐷
𝜙𝜙𝜙𝜙 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 ⎢ ⎥ (37)
⎢ 𝑓𝑓 2 ⎥
⎣ 16 � � − 1 ⎦
𝐷𝐷𝐷𝐷

Otro parámetro necesario para la geometría del colector es el radio de borde (rr) o
máximo valor de distancia entre el punto focal y el extremo del paraboloide, que se
calcula según la ecuación (38), esta permite determinar el valor de la distancia que
existe entre el punto focal y el plato para cualquier ángulo ϕ de la parábola.

2𝑓𝑓
𝑟𝑟𝑟𝑟 = (38)
1 + cos 𝜙𝜙𝜙𝜙

A partir de la relación f/Da y por medio de la ecuación (39) es posible establecer

una relación matemática para determinar la distancia óptima del foco:

𝐷𝐷𝐷𝐷
𝑓𝑓 =
𝜙𝜙𝜙𝜙 (39)
4 tan � �
2

La razón de concentración C se define como la relación entre el área de apertura

(Aa) y el área del recibidor (Ar), e indica el factor o proporción en que se intensifica
la concentración de los rayos solares en el recibidor. El incremento de este valor
significa que existe un aumento de la energía disponible. No obstante, esto también
implica mayor precisión en la calidad óptica (manufactura del disco) así como de su
sistema de posicionamiento.

𝐴𝐴𝐴𝐴
𝐶𝐶 = (40)
𝐴𝐴𝐴𝐴

𝐷𝐷𝑎𝑎2 𝜋𝜋
𝐴𝐴𝐴𝐴 = (41)
4

48
Para determinar el área del recibidor en caso que sea un cilindro se ha de tener en
cuenta el radio a del receptor y la zona de contacto con la radiación de la superficie
h.

𝐴𝐴𝐴𝐴 = 2𝜋𝜋𝜋𝜋ℎ (42)

La radiación solar incide con una amplitud angular de 0,53 °. Si se asume que la
radiación directa es normal a la apertura y que el concentrador es simétrico, en un
caso ideal el ancho de la imagen focal se ve afectada por el ángulo de borde (𝝓𝝓𝝓𝝓 )
y la amplitud angular de la radiación solar. Sin embargo, a nivel practico existe otro
factor que es necesario considerar, el ángulo de dispersión (δ).

EL ángulo de dispersión (δ) considera los errores angulares asociados a un

seguimiento solar inapropiado, rugosidad de la superficie reflectora y mal
conformado en la curvatura del disco. Estos errores permiten determinar el aumento
de la amplitud de la radiación reflejada por el concentrador y han sido obtenidos de
forma experimental haciendo uso de la estadística y se indican en la tabla 3. 20

20
STINE, W. Solar energy fundamentals and design with computer applications. Wiley, Estados
Unidos. 1985. Citado por DURÁN MEDINA, Pino; BARBOSA SALDAÑA, Juan Gabriel. Diseño y
construcción de un prototipo de concentrador solar parabólico de disco reflector para generación de
energía térmica . Obtención del grado de Maestro en ciencias en Ingeniería Mecánica.

49
Tabla 3. Errores típicos para concentradores parabólicos.

Fuente. DURÁN MEDINA, Pino; BARBOSA SALDAÑA, Juan Gabriel. Diseño y

construcción de un prototipo de concentrador solar parabólico de disco reflector para
generación de energía térmica [Tabla]. Instituto Politécnico Nacional. Escuela
superior de Ingeniería Mecánica y Eléctrica, 2012. 54 p. [Consultado: 15 de mayo
del 2018]. Disponible en http://tesis.ipn.mx/jspui/handle/123456789/17648

El valor de la dispersión depende de los errores unidimensionales (σ1d) y

bidimensionales (σ2d). Los errores unidimensionales están medidos en mili
radianes (mrad) son los que influyen en la radiación directa en el plano de curvatura.
Y los errores bidimensionales son los asociados a la parábola cuando los rayos
incidentes no llegan al plano de curvatura, al ángulo de apertura del Sol y a los rayos
reflejados especulares. De tal modo ambos errores se determinan por medio de las
ecuaciones (43 y 44) y la tabla 3.

1
𝜎𝜎1𝑑𝑑 = [(2𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎)2 + (𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎)2 + (𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎)2 + (𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎)2 ]2 (43)

1
(𝜎𝜎𝜎𝜎𝜎𝜎𝑙𝑙2 + 𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝑙𝑙2 )2 (44)
𝜎𝜎2𝑑𝑑 =
cos 𝜃𝜃

50
A partir del resultado de los errores unidimensionales y bidimensionales se
determina el error total con la ecuación (45). Este valor que se obtiene es en
radianes y es considerado como la dispersión que ocurre durante la reflexión en la
parábola.

1
𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎 = (𝜎𝜎1𝑑𝑑 2 + 𝜎𝜎2𝑑𝑑2 )2 (45)

Por ello si se desea el valor en grados se hace uso de la ecuación (46):

180
(𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎) ≈ 𝛿𝛿 (46)
𝜋𝜋

Los valores mostrados de desviación estándar mostrados en la tabla 3 son válido

para una magnitud (1σ), para obtener una mayor precisión del efecto generado por
la desviación angular es posible considerar tres desviaciones estándares (3σ), de
tal modo:

180
(3)(𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎) ≈ 𝛿𝛿 (47)
𝜋𝜋

Una vez que se conoce el valor de la dispersión angular (δ), el ángulo de borde (ϕr)
y la amplitud angular de la radiación del Sol es posible determinar el ancho mínimo
(W) de la imagen focal por medio de la ecuación (48) o el diámetro mínimo (D) por
medio de la ecuación (49) que debe de tener un cilindro 21:

𝛿𝛿
2𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 �0.267 + � ��
2
𝑊𝑊 = (48)
𝛿𝛿
cos �𝜙𝜙𝜙𝜙 + 0.267 + � ��
2

21
BECKMAN, William A.; DUFFIR, John A. Solar Engineering of Thermal Processes. 4 ed. New
Jersey: John Wiley & Sons, Inc, 2013. 359 p. ISBN 978-0-470-87366-3.

51
 𝛿𝛿
𝐷𝐷 = 2 ∗ 𝑟𝑟𝑟𝑟 ∗ 𝑠𝑠𝑠𝑠𝑠𝑠 �0.267 + � �� (49)
2

Figura 8. Dimensión de la imagen generada en el punto focal de un colector de

concentración.

El ángulo de dispersión (δ), el ángulo de borde (ϕr) y el ángulo con el que incide el
Sol en la Tierra influyen en el cálculo de la razón de concentración, permitiéndose
calcular por medio de la ecuación (50) la razón de concentración máxima que se
puede obtener por medio del colector parabólico 22.

𝛿𝛿
𝑠𝑠𝑠𝑠𝑛𝑛2 𝜙𝜙𝜙𝜙 cos 2 �𝜙𝜙𝜙𝜙 + 0.267 + �
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 = 2 −1 (50)
𝛿𝛿
4𝑠𝑠𝑠𝑠𝑛𝑛2 �0.267 + �
2

22
Ibid., p. 367.

52
La eficiencia óptica del colector depende de la geometría del sistema, factor de
concentración y los materiales utilizados en la construcción:

𝑛𝑛𝑛𝑛 = 𝛼𝛼𝛼𝛼𝛼𝛼𝛼𝛼
(51)
Tabla 4. Simbología de ecuación para la eficiencia óptica

Abreviatura Variable
α Absortancia del receptor
ρ Reflectividad del concentrador
τ Transmitancia que tiene el material del receptor
F Factor de forma

A partir de los parámetros ópticos y conociendo la radiación directa promedio (I) en

la zona donde se ubica el concentrador se puede determinar la energía absorbida
por el receptor por medio de la ecuación (52):

𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄 = 𝐴𝐴𝐴𝐴 ∗ 𝜌𝜌 ∗ 𝜏𝜏 ∗ 𝛼𝛼 ∗ 𝐹𝐹 ∗ 𝐼𝐼
(52)

Tabla 5. Simbología de ecuación para la energía absorbida.

Abreviatura Variable
Α Absortancia del receptor
Ρ Reflectividad del concentrador
Transmitancia que tiene el
Τ
material del receptor
F Factor de forma
Aa Área del colector
I Radiación directa promedio

53
 6. METODOLOGÍA
Diagrama 1. Metodología

1. ESTADO DE ARTE DEL MOTOR STIRLING SOLAR

- Condiciones
iniciales de
2. DESARROLLO DE LOS PARÁMETROS operación
INICIALES DEL DISEÑO DEL MOTOR STILING - Volumen
desplazado
Y EL COLECTOR SOLAR - Diámetro del
plato

- Dimensionamiento
del motor
- Selección de
3. DISEÑO DEL MOTOR STIRLING SOLAR materiales
- Estudio de
movimiento en
software CAD
- Diseño del colector

- Selección de
4. CONSTRUCCIÓN DEL MOTOR los procesos de
STIRLING SOLAR manufactura

Análisis por
5. VALIDACIÓN CON EL MODELO MATEMÁTICO el método
de Schmidt

54
 7. DESARROLLO DE PARÁMETROS INICIALES DEL DISEÑO DEL
MOTOR STIRLING Y EL COLECTOR SOLAR

En este capítulo se presenta la metodología de diseño del motor Stirling tipo alfa y
del colector solar teniendo en cuenta los parámetros iniciales y métodos de cálculo
para ello.

7.1 DISEÑO DEL MOTOR STIRLING TIPO ALFA

Diagrama 2. Diseño del motor

7.1.1 Cálculo para el volumen barrido

El diseño de un motor Stirling parte de la obtención del volumen que va a ser

desplazado por el pistón para determinar las dimensiones. Como se menciona en
el capítulo anterior, existen modelos que permiten validar los parámetros teóricos
como lo es la teoría de Schmidt. Sin embargo, esta requiere de geometrías
existentes y por consiguiente no permite realizar un proceso inverso de diseño, por
lo tanto, la ecuación de Beale es la indicada para iniciar con el diseño ya que
establece el volumen de barrido por el pistón a partir de la potencia deseada, la
presión media y la frecuencia.

Se establecen los siguientes parámetros iniciales, considerando valores a partir de

la investigación realizada en el capítulo 4 “Estado del arte”:

• Potencia deseada: 14 W
• Presión media: 1,2 bar
• Frecuencia: 60 Hz

55
• Temperatura en la zona de
expansión deseada 575K

Con base al gráfico 1 y la temperatura deseada en la zona de expansión, se

determina el número adimensional de Beale que aplica para este tipo de motor el
cual es de: 0.0025, ya que el motor Stirling tiende a ser un motor pequeño, de baja
eficiencia.

A partir de la ecuación (11) se determina el volumen barrido por el pistón dando

como resultado:

14 𝑊𝑊
𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 = = 𝟕𝟕𝟕𝟕. 𝟒𝟒 𝒄𝒄𝒎𝒎𝟑𝟑
1.19 𝑏𝑏𝑏𝑏𝑏𝑏 ∗ 60 𝐻𝐻𝐻𝐻 ∗ 0.0025

7.1.2 Cálculo de diámetro de los pistones

Para el diseño se utilizó un cigüeñal del cabezal de un compresor de aire, debido a

que este elemento mecánico ya tiene sus dimensiones balanceadas correctamente.
Por lo tanto, con base a este se establece el recorrido barrido por los pistones dentro
de los cilindros para finalmente darle las dimensiones al motor.

La carrera del cigüeñal utilizado es de 4 cm por consiguiente el valor de los

diámetros son:

𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 × 4 78.5 𝑐𝑐𝑚𝑚3 × 4

𝑑𝑑 = � = � = 𝟒𝟒. 𝟗𝟗𝟗𝟗 𝒄𝒄𝒄𝒄 ~ 𝟓𝟓 𝒄𝒄𝒄𝒄
𝜋𝜋 × ℎ 𝜋𝜋 × 4 𝑐𝑐𝑐𝑐

7.1.3 Diseño de las camisas de los pistones

A continuación, en la tabla 6 se presenta la simbología de las cotas desarrolladas

en la figura 9 con sus respectivas medidas de diseño.

56
Tabla 6. Simbología de las dimensiones de los cilindros de las zonas de expansión
y compresión

Parámetros Símbolo Valor Unidades

Diámetro del cilindro de expansión 𝐷𝐷𝐷𝐷 50 mm

Diámetro del cilindro de
𝐷𝐷𝐷𝐷 50 mm
compresión
Diámetro del regenerador 𝐷𝐷𝐷𝐷 9.5 mm
Longitud total del cilindro de
𝐿𝐿𝐿𝐿𝐿𝐿 62 mm
expansión
Longitud total del cilindro de
𝐿𝐿𝐿𝐿𝐿𝐿 62 mm
compresión
Recorrido barrido por el pistón de
𝐿𝐿𝐿𝐿𝐿𝐿 40 mm
expansión
Recorrido barrido por el pistón de
𝐿𝐿𝐿𝐿𝐿𝐿 40 mm
compresión
Longitud inicial del pistón de
𝐿𝐿𝐿𝐿𝐿𝐿 12 mm
expansión
Longitud inicial del pistón de
𝐿𝐿𝐿𝐿𝐿𝐿 35 mm
compresión
Longitud de arco del regenerador 𝐿𝐿𝐿𝐿 190 mm

57
Figura 9. Parámetros de las dimensiones del motor Stirling solar tipo Alfa

58
Figura 10. Dimensiones del cilindro caliente (cotas en mm)

Figura 11. Dimensiones del cilindro frio (cotas en mm)

59
7.2 DISEÑO DEL COLECTOR SOLAR

Como se explica con anterioridad en el “Marco teórico” uno de los parámetros

importantes es la relación de concentración del colector, ya que determina la
capacidad de intensificar la radiación solar que capta el colector en un solo punto.
Para determinar este parámetro se utiliza la ecuación (501), la cual requiere que se
conozca el ángulo de borde (ϕr) y el ángulo de dispersión (δ).

De la tabla 3 se determina el ángulo de dispersión, donde se establecen los valores

de los errores típicos σ1d y σ2d. De tal modo:

𝜎𝜎1𝑑𝑑 = 6.1 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚

𝜎𝜎2𝑑𝑑 = 2.8 𝑚𝑚𝑟𝑟𝑟𝑟𝑟𝑟

A partir de la ecuación (45) se determina el σtotal:

1
𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎𝜎 = (6.12 + 2.82 )2 = 6.711 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚

Sustituyendo el valor obtenido en la ecuación (47), se obtiene con gran precisión el

efecto generado por la desviación angular (δ):

180°
𝛿𝛿 = ∗ 3 ∗ 6.711 ∗ 10−3 ≈ 1.15 °
𝜋𝜋

Una vez que se obtiene el valor de la desviación angular (δ) se determina el ángulo
de borde (ϕr) que debe de tener el plato, para ello se establecen diferentes valores
comprendidos entre cero y noventa grados (0°≤ϕr≤ 90°). De tal modo, a partir de la
ecuación (50) se busca cual es el valor que ofrece una mayor relación de
concentración.

60
Figura 12. Relación de concentración dependiendo el ángulo de borde (ϕr).

De la figura 12 se puede determinar que cuando el ángulo de borde (ϕr) es igual a

45° la relación de concentración es máxima, siendo:

𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 = 𝟐𝟐𝟐𝟐𝟐𝟐. 𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟗

Siendo el ángulo de borde de 45° el que ofrece dicha relación de concentración, se

establece como el ángulo de borde del colector, a partir de este se determinan los
diferentes valores del diámetro de apertura (Da) para analizar como varían los
parámetros de diseño como lo son: la distancia focal (f), área de apertura (Aa), radio
de borde (rr), altura del disco (y), diámetro mínimo del cilindro receptor (D) y el ancho
mínimo de la imagen focal (W).

Para ello primero se calcula el valor de la distancia focal utilizando la ecuación (39)
y a partir del valor que se obtiene, se sustituye en las ecuaciones (38), (48) y (49)
para determinar los valores de radio de borde (rr), diámetro mínimo del receptor (D)
y el ancho mínimo de la imagen focal (W).

61
Para determinar la altura del disco (y) se utiliza la ecuación (36), obteniendo la
ecuación (53):

𝑥𝑥 2
𝑦𝑦 =
4 ∗ 𝑓𝑓 (53)

Los valores de los parámetros anteriormente anunciados se pueden observar en la

tabla 7.

Tabla 7. Parámetros de diseño del colector según los diámetros de apertura (Da)

φr (°) Da (m) f (m) Aa (m2) rr (m) y (m) D (m) W (m)

45 0,66 0,39 0,34 0,46 0,070 0,013 0,019

45 1 0,60 0,79 0,70 0,104 0,020 0,029
45 1,2 0,72 1,13 0,84 0,125 0,025 0,035
45 1,5 0,90 1,77 1,05 0,156 0,030 0,044
45 1,8 1,08 2,54 1,26 0,188 0,037 0,053
45 2,1 1,26 3,46 1,47 0,219 0,043 0,062

Por medio de la tabla 7 se considera que el diámetro de apertura óptimo del colector
parabólico es de 0,66 m, el cual requiere de una distancia focal (f) de 0,39 m, radio
de borde (rr) de 0,46 m y una altura de disco de 0.07 m. Respecto a la imagen del
foco se determina que el diámetro mínimo (D) y ancho del foco (W) son de 0,013 y
0,019 m respectivamente.

Por otra parte, se puede observar en la tabla 7 que los parámetros de diseño son
directamente proporcionales al diámetro de apertura (Da), ya que al aumentar el
diámetro estos también aumentan. Sin embargo, en el caso de este trabajo no es
necesario un colector de gran tamaño al ser un motor de pequeño porte el cual tiene
como diámetro del cilindro del foco cerca de 0,05 m.

62
 8. DISEÑO DEL MOTOR STILING SOLAR

8.1 COMPONENTES DEL MOTOR STIRLING

A continuación, se especifican los componentes del Motor Stirling Solar con sus
respectivas descripciones de dimensiones y materiales, con el fin de dar a conocer
que consideraciones son importantes a la hora del diseño de un motor Stirling tipo
Alfa y que materiales son idóneos para la función que va a desempeñar cada parte.

• Cigüeñal:

Como anteriormente fue mencionado, el cigüeñal del motor se diseña con respecto
al cigüeñal de un cabezal de un compresor de aire debido a que este ya cuenta con
las dimensiones correctamente balanceadas.

Figura 13. Diseño cigüeñal en CAD

63
Figura 14. Dimensiones del cigüeñal (cotas en mm)

• Descripción:

Cigüeñal con dimensiones de 15cm de largo x 72 cm alto, aproximadamente.

Cuenta con dos levas que sirven para aumentar el momento de inercia. La carrera
que este cigüeñal genera es de 4 cm. Se establecen los cambios de diámetro para
contrarrestar las fuerzas axiales que generan los rodamientos.

• Material:

El cigüeñal del cabezal del compresor se usa como modelo ya que el material en el
cual está elaborado es de acero al carbono y para este motor Stirling se requieren
que sus partes tengan el menor peso posible. Por lo tanto, se establece que el
material idóneo para este sistema es duraluminio por medio del proceso de
fundición.

64
• Bielas:

Las bielas del Motor Stirling Solar se dimensionan con respecto a la carrera del
cigüeñal. Su estructura cuenta con dos partes para tener un ensamble simple y
están unidos por un tornillo de 5mm.

Figura 15. Diseño de las bielas en CAD

Figura 16. Secciones de las bielas

65
Figura 17. Dimensiones de las bielas (cotas en mm)

• Descripción:

Bielas de 10 cm de largo por 5 cm máximo de ancho, aproximadamente. Estas se

componen de dos partes ya que por facilidad de ensamble se busca que las dos
partes sean acopladas por un perno. El diseño cuenta con un vaciado para disminuir
el peso.

• Material:

Las bielas de los motores de combustión interna, por lo general, se fabrican en acero
templado, sin embargo, como se requiere que el Moto Stirling Solar sea ligero, se
establece que el material idóneo es el duraluminio. Este material cuenta con
importantes propiedades mecánicas y buena dureza. Con este material se fabrican
partes de maquinaria que requieren de gran precisión y de uso liviano. Adicional a
ello, este material permite que el proceso de elaboración sea sencillo.

• Pistones:

Los pistones del motor se diseñan con respecto a los pistones comerciales, sin
embargo, tienen en su estructura modificaciones para disminuir el peso.

66
Figura 18. Diseño de pistones en CAD

Figura 19. Dimensiones de los pistones (cotas en mm)

• Descripción:

Pistones de 5 cm de diámetro y 3 cm de alto. A diferencia de los pistones

comerciales, a estos se diseñan en el cielo (parte superior) una superficie esférica,
debido a que científicos que llevan a cabo estudios sobre los motores Stirling han
observado que al darle una redondez a la superficie del pistón que se encuentra en
67
contacto con el gas de trabajo, el funcionamiento del sistema es más óptimo. Esto
puede ser posible debido a que los gases en movimiento al tener contacto con
superficies esféricas tienden a tener mayor efectividad por la aerodinámica.

Los pistones no se diseñan con las dimensiones comerciales debido a que se busca
que tengan el menor peso posible, por lo que la parte de la falta se diseñan dos
semicírculos para eliminar ese material.

Adicional a ello, se diseñan dos ranuras en las cuales van ubicados los o’rings; estos
tendrán la función principal de confirmar el gas de trabajo y evitar fugas del
lubricante. Los o’rings a usar son los de referencia OM080.

• Material:

Los pistones deben ser diseñados de tal forma que permita tener una buena
propagación del calor y evitar altas tensiones moleculares para evitar daños en el
material. Generalmente los pistones están hechos en fundición de grano fino, pero
cuando se requieren que tengan dos partes, se usa el fondo de acero fundido para
resistir mejor las tensiones. Los materiales comunes son: hierro fundido, aleación
de níquel, aleación de aceros y aleación de aluminio.

Para este caso los pistones se realizan en duraluminio 7075 ya que es una aleación
de Zinc, Magnesio y Cobre que permite tener una buena maquinabilidad, alta
dureza, alto rango de fusión y buena conductividad térmica. 23

• Cilindro de la zona caliente:

El cilindro caliente es la camisa del pistón que se encuentra en la zona caliente, es

decir, en la zona de alta temperatura del Motor Stirling. El objetivo es transferirle la
mayor cantidad de calor al gas de trabajo proveniente del disco parabólico.

23
ALUMINIOS Y METALES UNICORNIO. Aluminio 7075 [en línea]. Aluminios y metales unicornio
[Consultado el 5 julio del 2018]. Disponible en:
http://www.aluminiosymetalesunicornio.com.mx/7075.html

68
Figura 20. Diseño del cilindro de la zona caliente en CAD

Figura 21. Dimensiones del cilindro de la zona caliente (cotas en mm)

• Descripción:

Camisa del pistón de la zona de alta temperatura de 6 cm de alto con un diámetro

de 5 cm, aproximadamente. La dimensión del diámetro se establece en el numeral
(7.1.3). En la parte superior se ubica la rosca del racor que une al cilindro con el
regenerador; el tipo de rosca es NTP 3/8 ya que esta cuenta con las normas de
hermeticidad.

69
• Material:

El material que se requiere para esta zona debe de tener propiedades de alta
absortividad y alta conductividad térmica, ya que, dependiendo de la capacidad del
material, así mismo es la cantidad de calor que se aprovecha proveniente del disco
parabólico. El paso posterior a la elección del material, es calcular la cantidad de
calor que se requiere, para ello se utiliza la ecuación (52).

Teniendo en cuenta que el Motor Stirling Solar está ubicado en la ciudad de

Santiago de Cali, Colombia, se establece según el IDEAM que la radiación promedio
en esta zona es entre 900 − 1000 𝑊𝑊/𝑚𝑚2 .

El material seleccionado en el cilindro de la zona caliente es el latón ya que, a

comparación de otros metales, este cuenta con un valor elevado de absortividad
(véase anexo B) y a su vez tiene un alto grado de conductividad térmica (véase
anexo C). Adicional a ello, el latón cuenta con otras propiedades que permiten que
su proceso de manufactura sea sencillo. A continuación, se especifican los
parámetros que se requieren para calcular el calor del receptor.

Tabla 8. Valores para determinar el calor del receptor (Qrecep).

Abreviatura Valor
Α 0,65
Ρ 0,97
Τ 1
F 0,9
Aa 0,31
I 900 W/m2
𝑸𝑸𝑸𝑸𝑸𝑸𝑸𝑸𝑸𝑸𝑸𝑸 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟖𝟖𝟖𝟖 𝑾𝑾

• Cilindro de la zona fría:

El cilindro frio es la camisa del pistón que se encuentra en la zona fría, es decir, de
la zona de baja temperatura del Motor Stirling Solar. Su objetivo es disipar el calor
con el cual llega a esta zona el gas, para generar un delta de temperatura ∆𝑇𝑇 mayor
y así mejorar su eficiencia. En este diseño se tiene en cuenta la cantidad de calor
suministrada por el disco parabólico, como se muestra en la tabla 8, para poder

70
calcular la cantidad de calor que se requiere disipar. El método por el cual se va a
disipar el calor en este Motor Stirling Solar es por medio de aletas. El número de
aletas también es calculado con respecto a la cantidad de calor a disipar.

Figura 22. Diseño del cilindro de la zona fría en CAD

71
Figura 23. Dimensiones del cilindro de la zona fría

• Descripción:

Camisa del pistón que se encuentra en la zona de baja temperatura con un diámetro
de 5 cm y una altura de 6 cm, aproximadamente. En la parte superior se encuentra
la rosca hermética la cual une el regenerador con el cilindro de la zona fría; el tipo
de rosca que se requiere es NPS 3/8. La dimensión del diámetro está determinada
en el numeral (7.1.3).

Teniendo en cuenta que en esta zona se requiere que la temperatura disminuya, se

diseñan aletas las cuales van a hacer el medio de disipación. Para su diseño se
requiere establecer el diámetro, el espesor y temperaturas de las fronteras (figura
24) para determinar la cantidad de aletas. A continuación, se especifica las
ecuaciones que determinan lo anteriormente descrito:

72
 Figura 24. Perfil de una aleta (cotas en metros)

𝑡𝑡
𝑟𝑟2𝑐𝑐 = 𝑟𝑟2 +
2 (54)

Longitudes corregidas 𝑡𝑡
𝐿𝐿𝐿𝐿 = 𝐿𝐿 +
2 (55)

𝐴𝐴𝐴𝐴 = 𝐿𝐿𝐿𝐿 ∗ 𝑡𝑡 (56)

Área de la aleta total 𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 = 𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 + 𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 (57)

𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 = 2 ∗ 𝜋𝜋 ∗ (𝑟𝑟2𝑐𝑐 2 − 𝑟𝑟12 ) (58)

𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 = 2 ∗ 𝜋𝜋 ∗ 𝑟𝑟2 ∗ 𝑡𝑡 (59)

Calor de la aleta 𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄 = 𝜂𝜂 ∗ 𝐴𝐴𝑎𝑎𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 ∗ ℎ ∗ (𝑇𝑇𝑎𝑎𝑎𝑎𝑎𝑎 − 𝑇𝑇𝑠𝑠 ) (60)

73
Área de la zona no
𝐴𝐴𝐴𝐴𝐴𝐴 = 2 ∗ 𝜋𝜋 ∗ 𝑟𝑟1 ∗ 𝑡𝑡 (61)
aleteada

Calor de la zona no
𝑄𝑄𝑄𝑄𝑄𝑄 = 𝐴𝐴𝐴𝐴𝐴𝐴 ∗ ℎ ∗ (𝑇𝑇𝑎𝑎𝑎𝑎𝑎𝑎 − 𝑇𝑇𝑠𝑠 ) (62)
aleteada

Calor total disipado

𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄 = (𝑄𝑄𝑄𝑄𝑄𝑄 + 𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄) ∗ 𝑛𝑛°𝑎𝑎 (63)
por las aletas

Eficiencia para 𝑄𝑄𝑄𝑄𝑄𝑄−𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄

𝜂𝜂 = 0,5 ∗ ( )
maquinas térmicas 𝑄𝑄𝑄𝑄𝑄𝑄 (64)

𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄
Número de aletas 𝑛𝑛°𝑎𝑎 =
𝑄𝑄𝑄𝑄𝑄𝑄 + 𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄 (65)

A partir del autor Walker Graham 24 se determina que la eficiencia del motor Stirling
es cerca del 40% y 50% de la eficiencia de ideal de Carnot como se establece en la
ecuación 64. Además, partiendo de la experiencia de autores como Martini 25 y
docentes de la universidad 26 en la elaboración de motores Stirling, han obtenido
eficiencias dentro de un rango del 15% al 30%. De tal modo, en este proyecto se
trabaja con una eficiencia del motor del 23% que corresponde a motores con
sistemas solares como fuente de energía.

Teniendo en cuenta el calor teórico brindado por el colector (tabla 8), la eficiencia
establecida en el párrafo anterior y la ecuación 64, se determina que el calor
necesario que se debe disipar es:

24
WALKER Graham. Stirling engine [en línea]. Lugar de publicación: Oxford University Press. Cap
5: Preliminary engine design, Thermal efficiency. New York, 1980. [Consultado: 10 de marzo de
2018]. Disponible en Internet: https://es.scribd.com/doc/45062130/Stirling-Engines-G-Walker-
Oxford-1980-WW

25 MARTINI, William R. Stirling engine design manual. [en linea] Segunda edición. EEUU: NASA,
1983. [Consultado: 05 de marzo del 2018]. Disponible en Internet:
https://ntrs.nasa.gov/search.jsp?R=19830022057

26
MENDOZA, Sebastián, VIDAL, Juan Ricardo, MELIÁN Cobas, VLADIMIR, Silva, ELECTO
Eduardo, ULIANOV LÓPEZ, Yuri. Proyección de un sistema Dish/Stirling para generación distribuida.
[en linea] En: El Hombre y la Máquina Julio-Diciembre. 2014, [consultado: 15 de marzo de septiembre
de 2018] Disponible en: http://4www.redalyc.org/articulo.oa?id=47838946008> ISSN 0121-0777

74
 0,23
𝑄𝑄𝑄𝑄𝑄𝑄𝑄𝑄 = � ∗ 158,8� − 158,8 = −𝟖𝟖𝟖𝟖, 𝟕𝟕𝟕𝟕𝟕𝟕
0,5

Por otro lado, partiendo de la temperatura de trabajo supuesta y la del ambiente,

300°C y 27°C respectivamente, se realiza un análisis de convección natural para
determinar el coeficiente de transferencia de calor (h), siendo de 9 W/°C, el cual se
utiliza en las ecuaciones 60 y 62.

A partir de la geometría del cilindro de la zona de compresión, propuesta en el anexo

E (planos), se determinan los valores de las longitudes corregidas (ecuaciones 54,
55 y 56) siendo los que se requieren dentro del gráfico (ver anexo D) para determinar
una eficiencia del 97%. Además, estos valores son utilizados para determinar el
área de la aleta por medio de las ecuaciones 57, 58 y 59, siendo de:

𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 = �2 ∗ 𝜋𝜋 ∗ (0,03952 − 0,02852 )� + (2 ∗ 𝜋𝜋 ∗ 0,0385 ∗ 0,002) =

𝟓𝟓, 𝟗𝟗𝟗𝟗 ∗ 𝟏𝟏𝟏𝟏−𝟑𝟑 𝒎𝒎𝟐𝟐

Una vez que se conocen los valores anteriores, se obtiene el número de aletas
mínimo que se requieren para disipar la cantidad de calor establecida anteriormente
a partir de la ecuación 65, obteniéndose un valor mínimo de 6 aletas, sin embargo,
para garantizar que el sistema disipe la cantidad deseada se establece un número
de 10 aletas.

• Material:

Como el principal objetivo del cilindro de la zona fría es disipar el calor, el material
idóneo para este requerimiento es el aluminio, ya que tiene una alta conductividad
térmica. En este caso se establece el duraluminio 7075 debido a que el proceso de
manufactura que se requiere para fabricarlo, necesita que el material no contenga
poros, lo cual una fundición de aluminio tiene probabilidades que queden poros en
su interior.

• Regenerador:

75
Esta es una parte muy importante del motor Stirling. El regenerador absorbe y
entrega calor al fluido de trabajo compensando una parte del calor perdido por el
motor, generando que la potencia y velocidad del motor se incrementen. Esto
sucede debido a que cuando trabaja el regenerador el gas necesitaría absorber
menor calor en cada ciclo, con lo cual hace que el ciclo necesite de menos tiempo
para realizarse y también se consuma menos cantidad de energía.

Figura 25. Diseño del regenerador en CAD

• Descripción:

El regenerador a utilizar en el Motor Stirling Solar se diseña con diámetro de 3/8’’ y

su longitud de arco es de aproximadamente 190mm. La efectividad de este
regenerador depende del tipo de material ya que al tener menor difusividad térmica
el material puede absorber con facilidad el calor.

• Material:

Teniendo en cuenta lo anterior, el material que mejor cumple con la especificación

es el acero inoxidable ya que tiene una difusividad térmica de 3.95x106 m2/s.

• Diseño de estructura:

76
La estructura del Motor Stirling Solar se diseña con respecto a las medidas del
cigüeñal, la carrera de los pistones y dimensiones de los cilindros. Esta consta de
dos partes, caja donde se encuentran las partes mecánicas y una tapa.

Figura 26. Diseño de la estructura del Motor Stirling Solar en CAD

Figura 27. Dimensiones de la estructura del Motor Stirling Solar

77
• Descripción:

La estructura es de 14 cm de ancho por 12 cm de alto y 9 cm de profundidad

aproximadamente. Cuenta con las perforaciones de los cilindros y del eje del
cigüeñal con las dimensiones a partir de la ubicación de las piezas.

• Material:

La estructura requiere que sea estable para dar firmeza al mecanismo, es por esto
que se diseña en láminas de acero inoxidable de 6 mm de espesor, debido a que el
motor estará en la intemperie y las propiedades del material no puede variar.

• Volante de inercia:

El volante de inercia es importante en todo sistema mecánico que tenga velocidad

angular ya que este permite almacenar energía y mejorar su eficiencia. Para el
diseño del volante del Motor Stirling Solar se consideran las siguientes ecuaciones:

Figura 28. Diseño del volante de inercia en CAD

78
Figura 29. Diámetro del volante de inercia

• Material:

El volante se diseña de aluminio 1060 debido a que es un material liviano y fácil de

mecanizar.

8.2 COMPONENTES DEL COLECTOR SOLAR

A partir de los parámetros de diseño que se obtienen del colector solar en el capítulo
anterior, se diseña el disco en el software SolidWorks como se puede observar en
la figura 30.

79
Figura 30. Croquis del diseño del colector solar con parámetros calculados

A partir del croquis de la figura 30, se obtuvo el diseño del colector (figura 31):

Figura 31. Diseño del colector solar en CAD

Teniendo en cuenta las dimensiones que se determinan del disco parabólico, se

trabaja con un disco de una parabólica de televisión que posee unas dimensiones
semejantes y conservan el ángulo de borde. Este tipo de discos no son totalmente
esféricos más bien ovalados, teniendo un diámetro mayor de 0,66 m y un diámetro
menor de 0,6 m, están elaborados de acero inoxidable teniendo en cuenta que son
para uso en intemperie.

80
Debido a la escasa reflectividad que tiene la superficie del disco se cubre con un
film Mylar con acabado de platino, que posee una reflectividad de cerca del 95 %.

A partir de la ecuación (41) se determina el área de un colector totalmente circular,

sin embargo, para este caso se determina el área teniendo la ecuación de una elipse
(ecuación 66), la cual tiene en cuenta ambos diámetros (diámetro mayor dmáx y
diámetro mínimo dmin).

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑
𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 = 𝜋𝜋 ∗ ∗
2 2 (66)

De tal modo sustituyendo el diámetro en la ecuación (66) el área del colector es:

𝐴𝐴𝐴𝐴 = 0,312 𝑚𝑚2

81
 9. CONSTRUCCIÓN DEL MOTOR STIRLING SOLAR

A continuación, se exponen los procesos de manufactura que se requirieron para

desarrollar la construcción tal como se ha planteado en el diseño.

9.1 FABRICACIÓN DE COMPONENTES DEL MOTOR STIRLING

• Cigüeñal:

Figura 32. Cigüeñal del mecanismo

El proceso de fabricación del cigüeñal se llevó a cabo por medio de una fundición
de duraluminio, en la cual fue necesario utilizar el cigüeñal que se usó en el diseño
como modelo, dejando un material en bruto que posteriormente se mecanizó y
rectificó dejándolo con dimensiones de planos para garantizar su balance.

El cambio de diámetro del eje se establece en función de los rodamientos

seleccionados para el mecanismo. (Ref. 61804)

82
Figura 33. Rodamientos seleccionados para el mecanismo

• Bielas:

Figura 34. Bielas del mecanismo

El proceso de manufactura de las bielas fue desarrollado por medio del software
MasterCam, debido a que se fabricaron en una fresadora CNC (Control Numérico
por Computadora) por su alta precisión.

La tolerancia radial que deben manejar estos sistemas debe ser aproximadamente
+-0.05mm con respecto al diámetro del cigüeñal y al diámetro del buje del pistón,
para así garantizar un movimiento estable y libre.

83
• Pistones:

Figura 35. Pistones del mecanismo

Los pistones se fabricaron por medio de un proceso de electro-erosionado, el cual

le da la forma final en su interior. En el exterior se utilizó mecanizado en torno para
dar las medidas precisas donde van ubicados los anillos.

Inicialmente se prueba el pistón dentro de la camisa para observar el movimiento

que debe ser libre y hermético con los o’rings seleccionados en el capítulo de
diseño, sin embargo, se observa que el material del cual está compuesto el o’ring
no permite generar un movimiento armónico, por lo tanto, se decide utilizar anillos
de teflón grafito ya que estos cuentan con auto lubricación que permiten generar el
movimiento libre del pistón y garantiza la hermeticidad.

• Cilindros:

Figura 36. Cilindros de la zona fría y caliente del mecanismo

84
Los cilindros se fabrican por medio del mecanizado en torno. El cilindro frio en su
parte superior va conectado un racor que tiene rosca NPT, la cual fue hecha en su
interior para el acople. Del mismo modo en la parte lateral del cilindro caliente se
realiza la rosca NPT. En el interior del cilindro frio se inserta un buje de latón debido
a que el coeficiente de rozamiento entre duraluminio y el latón es menor.

Figura 37. Buje de latón para el cilindro de la zona fría

• Regenerador:

Figura 38. Regenerador del mecanismo

Teniendo en cuenta las condiciones de diseño de acuerdo al material, se encontró

que las mangueras flexo metálicas que existen en la industria están compuestas de
acero inoxidable y en su interior constan de una superficie corrugada la cual es
idónea para generar turbulencia y aumentar la transferencia de calor. Adicional a
ello, se le inserta esponjilla metálica y balines de acero inoxidable para disminuir el
volumen muerto y aumentar la capacidad de almacenamiento del calor.

85
• Caja:

Figura 39. Caja del mecanismo

El proceso de manufactura de la caja del motor consta de dos partes, la primera es

el corte de las láminas de 6mm de espesor con láser de acuerdo con las
dimensiones establecidas, y la segunda es la soldadura de las partes mediante el
proceso TIG. En la parte interior se ubican cuatro (4) laminas pequeñas que
aseguran la tapa de la caja.

86
9.2 FABRICACIÓN DE COMPONENTES DEL DISCO PARABÓLICO

• Plato Solar:

Figura 40. Plato solar del motor Stirling

Como se estableció anteriormente en el capítulo de diseño, el plato solar se

construye utilizando una antena parabólica y en ella se ubica el papel reflectivo
(Space Paper Platinum Reflective Mylar Film). Para instalarlo primero se pule la
superficie con lija y se pega el papel con un aerosol PL 285, el cual es un pegante
especial que permite adherir a cualquier superficie de manera rápida y permanente.

• Base del motor:

La base del motor es de acero 1060 de 3mm de espesor el cual es diseñado para
darle estabilidad al motor por medio de prisioneros. Esta base está instalada en un
brazo que conecta al disco parabólico a la distancia requerida del foco.

87
9.3 ENSAMBLE DEL MOTOR STIRLING SOLAR

Figura 41. Motor Stirling Solar

Figura 42. Ensamble del motor Stirling con el plato solar

88
9.4 COSTOS DE ELABORACIÓN DEL MOTOR STIRLING SOLAR:

A continuación, se presentan los costos de fabricación del prototipo del Motor

Stirling Solar.

Tabla 9. Costos de fabricación

N° Descripción Cantidad Valor

1 Antena parabólica 1 $ 28.000
2 Volante (polea) 1 $ 12.000
3 Fundición duraluminio 1 $ 50.000
4 Racores 4 $ 43.000
5 Corte de laminas 11 $ 265.370
Fabricación de piezas:
(2) pistones
6 (1) cilindro latón interno 5 $ 1.118.600
(1) cilindro zona fría
(1) cilindro zona caliente
7 Maquinado de cigüeñal 1 $ 285.600
8 Orings y anillos 4 $ 106.000
9 Soldadura TIG 1 $ 350.000
10 Tornillería 20 $ 40.000
11 Empaques 2 $ 4.000
12 Codal 1 $ 10.000
13 Base motor 1 $ 120.000
14 Rodamientos 2 $ 26.000
15 Material de bielas 2 $ 70.000
16 Manguera flexo-metálica 1 $ 50.000
17 Hielo seco 2 $ 15.000
18 Manómetro 1 $ 15.000
TOTAL (COP) $ 2’608.500

Para el desarrollo de este proyecto se contó con el apoyo del Semillero de

Investigación de Sistemas Térmicos quien financió el 60% del proyecto.

89
 10. RESULTADOS Y VALIDACIÓN CON EL MODELO MATEMÁTICO

10.1 RESULTADOS

En esta etapa del proyecto se evalúan los datos obtenidos a partir de las pruebas
experimentales y se realiza la validación con el modelo matemático con base a las
temperaturas alcanzadas, para comparar las potencias y eficiencias obtenidas
teóricamente vs experimentalmente.

Pruebas realizadas:

• Temperatura del colector experimental

• Prueba en seco del Motor Stirling

• Prueba del funcionamiento del Motor Stirling con el Colector Solar

• Prueba del motor Stirling con diferentes fuentes de calor

La toma de datos de las temperaturas fue por medio de un sensor de temperatura

infrarrojo.

• Temperatura del colector experimental:

Se realizaron 4 pruebas en diferentes días, entre la 1 pm y 3 pm, en las cuales el

plato solar alcanzó una temperatura promedio de 88°C, como se muestra en la
Tabla 10.

Tabla 10. Valores de temperaturas registrados en cada prueba

Prueba Hora Temperatura

1 13:30 92°C
2 15:00 89°C
3 14:00 85°C
4 14:30 82°C

Con base a la temperatura promedio alcanzada del plato solar, se estima la cantidad
de calor experimental absorbido por el sistema teniendo en cuenta los valores del
calor específico (véase anexo C), el tiempo que tardo en estabilizarse la temperatura
y la masa de la superficie del material. Para ello se utiliza la ecuación 67:

90
 Cantidad de calor 𝑸𝑸𝑸𝑸𝑸𝑸𝑸𝑸 = 𝒎𝒎 ∗ 𝒄𝒄𝒄𝒄 ∗ 𝜟𝜟𝜟𝜟
experimental (67)

𝑄𝑄 = 0,04 ∗ 8522 ∗ (88 − 28) = 20452,8 𝐽𝐽

𝑡𝑡 = 180 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑒𝑒𝑒𝑒 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

20452,8
𝑸𝑸𝑸𝑸𝑸𝑸𝑸𝑸 = = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟔𝟔 𝑾𝑾
180

Cantidad de calor teórica 158,83 W

Cantidad de calor experimental 113,6 W

Los resultados obtenidos de la cantidad de calor teórico y experimental son

diferentes debido a que los cálculos realizados suponen un promedio de radiación
solar de 900 W/m2 donde requiere que el cielo este totalmente despejado, por lo
tanto, varia con respecto al experimental ya que en los días de prueba se contada
con nubosidad.

• Prueba en seco del Motor Stirling:

Se realiza una prueba en seco la cual consiste en mover el motor Stirling conectado
a un motor eléctrico, el cual le transmitirá potencia para realizar su movimiento a
varias revoluciones con el fin de observar las temperaturas que se alcanzan en los
cilindros. Para esta prueba las temperaturas deben se invertidas, es decir, el cilindro
caliente debe estar a menor temperatura y el frio a mayor temperatura, para así
observar si el motor se encuentra óptimo para la operación. Lo ideal en esta prueba
es que el cilindro de la zona de expansión alcance una temperatura tan baja que
permita mostrar escarcha.

91
Figura 43. Prueba en seco del motor Stirling con motor trifásico

En la prueba se conectó el eje del Motor Stirling Solar con el eje de un motor trifásico
que contaba con un variador de frecuencia para ir aumentando las revoluciones y
observar su comportamiento, como se muestra en la figura (43). A continuación, se
muestra en la tabla (11) los valores de revoluciones sometidas al motor Stirling y
sus respectivas temperaturas alcanzadas.

Tabla 11. Datos obtenidos de la prueba en seco

Revoluciones (rpm) Temp. zona fría Temp. zona caliente

800 32°C 29°C
900 36°C 30°C
1300 41°C 33°C
1500 43°C 36°C

Los valores registrados muestran un comportamiento aproximado de las

temperaturas en los cilindros. Sin embargo, se esperaba que la diferencia de
temperaturas fuera mayor.

En esta prueba se puede evidenciar que las piezas del motor Stirling y el mecanismo
cuentan con resistencia e idoneidad ya que se alcanzaron 1500 RPM y no se vieron
afectadas.

92
• Prueba del funcionamiento del Motor Stirling con el Colector Solar:

Se desarrollan las pruebas del motor Stirling con el colector solar en donde se
registran las siguientes temperaturas en los cilindros:

Cilindro de la zona de expansión: 85°C

Cilindro de la zona de compresión: 28°C

En las pruebas no se registra un movimiento del cigüeñal, sin embargo, al realizarse

manualmente un empuje al volante de inercia se siente el mecanismo con presión
interna, el cual valida que efectivamente el proceso de expansión del gas se está
realizando.

En la figura 44 se puede observar el motor Stirling con el Colector incidiendo la

radiación solar en la zona de expansión del gas de trabajo.

Figura 44. Prueba de temperatura de los cilindros con el plato solar

93
• Prueba del motor Stirling con diferentes fuentes de calor:

Al no evidenciar movimiento del cigüeñal con la temperatura alcanzada por el

colector solar, se procede a realizar pruebas con mechero y soplete para elevar la
temperatura en la zona de expansión, y en la zona de compresión se ubica hielo
seco para bajar la temperatura. Los valores de temperaturas alcanzadas con ellos
se pueden observar en la tabla 12.

Tabla 12. Valores de temperaturas en las zonas de expansión y compresión con

diferentes fuentes de calor

Fuente de calor Temperatura en la zona Temperatura en la zona

de expansión de compresión
Soplete 500 °C -
Mechero 180°C 12 °C

Realizadas estas pruebas no se evidencia movimiento por parte del cigüeñal del
motor Stirling. Se procedió a ubicar el motor eléctrico con variador de frecuencia y
realizar de nuevo las pruebas para darle un arranque inicial al motor Stirling. Durante
la prueba el motor daba indicios de arranque, pero no se efectuó.

Figura 45. Prueba del motor Stirling con soplete

94
Durante la prueba con el soplete, los anillos de teflón grafito ubicados en el pistón
de la zona caliente alcanzaron temperaturas elevadas los cuales están diseñados
para soportar máximo 260ºC, provocando que rayaran la superficie del cilindro y por
allí probablemente se generaran fugas del gas de trabajo.

En el proceso de poner a funcionar el motor Stirling, se comenzaron a detectar

comportamientos en las temperaturas de los cilindros no adecuados cuando se
presentaban las siguientes condiciones:

• Ubicar el mechero y quitar el hielo seco

• Ubicar el hielo seco y quitar el mechero

Al realizar estas acciones los cilindros que no tenían la fuente externa tendían a
adquirir la temperatura de la zona que si la tenía (ya sea fría o caliente), por lo tanto
se pudo considerar que el regenerador no estaba realizando la función
correspondiente ya que, el principal objetivo de este es poder almacenar la mayor
cantidad posible de calor que lleva el gas de trabajo de la zona de expansión a la
zona de compresión y al observar que las temperaturas pasaban al otro cilindro, el
regenerador no tenía la suficiente capacidad para almacenar el calor transferido.

10.2 VALIDACIÒN CON EL MODELO MATEMÀTICO

Para validar el funcionamiento del motor Stirling se cuentan con diferentes modelos
termodinámicos que permiten estimar la potencia y eficiencia que se obtendrá
teóricamente a partir de unas dimensiones estipuladas. El modelo termodinámico
se basa en el análisis que realizó el ingeniero Israel Urielli 27 el cual relaciona las
ecuaciones y parámetros iniciales basado en las ecuaciones del modelo de Schmidt.
Los parámetros iniciales de operación se describen en la tabla 13.

27
URIELI Israel; Berchowitz David M. Stirling Cycle Machine Analysis [en linea]. Bristol: A. Hilger,
1984 [Consultado: 10 de abril de 2018]. Disponible en Internet:
https://www.ohio.edu/mechanical/stirling/isothermal/Schm_summary.html

95
Tabla 13. Parámetros de operación iniciales para el desarrollo del modelo
termodinámico del motor Stirling solar

Parámetros Símbolo Valor Unidades

Presión inicial 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 101.3 KPa

Gas de trabajo − Aire -

Temperatura foco caliente 𝑇𝑇ℎ 358 K

Temperatura foco frio 𝑇𝑇𝑇𝑇 301 K

𝐽𝐽
Constante de gas ideal 𝑅𝑅 287
𝐾𝐾𝐾𝐾. 𝐾𝐾
Ángulo del ciclo 𝑥𝑥 0-360 °

Fase 𝑑𝑑𝑑𝑑 90 °

Revoluciones del motor 𝑛𝑛 900 Rpm

Los valores de las temperaturas utilizadas para el modelo matemático son

equivalentes a las temperaturas experimentales alcanzadas en los cilindros de
expansión y de compresión en la prueba del funcionamiento del Motor Stirling con
el Colector.

Inicialmente el modelo de Schmidt requiere introducirle datos de: parámetros de

operación iniciales y las magnitudes de volúmenes que son hallados a partir del
diseño del motor Stirling. En la tabla 14 se expresan los volúmenes con sus
respectivas magnitudes y unidades.

Estos datos se ingresan al software Matlab con el fin de realizar la iteración de forma
rápida y precisa.

96
Tabla 14. Magnitud de volúmenes con sus respectivas unidades

Nombre Símbolo Magnitud Unidades

Volumen de amplitud de
𝑉𝑉𝑉𝑉𝑉𝑉 7.85𝑥𝑥10−5 𝑚𝑚3
compresión

Volumen de amplitud de
𝑉𝑉𝑉𝑉𝑉𝑉 7.85𝑥𝑥10−5 𝑚𝑚3
expansión

Volumen total de
𝑉𝑉𝑉𝑉𝑉𝑉 1.02𝑥𝑥10−4 𝑚𝑚3
compresión

Volumen total de expansión 𝑉𝑉𝑉𝑉𝑉𝑉 1.04𝑥𝑥10−4 𝑚𝑚3

Volumen barrido
𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 1.57𝑥𝑥10−4 𝑚𝑚3
compresión
Volumen barrido expansión 𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 1.57𝑥𝑥10−4 𝑚𝑚3

Volumen inicial de
𝑉𝑉𝑉𝑉 6.87𝑥𝑥10−5 𝑚𝑚3
compresión
Volumen inicial de
𝑉𝑉ℎ 2.35𝑥𝑥10−5 𝑚𝑚3
expansión
Volumen muerto de
𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 2.35𝑥𝑥10−5 𝑚𝑚3
compresión (Vco - Vca)
Volumen muerto de
𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 2.5510−5 𝑚𝑚3
expansión (Veo - Vea)

Volumen del regenerador 𝑉𝑉𝑉𝑉 1.346𝑥𝑥10−5 𝑚𝑚3

Las ecuaciones que se utilizan para el desarrollo de la teoría de Schmidt se

encuentran desde la ecuación 12 hasta ecuación 35, descritas en el numeral
(5.1.6.1).

• Volúmenes actuales de compresión (𝑉𝑉𝑉𝑉) y expansión (𝑉𝑉𝑉𝑉):

Los valores de estos volúmenes varían con respecto a la posición del cigüeñal (0º -
360°), por lo tanto, se obtienen 361 resultados los cuales varían entre:
97
 [0.0236E-3 m3 - 0.1806E-3 m3].
• Volumen total:

El volumen total como se explica anteriormente en el numeral (5.1.6) es la sumatoria

de los volúmenes actuales de compresión, los volúmenes actuales de expansión y
el volumen del regenerador.

• Relación de volúmenes y temperaturas:

𝑡𝑡 = 𝟎𝟎. 𝟖𝟖𝟖𝟖𝟖𝟖𝟖𝟖

𝑣𝑣 = 𝟏𝟏

𝑋𝑋𝑋𝑋𝑋𝑋 = 𝟎𝟎. 𝟑𝟑

𝑋𝑋𝑋𝑋𝑋𝑋 = 𝟎𝟎. 𝟑𝟑

𝑋𝑋𝑋𝑋 = 𝟎𝟎. 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏

• Coeficientes:

Las operaciones matemáticas que se realizan para encontrar la ecuación que

determina la potencia y la presión del ciclo son extensas, lo cual se simplifican
determinando los siguientes coeficientes:

𝑎𝑎 = 𝟒𝟒𝟒𝟒. 𝟗𝟗𝟗𝟗 °

𝑆𝑆 = 𝟑𝟑. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐

𝐵𝐵 = 𝟏𝟏. 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑

• Masa Total:

98
 𝑐𝑐 = 𝟎𝟎. 𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒

Para determinar la masa se asume la presión como a la que el sistema va a estar

inicialmente, en este caso, a presión atmosférica. Adicional a ello, se requiere
asumir que la masa total del gas en el sistema es constante.

𝑀𝑀 = 𝟏𝟏. 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏−𝟒𝟒 𝑲𝑲𝑲𝑲

• Presión:

Para el cálculo de la presión se obtienen diferentes valores ya que depende de la

variación del ángulo que recorre el cigüeñal. Sin embargo, la presión media se
obtiene a partir de la ecuación (30) y establece la presión promedio del ciclo Stirling.

1
𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 = ∮ 𝑃𝑃𝑃𝑃𝑃𝑃 = 𝟏𝟏𝟏𝟏𝟏𝟏 𝑲𝑲𝑲𝑲𝑲𝑲
2𝜋𝜋

• Trabajo en la zona de expansión y compresión:

Teniendo en cuenta que el modelo isotérmico para los espacios de compresión y

expansión se tiene que la temperatura es un valor constante, al igual que un flujo
másico constante. La ecuación de energía se simplifica a:

𝑑𝑑𝑑𝑑 = 𝑑𝑑𝑑𝑑 − 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅

𝑄𝑄𝑄𝑄 = 𝑊𝑊𝑊𝑊
𝑄𝑄𝑄𝑄 = 𝑊𝑊𝑊𝑊

Por lo tanto, las ecuaciones (31) y (32) establecidas con los parámetros de Schmidt
dan como resultado el trabajo realizado en la zona de compresión como en la zona
de expansión, siendo estos equivalentes al calor requerido para la expansión y el
calor que se debe disipar.

99
 𝑊𝑊𝑊𝑊 = 𝟔𝟔. 𝟕𝟕𝟕𝟕𝟕𝟕𝟕𝟕 𝑱𝑱

𝑊𝑊𝑊𝑊 = −𝟓𝟓. 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔 𝑱𝑱

Se obtiene un trabajo total del ciclo:

𝑊𝑊𝑊𝑊 = 5.2978 + (−2.9864) = 𝟏𝟏. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 𝑱𝑱

• Potencia y eficiencia del motor Stirling solar tipo Alfa:

Para estimar la potencia se asume que el motor alcanza 900 RPM, sin embargo,
este valor se verifica en las pruebas que se le realizan al equipo en funcionamiento
para determinar las revoluciones que el motor puede mantener constante, ya que al
ser un dispositivo que depende de la energía solar tiene fluctuaciones en sus
revoluciones.

𝐿𝐿𝐿𝐿 = 𝟏𝟏𝟏𝟏. 𝟎𝟎𝟎𝟎 𝑾𝑾

𝑒𝑒 = 𝟎𝟎. 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 = 𝟏𝟏𝟏𝟏. 𝟗𝟗%

El valor de la eficacia depende de los resultados obtenidos del trabajo del sistema,
es por esto que puede variar dependiendo de la temperatura alcanzada en el foco
y la cantidad de calor que el cilindro frio alcance a disipar.

En la tabla 15 se relacionan las magnitudes obtenidas de presión y volumen para

cada ángulo del cigüeñal (0° a 360°) las cuales forman el diagrama P-V del ciclo
Stirling

100
Tabla 15. Valores de Presión y Volumen en función del ángulo del cigüeñal

Ángulo del cigüeñal

Presión (KPa) Volumen (cm3)
(grados)
0 198.6 139.1

45 250.3 106.6

90 213.9 139.1

135 149.8 217.7

180 113.3 296.2

225 101.4 328.7

270 108.9 296.2

315 139.2 217.7

360 198.6 139.1

Como se puede observar en la tabla 15, la presión máxima alcanzada es cuando el

cigüeñal se encuentra aproximadamente a 45° y el volumen máximo cuando se
encuentra a 225°.

101
Figura 46. Presión vs Volumen del motor Stirling solar tipo alfa

En la figura 46 se puede evidenciar que el análisis que se realiza por el modelo de

Schmidt arroja una presión máxima aproximadamente de 250 KPa, partiendo de la
atmosférica que es la presión a la que está el gas inicialmente. Una presión media
de 160 KPa aproximadamente.

102
Figura 47. Variación del volumen de expansión vs el ángulo del cigüeñal

Figura 48. Variación del volumen de expansión vs el ángulo del cigüeñal

103
Figura 49. Variación de la presión vs el ángulo del cigüeñal

Figura 50. Variación de potencia con respecto al diámetro de diseño

104
Según el modelo de Schmidt, el motor Stirling solar tipo Alfa diseñado alcanzaría
una potencia aproximadamente de 16 W.

Como se puede observar en la Figura 50 la variación de la potencia neta que puede

entregar el motor Stirling depende del diámetro de los cilindros. Si el diámetro
incrementa, aumenta el volumen en las zonas de expansión y compresión, el motor
entregaría mayor potencia. Sin embargo, los motores Stirling al requerir una fuente
de calor externa, el volumen de gas a calentar no puede excederse en tamaño,
debido a que se necesitaría que la fuente de calor entregara mayor calor. Para este
caso en donde la fuente de calor es proveniente de un disco parabólico es difícil
diseñar un motor Stirling con diámetros de 10 o 15 cm ya que se necesitaría un
disco con diámetro mínimo de 8 a 10m. Estos sistemas en los cuales se busca tener
una potencia elevada son los que actualmente se encuentran instalados en parques
solares de motores Stirling y cuentan con discos de 11 m.

Otro factor que influye en la potencia entregada del motor Stirling es el ∆𝑇𝑇, es decir,
la diferencia entre las zonas de alta y de baja temperatura ya que al tener una
diferencia de temperaturas considerable se entiende que el proceso del
intercambiador de calor frio (aletas), el proceso de expansión del gas y el
regenerador están siendo efectivos. Como se puede observar en la siguiente tabla,
los valores de eficiencia del ciclo aumentan al tener un mayor delta de temperatura.

Tabla 16. Delta de temperatura vs eficiencia del ciclo Stirling

∆𝑇𝑇 Eficiencia
10 0,03
50 0,134
100 0,2364
150 0,3171
200 0,3824
250 0,4363
300 0,4815

105
• Comparativo de valores teóricos vs experimentales:

Teniendo en cuenta que el motor Stirling no efectuó ningún movimiento, no es

posible realizar un comparativo con los valores de potencia y eficiencia, sin
embargo, es posible comparar otros datos relevantes que se utilizaron en el diseño
y si tienen una variable experimental para ser contrarrestada como:

• Calor del receptor teórica vs calor del receptor experimental:

Tabla 17. Calor captado por el receptor teórica vs experimental

Cantidad de calor teórica 158,83 W

Cantidad de calor experimental 113,6 W

• Potencia deseada vs potencia teórica:

Se toma la potencia de diseño que se establece en la ecuación de Beale y la

potencia hallada por la teoría de Schmidt.

Tabla 18. Potencia alcanzada teóricamente por medio de Beale y Schmidt

Potencia de diseño (300ºC) 14 W

Potencia Schmidt (88ºC) 16 W
Potencia Schmidt (300ºC) 70 W

• Eficiencias con valores de temperaturas alcanzadas:

Con el modelo matemático se halla la eficiencia que se obtendría en el motor Stirling

a partir de las temperaturas alcanzadas con las diferentes fuentes de calor.

Tabla 19. Eficiencia teórica con las temperaturas alcanzadas

Temperaturas
Fuentes de
Calor Zona
Zona fría Eficiencia
Caliente
Plato Solar 28°C 88°C 15,9%
Mechero 12°C 150°C 32,6%
106
 11. CONCLUSIONES

Los valores hallados de la cantidad de calor del receptor teóricos y experimentales

tuvieron un error del 40% debido a que el parámetro de radiación con el que se
diseño fue en condiciones atmosféricas ideales, y las condiciones en el momento
de la prueba no fueron las óptimas.

Es importante observar que la potencia de diseño utilizada en la ecuación de Beale

y la potencia obtenida por el modelo de Schmidt difieren en gran medida, ya que en
la teoría de Schmidt se considera la temperatura de la zona de compresión, es decir
existe un delta de temperatura que va directamente relacionado con la eficiencia
que se obtendrá, mientras que en la teoría de Beale solo se parte de la temperatura
en la zona de expansión.

El diseño de aletas depende en gran magnitud del gradiente de temperatura que

exista dentro del cilindro y la temperatura ambiente, de tal modo al no alcanzar a la
temperatura supuesta la efectividad de las aletas se ve directamente afectada y por
consiguiente el calor disipado será menor.

A partir el estado del arte se pudo evidenciar que por medio de un colector
parabólico con las mismas dimensiones se pueden alcanzar temperaturas
aproximadas de 300ºC sin embargo en este proyecto se alcanzó un máximo de
90ºC lo cual indica que el papel Mylar utilizado no tiene una reflectividad del 0,97
como indica el fabricante. Por lo tanto, se requiere de un proceso de manufactura
especial para construir el plato donde se utilicen materiales de alta reflectividad.

Las pruebas en seco permitieron evidenciar que tanto el mecanismo de biela

manivela como los materiales seleccionados fueron adecuados, ya que al ser
sometido el sistema a 1500 RPM no se tuvo problemas con daños en el mecanismo.
Se realizaron varias pruebas al Motor Stirling Solar con el fin de observar su
comportamiento en donde no se logró efectuar un movimiento en el motor. Es
posible considerar que las causas de ellos hayan sido producto del regenerador, al
observar el comportamiento de las temperaturas como se expresó en el capítulo 9.

También se puede considerar que el aumento de temperatura con el soplete en la

zona de expansión haya generado fugas del gas de trabajo, ya que al tener rayones
en la zona se procedió a rectificar el cilindro interno lo cual modifico su medida.

107
A pesar no haber logrado un resultado de generación en el prototipo, este proyecto
da como primera instancia para continuar con la investigación del motor Stirling, ya
que si bien es un motor que puede llegar a revolucionar la generación de energía y
ser una posible solución a las zonas no interconectadas del país. Pero como todo
en la ciencia se requiere de más estudios para llevar a cabo un diseño eficiente es
por esto que este proyecto estará como referencia en el semillero de investigación
de sistemas térmicos de la Universidad Autónoma de Occidente para futuras
investigaciones.

108
 12. RECOMENDACIONES

• Para tener un óptimo tiempo de calentamiento desde el plato solar, se

recomienda realizar un sistema de seguimiento solar por medio de sensores, así se
podrá garantizar que el plato alimente continuamente el sistema.

• Los anillos que se utilizan en los pistones, deben ser fabricados con gran
exactitud para evitar el rozamiento con la camisa sin perder la hermeticidad y no
llevar el sistema a mayor de 260°C.

• Utilizar otro material para realizar el regenerador o diseñarlo de tal modo que
garantice la conservación de energía y no permita el calentamiento en el cilindro frio
ni el enfriamiento del cilindro caliente.

• No realizar las perforaciones de las placas de la base del cigüeñal por separado,
realizarlas de forma simultánea para garantizar la concentricidad. Esto puede
generar que los rodamientos se encuentren sobre ajustados y no permitan un libre
movimiento.

• Para alcanzar un valor de temperatura mayor proveniente del plato solar, es

necesario diseñarlo con un diámetro mayor para captar mayor radiación solar.

• Es recomendable ubicar en el cilindro caliente un cono con mayor diámetro para

poder captar la mayor radiación proveniente del foco, ya que esta se puede llegar a
propagar de manera difusa.

• En el cilindro frio es idóneo diseñarle un sistema de intercambiador de calor

adicional a las aletas, ya que ellas pueden ayudar a disipar calor, pero es
recomendable disminuir la temperatura por debajo de la temperatura ambiente.

109
 BIBLIOGRAFIA

AGREDA Víctor, MARASCA Tony. Diseño y fabricación a escala de un motor Stirling

accionado por energía solar [en línea]. Trabajo de investigación Ingeniero Mecánico.
Barbula, Venezuela.Universidad de Carabobo. Facultad de Inegniería. Escuela de
Ingeniería Mecánica. Mayo, 2013. [Consultado en 15 de abril de 2018]. Disponible
en internet: https://es.scribd.com/document/207121829/DISENO-Y-
FABRICACION-A-ESCALA-DE-UN-MOTOR-STIRLING-ACCIONADO-POR-
ENERGIA-SOLAR

AGÜERO ZAMORA, Víctor Raúl. Diseño y construcción de un motor Stirling para la

generación de energía eléctrica [en línea]. Proyecto de grado Ingeniero Mecánico.
Lima Perú. Universidad Nacional de Ingeniería. Facultad de Ingeniería Mecánica.
2011. p. 42. Disponible en Internet: http://cybertesis.uni.edu.pe/handle/uni/756

AHUMADA Omar. Recomiendan al Gobierno racionamiento programado de energía

del 5 %. 2016 [en línea]. En: El Tiempo, marzo 6 2016. [Consultado: 5 de mayo del
2018]. Disponible en http://www.eltiempo.com/archivo/documento/CMS-16531353

ALUMINIOS Y METALES UNICORNIO. Aluminio 7075 [en línea]. Aluminios y

metales unicornio [Consultado el 5 julio del 2018]. Disponible en internet:
http://www.aluminiosymetalesunicornio.com.mx/7075.html

BECKMAN, William A.; DUFFIR, John A. Solar Engineering of Thermal Processes.

4 ed. New Jersey: John Wiley & Sons, Inc, 2013. 359 p. ISBN 978-0-470-87366-3.

DURÁN MEDINA, Pino; BARBOSA SALDAÑA, Juan Gabriel. Diseño y construcción

de un prototipo de concentrador solar parabólico de disco reflector para generación
de energía térmica [En línea]. Instituto Politécnico Nacional. Escuela superior de
Ingeniería Mecánica y Eléctrica, 2012. 54 p. [Consultado: 15 de mayo del 2018].
Disponible en http://tesis.ipn.mx/jspui/handle/123456789/17648

Energías renovables, con gran potencial para crecer en el país [en línea]. En: El
Tiempo: abril 2018. [Consultado: 20 de mayo del 2018]. Disponible en Internet:
https://www.eltiempo.com/vida/medio-ambiente/las-energias-renovables-tienen-un-
gran-potencial-para-crecer-en-colombia-208116

110
GARCÍA José Luis. Diseño de un sistema prototipo de conversión de energía solar
en potencia mecánica compuesto por un concentrador solar y motor Stirling [en
línea]. Proyecto de grado Ingeniero Mecánico. Medellín. Universidad EAFIT.
Facultad de Ingeniería. Departamento de Ingeniería, 2010. [Consultado el 10 de
abril de 2018]. Disponible en internet:
https://repository.eafit.edu.co/handle/10784/4273

HERNANDEZ Miguel, SANCHEZ José Salvador. Diseño y prototipo de un motor

Stirling. [En línea]. Trabajo de investigación Ingeniero Mecánico. Universidad
Nacional Autónoma de México. Facultad de Ingeniería. 2013 [Consultado el 20 de
marzo de 2018]. Disponible en internet:
http://132.248.52.100:8080/xmlui/handle/132.248.52.100/4961

INSTITUTO DE PLANIFICACIÓN Y PROMOCIÓN DE SOLUCIONES

ENERGÉTICAS PARA LAS ZONAS NO INTERCONECTADAS (IPSE). Informe de
gestión 2017 [En línea]. p. 5. [Consultado: 20 de junio del 2018]. Disponible en
Internet: http://www.ipse.gov.co/ipse/informes-de-gestion

MARTINI, William R. Stirling engine design manual. [En linea] Segunda edición.
EEUU: NASA, 1983. [Consultado: 05 de marzo del 2018]. Disponible en Internet:
https://ntrs.nasa.gov/search.jsp?R=19830022057

MENDOZA, Sebastián, VIDAL, e tal. Proyección de un sistema Dish/Stirling para

generación distribuida. [en linea]En: El Hombre y la Máquina 2014, (Julio-
diciembre): [Fecha de consulta: 15 de marzo de septiembre de 2018] Disponible
en:<http://4www.redalyc.org/articulo.oa?id=47838946008> ISSN 0121-0777

PASCUAL TOLEDANO, Juan Ignacio. Configuraciones Alfa, Beta y Gamma

[Imagen]. Simulación de un motor Stirling mediante modelos de primer orden [en
línea]. Trabajo de investigación para Ingeniero Técnico Mecánico. España.
Universidad Carlos III de Madrid. Departamento de ingeniería térmica y fluidos,
2012. p. 24. [Consultado: 10 de febrero del 2018]. Disponible en Internet: https://e-
archivo.uc3m.es/handle/10016/18099#preview

SAAVEDRA, Mónica Andrea. Generación eléctrica en Colombia está holgada ante

la demanda [en línea]. EL MUNDO. Agosto, 2017. [Consultado el 15 de mayo de
2018] Disponible en Internet: http://www.elmundo.com/noticia/Generacion-
electrica-en-Colombia-esta-holgada-ante-la-demanda/357381

111
SÁENZ, Jorge. En penumbras: 2,5 millones de colombianos no tienen energía
eléctrica [En línea]. En: El Espectador. Agosto. 2017. [Consultado: 19 de abril del
2018]. Disponible en Internet: https://www.elespectador.com/economia/en-
penumbras-25-millones-de-colombianos-no-tienen-energia-electrica-articulo-
706892

STINE, W. Solar energy fundamentals and design with computer applications. Wiley,
Estados Unidos. 1985. Citado por DURÁN MEDINA, Pino; BARBOSA SALDAÑA,
Juan Gabriel. Diseño y construcción de un prototipo de concentrador solar
parabólico de disco reflector para generación de energía térmica [En línea].
Obtención del grado de Maestro en ciencias en Ingeniería Mecánica. [Consultado:
19 de abril del 2018]. Disponible en internet:
http://tesis.ipn.mx/jspui/handle/123456789/17648

RODRÍGUEZ José Arturo. Diseño y construcción de un motor Stirling tipo beta

acoplado a un concentrador solar parabólico cóncavo. [En línea]. Trabajo de
investigación Ingeniero Mecánico. Universidad de San Carlos de Guatemala.
Facultad de Ingeniería. Escuela de Ingeniería mecánica, 2015. [Consultado el 5 de
abril de 2018]. Disponible en internet: http://www.repositorio.usac.edu.gt/3173

UNIVERSIDAD NACIONAL DE COLOMBIA. Microrred generaría energía a bajo

costo en zonas no interconectadas [En línea]. Agenciadenoticias.unal.edu.co.
Medellín, 18 de abril de 2017. [Consultado: 20 de junio del 2018]. Disponible en
Internet: http://agenciadenoticias.unal.edu.co/detalle/article/microrred-generaria-
energia-a-bajo-costo-en-zonas-no-interconectadas.html

UNIDAD DE PLANEACIÓN MINERO ENERGÉTICA. Informe de rendición de

cuentas UPME 2016 2017 [en línea]. Unidad de Planeación Minero Energética:
Octubre, 2017, p.29. [Consultado: 15 de mayo de 2018]. Disponible en Internet:
http://www1.upme.gov.co/InformesGestion/Memorias_Rendicion_2017.pdf. ISSN
2463-123X

UPME. Energía eléctrica SIN (Sistema Interconectado Nacional) [En línea]. Unidad
de Planeación Minero Energética. [Consultado: 20 de junio del 2018]. Disponible en
Internet: http://www1.upme.gov.co/InformacionCifras/Paginas/PETROLEO.aspx

112
URIELI Israel; Berchowitz David M. Stirling Cycle Machine Analysis [en linea].
Bristol: A. Hilger, 1984 [Consultado: 10 de abril de 2018]. Disponible en Internet:
https://www.ohio.edu/mechanical/stirling/isothermal/Schm_summary.html

WALKER Graham. Stirling engine [en línea]. Lugar de publicación: New York, 1980.
Oxford University Press. Cap 5: Preliminary engine design, Thermal efficiency.
[Consultado: 10 de marzo de 2018]. Disponible en Internet:
https://es.scribd.com/doc/45062130/Stirling-Engines-G-Walker-Oxford-1980-WW

113
 ANEXOS

Anexo A. Tabla de emisividad normal de los metales

114
Anexo B. Tabla de absortividad en superficies de los metales

115
Anexo C. Tabla de propiedades térmicas de algunos elementos metálicos

116
Anexo D. Eficiencia de aletas circulares de espesor constante

117
Anexo E. Planos

PLANOS DEL SISTEMA MOTOR STIRLING Y COLECTOR

118
120
122
123
125
126
128
129
130
135