análisis termodinámico-mecánico de un prototipo de motor stirling

ANÁLISIS TERMODINÁMICO-MECÁNICO DE UN PROTOTIPO DE MOTOR
STIRLING DE CONFIGURACIÓN GAMMA DE BAJA POTENCIA
MARIO SUÁREZ LÓPEZ
UNIVERSIDAD AUTÓNOMA DE OCCIDENTE
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE ENERGÉTICA Y MECÁNICA
PROGRAMA DE INGENIERÍA MECANICA
SANTIAGO DE CALI
2015
ANÁLISIS TERMODINÁMICO-MECÁNICO DE UN PROTOTIPO DE MOTOR
STIRLING DE CONFIGURACIÓN GAMMA DE BAJA POTENCIA
Trabajo de grado para optar por el título de
Ingeniero Mecánico
Director
ELVER MAURICIO BARRERA
Ingeniero Mecánico
UNIVERSIDAD AUTÓNOMA DE OCCIDENTE
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE ENERGÉTICA Y MECÁNICA
PROGRAMA DE INGENIERÍA MECANICA
SANTIAGO DE CALI
2015
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
DUCARDO LEÓN MOLINA
JURADO
JUAN RICARDO VIDAL
JURADO
Santiago de Cali, Junio 13 de 2013
3
MIS MÁS SINCEROS AGRADECIMIENTOS A:
Dios que ha sido fuente inagotable de amor, fortaleza y sabiduría.
Mi madre por su incansable e incondicional colaboración y comprensión.
Mi Padre y hermanos por su valioso apoyo en todo momento.
Esa amiga por su sincera amistad y ayuda oportuna.
4
CONTENIDO
pág.
RESUMEN ............................................................................................................. 13
INTRODUCCIÓN ................................................................................................... 14
1. ANTECEDENTES ........................................................................................... 15
1.1
DESCRIPCION DEL PROTOTIPO EXISTENTE ......................................... 16
2. JUSTIFICACIÓN ............................................................................................. 18
3. OBJETIVOS ..................................................................................................... 19
3.1
OBJETIVO GENERAL ................................................................................. 19
3.2
OBJETIVOS ESPECÍFICOS ........................................................................ 19
4. MARCO TEÓRICO .......................................................................................... 20
4.1
CICLO STIRLING ........................................................................................ 20
4.2
MOTOR STIRLING TIPO ALFA .................................................................. 24
4.3
MOTOR STIRLING TIPO BETA .................................................................. 25
4.4
MOTOR STIRLING TIPO GAMMA .............................................................. 26
4.5
MECANISMO BIELA MANIVELA ................................................................ 26
4.5.1
Estudio cinemático .............................................................................. 27
5. METODOLOGÍA .............................................................................................. 29
6. FORMULACIÓN DEL ANÁLISIS MECÁNICO ................................................ 32
5
6.1
PARÁMETROS CINEMÁTICOS .................................................................. 32
6.2
CINEMÁTICA DEL MECANISMO ............................................................... 34
6.3
ANGULO DE TRANSMISIÓN ...................................................................... 37
6.4
MASAS EQUIVALENTES............................................................................ 39
6.4.1
Modelo dinámicamente equivalente biela. ......................................... 39
6.4.2
Modelo estáticamente equivalente de manivela................................ 40
6.5 FUERZAS DEINERCIA Y DE SACUDIMIENTO .............................................. 43
7. FORMULACIÓN DEL ANALISIS TERMODINÁMICO ..................................... 46
7.1
VARIACIÓN DE VOLÚMENES DE COMPRESIÓN Y EXPANSIÓN ........... 48
7.2
EXPRESIÓN PARA LA PRESIÓN .............................................................. 53
8. PARES DE TORSIÓN ...................................................................................... 55
8.1
PAR DE TORSIÓN DEL GAS ...................................................................... 55
8.2
PARES DE TORSIÓN DE INERCIA Y SACUDIMIENTO ............................ 58
8.3
TORQUE DE LA FUERZA DEL GAS Y EFECTO DEL DESCENTRADO .. 59
8.4
TORQUETOTAL DEL MOTOR .................................................................... 59
9. IMPLEMENTACIÓN HOJA DE CÁLCULO .................................................... 61
9.1
VARIABLES ................................................................................................. 61
9.2
DESCRIPCIÓN DE HOJA DE CÁLCULO .................................................. 63
10. ANÁLISIS DE RESULTADOS ......................................................................... 70
11. CONCLUSIONES ............................................................................................ 83
6
BIBLIOGRAFÍA ...................................................................................................... 85
7
LISTA DE FIGURAS
pág.
Figura 1. Fotografía del prototipo de motor Stirling existente
17
Figura 2. Dos pistones y un cilindro del motor prototipo de motor Stirling 17
Figura 3. Motor Original de Robert Stirling
21
Figura 4 Diagrama P-V y T-S del ciclo Stirling
22
Figura 5. Procesos del ciclo Stirling
23
Figura 6. Motor tipo alfa con disposición en V
24
Figura 7. Esquema del motor Stirling tipo alfa
25
Figura 8. Esquema del Motor Stirling tipo Beta
25
Figura 9. Esquema del Motor Stirling tipo gamma
26
Figura 10. Mecanismo biela manivela
27
Figura 11. Diagrama de flujo para la metodología
31
Figura 12 Parámetros geométricos del mecanismo manivela -biela- pistón 33
Figura 13. Parámetros geométricos del mecanismo manivela - biela desplazador
33
Figura 14. Geometría del mecanismo
35
Figura 15. Esquema para ángulo de transmisión µ en un mecanismo bielamanivela
38
Figura 16. Modelos dinámicos de masa concentrada de la biela
41
Figura 17. Dinámico masa concentrada estáticamente equivalente de
manivela
41
Figura 18. Masa concentrada de manivela-biela para pistón y desplazador 43
8
Figura 19. . Modelo dinámico masa concentrada de mecanismo bielamanivela con h=0
44
Figura 20. Diagrama de cuerpo libre de un mecanismo biela manivela
con h=0
45
Figura 21. Parámetros geométricos y termodinámicos
47
Figura 22. Volúmenes de trabajo para el análisis de termodinámico
49
Figura 23. Diagrama de fasores para desplazamientos del pistón y
desplazador
51
Figura 24. Diagrama fasores para volúmenes de desplazamiento de
pistón y desplazador
52
Figura 25. Presión y fuerza del gas sobre el pistón
56
Figura 26. Diagrama cuerpo libre análisis de fuerza del gas mecanismo
manivela-corredera
56
Figura 27. Pantalla principal hoja de calculo
62
Figura 28 Bloque de variable dimensional para los mecanismos
manivela biela del pistón y desplazador.
64
Figura 29. Bloque de variable dimensional para los mecanismos
manivela biela del pistón y desplazador
65
Figura 30. Bloque correspondiente a las dimensiones del pistón,
desplazador y sus volúmenes muertos.
65
Figura 31 Asignación de variables en celdas de hoja de cálculo para
espacios de pistón y desplazador.
66
Figura 32. Bloque de propiedades del gas de trabajo
67
Figura 33. Bloque para variables de propiedades del gas de trabajo
67
Figura 34. Ubicación de las componentes de las masas equivalentes para
los puntos A, A’, B, B’ en el modelo del mecanismo manivela pistón.
68
Figura 35. Valores iníciales para análisis de configuraciones del motor
9
70
Figura 36. Presión del gas para la primera configuración del motor con
h = 0 m.
71
Figura 37. Valores para el torque en la primera configuración del motor
con h = 0
71
Figura 38. Valores para la segunda configuración h=0,01 m
72
Figura 39. Curva del gas de trabajo para la segunda disposición
h = 0,01 m.
72
Figura 40. Curva del torque para segunda disposición del motor con
h = 0,01 m
73
Figura 41. Comportamiento de la presión para h= 0,02m
73
Figura 42. Comportamiento del torque para h= 0,02m
74
Figura 43. Comportamiento de la presión para h=0,04 m
74
Figura 44. Comportamiento del torque total para h=0,04 m
75
Figura 45. Parámetros geométrico para la quita configuración
75
Figura 46. Comportamiento de la presión para h=0,055 m
76
Figura 47. Comportamiento del torque total para h=0,055 m
76
Figura 48. Curva de presión para el sexto arreglo con biela = 0,14
77
Figura 49. Comportamiento del torque total para h=0,14 m
77
Figura 50. Variables para séptimo arreglo a comparar diámetro
desplazador = 0,1 m
78
Figura 51. Curva de presión séptimo arreglo diámetro desplazador = 0,1 m 78
Figura 52. Curva de torque séptimo arreglo diámetro desplazador = 0,1 m 78
Figura 53. Variables para el octavo arreglo del motor Stirling a comparar
79
Figura 54. Curva de presión para un desplazador < pistón
79
Figura 55. Curva del torque para un desplazador < pistón
80
10
Figura 56. Valores que caracterizan la novena configuración del motor.
80
Figura 57. Curva de presión para la novena configuración
81
Figura 58. Curva de torque total para la novena configuración
81
Figura 59. Parámetros que caracterizan la décima configuración
82
Figura 60. Curva de presión para densidad del aire = 1,7 Kg/m3
82
Figura 61. Curva de torque para densidad del aire = 1,7 Kg/m3
82
11
LISTA DE ANEXOS
Anexo A. Análisis termodinámico-mecánico de un prototipo de motor Stirling de
configuración gamma y de baja potencia (ver en cd)
12
RESUMEN
El presente documento de grado, ha sido desarrollado con la pretensión de poner a
consideración de la academia, un análisis termodinámico-mecánico de un prototipo
de motor stirling de configuración gamma de baja potencia, procedimiento llevado a
cabo por el autor de este proyecto. En consecuencia, la primera actividad que se
realizó, fue la de describir del prototipo existente, para exhibir los diferentes
compuestos y partes de dicho prototipo en aras de contextualizar al lector.
Tomando como referente que el propósito capital de este análisis era brindar una
interpretación de la condición de funcionamiento actual de un prototipo de motor
Stirling tipo Gamma de baja potencia construido por el autor, a partir de un análisis de
ingeniería, se formuló un modelo termodinámico/mecánico de apoyo sobre el
funcionamiento del prototipo del motor Stirling objeto de este proyecto. Para tal final,
se hizo una selección de las variables relevantes en el funcionamiento del motor y
también del rango de sus valores
En la parte final del documento, se propone un plan sistemático de simulaciones
utilizando el modelo matemático formulado, orientado a interpretar el funcionamiento
del prototipo.
13
INTRODUCCIÓN
El consumo de energía es necesario para el desarrollo económico y social. Sin
embargo, en los últimos años, las concentraciones de gases con efecto invernadero
están creciendo rápidamente en el planeta como consecuencia del uso generalizado
de los combustibles fósiles, provocando cambios drásticos en el clima mundial y
haciéndolo cada vez más impredecible.1
El desarrollo de fuentes energéticas renovables y el perfeccionamiento de tecnologías
que permitan su utilización en los diferentes escenarios industriales pueden marcar la
pauta en el camino hacia un uso sostenible de los recursos del medio ambiente. Por
ello hay un interés a nivel mundial por el desarrollo de nuevas tecnologías que
permitan la generación de potencia en forma limpia y económica. En este contexto
emergen las posibilidades que ofrece el motor tipo Stirling.
Como éste necesita solamente una fuente de calor externa es posible usar una gran
variedad de fuentes energéticas (energía solar, biomasa, energía geotérmica,
etcétera), lo que ha motivado el estudio de su desempeño a través de modelos
matemáticos.
Dado que el ciclo Stirling puede adaptarse a una variedad de configuraciones
mecánicas para obtener trabajo útil, cualquier intención de diseño debe estar
respaldada por un análisis termodinámico-mecánico, aplicado a algún caso particular.
Así, este trabajo propone realizar un análisis tal a un prototipo de motor Stirling tipo
Gamma de baja potencia, construido por el autor de este trabajo de grado, del cual
puedan obtenerse conclusiones sobre su puesta en marcha.
Es tiempo de cambiar sustentabilidad a partir de energías limpias [en línea]. En: revista de
humanidades. 2009 [consultado 20 de marzo de 2015]. Disponible en Internet:
http://www.up.edu.mx/files_uploads/20150_Revista_Humanidades_UP_2009.pdf
1
14
1. ANTECEDENTES
El gran desarrollo de los motores de combustión interna a partir de la mitad del siglo
XIX y la mejora experimentada en el refinamiento de los derivados del petróleo, colocó
a los motores alternativos de combustión externa en gran desventaja. Este hecho,
acompañado de la invención de los motores eléctricos, consiguió que desde principios
del siglo XX, la máquina de vapor y los motores Stirling y Ericsson fueran dejados de
lado en la carrera por la industrialización.[A]
Hacia mediados del siglo XX aparece un renovado interés en estos dispositivos para
nuevas aplicaciones como refrigeración, calefacción y generación eléctrica – incluso
automoción-, a partir de fuentes de calor alternativas a los combustibles fósiles. Martini
[A] preparó un documento donde se sistematiza la información al respecto, y ante
todo, las estrategias de modelado que marcan la pauta en los estudios más recientes.
Este trabajo ha sido incluso reeditado [B];
Por la manera en que está estructurado el trabajo de Martini[A] y otros trabajos que
se citan como referencia central [H] puede decirse que ofrece una estrategia plausible
para lo que se propone mediante este trabajo.
Una muestra de la literatura reciente con respecto al diseño de motores Stirling [C, D,
E, F] confirma que, un principio subyacente a cualquier esfuerzo sistemático por lograr
su funcionamiento confiable, es trabajar sobre la base de algún modelo que tenga
forma matemática.
Como motivación importante para este documento se considera el trabajo realizado
por Robson y colaboradores [G], donde se toma un prototipo de motor Stirling como
un sistema dinámico conformado por unidades de compresión, expansión, movimiento
del pistón de desplazamiento libre y pistón de potencia de salida. Las predicciones
logradas fueron consideradas satisfactorias por sus autores.
Por otra parte, alrededor del año 2008, Al cursar la materia “Proyecto integrador” en
la Universidad Autónoma de Occidente, el objetivo general de dicho curso fue la
construcción de un motor Stirling. En este momento el autor elaboró un prototipo de
dicha máquina a partir de la teoría del ciclo termodinámico del mismo nombre. En
algún momento de la experimentación con el modelo, éste funcionó por un tiempo
corto y desde entonces no ha sido posible reanudar su funcionamiento, a pesar de
una serie de variaciones realizadas.
15
Finalmente, y en virtud de la consulta bibliográfica llevada a cabo, se piensa que
debería formularse un análisis de tipo termodinámico-mecánico, por el que pueda
darse una interpretación lógica a las dificultades observadas con este prototipo.
A continuación se presenta una descripción del prototipo objeto de estudio.
1.1
DESCRIPCION DEL PROTOTIPO EXISTENTE
El modelo de motor Stirling construido por el autor, tiene una configuración gamma
que dispone de dos pistones de movimiento horizontal y paralelo unidos a un cigüeñal
por dos bielas, ver figura 1.
Los dos pistones tienen un diámetro de 15mm y al igual que los cilindros son de vidrio;
la longitud del pistón del lado caliente es de 64mm mientras que. El pistón lado frío es
de 40mm y se muestran en la figura 2. Las bielas son fabricadas en aluminio y tienen
una serie de perforaciones para reducir su peso.
El eje del cigüeñal es de acero inoxidable con 6mm de diámetro y 73mm de longitud
y rueda en dos cojinetes de bolas. Adicionalmente el prototipo cuenta con un volante
fabricado en bronce con un peso de 80 gr.
La recámara caliente está compuesta por dos elementos: el elemento que recibe el
calor, fabricado en cobre, y su parte posterior, en teflón, que es el dispositivo
encargado de alojar el cilindro de vidrio. La recámara del lado frío es una sola pieza,
construida en cobre y provista de aletas para ayudar a disipar de calor.
El chasis de este prototipo de motor Stirling es un bloque realizado en láminas de
acero inoxidable atornillado a una base de madera. Estas láminas están separadas
por un bloque de empack D usado en este caso como un aislante térmico. La fuente
de calor es un mechero que usa como combustible alcohol industrial.
Este motor opera en un ciclo cerrado que usa aire como fluido de trabajo,
desplazándose alternativamente entre el lado frío y el lado caliente del motor. Al
hacerlo pasa por un regenerador fabricado en lámina de cobre enrollada dentro de un
tubo -también de cobre- de 10mm de diámetro y 17mm de longitud.
16
Figura 1. Fotografía del prototipo de motor Stirling existente
Figura 2. Dos pistones y un cilindro del motor prototipo de motor Stirling
17
2. JUSTIFICACIÓN
Hoy día el agotamiento de las fuentes de energía no renovables y el alto nivel de
contaminación del medio ambiente en el planeta, reclaman procesos industriales cada
vez más eficientes y limpios. En este contexto los motores de combustión eficiente
pueden alcanzar grandes ahorros energéticos y reducciones en las emisiones de CO 2.
Dadas las características de su funcionamiento, los motores Stirling pueden cumplir
con tales requerimientos.
Por este motivo se ha despertado un interés en la industria y la academia de todo el
mundo para desarrollar programas de cara al diseño y experimentación de motores
Stirling suficientemente eficientes para competir con las soluciones actuales a la crisis
energética.
Dado que los países latinoamericanos no se escapan a los problemas energéticos del
planeta, es necesario promover también la investigación acerca de estos motores. Ello
implica, en primer lugar, tener un prototipo respaldado por un análisis desde el punto
de vista de las ciencias de la ingeniería (mecánica y termodinámica), que ayude a
participar de manera activa en las discusiones sobre el tema.
18
3.
OBJETIVOS
3.1 OBJETIVO GENERAL
Dar una interpretación de la condición de funcionamiento actual de un prototipo de
motor Stirling tipo Gamma de baja potencia construido por el autor, a partir de un
análisis de ingeniería.
3.2 OBJETIVOS ESPECÍFICOS
 Formular un modelo termodinámico/mecánico de apoyo sobre el funcionamiento
del prototipo del motor Stirling objeto de este proyecto.
 Hacer una selección de las variables relevantes en el funcionamiento del motor y
también del rango de sus valores
 Proponer y llevar a cabo un plan sistemático de simulaciones utilizando el modelo
matemático formulado, orientado a interpretar el funcionamiento del prototipo.
19
4. MARCO TEÓRICO
En 1816 Robert Stirling, reverendo de origen escocés proyecta el motor Stirling con el
objetivo de lograr un motor menos peligroso que la máquina de vapor que trabajaban
a altas presiones.
El motor Stirling es un dispositivo que opera entre una fuente a temperatura alta y un
sumidero a baja temperatura. Es decir, es necesaria la presencia de una diferencia de
temperaturas entre dos focos, lo cual lo caracteriza como un motor térmico. El calor
que toma de la fuente a alta temperatura lo convierte en trabajo, usando como medio
de trabajo un gas en un ciclo termodinámico cerrado conocido como el ciclo Stirling.
Este motor continúa siendo investigado debido a la gran variedad de fuentes de
energía que pueden ser utilizadas para su funcionamiento como la energía solar,
geotérmica, biomasa, etc.
La figura 3 muestra el esquema de un motor Stirling original y la empleamos para
explicar su funcionamiento. Un cilindro vertical es calentado en su parte superior por
el flujo de gases calientes provenientes de la caldera. En el interior del cilindro se
encuentra un pistón de potencia y un desplazador. El desplazador es liviano y no es
buen conductor de calor; al medio del desplazador existe un anillo de material capaz
de absorber y ceder calor que es el regenerador. Cuando el desplazador se mueve
hacia abajo, la mayor parte del aire que se encuentra dentro del cilindro queda en la
zona caliente y se expande, empujando el pistón de trabajo hacia abajo. Aquí se
entrega trabajo al exterior y gira el volante. Al suceder esto, una serie de bielas
mueven el desplazador hacia arriba, desplazando la mayor parte del aire a través del
regenerador hacia la zona fría donde se comprime el aire y se inicia nuevamente el
ciclo.
4.1 CICLO STIRLING
El ciclo Stirling ideal consiste de dos procesos isotérmicos totalmente reversibles y
dos procesos reversibles a volumen constante. La figura 4. Muestra los procesos del
ciclo en una curva p-v y en un cura T-s.
Durante el proceso de expansión 3-4 el calor es suministrado a temperatura constante
TH y durante el proceso 1-2 el calor es rechazado a temperatura constante TL. De 2-3
el calor interactúa a volumen constante al igual que de 4-1. La cantidad de calor en
estos dos procesos es esencialmente igual, pero en dirección opuesta; este proceso
de intercambio necesita ser realizado por un regenerador. La función del regenerador
20
es actuar como un reservorio temporal, capaz de absorber calor durante el proceso
de 4-1 y entregar idealmente la misma cantidad de calor durante el proceso de 2-3.
Figura 3. Motor Original de Robert Stirling
Fuente: El motor sterling [en línea]. Chile: Universidad de Chile, 1998 [consultado 20
de marzo de 2015]. Disponible en Internet:
http://www.cec.uchile.cl/~roroman/cap_10/STRLNG2.HTM
21
Figura 4 Diagrama P-V y T-S del ciclo Stirling
Fuente: SAAD, Michael. Thermodynamics principles and practice. México: Prentice
Hall, 1997. capítulo 7.
En la figura 5 se representa un sistema que puede ejecutar los diferentes procesos
del ciclo Stirling.
Al fluido de trabajo se le añade calor isotérmicamente de una fuente externa de
temperatura TH durante el proceso 1-2, y se rechaza también isotérmicamente en un
sumidero externo a temperatura TL durante el proceso 3-4. En un proceso isotérmico
reversible, la transferencia de calor se relaciona con el cambio de entropía.
En una carrera permanente por realizar mejoras al motor Stirling se han generado
innumerables arreglos los cuales se clasifican dentro de tres grupos principales que
son los motores tipo alfa, beta y gamma.
22
Figura 5. Procesos del ciclo Stirling
Fuente: SAAD, Michael. Thermodynamics principles and practice. Mèxico: Prentice
Hall, 1997. capítulo 7.
23
4.2 MOTOR STIRLING TIPO ALFA
Es posible que esta sea la configuración más sencilla de las tres, siendo conformada
por dos pistones en cilindros diferentes conectados en serie mediante un regenerador.
Uno de los cilindros está conectado a un calentador y el otro a un enfriador. Richard
Wheeler2 desarrolló un motor tipo alfa con disposición en V, este se muestra en la
figura 6
Figura 6. Motor tipo alfa con disposición en V
Fuente: Stirling engines - mechanical configurations [en línea]. Ohio: Universidad de
Ohio [consultado 20 de marzo de 2015]. Disponible en Internet:
http://www.ohio.edu/mechanical/stirling/engines/engines.html
Un diagrama esquemático de esta configuración se muestra en la figura 7, donde se
observa la disposición de los elementos termodinámicos de la configuración alfa, es
de notar el regenerador entre el intercambiador de calor caliente y el frío.
2
Stirling engines - mechanical configurations [en línea]. Ohio: Universidad de Ohio [consultado 20 de
marzo de 2015]. Disponible en Internet: http://www.ohio.edu/mechanical/stirling/engines/engines.html
24
Figura 7. Esquema del motor Stirling tipo alfa
Fuente: Stirling engines - mechanical configurations [en línea]. Ohio: Universidad de
Ohio [consultado 20 de marzo de 2015]. Disponible en Internet:
http://www.ohio.edu/mechanical/stirling/engines/engines.html
4.3 MOTOR STIRLING TIPO BETA
El motor Stirling original de acuerdo a la patente del dibujo de 1816, se muestra un
tipo beta. En esta configuración el pistón y el desplazador se encuentran en un mismo
cilindro el cual tiene en un extremo la zona caliente y en el otro la zona fría. El oficio
del desplazador es impulsar el gas de la zona fría a la caliente y viceversa. Un
esquema de esta configuración lo vemos en la figura 8.
Figura 8. Esquema del Motor Stirling tipo Beta
Fuente: Stirling engines - mechanical configurations [en línea]. Ohio: Universidad de
Ohio [consultado 20 de marzo de 2015]. Disponible en Internet:
http://www.ohio.edu/mechanical/stirling/engines/engines.html
25
4.4 MOTOR STIRLING TIPO GAMMA
El motor tipo Gamma tiene una configuración similar a la máquina del tipo beta, sin
embargo el desplazador y el pistón se encuentran en distintos cilindros. Esto permite
una separación completa entre el cilindro del desplazador y el cilindro que contiene el
pistón, lo cual hace que los intercambiadores de calor estén separados el uno del otro,
por lo tanto, tienden a tener un poco más grandes los volúmenes de las zonas muertas
que cualquiera de los motores alfa o beta. Ver figura 9.
Figura 9. Esquema del Motor Stirling tipo gamma
Fuente: Stirling engines - mechanical configurations [en línea]. Ohio: Universidad de
Ohio [consultado 20 de marzo de 2015]. Disponible en Internet:
http://www.ohio.edu/mechanical/stirling/engines/engines.html
4.5 MECANISMO BIELA MANIVELA
El mecanismo de biela - manivela es un mecanismo que transforma un movimiento
circular en un movimiento de traslación, o viceversa. El ejemplo actual más común se
encuentra en el motor de combustión interna de pistones, en el cual el movimiento
lineal del pistón producido por la explosión del combustible, se trasmite a la biela y se
convierte en movimiento circular en el cigüeñal.
26
4.5.1 Estudio cinemático. Definimos la dirección de la coordenada X de manera
paralela a la dirección del movimiento cruceta - pistón y definimos el origen de
coordenadas correspondiente en el centro de rotación de la manivela; el problema
consiste en hallar las funciones que determinan la posición, la velocidad lineal y la
aceleración también lineal de la cruceta en función de la velocidad angular ω de la
manivela y del ángulo α de rotación de la misma.
Siendo:
R = Longitud de la manivela
L = Longitud de la biela
Ө = Ángulo de rotación descripto por la manivela
Ø = Ángulo de rotación descripto por la biela alrededor de su articulación con la cruceta
o pistón.
Las ecuaciones para el desplazamiento, velocidad y aceleración se obtienen a partir
de la figura 10.
Figura 10. Mecanismo biela manivela
Fuente: NORTON, Robert. Diseño de maquinaria. Cuarta edición México: Prentice
Hall, 2011. 584 p.
X = R + L – RcosӨ – Lcos Ø
(3.1)
X = R( 1 – cosӨ) + L(1 - cos Ø)
(3.2)
2
𝑋 = 𝑅(1 − 𝑐𝑜𝑠𝜃) + 𝐿 (1 − √1 − (𝑅⁄𝐿) 𝑠𝑒𝑛2 𝜃)
27
(3.3)
Si tiene en cuenta que
2
2
√1 − (𝑅⁄ ) 𝑠𝑒𝑛2 𝜃 = 1 − 1 (𝑅) 𝑠𝑒𝑛2 𝜃
𝐿
2 𝐿
(3.4)
La ecuación (15) se puede simplificar:
𝑅2
(3.5)
𝑋 = 𝑅(1 − 𝑐𝑜𝑠𝜃) + 2𝐿 𝑠𝑒𝑛2 𝜃
Donde
serán:
𝑉=
𝐴=
𝑑𝑥
𝑑𝑡
debido a que se toma
como constante. La velocidad y la aceleración
𝑅
(3.6)
= 𝑅𝑤 (𝑠𝑒𝑛𝜃 + 2𝐿 𝑠𝑒𝑛2𝜃)
𝑑2 𝑥
𝑑𝑡 2
𝑅
(3.7)
= 𝑅𝑤 2 (𝑐𝑜𝑠𝜃 + 𝐿 𝑐𝑜𝑠2𝜃)
28
5. METODOLOGÍA
La formulación del modelo termodinámico-mecánico de un prototipo de motor Stirling
tipo gamma inició por la parte mecánica, tomando como patrón el análisis dinámico
de un mecanismo manivela-deslizador que se presenta en el texto “Diseño de
Maquinaria” de Norton, indicado en la bibliografía al final del informe. Dicho análisis
viene precedido por un análisis cinemático, en el cual se establecen expresiones para
la posición, velocidad y aceleración del deslizador dada una posición angular de la
manivela.
A diferencia del estudio presentado en el texto del profesor Norton, el análisis que se
desea hacer al prototipo incluye un desfase u “Offset”, que es la longitud del segmento
que parte del eje de la manivela y es perpendicular a la línea de eje del deslizador,
cuando se mira al mecanismo perpendicular al plano descrito por la rotación de la
manivela. Luego, la obtención de una expresión acorde con las simplificaciones que
presenta el profesor Norton en su texto (omisión de términos armónicos de alto orden)
obligó a una reformulación propia.
El análisis dinámico, en cuanto a las expresiones para el torque de gas y el torque de
inercia, no presenta variaciones con respecto de la referencia utilizada, a excepción
del valor de la presión del gas, que fue el objeto del análisis termodinámico.
El análisis termodinámico es el que se presenta bajo el nombre de Gustav Schmidt,
tal como se expone en el portal “Stirling cycle machine analysis” de la universidad de
Ohio (EEUU), cuyo enlace se indica en la bibliografía al final de este documento. Se
trata de un análisis donde se asume un régimen isotérmico para el gas en los espacios
de compresión y expansión tanto barridos por el desplazador y el pistón, como los
correspondientes a volúmenes muertos. Así mismo, se asume un perfil lineal de
temperatura a través del regenerador, o zona que media entre los volúmenes de aire
frío y caliente.
El análisis Schmidt produce como resultado un perfil de presión de trabajo en función
de la coordenada angular de la manivela. Este resulta altamente dependiente de la
fase que se tenga entre el movimiento del pistón y el del desplazador, que a su vez va
a incidir en la variación periódica del volumen de compresión. Esta relación se hace
evidente a través del denominado “diagrama fasor”, que aunque está disponible en la
referencia indicada, hubo necesidad de reformularlo para hacer evidente los términos
que representan los parámetros de diseño en el presente trabajo. Esta acción
simplifica la implementación de las expresiones resultantes en hoja de cálculo y
29
también la explicación a quien no estuviera familiarizado con el funcionamiento del
motor Stirling.
La expresión para la presión del gas, que resulta del análisis termodinámico, se
reemplaza en su lugar correspondiente en la del torque de gas; esta última sumada al
torque de inercia resulta en la curva de torque total. Dichas expresiones muestran las
siguientes cantidades como parámetros de diseño: longitud de manivela y de biela del
lado del pistón de potencia y del desplazador, masas equivalentes en los mecanismos
de manivela-deslizador del lado del pistón de potencia y del desplazador, diámetros
del pistón de potencia y del desplazador, longitudes de los tramos muertos de
compresión y de expansión, longitud y diámetro del regenerador, temperaturas de lado
frío y de lado caliente, y densidad y presión del aire en la cámara a temperatura
ambiente (condición de inicio del motor)
La curva de presión de gas y la de torque total, se pueden utilizar para estudiar el
efecto de estos parámetros en el tamaño del volante requerido y en el trabajo neto de
salida del motor.
Con el fin de sintetizar la metodología expuesta en la figura 11 se presenta un
diagrama de flujo con los pasos que se sigue en este proyecto para lograr los
objetivos.
30
Figura 11. Diagrama de flujo para la metodología
Mecanismo a considerar
Análisis cinemático: posición,
velocidad y aceleración
x , r , l , h 
x  , r , l , h 
x , r , l , h 
ANALISIS
MECÁNICO
TTOTAL  TGAS  TINERCIA
Torque debido
a la acción del
gas sobre el
pistón
Modelo dinámico de
masas concentradas
TGAS
TINERCIA
Torque debido
a las fuerzas de
inercia del
mecanismo
ANALISIS
TERMODINÁMICO
Motor Stirling
tipo Gamma
31
Presión al
interior de la
cámara
p 
6. FORMULACIÓN DEL ANÁLISIS MECÁNICO
6.1 PARÁMETROS CINEMÁTICOS
El mecanismo de un motor Stirling de configuración gamma está constituido por un
pistón y un desplazador, cada uno de estos elementos posee una biela y una manivela
las cuales los unen a un mismo eje. Para realizar el análisis cinemático del motor se
asume que:
-Las bielas del pistón y del desplazador tienen las mismas longitudes.
-Las manivelas del pistón y del desplazador tienen las mismas longitudes.
-La presión del gas de trabajo se ejerce solamente sobre el pistón.
Para este estudio se eligió una configuración de manivela-biela descentrada, es decir
que el eje de la corredera no está centrado con el pivote de la biela, esto con el fin de
abarcar una mayor cantidad de configuraciones para dicho motor.
En las figuras 12 y 13 se muestran los mecanismos manivela - biela para el pistón y
el desplazador en la configuración propuesta, donde r es la longitud de la manivela, l
la longitud de la biela y h es distancia entre el eje de la corredera y el eje de la
manivela.
Las longitudes r, l y h son variables de entrada para el análisis y es necesario que al
darle valores se configure un mecanismo que al girar la manivela no presente posición
de agarrotamiento y sea capaz de dar una revolución completa alrededor de su eje
32
Figura 12 Parámetros geométricos del mecanismo manivela -biela- pistón
Figura 13. Parámetros geométricos del mecanismo manivela - biela desplazador
La condicione básica para que la ecuación de desplazamiento de la corredera (esta
se verá en el siguiente numeral) no quede indeterminada es que los valores dados a
las variables r, hy l cumplan con:
Si h es diferente de cero:
𝑟+ℎ
𝑙
(6.1)
<1
Si h es igual a cero:
𝑟
𝑙
(6.2)
<1
33
6.2 CINEMÁTICA DEL MECANISMO
Por parámetros cinemáticos se entienden los que resultan del análisis del mecanismo
cuando el propósito es obtener expresiones para la posición, velocidad y aceleración
del deslizador en función de la posición angular de la manivela.
Un mecanismo de biela y corredera puede ser analizado empleando los métodos
usados en el análisis de mecanismos de cuatro barras con juntas de pasadores. En la
figura 14 se observa la geometría de un mecanismo de corredera manivela de cuatro
barras descentrado lo cual significa que el eje de la corredera no pasa por el punto de
pivote de la manivela.
Este mecanismo se puede representar por cuatro vectores de posición R1, R2, R3, R4,
donde R1 equivale a la posición de la corredera con magnitud X. R4 corresponde a
una magnitud constante h la cual define el descentre entre el eje de la corredera y el
pivote de la biela. El vector R2 define la longitud de la manivela r y el vector R3
representa la magnitud L dela biela que es elemento acoplador entre la manivela y el
pistón.
Para el análisis se crea un eje de coordenadas paralelo al eje de deslizamiento de la
corredera. El vector R3 que representa el acoplador se coloca con su origen en la
corredera.
El ángulo ɵ es el ángulo de la manivela y para una velocidad angular constante w este
ángulo está dado por ɵ= wt. El ángulo formado por la biela con el eje x es Ø
De la figura 14 se construye el triángulo rsq y el triángulo ul(q-h)y por geometría se
tiene que:
(6.3)
𝑞 = 𝑟 sin 𝜃 = (𝑙 sin ∅) + ℎ
De (6.3) se despeja senØ y se tiene:
sin ∅ =
𝑟 sin 𝑤𝑡−ℎ
(6.4)
𝑙
34
Figura 14. Geometría del mecanismo
De Figura 14 también se tiene que:
𝑠 = 𝑟𝑐𝑜𝑠𝑤𝑡 y 𝑢 = 𝑙 cos ∅ y𝑋 = 𝑆 + 𝑈
(6.5)
𝑥 = 𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙𝑐𝑜𝑠∅
De la identidad trigonométrica se tiene que:
(6.6)
𝑐𝑜𝑠∅ = √𝑙 − 𝑠𝑒𝑛2 ∅
Reemplazando (6.4) en (6.6)
2
2
𝑟𝑠𝑒𝑛𝑤𝑡−ℎ
𝑐𝑜𝑠∅ = √1 − (
)
𝑙
(6.7)
Reemplazando (6.7) en (6.5)
𝑥 = 𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 √1 − (
𝑟𝑠𝑒𝑛𝑤𝑡−ℎ 2
𝑙
(6.8)
)
El radical de la ecuación (6.8) es
2
√1 − (𝑟𝑠𝑒𝑛𝑤𝑡−ℎ)
𝑙
(6.9)
35
La ecuación (6.9) se expresa como:
(1 − (
1/2
𝑟𝑠𝑒𝑛𝑤𝑡−ℎ 2
𝑙
(6.10)
) )
1/2
ℎ 2
𝑟
(6.11)
= (1 − ( 𝑙 𝑟𝑠𝑒𝑛𝑤𝑡 − 𝑙 ) )
Se emplea el teorema binomial para expandir el radical de la ecuación (6.8) que nos
da la posición pistón. La forma general del teorema binomial ó del también llamado
binomio de Newton expresa la enésima potencia de un binomio como un polinomio:
(𝑎 + 𝑏)𝑛 = 𝑎𝑛 + 𝑛𝑎𝑛−1 𝑏 +
𝑛(𝑛 − 1) 𝑛−2 2 𝑛(𝑛 − 1)(n − 2) 𝑛−3 3
𝑎 𝑏 +
𝑎 𝑏 +⋯
2!
3!
Para la expansión binomial de (6.11) se tiene que:
ℎ 2
𝑟
a=1, 𝑏 = − ( 𝑙 𝑠𝑒𝑛𝑤𝑡 − 𝑙 ) , n= ½
Entonces (6.11) se expande a:
1 𝑟
ℎ 2
1 − 2 ( 𝑙 𝑠𝑒𝑛𝑤𝑡 − 𝑙 ) +
ℎ 4
1 𝑟
1
𝑟
ℎ 6
( 𝑠𝑒𝑛𝑤𝑡 − 𝑙 ) − 16 ( 𝑙 𝑠𝑒𝑛𝑤𝑡 − 𝑙 ) + ⋯
8 𝑙
(6.12)
Desarrollando los cuatro primeros términos de la expresión binomial se tiene que:
1 𝑟2
1 − 2 ( 𝑙2 𝑠𝑒𝑛2 𝑤𝑡 −
2𝑟
𝑙
𝑠𝑒𝑛𝑤𝑡 +
1
ℎ2
𝑙2
1
) + 𝟖𝒍𝟒 (𝑟 4 𝑠𝑒𝑛4 𝑤𝑡 − 4𝑟 3 ℎ𝑠𝑒𝑛3 𝑤𝑡 + 6𝑟 2 ℎ2 𝑠𝑒𝑛2 𝑤𝑡 −
4𝑟ℎ3 𝑠𝑒𝑛𝑤𝑡 + ℎ4 ) − 𝟏𝟔𝒍𝟔 (𝑟 6 𝑠𝑒𝑛6 𝑤𝑡 − 6𝑟 5 ℎ𝑠𝑒𝑛5 𝑤𝑡 + 15𝑟 4 ℎ2 𝑠𝑒𝑛4 𝑤𝑡 − 20𝑟 3 ℎ3 𝑠𝑒𝑛3 𝑤𝑡 +
15𝑟 2 ℎ4 𝑠𝑒𝑛2 𝑤𝑡 − 6𝑟ℎ5 𝑠𝑒𝑛𝑤𝑡 + ℎ6 )
(6.13)
Si se tiene en cuenta que l siempre será mayor que r, es decir que la relación 𝑟⁄𝑙 será
menor que uno. Entonces para (6.13) los elementos a partir del tercer término se
harán cada vez muy pequeños por lo que la podemos simplificar a:
1 𝑟2
1 − 2 ( 𝑙2 𝑠𝑒𝑛2 𝑤𝑡 −
2𝑟
𝑙
𝑠𝑒𝑛𝑤𝑡 +
ℎ2
𝑙2
(6.14)
)
Se reemplaza la expresión (6.14) en la ecuación (6.8) y se tiene que:
36
1 𝑟2
𝑥 = 𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 (1 − 2 ( 𝑙2 𝑠𝑒𝑛2 𝑤𝑡 −
2𝑟
𝑙
𝑠𝑒𝑛𝑤𝑡 +
ℎ2
𝑙2
))
(6.15)
En la ecuación (6.15) se sustituye:
𝑠𝑒𝑛2 𝑤𝑡 = 1−𝑐𝑜𝑠2𝑤𝑡
2
(6.16)
𝑟2
𝑟2
ℎ2
𝑥 = 𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 − 4𝑙 + 4𝑙 𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − 2𝑙
(6.17)
La ecuación (6.17) es un modelo matemático bastante aproximado para describir el
desplazamiento del pistón en la corredera.
Una expresión para la velocidad del pistón se encuentra derivando la ecuación (6.17)
𝑥 ′ = −𝑟𝑤𝑠𝑒𝑛𝑤𝑡 −
𝑟2𝑤
2𝑙
(6.18)
𝑠𝑒𝑛2𝑤𝑡 + 𝑟𝑤𝑐𝑜𝑠𝑤𝑡
La aceleración de un punto en x se obtiene derivando la ecuación de velocidad (6.18)
𝑥 ′′ = −𝑟𝑤 2 𝑐𝑜𝑠𝑤𝑡 −
𝑟2𝑤2
𝑙
𝑐𝑜𝑠2𝑤𝑡 − 𝑟𝑤 2 𝑠𝑒𝑛𝑤𝑡
(6.19)
Del análisis presentado es claro que además de la posición angular de la manivela
reflejada por el término t , se tiene el largo de la manivela r y el de la biela l como
parámetros cinemáticos.
6.3 ANGULO DE TRANSMISIÓN
El ángulo de transmisión definido como el ángulo formado entre los eslabones
acoplador y de salida en un mecanismo de cuatro barras, es una forma rápida de
apreciar la eficiencia con la que la potencia es transmitida. Este ángulo varia
continuamente entre un valor mínimo y un valor máximo conforme el eslabonamiento
pasa por su intervalo en movimiento. En la figura 15 se muestra un mecanismo de
cuatro barras donde la barra 4 es la motriz, la barra 3 es el elemento acoplador y la
barra 2 es el eslabón de salida. Por lo tanto el ángulo de transmisión es el formado
entre las barras 2 y 3.
37
Figura 15. Esquema para ángulo de transmisión µ en un mecanismo bielamanivela
De la ecuación (6.4) se despeja Ø para obtener:
𝑟 sin 𝑤𝑡−ℎ
∅ = 𝑎𝑟𝑐𝑠𝑒𝑛 (
𝑙
(6.20)
)
Ø corresponde al ángulo menor que se pueda formar entre el elemento acoplador y la
horizontal en el origen del vector R4 que corresponde al pistón.
El ángulo entre la manivela y la biela está dado por:
(6.21)
𝜇 = 180 − (𝜃2 − ∅31 )
Por lo que el ángulo de transmisión está dado por:
𝑟 sin 𝑤𝑡−ℎ
𝜇 = 180 − (𝜃2 − (𝑎𝑟𝑐𝑠𝑒𝑛 (
𝑙
(6.22)
)))
Una implementación en hoja de cálculo permitirá visualizar aquellas configuraciones
del mecanismo (a través de la posición angular de la manivela) donde la cercanía de
µ a 0° o 180° indica una transmisión poco eficiente del movimiento.
38
6.4 MASAS EQUIVALENTES
La biela o acoplador del mecanismo manivela-deslizador tiene un movimiento plano,
lo que asociado a su naturaleza de cuerpo rígido, hace complicado considerar un
modelo en términos de su momento de inercia y del movimiento de su centro de
gravedad. Por lo anterior se considera la representación del mecanismo como un
sistema de masas puntuales equivalentes. A continuación se establecen los detalles
de dicha formulación y se concluye señalando los parámetros inerciales de diseño del
motor.
En el caso del mecanismo de manivela-corredera, la manivela está en rotación pura y
el pistón en traslación pura y sus movimientos cinemáticos son fáciles de determinar.
La biela está en movimiento complejo y para realizar un análisis dinámico se debe
determinar la aceleración lineal de su CG en todas las posiciones para lo cual se
requiere un modelo simplificado de esta biela. Los requerimientos para obtener un
modelo dinámicamente equivalente son:
 La masa del modelo debe ser igual a la masa del cuerpo original.
 El centro de gravedad debe estar en el mismo lugar que el del cuerpo original.
 El momento de inercia de masa debe ser igual al del cuerpo original.
6.4.1 Modelo dinámicamente equivalente biela. En la figura 16.se muestra un
modelo genérico de dos masas de la biela. Una masa mt se localiza a una distancia lt
del CG de la biela original y la segunda masa mp a una distancia lp del CG. La masa
de la pieza original es m3 y su momento de inercia con respecto a su CG es IG3.
Una expresión matemática para los requisitos de la equivalencia dinámica en función
de estas variables es:
𝑚𝑝 + 𝑚𝑡 = 𝑚3
(6.23)
𝑚𝑝 𝑙 𝑝 = 𝑚𝑡 𝑙 𝑡
(6.24)
𝑚𝑝 𝑙𝑝2 + 𝑚𝑡 𝑙𝑡2 = 𝐼𝐺3
(6.25)
39
Se elige la distancia lt igual a la distancia del pasador del pistón lb como se muestra
en la figura 16.c. Al resolver las ecuaciones (6.23), (6.24), (6.25) se obtiene una
expresión para las dos masas concentradas.
𝑚𝑝 = 𝑚3 𝑙
𝑙𝑏
(6.26)
𝑝 +𝑙𝑏
𝑚𝑏 = 𝑚3 𝑙
𝑙𝑝
(6.27)
𝑝 +𝑙𝑏
6.4.2 Modelo estáticamente equivalente de manivela. Puede generarse un
modelo para la manivela, su CG se localiza a cierta distancia rG2 del pivote O2 sobre
la línea del muñón de la manivela A como se muestra en la figura 17. Para facilitar el
análisis se plantea un modelo como una masa concentrada en A unida a una masa
en el pivote O2 unidas por una barra sin masa.
40
Figura 16. Modelos dinámicos de masa concentrada de la biela
Fuente: NORTON, Robert. Diseño de maquinaria. Cuarta edición México: Prentice
Hall, 201. p. 584.
Figura 17. Dinámico masa concentrada estáticamente equivalente de manivela
Fuente: NORTON, Robert. Diseño de maquinaria. Cuarta edición México: Prentice
Hall, 2011. p. 584.
41
Dado que se asume que la velocidad angular w de la manivela como constante, es
posible emplear un modelo estáticamente equivalente debido a que la aceleración es
cero y las condiciones que debe cumplir son:

La masa del modelo debe ser igual que la del cuerpo original

El centro de gravedad debe estar en el mismo lugar que el del cuerpo original
Esto se puede escribir de la siguiente manera.
𝑚2 = 𝑚2 + 𝑚2𝑂2
(6.28)
𝑚2𝑎 𝑟 = 𝑚2 𝑟𝐺2
(6.29)
𝑚2𝑎 = 𝑚2
𝑟𝐺2
(6.30)
𝑟
El muñón de la manivela tiene dos masas concentradas en el punto A, la masa
equivalente a la manivela 2 esm2A y la masa de la biela 3 esm3A. El total de estas
masas es:
(6.31)
𝑚𝐴 = 𝑚2𝐴 + 𝑚3𝐴
En el punto B se concentran las masas del pistón 4 es m4y la masa del modelo de la
biela 3 es m3B. La masa total en este punto es
(6.32)
𝑚𝐵 = 𝑚3𝐵 + 𝑚4
En la figura 18 se representan las masas totales en los puntos A y B.
Nombrando la manivela del desplazador como barra 5, la biela del desplazador como
barra 6 y el desplazador como barra 7. Se tiene que el muñón entre la manivela y la
biela del desplazador sea el punto A’ y el muñón entre la biela y la extensión del
desplazador sea el punto B’. En A’ se tiene que la masa equivalente de la manivela
del desplazador es m5A’ y la masa equivalente de la biela es m6A’.
La masa total en este punto es:
(6.31’)
𝑚𝐴′ = 𝑚5𝐴′ + 𝑚6𝐴′
42
En B’ la masa equivalente a la biela es m6B’ y la masa equivalente al desplazador y su
extensión es m7. La masa total en este puto es:
(6.32’)
𝑚𝐵′ = 𝑚6𝐵′ + 𝑚7
En la figura 18 se exponen las masa equivalentes en mA, mB, mA’ y mB’
Figura 18. Masa concentrada de manivela-biela para pistón y desplazador
De la discusión en esta sección del informe, queda claro que los parámetros inerciales
de diseño son las masas puntuales m2a, m3a, m3b, m4, m5a, m6a’, m6b’y m7 que junto con
los parámetros cinemáticos r y l se pueden utilizar para establecer una forma física
para la manivela y la biela; en este punto las limitaciones vienen por cuenta de la
facilidad de fabricación y por los esfuerzos que tengan lugar al interior de las mismas,
que ya no son el tema de este proyecto.
6.5 FUERZAS DEINERCIA Y DE SACUDIMIENTO
Utilizando el modelo simplificado de masa concentrada del mecanismo manivela
corredera para h=0 como el de la figura 19 es posible desarrollar una expresión para
las fuerzas y pares de torsión generados por las aceleraciones de las masas presentes
en el sistema. El método de d’Alembert permite visualizar los efectos de estas masa
en movimiento en el sistema del plano de la bancada.
43
Figura 19. . Modelo dinámico masa concentrada de mecanismo biela-manivela
con h=0
Fuente: NORTON, Robert. Diseño de maquinaria. Cuarta edición México: Prentice
Hall, 2011. p. 523
Los diagramas de cuerpo libre de la figura 20 muestran las fuerzas de inercia de
d’Adembert que actúan en las masas localizadas en los puntos A y B. Se ignoran las
fricciones.
La aceleración del punto A en rotación se obtiene al diferenciar dos veces el vector de
posición RA, con w constante se tiene:
𝑅𝐴 = 𝑟𝑐𝑜𝑠𝑤𝑡𝒊 + 𝑟𝑠𝑒𝑛𝑤𝑡𝒋
(6.33)
𝑎𝐴 = −𝑟𝑤 2 𝑐𝑜𝑠𝑤𝑡𝒊 − 𝑟𝑤 2 𝑠𝑒𝑛𝑤𝑡𝒋
(6.34)
La fuerza de inercia total Fi es la suma de la fuerza centrífuga en el punto A y la fuerza
de inercia en el punto B
(6.35)
𝐹𝑖 = −𝑚𝐴 𝑎𝐴 − 𝑚𝐵 𝑎𝐵
Sustituyendo la ecuación (6.34) en (6.35) se tiene que:
(6.36)
𝐹𝑖 = −𝑚𝐴 (−𝑟𝑤 2 𝑐𝑜𝑠𝑤𝑡) − 𝑚𝐵 𝑎𝐵
44
Figura 20. Diagrama de cuerpo libre de un mecanismo biela manivela con h=0
Fuente: NORTON, Robert. Diseño de maquinaria. Cuarta edición México: Prentice
Hall, 2011. p. 523
𝑎𝐵 es la aceleración del pistón que está dada por la ecuación (6.19), reemplazándola
en (6.36):
𝐹𝑖 = −𝑚𝐴 (−𝑟𝑤 2 𝑐𝑜𝑠𝑤𝑡) − 𝑚𝐵 (−𝑟𝑤 2 𝑐𝑜𝑠𝑤𝑡 −
𝑟2𝑤2
𝑙
𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑤 2 𝑠𝑒𝑛𝑤𝑡)
(6.37)
Esta es la ecuación para la fuerza de sacudimiento que se define como la suma de
las fuerzas que actúan en el plano de la bancada.
45
7. FORMULACIÓN DEL ANALISIS TERMODINÁMICO
El motor Stirling tipo gamma tiene un pistón desplazador y un pistón de potencia en
cilindros diferentes, lo que le permite una separación completa entre el intercambiador
de calor de la zona caliente y el intercambiador de calor de la zona fría. El pistón
desplazador se mueve en el volumen asociado con el lado caliente y el pistón de
potencia se mueve en el volumen de compresión asociado con el lado frio.
El volumen que se encuentra en el cilindro caliente y delante del desplazador es
considerado como volumen de expansión y se suponen a una misma temperatura T c.
El volumen que comprende desde el cilindro caliente por la parte posterior del
desplazador hasta el volumen barrido por el pistón se considera volumen de
compresión y se suponen a una misma temperatura Tf. El cuerpo del desplazador se
considera un elemento poroso que hace las veces de regenerador, esta porosidad
puede contener una cantidad de gas que corresponde a un porcentaje del volumen
total del cuerpo del desplazador.
La fórmula para la variación de la presión con la posición angular de la manivela,
resultado del análisis de Schmidt3, indica que dicha dependencia se da a través de la
variación de los volúmenes de compresión Vcy de expansión Ve; luego estos
volúmenes pasan a considerarse como parámetros de diseño surgidos del análisis
termodinámico del motor Stirling tipo gamma.
Antes de entrar a detallar las formulaciones correspondientes que se obtienen
mediante el denominado diagrama de fasor, es necesario esclarecer los parámetros
geométricos y termodinámicos que se utilizarán para caracterizar el motor y su fluido
de trabajo (aire).
En estos parámetros se citamos las variables que tienen que ver con el análisis
termodinámico del motor y que son datos de entrada como son las temperaturas y las
dimensiones de los distintos elementos que conforman las zonas donde se realiza el
intercambio de calor como se muestra en la figura 21.
Dd = diámetro del elemento desplazador
Dp = diámetro del pistón en el lado frío
3 Sinusoidal
volume variations [en línea]. Ohio: Universidad de Ohio [consultado 20 de marzo de 2015].
Disponible en internet: http://www.ohio.edu/mechanical/stirling/engines/volumes.html
46
Ddt = diámetro del ducto que comunica los cilindros
Ld = longitud del desplazador
Lp = longitud del pistón
Figura 21. Parámetros geométricos y termodinámicos
Ldt = longitud del ducto que une los cilindros
Lcled = longitud del volumen muerto en la zona de expansión o lado caliente del
desplazador
Lclcd = longitud del volumen muerto en la zona de compresión al lado del
desplazadoro lado frío del desplazador
Lclcp = longitud del volumen muerto en la zona del pistón
Tf = temperatura del lado frío
47
Tc = temperatura del lado caliente
Tr = temperatura del regenerador que hace las veces de elemento desplazador
%Vd = porcentaje del volumen del desplazador ocupado por el regenerador
7.1 VARIACIÓN DE VOLÚMENES DE COMPRESIÓN Y EXPANSIÓN
Para analizar el motor Stirling tipo gamma es necesario hallar una expresión que
ayude a determinar la variación de sus volúmenes con respecto al ángulo ɵ, para lo
cual nos basaremos en un estudio realizado por Gustav Schmidt del Instituto Alemán
Politécnico de Praga en el año de 18714.
Dicho razonamiento parte de la suposición, de que, la temperatura en todo el
intercambiador de calor del lado caliente, se encuentra a una misma temperatura
superior y que todo el intercambiador de calor del lado frio se encuentra a una misma
temperatura inferior. También se supone que la masa de aire contenida en los
volúmenes es constante durante su funcionamiento. En esta configuración los
espacios de trabajo del motor en cuestión son los volúmenes donde se desplazan el
pistón, el desplazador y gas que en este caso es aire.
Estos espacios o volúmenes se muestran en la figura 22 y se nombran según las
características del aire contenido en él, así:
Vcled = volumen muerto en el lado caliente
Vda = mitad del volumen barrido por el desplazador en el lado caliente
Vclcd = volumen muerto en el lado frío correspondiente al desplazador
4
Ibíd., http://www.ohio.edu/mechanical/stirling/engines/volumes.html
48
Figura 22. Volúmenes de trabajo para el análisis de termodinámico
Vclcp = volumen muerto correspondiente al lado del pistón
Vpa = mitad del volumen barrido por el pistón
Xdmin = valor mínimo que alcanza el desplazador en la dirección X
Xdmax = valor máximo que alcanza el desplazador en la dirección X
Xpmin = valor mínimo que alcanza el pistón en la dirección X
Xpmax = valor máximo que alcanza el pistón en la dirección X
Xpa = longitud efectiva del volumen de comprimido por el pistón
𝑋𝑝𝑎 =
𝑋𝑝𝑚𝑎𝑥 −𝑋𝑝𝑚𝑖𝑛
2
−
𝐿𝑝
(7.1)
2
Xda = longitud efectiva del volumen recorrido por el desplazador en la zona de
expansión:
𝑋𝑑𝑎 =
𝑋𝑑𝑚𝑎𝑥 −𝑋𝑑𝑚𝑖𝑛
2
−
𝐿𝑑
(7.2)
2
49
Ap= área de la sección del pistón y es igual a:
𝜋
(7.3)
𝐴𝑝 = 4 𝐷𝑝2
Ad = área de la sección del desplazador y está dada por:
𝜋
(7.4)
𝐴𝑑 = 4 𝐷𝑑2
Adt = área de la sección del ducto que comunica los cilindros:
𝜋
(7.5)
2
𝐴𝑑𝑡 = 4 𝐷𝑑𝑡
El volumen barrido por el desplazador en la zona de expansión es:
(7.6)
𝑉𝑑𝑎 = 𝐴𝑑 𝐿𝑑𝑎
El volumen barrido por el pistón a partir de un punto medio hasta el límite superior de
la carrera es:
(7.7)
𝑉𝑝𝑎 = 𝐴𝑝 𝐿𝑝𝑎
Volumen muerto en el lado caliente o volumen muerto en la zona de expansión:
(7.8)
𝑉𝑐𝑙𝑒𝑑 = 𝐴𝑑 𝐿𝑐𝑙𝑒𝑑
Volumen muerto en el lado frio correspondiente al cilindro del desplazador:
(7.9)
𝑉𝑐𝑙𝑐𝑑 = 𝐴𝑑 𝐿𝑐𝑙𝑐𝑑
Volumen muerto en el lado del pistón:
(7.10)
𝑉𝑐𝑙𝑐𝑝 = 𝐴𝑝 𝐿𝑐𝑙𝑐𝑝
Volumen del ducto que comunica los cilindros
(7.11)
𝑉𝑑𝑡 = 𝐴𝑑𝑡 𝐿𝑑𝑡
Volumen total del cuerpo del desplazador
(7.12)
𝑉𝑑 = 𝐴𝑑 𝐿𝑑
Volumen del regenerador:
50
𝑉𝑟 =
(%𝑉𝑑 )𝑉𝑑
(7.13)
100
Vcoes es la suma de volúmenes de compresión comprendidos entre el lado del pistón
y del lado caliente en la parte inferior del desplazador, esto incluye los volúmenes
muertos.
(7.14)
𝑉𝑐𝑜 = 𝑉𝑑𝑎 + 𝑉𝑐𝑙𝑐𝑑 + 𝑉𝑑𝑡 + 𝑉𝑝𝑎 + 𝑉𝑐𝑙𝑐𝑝
Veo es la suma del volúmenes de expansión en la zona caliente incluyendo el volumen
muerto.
(7.15)
𝑉𝑒𝑜 = 𝑉𝑑𝑎 + 𝑉𝑐𝑙𝑒𝑑
Se designa φ como el ángulo de avance del desplazador con respecto al pistón. En la
figura 23 se presenta el diagrama de fasores para los vectores Xpa que representa el
desplazamiento del pistón y el vector Xda que representa el desplazamiento del
elemento desplazador y entre ellos el ángulo φ marcando el desfase entre los
vectores.
Figura 23. Diagrama de fasores para desplazamientos del pistón y desplazador
En la figura 24 se observa que el comportamiento volumen Vpa es inversamente
proporcional al desplazamiento del pistón; es decir que cuando el pistón se desplaza
en sentido positivo el Vpa disminuye. Este comportamiento ubica los vectores Xpa y
Vpa con un desfase de 180 grados.
De igual manera, se observa que el desplazamiento positivo del desplazador
proporciona un incremento del Vda en el lado frío lo que hace que tengan un
comportamiento lineal sin desfase con el vector Xda. Cuando el volumen Vda del lado
51
frío se incrementa, el volumen Vda del lado caliente disminuye por lo que existe un
desfase de 180 grados entre los vectores de Vda del lado frío y Vda del lado caliente.
Estas condiciones las podemos expresaren un diagrama de fasores como se muestra
en la figura 24, donde Vca es la suma de los volúmenes o espacios en los que se
desplazan el pistón y el desplazador, esto no incluye los volúmenes muertos. Este
volumen se encuentra desplazado con respecto al pistón un ángulo δ.
El volumen de comprensión Vc en un instante dado en función del ángulo 𝜃 está dado
por:
(7.16)
𝑉𝑐 = 𝑉𝑐𝑜 − 𝑉𝑝𝑎 𝑐𝑜𝑠𝜃 + 𝑉𝑑𝑎 𝑐𝑜𝑠(𝜃 + 𝜑)
Figura 24. Diagrama fasores para volúmenes de desplazamiento de pistón y
desplazador
Expandiendo la ecuación (7.16)
𝑉𝑐 = 𝑉𝑐𝑜 + (𝑉𝑑𝑎 𝑐𝑜𝑠𝜑 − 𝑉𝑝𝑎 )𝑐𝑜𝑠𝜃 − 𝑉𝑑𝑎 𝑠𝑒𝑛𝜑𝑠𝑒𝑛𝜃
(7.17)
Otra forma de expresar Vc es:
(7.18)
𝑉𝑐 = 𝑉𝑐𝑜 + 𝑉𝑐𝑎 𝑐𝑜𝑠(𝜃 + 𝛿)
Expandiendo (7.18) se tiene:
(7.19)
𝑉𝑐 = 𝑉𝑐𝑜 + 𝑉𝑐𝑎 (𝑐𝑜𝑠𝜃𝑐𝑜𝑠𝛿 − 𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝛿)
52
Igualando términos para las ecuaciones (7. 17) y (7.19) se tienen las siguientes
igualdades:
𝑉𝑐𝑎 𝑐𝑜𝑠𝛿 = 𝑉𝑑𝑎 𝑐𝑜𝑠𝜑 − 𝑉𝑝𝑎
(7.20)
𝑉𝑐𝑎 𝑠𝑒𝑛𝛿 = 𝑉𝑑𝑎 𝑠𝑒𝑛𝜑
(7.21)
Donde podemos despejar 𝛿:
𝛿 = 𝑇𝑎𝑛−1 (𝑉
𝑉𝑑𝑎 𝑠𝑒𝑛𝜑
𝑑𝑎 𝑐𝑜𝑠𝜑−𝑉𝑝𝑎
(7.22)
)
(7.23)
2 − 2𝑉 𝑉 𝑐𝑜𝑠𝜑 + 𝑉 2
𝑉𝑐𝑎 = √𝑉𝑝𝑎
𝑝𝑎 𝑑𝑎
𝑑𝑎
1⁄
2
2
2
− 2𝑉𝑝𝑎 𝑉𝑑𝑎 + 𝑉𝑑𝑎
𝑉𝑐 = (𝑉𝑐𝑙𝑐𝑑 + 𝑉𝑐𝑙𝑐𝑟 + 𝑉𝑐𝑙𝑐𝑝 ) + (𝑉𝑝𝑎
)
(1 + 𝑐𝑜𝑠(𝜃 + 𝛿))
(7.24)
El volumen de expansión Ve está dado por:
(7.25)
𝑉𝑒 = 𝑉𝑐𝑙𝑒𝑑 + 𝑉𝑑𝑎 (1 + 𝑐𝑜𝑠(𝜃 + 𝜑 + 𝜋))
7.2 EXPRESIÓN PARA LA PRESIÓN
La suposición de un sistema isotérmico hace que sea posible generar una expresión
sencilla para la presión del gas de trabajo en función de las variaciones de volumen.
De la ley de los gases ideales se tiene que:
𝑃𝑉
(7.26
𝑀 = 𝑅𝑇
Para cada uno de los volúmenes de trabajo se tiene
𝑃 𝑉𝑐𝑙𝑒𝑑
𝑀 = 𝑅(
𝑇𝑐
+
𝑉𝑑𝑎
𝑇𝑐
𝑃 𝑉𝑐𝑙𝑒𝑑 +𝑉𝑑𝑎
𝑀 = 𝑅(
𝑇𝑐
𝑉
+ 𝑇𝑟 +
𝑟
𝑉
+ 𝑇𝑟 +
𝑟
𝑉𝑑𝑎
𝑇𝑓
+
𝑉𝑐𝑙𝑐𝑑
𝑇𝑓
+
𝑉𝑐𝑙𝑐𝑝
𝑇𝑓
𝑉𝑑𝑎 +𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑝 +𝑉𝑝𝑎
𝑇𝑓
+
𝑉𝑝𝑎
𝑇𝑓
)
(7.27)
(7.28)
)
El regenerador se encuentra entre el lado caliente y el lado frío del motor, se supone
que éste tiene un comportamiento ideal y que tiene un perfil de temperatura lineal
53
entre la temperatura Tf frío y la temperatura Tc caliente. Por estas suposiciones la
media de la temperatura en el regenerador está dada por:
𝑇𝑐 −𝑇
(7.29)
𝑇𝑟 = 𝑙𝑛(𝑇 ⁄𝑇𝑓
𝑐
𝑓)
Se reemplaza (7.24) en (7.23)
𝑃
𝑉𝑐𝑙𝑒𝑑 +𝑉𝑑𝑎
𝑅
𝑇𝑐
𝑀= (
𝑉𝑟
𝑇𝑐 −𝑇𝑓
+
+
𝑉𝑑𝑎 +𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑝 +𝑉𝑝𝑎
𝑇𝑓
(7.30)
)
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
Despejando P:
𝑀𝑅
𝑃=
𝑉𝑐𝑙𝑒𝑑 +𝑉𝑑𝑎
+
𝑇𝑐
𝑉𝑟
𝑇𝑐 −𝑇𝑓
(7.31)
𝑉𝑑𝑎 +𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑝 +𝑉𝑝𝑎
𝑇𝑓
+
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
(
)
𝑀𝑅
𝑃=
𝑉𝑐
𝑇𝑐
(
+
𝑉𝑟
𝑇𝑐 −𝑇𝑓
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
(7.32)
𝑉
+ 𝑒
𝑇𝑓
)
Se reemplazan las ecuaciones (7.24) y (7.26) en (7.32) para obtener una expresión
que describe la presión en términos de la variación de los volúmenes.
𝑀𝑅
𝑃=
1⁄
2
2
(𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑟 +𝑉𝑐𝑙𝑐𝑝 )+(𝑉2
(1+𝑐𝑜𝑠(𝜃+𝛿))
𝑝𝑎 −2𝑉𝑝𝑎 𝑉𝑑𝑎 +𝑉𝑑𝑎 )
+
𝑇𝑐
(7.28)
𝑉𝑟
𝑇𝑐 −𝑇𝑓
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
(
54
𝑉
+𝑉 (1+𝑐𝑜𝑠(𝜃+∅+𝜋))
+ 𝑐𝑙𝑒 𝑑𝑎
𝑇𝑓
)
8. PARES DE TORSIÓN
8.1 PAR DE TORSIÓN DEL GAS
El torque de gas es el torque de salida (en el eje de la manivela) debido a la acción de
la presión del gas sobre el pistón; la fuerza del gas se debe a la presión que ejerce
este sobre el pistón en el proceso de expansión. Sea Fg la fuerza del gas, Pg la presión
del gas, Ap el área del pistón y Dp el diámetro interior del cilindro que en este caso es
igual al diámetro del pistón, por lo que el área del pistón es:
𝐴𝑔 =
𝜋
(8.1)
𝐷𝑝2
4
La fuerza ejercida sobre esta área es:
(8.2)
Fg = −Pg Ag
Al reemplazar el área se obtiene
𝜋
(8.3)
𝐹𝑔 = − 4 𝑃𝑔 𝐷𝑝2
Se reemplaza la ecuación (7.28) en (8.3) para obtener:
𝜋𝑀𝑅𝐷𝑝2
𝐹𝑔 = −
2
(𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑟 +𝑉𝑐𝑙𝑐𝑝 )+(𝑉2
𝑝𝑎 −2𝑉𝑝𝑎 𝑉𝑑𝑎 +𝑉𝑑𝑎 )
𝑇𝑐
4
1⁄
2
(1+𝑐𝑜𝑠(𝜃+𝛿))
(8.4)
𝑉
+𝑉 (1+𝑐𝑜𝑠(𝜃+∅+𝜋))
𝑉𝑟
+ 𝑇𝑐−𝑇
+ 𝑐𝑙𝑒 𝑑𝑎
𝑇𝑓
𝑓
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
(
)
La fuerza del gas distribuida sobre la superficie del pistón se convierte en una sola
fuerza Fg que actúa a través del centro de masa del pistón, el signo negativo se debe
a la orientación de la fuerza contraria a la dirección positiva de x en la figura 25.
Este par de torsión se debe a la fuerza del gas Fg que actúa en un brazo de momento
con respecto al centro de la manivela O2. El brazo varía desde cero
55
Figura 25. Presión y fuerza del gas sobre el pistón
Hasta un máximo conforme voltea la manivela. Un análisis al sistema de fuerzas que
concurren en el punto B de la figura 26 determina que:
𝐹𝑔14 = 𝐹𝑔 𝑡𝑎𝑛∅
(8.5)
𝐹𝑔34 = −𝐹𝑔 𝒊 − 𝐹𝑔 𝑡𝑎𝑛∅𝒋
(8.6)
Figura 26. Diagrama cuerpo libre análisis de fuerza del gas mecanismo
manivela-corredera
Fuente: NORTON, Robert. Diseño de maquinaria. Cuarta edición México: Prentice
Hall, 2011. p. 523
56
El par de torsión motriz Tg21 generado por la fuerza del gas puede determinarse a partir
del producto de la fuerza Fg14 o Fg41 y la distancia x que es el brazo de momento
instantáneo con respecto a O2.
Tg21 = Fg41 . xk
(8.6)
Al sustituir la ecuación (6.17) y (8.2) en la ecuación (8.3) se tiene que:
𝑟2
𝑟2
ℎ2
(8.7)
𝑇𝑔21 = (𝐹𝑔 𝑡𝑎𝑛∅) (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 − 4𝑙 + 4𝑙 𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − 2𝑙 )
Es posible expresar 𝑡𝑎𝑛∅ en términos de wt
𝑡𝑎𝑛∅ =
𝑞−ℎ
𝑢
=
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
(8.8)
𝑙𝑐𝑜𝑠∅
Reemplazando (6.7) en la ecuación (8.8) se tiene:
𝑡𝑎𝑛∅ =
𝑞−ℎ
𝑢
=
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
2
(8.10)
𝑟𝑠𝑒𝑛𝑤𝑡−ℎ 2
) )
𝑙
𝑙( √1−(
El radical de la ecuación (8.10) se reemplaza por la expresión (6.14)
𝑡𝑎𝑛∅ =
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
(8.10)
1 𝑟2
2𝑟
ℎ2
𝑙(1− ( 2 𝑠𝑒𝑛2 𝑤𝑡− 𝑠𝑒𝑛𝑤𝑡+ 2 ))
2 𝑙
𝑙
𝑙
Desarrollando el denominador
𝑡𝑎𝑛∅ =
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
(8.11)
𝑟2
ℎ2
𝑙− 𝑠𝑒𝑛2 𝑤𝑡+𝑟𝑠𝑒𝑛𝑤𝑡+
2𝑙
𝑙
Reemplazando (8.11) en (8.7) se obtiene: la ecuación para el par de torsión motriz
Tg21 en términos de wt
𝑇𝑔21 = 𝐹𝑔 (
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
𝑟2
ℎ2
𝑙− 𝑠𝑒𝑛2 𝑤𝑡+ 𝑟𝑠𝑒𝑛𝑤𝑡−
2𝑙
2𝑙
𝑟2
𝑟2
ℎ2
) (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 − 4𝑙 + 4𝑙 𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − 2𝑙 )
(8.12)
En la ecuación (8.12) se reemplaza la expresión para la fuerza (8.4) y se obtiene la
ecuación para el par de torsión motriz debido al gas Tg21 en términos de wt
57
𝜋𝑀𝑅𝐷𝑝2
𝑇𝑔21 = −
4
(
( (
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
𝑟2
2𝑙
∗
1⁄
2
2
(𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑟 +𝑉𝑐𝑙𝑐𝑝 )+(𝑉2
(1+𝑐𝑜𝑠(𝜃+𝛿))
𝑝𝑎 −2𝑉𝑝𝑎 𝑉𝑑𝑎 +𝑉𝑑𝑎 )
+
𝑇𝑐
𝑙− 𝑠𝑒𝑛2 𝑤𝑡+ 𝑟𝑠𝑒𝑛𝑤𝑡−
𝑉
+𝑉 (1+𝑐𝑜𝑠(𝜃+∅+𝜋))
+ 𝑐𝑙𝑒 𝑑𝑎
𝑉𝑟
𝑇𝑐 −𝑇𝑓
𝑇𝑓
))
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
ℎ2
2𝑙
) (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 −
𝑟2
4𝑙
+
𝑟2
𝑐𝑜𝑠2𝑤𝑡
4𝑙
+ 𝑟𝑠𝑒𝑛𝑤𝑡 −
ℎ2
)
2𝑙
(8.13)
8.2 PARES DE TORSIÓN DE INERCIA Y SACUDIMIENTO
El torque de inercia también llamado de sacudimiento es el torque debido a las fuerzas
inerciales de los eslabones en movimiento. Su contribución neta al movimiento del
mecanismo es cero, pero ayudan a que el mecanismo pueda pasar a través de sus
puntos muertos (aquellas configuraciones en que el ángulo de transmisión se acerca
o es igual a 0° o 180°). La expresión para el torque de inercia se aplica dos veces, ya
que tanto el pistón como el desplazador juegan un papel en este sentido.
En el punto A la componente tangencial tiene un brazo de momento de radio de
manivela r, Si su w es constante la masa en A no aporta el par de torsión de inercia.
El pistón tiene una componente de inercia en la dirección x, excepto cuando este se
encuentra en PMS o PMI. De la figura 20 el par de torsión generado por esta fuerza
puede expresarse como:
𝐹𝑖41
⟹ 𝐹𝑖41 = 𝑡𝑔𝜙𝑚𝐵 𝑥̈
(8.14)
𝑇𝑖21 = (𝐹𝑖41 . 𝑥)𝒌 = (−𝐹𝑖14 . 𝑥)𝒌
(8.15)
𝑡𝑔∅ =
𝑚𝐵 𝑥̈
En (8.15) se sustituye la expresión (6.17) para x
𝑟2
𝑟2
ℎ2
𝑇𝑖21 = −(−𝑡𝑔𝜙𝑚𝐵 𝑥̈ ) (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 − 4𝑙 + 4𝑙 𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − 2𝑙 ) 𝒌
Al reemplazar la ecuación (8.11) de 𝑡𝑔𝜙 en (8.16) se tiene:
58
(8.16)
𝑇𝑖21 = −(−𝑚𝐵 𝑥̈ ) [
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
ℎ2
𝑟2
𝑙− 𝑠𝑒𝑛2 𝑤𝑡+𝑟𝑠𝑒𝑛𝑤𝑡+
𝑙
2𝑙
(𝑟𝑠𝑒𝑛𝑤𝑡)−ℎ
𝑇𝑖21 = 𝑚𝐵 [
𝑟2
ℎ2
𝑙− 𝑠𝑒𝑛2 𝑤𝑡+𝑟𝑠𝑒𝑛𝑤𝑡+
2𝑙
𝑙
𝑟2𝑤2
(−𝑟𝑤 2 𝑐𝑜𝑠𝑤𝑡 −
𝑙
] (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 −
𝑟2
4𝑙
+
𝑟2
𝑟2
4𝑙
𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 −
ℎ2
2𝑙
𝑟2
(8.17)
)𝒌
ℎ2
] (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 − 4𝑙 + 4𝑙 𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − 2𝑙 ) ∗
(8.18)
𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑤 2 𝑠𝑒𝑛𝑤𝑡) 𝒌
8.3 TORQUE DE LA FUERZA DEL GAS Y EFECTO DEL DESCENTRADO
La fuerza del gas que es perpendicular al área del pistón y paralela al eje de la
corredera, actúa con el brazo h para generar un torque respecto al eje de la manivela
como se aprecia en la figura (25), por lo que el torque está dado por:
(8.19)
𝑇𝑔ℎ = 𝐹𝑔 ℎ
Reemplazando en (8.19) la expresión para Fg de (8.4) se tiene:
ℎ𝜋𝑀𝑅𝐷𝑝2
𝑇𝑔ℎ = −
4
1⁄
2
2
(𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑟 +𝑉𝑐𝑙𝑐𝑝 )+(𝑉2
(1+𝑐𝑜𝑠(𝜃+𝛿))
𝑝𝑎 −2𝑉𝑝𝑎 𝑉𝑑𝑎 +𝑉𝑑𝑎 )
+
𝑇𝑐
𝑉𝑟
𝑇𝑐 −𝑇𝑓
𝑉
+𝑉 (1+𝑐𝑜𝑠(𝜃+∅+𝜋))
+ 𝑐𝑙𝑒 𝑑𝑎
𝑇𝑓
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
(
)
(8.20)
8.4 TORQUETOTAL DEL MOTOR
El par de torsión absoluto que se genera en el eje de la manivela es la suma de los
pares debido a las fuerzas de inercia y los pares generados por la fuerza del gas:
(8.21)
𝑇𝑡 = 𝑇𝑔 + 𝑇𝑖 + 𝑇𝑔ℎ
Donde:
59
ℎ𝜋𝑀𝑅𝐷𝑝2
𝑻𝒈𝒉 = −
1
4(
2 −2𝑉 𝑉 +𝑉 2 ) ⁄2 (1+𝑐𝑜𝑠(𝜃+𝛿))
(𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑟 +𝑉𝑐𝑙𝑐𝑝 )+(𝑉𝑝𝑎
𝑝𝑎 𝑑𝑎
𝑑𝑎
𝑇𝑐
+
𝑉𝑟
𝑇𝑐 −𝑇𝑓
+
𝑉𝑐𝑙𝑒 +𝑉𝑑𝑎 (1+𝑐𝑜𝑠(𝜃+∅+𝜋))
𝑇𝑓
)
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
𝑻𝒊𝟐𝟏 = 𝑚𝐵 [
(𝑟𝑠𝑒𝑛𝑤𝑡) − ℎ
] (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 −
ℎ2
𝑟2
𝑟2 𝑟2
ℎ2
+ 𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − )
4𝑙 4𝑙
2𝑙
𝑙 − 2𝑙 𝑠𝑒𝑛2 𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 +
𝑙
2 2
𝑟
𝑤
𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑤 2 𝑠𝑒𝑛𝑤𝑡) 𝒌
∗ (−𝑟𝑤 2 𝑐𝑜𝑠𝑤𝑡 −
𝑙
𝜋𝑀𝑅𝐷𝑝2
𝑻𝒈 = −
1
4(
2 −2𝑉 𝑉 +𝑉 2 ) ⁄2 (1+𝑐𝑜𝑠(𝜃+𝛿))
(𝑉𝑐𝑙𝑐𝑑 +𝑉𝑐𝑙𝑐𝑟 +𝑉𝑐𝑙𝑐𝑝 )+(𝑉𝑝𝑎
𝑝𝑎 𝑑𝑎
𝑑𝑎
𝑇𝑐
(
+
𝑉𝑟
𝑇𝑐 −𝑇𝑓
𝑙𝑛(𝑇𝑐 ⁄𝑇𝑓 )
∗(
(𝑟𝑠𝑒𝑛𝑤𝑡) − ℎ
𝑟2
) ∗ (𝑟𝑐𝑜𝑠𝑤𝑡 + 𝑙 −
ℎ2
𝑙 − 2𝑙 𝑠𝑒𝑛2 𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − 2𝑙
60
+
𝑉𝑐𝑙𝑒 +𝑉𝑑𝑎 (1+𝑐𝑜𝑠(𝜃+∅+𝜋))
𝑇𝑓
)
)
𝑟2 𝑟2
ℎ2
+ 𝑐𝑜𝑠2𝑤𝑡 + 𝑟𝑠𝑒𝑛𝑤𝑡 − )
4𝑙 4𝑙
2𝑙
9. IMPLEMENTACIÓN HOJA DE CÁLCULO
En la preparación de este análisis se han logrado expresiones con la finalidad de
obtener un modelo matemático que describa aproximada y teóricamente el
comportamiento de un motor Stirling de configuración gamma de baja potencia. Las
ecuaciones que conforman este modelo están encaminadas a hallar la presión del gas
(7.28)y el torque total que se ejerce sobre el eje de la manivela (8.1), por lo tanto estos
serían el punto de análisis y de comparación para las diferentes arquitecturas que
generan cada posible conjunto de variables.
Dado que no es fácil fabricar una serie de motores Stirling con diferentes
configuraciones para su respectivo análisis, se advierte que este modelo matemático
ayuda a estudiar el comportamiento o las tendencias de las variaciones de los
parámetros del motor. Montar y desarrollar estas ecuaciones en una hoja de cálculo
como Microsoft office Excel 2007 permite tener una herramienta que se convierte en
un simulador matemático con el cual es posible analizar y optimizar la presión del gas
y torque total de un motor Stirling tipo gamma.
En esta hoja de cálculo las celdas donde se digita el valor para las variables de entrada
tienen fondo amarillo y se emplea el sistema internacional de unidades SI. En la figura
27 se ilustra la pantalla principal de esta hoja de cálculo.
9.1 VARIABLES
Las variables creadas para esta hoja de cálculo son:
r = longitud de la manivela
Lb = longitud de la biela
h = distancia entre el eje de la corredera y el eje de la manivela.
Dd = es el diámetro del elemento desplazador
Dp = es el diámetro del pistón en el lado frío
Ddt = es el diámetro del ducto que comunica los cilindros
61
Figura 27. Pantalla principal hoja de calculo
Ld = longitud del desplazador
Lp = longitud del pistón
Ldt = longitud del ducto que une los cilindros
Lclde = longitud del volumen muerto en la zona de expansión o lado caliente del
desplazador
Lcldc = longitud del volumen muerto en la zona de compresión al lado del
desplazadoro lado frío del desplazador
Lclp = longitud del volumen muerto en la zona del pistón
Lclp = longitud zona muerta lado del pistón
Tf = temperatura del lado frío
Tc = temperatura del lado caliente
Tr = temperatura del gas en el regenerador.
62
%Vd = porcentaje del volumen del desplazador por el gas en el regenerador
R = constante de los gases ideales
m2a = sección de masa equivalente de la manivela-pistón en el punto A
m3a = sección de masa equivalente de la biela-pistón en el punto A
m3b =sección de masa equivalente de la biela-pistón en el punto B
m4 =masa del pistón en el punto B
m5a’= sección de masa equivalente de la manivela-desplazador en el punto A
m6a’ = sección de masa equivalente de biela - desplazador en el punto A
m6b’ = sección de masa equivalente de biela – desplazador en el punto B
m7 = masa del desplazador
φ = ángulo de desfase del desplazador
Ƿ = densidad del gas de trabajo
9.2 DESCRIPCIÓN DE HOJA DE CÁLCULO
La hoja de cálculo consta de los siguientes bloques:
 Variables dimensionales de los mecanismos manivela biela. Se recuerda que estas
dimensiones son idénticas para el pistón y el desplazador.
 Variables dimensionales del pistón, desplazador, y volúmenes muertos.
 Propiedades del gas de trabajo.
 Masas equivalentes para los mecanismos manivela biela del pistón y desplazador.
 Alertas por inconsistencias en valores asignados.
63
 Curva de presión del gas de trabajo y del torque toral.
 Cálculos generales.
Variables dimensionales de los mecanismos manivela biela. La longitud de la
manivela se digita en la celda C7, la altura de descentrado h en la celda C8 y la
longitud de la biela en C9. La celda E7 contiene la velocidad angular y la casilla E8
contiene el ángulo de desplazamiento entre la manivela del pistón y la manivela del
desplazador. En la celda E9 se realiza el chequeo de la restricción para las valores de
r, L y h. Ver figura 28. En este bloque no es correcto digitar valores negativos.
En la figura 28 se ilustra la correspondencia de estas variables en la hoja de cálculo
con el modelo del mecanismo manivela biela de pistón y desplazador. Es de notar que
las variables r, h y L solo se solicitan una vez dada la suposición de que son idénticas
para el mecanismo del pistón y del desplazador.
Figura 28 Bloque de variable dimensional para los mecanismos manivela biela
del pistón y desplazador.
64
Figura 29. Bloque de variable dimensional para los mecanismos manivela biela
del pistón y desplazador
Variables dimensionales del pistón, desplazador, y volúmenes muertos. En este
bloque se solicitan las dimensiones del pistón y desplazador, de los volúmenes
muertos y del elemento que los comunica. En la figura 30 se presenta una ilustración
de este bloque.
Figura 30. Bloque correspondiente a las dimensiones del pistón, desplazador y sus
volúmenes muertos.
En H7 y H8 se digitan los valores para la longitud del pistón y del desplazador. J7 y J8
contienen los diámetros del pistón y del desplazador. Las celdas H9 y J9 son la
longitud y diámetro respectivamente del conducto que comunica el volumen de
compresión del lado desplazador con el lado del pistón. L7 contiene la longitud del
volumen muerto de compresión del lado pistón, L8 la longitud del volumen muerto de
compresión lado desplazador y L9 es la longitud del volumen muerto del expansión
65
lado desplazador. La figura 31 ilustraciones la asignación de variables en celdas de
la hoja de cálculo según corresponde en los volúmenes del modelo planteado para el
motor Stirling analizado en el presente informe. En este bloque no es correcto digitar
valores negativos.
Figura 31 Asignación de variables en celdas de hoja de cálculo para espacios de pistón
y desplazador.
Propiedades del gas de trabajo. En este bloque se consignan las propiedades del
gas de trabajo. La celda O7 contiene la temperatura del lado frio que corresponde al
pistón y O8 la temperatura del lado caliente relacionada al desplazador. Q7 es la
densidad del gas de trabajo y Q8 es la constante para los gases ideales. O9 almacena
el porcentaje de gas que puede albergar el cuerpo del desplazador que a su vez es el
regenerador. Ver figura 32.
66
Figura 32. Bloque de propiedades del gas de trabajo
Masas equivalentes. En este bloque se introducen las masas equivalentes para los
mecanismos manivela biela del pistón y desplazador. En la hoja de cálculo aparece
como de muestra en la figura 33.
Figura 33. Bloque para variables de propiedades del gas de trabajo
La figura 33 presenta un esquema de la ubicación de las componentes de las masas
equivalentes para los puntos A, A’, B, B’ que son los puntos tenido en cuenta en esta
hoja de cálculo.
m2A masa de la manivela 2 del pistón concentrado en el punto A.
m3A corresponde a la masa de biela 3 concentrada en el punto A.
m3B masa de biela concentrada en el punto B equivalente al pistón.
m4 masa de pistón concentrada en el eje de B.
m5A’ masa de la manivela de desplazador concentrada en el pasador A’
67
m6A’ masa de la biela del desplazador concentrada en el pasador A’
m6B’ masa de la biela del desplazador concentrada en el pasador B’
m7 masa del desplazador concentrada en el eje 7.
Figura 34. Ubicación de las componentes de las masas equivalentes para los
puntos A, A’, B, B’ en el modelo del mecanismo manivela pistón.
Alertas. Estas son elementos que se activan en el momento que se introducen
algunos valores inconsistentes y se manifiestan colocando en rojo la casilla en la que
se digito el valor errado, además se presenta un rótulo pidiendo verificar el error
cometido. Las alertas que activa esta hoja de cálculo son:
Verificar que (r+h)/L< 1
Que el cociente (r+h)/L sea menor que uno es una condición necesaria para la
solución de la ecuación (6.8). Esta alarma advierte que se han digitado una
configuración que no cumple esta condición y el programa detiene los cálculos. Se
debe corregir esta situación cambiando valores de r, h y L hasta que se cumpla la
condición (r+h)/L< 1
Verificar que Tf < Tc
68
Esta alerta verifica que la temperatura relacionada con el desplazador sea siempre la
caliente o la mayor de las dos y que la temperatura relacionada con el pistón sea
siempre la menor o la temperatura fría.
Curva de presión del gas de trabajo y del torque total. La hoja de cálculo presenta
la curva del comportamiento de la presión del gas de trabajo con respecto a la posición
de la manivela en la gráfica izquierda de la pantalla y el comportamiento del torque
total del motor durante una revolución en el eje de la manivela en la gráfica derecha
de la pantalla.
69
10. ANÁLISIS DE RESULTADOS
Variando los parámetros geométricos y termodinámicos es posible obtener una
variedad de configuraciones para el motor Stirling estudiado en el presente informe.
La hoja de cálculo permite analizar y comparar los efectos que causan en la presión
del gas y el torque total para cada una de estas arquitecturas del motor Stirling.
Manteniendo constantes las variables termodinámicas se analiza el comportamiento
de la presión del gas y del torque según se varía la altura h. La temperatura del lado
frío permanece en 28oC y la temperatura del lado caliente en 200oC, la densidad del
aire está en 1,184 kg/cm3. Los valores para las masas equivalentes y los espesores
del pistón y desplazador son como aparecen en la figura 35.
Figura 35. Valores iníciales para análisis de configuraciones del motor
Primera configuración h = 0
Este arreglo de motor Stirling está caracterizado por las variables que se presentan
en la figura 35. La presión del gas varia aproximadamente entre 114100 y 112800
Pas, ver figura 36. Esta presión alcanza su máximo valor a los 275 grados de giro de
la manivela aproximadamente.
El torque para esta primera configuración toma valores entre -0,8 y 0,7 Nm y su
máximo valor se alcanza a los 90 grados de giro para la manivela como se muestra
en la figura 37.
70
Figura 36. Presión del gas para la primera configuración del motor con h = 0 m.
Figura 37. Valores para el torque en la primera configuración del motor con h =
0
Segunda configuración. h < r
Sin variar los parámetros termodinámicos solamente se modifica la longitud de h a
0,01m que es la segunda configuración, los otros valores se muestran en la figura 38.
Para esta segunda disposición la curva de presión se muestra en la figura 39 donde
se nota que la presión del gas de trabajo no sufre cambios notables en comparación
con la anterior configuración.
71
Figura 38. Valores para la segunda configuración h=0,01 m
Los resultados para valores del torque se presentan el al figura 40 y se observa que
estos oscilan entre -0,8 y 0,7 Nm aproximadamente. El valor máximo del torque se a
los 90 grados de giro de la manivela. En comparación con la disposición anterior, se
nota que al incrementar la longitud de la altura h de 0 m a 0,01m la presión en el gas
de trabajo y el torque permanece relativamente constante.
Figura 39. Curva del gas de trabajo para la segunda disposición h = 0,01 m.
72
Figura 40. Curva del torque para segunda disposición del motor con h = 0,01 m
Tercera configuración. h = r
La tercer configuración del motor Stirling se logra haciendo la longitud de h = 0,02m.
Para esta disposición la curva de presión se presenta en la gráfica 10.7 y no se
observa mayor variación en estos valores. Los valores del torque son muy similares
oscilando entre -0,7 Nm y 0,7 Nm, la diferencia es que el área bajo la curva de valores
positivos de toque se hace mayor, se puede concluir que en esta configuración el
motor puede funcionar con mayor facilidad al lograr producir un torque mayor que el
consumido para su funcionamiento.
Figura 41. Comportamiento de la presión para h= 0,02m
73
Figura 42. Comportamiento del torque para h= 0,02m
Cuarta configuración. h > r
Una cuarta configuración del motor Stirling se establece con h=0,04m sin variar los
otros parámetros geométricos. La figura 43 enseña la curva de presión para esta
configuración donde se advierte un comportamiento muy similar a las anteriores
comparaciones, lo que lleva a concluir que la presión del gas no es sensible a los
cambios en la longitud de la altura h. La figura 44 exhibe la curva del toque total para
esta configuración, los valores oscila entre -0,2 Nm y 0,9Nm y el punto máximo lo
alcanza alrededor de los 110 grados de giro de la manivela. En esta curva de torque
total el área bajo la curva de valores positivos se incrementó con respecto a la
configuración anterior.
Figura 43. Comportamiento de la presión para h=0,04 m
74
Figura 44. Comportamiento del torque total para h=0,04 m
Quinta configuración. h > r
Una quinta configuración se logra con h=0,055 m sin variar los demás parámetros
geométricos como se muestra en la figura 45 y en la figura 46 se enseña la curva de
presión en la que se observa un comportamiento muy similar a las anteriores
configuraciones. La figura 47 exhibe la curva del toque total para esta disposición, los
valores oscila entre 0,5 Nm y 1,3 Nm y el punto máximo lo alcanza entre los 110 y 200
grados de giro de la manivela. En esta curva no existen valores negativos por lo que
el área bajo la curva se incrementa en gran proporción por lo tanto en este arreglo el
motor puede brindar un mayor torque de trabajo. Los parámetros del mecanismo
manivela corredera para esta configuración que hasta el momento ha sido la de
mejores resultados en la generación de torque.
Figura 45. Parámetros geométrico para la quita configuración
75
Figura 46. Comportamiento de la presión para h=0,055 m
Figura 47. Comportamiento del torque total para h=0,055 m
Sexta configuración biela = 0,14 m
Si la biela se hace menor a 0,08 m no cumple la restricción de la ecuación (6.1) Para
conformar una nueva disposición se le da un valor de 0,14 m a la biela. Las curvas de
presión y torque para este arreglo se muestran en las figuras 10.14, y 10.15
respectivamente. En la curva de torque de esta configuración se presentan valores
negativos de hasta -0,2 Nm aproximadamente el área bajo la curva de valores
positivos se hace más pequeña, lo que indica que una biela muy larga no beneficia el
rendimiento del motor. Hasta el momento se observa que la presión no es función
directa de las longitudes del mecanismo manivela corredera.
76
Figura 48. Curva de presión para el sexto arreglo con biela = 0,14
Figura 49. Comportamiento del torque total para h=0,14 m
Séptima configuración Diámetro desplazador > diámetro pistón
A partir de los valores de la quinta configuración que ha sido la mejor de las seis
anteriores se procede a varia los diámetros de pistón y desplazador. El séptimo arreglo
a comparar será con un desplazador de 0,1 m al diámetro mayor que el diámetro del
pistón, las variables restantes se muestran en la figura 50. Las curvas de presión y el
torque se comportan como lo muestran las figuras 51 y 52 respectivamente. El área
bajo estas dos curvas disminuye levemente, esto indica que un desplazador de mayor
diámetro que el pistón no mejora la eficiencia del motor Stirling estudiado en el
presente informe
77
Figura 50. Variables para séptimo arreglo a comparar diámetro desplazador =
0,1 m
Figura 51. Curva de presión séptimo arreglo diámetro desplazador = 0,1 m
Figura 52. Curva de torque séptimo arreglo diámetro desplazador = 0,1 m
78
Octava configuración Diámetro desplazador < diámetro pistón
Un octavo arreglo a comparar será con un diámetro de desplazador menor que el
diámetro del pistón con los parámetros expresados en la figura 53. Las curvas de
presión y el torque se comportan como lo muestran las figuras 54 y 55
respectivamente. El área bajo estas dos curvas sufre un notable incremento, esto
indica que un diámetro en el desplazador menor que el diámetro del pistón
incrementa la eficiencia del motor Stirling estudiado en el presente informe.
Figura 53. Variables para el octavo arreglo del motor Stirling a comparar
Figura 54. Curva de presión para un desplazador < pistón
Dado que al presentarse un incremento notable en los valores de presión hasta
138000 Pa, un prototipo construido bajo estas especificaciones debe tener un sistema
de sellado más eficiente para evitar fugas del gas de trabajo.
79
Figura 55. Curva del torque para un desplazador < pistón
Novena configuración. Densidad aire < 1,184
La novena configuración a revisarse establece con una densidad de aire de menor
valor, igual a 0,8 Kg/m3 y con los parámetros expresados en la figura 56. Las curvas
de presión y de torque se muestran las figuras 57 y 10.24correspondientemente.
El área bajo estas dos curvas sufre una notable reducción, esto indica que una
densidad menor a 1,184 Kg/m3disminuye la eficiencia del motor Stirling estudiado
en el presente informe.
Figura 56. Valores que caracterizan la novena configuración del motor.
80
Figura 57. Curva de presión para la novena configuración
Figura 58. Curva de torque total para la novena configuración
Decima configuración. Densidad del aire > 1,184
La décima configuración a revisar se establece con una densidad de aire de mayor
valor igual a 1,7 Kg/m3 y con los parámetros expresados en la figura 59. Las curvas
de presión y de torque se muestran las figuras 60 y 61 correspondientemente. El área
bajo estas dos curvas se incrementa claramente, esto indica que una densidad mayor
en el gas de trabajo mejora la eficiencia del motor Stirling.
81
Figura 59. Parámetros que caracterizan la décima configuración
Figura 60. Curva de presión para densidad del aire = 1,7 Kg/m3
Figura 61. Curva de torque para densidad del aire = 1,7 Kg/m3
82
11. CONCLUSIONES
Existe un interés en la industria y la academia en desarrollar programas de cara al
diseño y experimentación de motores Stirling que compitan con las soluciones
actuales a la crisis energética. Tener un prototipo respaldado por un análisis mecánico
y termodinámico es determinante para participar de manera activa en las discusiones
sobre el tema
En este trabajo, la formulación del modelo termodinámico-mecánico de un prototipo
de motor Stirling tipo gamma inició por la parte mecánica valiéndose del análisis
dinámico de un mecanismo manivela-deslizador. Se partió de un modelo cinemático,
en el cual se establecen expresiones para la posición, velocidad y aceleración del
deslizador dada una posición angular de la manivela. La consideración del desfase u
“Offset”, que es la longitud del segmento que parte del eje de la manivela y es
perpendicular a la línea de eje del deslizador, cuando se mira al mecanismo
perpendicular al plano descrito por la rotación de la manivela, obligó a una
reformulación propia de la simplificación que por medio de la teoría del binomio se
aplica a la expresión para la posición del deslizador.
El análisis de fuerzas se realizó sobre la base del concepto de masas puntuales para
representar el mecanismo. Se produjeron expresiones para el torque de gas y el
torque inercial o dinámico, dejando como entrada una expresión para la presión del
gas, que se obtuvo como resultado de un análisis termodinámico. Dicho análisis es el
que se presenta bajo el nombre de Gustav Schmidt, donde se asume un régimen
isotérmico para el gas en los espacios de compresión y expansión tanto barridos por
el desplazador y el pistón, como los correspondientes a volúmenes muertos. Así
mismo, se asume un perfil lineal de temperatura a través del regenerador, o zona que
media entre los volúmenes de aire frío y caliente. Se produce como resultado un perfil
de presión de trabajo en función de la coordenada angular de la manivela, que resulta
ser altamente dependiente de la fase que se tenga entre el movimiento del pistón y el
del desplazador, que a su vez va a incidir en la variación periódica del volumen de
compresión.
Se utilizó el modelo implementado en hoja de cálculo de Excel para examinar el efecto
en las curvas de fuerza de gas y de torque total producido por las siguientes variables:
magnitudes de manivela biela y Offset, diámetro del desplazador, diámetro del pistón,
longitud del desplazador, longitud del pistón, volumen del lado caliente, volumen del
lado frio, temperatura del lado caliente, temperatura del lado frio, densidad del gas de
trabajo. Se encontró que los factores de mayor incidencia fueron:
- La relación entre el tamaño de la manivela y la altura de descentrado offset, con
una manivela de menor tamaño que el offset se mejora el toque de salida del motor.
83
- Un pistón de mayor diámetro que el desplazador mejora la eficiencia del motor.
- El toque total y la presión del gas son directamente proporcionales a la densidad
del gas de trabajo.
También se observó que temas como la longitud de las bielas, los volúmenes muertos
de gas tuvieron poco efecto en los datos de salida.
El modelo propuesto no considera procesos transigentes de transferencia de calor así
como fenómenos de flujo del gas al interior de la cámara, perdidas por fricción,
balanceo de volantes entre otros. Para llegar a conclusiones más aproximada sal
comportamiento real del motor Stirling analizado en el presente trabajo se hace
necesario incorporarle al modelo matemático obtenido ecuaciones que procesen estos
temas.
Los resultados producidos sugieren que sería interesante realizar un montaje donde
se controlen las variables como el offset, diámetros de pistón y desplazador pues en
el modelo matemático resultan ser determinantes en la producción de torque de salida
en el motor Stirling considerado en este trabajo.
84
BIBLIOGRAFÍA
DEUTSCHMAND. Aaron. MICHELS Walter, WILSON Ch. Diseño de máquinas teoría
y práctica. Editorial continental S.A. México. Cuarta impresión 1991. 157 p.
BOLES, Michael A. Termodinámica. Cuarta edición: Mc Graw Hill 2003. 117 p.
CALLEJÓN AGRAMUNT, Ismael. Máquinas térmicas motoras. México: Alfaomega,
2005. 533 p.
CHIN HSIANG Cheng, YING JU Yu. Combining dynamic and thermodynamic models
for dynamic simulation of a beta type stirling engine En: Renewable energy. 2012. No.
37, p. 161-173.
CONFERENCIA LATINO AMERICANA DE ENERGÍA SOLAR. (4: 2010: Cusco)
ISES_CLA. Chile: Universidad de Chile, Departamento de Ingeniería Mecánica, 2010.
HAMILTON H. Mabie y OCVIRK, Fred. Mecanismos y dinámica de maquinaria.
segunda edición México: Editorial Limusa,1981 530 p.
HIN HSIANG Cheng, YING JU Yu. Dynamic simulation of a beta type stirling engine
with cam drive mechanism via the combination of the thermodynamic. En: Renewable
energy. 2011. No. 36, p. 714-725.
MIRANDA MORALES José. Análisis termodinámico de un motor de ciclo Stirling.
Tesis para obtener el grado de Maestro en Ciencias. México: Instituto Politécnico
Nacional, 2008. 91 p.
MORAN, Michael J. Fundamentals of engineering thermodynamics. New York: John
Wiley sons, 1996. 859 p.
NORTON, Robert. Diseño de maquinaria. Cuarta edición México: Prentice Hall, 2011.
856 p.
PUECH, Pascal y TISHKOVA, Victoria. Thermodynamic analysis of a stirling engine
including regenerator dead volume. En: Energy Renewable. 2011. No. 36 p. 872-878.
85
RODGER W. Dyson y SCOTT D., Wilson. Fast whole-engine stirling analysis. En:
Nasa/TM-2005-213960, AIAA Paper 2005-5558, E-15276. 2005. 31 p.
ROGDAKISE, D. Bormpilas. Athermodynamic study for the optimization of stable
operation of free piston stirling engines. En: Energy conversion and management.
2004. No. 45, p. 575-593.
ROBSON A, GRASSIE T y KUBIE J. Modelling of a low temperature differential stirling
engine. En: Mechanical engineering science proc. Imeche. 2008. Vol 22, no. 8. p. 927.
SAAD, Michael A. Thermodynamics principles and practice. New Jersey: PrenticeHall, 1997. 935p.
Stirling engines - mechanical configurations [en línea]. Ohio: Universidad de Ohio
[consultado
20
de
marzo
de
2015].
Disponible
en
Internet:
http://www.ohio.edu/mechanical/stirling/engines/engines.html
Sinusoidal volume variations [en línea]. ohio: universidad de ohio [consultado 20 de
marzo
de
2015].
Disponible
en
internet:
http://www.ohio.edu/mechanical/stirling/engines/volumes.html
86