CARGAS SOBRE EL ROTOR DE UNA TURBINA
1. TIPOS DE CARGAS
El rotor de un aerogenerador de eje horizontal esta expuesto a numerosas cargas que se clasifican
segn su forma de actuacin. En este tema se revisarn las distintas clases de cargas sobre un
aerogenerador que se clasificarn segn su origen y segn su variacin en el tiempo. Las
convenciones y la nomenclatura utilizada en este documento se pueden apreciar en el ejemplo de un
aerogenerador bipala a sotavento que muestra la figura 1.
Ilustracin 1:Solicitaciones
Para generar la suficiente energa a lo largo de su vida til, se debe sobrevivir las adversidades a las
que va a estar sometido. Los prototipos anteriores a lo largo de los aos no pudieron cumplir con este
simple pero necesario objetivo, ya que no se tena en cuenta por diversas razones. A lo largo de los
aos se ha empezado a tener ms en cuenta para alargar la vida del rotor del aerogenerador.
Si consideramos el aerogenerador como una construccin, no debemos pensar en ste como un
sistema cerrado que nicamente esta cargado exteriormente por el viento, principalmente por el efecto
aerodinmico sobre las palas en rotacin. No debemos olvidar la fuerza de gravedad, que cargar las
partes msicas del sistema. Tambin aparecen fuerzas y momentos inerciales ciando las partes del
aerogenerador, principalmente las palas, se aceleran bajo influencia de las cargas aerodinmicas y
msicas antes descritas. La rotacin de las palas genera un cargas inerciales muy importantes: las
cargas centrfugas. Finalmente hay que considerar que todas las fuerzas y momentos deben estar en
equilibrio con las fuerzas reactivas internas de la estructura. En resumen, las cargas en los
aerogeneradores se dividen segn su origen en las siguientes categoras.
cargas aerodinmicas (aerodymamic loading)
cargas gravitatorias (gravity loading)
cargas inerciales (inertial Loads)
cargas estructurales (structural loads)
Como se puede ver, en ingles se diferencia, siguiendo una convencin, entre loads y loading donde
loading denota una accin externa sobre la estructura, mientras que load se refiere a la reaccin
interna inercial y elstica. El equilibrio final de las cargas se realiza en dos contornos del sistema:
En el generador
y a nivel de suelo en la cimentacin
1.1 DISEO GENERAL
Generalmente el diseo es un proceso cclico dividido en tres fases:
Concepto
Prediseo
Diseo
CONCEPTO
En la fase de concepto se debern elegir unas opciones de diseo y las dimensiones principales
debern ser decididas. Para un aerogenerador conectado a red con un rgimen de viento especificado
las opciones principales son:
1. nmero de palas
2. radio de pala
3. perfil de la pala
4. distribucin de cuerda y de torsin de la pala
5. potencia nominal
6. altura del buje
7. barlovento o sotavento
8. control de potencia por perdida aerodinmica o por control de paso
9. control de potencia activo o pasivo
10. buje fijo, teeter o buje articulado
11. vueltas fijas o variables
12. torre rgida o flexible
13. multiplicadora o conexin directa al generador
PREDISEO
En la fase de prediseo se realizan los clculos de cargas estticas, que aparecen como niveles de
tensiones en un nmero limitado de puntos, que normalmente suelen ser los que mayores cargas
sufren. Asimismo, la viabilidad econmica del concepto elegido es juzgada mediante la estimacin del
coste de produccin de una unidad de energa (kWh) . El resultado del prediseo depender de si el
concepto elegido es econmica y a su vez tcnicamente viable. En el caso de no ser as es necesario
volver a tras y realizar un concepto ms realista.
DISEO
En la siguiente fase se genera el diseo real en varios ciclos. En cada ciclo se calcula la media y las
fluctuaciones de las cargas, as como las deflexiones, incluyendo la parte dinmica en el modelo de la
turbina. La resistencia esttica y la resistencia a fatiga se determina para un gran nmero de puntos de
la construccin. El uso de un buen cdigo aeroelstico computerizado es indispensable en esta fase de
diseo.
Ms adelante se identificarn y se cuantificarn las cargas principales de una aeroturbina, asumiendo
algunas simplificaciones.
1.1.1 La variacin en el tiempo de las cargas
Cada tipo de carga puede consistir de una parte estacionaria y una parte variable, generando cargas
estticas y dinmicas en las diferentes parte de la aeroturbina. Las cargas dinmicas se subdividen en
deterministas y estocsticas. Una carga determinista es peridica en el tiempo, de tal manera que si se
obtienen las cargas para un periodo son predecibles ms adelante. Las cargas estocsticas no son ni
predecibles ni peridicas. Si estn dadas para un cierto rango de tiempo, solamente se podrn predecir
cantidades estadsticas como son la media, la desviacin tpica, la varianza y la densidad espectral
para un futuro rango de tiempo.
Ilustracin 2: Tringulo de velocidades
CARGAS SOBRE EL ROTOR DE UNA TURBINA
2. CARGAS AERODINMICAS SOBRE LA PALA
Las cargas aerodinmicas sobre la pala son aquellas causadas por la corriente del aire y su interaccin
con las palas en rotacin del aerogenerador. El tringulo de velocidades para una seccin de pala y la
nomenclatura usada en esta seccin se pueden apreciar en la figura 2. Ut es la corriente debida a la
rotacin de las palas, incluyendo el coeficiente de velocidad inducida angular. Up es la velocidad del
viento perpendicular al plano del rotor, tambin corregido con el coeficiente de velocidad inducida axial.
La carga ms importante que sufre la pala es la fuerza de empuje Fax, apuntando en la direccin
perpendicular al plano del rotor.
(1)
(2)
(3)
(4)
Fax= fuerza de empuje
r = posicin del elemento de pala a lo largo del eje de la pala
= densidad del aire
U = velocidad de viento resultante en el elemento de pala
c = longitud de cuerda
Cl = coeficiente de sustentacin
Cd= coeficiente de resistencia
= ngulo de ataque
Ut = velocidad de viento tangencial al plano del rotor
Up = velocidad de viento perpendicular al plano del rotor
a = coeficiente de velocidad inducida axial
a = coeficiente de velocidad inducida angular
wr = velocidad de rotacin
V = velocidad de viento aguas arriba
La fuerza total de empuje Fax en direccin del eje del rotor se puede obtener integrando a lo largo de la
pala y sumndola para cada una de las palas.
2.1 EL EMPUJE ESTTICO FRENTE A LA VELOCIDAD DE VIENTO
La relacin entre la fuerza de empuje y la velocidad de viento a la altura del buje depende del sistema
de limitacin de potencia. Los dos sistemas ms comunes de regulacin de potencia son:
* control pasivo por perdida aerodinmica manteniendo la velocidad de rotacin constante en todo el
rango operacional de velocidades de viento
* control activo por cambio de paso a una velocidad de viento igual o mayor a la velocidad de viento
nominal. La velocidad de viento nominal es la velocidad de viento a la que el generador alcanza la
potencia nominal de diseo
Esto queda lastrado en las medidas tomadas de la fuerza axial como funcin de la velocidad de viento
en tres casos diferentes:
* En el rotor 20 WPX de perdida aerodinmica, (dimetro 21 .5m) medida a velocidad de giro constante
de 47 r.p.m.
* En el mismo rotor pero ahora equipado con un control pasivo de cambio de paso de punta de pala
(punta de pala de 2.5 m) y operando a velocidades de rotacin variables
* En el NEWECS 45 con control activo por cambio de paso a velocidad de giro variable En el caso del
control por perdida aerodinmica (20 WPX) Fax aumenta de forma montona en todo el rango de
velocidades de viento. Para los aerogeneradores con control por cambio de paso (20 WPX-PRT1 y
NEWECS45) la fuerza de empuje alcanza su mximo en la velocidad de viento nominal.
La fuerza de empuje en cada pala genera en la raz de pala un momento flector (fuera del plano del
rotor o momento flector en flap) . La parte esttica puede ser compensada por la fuerza centrifuga , si
se introduce conicidad (vase seccin 2.4). La parte fluctuante junto al momento flector en el plano del
rotor causado principalmente por la masa de las palas en rotacin, determinan las tensiones fluctuantes
en la raz de pala. La dimensin de la raz de pala est determinada por estas cargas.
2.2 FLUCTUACIN DEL EMPUJE
Para poder entender la parte fluctuante del empuje la ecuacin (1) es aproximada por
(5)
asumiendo lo siguiente
* = wr/ V > 5, de tal manera que el perfil opere antes de entrar en prdida, con wr velocidad de
rotacin y V = velocidad de viento aguas arriba
* la componente en el plano del rotor de la velocidad de viento Up>> Ut
* cos =1
* Cl>>Cd
Para un aerogenerador operando a velocidad de giro constante, Ut es constante ( a parte de cambios
en el coeficiente de velocidad inducida angular) y si se excluyen cambios de ngulo de paso (operacin
antes de llegar a potencia nominal), la variacin del ngulo de ataque es
(6)
Entonces la variacin en Cl es
(7)
Una aproximacin de la variacin en el empuje axial por unidad de longitud de palas se obtiene
combinando las ecuaciones (5) y (7)
(8)
Para la fuerza tangencial (la fuerza en el plano de rotacin) es valida una expresin similar a la
ecuacin (1) (comparar con figura 2):
(9)
Aplicando las mismas aproximaciones que en ecuacin (1) nos lleva a
(10)
Asumiendo funcionamiento antes de entrada en perdida, aerodinmica lineal, velocidad constante y un
alto, las fuerzas tanto axiales como tangenciales varan linealmente con las variaciones en la
velocidad de viento de axial.
2.2.1 Sombra de torre
La estructura de soporte de la aeroturbina influencia la corriente que pasa a travs del rotor, tanto a
sotavento como a barlovento para una torre tubular. Esta carga aerodinmica determinista tiene
carcter impulsivo si la frecuencia es igual a la velocidad de rotacin. En la literatura la frecuencia del
rotor se conoce como 1 P. Una descomposicin de Fourier en el dominio de la frecuencia de las cargas
de la pala nos muestran picos de intensidad en 1 P as como en los armnicos ms altos 2P, 3P, 4P
debido a el caracter impulsivo. Con acuerdo a las frmulas (9) y (10) esto tambin producira
componentes armnicos y mltiplos en la fuerza de empuje y tangencial.
2.2.2 Cortadura del viento
La velocidad de viento aumenta con la altura sobre el nivel del mar. As cada vez que una pala se
encuentra en la parte superior del rotor las cargas son ms altas y cuando est en la parte inferior del
rotor las cargas son ms bajas. Despreciando los trminos de alto orden, el perfil de cortadura de
viento a la altura de buje se puede describir como:
(11)
z posicin vertical de la punta de pala con H altura del buje U(H) velocidad del viento a la altura H z0
rugosidad del terreno Por ejemplo, para un aerogenerador con H = 30 m , R - 10 m, situado en terreno
llano con z0=0.05 m la cortadura del viento variar a lo largo del dimetro de la pala un 10.4% (2 x
0.156 x 0.333 = 0.104) de la velocidad de viento a la altura del buje. Siguiendo la ecuacin (8), la
variacin de velocidad de viento genera una variacin en la fuerza de empuje. La variacin relativa de
la fuerza de empuje se consigue dividiendo ecuacin (8) entre ecuacin (5):
(12)
Para un perfil aerodinmico NACA 4415 con = 4 grados, Cl = 0.8, dCl/d= 5.73 rad-1 y =Ut / Up =
6, la variacin en la carga es un coseno puro en 1P, si consideramos aerodinmica lineal, as como una
aproximacin lineal de la cortadura del viento. En realidad tambin se generan armnicos ms altos,
aunque con intensidades comparadas a 1 P mucho ms pequeas que en el caso de las cargas por
sombra de torre.
2.2.3 Desalineacin del rotor
Un aerogenerador debe operar con el eje del rotor alineado a la direccin del viento . En la prctica
esta alineacin nunca es perfecta, especialmente cuando los cambios en la direccin del viento son
muy rpidos y los momentos inerciales en el eje de orientacin son grandes. Esta desalineacin genera
una modulacin sinusoidal de la velocidad relativa en el plano del rotor:
(13)
y la variacin del empuje ahora es originada por una variacin de Ut. No debemos olvidar que Cl
depende de Ut a travs del ngulo de ataque, de tal manera que
(14)
Asumiendo las mismas simplificaciones que arriba, la variacin en Cl se escribe como:
(15)
Combinando ecuacin (14) y (15), la variacin aproximada del empuje se convierte en:
(16)
As, en el caso de desalineacin del rotor, la variacin del empuje contiene dos trminos , el primero
directamente de la variacin de Ui, y el segundo a travs de la variacin de Cl. La relativa importancia
del segundo trmino depende de las condiciones aerodinmicas: Para un perfil aerodinmico NACA
4415 con = 4 grados, Cl = 0.8, dCl/d = 5.73 rad-1 y = Ut / Up = 6 el segundo trmino es 0.73
veces el primero.
2.2.4 ngulo de tilt
En la mayora de aerogeneradores el eje del rotor forma un pequeo ngulo (> 5 grados) con el plano
horizontal (el ngulo de tilt), para aumentar la separacin entre palas y torre y evitar su interaccin. Con
respecto a las cargas aerodinmicas el tilt es equivalente a la desalineacin del rotor, dado que el
plano del rotor no esta perpendicular a la velocidad del viento V.
2.3 CARGAS GRAVITACIONALES
Para un aerogenerador sin ngulo de tilt y sin conicidad, la gravedad genera un momento flector en la
raz de pala con una frecuencia 1P y con una amplitud MgZ [Nm], con M como masa de la pala, Z
como centro de gravedad de la pala y g como aceleracin debido a la gravedad. Especialmente en los
grandes aerogeneradores esta es la carga dominante en el plano del rotor. La introduccin del ngulo
de tilt causa un momento flector de flap adicional constante en la raz de pala proporcional al seno del
ngulo de tilt. La conicidad aade un momento flector de flap variable adicional con una frecuencia 1p,
proporcional al seno del ngulo de conicidad. Para un ngulo de conicidad de 5 grados el momento
flector de flap se aumenta en un 8.7% de las cargas gravitatorias en el plano de rotor.
Debido a la menor resistencia frente a la flexin en flap que a la flexin en lag, la tensin resultante en
una seccin de la pala debido a las cargas gravitacionales en un aerogenerador con conicidad no debe
de ser menospreciada.
2.4 CARGAS INERCIALES
La principal carga inercial en un aerogenerador es la fuerza centrifuga que generan las palas en
rotacin. Algunas veces se aplica un ngulo de conicidad para compensar la parte esttica del empuje
aerodinmico. Aunque la compensacin del empuje a travs de la conicidad tiene sus desventajas:
* las cargas fluctuantes alrededor del nivel cero de tensin con cambios de signo son muchas veces
ms dainas que las fluctuaciones alrededor de una media sin cambios de signo ver, ms adelante, las
propiedades a fatiga de los materiales compuestos)
* en el caso de embalamiento, la sobrecompensacion agrava esta condicin de fallo
Los momentos giroscpicos ocurren cuando cambia de direccin la orientacin del vector de velocidad
angular. La orientacin al viento de un aerogenerador en operacin es una situacin as. Aunque
cuando se orienta activamente como es el caso de la mayora de aerogeneradores a barlovento, las
velocidades de orientacin son tan bajas (1 a 3 grados/s) que el momento de par giroscpico resultante
es despreciable. Los aerogeneradores a sotavento suelen orientarse libremente. En ese caso las
velocidades de orientacin son mucho mayores (10 a 20 grados /s) y el momento giroscpico,
orientado perpendicularmente al eje del rotor, debe ser tenido en cuenta.
2.5 TURBULENCIA
En la discusin sobre cargas aerodinmicas, el viento era considerado como estacionario. En realidad
el viento es turbulento como se ha descrito en la parte de recursos (ver captulo referente a recursos
elicos). El efecto en las cargas aerodinmicas es doble:
* Un contribucin directa. Las fluctuaciones estocsticas son superpuestas a las cargas aerodinmicas
directamente a travs de las fluctuaciones estocsticas en la velocidad de viento relativo Up (axial)
(comparar figura 2 y ecuacin (5) para la fuerza de empuje y la aerodinmica lineal). Las frecuencias
de estas fluctuaciones son bajas, es decir , menores de 0.2P (1 P = frecuencia de rotacin)
* "Rotational sampling": Cuando una pala rota a travs de un campo de viento turbulento esta expuesto
a un espectro de intensidad de viento diferente al de un observador estacionado a la altura del buje. La
naturaleza de esta transformacin provoca que las fluctuaciones de baja frecuencia sean transportadas
a frecuencias ms altas alrededor de 1 P, 2P, 3P, 4P de tal manera que la varianza y el contenido
energtico de la seal quedan inalterados.
La importancia relativa de la fluctuaciones estocsticas con respecto a las fluctuaciones determinantes
aumenta con el dimetro del rotor y con la intensidad de turbulencia en el emplazamiento del
aerogenerador. La intensidad de turbulencia 1. esta definida por:
(17)
con v = desviacin estndar de la velocidad del viento
Vave= velocidad media del viento
Para un aerogenerador de 25 m de dimetro y una intensidad de turbulencia de 10%, la varianza ( el
cuadrado de la desviacin estndar) del momento flector en flap en la raz de pala es alrededor del
50% estocstico y del 50% determinista.
2.6 RESPUESTA DINMICA
Hasta ahora hemos discutido las cargas sobre el aerogenerador, olvidando la respuesta dinmica de la
turbina a las cargas. En principio el aerogenerador completo (palas, buje, multiplicadora, generador,
torre y cimentacin) define la respuesta dinmica a las cargas. Como primer paso de diseo
consideraremos las palas y el rotor separado de la torre, una aproximacin que solo es vlida para una
torre que sea mucho ms rgida comparada con las palas. Concretamente, nos deberamos restringir a
las frecuencias de vibracin en flap (fuera del plano del rotor) de las palas individualmente y de las
frecuencias de torre en direccin axial.
2.6.1 Frecuencias propias
Cuando un sistema, como el de un aerogenerador es excitado con cargas fluctuantes , la pregunta
clave que hay que hacerse es: donde estn las frecuencias propias ms bajas ? Una complicacin
posterior reside en que las frecuencias propias de las palas dependen de la velocidad de rotacin. Este
fenmeno se llama rigidez centrfuga ("centrifugal stiffness"). En la figura 10 se puede ver un elemento
de pala, que es empujado fuera del plano de rotacin. La fuerza centrfuga intenta devolver a su sitio al
elemento de pata, de tal manera que se genera una rigidez adicional hacia la flexin fuera del plano,
dependiendo de la velocidad de rotacin. El resultado es que las frecuencias propias en flap no son
constantes sino dependientes de la velocidad de rotacin.
Sin olvidar que las cargas fluctuantes externas estn concentradas alrededor de las frecuencias OP, 1
P, 2P, 3P ,..*. es muy importante que las frecuencias propias del sistema no coincidan (20 %) con los
mltiplos de la velocidad de rotacin. Un control (ver figura 11) en forma de diagrama de resonancia o
tambin conocido como diagrama de Campbell, es indispensable en todas las etapas del diseo. En el
diagrama las frecuencias propias ms importantes del sistema estn representadas frente a la
velocidad de rotacin junto a las lneas 1P, 2P, 3P,... de cargas fluctuantes. En la interseccin entre un
lnea nP y una frecuencia propia se produce resonancia para esa velocidad de rotacin en particular. El
diseador evitar a toda costa esta condicin de operacin de fallo.
2.6.2 Sistemas rgidos frente a sistemas flexibles
Un sistema se llama rgido ("stiff"), cuando las frecuencias de una fuerza fluctuante de entrada son
apreciablemente ms pequeas que las frecuencias propias ms bajas del sistema. Un sistema se
llama flexible ("soft"), cuando las frecuencias de entrada estn bastante por encima de las frecuencias
propias ms bajas del sistema. Consideremos ahora por ejemplo, la flexin de la torre causada por la
sombra de forre. Cada pala es cargada cuando pasa por la torre, pero la torre misma no distingue
individualmente entre las palas y as carga a las frecuencias: BP, 2BP, 3BP....con B como nmero de
palas. As se considerar una torre flexible ("soft") cuando la frecuencia con la que pasan las palas
(BP) es ms alta que la frecuencia propia de flexin de la torre. Algunas veces tambin se considera
una torre "soft-soft" en el diseo, es decir, cuando tambin 1P es ms alto que la frecuencia propia de
flexin de la torre. Las cargas de la torre a frecuencia 1 P son bajas y son originadas por asimetra
entre palas y por interaccin entre palas y torre. Para que las palas individualmente sean consideradas
como "soft" a la flexin en flap es necesario que 1P sea ms alto que la primera frecuencia propia. En
la prctica pocas palas "soft" son fabricadas. Para aerogeneradores bipala se utiliza frecuentemente un
buje articulado "teeter", es decirlas palas se conectan al buje a travs de una articulacin perpendicular
al eje principal.
