1 UNIVERSIDAD DE CÓRDOBA INGENIERÍA Y ARQUITECTURA Escuela Politécnica Superior de Córdoba MODELADO Y SIMULACIÓN DE GENERADORES EÓLICOS Milton Fabricio Pérez Gutiérrez Trabajo fin de Máster en Control de Procesos Industriales. Córdoba, Septiembre 2011
53
Embed
UNIVERSIDAD DE CÓRDOBA - …repositorio.educacionsuperior.gob.ec/bitstream/28000/274/1/T... · los días despejados ... generan dióxido de carbono y otros contaminantes atmosféricos,
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
1
UNIVERSIDAD DE CÓRDOBA
INGENIERÍA Y ARQUITECTURA
Escuela Politécnica Superior de Córdoba
MODELADO Y SIMULACIÓN DE
GENERADORES EÓLICOS
Milton Fabricio Pérez Gutiérrez
Trabajo fin de Máster en Control de Procesos Industriales.
Córdoba, Septiembre 2011
2
D. Jorge Eugenio Jiménez Hornero, profesor de la Escuela Politécnica Superior de la Universidad de Córdoba, perteneciente al Área de Ingeniería de Sistemas y Automática del departamento de Informática y Análisis Numérico. INFORMA Que el alumno Milton Fabricio Pérez Gutiérrez ha desarrollado, bajo su dirección, el Trabajo Fin de Máster titulado “Modelado y simulación de generadores eólicos”, habiendo superado, a su juicio, todos los objetivos inicialmente propuestos. Y para que así conste
Córdoba, septiembre de 2011.
Fdo: Prof. Dr. Jorge Eugenio Jiménez Hornero
3
AGRADECIMIENTOS
Expreso mi agradecimiento a todos los profesores, compañeros y amigos, del
Máster de Control de Procesos Industriales y de la Universidad de Córdoba, que
con su ayuda continua tanto a nivel científico como personal han colaborado en
la realización de este trabajo, especialmente a Antonio que me brindo su amistad
al estar tan lejos de mi familia.
Agradezco a mi familia, fundamentalmente a mi padre, a mi madre, hermano,
hermana, mi esposa y mis hijos que con su ánimo desde el otro lado del mundo
me apoyaron en los buenos y en los malos momentos.
Finalmente, mi más especial agradecimiento al director de este trabajo, Dr. D.
Jorge Eugenio Jiménez Hornero, por su rigurosidad científica, experiencia
investigadora y, sobre todo, por su paciencia, dedicación y apoyo diario. Sin
todo ello, habría sido imposible la realización de este trabajo.
3. REVISIÓN DEL ESTADO DEL ARTE ..............................................................................8
3.1 Los aerogeneradores y su impacto en el medio ambiente. ................................................8
3.2 Potencialidad del programa FAST utilizado en el cálculo y simulación del funcionamiento del aerogenerador. ..............................................................................................9
3.3 Revisión bibliográfica con relación al recurso eólico. .................................................... 12
3.4 Modelización Matemática de la Eficiencia Aerodinámica de Aerogeneradores de Eje Horizontal. ................................................................................................................................. 14
Figura 1.Sustentación del ala. ..................................................................................................... 16Figura 2. Efecto capa límite ....................................................................................................... 16Figura 3. Efecto Coanda ............................................................................................................. 17Figura 4. Resistencia Aerodinámica según la forma ................................................................... 19Figura 5: Esquema de programas usados en FAST .................................................................... 21Figura 6. Balanceo del Rotor con viento uniforme. .................................................................... 25Figura 7. Balanceo del rotor con turbulencia .............................................................................. 25Figura 8. Momento por el movimiento de la torre viento constante ........................................... 27Figura 9. Momento por el movimiento de la torre viento turbulento .......................................... 27Figura 10.Par en LSS con viento uniforme ................................................................................. 29Figura 11. Par en LSS con turbulencia ....................................................................................... 29Figura 12. Dirección de la Góndola con viento uniforme ........................................................... 31Figura 13. Dirección de la Góndola con turbulencias. ................................................................ 31Figura 14. Aceleración angular en LSS viento uniforme ............................................................ 35Figura 15. Aceleración angular en LSS viento turbulento. ......................................................... 35Figura 16. Error estimado de la orientación de la Góndola con viento uniforme ........................ 37Figura 17. Error estimado de la orientación de la Góndola con viento turbulento. ..................... 37Figura 18. Momento de giro de la corona de orientación viento uniforme ................................. 39Figura 19. Momento de giro de la corona de orientación viento turbulento ............................... 39
6
1. INTRODUCCIÓN
El aprovechamiento por el hombre de las fuentes de energía renovables, entre ellas
la energía solar, eólica e hidráulica, es muy antiguo; desde muchos siglos antes de
nuestra era ya se utilizaban y su empleo continuó durante toda la historia hasta la
llegada de la "Revolución Industrial", en la que, debido al bajo precio del petróleo,
fueron abandonadas.
Durante los últimos años, debido al incremento del coste de los combustibles fósiles
y los problemas medioambientales derivados de su explotación, estamos asistiendo a un
renacer de las energías renovables.
Es una de las fuentes más baratas, puede competir en rentabilidad con otras fuentes
energéticas tradicionales como las centrales térmicas de carbón (considerado
tradicionalmente como el combustible más barato), las centrales de combustible e
incluso con la energía nuclear, si se consideran los costes de reparar los daños
medioambientales.
El generar energía eléctrica sin que exista un proceso de combustión o una etapa de
transformación térmica supone, desde el punto de vista medioambiental, un
procedimiento muy favorable por ser limpio, exento de problemas de contaminación,
etc. Se suprimen radicalmente los impactos originados por los combustibles durante su
extracción, transformación, transporte y combustión, lo que beneficia la atmósfera, el
suelo, el agua, la fauna, la vegetación, etc.
Las energías renovables son inagotables, limpias y se pueden utilizar de forma
autogestionada (ya que se pueden aprovechar en el mismo lugar en que se producen).
Además tienen la ventaja adicional de complementarse entre sí, favoreciendo la
integración entre ellas. Por ejemplo, la energía solar fotovoltaica suministra electricidad
los días despejados (por lo general con poco viento), mientras que en los días fríos y
ventosos, frecuentemente nublados, son los aerogeneradores los que pueden producir
mayor energía eléctrica.
Al contrario de lo que puede ocurrir con las energías convencionales, la energía
eólica no produce ningún tipo de alteración sobre los acuíferos ni por consumo, ni por
7
contaminación por residuos o vertidos. La generación de electricidad a partir del viento
no produce gases tóxicos, ni contribuye al efecto invernadero, ni destruye la capa de
ozono, tampoco crea lluvia ácida. No origina productos secundarios peligrosos ni
residuos contaminantes.
La electricidad producida por un aerogenerador evita que se quemen diariamente
miles de litros de petróleo y miles de kilogramos de lignito negro en las centrales
térmicas. Ese mismo generador produce idéntica cantidad de energía que la obtenida
por quemar diariamente 1.000 Kg. de petróleo. Al no quemarse esos Kg. de carbón, se
evita la emisión de 4.109 Kg. de CO2, lográndose un efecto similar al producido por
200 árboles. Se impide la emisión de 66 Kg. de dióxido de azufre -SO2- y de 10 Kg. de
óxido de nitrógeno -NOx- principales causantes de la lluvia ácida.
La energía eólica es independiente de cualquier política o relación comercial, se
obtiene en forma mecánica y por tanto es directamente utilizable.
Al finalizar la vida útil de la instalación, el desmantelamiento no deja huellas.
Un Parque de 10 MW:
• Evita: 28.480 Tn. Al año de CO2.
• Sustituye: 2.447 Tep. (Toneladas equivalentes de petróleo).
• Aporta: Trabajo a 130 personas al año durante el diseño y la construcción.
• Proporciona: Industria y desarrollo de tecnología.
• Genera: Energía eléctrica para 11.000 familias. (Energía Renovables:
ventajas y desentajas de la energía eólica, 2006)
8
2. OBJETIVOS • Obtener el modelo matemático de un aerogenerador.
• Simular el modelo matemático de un aerogenerador con el Programa FAST.
3. REVISIÓN DEL ESTADO DEL ARTE
3.1 Los aerogeneradores y su impacto en el medio ambiente.
La generación eléctrica por medio de aerogeneradores aumenta cada año en
el mundo tomando en consideración, que es una energía renovable que no
contamina el medio ambiente con emisiones de efecto invernadero, como lo que
ocurre con las centrales de energía termoeléctricas del tipo convencional, que
generan dióxido de carbono y otros contaminantes atmosféricos, debido al
proceso de combustión del aire con petróleo, carbón o gas. La Comisión
Nacional Europea estima que con la producción de 1 Kw de energía eléctrica
generada por el viento se dejaría de producir la emisión de 2 toneladas de
dióxido de carbono al año, además, se generarían ahorros considerables por
concepto de combustibles fósiles no quemados en un proceso de combustión.
(Energía Renovables: ventajas y desentajas de la energía eólica, 2006)
Los impactos negativos de utilizar la generación eólica están representados
en la fase de construcción, lo que repercute en impactos directos e indirectos en
la vegetación. Sin embargo, los aerogeneradores de 40 a 50 metros de diámetro
(SÁNCHEZ & VIEDMA, 2003), solo afectan alrededor del 5% del terreno
donde se instalan y al resto del mismo, es factible darle su utilización inicial.
Otro impacto negativo que se le atribuye a la energía eólica para la producción
de electricidad, es la contaminación acústica. En la actualidad este impacto
negativo se encuentra prácticamente superado, debido a que los aerogeneradores
se instalan a una distancia tal que los habitantes no están expuestos a los ruidos
que superen los 45 db, por lo general la distancia que cumple con esta condición
es inferior a 1000 metros.
El impacto visual, es un factor que a veces se considera negativo, sin
embargo, actualmente los aerogeneradores tienen una aceptación que está por
sobre el 80% en la mayoría de los países donde se encuentran instaladas granjas
de aerogeneradores. Esta estadística abarca en forma especial a los habitantes
9
que tienen su residencia cercana a la central, quienes la visualizan como una
tecnología atractiva y generadora de empleo.
Otro efecto atribuible a las Granjas eólicas está relacionado con la
mortalidad de las aves. En la práctica se ha demostrado que aún cuando las
palas aerodinámicas tienen una gran envergadura dependiendo de su potencia
nominal, estudios realizados en países como Alemania, Dinamarca, Reino
Unido y otros, han concluido que los aerogeneradores no representan una
amenaza para las aves; el mayor problema que tienen éstas se relaciona con las
líneas de transmisión de energía eléctrica.
La Asociación Europea de Energía Eólica plantea que los aerogeneradores a
pesar de su tamaño y sus palas en movimiento no representan un problema
particular para la vida de las aves (SÁNCHEZ & VIEDMA, 2003). En general,
las granjas de turbinas de viento tienen beneficios que las hacen atractivas en la
actualidad por la contribución a la disminución del efecto invernadero y a la
contaminación atmosférica, ellas generan un ambiente positivo frente al cambio
climático global. Por otra parte, el aumento de estas centrales en el mundo
provoca la menor dependencia del consumo de combustibles fósiles.
3.2 Potencialidad del programa FAST utilizado en el cálculo y
simulación del funcionamiento del aerogenerador.
Teniendo en consideración el avance tecnológico de los sistemas eólicos, ha
sido necesario desarrollar herramientas computacionales para tomar datos en
campo de los vientos existentes en un lugar determinado donde se considere
necesario llevar adelante un emplazamiento eólico, o híbridos relacionados con
esta vital energía del viento. Para el análisis de aerogeneradores y sistemas
eólicos, existen varios programa que permiten analizar el comportamiento y la
producción de potencia y energía de dichos sistemas, en este caso particular se
estudiará el programa FAST (Fatigue, Aerodynamics, Structures, and
Turbulence) por sus siglas en inglés.
“El código de FAST es un simulador integral, capaz de predecir tanto la
carga extrema y carga de fatiga de dos y tres palas de aerogeneradores de eje
horizontal (HAWTs). El código FAST es el resultado de la unión de tres
códigos distintos, el Código FAST2 de dos palas HAWTs, el Código FAST3 de
tres palas HAWTs, y el Aerodyn subrutinas aerodinámica para HAWTs
10
(LAINO & HANSEN, 2001). Una versión intermedia de FAST, llamada
FAST_AD, utiliza diferentes archivos ejecutables de turbina de dos y tres palas.
La versión de FAST que fue desarrollado en 2002, utiliza un único archivo
ejecutable para ambos tipos de turbinas. Estos cambios dieron lugar a que el
código se ejecutara muy rápidamente, por lo que el código es, en efecto, rápido.
En 2003, las características adicionales fueron incluidos en el Código de
FAST, incluyendo la capacidad de desarrollar matrices de estado periódicas
linealizadas para el diseño de los controles y la capacidad de utilizar FAST
como un preprocesador para la generación de bases de datos ADAMS ®
(LAINO & CRAIG HANSEN, 2001) de modelos de aerogeneradores.
Algoritmos de predicción de ruidos aeroacústicos también se han introducido.
Las características adicionales se agregaron de nuevo al Código FAST en
2004. El nuevo modelo añade nuevas características incluido el ángulo de
desplazamiento lateral y la inclinación del eje del rotor, balanceo del rotor,
balanceo de la cola, la inercia y la aerodinámica de cola, control de orientación,
y control de freno del eje de alta velocidad (HSS). Una interfaz se ha
desarrollado entre FAST y un controlador maestro implementado como un
proceso dynamic-link-library (librerías de unión dinámicas) (DLL) en el estilo
de paquete de programa Bladed Garrad Hassan de aerogeneradores (Bossanyi E.
A., 2003). Una interfaz también ha sido desarrollado entre FAST y Simulink ®
con MATLAB ®, lo que permite a los usuarios implementar controles
avanzados de turbinas en Simulink con la conveniente forma de diagrama de
bloques. En 2005, FAST y Adams con Aerodyn fueron evaluados por el
Germanischer Lloyd Windenergie y considerado adecuado para "el cálculo,
diseño y certificación de las cargas de turbinas eólicas en tierra" (MANJOCK,
2005). Estas características adicionales se sumaron también a los Códigos
anteriormente mencionados. Estos incluyen nuevas unidades de medición
inercial en la góndola y galgas extensiométricas en las salidas de la torre;
mejoras en el modelo de control simple de velocidad variable, y el movimiento
de la nueva plataforma de apoyo y la funcionalidad de carga. A pesar de la
incorporación en la plataforma de seis nuevos grados de libertad (DOF), el
Código fue también optimizado para que se ejecute un 15% más rápido que las
versiones anteriores (dependiendo de las opciones que está siendo modelado).
Este manual es un subconjunto actualizada de uno escrito originalmente de
Oregon State University (OSU) (WILSON, WALKER, & HEH, 1999). El
11
manual original incluyó una discusión detallada de la teoría detrás de
FAST_AD (una encarnación anterior de FAST) y la validación del código. Para
estos dos temas, por favor refiérase a la original (WILSON, WALKER, & HEH,
1999). También está disponible Buhl y otros (BUHL, WRIGHT, & PIERCE,
2001), que es una verificación estructural de FAST_AD vs Adams. FAST y
ADAMS utilizan el conjunto de subrutinas Aerodyn, por lo que el estudio de
verificación estructural, no aportó ninguna verificación de la aerodinámica del
FAST_AD. Una verificación más reciente de FAST a través de ADAMS se
presenta en Jonkman y Buhl (Jonkman & Buhl, 2004).
Los modos de operación del capítulo describe los diferentes tipos de análisis
disponibles en el programa FAST y una breve descripción sobre la forma de
ejecutar el código se presenta en el capítulo (RUNNING FAST) Ejecución
FAST. Si desea volver a compilar FAST, usted puede encontrar la información
que necesita en el capítulo (COMPILING FAST) Compilación FAST. El
capítulo (MODEL DESCRIPTION) Descripción del Modelo se analizan los
grados de libertad para aerogeneradores de dos y tres palas HAWTs. El capítulo
(CONTROLS) Controladores describe los métodos para controlar activamente
en muchos aspectos la operación de la turbina durante la simulación.
El Control activo también puede ser implementado en Simulink como se
describe en el capítulo (SIMULINK INTERFACE) Interfaz de Simulink. El
La fuerza constante ejercida por el viento uniforme hace que la aceleración vaya
incrementándose a medida que las palas van aumentando su velocidad y cuando estas
han llegado a su valor máximo la aceleración prácticamente es cero ya que su velocidad
final es constante (figura 13) existen dos picos porque pasan por dos cambios de
velocidad pero las dos son uniformes por tanto la aceleración final es cero, en el
segundo caso (figura 14) como es viento turbulento las fuerzas son cambiantes la
aceleración producida para compensar ese nuevo cambio hace que varíe
constantemente.
37
Error estimado de la orientación de la Góndola (NacYawErr) (grados).
Test12
Figura 16. Error estimado de la orientación de la Góndola con viento uniforme
Test 13
Figura 17. Error estimado de la orientación de la Góndola con viento turbulento.
0 1000 2000 3000 4000 5000 6000 7000 8000 9000-70
-60
-50
-40
-30
-20
-10
0
0 1000 2000 3000 4000 5000 6000 7000 8000 9000-20
-15
-10
-5
0
5
10
15
20
38
El error estimado para la orientación de la góndola es la diferencia entre la
referencia (ubicación actual) y el valor de la dirección del viento, por tanto el error
estimado (figura 17) es de valor negativo porque la dirección ejercida por el viento
es contrario a la referencia, pero cuando es turbulento la dirección cambia
constantemente en valores positivos y negativos ya que la góndola trata de estar
siempre en forma perpendicular a la dirección del viento (figura 18).
39
Rotación de la parte superior de la torre (con góndola) / Momento de giro de la
corona de orientación (YawBrMxn) (kN·m)
Test12
Figura 18. Momento de giro de la corona de orientación viento uniforme
Test13
Figura 19. Momento de giro de la corona de orientación viento turbulento
0 1000 2000 3000 4000 5000 6000 7000 8000 90000
100
200
300
400
500
600
700
800
900
0 1000 2000 3000 4000 5000 6000 7000 8000 90000
200
400
600
800
1000
1200
40
La orientación de la corona depende de la dirección del viento, pero el giro del rotor
depende de la fuerza que tenga en ese momento el viento, por tanto esta fuerza produce
un Momento que hace que la góndola completa se mueva horizontalmente de un lado a
otro (figura 21) este movimiento aumenta hasta su valor característico y en un momento
determinado se vuelve oscilatorio y constante, pero cuando cambiamos a un viento
turbulento las fuerzas ejercidas así como la dirección varían en el tiempo y su Momento
variar de acuerdo a estas fuerzas producidas por la variación del viento (figura 22).
41
5. CONCLUSIONES
A lo largo del presente trabajo se analizó los aspectos físicos que rodean al
aerogenerador y su movimiento y se encontró distintas ecuaciones matemáticas que
representaban a cada uno de estos principios físicos para generar un modelo
matemático, llegando a la conclusión después de un exhaustiva investigación
bibliográfica que las ecuaciones a resolver son las Ecuaciones de Navier-Stokes para
obtener resultados de las distintas partes del aerogenerador en cualquier
condición, pero estas ecuaciones (modelo matemático) se deben resolver en
forma numérica, es decir con la ayuda de un ordenador y técnicas CFD (Cálculo
de Dinámica de Fluidos), como es el caso del programa FAST y Simulink de
Matlab este último tiene los algoritmos necesarios para realizar los cálculos
numéricos en la resolución del modelo matemático. El programa FAST es un programa libre del gobierno de EEUU,
Departamento de Energía de la Oficina de Eficiencia Energética y Energía
Renovable (NREL National Renewable Energy Laboratory) que permite obtener
resultados a través del programa Simulink de Matlab, estos resultados son
valores numéricos que se pueden observar en cualquier programa que pueda
dibujar valores en tablas, se puede cambiar varios parámetros entre los que se
encuentran las características físicas del aerogenerador como peso, diámetro de
las palas, etc.
Existe un problema con este programa que el momento que no se
configura los parámetros de salida correctamente genera un error y el programa
deja de ejecutarse, por tanto debe conocerse más a profundidad el programa
como para poder realizar cambios en sus parámetros sin producir un error.
Otro problema es la cantidad de archivos adjuntos que tiene FAST
necesarios para ejecutar la simulación, se debe tener mucha precaución el
momento de cambiar los parámetros en alguno de estos archivos incluidos los
más básicos indicados en los resultados, ya que un simple cambio producirá un
error en la simulación y no se obtendrá ningún resultado.
Para demostrar el funcionamiento de FAST se realizo dos ejemplos con
dos aerogeneradores diferentes y cada uno de ellos variando la forma del viento,
se obtuvo las siguientes conclusiones:
42
Para el ejemplo 1 el viento uniforme provoca fuerzas constantes sobre el
aerogenerador, estas producirán movimiento en el rotor por medio de las palas y
estas generaran movimiento oscilatorios ya que no es un sistema balanceado,
tomando referencia esto vemos que en el Balance, el Momento en la torre y el
Par en LSS llegan a un valor determinado constate pero con oscilaciones propias
de un sistema desbalanceado, incluso al observar la orientación de la góndola
tenemos estas oscilaciones presentes, pero el comportamiento varía cuando
cambiamos a un viento turbulento donde la velocidad y la dirección son siempre
cambiantes y esto hace que el comportamiento del aerogenerador sea variable,
en muchos casos es tan errático que se debería detener el aerogenerador para que
no sufra desperfectos.
Para el ejemplo 2, como en el caso anterior cuando existe viento que
ejerce una fuerza las oscilaciones constantes son producidas por el movimiento
del rotor por el viento uniforme, en el caso de la aceleración, irá en aumento
hasta que llegue a su valor máximo de rotación y en ese momento esta
aceleración será cero porque la velocidad es uniforme, en el caso del error de la
orientación de la góndola se encontrará por la diferencia entre su referencia y su
valor que en ese momento tenga, como el viento es uniforme llega a su dirección
máxima, y el momento de giro de la corona de orientación se producirá por las
fuerzas ejercidas por las palas que se explico anteriormente y por la fuerza
ejercida del viento al cambiar su posición para encontrar su máximo valor que es
en sentido perpendicular, todas estas tienen una oscilación constante. El mismo
aerogenerador pero con viento turbulento se observa que existe variaciones
aleatorias, cambios constantes de velocidad y dirección, y estos deben adaptarse
a estos nuevos cambios y producen movimientos erráticos en la orientación y la
aceleración porque a cada momento tiene que cambiar a una nueva velocidad y
esto hace que haya aceleraciones muy pequeñas hasta que exista un nuevo
cambio, el error producido es más grande ya que sus cambios son constantes
para llegar a una nueva posición, en todas ellas los cambios pasan de un valor a
otro muy rápido, como en el caso anterior hay que analizar la situación y detener
cuando se necesite hacerlo.
43
6. PLANTEAMIENTO DE TRABAJO FUTURO
Las investigaciones que podrían desarrollarse a partir de este trabajo pueden ser:
• Estudio de un cambio de perfiles aerodinámicos en AeroDyn y
geometrías de pala que permitan conocer mejor aún las influencias de
determinados parámetros en la estructura del aerogenerador.
• Comparativa de resultados mediante estudios numéricos de gran
precisión. En este proyecto se ha utilizado los modelos de FAST para
generar resultados con ayuda del programa de desarrollo matemático
SIMULINK de MATLAB. Sin embargo un estudio 3D de la estructura
sometiéndola también a una distribución de cargas generadas por el
viento permitiría observar la respuesta numérica dinámica del sistema,
obteniendo resultados mucho más precisos.
• Muchas de las mejoras pertenecen al campo de la mecánica.
Considerando este tema complejo y relevante, lo más oportuno es que la
optimización en esta área la lleve a cabo un estudiante de ingeniería
mecánica.
• El sistema de control propuesto es básico. No obstante, se pueden
realizar estudios de algoritmos de control como control experto, óptimo,
robusto, adaptativo, neuronal, o la variación de ellos
44
7. REFERENCIAS BIBLIOGRÁFICAS
Adsuar, J. C. (2006). Principio de Vuelo. En J. C. Adsuar, Principio de Vuelo (págs. 43 - 44). Madrid: Thomson Paraninfo.
Anaya-Lara, O., Jenkins, N., Ekanayake, J., Cartwright, P., & Hughes, M. (2009). WIND ENERGY GENERATION Modelling an Control. West Sussex: John Wiley & Sons. Ltd.
BECHRAKIS, D., & DEANE, J. (2004). Wind Resourse Assessment of and Area Using Short Term Data Correlated to a Long Term Data Set. Solar Energy, Volumen 76, Elsevier Ltd. , 725-732.
BERRIOS, R., & CERROLAZA, M. (1996). Tridimensional Numerical Model for the Wind Energy Autonomy Evaluation Over Some Isolated Venezuela Coast Regions. Proceeding of the 1996 3rd International Congress on Numerical , 87-94.
Bossanyi, E. (2003). GH BLADED VERSION 3.6 USER MANUAL . Garrad Hassan asn Partners.
Bossanyi, E. A. (December de 2003). GH Bladed Theory Manual. Bristol, United Kingdom: Garrad Hassan and Partners.
Buhl, J. M. (29 de April de 2002). NWTC Design Ccodes (Modes) .
BUHL, J. M., WRIGHT, A. D., & PIERCE, K. G. (2001). The 20th ASME Wind Energy Symposium Reno, Nevada, January 8-11, 2001. NREL/CP-500-28848. FAST_AD Code Verificaction: A Comparison to ADAMS. Reno, Nevada: National Renewable Energy Laboratory.
Buhl, J. M., Wright, A., & Pierce, K. (8 - 11 de January de 2001). The 20th ASME Wind Energy Symposium, Reno, Nevada. FAST_AD Code Verification: A Comparison to ADAMS . Golden, Colorado, USA: National Renewable Energy Laboratory.
Buhl, Jr. M.L. (5 de April de 2002). Installing NWTC Design Codes on PCs Runnign Windows NT®. Golden, Colorado, USA: National Renewable Energy Laboratory.
Center, N. W. (18 de March de 2010). NWTC Design Codes. Recuperado el 2 de August de 2011, de http://wind.nrel.gov/designcodes/preprocessors/iecwind/
Crespo Martínez, A. (2006). Mecánica de Fluidos. Madrid: INDUSTRIALES Paraninfo.
ed., I. 6.-1. (2001). Wind Turbine Generator Systems.
Energía Renovables: ventajas y desentajas de la energía eólica. (2006). Recuperado el 2011, de http://www.revistafuturos.info/futuros14/energia_eolica.htm
F1 Revolution, una mirada diferente. (19 de February de 2009). Recuperado el 20 de August de 2011, de http://f1revolution.com/2009/02/capitulo-5-fundamentos-basicos-de-aerodinamica/
GASTÓN, M. (16 de April de 2008). Wind resources map of Spain at mesoscale. Methodology and validation. Recuperado el 2 de August de 2011, de http://secure.cener.com/documentos/wind-resources-map-mesoscale-PaperEwec08.pdf
45
Hewitt, P. G. (2004). Física Conceptual. En P. G. Hewitt, Física Conceptual (págs. 276 - 281). Addison Wesley.
IMAMURA, H., TSUMAMUNA, T., & KUROKAWA, J. (2004). Study on the Wind Measurements and Performance Evaluation of a WTGS in Complex Terrain. Japan Society of Mechanical Engineers, Volumen 70, Tokyo, Japan , 1230-1236.
Jason M. Jonkman, M. L. (August de 2005). FAST User’s Guide . 1617 Cole Boulevard, Golden, Colorado 80401-3393, Colorado, USA.
JONKMAN, B. J. (September de 2009). TurbSim User's Guide. Golden, Colorado: National Renewable Energy Laboratory.
Jonkman, B. J., & Buhl, J. M. (26 de May de 2005). TurbSim User's Guide. Recuperado el 11 de August de 2005, de http://wind.nrel.gov/designcodes/preprocessors/turbsim/turbsim.pdf
Jonkman, J. (December de 2003). Modeling of the UAE Wind Turbine for Refinement of FAST_AD . 1617 Cole Boulevard Golden, Colorado, USA.
Jonkman, J. M., & Buhl, J. M. (2004). New Development's for the NWTC's FAST Aerroelastic HAWT Simulator. The 23rd ASME Wind Energy Symposium, Reno, Nevada January 5-8, 2004 . Golden, Colorado, USA: National Renewable Energy Laboratory.
Kane, J. W., & Sterheim, M. M. (2007). Física. Barcelona: REVERTÉ.
KOSE, R., ARIF, M., & ERBUS, O. (2004). The Analysis of Wind Data and Wind Energy Potencial in Kutaya, Turkey. Renowable and Sustainable Energy Reviews, Volumen 8 , 277-288.
Laino, D. (12 de May de 2000). Recuperado el 9 de July de 2002, de hhtp://wind.nrel.gov/designcodes/preprocessors/windmaker/
Laino, D. (19 de June de 2001). Recuperado el 9 de July de 2011, de http://wind.nrel.gov/designcodes/preprocessors/iecwind/
LAINO, D. J., & CRAIG HANSEN, A. (06 de September de 2001). AeroDyn Interface for ADAMS®. USER'S GUIDE to the Computer Software Routines . Salta Lake City, Utah: National Renewable Energy Laboratory.
Laino, D. J., & Craig Hasen, A. (06 de Sep de 2001). USER'S GUIDE to the Computer Software Routines AeroDyn Interface for ADAMS®. Salt Lake City, Utah, USA: National Renewable Energy Laboratory - Windward Engineering.
LAINO, D. J., & HANSEN, A. C. (August de 2001). User's Guide to the Wind Turbine Dynamics Computer Software AeroDyn . Salt Lake City, Utah: Windward Engineering, LC.
Laino, David J.; Hansen, A. Craig. (24 de August de 2001). AERODYN USER'S GUIDE to the Wind Turbine Aerodynamics Computer Software . Salt Lake City, Utah, USA: Utah: Windward Engineering, LC.
Larsen, T. H. (December de 2007). How 2 HAWC2, the user's manual. Roskilde, Denmark: Risø National Laboratory, Technical University of Denmark.
Manjock, A. (26 de May de 2005). Evaluation Report: Design Codes FAST and ADAMS® for Load Calculations of Onshore Wind Turbines . Humburg, Germany: Germanischer Lloyd WindEnergie GmbH.
46
MANJOCK, A. (26 de May de 2005). Evaluation Report: Design Codes FAST and ADAMS® for Load Calculations of Onshore Wind Turbines. Report Nº. 72042 . Humburg, Germany: Germanischer Lloyd WindEnergie GmbH.
Moriarty, P. J., & Butterfield, S. B. (2009). Wind Turbine Modeling Overview for Control Engineers. 6.
Munteanu, I., Iuliana Bratcu, A., Cutululis, A., & Ceanga, E. (2008). Optimal Control of Wind Energy Systems. London: Springer-Verlag London Limited.
MURAKAMI, S., MOCHIDU, A., & KATO, S. (2003). Development of Local Area Wind Prediction System for Selecting Suitable Site for Windmill. Journal of Wind Engineering and Industrial Aerodynamics, Volumen 91, Elsevier Ltd. , 1759-1776.
OYARZO PEREZ, H. (2008). Modelización matemática para evaluar energías para sistemas eólicos e híbridos eólico-diesel , 19-20. Madrid.
Øye, S. (1991). Simulations of Loads on a Wind Turbine including Turbulence. Stockholm, Sweden: IEA Meeting.
POLLABAZER, R. (1995). Evaluation of Wind-Generator Potentiality. Solar Energy. Volumen 55, Pergamon Press Inc. Tarrytowm. NY. USA , 49-59.
Rodríguez Amenedo, J. L., Burgos Díaz, J. C., & Arnalte Gómez, S. (2003). Sistemas Eólicos de Producción de Energía Eléctrica. Madrid: Rueda S.L.
SÁNCHEZ, A., & VIEDMA, A. (2003). Energía Eólica. España: Editorial Horacio Escarabajal.
Sanchis Sabater, A. (1999). Fundamentos Físicos para Ingenieros. En A. Sanchis Sabater, Fundamentos Físicos para Ingenieros (págs. 368 - 370). Valencia: REPROVAL S. L.
TAYLOR, P., & WALMSLEY, J. (1986). Guidelines and Models for Estimating Wind Speed at WECS in Complex Terrain. New York, NY, USA: Pergamon Press.
van Dam, J., Link, M., Meadors, M., & Bianchi, J. (June de 2002). Wind Turbine Generator System Safety and Function Test. Golden, Colorado, USA: National Renewable Energy Laboratory.
Wilson, R. E., Walker, S. N., & Heh, P. (May de 1999). Technical and User's Manual for the FAST_AD Advanced Dynamics Code . Oregon, USA: Oregon State University.
WILSON, R. E., WALKER, S. N., & HEH, P. (May de 1999). Technical and User's Manual for the FAST_AD Advanced Dynamics Code. Oregon: OSU/NREL Report 99-01. Corvallis.
YANG, Z., & XUE, H. (2004). Study on the Estimation of Wind Energy over Complex Terrain. Science Press, Begin, 100085, China.Volumen 25 , 744-749.
47
ANEXO
48
PROCEDIMIENTO DE USO DEL PROGRAMA FAST
Este programa se puede ejecutar de dos formas: bajo entorno Windows o bajo el entorno Matlab, a continuación se procederá a explicar cada una de ellas.
ENTORNO WINDOWS.
Introducción
Instalar y ejecutar las aplicaciones de consola como es el código diseñado por National Wind Technology Center (NWTC) es muy diferente a las aplicaciones normales de Windows de 32 o 64 bits. Por lo general, el usuario tendrá que instalar los códigos de tal forma que se ejecuta desde una ventana de comandos y de cualquier carpeta. Este documento será una guía de la instalación de los códigos para el sistema operativo Windows 7®. En este tutorial el código que será instalado es FAST, pero estas instrucciones deben trabajar para cualquiera de los otros códigos que distribuye NWTC.
Inicialmente debe descargar el archivo de un código del servidor web http://wind.nrel.gov/designcodes/. El archivo almacenado debe tener un nombre como FAST_v70001a-bjj.exe. Crear una carpeta llamada FAST en algún lugar de la unidad C y poner los archivos entro de ella. Para este ejemplo, asumiremos que el nombre de la carpeta es C:\DesignCodes\Fast\. Puede hacer doble clic en el archivo desde el Explorador de Windows o ingresar FAST_v70001a-bjj en la ruta de archivos de la carpeta Fast como directorio actual. Esto creará algunos archivos y carpetas. Incluido en la carpeta Fast estará un archivo llamado Fast.exe. Este es el archivo ejecutable para FAST.
Configurar una carpeta Batch para los archivos del programa FAST.
Esto es algo que necesita hacer sólo una vez. Después de que haya hecho para este código, no será necesario repetir el proceso para otros códigos. Crear una carpeta en algún lugar de su sistema de archivos. Este tutorial se utiliza C:\Batch. Ahora, agregue la carpeta a la ruta de búsqueda, debe tener privilegios de administrador en el PC para hacer esto. En Windows 7, haga clic en Inicio, Panel de control, Sistema y Seguridad, Sistema, Configuración avanzada del sistema. En la pestaña Opciones Avanzadas dar un clic sobre Variables de entorno en la parte inferior derecha del cuadro de dialogo. Esto abrirá un cuadro de dialogo con dos ventanas una de variables de usuario y otra de variables del sistema. Seleccione Ruta de la lista Variables del sistema y haga clic en Editar.
La palabra Path se mostrará en el cuadro Nombre de variable y la ruta de búsqueda actual se situará por debajo en el cuadro Valor de variable. Escriba "C:\Batch," al inicio de la ruta en el cuadro Valor de variable. No olvide poner en el separador punto y coma. Haga clic en el botón Aceptar para cerrar todas las ventanas que usa para cambiar la ruta y poner el cambio en vigor.
El cambio que acaba de hacer será permanente y se aplicará a todas los nuevos comandos de Windows que usted abre. No es necesario reiniciar el PC. Las rutas de búsqueda para Windows que ya estaban abiertas cuando se hace el cambio no se actualizan, por lo que no se puede utilizar para ejecutar el código.
El siguiente paso es crear un archivo para traer el código. El uso de un editor de texto nos ayudará a crear un archivo llamado fast.bat en la carpeta C:\Batch. Sólo es necesaria una línea en el archivo fast.bat. Su contenido debe ser algo similar a esto:
C:\DesignCodes\Fast\Fast.exe% *
Esto le dice al sistema operativo para encontrar un archivo ejecutable llamado Fast.exe en la carpeta Fast. El %* indica al archivo para pasar de argumento en la línea de comandos para Fast.exe.
Invocando el Código
Para ejecutar un código, necesita inicialmente copiar los tres archivos que se muestran a continuación Test##.fst, Test##.cru y Test##_AD.ipt (##: Número de archivo a ser ejecutado, estos archivos se encuentran en la carpeta CertTest.) en la carpeta donde se encuentra Fast.exe
Se escribe la siguiente línea de código:
Fast Test##.fst
Recomendación: El momento de ejecutar Fast.exe con el archivo a analizar, si existe algún error el programa va indicando en pantalla que error se produce, principalmente los archivos adicionales que necesita para ejecutar, el momento que suceda esto, se debe buscar el archivo y copiarlo en la carpeta donde se encuentra FAST.exe, esto puede suceder varias veces, por eso es preferible crear una carpeta diferente donde se vaya copiando todos los archivos que necesita para realizar una correcta ejecución.
Cuando el archivo se ejecuta correctamente debe salir la siguiente ventana (figura 1) y dentro de la carpeta existe un archivo con nombre test##.out donde se encuentran valores numéricos para graficar en cualquier programa que permita su visualización como por ejemplo Excel.
Figura 20. Ejecución correcta
50
Entorno SIMULINK
Descripción General
SIMULINK es una herramienta de simulación popular para el diseño de controles que vende Mathworks, Inc., en un solo paquete con MATLAB.
Simulink tiene la capacidad de incorporar rutinas personalizadas Fortran en un bloque llamado S-function. Las subrutinas de FAST han sido unidas con una subrutina estándar de MATLAB con el fin de utilizar las ecuaciones de movimiento de FAST en un S-función que se puede incorporar en un modelo Simulink. Esto introduce una gran flexibilidad en los controles de turbina de viento implementados durante la simulación.
Los módulos control del generador de par, control de orientación de la góndola, y los módulos de control de inclinación pueden ser diseñados en el entorno de Simulink y simular mientras hacen uso de las ecuaciones completas de movimiento no lineales aeroelásticas de la turbina de viento disponible en FAST.
El bloque de la turbina de viento, como se muestra en la Figura 2, contiene el bloque S-Function con las ecuaciones de movimiento de FAST. También contiene los bloques que integran los grados de libertas (DOF por sus siglas en inglés) como el desplazamiento, velocidad y aceleración. Así, las ecuaciones de movimiento son formuladas en la S-Function de FAST, pero resuelve usando los solvers de Simulink.
Figura 21. S-Function bloque turbina de viento.
La interfaz entre FAST y Simulink es muy similar a la interfaz desarrollada para el código Dinámica Simbólica (SymDyn), que es una herramienta de control orientada a análisis HAWT desarrollado por investigadores del NREL. El modelo estructural de FAST, sin embargo, es de mayor fidelidad que la de SymDyn.
Primeros pasos
Para construir un modelo Simulink que usa turbina de viento dinámica en FAST en S-Function, usted debe comprar el software comercial MATLAB con el paquete adicional de Simulink. MATLAB está disponible en The Mathworks, Inc. (http://www.mathworks.com/). El conocimiento práctico de desarrollo de modelos de Simulink es también esencial.
El archivo FAST contiene varios archivos que son pertinentes a la interfaz FAST con Simulink, como se describe a continuación:
FAST_SFunc.mex32 La compilación S-Función de FAST tiene un vínculo dinámico librería (DLL) de Matlab llamada mexw32. Esta DLL contiene las rutinas dinámicas estructurales de FAST, las rutinas de aerodinámica de Aerodyn, y las interfaces con Simulink.
Simsetup.m Este archivo de MATLAB solicita al usuario el nombre del archivo a ser analizado por FAST y llama a Read_FAST_Input.m, que inicializa las variables del modelo. Debe ser llamado desde el espacio de trabajo MATLAB antes de ejecutar un modelo de Simulink con la S-Function de FAST.
Read_FAST_Input.m Este archivo de MATLAB es llamado por Simsetup.m y lee los archivos de entrada FAST para inicializar los parámetros de un modelo Simulink. Los usuarios no deberían modificar este archivo.
OpenLoop.mdl Un ejemplo de modelo Simulink que contiene el bloque de S-Function FAST, estos bloques que integran los grados de libertad (DOFs) y los bloques de entrada de lazo abierto constante.
Test01_SIG.mdl U ejemplo de un modelo Simulink que contiene el bloque S-Function de FAST, los bloques que integran los grados de libertad (DOFs), y el simple modelo de generador de inducción para FAST con certificación prueba # 01 implementado en Simulink.
Para ejecutar un modelo de FAST en Simulink se debe usar solo plataformas Windows y Matlab de 32 bits, ya que en 64 bits da errores de versión por la librería DLL.
1) Se crea una carpeta donde se realizará las simulaciones 2) Se copia los archivos Simsetup.m, Read_FAST_Input.m, FAST_SFunc.mex32 y
OpenLoop.mdl que se encuentran dentro de la carpeta FAST (lea apartado ENTORNO WINDOWS)
3) Copiar todos los archivos que se encuentras dentro de la carpeta CertTest en la carpeta generada para la simulación.
4) Se ejecuta Matlab, luego se busca en la ventana Directorio Actual (Current Directory) la carpeta donde se encuentran los archivos anteriormente copiados.
5) En la ventana de comandos (Command Window) escriba Simsetup y presione Entrar (Enter) como se muestra en la figura 3.
52
Figura 22. Ejecución Simsetup
6) El siguiente paso es escribir el archivo a ser analizado Test##.fst 7) Si es correcta la ejecución tendrá una ventana como se muestra en la figura 4.
Figura 23. Ejecución archivo Test##.fst
8) Abrimos el modelo OpenLoop.mdl, este modelo debe ser como la Figura 5.
Figura 24. Modelo Simulink OpenLoop.mdl
53
9) Por último, haga clic en el Play (►) en la ventana de Simulink para ejecutar la simulación. Si la simulación es correcta debe salir una ventana semejante a la figura 6.
Figura 25. Ejecución correcta del modelo en Simulink
La S-Function de FAST generará los mismos ficheros ASCII de salida como se genera durante una simulación rápida normal. Estos archivos de salida usar el nombre de la raíz del archivo de entrada principal pero con extensión OUT. Por ejemplo, si el archivo de entrada principal fue nombrado Test01t.fst, el archivo de salida principal de salida Test01_SFunc.out mientras que el ejecutable de FAST generaría Test01.out. Se puede realizar cambios de parámetros iniciales de ejecución de Simulink, cambios estructurales del aerogenerador al abrir los tres archivos principales y buscar en cada uno de ellos los valores adecuados del aerogenerador a utilizar, FAST tiene varios archivos para su análisis (como se indica en el manual de usuario del programa FAST). Si el usuario es inexperto se recomienda no abrir esto tres archivos antes indicados o cambiar sus parámetros hasta que se familiarice completamente con este programa.