Diseño de Turbina Eólica de Baja Velocidad
Diseño de Turbina Eólica de Baja Velocidad
1
2
AGRADECIMIENTOS
Por último agradezco a mis padres y a mi hermano por ayudarme y apoyarme durante
todo el proyecto. Gracias a ellos soy el ser humano que soy.
3
ÍNDICE
AGRADECIMIENTOS .................................................................................................... 3
LISTA DE FIGURAS ...................................................................................................... 6
NOMENCLATURA....................................................... ¡Error! Marcador no definido.
1. INTRODUCCIÓN ........................................................................................................ 8
2. OBJETIVOS ............................................................................................................... 10
2.1. Objetivo general ................................................................................................... 10
2.2. Objetivos específicos ........................................................................................... 10
3. MARCO TEÓRICO ................................................................................................... 11
3.1. El recurso eólico .................................................................................................. 11
3.1.1. Características del viento .................................................................................. 11
3.1.2. Potencia del viento ............................................................................................ 12
3.1.3. Turbinas eólicas ................................................................................................ 12
3.1.4. Tipos de turbinas eólicas................................................................................... 14
3.2. Fundamentos aerodinámicos................................................................................ 16
3.2.1. Teoría del disco actuador .................................................................................. 16
3.2.2. Coeficientes aerodinámicos .............................................................................. 17
3.3. Perfiles aerodinámicos ......................................................................................... 18
3.3.1. Fuerzas sobre el perfil ....................................................................................... 18
3.3.2. Coeficientes de sustentación y arrastre ............................................................. 19
3.3.3. Velocidades y ángulos de estudio en un perfil ................................................. 19
3.4. Metodología de diseño de rotores eólicos............................................................ 21
3.4.1. Condición de extracción de potencia ................................................................ 21
3.4.2. Diseño simplificado de rotores ......................................................................... 21
3.4.3. Diseño con arrastre y pérdidas en punta de rotores .......................................... 22
4. DISEÑO DEL ROTOR EÓLICO............................................................................... 24
4.1 Definición de parámetros de diseño ...................................................................... 24
4.2 Selección del perfil aerodinámico ......................................................................... 25
4.3. Geometría óptima del rotor .................................................................................. 25
5. MANUFACTURA DEL ROTOR .............................................................................. 27
5.1. Modelamiento del aspa ........................................................................................ 27
5.2. Modelamiento de la nariz y cubo ......................................................................... 27
5.3. Selección de torre y generador ............................................................................. 28
4
5.4. Diseño del acople eje – cubo – nariz ................................................................... 30
5.5. Configuraciones del rotor .................................................................................... 30
5.6. Construcción del rotor ......................................................................................... 32
6. MONTAJE EXPERIMENTAL .................................................................................. 34
6.1. Medición de la densidad del aire ......................................................................... 34
6.2 Medición de la velocidad del viento ..................................................................... 35
6.3. Medición del momento de Inercia ....................................................................... 36
6.3.1. Metodología péndulo trifilar ............................................................................. 36
6.3.2. Metodología Inventor ....................................................................................... 37
6.4. Descripción del experimento ............................................................................... 38
6.5. Banco de pruebas ................................................................................................. 39
7. RESULTADOS .......................................................................................................... 40
7.1. Curvas de rendimiento ......................................................................................... 40
7.2. Análisis general .................................................................................................... 42
7.3. Limitaciones del experimento .............................................................................. 44
8. CONCLUSIONES Y RECOMENDACIONES ......................................................... 45
8.1. Conclusiones ........................................................................................................ 45
8.2. Recomendaciones ................................................................................................ 45
9. BIBLIOGRAFÍA ........................................................................................................ 46
10. ANEXOS .................................................................................................................. 47
10.1. Resultados geometría óptima de diseño............................................................. 47
10.2. Plano aspa .......................................................................................................... 48
10.3. Plano de la nariz y el cubo ................................................................................. 49
10.4. Plano del acople eje – cubo – nariz .................................................................... 51
10.5. Medición momento de Inercia – Configuración 4 Palas.................................... 52
10.6. Medición momento de Inercia – Configuración 8 Palas.................................... 53
10.7. Medición momento de Inercia – Configuración 16 Palas.................................. 54
5
LISTA DE FIGURAS
Figura 1. Patrones globales de viento (Ramírez, 2006). ................................................. 11
Figura 2. Acoplamiento entre el rotor y la góndola ([Link], 2018). .... 13
Figura 3. Elementos pertenecientes a la góndola de una turbina de aerogeneración
([Link], 2018). ...................................................................................... 14
Figura 4. Turbinas de eje (a) horizontal y (b) vertical (Battista, 2000). ......................... 15
Figura 5. Tubo de corriente y disco actuador (Burton, Sharpe, & Jenckins, 2000). ...... 16
Figura 6. Fuerza de sustentación (FL) y fuerza de arrastre (FD) sobre un perfil alar
(Battista, 2000). .............................................................................................................. 18
Figura 7. Coeficientes de sustentación y arrastre de un perfil a diferentes ángulos de
ataque (Pintea, Popescu, & Borne, 2010). ...................................................................... 19
Figura 8. Diagrama de velocidades (Pinilla, Notas de lectura - Curso electivo en energía
eólica, 2017). .................................................................................................................. 20
Figura 9. Velocidades y fuerzas sobre el elemento de aspa (Pinilla, Curso Electivo de
Aerodinámica Básica, 2018)........................................................................................... 21
Figura 10. Características perfil Placa Curvada con 5% de flecha. ................................ 25
Figura 11. Descripción de la geometría del aspa . .......................................................... 26
Figura 12. Forma del aspa. ............................................................................................. 26
Figura 13. Modelo CAD del aspa. .................................................................................. 27
Figura 14. Modelo CAD de la nariz y su tapa. ............................................................... 28
Figura 15. Modelo CAD del cubo. ................................................................................. 28
Figura 16. Comparación generadores. ............................................................................ 29
Figura 17. Caracterización motor 3 Phase Generator. .................................................... 29
Figura 18. Modelo CAD del acople................................................................................ 30
Figura 19. Rotor bajo la configuración 4 palas. ............................................................. 30
Figura 20. Rotor bajo la configuración 8 palas. ............................................................. 31
Figura 21. Rotor bajo la configuración 16 palas. ........................................................... 31
Figura 22. Proceso de manufactura de las aspas. ........................................................... 32
Figura 23. Aspa tras el proceso de manufactura. ............................................................ 32
Figura 24. Nariz prototipada. .......................................................................................... 33
Figura 25. Rotor manufacturado para configuración 8 palas y 16 paslas....................... 33
Figura 26. Montaje del tubo de Pitot para la medición de velocidad del viento ............ 35
Figura 27. Esquema de conexiones del instrumento de medición VAISALA PTB330TS.
........................................................................................................................................ 35
Figura 28. Medición del momento de inercia para el rotor de 8 aspas mediante el
péndulo trifilar. ............................................................................................................... 37
Figura 29. Opción Iproperties y medición del momento de Inercia. .............................. 37
Figura 30. Turbina eólica manufacturada y montada en el banco de pruebas. ............... 39
Figura 31. Vaisala PTB330TS (Izquierda) y Tarjeta de adquisición NI9215 (Derecha).
........................................................................................................................................ 39
Figura 32. Curva de rendimiento del rotor vs velocidad específica ............................... 40
Figura 33. Curva de momento par vs velocidad específica ........................................... 40
6
Figura 34. Comparación coeficiente de rendimiento para las 3 configuraciones de rotor
........................................................................................................................................ 41
Figura 35. Comparación coeficiente de momento-par para las 3 configuraciones de rotor
........................................................................................................................................ 41
Figura 36. Comparación coeficiente de momento-par experimental vs teórico ............. 42
Figura 37. Comparación coeficientes de rendimiento para diferentes rotores (UCLM,
2012) ............................................................................................................................... 43
7
1. INTRODUCCIÓN
La energía eólica está creciendo en todo el mundo como una fuente importante de
energía renovable (DaSolar, 2018). Para satisfacer la demanda de fuentes de energía
alternativas, un gran número de generadores eólicos o aerogeneradores se están
implementando en parques eólicos tanto terrestres como marinos. La energía eólica
doméstica también está creciendo rápidamente a medida que los generadores de energía
eólica residencial son cada vez más eficientes y prácticos. Sin embargo, existen otras
aplicaciones además de la extracción de energía eléctrica, especificamente, el
aprovechamiento de la energía que posee el viento para bombear agua o aerobombeo.
8
Tabla 1. Potencial eólico (UPME, 2015).
La ausencia de proyectos eólicos hoy en día más allá del parque eólico de Jepirachi
(19,5 MW de capacidad nominal) responde a la existencia de barreras que impiden o
complican la viabilidad de estos proyectos por factores como son, en el caso de La
Guajira, la falta de la infraestructura eléctrica necesaria para desalojar la energía
producida hacia el interior del país, la complejidad de los procesos de negociación con
las comunidades que habitan la región y, en términos generales, la ausencia de un marco
normativo y regulatorio que viabilice la participación de este tipo de energía de carácter
variable en el mercado eléctrico nacional y que valore el aporte que esta le haría al
sistema en términos de complementariedad hídrica (UPME, 2015).
Para tal fin, en este documento se desarrolla una pequeña introducción a la energía
eólica, sus diferentes aplicaciones y sus principales diferencias; posteriormente se
estudiarán las diferentes características de diseño de los rotores eólicos, terminando con
la explicación del diseño implementado; luego se describirá el proceso de manufactura
del rotor eólico mediante la utilización de un proceso de manufactura no convencional,
conocido como manufactura aditiva, especificamente, impresión 3D. Paso siguiente, se
describirá el banco de pruebas utilizado para la toma de datos, así como la metodología
para la obtención y procesamiento de los mismo; luego, se muestran los resultados
obtenidos de las pruebas con su respectivo análisis. Por último, se desarrollan las
conclusiones obtenidas a partir del desarrollo del proyecto así como posibles cambios o
modificaciones para tener en cuenta en trabajos futuros en la misma área del
conocimiento.
9
2. OBJETIVOS
2.1. Objetivo general
Diseñar y realizar un montaje experimental para el rotor con el fin de encontrar sus
curvas de rendimiento (Cp y Ct) y con ello entender su comportamiento.
10
3. MARCO TEÓRICO
Las masas de aire caliente de la zona ecuatorial ascienden generando así la Zona de
Confluencia Intertropical (ZCIT) ubicada paralela al ecuador. También, los vientos de la
superficie se devuelven hacia la zona ecuatorial, creando así los vientos alisios del
noreste y sureste. Por último, existen los vientos del oeste ubicados en zonas templadas
en latitudes entre 30 y 60 grados de ambos hemisferios.
11
3.1.2. Potencia del viento
1 3
𝑃𝑣𝑖𝑒𝑛𝑡𝑜 = 𝜌𝑈∞ 𝐴 [𝑊] (1)
2
En cuanto a la velocidad del viento, es el componente que más afecta la potencia debido
a la relación cúbica, lo que genera una gran amplificación en la potencia cuando la
velocidad del viento aumenta. En términos específicos, si la velocidad de viento se
duplica, por ejemplo de 2 a 4 m/s, la potencia disponible se incrementa en 8 veces. Es
por esto que localizaciones con una alta velocidad promedio de viento son preferidos
para la implementación de parques eólicos.
Por último, al área por el que el viento fluye es un factor también importante en la
estimación de la potencia del viento, sin embargo, esta área está definida por el
diseñador del rotor eólico al momento de determinar el diámetro del rotor.
Las turbinas eólicas son máquinas accionadas por la energía eólica, sobre las cuáles se
intercambia cantidad de movimiento con el viento mediante el accionamiento de un
rotor. La energía mecánica del eje del rotor puede ser aprovechada en diversas
aplicaciones como bombear agua (aerobombas), moler (molinos de viento) o para la
generación de energía eléctrica (aerogeneradores).
Aunque las turbinas eólicas cuentan con diferentes fines o usos, independientemente del
fin, estas están compuestas principalmente por los mismos componentes, los cuáles son
(Figueredo, 2010):
12
Torre y cimiento: Los cimientos se construyen para garantizar la estabilidad de la
turbina eólica, estos pueden ser superficiales o profundos. El tipo de cimentación
depende tanto del tipo de suelo donde se implementará la turbina eólica, así como del
peso de la turbina. La torre se construye tanto para resistir el peso de la góndola y de los
álabes o aspas del rotor, como para absorber las cargas cuasadas por la variación de la
potencia del viento.
Existen varios tipos de torres sin embargo las más utilizadas son las torres tubulares de
acero, las torres de concreto y las torres híbridas (combinación entre acero y concreto).
Los criterios de selección de torre más utilizados son la altura de la turbina, de manera
que con torres de acero es más fácil alcanzar mayores alturas que con torres de cemento.
El tipo de aplicación, de manera que las turbinas para aerobombeo tienen un diseño de
torre diferente a turbinas para aerogeneración y la facilidad de acceso al sitio de
operación de la turbina, puesto que las torres de concreto se construyen en el lugar
evitando el problema de transportación.
Rotor: El rotor es el componente en el cual las palas convierten la energía del viento en
movimiento mecánico rotacional. El rotor está compuesto por las aspas y el buje o cubo,
en el cual se unen las aspas, el cono central y el eje del rotor, este último permite una
conexión entre el rotor y la góndola (Figura 2).
La disposición del eje del rotor permite la clasificación de turbinas eólicas entre turbinas
de eje vertical y horizontal. Actualmente la mayoría de rotores empleados en la
generación de energía tienen tres palas, un eje horizontal y un diámetro entre 40 y 90
metros. Según Conrado Moreno Figueredo, doctor en ciencias técnicos y miembro de la
Academia de Ciencias de Cuba, la experiencia ha demostrado que el rotor de tres palas
es más eficiente para la generación de energía en las grandes turbinas, además los
rotores de tres palas tienen una mejor distribución de masa, lo que permite una rotación
más estable (Figueredo, 2010).
13
Las aspas del rotor son fabricados principalmente en fibra de vidrio o fibra de carbón
reforzado con plástico para aplicaciones relacionadas con la aerogeneración, mientras
que en otras aplicaciones como el aerobombeo las aspas son fabricadas principalmente
en hojas de acero o aluminio.
14
las palas están situadas a una elevada altura, con lo cual la velocidad media del viento es
mayor, y la intensidad de turbulencia es menor, que a nivel del suelo.
Por último, las turbinas verticales más utilizadas actualmente son las turbinas Darrieus y
Savonious; la primera para aerogeneración mientras que la segunda para aerobombeo.
Aunque actualmente se utilicen en mayor proporción las turbinas horizontales, las
turbinas verticales también son ampliamente estudiadas debido a su facilidad de
construcción y controlabilidad.
15
3.2. Fundamentos aerodinámicos
Esta teoría recurre a un modelo simplificado del rotor eólico para explicar cómo y por
qué se produce un intercambio de energía cuando el viento atraviesa el aerogenerador.
Se considera una turbina ideal, con un número infinito de palas y representada por un
disco de espesor nulo y área igual a la barrida por las palas al girar, conocido como
disco actuador (Figura 5)
Figura 5. Tubo de corriente y disco actuador (Burton, Sharpe, & Jenckins, 2000).
Por lo cual, con lo expresado en la figura anterior es posible definir 3 estados diferentes,
el primer estado en el que el viento tiene una velocidad no perturbada (𝑈∞ ) y se
encuentra antes del disco actuador, el segundo estado en el que el viento interactua con
el disco actuador y por último, el tercer estado en el que el viento está después del disco
actuador con una velocidad aguas abajo (𝑉2).
𝑈∞ + 𝑉2 𝑚
𝑉𝑑𝑖𝑠𝑐𝑜 = [ ] (2)
2 𝑠
Adicionalmente, la potencia extraída por el disco está relacionada por la energía
contenida en el volumen de viento que siente el cambio en la velocidad, y la cantidad de
volumen que pasa a través del rotor en un determinado tiempo:
1 2
𝑃 = 𝜌(𝑈∞ − 𝑉22 )𝑉𝑑𝑖𝑠𝑐𝑜 𝐴 [𝑊] (3)
2
La anterior relación obtenida, puede ser optimizada, con el objetivo de obtener la
máxima potencia que se puede extraer en términos de la velocidad:
16
𝑑𝑃 𝑈∞ 𝑚 2𝑈∞ 𝑚
=0 → 𝑉2 = [ ] 𝑉𝑑𝑖𝑠𝑐𝑜 = [ ] (4)
𝑑𝑉1 3 𝑠 3 𝑠
Por lo cual, la máxima potencia extraída por el disco actuador se da cuando este
experimenta una velocidad correspondiente al 66.6% de la velocidad del viento antes de
ser perturbada, y la velocidad del viento después de interactuar con el disco es 33.3% de
la velocidad del viento sin perturbar. Replanteando la ecuación de potencia en la
condición de máxima extracción se obtiene:
16 1 3
𝑃𝑚𝑎𝑥 = ( 𝜌𝑈 𝐴) [𝑊] (5)
27 2 ∞
Este factor, el cual multiplica la potencia del viento antes de ser perturbada (16/27), es
conocido como el Limite de Betz (59,26%) y representa la potencia máxima que puede
ser extraída del viento. Este límite significa entonces que el disco actuador no puede
extraer el total de la energía del viento, ya que debe dejar fluir el viento.
𝑈∞ − 𝑉𝑑𝑖𝑠𝑐𝑜
𝑎= (6)
𝑈∞
Velocidad Específica (𝝀): Este coeficiente relaciona la velocidad lineal de punta del
rotor con la velocidad no perturbada del viento.
Ω𝑅
𝜆= (7)
𝑈∞
𝑃
𝐶𝑝 = (8)
1 3
2 𝜌𝑈∞ 𝐴
17
Coeficiente de Momento-Par (𝑪𝑻 ): Este coeficiente surge de la realidad que se tiene a
la hora de analizar rotores eólicos, ya que estos no pueden tener físicamente un número
infinito de aspas ni una velocidad de rotación infinita; por lo que se implica que parte de
la generación de potencia se ve manifestada como un momento-par en el eje de
rotación. En consecuencia, la potencia extraída por un rotor eólico es expresada
por: 𝑃 = 𝜔𝑇. Aplicando la relación anterior, en conjunto con las definiciones de
coeficiente de rendimiento y velocidad específica; es posible econtrar una relación para
el coeficiente de momento par:
𝑇
𝐶𝑇 = (9)
1 2
2 𝜌𝑈∞ 𝐴𝑅
Figura 6. Fuerza de sustentación (FL) y fuerza de arrastre (FD) sobre un perfil alar (Battista, 2000).
18
representados en cambios de presión alrededor del cuerpo y estas diferencias de presión
es lo que produce la fuerza aerodinámica.
𝑙
𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑆𝑢𝑠𝑡𝑒𝑛𝑡𝑎𝑐𝑖ó𝑛 → 𝐶𝑙 = (10)
1 2
2 𝜌𝑈∞ 𝑐
𝑑
𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝐴𝑟𝑟𝑎𝑠𝑡𝑟𝑒 → 𝐶𝑑 = (11)
1 2
2 𝜌𝑈∞ 𝑐
𝛼
Figura 7. Coeficientes de sustentación y arrastre de un perfil a diferentes ángulos de ataque (Pintea, Popescu, &
Borne, 2010).
Adicional a la velocidad no perturbada del viento y la velocidad lineal del aspa, existe
una velocidad relativa que experimenta el aspa, también llamada velocidad efectiva
(𝑊). Su principal diferencia con respecto a la velocidad geométrica es que la velocidad
efectiva tiene dentro de sus componentes los factores de inducción axial y tangencial, de
manera que:
19
𝑚
𝑉𝑔𝑒𝑜 = √(Ω𝑟)2 + 𝑈∞
2 [ ] (12)
𝑠
2
𝑚
𝑊 = √(1 + 𝑎′ )2 (Ω𝑟)2 + (1 − 𝑎)2 𝑈∞ [ ] (13)
𝑠
Por lo cual, es posible crear un diagrama de velocidades en el cual se observen todas las
velocidades implicadas en el estudio de perfiles.
Figura 8. Diagrama de velocidades (Pinilla, Notas de lectura - Curso electivo en energía eólica, 2017).
Velocidades: Ángulos:
20
Figura 9. Velocidades y fuerzas sobre el elemento de aspa (Pinilla, Curso Electivo de Aerodinámica Básica, 2018).
El momento-par y empuje diferenciales pueden ser expresados como una función de las
características aerodinámicas del perfil, es específico, entre menor el arrastre del perfil,
mayor el momento-par producido en el elemento de aspa.
Por tal motivo, el perfil debe contar con una furza de sustentación alta y una fuerza de
arrastre lo más pequeña posible para obtener una máxima extracción de potencia. Por
ello, el perfil a elegir deberá tener una alta relación entre sustentación y arrastre (𝐸), lo
cual determina un ángulo de ataque óptimo (𝛼𝑜𝑝𝑡 ).
Teniendo esto en cuenta, en la condición de diseño óptimo del rotor, el ángulo de ataque
debe permanecer constante desde la raíz hasta la punta, esto genera que el ángulo de
calaje y la cuerda varíe entre la raiz y la punta. En específico, existe una relación entre
el ángulo 𝜑, el ángulo de calaje 𝛽 y el ángulo de ataque 𝛼:
𝜑 =𝛽+𝛼 (14)
Por último, para tener definido el diseño del rotor es necesario conocer en todo lugar
entre la raíz y la punta el valor de la cuerda y el ángulo de calaje que preserve el ángulo
de ataque óptimo para así garantizar la máxima extracción de potencia.
21
Por lo cual el ángulo que optimiza el momento-par a lo largo del aspa para cada
posición radial es 𝜑 = 2𝜖. Teniendo esto en cuenta, el valor de la cuerda y el ángulo de
calaje a lo largo del aspa son:
2𝜋𝑟
𝑐(𝑟) = ∗ 4(1 − 𝑐𝑜𝑠𝜑) [𝑚] (15)
𝐶𝑙−𝑜𝑝𝑡 ∗ 𝐵
2 𝑅
𝜑(𝑟) = tan−1 [ ] [°] (16)
3 𝜆𝑑𝑖𝑠𝑒ñ𝑜 𝑟
De esta manera, a partir de definir una velocidad diseño, un diámetro del rotor y un
número de palas es posible diseñar el rotor.
Gracias a esto, es necesario encontrar para cada posición radial debe cumplirse la
relación entre factores de inducción (Ecuación 19). Adicionalmente la relación entre
coeficientes de sustentación y arrastre puede re-escribirse como (Ecuación 18):
𝐶𝑙 𝜆𝑟 (𝑎′ + 𝑎) 𝜆𝑟 𝑑𝑎′ 𝑎′
= 𝑐𝑜𝑛 𝜆 𝑟 = (18) = (19)
𝐶𝑑 𝑎(1 − 𝑎) − 𝜆2𝑟 𝑎′(1 + 𝑎′ ) 𝑅 𝑑𝑎 1 − 𝑎
Por lo cual al combinar las ecuaciones 18 y 19 es posible encontrar dos ecuaciones que
parten de la relación sustentación-arrastre (𝐸) y la velocidad específica de diseño en
cada posición radial (𝜆𝑟 ) para resolver 𝑎 y 𝑎′:
22
1−𝑎 1
𝛽(𝑟) = tan−1 [ ∗ ] − 𝛼𝑜𝑝𝑡 [°] (23)
1 + 𝑎 ′ 𝜆𝑟
De esta manera, a partir de definir una velocidad diseño, un diámetro del rotor y un
número de palas es posible diseñar el rotor.
23
4. DISEÑO DEL ROTOR EÓLICO
Los parámetros iniciales de diseño estan contenidos en la Tabla 2, tal como se muestra
más adelante:
Tabla 2. Parámetros de diseño seleccionados.
Parámetro Valor
Diámetro (D) 0,3 m
Velocidad específica (λ) 0,5
Número de aspas (B) 16
Densidad (ρ) 0,9 kg/m3
Por último, para definir el número de aspas es necesario tener en cuenta que existe una
relación entre el número de aspas y el tamaño del aspa, de manera que entre mayor sea
el número de aspas, menor será el tamaño del aspa. Esto es importante, puesto que al ser
una turbina diseñada para aerobombeo, esta debe contar con un número alto de aspas
para así contar con una alta solidez. Por lo cual, se definío el número de aspas como 16
tras realizar varias iteraciones al momento del diseño, variando el número de aspas hasta
encontrar una geometría del aspa que permita conectar el número de aspas definido en
el diseño.
24
4.2 Selección del perfil aerodinámico
Para la selección del perfil aerodinámico se revisó en primer lugar literatura acerca de
perfiles aerodinámicos, enfocándose en perfiles para uso en condiciones de bajo número
de Reynolds (10000-100000), con curvas de rendimiento encontradas a partir de
experimentos en el túnel de viento y con información detallada en un rango extenso de
ángulos de ataque (0° y 90°).
Teniendo todo esto en cuenta, se escoge el perfil de placa curvada con 5% de flecha
debido no solo gracias a sus características aerodinámicas sino también por su
geometría simple, lo cual facilitará el procedimiento de manufactura del rotor. A
continuación se presenta las características aerodinámicas del perfil para el rango de
ángulos de ataque especificado anteriormente:
Para encontrar la geometría óptima del rotor es necesario primero definir el tipo de
metodología de diseño. En este caso, se selecciona la metodología de diseño incluyendo
arrastre y pérdidas en punta del aspa. Así pues, la geometría resultante es expresada a
partir del ángulo de calaje 𝜷(𝒓) y el tamaño de la cuerda 𝒄(𝒓) como función de la
posición radial, tal como se muestra en la Figura 11:
25
Figura 11. Descripción de la geometría del aspa .
A partir de la Figura 11, es posible concluir que la cuerda varía de 0,0095 m a 0,088 m
para volver a disminuir a 0 m, mientras que el ángulo de calaje disminuye de 56° en la
raíz hasta 38° en la punta. También, si se considera que el punto de cuarta cuerda debe
ser una línea recta, es posible establecer la forma del aspa, la cual se muestra en la
Figura 12 a continuación:
0,2
0,1
0
0 0,02 0,04 0,06 0,08 0,1 0,12 0,14 0,16
Cuerda/R
-0,1
-0,2
-0,3
-0,4
-0,5
Radio (m)
Por último, en el Anexo 1 es posible encontrar la hoja de cálculo del diseño final
utilizado.
26
5. MANUFACTURA DEL ROTOR
5.1. Modelamiento del aspa
Para empezar con la manufactura, fue necesario la realización de un modelo CAD del
aspa para así poder manufacturar la misma a partir de la máquina de prototipado rápido.
En la construcción del modelo CAD se respetaron tanto los ángulos de calaje como la
cuerda para cada posición radial. A continuación se muestra el modelo CAD del aspa:
En la realización del CAD para las aspas, también se realizaron 2 huecos en la base de
cada aspa para así después conectar todas las aspas y la nariz en el cubo. El material
escogido para la construcción fue PLA, después de comprobar la factibilidad de las
propiedades mecanicas para la aplicación. Por último, también se aumentó el área de
contacto entre la base y la raíz del aspa para mejorar la resistencia del aspa puesto que
esta zona es un concentrador de esfuerzos. En el Anexo 2 es posible encontrar el plano
del aspa.
27
Figura 14. Modelo CAD de la nariz y su tapa.
Por otra parte, el cubo se diseña para conectar de manera fácil mediante prisioneros M3
de 10 mm de longitud a todas las aspas. El cubo cuenta con un diámetro total de 48 mm,
sin embargo, tiene un hueco de 12.5 mm o ½ pulgada donde ira conectado el eje de
transmisión de potencia del rotor. A continuación se presenta el modelo CAD del cubo:
Por último, tanto en el Anexo 3 como en el Anexo 4 es posible encontrar tanto el plano
de la nariz y el cubo.
28
Voltaje sin resistencia vs RPM
18
16
14
Voltaje (V) 12
10
8 3 Phase Generator
6
DunkenMotoren G30.2
4
2
0
0 200 400 600 800
RPM
Corriente vs RPM
0,07
0,06
Corriente (A)
0,05
0,04 5 Ohms
0,03 10 Ohms
0,02
50 Ohms
0,01
210 Ohms
0
0 200 400 600 800 1000
RPM
Voltaje vs RPM
0,6
0,5
0,4 5 Ohms
Voltaje (V)
0,3 10 Ohms
0,2 50 Ohms
210 Ohms
0,1
Cortocircuito
0
0 200 400 600 800 1000
RPM
29
5.4. Diseño del acople eje – cubo – nariz
Para unir la nariz, el cubo, las aspas, el eje de rotación y el eje del generador es
necesario diseñar un acople, el cual pueda ensamblar el cubo, la nariz y las aspas y al
mismo tiempo permitir una conexión entre el eje de rotación y el eje del generador. A
continuación se presenta el modelo CAD del acople:
Por último, antes de empezar con la construcción del rotor, es necesario definir las tres
distintas configuraciones en las que más adelante será puesto a prueba el rotor, en este
sentido, cada configuración varia de las demás en el número de palas.
30
Configuración 8 palas: En esta configuración el rotor cuenta con una velocidad de
arranque de 3,2 m/s. A continuación se presenta el modelo CAD de la configuración:
Por último, se espera que la configuración que presente un mejor rendimiento sea la
configuración completa o 16 aspas puesto que tiene la menor velocidad de arranque,a
adicionalmente el rotor fue diseñado para contar con un número de aspas de 16, por lo
cual se espera que las condiciones de diseño se alcancen bajo dicha configuración.
31
5.6. Construcción del rotor
Para la construcción del rotor fue necesario diseñar la manera de como permitir que el
peso del aspa al momento de su construcción no afectará la base o los agujeros
diseñados para los prisioneros M3. En adición, debido a las limitaciones de espacio de
la máquina de prototipado rápido no fue posible imprimir todo el rotor, sino que este
tuvo que ser impreso por partes. Por tal motivo, se dividío el rotor en 8 pares de aspas.
A continuación se aprecia el proceso de manufactura aditiva así como la realización de
un par de aspas.
32
A continuación, se fabricó la nariz para el rotor a partir de resina en una máquina de
prototipado rápido, siguiendo así el diseño encontrado en el numeral 5.2. Es importante
recordar que esta geometría es la óptima para dirigir el viento hacia el rotor puesto que
minimiza las perturbaciones en el viento incidente.
33
6. MONTAJE EXPERIMENTAL
Para obtener un buen nivel de precisión al momento de medir la densidad del viento, se
realiza una medición indirecta de la misma teniendo en cuenta el cambio de las
condiciones experimentales manifestada en características como la temperatura, la
presión y la humedad relativa.
En primera instancia, se parte desde la ecuación de estado del aire como un gas real:
𝜌𝑍𝑅𝑇 = 𝑝𝑀 (27)
Ahora, se tiene en cuenta que el aire es aire húmedo, es decir, que posee un porcentaje
significativo de masa molar representado como vapor de agua (𝑥𝑣 ) acompañado de un
restante de aire seco. Esta relación de molaridad se representa de la siguiente manera:
𝑀𝑣 𝑘𝑔
𝑀 = 𝑀𝑎 [1 − 𝑥𝑣 (1 − )] [ ] (28)
𝑀𝑎 𝑚𝑜𝑙
Donde 𝑀𝑎 corresponde a la masa molar del aire húmedo (0.02896 𝑘𝑔/𝑚𝑜𝑙 ) y 𝑀𝑣
corresponde a la masa molar del agua (0.01812 𝑘𝑔/𝑚𝑜𝑙 ). Con la información anterior,
se hace necesario encontrar el valor correspondiente a la fracción de masa molar del
vapor de agua; este, se puede estimar a partir de la humedad relativa del aire relacionado
con otros factores característicos del vapor de agua:
𝑃𝑠𝑣(𝑇)
𝑥𝑣 = ℎ𝑓(𝑝, 𝑇) (29)
𝑝
𝑝𝑀𝑎 𝑀𝑣
𝜌 = [1 − 𝑥𝑣 (1 − )] (30)
𝑍𝑅𝑇 𝑀𝑎
34
6.2 Medición de la velocidad del viento
La velocidad de viento en la sección de pruebas del túnel de viento se mide por medio
de la implementación de un tubo de Pitot y un instrumento meteorológico llamado
Vaisala PTB330TS. A continuación, se muestra la configuración del tubo de Pitot en la
sección de pruebas:
Figura 26. Montaje del tubo de Pitot para la medición de velocidad del viento
A partir del tubo de Pitot se obtienen los datos relacionados a la presión total y la
presión estática del aire en la sección del túnel de viento gracias a la implementación del
Vaisala, el cual se conecta al tubo de Pitot por medio de mangueras elásticas para
obtener la medición correspondiente a las presiones de manera digital.
Adicionalmente, para poder calcular de manera precisa la densidad del aire en la sección
de pruebas, se utilizan un par de accesorios adicionales del instrumento Vaisala, que
permiten medir la temperatura y humedad relativa en la sección de pruebas. A
continuación, se muestran los instrumentos y sensores relacionados al Vaisala y la
interfaz gráfica de la obtención de datos:
35
Una vez obtenidas las mediciones relacionadas con la densidad y las diferentes
presiones del fluido, se encuentra la velocidad del viento con la siguiente fórmula:
2(𝑃∞ − 𝑃0 )
𝑉∞ = √ (31)
𝜌
Donde 𝑃∞ es la presión total del viento y 𝑃0 corresponde a la presión estática del viento;
la resta de estos dos valores corresponde a la presión dinámica del fluido.
Para la obtención de las curvas de rendimiento del rotor para cada configuración, es
necesario medir el momento de inercia para cada configuración. Para medir el momento
de inercia se utilizan dos metodologías diferentes, la primera el péndulo trifilar y la
segunda mediante Autopdesk Inventor.
1. Colocar el rotor en una posición en la que la longitud de los 3 hilos sea L. Los
hilos han de estar paralelos (¡Error! No se encuentra el origen de la
referencia.).
4. Repetir el paso anterior 5 veces, para así poder realizar un promedio de los
periodos obtenidos.
36
16 Palas 0,003627
Figura 28. Medición del momento de inercia para el rotor de 8 aspas mediante el péndulo trifilar.
37
Tabla 4. Resultados metodología inventor
Metodología Inventor
Configuración I [kg m2]
4 Palas 0,000921
8 Palas 0,001821
16 Palas 0,003620
Por último, del Anexo 6 al Anexo 8 se encuentran todas las mediciones para las
distintas configuraciones del rotor mediante el Software Inventor.
1. Medición del voltaje del generador acoplado al motor como función de tiempo,
cuando esta sujeto a una velocidad de viento constante, desde su condición de
parada hasta que el rotor alcance su máxima velocidad angular.
Para garantizar que la toma de datos sea ordenada y precisa, para cada configuración se
realizará el experimento para 3 diferentes velocidades: la velocidad de arranque, la cual
varia entre las diferentes configuraciones, 4.3 m/s y 6.4 m/s. Las otras dos velocidades
fueron seleccionadas puesto que a velocidades superiores de 7 m/s la torre empieza a
vibrar, dañando la toma de datos.
38
6.5. Banco de pruebas
El banco de pruebas para este proyecto será la sección de pruebas del túnel de viento
TVIM-49-60-1X1 del Departamento de Ingeniería Mecánica de la Universidad de los
Andes. El túnel de viento permitirá generar una corriente de viento con las velocidades
establecidas en la descripción del experimento. Adicional al túnel de viento, es
necesario los siguientes instrumentos:
Turbina eólica: Tal como lo muestra la Figura 30, la turbina eólica esta compuesta por
el rotor, el generador y la torre.
39
7. RESULTADOS
Cp
0,18
0,16
0,14
0,12
0,1
Cp
0,08
0,06
0,04
0,02
0
-0,2 0 0,2 0,4 0,6 0,8 1 1,2 1,4
-0,02
Vel esp
Ct
0,45
0,4
0,35
0,3
0,25
Cp
0,2
0,15
0,1
0,05
0
-0,2 0 0,2 0,4 0,6 0,8 1 1,2 1,4
-0,05
Vel esp
40
Adicionalmente se presenta comparaciones para el coeficiente de rendimiento y
coeficiente de momento par para las tres configuraciones del rotor:
0,4
0,35
0,3
0,25
Ct
0,2
0,15
0,1
0,05
0
0 0,2 0,4 0,6 0,8 1 1,2
Velocidad específica
41
Por último, se realiza una comprobación teórica de la curva de momento-par,
especificamente para el rotor completo (16 aspas) y a 2.3 m/s:
0,4
0,35
0,3
0,25
Ct
0,2
0,15
0,1
0,05
0
-0,2 0 0,2 0,4 0,6 0,8 1 1,2
Vel esp
En primer lugar, las curvas obtenidas se comportan de manera similar a las gráficas
expuestas por (Pinilla, 2017). En este sentido, el coeficiente de rendimiento relaciona la
potencia extraída por el rotor a tráves del momento-par aerodinámico 𝑻, transmitido en
el eje a una velocidad angular 𝛀, y la potencia que contiene el viento a una velocidad
no-perturbada 𝑼∞ enfrente del rotor, con un área de referencia 𝑨 correspondiente al
área barrida por las aspas en un giro.
Tal como se pueden observar en la Figura 32 y Figura 33, existe una conexión entre los
coeficientes de rendimiento, momento par y velocidad específica. En específico, para
velocidades específicas altas es previsible valores bajos de momento-par,
correspondiente a baja rotación de la estela y finalmente a un coeficiente de rendimiento
alto. Este tipo de comportamiento es el de las turbinas de aerogeneración donde se
requieren altas velocidades de rotación en el eje y un bajo momento-par de alimentación
a los generadores. Por el contrario, para velocidades específicas bajas es previsible
valores altos de momento-par, correspondiente a alta rotación de la estela y en
consecuencia un coeficiente de rendimiento bajo. Este tipo de comportamiento es el de
las turbinas de aerobombeo, en los cuáles se busca un alto momento-par capaz de
accionar las bombas.
42
específica de diseño. Otra condición que se puede determinar es el valor del coeficiente
de momento-par de arranque, cuando 𝜆 = 0. En turbinas de aerogeneración es de
esperar que este valor sea menor que el coeficiento máximo de momento-par, sin
embargo como lo muestra la Figura 33, en turbinas de aerobombeo es posible que el
momento-par máximo se encuentre en la condición de arranque, lo cual beneficia a la
bomba conectada a la turbina. Además, es posible observar que cuando 𝜆 = 1, el valor
de los coeficientes son iguales 𝑪𝑷 = 𝑪𝑻 .
Figura 37. Comparación coeficientes de rendimiento para diferentes rotores (UCLM, 2012)
43
7.3. Limitaciones del experimento
Durante el desarrollo de los experimentos fue posible observar ciertas limitaciones que
no permiten obtener con claridad resultados en algunas ocasiones. En específico, el
túnel de viento se encuentra conectado al sistema de electricidad de la Universidad, en
este sentido, mantener una velocidad constante en el túnel es muy díficil puesto que si
en algún otro laboratorio prenden alguna máquina como por ejemplo una bomba, la
energía eléctrica que recibe el túnel es menor, generando una pérdida en la velocidad
del viento adentro del mismo.
Del mismo modo, por la naturaleza del experimento de partir con el rotor sin
movimiento, no fue posible realizar pruebas con la sección cerrada del túnel de viento,
esto genera pérdidas adicionales en el rendimiento del rotor.
Por último, el método de medición no permite conocer con claridad ambas curvas a
velocidades altas, puesto que entre mayor sea la velocidad, menor tiempo se demorará
el rotor de estar quieto a alcanzar su velocidad ángular nominal, por tal motivo,
experimentos con velocidades mayores de 7 m/s no permiten un adecuado tratamiento
de los datos. Adicionalmente, desde esta misma velocidad con cualquiera de las
configuraciones del rotor es imposible mantener la velocidad angular constante puesto
que la torre comienza a vibrar, dañando así la toma de datos.
44
8. CONCLUSIONES Y RECOMENDACIONES
8.1. Conclusiones
8.2. Recomendaciones
45
9. BIBLIOGRAFÍA
2. Burton, T., Sharpe, D., & Jenckins, N. (2000). Wind Energy Handbook. Wiley.
10. Pintea, A., Popescu, D., & Borne, P. (2010). Modelling and control of wind
turbines. IFAC Proceedings Volumes, 251-256.
46
10. ANEXOS
r/R r(m) 𝜶𝒐𝒑 (°) a a' Eq1 Eq2 Error Lambda r 𝝋 (°) 𝜷 (°) 𝝐(°) f F cuerda (m) Cuerda (cm)
0,05 0,0075 0,2639 16,4805 -0,594 -0,594 7,46E-09 0,025 59,30 56,91 29,65 176,77 1 0,01 0,91
0,1 0,015 0,2671 8,0413 -0,581 -0,581 6,73E-09 0,05 58,33 55,94 29,17 84,59 1 0,02 1,76
0,15 0,0225 0,2702 5,2305 -0,577 -0,577 6,00E-09 0,075 57,37 54,98 28,68 53,83 1 0,03 2,55
0,2 0,03 0,2732 3,8278 -0,569 -0,569 5,30E-09 0,1 56,41 54,01 28,20 38,42 1 0,03 3,30
0,25 0,0375 0,2760 2,9880 -0,565 -0,565 4,63E-09 0,125 55,45 53,06 27,72 29,14 1 0,04 3,99
0,3 0,045 0,2788 2,4296 -0,562 -0,562 4,00E-09 0,15 54,50 52,11 27,25 22,93 1 0,05 4,63
0,35 0,0525 0,2815 2,0320 -0,562 -0,562 3,41E-09 0,175 53,56 51,16 26,78 18,47 1 0,05 5,23
0,4 0,06 0,2840 1,7351 -0,558 -0,558 2,88E-09 0,2 52,62 50,23 26,31 15,10 1 0,06 5,78
0,45 0,0675 0,2865 1,5051 -0,561 -0,561 2,40E-09 0,225 51,69 49,30 25,85 12,46 1 0,06 6,30
0,5 0,075 0,2888 1,3221 -0,560 -0,560 1,97E-09 0,25 50,78 48,38 25,39 10,33 1 0,07 6,76
0,55 0,0825 0,2910 1,1732 -0,561 -0,561 1,60E-09 0,275 49,87 47,48 24,94 8,56 1 0,07 7,19
2,39
0,6 0,09 0,2932 1,0498 -0,562 -0,562 1,28E-09 0,3 48,98 46,58 24,49 7,07 1 0,08 7,58
0,65 0,0975 0,2952 0,9461 -0,564 -0,564 1,01E-09 0,325 48,09 45,70 24,05 5,79 1 0,08 7,93
0,7 0,105 0,2973 0,8577 -0,572 -0,572 7,77E-10 0,35 47,22 44,83 23,61 4,67 0,99 0,08 8,23
0,75 0,1125 0,2991 0,7818 -0,574 -0,574 5,90E-10 0,375 46,37 43,98 23,18 3,68 0,98 0,08 8,47
0,8 0,12 0,3010 0,7158 -0,582 -0,582 -2,48E-08 0,4 45,53 43,13 22,76 2,80 0,96 0,09 8,61
0,85 0,1275 0,3027 0,6581 -0,585 -0,585 -1,87E-08 0,425 44,70 42,31 22,35 2,01 0,91 0,09 8,57
0,9 0,135 0,3043 0,6073 -0,591 -0,591 -1,37E-08 0,45 43,89 41,49 21,94 1,28 0,82 0,08 8,19
0,95 0,1425 0,3059 0,5622 -0,599 -0,599 -9,72E-09 0,475 43,09 40,70 21,54 0,62 0,64 0,07 6,97
0,97 0,1455 0,3021 0,5504 -0,308 -0,308 4,59E-09 0,485 42,86 40,47 21,43 0,36 0,51 0,06 5,78
0,99 0,1485 0,2987 0,5388 -0,037 -0,037 2,46E-08 0,495 42,64 40,24 21,32 0,12 0,30 0,04 3,66
1 0,15 0,3075 0,5219 -0,609 -0,609 -6,61E-09 0,5 42,30 39,91 21,15 0,00 0 0,00 0,00
Anexo 1. Resultados metodología de diseño perfil aspa.
47
10.2. Plano aspa
48
10.3. Plano de la nariz y el cubo
49
Anexo 4. Plano Cubo.
50
10.4. Plano del acople eje – cubo – nariz
51
10.5. Medición momento de Inercia – Configuración 4 Palas
52
10.6. Medición momento de Inercia – Configuración 8 Palas
53
10.7. Medición momento de Inercia – Configuración 16 Palas
54